Multiple Environment Single System Quantum Mechanical/Molecular Mechanical (MESS-QM/MM) Calculations. 1. Estimation of Polarization Energies

In combined quantum mechanical/molecular mechanical (QM/MM) free energy calculations, it is often advantageous to have a frozen geometry for the quantum mechanical (QM) region. For such multiple-environment single-system (MESS) cases, two schemes are proposed here for estimating the polarization energy: the first scheme, termed MESS-E, involves a Roothaan step extrapolation of the self-consistent field (SCF) energy; whereas the other scheme, termed MESS-H, employs a Newton–Raphson correction using an approximate inverse electronic Hessian of the QM region (which is constructed only once). Both schemes are extremely efficient, because the expensive Fock updates and SCF iterations in standard QM/MM calculations are completely avoided at each configuration. They produce reasonably accurate QM/MM polarization energies: MESS-E can predict the polarization energy within 0.25 kcal/mol in terms of the mean signed error for two of our test cases, solvated methanol and solvated β-alanine, using the M06-2X or ωB97X-D functionals; MESS-H can reproduce the polarization energy within 0.2 kcal/mol for these two cases and for the oxyluciferin–luciferase complex, if the approximate inverse electronic Hessians are constructed with sufficient accuracy.

[1]  So Hirata,et al.  Time-dependent density functional theory for radicals: An improved description of excited states with substantial double excitation character , 1999 .

[2]  Rustam Z. Khaliullin,et al.  An efficient self-consistent field method for large systems of weakly interacting components. , 2006, The Journal of chemical physics.

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

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

[5]  Michiel Sprik,et al.  Free energy from constrained molecular dynamics , 1998 .

[6]  T. P. Straatsma,et al.  Treatment of rotational isomers in free energy evaluations. Analysis of the evaluation of free energy differences by molecular dynamics simulations of systems with rotational isomeric states , 1989 .

[7]  P. Senet.  Kohn-Sham orbital formulation of the chemical electronic responses, including the hardness , 1997 .

[8]  P. Hohenberg,et al.  Inhomogeneous Electron Gas , 1964 .

[9]  Y. M. Rhee,et al.  Dynamics on the electronically excited state surface of the bioluminescent firefly luciferase-oxyluciferin system. , 2011, Journal of the American Chemical Society.

[10]  Yihan Shao,et al.  Dual-basis second-order Moller-Plesset perturbation theory: A reduced-cost reference for correlation calculations. , 2006, The Journal of chemical physics.

[11]  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.

[12]  Stefan Boresch,et al.  THE ROLE OF BONDED TERMS IN FREE ENERGY SIMULATIONS : 1. THEORETICAL ANALYSIS , 1999 .

[13]  T. P. Straatsma,et al.  Treatment of rotational isomers in free energy calculations. II. Molecular dynamics simulation study of 18‐crown‐6 in aqueous solution as an example of systems with large numbers of rotational isomeric states , 1989 .

[14]  D. Truhlar,et al.  The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functionals , 2008 .

[15]  S. Kato,et al.  Ab Initio Molecular Orbital Theory on Intramolecular Charge Polarization: Effect of Hydrogen Abstraction on the Charge Sensitivity of Aromatic and Nonaromatic Species , 1997 .

[16]  C. Cramer,et al.  A universal approach to solvation modeling. , 2008, Accounts of chemical research.

[17]  K. Burke.  Perspective on density functional theory. , 2012, The Journal of chemical physics.

[18]  E. Davidson.  The iterative calculation of a few of the lowest eigenvalues and corresponding eigenvectors of large real-symmetric matrices , 1975 .

[19]  P. Ayers,et al.  Potentialphilicity and potentialphobicity: Reactivity indicators for external potential changes from density functional reactivity theory. , 2009, The Journal of chemical physics.

[20]  Roberto Car,et al.  Free energy profile along a discretized reaction path via the hyperplane constraint force and torque. , 2005, The Journal of chemical physics.

