Further along the Road Less Traveled: AMBER ff15ipq, an Original Protein Force Field Built on a Self-Consistent Physical Model

We present the AMBER ff15ipq force field for proteins, the second-generation force field developed using the Implicitly Polarized Q (IPolQ) scheme for deriving implicitly polarized atomic charges in the presence of explicit solvent. The ff15ipq force field is a complete rederivation including more than 300 unique atomic charges, 900 unique torsion terms, 60 new angle parameters, and new atomic radii for polar hydrogens. The atomic charges were derived in the context of the SPC/Eb water model, which yields more-accurate rotational diffusion of proteins and enables direct calculation of nuclear magnetic resonance (NMR) relaxation parameters from molecular dynamics simulations. The atomic radii improve the accuracy of modeling salt bridge interactions relative to contemporary fixed-charge force fields, rectifying a limitation of ff14ipq that resulted from its use of pair-specific Lennard-Jones radii. In addition, ff15ipq reproduces penta-alanine J-coupling constants exceptionally well, gives reasonable agreement with NMR relaxation rates, and maintains the expected conformational propensities of structured proteins/peptides, as well as disordered peptides—all on the microsecond (μs) time scale, which is a critical regime for drug design applications. These encouraging results demonstrate the power and robustness of our automated methods for deriving new force fields. All parameters described here and the mdgx program used to fit them are included in the AmberTools16 distribution.

[1]  C. Tanford The Association of Acetate with Ammonium and Guanidinium Ions , 1954 .

[2]  F. Richards,et al.  The preparation of subtilisn-modified ribonuclease and the separation of the peptide and protein components. , 1959, The Journal of biological chemistry.

[3]  Paul Haake,et al.  Equilibrium constants for association of guanidinium and ammonium ions with oxyanions: The effect of changing basicity of the oxyanion , 1977 .

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

[5]  G J Williams,et al.  The Protein Data Bank: a computer-based archival file for macromolecular structures. , 1978, Archives of biochemistry and biophysics.

[6]  W. Kabsch,et al.  Dictionary of protein secondary structure: Pattern recognition of hydrogen‐bonded and geometrical features , 1983, Biopolymers.

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

[8]  R. Huber,et al.  Structure of bovine pancreatic trypsin inhibitor. Results of joint neutron and X-ray refinement of crystal form II. , 1984, Journal of molecular biology.

[9]  C. Bugg,et al.  Structure of ubiquitin refined at 1.8 A resolution. , 1987, Journal of molecular biology.

[10]  T. H. Dunning Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen , 1989 .

[11]  H. G. Petersen,et al.  Error estimates on averages of correlated data , 1989 .

[12]  P. Kollman,et al.  Settle: An analytical version of the SHAKE and RATTLE algorithm for rigid water models , 1992 .

[13]  F M Richards,et al.  Refinement of the crystal structure of ribonuclease S. Comparison with and between the various ribonuclease A structures. , 1994, Biochemistry.

[14]  Peter A. Kollman,et al.  How transferable are hydrogen parameters in molecular mechanics calculations? , 1992 .

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

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

[17]  E. Stellwagen,et al.  Distribution of Helicity within the Model Peptide Acetyl(AAQAA)3amide , 1994 .

[18]  L. Serrano,et al.  A short linear peptide that folds into a native stable β-hairpin in aqueous solution , 1994, Nature Structural Biology.

[19]  A. Bax,et al.  Rotational diffusion anisotropy of human ubiquitin from 15N NMR relaxation , 1995 .

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

[21]  Shankar Kumar,et al.  Multidimensional free‐energy calculations using the weighted histogram analysis method , 1995, J. Comput. Chem..

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

[23]  A. Levine,et al.  Structure of the MDM2 Oncoprotein Bound to the p53 Tumor Suppressor Transactivation Domain , 1996, Science.

[24]  G. Huber,et al.  Weighted-ensemble Brownian dynamics simulations for protein association reactions. , 1996, Biophysical journal.

[25]  U. Hansmann Parallel tempering algorithm for conformational studies of biological molecules , 1997, physics/9710041.

[26]  V. Muñoz,et al.  Folding dynamics and mechanism of β-hairpin formation , 1997, Nature.

[27]  W. Jahnke,et al.  15N backbone dynamics of the S‐peptide from ribonuclease A in its free and S‐protein bound forms: Toward a site‐specific analysis of entropy changes upon folding , 1998, Protein science : a publication of the Protein Society.

