Development of polarizable models for molecular mechanical calculations. 4. van der Waals parametrization.

In the previous publications of this series, we presented a set of Thole induced dipole interaction models using four types of screening functions. In this work, we document our effort to refine the van der Waals parameters for the Thole polarizable models. Following the philosophy of AMBER force field development, the van der Waals (vdW) parameters were tuned for the Thole model with linear screening function to reproduce both the ab initio interaction energies and the experimental densities of pure liquids. An in-house genetic algorithm was applied to maximize the fitness of "chromosomes" which is a function of the root-mean-square errors (RMSE) of interaction energy and liquid density. To efficiently explore the vdW parameter space, a novel approach was developed to estimate the liquid densities for a given vdW parameter set using the mean residue-residue interaction energies through interpolation/extrapolation. This approach allowed the costly molecular dynamics simulations be performed at the end of each optimization cycle only and eliminated the simulations during the cycle. Test results show notable improvements over the original AMBER FF99 vdW parameter set, as indicated by the reduction in errors of the calculated pure liquid densities (d), heats of vaporization (H(vap)), and hydration energies. The average percent error (APE) of the densities of 59 pure liquids was reduced from 5.33 to 2.97%; the RMSE of H(vap) was reduced from 1.98 to 1.38 kcal/mol; the RMSE of solvation free energies of 15 compounds was reduced from 1.56 to 1.38 kcal/mol. For the interaction energies of 1639 dimers, the overall performance of the optimized vdW set is slightly better than the original FF99 vdW set (RMSE of 1.56 versus 1.63 kcal/mol). The optimized vdW parameter set was also evaluated for the exponential screening function used in the Amoeba force field to assess its applicability for different types of screening functions. Encouragingly, comparable performance was observed when the optimized vdW set was combined with the Thole Amoeba-like polarizable model, particularly for the interaction energy and liquid density calculations. Thus, the optimized vdW set is applicable to both types of Thole models with either linear or Amoeba-like screening functions.

[1]  P. Kollman,et al.  How well does a restrained electrostatic potential (RESP) model perform in calculating conformational energies of organic and biological molecules? , 2000 .

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

[3]  Alexander D. MacKerell,et al.  Combined ab initio/empirical approach for optimization of Lennard‐Jones parameters for polar‐neutral compounds , 2002, J. Comput. Chem..

[4]  Jie Li,et al.  Development of polarizable models for molecular mechanical calculations II: induced dipole models significantly improve accuracy of intermolecular interaction energies. , 2011, The journal of physical chemistry. B.

[5]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

[6]  S. F. Boys,et al.  The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors , 1970 .

[7]  H. Sun,et al.  COMPASS: An ab Initio Force-Field Optimized for Condensed-Phase ApplicationsOverview with Details on Alkane and Benzene Compounds , 1998 .

[8]  George A. Kaminski,et al.  Development of an Accurate and Robust Polarizable Molecular Mechanics Force Field from ab Initio Quantum Chemistry , 2004 .

[9]  Junmei Wang,et al.  Genetic Algorithm-Optimized QSPR Models for Bioavailability, Protein Binding, and Urinary Excretion , 2006, J. Chem. Inf. Model..

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

[11]  J. Šponer,et al.  Comparison of intrinsic stacking energies of ten unique dinucleotide steps in A-RNA and B-DNA duplexes. Can we determine correct order of stability by quantum-chemical calculations? , 2010, The journal of physical chemistry. B.

[12]  Junmei Wang,et al.  Application of Molecular Dynamics Simulations in Molecular Property Prediction I: Density and Heat of Vaporization. , 2011, Journal of chemical theory and computation.

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

[14]  Marian Anghel,et al.  Synchronization of trajectories in canonical molecular-dynamics simulations: observation, explanation, and exploitation. , 2004, The Journal of chemical physics.

[15]  Daniel Svozil,et al.  Reference MP2/CBS and CCSD(T) quantum-chemical calculations on stacked adenine dimers. Comparison with DFT-D, MP2.5, SCS(MI)-MP2, M06-2X, CBS(SCS-D) and force field descriptions. , 2010, Physical chemistry chemical physics : PCCP.

[16]  B. Brooks,et al.  Langevin dynamics of peptides: The frictional dependence of isomerization rates of N‐acetylalanyl‐N′‐methylamide , 1992, Biopolymers.

[17]  Kazutoshi Tanabe,et al.  Refinement of Nonbonding Interaction Potential Parameters for Methane on the Basis of the Pair Potential Obtained by MP3/6-311G(3d,3p)-Level ab Initio Molecular Orbital Calculations: The Anisotropy of H/H Interaction , 1994 .

[18]  P. Kollman,et al.  Solvation Model Based on Weighted Solvent Accessible Surface Area , 2001 .

