Evidence of protein collective motions on the picosecond timescale

We investigate the presence of structural collective motions on a picosecond timescale for the heme protein, cytochrome c, as a function of oxidation and hydration, using terahertz (THz) time domain spectroscopy and molecular dynamics simulations. The THz response dramatically increases with oxidation, with the largest increase for lowest hydrations, and highest frequencies. For both oxidation states the THz response rapidly increases with hydration saturating above ∼25% (g H₂O/g protein). Quasiharmonic vibrational modes and dipole-dipole correlation functions were calculated from molecular dynamics trajectories. The collective mode density of states alone reproduces the measured hydration dependence, providing strong evidence of the existence of these motions. The large oxidation dependence is reproduced only by the dipole-dipole correlation function, indicating the contrast arises from diffusive motions consistent with structural changes occurring in the vicinity of buried internal water molecules. This source for the observed oxidation dependence is consistent with the lack of an oxidation dependence in nuclear resonant vibrational spectroscopy measurements.     

Dielectric Relaxation Dynamics of Water in Model Membranes Probed by Terahertz Spectroscopy

We study hydrated model membranes, consisting of stacked bilayers of 1,2-dioleoyl-sn-glycero-3-phosphocholine lipids, using terahertz time-domain spectroscopy and infrared spectroscopy. Terahertz spectroscopy enables the investigation of water dynamics, owing to its sensitivity to dielectric relaxation processes associated with water reorientation. By controlling the number of water molecules per lipid molecule in the system, we elucidate how the interplay between the model membrane and water molecules results in different water dynamics. For decreasing hydration levels, we observe the appearance of new types of water dynamics: the collective bulklike dynamics become less pronounced, whereas an increased amount of both very slowly reorienting (i.e., irrotational) and very rapidly reorienting (i.e., fast) water molecules appear. Temperature-dependent measurements reveal the interconversion between the three distinct types of water present in the system.     

Hydration dynamics of hyaluronan and dextran

Hyaluronan is a polysaccharide, which is ubiquitous in vertebrates and has been reported to be strongly hydrated in a biological environment. We study the hydration of hyaluronan in solution using the rotational dynamics of water as a probe. We measure these dynamics with polarization-resolved femtosecond-infrared and terahertz time-domain spectroscopies. Both experiments reveal that a subensemble of water molecules is slowed down in aqueous solutions of hyaluronan amounting to ～15 water molecules per disaccharide unit. This quantity is consistent with what would be expected for the first hydration shell. Comparison of these results to the water dynamics in aqueous dextran solution, a structurally similar polysaccharide, yields remarkably similar results. This suggests that the observed interaction with water is a common feature for hydrophilic polysaccharides and is not specific to hyaluronan.     

Hydration affects both harmonic and anharmonic nature of protein dynamics

To understand the effect of hydration on protein dynamics, inelastic neutron-scattering experiments were performed on staphylococcal nuclease samples at differing hydration levels: dehydrated, partially hydrated, and hydrated. At cryogenic temperatures, hydration affected the collective motions with energies lower than 5 meV, whereas the high-energy localized motions were independent of hydration. The prominent change was a shift of boson peak toward higher energy by hydration, suggesting a hardening of harmonic potential at local minima on the energy landscape. The 240 K transition was observed only for the hydrated protein. Significant quasielastic scattering at 300 K was observed only for the hydrated sample, indicating that the origin of the transition is the motion activated by hydration water. The neutron-scattering profile of the partially hydrated sample was quite similar to that of the hydrated sample at 100 K and 200 K, whereas it was close to the dehydrated sample at 300 K, indicating that partial hydration is sufficient to affect the harmonic nature of protein dynamics, and that there is a threshold hydration level to activate anharmonic motions. Thus, hydration water controls both harmonic and anharmonic protein dynamics by differing means.     

Dynamics of tRNA at Different Levels of Hydration

