Coincidence of dynamical transitions in a soluble protein and its hydration water: direct measurements by neutron scattering and MD simulations

Computational study on the impact of linkage sequence on the structure and dynamics of lignin

Lignin, one of the most abundant biopolymers on Earth, is of great research interest due to its industrial applications including biofuel production and materials science. The structural composition of lignin plays an important role in shaping its properties and functionalities. Notably, lignin exhibits substantial compositional diversity, which varies not only between different plant species but even within the same plant. Currently, it is unclear to what extent this compositional diversity plays on the overall structure and dynamics of lignin. To address this question, this paper reports on the development of two models of lignin containing all guaiacyl (G) subunits with varied linkage sequences and makes use of all-atom molecular dynamics simulations to examine the impact of linkage sequence alone on the lignin's structure and dynamics. This work demonstrates that the structure of the lignin polymer depends on its linkage sequence at temperatures above and below the glass transition temperature ( T g ), but the polymers exhibit similar structural properties as it is approaching the viscous flow state (480 K). At low temperatures, both of lignin models have a local dynamics confined in a cage, but the size of cages varies depending on structural differences. Interestingly, at temperatures higher than T g , the different linkage sequence leads to the subtle dynamical difference which diminishes at 480 K.

Decoupling of the onset of anharmonicity between a protein and its surface water around 200 K

The protein dynamical transition at ~200 K, where the biomolecule transforms from a harmonic, non-functional form to an anharmonic, functional state, has been thought to be slaved to the thermal activation of dynamics in its surface hydration water. Here, by selectively probing the dynamics of protein and hydration water using elastic neutron scattering and isotopic labeling, we found that the onset of anharmonicity in the two components around 200 K is decoupled. The one in protein is an intrinsic transition, whose characteristic temperature is independent of the instrumental resolution time, but varies with the biomolecular structure and the amount of hydration, while the one of water is merely a resolution effect.

Structure and slow dynamics of protein hydration water with cryopreserving DMSO and trehalose upon cooling

We study, through molecular dynamics simulations, three aqueous solutions with one lysozyme protein and three different concentrations of trehalose and dimethyl sulfoxide (DMSO). We analyze the structural and dynamical properties of the protein hydration water upon cooling. We find that trehalose plays a major role in modifying the structure of the network of HBs between water molecules in the hydration layer of the protein. The dynamics of hydration water presents, in addition to the α-relaxation, typical of glass formers, a slower long-time relaxation process, which greatly slows down the dynamics of water, particularly in the systems with trehalose, where it becomes dominant at low temperatures. In all the solutions, we observe, from the behavior of the α-relaxation times, a shift of the Mode Coupling Theory crossover temperature and the fragile-to-strong crossover temperature toward higher values with respect to bulk water. We also observe a strong-to-strong crossover from the temperature behavior of the long-relaxation times. In the aqueous solution with only DMSO, the transition shifts to a lower temperature than in the case with only lysozyme reported in the literature. We observe that the addition of trehalose to the mixture has the opposite effect of restoring the original location of the strong-to-strong crossover. In all the solutions analyzed in this work, the observed temperature of the protein dynamical transition is slightly shifted at lower temperatures than that of the strong-to-strong crossover, but their relative order is the same, showing a correlation between the motion of the protein and that of the hydration water.

Water slowing down drives the occurrence of the low temperature dynamical transition in microgels

The protein dynamical transition marks an increase in atomic mobility and the onset of anharmonic motions at a critical temperature (T d), which is considered relevant for protein functionality. This phenomenon is ubiquitous, regardless of protein composition, structure and biological function and typically occurs at large protein content, to avoid water crystallization. Recently, a dynamical transition has also been reported in non-biological macromolecules, such as poly(N-isopropyl acrylamide) (PNIPAM) microgels, bearing many similarities to proteins. While the generality of this phenomenon is well-established, the role of water in the transition remains a subject of debate. In this study, we use atomistic molecular dynamics (MD) simulations and elastic incoherent neutron scattering (EINS) experiments with selective deuteration to investigate the microscopic origin of the dynamical transition and distinguish water and PNIPAM roles. While a standard analysis of EINS experiments would suggest that the dynamical transition occurs in PNIPAM and water at a similar temperature, simulations reveal a different perspective, also qualitatively supported by experiments. From room temperature down to about 180 K, PNIPAM exhibits only modest changes of dynamics, while water, being mainly hydration water under the probed extreme confinement, significantly slows down and undergoes a mode-coupling transition from diffusive to activated. Our findings therefore challenge the traditional view of the dynamical transition, demonstrating that it occurs in proximity of the water mode-coupling transition, shedding light on the intricate interplay between polymer and water dynamics.

An experimentally representative in-silico protocol for dynamical studies of lyophilised and weakly hydrated amorphous proteins

Characterization of biopolymers in both dry and weakly hydrated amorphous states has implications for the pharmaceutical industry since it provides understanding of the effect of lyophilisation on stability and biological activity. Atomistic Molecular Dynamics (MD) simulations probe structural and dynamical features related to system functionality. However, while simulations in homogenous aqueous environments are routine, dehydrated model assemblies are a challenge with systems investigated in-silico needing careful consideration; simulated systems potentially differing markedly despite seemingly negligible changes in procedure. Here we propose an in-silico protocol to model proteins in lyophilised and weakly hydrated amorphous states that is both more experimentally representative and routinely applicable. Since the outputs from MD align directly with those accessed by neutron scattering, the efficacy of the simulation protocol proposed is shown by validating against experimental neutron data for apoferritin and insulin. This work also highlights that without cooperative experimental and simulative data, development of simulative procedures using MD alone would prove most challenging.