[19]  David L Mobley,et al.  Small molecule hydration free energies in explicit solvent: An extensive test of fixed-charge atomistic simulations. , 2009, Journal of chemical theory and computation.

[20]  William L. Jorgensen,et al.  Gas‐phase and liquid‐state properties of esters, nitriles, and nitro compounds with the OPLS‐AA force field , 2001, J. Comput. Chem..

[21]  Alexander D. MacKerell,et al.  An Improved Empirical Potential Energy Function for Molecular Simulations of Phospholipids , 2000 .

[22]  H. V. Kehiaian,et al.  Enthalpies of vaporization of organic compounds : a critical review and data compilation , 1985 .

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

[24]  Yang Yang,et al.  An extensible and systematic force field, ESFF, for molecular modeling of organic, inorganic, and organometallic systems , 2003, J. Comput. Chem..

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

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

[27]  Donald G. Truhlar,et al.  Basis-set extrapolation , 1998 .

[28]  W. M. Haynes.  CRC Handbook of Chemistry and Physics , 1990 .

[29]  W. Goddard,et al.  UFF, a full periodic table force field for molecular mechanics and molecular dynamics simulations , 1992 .

[30]  Wei Zhang,et al.  Strike a balance: Optimization of backbone torsion parameters of AMBER polarizable force field for simulations of proteins and peptides , 2006, J. Comput. Chem..

[31]  David L Mobley,et al.  Nonlinear scaling schemes for Lennard-Jones interactions in free energy calculations. , 2007, The Journal of chemical physics.

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

[33]  Y. Duan,et al.  Accurate ab Initio Study on the Hydrogen-Bond Pairs in Protein Secondary Structures. , 2007, Journal of chemical theory and computation.

[34]  W. V. Gunsteren,et al.  Validation of the 53A6 GROMOS force field , 2005, European Biophysics Journal.

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

[36]  R. Mannella,et al.  Langevin stabilization of molecular-dynamics simulations of polymers by means of quasisymplectic algorithms. , 2007, The Journal of chemical physics.

[37]  William L. Jorgensen,et al.  Development of an All-Atom Force Field for Heterocycles. Properties of Liquid Pyrrole, Furan, Diazoles, and Oxazoles , 1998 .

[38]  D. Truhlar,et al.  Infinite basis limits in electronic structure theory , 1999 .

[39]  Pavel Hobza,et al.  Accurate interaction energies of hydrogen-bonded nucleic acid base pairs. , 2004, Journal of the American Chemical Society.

[40]  T Darden,et al.  New tricks for modelers from the crystallography toolkit: the particle mesh Ewald algorithm and its use in nucleic acid simulations. , 1999, Structure.

[41]  Piotr Cieplak,et al.  Polarization effects in molecular mechanical force fields , 2009, Journal of physics. Condensed matter : an Institute of Physics journal.

[42]  T. Halgren.  Merck molecular force field. II. MMFF94 van der Waals and electrostatic parameters for intermolecular interactions , 1996 .

[43]  David L Mobley,et al.  Accurate and efficient corrections for missing dispersion interactions in molecular simulations. , 2007, The journal of physical chemistry. B.

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

[45]  Alexander D. MacKerell,et al.  Combined ab initio/empirical approach for optimization of Lennard–Jones parameters , 1998 .

[46]  J. Åqvist,et al.  Ion-water interaction potentials derived from free energy perturbation simulations , 1990 .

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

[48]  J. Šponer,et al.  Balance of Attraction and Repulsion in Nucleic-Acid Base Stacking: CCSD(T)/Complete-Basis-Set-Limit Calculations on Uracil Dimer and a Comparison with the Force-Field Description. , 2009, Journal of chemical theory and computation.

[49]  Tingjun Hou,et al.  Automated docking of peptides and proteins by using a genetic algorithm combined with a tabu search. , 1999, Protein engineering.

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

[51]  P. Kollman,et al.  Automatic atom type and bond type perception in molecular mechanical calculations. , 2006, Journal of molecular graphics & modelling.

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

[53]  Celeste Sagui,et al.  Towards an accurate representation of electrostatics in classical force fields: efficient implementation of multipolar interactions in biomolecular simulations. , 2004, The Journal of chemical physics.

[54]  Jie Li,et al.  Development of polarizable models for molecular mechanical calculations I: parameterization of atomic polarizability. , 2011, The journal of physical chemistry. B.

[55]  R. Skeel,et al.  Langevin stabilization of molecular dynamics , 2001 .

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

[57]  Miquel Duran,et al.  How does basis set superposition error change the potential surfaces for hydrogen-bonded dimers? , 1996 .

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

[59]  Pavel Hobza,et al.  S66: A Well-balanced Database of Benchmark Interaction Energies Relevant to Biomolecular Structures , 2011, Journal of chemical theory and computation.