Hydration affects both harmonic and anharmonic nature of protein dynamics

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.

Response of Terahertz Protein Vibrations to Ligand Binding: Calmodulin-Peptide Complexes as a Case Study

Terahertz vibrations are sensitive reporters of the structure and interactions of proteins. Ligand binding alters the nature and distribution of these collective vibrations. The ligand-induced changes in the terahertz protein vibrations contribute to the binding entropy and to the overall thermodynamic stability of the resultant protein-ligand complexes. Here, we have examined the response of the low-frequency (below 6 terahertz) collective vibrations of the calcium-loaded calmodulin (CaM) to binding to five different ligands, both in the presence and absence of water, using normal-mode analysis and molecular dynamics simulations. A comparison of the vibrational spectra of hydrated and dry systems reveals that protein-solvent interactions stiffen the terahertz protein vibrations and that these solvent-coupled collective vibrations contribute significantly to the hydration-sensitive variation in the vibrational entropy of CaM. In the absence of water, the low-frequency vibrations of CaM are stiffened by ligand binding. On the contrary, the number and the cumulative vibrational entropy of low-frequency vibrational modes (ω < 200 cm^-1) of the hydrated CaM are increased noticeably after binding to the peptides, indicating binding-induced softening of collective vibrations of the protein. Although the calculated and experimental binding affinities of the chosen complexes correlated reasonably well, no systematic correlation was observed between the protein vibrational entropy and the binding affinity. The results underscored the importance of the interplay of protein-ligand and solvent interactions in modulating the low-frequency vibrations of proteins.

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.

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.

Freezable and Unfreezable Hydration Water: Distinct Contributions to Protein Dynamics Revealed by Neutron Scattering

Hydration water plays a crucial role for activating the protein dynamics required for functional expression. Yet, the details are not understood about how hydration water couples with protein dynamics. A temperature hysteresis of the ice formation of hydration water is a key phenomenon to understand which type of hydration water, unfreezable or freezable hydration water, is crucial for the activation of protein dynamics. Using neutron scattering, we observed a temperature-hysteresis phenomenon in the diffraction peaks of the ice of freezable hydration water, whereas protein dynamics did not show any temperature hysteresis. These results show that the protein dynamics is not coupled with freezable hydration water dynamics, and unfreezable hydration water is essential for the activation of protein dynamics. Decoupling of the dynamics between unfreezable and freezable hydration water could be the cause of the distinct contributions of hydration water to protein dynamics.

Molecular Recognition and Self-Organization in Life Phenomena Studied by a Statistical Mechanics of Molecular Liquids, the RISM/3D-RISM Theory

There are two molecular processes that are essential for living bodies to maintain their life: the molecular recognition, and the self-organization or self-assembly. Binding of a substrate by an enzyme is an example of the molecular recognition, while the protein folding is a good example of the self-organization process. The two processes are further governed by the other two physicochemical processes: solvation and the structural fluctuation. In the present article, the studies concerning the two molecular processes carried out by Hirata and his coworkers, based on the statistical mechanics of molecular liquids or the RISM/3D-RISM theory, are reviewed.

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.

Cancellation between auto- and mutual correlation contributions of protein/water dynamics in terahertz time-domain spectra

Terahertz time-domain spectra (THz-TDS) were investigated using the results of molecular dynamics (MD) simulations of Staphylococcal nuclease at two hydration states in the temperature range between 100 and 300 K. The temperature dependence of THz-TDS was found to differ significantly from that of the incoherent neutron scattering spectra (INSS) calculated from the same MD simulation results. We further examined contributions of the mutual and auto-correlations of the atomic fluctuations to THz-TDS and found that the negative value of the former contribution nearly canceled out the positive value of the latter, resulting in a monotonic increase of the reduced absorption cross section. Because of this cancellation, no distinct broad peak was observed in the absorption lineshape function of THz-TDS, whereas the protein boson peak was observed in INSS. The contribution of water molecules to THz-TDS was extremely large for the hydrated protein at temperatures above 200 K, in which large-amplitude motions of water were excited. The combination of THz-TDS, INSS and MD simulations has the potential to extract function-relevant protein dynamics occurring on the picosecond to nanosecond timescale.

