Distinguishing dynamical features of water inside protein hydration layer: Distribution reveals what is hidden behind the average

Structural and aggregation changes of silver carp myosin induced with alcohols: Effects of ethanol, 1,2-propanediol, and glycerol

This study investigated the effects of ethanol, 1,2-propanediol, and glycerol on the structure and aggregation behavior of silver carp (Hypophthalmichthys molitrix) myosin. All alcohols induced extensive alteration in the tertiary structure of myosin. Both ethanol and 1,2-propanediol further promoted an increase in the content of β-sheets in myosin and induced myosin aggregation. While glycerol had almost no impact on the secondary structure of myosin. Molecular dynamics simulations revealed that increasing the concentration of ethanol and 1,2-propanediol affected the overall structural changes in the myosin heavy chain (MHC), while glycerol exerted a more pronounced effect on the MHC tail when compared to the MHC head. Disruption of the hydration layers induced by ethanol and 1,2-propanediol contributed to local structural changes in myosin. Glycerol at a concentration of 20% induced the formation of a larger hydration layer around the MHC tail, which facilitated the stabilization of the protein structure.

Probing the structure of water in individual living cells

Water regulates or even governs a wide range of biological processes. Despite its fundamental importance, surprisingly little is known about the structure of intracellular water. Herein we employ a Raman micro-spectroscopy technique to uncover the composition, abundance and vibrational spectra of intracellular water in individual living cells. In three different cell types, we show a small but consistent population (~3%) of non-bulk-like water. It exhibits a weakened hydrogen-bonded network and a more disordered tetrahedral structure. We attribute this population to biointerfacial water located in the vicinity of biomolecules. Moreover, our whole-cell modeling suggests that all soluble (globular) proteins inside cells are surrounded by, on average, one full molecular layer (about 2.6 Angstrom) of biointerfacial water. Furthermore, relative invariance of biointerfacial water is observed among different single cells. Overall, our study not only opens up experimental possibilities of interrogating water structure in vivo but also provides insights into water in life.

Glassy dynamics in a liquid of anisotropic molecules: Bifurcation of relaxation spectrum

In experimental and theoretical studies of glass transition phenomena, one often finds a sharp crossover in dynamical properties at a temperature Tcr. A bifurcation of a relaxation spectrum is also observed at a temperature TB≈Tcr; both lie significantly above the glass transition temperature. In order to better understand these phenomena, we introduce a new model of glass-forming liquids, a binary mixture of prolate and oblate ellipsoids. This model system exhibits sharp thermodynamic and dynamic anomalies, such as the specific heat jump during heating and a sharp variation in the thermal expansion coefficient around a temperature identified as the glass transition temperature, Tg. The same temperature is obtained from the fit of the calculated relaxation times to the Vogel-Fulcher-Tammann (VFT) form. As the temperature is lowered, the calculated single peak rotational relaxation spectrum splits into two peaks at TB above the estimated Tg. Similar bifurcation is also observed in the distribution of short-to-intermediate time translational diffusion. Interrogation of the two peaks reveals a lower extent of dynamic heterogeneity in the population of the faster mode. We observe an unexpected appearance of a sharp peak in the product of rotational relaxation time τ2 and diffusion constant D at a temperature Tcr, close to TB, but above the glass transition temperature. Additionally, we coarse-grain the system into cubic boxes, each containing, on average, ∼62 particles, to study the average dynamical properties. Clear evidence of large-scale sudden changes in the diffusion coefficient and rotational correlation time signals first-order transitions between low and high-mobility domains.

Significance of the Disulfide Bridge in the Structure and Stability of Metalloprotein Azurin

Metalloproteins make up a class of proteins that incorporate metal ions into their structures, enabling them to perform essential functions in biological systems, such as catalysis and electron transport. Azurin is one such metalloprotein with copper cofactor, having a β-barrel structure with exceptional thermal stability. The copper metal ion is coordinated at one end of the β-barrel structure, and there is a disulfide bond at the opposite end. In this study, we explore the effect of this disulfide bond in the high thermal stability of azurin by analyzing both the native S-S bonded and S-S nonbonded (S-S open) forms using temperature replica exchange molecular dynamics (REMD). Similar to experimental observations, we find a 35 K decrease in denaturation temperature for S-S open azurin compared to that of the native holo form (420 K). As observed in the case of native holo azurin, the unfolding process of the S-S open form also started with disruptions of the α-helix. The free energy surfaces of the unfolding process revealed that the denaturation event of the S-S open form progresses through different sets of conformational ensembles. Subsequently, we compared the stabilities of individual β-sheet strands of both the S-S bonded and the S-S nonbonded forms of azurin. Further, we examined the contacts between individual residues for the central structures from the free energy surfaces of the S-S nonbonded form. The microscopic origin of the lowering in the denaturation temperature is further supplemented by thermodynamic analysis.

