Molecular mechanism of chorismate mutase activity of promiscuos MbtI

Best Practices on QM/MM Simulations of Biological Systems

During the second half of the 20th century, following structural biology hallmark works on DNA and proteins, biochemists shifted their questions from "what does this molecule look like?" to "how does this process work?". Prompted by the theoretical and practical developments in computational chemistry, this led to the emergence of biomolecular simulations and, along with the 2013 Nobel Prize in Chemistry, to the development of hybrid QM/MM methods. QM/MM methods are necessary whenever the problem we want to address involves chemical reactivity and/or a change in the system's electronic structure, with archetypal examples being the studies of an enzyme's reaction mechanism and a metalloprotein's active site. In the last decades QM/MM methods have seen an increasing adoption driven by their incorporation in widely used biomolecular simulation software. However, properly setting up a QM/MM simulation is not an easy task, and several issues need to be properly addressed to obtain meaningful results. In the present work, we describe both the theoretical concepts and practical issues that need to be considered when performing QM/MM simulations. We start with a brief historical perspective on the development of these methods and describe when and why QM/MM methods are mandatory. Then we show how to properly select and analyze the performance of the QM level of theory, the QM system size, and the position and type of the boundaries. We show the relevance of performing prior QM model system (or QM cluster) calculations in a vacuum and how to use the corresponding results to adequately calibrate those derived from QM/MM. We also discuss how to prepare the starting structure and how to select an adequate simulation strategy, including those based on geometry optimizations as well as free energy methods. In particular, we focus on the determination of free energy profiles using multiple steered molecular dynamics (MSMD) combined with Jarzynski's equation. Finally, we describe the results for two illustrative and complementary examples: the reaction performed by chorismate mutase and the study of ligand binding to hemoglobins. Overall, we provide many practical recommendations (or shortcuts) together with important conceptualizations that we hope will encourage more and more researchers to incorporate QM/MM studies into their research projects.

Shedding X-ray Light on the Role of Magnesium in the Activity of Mycobacterium tuberculosis Salicylate Synthase (MbtI) for Drug Design

The Mg²⁺-dependent Mycobacterium tuberculosis salicylate synthase (MbtI) is a key enzyme involved in the biosynthesis of siderophores. Because iron is essential for the survival and pathogenicity of the microorganism, this protein constitutes an attractive target for antitubercular therapy, also considering the absence of homologous enzymes in mammals. An extension of the structure–activity relationships of our furan-based candidates allowed us to disclose the most potent competitive inhibitor known to date (10, K_i = 4 μM), which also proved effective on mycobacterial cultures. By structural studies, we characterized its unexpected Mg²⁺-independent binding mode. We also investigated the role of the Mg²⁺ cofactor in catalysis, analyzing the first crystal structure of the MbtI–Mg²⁺–salicylate ternary complex. Overall, these results pave the way for the development of novel antituberculars through the rational design of improved MbtI inhibitors.

Tuning the Phosphoryl Donor Specificity of Dihydroxyacetone Kinase from ATP to Inorganic Polyphosphate. An Insight from Computational Studies

Dihydroxyacetone (DHA) kinase from Citrobacter freundii provides an easy entry for the preparation of DHA phosphate; a very important C3 building block in nature. To modify the phosphoryl donor specificity of this enzyme from ATP to inorganic polyphosphate (poly-P); a directed evolution program has been initiated. In the first cycle of evolution, the native enzyme was subjected to one round of error-prone PCR (EP-PCR) followed directly (without selection) by a round of DNA shuffling. Although the wild-type DHAK did not show activity with poly-P, after screening, sixteen mutant clones showed an activity with poly-phosphate as phosphoryl donor statistically significant. The most active mutant presented a single mutation (Glu526Lys) located in a flexible loop near of the active center. Interestingly, our theoretical studies, based on molecular dynamics simulations and hybrid Quantum Mechanics/Molecular Mechanics (QM/MM) optimizations, suggest that this mutation has an effect on the binding of the poly-P favoring a more adequate position in the active center for the reaction to take place.