Strike a balance: Optimization of backbone torsion parameters of AMBER polarizable force field for simulations of proteins and peptides

Based on the AMBER polarizable model (ff02), we have reoptimized the parameters related to the main‐chain (Φ, Ψ) torsion angles by fitting to the Boltzmann‐weighted average quantum mechanical (QM) energies of the important regions (i.e., β, PII, αR, and αL regions). Following the naming convention of the AMBER force field series, this release will be called ff02pol.rl The force field has been assessed both by energetic comparison against the QM data and by the replica exchange molecular dynamics simulations of short alanine peptides in water. For Ace‐Ala‐Nme, the simulated populations in the β, PII and αR regions were approximately 30, 43, and 26%, respectively. For Ace‐(Ala)7‐Nme, the populations in these three regions were approximately 24, 49, and 26%. Both were in qualitative agreement with the NMR and CD experimental conclusions. In comparison with the previous force field, ff02pol.rl demonstrated good balance among these three important regions. The optimized torsion parameters, together with those in ff02, allow us to carry out simulations on proteins and peptides with the consideration of polarization. © 2006 Wiley Periodicals, Inc. J Comput Chem 27: 781–790, 2006

[1]  Ian W. Davis,et al.  Structure validation by Cα geometry: ϕ,ψ and Cβ deviation , 2003, Proteins.

[2]  K. Sanbonmatsu,et al.  α-Helical stabilization by side chain shielding of backbone hydrogen bonds , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[3]  Peter A. Kollman,et al.  Ion solvation in polarizable water: molecular dynamics simulations , 1991 .

[4]  Yong Duan,et al.  Computational protein folding: From lattice to all-atom , 2001, IBM Syst. J..

[5]  Charles L Brooks,et al.  Protein and peptide folding explored with molecular simulations. , 2002, Accounts of chemical research.

[6]  Robert W Woody,et al.  Solvent dependence of PII conformation in model alanine peptides. , 2004, Journal of the American Chemical Society.

[7]  Xiaolin Cao,et al.  Tripeptides adopt stable structures in water. A combined polarized visible Raman, FTIR, and VCD spectroscopy study. , 2002, Journal of the American Chemical Society.

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

[9]  Peter A. Kollman,et al.  Implementation of nonadditive intermolecular potentials by use of molecular dynamics: development of a water-water potential and water-ion cluster interactions , 1990 .

[10]  P. Kollman,et al.  Structure and Properties of Neat Liquids Using Nonadditive Molecular Dynamics: Water, Methanol, and N-Methylacetamide , 1995 .

[11]  Alexander D. MacKerell,et al.  Computational Biochemistry and Biophysics , 2001 .

[12]  Charles L. Brooks,et al.  CHARMM fluctuating charge force field for proteins: I parameterization and application to bulk organic liquid simulations , 2004, J. Comput. Chem..

[13]  Harry A. Stern,et al.  Fluctuating Charge, Polarizable Dipole, and Combined Models: Parameterization from ab Initio Quantum Chemistry , 1999 .

[14]  Pengyu Y. Ren,et al.  Consistent treatment of inter‐ and intramolecular polarization in molecular mechanics calculations , 2002, J. Comput. Chem..

[15]  Michael L. McKee,et al.  Theoretical investigation of the thermal ring opening of bicyclobutane to butadiene. Evidence for a nonsynchronous process , 1988 .

[16]  C. R. Mann.  The Theory of Optics , 1903, Nature.

[17]  Alexander D. MacKerell,et al.  Improved treatment of the protein backbone in empirical force fields. , 2004, Journal of the American Chemical Society.

[18]  Peter Hamm,et al.  Structure Determination of Trialanine in Water Using Polarization Sensitive Two-Dimensional Vibrational Spectroscopy , 2000 .

[19]  B. Berne,et al.  The free energy landscape for β hairpin folding in explicit water , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[20]  R. Friesner,et al.  Evaluation and Reparametrization of the OPLS-AA Force Field for Proteins via Comparison with Accurate Quantum Chemical Calculations on Peptides† , 2001 .

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

[22]  K. Kopple,et al.  Solvent-dependent conformational distributions of some dipeptides , 1980 .

[23]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[24]  Yoshitake Sakae,et al.  Optimization of protein force-field parameters with the Protein Data Bank , 2003, cond-mat/0309110.

[25]  Yong Duan,et al.  Solvation effects on alanine dipeptide: A MP2/cc‐pVTZ//MP2/6‐31G** study of (Φ, Ψ) energy maps and conformers in the gas phase, ether, and water , 2004, J. Comput. Chem..

[26]  Harry A. Stern,et al.  Development of a polarizable force field for proteins via ab initio quantum chemistry: First generation model and gas phase tests , 2002, J. Comput. Chem..

[27]  Chris Oostenbrink,et al.  A biomolecular force field based on the free enthalpy of hydration and solvation: The GROMOS force‐field parameter sets 53A5 and 53A6 , 2004, J. Comput. Chem..

[28]  Peter A. Kollman,et al.  Cation-.pi. Interactions: Nonadditive Effects Are Critical in Their Accurate Representation , 1995 .

