A point‐charge force field for molecular mechanics simulations of proteins based on condensed‐phase quantum mechanical calculations

Molecular mechanics models have been applied extensively to study the dynamics of proteins and nucleic acids. Here we report the development of a third‐generation point‐charge all‐atom force field for proteins. Following the earlier approach of Cornell et al., the charge set was obtained by fitting to the electrostatic potentials of dipeptides calculated using B3LYP/cc‐pVTZ//HF/6‐31G** quantum mechanical methods. The main‐chain torsion parameters were obtained by fitting to the energy profiles of Ace‐Ala‐Nme and Ace‐Gly‐Nme di‐peptides calculated using MP2/cc‐pVTZ//HF/6‐31G** quantum mechanical methods. All other parameters were taken from the existing AMBER data base. The major departure from previous force fields is that all quantum mechanical calculations were done in the condensed phase with continuum solvent models and an effective dielectric constant of ε = 4. We anticipate that this force field parameter set will address certain critical short comings of previous force fields in condensed‐phase simulations of proteins. Initial tests on peptides demonstrated a high‐degree of similarity between the calculated and the statistically measured Ramanchandran maps for both Ace‐Gly‐Nme and Ace‐Ala‐Nme di‐peptides. Some highlights of our results include (1) well‐preserved balance between the extended and helical region distributions, and (2) favorable type‐II poly‐proline helical region in agreement with recent experiments. Backward compatibility between the new and Cornell et al. charge sets, as judged by overall agreement between dipole moments, allows a smooth transition to the new force field in the area of ligand‐binding calculations. Test simulations on a large set of proteins are also discussed. © 2003 Wiley Periodicals, Inc. J Comput Chem 24: 1999–2012, 2003

[1]  Thomas A. Darden,et al.  Adventures in Improving the Scaling and Accuracy of a Parallel Molecular Dynamics Program , 1997, The Journal of Supercomputing.

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

[3]  Valerie Daggett,et al.  The role of α‐, 310‐, and π‐helix in helix→coil transitions , 2003 .

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

[5]  M. Levitt,et al.  Energy functions that discriminate X-ray and near native folds from well-constructed decoys. , 1996, Journal of molecular biology.

[6]  A. Voter.  Parallel replica method for dynamics of infrequent events , 1998 .

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

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

[9]  P. A. Bash,et al.  Free energy calculations by computer simulation. , 1987, Science.

[10]  Charles L. Brooks,et al.  Molecular picture of folding of a small α/β protein , 1998 .

[11]  Shibasish Chowdhury,et al.  Ab initio folding simulation of the Trp-cage mini-protein approaches NMR resolution. , 2003, Journal of molecular biology.

[12]  J. Skolnick,et al.  A distance‐dependent atomic knowledge‐based potential for improved protein structure selection , 2001, Proteins.

[13]  Parr,et al.  Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. , 1988, Physical review. B, Condensed matter.

[14]  S C Harvey,et al.  Dielectric relaxation spectra of water adsorbed on lysozyme. , 1972, The Journal of physical chemistry.

[15]  E. Lattman,et al.  High apparent dielectric constants in the interior of a protein reflect water penetration. , 2000, Biophysical journal.

[16]  W A Shirley,et al.  Curious structure in “canonical” alanine‐based peptides , 1997, Proteins.

[17]  E. Shakhnovich.  Theoretical studies of protein-folding thermodynamics and kinetics. , 1997, Current opinion in structural biology.

[18]  R. L. Baldwin,et al.  Comparison of NH exchange and circular dichroism as techniques for measuring the parameters of the helix-coil transition in peptides. , 1997, Biochemistry.

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

[20]  Alexander D. MacKerell,et al.  Force field influence on the observation of π-helical protein structures in molecular dynamics simulations , 2003 .

[21]  Chun Wu,et al.  Breaking non-native hydrophobic clusters is the rate-limiting step in the folding of an alanine-based peptide. , 2003, Biopolymers.

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

