Hydration thermodynamic properties of amino acid analogues: a systematic comparison of biomolecular force fields and water models

We present an extensive study on hydration thermodynamic properties of analogues of 13 amino acid side chains at 298 K and 1 atm. The hydration free energies DeltaG, entropies DeltaS, enthalpies DeltaH, and heat capacities Deltac(P)() were determined for 10 combinations of force fields and water models. The statistical sampling was extended such that precisions of 0.3, 0.8, 0.8 kJ/mol and 25 J/(mol K) were reached for DeltaG, TDeltaS, DeltaH, and Deltac(P)(), respectively. The three force fields used in this study are AMBER99, GROMOS 53A6, and OPLS-AA; the five water models are SPC, SPC/E, TIP3P, TIP4P, and TIP4P-Ew. We found that the choice of water model strongly influences the accuracy of the calculated hydration entropies, enthalpies, and heat capacities, while differences in accuracy between the force fields are small. On the basis of an analysis of the hydrophobic analogues of the amino acid side chains, we discuss what properties of the water models are responsible for the observed discrepancies between computed and experimental values. The SPC/E water model performs best with all three biomolecular force fields.     

Prediction of self-diffusion coefficient and shear viscosity of water and its binary mixtures with methanol and ethanol by molecular simulation

Density, self-diffusion coefficient, and shear viscosity of pure liquid water are predicted for temperatures between 280 and 373 K by molecular dynamics simulation and the Green-Kubo method. Four different rigid nonpolarizable water models are assessed: SPC, SPC/E, TIP4P, and TIP4P/2005. The pressure dependence of the self-diffusion coefficient and the shear viscosity for pure liquid water is also calculated and the anomalous behavior of these properties is qualitatively well predicted. Furthermore, transport properties as well as excess volume and excess enthalpy of aqueous binary mixtures containing methanol or ethanol, based on the SPC/E and TIP4P/2005 water models, are calculated. Under the tested conditions, the TIP4P/2005 model gives the best quantitative and qualitative agreement with experiments for the regarded transport properties. The deviations from experimental data are of 5% to 15% for pure liquid water and 5% to 20% for the water + alcohol mixtures. Moreover, the center of mass power spectrum of water as well as the investigated mixtures are analyzed and the hydrogen-bonding structure is discussed for different states.     

Influence of bond flexibility on the vapor-liquid phase equilibria of water

The authors performed Gibbs ensemble simulations on the vapor-liquid equilibrium of water to investigate the influence of incorporating intramolecular degrees of freedom in the simple point charge (SPC) water model. Results for vapor pressures, saturation densities, heats of vaporization, and the critical point for two different flexible models are compared with data for the corresponding rigid SPC and SPC/E models. They found that the introduction of internal vibrations, and also their parametrization, has an observable effect on the prediction of the vapor-liquid coexistence curve. The flexible SPC/Fw model, although optimized to describe bulk diffusion and dielectric constants at ambient conditions, gives the best prediction of saturation densities and the critical point of the examined models.     

GEMC和GDI方法计算流体气液相平衡的比较

考察采用TraPPE联合原子和OPLS全原子力场两种分子力场, Gibbs系综蒙特卡罗(GEMC)方法和Gibbs-Duhem积分(GDI)方法计算流体气液相平衡的适用性、计算速度、计算精度等问题. 结果表明, 在采用全原子力场情况下, GDI方法比GEMC方法极大地节省了计算时间. 从计算结果来看, 两种方法各有适用范围, 在使用时可互为补充. 在给定力场的前提下, 两种方法所得到的液相密度、蒸发焓、临界温度和临界密度相差不大, 而当力场中的缺陷导致蒸发焓的计算不够准确时, 两种计算方法得到的气体的压力和密度明显不同,进而导致预测的临界压力也明显不同.     

Experimental and Molecular Simulation of Volumetric Properties of Methyl Nonanoate,n-Dodecane,and Their Binary Mixtures

