Water dynamics in protein hydration shells: the molecular origins of the dynamical perturbation

Picosecond Solvation Dynamics in Nanoconfinement: Role of Water and Host-Guest Complexation

The dynamics of solvation of an excited chromophore, 5-(4″-dimethylaminophenyl)-2-(4'-sulfophenyl)oxazole, sodium salt (DMO), has been explored in confined nanoscopic environments of β-cyclodextrin (βCD) and heptakis(2,6-di- O-methyl)-β-cyclodextrin (DIMEB). Solvation occurs on a distinctly slower time scale (τ_S³ ∼ 47 ps, τ_S⁴ ∼ 517 ps) in the host cavity of DIMEB than in that of βCD (τ_S³ ∼ 20 ps, τ_S⁴ ∼ 174 ps). The calculated equilibrium solvation response of DMO was characterized by four relaxation components (τ_S¹ ∼ 0.46-0.48 ps, τ_S² ∼ 3.2-3.4 ps, τ_S³ ∼ 32.3-37.7 ps, and τ_S⁴ ∼ 232-485 ps), of which the longer ones (τ_S³, τ_S⁴) are well-consistent with experiments, whereas the ultrafast components (τ_S¹, τ_S²) are unresolved. The observed time constant (τ_S³) within the ∼20-47 ps range arises from slow water molecules in the primary hydration layers of the host CDs and is slower for DIMEB than for βCD presumably due to longer-lived and stronger hydrogen bonds that water forms with the former host relative to the latter. Decomposition of the calculated solvation response of DMO has revealed that conformational fluctuations of the macrocyclic hosts give rise to the observed long-time relaxation component (τ_S⁴), which is much slower for the inclusion complexes with DIMEB than for those with βCD because of slower conformational dynamics of the former host than that of the latter.

Hydrophobic Collapse in N-Methylacetamide-Water Mixtures

Aqueous N-methylacetamide solutions were investigated by polarization-resolved pump-probe and 2D infrared spectroscopy (2D IR), using the amide I mode as a reporter. The 2D IR results are compared with molecular dynamics simulations and spectral calculations to gain insight into the molecular structures in the mixture. N-Methylacetamide and water molecules tend to form clusters with "frozen" amide I dynamics. This is driven by a hydrophobic collapse as the methyl groups of the N-methylacetamide molecules cluster in the presence of water. Since the studied system can be considered as a simplified model for the backbone of proteins, the present study forms a convenient basis for understanding the structural and vibrational dynamics in proteins. It is particularly interesting to find out that a hydrophobic collapse as the one driving protein folding is observed in such a simple system.

HullRad: Fast Calculations of Folded and Disordered Protein and Nucleic Acid Hydrodynamic Properties

Hydrodynamic properties are useful parameters for estimating the size and shape of proteins and nucleic acids in solution. The calculation of such properties from structural models informs on the solution properties of these molecules and complements corresponding structural studies. Here we report, to our knowledge, a new method to accurately predict the hydrodynamic properties of molecular structures. This method uses a convex hull model to estimate the hydrodynamic volume of the molecule and is orders of magnitude faster than common methods. It works well for both folded proteins and ensembles of conformationally heterogeneous proteins and for nucleic acids. Because of its simplicity and speed, the method should be useful for the modification of computer-generated, intrinsically disordered protein ensembles and ensembles of flexible, but folded, molecules in which rapid calculation of experimental parameters is needed. The convex hull method is implemented in a Python script called HullRad. The use of the method is facilitated by a web server and the code is freely available for batch applications.

The effect of Hofmeister anions on water structure at protein surfaces

To understand the effects of specific ions on protein-water interactions and the thermodynamic stability of proteins in salt solutions, we use a molecular dynamics (MD) simulation to examine the water structure, orientational distribution, and dynamics near the surface of ubiquitin. In particular, we consider NaCl, NaBF₄, NaSCN, and NaClO₄ salt solutions containing ubiquitin, where the anions of the latter three salts are well-known chaotropic ions in the Hofmiester anion series. The number of hydrogen bonds (H-bonds) per water molecule is found to decrease significantly at the ubiquitin-water interface, indicating a significant disruption of the water H-bonding network. The distribution of the water H-bond numbers near the protein surface is modulated by dissolved ions, and the extent of the ion effect on the H-bonding network structure follows the order of the Hofmeister anion series, while there are no specific ion effects on water properties at distances larger than 5 Å from the protein surface. From detailed analyses of the surface area, volume, and root-mean-square deviation (RMSD) of ubiquitin, we show that changes in the properties of the protein could originate from the disruption of the water H-bond network induced by ions with a higher affinity for the protein surface instead of direct protein residue-ion interactions. An interesting observation made here is that the orientational distribution of water molecules at the protein-water interface is close to random, but there is a slight preference for interfacial water molecules with a straddle structure within 2.5 Å of the protein surface, where one of the two OH groups points away from the protein surface and the other points toward the surface. In addition, comparing the MD simulation results for ubiquitin solutions with dissolved NaSCN and KSCN, we show that Na⁺ affects the water H-bonding structure at the protein surface more than K⁺. It is clear that the H-bonding network structure of water more than one water layer away from the protein surface is not distinguishably different from that of neat water. We thus anticipate that the present work will provide insights into the scale of specific ion effects on the H-bonding structure and orientational distribution of water in the vicinity of protein surfaces in aqueous solutions.

Effect of Temperature and Hydration Level on Purple Membrane Dynamics Studied Using Broadband Dielectric Spectroscopy from Sub-GHz to THz Regions

To investigate the effects of temperature and hydration on the dynamics of purple membrane (PM), we measured the broadband complex dielectric spectra from 0.5 GHz to 2.3 THz using a vector network analyzer and terahertz time-domain spectroscopy from 233 to 293 K. In the lower temperature region down to 83 K, the complex dielectric spectra in the THz region were also obtained. The complex dielectric spectra were analyzed through curve fitting using several model functions. We found that the hydrated states of one relaxational mode, which was assigned as the coupled motion of water molecules with the PM surface, began to overlap with the THz region at approximately 230 K. On the other hand, the relaxational mode was not observed for the dehydrated state. On the basis of this result, we conclude that the protein-dynamical-transition-like behavior in the THz region is due to the onset of the overlap of the relaxational mode with the THz region. Temperature hysteresis was observed in the dielectric spectrum at 263 K when the hydration level was high. It is suggested that the hydration water behaves similarly to supercooled liquid at that temperature. The third hydration layer may be partly formed to observe such a phenomenon. We also found that the relaxation time is slower than that of a globular protein, lysozyme, and the microscopic environment in the vicinity of the PM surface is suggested to be more heterogeneous than lysozyme. It is proposed that the spectral overlap of the relaxational mode and the low-frequency vibrational mode is necessary for the large conformational change of protein.

Anomalous water dynamics at surfaces and interfaces: synergistic effects of confinement and surface interactions

In nature, water is often found in contact with surfaces that are extended on the scale of molecule size but small on a macroscopic scale. Examples include lipid bilayers and reverse micelles as well as biomolecules like proteins, DNA and zeolites, to name a few. While the presence of surfaces and interfaces interrupts the continuous hydrogen bond network of liquid water, confinement on a mesoscopic scale introduces new features. Even when extended on a molecular scale, natural and biological surfaces often have features (like charge, hydrophobicity) that vary on the scale of the molecular diameter of water. As a result, many new and exotic features, which are not seen in the bulk, appear in the dynamics of water close to the surface. These different behaviors bear the signature of both water-surface interactions and of confinement. In other words, the altered properties are the result of the synergistic effects of surface-water interactions and confinement. Ultrafast spectroscopy, theoretical modeling and computer simulations together form powerful synergistic approaches towards an understanding of the properties of confined water in such systems as nanocavities, reverse micelles (RMs), water inside and outside biomolecules like proteins and DNA, and also between two hydrophobic walls. We shall review the experimental results and place them in the context of theory and simulations. For water confined within RMs, we discuss the possible interference effects propagating from opposite surfaces. Similar interference is found to give rise to an effective attractive force between two hydrophobic surfaces immersed and kept fixed at a separation of d, with the force showing an exponential dependence on this distance. For protein and DNA hydration, we shall examine a multitude of timescales that arise from frustration effects due to the inherent heterogeneity of these surfaces. We pay particular attention to the role of orientational correlations and modification of the same due to interaction with the surfaces.

Insights into the role of electrostatics in temperature adaptation: a comparative study of psychrophilic, mesophilic, and thermophilic subtilisin-like serine proteases

To investigate the role of electrostatics in different temperature adaptations, we performed a comparative study on subtilisin-like serine proteases from psychrophilic Vibrio sp. PA-44 (VPR), mesophilic Engyodontium album (Tritirachium album) (PRK), and thermophilic Thermus aquaticus (AQN) using multiple-replica molecular dynamics (MD) simulations combined with continuum electrostatics calculations. The results reveal that although salt bridges are not a crucial factor in determining the overall thermostability of these three proteases, they on average provide the greatest, moderate, and least electrostatic stabilization to AQN, PRK, and VPR, respectively, at the respective organism growth temperatures. Most salt bridges in AQN are effectively stabilizing and thus contribute to maintaining the overall structural stability at 343 K, while nearly half of the salt bridges in VPR interconvert between being stabilizing and being destabilizing, likely aiding in enhancing the local conformational flexibility at 283 K. The individual salt bridges, salt-bridge networks, and calcium ions contribute differentially to local stability and flexibility of these three enzyme structures, depending on their spatial distributions and electrostatic strengths. The shared negatively charged surface potential at the active center of the three enzymes may provide the active-center flexibility necessary for nucleophilic attack and proton transfer. The differences in distributions of the electro-negative, electro-positive, and electro-neutral potentials, particularly over the back surfaces of the three proteases, may modulate/affect not only protein solubility and thermostability but also structural stability and flexibility/rigidity. These results demonstrate that electrostatics contributes to both heat and cold adaptation of subtilisin-like serine proteases through fine-tuning, either globally or locally, the structural stability and conformational flexibility/rigidity, thus providing a foundation for further engineering and mutagenesis studies.

Differently charged surface patches contribute to temperature adaptation of subtilisin-like serine proteases through affecting/modulating the protein solubility and thermostability and the structural flexibility/rigidity/stability.

Revival of collective water structure and dynamics in reverse micelles brought about by protein encapsulation

A novel mechanism of depolarization in reverse micelles with zwitterionic surfactants and containing polar species but lacking ions is reported.

Macromolecular crowding and the importance of proper hydration for the structure and dynamics of protein solutions

Extensive computational studies of ubiquitin crowding with a special focus on protein hydration directly visible in dielectric spectra.

Spatially Heterogeneous Surface Water Diffusivity around Structured Protein Surfaces at Equilibrium

Hydration water on the surface of a protein is thought to mediate the thermodynamics of protein-ligand interactions. For hydration water to play a role beyond modulating global protein solubility or stability, the thermodynamic properties of hydration water must reflect on the properties of the heterogeneous protein surface and thus spatially vary over the protein surface. A potent read-out of local variations in thermodynamic properties of hydration water is its equilibrium dynamics spanning picosecond to nanosecond time scales. In this study, we employ Overhauser dynamic nuclear polarization (ODNP) to probe the equilibrium hydration water dynamics at select sites on the surface of Chemotaxis Y (CheY) in dilute solution. ODNP reports on site-specific hydration water dynamics within 5-10 Å of a label tethered to the biomolecular surface on two separate time scales of motion, corresponding to diffusive water (DW) and protein-water coupled motions, referred to as bound water (BW). We find DW dynamics to be highly heterogeneous across the surface of CheY. We identify a significant correlation between DW dynamics and the local hydropathy of the CheY protein surface, empirically determined by molecular dynamics (MD) simulations, and find the more hydrophobic sites to be hydrated with slower diffusing water. Furthermore, we compare the hydration water dynamics on different polypeptides and liposome surfaces and find the DW dynamics on globular proteins to be significantly more heterogeneous than on intrinsically disordered proteins (IDPs), peptides, and liposomes. The heterogeneity in the hydration water dynamics suggests that structured proteins have the capacity to encode information into the surrounding hydration shell.

Picosecond orientational dynamics of water in living cells

Cells are extremely crowded, and a central question in biology is how this affects the intracellular water. Here, we use ultrafast vibrational spectroscopy and dielectric-relaxation spectroscopy to observe the random orientational motion of water molecules inside living cells of three prototypical organisms: Escherichia coli, Saccharomyces cerevisiae (yeast), and spores of Bacillus subtilis. In all three organisms, most of the intracellular water exhibits the same random orientational motion as neat water (characteristic time constants ~9 and ~2 ps for the first-order and second-order orientational correlation functions), whereas a smaller fraction exhibits slower orientational dynamics. The fraction of slow intracellular water varies between organisms, ranging from ~20% in E. coli to ~45% in B. subtilis spores. Comparison with the water dynamics observed in solutions mimicking the chemical composition of (parts of) the cytosol shows that the slow water is bound mostly to proteins, and to a lesser extent to other biomolecules and ions.

The cytoplasm’s crowdedness leads one to expect that cell water is different from bulk water. By measuring the rotational motion of water molecules in living cells, Tros et al. find that apart from a small fraction of water solvating biomolecules, cell water has the same dynamics as bulk water.

Biopreservation of Myoglobin in Crowded Environment: A Comparison between Gelatin and Trehalose Matrixes

Biopreservation by sugar and/or polymeric matrixes is a thoroughly studied research topic with wide technological relevance. Ternary amorphous systems containing both saccharides and proteins are extensively exploited to model the in vivo biopreservation process. With the aim of disentangling the effect of saccharides and polypeptidic crowders (such as gelatin) on the preservation of a model protein, we present here a combined differential scanning calorimetry and UV-vis spectrophotometry study on samples of myoglobin embedded in amorphous gelatin and trehalose + gelatin matrixes at different hydrations, and compare them with amorphous myoglobin-only and myoglobin-trehalose samples. The results point out the different effects of gelatin, which acts mainly as a crowding agent, and trehalose, which acts mainly by direct interaction. Gelatin is able to improve effectively the protein thermal stability at very low hydration; however, it has small effects at medium to high hydration. Consistently, gelatin appears to be more effective than trehalose against massive denaturation in the long time range, while the mixed trehalose + collagen matrix is most effective in preserving protein functionality, outdoing both gelatin-only and trehalose-only matrixes.

The dynamics of peptide-water interactions in dialanine: An ultrafast amide I 2D IR and computational spectroscopy study

We present a joint experimental and computational study of the dynamic interactions of dialanine (Ala-Ala) with water, comparing the results of ultrafast 2D IR and infrared transient absorption spectroscopy of its amide I vibration with spectra modeled from molecular dynamics (MD) simulations. The experimental data are analyzed to describe vibrational frequency fluctuations, vibrational energy relaxation, and chemical exchange processes. The origin of these processes in the same underlying fluctuating forces allows a common description in terms of the fluctuations and conformational dynamics of the peptide and associated solvent. By comparing computational spectroscopy from MD simulations with multiple force fields and water models, we describe how the dynamics of water hydrogen bond fluctuations and switching processes act as a source of friction that governs the dephasing and vibrational relaxation, and provide a description of coupled water and peptide motions that give rise to spectroscopic exchange processes.

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.

Computational Studies on Acetylcholinesterases

Functions of biomolecules, in particular enzymes, are usually modulated by structural fluctuations. This is especially the case in a gated diffusion-controlled reaction catalyzed by an enzyme such as acetylcholinesterase. The catalytic triad of acetylcholinesterase is located at the bottom of a long and narrow gorge, but it catalyzes the extremely rapid hydrolysis of the neurotransmitter, acetylcholine, with a reaction rate close to the diffusion-controlled limit. Computational modeling and simulation have produced considerable advances in exploring the dynamical and conformational properties of biomolecules, not only aiding in interpreting the experimental data, but also providing insights into the internal motions of the biomolecule at the atomic level. Given the remarkably high catalytic efficiency and the importance of acetylcholinesterase in drug development, great efforts have been made to understand the dynamics associated with its functions by use of various computational methods. Here, we present a comprehensive overview of recent computational studies on acetylcholinesterase, expanding our views of the enzyme from a microstate of a single structure to conformational ensembles, strengthening our understanding of the integration of structure, dynamics and function associated with the enzyme, and promoting the structure-based and/or mechanism-based design of new inhibitors for it.

Revival of the Intermolecular Nuclear Overhauser Effect for Mapping Local Protein Hydration Dynamics

The highly heterogeneous hydration dynamics of protein-water interfaces is considered important for protein stability and dynamics, protein folding, enzymatic activity, and even drug design. The nuclear Overhauser effect (NOE) between protein and water protons is the only experimental observable which, in principle, can provide a map of locally resolved hydration dynamics. However, its utility was questioned in various theoretical studies that emphasized the contributions of long-range NOE interactions. We show by a detailed analysis based on molecular dynamics simulations that, contrary to recent claims, the protein-water NOE is an excellent observable to map local hydration dynamics at the protein surface.

In Silico Studies of Small Molecule Interactions with Enzymes Reveal Aspects of Catalytic Function

Small molecules, such as solvent, substrate, and cofactor molecules, are key players in enzyme catalysis. Computational methods are powerful tools for exploring the dynamics and thermodynamics of these small molecules as they participate in or contribute to enzymatic processes. In-depth knowledge of how small molecule interactions and dynamics influence protein conformational dynamics and function is critical for progress in the field of enzyme catalysis. Although numerous computational studies have focused on enzyme-substrate complexes to gain insight into catalytic mechanisms, transition states and reaction rates, the dynamics of solvents, substrates, and cofactors are generally less well studied. Also, solvent dynamics within the biomolecular solvation layer play an important part in enzyme catalysis, but a full understanding of its role is hampered by its complexity. Moreover, passive substrate transport has been identified in certain enzymes, and the underlying principles of molecular recognition are an area of active investigation. Enzymes are highly dynamic entities that undergo different conformational changes, which range from side chain rearrangement of a residue to larger-scale conformational dynamics involving domains. These events may happen nearby or far away from the catalytic site, and may occur on different time scales, yet many are related to biological and catalytic function. Computational studies, primarily molecular dynamics (MD) simulations, provide atomistic-level insight and site-specific information on small molecule interactions, and their role in conformational pre-reorganization and dynamics in enzyme catalysis. The review is focused on MD simulation studies of small molecule interactions and dynamics to characterize and comprehend protein dynamics and function in catalyzed reactions. Experimental and theoretical methods available to complement and expand insight from MD simulations are discussed briefly.

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.

Hydrogen Bond Network of Water around Protein Investigated with Terahertz and Infrared Spectroscopy

The dynamical and structural properties of water at protein interfaces were characterized on the basis of the broadband complex dielectric constant (0.25 to 400 THz) of albumin aqueous solutions. Our analysis of the dielectric responses between 0.25 and 12 THz first revealed hydration water with retarded reorientational dynamics extending ∼8.5 Å (corresponding to three to four layers) out from the albumin surface. Second, the number of nonhydrogen-bonded water was decreased in the presence of the albumin solute, indicating protein inhibits the fragmentation of the water hydrogen-bond network. Finally, water molecules at the albumin interface were found to form a distorted hydrogen-bond structure due to topological and energetic disorder of the protein surface. In addition, the intramolecular O-H stretching vibration of water (∼100 THz), which is sensitive to hydrogen-bond environment, pointed to a trend that hydration water has a larger population of strongly hydrogen-bonded water molecules compared with that of bulk water. From these experimental results, we concluded that the "strengthened" water hydrogen bonds at the protein interface dynamically slow down the reorientational motion of water and form the less-defective hydrogen-bond network by inhibiting the fragmentation of water-water hydrogen bonds. Nevertheless, such a strengthened water hydrogen-bond network is composed of heterogeneous hydrogen-bond distances and angles, and thus characterized as structurally "distorted."

Water is an active matrix of life for cell and molecular biology

Szent-Győrgi called water the "matrix of life" and claimed that there was no life without it. This statement is true, as far as we know, on our planet, but it is not clear whether it must hold throughout the cosmos. To evaluate that question requires a close consideration of the many varied and subtle roles that water plays in living cells-a consideration that must be free of both an assumed essentialism that gives water an almost mystical life-giving agency and a traditional tendency to see it as a merely passive solvent. Water is a participant in the "life of the cell," and here I describe some of the features of that active agency. Water's value for molecular biology comes from both the structural and dynamic characteristics of its status as a complex, structured liquid as well as its nature as a polar, protic, and amphoteric reagent. Any discussion of water as life's matrix must, however, begin with an acknowledgment that our understanding of it as both a liquid and a solvent is still incomplete.

Photoswitchable Sn-Cyt c Solid-State Devices

Electron transfer across proteins plays an important role in many biological processes, including those relevant for the conversion of solar photons to chemical energy. Previous studies demonstrated the generation of photocurrents upon light irradiation in a number of photoactive proteins, such as photosystem I or bacteriorhodopsin. Here, it is shown that Sn-cytochrome c layers act as reversible and efficient photoelectrochemical switches upon integration into large-area solid-state junctions. Photocurrents are observed both in the Soret band (λ = 405 nm) and in the Q band (λ = 535 nm), with current on/off ratios reaching values of up to 25. The underlying modulation in charge-transfer rate is attributed to a hole-transport channel created by the photoexcitation of the Sn-porphyrin.

Hydration-Shell Transformation of Thermosensitive Aqueous Polymers

Although water plays a key role in the coil-globule transition of polymers and biomolecules, it is not clear whether a change in water structure drives or follows polymer collapse. Here, we address this question by using Raman multivariate curve resolution (Raman-MCR) spectroscopy to investigate the hydration shell structure around poly(N-isopropylacrylamide) (PNIPAM) and poly(propylene oxide) (PPO), both below and above the cloud point temperature at which the polymers collapse and form mesoscopic polymer-rich aggregates. We find that, upon clouding, the water surrounding long PNIPAM chains transforms to a less ordered and more weakly hydrogen bonded structure, while the water surrounding short PNIPAM and PPO chains remains similar above and below the cloud point. Furthermore, microfluidic temperature jump studies demonstrate that the onset of clouding precedes the hydration-shell structural transformation, and thus the observed water structural transformation is associated with ripening of aggregates composed of long-chain polymers, on a time scale that is long compared to the onset of clouding.

How many water molecules are detected in X-ray protein crystal structures?

The positions of several water molecules can be determined in protein crystallography, either buried in internal cavities or at the protein surface. It is important to be able to estimate the expected number of these water molecules to facilitate crystal structure determination. Here, a multiple Poisson regression model implemented on nearly 10 000 protein crystal structures shows that the number of detectable water molecules depends on eight variables: crystallographic resolution, R factor, percentage of solvent in the crystal, average B factor of the protein atoms, percentage of amino acid residues in loops, average solvent-accessible surface area of the amino acid residues, grand average of hydropathy of the protein(s) in the asymmetric unit and normalized number of heteroatoms that are not water molecules. Furthermore, a secondary analysis tested the effect of different software packages. Given the values of these eight variables, it is possible to compute the expected number of water molecules detectable in electron-density maps with reasonable accuracy (as suggested by an external validation of the model).

Intracellular water – an overlooked drug target? Cisplatin impact in cancer cells probed by neutrons

Intracellular water as a secondary pharmacological target?

Towards a complete characterization of the δ-dispersion in dielectric spectroscopy of protein–water systems

The δ-process in dielectric spectroscopy of protein–water systems is computationally analyzed in great detail, in relation to other experiments.

On the coupling between the dynamics of protein and water

The solvation entropy of flexible protein regions is higher than that of rigid regions and contributes differently to the overall thermodynamic stability.

Local water sensing: water exchange in bacterial photosynthetic reaction centers embedded in a trehalose glass studied using multiresonance EPR

Pulsed EPR spectroscopies and isotope labeled water are applied to detect and quantify the local water in a bacterial reaction center embedded into a trehalose glass.

Description of Hydration Water in Protein (Green Fluorescent Protein) Solution

The structurally and dynamically perturbed hydration shells that surround proteins and biomolecules have a substantial influence upon their function and stability. This makes the extent and degree of water perturbation of practical interest for general biological study and industrial formulation. We present an experimental description of the dynamical perturbation of hydration water around green fluorescent protein in solution. Less than two shells (∼5.5 Å) were perturbed, with dynamics a factor of 2-10 times slower than bulk water, depending on their distance from the protein surface and the probe length of the measurement. This dependence on probe length demonstrates that hydration water undergoes subdiffusive motions (τ ∝ q^-2.5 for the first hydration shell, τ ∝ q^-2.3 for perturbed water in the second shell), an important difference with neat water, which demonstrates diffusive behavior (τ ∝ q^-2). These results help clarify the seemingly conflicting range of values reported for hydration water retardation as a logical consequence of the different length scales probed by the analytical techniques used.

Assessment of Hydration Thermodynamics at Protein Interfaces with Grid Cell Theory

Molecular dynamics simulations have been analyzed with the Grid Cell Theory (GCT) method to spatially resolve the binding enthalpies and entropies of water molecules at the interface of 17 structurally diverse proteins. Correlations between computed energetics and structural descriptors have been sought to facilitate the development of simple models of protein hydration. Little correlation was found between GCT-computed binding enthalpies and continuum electrostatics calculations. A simple count of contacts with functional groups in charged amino acids correlates well with enhanced water stabilization, but the stability of water near hydrophobic and polar residues depends markedly on its coordination environment. The positions of X-ray-resolved water molecules correlate with computed high-density hydration sites, but many unresolved waters are significantly stabilized at the protein surfaces. A defining characteristic of ligand-binding pockets compared to nonbinding pockets was a greater solvent-accessible volume, but average water thermodynamic properties were not distinctive from other interfacial regions. Interfacial water molecules are frequently stabilized by enthalpy and destabilized entropy with respect to bulk, but counter-examples occasionally occur. Overall detailed inspection of the local coordinating environment appears necessary to gauge the thermodynamic stability of water in protein structures.

In silico characterization of protein partial molecular volumes and hydration shells

In this paper we present a computational approach, based on NVT molecular dynamics trajectories, that allows the direct evaluation of the protein partial molecular volume. The results obtained for five different globular proteins demonstrate the accuracy of this computational procedure in reproducing protein partial molecular volumes, providing quantitative characterization of the hydration shell in terms of the protein excluded volume, hydration shell ellipsoidal volume and related solvent density. Remarkably, our data indicate for the hydration shell a ≈10% solvent density increase with respect to the liquid water bulk density, in excellent agreement with the available experimental data.

Hydration Dynamics of a Peripheral Membrane Protein

Water dynamics in the hydration shell of the peripheral membrane protein annexin B12 were studied using MD simulations and Overhauser DNP-enhanced NMR. We show that retardation of water motions near phospholipid bilayers is extended by the presence of a membrane-bound protein, up to around 10 Å above that protein. Near the membrane surface, electrostatic interactions with the lipid head groups strongly slow down water dynamics, whereas protein-induced water retardation is weaker and dominates only at distances beyond 10 Å from the membrane surface. The results can be understood from a simple model based on additive contributions from the membrane and the protein to the activation free energy barriers of water diffusion next to the biomolecular surfaces. Furthermore, analysis of the intermolecular vibrations of the water network reveals that retarded water motions near the membrane shift the vibrational modes to higher frequencies, which we used to identify an entropy gradient from the membrane surface toward the bulk water. Our results have implications for processes that take place at lipid membrane surfaces, including molecular recognition, binding, and protein-protein interactions.

Rotational dynamics of water molecules near biological surfaces with implications for nuclear quadrupole relaxation

Based on Molecular Dynamics simulations of two different systems, the protein ubiquitin dissolved in water and an AOT reverse micelle, we present a broad analysis of the single particle rotational dynamics of water. A comprehensive connection to NQR, which is a prominent experimental method in this field, is developed, based on a reformulation of its theoretical framework. Interpretation of experimental NQR results requires a model which usually assumes that the NQR experiences retardation only in the first hydration shell. Indeed, the present study shows that this first-shell model is correct. Moreover, previous experimental retardation factors are quantitatively reproduced. All of this is seemingly contradicted by results of other methods, e.g., dielectric spectroscopy, responsible for a long-standing debate in this field. Our detailed analysis shows that NQR omits important information contained in overall water dynamics, most notably, the retardation of the water dipole axis in the electric field exerted by a biological surface.

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.

Proteins Take up Water Before Unfolding

Proteins perform specific biological functions that strongly depend on their three-dimensional structure. This three-dimensional structure, i.e. the way the protein folds, is strongly determined by the interaction between the protein and the water solvent. We study the dynamics of water in aqueous solutions of several globular proteins at different degrees of urea-induced unfolding, using polarization-resolved femtosecond infrared spectroscopy. We observe that a fraction of the water molecules is strongly slowed down by their interaction with the protein surface. By monitoring the slow water fraction we can directly probe the amount of water-exposed protein surface. We find that at mild denaturing conditions, the water-exposed surface increases by almost 50%, while the secondary structure is still intact. This finding indicates that protein unfolding starts with the protein structure becoming less tight, thereby allowing water to enter.