Inelastic neutron scattering and spectroscopy methods to characterize dynamics in colloidal and soft matter systems

Colloidal systems and soft materials are well suited to neutron scattering, and the community has readily adopted elastic scattering techniques to investigate their structure. Due to their unique properties, neutrons may also be used to characterize the dynamics of soft materials over a wide range of length and time scales in situ. Both static structures and an understanding of how molecules move about their equilibrium positions is essential if we are to deliver on the promise of rationally designing soft materials. In this review we introduce the basics of neutron spectroscopy and explore the ways in which inelastic neutron scattering can be used to study colloidal and soft materials. Illustrative examples are chosen that highlight the phenomena suitable for investigation using this suite of techniques.

Direct Observation of the Mutual Coupling Effect in the Protein-Water-Glycerol Mixture by Combining Neutron Scattering and Selective Deuteration

Numerous studies have discussed the impact of cosolvents on the structure, dynamics, and stability of proteins in aqueous solutions. However, the dynamics of cosolvents in the protein-water-cosolvent ternary system is largely unexplored in experiments due to technical difficulty. Consequently, a comprehensive understanding of the interplay among proteins, water, and cosolvents is still lacking. Here, we employed selective deuteration and neutron scattering techniques to characterize the individual motions of each component in the protein/water/glycerol (GLY) mixture across various temperatures. The consistent dynamic onset temperatures and the correlation between the MSD of the protein and the viscosity of solvents revealed the mutual coupling effects among the three components. Furthermore, our experimental and simulation results showed that the hydrogen bond relaxation energy barrier in the ternary system is ∼43 kJ/mol, whereas in the protein-water binary system it is merely ∼35 kJ/mol. Therefore, we suggest that GLY can enhance hydrogen bond interactions in the ternary system through the mutual coupling effect, thereby serving as one of the protective mechanisms of protein preservation by GLY.

Two distinct mechanisms of flavoprotein spectral tuning revealed by low-temperature and time-dependent spectroscopy

Flavins such as flavin mononucleotide or flavin adenine dinucleotide are bound by diverse proteins, yet have very similar spectra when in the oxidized state. Recently, we developed new variants of flavin-binding protein CagFbFP exhibiting notable blue (Q148V) or red (I52V A85Q) shifts of fluorescence emission maxima. Here, we use time-resolved and low-temperature spectroscopy to show that whereas the chromophore environment is static in Q148V, an additional protein-flavin hydrogen bond is formed upon photoexcitation in the I52V A85Q variant. Consequently, in Q148V, excitation, emission, and phosphorescence spectra are shifted, whereas in I52V A85Q, excitation and low-temperature phosphorescence spectra are relatively unchanged, while emission spectrum is altered. We also determine the x-ray structures of the two variants to reveal the flavin environment and complement the spectroscopy data. Our findings illustrate two distinct color-tuning mechanisms of flavin-binding proteins and could be helpful for the engineering of new variants with improved optical properties.

Frequency-Selective Anharmonic Mode Analysis of Thermally Excited Vibrations in Proteins

Low-frequency molecular vibrations at far-infrared frequencies are thermally excited at room temperature. As a consequence, thermal fluctuations are not limited to the immediate vicinity of local minima on the potential energy surface, and anharmonic properties cannot be ignored. The latter is particularly relevant in molecules with multiple conformations, such as proteins and other biomolecules. However, existing theoretical and computational frameworks for the analysis of molecular vibrations have so far been limited by harmonic or quasi-harmonic approximations, which are ill-suited to describe anharmonic low-frequency vibrations. Here, we introduce a fully anharmonic analysis of molecular vibrations based on a time correlation formalism that eliminates the need for harmonic or quasi-harmonic approximations. We use molecular dynamics simulations of a small protein to demonstrate that this new approach, in contrast to harmonic and quasi-harmonic normal modes, correctly identifies the collective degrees of freedom associated with molecular vibrations at any given frequency. This allows us to unambiguously characterize the anharmonic character of low-frequency vibrations in the far-infrared spectrum.

Incoherent Neutron Scattering and Terahertz Time-Domain Spectroscopy on Protein and Hydration Water

Incoherent inelastic and quasi-elastic neutron scattering (INS) and terahertz time-domain spectroscopy (THz-TDS) are spectroscopy methods that directly detect molecular dynamics, with an overlap in the measured energy regions of each method. Due to the different characteristics of their probes (i.e., neutron and light), the information obtained and the sample conditions suitable for each method differ. In this review, we introduce the differences in the quantum beam properties of the two methods and their associated advantages and disadvantages in molecular spectroscopy. Neutrons are scattered via interaction with nuclei; one characteristic of neutron scattering is a large incoherent scattering cross-section of a hydrogen atom. INS records the auto-correlation functions of atomic positions. By using the difference in neutron scattering cross-sections of isotopes in multi-component systems, some molecules can be selectively observed. In contrast, THz-TDS observes the cross-correlation function of dipole moments. In water-containing biomolecular samples, the absorption of water molecules is particularly large. While INS requires large-scale experimental facilities, such as accelerators and nuclear reactors, THz-TDS can be performed at the laboratory level. In the analysis of water molecule dynamics, INS is primarily sensitive to translational diffusion motion, while THz-TDS observes rotational motion in the spectrum. The two techniques are complementary in many respects, and a combination of the two is very useful in analyzing the dynamics of biomolecules and hydration water.

