Modelagem molecular no estudo das interações receptor-ligante e no desenho racional de inibidores da biossíntese de petrobactina em Bacillus Anthracis deidroshikimato desidratase como alvo de novas terapias anti-antraz

Abstract

O antraz é uma doença infecciosa aguda grave, com uma taxa de mortalidade superior a 90% em sua forma respiratória, causada pelo Bacillus anthracis, uma bactéria altamente virulenta, que está desenvolvendo resistência e que tem potencial aplicação como arma biológica e agente de bioterrorismo. Nesse trabalho, a inibição de deidroshikimato desidratase do B. anthracis foi estuda por meio docking, dinâmica molecular e ensemble docking. Essa enzima é responsável por uma etapa chave na biossíntese de petrobactina, molécula através da qual o B. anthracis adquire ferro – micronutriente essencial para seu desenvolvimento e proliferação no hospedeiro. O docking de 25 compostos com ação inibitória conhecida na estrutura cristalográfica da enzima indicou interações importantes com os resíduos His144, His175, Phe211, Tyr217 (ligações de hidrogênio), Arg102 (ponte salina), His144 e Phe255 (interações π-π). Ligantes estruturalmente semelhantes ao cristalográfico (3,4-DHBA) foram docados adequadamente no sítio ativo, enquanto ligantes mais volumosos foram docados na entrada do sítio, resultando em baixa correlação entre as energias livres de ligação experimentais e os escores de docking (R² = 0,1295; R-Pearson = 0,360) e desvios de 23%, em média, frente ao experimental. Simulações de dinâmica molecular mostraram que essa proteína apresenta uma grande rigidez estrutural intrínseca, porém porções do seu sítio ativo, sobretudo da estrutura em forma de laço que o recobre, apresentaram flexibilidade significativa. A presença de ligantes induz a alterações conformacionais que proporcionam o alargamento do sítio e permitem a entrada de ligantes mais volumosos, indicando que o sítio cristalográfico era, de fato, muito restrito. A atividade inibitória aparenta estar relacionada com a formação de uma rede de ligações de hidrogênio entre os ligantes e resíduos do sítio ativo, sendo as principais entre grupos 3-OH do anel aromático dos ligantes e a His175; entre o grupo carboxílico e a Arg102 (ponte salina); entre o grupo 4-OH e a Phe211 e principalmente entre o grupo 5-OH e a His144, um resíduo importante no mecanismo enzimático. O ensemble docking em três estruturas extraídas das simulações de dinâmica molecular permitiu a aprimorar a correlação entre os escores de docking e atividade inibitória experimental, com R² = 0,363 e R-Pearson = 0,602 considerando a totalidade dos ligantes ou com R² = 0,8157 e R-Pearson = 0,903 considerando-se os dez ligantes mais potentes (contra R² = 0,5683 e R-Pearson = 0,754 na estrutura cristalográfica), evidenciando a necessidade de se considerar a flexibilidade do receptor para o docking adequado. Esse modelo linear juntamente com essa compreensão mais profunda dos mecanismos relacionados com a inibição dessa enzima permitirão o desenho e a triagem in silico de novas moléculas com potência e seletividade aprimoradas e com potencial aplicação como uma nova terapia contra o Bacillus anthracis.

Anthrax is a serious acute infectious disease with a mortality rate higher than 90% in its inhalational form. This disease is caused by Bacillus anthracis, a highly virulent bacterium that is developing resistance and which has potential application as a biological weapon and bioterrorism agent. In this work, the inhibition of dehydroshikimate dehydratase from B. anthracis was studied through docking, molecular dynamics and ensemble docking. This enzyme is responsible for a key step in the biosynthetic pathway of petrobactin, a molecule released by B. anthracis to acquire iron, an essential micronutrient for its development and proliferation within the host. Molecular dockings of 25 compounds with known inhibitory activity against dehydroshikimate dehydratase in the crystallographic structure of this enzyme indicated important interactions with the residues His144, His175, Phe211, Tyr217 (hydrogen bonds), Arg102 (salt bridge), His144 and Phe255 (π-π interactions). Ligands structurally similar to the crystallographic (3,4-DHBA) were appropriately docked within the active site, while bulkier ligands were docked at the site's entrance, resulting in a low correlation between the experimental binding free energies and the docking scores (R² = 0,1295; R-Pearson = 0,360), as well as a deviation of 23%, on average, compared to the experimental data. Molecular dynamics simulations showed that this protein has a high structural rigidity, however portions of its active site, especially the loop-like structure that covers it, showed a significant mobility. The presence of ligands induced conformational changes that lead to the widening of the site and allowed bulkier ligands to enter it, what indicates the crystallographic site was, in fact, very restricted. The inhibitory activity appears to be related with the formation of a network of hydrogen bonds between ligands and active site residues, mainly between the 3-OH moiety in the aromatic ring of ligands and His175; between the carboxylic group and Arg102 (salt bridge); between the 4-OH moiety and Phe211 and specially between the 5-OH group and His144, a residue with an important role in the enzymatic mechanism. Ensemble docking with three structures extracted from molecular dynamics simulations allowed to improve the correlation between docking scores and experimental inhibitory activity, with R² = 0,363 and R-Pearson = 0,602, when considering all ligands, and R² = 0,8157 and R-Pearson = 0,903 when considering the ten ligands of higher activity (against the values of R² = 0,5683 and R-Pearson = 0,754 for their docking in the crystallographic structure). This point out the need to account for receptor's flexibility for an appropriate docking. This linear model coupled with this deeper understanding about the mechanisms related with enzymatic inhibition will allow the in silico drug design and screening of new molecules with improved potency and selectivity and with potential application as a new therapy against Bacillus anthracis.