Densities of methyl nonanoate,n-dodecane,and their binary mixtures were investigated to provide the necessary data for their engineering applications as promising fuels and fuel additives.In the present work,densities were measured under atmospheric pressure at 293.15—463.15 K.The density data for the binary mixtures were fitted into a forai of excess molar volume.The excess molar volumes were mostly positive,and the maximum value was obtained at molar fractions of n-dodecane between 0.5 and 0.6.Molecular simulations of specified systems were carried out by using four kinds of force fields,and the suitable force fields for describing the volume properties of the system were AMBER96 and OPLS-AA.The relative deviations for these two force fields between the simulated and the experimental data were well within±4%,which meets the general engineering requirement.     

Solvation free energies of amino acid side chain analogs for common molecular mechanics water models

Quantitative free energy computation involves both using a model that is sufficiently faithful to the experimental system under study (accuracy) and establishing statistically meaningful measures of the uncertainties resulting from finite sampling (precision). In order to examine the accuracy of a range of common water models used for protein simulation for their solute/solvent properties, we calculate the free energy of hydration of 15 amino acid side chain analogs derived from the OPLS-AA parameter set with the TIP3P, TIP4P, SPC, SPC/E, TIP3P-MOD, and TIP4P-Ew water models. We achieve a high degree of statistical precision in our simulations, obtaining uncertainties for the free energy of hydration of 0.02-0.06 kcal/mol, equivalent to that obtained in experimental hydration free energy measurements of the same molecules. We find that TIP3P-MOD, a model designed to give improved free energy of hydration for methane, gives uniformly the closest match to experiment; we also find that the ability to accurately model pure water properties does not necessarily predict ability to predict solute/solvent behavior. We also evaluate the free energies of a number of novel modifications of TIP3P designed as a proof of concept that it is possible to obtain much better solute/solvent free energetic behavior without substantially negatively affecting pure water properties. We decrease the average error to zero while reducing the root mean square error below that of any of the published water models, with measured liquid water properties remaining almost constant with respect to our perturbations. This demonstrates there is still both room for improvement within current fixed-charge biomolecular force fields and significant parameter flexibility to make these improvements. Recent research in computational efficiency of free energy methods allows us to perform simulations on a local cluster that previously required large scale distributed computing, performing four times as much computational work in approximately a tenth of the computer time as a similar study a year ago.     

Assessment of biomolecular force fields for molecular dynamics simulations in a protein crystal

Different biomolecular force fields (OPLS‐AA, AMBER03, and GROMOS96) in conjunction with SPC, SPC/E and TIP3P water models are assessed for molecular dynamics simulations in a tetragonal lysozyme crystal. The root mean square deviations for the C_a atoms of lysozymes are about 0.1 to 0.2 nm from OPLS‐AA and AMBER03, smaller than 0.4 nm from GROMOS96. All force fields exhibit similar pattern in B‐factors, whereas OPLS‐AA and AMBER03 accurately reproduce experimental measurements. Despite slight variations, the primary secondary structures are well conserved using different force fields. Water diffusion in the crystal is approximately ten‐fold slower than in bulk phase. The directional and average water diffusivities from OPLS‐AA and AMBER03 along with SPC/E model match fairly well with experimental data. Compared to GROMOS96, OPLS‐AA and AMBER03 predict larger hydrophilic solvent‐accessible surface area of lysozyme, more hydrogen bonds between lysozyme and water, and higher percentage of water in hydration shell. SPC, SPC/E and TIP3P water models have similar performance in most energetic and structural properties, but SPC/E outperforms in water diffusion. While all force fields overestimate the mobility and electrical conductivity of NaCl, a combination of OPLS‐AA for lysozyme and the Kirkwood‐Buff model for ions is superior to others. As attributed to the steric restraints and surface interactions, the mobility and conductivity in the crystal are reduced by one to two orders of magnitude from aqueous solution. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Detailed considerations for a balanced and broadly applicable force field: a study of substituted benzenes modeled with OPLS-AA

