OPLS3: A Force Field Providing Broad Coverage of Drug-like Small Molecules and Proteins.

The parametrization and validation of the OPLS3 force field for small molecules and proteins are reported. Enhancements with respect to the previous version (OPLS2.1) include the addition of off-atom charge sites to represent halogen bonding and aryl nitrogen lone pairs as well as a complete refit of peptide dihedral parameters to better model the native structure of proteins. To adequately cover medicinal chemical space, OPLS3 employs over an order of magnitude more reference data and associated parameter types relative to other commonly used small molecule force fields (e.g., MMFF and OPLS_2005). As a consequence, OPLS3 achieves a high level of accuracy across performance benchmarks that assess small molecule conformational propensities and solvation. The newly fitted peptide dihedrals lead to significant improvements in the representation of secondary structure elements in simulated peptides and native structure stability over a number of proteins. Together, the improvements made to both the small molecule and protein force field lead to a high level of accuracy in predicting protein-ligand binding measured over a wide range of targets and ligands (less than 1 kcal/mol RMS error) representing a 30% improvement over earlier variants of the OPLS force field.

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

[2]  K. Dill,et al.  Predicting absolute ligand binding free energies to a simple model site. , 2007, Journal of molecular biology.

[3]  Eric T. Kim,et al.  How does a drug molecule find its target binding site? , 2011, Journal of the American Chemical Society.

[4]  A. Barducci,et al.  Accuracy of current all-atom force-fields in modeling protein disordered states. , 2015, Journal of chemical theory and computation.

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

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

[7]  A. Moretto,et al.  Structure-based optimization of protein tyrosine phosphatase 1B inhibitors: from the active site to the second phosphotyrosine binding site. , 2007, Journal of medicinal chemistry.

[8]  William L. Jorgensen,et al.  Characterization of Biaryl Torsional Energetics and its Treatment in OPLS All-Atom Force Fields , 2013, J. Chem. Inf. Model..

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

[10]  V S Pande,et al.  Molecular dynamics simulations of unfolding and refolding of a beta-hairpin fragment of protein G. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[11]  W. L. Jorgensen,et al.  The OPLS [optimized potentials for liquid simulations] potential functions for proteins, energy minimizations for crystals of cyclic peptides and crambin. , 1988, Journal of the American Chemical Society.

[12]  William L Jorgensen,et al.  Halide, Ammonium, and Alkali Metal Ion Parameters for Modeling Aqueous Solutions. , 2006, Journal of chemical theory and computation.

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

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

[15]  P. Kovar,et al.  Synthesis and biological evaluation of 3-ethylidene-1,3-dihydro-indol-2-ones as novel checkpoint 1 inhibitors. , 2006, Bioorganic & medicinal chemistry letters.

[16]  Yuzhu Chen,et al.  N2-substituted O6-cyclohexylmethylguanine derivatives: potent inhibitors of cyclin-dependent kinases 1 and 2. , 2004, Journal of medicinal chemistry.

[17]  Jennifer L. Knight,et al.  Modeling Local Structural Rearrangements Using FEP/REST: Application to Relative Binding Affinity Predictions of CDK2 Inhibitors. , 2013, Journal of chemical theory and computation.

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

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

[20]  Wei Yang,et al.  The Effects of Computational Modeling Errors on the Estimation of Statistical Mechanical Variables. , 2012, Journal of chemical theory and computation.

[21]  David L. Mobley,et al.  Alchemical prediction of hydration free energies for SAMPL , 2012, Journal of Computer-Aided Molecular Design.

[22]  W. Sherman,et al.  Prediction of Absolute Solvation Free Energies using Molecular Dynamics Free Energy Perturbation and the OPLS Force Field. , 2010, Journal of chemical theory and computation.

[23]  Brian A. Chauder,et al.  Discovery of potent myeloid cell leukemia 1 (Mcl-1) inhibitors using fragment-based methods and structure-based design. , 2013, Journal of medicinal chemistry.

[24]  Chang Park,et al.  Structure-based design, synthesis, and biological evaluation of potent and selective macrocyclic checkpoint kinase 1 inhibitors. , 2007, Journal of medicinal chemistry.

[25]  Michael Reilly,et al.  Development of 6-substituted indolylquinolinones as potent Chek1 kinase inhibitors. , 2006, Bioorganic & medicinal chemistry letters.

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

[27]  R. L. Baldwin,et al.  N‐ and C‐capping preferences for all 20 amino acids in α‐helical peptides , 1995, Protein science : a publication of the Protein Society.

[28]  Pavel Hobza,et al.  On Extension of the Current Biomolecular Empirical Force Field for the Description of Halogen Bonds. , 2012, Journal of chemical theory and computation.

[29]  Richard A Friesner,et al.  Accurate Force Field Development for Modeling Conjugated Polymers. , 2012, Journal of chemical theory and computation.

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

[31]  Michael R. Shirts,et al.  Direct calculation of the binding free energies of FKBP ligands. , 2005, The Journal of chemical physics.

[32]  Jennifer L. Knight,et al.  Accurate and reliable prediction of relative ligand binding potency in prospective drug discovery by way of a modern free-energy calculation protocol and force field. , 2015, Journal of the American Chemical Society.

[33]  Scott P. Brown,et al.  Healthy skepticism: assessing realistic model performance. , 2009, Drug discovery today.

[34]  Timothy Clark,et al.  Halogen bonding: an electrostatically-driven highly directional noncovalent interaction. , 2010, Physical chemistry chemical physics : PCCP.

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

