Multi-level free energy simulation with a staged transformation approach.

Combining multiple levels of theory in free energy simulations to balance computational accuracy and efficiency is a promising approach for studying processes in the condensed phase. While the basic idea has been proposed and explored for quite some time, it remains challenging to achieve convergence for such multi-level free energy simulations as it requires a favorable distribution overlap between different levels of theory. Previous efforts focused on improving the distribution overlap by either altering the low-level of theory for the specific system of interest or ignoring certain degrees of freedom. Here, we propose an alternative strategy that first identifies the degrees of freedom that lead to gaps in the distributions of different levels of theory and then treats them separately with either constraints or restraints or by introducing an intermediate model that better connects the low and high levels of theory. As a result, the conversion from the low level to the high level model is done in a staged fashion that ensures a favorable distribution overlap along the way. Free energy components associated with different steps are mostly evaluated explicitly, and thus, the final result can be meaningfully compared to the rigorous free energy difference between the two levels of theory with limited and well-defined approximations. The additional free energy component calculations involve simulations at the low level of theory and therefore do not incur high computational costs. The approach is illustrated with two simple but non-trivial solution examples, and factors that dictate the reliability of the result are discussed.

[1]  Anders S. Christensen,et al.  Improving intermolecular interactions in DFTB3 using extended polarization from chemical-potential equalization. , 2015, Journal of Chemical Physics.

[2]  Walter Thiel,et al.  QM/MM studies of enzymes. , 2007, Current opinion in chemical biology.

[3]  Robert H. Wood,et al.  Systematic errors in free energy perturbation calculations due to a finite sample of configuration space: sample-size hysteresis , 1991 .

[4]  Straub,et al.  Generalized simulated annealing algorithms using Tsallis statistics: Application to conformational optimization of a tetrapeptide. , 1996, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[5]  Bernard R Brooks,et al.  Correcting for the free energy costs of bond or angle constraints in molecular dynamics simulations. , 2015, Biochimica et biophysica acta.

[6]  Steve Kaminski,et al.  Extended polarization in third-order SCC-DFTB from chemical-potential equalization. , 2012, The journal of physical chemistry. A.

[7]  Tai-Sung Lee,et al.  A variational linear-scaling framework to build practical, efficient next-generation orbital-based quantum force fields. , 2013, Journal of chemical theory and computation.

[8]  W. Thompson.  Solvation dynamics and proton transfer in nanoconfined liquids. , 2011, Annual review of physical chemistry.

[9]  Y. Shao,et al.  Efficient Computation of Free Energy Surfaces of Diels–Alder Reactions in Explicit Solvent at Ab Initio QM/MM Level , 2018, Molecules.

[10]  Thomas Simonson,et al.  Conformational substrates and uncertainty in macromolecular free energy calculations , 1993 .

[11]  Timothy J. Giese,et al.  Charge-dependent model for many-body polarization, exchange, and dispersion interactions in hybrid quantum mechanical/molecular mechanical calculations. , 2007, The Journal of chemical physics.

[12]  G. Hummer.  Fast-growth thermodynamic integration: Error and efficiency analysis , 2001 .

[13]  Noam Bernstein,et al.  Hybrid atomistic simulation methods for materials systems , 2009 .

[14]  D. Truhlar,et al.  Quantum mechanical methods for enzyme kinetics. , 2003, Annual review of physical chemistry.

[15]  Bernard R Brooks,et al.  On the convergence of multi-scale free energy simulations , 2018, Molecular simulation.

[16]  H. Woodcock,et al.  On the use of interaction energies in QM/MM free energy simulations. , 2019, Journal of chemical theory and computation.

[17]  Stefan Bruckner,et al.  Efficiency of alchemical free energy simulations. I. A practical comparison of the exponential formula, thermodynamic integration, and Bennett's acceptance ratio method , 2011, J. Comput. Chem..

[18]  Host-Guest Relative Binding Affinities at Density-Functional Theory Level from Semiempirical Molecular Dynamics Simulations. , 2019, Journal of chemical theory and computation.

[19]  Christophe Chipot,et al.  Good practices in free-energy calculations. , 2010, The journal of physical chemistry. B.

[20]  Anders S. Christensen,et al.  Semiempirical Quantum Mechanical Methods for Noncovalent Interactions for Chemical and Biochemical Applications , 2016, Chemical reviews.

[21]  Timothy J. Giese,et al.  Development of a Robust Indirect Approach for MM→QM Free Energy Calculations that Combines Force-matched Reference Potential and Bennett's Acceptance Ratio Methods. , 2019, Journal of chemical theory and computation.

[22]  M. Maroncelli,et al.  Nonreactive Dynamics in Solution: The Emerging Molecular View of Solvation Dynamics and Vibrational Relaxation , 1996 .

