Application of Molecular Dynamics Simulations in ... - ACS Publications

May 30, 2011 - similar MD protocol to the above liquid MD simulation was applied to do the ..... minimization. Fourth, protocol 4, with the weak coupl...
1 downloads 0 Views 3MB Size
ARTICLE pubs.acs.org/JCTC

Application of Molecular Dynamics Simulations in Molecular Property Prediction. 1. Density and Heat of Vaporization Junmei Wang*,† and Tingjun Hou*,‡ †

Department of Pharmacology, University of Texas Southwestern Medical Center at Dallas, 5323 Harry Hines Boulevard, Dallas, Texas 75390-9050, United States ‡ Institute of Nano & Soft Materials (FUNSOM) and Jiangsu Key Laboratory for Carbon-Based Functional Materials & Devices, Soochow University, Suzhou, Jiangsu 215123, P. R. China

bS Supporting Information ABSTRACT: Molecular mechanical force field (FF) methods are useful in studying condensed phase properties. They are complementary to experiments and can often go beyond experiments in atomic details. Even if 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 the 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) parametrization has been discussed. An encouraging finding in vdW parametrization 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 parametrization 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 parametrization, which will be reported in a separate paper.

1. INTRODUCTION Molecular dynamics simulation has become increasingly important in studying the structures, dynamics, and functions of biomolecules. 1 Nevertheless, it is also an important tool for predicting a variety of molecular properties. 2 The power of MD simulations lies in that it can provide atomic details on the dynamics of a simulation system under conditions that are beyond the scope of experiments. Even for some specific force fields, such as CHARMM, 3,4 AMBER,5,6 OPLS, 7 and GROMOS, 8,9 which were developed mainly to study the biomolecules, it is often critical for them to reproduce the bulk properties of some small molecules that represent the moieties of biomolecules. With the continually increased demand of “ideal” molecular mechanical models to study ligandprotein and ligandnucleic acid interactions, it is important to develop a general purpose force field that is totally compatible with one or a set of biomolecular force fields for studying the interactions between biomolecules and organic molecules. The General AMBER force field (GAFF)10 was developed in this spirit to study biomoleculeligand interactions in conjunction with the AMBER force r 2011 American Chemical Society

fields. 5,6 Vpotential ¼

Kr ðr  req Þ2 þ ∑ Kθ ðθ  θeq Þ2 ∑ bonds angles





Vn ½1 þ cosðnϕ  γÞþ þ i