[23]  J M Rosenberg,et al.  Molecular dynamics simulation study of DNA dodecamer d(CGCGAATTCGCG) in solution: conformation and hydration. , 1997, Journal of molecular biology.

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

[25]  Eric J. Sorin,et al.  β-hairpin folding simulations in atomistic detail using an implicit solvent model1 , 2001 .

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

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

[28]  Charles L. Brooks,et al.  Identifying native‐like protein structures using physics‐based potentials , 2002, J. Comput. Chem..

[29]  Christian Silvio Pomelli,et al.  An improved iterative solution to solve the electrostatic problem in the polarizable continuum model , 2001 .

[30]  E. Lattman,et al.  Experimental measurement of the effective dielectric in the hydrophobic core of a protein. , 1997, Biophysical chemistry.

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

[32]  Junmei Wang,et al.  Automatic parameterization of force field by systematic search and genetic algorithms , 2001, J. Comput. Chem..

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

[34]  T. Dunning,et al.  Electron affinities of the first‐row atoms revisited. Systematic basis sets and wave functions , 1992 .

[35]  김삼묘,et al.  “Bioinformatics” 특집을 내면서 , 2000 .

[36]  Guoli Wang,et al.  PISCES: a protein sequence culling server , 2003, Bioinform..

[37]  L Wang,et al.  The early stage of folding of villin headpiece subdomain observed in a 200-nanosecond fully solvated molecular dynamics simulation. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[38]  F. Young.  Biochemistry , 1955, The Indian Medical Gazette.

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

[40]  D. Case,et al.  Molecular Dynamics Simulations of Nucleic Acids with a Generalized Born Solvation Model , 2000 .

[41]  P A Kollman,et al.  Observation of the A-DNA to B-DNA transition during unrestrained molecular dynamics in aqueous solution. , 1996, Journal of molecular biology.

[42]  D. Case,et al.  Generalized born models of macromolecular solvation effects. , 2000, Annual review of physical chemistry.

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

[44]  J. Skolnick,et al.  TOUCHSTONE: An ab initio protein structure prediction method that uses threading-based tertiary restraints , 2001, Proceedings of the National Academy of Sciences of the United States of America.

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

[46]  A. Roitberg,et al.  All-atom structure prediction and folding simulations of a stable protein. , 2002, Journal of the American Chemical Society.

[47]  P. Kollman,et al.  COMBINED LOCALLY ENHANCED SAMPLING AND PARTICLE MESH EWALD AS A STRATEGY TO LOCATE THE EXPERIMENTAL STRUCTURE OF A NONHELICAL NUCLEIC ACID , 1998 .

[48]  Anthony K. Felts,et al.  Distinguishing native conformations of proteins from decoys with an effective free energy estimator based on the OPLS all‐atom force field and the surface generalized born solvent model , 2002, Proteins.

[49]  Andreas Klamt,et al.  Treatment of the outlying charge in continuum solvation models , 1996 .

[50]  R Elber,et al.  Computer determination of peptide conformations in water: different roads to structure. , 1995, Proceedings of the National Academy of Sciences of the United States of America.

[51]  Xiongwu Wu,et al.  Molecular dynamics simulations of synthetic peptide folding , 1996, Proteins.

[52]  V. Pande,et al.  Absolute comparison of simulated and experimental protein-folding dynamics , 2002, Nature.

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

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

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

[56]  J. Apostolakis,et al.  Thermodynamics and Kinetics of Folding of Two Model Peptides Investigated by Molecular Dynamics Simulations , 2000 .

[57]  Hermann Stoll,et al.  Results obtained with the correlation energy density functionals of becke and Lee, Yang and Parr , 1989 .

[58]  Andreas Klamt,et al.  Incorporation of solvent effects into density functional calculations of molecular energies and geometries , 1995 .

[59]  V. Pande,et al.  The Trp cage: folding kinetics and unfolded state topology via molecular dynamics simulations. , 2002, Journal of the American Chemical Society.

[60]  A. Klamt.  Conductor-like Screening Model for Real Solvents: A New Approach to the Quantitative Calculation of Solvation Phenomena , 1995 .