Molecular view of water dynamics near model peptides

A peptide's perspective of water dynamics

This Perspective is focused on amide groups of peptides interacting with water. The 2D IR spectroscopy has already enabled structural aspects of the peptide backbone to be determined through its ability to measure the coupling between different amide-I modes. Here we describe why nonlinear IR is emerging as the method of choice to examine the fast components of the water dynamics near peptides and how isotopically edited peptide links can be used to probe the local water at a residue level in proteins. This type of research necessarily involves an intimate mix of theory and experiment. The description of the results is underpinned by relatively well established quantum-statistical theories that describe the important manifestations of peptide vibrational frequency fluctuations.

Vibrational spectroscopy of water in hydrated lipid multi-bilayers. I. Infrared spectra and ultrafast pump-probe observables

The vibrational spectroscopy of hydration water in dilauroylphosphatidylcholine lipid multi-bilayers is investigated using molecular dynamics simulations and a mixed quantum/classical model for the OD stretch spectroscopy of dilute HDO in H(2)O. FTIR absorption spectra, and isotropic and anisotropic pump-probe decay curves have been measured experimentally as a function of the hydration level of the lipid multi-bilayer, and our goal is to make connection with these experiments. To this end, we use third-order response functions, which allow us to include non-Gaussian frequency fluctuations, non-Condon effects, molecular rotations, and a fluctuating vibrational lifetime, all of which we believe are important for this system. We calculate the response functions using existing transition frequency and dipole maps. From the experiments it appears that there are two distinct vibrational lifetimes corresponding to HDO molecules in different molecular environments. In order to obtain these lifetimes, we consider a simple two-population model for hydration water hydrogen bonds. Assuming a different lifetime for each population, we then calculate the isotropic pump-probe decay, fitting to experiment to obtain the two lifetimes for each hydration level. With these lifetimes in hand, we then calculate FTIR spectra and pump-probe anisotropy decay as a function of hydration. This approach, therefore, permits a consistent calculation of all observables within a unified computational scheme. Our theoretical results are all in qualitative agreement with experiment. The vibrational lifetime of lipid-associated OD groups is found to be systematically shorter than that of the water-associated population, and the lifetimes of each population increase with decreasing hydration, in agreement with previous analysis. Our theoretical FTIR absorption spectra successfully reproduce the experimentally observed red-shift with decreasing lipid hydration, and we confirm a previous interpretation that this shift results from the hydrogen bonding of water to the lipid phosphate group. From the pump-probe anisotropy decay, we confirm that the reorientational motions of water molecules slow significantly as hydration decreases, with water bound in the lipid carbonyl region undergoing the slowest rotations.

Molecular evolution of rbcL in three gymnosperm families: identifying adaptive and coevolutionary patterns

Background

The chloroplast-localized ribulose-1, 5-biphosphate carboxylase/oxygenase (Rubisco), the primary enzyme responsible for autotrophy, is instrumental in the continual adaptation of plants to variations in the concentrations of CO₂. The large subunit (LSU) of Rubisco is encoded by the chloroplast rbcL gene. Although adaptive processes have been previously identified at this gene, characterizing the relationships between the mutational dynamics at the protein level may yield clues on the biological meaning of such adaptive processes. The role of such coevolutionary dynamics in the continual fine-tuning of RbcL remains obscure.

Results

We used the timescale and phylogenetic analyses to investigate and search for processes of adaptive evolution in rbcL gene in three gymnosperm families, namely Podocarpaceae, Taxaceae and Cephalotaxaceae. To understand the relationships between regions identified as having evolved under adaptive evolution, we performed coevolutionary analyses using the software CAPS. Importantly, adaptive processes were identified at amino acid sites located on the contact regions among the Rubisco subunits and on the interface between Rubisco and its activase. Adaptive amino acid replacements at these regions may have optimized the holoenzyme activity. This hypothesis was pinpointed by evidence originated from our analysis of coevolution that supported the correlated evolution between Rubisco and its activase. Interestingly, the correlated adaptive processes between both these proteins have paralleled the geological variation history of the concentration of atmospheric CO₂.

Conclusions

The gene rbcL has experienced bursts of adaptations in response to the changing concentration of CO₂ in the atmosphere. These adaptations have emerged as a result of a continuous dynamic of mutations, many of which may have involved innovation of functional Rubisco features. Analysis of the protein structure and the functional implications of such mutations put forward the conclusion that this evolutionary scenario has been possible through a complex interplay between adaptive mutations, often structurally destabilizing, and compensatory mutations. Our results unearth patterns of evolution that have likely optimized the Rubisco activity and uncover mutational dynamics useful in the molecular engineering of enzymatic activities.

