Lennard–Jones parameters for the combined QM/MM method using the B3LYP/6‐31G*/AMBER potential

A combined DFT quantum mechanical and AMBER molecular mechanical potential (QM/MM) is presented for use in molecular modeling and molecular simulations of large biological systems. In our approach we evaluate Lennard–Jones parameters describing the interaction between the quantum mechanical (QM) part of a system, which is described at the B3LYP/6‐31+G* level of theory, and the molecular mechanical (MM) part of the system, described by the AMBER force field. The Lennard–Jones parameters for this potential are obtained by calculating hydrogen bond energies and hydrogen bond geometries for a large set of bimolecular systems, in which one hydrogen bond monomer is described quantum mechanically and the other is treated molecular mechanically. We have investigated more than 100 different bimolecular systems, finding very good agreement between hydrogen bond energies and geometries obtained from the combined QM/MM calculations and results obtained at the QM level of theory, especially with respect to geometry. Therefore, based on the Lennard–Jones parameters obtained in our study, we anticipate that the B3LYP/6‐31+G*/AMBER potential will be a precise tool to explore intermolecular interactions inside a protein environment. © 2005 Wiley Periodicals, Inc. J Comput Chem 26: 1270–1278, 2005

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

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

[3]  Y. Mo,et al.  Ab initio QM/MM simulations with a molecular orbital‐valence bond (MOVB) method: application to an SN2 reaction in water , 2000 .

[4]  Kenneth M. Merz,et al.  Calculation of solvation free energies using a density functional/molecular dynamics coupled potential , 1993 .

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

[6]  W. Thiel,et al.  Semiempirical treatment of electrostatic potentials and partial charges in combined quantum mechanical and molecular mechanical approaches , 1996, J. Comput. Chem..

[7]  Claude Millot,et al.  A coupled density functional‐molecular mechanics Monte Carlo simulation method: The water molecule in liquid water , 1996 .

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

[9]  Sason Shaik,et al.  The elusive oxidant species of cytochrome P450 enzymes: characterization by combined quantum mechanical/molecular mechanical (QM/MM) calculations. , 2002, Journal of the American Chemical Society.

[10]  P. Kozlowski.  Quantum chemical modeling of Co--C bond activation in B(12)-dependent enzymes. , 2001, Current opinion in chemical biology.

[11]  D. Truhlar,et al.  Quantum Mechanical Dynamical Effects in an Enzyme-Catalyzed Proton Transfer Reaction , 1999 .

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

[13]  Kenneth M. Merz,et al.  Dynamic Force Field Models: Molecular Dynamics Simulations of Human Carbonic Anhydrase II Using a Quantum Mechanical/Molecular Mechanical Coupled Potential , 1995 .

[14]  F. Grozema,et al.  Combined Quantum Mechanical and Molecular Mechanical Methods , 1999 .

[15]  D. Hartsough,et al.  Potential of mean force calculations on the SN1 fragmentation of tert-butyl chloride , 1995 .

[16]  Jiali Gao,et al.  Computation of Intermolecular Interactions with a Combined Quantum Mechanical and Classical Approach , 1994 .

[17]  B. Brooks,et al.  HIV-1 Protease Cleavage Mechanism Elucidated with Molecular Dynamics Simulation , 1995 .

[18]  Jiali Gao,et al.  Optimization of the Lennard‐Jones parameters for a combined ab initio quantum mechanical and molecular mechanical potential using the 3‐21G basis set , 1996 .

[19]  Ramkumar Rajamani,et al.  Combined QM/MM study of the opsin shift in bacteriorhodopsin , 2002, J. Comput. Chem..

[20]  David Martin,et al.  Computational Molecular Biology: An Algorithmic Approach , 2001 .

[21]  Feliu Maseras,et al.  IMOMM: A new integrated ab initio + molecular mechanics geometry optimization scheme of equilibrium structures and transition states , 1995, J. Comput. Chem..

[22]  Martin J. Field,et al.  Free energy perturbation method for chemical reactions in the condensed phase: a dynamic approach based on a combined quantum and molecular mechanics potential , 1987 .

[23]  Sason Shaik,et al.  Quantum mechanical/molecular mechanical investigation of the mechanism of C-H hydroxylation of camphor by cytochrome P450cam: theory supports a two-state rebound mechanism. , 2004, Journal of the American Chemical Society.

[24]  R. Cramer,et al.  Validation of the general purpose tripos 5.2 force field , 1989 .

[25]  A. Lesk.  COMPUTATIONAL MOLECULAR BIOLOGY , 1988, Proceeding of Data For Discovery.

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

[27]  Martin J. Field,et al.  Chapter 1 – Hybrid Potentials for Large Molecular Systems , 1999 .

[28]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[29]  R. Ornstein,et al.  A combined density functional theory/molecular mechanics formalism and its application to small water clusters , 1997 .

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

[31]  Douglas A. Smith.  Modeling the hydrogen bond , 1994 .

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

[33]  Jon Baker,et al.  Q‐Chem 2.0: a high‐performance ab initio electronic structure program package , 2000, J. Comput. Chem..

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

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

[36]  W. L. Jorgensen,et al.  The OPLS [optimized potentials for liquid simulations] potential functions for proteins, energy minimizations for crystals of cyclic peptides and crambin. , 1988, Journal of the American Chemical Society.

[37]  Jean-Louis Rivail,et al.  Specific force field parameters determination for the hybrid ab initio QM/MM LSCF method , 2002, J. Comput. Chem..

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

[39]  Richard A. Friesner,et al.  Mixed ab initio QM/MM modeling using frozen orbitals and tests with alanine dipeptide and tetrapeptide , 1999 .

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

[41]  U. Rothlisberger,et al.  The Role of p-p Stacking Interactions in Square Planar Palladium Complexes. Combined Quantum Mechanics/Molecular Mechanics QM/MM Studies , 2001 .

[42]  H Liu,et al.  A combined quantum/classical molecular dynamics study of the catalytic mechanism of HIV protease. , 1996, Journal of molecular biology.

[43]  M. Freindorf,et al.  A combined density functional theory and molecular mechanics study of the relationship between the structure of coenzyme B12 and its binding to methylmalonyl-CoA mutase. , 2004, Journal of the American Chemical Society.

[44]  Jiali Gao,et al.  Walden-Inversion-Enforced Transition-State Stabilization in a Protein Tyrosine Phosphatase , 1998 .

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

[46]  J. Gao,et al.  A priori evaluation of aqueous polarization effects through Monte Carlo QM-MM simulations. , 1992, Science.

[47]  J. Rivail,et al.  The Local Self-Consistent Field Principles and Applications to Combined Quantum Mechanical-Molecular Mechanical Computations on Biomacromolecular Systems , 1998 .

[48]  Jiali Gao,et al.  Hybrid Quantum and Molecular Mechanical Simulations: An Alternative Avenue to Solvent Effects in Organic Chemistry , 1996 .