Atomistic details of protein dynamics and the role of hydration water

Temporal Resolution of Activity-Related Solvation Dynamics in the TIM Barrel Enzyme Murine Adenosine Deaminase

Murine adenosine deaminase (mADA) is a prototypic system for studying the thermal activation of active site chemistry within the TIM barrel family of enzyme reactions. Previous temperature-dependent hydrogen deuterium exchange studies under various conditions have identified interconnected thermal networks for heat transfer from opposing protein-solvent interfaces to active site residues in mADA. One of these interfaces contains a solvent exposed helix-loop-helix moiety that presents the hydrophobic face of its long α-helix to the backside of bound substrate. Herein we pursue the time and temperature dependence of solvation dynamics at the surface of mADA, for comparison to established kinetic parameters that represent active site chemistry. We first created a modified protein devoid of native tryptophans with close to native kinetic behavior. Single site-specific tryptophan mutants were back inserted into each of the four positions where native tryptophans reside. Measurements of nanosecond fluorescence relaxation lifetimes and Stokes shift decays, that reflect time dependent environmental reoroganization around the photo-excited state of Trp*, display minimal temperature dependences. These regions serve as controls for the behavior of a new single tryptophan inserted into a solvent exposed region near the helix-loop-helix moiety located behind the bound substrate, Lys54Trp. This installed Trp displays a significantly elevated value for Ea ( k Stokes shift ) ; further, when Phe61 within the long helix positioned behind bound substrate is replaced by a series of aliphatic hydrophobic side chains, the trends in Ea ( k Stokes shift ) mirror the earlier reported impact of the same series of function-altering hydrophobic side chains on the activation energy of catalysis, Ea ( k cat ) .The reported experimental findings implicate a solvent initiated and rapid (>ns) protein restructuring that contributes to the enthalpic activation barrier to catalysis in mADA.

Water molecule ordering on the surface of an intrinsically disordered protein

The dynamics and local structure of the hydration water on surfaces of folded proteins have been extensively investigated. However, our knowledge of the hydration of intrinsically disordered proteins (IDPs) is more limited. Here, we compare the local structure of water molecules hydrating a globular protein, lysozyme, and the intrinsically disordered N-terminal of c-Src kinase (SH4UD) using molecular dynamics simulation. The radial distributions from the protein surface of the first and the second hydration shells are similar for the folded protein and the IDP. However, water molecules in the first hydration shell of both the folded protein and the IDP are perturbed from the bulk. This perturbation involves a loss of tetrahedrality, which is, however, significantly more marked for the folded protein than the IDP. This difference arises from an increase in the first hydration shell of the IDP of the fraction of hydration water molecules interacting with oxygen. The water ordering is independent of the compactness of the IDP. In contrast, the lifetimes of water molecules in the first hydration shell increase with IDP compactness, indicating a significant impact of IDP configuration on water surface pocket kinetics, which here is linked to differential pocket volumes and polarities.

Nonthermal acceleration of protein hydration by sub-terahertz irradiation

The collective intermolecular dynamics of protein and water molecules, which overlap in the sub-terahertz (THz) frequency region, are relevant for expressing protein functions but remain largely unknown. This study used dielectric relaxation (DR) measurements to investigate how externally applied sub-THz electromagnetic fields perturb the rapid collective dynamics and influence the considerably slower chemical processes in protein-water systems. We analyzed an aqueous lysozyme solution, whose hydration is not thermally equilibrated. By detecting time-lapse differences in microwave DR, we demonstrated that sub-THz irradiation gradually decreases the dielectric permittivity of the lysozyme solution by reducing the orientational polarization of water molecules. Comprehensive analysis combining THz and nuclear magnetic resonance spectroscopies suggested that the gradual decrease in the dielectric permittivity is not induced by heating but is due to a slow shift toward the hydrophobic hydration structure in lysozyme. Our findings can be used to investigate hydration-mediated protein functions based on sub-THz irradiation.