[28]  V S Lamzin,et al.  Refinement of triclinic hen egg-white lysozyme at atomic resolution. , 1998, Acta crystallographica. Section D, Biological crystallography.

[29]  D. Nolde,et al.  Three‐dimensional structure of binase in solution , 1998, FEBS letters.

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

[31]  Laxmikant V. Kale,et al.  NAMD2: Greater Scalability for Parallel Molecular Dynamics , 1999 .

[32]  David,et al.  Gaussian basis sets for use in correlated molecular calculations . Ill . The atoms aluminum through argon , 1999 .

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

[34]  D. Case,et al.  Static and Dynamic Effects on Vicinal Scalar J Couplings in Proteins and Peptides: A MD/DFT Analysis , 2000 .

[35]  H. Schwalbe,et al.  Determination of ψ Torsion Angle Restraints from 3J(Cα,Cα) and 3J(HN,Cα) Coupling Constants in Proteins. , 2000 .

[36]  A Mitsutake,et al.  Generalized-ensemble algorithms for molecular simulations of biopolymers. , 2000, Biopolymers.

[37]  Angela K. Wilson,et al.  Gaussian basis sets for use in correlated molecular calculations. X. The atoms aluminum through argon revisited , 2001 .

[38]  S. Grzesiek,et al.  Observation of the closing of individual hydrogen bonds during TFE–induced helix formation in a peptide , 2001, Protein science : a publication of the Protein Society.

[39]  A. Lebedev,et al.  The structure of substrate-free microbial ribonuclease binase and of its complexes with 3'GMP and sulfate ions. , 2002, Acta crystallographica. Section D, Biological crystallography.

[40]  Rafael Brüschweiler,et al.  General framework for studying the dynamics of folded and nonfolded proteins by NMR relaxation spectroscopy and MD simulation. , 2002, Journal of the American Chemical Society.

[41]  J. W. Neidigh,et al.  Designing a 20-residue protein , 2002, Nature Structural Biology.

[42]  H. Schwalbe,et al.  Angular dependence of 1J(Ni,Cα i) and 2J(Ni,Cα (i − 1)) coupling constants measured in J-modulated HSQCs , 2002, Journal of biomolecular NMR.

[43]  E. Zuiderweg,et al.  Temperature dependence of anisotropic protein backbone dynamics. , 2003, Journal of the American Chemical Society.

[44]  Ad Bax,et al.  Evaluation of backbone proton positions and dynamics in a small protein by liquid crystal NMR spectroscopy. , 2003, Journal of the American Chemical Society.

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

[46]  D. Case,et al.  Exploring protein native states and large‐scale conformational changes with a modified generalized born model , 2004, Proteins.

[47]  A. Gronenborn,et al.  Protein Backbone 1H(N)-13Calpha and 15N-13Calpha residual dipolar and J couplings: new constraints for NMR structure determination. , 2004, Journal of the American Chemical Society.

[48]  Greg L. Hura,et al.  Development of an improved four-site water model for biomolecular simulations: TIP4P-Ew. , 2004, The Journal of chemical physics.

[49]  C. Vega,et al.  Phase diagram of water from computer simulation. , 2004, Physical review letters.

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

[51]  Shinya Honda,et al.  10 residue folded peptide designed by segment statistics. , 2004, Structure.

[52]  C. Vega,et al.  The melting temperature of the most common models of water. , 2005, The Journal of chemical physics.

[53]  K. Lindorff-Larsen,et al.  Interpreting Dynamically-Averaged Scalar Couplings in Proteins , 2005, Journal of biomolecular NMR.

[54]  Bosco K. Ho,et al.  The Ramachandran plots of glycine and pre-proline , 2005, BMC Structural Biology.

[55]  J. Hofrichter,et al.  Sub-microsecond protein folding. , 2006, Journal of molecular biology.

[56]  Neal J. Zondlo,et al.  Determinants of specificity of MDM2 for the activation domains of p53 and p65: proline27 disrupts the MDM2-binding motif of p53. , 2006, Biochemistry.

[57]  Virginia A. Jarymowycz,et al.  Fast time scale dynamics of protein backbones: NMR relaxation methods, applications, and functional consequences. , 2006, Chemical reviews.