[23]  D. Kofke.  Free energy methods in molecular simulation , 2005 .

[24]  C. Jarzynski.  Nonequilibrium Equality for Free Energy Differences , 1996, cond-mat/9610209.

[25]  Donald G. Truhlar,et al.  Parameterization of NDDO wavefunctions using genetic algorithms. An evolutionary approach to parameterizing potential energy surfaces and direct dynamics calculations for organic reactions , 1995 .

[26]  G. Crooks.  Entropy production fluctuation theorem and the nonequilibrium work relation for free energy differences. , 1999, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[27]  Martin Karplus,et al.  Molecular properties from combined QM/MM methods. I. Analytical second derivative and vibrational calculations , 2000 .

[28]  H. Woodcock,et al.  Convergence of single-step free energy perturbation , 2017 .

[29]  B. Roux,et al.  Computations of standard binding free energies with molecular dynamics simulations. , 2009, The journal of physical chemistry. B.

[30]  Pär Söderhjelm,et al.  Converging ligand‐binding free energies obtained with free‐energy perturbations at the quantum mechanical level , 2016, J. Comput. Chem..

[31]  Michael R. Shirts,et al.  Comparison of efficiency and bias of free energies computed by exponential averaging, the Bennett acceptance ratio, and thermodynamic integration. , 2005, The Journal of chemical physics.

[32]  Ulf Ryde,et al.  Convergence of QM/MM free-energy perturbations based on molecular-mechanics or semiempirical simulations. , 2012, Physical chemistry chemical physics : PCCP.

[33]  Quantum mechanics/molecular mechanics dual Hamiltonian free energy perturbation. , 2013, The Journal of chemical physics.

[34]  David A. Kofke,et al.  Accuracy of free-energy perturbation calculations in molecular simulation. I. Modeling , 2001 .

[35]  Yong Wang,et al.  P450 enzymes: their structure, reactivity, and selectivity-modeled by QM/MM calculations. , 2010, Chemical reviews.

[36]  M. Elstner,et al.  Parametrization and Benchmark of DFTB3 for Organic Molecules. , 2013, Journal of chemical theory and computation.

[37]  U. Rothlisberger,et al.  Mixed Quantum Mechanical/Molecular Mechanical Molecular Dynamics Simulations of Biological Systems in Ground and Electronically Excited States. , 2015, Chemical reviews.

[38]  H Lee Woodcock,et al.  Accelerating QM/MM Free Energy Computations via Intramolecular Force Matching. , 2018, Journal of chemical theory and computation.

[39]  Ye Mei,et al.  Accelerated Computation of Free Energy Profile at ab Initio Quantum Mechanical/Molecular Mechanics Accuracy via a Semi-Empirical Reference Potential. I. Weighted Thermodynamics Perturbation. , 2018, Journal of chemical theory and computation.

[40]  M. Karplus,et al.  A combined quantum mechanical and molecular mechanical potential for molecular dynamics simulations , 1990 .

[41]  Jiali Gao,et al.  Absolute free energy of solvation from Monte Carlo simulations using combined quantum and molecular mechanical potentials , 1992 .

[42]  U. Ryde.  How Many Conformations Need To Be Sampled To Obtain Converged QM/MM Energies? The Curse of Exponential Averaging. , 2017, Journal of chemical theory and computation.

[43]  Alexander D. MacKerell,et al.  Optimization of the additive CHARMM all-atom protein force field targeting improved sampling of the backbone φ, ψ and side-chain χ(1) and χ(2) dihedral angles. , 2012, Journal of chemical theory and computation.

[44]  M. Retegan,et al.  Free energy calculations using dual-level Born-Oppenheimer molecular dynamics. , 2010, The Journal of chemical physics.

[45]  G. Torrie,et al.  Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling , 1977 .

[46]  Phillip S. Hudson,et al.  Computing converged free energy differences between levels of theory via nonequilibrium work methods: Challenges and opportunities , 2017, J. Comput. Chem..

[47]  Michael Gaus,et al.  DFTB3: Extension of the self-consistent-charge density-functional tight-binding method (SCC-DFTB). , 2011, Journal of chemical theory and computation.

[48]  Arieh Warshel,et al.  Microscopic models for quantum mechanical calculations of chemical processes in solutions: LD/AMPAC and SCAAS/AMPAC calculations of solvation energies , 1992 .

[49]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[50]  Arieh Warshel,et al.  Towards accurate ab initio QM/MM calculations of free-energy profiles of enzymatic reactions. , 2006, The journal of physical chemistry. B.