Universality and Structural Implications of the Boson Peak in Proteins

The softness and rigidity of proteins are reflected in the structural dynamics, which are in turn affected by the environment. The characteristic low-frequency vibrational spectrum of a protein, known as boson peak, is an indication of the structural rigidity of the protein at a cryogenic temperature or dehydrated conditions. In this article, the effect of hydration, temperature, and pressure on the boson peak and volumetric properties of a globular protein are evaluated by using inelastic neutron scattering and molecular dynamics simulation. Hydration, pressurization, and cooling shift the boson peak position to higher energy and depress the peak intensity and decreases the protein and cavity volumes. We found the correlation between the boson peak and cavity volume in a protein. A decrease of cavity volume means the increase of rigidity, which is the origin of the boson peak shift. Boson peak is the universal property of a protein, which is rationalized by the correlation.

The unusual internal motion of the villin headpiece subdomain

The thermostable 36-residue subdomain of the villin headpiece (HP36) is the smallest known cooperatively folding protein. Although the folding and internal dynamics of HP36 and close variants have been extensively studied, there has not been a comprehensive investigation of side-chain motion in this protein. Here, the fast motion of methyl-bearing amino acid side chains is explored over a range of temperatures using site-resolved solution nuclear magnetic resonance deuterium relaxation. The squared generalized order parameters of methyl groups extensively spatially segregate according to motional classes. This has not been observed before in any protein studied using this methodology. The class segregation is preserved from 275 to 305 K. Motions detected in Helix 3 suggest a fast timescale of conformational heterogeneity that has not been previously observed but is consistent with a range of folding and dynamics studies. Finally, a comparison between the order parameters in solution with previous results based on solid-state nuclear magnetic resonance deuterium line shape analysis of HP36 in partially hydrated powders shows a clear disagreement for half of the sites. This result has significant implications for the interpretation of data derived from a variety of approaches that rely on partially hydrated protein samples.

Rigidity, secondary structure, and the universality of the boson peak in proteins

Complementary neutron- and light-scattering results on nine proteins and amino acids reveal the role of rigidity and secondary structure in determining the time- and lengthscales of low-frequency collective vibrational dynamics in proteins. These dynamics manifest in a spectral feature, known as the boson peak (BP), which is common to all disordered materials. We demonstrate that BP position scales systematically with structural motifs, reflecting local rigidity: disordered proteins appear softer than α-helical proteins; which are softer than β-sheet proteins. Our analysis also reveals a universal spectral shape of the BP in proteins and amino acid mixtures; superimposable on the shape observed in typical glasses. Uniformity in the underlying physical mechanism, independent of the specific chemical composition, connects the BP vibrations to nanometer-scale heterogeneities, providing an experimental benchmark for coarse-grained simulations, structure/rigidity relationships, and engineering of proteins for novel applications.

Local dynamics coupled to hydration water determines DNA-sequence-dependent deformability

Molecular dynamics (MD) simulations and quasielastic neutron scattering (QENS) experiments were conducted on two hydrated DNA dodecamers with distinct deformability: 5'CGCG[under AATT]̲CGCG3' and 5'CGCG[under TTAA]̲CGCG3'. The former is known to be rigid and the latter to be flexible. The mean-square displacements of DNA dodecamers exhibit so-called dynamical transition around 200-240 K for both sequences. To investigate the DNA-sequence-dependent dynamics, the dynamics of DNA and hydration water above the transition temperature were examined using both MD simulations and QENS experiments. The fluctuation amplitude of the AATT central tetramer is smaller, and its relaxation time is longer, than that observed in TTAA, suggesting that the AT step is kinetically more stable than TA. The sequence-dependent local base pair step dynamics correlates with the kinetics of breaking the hydrogen bond between DNA and hydration water. The sequence-dependent DNA base pair step fluctuations appear above the dynamical transition temperature. Together with these results, we conclude that DNA deformability is related to the local dynamics of the base pair steps, themselves coupled to hydration water in the minor groove.