The influence of hydration on the nanosecond timescale dynamics of tRNA is investigated using neutron scattering spectroscopy. Unlike protein dynamics, the dynamics of tRNA is not affected by methyl group rotation. This allows for a simpler analysis of the influence of hydration on the conformational motions in RNA. We find that hydration affects the dynamics of tRNA significantly more than that of lysozyme. Both the characteristic length scale and the timescale of the conformational motions in tRNA depend strongly on hydration. Even the characteristic temperature of the so-called “dynamical transition” appears to be hydration-dependent in tRNA. The amplitude of the conformational motions in fully hydrated tRNA is almost twice as large as in hydrated lysozyme. We ascribe these differences to a more open and flexible structure of hydrated RNA, and to a larger fraction and different nature of hydrophilic sites. The latter leads to a higher density of water that makes the biomolecule more flexible. All-atom molecular-dynamics simulations are used to show that the extent of hydration is greater in tRNA than in lysozyme. We propose that water acts as a “lubricant” in facilitating enhanced motion in solvated RNA molecules.     

Do hydration dynamics follow the structural perturbation during thermal denaturation of a protein: a terahertz absorption study

We investigate the thermal denaturation of human serum albumin and the associated solvation using terahertz (THz) spectroscopy in aqueous buffer solution. Far- and near-ultraviolet circular dichroism spectroscopy reveal that the protein undergoes a native (N) to extended (E) state transition at temperature ≤55°C with a marginal change in the secondary and tertiary structure. At 70°C, the protein transforms into an unfolded (U) state with significant irreversible disruption of its structures. We measure the concentration- and temperature-dependent THz absorption coefficient (α) of the protein solution using a p-Ge THz difference spectrometer (2.1–2.8 THz frequency range), thereby probing the collective protein-water network dynamics. When the solvated protein is heated up to 55°C and cooled down again, a reversible change in THz absorption is observed. When increasing the temperature up to 70°C, we find a dramatic irreversible change of THz absorption. The increase in THz absorption compared to bulk water is attributed to a blue shift in the spectrum of the solvated protein compared to bulk water. This is supported by measurements of THz absorption coefficients using THz time-domain spectroscopy (0.1–1.2 THz frequency range). We also use picosecond-resolved fluorescence spectroscopy of the tryptophan 214 moiety of human serum albumin. All experimental observations can be explained by a change in the hydration dynamics of the solvated protein due to the additional exposure of hydrophobic residues upon unfolding.     

Site-Specific Peptide Probes Detect Buried Water in a Lipid Membrane

Transmembrane peptides contain polar residues in the interior of the membrane, which may alter the electrostatic environment and favor hydration in the otherwise nonpolar environment of the membrane core. Here, we demonstrate a general, nonperturbative strategy to probe hydration of the peptide backbone at specific depths within the bilayer using a combination of site-specific isotope labels, ultrafast two-dimensional infrared spectroscopy, and spectral modeling based on molecular dynamics simulations. Our results show that the amphiphilic pH-low insertion peptide supports a highly heterogeneous environment, with significant backbone hydration of nonpolar residues neighboring charged residues. For example, a leucine residue located as far as 1 nm into the hydrophobic bulk reports hydrogen-bonded populations as high as ∼20%. These findings indicate that the polar nature of these residues may facilitate the transport of water molecules into the hydrophobic core of the membrane.     

Hydration of short DNA, RNA and 2'-OMe oligonucleotides determined by osmotic stressing

Studies on hydration are important for better understanding of structure and function of nucleic acids. We compared the hydration of self-complementary DNA, RNA and 2′-O-methyl (2′-OMe) oligonucleotides GCGAAUUCGC, (UA)₆ and (CG)₃ using the osmotic stressing method. The number of water molecules released upon melting of oligonucleotide duplexes, Δn_W, was calculated from the dependence of melting temperature on water activity and the enthalpy, both measured with UV thermal melting experiments. The water activity was changed by addition of ethylene glycol, glycerol and acetamide as small organic co-solutes. The Δn_W was 3–4 for RNA duplexes and 2–3 for DNA and 2′-OMe duplexes. Thus, the RNA duplexes were hydrated more than the DNA and the 2′-OMe oligonucleotide duplexes by approximately one to two water molecules depending on the sequence. Consistent with previous studies, GC base pairs were hydrated more than AU pairs in RNA, whereas in DNA and 2′-OMe oligonucleotides the difference in hydration between these two base pairs was relatively small. Our data suggest that the better hydration of RNA contributes to the increased enthalpic stability of RNA duplexes compared with DNA duplexes.     

Concentration-dependent effects on fully hydrated DNA at terahertz frequencies