[36]  Kyle A Beauchamp,et al.  Toward Automated Benchmarking of Atomistic Force Fields: Neat Liquid Densities and Static Dielectric Constants from the ThermoML Data Archive. , 2015, The journal of physical chemistry. B.

[37]  Michael R. Shirts,et al.  Simulation of folding of a small alpha-helical protein in atomistic detail using worldwide-distributed computing. , 2002, Journal of molecular biology.

[38]  Benoît Roux,et al.  Atomic Level Anisotropy in the Electrostatic Modeling of Lone Pairs for a Polarizable Force Field Based on the Classical Drude Oscillator. , 2006, Journal of chemical theory and computation.

[39]  William L. Jorgensen,et al.  OPLS ALL-ATOM MODEL FOR AMINES : RESOLUTION OF THE AMINE HYDRATION PROBLEM , 1999 .

[40]  Woody Sherman,et al.  Improving the Prediction of Absolute Solvation Free Energies Using the Next Generation OPLS Force Field. , 2012, Journal of chemical theory and computation.

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

[42]  B. Graves,et al.  Discovery of 6-(2,4-difluorophenoxy)-2-[3-hydroxy-1-(2-hydroxyethyl)propylamino]-8-methyl-8H-pyrido[2,3-d]pyrimidin-7-one (pamapimod) and 6-(2,4-difluorophenoxy)-8-methyl-2-(tetrahydro-2H-pyran-4-ylamino)pyrido[2,3-d]pyrimidin-7(8H)-one (R1487) as orally bioavailable and highly selective inhibitors , 2011, Journal of medicinal chemistry.

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

[44]  David L Mobley,et al.  Calculating the sensitivity and robustness of binding free energy calculations to force field parameters. , 2013, Journal of chemical theory and computation.

[45]  R. Dror,et al.  How Fast-Folding Proteins Fold , 2011, Science.

[46]  J. Breed,et al.  Discovery of checkpoint kinase inhibitor (S)-5-(3-fluorophenyl)-N-(piperidin-3-yl)-3-ureidothiophene-2-carboxamide (AZD7762) by structure-based design and optimization of thiophenecarboxamide ureas. , 2012, Journal of medicinal chemistry.

[47]  Adam R. Johnson,et al.  Lead optimization of a 4-aminopyridine benzamide scaffold to identify potent, selective, and orally bioavailable TYK2 inhibitors. , 2013, Journal of medicinal chemistry.

[48]  Lynn A. Hyde,et al.  Structure based design of iminohydantoin BACE1 inhibitors: identification of an orally available, centrally active BACE1 inhibitor. , 2012, Bioorganic & medicinal chemistry letters.

[49]  Chaohong Sun,et al.  Aminopyridine-based c-Jun N-terminal kinase inhibitors with cellular activity and minimal cross-kinase activity. , 2006, Journal of medicinal chemistry.

[50]  Mahmoud A. A. Ibrahim,et al.  Molecular mechanical study of halogen bonding in drug discovery , 2011, J. Comput. Chem..

[51]  Adam R. Johnson,et al.  Lead identification of novel and selective TYK2 inhibitors. , 2013, European journal of medicinal chemistry.

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

[53]  R. Dror,et al.  Systematic Validation of Protein Force Fields against Experimental Data , 2012, PloS one.

[54]  Richard A. Friesner,et al.  Pseudospectral localized Mo/ller–Plesset methods: Theory and calculation of conformational energies , 1995 .

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

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

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

[58]  William L. Jorgensen,et al.  Development of an all-atom force field for heterocycles. Properties of liquid pyridine and diazenes , 1998 .

[59]  G. Klebe,et al.  More than a simple lipophilic contact: a detailed thermodynamic analysis of nonbasic residues in the s1 pocket of thrombin. , 2009, Journal of molecular biology.

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

[61]  Wilfred F. van Gunsteren,et al.  An improved GROMOS96 force field for aliphatic hydrocarbons in the condensed phase , 2001, J. Comput. Chem..

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

[63]  Jing Zhang,et al.  Jaguar: A high-performance quantum chemistry software program with strengths in life and materials sciences , 2013 .

[64]  David Parry,et al.  Design, synthesis and SAR of thienopyridines as potent CHK1 inhibitors. , 2010, Bioorganic & medicinal chemistry letters.

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

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

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

[68]  David J. Giesen,et al.  Class IV charge models: A new semiempirical approach in quantum chemistry , 1995, J. Comput. Aided Mol. Des..

[69]  Thierry O Fischmann,et al.  Discovery of a Novel Series of CHK1 Kinase Inhibitors with a Distinctive Hinge Binding Mode. , 2012, ACS medicinal chemistry letters.

[70]  Stefan Senger,et al.  Factor Xa inhibitors: S1 binding interactions of a series of N-{(3S)-1-[(1S)-1-methyl-2-morpholin-4-yl-2-oxoethyl]-2-oxopyrrolidin-3-yl}sulfonamides. , 2007, Journal of medicinal chemistry.

[71]  David L. Mobley,et al.  FreeSolv: a database of experimental and calculated hydration free energies, with input files , 2014, Journal of Computer-Aided Molecular Design.

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

[73]  Rhiju Das,et al.  Are Protein Force Fields Getting Better? A Systematic Benchmark on 524 Diverse NMR Measurements. , 2012, Journal of chemical theory and computation.