Accelerating Ab Initio Quantum Mechanical and Molecular Mechanical (QM/MM) Molecular Dynamics Simulations with Multiple Time Step Integration and a Recalibrated Semiempirical QM/MM Hamiltonian

The Effects of Conformational Sampling and QM Region Size in QM/MM Simulations: An Adaptive QM/MM Study With Model Systems

Molecular properties in combined quantum mechanics and molecular mechanics (QM/MM) simulations have been shown to be dependent on the size of the quantum mechanical (QM) region and the amount of conformational sampling. Previous studies have largely focused on enzymatic systems, which have made it difficult to distinguish the effects of QM region size and conformational sampling from other factors including QM-MM boundary artifacts and the boundary effects. This study uses the difference-based adaptive solvation QM/MM method to investigate the tautomerization reactions of alanine and aspartate in explicit solvent. The choice of computationally tractable systems enables the decoupling of QM region size effects from other factors and a direct comparison of free energy surfaces with potential energy surfaces (PESs). The results show that (1) it is crucial to properly account for thermal fluctuations along the reaction pathways, and (2) free energy surfaces converge rapidly with increasing QM region size, whereas charge transfer requires a slightly larger QM region to achieve convergence. These findings are expected to guide future studies of enzymatic systems and other complex systems where QM/MM methods are applied.

Assessment of Free Energies From Electrostatic Embedding Density Functional Tight Binding-Based/Molecular Mechanics in Periodic Boundary Conditions

Electrostatic embedding quantum mechanics/molecular mechanics (QM/MM) methods in periodic boundary conditions (PBC) can successfully describe the condensed phase reactivity of a fragment treated at the QM level with an atomistic description of an electrostatic environment treated at the MM level. The computational cost of ab initio QM methods limits the phase space sampling, thus affecting statistical quantities like free energies. Here, we describe the implementation of a PBC-adapted QM/MM model based on the semi-empirical density-functional based tight-binding (DFTB) method within the GAMESS-US quantum package interfaced with Tinker. Further, we take advantage of the free energy methods provided by a newly developed interface with the PLUMED plugin. The versatility of the implementation is illustrated by the prediction of the free energy profile for three different families of reactions in solution. Overall, using the DFTB/MM, it has been possible to obtain results that are at least in a qualitatively agreement with respect to the experimental data or high-level ab initio simulations.

High-performance QM/MM Enhanced Sampling Molecular Dynamics Simulations with GENESIS SPDYN and QSimulate-QM

A new module for quantum mechanical/molecular mechanical (QM/MM) calculations is implemented in a molecular dynamics (MD) program, SPDYN in GENESIS, interfaced with an electronic structure program, QSimulate-QM. The periodic boundary condition (PBC) in QM/MM simulation is incorporated via QM calculation in real space with duplicated MM charges and particle mesh Ewald (PME) calculation with QM and MM charges. A highly optimized code is implemented in QSimulate-QM, particularly for the density functional tight-binding (DFTB) method, where the interaction between the QM and MM regions is computed utilizing multipole expansions. Together with highly parallelized algorithms in SPDYN, the developed program performs MD simulations based on DFTB in the QM size of ∼100 atoms and MM of ∼100,000 atoms with a better performance than 1 ns/day using one computer node. This feature paves the way for QM/MM-MD enhanced sampling simulations. Various enhanced sampling methods in GENESIS, namely, generalized replica exchange solute tempering (gREST), replica-exchange umbrella sampling (REUS), and path sampling with the string method, are demonstrated at the QM/MM level to compute the Ramachandran plot of alanine dipeptide and the potential of mean force (PMF) of a proton transfer reaction in an enzyme.

Machine Learning Quantum Mechanical/Molecular Mechanical Potentials: Evaluating Transferability in Dihydrofolate Reductase-Catalyzed Reactions