Using terahertz time-domain spectroscopy (THz-TDS), the frequency-dependent dielectric constant of deoxyribonucleic acid (DNA) in solution was measured. The response of the buffer solution is dominated by two Debye modes in this frequency range, and, from an analysis of the concentration dependence, the presence of the DNA increases the main relaxation time and dielectric constant. This reflects the fact that the water in the hydration layer is more tightly bound under the influence of the DNA molecule in comparison to bulk water. This dynamical slowing down with increasing DNA concentration is similar to what is observed with purine nucleotides, but opposite to the behavior of pyrimidine nucleotides. In addition, a suspension model was used with the concentration-dependent data to isolate the dielectric response of the hydrated DNA molecule. The data for the hydrated DNA molecule is still dominated by a Debye response. It is also possible to determine the thickness of the hydration layer, and the DNA molecule influences the surrounding water out to 16 or 17 Å, which corresponds to about six effective hydration layers.     

Dynamics of Biological Macromolecules: Not a Simple Slaving by Hydration Water

We studied the dynamics of hydrated tRNA using neutron and dielectric spectroscopy techniques. A comparison of our results with earlier data reveals that the dynamics of hydrated tRNA is slower and varies more strongly with temperature than the dynamics of hydrated proteins. At the same time, tRNA appears to have faster dynamics than DNA. We demonstrate that a similar difference appears in the dynamics of hydration water for these biomolecules. The results and analysis contradict the traditional view of slaved dynamics, which assumes that the dynamics of biological macromolecules just follows the dynamics of hydration water. Our results demonstrate that the dynamics of biological macromolecules and their hydration water depends strongly on the chemical and three-dimensional structures of the biomolecules. We conclude that the whole concept of slaving dynamics should be reconsidered, and that the mutual influence of biomolecules and their hydration water must be taken into account.     

Collective vibrational modes in biological molecules investigated by terahertz time-domain spectroscopy

We present well-resolved absorption spectra of biological molecules in the far-IR (FIR) spectral region recorded by terahertz time-domain spectroscopy (THz-TDS). As an illustrative example we discuss the absorption spectra of benzoic acid, its monosubstitutes salicylic acid (2-hydroxy-benzoic acid), 3- and 4-hydroxybenzoic acid, and aspirin (acetylsalicylic acid) in the spectral region between 18 and 150 cm(-1). The spectra exhibit distinct features originating from low-frequency vibrational modes caused by intra- or intermolecular collective motion and lattice modes. Due to the collective origin of the observed modes the absorption spectra are highly sensitive to the overall structure and configuration of the molecules, as well as their environment. The THz-TDS procedure can provide a direct fingerprint of the molecular structure or conformational state of a compound.     

Probing the collective vibrational dynamics of a protein in liquid water by terahertz absorption spectroscopy

Biological polymers are expected to exhibit functionally relevant, global, and subglobal collective modes in the terahertz (THz) frequency range (i.e., picosecond timescale). In an effort to monitor these collective motions, we have experimentally determined the absorption spectrum of solvated bovine serum albumin (BSA) from 0.3 to 3.72 THz (10-124 cm(-1)). We successfully extract the terahertz molar absorption of the solvated BSA from the much stronger attenuation of water and observe in the solvated protein a dense, overlapping spectrum of vibrational modes that increases monotonically with increasing frequency. We see no evidence of distinct, strong, spectral features, suggesting that no specific collective vibrations dominate the protein's spectrum of motions, consistent with the predictions of molecular dynamics simulations and normal mode analyses of a range of small proteins. The shape of the observed spectrum resembles the ideal quadratic spectral density expected for a disordered ionic solid, indicating that the terahertz normal mode density of the solvated BSA may be modeled, to first order, as that of a three-dimensional elastic nanoparticle with an aperiodic charge distribution. Nevertheless, there are important detailed departures from that of a disordered inorganic solid or the normal mode densities predicted for several smaller proteins. These departures are presumably the spectral features arising from the unique molecular details of the solvated BSA. The techniques used here and measurements have the potential to experimentally confront theoretical calculations on a frequency scale that is important for macromolecular motions in a biologically relevant water environment.     

Low-frequency vibrational modes and infrared absorbance of red, blue and green opsin