Triplet States Reveal Slow Local Dynamics in the Solvation Shell of Biomolecules

Protein hydration shell dynamics plays a pivotal role in biochemical processes such as protein folding, enzyme function, molecular recognition and interaction with biological membranes. Thus, it is crucial to understand the mobility of the solvation shell at the surface of biomolecules. Triplet state solvation dynamics can reveal the slow dynamics of the solvation shell. This is done in the present work without adding separate dye molecules but instead by using a phosphorescent subgroup of the biomolecule itself. In particular, we study a small heptapeptide in a glycerol-water mixture under cryoconservation conditions so that the system can be supercooled without crystallization. We find a significant slowing of molecules in the solvation shell in the millisecond range compared to the bulk. This opens up the possibility to unravel the nature of relaxation processes in the solvation shell usually overlapping at room temperature.

Hydration Dynamics of Model Peptides with Different Hydrophobic Character

The multi-scale dynamics of aqueous solutions of the hydrophilic peptide N-acetyl-glycine-methylamide (NAGMA) have been investigated through extended frequency-range depolarized light scattering (EDLS), which enables the broad-band detection of collective polarizability anisotropy fluctuations. The results have been compared to those obtained for N-acetyl-leucinemethylamide (NALMA), an amphiphilic peptide which shares with NAGMA the same polar backbone, but also contains an apolar group. Our study indicates that the two model peptides induce similar effects on the fast translational dynamics of surrounding water. Both systems slow down the mobility of solvating water molecules by a factor 6–8, with respect to the bulk. Moreover, the two peptides cause a comparable far-reaching spatial perturbation extending to more than two hydration layers in diluted conditions. The observed concentration dependence of the hydration number is explained considering the random superposition of different hydration shells, while no indication of solute aggregation phenomena has been found. The results indicate that the effect on the dynamics of water solvating the amphiphilic peptide is dominated by the hydrophilic backbone. The minor impact of the hydrophobic moiety on hydration features is consistent with structural findings derived by Fourier transform infrared (FTIR) measurements, performed in attenuated total reflectance (ATR) configuration. Additionally, we give evidence that, for both systems, the relaxation mode in the GHz frequency range probed by EDLS is related to solute rotational dynamics. The rotation of NALMA occurs at higher timescales, with respect to the rotation of NAGMA; both processes are significantly slower than the structural dynamics of hydration water, suggesting that solute and solvent motions are uncoupled. Finally, our results do not indicate the presence of super-slow water (relaxation times in the order of tens of picoseconds) around the peptides investigated.

Recognition of N6-Methyladenosine by the YTHDC1 YTH Domain Studied by Molecular Dynamics and NMR Spectroscopy: The Role of Hydration

The YTH domain of YTHDC1 belongs to a class of protein "readers", recognizing the N6-methyladenosine (m⁶A) chemical modification in mRNA. Static ensemble-averaged structures revealed details of N6-methyl recognition via a conserved aromatic cage. Here, we performed molecular dynamics (MD) simulations along with nuclear magnetic resonance (NMR) and isothermal titration calorimetry (ITC) to examine how dynamics and solvent interactions contribute to the m⁶A recognition and negative selectivity toward an unmethylated substrate. The structured water molecules surrounding the bound RNA and the methylated substrate's ability to exclude bulk water molecules contribute to the YTH domain's preference for m⁶A. Intrusions of bulk water deep into the binding pocket disrupt binding of unmethylated adenosine. The YTHDC1's preference for the 5'-Gm⁶A-3' motif is partially facilitated by a network of water-mediated interactions between the 2-amino group of the guanosine and residues in the m⁶A binding pocket. The 5'-Im⁶A-3' (where I is inosine) motif can be recognized too, but disruption of the water network lowers affinity. The D479A mutant also disrupts the water network and destabilizes m⁶A binding. Our interdisciplinary study of the YTHDC1 protein-RNA complex reveals an unusual physical mechanism by which solvent interactions contribute toward m⁶A recognition.

