Dynamics and Energetics of a Bromodomain in complex with bromosporine from Leishmania donovani

ABSTRACT

Leishmaniasis continues to be a neglected tropical disease, affecting people and animals and causing significant economic losses. Therefore, there is interest in the study and evaluation of new drug targets. In fact, it has been shown that by interfering with lysine-reading proteins such as bromodomain (BMD) there is a decrease in parasite survival. In this study, we researched the dynamics and energetics of the Leishmania donovani BMD in complex with bromosporin, which is considered to be a pan-inhibitor of BMDs, with the aim of understanding the molecular recognition mechanism. Molecular dynamics (MD) and non-equilibrium free energy calculation guided by steered molecular dynamics (SMD) simulations showed that the BMD has three flexible amino acid regions and bromosporin exhibiting various recognition states during the interaction. These results corroborate the promiscuity of bromosporin for energetically favourable sites, with the possibility of expanding its inhibition to other bromodomains. Furthermore, these results suggest that Van der Waals interactions have more relevance for complex recognition and residues ASN-87 and TRP-93 are key in forming hydrophobic and H-bond interactions, respectively. This research provides new insights for understanding the recognition mechanism, dynamics and energetics of the complex for the development of new therapeutic strategies.

RESUMEN

La leishmaniasis sigue siendo una enfermedad tropical desatendida, que afecta a personas y animales y causa importantes pérdidas económicas. De ahí el interés por estudiar y evaluar nuevas dianas farmacológicas. De hecho, se ha demostrado que al interferir con proteínas lectoras de lisina como el bromodominio ("bromodomain", BMD) se produce una disminución de la supervivencia del parásito. En este artículo estudiamos la dinámica y la energética del BMD de Leishmania donovani en complejo con bromosporina, que se considera un pan-inhibidor de BMD, con el objetivo de comprender el mecanismo de reconocimiento molecular. Las simulaciones de dinámica molecular (DM) y el cálculo de energía libre de no-equilibrio guiado por dinámica molecular de estiramiento (DMS) mostraron que BMD tiene tres regiones de aminoácidos flexibles y la bromosporina presenta varios estados de reconocimiento durante la interacción. Estos resultados corroboran la promiscuidad de la bromosporina por sitios energéticamente favorables, siendo posible expandir su inhibición a otros bromodominios. Además, los resultados sugieren que las interacciones de Van der Waals tienen más relevancia para el reconocimiento del complejo y los residuos ASN-87 y TRP-93 son clave en la formación de interacciones hidrofóbicas y de puentes de hidrógeno, respectivamente. Esta investigación proporciona nuevos conocimientos para comprender el mecanismo de reconocimiento molecular, la dinámica y la energética del complejo para el desarrollo de nuevas estrategias terapéuticas.

RESUMO

A leishmaniose continua a ser uma doença tropical negligenciada, que afeta os seres humanos e os animais e causa perdas econômicas significativas. Daí o interesse em estudar e avaliar novos alvos de medicamentos. De fato, a interferência com proteínas leitoras de lisina, como o bromo domínio ("bromodomain", BMD), tem demonstrado diminuir a sobrevivência do parasita. Neste trabalho, estudamos a dinâmica e a energética do BMD de Leish-mania donovani em complexo com a bro-mosporina, considerada um pan-inibidor da BMDs, com o objetivo de compreender o mecanismo de reconhecimento molecular. As simulações de dinâmica molecular (MD) e cálculo de energia livre de não-equilíbrio guiada por dinâmica molecular esticamento (MDS) mostraram que o BMD tem três regiões de aminoácidos flexíveis e que a bromosporina apresenta vários estados de reconhecimento durante a interação. Esses resultados corroboram a promiscuidade da bromosporina para locais energeticamente favoráveis, possibilitando a expansão de sua inibição para outros bromodomínios. Além disso, os resultados sugerem que as interações de Van der Waals são mais relevantes no momento do reconhecimento do complexo e os resíduos ASN-87 e TRP-93 são fundamentais na formação de interações hidrofóbicas e de ligações de hidrogênio, respectivamente. Essa pesquisa fornece novos insights para compreender o mecanismo de reconhecimento, a dinâmica e a energética do complexo para o desenvolvimento de novas estratégias terapêuticas.