On possible pitfalls in ab initio quantum mechanics/molecular mechanics minimization approaches for studies of enzymatic reactions.

Reliable studies of enzymatic reactions by combined quantum mechanics/molecular mechanics (QM/MM) approaches, with an ab initio description of the quantum region, presents a major challenge to computational chemists. The main problem is the need for a very large computer time for the evaluation of the QM energy, which in turn makes it extremely challenging to perform proper configurational sampling. A seemingly reasonable alternative is to perform energy minimization studies of the type used in gas-phase ab initio studies. However, it is hard to see why such an approach should give reliable results in protein active sites. To examine the problems with energy minimization QM/MM approaches, we chose the hypothetical reaction of a metaphosphate ion with water in the Ras.GAP complex. This hypothetical reaction served as a simple benchmark reaction. The possible problems with the QM/MM minimization were explored by generating several protein configurations from long MD simulations and using energy minimization and scanning of the reaction coordinates to evaluate the corresponding potential energy surfaces of the reaction for each of these different protein configurations. Comparing these potential energy surfaces, we found major variations of the corresponding minima. Furthermore, the reaction energies and activation energies also varied significantly even for similar protein configurations. The specific coordination of a magnesium ion, present in the active center of the protein complex, turned out to influence the energetics of the reaction in a major way, where a direct coordination to the reactant leads to an increase of the activation energy by 17 kcal/mol. Apparently, using energy minimization to generate potential surfaces for an enzymatic reaction, while starting from a single protein structure, could lead to major errors in calculations of activation free energies and binding free energies. Thus we believe that extensive samplings of the configurational space of the protein are essential for meaningful determination of the energetics of enzymatic reactions. The possible relevance of our conclusion with regard to a recent study of the RasGAP reaction is discussed.

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

[2]  W. Kabsch,et al.  The Ras-RasGAP complex: structural basis for GTPase activation and its loss in oncogenic Ras mutants. , 1997, Science.

[3]  J. M. Lluch,et al.  A QM/MM study of the racemization of vinylglycolate catalyzed by mandelate racemase enzyme. , 2001, Journal of the American Chemical Society.

[4]  Arieh Warshel,et al.  A local reaction field method for fast evaluation of long‐range electrostatic interactions in molecular simulations , 1992 .

[5]  John A Pople.  Quantum Chemical Models (Nobel Lecture). , 1999, Angewandte Chemie.

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

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

[8]  Arieh Warshel,et al.  Computer simulations of protein functions: searching for the molecular origin of the replication fidelity of DNA polymerases. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

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

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

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

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

[13]  K. Gerwert,et al.  Theoretical IR spectroscopy based on QM/MM calculations provides changes in charge distribution, bond lengths, and bond angles of the GTP ligand induced by the Ras-protein. , 2005, Biophysical journal.

[14]  U. Ryde,et al.  Binding of Benzylpenicillin to Metallo-β-lactamase: A QM/MM Study , 2004 .

[15]  A. Klamt,et al.  Fast solvent screening via quantum chemistry: COSMO‐RS approach , 2002 .

[16]  Arieh Warshel,et al.  Calculations of ground- and excited-state potential surfaces for conjugated heteroatomic molecules , 1981 .

[17]  A. Warshel,et al.  Phosphate Ester Hydrolysis in Aqueous Solution: Associative versus Dissociative Mechanisms , 1998 .

[18]  A. Warshel,et al.  Constraining the electron densities in DFT method as an effective way for ab initio studies of metal‐catalyzed reactions , 2000 .

[19]  Arieh Warshel,et al.  Langevin Dipoles Model for ab Initio Calculations of Chemical Processes in Solution: Parametrization and Application to Hydration Free Energies of Neutral and Ionic Solutes and Conformational Analysis in Aqueous Solution , 1997 .

[20]  M. Rodnina,et al.  Essential role of histidine 84 in elongation factor Tu for the chemical step of GTP hydrolysis on the ribosome. , 2003, Journal of molecular biology.

[21]  Bernard R. Brooks,et al.  ENZYME MECHANISMS WITH HYBRID QUANTUM AND MOLECULAR MECHANICAL POTENTIALS.I. THEORETICAL CONSIDERATIONS , 1996 .

[22]  Arieh Warshel,et al.  Microscopic and semimicroscopic calculations of electrostatic energies in proteins by the POLARIS and ENZYMIX programs , 1993, J. Comput. Chem..

[23]  Arieh Warshel,et al.  Frozen density functional free energy simulations of redox proteins: computational studies of the reduction potential of plastocyanin and rusticyanin. , 2003, Journal of the American Chemical Society.

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

[25]  R A Friesner,et al.  Quantum mechanical calculations on biological systems. , 1998, Current opinion in structural biology.

[26]  R. Cachau,et al.  QM/MM modeling the Ras–GAP catalyzed hydrolysis of guanosine triphosphate , 2005, Proteins.

[27]  Martin J. Field,et al.  Simulating enzyme reactions: Challenges and perspectives , 2002, J. Comput. Chem..

[28]  Arieh Warshel,et al.  Computer Simulations of the Action of Metalloenzymes , 1997 .

[29]  B. Grigorenko,et al.  Modeling of serine protease prototype reactions with the flexible effective fragment potential quantum mechanical/molecular mechanical method , 2004 .

[30]  Paul Tavan,et al.  IR Spectra of Phosphate Ions in Aqueous Solution: Predictions of a DFT/MM Approach Compared with Observations , 2004 .

[31]  Arieh Warshel,et al.  Ab Initio QM/MM Simulation with Proper Sampling: “First Principle” Calculations of the Free Energy of the Autodissociation of Water in Aqueous Solution , 2002 .

[32]  Arieh Warshel,et al.  Converting conformational changes to electrostatic energy in molecular motors: The energetics of ATP synthase , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[33]  Efthimios Kaxiras,et al.  A QM/MM Implementation of the Self-Consistent Charge Density Functional Tight Binding (SCC-DFTB) Method , 2001 .

[34]  A. Warshel,et al.  Why does the Ras switch “break” by oncogenic mutations? , 2004, Proteins.

[35]  A. Warshel,et al.  Structure/function correlations of proteins using MM, QM/MM, and related approaches: methods, concepts, pitfalls, and current progress. , 2003, Advances in protein chemistry.

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

[37]  R. Cachau,et al.  Quantum chemical modeling of the GTP hydrolysis by the RAS-GAP protein complex. , 2004, Biochimica et biophysica acta.

[38]  H. Bernhard Schlegel,et al.  Reaction Path Following in Mass-Weighted Internal Coordinates , 1990 .

[39]  J. Bertrán,et al.  Transition structure selectivity in enzyme catalysis: a QM/MM study of chorismate mutase , 2001 .

[40]  Arieh Warshel,et al.  Ab Initio Calculations of Free Energy Barriers for Chemical Reactions in Solution , 1995 .

[41]  H. Schlegel,et al.  A combined method for determining reaction paths, minima, and transition state geometries , 1997 .