Accurate Calculation of Hydration Free Energies using Pair-Specific Lennard-Jones Parameters in the CHARMM Drude Polarizable Force Field

Abstract

Lennard-Jones (LJ) parameters for a variety of model compounds have previously been optimized within the CHARMM Drude polarizable force field to reproduce accurately pure liquid phase thermodynamic properties as well as additional target data. While the polarizable force field resulting from this optimization procedure has been shown to satisfactorily reproduce a wide range of experimental reference data across numerous series of small molecules, a slight but systematic overestimate of the hydration free energies has also been noted. Here, the reproduction of experimental hydration free energies is greatly improved by the introduction of pair-specific LJ parameters between solute heavy atoms and water oxygen atoms that override the standard LJ parameters obtained from combining rules. The changes are small and a systematic protocol is developed for the optimization of pair-specific LJ parameters and applied to the development of pair-specific LJ parameters for alkanes, alcohols and ethers. The resulting parameters not only yield hydration free energies in good agreement with experimental values, but also provide a framework upon which other pair-specific LJ parameters can be added as new compounds are parametrized within the CHARMM Drude polarizable force field. Detailed analysis of the contributions to the hydration free energies reveals that the dispersion interaction is the main source of the systematic errors in the hydration free energies. This information suggests that the systematic error may result from problems with the LJ combining rules and is combined with analysis of the pair-specific LJ parameters obtained in this work to identify a preliminary improved combining rule.

A GROMOS-Compatible Force Field for Small Organic Molecules in the Condensed Phase: The 2016H66 Parameter Set

“…Condensed-phase force fields such as CHARMM [5][6][7][8][9][10][11][12][13][14][15] , AMBER [16][17][18][19][20][21][22][23][24][25][26][27][28][29][30][31] , OPLS 32-44 , TraPPE [45][46][47][48][49][50][51][52] and GROMOS [53][54][55][56][57][58][59][60][61][62][63][64][65][66][67][68] mainly aim at the description of solids, liquids and solutions, including solvated biomolecules as a particularly relevant special case. ...…”

Evaluation of solvation free energies for small molecules with the AMOEBA polarizable force field

“…Commonly, solvation free energy calculation approaches20, 21, 22, 23, 24, 25, 26 have been performed to assess force field properties. This is thanks to the availability of high accuracy experimental data, and the straightforward computational methodologies for free energy prediction.…”

The polarizing forces of water

“…However, it should be noted that the goal has been to efficiently introduce polarization into large biological simulations, for example, proteins, which are typically carried out at ambient conditions. The models now encompass not simply water, but other molecules like alkanes and aromatic compounds [48][49][50].…”