Detecting the First Hydration Shell Structure around Biomolecules at Interfaces

Understanding the role of water in biological processes remains a central challenge in the life sciences. Water structures in hydration shells of biomolecules are difficult to study in situ due to overwhelming background from aqueous environments. Biological interfaces introduce additional complexity because biomolecular hydration differs at interfaces compared to bulk solution. Here, we perform experimental and computational studies of chiral sum frequency generation (chiral SFG) spectroscopy to probe chirality transfer from a protein to the surrounding water molecules. This work reveals that chiral SFG probes the first hydration shell around the protein almost exclusively. We explain the selectivity to the first hydration shell in terms of the asymmetry induced by the protein structure and specific protein-water hydrogen-bonding interactions. This work establishes chiral SFG as a powerful technique for studying hydration shell structures around biomolecules at interfaces, presenting new possibilities to address grand research challenges in biology, including the molecular origins of life.

Structure and dynamics of supercooled water in the hydration layer of poly(ethylene glycol)

The statics and dynamics of supercooled water in the hydration layer of poly(ethylene glycol) (PEG) were studied by a combination of quasi-elastic neutron scattering (QENS) and molecular dynamics (MD) simulations. Two samples, that is, hydrogenated PEG/deuterated water (h-PEG/D₂O) and fully deuterated PEG/hydrogenated water (d-PEG/H₂O) with the same molar ratio of ethylene glycol (EG) monomer to water, 1:1, are compared. The QENS data of h-PEG/D₂O show the dynamics of PEG, and that of d-PEG/H₂O reveals the motion of water. The temperature-dependent elastic scattering intensity of both samples has shown transitions at supercooled temperature, and these transition temperatures depend on the energy resolution of the instruments. Therefore, neither one is a phase transition, but undergoes dynamic process. The dynamic of water can be described as an Arrhenius to super-Arrhenius transition, and it reveals the hydrogen bonding network relaxation of hydration water around PEG at supercooled temperature. Since the PEG-water hydrogen bond structural relaxation time from MD is in good agreement with the average relaxation time from QENS (d-PEG/H₂O), MD may further reveal the atomic pictures of the supercooled hydration water. It shows that hydration water molecules form a series of pools around the hydrophilic oxygen atom of PEG. At supercooled temperature, they have a more bond ordered structure than bulk water, proceed a trapping sites diffusion on the PEG surface, and facilitate the structural relaxation of PEG backbone.

Sub-Nanosecond Dynamics of Pathologically Relevant Bio-Macromolecules Observed by Incoherent Neutron Scattering

Incoherent neutron scattering (iNS) is one of the most powerful techniques to study the dynamical behavior of bio-macromolecules such as proteins and lipid molecules or whole cells. This technique has widely been used to elucidate the fundamental aspects of molecular motions that manifest in the bio-macromolecules in relation to their intrinsic molecular properties and biological functions. Furthermore, in the last decade, iNS studies focusing on a possible relationship between molecular dynamics and biological malfunctions, i.e., human diseases and disorders, have gained importance. In this review, we summarize recent iNS studies on pathologically relevant proteins and lipids and discuss how the findings are of importance to elucidate the molecular mechanisms of human diseases and disorders that each study targets. Since some diseases such as amyloidosis have become more relevant in the aging society, research in this field will continue to develop further and be more important in the current increasing trend for longevity worldwide.

Water Networks Repopulate Protein-Ligand Interfaces with Temperature

High-resolution crystal structures highlight the importance of water networks in protein-ligand interactions. However, as these are typically determined at cryogenic temperature, resulting insights may be structurally precise but not biologically accurate. By collecting 10 matched room-temperature and cryogenic datasets of the biomedical target Hsp90α, we identified changes in water networks that impact protein conformations at the ligand binding interface. Water repositioning with temperature repopulates protein ensembles and ligand interactions. We introduce Flipper conformational barcodes to identify temperature-sensitive regions in electron density maps. This revealed that temperature-responsive states coincide with ligand-responsive regions and capture unique binding signatures that disappear upon cryo-cooling. Our results have implications for discovering Hsp90 selective ligands, and, more generally, for the utility of hidden protein and water conformations in drug discovery.

Joint neutron/molecular dynamics vibrational spectroscopy reveals softening of HIV-1 protease upon binding of a tight inhibitor

Biomacromolecules are inherently dynamic, and their dynamics are interwoven into function. The fast collective vibrational dynamics in proteins occurs in the low picosecond timescale corresponding to frequencies of ∼5-50 cm^-1. This sub-to-low THz frequency regime covers the low-amplitude collective breathing motions of a whole protein and vibrations of the constituent secondary structure elements, such as α-helices, β-sheets and loops. We have used inelastic neutron scattering experiments in combination with molecular dynamics simulations to demonstrate the vibrational dynamics softening of HIV-1 protease, a target of HIV/AIDS antivirals, upon binding of a tight clinical inhibitor darunavir. Changes in the vibrational density of states of matching structural elements in the two monomers of the homodimeric protein are not identical, indicating asymmetric effects of darunavir on the vibrational dynamics. Three of the 11 major secondary structure elements contribute over 40% to the overall changes in the vibrational density of states upon darunavir binding. Molecular dynamics simulations informed by experiments allowed us to estimate that the altered vibrational dynamics of the protease would contribute -3.6 kcal mol^-1 at 300 K, or 25%, to the free energy of darunavir binding. As HIV-1 protease drug resistance remains a concern, our results open a new avenue to help establish a direct quantitative link between protein vibrational dynamics and drug resistance.