Vibrational excitations of low-frequency collective modes are essential for functionally important conformational transitions in proteins. We carried out an analysis of the low-frequency modes in the G protein coupled receptors (GPCR) family of cone opsins based on both normal-mode analysis and molecular dynamics (MD) simulations. Power spectra obtained by MD can be compared directly with normal modes. In agreement with existing experimental evidence related to transmembrane proteins, cone opsins have functionally important transitions that correspond to approximately 950 modes and are found below 80 cm⁻¹. This is in contrast to bacteriorhodopsin and rhodopsin, where the important low-frequency transition modes are below 50 cm⁻¹. We find that the density of states (DOS) profile of blue opsin in a solvent (e.g. water) has increased populations in the very lowest frequency modes (<15 cm⁻¹); this is indicative of the increased thermostability of blue opsin. From our work we found that, although light absorption behaves differently in blue, green and red opsins, their low-frequency vibrational motions are similar. The similarities and differences in the domain motions of blue, red and green opsins are discussed for several representative modes. In addition, the influence of the presence of a solvent is reported and compared with vacuum spectra. We thus demonstrate that terahertz spectroscopy of low-frequency modes might be relevant for identifying those vibrational degrees of freedom that correlate to known conformational changes in opsins. An erratum to this article can be found at     

Proteins maintain hydration at high [KCl] concentration regardless of content in acidic amino acids

Proteins of halophilic organisms, which accumulate molar concentrations of KCl in their cytoplasm, have a much higher content in acidic amino acids than proteins of mesophilic organisms. It has been proposed that this excess is necessary to maintain proteins hydrated in an environment with low water activity, either via direct interactions between water and the carboxylate groups of acidic amino acids or via cooperative interactions between acidic amino acids and hydrated cations. Our simulation study of five halophilic proteins and five mesophilic counterparts does not support either possibility. The simulations use the AMBER ff14SB force field with newly optimized Lennard-Jones parameters for the interactions between carboxylate groups and potassium ions. We find that proteins with a larger fraction of acidic amino acids indeed have higher hydration levels, as measured by the concentration of water in their hydration shell and the number of water/protein hydrogen bonds. However, the hydration level of each protein is identical at low (b_KCl = 0.15 mol/kg) and high (b_KCl = 2 mol/kg) KCl concentrations; excess acidic amino acids are clearly not necessary to maintain proteins hydrated at high salt concentration. It has also been proposed that cooperative interactions between acidic amino acids in halophilic proteins and hydrated cations stabilize the folded protein structure and would lead to slower dynamics of the solvation shell. We find that the translational dynamics of the solvation shell is barely distinguishable between halophilic and mesophilic proteins; if such a cooperative effect exists, it does not have that entropic signature.     

Dynamics of uncrystallized water and protein in hydrated elastin studied by thermal and dielectric techniques

Dynamics of uncrystallized water and protein was studied in hydrated pellets of the fibrous protein elastin in a wide hydration range (0 to 23 wt.%), by differential scanning calorimetry (DSC), thermally stimulated depolarization current technique (TSDC) and dielectric relaxation spectroscopy (DRS). Additionally, water equilibrium sorption–desorption measurements (ESI) were performed at room temperature. The glass transition of the system was studied by DSC and its complex dependence on hydration water was verified. A critical water fraction of about 18 wt.% was found, associated with a reorganization of water in the material. Three dielectric relaxations, associated to dynamics related to distinct uncrystallized water populations, were recorded by TSDC and DRS. The low temperature secondary relaxation of hydrophilic polar groups on the protein surface triggered by hydration water for almost dry samples contains contributions from water molecules themselves at higher water fractions (ν relaxation). This particular relaxation is attributed to water molecules in the primary and secondary hydration shells of the protein fibers. At higher temperatures and for water fraction values equal to or higher than 10 wt.%, a local relaxation of water molecules condensed within small openings in the interior of the protein fibers was recorded. The evolution of this relaxation (w relaxation) with hydration level results in enhanced cooperativity at high water fraction values, implying the existence of “internal” water confined within the protein structure. At higher temperatures a relaxation associated with water dynamics within clusters between fibers (p relaxation) was also recorded, in the same hydration range.     

Thermal fluctuations of DNA enclosed by glycerol–water glassy matrices: an elastic neutron scattering investigation

