Analysis of polarization in QM/MM modelling of biologically relevant hydrogen bonds

Combined quantum mechanics/molecular mechanics (QM/MM) methods are increasingly important for the study of chemical reactions and systems in condensed phases. Here, we have tested the accuracy of a density functional theory-based QM/MM implementation (B3LYP/6-311+G(d,p)/CHARMM27) on a set of biologically relevant interactions by comparison with full QM calculations. Intermolecular charge transfer due to hydrogen bond formation is studied to assess the severity of spurious polarization of QM atoms by MM point charges close to the QM/MM boundary. The changes in total electron density and natural bond orbital atomic charges due to hydrogen bond formation in selected complexes obtained at the QM/MM level are compared with full QM results. It is found that charge leakage from the QM atoms to MM atomic point charges close to the QM/MM boundary is not a serious problem, at least with limited basis sets. The results are encouraging in showing that important properties of key biomolecular interactions can be treated well at the QM/MM level employing good-quality levels of QM theory.

[1]  W. Thiel,et al.  Hybrid Models for Combined Quantum Mechanical and Molecular Mechanical Approaches , 1996 .

[2]  Keiji Morokuma,et al.  The IMOMO method: Integration of different levels of molecular orbital approximations for geometry optimization of large systems: Test for n‐butane conformation and SN2 reaction: RCl+Cl− , 1996 .

[3]  Thomas Heine,et al.  Description of electron delocalization via the analysis of molecular fields. , 2005, Chemical reviews.

[4]  Lars Ridder,et al.  Transition state stabilization and substrate strain in enzyme catalysis: ab initio QM/MM modelling of the chorismate mutase reaction. , 2004, Organic & biomolecular chemistry.

[5]  Dennis R. Salahub,et al.  Gaussian density functional calculations on hydrogen-bonded systems , 1992 .

[6]  P K Biswas,et al.  A regularized and renormalized electrostatic coupling Hamiltonian for hybrid quantum-mechanical-molecular-mechanical calculations. , 2005, The Journal of chemical physics.

[7]  K. Morokuma,et al.  Molecular orbital studies of hydrogen bonds. IX. Electron distribution analysis , 1975 .

[8]  Kuan-Jiuh Lin,et al.  Charge Density Distributions, Interaction Energies, and Electrostatic Potentials of Hydrogen Bonds , 1994 .

[9]  Richard A. Friesner,et al.  A mixed quantum mechanics/molecular mechanics (QM/MM) method for large-scale modeling of chemistry in protein environments , 2000, J. Comput. Chem..

[10]  Martin J. Field,et al.  HYBRID QUANTUM MECHANICAL/MOLECULAR MECHANICAL FLUCTUATING CHARGE MODELS FOR CONDENSED PHASE SIMULATIONS , 1997 .

[11]  D. Peeters,et al.  Insight into metal-phosphorus bonding from analysis of the electronic structure of redox pairs of metal-phosphine complexes , 2005 .

[12]  Alexander D. MacKerell,et al.  Development and current status of the CHARMM force field for nucleic acids , 2000, Biopolymers.

[13]  I. H. Hillier,et al.  Combined quantum mechanical–molecular mechanical study of catalysis by the enzyme phospholipase A2: an investigation of the potential energy surface for amide hydrolysis , 1991 .

[14]  Donald G. Truhlar,et al.  Generalized Hybrid Orbital (GHO) Method for Combining Ab Initio Hartree-Fock Wave Functions with Molecular Mechanics , 2004 .

[15]  Iris Antes,et al.  Adjusted Connection Atoms for Combined Quantum Mechanical and Molecular Mechanical Methods , 1999 .

[16]  Adrian J Mulholland,et al.  Differential transition-state stabilization in enzyme catalysis: quantum chemical analysis of interactions in the chorismate mutase reaction and prediction of the optimal catalytic field. , 2004, Journal of the American Chemical Society.

[17]  M. Dulak,et al.  On the electron leak problem in orbital-free embedding calculations. , 2006, The Journal of chemical physics.

[18]  Jing Kong,et al.  Lennard–Jones parameters for the combined QM/MM method using the B3LYP/6‐31G*/AMBER potential , 2005, J. Comput. Chem..

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

[20]  U. Singh,et al.  A combined ab initio quantum mechanical and molecular mechanical method for carrying out simulations on complex molecular systems: Applications to the CH3Cl + Cl− exchange reaction and gas phase protonation of polyethers , 1986 .

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

[22]  A. Mulholland,et al.  The Fe-CO bond energy in myoglobin: a QM/MM study of the effect of tertiary structure. , 2006, Biophysical journal.

[23]  Donald G Truhlar,et al.  Generalized hybrid-orbital method for combining density functional theory with molecular mechanicals. , 2005, Chemphyschem : a European journal of chemical physics and physical chemistry.

[24]  F. Javier Luque,et al.  Energy decomposition in molecular complexes: Implications for the treatment of polarization in molecular simulations , 2003, J. Comput. Chem..

[25]  György G. Ferenczy,et al.  Quantum mechanical computations on very large molecular systems: The local self‐consistent field method , 1994, J. Comput. Chem..

[26]  L. Curtiss,et al.  Intermolecular interactions from a natural bond orbital, donor-acceptor viewpoint , 1988 .