Universal dynamical onset in water at distinct material interfaces

Interfacial water remains liquid and mobile much below 0 °C, imparting flexibility to the encapsulated materials to ensure their diverse functions at subzero temperatures. However, a united picture that can describe the dynamical differences of interfacial water on different materials and its role in imparting system-specific flexibility to distinct materials is lacking. By combining neutron spectroscopy and isotope labeling, we explored the dynamics of water and the underlying substrates independently below 0 °C across a broad range of materials. Surprisingly, while the function-related anharmonic dynamical onset in the materials exhibits diverse activation temperatures, the surface water presents a universal onset at a common temperature. Further analysis of the neutron experiment and simulation results revealed that the universal onset of water results from an intrinsic surface-independent relaxation: switching of hydrogen bonds between neighboring water molecules with a common energy barrier of ∼35 kJ mol⁻¹.

We demonstrated that the dynamical onset of interfacial water is an intrinsic property of water itself, resulting from a surface independent relaxation process in water with an approximately universal energy barrier of ∼35 kJ mol⁻¹.

Hydration of methemoglobin studied by in silico modeling and dielectric spectroscopy

The hemoglobin concentration of 35 g/dl of human red blood cells is close to the solubility threshold. Using microwave dielectric spectroscopy, we have assessed the amount of water associated with hydration shells of methemoglobin as a function of its concentration in the presence or absence of ions. We estimated water-hemoglobin interactions to interpret the obtained data. Within the concentration range of 5-10 g/dl of methemoglobin, ions play an important role in defining the free-to-bound water ratio competing with hemoglobin to recruit water molecules for the hydration shell. At higher concentrations, hemoglobin is a major contributor to the recruitment of water to its hydration shell. Furthermore, the amount of bound water does not change as the hemoglobin concentration is increased from 15 to 30 g/dl, remaining at the level of ∼20% of the total intracellular water pool. The theoretical evaluation of the ratio of free and bound water for the hemoglobin concentration in the absence of ions corresponds with the experimental results and shows that the methemoglobin molecule binds about 1400 water molecules. These observations suggest that within the concentration range close to the physiological one, hemoglobin molecules are so close to each other that their hydration shells interact. In this case, the orientation of the hemoglobin molecules is most likely not stochastic, but rather supports partial neutralization of positive and negative charges at the protein surface. Furthermore, deformation of the red blood cell shape results in the rearrangement of these structures.

Dehydration induced dynamical heterogeneity and ordering mechanism of lipid bilayers

Understanding the influence of dehydration on the membrane structure is crucial to control membrane functionality related to domain formation and cell fusion under anhydrobiosis conditions. To this end, we perform all-atom molecular dynamic simulations of 1,2-dimyristoyl-sn-glycero-3-phosphocholine dimyristoylphosphatidylcholine lipid membranes at different hydration levels at 308 K. As dehydration increases, the lipid area per head group decreases with an increase in bilayer thickness and lipid order parameters indicating bilayer ordering. Concurrently, translational and rotational dynamics of interfacial water (IW) molecules near membranes slow down. On the onset of bilayer ordering, the IW molecules exhibit prominent features of dynamical heterogeneity evident from non-Gaussian parameters and one-dimensional van Hove correlation functions. At a fully hydrated state, diffusion constants (D) of the IW follow a scaling relation, D∼τ_α ^-1, where the α relaxation time (τ_α) is obtained from self-intermediate scattering functions. However, upon dehydration, the relation breaks and the D of the IW follows a power law behavior as D∼τ_α ^-0.57, showing the signature of glass dynamics. τ_α and hydrogen bond lifetime calculated from intermittent hydrogen bond auto-correlation functions undergo a similar crossover in association with bilayer ordering on dehydration. The bilayer ordering is accompanied with an increase in fraction of caged lipids spanned over the bilayer surface and a decrease in fraction of mobile lipids due to the non-diffusive dynamics. Our analyses reveal that the microscopic mechanism of lipid ordering by dehydration is governed by dynamical heterogeneity. The fundamental understanding from this study can be applied to complex bio-membranes to trap functionally relevant gel-like domains at room temperature.

Photo-Switching of Protein Dynamical Collectivity

We examine changes in the picosecond structural dynamics with irreversible photobleaching of red fluorescent proteins (RFP) mCherry, mOrange2 and TagRFP-T. Measurements of the protein dynamical transition using terahertz time-domain spectroscopy show in all cases an increase in the turn-on temperature in the bleached state. The result is surprising given that there is little change in the protein surface, and thus, the solvent dynamics held responsible for the transition should not change. A spectral analysis of the measurements guided by quasiharmonic calculations of the protein absorbance reveals that indeed the solvent dynamical turn-on temperature is independent of the thermal stability/photostate however the protein dynamical turn-on temperature shifts to higher temperatures. This is the first demonstration of switching the protein dynamical turn-on temperature with protein functional state. The observed shift in protein dynamical turn-on temperature relative to the solvent indicates an increase in the required mobile waters necessary for the protein picosecond motions, that is, these motions are more collective. Melting-point measurements reveal that the photobleached state is more thermally stable, and structural analysis of related RFP’s shows that there is an increase in internal water channels as well as a more uniform atomic root mean squared displacement. These observations are consistent with previous suggestions that water channels form with extended light excitation providing O2 access to the chromophore and subsequent fluorescence loss. We report that these same channels increase internal coupling enhancing thermal stability and collectivity of the picosecond protein motions. The terahertz spectroscopic characterization of the protein and solvent dynamical onsets can be applied generally to measure changes in collectivity of protein motions.