Reviewers

This article was reviewed by Prof. Christian Blouin (nominated by Dr W Ford Doolittle), Dr Endre Barta (nominated by Dr Sandor Pongor), and Dr Nicolas Galtier.

Probing deuterium isotope effect on structure and solvation dynamics of human serum albumin

The deuterium isotopic effect on the structure and solvation dynamics of the protein, human serum albumin (HSA), has been studied by using circular dichroism (CD), femtosecond up-conversion, FRET, and single-molecule spectroscopy. The CD spectra suggest that D(2)O affects the structure of HSA, leading to a 20% decrease in the helical structure. The FRET study indicates that the distance of C153 from the lone tryptophan residue of HSA is quite similar (≈21 Å) in H(2)O and D(2)O, and hence, the location of the probe in the protein remains the same in the two solvents. The single-molecule study suggests that coumarin 153 (C153) binds almost exclusively (>96%) to one site of HSA. Solvation dynamics of C153 in HSA is found to be markedly retarded in D(2)O compared with H(2)O. In H(2)O, the solvation of C153 bound to HSA is found to be biexponential with one component of 7 ps (30%) and a long component of 350 ps (70%). In D(2)O, we detected a short component of 4 ps (41%) and a long component of 950 ps (59%). Thus, the ultraslow component of the solvation dynamics of C153 bound to HSA in D(2)O (950 ps) is 2.5-fold slower than that in H(2)O (350 ps). The marked deuterium isotope effect has been ascribed to water molecules confined in the protein environment and to a lesser extent to the structural modification of protein by D(2)O.

Hydration dynamics at fluorinated protein surfaces

Water-protein interactions dictate many processes crucial to protein function including folding, dynamics, interactions with other biomolecules, and enzymatic catalysis. Here we examine the effect of surface fluorination on water-protein interactions. Modification of designed coiled-coil proteins by incorporation of 5,5,5-trifluoroleucine or (4S)-2-amino-4-methylhexanoic acid enables systematic examination of the effects of side-chain volume and fluorination on solvation dynamics. Using ultrafast fluorescence spectroscopy, we find that fluorinated side chains exert electrostatic drag on neighboring water molecules, slowing water motion at the protein surface.

Evidence of functional protein dynamics from X-ray crystallographic ensembles

It is widely recognized that representing a protein as a single static conformation is inadequate to describe the dynamics essential to the performance of its biological function. We contrast the amino acid displacements below and above the protein dynamical transition temperature, T_D∼215K, of hen egg white lysozyme using X-ray crystallography ensembles that are analyzed by molecular dynamics simulations as a function of temperature. We show that measuring structural variations across an ensemble of X-ray derived models captures the activation of conformational states that are of functional importance just above T_D, and they remain virtually identical to structural motions measured at 300K. Our results highlight the ability to observe functional structural variations across an ensemble of X-ray crystallographic data, and that residue fluctuations measured in MD simulations at room temperature are in quantitative agreement with the experimental observable.

There is a well-recognized gap between the dynamical motions of proteins required to execute function and the experimental techniques capable of capturing that motion at the atomic level. We show that much experimental detail of dynamical motion is already present in X-ray crystallographic data, which arises from being solved by different research groups using different methodologies under different crystallization conditions, which then capture an ensemble of structures whose variations can be quantified on a residue-by-residue level using local density correlations. We contrast the amino acid displacements below and above the protein dynamical transition temperature, T_D∼215K, of hen egg white lysozyme by comparing the X-ray ensemble to MD ensembles as a function of temperature. We show that measuring structural variations across an ensemble of X-ray derived models captures the activation of conformational states that are of functional importance just above T_D and they remain virtually identical to structural motions measured at 300K. It provides a novel analysis of large X-ray ensemble data that is useful for the broader structural biology community.

Dynamics in Biological Systems as seen by QENS

Quasielastic incoherent neutron scattering is a well suited and established experimental method to study protein and water dynamics in the picosecond to nanosecond time- and Ångstrom length-scale. Using deuterium labelling either protein or water motions can be selected and brought into focus. Protein and cell water dynamics were separately studied in red blood cells. A consistent picture of cytoplasmic water and protein dynamics in whole cells is emerging from recent experimental results.

Protein hydration dynamics and molecular mechanism of coupled water-protein fluctuations