Mode coupling points to functionally important residues in myosin II

Relevance of mode coupling to energy/information transfer during protein function, particularly in the context of allosteric interactions is widely accepted. However, existing evidence in favor of this hypothesis comes essentially from model systems. We here report a novel formal analysis of the near-native dynamics of myosin II, which allows us to explore the impact of the interaction between possibly non-Gaussian vibrational modes on fluctutational dynamics. We show that an information-theoretic measure based on mode coupling alone yields a ranking of residues with a statistically significant bias favoring the functionally critical locations identified by experiments on myosin II.

Specific ions modulate diffusion dynamics of hydration water on lipid membrane surfaces

Effects of specific ions on the local translational diffusion of water near large hydrophilic lipid vesicle surfaces were measured by Overhauser dynamic nuclear polarization (ODNP). ODNP relies on an unpaired electron spin-containing probe located at molecular or surface sites to report on the dynamics of water protons within ∼10 Å from the spin probe, which give rise to spectral densities for electron–proton cross-relaxation processes in the 10 GHz regime. This pushes nuclear magnetic resonance relaxometry to more than an order of magnitude higher frequencies than conventionally feasible, permitting the measurement of water moving with picosecond to subnanosecond correlation times. Diffusion of water within ∼10 Å of, i.e., up to ∼3 water layers around the spin probes located on hydrophilic lipid vesicle surfaces is ∼5 times retarded compared to the bulk water translational diffusion. This directly reflects on the activation barrier for surface water diffusion, i.e., how tightly water is bound to the hydrophilic surface and surrounding waters. We find this value to be modulated by the presence of specific ions in solution, with its order following the known Hofmeister series. While a molecular description of how ions affect the hydration structure at the hydrophilic surface remains to be answered, the finding that Hofmeister ions directly modulate the surface water diffusivity implies that the strength of the hydrogen bond network of surface hydration water is directly modulated on hydrophilic surfaces.

Protein surface and core dynamics show concerted hydration-dependent activation

By specifically labeling leucine/valine methyl groups and lysine side chains "inside" and "outside" dynamics of proteins on the nanosecond timescale are compared using neutron scattering. Surprisingly, both groups display similar dynamics as a function of temperature, and the buried hydrophobic core is sensitive to hydration and undergoes a dynamical transition.

Dynamic behavior of oligomeric inorganic pyrophosphatase explored by quasielastic neutron scattering

The purpose of this investigation is to determine whether a large oligomeric protein, inorganic pyrophosphatase (IPPase) from Thermococcus thioreducens with quaternary structural complexity, would have distinguishable dynamic characteristics compared to those of the small simple monomeric model protein, lysozyme. In this study, the β-relaxational dynamics of the two proteins, IPPase and lysozyme, are compared in the 10 ps to 0.5 ns time interval using quasi-elastic neutron scattering (QENS). Both of the protein dynamics show a characteristic logarithmic-like decay in the intermediate scattering function (ISF) of the hydrogen atoms. Distinguishable dynamical behavior found between two proteins reveals local flexibility and conformational substates unique to oligomeric structures. Moreover, the temperature dependence of the mean square displacement (MSD) of the hydrogen atoms in protein molecules, which is a traditional way to determine the "softness" of the protein molecule, is measured and shows no difference for the two proteins.

Thermal fluctuations of haemoglobin from different species: adaptation to temperature via conformational dynamics