Nonthermal excitation effects mediated by sub-terahertz radiation on hydrogen exchange in ubiquitin

Water dynamics in the hydration layers of biomolecules play crucial roles in a wide range of biological functions. A hydrated protein contains multiple components of diffusional and vibrational dynamics of water and protein, which may be coupled at ∼0.1-THz frequency (10-ps timescale) at room temperature. However, the microscopic description of biomolecular functions based on various modes of protein-water-coupled motions remains elusive. A novel approach for perturbing the hydration dynamics in the subterahertz frequency range and probing them at the atomic level is therefore warranted. In this study, we investigated the effect of klystron-based, intense 0.1-THz excitation on the slow dynamics of ubiquitin using NMR-based measurements of hydrogen-deuterium exchange. We demonstrated that the subterahertz irradiation accelerated the hydrogen-deuterium exchange of the amides located in the interior of the protein and hydrophobic surfaces while decelerating this exchange in the amides located in the surface loop and short 3₁₀ helix regions. This subterahertz-radiation-induced effect was qualitatively contradictory to the increased-temperature-induced effect. Our results suggest that the heterogeneous water dynamics occurring at the protein-water interface include components that are nonthermally excited by the subterahertz radiation. Such subterahertz-excited components may be linked to the slow function-related dynamics of the protein.

Similarities and Differences among Protein Dynamics Studied by Variable Temperature Nuclear Magnetic Resonance Relaxation

Understanding and describing the dynamics of proteins is one of the major challenges in biology. Here, we use multifield variable-temperature NMR longitudinal relaxation (R₁) measurements to determine the hierarchical activation energies of motions of four different proteins: two small globular proteins (GB1 and the SH3 domain of α-spectrin), an intrinsically disordered protein (the C-terminus of the nucleoprotein of the Sendai virus, Sendai Ntail), and an outer membrane protein (OmpG). The activation energies map the motions occurring in the side chains, in the backbone, and in the hydration shells of the proteins. We were able to identify similarities and differences in the average motions of the proteins. We find that the NMR relaxation properties of the four proteins do share similar features. The data characterizing average backbone motions are found to be very similar, the same for methyl group rotations, and similar activation energies are measured. The main observed difference occurs for the intrinsically disordered Sendai Ntail, where we observe much lower energy of activation for motions of protons associated with the protein-solvent interface as compared to the others. We also observe variability between the proteins regarding side chain ¹⁵N relaxation of lysine residues, with a higher activation energy observed in OmpG. This hints at strong interactions with negatively charged lipids in the bilayer and provides a possible mechanistic clue for the "positive-inside" rule for helical membrane proteins. Overall, these observations refine the understanding of the similarities and differences between hierarchical dynamics in proteins.

Tuning Local Hydration Enables a Deeper Understanding of Protein-Ligand Binding: The PP1-Src Kinase Case

Water plays a key role in biomolecular recognition and binding. Despite the development of several computational and experimental approaches, it is still challenging to comprehensively characterize water-mediated effects on the binding process. Here, we investigate how water affects the binding of Src kinase to one of its inhibitors, PP1. Src kinase is a target for treating several diseases, including cancer. We use biased molecular dynamics simulations, where the hydration of predetermined regions is tuned at will. This computational technique efficiently accelerates the SRC-PP1 binding simulation and allows us to identify several key and yet unexplored aspects of the solvent's role. This study provides a further perspective on the binding phenomenon, which may advance the current drug design approaches for the development of new kinase inhibitors.

Quasielastic neutron scattering studies on couplings of protein and water dynamics in hydrated elastin

