On the interpretation of the observed linear free energy relationship in phosphate hydrolysis: a thorough computational study of phosphate diester hydrolysis in solution.

The hydrolysis of phosphate esters is crucially important to biological systems, being involved in, among other things, signaling, energy transduction, biosynthesis, and the regulation of protein function. Despite this, there are many questions that remain unanswered in this important field, particularly with regard to the preferred mechanism of hydrolysis of phosphate esters, which can proceed through any of multiple pathways that are either associative or dissociative in nature. Previous comparisons of calculated and observed linear free energy relationships (LFERs) for phosphate monoester dianions with different leaving groups showed that the TS character gradually changes from associative to dissociative with the increasing acidity of the leaving group, while reproducing the experimental LFER. Here, we have generated ab initio potential energy surfaces for the hydrolysis of phosphate diesters in solution, with a variety of leaving groups. Once again, the reaction changes from a compact concerted pathway to one that is more expansive in character when the acidity of the leaving group increases. When such systems are examined in solution, it is essential to take into consideration the contribution of solute to the overall activation entropy, which remains a major computational challenge. The popular method of calculating the entropy using a quasi-harmonic approximation appears to markedly overestimate the configurational entropy for systems with multiple occupied energy wells. We introduce an improved restraint release approach for evaluating configurational entropies and apply this approach to our systems. We demonstrate that when this factor is taken into account, it is possible to reproduce the experimental LFER for this system with reasonable accuracy.

[1]  D. Herschlag,et al.  Do electrostatic interactions with positively charged active site groups tighten the transition state for enzymatic phosphoryl transfer? , 2004, Journal of the American Chemical Society.

[2]  W. Jencks.  Catalysis in chemistry and enzymology , 1969 .

[3]  R. Cachau,et al.  Computational study of a transition state analog of phosphoryl transfer in the Ras–RasGAP complex: AlFx versus MgF3– , 2005, Journal of molecular modeling.

[4]  S. Sprang,et al.  Kinetic isotope effects in Ras-catalyzed GTP hydrolysis: Evidence for a loose transition state , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[5]  Arieh Warshel,et al.  Computer simulation of the chemical catalysis of DNA polymerases: discriminating between alternative nucleotide insertion mechanisms for T7 DNA polymerase. , 2003, Journal of the American Chemical Society.

[6]  A. Warshel,et al.  Mechanistic alternatives in phosphate monoester hydrolysis: what conclusions can be drawn from available experimental data? , 1999, Chemistry & biology.

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

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

[9]  Richard A. Friesner,et al.  Quasi-harmonic method for calculating vibrational spectra from classical simulations on multi-dimensional anharmonic potential surfaces , 1984 .

[10]  A. Warshel,et al.  How important are entropic contributions to enzyme catalysis? , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[11]  I R Vetter,et al.  Nucleoside triphosphate-binding proteins: different scaffolds to achieve phosphoryl transfer , 1999, Quarterly Reviews of Biophysics.

[12]  P. Reinemer,et al.  Crystal structure of the catalytic subunit of human protein phosphatase 1 and its complex with tungstate. , 1995, Journal of molecular biology.

[13]  R. Levy,et al.  Corrections to the quasiharmonic approximation for evaluating molecular entropies , 1986 .

[14]  Arieh Warshel,et al.  Calculations of Activation Entropies of Chemical Reactions in Solution , 2000 .

[15]  O. Ramsay,et al.  Mechanisms of Nucleophilic Substitution in Phosphate Esters , 1964 .

[16]  Peter Pulay,et al.  Geometry optimization by direct inversion in the iterative subspace , 1984 .

[17]  A. Warshel,et al.  Why have mutagenesis studies not located the general base in ras p21 , 1994, Nature Structural Biology.

[18]  K. Taira,et al.  The Hydrolysis of RNA: From Theoretical Calculations to the Hammerhead Ribozyme-Mediated Cleavage of RNA. , 1998, Chemical reviews.