[21]  Wilfred F. van Gunsteren,et al.  Basic ingredients of free energy calculations: A review , 2009, J. Comput. Chem..

[22]  M. Berkowitz,et al.  Molecular hardness and softness, local hardness and softness, hardness and softness kernels, and relations among these quantities , 1988 .

[23]  S. Fias,et al.  Evaluating and interpreting the chemical relevance of the linear response kernel for atoms II: open shell. , 2014, Physical chemistry chemical physics : PCCP.

[24]  Anna I. Krylov,et al.  Q‐Chem: an engine for innovation , 2013 .

[25]  C. Sherrill.  Frontiers in electronic structure theory. , 2010, The Journal of chemical physics.

[26]  Arieh Warshel,et al.  Hybrid ab Initio Quantum Mechanics/Molecular Mechanics Calculations of Free Energy Surfaces for Enzymatic Reactions: The Nucleophilic Attack in Subtilisin , 1998 .

[27]  M. Head‐Gordon,et al.  Dual-basis analytic gradients. 1. Self-consistent field theory. , 2006, The journal of physical chemistry. A.

[28]  Trygve Helgaker,et al.  Quantitative quantum chemistry , 2008 .

[29]  Wilfred F. van Gunsteren,et al.  Computation of free energy , 2002 .

[30]  R. Parr.  Density-functional theory of atoms and molecules , 1989 .

[31]  Walter Thiel,et al.  QM/MM methods for biomolecular systems. , 2009, Angewandte Chemie.

[32]  Andrew T. B. Gilbert,et al.  Density functional triple jumping. , 2010, Physical chemistry chemical physics : PCCP.

[33]  Martin Head-Gordon,et al.  Quantum chemistry and molecular processes , 1996 .

[34]  Adrian J Mulholland,et al.  Computational enzymology. , 2010, Chemical communications.

[35]  Peter Pulay,et al.  Ultrafast Quantum Mechanics/Molecular Mechanics Monte Carlo simulations using generalized multipole polarizabilities , 2012 .

[36]  A. Shi.  Self-Consistent Field Theory , 2013 .

[37]  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.

[38]  H. Bernhard Schlegel,et al.  Geometry optimization , 2011 .

[39]  Shawn T. Brown,et al.  Advances in methods and algorithms in a modern quantum chemistry program package. , 2006, Physical chemistry chemical physics : PCCP.

[40]  Tamar Schlick,et al.  Geometry Optimization: 2 , 2002 .

[41]  P. Senet.  Nonlinear electronic responses, Fukui functions and hardnesses as functionals of the ground‐state electronic density , 1996 .

[42]  M. Karplus,et al.  The Jacobian factor in free energy simulations , 1996 .

[43]  Clemens C. J. Roothaan,et al.  New Developments in Molecular Orbital Theory , 1951 .

[44]  Pär Söderhjelm,et al.  On the Convergence of QM/MM Energies. , 2011, Journal of chemical theory and computation.

[45]  P. Pulay,et al.  Efficient calculation of the energy of a molecule in an arbitrary electric field , 2009 .

[46]  J. Pople,et al.  Self‐consistent molecular orbital methods. XX. A basis set for correlated wave functions , 1980 .

[47]  Masataka Nagaoka,et al.  Structure optimization via free energy gradient method: Application to glycine zwitterion in aqueous solution , 2000 .

[48]  Bernard R Brooks,et al.  Maintain rigid structures in Verlet based cartesian molecular dynamics simulations. , 2012, The Journal of chemical physics.

[49]  A. Becke,et al.  Density-functional exchange-energy approximation with correct asymptotic behavior. , 1988, Physical review. A, General physics.

[50]  M. V. Ganduglia-Pirovano,et al.  Reactivity kernels, the normal modes of chemical reactivity, and the hardness and softness spectra , 1967 .

[51]  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..

[52]  Adrian Saldanha,et al.  Structural basis for the spectral difference in luciferase bioluminescence , 2006, Nature.