Modern classical force fields have been traditionally parameterized by attempting to maximize agreement to any number of experimental and/or quantum mechanical target properties. As these force fields are pushed towards obtaining quantitative estimates of often subtle energetic differences, stringent and consistent parameterization criteria, particularly in regard to charge distributions, are required to ensure that systematic errors cancel, that parameters are transferable between molecules, and that performance does not significantly deteriorate when using more approximate methods, such as with continuum solvent models. Relative free energies of hydration are presented here for 40 mono- and disubstituted benzenes modeled with the OPLS-AA force field; heats of vaporization and pure liquid densities at standard conditions are presented when experimental data is available. Overall agreement between OPLS-AA and experiment is remarkable (average error = 0.5 kcal/mol for DeltaDeltaG(hydration), 1.0 kcal/mol for DeltaH(vap) (0), 0.02 g/mL for densities), yet several functional groups are identified as having consistent and correctable errors (alkyl-, nitro-, and thiobenzenes). Relative free energies of hydration obtained with rigorous free energy perturbations using explicit solvent are also compared with energies from minimizations using a generalized Born model (GB). There is high correlation between these estimates (R = 0.99), and as demonstrated here, reparameterization of the aforementioned groups can be guided with rapid GB calculations.     

An advanced Gibbs-Duhem integration method: theory and applications

The conventional Gibbs-Duhem integration method is very convenient for the prediction of phase equilibria of both pure components and mixtures. However, it turns out to be inefficient. The method requires a number of lengthy simulations to predict the state conditions at which phase coexistence occurs. This number is not known from the outset of the numerical integration process. Furthermore, the molecular configurations generated during the simulations are merely used to predict the coexistence condition and not the liquid- and vapor-phase densities and mole fractions at coexistence. In this publication, an advanced Gibbs-Duhem integration method is presented that overcomes above-mentioned disadvantage and inefficiency. The advanced method is a combination of Gibbs-Duhem integration and multiple-histogram reweighting. Application of multiple-histogram reweighting enables the substitution of the unknown number of simulations by a fixed and predetermined number. The advanced method has a retroactive nature; a current simulation improves the predictions of previously computed coexistence points as well. The advanced Gibbs-Duhem integration method has been applied for the prediction of vapor-liquid equilibria of a number of binary mixtures. The method turned out to be very convenient, much faster than the conventional method, and provided smooth simulation results. As the employed force fields perfectly predict pure-component vapor-liquid equilibria, the binary simulations were very well suitable for testing the performance of different sets of combining rules. Employing Lorentz-Hudson-McCoubrey combining rules for interactions between unlike molecules, as opposed to Lorentz-Berthelot combining rules for all interactions, considerably improved the agreement between experimental and simulated data.     

A general perturbation approach for equation of state development: applications to simple fluids, ab initio potentials, and fullerenes

A new perturbation scheme based on the Barker-Henderson perturbation theory [J. Chem. Phys. 47, 4714 (1967)] is proposed to predict the thermodynamic properties of spherical molecules. Accurate predictions of second virial coefficients and vapor-liquid coexistence properties are obtained for a large variety of potential functions (square well, Yukawa, Sutherland, Lennard-Jones, Buckingham, Girifalco). New Gibbs ensemble Monte Carlo simulations of the generalized exp-m Buckingham potential are reported. An extension of the perturbation approach to mixtures is proposed, and excellent predictions of vapor-liquid equilibria are obtained for Lennard-Jones mixtures. The perturbation scheme can be applied to complex potential functions fitted to ab initio data to predict the properties of real molecules such as neon. The new approach can also be used as an auxiliary tool in molecular simulation studies, to efficiently optimize an intermolecular potential on macroscopic properties or match force fields based on different potential functions.     

Effects of water model and simulation box size on protein diffusional motions