[19]  Jacopo Tomasi,et al.  A new integral equation formalism for the polarizable continuum model: Theoretical background and applications to isotropic and anisotropic dielectrics , 1997 .

[20]  Michael K Gilson,et al.  Evaluating the Accuracy of the Quasiharmonic Approximation. , 2005, Journal of chemical theory and computation.

[21]  D. Herschlag,et al.  Alkaline phosphatase mono- and diesterase reactions: comparative transition state analysis. , 2006, Journal of the American Chemical Society.

[22]  D. Herschlag,et al.  Ras-catalyzed hydrolysis of GTP: a new perspective from model studies. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[23]  A. Warshel,et al.  On the Reactivity of Phosphate Monoester Dianions in Aqueous Solution: Brønsted Linear Free-Energy Relationships Do Not Have an Unique Mechanistic Interpretation , 1998 .

[24]  D. Herschlag,et al.  Mapping the transition state for ATP hydrolysis: implications for enzymatic catalysis. , 1995, Chemistry & biology.

[25]  E. E. Kim,et al.  Reaction mechanism of alkaline phosphatase based on crystal structures. Two-metal ion catalysis. , 1991, Journal of molecular biology.

[26]  A. Hengge,et al.  Examination of P-OR bridging bond orders in phosphate monoesters using (18)O isotope shifts in 31P NMR. , 2005, The Journal of organic chemistry.

[27]  Lattice dynamics of solid cubane within the quasi-harmonic approximation , 2002, cond-mat/0301564.

[28]  D. Herschlag,et al.  The nature of the transition state for enzyme-catalyzed phosphoryl transfer. Hydrolysis of O-aryl phosphorothioates by alkaline phosphatase. , 1995, Biochemistry.

[29]  J. Springer,et al.  Structural analysis of inositol monophosphatase complexes with substrates. , 1994, Biochemistry.

[30]  D. Herschlag,et al.  Kinetic isotope effects for alkaline phosphatase reactions: implications for the role of active-site metal ions in catalysis. , 2007, Journal of the American Chemical Society.

[31]  K. Sharp,et al.  On the calculation of absolute macromolecular binding free energies , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[32]  Philippe Y. Ayala,et al.  Identification and treatment of internal rotation in normal mode vibrational analysis , 1998 .

[33]  D. Herschlag,et al.  Functional interrelationships in the alkaline phosphatase superfamily: phosphodiesterase activity of Escherichia coli alkaline phosphatase. , 2001, Biochemistry.

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

[35]  R. Goody,et al.  The mechanism of guanosine nucleotide hydrolysis by p21 c-Ha-ras. The stereochemical course of the GTPase reaction. , 1989, The Journal of biological chemistry.

[36]  V. Barone,et al.  Quantum Calculation of Molecular Energies and Energy Gradients in Solution by a Conductor Solvent Model , 1998 .

[37]  A. Warshel,et al.  On the mechanism of guanosine triphosphate hydrolysis in ras p21 proteins. , 1992, Biochemistry.

[38]  Richard H. Henchman,et al.  Revisiting free energy calculations: a theoretical connection to MM/PBSA and direct calculation of the association free energy. , 2004, Biophysical journal.

[39]  Lu Wang,et al.  Inclusion of Loss of Translational and Rotational Freedom in Theoretical Estimates of Free Energies of Binding. Application to a Complex of Benzene and Mutant T4 Lysozyme , 1997 .

[40]  Samuel H. Wilson,et al.  Energy analysis of chemistry for correct insertion by DNA polymerase β , 2006, Proceedings of the National Academy of Sciences.

[41]  A. Warshel,et al.  How does GAP catalyze the GTPase reaction of Ras? A computer simulation study. , 2000, Biochemistry.

[42]  Jacopo Tomasi,et al.  Evaluation of Solvent Effects in Isotropic and Anisotropic Dielectrics and in Ionic Solutions with a Unified Integral Equation Method: Theoretical Bases, Computational Implementation, and Numerical Applications , 1997 .

