Best Practices on QM/MM Simulations of Biological Systems

• Published in: Journal of chemical information and modeling Vol. 63; no. 9; pp. 2609 - 2627
• Main Authors: Clemente, Camila M., Capece, Luciana, Martí, Marcelo A.
• Format: Journal Article
• Language: English
• Published: United States American Chemical Society 08.05.2023
• Subjects: Best practice; Chorismate Mutase; Computational chemistry; Electronic structure; Energy methods; Entropy; Free energy; Hemoglobin; Models, Biological; Molecular dynamics; Molecular Dynamics Simulation; Proteins - chemistry; Quantum Theory; Reaction mechanisms; Research projects; Simulation
• ISSN: 1549-9596; 1549-960X
• DOI: 10.1021/acs.jcim.2c01522

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.