[58]  V. Hornak,et al.  Comparison of multiple Amber force fields and development of improved protein backbone parameters , 2006, Proteins.

[59]  Hao Hu,et al.  Fitting Molecular Electrostatic Potentials from Quantum Mechanical Calculations. , 2007, Journal of chemical theory and computation.

[60]  H. Schwalbe,et al.  Structure and dynamics of the homologous series of alanine peptides: a joint molecular dynamics/NMR study. , 2007, Journal of the American Chemical Society.

[61]  Shinya Honda,et al.  Crystal structure of a ten-amino acid protein. , 2008, Journal of the American Chemical Society.

[62]  C. Simmerling,et al.  Structural insights for designed alanine‐rich helices: Comparing NMR helicity measures and conformational ensembles from molecular dynamics simulation , 2008, Biopolymers.

[63]  D. Case,et al.  Evaluating rotational diffusion from protein MD simulations. , 2008, The journal of physical chemistry. B.

[64]  G. Hummer,et al.  Are current molecular dynamics force fields too helical? , 2008, Biophysical journal.

[65]  R. Dror,et al.  Long-timescale molecular dynamics simulations of protein structure and function. , 2009, Current opinion in structural biology.

[66]  G. Hummer,et al.  Optimized molecular dynamics force fields applied to the helix-coil transition of polypeptides. , 2009, The journal of physical chemistry. B.

[67]  Michael I. Jordan,et al.  Neighbor-Dependent Ramachandran Probability Distributions of Amino Acids Developed from a Hierarchical Dirichlet Process Model , 2010, PLoS Comput. Biol..

[68]  R. Best,et al.  Protein simulations with an optimized water model: cooperative helix formation and temperature-induced unfolded state collapse. , 2010, The journal of physical chemistry. B.

[69]  Amparo Galindo,et al.  The Use of Anisotropic Potentials in Modeling Water and Free Energies of Hydration. , 2010, Journal of chemical theory and computation.

[70]  Klaus Schulten,et al.  GPU-accelerated molecular modeling coming of age. , 2010, Journal of molecular graphics & modelling.

[71]  R. Brüschweiler,et al.  NMR-based protein potentials. , 2010, Angewandte Chemie.

[72]  R. Dror,et al.  Improved side-chain torsion potentials for the Amber ff99SB protein force field , 2010, Proteins.

[73]  Klaus Schulten,et al.  High-performance scalable molecular dynamics simulations of a polarizable force field based on classical Drude oscillators in NAMD. , 2011, The journal of physical chemistry letters.

[74]  A. Garcia,et al.  Influence of water-protein hydrogen bonding on the stability of Trp-cage miniprotein. A comparison between the TIP3P and TIP4P-Ew water models. , 2011, Physical chemistry chemical physics : PCCP.

[75]  K. Lindorff-Larsen,et al.  How robust are protein folding simulations with respect to force field parameterization? , 2011, Biophysical journal.

[76]  R. Brüschweiler,et al.  Iterative Optimization of Molecular Mechanics Force Fields from NMR Data of Full-Length Proteins. , 2011, Journal of chemical theory and computation.

[77]  Pengyu Ren,et al.  Automation of AMOEBA polarizable force field parameterization for small molecules , 2012, Theoretical Chemistry Accounts.

[78]  A. Kitao,et al.  Water model tuning for improved reproduction of rotational diffusion and NMR spectral density. , 2012, The journal of physical chemistry. B.

[79]  Patric Schyman,et al.  Treatment of Halogen Bonding in the OPLS-AA Force Field; Application to Potent Anti-HIV Agents. , 2012, Journal of chemical theory and computation.

[80]  Alexander D. MacKerell,et al.  Optimization of the additive CHARMM all-atom protein force field targeting improved sampling of the backbone φ, ψ and side-chain χ(1) and χ(2) dihedral angles. , 2012, Journal of chemical theory and computation.

[81]  Marcus D. Hanwell,et al.  Avogadro: an advanced semantic chemical editor, visualization, and analysis platform , 2012, Journal of Cheminformatics.

[82]  Duncan Poole,et al.  Routine Microsecond Molecular Dynamics Simulations with AMBER on GPUs. 1. Generalized Born , 2012, Journal of chemical theory and computation.

[83]  D. Case,et al.  Derivation of fixed partial charges for amino acids accommodating a specific water model and implicit polarization. , 2013, The journal of physical chemistry. B.