[43]  M. Karplus,et al.  Phosphate ester hydrolysis: calculation of gas-phase reaction paths and solvation effects , 1994 .

[44]  H. Kalbitzer,et al.  Substrate-assisted catalysis as a mechanism for GTP hydrolysis of p21ras and other GTP-binding proteins , 1995, Nature Structural Biology.

[45]  Michael Y. Galperin,et al.  A superfamily of metalloenzymes unifies phosphopentomutase and cofactor‐independent phosphoglycerate mutase with alkaline phosphatases and sulfatases , 1998, Protein science : a publication of the Protein Society.

[46]  Giovanni Scalmani,et al.  New developments in the polarizable continuum model for quantum mechanical and classical calculations on molecules in solution , 2002 .

[47]  Miquel Duran,et al.  How does basis set superposition error change the potential surfaces for hydrogen-bonded dimers? , 1996 .

[48]  M Geyer,et al.  Linear free energy relationships in the intrinsic and GTPase activating protein-stimulated guanosine 5'-triphosphate hydrolysis of p21ras. , 1996, Biochemistry.

[49]  M. Karplus,et al.  Evaluation of the configurational entropy for proteins: application to molecular dynamics simulations of an α-helix , 1984 .

[50]  M. Ruiz‐López,et al.  Theoretical evaluation of the substrate-assisted catalysis mechanism for the hydrolysis of phosphate monoester dianions. , 2007, Chemistry.

[51]  D. Herschlag,et al.  Phosphoryl Transfer to Anionic Oxygen Nucleophiles. Nature of the Transition State and Electrostatic Repulsion , 1989 .

[52]  J. Wilkie,et al.  Metal ions in the mechanism of enzyme-catalysed phosphate monoester hydrolyses , 1997 .

[53]  T. Brinck,et al.  Theoretical Studies of the Hydrolysis of the Methyl Phosphate Anion , 1999 .

[54]  S. Benkovic,et al.  6 Chemical Basis of Biological Phosphoryl Transfer , 1973 .

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

[56]  Hua Guo,et al.  Supermolecule density functional calculations suggest a key role for solvent in alkaline hydrolysis of p-nitrophenyl phosphate. , 2007, Chemical communications.

[57]  A. Warshel,et al.  What are the roles of substrate-assisted catalysis and proximity effects in peptide bond formation by the ribosome? , 2005, Biochemistry.

[58]  G. Maga,et al.  Eukaryotic DNA polymerases. , 2002, Annual review of biochemistry.

[59]  W. Cleland,et al.  Enzymatic Mechanisms of Phosphate and Sulfate Transfer , 2006 .

[60]  J. Liu,et al.  Insight into the catalytic mechanism of DNA polymerase beta: structures of intermediate complexes. , 2001, Biochemistry.

[61]  R. Wolfenden,et al.  The time required for water attack at the phosphorus atom of simple phosphodiesters and of DNA. , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[62]  M. Bollen,et al.  Structural and Catalytic Similarities between Nucleotide Pyrophosphatases/Phosphodiesterases and Alkaline Phosphatases* , 2001, The Journal of Biological Chemistry.

[63]  William P. Jencks,et al.  A primer for the Bema Hapothle. An empirical approach to the characterization of changing transition-state structures , 1985 .

[64]  J. Wilkie,et al.  Comparison of inline and non-inline associative and dissociative reaction pathways for model reactions of phosphate monoester hydrolysis , 1996 .

[65]  Samuel H. Wilson,et al.  Structures of ternary complexes of rat DNA polymerase beta, a DNA template-primer, and ddCTP. , 1994, Science.

[66]  Shina C L Kamerlin,et al.  The role of metal ions in phosphate ester hydrolysis. , 2007, Organic & biomolecular chemistry.

[67]  A. Becke.  Density-functional thermochemistry. III. The role of exact exchange , 1993 .

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