Performing quasielastic neutron scattering measurements and analyzing both elastic and quasielasic contributions, we study protein and water dynamics of hydrated elastin. At low temperatures, hydration-independent methyl group rotation dominates the findings. It is characterized by a Gaussian distribution of activation energies centered at about E_m = 0.17 eV. At ∼195 K, coupled protein-water motion sets in. The hydration water shows diffusive motion, which is described by a Gaussian distribution of activation energies with E_m = 0.57 eV. This Arrhenius behavior of water diffusion is consistent with previous results for water reorientation, but at variance with a fragile-to-strong crossover at ∼225 K. The hydration-related elastin backbone motion is localized and can be attributed to the cage rattling motion. We speculate that its onset at ∼195 K is related to a secondary glass transition, which occurs when a β relaxation of the protein has a correlation time of τ_β ∼ 100 s. Moreover, we show that its temperature-dependent amplitude has a crossover at the regular glass transition T_g = 320 K of hydrated elastin, where the α relaxation of the protein obeys τ_α ∼ 100 s. By contrast, we do not observe a protein dynamical transition when water dynamics enters the experimental time window at ∼240 K.

On the microscopic origins of relaxation processes in aqueous peptide solutions undergoing a glass transition

We combine broadband dielectric spectroscopy (BDS) with ¹H and ²H nuclear magnetic resonance (NMR) to study molecular dynamics in mixtures of ε-polylysine with H₂O or D₂O. In BDS, four relaxation processes can be attributed to molecular dynamics. While the fastest process P1 obeys the Arrhenius law, the slowest process P4 shows prominent non-Arrhenius behavior typical of structural α relaxation. For the intermediate processes P2 and P3, the temperature dependence changes at the glass transition temperature T_g. The ¹H and ²H NMR results yield insights into the molecular origins of these relaxation phenomena. In these NMR analyses, we exploit, in addition to the isotope selectivity of the method, the possibility to distinguish between various types of motion based on their respective line-shape effects and the capability to single out specific molecular moieties based on different spin-lattice relaxation behaviors. In this way, we reveal that process P1 results from the rotation of side and end groups of the peptide, while process P2 is caused by a reorientation of essentially all water molecules, which are quasi-isotropic and survive well below T_g. As for the peptide backbone dynamics, we find evidence that rotational motion of polar groups is involved in process P3 and that nonpolar regions show a dynamical process, which is located between P3 and P4. Thus, the NMR analyses do not yield evidence for coexisting fast peptide-decoupled and slow peptide-coupled water species, which contribute to BDS processes P2 and P3, respectively, but minor bimodality of water motion may remain undetected. Finally, it is demonstrated that the proton/deuteron exchange needs to be considered when interpreting experimental results for molecular dynamics in aqueous peptide solutions.

Length-scale dependence of protein hydration-shell density

Here we present a computational approach based on molecular dynamics (MD) simulation to study the dependence of the protein hydration-shell density on the size of the protein molecule. The hydration-shell density of eighteen different proteins, differing in size, shape and function (eight of them are antifreeze proteins), is calculated. The results obtained show that an increase in the hydration-shell density, relative to that of the bulk, is observed (in the range of 4-14%) for all studied proteins and that this increment strongly correlates with the protein size. In particular, a decrease in the density increment is observed for decreasing protein size. A simple model is proposed in which the basic idea is to approximate the protein molecule as an effective ellipsoid and to partition the relevant parameters, i.e. the solvent-accessible volume and the corresponding solvent density, into two regions: inside and outside the effective protein ellipsoid. It is found that, within the model developed here, almost all of the hydration-density increase is located inside the protein ellipsoid, basically corresponding to pockets within, or at the surface of the protein molecule. The observed decrease in the density increment is caused by the protein size only and no difference is found between antifreeze and non-antifreeze proteins.

Water in protein hydration and ligand recognition

This review describes selected basics of water in biomolecular recognition. We focus on a qualitative understanding of the most important physical aspects, how these change in magnitude between bulk water and protein environment, and how the roles that water plays for proteins arise from them. These roles include mechanical support, thermal coupling, dielectric screening, mass and charge transport, and the competition with a ligand for the occupation of a binding site. The presence or absence of water has ramifications that range from the thermodynamic binding signature of a single ligand up to cellular survival. The large inhomogeneity in water density, polarity and mobility around a solute is hard to assess in experiment. This is a source of many difficulties in the solvation of protein models and computational studies that attempt to elucidate or predict ligand recognition. The influence of water in a protein binding site on the experimental enthalpic and entropic signature of ligand binding is still a point of much debate. The strong water‐water interaction in enthalpic terms is counteracted by a water molecule's high mobility in entropic terms. The complete arrest of a water molecule's mobility sets a limit on the entropic contribution of a water displacement process, while the solvent environment sets limits on ligand reactivity.

Surface Topography Effects of Globular Biomolecules on Hydration Water

Hydration water serves as a microscopic manifestation of structural stability and functions of biomolecules. To develop bio-nanomaterials in applications, it is important to study how the surface topography and heterogeneity of biomolecules result in their diversity of the hydration dynamics and energetics. We here performed molecular dynamics simulations combined with the steered molecular dynamics and umbrella sampling to investigate the dynamics and escape process associated with the free energy change of water molecules close to a globular biomolecule, i.e., hemoglobin (Hb) and G-quadruplex DNA (GDNA). The residence time, power of long-time tail, and dipole relaxation time were found to display drastic changes within the averaged hydration shell of 3.0-5.0 Å. Compared with bulk water, in the inner hydration shell, the water dipole moment displays a slower relaxation process and is more oriented toward GDNA than toward Hb, forming a hedgehog-like structure when it surrounds GDNA. In particular, a spine water structure is observed in the GDNA narrow groove. The water isotope effect not only prolongs the dynamic time scales of libration motion in the inner hydration shell and the dipole relaxation processes in the bulk but also strengthens the DNA spine water structure. The potential of the mean force profile reflects the integrity of the hydration shell structure and enables us to obtain detailed insights into the structures formed by water, such as the caged H-bond network and the edge bridge structures; it also reveals that local hydration shell free energy (LHSFE) depends on H-bond rupture processes and ranges from 0.2 to 4.2 kcal/mol. Our results demonstrate that the surface topography of a biomolecule influences the integrity of the hydration shell structure and LHSFE. Our studies are able to identify various further applications in the areas of microfluid devices and nano-dewetting on bioinspired surfaces.

Observation of high-temperature macromolecular confinement in lyophilised protein formulations using terahertz spectroscopy

•

Structural dynamics in lyophilised protein formulations can be probed with terahertz spectroscopy and two glass transition processes, Tg,α and Tg,β, are observed.

•

Vibrational confinement upon thermal activation is observed resulting in no detectable changes in secondary structure but strongly reduced the molecular mobility at temperatures above Tg,α.

•

The confinement was found to be strongly dependent on the formulation.

•

We hypothesise that confinement is linked to conformational states with potential effects on physical and chemical stability of the biomolecule during storage.

Characterising the structural dynamics of proteins and the effects of excipients are critical for optimising the design of formulations. In this work we investigated four lyophilised formulations containing bovine serum albumin (BSA) and three formulations containing a monoclonal antibody (mAb, here mAb1), and explored the role of the excipients polysorbate 80, sucrose, trehalose, and arginine on stabilising proteins. By performing temperature variable terahertz time-domain spectroscopy (THz-TDS) experiments it is possible to study the vibrational dynamics of these formulations. The THz-TDS measurements reveal two distinct glass transition processes in all tested formulations. The lower temperature transition, Tg,β, is associated with the onset of local motion due to the secondary relaxation whilst the higher temperature transition, Tg,α, marks the onset of the α-relaxation. For some of the formulations, containing globular BSA as well as mAb1, the absorption at terahertz frequencies does not increase further at temperatures above Tg,α. Such behaviour is in contrast to our previous observations for small organic molecules as well as linear polymers where absorption is always observed to steadily increase with temperature due to the stronger absorption of terahertz radiation by more mobile dipoles. The absence of such further increase in absorption with higher temperatures therefore suggests a localised confinement of the protein/excipient matrix at high temperatures that hinders any further increase in mobility. We found that subtle changes in excipient composition had an effect on the transition temperatures Tg,α and Tg,β as well as the vibrational confinement in the solid state. Further work is required to establish the potential significance of the vibrational confinement in the solid state on formulation stability and chemical degradation as well as what role the excipients play in achieving such confinement.