Through elastic neutron scattering measurements, we investigated the thermal fluctuations of DNA enclosed by glycerol–water glassy matrices, at different levels of hydration, over the wide temperature range from 20 to 300 K. For all the samples, the extracted hydrogen mean square displacements (MSD) show a purely vibrational harmonic trend at very low temperatures, and a first onset of anharmonic dynamics above ∼100 K. Such onset is consistent with the activation of DNA methyl group rotational motions. Then, at a certain temperature T _d, the MSD show a second onset of anharmonicity, which corresponds to the DNA dynamical transition. The T _d values vary as a function of the hydration degree of the environment. The crucial role of the solvent mobility to activate the DNA thermal fluctuations is proposed, together with a preferential hydration effect of the DNA phosphate groups. Finally, a comparison between the average mobility of homologous samples of DNA and the lysozyme protein is considered. Advanced neutron scattering and complementary techniques to study biological systems. Contributions from the meetings, “Neutrons in Biology”, STFC Rutherford Appleton Laboratory, Didcot, UK, 11–13 July and “Proteins At Work 2007”, Perugia, Italy, 28–30 May 2007.     

Effect of conformational states on protein dynamical transition

In order to examine the properties specific to the folded protein, the effect of the conformational states on protein dynamical transition was studied by incoherent elastic neutron scattering for both wild type and a deletion mutant of staphylococcal nuclease. The deletion mutant of SNase which lacks C-terminal 13 residues takes a compact denatured structure, and can be regarded as a model of intrinsic unstructured protein. Incoherent elastic neutron scattering experiments were carried out at various temperature between 10 K and 300 K on IN10 and IN13 installed at ILL. Temperature dependence of mean-square displacements was obtained by the q-dependence of elastic scattering intensity. The measurements were performed on dried and hydrated powder samples. No significant differences were observed between wild type and the mutant for the hydrated samples, while significant differences were observed for the dried samples. A dynamical transition at ∼ 140 K observed for both dried and hydrated samples. The slopes of the temperature dependence of MSD before transition and after transition are different between wild type and the mutant, indicating the folding induces hardening. The hydration water activates a further transition at ∼ 240 K. The behavior of the temperature dependence of MSD is indistinguishable for wild type and the mutant, indicating that hydration water dynamics dominate the dynamical properties.     

Incoherent neutron scattering of copper azurin: a comparison with molecular dynamics simulation results

The low-frequency dynamics of copper azurin has been studied at different temperatures for a dry and deuterium hydrated sample by incoherent neutron scattering and the experimental results have been compared with molecular dynamics (MD) simulations carried out in the same temperature range. Experimental Debye-Waller factors are consistent with a dynamical transition at approximately 200 K which appears partially suppressed in the dry sample. Inelastic and quasielastic scattering indicate that hydration water modulates both vibrational and diffusive motions. The low-temperature experimental dynamical structure factor of the hydrated protein shows an excess of inelastic scattering peaking at about 3 meV and whose position is slightly shifted downwards in the dry sample. Such an excess is reminiscent of the “boson peak” observed in glass-like materials. This vibrational peak is quite well reproduced by MD simulations, although at a lower energy. The experimental quasielastic scattering of the two samples at 300 K shows a two-step relaxation behaviour with similar characteristic times, while the corresponding intensities differ only by a scale factor. Also, MD simulations confirm the two-step diffusive trend, but the slow process seems to be characterized by a decay faster than the experimental one. Comparison with incoherent neutron scattering studies carried out on proteins having different structure indicates that globular proteins display common elastic, quasielastic and inelastic features, with an almost similar hydration dependence, irrespective of their secondary and tertiary structure. Received: 12 October 1998 / Revised version: 19 February 1999 / Accepted: 1 March 1999     

A Raman scattering study of the interactions of DNA with its water of hydration

Raman spectroscopy is used to probe the nature of the hydrogen bonds which hold the water of hydration to DNA. The ～ 3450?cm^?1 molecular O–H stretching mode shows that the first six water molecules per base pair of the primary hydration shell are very strongly bound to the DNA. The observed shift in the peak position of this mode permits a determination of the length of the hydrogen bonds for these water molecules. These hydrogen bonds appear to be about 0.3?Å shorter than the hydrogen bonds in bulk water. The linewidth of this mode shows no significant changes above water contents of about 15 water molecules per base pair. This technique of using a vibrational spectroscopy to obtain structural information about the hydration shells of DNA could be used to study the hydration shells of other biomolecules.