[69]  Nicholas H Williams,et al.  Reactivity of Phosphate Diesters Doubly Coordinated to a Dinuclear Cobalt(III) Complex: Dependence of the Reactivity on the Basicity of the Leaving Group , 1998 .

[70]  H. Bernhard Schlegel,et al.  Methods for optimizing large molecules. II. Quadratic search , 1999 .

[71]  A. Mildvan.  The role of metals in enzyme-catalyzed substitutions at each of the phosphorus atoms of ATP. , 1979, Advances in enzymology and related areas of molecular biology.

[72]  D. Herschlag,et al.  The Substrate-Assisted General Base Catalysis Model for Phosphate Monoester Hydrolysis: Evaluation Using Reactivity Comparisons , 2000 .

[73]  A. Warshel,et al.  A Fundamental Assumption about OH- Attack in Phosphate Ester Hydrolysis Is Not Fully Justified , 1997 .

[74]  J. Wilkie,et al.  Stereochemical, mechanistic, and structural features of enzyme-catalysed phosphate monoester hydrolyses , 1995 .

[75]  A. Warshel,et al.  Mechanistic analysis of the observed linear free energy relationships in p21ras and related systems. , 1996, Biochemistry.

[76]  R. A. More O'Ferrall,et al.  Relationships between E2 and E1cB mechanisms of β-elimination , 1970 .

[77]  M. Karplus,et al.  Method for estimating the configurational entropy of macromolecules , 1981 .

[78]  Arieh Warshel,et al.  On the mechanism of hydrolysis of phosphate monoesters dianions in solutions and proteins. , 2006, Journal of the American Chemical Society.

[79]  D. Rinaldi,et al.  Theoretical studies of the hydroxide-catalyzed P-O cleavage reactions of neutral phosphate triesters and diesters in aqueous solution: examination of the changes induced by H/Me substitution. , 2005, The journal of physical chemistry. B.

[80]  J. Coleman,et al.  Structure and mechanism of alkaline phosphatase. , 1992, Annual review of biophysics and biomolecular structure.

[81]  Suse Broyde,et al.  A water-mediated and substrate-assisted catalytic mechanism for Sulfolobus solfataricus DNA polymerase IV. , 2007, Journal of the American Chemical Society.

[82]  T. N. Bhat,et al.  The Protein Data Bank , 2000, Nucleic Acids Res..

[83]  Andrew Williams.  Effective charge and Leffler's index as mechanistic tools for reactions in solution , 1984 .

[84]  Michael Y. Galperin,et al.  Conserved core structure and active site residues in alkaline phosphatase superfamily enzymes , 2001, Proteins.

[85]  J. Tomasi,et al.  The IEF version of the PCM solvation method: an overview of a new method addressed to study molecular solutes at the QM ab initio level , 1999 .

[86]  A. Klamt,et al.  COSMO : a new approach to dielectric screening in solvents with explicit expressions for the screening energy and its gradient , 1993 .

[87]  S. F. Boys,et al.  The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors , 1970 .

[88]  Arieh Warshel,et al.  Computer Modeling of Chemical Reactions in Enzymes and Solutions , 1991 .

[89]  Jacopo Tomasi,et al.  Continuum solvation models: A new approach to the problem of solute’s charge distribution and cavity boundaries , 1997 .

[90]  A. Warshel,et al.  Simulating the effect of DNA polymerase mutations on transition-state energetics and fidelity: evaluating amino acid group contribution and allosteric coupling for ionized residues in human pol beta. , 2006, Biochemistry.

[91]  W. Jencks.  General acid-base catalysis of complex reactions in water , 1972 .

[92]  J E Leffler,et al.  Parameters for the Description of Transition States. , 1953, Science.

[93]  P. Schultz,et al.  Probing the structure and mechanism of Ras protein with an expanded genetic code. , 1993, Science.

[94]  N. Allan,et al.  Ionic solids at elevated temperatures and/or high pressures: lattice dynamics, molecular dynamics, Monte Carlo and ab initio studies , 2000 .