Understanding the intricacy of protein in hydrated deep eutectic solvent: Solvation dynamics, conformational fluctuation dynamics, and stability

Deep eutectic solvents (DESs) are potential biocatalytic media due to their easy preparation, fine-tuneability, biocompatibility, and most importantly, due to their ability to keep protein stable and active. However, there are many unanswered questions and gaps in our knowledge about how proteins behave in these alternate media. Herein, we investigated solvation dynamics, conformational fluctuation dynamics, and stability of human serum albumin (HSA) in 0.5 Acetamide/0.3 Urea/0.2 Sorbitol (0.5Ac/0.3Ur/0.2Sor) DES of varying concentrations to understand the intricacy of protein behaviour in DES. Our result revealed a gradual decrease in the side-chain flexibility and thermal stability of HSA beyond 30 % DES. On the other hand, the associated water dynamics around domain-I of HSA decelerate only marginally with increasing DES content, although viscosity rises considerably. We propose that even though macroscopic solvent properties are altered, a protein feels only an aqueous type of environment in the presence of DES. This is probably the first experimental study to delineate the role of the associated water structure of the enzyme for maintaining its stability inside DES. Although considerable effort is necessary to generalize such claims, it might serve as the basis for understanding why proteins remain stable and active in DES.

Osmolyte induced protein stabilization: modulation of associated water dynamics might be a key factor

The mechanism of protein stabilization by osmolytes remains one of the most important and long-standing puzzles. The traditional explanation of osmolyte-induced stability through the preferential exclusion of osmolytes from the protein surface has been seriously challenged by the observations like the concentration-dependent reversal of osmolyte-induced stabilization/destabilization. The more modern explanation of protein stabilization/destabilization by osmolytes considers an indirect effect due to osmolyte-induced distortion of the water structure. It provides a general mechanism, but there are numerous examples of protein-specific effects, i.e., a particular osmolyte might stabilize one protein, but destabilize the other, that could not be rationalized through such an explanation. Herein, we hypothesized that osmolyte-induced modulation of associated water might be a critical factor in controlling protein stability in such a medium. Taking different osmolytes and papain as a protein, we proved that our proposal could explain protein stability in osmolyte media. Stabilizing osmolytes rigidify associated water structures around the protein, whereas destabilizing osmolytes make them flexible. The strong correlation between the stability and the associated water dynamics, and the fact that such dynamics are very much protein specific, established the importance of considering the modulation of associated water structures in explaining the osmolyte-induced stabilization/destabilization of proteins. More interestingly, we took another protein, bromelain, for which a traditionally stabilizing osmolyte, sucrose, acts as a stabilizer at higher concentrations but as a destabilizer at lower concentrations. Our proposal successfully explains such observations, which is probably impossible by any known mechanisms. We believe this report will trigger much research in this area.

Dynamic Electric Field Alignment Determines the Water Rotational Motion around Protein

Water rotational dynamics in biomolecular solution is crucial to evaluating and controlling biomolecule stability. In this molecular dynamics simulation (MD) study on lysozyme solutions, we present how the exerted internal electric field determines water rotational dynamics. We find that the relaxation time of water rotation is equivalent to that of the reorientation of the exerted overall electric field for every single water molecule, regardless of its translation mode. Namely, water molecular rotation synchronizes with the exerted field reorientation. We also map the reorientation process of the electric field at fixed points relative to protein in the solution, which displays the local hydration dynamics commensurate with the reported time-dependent fluorescence Stokes shift (TDFSS) measurements. Comparing the spatial distribution of local field reorientation relaxation time with that of rotational relaxation time, we further suggest that water rotation dynamics are subject to the reorientation of the local overall field within the hydration layer. While outside the hydration layer, the relaxation time of the local electric field reorientation is short enough (subpicosecond) to assume the δ function, showing the electric force with randomly changing orientation is applied to each water molecule.

Stochastic Analysis of Molecular Dynamics Reveals the Rotation Dynamics Distribution of Water around Lysozyme