[53]  Ewald mesh method for quantum mechanical calculations. , 2012, The Journal of chemical physics.

[54]  Bernard R. Brooks,et al.  Comparing normal modes across different models and scales: Hessian reduction versus coarse‐graining , 2012, J. Comput. Chem..

[55]  W. Briels,et al.  THE CALCULATION OF FREE-ENERGY DIFFERENCES BY CONSTRAINED MOLECULAR-DYNAMICS SIMULATIONS , 1998 .

[56]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[57]  Takeshi Yamamoto.  Variational and perturbative formulations of quantum mechanical/molecular mechanical free energy with mean-field embedding and its analytical gradients. , 2008, The Journal of chemical physics.

[58]  Weitao Yang,et al.  Challenges for density functional theory. , 2012, Chemical reviews.

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

[60]  A. Morita,et al.  The Charge Response Kernel with Modified Electrostatic Potential Charge Model , 2002 .

[61]  Martin Head-Gordon,et al.  Approaching the Basis Set Limit in Density Functional Theory Calculations Using Dual Basis Sets without Diagonalization , 2004 .

[62]  Bernard R. Brooks,et al.  Computing the Free Energy along a Reaction Coordinate Using Rigid Body Dynamics , 2014, Journal of chemical theory and computation.

[63]  P. Geerlings,et al.  Evaluating and Interpreting the Chemical Relevance of the Linear Response Kernel for Atoms. , 2013, Journal of chemical theory and computation.

[64]  Bernard R. Brooks,et al.  Interfacing Q‐Chem and CHARMM to perform QM/MM reaction path calculations , 2007, J. Comput. Chem..

[65]  Constrained reaction coordinate dynamics for systems with constraints , 2003 .

[66]  J. Pople,et al.  Self—Consistent Molecular Orbital Methods. XII. Further Extensions of Gaussian—Type Basis Sets for Use in Molecular Orbital Studies of Organic Molecules , 1972 .

[67]  Parr,et al.  Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. , 1988, Physical review. B, Condensed matter.

[68]  H. Jónsson,et al.  Optimization of hyperplanar transition states , 2001 .

[69]  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.

[70]  Ulf Ryde,et al.  Accurate QM/MM Free Energy Calculations of Enzyme Reactions:  Methylation by Catechol O-Methyltransferase. , 2005, Journal of chemical theory and computation.

[71]  J. C. D. Silva,et al.  Computational Studies of the Luciferase Light-Emitting Product: Oxyluciferin. , 2011, Journal of chemical theory and computation.

[72]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[73]  Wei Yang,et al.  Catalytic mechanism of RNA backbone cleavage by ribonuclease H from quantum mechanics/molecular mechanics simulations. , 2011, Journal of the American Chemical Society.

[74]  Masataka Nagaoka,et al.  Hydrated structure of ammonia–water molecule pair via the free energy gradient method: Realization of zero gradient and force balance on free energy surfaces , 2003 .

[75]  Marco De Vivo,et al.  The increasing role of QM/MM in drug discovery. , 2012, Advances in protein chemistry and structural biology.

[76]  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.

[77]  Yihan Shao,et al.  Efficient evaluation of the Coulomb force in density-functional theory calculations , 2001 .

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

[79]  Ross C. Walker,et al.  The implementation of a fast and accurate QM/MM potential method in Amber , 2008, J. Comput. Chem..

[80]  Weitao Yang,et al.  Free energy calculation on enzyme reactions with an efficient iterative procedure to determine minimum energy paths on a combined ab initio QM/MM potential energy surface , 2000 .

[81]  M. Head‐Gordon,et al.  Long-range corrected hybrid density functionals with damped atom-atom dispersion corrections. , 2008, Physical chemistry chemical physics : PCCP.

[82]  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.

[83]  G. Ackland.  Calculation of free energies from ab initio calculation , 2002 .