[51]  David A. Kofke,et al.  Accuracy of free-energy perturbation calculations in molecular simulation. II. Heuristics , 2001 .

[52]  Adrian J Mulholland,et al.  High-accuracy computation of reaction barriers in enzymes. , 2006, Angewandte Chemie.

[53]  Di Wu,et al.  Phase-space overlap measures. II. Design and implementation of staging methods for free-energy calculations. , 2005, The Journal of chemical physics.

[54]  R. Friesner,et al.  Ab initio quantum chemical and mixed quantum mechanics/molecular mechanics (QM/MM) methods for studying enzymatic catalysis. , 2005, Annual review of physical chemistry.

[55]  M. Levitt,et al.  Theoretical studies of enzymic reactions: dielectric, electrostatic and steric stabilization of the carbonium ion in the reaction of lysozyme. , 1976, Journal of molecular biology.

[56]  Charles H. Bennett,et al.  Efficient estimation of free energy differences from Monte Carlo data , 1976 .

[57]  Alexander D. MacKerell,et al.  CHARMM general force field: A force field for drug‐like molecules compatible with the CHARMM all‐atom additive biological force fields , 2009, J. Comput. Chem..

[58]  Florian Sittel,et al.  Perspective: Identification of collective variables and metastable states of protein dynamics. , 2018, The Journal of chemical physics.

[59]  Q. Cui.  Perspective: Quantum mechanical methods in biochemistry and biophysics. , 2016, The Journal of chemical physics.

[60]  Hao Hu,et al.  Free energies of chemical reactions in solution and in enzymes with ab initio quantum mechanics/molecular mechanics methods. , 2008, Annual review of physical chemistry.

[61]  D. Truhlar,et al.  Mechanisms and free energies of enzymatic reactions. , 2006, Chemical reviews.

[62]  Jianpeng Ma,et al.  CHARMM: The biomolecular simulation program , 2009, J. Comput. Chem..

[63]  Michael Gaus,et al.  Density functional tight binding: application to organic and biological molecules , 2014 .

[64]  Yihang Wang,et al.  Machine learning approaches for analyzing and enhancing molecular dynamics simulations. , 2019, Current opinion in structural biology.

[65]  Sergio Martí,et al.  Improving the QM/MM Description of Chemical Processes:  A Dual Level Strategy To Explore the Potential Energy Surface in Very Large Systems. , 2005, Journal of chemical theory and computation.

[66]  Di Wu,et al.  Phase-space overlap measures. I. Fail-safe bias detection in free energies calculated by molecular simulation. , 2005, The Journal of chemical physics.

[67]  Arieh Warshel,et al.  Paradynamics: an effective and reliable model for ab initio QM/MM free-energy calculations and related tasks. , 2011, The journal of physical chemistry. B.

[68]  Guohui Li,et al.  Development of effective quantum mechanical/molecular mechanical (QM/MM) methods for complex biological processes. , 2006, The journal of physical chemistry. B.

[69]  Yuko Okamoto,et al.  QM/MM free energy simulations: recent progress and challenges , 2016, Molecular simulation.

[70]  A. Mulholland,et al.  QM/MM modelling of drug-metabolizing enzymes. , 2014, Current topics in medicinal chemistry.

[71]  Donald G Truhlar,et al.  Multidimensional tunneling, recrossing, and the transmission coefficient for enzymatic reactions. , 2006, Chemical reviews.

[72]  Bernard R Brooks,et al.  Vibrational subsystem analysis: A method for probing free energies and correlations in the harmonic limit. , 2008, The Journal of chemical physics.

[73]  David A. Kofke,et al.  Appropriate methods to combine forward and reverse free-energy perturbation averages , 2003 .

[74]  David A. Kofke,et al.  On the sampling requirements for exponential-work free-energy calculations , 2006 .

[75]  Fiona L. Kearns,et al.  The Good, the Bad, and the Ugly: “HiPen”, a New Dataset for Validating (S)QM/MM Free Energy Simulations , 2019, Molecules.

[76]  Weitao Yang,et al.  Molecular Dynamics Simulations with Quantum Mechanics/Molecular Mechanics and Adaptive Neural Networks. , 2018, Journal of chemical theory and computation.

[77]  K. Merz,et al.  Combined Quantum Mechanical/Molecular Mechanical Methodologies Applied to Biomolecular Systems , 1999 .

[78]  Thomas B Woolf,et al.  Theory of a systematic computational error in free energy differences. , 2002, Physical review letters.

[79]  Gerhard König,et al.  Multiscale Free Energy Simulations: An Efficient Method for Connecting Classical MD Simulations to QM or QM/MM Free Energies Using Non-Boltzmann Bennett Reweighting Schemes , 2014, Journal of chemical theory and computation.