Effect of trehalose on protein cryoprotection: Insights into the mechanism of slowing down of hydration water

We study, with molecular dynamics simulations, a lysozyme protein immersed in a water-trehalose solution upon cooling. The aim is to understand the cryoprotectant role played by this disaccharide through the modifications that it induces on the slow dynamics of protein hydration water with its presence. The α-relaxation shows a fragile to strong crossover about 20° higher than that in the bulk water phase and 15° higher than that in lysozyme hydration water without trehalose. The protein hydration water without trehalose was found to show a second slower relaxation exhibiting a strong to strong crossover coupled with the protein dynamical transition. This slower relaxation time importantly appears enormously slowed down in our cryoprotectant solution. On the other hand, this long-relaxation in the presence of trehalose is also connected with a stronger damping of the protein structural fluctuations than that found when the protein is in contact with the pure hydration water. Therefore, this appears to be the mechanism through which trehalose manifests its cryoprotecting function.

Zinc determines dynamical properties and aggregation kinetics of human insulin

Protein aggregation is a widespread process leading to deleterious consequences in the organism, with amyloid aggregates being important not only in biology but also for drug design and biomaterial production. Insulin is a protein largely used in diabetes treatment, and its amyloid aggregation is at the basis of the so-called insulin-derived amyloidosis. Here, we uncover the major role of zinc in both insulin dynamics and aggregation kinetics at low pH, in which the formation of different amyloid superstructures (fibrils and spherulites) can be thermally induced. Amyloid aggregation is accompanied by zinc release and the suppression of water-sustained insulin dynamics, as shown by particle-induced x-ray emission and x-ray absorption spectroscopy and by neutron spectroscopy, respectively. Our study shows that zinc binding stabilizes the native form of insulin by facilitating hydration of this hydrophobic protein and suggests that introducing new binding sites for zinc can improve insulin stability and tune its aggregation propensity.

The Dynamics of Hydrated Proteins Are the Same as Those of Highly Asymmetric Mixtures of Two Glass-Formers

Customarily, the studies of dynamics of hydrated proteins are focused on the fast hydration water ν-relaxation, the slow structural α-relaxation responsible for a single glass transition, and the protein dynamic transition (PDT). Guided by the analogy with the dynamics of highly asymmetric mixtures of molecular glass-formers, we explore the possibility that the dynamics of hydrated proteins are richer than presently known. By providing neutron scattering, dielectric relaxation, calorimetry, and deuteron NMR data in two hydrated globular proteins, myoglobin and BSA, and the fibrous elastin, we show the presence of two structural α-relaxations, α1 and α2, and the hydration water ν-relaxation, all coupled together with interconnecting properties. There are two glass transition temperatures T _g ^α1and T _g ^α2 corresponding to vitrification of the α1 and α2 processes. Relaxation time τ_α2(T) of the α2-relaxation changes its Arrhenius temperature dependence to super-Arrhenius on crossing T _g ^α1 from below. The ν-relaxation responds to the two vitrifications by changing the T-dependence of its relaxation time τ_ν(T) on crossing consecutively T _g ^α2 and T _g ^α1. It generates the PDT at T _d where τ_ν(T _d) matches about five times the experimental instrument timescale τ_exp, provided that T _d > T _g ^α1. This condition is satisfied by the hydrated globular proteins considered in this paper, and the ν-relaxation is in the liquid state with τ_ν(T) having the super-Arrhenius temperature dependence. However, if T _d < T _g ^α1, the ν-relaxation fails to generate the PDT because it is in the glassy state and τ_ν(T) has Arrhenius T-dependence, as in the case of hydrated elastin. Overall, the dynamics of hydrated proteins are the same as the dynamics of highly asymmetric mixtures of glass-formers. The results from this study have expanded the knowledge of the dynamic processes and their properties in hydrated proteins, and impact on research in this area is expected.

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.

How can we derive hydration water dynamics with incoherent neutron scattering and molecular dynamics simulation?

Incoherent neutron scattering (INS) is one of the useful experimental methods for studying protein dynamics at the pico-nanosecond timescale. At this timescale, protein dynamics is highly coupled with hydration, which is observed as protein dynamical transition (PDT). INS is very sensitive to hydrogen atomic dynamics because of the large incoherent scattering cross section of hydrogen atom, and thus, the INS of a hydrated protein provides overall dynamic information about the protein, including hydration water. Separation of hydration water dynamics is essential for understanding hydration-related protein dynamics. H₂O/D₂O exchange is an effective method in the context of INS experiments for observing the dynamics of protein and hydration water separately. Neutron scattering is directly related to the van Hove space-time correlation function, which can be calculated quantitatively by performing molecular dynamics (MD) simulations. Diffusion and hydrogen bond dynamics of hydration water can be analyzed by performing MD simulation. MD simulation is useful for analyzing the dynamic coupling mechanism in hydration-related protein dynamics from the viewpoint of interpreting INS data because PDT is induced by hydration. In the present work, we demonstrate the methodological advantages of the H₂O/D₂O exchange technique in INS and the compatibility of INS and MD simulation as tools for studying protein dynamics and hydration water.

Probing the dynamics of biological matter by elastic, quasi-elastic, and inelastic neutron scattering