We performed molecular dynamics simulations of ubiquitin with the distinct water models, TIP3P, SPC, SPC/E, and SPC/Fw, in different system sizes with different box shapes. The translational diffusion constants of pure water linearly depend on the effective box length, which is known as finite size effect, whereas the first and second rotational times of pure water are nearly constant. We then observed that both the overall translation and rotational motions of the protein are linearly correlated to the viscosity of pure water. As expected from the finite size effect in pure water, the translational diffusion of the protein is significantly affected by the system size, and rotational diffusion is nearly size-independent. After correction for the finite size effect, the SPC/E and SPC/Fw models reproduce both the translational and rotational motion of the protein relatively well. Thus, water models that reproduce the experimentally derived diffusional properties of pure water accurately are expected to also be suitable for simulating protein diffusion quantitatively.     

Interatomic Lennard-Jones potentials of linear and branched alkanes calibrated by Gibbs ensemble simulations for vapor-liquid equilibria

We propose Lennard-Jones potential parameters for interatomic interactions of linear and branched alkanes based on matching the results of Gibbs ensemble simulations of vapor-liquid equilibria to experimental data. The alkane model is similar as in the OPLS-AA, but multiple atom types for carbon based on the number of covalently bonded hydrogen atoms are necessary to accurately reproduce liquid densities and enthalpies of vaporization with the errors of 2.1% and 3.3%, respectively, for hydrocarbons of various chain lengths and structures. We find that the attraction energies of the carbon atoms are almost proportional to the number of covalent hydrogen atoms with each increasing the carbon energy parameter by approximately 0.033 kcal/mol. Though the present force field outperforms the OPLS-AA force field for alkanes we studied, systematic deviations for vapor pressures are still observed with errors of 15%-30%, and critical temperatures are slightly underestimated. We think that these shortcomings are probably due to the inadequacy of the two-parameter Lennard-Jones potential, and especially its behavior at short distances.     

Simulation of phase separation in alcohol/water mixtures using two-body force field and standard molecular dynamics

Standard molecular dynamics simulations have been carried out on pure alcohols and alcohol/water mixtures. A simple atom-atom force field consisting of Lennard-Jones potentials plus coulombic terms over atomic point charges, but without explicit polarization terms, has been specifically fitted to reproduce several experimental properties of the pure alcohols, and has been used for mixtures by developing combination rules with the TIP3P water model. Densities, enthalpies of vaporization, radial distribution functions, self-diffusion coefficients, and rotational correlation functions of the pure alcohols are well reproduced and compare favorably with those from more sophisticated force fields. Some key aspects of the phase behaviour are correctly reproduced by the molecular dynamics simulation, showing a distinct demixing process for the n-butanol/water mixture as opposed to the stability of the t-butanol/water mixtures. The results demonstrate the ability of a molecular dynamics simulation, even in its standard form and with easily accessible time ranges, but with a carefully optimized force field, to simulate and, to a certain extent, predict the properties of binary mixtures.     

Molecular dynamics simulations of methane hydrate using polarizable force fields

Molecular dynamics simulations of methane hydrate have been carried out using the polarizable AMOEBA and COS/G2 force fields. Properties calculated include the temperature dependence of the lattice constant, the OC and OO radial distribution functions, and the vibrational spectra. Both the AMOEBA and COS/G2 force fields are found to successfully account for the available experimental data, with overall somewhat better agreement with experiment being found for the AMOEBA model. Comparison is made with previous results obtained using TIP4P and SPC/E effective two-body force fields and the polarizable TIP4P-FQ force field, which allows for in-plane polarization only. Significant differences are found between the properties calculated using the TIP4P-FQ model and those obtained using the other models, indicating an inadequacy of restricting explicit polarization to in-plane only.     

Dynamical properties of the soft sticky dipole-quadrupole-octupole water model: a molecular dynamics study