Protein surface hydration is fundamental to its structural stability and flexibility, and water-protein fluctuations are essential to biological function. Here, we report a systematic global mapping of water motions in the hydration layer around a model protein of apomyoglobin in both native and molten globule states. With site-directed mutagenesis, we use intrinsic tryptophan as a local optical probe to scan the protein surface one at a time with single-site specificity. With femtosecond resolution, we examined 16 mutants in two states and observed two types of water-network relaxation with distinct energy and time distributions. The first water motion results from the local collective hydrogen-bond network relaxation and occurs in a few picoseconds. The initial hindered motions, observed in bulk water in femtoseconds, are highly suppressed and drastically slow down due to structured water-network collectivity in the layer. The second water-network relaxation unambiguously results from the lateral cooperative rearrangements in the inner hydration shell and occurs in tens to hundreds of picoseconds. Significantly, this longtime dynamics is the coupled interfacial water-protein motions and is the direct measurement of such cooperative fluctuations. These local protein motions, although highly constrained, are necessary to assist the longtime water-network relaxation. A series of correlations of hydrating water dynamics and coupled fluctuations with local protein's chemical and structural properties were observed. These results are significant and reveal various water behaviors in the hydration layer with wide heterogeneity. We defined a solvation speed and an angular speed to quantify the water-network rigidity and local protein flexibility, respectively. We also observed that the dynamic hydration layer extends to more than 10 A. Finally, from native to molten globule states, the hydration water networks loosen up, and the protein locally becomes more flexible with larger global plasticity and partial unfolding.

Effects of co-solvents on peptide hydration water structure and dynamics

We evaluate the molecular response of hydration water as a function of temperature and proximity to the surface of the peptide N-acetyl-leucine-methylamide (NALMA) when in the presence of the kosmotrope co-solvent glycerol or the chaotrope co-solvent dimethyl sulfoxide (DMSO), using molecular dynamics simulation with a polarizable force field. These detailed microscopic studies complement established thermodynamic analysis on the role of co-solvents in shifting the equilibrium for proteins away from or towards the native folded state. We find that the structure of the water at the peptide interfaces reflects an increase in hydration number in the glycerol solution and a decrease in hydration numbers in the DMSO solution. While the water dynamics around NALMA in the presence of both co-solvents is slower than that observed with the water solvent alone, in the DMSO mixture we no longer measure a separation in water motion time scales at low temperatures as is seen in the pure water solvent, but rather one single relaxation time. In the glycerol, however, we do observe a separation of time scales at low temperatures, supporting the hypothesis that hydration water near a hydrophobic solute evolves on a separate time scale than the extensive hydrogen-bonding network of more bulk-like water. Our simulation studies highlight the differences in the two co-solvent solutions due to the relative frequency of water contacts with the hydrophobic vs. hydrophilic peptide surface, and direct water interactions with the co-solvents.

Hydration water dynamics near biological interfaces

We performed classical molecular dynamics simulations using both fixed-charge and polarizable water and protein force fields to contrast the hydration dynamics near hydrophilic and amphiphilic peptides as a function of temperature. The high peptide concentrations we use serve as a model for the surface of folded proteins where hydration layers around each residue overlap significantly. Through simulation we determine that there are notable differences in the water dynamics analyzed from the outer and inner hydration layer regions of the amphiphilic peptide solution that explains the experimentally observed presence of two translational relaxations, while the hydrophilic peptide solution shows only a single non-Arrhenius translational process with no distinction between hydration layers. Given that water dynamics for the amphiphilic peptide system reproduces all known rotational and translational hydration dynamical anomalies exhibited by hydration water near protein surfaces, our analysis provides strong evidence that dynamical signatures near biological interfaces arises because of frustration in the hydration dynamics induced by chemical heterogeneity, as opposed to just topological roughness, of the protein surface.

Time scales of water dynamics at biological interfaces: peptides, proteins and cells

Water 2H and 17O spin relaxation is used to study water dynamics in the hydration layers of two small peptides, two globular proteins and in living cells of two microorganisms. The dynamical heterogeneity of hydration water is characterized by performing relaxation measurements over a wide temperature range, extending deeply into the supercooled regime, or by covering a wide frequency range. Protein hydration layers can be described by a power-law distribution of rotational correlation times with an exponent close to 2. This distribution comprises a small fraction of protein-specific hydration sites, where water rotation is strongly retarded, and a dominant fraction of generic hydration sites, where water rotation is as fast as in the hydration shells of small peptides. The generic dynamic perturbation factor is less than 2 at room temperature and exhibits a maximum near 260 K. The dynamic perturbation is induced by H-bond constraints that interfere with the cooperative mechanism that facilitates rotation in bulk water. Because these constraints are temperature-independent, hydration water does not follow the super-Arrhenius temperature dependence of bulk water. Water in living cells behaves as expected from studies of simpler model systems, the only difference being a larger fraction of secluded (strongly perturbed) hydration sites associated with the supramolecular organization in the cell. Intracellular water that is not in direct contact with biopolymers has essentially the same dynamics as bulk water. There is no significant difference in cell water dynamics between mesophilic and halophilic organisms, despite the high K+ and Na+ concentrations in the latter.

