Dataset Information

Assimilating Radial Distribution Functions To Build Water Models with Improved Structural Properties.

ABSTRACT: The structural properties of three- and four-site water models are improved by extending the ForceBalance parametrization code to include a new methodology allowing for the targeting of any radial distribution function (RDF) during the parametrization of a force field. The mean squared difference (MSD) between the experimental and simulated RDFs contributes to an objective function, allowing for the systematic optimization of force field parameters to reach closer overall agreement with experiment. RDF fitting is applied to develop modified versions of the TIP3P and TIP4P/2005 water models in which the Lennard-Jones potential is replaced by a Buckingham potential. The optimized TIP3P-Buckingham and TIP4P-Buckingham potentials feature 93 and 98% lower MSDs in the OO RDF compared to the TIP3P and TIP4P/2005 models respectively, with marked decreases in the height of the first peak. Additionally, these Buckingham models predict the entropy of water more accurately, reducing the error in the entropy of TIP3P from 11 to 3% and the error in the entropy of TIP4P/2005 from 11 to 2%. These new Buckingham models have improved predictive power for many nonfitted properties particularly in the case of TIP3P. Our work directly demonstrates how the Buckingham potential can improve the description of water's structural properties beyond the Lennard-Jones potential. Moreover, adding a Buckingham potential is a favorable alternative to adding interaction sites in terms of computational speed on modern GPU hardware.

SUBMITTER: Wade AD

PROVIDER: S-EPMC6175054 | biostudies-literature | 2018 Sep

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Assimilating Radial Distribution Functions To Build Water Models with Improved Structural Properties.

Wade Alexander D AD Wang Lee-Ping LP Huggins David J DJ

Journal of chemical information and modeling 20180830 9

The structural properties of three- and four-site water models are improved by extending the ForceBalance parametrization code to include a new methodology allowing for the targeting of any radial distribution function (RDF) during the parametrization of a force field. The mean squared difference (MSD) between the experimental and simulated RDFs contributes to an objective function, allowing for the systematic optimization of force field parameters to reach closer overall agreement with experime ...[more]

PMID: 30113842

Similar Datasets

Project description:Water plays a critical role in the structure and function of proteins, although the experimental properties of water around protein structures are not well understood. The water can be classified by the separation from the protein surface into bulk water and hydration water. Hydration water interacts closely with the protein and contributes to protein folding, stability, and dynamics, as well as interacting with the bulk water. Water potential functions are often parametrized to fit bulk water properties because of the limited experimental data for hydration water. Therefore, the structural and energetic properties of the hydration water were assessed for 105 atomic resolution (<or=1.0 A) protein crystal structures with a high level of hydration water by calculating the experimental water-protein radial distribution function or surface distribution function (SDF) and water radial distribution function (RDF). Two maxima are observed in SDF: the first maximum at a radius of 2.75 A reflects first shell and hydrogen bond interactions between protein and water, and the second maximum at 3.65 A reflects second shell and van der Waals interactions between water and nonpolar atoms of protein-forming clathrate-hydrate-like structures. Thus, the two shells do not overlap. The RDF showed the features of liquid water rather than solid ice. The first and second maxima of RDF at 2.75 and 4.5 A, respectively, are the same as for bulk water, but the peaks are sharper, indicating hydration water is more stable than bulk water. Both distribution functions are inversely correlated with the distribution of B factors (atomic thermal factors) for the waters, suggesting that the maxima reflect stable positions. Therefore, the average water structure near the protein surface has experimentally observable differences from bulk water. This analysis will help improve the accuracy for models of water on the protein surface by providing rigorous data for the effects of the apparent chemical potential of the water near a protein surface.

Dataset Information

Assimilating Radial Distribution Functions To Build Water Models with Improved Structural Properties.

Publications

Assimilating Radial Distribution Functions To Build Water Models with Improved Structural Properties.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets