High-density hydration layer of lysozymes: molecular dynamics decomposition of solution scattering data

A characterization of the physical properties of protein hydration water is critical for understanding protein structure and function. Recent small-angle X-ray and neutron scattering data indicate that the density of water on the surface of lysozyme is significantly higher than in bulk water. Here, we provide an interpretation of the scattering results using a molecular dynamics simulation, which allows us to make quantitative predictions about density variations in the first hydration shell. The perturbation relative to bulk water involves statistically significant changes in the average water structure in the first hydration layer. The water density in the first hydration shell is increased by 5% with respect to the bulk. In regions of higher water density, the water dipoles align more parallel to each other and the number of hydrogen bonds per water molecule is higher. Increased water density is found for water molecules interacting with hydrogen and carbon atoms in the backbone or with nonpolar or negatively charged side-chain groups.     

Correlation between the dynamics of hydrogen bonds and the local density reorganization in the protein hydration layer

An atomistic molecular dynamics simulation of the protein villin headpiece subdomain or HP-36 has been carried out with explicit water to explore the microscopic inhomogeneity of local density reorganization of the hydration layers of the three alpha-helical segments of the protein. The density reorganization of the hydration layer of helix-3 is found to occur faster than that for the hydration layers of the other two helices. It is noticed that such inhomogeneous density reorganization at the surface of different secondary structures exhibits excellent correlation with the microscopic dynamics of hydrogen bonds between the protein residues and the hydration water. Further, it is observed that the reorientation of water molecules involved in the formation and breaking of protein-water or water-water hydrogen bonds plays an important role in determining the dynamics of local density of the hydration layer. The faster density reorganization of the hydration layer of helix-3 is also consistent with the functionality of HP-36, as helix-3 contains several active site residues.     

Dielectric relaxation of aqueous solutions of hydrophilic versus amphiphilic peptides

We report on molecular dynamics simulations of the frequency-dependent dielectric relaxation spectra at room temperature for aqueous solutions of a hydrophilic peptide and an amphiphilic peptide at two concentrations. We find that only the high-concentration amphiphilic peptide solution exhibits an anomalous dielectric increment over that of pure water, while the hydrophilic peptide exhibits a significant dielectric decrement. The dielectric component analysis carried out by decomposing these peptide solutions into peptide, hydration layer, and outer layer(s) of water clearly shows the presence of a unique dipolar component with a relaxation time scale on the order of approximately 25 ps (compared to the bulk water time scale of approximately 11 ps) that originates from the interaction between the hydration layer water and the outer layer(s) of water. Results obtained from the dielectric component analysis further show the emergence of a distinct and much lower frequency relaxation process for the high-concentration amphiphilic peptide compared to the hydrophilic peptide due to strong peptide dipolar couplings to all constituents, accompanied by a slowing of the structural relaxation in all water layers, giving rise to time scales close to approximately 1 ns. We suggest that the molecular origin of the dielectric relaxation anomalies is due to frustration in the water network arising from the amphiphilic chemistry of the peptide that does not allow it to reorient on the picosecond time scale of bulk water motions. This explanation is consistent with the idea of the "slaving" of residue side chain motions to protein surface water, and furthermore offers the possibility that the anomalous dynamics observed from a number of spectroscopies arises at the interface of hydrophobic and hydrophilic domains on the protein surface.     

Secondary structure sensitivity of hydrogen bond lifetime dynamics in the protein hydration layer

The heterogeneous nature of a protein surface plays an essential role in its biological activity and molecular recognition, and this role is mediated at least partly through the surrounding water molecules. We have performed atomistic molecular dynamics simulations of an aqueous solution of HP-36 to investigate the correlation between the dynamics of the hydration layer water molecules and the lifetimes of protein-water hydrogen bonds. The nonexponential hydrogen bond lifetime correlation functions have been analyzed by using the formalism of Luzar and Chandler, which allowed identification of the quasi-bound states in the surface and quantification of the dynamic equilibrium between quasi-bound and free water molecules in terms of time-dependent rate of interconversion. It is noticed that, irrespective of the structural heterogeneity of different segments of the protein, namely the three alpha-helices, the positively charged amino acid residues form longer-lived hydrogen bonds with water. The overall relaxation behavior of protein-water hydrogen bonds is found to differ significantly among the three helices of the protein. Study of water number density fluctuation reveals that the hydration layer of helix-3 is much less rigid, which can be correlated with faster structural relaxation of the hydrogen bonds between its residues and water. This also agrees excellently with faster translational and rotational motions of water near helix-3, and hence the lower rigidity of its hydration layer. The lower rigidity of the helix-3 hydration layer also correlates well with the biological activity of the protein, as several of the active-site residues of HP-36 are located in helix-3.     