Dynamics of proteins in solution

The dynamics of proteins in solution includes a variety of processes, such as backbone and side-chain fluctuations, interdomain motions, as well as global rotational and translational (i.e. center of mass) diffusion. Since protein dynamics is related to protein function and essential transport processes, a detailed mechanistic understanding and monitoring of protein dynamics in solution is highly desirable. The hierarchical character of protein dynamics requires experimental tools addressing a broad range of time- and length scales. We discuss how different techniques contribute to a comprehensive picture of protein dynamics, and focus in particular on results from neutron spectroscopy. We outline the underlying principles and review available instrumentation as well as related analysis frameworks.

Protein Solvent Shell Structure Provides Rapid Analysis of Hydration Dynamics

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.

An isolate of Potato Virus X capsid protein from N. benthamiana: Insights from homology modeling and molecular dynamics simulation

Since Potato Virus X (PVX) is easily transmitted mechanically between their hosts, its control is difficult. We have previously reported new isolate of this virus (PVX-Iran, GenBank Accession number FJ461343). However, the molecular basis of resistance breaking activity and its relation to capsid protein structure are still not well-understood. SDS-PAGE, ELISA, Western blot and RT-PCR molecular examinations were performed on the inoculated plants Nicotiana benthamiana. The pathological symptoms were related to the PVX isolate. The capsid protein (CP) structure were modeled based on homology and subjected to three independent 80 ns molecular dynamics minimization (GROMACS, OPLS force field) in the SPC water box. The RMSD, RMSF, SASA, and electrostatic properties were retrieved from the trajectories. Flexibility and hydrophilic nature of the N-terminal residues (1-34) of solvated CP could be observed in conformational changes upon minimization. The obtained structure was then docked with NbPCIP1 using ClusPro 2.0. The strong binding affinity of these two proteins (≈-16.0 Kcal mol-1) represents the formation of inclusion body and hence appearance of the symptoms.

Slow Dynamics around a Protein and Its Coupling to Solvent

Solvent is essential for protein dynamics and function, but its role in regulating the dynamics remains debated. Here, we employ saturation transfer electron spin resonance (ST-ESR) to explore the issue and characterize the dynamics on a longer (from μs to s) time scale than has been extensively studied. We first demonstrate the reliability of ST-ESR by showing that the dynamical changeovers revealed in the spectra agree to liquid-liquid transition (LLT) in the state diagram of the glycerol/water system. Then, we utilize ST-ESR with four different probes to systematically map out the variation in local (site-specific) dynamics around a protein surface at subfreezing temperatures (180-240 K) in 10 mol % glycerol/water mixtures. At highly exposed sites, protein and solvent dynamics are coupled, whereas they deviate from each other when temperature is greater than LLT temperature (∼190 K) of the solvent. At less exposed sites, protein however exhibits a dynamic, which is distinct from the bulk solvent, throughout the temperature range studied. Dominant dynamic components are thus revealed, showing that (from low to high temperatures) the overall structural fluctuation, rotamer dynamics, and internal side-chain dynamics, in turn, dominate the temperature dependence of spin-label motions. The structural fluctuation component is relatively slow, collective, and independent of protein structural segments, which is thus inferred to a fundamental dynamic component intrinsic to protein. This study corroborates that bulk solvent plasticizes protein and facilitates rather than slaves protein dynamics.

Effect of glycerol on photobleaching of cytochrome Raman lines in frozen yeast cells