[84]  G. Ciccotti,et al.  Rare events by constrained molecular dynamics , 2000 .

[85]  J. Herbert,et al.  Periodic boundary conditions for QM/MM calculations: Ewald summation for extended Gaussian basis sets. , 2013, The Journal of chemical physics.

[86]  Arieh Warshel,et al.  Towards Quantitative Computer‐Aided Studies of Enzymatic Enantioselectivity: The Case of Candida antarctica Lipase A , 2012, Chembiochem : a European journal of chemical biology.

[87]  Yoshiyuki Koyano,et al.  Transition-state characterization of the ammonia ionization process in aqueous solution via the free-energy gradient method. , 2006, The journal of physical chemistry. A.

[88]  P. Ayers.  Strategies for computing chemical reactivity indices , 2001 .

[89]  Davide Branduardi,et al.  String method for calculation of minimum free-energy paths in Cartesian space in freely-tumbling systems. , 2013, Journal of chemical theory and computation.

[90]  A. Becke.  Perspective: Fifty years of density-functional theory in chemical physics. , 2014, The Journal of chemical physics.

[91]  H. Bernhard Schlegel,et al.  Exploring potential energy surfaces for chemical reactions: An overview of some practical methods , 2003, J. Comput. Chem..

[92]  Zhong-Zhi Yang,et al.  Calculation of the linear response function by the atom-bond electronegativity equalization method (ABEEM) , 2000 .

[93]  D. Herschbach,et al.  Molecular Partition Functions in Terms of Local Properties , 1959 .

[94]  Darrin M York,et al.  An Efficient Linear-Scaling Ewald Method for Long-Range Electrostatic Interactions in Combined QM/MM Calculations. , 2005, Journal of chemical theory and computation.

[95]  S. Goedecker.  Linear scaling electronic structure methods , 1999 .

[96]  Manuel A. Aguilar,et al.  A new method to locate saddle points for reactions in solution by using the free‐energy gradient method and the mean field approximation , 2004, J. Comput. Chem..

[97]  Lochana C. Menikarachchi,et al.  QM/MM approaches in medicinal chemistry research. , 2010, Current topics in medicinal chemistry.

[98]  Christophe Chipot,et al.  Free Energy Calculations. The Long and Winding Gilded Road , 2002 .

[99]  B. Brooks,et al.  Reaction Path Optimization and Sampling Methods and Their Applications for Rare Events , 2012 .

[100]  Ulf Ryde,et al.  Quantum mechanical free energy barrier for an enzymatic reaction. , 2005, Physical review letters.

[101]  Eric Vanden-Eijnden,et al.  Some recent techniques for free energy calculations , 2009, J. Comput. Chem..

[102]  Rustam Z. Khaliullin,et al.  Analysis of charge transfer effects in molecular complexes based on absolutely localized molecular orbitals. , 2008, The Journal of chemical physics.

[103]  A. Becke.  A New Mixing of Hartree-Fock and Local Density-Functional Theories , 1993 .

[104]  Yihan Shao,et al.  An improved J matrix engine for density functional theory calculations , 2000 .

[105]  P W Ayers,et al.  Variational principles for describing chemical reactions. Reactivity indices based on the external potential. , 2001, Journal of the American Chemical Society.

[106]  J. Kirkwood.  Statistical Mechanics of Fluid Mixtures , 1935 .

[107]  A. Warshel.  Computer simulations of enzyme catalysis: methods, progress, and insights. , 2003, Annual review of biophysics and biomolecular structure.

[108]  W. Kohn,et al.  Self-Consistent Equations Including Exchange and Correlation Effects , 1965 .

[109]  Free energy from molecular dynamics with multiple constraints , 2000 .

[110]  H. C. Andersen.  Rattle: A “velocity” version of the shake algorithm for molecular dynamics calculations , 1983 .

[111]  John A. Pople,et al.  Self‐consistent molecular orbital methods. XVIII. Constraints and stability in Hartree–Fock theory , 1977 .