Thermodynamic stability, configurational motions and internal forces of haemoglobin (Hb) of three endotherms (platypus, Ornithorhynchus anatinus; domestic chicken, Gallus gallus domesticus and human, Homo sapiens) and an ectotherm (salt water crocodile, Crocodylus porosus) were investigated using circular dichroism, incoherent elastic neutron scattering and coarse-grained Brownian dynamics simulations. The experimental results from Hb solutions revealed a direct correlation between protein resilience, melting temperature and average body temperature of the different species on the 0.1 ns time scale. Molecular forces appeared to be adapted to permit conformational fluctuations with a root mean square displacement close to 1.2 Å at the corresponding average body temperature of the endotherms. Strong forces within crocodile Hb maintain the amplitudes of motion within a narrow limit over the entire temperature range in which the animal lives. In fully hydrated powder samples of human and chicken, Hb mean square displacements and effective force constants on the 1 ns time scale showed no differences over the whole temperature range from 10 to 300 K, in contrast to the solution case. A complementary result of the study, therefore, is that one hydration layer is not sufficient to activate all conformational fluctuations of Hb in the pico- to nanosecond time scale which might be relevant for biological function. Coarse-grained Brownian dynamics simulations permitted to explore residue-specific effects. They indicated that temperature sensing of human and chicken Hb occurs mainly at residues lining internal cavities in the β-subunits.

Protein dynamical transition at 110 K

Proteins are known to undergo a dynamical transition at around 200 K but the underlying mechanism, physical origin, and relationship to water are controversial. Here we report an observation of a protein dynamical transition as low as 110 K. This unexpected protein dynamical transition precisely correlated with the cryogenic phase transition of water from a high-density amorphous to a low-density amorphous state. The results suggest that the cryogenic protein dynamical transition might be directly related to the two liquid forms of water proposed at cryogenic temperatures.

Quasielastic neutron scattering study of hydrogen motions in an aqueous poly(vinyl methyl ether) solution

We present a quasielastic neutron scattering (QENS) investigation of the component dynamics in an aqueous Poly(vinyl methyl ether) (PVME) solution (30% water content in weight). In the glassy state, an important shift in the Boson peak of PVME is found upon hydration. At higher temperatures, the diffusive-like motions of the components take place with very different characteristic times, revealing a strong dynamic asymmetry that increases with decreasing T. For both components, we observe stretching of the scattering functions with respect to those in the bulk and non-Gaussian behavior in the whole momentum transfer range investigated. To explain these observations we invoke a distribution of mobilities for both components, probably originated from structural heterogeneities. The diffusive-like motion of PVME in solution takes place faster and apparently in a more continuous way than in bulk. We find that the T-dependence of the characteristic relaxation time of water changes at T ≲ 225 K, near the temperature where a crossover from a low temperature Arrhenius to a high temperature cooperative behavior has been observed by broadband dielectric spectroscopy (BDS) [S. Cerveny, J. Colmenero and A. Alegría, Macromolecules, 38, 7056 (2005)]. This observation might be a signature of the onset of confined dynamics of water due to the freezing of the PVME dynamics, that has been selectively followed by these QENS experiments. On the other hand, revisiting the BDS results on this system we could identify an additional "fast" process that can be attributed to water motions coupled with PVME local relaxations that could strongly affect the QENS results. Both kinds of interpretations, confinement effects due to the increasing dynamic asymmetry and influence of localized motions, could provide alternative scenarios to the invoked "strong-to-fragile" transition.

Redox-promoting protein motions in rubredoxin

Proteins are dynamic objects, constantly undergoing conformational fluctuations, yet the linkage between internal protein motion and function is widely debated. This study reports on the characterization of temperature-activated collective and individual atomic motions of oxidized rubredoxin, a small 53 residue protein from thermophilic Pyrococcus furiosus (RdPf). Computational modeling allows detailed investigations of protein motions as a function of temperature, and neutron scattering experiments are used to compare to computational results. Just above the dynamical transition temperature which marks the onset of significant anharmonic motions of the protein, the computational simulations show both a significant reorientation of the average electrostatic force experienced by the coordinated Fe(3+) ion and a dramatic rise in its strength. At higher temperatures, additional anharmonic modes become activated and dominate the electrostatic fluctuations experienced by the ion. At 360 K, close to the optimal growth temperature of P. furiosus, simulations show that three anharmonic modes including motions of two conserved residues located at the protein active site (Ile7 and Ile40) give rise to the majority of the electrostatic fluctuations experienced by the Fe(3+) ion. The motions of these residues undergo displacements which may facilitate solvent access to the ion.