Cytoplasmic water and hydration layer dynamics in human red blood cells

The dynamics of water in human red blood cells was measured with quasielastic incoherent neutron scattering in the temperature range between 290 and 320 K. Neutron spectrometers with time resolutions of 40, 13, and 7 ps were combined to cover time scales of bulk water dynamics to reduced mobility interfacial water motions. A major fraction of approximately 90% of cell water is characterized by a translational diffusion coefficient similar to bulk water. A minor fraction of approximately 10% of cellular water exhibits reduced dynamics. This slow water fraction was attributed to dynamically bound water on the surface of hemoglobin which accounts for approximately half of the hydration layer.

Hydration profiles of amyloidogenic molecular structures

Hydration shells of normal proteins display regions of highly structured water as well as patches of less structured bulk-like water. Recent studies suggest that isomers with larger surface densities of patches of bulk-like water have an increased propensity to aggregate. These aggregates are toxic to the cellular environment. Hence, the early detection of these toxic deposits is of paramount medical importance. We show that various morphological states of association of such isomers can be differentiated from the normal protein background based on the characteristic partition between bulk, caged, and surface hydration water and the magnetic resonance (MR) signals of this water. We derive simple mathematical equations relating the compartmentalization of water to the local hydration fraction and the packing density of the newly formed molecular assemblies. Then, we employ these equations to predict the MR response of water constrained by protein aggregation. Our results indicate that single units and compact aggregates that contain no water between constituents induce a shift of the MR signal from normal protein background to values in the hyperintensity domain (bright spots), corresponding to bulk water. In contrast, large plaques that cage significant amounts of water between constituents are likely to generate MR responses in the hypointensity domain (dark spots), typical for strongly correlated water. The implication of these results is that amyloids can display both dark and bright spots when compared to the normal gray background tissue on MR images. In addition, our findings predict that the bright spots are more likely to correspond to amyloids in their early stage of development. The results help explain the MR contrast patterns of amyloids and suggest a new approach for identifying unusual protein aggregation related to disease.

Water hydrogen bond analysis on hydrophilic and hydrophobic biomolecule sites

Elastic and quasielastic neutron scattering experiments have been used to investigate the hydrogen bonding network dynamics of hydration water on hydrophilic and hydrophobic sites. To this end the evolution of hydration water dynamics of a prototypical hydrophobic amino acid with polar backbone, N-acetyl-leucine-methylamide (NALMA), and hydrophilic amino acid, N-acetyl-glycine-methylamide (NAGMA), has been investigated as a function of the molecular ratio water : peptide. The results suggest that the dynamical contribution of the intrinsic and low hydration molecules of water is characteristic of pure librational/rotational movement. The water molecule remains attached to the hydrophilic site with only the possibility of hindered rotations that eventually break the bond with the peptide and reform it immediately after. A gradual evolution from librational motions to hindered rotations is observed as a function of temperature. When the hydration increases, we observe (together with the hindered rotations of hydrogen bonds) a slow diffusion of water molecules on the surface of the peptides.

Thickness of the hydration layer of a protein from molecular dynamics simulation

Water molecules around a protein exhibit slow dynamics with respect to that of pure bulk water. One important issue in protein hydration is the thickness of the hydration layer (i.e., the distance from the protein surface up to which the water dynamics is influenced by the protein). Estimation of thickness is crucial to understand better the properties of "biological water" and the role that it plays in guiding the protein's function. We have performed an atomistic molecular dynamics simulation of an aqueous solution of the protein villin headpiece subdomain or HP-36 to estimate the thickness of its hydration water. In particular, several dynamical properties of water around different segments (three alpha-helices) of the protein have been calculated by varying the thickness of the hydration layers. It is found that in general the influence of the helices on water properties extends beyond the first hydration layer. However, the heterogeneous nature of water among the first hydration layers of the three helices diminishes as the thickness is increased. It indicates that, for a small protein such as HP-36, the thickness of "biological water" is uniform for different segments of the protein.