[27]  O. Kapuy,et al.  Analysis of weakly bound structures: Hydrogen bond and the electron density in a water dimer , 2003 .

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

[29]  P. Geerlings,et al.  Conceptual density functional theory. , 2003, Chemical reviews.

[30]  Adrian J Mulholland,et al.  Modelling enzyme reaction mechanisms, specificity and catalysis. , 2005, Drug discovery today.

[31]  Kenneth M. Merz,et al.  Combined Quantum Mechanical/Molecular Mechanical Methodologies Applied to Biomolecular Systems , 1999 .

[32]  J. Ponder,et al.  Force fields for protein simulations. , 2003, Advances in protein chemistry.

[33]  William L. Jorgensen,et al.  Ab Initio Study Of Hydrogen-Bonded Complexes Of Small Organic Molecules With Water , 1998 .

[34]  Adrian J Mulholland,et al.  Multiple high-level QM/MM reaction paths demonstrate transition-state stabilization in chorismate mutase: correlation of barrier height with transition-state stabilization. , 2005, Chemical communications.

[35]  A. Mulholland,et al.  Acetyl‐CoA enolization in citrate synthase: A quantum mechanical/molecular mechanical (QM/MM) study , 1997, Proteins.

[36]  Arieh Warshel,et al.  Simulation of enzyme reactions using valence bond force fields and other hybrid quantum/classical approaches , 1993 .

[37]  Qiang Cui,et al.  Importance of van der Waals Interactions in QM/MM Simulations. , 2004, The journal of physical chemistry. B.

[38]  Kazuo Kitaura,et al.  A new energy decomposition scheme for molecular interactions within the Hartree‐Fock approximation , 1976 .

[39]  Nathalie Reuter,et al.  Frontier Bonds in QM/MM Methods: A Comparison of Different Approaches , 2000 .

[40]  Jiali Gao,et al.  Combined Quantum Mechanical and Molecular Mechanical Methods , 1999 .

[41]  Manuel F. Ruiz-López,et al.  Calibration of the Quantum/Classical Hamiltonian in Semiempirical QM/MM AM1 and PM3 Methods , 2000 .

[42]  Manuel F. Ruiz-López,et al.  An iterative procedure to determine Lennard-Jones parameters for their use in quantum mechanics/molecular mechanics liquid state simulations , 2002 .

[43]  Alessandro Laio,et al.  A Hamiltonian electrostatic coupling scheme for hybrid Car-Parrinello molecular dynamics simulations , 2002 .

[44]  Alexander D. MacKerell.  Chapter 7 Empirical Force Fields for Proteins: Current Status and Future Directions , 2005 .

[45]  Richard A. Friesner,et al.  A mixed quantum mechanics/molecular mechanics (QM/MM) method for large‐scale modeling of chemistry in protein environments , 2000, J. Comput. Chem..

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

[47]  Bernard R. Brooks,et al.  Optimization of quantum mechanical molecular mechanical partitioning schemes: Gaussian delocalization of molecular mechanical charges and the double link atom method , 2002 .

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

[49]  M. Field,et al.  A Generalized Hybrid Orbital (GHO) Method for the Treatment of Boundary Atoms in Combined QM/MM Calculations , 1998 .

[50]  D. Feil,et al.  A local density‐functional study of the electron density distribution in the H2O dimer , 1988 .

[51]  A. Rappé,et al.  Ab Initio Calculation of Nonbonded Interactions: Are We There Yet? , 2000 .

[52]  Jean-Louis Rivail,et al.  Quantum chemical computations on parts of large molecules: the ab initio local self consistent field method , 1996 .

[53]  K. Morokuma,et al.  Molecular orbital studies of hydrogen bonds , 1977 .

[54]  Christine M. Bathelt,et al.  QM/MM studies of the electronic structure of the compound I intermediate in cytochrome c peroxidase and ascorbate peroxidase. , 2005, Dalton transactions.

[55]  B. Silvi,et al.  Structure and vibrational properties of water hydrogen halide complexes , 1991 .

[56]  M. Swart.  AddRemove: A new link model for use in QM/MM studies , 2003 .

[57]  Adrian J Mulholland,et al.  Electronic structure of compound I in human isoforms of cytochrome P450 from QM/MM modeling. , 2005, Journal of the American Chemical Society.

[58]  Adrian J. Mulholland,et al.  Ab initio QM/MM study of the citrate synthase mechanism. A low-barrier hydrogen bond is not involved , 2000 .

[59]  J. Harvey,et al.  Spin-forbidden CO ligand recombination in myoglobin. , 2004, Faraday discussions.

[60]  Adrian J Mulholland,et al.  Chemical accuracy in QM/MM calculations on enzyme-catalysed reactions , 2007, Chemistry Central journal.

[61]  Victor Guallar,et al.  Peripheral heme substituents control the hydrogen-atom abstraction chemistry in cytochromes P450 , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[62]  V. Vasilyev,et al.  Tetrahedral intermediate formation in the acylation step of acetylcholinesterases. A combined quantum chemical and molecular mechanical model , 1994 .

[63]  Paul Tavan,et al.  A hybrid method for solutes in complex solvents: Density functional theory combined with empirical force fields , 1999 .