[29]  William L. Jorgensen,et al.  OPLS all‐atom force field for carbohydrates , 1997 .

[30]  U. Singh,et al.  A NEW FORCE FIELD FOR MOLECULAR MECHANICAL SIMULATION OF NUCLEIC ACIDS AND PROTEINS , 1984 .

[31]  Karl-Heinz Ott,et al.  Parametrization of GROMOS force field for oligosaccharides and assessment of efficiency of molecular dynamics simulations , 1996, J. Comput. Chem..

[32]  Y. Sugita,et al.  Replica-exchange molecular dynamics method for protein folding , 1999 .

[33]  Jenn-Huei Lii,et al.  An improved force field (MM4) for saturated hydrocarbons , 1996, J. Comput. Chem..

[34]  P. Kollman,et al.  Pathways to a protein folding intermediate observed in a 1-microsecond simulation in aqueous solution. , 1998, Science.

[35]  Alexander D. MacKerell,et al.  Extending the treatment of backbone energetics in protein force fields: Limitations of gas‐phase quantum mechanics in reproducing protein conformational distributions in molecular dynamics simulations , 2004, J. Comput. Chem..

[36]  W. L. Jorgensen,et al.  Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids , 1996 .

[37]  Vijay S. Pande,et al.  Empirical force‐field assessment: The interplay between backbone torsions and noncovalent term scaling , 2005, J. Comput. Chem..

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

[39]  Norman L. Allinger,et al.  Molecular mechanics. The MM3 force field for hydrocarbons. 1 , 1989 .

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

[41]  Scott R. Wilson,et al.  Hyperconjugative effects on carbon—Carbon bond lengths in molecular mechanics (MM4) , 1996 .

[42]  P. Hamm,et al.  Subpicosecond conformational dynamics of small peptides probed by two-dimensional vibrational spectroscopy , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[43]  Chun Wu,et al.  Convergence of replica exchange molecular dynamics. , 2005, The Journal of chemical physics.

[44]  Mihaly Mezei,et al.  A novel method reveals that solvent water favors polyproline II over β‐strand conformation in peptides and unfolded proteins: conditional hydrophobic accessible surface area (CHASA) , 2005, Protein science : a publication of the Protein Society.

[45]  Neysa Nevins,et al.  Molecular mechanics (MM4) calculations on conjugated hydrocarbons , 1996, J. Comput. Chem..

[46]  Peter A. Kollman,et al.  AMBER: Assisted model building with energy refinement. A general program for modeling molecules and their interactions , 1981 .

[47]  Alexander D. MacKerell,et al.  An all-atom empirical energy function for the simulation of nucleic acids , 1995 .

[48]  Steven J. Stuart,et al.  Effects of Polarizability on the Hydration of the Chloride Ion , 1996 .

[49]  Alexander D. MacKerell,et al.  CHARMM fluctuating charge force field for proteins: II Protein/solvent properties from molecular dynamics simulations using a nonadditive electrostatic model , 2004, J. Comput. Chem..

[50]  Steven J. Stuart,et al.  Dynamical fluctuating charge force fields: Application to liquid water , 1994 .

[51]  Richard A. Friesner,et al.  Accurate ab Initio Quantum Chemical Determination of the Relative Energetics of Peptide Conformations and Assessment of Empirical Force Fields , 1997 .

[52]  Wei Zhang,et al.  A point‐charge force field for molecular mechanics simulations of proteins based on condensed‐phase quantum mechanical calculations , 2003, J. Comput. Chem..

[53]  George D Rose,et al.  Polyproline II structure in a sequence of seven alanine residues , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[54]  H. Berendsen,et al.  A consistent empirical potential for water–protein interactions , 1984 .

[55]  Gerhard Stock,et al.  Peptide conformational heterogeneity revealed from nonlinear vibrational spectroscopy and molecular dynamics simulations , 2002 .

[56]  Gerhard Stock,et al.  Conformational dynamics of trialanine in water. 2. Comparison of AMBER, CHARMM, GROMOS, and OPLS force fields to NMR and infrared experiments , 2003 .

[57]  J. Hermans,et al.  Comparison of a QM/MM force field and molecular mechanics force fields in simulations of alanine and glycine “dipeptides” (Ace‐Ala‐Nme and Ace‐Gly‐Nme) in water in relation to the problem of modeling the unfolded peptide backbone in solution , 2003, Proteins.

[58]  P. Kollman,et al.  An all atom force field for simulations of proteins and nucleic acids , 1986, Journal of computational chemistry.

[59]  Pengyu Y. Ren,et al.  Polarizable Atomic Multipole Water Model for Molecular Mechanics Simulation , 2003 .

[60]  Martin Karplus,et al.  Molecular dynamics simulations of biomolecules. , 2002, Nature structural biology.

[61]  Charles L. Brooks,et al.  Modern protein force fields behave comparably in molecular dynamics simulations , 2002, J. Comput. Chem..

[62]  Robert W Woody,et al.  Single peptide bonds exhibit poly(pro)II ("random coil") circular dichroism spectra. , 2005, Journal of the American Chemical Society.

[63]  S. Gnanakaran,et al.  Validation of an all-atom protein force field: From dipeptides to larger peptides , 2003 .