Revised AMBER parameters for bioorganic phosphates.

We report AMBER force field parameters for biological simulations involving phosphorylation of serine, threonine or tyrosine. The initial parameters used RESP fitting for the atomic partial charges and standard values for all other parameters such as Lennard-Jones coefficients. These were refined with the aid of a thermodynamic cycle consisting of experimentally determined pKa values, solvation energies from molecular dynamics free energy simulations, and gas phase basicities from QM calculations. A polarization energy term was included to account for the charge density change between the gas-phase and solution, and solvation free energies were determined using thermodynamic integration. Parameter adjustment is required to obtain consistent thermodynamic results with better balanced electrostatic interactions between water and the phosphate oxygens. To achieve this we modified the phosphate oxygen radii. A thermodynamically consistent parameter set can be derived for monoanions and requires an increase of the van der Waals phosphate oxygen radii of approximately 0.09 Å. Larger, residue-specific radii appear to be needed for dianions. The revised parameters developed here should be of particular interest for environments where simulations of multiple protonation states may be of interest.

[1]  P. Kollman,et al.  A well-behaved electrostatic potential-based method using charge restraints for deriving atomic char , 1993 .

[2]  Michiel Sprik,et al.  Solvent polarization and hydration of the chlorine anion , 1990 .

[3]  Ralph G. Pearson,et al.  Ionization potentials and electron affinities in aqueous solution , 1986 .

[4]  D. R. Llewellyn,et al.  716. The reactions of organic phosphates. Part I. The hydrolysis of methyl dihydrogen phosphate , 1958 .

[5]  B. Roux,et al.  Absolute hydration free energy scale for alkali and halide ions established from simulations with a polarizable force field. , 2006, The journal of physical chemistry. B.

[6]  J. Eiler,et al.  The Acid Strength of Mono and Diesters of Phosphoric Acid. The n-Alkyl Esters from Methyl to Butyl, the Esters of Biological Importance, and the Natural Guanidine Phosphoric Acids , 1943 .

[7]  Alexander D. MacKerell,et al.  All‐atom empirical force field for nucleic acids: II. Application to molecular dynamics simulations of DNA and RNA in solution , 2000 .

[8]  Yizhak Marcus,et al.  Thermodynamics of solvation of ions. Part 5.—Gibbs free energy of hydration at 298.15 K , 1991 .

[9]  M. Trachtman,et al.  "Squiggle-H2O". An enquiry into the importance of solvation effects in phosphate ester and anhydride reactions. , 1970, Biochimica et biophysica acta.

[10]  Richard M. Noyes,et al.  Thermodynamics of Ion Hydration as a Measure of Effective Dielectric Properties of Water , 1962 .

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

[12]  D. York,et al.  Benchmark calculations of proton affinities and gas-phase basicities of molecules important in the study of biological phosphoryl transfer. , 2005, Physical chemistry chemical physics : PCCP.

[13]  J. Chanley,et al.  The Mechanism of the Hydrolysis of Organic Phosphates. III. Aromatic Phosphates1 , 1955 .

[14]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[15]  J. Max,et al.  Titrage infrarouge de l'acide phosphorique , 2000 .

[16]  R. Wolfenden,et al.  AFFINITIES OF PHOSPHORIC ACIDS, ESTERS, AND AMIDES FOR SOLVENT WATER. , 1983 .

[17]  Anselm H. C. Horn,et al.  AMBER force-field parameters for phosphorylated amino acids in different protonation states: phosphoserine, phosphothreonine, phosphotyrosine, and phosphohistidine , 2006, Journal of molecular modeling.

[18]  B. Neel,et al.  From Form to Function: Signaling by Protein Tyrosine Phosphatases , 1996, Cell.

[19]  J. Antosiewicz,et al.  Prediction of secondary ionization of the phosphate group in phosphotyrosine peptides. , 2003, Biophysical journal.

[20]  N. Ahn.  Introduction: Protein Phosphorylation and Signaling , 2001 .

[21]  T. Gojobori,et al.  Cancer‐related mutations in BRCA1‐BRCT cause long‐range structural changes in protein–protein binding sites: A molecular dynamics study , 2006, Proteins.

[22]  S. Pantano,et al.  The role of phosphorylation on the structure and dynamics of phospholamban: A model from molecular simulations , 2006, Proteins.

[23]  L. Curtiss,et al.  Gaussian-3 (G3) theory for molecules containing first and second-row atoms , 1998 .

[24]  T. Cheatham,et al.  Determination of Alkali and Halide Monovalent Ion Parameters for Use in Explicitly Solvated Biomolecular Simulations , 2008, The journal of physical chemistry. B.

[25]  T. Straatsma,et al.  THE MISSING TERM IN EFFECTIVE PAIR POTENTIALS , 1987 .

[26]  T. Hunter,et al.  Transcriptional control by protein phosphorylation: signal transmission from the cell surface to the nucleus , 1995, Current Biology.

[27]  Hans W. Horn,et al.  Accounting for polarization cost when using fixed charge force fields. I. Method for computing energy. , 2010, The journal of physical chemistry. B.

[28]  Tony Pawson,et al.  Specificity in Signal Transduction From Phosphotyrosine-SH2 Domain Interactions to Complex Cellular Systems , 2004, Cell.

[29]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[30]  Thomas Steinbrecher,et al.  A multistep approach to structure-based drug design: studying ligand binding at the human neutrophil elastase. , 2006, Journal of medicinal chemistry.

[31]  Alexander D. MacKerell,et al.  All‐atom empirical force field for nucleic acids: I. Parameter optimization based on small molecule and condensed phase macromolecular target data , 2000 .

[32]  David A. Case,et al.  Soft‐core potentials in thermodynamic integration: Comparing one‐ and two‐step transformations , 2011, J. Comput. Chem..

[33]  M. Tissandier,et al.  The Proton's Absolute Aqueous Enthalpy and Gibbs Free Energy of Solvation from Cluster-Ion Solvation Data , 1998 .

[34]  Gregory D. Hawkins,et al.  Extension of the platform of applicability of the SM5.42R universal solvation model , 1999 .

[35]  Shulin Zhuang,et al.  Insights into unbinding mechanisms upon two mutations investigated by molecular dynamics study of GSK3β–axin complex: Role of packing hydrophobic residues , 2007, Proteins.

[36]  Pengyu Y. Ren,et al.  Ion solvation thermodynamics from simulation with a polarizable force field. , 2003, Journal of the American Chemical Society.

[37]  Tony Hunter,et al.  The regulation of transcription by phosphorylation , 1992, Cell.