The so-called function-structure-dynamics paradigm established that a close relationship links the way biological molecules work (function), their 3-dimensional organization (structure) and the changes of this organization in time (dynamics), which characterize biomolecules as highly dynamic objects. A typical example of protein dynamics is provided by protein reactions with substrates: equilibrium thermal fluctuations of protein structure are necessary to allow the access of substrates to the active site, where the functional reaction occurs. Neutron scattering is a powerful technique to study equilibrium protein structural dynamics. The incoherent structure factor, which is dominant in neutron scattering from biological matter, is related to the time-position self correlation function of protein/solvent nuclei. Here the basic theory of neutron scattering and the principles of the technologies used to measure it are described. Some selected applications of neutron scattering for investigating the structural dynamics of biological molecules are also reviewed.

Role of hydration water in the onset of protein structural dynamics

Proteins are the molecular workhorses in a living organism. Their 3D structures are animated by a multitude of equilibrium fluctuations and specific out-of-equilibrium motions that are required for proteins to be biologically active. When studied as a function of temperature, functionally relevant dynamics are observed at and above the so-called protein dynamical transition (~240 K) in hydrated, but not in dry proteins. In this review we present and discuss the main experimental and computational results that provided evidence for the dynamical transition, with a focus on the role of hydration water dynamics in sustaining functional protein dynamics. The coupling and mutual influence of hydration water dynamics and protein dynamics are discussed and the hypotheses illustrated that have been put forward to explain the physical origin of their onsets.

The Hydration Shell of Monomeric and Dimeric Insulin Studied by Terahertz Time-Domain Spectroscopy

Protein aggregation is believed to be a significant biological mechanism related to neurodegenerative disease, which makes the early-stage detection of aggregates a major concern. We demonstrated the use of terahertz (THz) time-domain spectroscopy to study protein-water interaction of monomeric and dimeric bovine insulin in aqueous samples. Regulated by changing pH and verified by size-exclusion chromatography and dynamic light scattering, we then measured their concentration-dependent changes in THz absorption between 0.5 and 3.0 THz and quantitatively deduced the extended hydration shell thickness by cubic distribution model and random distribution model. Under a random distribution assumption, the extended hydration thickness is 15.4 ± 0.4 Å for monomeric insulin and 17.5 ± 0.5 Å for dimeric insulin, with the hydration number of 6700 and 11,000, respectively. The hydration number of dimeric insulin is not twice but 1.64 times that of monomeric insulin, further supported by the ratio of solvent-accessible surface area. This "1.64-times" relation probably originates from the structural and conformational changes accompanied with dimerization. Combined with the investigations on insulin samples with different single amino acid mutations, residue B24 is believed to play an important role in the dimerization process. It is demonstrated that THz time-domain spectroscopy is a useful tool and has the sensitivity to provide the hydration information of different protein aggregates at an early stage.

Characterization of hydration water in supercooled water-trehalose solutions: The role of the hydrogen bonds network

The structural and dynamical properties of hydration water in aqueous solutions of trehalose are studied with molecular dynamics simulation. We simulate the systems in the supercooled region to investigate how the interaction with the trehalose molecules modifies the hydrogen bond network, the structural relaxation, and the diffusion properties of hydration water. The analysis is performed by considering the radial distribution functions, the residence time of water molecules in the hydration shell, the two body excess entropy, and the hydrogen bond water-water and water-trehalose correlations of the hydration water. The study of the two body excess entropy shows the presence of a fragile to strong crossover in supercooled hydration water also found in the relaxation time of the water-water hydrogen bond correlation function, and this is in agreement with predictions of the mode coupling theory and of previous studies of the oxygen-oxygen density correlators [A. Iorio et al., J. Mol. Liq. 282, 617 (2019); Sci. China: Phys., Mech. Astron. 62, 107011 (2019)]. The water-trehalose hydrogen bond correlation function instead evidences a strong to strong crossover in the relaxation time, and this crossover is related to a trehalose dynamical transition. This signals the role that the strong interplay between the soluted molecules and the surrounding solvent has in determining the dynamical transition common to both components of the system that happens upon cooling and that is similar to the well known protein dynamical transition. We connect our results with the cryoprotecting role of trehalose molecules.

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.

Temperature Induced Dynamical Transition of Biomolecules in Polarizable and Nonpolarizable TIP3P Water

Temperature induced dynamical transition (DT), associated with a sharp rise in molecular flexibility, is well-known to be exhibited between 270 and 280 K in glycerol to 200-230 K in hydrated biomolecules and is controlled by diffusivity (viscosity) of the solvation layer. In the molecular dynamics (MD) community, especially for water as a solvent, this has been an intense area of research despite decades of investigations. However, in general, water in these studies is described by empirical nonpolarizable force fields in which electronic polarizability is treated implicitly with effective charges and related parameters. This might have led to the present trait of discovery that DTs of biomolecules, irrespective of the potential functions for water models used, occur within a narrow band of temperature variation (30-40 K). Whereas a water molecule in a biomolecular surface and one in bulk are polarized differently, therefore explicit treatment of water polarizability would be a powerful approach toward the treatment of hydration water, believed to cause the DT manifestation. Using MD simulations, we investigated the effects of polarizable water on the DT of biomolecules and the dynamic properties of hydration water. We chose two types of solutes: globular protein (lysozyme) and more open and flexible RNAs (a hairpin and a riboswitch) with different natures of hydrophilic sites than proteins in general. We found that the characteristic temperature of DT ( T_DT) for the solutes in polarizable water is always higher than that in its nonpolarizable counterpart. In particular, for RNAs, the variations are found to be ∼45 K between the two water models, whereas for the more compact lysozyme, it is only ∼4 K. The results are discussed in light of the enormous increase in relaxation times of a liquid upon cooling in the paradigm of dynamic switchover in hydration water with liquid-liquid phase transition, derived from the existence of the second critical point. Our result supports the idea that structures of biomolecules and their interactions with the hydration water determines T_DT and provides evidence for the decisive role of polarizable water on the onset of DT, which has been hitherto ignored.