Dynamics of water in the hydration layer of a partially unfolded structure of the protein HP-36

Atomistic molecular dynamics simulations of the folded native structure and a partially unfolded molten globule structure of the protein villin headpiece subdomain or HP-36 have been carried out with explicit solvent to explore the effects of unfolding on the dynamical behavior of water present in the hydration layers of different segments (three alpha-helices) of the protein. The calculations revealed that the unfolding of helix-2 influences the translational and rotational motions of water present in the hydration layers of the three helices in a heterogeneous manner. It is observed that a correlation exists between the unfolding of helix-2 and the microscopic kinetics of protein-water hydrogen bonds formed by its residues. This in turn has an influence on the rigidity of the hydration layers of the helices in the unfolded structure versus that in the folded native structure. These results should provide a microscopic explanation to recent solvation dynamics experiments on folded native and unfolded structures of proteins.

How protein surfaces induce anomalous dynamics of hydration water

Water around biomolecules slows down with respect to pure water, and both rotation and translation exhibit anomalous time dependence in the hydration shell. The origin of such behavior remains elusive. We use molecular dynamics simulations of water dynamics around several designed protein models to establish the connection between the appearance of the anomalous dynamics and water-protein interactions. For the first time we quantify the separate effect of protein topological and energetic disorder on the hydration water dynamics. When a static protein structure is simulated, we show that both types of disorder contribute to slow down water diffusion, and that allowing for protein motion, increasing the spatial dimensionality of the interface, reduces the anomalous character of hydration water. The rotation of water is, instead, altered by the energetic disorder only; indeed, when electrostatic interactions between the protein and water are switched off, water reorients even faster than in the bulk. The dynamics of water is also related to the collective structure--à voir the hydrogen bond (H-bond) network--formed by the solvent enclosing the protein surface. We show that, as expected for a full hydrated protein, when the protein surface offers pinning sites (charged or polar sites), the superficial water-water H-bond network percolates throughout the whole surface, hindering the water diffusion, whereas it does not when the protein surface lacks electrostatic interactions with water and the water diffusion is enhanced.

Single particle and collective hydration dynamics for hydrophobic and hydrophilic peptides

We have conducted extensive molecular dynamics simulations to study the single particle and collective dynamics of water in solutions of N-acetyl-glycine-methylamide, a model hydrophilic protein backbone, and N-acetyl-leucine-methylamide, a model (amphiphilic) hydrophobic peptide, as a function of peptide concentration. Various analytical models commonly used in the analysis of incoherent quasielastic neutron scattering (QENS), are tested against the translational and rotational intermediate scattering function, the mean square displacement of the water molecule center of mass, and fits to the second-order rotational correlation function of water evaluated directly from the simulation data. We find that while the agreement between the model-free analysis and analytical QENS models is quantitatively poor, the qualitative feature of dynamical heterogeneity due to caging is captured well by all approaches. The center of mass collective and single particle intermediate scattering functions of water calculated for these peptide solutions show that the crossover from collective to single particle-dominated motions occurs at a higher value of Q for high concentration solutions relative to low concentration because of the greater restriction in movement of water molecules due to confinement. Finally, we have shown that at the same level of confinement of the two peptides, the aqueous amphiphilic amino acid solution shows the strongest deviation between single particle and collective dynamics relative to the hydrophilic amino acid, indicating that chemical heterogeneity induces even greater spatial heterogeneity in the water dynamics.

Effects of hydration water on protein methyl group dynamics in solution

Elastic and quasielastic neutron scattering experiments have been used to investigate the dynamics of methyl groups in a protein-model hydrophobic peptide in solution. The results suggest that, when the hydrophobic side chains are hydrated by a single hydration water layer, the only allowed motions are confined and attributed to librational and rotational movement associated with the methyl groups. They provide unique experimental evidence that the structural and dynamical properties of the interfacial water strongly influence the side-chain dynamics and the activation of diffusive motion.

Inhomogeneity effects on the structure and dynamics of water at the surface of a membrane: A computer simulation study

The authors report the structural and dynamical properties of water interacting with the surface of a lipid bilayer. Three regions have been identified, which show different dynamical regimes of water: a region of strong water-solute interaction, a transition region, and the bulk water region. The dynamics of the strong-interacting water is dominated by caging effects, as shown by the analysis of the self-intermediate scattering function, and by the disrupture of water's hydrogen bond network, while the smooth transition to bulk water is traced back to the roughness of the bilayer surface.