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

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.