Low-frequency vibrational spectrum of water in the hydration layer of a protein: a molecular dynamics simulation study

An atomistic molecular dynamics simulation has been carried out to understand the low-frequency intermolecular vibrational spectrum of water present in the hydration layer of the protein villin headpiece subdomain or HP-36. An attempt is made to explore how the heterogeneous rigidity of the hydration layers of different segments (three alpha helices) of the protein, strength of the protein-water hydrogen bonds, and their differential relaxation behavior influence the distribution of the intermolecular vibrational density of states of water in the hydration layers. The calculations revealed that compared to bulk water these bands are nonuniformly blue-shifted for water near the helices, the extent of shifts being more pronounced for water molecules hydrogen bonded to the protein residues. It is further noticed that the larger blue shift observed for the water molecules hydrogen bonded to helix 2 residues correlates excellently with the slowest structural relaxation of these hydrogen bonds. These results can be verified by suitable experimental measurements.     

Experimental estimates of compression heating and decompression cooling in ethylene glycol

The chemical shift difference, Δσ, between the methylene and hydroxyl protons in the high resolution ¹H nuclear magnetic resonance spectrum of ethylene glycol is shown to be pressure dependent. The equilibrium Δσ values for ethylene glycol are reported as a function of temperature and pressure between ambient conditions, 323 K and 2 kbar, respectively. This surface is used along with Δσ values measured in response to a rapid pressure increase to calculate a temperature rise that is used to infer a temperature change for water that is consistent with theoretical estimates. This work implies that compression heating and decompression cooling are not significant enough to interfere with pressure induced protein folding studies.     

Electric-field-controlled water and ion permeation of a hydrophobic nanopore

The permeation of hydrophobic, cylindrical nanopores by water molecules and ions is investigated under equilibrium and out-of-equilibrium conditions by extensive molecular-dynamics simulations. Neglecting the chemical structure of the confining pore surface, we focus on the effects of pore radius and electric field on permeation. The simulations confirm the intermittent filling of the pore by water, reported earlier under equilibrium conditions for pore radii larger than a critical radius R(c). Below this radius, water can still permeate the pore under the action of a strong electric field generated by an ion concentration imbalance at both ends of the pore embedded in a structureless membrane. The water driven into the channel undergoes considerable electrostriction characterized by a mean density up to twice the bulk density and by a dramatic drop in dielectric permittivity which can be traced back to a considerable distortion of the hydrogen-bond network inside the pore. The free-energy barrier to ion permeation is estimated by a variant of umbrella sampling for Na(+), K(+), Ca(2+), and Cl(-) ions, and correlates well with known solvation free energies in bulk water. Starting from an initial imbalance in ion concentration, equilibrium is gradually restored by successive ion passages through the water-filled pore. At each passage the electric field across the pore drops, reducing the initial electrostriction, until the pore, of radius less than R(c), closes to water and hence to ion transport, thus providing a possible mechanism for voltage-dependent gating of hydrophobic pores.     

Thickness of the hydration layer of a protein from molecular dynamics simulation

Water molecules around a protein exhibit slow dynamics with respect to that of pure bulk water. One important issue in protein hydration is the thickness of the hydration layer (i.e., the distance from the protein surface up to which the water dynamics is influenced by the protein). Estimation of thickness is crucial to understand better the properties of "biological water" and the role that it plays in guiding the protein's function. We have performed an atomistic molecular dynamics simulation of an aqueous solution of the protein villin headpiece subdomain or HP-36 to estimate the thickness of its hydration water. In particular, several dynamical properties of water around different segments (three alpha-helices) of the protein have been calculated by varying the thickness of the hydration layers. It is found that in general the influence of the helices on water properties extends beyond the first hydration layer. However, the heterogeneous nature of water among the first hydration layers of the three helices diminishes as the thickness is increased. It indicates that, for a small protein such as HP-36, the thickness of "biological water" is uniform for different segments of the protein.     

Computer simulation of dissociative equilibrium in aqueous NaCl electrolyte with account for polarization and ion recharging. Ionization mechanism

The Monte Carlo method in a system with periodic boundary conditions was used within the model with explicit account for many-bod interactions to calculate ion-water correlation functions and the mean force ion-ion potential for extremely dilute aqueous electrolyte. Many-body interactions result in a decrease in the first coordination number of ions by approximately one molecule. The same effect is observed in the case of hydration in water vapors. Partial displacement of molecules from the lower layer into the higher hydrate layers occurs mainly by means of interactions of dipoles induced on molecules. Many-body interactions enhance the stability of unrecombined ion pairs separated by solvent molecules (SSIP states). The depth of the minimum in the dependence of the ion-ion mean force potential with account for many-body interaction forces is several times higher than in primitive interaction models. The value of effective relative dielectric permeability of the solvent at short distances from the ions grows faster than 1/R. Due to solvent polarization, counterions are strongly repelled at distances corresponding to overlapping of their hydrate shells and are weakly attracted at large distances. Stability of ion pair SSIP states in liquid electrolyte is due to rearrangement of the molecular structure of the solvent in the interion space and is an entropy effect. This mechanism differs qualitatively from that observed under hydration in water vapor and the depth of the minimum corresponding to SSIP states is by an order of magnitude lower in liquid electrolyte as compared to that in saturated water vapor.     

Extracting hydration sites around proteins from explicit water simulations

Two new methods are assessed for determining the location of hydration sites around proteins from computer simulation. Current methods extract hydration sites from peaks in the water density constructed in the protein frame. However, the dynamic nature of the water molecules, the nearby protein residues, and the protein reference frame as a whole tend to smear out the water density, making it more difficult to resolve sites. Two techniques are introduced to better resolve the water density. The first is to construct the water density from the time-averaged position of each water molecule in the protein frame while the water remains within a given distance of this averaged position. The second technique is to construct the water density from the time-averaged position of each water in the reference frame only of the nearby residues. Criteria for determining hydration sites from the water density are examined. Both techniques are found to significantly improve the detail in the water density and the number of hydration sites detected.     

Investigation of the hydration of nonfouling material poly(sulfobetaine methacrylate) by low-field nuclear magnetic resonance

The strong surface hydration layer of nonfouling materials plays a key role in their resistance to nonspecific protein adsorption. Poly(sulfobetaine methacrylate) (polySBMA) is an effective material that can resist nonspecific protein adsorption and cell adhesion. About eight water molecules are tightly bound with one sulfobetaine (SB) unit, and additional water molecules over 8:1 ratio mainly swell the polySBMA matrix, which is obtained through the measurement of T(2) relaxation time by low-field nuclear magnetic resonance (LF-NMR). This result was also supported by the endothermic behavior of water/polySBMA mixtures measured by differential scanning calorimetry (DSC). Furthermore, by comparing both results of polySBMA and poly(ethylene glycol) (PEG), it is found that (1) the hydrated water molecules on the SB unit are more tightly bound than on the ethylene glycol (EG) unit before saturation, and (2) the additional water molecules after forming the hydration layer in polySBMA solutions show higher freedom than those in PEG. These results might illustrate the reason for higher resistance of zwitterionic materials to nonspecific protein adsorptions compared to that of PEGs.     

Faster proton transfer dynamics of water on SnO2 compared to TiO2

Proton jump processes in the hydration layer on the iso-structural TiO(2) rutile (110) and SnO(2) cassiterite (110) surfaces were studied with density functional theory molecular dynamics. We find that the proton jump rate is more than three times faster on cassiterite compared with rutile. A local analysis based on the correlation between the stretching band of the O-H vibrations and the strength of H-bonds indicates that the faster proton jump activity on cassiterite is produced by a stronger H-bond formation between the surface and the hydration layer above the surface. The origin of the increased H-bond strength on cassiterite is a combined effect of stronger covalent bonding and stronger electrostatic interactions due to differences of its electronic structure. The bridging oxygens form the strongest H-bonds between the surface and the hydration layer. This higher proton jump rate is likely to affect reactivity and catalytic activity on the surface. A better understanding of its origins will enable methods to control these rates.     

Hydrogen bond breaking mechanism and water reorientational dynamics in the hydration layer of lysozyme

The mechanism and the rate of hydrogen bond-breaking in the hydration layer surrounding an aqueous protein are important ingredients required to understand the various aspects of protein dynamics, its function, and stability. Here, we use computer simulation and a time correlation function technique to understand these aspects in the hydration layer of lysozyme. Water molecules in the layer are found to exhibit three distinct bond-breaking mechanisms. A large angle orientational jump of the donor water molecule is common among all of them. In the most common ( approximately 80%) bond-breaking event in the layer, the new acceptor water molecule comes from the first coordination shell (initially within 3.5 A of the donor), and the old acceptor water molecule remains within the first coordination shell, even after the bond-breaking. This is in contrast to that in bulk water, in which both of the acceptor molecules involve the second coordination shell. Additionally, the motion of the incoming and the outgoing acceptor molecules involved is not diffusive in the hydration layer, in contrast to their observed diffusive motion in the bulk. The difference in rotational dynamics between the bulk and the hydration layer water molecules is clearly manifested in the calculated time-dependent angular van Hove self-correlation function ( G(theta, t)) which has a pronounced two-peak structure in the layer, and this can be traced to the constrained translational motion in the layer. The longevity of the surrounding hydrogen bond network is found to be significantly enhanced near a hydrophilic residue.     

Multiple relaxation processes versus the fragile-to-strong transition in confined water

Broadband dielectric spectroscopy data on water confined in three different environments, namely at the surface of a globular protein or inside the small pores of two silica substrates, in the temperature range 140 K ≤ T ≤ 300 K, are presented and discussed in comparison with previous results from different techniques. It is found that all samples show a fast relaxation process, independently of the hydration level and confinement size. This relaxation is well known in the literature and its cross-over from Arrhenius to non-Arrhenius temperature behavior is the object of vivid debate, given its claimed relation to the existence of a second critical point of water. We find such a cross-over at a temperature of ~180 K, and assign the relaxation process to the layer of molecules adjacent and strongly interacting with the substrate surface. This is the water layer known to have the highest density and slowest translational dynamics compared to the average: its apparent cross-over may be due to the freezing of some degree of freedom and survival of very localized motions alone, to the onset of finite size effects, or to the presence of a calorimetric glass transition of the hydration shell at ~170 K. Another relaxation process is visible in water confined in the silica matrices: this is slower than the previous one and has distinct temperature behaviors, depending on the size of the confining volume and consequent ice nucleation.     

Spatial hydration structures and dynamics of phenol in sub- and supercritical water

The hydration structures and dynamics of phenol in aqueous solution at infinite dilution are investigated using molecular-dynamics simulation technique. The simulations are performed at several temperatures along the coexistence curve of water up to the critical point, and above the critical point with density fixed at 0.3 g/cm3. The hydration structures of phenol are characterized using the radial, cylindrical, and spatial distribution functions. In particular, full spatial maps of local atomic (solvent) density around a solute molecule are presented. It is demonstrated that in addition to normal H bonds with hydroxyl group of phenol, water forms pi-type complexes with the center of the benzene ring, in which H2O molecules act as H-bond donor. At ambient conditions phenol is solvated by 38 water molecules, which make up a large hydrophobic cavity, and forms on average 2.39 H bonds (1.55 of which are due to the hydroxyl group-water interactions and 0.84 are due to the pi complex) with its hydration shell. As temperature increases, the hydration structure of phenol undergoes significant changes. The disappearance of the pi-type H bonding is observed near the critical point. Self-diffusion coefficients of water and phenol are also calculated. Dramatic increase in the diffusivity of phenol in aqueous solution is observed near the critical point of simple point-charge-extended water and is related to the changes in water structure at these conditions.     

Nature of the water molecules in the palisade layer of a triton X-100 micelle in the presence of added salts: A solvation dynamics study

The effect of added electrolytes on the nature of water molecules in the palisade layer of a Triton X-100 (TX-100) micelle has been investigated using solvation dynamics studies of C153 dye in the presence of different concentrations of NaCl, KCl, and CsCl salts. In all of the cases, the solvation dynamics is found to be biexponential in nature. It is seen that in the presence of added salts the solvation dynamics becomes slower. As previously reported (Charlton et al. J. Phys. Chem. B 2000, 104, 8327; Molina-Bolivar et al. J. Phys. Chem. B 2002, 106, 870), the presence of salt increases micellar hydration (and also size) for TX-100, mainly due to enhancement in the mechanically trapped water content in the palisade layer. Under normal circumstances, increased micellar hydration was expected to cause faster solvation dynamics (Kumbhakar et al. J. Phys. Chem. B 2004, 108, 19246), though in the present work, a reverse trend is in fact observed with the added salts. In accordance with solvation dynamics results, fluorescence anisotropy studies also indicate an increase in microviscosity for the palisade layer of the TX-100 micelle with the added salts. The present results have been rationalized assuming that the ions reside in the palisade layer, and due to the hydration of the ions, especially the cations, the water molecules in the palisade layer undergo a kind of clustering, causing the microviscosity to in fact increase rather than decrease as expected due to increased micellar hydration. A partial collapse of the surfactant chains due to their dehydration as caused by the hydration of the ions in the palisade layer may also add to the increase in microviscosity and the consequent retardation in relaxation dynamics in the presence of salts.     

Influence of pyrogenic silica hydration on water diffusion in surface layers

The surface hydration of pyrogenic silica (aerosil) has been studied by slow neutron scattering. It was shown that the rate of diffusion changes with the thickness of the layer monotonically. The mean square displacement of water molecules from the equilibrium position in aerosil hydration shell is smaller than in bulk water, but twice larger than in ice.     

Preferential hydration of lysozyme in water/glycerol mixtures: a small-angle neutron scattering study

In solution small-angle neutron scattering has been used to study the solvation properties of lysozyme dissolved in water/glycerol mixtures. To detect the characteristics of the protein-solvent interface, 35 different experimental conditions (i.e., protein concentration, water/glycerol fraction in the solvent, content of deuterated compounds) have been considered and a suitable software has been developed to fit simultaneously the whole set of scattering data. The average composition of the solvent in the close vicinity of the protein surface at each experimental condition has been derived. In all the investigated conditions, glycerol resulted especially excluded from the protein surface, confirming that lysozyme is preferentially hydrated. By considering a thermodynamic hydration model based on an equilibrium exchange between water and glycerol from the solvation layer to the bulk, the preferential binding coefficient and the excess solvation number have been estimated. Results were compared with data previously derived for ribonuclease A in the same mixed solvent: even if the investigated solvent compositions were very different, the agreement between data is noticeable, suggesting that a unique mechanism presides over the preferential hydration process. Moreover, the curve describing the excess solvation number as a function of the solvent composition shows the occurrence of a region of maximal hydration, which probably accounts for the changes in protein stability detected in the presence of cosolvents.     

Hydration layer of a cationic micelle, C(10)TAB: structure, rigidity, slow reorientation, hydrogen bond lifetime, and solvation dynamics

We report a theoretical study of the structure and dynamics of the water layer (the hydration layer) present at the surface of the cationic micelle decyltrimethylammonium bromide (DeTAB) by using atomistic molecular dynamics simulations. The simulated micelle consisted of 47 surfactant molecules (and an equal number of bromide ions), in good agreement with the pioneering light scattering experiments by Debye which found an aggregation number of 50. In this micelle, three partially positively charged methyl groups of each surfactant headgroup face the surrounding water. The nature of the cationic micellar surface is found to play an important role in determining the arrangement of water which is quite different from that in the bulk or on the surface of an anionic micelle, like cesium perfluorooctanoate. Water molecules present in the hydration layer are found to be preferentially distributed in the region between the three partially charged methyl headgroups. It is found that both the translational and rotational motions of water exhibit appreciably slower dynamics in the layer than those in the bulk. The solvation time correlation function (TCF) of bromide ions exhibits a long time component which is found to originate primarily from the interaction of the probe with the micellar headgroups. Thus, the decay of the solvation TCF is controlled largely by the residence time of the probe in the surface. The residence time distribution of the water molecules also exhibits a slow time component. We also calculate the collective number density fluctuation in the layer and find a prominent slow component compared to the similar quantity in the bulk. This slow component demonstrates that water structure in the hydration layer is more rigid than that in the bulk. These results demonstrate that the slow dynamics of hydration layer water is generic to macromolecular surfaces of either polarity.