Integrating machine learning potentials (MLPs) with quantum mechanical/molecular mechanical (QM/MM) free energy simulations has emerged as a powerful approach for studying enzymatic catalysis. However, its practical application has been hindered by the time-consuming process of generating the necessary training, validation, and test data for MLP models through QM/MM simulations. Furthermore, the entire process needs to be repeated for each specific enzyme system and reaction. To overcome this bottleneck, it is required that trained MLPs exhibit transferability across different enzyme environments and reacting species, thereby eliminating the need for retraining with each new enzyme variant. In this study, we explore this potential by evaluating the transferability of a pretrained ΔMLP model across different enzyme mutations within the MM environment using the QM/MM-based ML architecture developed by Pan, X. J. Chem. Theory Comput. 2021, 17(9), 5745-5758. The study includes scenarios such as single point substitutions, a homologous enzyme from different species, and even a transition to an aqueous environment, where the last two systems have MM environment that is substantially different from that used in MLP training. The results show that the ΔMLP model effectively captures and predicts the effects of enzyme mutations on electrostatic interactions, producing reliable free energy profiles of enzyme-catalyzed reactions without the need for retraining. The study also identified notable limitations in transferability, particularly when transitioning from enzyme to water-rich MM environments. Overall, this study demonstrates the robustness of the Pan et al.'s QM/MM-based ML architecture for application to diverse enzyme systems, as well as the need for further research and the development of more sophisticated MLP models and training methods.

Exploring CRISPR-Cas9 HNH-Domain-Catalyzed DNA Cleavage Using Accelerated Quantum Mechanical Molecular Mechanical Free Energy Simulation

The target DNA (tDNA) cleavage catalyzed by the CRISPR Cas9 enzyme is a critical step in the Cas9-based genome editing technologies. Previously, the tDNA cleavage from an active SpyCas9 enzyme conformation was modeled by Palermo and co-workers (Nierzwicki et al., Nat. Catal. 2022 5, 912) using ab initio quantum mechanical molecular mechanical (ai-QM/MM) free energy simulations, where the free energy barrier was found to be more favorable than that from a pseudoactive enzyme conformation. In this work, we performed ai-QM/MM simulations based on another catalytically active conformation (PDB 7Z4J) of the Cas9 HNH domain from cryo-electron microscopy experiments. For the wildtype enzyme, we acquired a free energy profile for the tDNA cleavage that is largely consistent with the previous report. Furthermore, we explored the role of the active-site K866 residue on the catalytic efficiency by modeling the K866A mutant and found that the K866A mutation increased the reaction free energy barrier, which is consistent with the experimentally observed reduction in the enzyme activity.

GENESIS 2.1: High-Performance Molecular Dynamics Software for Enhanced Sampling and Free-Energy Calculations for Atomistic, Coarse-Grained, and Quantum Mechanics/Molecular Mechanics Models

GENeralized-Ensemble SImulation System (GENESIS) is a molecular dynamics (MD) software developed to simulate the conformational dynamics of a single biomolecule, as well as molecular interactions in large biomolecular assemblies and between multiple biomolecules in cellular environments. To achieve the latter purpose, the earlier versions of GENESIS emphasized high performance in atomistic MD simulations on massively parallel supercomputers, with or without graphics processing units (GPUs). Here, we implemented multiscale MD simulations that include atomistic, coarse-grained, and hybrid quantum mechanics/molecular mechanics (QM/MM) calculations. They demonstrate high performance and are integrated with enhanced conformational sampling algorithms and free-energy calculations without using external programs except for the QM programs. In this article, we review new functions, molecular models, and other essential features in GENESIS version 2.1 and discuss ongoing developments for future releases.

Multiple-time scale integration method based on an interpolated potential energy surface for ab initio path integral molecular dynamics

A new multiple-time scale integration method is presented that propagates ab initio path integral molecular dynamics (PIMD). This method uses a large time step to generate an approximate geometrical configuration whose energy and gradient are evaluated at the level of an ab initio method, and then, a more precise integration scheme, e.g., the Bulirsch-Stoer method or velocity Verlet integration with a smaller time step, is used to integrate from the previous step using the computationally efficient interpolated potential energy surface constructed from two consecutive points. This method makes the integration of PIMD more efficient and accurate compared with the velocity Verlet integration. A Nosé-Hoover chain thermostat combined with this new multiple-time scale method has good energy conservation even with a large time step, which is usually challenging in velocity Verlet integration for PIMD due to the very small chain mass when a large number of beads are used. The new method is used to calculate infrared spectra and free energy profiles to demonstrate its accuracy and capabilities.

Perspectives on Computational Enzyme Modeling: From Mechanisms to Design and Drug Development

Understanding enzyme mechanisms is essential for unraveling the complex molecular machinery of life. In this review, we survey the field of computational enzymology, highlighting key principles governing enzyme mechanisms and discussing ongoing challenges and promising advances. Over the years, computer simulations have become indispensable in the study of enzyme mechanisms, with the integration of experimental and computational exploration now established as a holistic approach to gain deep insights into enzymatic catalysis. Numerous studies have demonstrated the power of computer simulations in characterizing reaction pathways, transition states, substrate selectivity, product distribution, and dynamic conformational changes for various enzymes. Nevertheless, significant challenges remain in investigating the mechanisms of complex multistep reactions, large-scale conformational changes, and allosteric regulation. Beyond mechanistic studies, computational enzyme modeling has emerged as an essential tool for computer-aided enzyme design and the rational discovery of covalent drugs for targeted therapies. Overall, enzyme design/engineering and covalent drug development can greatly benefit from our understanding of the detailed mechanisms of enzymes, such as protein dynamics, entropy contributions, and allostery, as revealed by computational studies. Such a convergence of different research approaches is expected to continue, creating synergies in enzyme research. This review, by outlining the ever-expanding field of enzyme research, aims to provide guidance for future research directions and facilitate new developments in this important and evolving field.

Exploration of biochemical reactivity with a QM/MM growing string method

In this work, we have implemented the single-ended growing string method using a hybrid internal/Cartesian coordinate scheme within our in-house QM/MM package, QoMMMa, representing the first implementation of the growing string method in the QM/MM framework. The goal of the implementation was to facilitate generation of QM/MM reaction pathways with minimal user input, and also to improve the quality of the pathways generated as compared to the widely used adiabatic mapping approach. We have validated the algorithm against a reaction which has been studied extensively in previous computational investigations - the Claisen rearrangement catalysed by chorismate mutase. The nature of the transition state and the height of the barrier was predicted well using our algorithm, where more than 88% of the pathways generated were deemed to be of production quality. Directly compared to using adiabatic mapping, we found that while our QM/MM single-ended growing string method is slightly less efficient, it readily produces reaction pathways with fewer discontinuites and thus minimises the need for involved remapping of unsatisfactory energy profiles.

Accurate and Efficient Multilevel Free Energy Simulations with Neural Network-Assisted Enhanced Sampling

Free energy differences (ΔF) are essential to quantitative characterization and understanding of chemical and biological processes. Their direct estimation with an accurate quantum mechanical potential is of great interest and yet impractical due to high computational cost and incompatibility with typical alchemical free energy protocols. One promising solution is the multilevel free energy simulation in which the estimate of ΔF at an inexpensive low level of theory is combined with the correction toward a higher level of theory. The poor configurational overlap generally expected between the two levels of theory, however, presents a major challenge. We overcome this challenge by using a deep neural network model and enhanced sampling simulations. An adversarial autoencoder is used to identify a low-dimensional (latent) space that compactly represents the degrees of freedom that encode the distinct distributions at the two levels of theory. Enhanced sampling in this latent space is then used to drive the sampling of configurations that predominantly contribute to the free energy correction. Results for both gas phase and condensed phase systems demonstrate that this data-driven approach offers high accuracy and efficiency with great potential for scalability to complex systems.

Bridging semiempirical and ab initio QM/MM potentials by Gaussian process regression and its sparse variants for free energy simulation

Free energy simulations that employ combined quantum mechanical and molecular mechanical (QM/MM) potentials at ab initio QM (AI) levels are computationally highly demanding. Here, we present a machine-learning-facilitated approach for obtaining AI/MM-quality free energy profiles at the cost of efficient semiempirical QM/MM (SE/MM) methods. Specifically, we use Gaussian process regression (GPR) to learn the potential energy corrections needed for an SE/MM level to match an AI/MM target along the minimum free energy path (MFEP). Force modification using gradients of the GPR potential allows us to improve configurational sampling and update the MFEP. To adaptively train our model, we further employ the sparse variational GP (SVGP) and streaming sparse GPR (SSGPR) methods, which efficiently incorporate previous sample information without significantly increasing the training data size. We applied the QM-(SS)GPR/MM method to the solution-phase SN2 Menshutkin reaction, NH3+CH3Cl→CH3NH3++Cl-, using AM1/MM and B3LYP/6-31+G(d,p)/MM as the base and target levels, respectively. For 4000 configurations sampled along the MFEP, the iteratively optimized AM1-SSGPR-4/MM model reduces the energy error in AM1/MM from 18.2 to 4.4 kcal/mol. Although not explicitly fitting forces, our method also reduces the key internal force errors from 25.5 to 11.1 kcal/mol/Å and from 30.2 to 10.3 kcal/mol/Å for the N-C and C-Cl bonds, respectively. Compared to the uncorrected simulations, the AM1-SSGPR-4/MM method lowers the predicted free energy barrier from 28.7 to 11.7 kcal/mol and decreases the reaction free energy from -12.4 to -41.9 kcal/mol, bringing these results into closer agreement with their AI/MM and experimental benchmarks.