Protein-Polymer Dynamics as Affected by Polymer Coating and Interactions

We investigate the relaxation dynamics of protein-polymer conjugates by neutron scattering spectroscopy to understand to which extent the coating of a protein by a polymer can replace water in promoting thermal structural fluctuations. For this purpose, we compare the dynamics of protein-polymer mixtures to that of conjugates with a variable number of polymers covalently attached to the protein. Results show that the flexibility of the protein is larger in protein-polymer mixtures than in native protein or in conjugates, even in the dry state. Upon hydration, both the native protein and the conjugate show equivalent dynamics, suggesting that the polymer grafted on the protein surface adsorbs all water molecules.

Role of Mono- and Divalent Surface Cations on the Structure and Adsorption Behavior of Water on Mica Surface

Understanding solid-water(vapor) interfacial systems is relevant for both industrial and academic scenarios for their presence in wide areas ranging from tribology to geochemistry. Using grand canonical Monte Carlo simulations, we have investigated the role of monovalent (lithium, Li⁺; sodium, Na⁺; and potassium, K⁺) and divalent (magnesium, Mg²⁺; calcium, Ca²⁺) cations on the structure and adsorption behavior of water on mica surface. The water density adjacent to the surface exhibits (a) oscillations due to hydration of surface cations (interfacial layer), (b) followed by a thick liquidlike layer. The thickness of the interfacial layer is strongly dependent on the hydration shell size and hydration energy of surface ions. Water molecules immediately next to the surface (contact layers) adsorb on ditrigonal (hexagonal) cavities of mica surface and form an ordered structure. The Li⁺, Na⁺, Mg²⁺, and Ca²⁺ surface ions are coadsorbed with water molecules on the ditrigonal cavities due to their smaller hydration shell. Majority of water molecules in the second contact layer hydrate the surface ions and, together with the rest of the water molecules, form hydrogen bonds among themselves. The structure of the water molecules in the third and subsequent layer is random and more bulk liquidlike, except those molecules that hydrate the surface ions. The adsorption isotherm of water on all ion-exposed mica surface exhibits three regimes: (a) an initial rapid increase in water loading for relative vapor pressure ( p/ p₀) ≤0.2 due to hydration of surface ions; (b) followed by a linear increase between p/ p₀ = 0.2 and 0.7, where the hydrogen bond formation between the water molecules of the interfacial layer occurs; and (c) exponential growth beyond p/ p₀ = 0.7 due to thickening of the liquidlike layer. The water loading on divalent-ion-exposed mica surface is higher compared to the monovalent ions case. Although the divalent ions have higher hydration energy, the fraction of water molecules hydrating the surface ions is less compared to nonhydrating water molecules. We found that ion hydration energy and size of hydration shell play a crucial role in their structure adjacent to mica surface. At lower water loadings, the surface ions form a hydration shell with surface bridging oxygens, whereas at higher water content, the hydration preference is shifted toward mobile water molecules. The detailed understanding obtained from this work will be useful in identifying the role of ions in cloud formation, nanotribological studies, and activities of biological molecules and catalysts.

Water-Protein Interactions Coupled with Protein Conformational Transition

Conformational fluctuations of proteins are crucially important for their functions. However, changes in the location and dynamics of hydrated water in many proteins accompanied by the conformational transition have not been fully understood. Here, we used phase-modulated clean chemical exchange NMR approach to investigate pressure-induced changes in water-to-amide proton exchange occurring at sub-second time scale. With the transition of ubiquitin from its native conformation (N1) to an alternative conformation (N2) at 250 MPa, proton exchange rates of residues 32-35, 40-41, and 71, which are located at the C-terminal side of the protein, were significantly increased. These observations can be explained by the destabilization of the hydrogen bonds in the backbone and partial exposure of those amide groups to solvent in N2. We conclude that phase-modulated clean chemical exchange NMR approach coupled with pressure perturbation will be a useful tool for investigations of more open and hydrated protein structures.

Evidence of a low-temperature dynamical transition in concentrated microgels

A low-temperature dynamical transition has been reported in several proteins. We provide the first observation of a "protein-like" dynamical transition in nonbiological aqueous environments. To this aim, we exploit the popular colloidal system of poly-N-isopropylacrylamide (PNIPAM) microgels, extending their investigation to unprecedentedly high concentrations. Owing to the heterogeneous architecture of the microgels, water crystallization is avoided in concentrated samples, allowing us to monitor atomic dynamics at low temperatures. By elastic incoherent neutron scattering and molecular dynamics simulations, we find that a dynamical transition occurs at a temperature T _d ~ 250 K, independently from PNIPAM mass fraction. However, the transition is smeared out on approaching dry conditions. The quantitative agreement between experiments and simulations provides evidence that the transition occurs simultaneously for PNIPAM and water dynamics. The similarity of these results with hydrated protein powders suggests that the dynamical transition is a generic feature in complex macromolecular systems, independently from their biological function.

