Dataset Information

Acceptable protein and solvent behavior in primary hydration shell simulations of hen lysozyme.

ABSTRACT: The "primary hydration shell" method in molecular dynamics simulations uses a two- to three-layer thick shell of explicitly represented water molecules as the solvent around the protein of interest. We show that despite its simplicity, this computationally cheap model is capable of predicting acceptable water and protein behavior using the CHARMM22/CMAP potential function. For protein dynamics, comparisons are made with Lipari-Szabo order parameters. These have been derived from NMR relaxation parameters for pico-nano second motions of the NH groups in the main-chain and NH(2) groups in Asn/Gln side chains in hen lysozyme. It is also shown that an even simpler, and therefore faster, water-shell model leads to results in similarly good agreement with experiments, and also compared with simulations using a full box of water with periodic boundary conditions or with an implicit solvation model. Thus, the primary hydration shell method should be useful in making larger systems accessible to extensive simulations.

SUBMITTER: Hamaneh MB

PROVIDER: S-EPMC1864839 | biostudies-literature | 2007 Apr

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Acceptable protein and solvent behavior in primary hydration shell simulations of hen lysozyme.

Hamaneh Mehdi Bagheri MB Buck Matthias M

Biophysical journal 20070126 7

The "primary hydration shell" method in molecular dynamics simulations uses a two- to three-layer thick shell of explicitly represented water molecules as the solvent around the protein of interest. We show that despite its simplicity, this computationally cheap model is capable of predicting acceptable water and protein behavior using the CHARMM22/CMAP potential function. For protein dynamics, comparisons are made with Lipari-Szabo order parameters. These have been derived from NMR relaxation p ...[more]

PMID: 17259273

Similar Datasets

Project description:The solvation layer surrounding a protein is clearly an intrinsic part of protein structure-dynamics-function, and our understanding of how the hydration dynamics influences protein function is emerging. We have recently reported simulations indicating a correlation between regional hydration dynamics and the structure of the solvation layer around different regions of the enzyme Candida antarctica lipase B, wherein the radial distribution function (RDF) was used to calculate the pairwise entropy, providing a link between dynamics (diffusion) and thermodynamics (excess entropy) known as Rosenfeld scaling. Regions with higher RDF values/peaks in the hydration layer (the first peak, within 6 Å of the protein surface) have faster diffusion in the hydration layer. The finding thus hinted at a handle for rapid evaluation of hydration dynamics at different regions on the protein surface in molecular dynamics simulations. Such an approach may move the analysis of hydration dynamics from a specialized venture to routine analysis, enabling an informatics approach to evaluate the role of hydration dynamics in biomolecular function. This paper first confirms that the correlation between regional diffusive dynamics and hydration layer structure (via water center of mass around protein side-chain atom RDF) is observed as a general relationship across a set of proteins. Second, it seeks to devise an approach for rapid analysis of hydration dynamics, determining the minimum amount of information and computational effort required to get a reliable value of hydration dynamics from structural data in MD simulations based on the protein-water RDF. A linear regression model using the integral of the hydration layer in the water-protein RDF was found to provide statistically equivalent apparent diffusion coefficients at the 95% confidence level for a set of 92 regions within five different proteins. In summary, RDF analysis of 10 ns of data after simulation convergence is sufficient to accurately map regions of fast and slow hydration dynamics around a protein surface. Additionally, it is anticipated that a quick look at protein-water RDFs, comparing peak heights, will be useful to provide a qualitative ranking of regions of faster and slower hydration dynamics at the protein surface for rapid analysis when investigating the role of solvent dynamics in protein function.

Project description:We demonstrate that the effect of protein crowding is critically dependent on the stability of the protein's hydration shell, which can dramatically vary between different proteins. In the human eye lens, γS-crystallin (γS-WT) forms a densely packed transparent hydrogel with a high refractive index, making it an ideal system for studying the effects of protein crowding. A single point mutation generates the cataract-related variant γS-G18V, dramatically altering the optical properties of the eye lens. This system offers an opportunity to explore fundamental questions regarding the effect of protein crowding, using γS-WT and γS-G18V: (i) how do the diffusion dynamics of hydration water change as a function of protein crowding?; and (ii) upon hydrogel formation of γS-WT, has a dynamic transition occurred generating a single population of hydration water, or do populations of bulk and hydration water coexist? Using localized spin probes, we separately probe the local translational diffusivity of both surface hydration and interstitial water of γS-WT and γS-G18V in solution. Surprisingly, we find that under the influence of hydrogel formation at highly crowded γS-WT concentrations up to 500 mg/mL, the protein hydration shell remains remarkably dynamic, slowing by less than a factor of 2, if at all, compared to that in dilute protein solutions of ∼5 mg/mL. Upon self-crowding, the population of this robust surface hydration water increases, while a significant bulk-like water population coexists even at ∼500 mg/mL protein concentrations. In contrast, surface water of γS-G18V irreversibly dehydrates with moderate concentration increases or subtle alterations to the solution conditions, demonstrating that the effect of protein crowding is highly dependent on the stability of the protein-specific hydration shell. The core function of γS-crystallin in the eye lens may be precisely its capacity to preserve a robust hydration shell, whose stability is abolished by a single G18V mutation.

Dataset Information

Acceptable protein and solvent behavior in primary hydration shell simulations of hen lysozyme.

Publications

Acceptable protein and solvent behavior in primary hydration shell simulations of hen lysozyme.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets