Reply to "Comment on 'Hydration and mobility of trehalose in aqueous solution'"

Universality in the Structure and Dynamics of Water under Lipidic Mesophase Soft Nanoconfinement

Water under soft nanoconfinement features physical and chemical properties fundamentally different from bulk water; yet, the multitude and specificity of confining systems and geometries mask any of its potentially universal traits. Here, we advance in this quest by resorting to lipidic mesophases as an ideal nanoconfinement system, allowing inspecting the behavior of water under systematic changes in the topological and geometrical properties of the confining medium, without altering the chemical nature of the interfaces. By combining Terahertz absorption spectroscopy experiments and molecular dynamics simulations, we unveil the presence of universal laws governing the physics of nanoconfined water, recapitulating the data collected at varying levels of hydration and nanoconfinement topologies. This geometry-independent universality is evidenced by the existence of master curves characterizing both the structure and dynamics of simulated water as a function of the distance from the lipid-water interface. Based on our theoretical findings, we predict a parameter-free law describing the amount of interfacial water against the structural dimension of the system (i.e., the lattice parameter), which captures both the experimental and numerical results within the same curve, without any fitting. Our results offer insight into the fundamental physics of water under soft nanoconfinement and provide a practical tool for accurately estimating the amount of nonbulk water based on structural experimental data.

How proteins modify water dynamics

Much of biology happens at the protein-water interface, so all dynamical processes in this region are of fundamental importance. Local structural fluctuations in the hydration layer can be probed by ¹⁷O magnetic relaxation dispersion (MRD), which, at high frequencies, measures the integral of a biaxial rotational time correlation function (TCF)-the integral rotational correlation time. Numerous ¹⁷O MRD studies have demonstrated that this correlation time, when averaged over the first hydration shell, is longer than in bulk water by a factor 3-5. This rotational perturbation factor (RPF) has been corroborated by molecular dynamics simulations, which can also reveal the underlying molecular mechanisms. Here, we address several outstanding problems in this area by analyzing an extensive set of molecular dynamics data, including four globular proteins and three water models. The vexed issue of polarity versus topography as the primary determinant of hydration water dynamics is resolved by establishing a protein-invariant exponential dependence of the RPF on a simple confinement index. We conclude that the previously observed correlation of the RPF with surface polarity is a secondary effect of the correlation between polarity and confinement. Water rotation interpolates between a perturbed but bulk-like collective mechanism at low confinement and an exchange-mediated orientational randomization (EMOR) mechanism at high confinement. The EMOR process, which accounts for about half of the RPF, was not recognized in previous simulation studies, where only the early part of the TCF was examined. Based on the analysis of the experimentally relevant TCF over its full time course, we compare simulated and measured RPFs, finding a 30% discrepancy attributable to force field imperfections. We also compute the full ¹⁷O MRD profile, including the low-frequency dispersion produced by buried water molecules. Computing a local RPF for each hydration shell, we find that the perturbation decays exponentially with a decay "length" of 0.3 shells and that the second and higher shells account for a mere 3% of the total perturbation measured by ¹⁷O MRD. The only long-range effect is a weak water alignment in the electric field produced by an electroneutral protein (not screened by counterions), but this effect is negligibly small for ¹⁷O MRD. By contrast, we find that the ¹⁷O TCF is significantly more sensitive to the important short-range perturbations than the other two TCFs examined here.

The spatial range of protein hydration

Proteins interact with their aqueous surroundings, thereby modifying the physical properties of the solvent. The extent of this perturbation has been investigated by numerous methods in the past half-century, but a consensus has still not emerged regarding the spatial range of the perturbation. To a large extent, the disparate views found in the current literature can be traced to the lack of a rigorous definition of the perturbation range. Stating that a particular solvent property differs from its bulk value at a certain distance from the protein is not particularly helpful since such findings depend on the sensitivity and precision of the technique used to probe the system. What is needed is a well-defined decay length, an intrinsic property of the protein in a dilute aqueous solution, that specifies the length scale on which a given physical property approaches its bulk-water value. Based on molecular dynamics simulations of four small globular proteins, we present such an analysis of the structural and dynamic properties of the hydrogen-bonded solvent network. The results demonstrate unequivocally that the solvent perturbation is short-ranged, with all investigated properties having exponential decay lengths of less than one hydration shell. The short range of the perturbation is a consequence of the high energy density of bulk water, rendering this solvent highly resistant to structural perturbations. The electric field from the protein, which under certain conditions can be long-ranged, induces a weak alignment of water dipoles, which, however, is merely the linear dielectric response of bulk water and, therefore, should not be thought of as a structural perturbation. By decomposing the first hydration shell into polarity-based subsets, we find that the hydration structure of the nonpolar parts of the protein surface is similar to that of small nonpolar solutes. For all four examined proteins, the mean number of water-water hydrogen bonds in the nonpolar subset is within 1% of the value in bulk water, suggesting that the fragmentation and topography of the nonpolar protein-water interface has evolved to minimize the propensity for protein aggregation by reducing the unfavorable free energy of hydrophobic hydration.

On the dielectric decrement of electrolyte solutions: a dressed-ion theory analysis

Based on the dressed-ion theory and a simple physical argument regarding the conductivity of the solution, we derive a relation between the ionic strength and dielectric constant of an electrolyte solution. At its simplest, this model gives the dielectric constant at low ionic strength I as ε_r(I) = ε_r(0)(1 + αI)^-1, where α (the excess polarization) is directly related to the dressed-ion charge. One contribution to the origin of the dielectric decrement is thus seen to stem from the electrostatic screening of the ions in solution, with no solvent contributions necessary.

X-ray diffraction measurement of cosolvent accessible volume in rhombohedral insulin crystals

X-ray crystallographic measurement of the number of solvent electrons in the unit cell of a protein crystal equilibrated with aqueous solutions of different densities provides information about preferential hydration in the crystalline state. Room temperature and cryo-cooled rhombohedral insulin crystals were equilibrated with 1.2M trehalose to study the effect of lowered water activity. The native and trehalose soaked crystals were isomorphous and had similar structures. Including all the low resolution data, the amplitudes of the structure factors were put on an absolute scale (in units of electrons per asymmetric unit) by constraining the integrated number of electrons inside the envelope of the calculated protein density map to equal the number deduced from the atomic model. This procedure defines the value of F(000), the amplitude at the origin of the Fourier transform, which is equal to the total number of electrons in the asymmetric unit (i.e. protein plus solvent). Comparison of the F(000) values for three isomorphous pairs of room temperature insulin crystals, three with trehalose and three without trehalose, indicates that 75±12 electrons per asymmetric unit were added to the crystal solvent when soaked in 1.2M trehalose. If all the water in the crystal were available as solvent for the trehalose, 304 electrons would have been added. Thus, the co-solvent accessible volume is one quarter of the total water in the crystal. Determination of the total number of electrons in a protein crystal is an essential first step for mapping the average density distribution of the disordered solvent.

Water Dynamics in the Hydration Shells of Biomolecules

The structure and function of biomolecules are strongly influenced by their hydration shells. Structural fluctuations and molecular excitations of hydrating water molecules cover a broad range in space and time, from individual water molecules to larger pools and from femtosecond to microsecond time scales. Recent progress in theory and molecular dynamics simulations as well as in ultrafast vibrational spectroscopy has led to new and detailed insight into fluctuations of water structure, elementary water motions, electric fields at hydrated biointerfaces, and processes of vibrational relaxation and energy dissipation. Here, we review recent advances in both theory and experiment, focusing on hydrated DNA, proteins, and phospholipids, and compare dynamics in the hydration shells to bulk water.

Perspective: Structure and ultrafast dynamics of biomolecular hydration shells

The structure and function of biomolecules can be strongly influenced by their hydration shells. A key challenge is thus to determine the extent to which these shells differ from bulk water, since the structural fluctuations and molecular excitations of hydrating water molecules within these shells can cover a broad range in both space and time. Recent progress in theory, molecular dynamics simulations, and ultrafast vibrational spectroscopy has led to new and detailed insight into the fluctuations of water structure, elementary water motions, and electric fields at hydrated biointerfaces. Here, we discuss some central aspects of these advances, focusing on elementary molecular mechanisms and processes of hydration on a femto- to picosecond time scale, with some special attention given to several issues subject to debate.

Molecular properties of aqueous solutions: a focus on the collective dynamics of hydration water

When a solute is dissolved in water, their mutual interactions determine the molecular properties of the solute on one hand, and the structure and dynamics of the surrounding water particles (the so-called hydration water) on the other. The very existence of soft matter and its peculiar properties are largely due to the wide variety of possible water-solute interactions. In this context, water is not an inert medium but rather an active component, and hydration water plays a crucial role in determining the structure, stability, dynamics, and function of matter. This review focuses on the collective dynamics of hydration water in terms of retardation with respect to the bulk, and of the number of molecules whose dynamics is perturbed. Since water environments are in a dynamic equilibrium, with molecules continuously exchanging from around the solute towards the bulk and vice versa, we examine the ability of different techniques to measure the water dynamics on the basis of the explored time scales and exchange rates. Special emphasis is given to the collective dynamics probed by extended depolarized light scattering and we discuss whether and to what extent the results obtained in aqueous solutions of small molecules can be extrapolated to the case of large biomacromolecules. In fact, recent experiments performed on solutions of increasing complexity clearly indicate that a reductionist approach is not adequate to describe their collective dynamics. We conclude this review by presenting current ideas that are being developed to describe the dynamics of water interacting with macromolecules.

Hydration of proteins and nucleic acids: Advances in experiment and theory. A review

BACKGROUND

Most biological processes involve water, and the interactions of biomolecules with water affect their structure, function and dynamics.

SCOPE OF REVIEW

This review summarizes the current knowledge of protein and nucleic acid interactions with water, with a special focus on the biomolecular hydration layer. Recent developments in both experimental and computational methods that can be applied to the study of hydration structure and dynamics are reviewed, including software tools for the prediction and characterization of hydration layer properties.

MAJOR CONCLUSIONS

In the last decade, important advances have been made in our understanding of the factors that determine how biomolecules and their aqueous environment influence each other. Both experimental and computational methods contributed to the gradually emerging consensus picture of biomolecular hydration.

GENERAL SIGNIFICANCE

An improved knowledge of the structural and thermodynamic properties of the hydration layer will enable a detailed understanding of the various biological processes in which it is involved, with implications for a wide range of applications, including protein-structure prediction and structure-based drug design.

Water-mediated influence of a crowded environment on internal vibrations of a protein molecule

The influence of crowding on the protein inner dynamics is examined by putting a single protein molecule close to one or two neighboring protein molecules. The presence of additional molecules influences the amplitudes of protein fluctuations. Also, a weak dynamical coupling of collective velocities of surface atoms of proteins separated by a layer of water is detected. The possible mechanisms of these phenomena are described. The cross-correlation function of the collective velocities of surface atoms of two proteins was decomposed into the Fourier series. The amplitude spectrum displays a peak at low frequencies. Also, the results of principal component analysis suggest that the close presence of an additional protein molecule influences the high-amplitude, low-frequency modes in the most prominent way. This part of the spectrum covers biologically important protein motions. The neighbor-induced changes in the inner dynamics of the protein may be connected with the changes in the velocity power spectrum of interfacial water. The additional protein molecule changes the properties of solvation water and in this way it can influence the dynamics of the second protein. It is suggested that this phenomenon may be described, at first approximation, by a damped oscillator driven by an external random force. This model was successfully applied to conformationally rigid Choristoneura fumiferana antifreeze protein molecules.

Molecular dynamics simulations and NMR spectroscopy studies of trehalose-lipid bilayer systems

The disaccharide trehalose (TRH) strongly affects the physical properties of lipid bilayers. We investigate interactions between lipid membranes formed by 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and TRH using NMR spectroscopy and molecular dynamics (MD) computer simulations. We compare dipolar couplings derived from DMPC/TRH trajectories with those determined (i) experimentally in TRH using conventional high-resolution NMR in a weakly ordered solvent (bicelles), and (ii) by solid-state NMR in multilamellar vesicles (MLV) formed by DMPC. Analysis of the experimental and MD-derived couplings in DMPC indicated that the force field used in the simulations reasonably well describes the experimental results with the exception for the glycerol fragment that exhibits significant deviations. The signs of dipolar couplings, not available from the experiments on highly ordered systems, were determined from the trajectory analysis. The crucial step in the analysis of residual dipolar couplings (RDCs) in TRH determined in a bicelle-environment was access to the conformational distributions derived from the MD trajectory. Furthermore, the conformational behavior of TRH, investigated by J-couplings, in the ordered and isotropic phases is essentially identical, indicating that the general assumptions in the analyses of RDCs are well founded.

Comparative study of hydration shell dynamics around a hyperactive antifreeze protein and around ubiquitin

The hydration layer surrounding a protein plays an essential role in its biochemical function and consists of a heterogeneous ensemble of water molecules with different local environments and different dynamics. What determines the degree of dynamical heterogeneity within the hydration shell and how this changes with temperature remains unclear. Here, we combine molecular dynamics simulations and analytic modeling to study the hydration shell structure and dynamics of a typical globular protein, ubiquitin, and of the spruce budworm hyperactive antifreeze protein over the 230-300 K temperature range. Our results show that the average perturbation induced by both proteins on the reorientation dynamics of water remains moderate and changes weakly with temperature. The dynamical heterogeneity arises mostly from the distribution of protein surface topographies and is little affected by temperature. The ice-binding face of the antifreeze protein induces a short-ranged enhancement of water structure and a greater slowdown of water reorientation dynamics than the non-ice-binding faces whose effect is similar to that of ubiquitin. However, the hydration shell of the ice-binding face remains less tetrahedral than the bulk and is not "ice-like". We finally show that the hydrogen bonds between water and the ice-binding threonine residues are particularly strong due to a steric confinement effect, thereby contributing to the strong binding of the antifreeze protein on ice crystals.