Dynamic Neutron Scattering by Biological Systems

Dynamic neutron scattering directly probes motions in biological systems on femtosecond to microsecond timescales. When combined with molecular dynamics simulation and normal mode analysis, detailed descriptions of the forms and frequencies of motions can be derived. We examine vibrations in proteins, the temperature dependence of protein motions, and concepts describing the rich variety of motions detectable using neutrons in biological systems at physiological temperatures. New techniques for deriving information on collective motions using coherent scattering are also reviewed.

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.

Dynamical Transition of Collective Motions in Dry Proteins

Water is widely assumed to be essential for protein dynamics and function. In particular, the well-documented "dynamical" transition at ∼200  K, at which the protein changes from a rigid, nonfunctional form to a flexible, functional state, as detected in hydrogenated protein by incoherent neutron scattering, requires hydration. Here, we report on coherent neutron scattering experiments on perdeuterated proteins and reveal that a transition occurs in dry proteins at the same temperature resulting primarily from the collective heavy-atom motions. The dynamical transition discovered is intrinsic to the energy landscape of dry proteins.

Distinct Protein Hydration Water Species Defined by Spatially Resolved Spectra of Intermolecular Vibrations

In this molecular dynamics simulation study, we analyze intermolecular vibrations in the hydration shell of a solvated enyzme, the membrane type 1–matrix metalloproteinase, with high spatial resolution. Our approach allows us to characterize vibrational signatures of the local hydrogen bond network, the translational mobility of water molecules, as well as the molecular entropy, in specific local environments. Our study demonstrates the heterogeneity of water properties within the hydration shell of a complex biomolecule. We define a classification scheme based on the vibrational density of states that allows us to distinguish separate classes of hydration water species and facilitates the description of hydration water properties at distinct hydration sites. The results demonstrate that no single characteristic of the protein surface is sufficient to determine the properties of nearby water. The protein surface geometry, quantified here by the number of protein atoms in the vicinity of a hydration water molecule, as well as the chemical nature of a solvated protein functional group, influences dynamic and thermodynamic properties of solvating water molecules.

Thermal activation of 'allosteric-like' large-scale motions in a eukaryotic Lactate Dehydrogenase

Conformational changes occurring during the enzymatic turnover are essential for the regulation of protein functionality. Individuating the protein regions involved in these changes and the associated mechanical modes is still a challenge at both experimental and theoretical levels. We present here a detailed investigation of the thermal activation of the functional modes and conformational changes in a eukaryotic Lactate Dehydrogenase enzyme (LDH). Neutron Spin Echo spectroscopy and Molecular Dynamics simulations were used to uncover the characteristic length- and timescales of the LDH nanoscale motions in the apo state. The modes involving the catalytic loop and the mobile region around the binding site are activated at room temperature, and match the allosteric reorganisation of bacterial LDHs. In a temperature window of about 15 degrees, these modes render the protein flexible enough and capable of reorganising the active site toward reactive configurations. On the other hand an excess of thermal excitation leads to the distortion of the protein matrix with a possible anti-catalytic effect. Thus, the temperature activates eukaryotic LDHs via the same conformational changes observed in the allosteric bacterial LDHs. Our investigation provides an extended molecular picture of eukaryotic LDH's conformational landscape that enriches the static view based on crystallographic studies alone.

Direct Observation of Nanosecond Water Exchange Dynamics at a Protein Metal Site

Nanosecond ligand exchange dynamics at metal sites within proteins is essential in catalysis, metal ion transport, and regulatory metallobiochemistry. Herein we present direct observation of the exchange dynamics of water at a Cd²⁺ binding site within two de novo designed metalloprotein constructs using ^111mCd perturbed angular correlation (PAC) of γ-rays and ¹¹³Cd NMR spectroscopy. The residence time of the Cd²⁺-bound water molecule is tens of nanoseconds at 20 °C in both proteins. This constitutes the first direct experimental observation of the residence time of Cd²⁺ coordinated water in any system, including the simple aqua ion. A Leu to Ala amino acid substitution ∼10 Å from the Cd²⁺ site affects both the equilibrium constant and the residence time of water, while, surprisingly, the metal site structure, as probed by PAC spectroscopy, remains essentially unaltered. This implies that remote mutations may affect metal site dynamics, even when structure is conserved.

Neutrons describe ectoine effects on water H-bonding and hydration around a soluble protein and a cell membrane

Understanding adaptation to extreme environments remains a challenge of high biotechnological potential for fundamental molecular biology. The cytosol of many microorganisms, isolated from saline environments, reversibly accumulates molar concentrations of the osmolyte ectoine to counterbalance fluctuating external salt concentrations. Although they have been studied extensively by thermodynamic and spectroscopic methods, direct experimental structural data have, so far, been lacking on ectoine-water-protein interactions. In this paper, in vivo deuterium labeling, small angle neutron scattering, neutron membrane diffraction and inelastic scattering are combined with neutron liquids diffraction to characterize the extreme ectoine-containing solvent and its effects on purple membrane of H. salinarum and E. coli maltose binding protein. The data reveal that ectoine is excluded from the hydration layer at the membrane surface and does not affect membrane molecular dynamics, and prove a previous hypothesis that ectoine is excluded from a monolayer of dense hydration water around the soluble protein. Neutron liquids diffraction to atomic resolution shows how ectoine enhances the remarkable properties of H-bonds in water—properties that are essential for the proper organization, stabilization and dynamics of biological structures.

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.