Water dynamics is essential to biochemical processes by mediating all such reactions, including biomolecular degeneration in solutions. To disentangle the molecular-scale distribution of water dynamics around a solute biomolecule, we investigated here the rotational dynamics of water around lysozyme by combining molecular dynamics (MD) simulations and broadband dielectric spectroscopy (BDS). A statistical analysis using the relaxation times and trajectories of every single water molecule was proposed, and the two-dimensional probability distribution of water at a distance from the lysozyme surface with a rotational relaxation time was given. For the observed lysozyme solutions of 34-284 mg/mL, we discovered that the dielectric relaxation time obtained from this distribution agrees well with the measured γ relaxation time, which suggests that rotational self-correlation of water molecules underlies the gigahertz domain of the dielectric spectra. Regardless of protein concentration, water rotational relaxation time versus the distance from the lysozyme surface revealed that the water rotation is severely retarded within 3 Å from the lysozyme surface and is nearly comparable to pure water when farther than 10 Å. The dimension of the first hydration layer was subsequently identified in terms of the relationship between the acceleration of water rotation and the distance from the protein surface.

Dielectric Constant of Aqueous Solutions of Proteins and Organic Polymers from Molecular Dynamics Simulations

The dielectric constant of water/oligomer mixtures, spanning the range from pure water to pure oligomeric melts, is investigated using molecular dynamics (MD) simulations. As prototypical water-soluble organic substances we consider neutral poly-glycine, poly-ethylene glycol and charged monomeric propanic acid. As the water content is reduced, the dielectric constant decreases but does not follow an ideal mixing behavior. The deviations from ideal mixing originate primarily in the non-linear relation between the oligomer mass fraction and collective polarization effects. We find that the dielectric constant is dominated by water polarization, even if the oligomer mass fraction exceeds 50%. By a double extrapolation of the MD simulation results to the limit of vanishing water fraction and to the limit of infinite oligomeric chain length, we estimate the orientational contribution to the dielectric constant of the pure polymeric melts. By this procedure, we obtain ε = 17 {plus minus} 2 for polyglycine and ε = 1 {plus minus} 0.3 for polyethylene glycol. The large difference is rationalized by polarization correlations of glycine units. Interestingly, we find constant temperature simulations to outperform replica exchange simulations in terms of equilibration speed.

Spatially Resolved Hydration Thermodynamics in Biomolecular Systems

Water is essential for the structure, dynamics, energetics, and thus the function of biomolecules. It is a formidable challenge to elicit, in microscopic detail, the role of the solvation-related driving forces of biomolecular processes, such as the enthalpy and entropy contributions to the underlying free-energy landscape. In this Perspective, we discuss recent developments and applications of computational methods that provide a spatially resolved map of hydration thermodynamics in biomolecular systems and thus yield atomic-level insights to guide the interpretation of experimental observations. An emphasis is on the challenge of quantifying the hydration entropy, which requires characterization of both the motions of the biomolecules and of the water molecules in their surrounding.

Site-Specific Water Dynamics in the First Hydration Layer of an Anti-Freeze Glyco-Protein: A Simulation Study

Antifreeze glycoproteins (AFGPs) inhibit ice recrystallization by a mechanism remaining largely elusive. Dynamics of AFGPs’ hydration water and its involvement in the antifreeze activity, for instance, have not been identified...

Analytical solution of a linear nonlocal Poisson-Boltzmann equation with multiple charges in a spherical solute region surrounded by a water spherical shell

In this paper, an analytical solution of a linear nonlocal Poisson-Boltzmann equation (NPBE) test model with multiple charges in a spherical solute region surrounded by a water spherical shell is derived as a single series of Legendre polynomials and modified spherical Bessel functions. The classic Kirkwood ball model is then shown to be a special case of the NPBE test model so that its analytical solution is regained from a double series of associated Legendre polynomials (derived by Kirkwood in 1934) to a new single series of Legendre polynomials, sharply reducing its computational cost. As an application of these series solutions, a comparison study is done to demonstrate the differences between the Kirkwood and NPBE test models.

Structural Stability of Insulin Oligomers and Protein Association-Dissociation Processes: Free Energy Landscape and Universal Role of Water

Association and dissociation of proteins are important biochemical events. In this Feature Article, we analyze the available studies of these processes for insulin oligomers in aqueous solution. We focus on the solvation of the insulin monomer in water, stability and dissociation of its dimer, and structural integrity of the hexamer. The intricate role of water in solvation of the dimer- and hexamer-forming surfaces, in long-range interactions between the monomers and the stability of the oligomers, is discussed. Ten water molecules inside the central cavity stabilize the structure of the insulin hexamer. We discuss how different order parameters can be used to understand the dissociation of the insulin dimer. The calculation of the rate using a recently computed multidimensional free energy provides considerable insight into the interplay between protein and water dynamics.

Hydration of Simple Model Peptides in Aqueous Osmolyte Solutions

The biology and chemistry of proteins and peptides are inextricably linked with water as the solvent. The reason for the high stability of some proteins or uncontrolled aggregation of others may be hidden in the properties of their hydration water. In this study, we investigated the effect of stabilizing osmolyte–TMAO (trimethylamine N-oxide) and destabilizing osmolyte–urea on hydration shells of two short peptides, NAGMA (N-acetyl-glycine-methylamide) and diglycine, by means of FTIR spectroscopy and molecular dynamics simulations. We isolated the spectroscopic share of water molecules that are simultaneously under the influence of peptide and osmolyte and determined the structural and energetic properties of these water molecules. Our experimental and computational results revealed that the changes in the structure of water around peptides, caused by the presence of stabilizing or destabilizing osmolyte, are significantly different for both NAGMA and diglycine. The main factor determining the influence of osmolytes on peptides is the structural-energetic similarity of their hydration spheres. We showed that the chosen peptides can serve as models for various fragments of the protein surface: NAGMA for the protein backbone and diglycine for the protein surface with polar side chains.

Structural and dynamical heterogeneity of water trapped inside Na+-pumping KR2 rhodopsin in the dark state

Photoisomerization in the retinal leads to a channel opening in rhodopsins that triggers translocation or pumping of ions/protons. Crystal structures of rhodopsins contain several structurally conserved water molecules. It has been suggested that water plays an active role in facilitating the ion pumping/translocation process by acting as a lubricant in these systems. In this paper, we systematically investigate the localization, structure, dynamics, and energetics of the water molecules along the channel for the resting/dark state of KR2 rhodopsin. By employing several microseconds long atomistic molecular dynamics simulation of this trans-membrane protein system, we demonstrate the presence of five distinct water containing pockets/cavities separated by gateways controlled by protein side-chains. There exists a strong hydrogen bonded network involving these buried water molecules and functionally important key residues. We present evidence of significant structural and dynamical heterogeneity in the water molecules present in these cavities, with very rare exchange between them. The exchange time scale of such buried water with the bulk has an extremely wide range, from tens of nanoseconds to >1.5 µs. The translational and rotational dynamics of buried water are found to be strongly dependent on the protein cavity size and local interactions with a classic signature of trapped diffusion and rotational anisotropy.

Slowdown of Water Dynamics from the Top to the Bottom of the GroEL Cavity

The GroE molecular chaperone system is a critical protein machine that assists the folding of substrate proteins in its cavity. Water in the cavity is suspected to play a role in substrate protein folding, but the mechanism is currently unknown. Herein, we report measurements of water dynamics in the equatorial and apical domains of the GroEL cavity in the apo and football states, using site-specific tryptophanyl mutagenesis as an intrinsic optical probe with femtosecond resolution combined with molecular dynamics simulations. We observed clearly different water dynamics in the two domains with a slowdown of the cavity water from the apical to equatorial region in the football state. The results suggest that the GroEL cavity provides a unique water environment that may facilitate substrate protein folding.

Molecular Mechanism of Acceleration and Retardation of Collective Orientation Relaxation of Water Molecules in Aqueous Solutions

The collective orientation relaxation (COR) of water molecules in aqueous solutions is faster or slower with an increase in the concentration of the solutions than that in pure water; for example, acceleration (deceleration) of the COR is observed in a solution of sodium chloride (tetramethylammonium chloride) with increasing concentration. However, the molecular mechanism of the solution and concentration dependence of the relaxation time of the COR has not yet been clarified. We theoretically investigate the concentration dependence of the COR of water molecules in solutions of tetramethylammonium chloride (TMACl), guanidinium chloride (GdmCl), and sodium chloride (NaCl). Based on the Mori-Zwanzig equation, we identify two opposing factors that determine the COR of water molecules in any aqueous solution: the correlation of dipole moments and the single-molecule orientation relaxation. We reveal the molecular mechanism of the concentration dependence of the relaxation time of the COR in the TMACl, GdmCl, and NaCl solutions in terms of these two factors.

Hydrogen-Bond Dynamics and Energetics of Biological Water