Site-specific hydration dynamics in the nonpolar core of a molten globule by dynamic nuclear polarization of water

Water-protein interactions play a direct role in protein folding. The chain collapse that accompanies protein folding involves extrusion of water from the nonpolar core. For many proteins, including apomyoglobin (apoMb), hydrophobic interactions drive an initial collapse to an intermediate state before folding to the final structure. However, the debate continues as to whether the core of the collapsed intermediate state is hydrated and, if so, what the dynamic nature of this water is. A key challenge is that protein hydration dynamics is significantly heterogeneous, yet suitable experimental techniques for measuring hydration dynamics with site-specificity are lacking. Here, we introduce Overhauser dynamic nuclear polarization at 0.35 T via site-specific nitroxide spin labels as a unique tool to probe internal and surface protein hydration dynamics with site-specific resolution in the molten globular, native, and unfolded protein states. The (1)H NMR signal enhancement of water carries information about the local dynamics of the solvent within ∼10 Å of a spin label. EPR is used synergistically to gain insights on local polarity and mobility of the spin-labeled protein. Several buried and solvent-exposed sites of apoMb are examined, each bearing a covalently bound nitroxide spin label. We find that the nonpoloar core of the apoMb molten globule is hydrated with water bearing significant translational dynamics, only 4-6-fold slower than that of bulk water. The hydration dynamics of the native state is heterogeneous, while the acid-unfolded state bears fast-diffusing hydration water. This study provides a high-resolution glimpse at the folding-dependent nature of protein hydration dynamics.

Solvent-induced backbone fluctuations and the collective librational dynamics of lysozyme studied by terahertz spectroscopy

THz spectroscopy is used to investigate the dynamics of the globular protein hen egg white lysozyme under varying hydration and temperature conditions. An analysis of the experimental spectra has revealed that the amount of solvent in the hydration shell has a strong influence on the low-frequency protein conformational dynamics and also the arrangement of hydrogen bonds in the protein secondary structure. Furthermore at a hydration level >0.2 we identify collective backbone fluctuations in the protein secondary structure that are not present at low hydration. It is possible that these solvent induced modes are important for the biological function of the protein.

Mean-squared atomic displacements in hydrated lysozyme, native and denatured

We use elastic neutron scattering to demonstrate that a sharp increase in the mean-squared atomic displacements, commonly observed in hydrated proteins above 200 K and often referred to as the dynamical transition, is present in the hydrated state of both native and denatured lysozyme. A direct comparison of the native and denatured protein thus confirms that the presence of the transition in the mean-squared atomic displacements is not specific to biologically functional molecules.

Understanding dynamics of myoglobin in heterogeneous aqueous environments using coupled water fractions

This work presents an analysis of near environment of myoglobin (Mb) in different aqueous solutions (in the presence of NaCl, sucrose, trehalose, urea, and glycerol) using the coupled water fractions measured using a quartz crystal microbalance (QCM). The secondary structural features of the protein from circular dichroic (CD) spectroscopy and the coupled water fractions give important clues to the overall dynamics of the protein. Using time resolved fluorescence, these leads have been applied to understand the observed lifetime relaxations of Mb. Though the time scales of observation of coupled water and the lifetimes are very different, our study suggests that the trends in coupled water fraction seem to be good indicators for regulation of the relaxation dynamics of the protein. The relaxations generally show a triphasic distribution of time scales. The initial relaxation in the picoseconds time scale represents the local motions of coupled water followed by a slightly slower decay in hundreds of picoseconds attributable to coupled water-'quasi free' water interactions. The third nanosecond lifetime is due to changes in transitions in isomers of hydrated protein. The dynamics of coupled water in Mb with NaCl is the fastest (around 21 ps) and is slowest in glycerol (250 ps). The results strongly indicate that it is the resident times of water molecules that play a dominant role in the overall stability of protein in a particular hydrated isomer and not just always the number of such water molecules in the hydrated protein.