[84]  Pengyu Y. Ren,et al.  The Polarizable Atomic Multipole-based AMOEBA Force Field for Proteins. , 2013, Journal of chemical theory and computation.

[85]  Andreas W. Götz,et al.  SPFP: Speed without compromise - A mixed precision model for GPU accelerated molecular dynamics simulations , 2013, Comput. Phys. Commun..

[86]  M. Shirts,et al.  Effects of Temperature Control Algorithms on Transport Properties and Kinetics in Molecular Dynamics Simulations. , 2013, Journal of chemical theory and computation.

[87]  Daniel R Roe,et al.  PTRAJ and CPPTRAJ: Software for Processing and Analysis of Molecular Dynamics Trajectory Data. , 2013, Journal of chemical theory and computation.

[88]  Alexander D. MacKerell,et al.  Force Field for Peptides and Proteins based on the Classical Drude Oscillator. , 2013, Journal of chemical theory and computation.

[89]  Duncan Poole,et al.  Routine Microsecond Molecular Dynamics Simulations with AMBER on GPUs. 2. Explicit Solvent Particle Mesh Ewald. , 2013, Journal of chemical theory and computation.

[90]  D. Case,et al.  ff14ipq: A Self-Consistent Force Field for Condensed-Phase Simulations of Proteins , 2014, Journal of chemical theory and computation.

[91]  R. Best,et al.  Balanced Protein–Water Interactions Improve Properties of Disordered Proteins and Non-Specific Protein Association , 2014, Journal of chemical theory and computation.

[92]  Vijay S Pande,et al.  Building Force Fields: An Automatic, Systematic, and Reproducible Approach. , 2014, The journal of physical chemistry letters.

[93]  Adrian E. Roitberg,et al.  Fitting of Dihedral Terms in Classical Force Fields as an Analytic Linear Least-Squares Problem , 2014, J. Chem. Inf. Model..

[94]  L. Chong,et al.  Simultaneous Computation of Dynamical and Equilibrium Information Using a Weighted Ensemble of Trajectories , 2012, Journal of chemical theory and computation.

[95]  Karl T. Debiec,et al.  Evaluating the Strength of Salt Bridges: A Comparison of Current Biomolecular Force Fields , 2014, The journal of physical chemistry. B.

[96]  Fangqiang Zhu,et al.  Role of electrostatic interactions in binding of peptides and intrinsically disordered proteins to their folded targets. 1. NMR and MD characterization of the complex between the c-Crk N-SH3 domain and the peptide Sos. , 2014, Biochemistry.

[97]  F. Jiang,et al.  Residue-specific force field based on the protein coil library. RSFF1: modification of OPLS-AA/L. , 2014, The journal of physical chemistry. B.

[98]  Alex Dickson,et al.  WExplore: hierarchical exploration of high-dimensional spaces using the weighted ensemble algorithm. , 2014, The journal of physical chemistry. B.

[99]  W. L. Jorgensen,et al.  Improved Peptide and Protein Torsional Energetics with the OPLS-AA Force Field , 2015, Journal of chemical theory and computation.

[100]  F. Jiang,et al.  Residue-specific force field based on protein coil library. RSFF2: modification of AMBER ff99SB. , 2015, The journal of physical chemistry. B.

[101]  A. Roitberg,et al.  Long-Time-Step Molecular Dynamics through Hydrogen Mass Repartitioning. , 2015, Journal of chemical theory and computation.

[102]  C. Simmerling,et al.  ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from ff99SB. , 2015, Journal of chemical theory and computation.

[103]  Paul Robustelli,et al.  Water dispersion interactions strongly influence simulated structural properties of disordered protein states. , 2015, The journal of physical chemistry. B.

[104]  Alexander D. MacKerell,et al.  Robustness in the fitting of molecular mechanics parameters , 2015, J. Comput. Chem..

[105]  N. Skrynnikov,et al.  Role of Electrostatic Interactions in Binding of Peptides and Intrinsically Disordered Proteins to Their Folded Targets: 2. The Model of Encounter Complex Involving the Double Mutant of the c-Crk N-SH3 Domain and Peptide Sos. , 2016, Biochemistry.

[106]  Jennifer L. Knight,et al.  OPLS3: A Force Field Providing Broad Coverage of Drug-like Small Molecules and Proteins. , 2016, Journal of chemical theory and computation.