Application of Molecular Dynamics Simulations in Molecular Property Prediction I: Density and Heat of Vaporization.

Molecular mechanical force field (FF) methods are useful in studying condensed phase properties. They are complementary to experiment and can often go beyond experiment in atomic details. Even a FF is specific for studying structures, dynamics and functions of biomolecules, it is still important for the FF to accurately reproduce the experimental liquid properties of small molecules that represent the chemical moieties of biomolecules. Otherwise, the force field may not describe the structures and energies of macromolecules in aqueous solutions properly. In this work, we have carried out a systematic study to evaluate the General AMBER Force Field (GAFF) in studying densities and heats of vaporization for a large set of organic molecules that covers the most common chemical functional groups. The latest techniques, such as the particle mesh Ewald (PME) for calculating electrostatic energies, and Langevin dynamics for scaling temperatures, have been applied in the molecular dynamics (MD) simulations. For density, the average percent error (APE) of 71 organic compounds is 4.43% when compared to the experimental values. More encouragingly, the APE drops to 3.43% after the exclusion of two outliers and four other compounds for which the experimental densities have been measured with pressures higher than 1.0 atm. For heat of vaporization, several protocols have been investigated and the best one, P4/ntt0, achieves an average unsigned error (AUE) and a root-mean-square error (RMSE) of 0.93 and 1.20 kcal/mol, respectively. How to reduce the prediction errors through proper van der Waals (vdW) parameterization has been discussed. An encouraging finding in vdW parameterization is that both densities and heats of vaporization approach their "ideal" values in a synchronous fashion when vdW parameters are tuned. The following hydration free energy calculation using thermodynamic integration further justifies the vdW refinement. We conclude that simple vdW parameterization can significantly reduce the prediction errors. We believe that GAFF can greatly improve its performance in predicting liquid properties of organic molecules after a systematic vdW parameterization, which will be reported in a separate paper.

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

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

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

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

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

[6]  S. Harvey,et al.  The flying ice cube: Velocity rescaling in molecular dynamics leads to violation of energy equipartition , 1998, J. Comput. Chem..

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

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

[9]  Bernard R Brooks,et al.  Long-range Lennard-Jones and electrostatic interactions in interfaces: application of the isotropic periodic sum method. , 2007, The journal of physical chemistry. B.

[10]  Wilfred F. van Gunsteren,et al.  Validation of molecular dynamics simulation , 1998 .

[11]  Bernard R Brooks,et al.  Isotropic periodic sum: a method for the calculation of long-range interactions. , 2005, The Journal of chemical physics.

[12]  Charles L. Brooks,et al.  Detailed considerations for a balanced and broadly applicable force field: A study of substituted benzenes modeled with OPLS‐AA , 2005, J. Comput. Chem..

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

[14]  D. J. Price,et al.  A modified TIP3P water potential for simulation with Ewald summation. , 2004, The Journal of chemical physics.

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

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

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

[18]  Berend Smit,et al.  Simulating the critical behaviour of complex fluids , 1993, Nature.

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

[20]  Tingjun Hou,et al.  Aqueous Solubility Prediction Based on Weighted Atom Type Counts and Solvent Accessible Surface Areas , 2009, J. Chem. Inf. Model..

[21]  Jussi T. S. Heikkilä,et al.  Enthalpies of Formation of Cyclic Acetals. 1,3-Dioxolane, 2-Methyl-1,3-Dioxolane, and 2,4-Dimethyl-1,3-Dioxolanes. , 1969 .

[22]  B. Kuhn,et al.  Validation and use of the MM-PBSA approach for drug discovery. , 2005, Journal of medicinal chemistry.

[23]  J. Berg,et al.  Molecular dynamics simulations of biomolecules , 2002, Nature Structural Biology.

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

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

[26]  Tingjun Hou,et al.  Assessing the Performance of the MM/PBSA and MM/GBSA Methods. 1. The Accuracy of Binding Free Energy Calculations Based on Molecular Dynamics Simulations , 2011, J. Chem. Inf. Model..

[27]  H. C. Andersen.  Molecular dynamics simulations at constant pressure and/or temperature , 1980 .

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

[29]  Haiyan Liu,et al.  Molecular dynamics simulations of liquid methanol and methanol–water mixtures with polarizable models , 2006, J. Comput. Chem..

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

[31]  J. D. Kemp,et al.  Hindered Rotation of the Methyl Groups in Propane. The Heat Capacity, Vapor Pressure, Heats of Fusion and Vaporization of Propane. Entropy and Density of the Gas , 1938 .

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

[33]  Peter A. Kollman,et al.  Application of the RESP Methodology in the Parametrization of Organic Solvents , 1998 .

[34]  Alexander D. MacKerell,et al.  Re-evaluation of the reported experimental values of the heat of vaporization of N-methylacetamide. , 2008, Journal of chemical theory and computation.

[35]  Alexander D. MacKerell,et al.  All‐atom empirical force field for nucleic acids: I. Parameter optimization based on small molecule and condensed phase macromolecular target data , 2000 .

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

[37]  Araz Jakalian,et al.  Fast, efficient generation of high‐quality atomic charges. AM1‐BCC model: I. Method , 2000 .

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

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

[40]  Christopher I. Bayly,et al.  Fast, efficient generation of high‐quality atomic charges. AM1‐BCC model: II. Parameterization and validation , 2002, J. Comput. Chem..

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

[42]  A. Shimizu,et al.  Enthalpies of Vaporization of Organic Compounds. VII. Some Carboxylic Acids. , 1970 .

[43]  W. L. Jorgensen,et al.  Temperature dependence of TIP3P, SPC, and TIP4P water from NPT Monte Carlo simulations: Seeking temperatures of maximum density , 1998 .

[44]  Peter A. Kollman,et al.  Application of the multimolecule and multiconformational RESP methodology to biopolymers: Charge derivation for DNA, RNA, and proteins , 1995, J. Comput. Chem..

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

[46]  Gerhard Hummer,et al.  Calculation of free‐energy differences from computer simulations of initial and final states , 1996 .

[47]  Chris Oostenbrink,et al.  An improved nucleic acid parameter set for the GROMOS force field , 2005, J. Comput. Chem..

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

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

[50]  W. V. Gunsteren,et al.  Can the density maximum of water be found by computer simulation , 1994 .

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

[52]  L. Crocker,et al.  Thermochemical studies of carbonyl compounds. 5. Enthalpies of reduction of carbonyl groups , 1991 .

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

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

[55]  William C. Swope,et al.  The role of long ranged forces in determining the structure and properties of liquid water , 1983 .

[56]  Wilfred F. van Gunsteren,et al.  Development of a simple, self-consistent polarizable model for liquid water , 2003 .

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

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

[59]  Tetsu Narumi,et al.  Cutoff radius effect of isotropic periodic sum method for transport coefficients of Lennard-Jones liquid. , 2007, The Journal of chemical physics.

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