We applied a Raman spectroscopy approach to investigate the effect of a cryoprotectant on the redox state of cytochromes on freezing yeast cells. The redox activity of cytochromes was studied using time-resolved photobleaching of the resonance Raman lines. It is found that ice formation causes a drastic change in the redox state of cytochromes in cells frozen without cryoprotectant, whereas in the presence of glycerol the effects of ice formation are more gradual. The photobleaching rate of cells frozen in glycerol solution shows a gradual slowing with temperature decrease and an abrupt slowdown below - 48 °C. This abrupt decrease was interpreted as originating from changes in protein conformational dynamics. Our findings provide important new insights into the transition from active to inactive cytochrome states as cells undergo freezing in the presence and absence of cryoprotectant.

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.

Complex Dynamics of Water in Protein Confinement

This paper studies single-molecule and collective dynamics of water confined in protein powders by means of molecular dynamics simulations. The single-particle dynamics show a modest retardation compared to the bulk but become highly stretched in the powder, with the stretching exponent of ≃0.2. The collective dynamics of the total water dipole are affected by intermolecular correlations inside water and by cross-correlations between the water and the protein. The dielectric spectrum of water in the powder has two nearly equal-amplitude peaks: a Debye peak with ≃16 ps relaxation time and a highly stretched peak with the relaxation time of ≃13 ns and a stretching exponent of ≃0.12. The slower relaxation component is not seen in the single-molecule correlation functions and can be assigned to elastic protein motions displacing water in the powder. The loss spectrum of the intermediate scattering function reported by neutron-scattering experiments is also highly stretched, with the high-frequency wing scaling according to a power law. Translational dynamics can become much slower in the powder than in the bulk but are overshadowed by the rotational loss in the overall loss spectrum of neutron scattering.

Ergodicity breaking of iron displacement in heme proteins

We present a model of the dynamical transition of atomic displacements in proteins. Increased mean-square displacement at higher temperatures is caused by the softening of the force constant for atomic/molecular displacements by electrostatic and van der Waals forces from the protein-water thermal bath. Displacement softening passes through a nonergodic dynamical transition when the relaxation time of the force-force correlation function enters, with increasing temperature, the instrumental observation window. Two crossover temperatures are identified. The lower crossover, presently connected to the glass transition, is related to the dynamical unfreezing of rotations of water molecules within nanodomains polarized by charged surface residues of the protein. The higher crossover temperature, usually assigned to the dynamical transition, marks the onset of water translations. All crossovers are ergodicity breaking transitions depending on the corresponding observation windows. Allowing stretched exponential relaxation of the protein-water thermal bath significantly improves the theory-experiment agreement when applied to solid protein samples studied by Mössbauer spectroscopy.

WatAA: Atlas of Protein Hydration. Exploring synergies between data mining and ab initio calculations

Water molecules represent an integral part of proteins and a key determinant of protein structure, dynamics and function. WatAA is a newly developed, web-based atlas of amino-acid hydration in proteins. The atlas provides information about the ordered first hydration shell of the most populated amino-acid conformers in proteins. The data presented in the atlas are drawn from two sources: experimental data and ab initio quantum-mechanics calculations. The experimental part is based on a data-mining study of a large set of high-resolution protein crystal structures. The crystal-derived data include 3D maps of water distribution around amino-acids and probability of occurrence of each of the identified hydration sites. The quantum mechanics calculations validate and extend this primary description by optimizing the water position for each hydration site, by providing hydrogen atom positions and by quantifying the interaction energy that stabilizes the water molecule at the particular hydration site position. The calculations show that the majority of experimentally derived hydration sites are positioned near local energy minima for water, and the calculated interaction energies help to assess the preference of water for the individual hydration sites. We propose that the atlas can be used to validate water placement in electron density maps in crystallographic refinement, to locate water molecules mediating protein-ligand interactions in drug design, and to prepare and evaluate molecular dynamics simulations. WatAA: Atlas of Protein Hydration is freely available without login at .

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.