Water molecules in the immediate vicinity of biomacromolecules and biomimetic organized assemblies often exhibit a markedly distinct behavior from that of their bulk counterparts. The overall sluggish behavior of biological water substantially affects the stability and integrity of biomolecules, as well as the successful execution of various crucial water-mediated biochemical phenomena. In this Minireview, insights are provided into the features of truncated hydrogen-bond networks that grant biological water its unique characteristics. In particular, experimental results and theoretical investigations, based on chemical kinetics, are presented that have shed light on the dynamics and energetics governing such characteristics. It is emphasized how such details help us to understand the energetics of biological water, an aspect relatively less explored than its dynamics. For instance, when biological water at hydrophilic or charged protein surfaces was explored, the free energy of H-bond breakage was found to be of the order of 0.4 kcal mol^-1 higher than that of bulk water.

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.

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.

Mechanism of Solvent Control of Protein Dynamics

We find that the coupled interactions between protein and water polarization fluctuations play a dominant role in driving the configuration space random walk of solvated proteins. We perform atomistic molecular dynamics simulations on five proteins. Owing to a very low dielectric constant of protein, its dipolar groups experience forces from water along with local forces due to protein atoms. Energy fluctuations reveal a pronounced anticorrelation between protein and water contributions. The protein energy spectrum shows bimodal 1/f noise, which can be attributed to the influence of water on the dynamics of protein.

Computational discovery of chemically patterned surfaces that effect unique hydration water dynamics

The interactions of water with solid surfaces govern their apparent hydrophobicity/hydrophilicity, influenced at the molecular scale by surface coverage of chemical groups of varied nonpolar/polar character. Recently, it has become clear that the precise patterning of surface groups, and not simply average surface coverage, has a significant impact on the structure and thermodynamics of hydration layer water, and, in turn, on macroscopic interfacial properties. Here we show that patterning also controls the dynamics of hydration water, a behavior frequently thought to be leveraged by biomolecules to influence functional dynamics, but yet to be generalized. To uncover the role of surface heterogeneities, we couple a genetic algorithm to iterative molecular dynamics simulations to design the patterning of surface functional groups, at fixed coverage, to either minimize or maximize proximal water diffusivity. Optimized surface configurations reveal that clustering of hydrophilic groups increases hydration water mobility, while dispersing them decreases it, but only if hydrophilic moieties interact with water through directional, hydrogen-bonding interactions. Remarkably, we find that, across different surfaces, coverages, and patterns, both the chemical potential for inserting a methane-sized hydrophobe near the interface and, in particular, the hydration water orientational entropy serve as strong predictors for hydration water diffusivity, suggesting that these simple thermodynamic quantities encode the way surfaces control water dynamics. These results suggest a deep and intriguing connection between hydration water thermodynamics and dynamics, demonstrating that subnanometer chemical surface patterning is an important design parameter for engineering solid-water interfaces with applications spanning separations to catalysis.

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.

Origin of diverse time scales in the protein hydration layer solvation dynamics: A simulation study

In order to inquire the microscopic origin of observed multiple time scales in solvation dynamics, we carry out several computer experiments. We perform atomistic molecular dynamics simulations on three protein-water systems, namely, lysozyme, myoglobin, and sweet protein monellin. In these experiments, we mutate the charges of the neighbouring amino acid side chains of certain natural probes (tryptophan) and also freeze the side chain motions. In order to distinguish between different contributions, we decompose the total solvation energy response in terms of various components present in the system. This allows us to capture the interplay among different self- and cross-energy correlation terms. Freezing the protein motions removes the slowest component that results from side chain fluctuations, but a part of slowness remains. This leads to the conclusion that the slow component approximately in the 20-80 ps range arises from slow water molecules present in the hydration layer. While the more than 100 ps component has multiple origins, namely, adjacent charges in amino acid side chains, hydrogen bonded water molecules and a dynamically coupled motion between side chain and water. In addition, the charges enforce a structural ordering of nearby water molecules and helps to form a local long-lived hydrogen bonded network. Further separation of the spatial and temporal responses in solvation dynamics reveals different roles of hydration and bulk water. We find that the hydration layer water molecules are largely responsible for the slow component, whereas the initial ultrafast decay arises predominantly (approximately 80%) due to the bulk. This agrees with earlier theoretical observations. We also attempt to rationalise our results with the help of a molecular hydrodynamic theory that was developed using classical time dependent density functional theory in a semi-quantitative manner.

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.

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.