The dynamical properties of the soft sticky dipole-quadrupole-octupole (SSDQO) water model using SPC/E moments are calculated utilizing molecular dynamics simulations. This new potential for liquid water describes the water-water interactions by a Lennard-Jones term and a sticky potential, which is an approximate moment expansion with point dipole, quadrupole, and octupole moments, and reproduces radial distribution functions of pure liquid water using the moments of SPC/E [Ichiye and Tan, J. Chem. Phys. 124, 134504 (2006)]. The forces and torques of SSDQO water for the dipole-quadrupole, quadrupole-quadrupole, and dipole-octupole interactions are derived here. The simulations are carried out at 298 K in the microcanonical ensemble employing the Ewald method for the long-range dipole-dipole interactions. Here, various dynamical properties associated with translational and rotational motions of SSDQO water using the moments of SPC/E (SSDQO:SPC/E) water are compared with the results from SPC/E and also experiment. The self-diffusion coefficient of SSDQO:SPC/E water is found to be in excellent agreement with both SPC/E and experiment whereas the single particle orientational relaxation time for dipole vector is better than SPC/E water but it is somewhat smaller than experiment. The dielectric constant of SSDQO:SPC/E is essentially identical to SPC/E, and both are slightly lower than experiment. Also, molecular dynamics simulations of the SSDQO water model are found to be about twice as fast as three-site models such as SPC/E.     

A modified clausius equation of state for calculation of multicomponent refrigerant vapor-liquid equilibria

Gow, A.S., 1993. A modified Clausius equation of state for calculation of multicomponent refrigerant vapor-liquid equilibria. Fluid Phase Equilibria, 90: 219-249.

A modified Clausius equation of state with a single temperature dependent energy-volume parameter a(T) in the attractive term was designed to describe the vapor pressure vs. temperature relationship of 39 pure refrigerant fluids including elementary cryogenic materials (e.g. He, Ar, N₂, CO₂, CH₄, etc.), chlorofluorocarbons (CFCs), hydrofluorocarbons (HFCs), hydrochlorofluorocarbons (HCFCs), fluorocarbons (FCs), and various other simple cryogenic compounds. The equation developed represents the vapor-liquid coexistence dome, and the superheated vapor compressibility factor and enthalpy for pure refrigerants.

The vapor-liquid equilibrium for refrigerant mixtures is calculated using a “phi-phi” method with “one fluid” van der Waals mixing and combining rules for the equation of state parameters a_M(T), b_M and c_M. A single interaction constant k₁₂ is used to describe non-ideal behavior of each binary. The binary interaction constant, which is a strong function of temperature, and the sign of which signifies the type of deviations from Raoult's law, is obtained by correlating experimental bubble point data for isothermal binary refrigerant mixtures. The proposed equation of state generally describes binary P-x,y data more accurately the higher the temperature for a given system. The method presented is extended to predict vapor-liquid equilibria for the R14-R23-R13 ternary system at 198.75 K using binary interaction constants at this temperature for the three binaries involved.     

Environmentally Friendly Green Materials from Plant-Based Resources: Modification of Soy Protein using Gellan and Micro/Nano-Fibrillated Cellulose

Fully biodegradable micro/nano-composite resins were prepared by reinforcing soy protein concentrate (SPC) with micro/nano-fibrillated cellulose (MFC) and then blending with gellan. The composite resins showed excellent mechanical and physical properties under testing conditions. Due to the high aspect ratio of MFC, excellent mechanical properties of MFC and MFC/SPC interfacial properties, the SPC (100 parts) reinforced with glycerol (1.5 parts) and MFC (40 parts) showed fracture stress of 88.2 MPa and Young's modulus of about 4.1 GPa, which are higher than those of many conventional petroleum-based plastics. MFC reinforced SPC composite resins were then further modified by blending with gellan to obtain further improvement in fracture stress and Young's modulus. SPC resin containing glycerol (1.5 parts), gellan (40 parts) and MFC (40 parts) had fracture stress of over 122 MPa and Young's modulus of about 5.8 GPa. Although the moisture sensitivity of the specimens was high, they have the potential to replace petroleum-based materials in many fields, particularly for indoor applications.