Hydration in protein dynamics and function

Contributions of higher-order proximal distribution functions to solvent structure around proteins

The proximal distribution function (pDF) quantifies the probability of finding a solvent molecule in the vicinity of solutes. The approach constitutes a hierarchically organized theory for constructing approximate solvation structures around solutes. Given the assumption of universality of atom cluster-specific solvation, reconstruction of the solvent distribution around arbitrary molecules provides a computationally convenient route to solvation thermodynamics. Previously, such solvent reconstructions usually considered the contribution of the nearest-neighbor distribution only. We extend the pDF reconstruction algorithm to terms including next-nearest-neighbor contribution. As a test, small molecules (alanine and butane) are examined. The analysis is then extended to include the protein myoglobin in the P6 crystal unit cell. Molecular dynamics simulations are performed, and solvent density distributions around the solute molecules are compared with the results from different pDF reconstruction models. It is shown that the next-nearest-neighbor modification significantly improves the reconstruction of the solvent number density distribution in concave regions and between solute molecules. The probability densities are then used to calculate the solute-solvent non-bonded interaction energies including van der Waals and electrostatic, which are found to be in good agreement with the simulated values.

Surface Topography Effects of Globular Biomolecules on Hydration Water

Hydration water serves as a microscopic manifestation of structural stability and functions of biomolecules. To develop bio-nanomaterials in applications, it is important to study how the surface topography and heterogeneity of biomolecules result in their diversity of the hydration dynamics and energetics. We here performed molecular dynamics simulations combined with the steered molecular dynamics and umbrella sampling to investigate the dynamics and escape process associated with the free energy change of water molecules close to a globular biomolecule, i.e., hemoglobin (Hb) and G-quadruplex DNA (GDNA). The residence time, power of long-time tail, and dipole relaxation time were found to display drastic changes within the averaged hydration shell of 3.0-5.0 Å. Compared with bulk water, in the inner hydration shell, the water dipole moment displays a slower relaxation process and is more oriented toward GDNA than toward Hb, forming a hedgehog-like structure when it surrounds GDNA. In particular, a spine water structure is observed in the GDNA narrow groove. The water isotope effect not only prolongs the dynamic time scales of libration motion in the inner hydration shell and the dipole relaxation processes in the bulk but also strengthens the DNA spine water structure. The potential of the mean force profile reflects the integrity of the hydration shell structure and enables us to obtain detailed insights into the structures formed by water, such as the caged H-bond network and the edge bridge structures; it also reveals that local hydration shell free energy (LHSFE) depends on H-bond rupture processes and ranges from 0.2 to 4.2 kcal/mol. Our results demonstrate that the surface topography of a biomolecule influences the integrity of the hydration shell structure and LHSFE. Our studies are able to identify various further applications in the areas of microfluid devices and nano-dewetting on bioinspired surfaces.

Vibrational dynamics of confined supercooled water

The quest for a possible liquid-liquid coexistence line in supercooled water below its homogeneous nucleation temperature is faced by confining water within a porous silica substrate (MCM-41). This system is investigated by synchrotron radiation infrared spectroscopy, exploring both the intramolecular and the intermolecular vibrational dynamics, in the temperature range from ambient down to ∼120 K, along several isobaric paths between 0.7 kbar and 3.0 kbar. Upon lowering the temperature, the OH-stretching band shows that the intramolecular vibrational dynamics continuously evolves from predominantly liquidlike to predominantly icelike. An abrupt change in the line shape of the intermolecular vibrational band between 220 K and 240 K, depending on the pressure, is the signature of nucleation of ice within the MCM-41 pores. These findings do not support the presence of two liquid phases and provide evidence for the coexistence of liquid water and ice in water confined in MCM-41.

The adaptive nature of protein residue networks

Protein residue networks PRNs are used to describe proteins. These networks are usually based on an average structure for the protein. However, proteins are dynamic entities that are affected by their surroundings. In this work, we study the effect of temperatures above and below the protein dynamical transition temperature(≈200 K), on three important network parameters gleaned from weighted PRNs for the solvated β-lactamase inhibitory protein BLIP: the betweenness centrality B, the closeness centrality C, and the clustering coefficient CC. The B and C values will be extracted for each node from PRNs at six different temperatures: 150 K, 180 K, 200 K, 220 K, 250 K, and 310 K respectively. The average value for the CC for each PRN will also be calculated at each temperature, respectively. We find that at temperatures ≤200 K, the network nodes with the most significant B and C values tend to have lower relative solvent accessibility RSA values, and to fall within the protein secondary structure elements (α helices and β sheets). At temperatures >200 K, the significant nodes in terms of B and C tend to have larger RSA values, and to fall on the connecting loops in the protein. The average CC decreases in value for the PRNs up to 200 K, and then remains basically constant above 200 K. This clearly shows that any conclusions based on static PRNs should be handled with care. The dynamic nature of proteins and its coupling to the surrounding environment should be taken into consideration when using the PRN paradigm. Proteins 2017; 85:917-923. © 2016 Wiley Periodicals, Inc.

The fluctuating ribosome: thermal molecular dynamics characterized by neutron scattering

Conformational changes associated with ribosome function have been identified by X-ray crystallography and cryo-electron microscopy. These methods, however, inform poorly on timescales. Neutron scattering is well adapted for direct measurements of thermal molecular dynamics, the ‘lubricant’ for the conformational fluctuations required for biological activity. The method was applied to compare water dynamics and conformational fluctuations in the 30 S and 50 S ribosomal subunits from Haloarcula marismortui, under high salt, stable conditions. Similar free and hydration water diffusion parameters are found for both subunits. With respect to the 50 S subunit, the 30 S is characterized by a softer force constant and larger mean square displacements (MSD), which would facilitate conformational adjustments required for messenger and transfer RNA binding. It has been shown previously that systems from mesophiles and extremophiles are adapted to have similar MSD under their respective physiological conditions. This suggests that the results presented are not specific to halophiles in high salt but a general property of ribosome dynamics under corresponding, active conditions. The current study opens new perspectives for neutron scattering characterization of component functional molecular dynamics within the ribosome.

THz frequency spectrum of protein-solvent interaction energy using a recurrence plot-based Wiener-Khinchin method

The dynamics of a protein and the water surrounding it are coupled via nonbonded energy interactions. This coupling can exhibit a complex, nonlinear, and nonstationary nature. The THz frequency spectrum for this interaction energy characterizes both the vibration spectrum of the water hydrogen bond network, and the frequency range of large amplitude modes of proteins. We use a Recurrence Plot based Wiener-Khinchin method RPWK to calculate this spectrum, and the results are compared to those determined using the classical auto-covariance-based Wiener-Khinchin method WK. The frequency spectra for the total nonbonded interaction energy extracted from molecular dynamics simulations between the β-Lactamase Inhibitory Protein BLIP, and water molecules within a 10 Å distance from the protein surface, are calculated at 150, 200, 250, and 310 K, respectively. Similar calculations are also performed for the nonbonded interaction energy between the residues 49ASP, 53TYR, and 142PHE in BLIP, with water molecules within 10 Å from each residue respectively at 150, 200, 250, and 310 K. A comparison of the results shows that RPWK performs better than WK, and is able to detect some frequency data points that WK fails to detect. This points to the importance of using methods capable of taking the complex nature of the protein-solvent energy landscape into consideration, and not to rely on standard linear methods. In general, RPWK can be a valuable addition to the analysis tools for protein molecular dynamics simulations. Proteins 2016; 84:1549-1557. © 2016 Wiley Periodicals, Inc.

Shape dependence of the radial distribution function of hydration water around proteins

The characterization of hydration water in proteins is important to understand their structure, function and folding properties. A calculation of the shape parameters reveals considerable asymmetry in the shapes of globular proteins. The present study suggests a generalized approach for the calculation of radial distribution of hydration water by accounting for the shape asymmetry in proteins. The surface and radial distribution function is analyzed for three groups of high resolution globular proteins and nonglobular proteins. This generalized approach depicts a considerable difference in the hydration water distribution pattern around aspherical proteins as compared to the earlier method and this difference is more pronounced for nonglobular proteins. The peaks for normalized RDF are found to be sharper compared to those of bulk water. The normalized RDF of hydration water exactly coincides with that of the bulk beyond 8.0 Å. The radial distribution of hydration water as a function of the water-protein distance matches with the experimentally observed distribution of hydration water around myoglobin. The results reveal that the distribution of hydration water is dependent on the shapes of proteins and hence a generalized approach should be used for the calculation of hydration water distribution around proteins, especially for nonglobular proteins.

The glassy state of crambin and the THz time scale protein-solvent fluctuations possibly related to protein function

BACKGROUND

THz experiments have been used to characterize the picosecond time scale fluctuations taking place in the model, globular protein crambin.

RESULTS

Using both hydration and temperature as an experimental parameter, we have identified collective fluctuations (<= 200 cm(-1)) in the protein. Observation of the protein dynamics in the THz spectrum from both below and above the glass transition temperature (Tg) has provided unique insight into the microscopic interactions and modes that permit the solvent to effectively couple to the protein thermal fluctuations.

CONCLUSIONS

Our findings suggest that the solvent dynamics on the picosecond time scale not only contribute to protein flexibility but may also delineate the types of fluctuations that are able to form within the protein structure.

Local order and mobility of water molecules around ambivalent helices

Water on a protein surface plays a key role in determining the structure and dynamics of proteins. Compared to the properties of bulk water, many aspects of the structure and dynamics of the water surrounding the proteins are less understood. It is interesting therefore to explore how the properties of the water within the solvation shell around the peptide molecule depend on its specific secondary structure. In this work we investigate the orientational order and residence times of the water molecules to characterize the structure, energetics, and dynamics of the hydration shell water around ambivalent peptides. Ambivalent sequences are identical sequences which display multiple secondary structures in different proteins. Molecular dynamics simulations of representative proteins containing variable helix, variable nonhelix, and conserved helix are also used to explore the local structure and mobility of water molecules in their vicinity. The results, for the first time, depict a different water distribution pattern around the conserved and variable helices. The water molecules surrounding the helical segments in variable helices are found to possess a less locally ordered structure compared to those around their corresponding nonhelical counterparts and conserved helices. The long conserved helices exhibit extremely high local residence times compared to the helical conformations of the variable helices, whereas the residence times of the nonhelical conformations of the variable helices are comparable to those of the short conserved helices. This differential pattern of the structure and dynamics of water molecules in the vicinity of conserved/variable helices may lend valuable insights for understanding the role of solvent effects in determining sequence ambivalency and help in improving the accuracy of water models used in the simulations of proteins.

Dynamical properties of the hydration shell of fully deuterated myoglobin

Freeze-dried perdeuterated sperm whale myoglobin was kept in a water-saturated atmosphere in order to obtain a hydration degree of 335 H(2)O molecules per one myoglobin molecule. Incoherent neutron scattering was performed at the neutron spectrometer TOFTOF at the FRM II in an angular range of q from 0.6 to 1.8 Å(-1) and a temperature range from 4 to 297 K. We used neutrons with a wavelength of λ αE 6 Å and an energy resolution of about 65 μeV corresponding to motions faster than 10 ps. At temperatures above 225 K, broad lines appear in the spectra caused by quasielastic scattering. For an explanation of these lines, we assumed that there are only two types of protons, those that are part of the hydration water (72%) and those that belong to the protein (28%). The protons of the hydration water were analyzed with the diffusion model of Singwi and Sjölander [Phys. Rev. 119, 863 (1960)]. In this model, a water molecule stays for a time τ(0) in a bound state performing oscillatory motions. Thereafter, the molecule performs free diffusion for the time τ(1) in a nonbound state followed again by the oscillatory motions for τ(0) and so forth. We used the general formulation with no simplifications as τ(0)≫τ(1) or τ(1)≫τ(0). At room temperature, we obtained τ(0) αE 104 ps and τ(1) αE 37 ps. For the protein bound hydrogen, the dynamics is described by a Brownian oscillator where the protons perform overdamped motions in limited space.

Structure of hydration water in proteins: a comparison of molecular dynamics simulations and database analysis

Hydration layer water molecules play important structural and functional roles in proteins. Despite being a critical component in biomolecular systems, characterizing the properties of hydration water poses a challenge for both experiments and simulations. In this context we investigate the local structure of hydration water molecules as a function of the distance from the protein and water molecules respectively in 188 high resolution protein structures and compare it with those obtained from molecular dynamics simulations. Tetrahedral order parameter of water in proteins calculated from previous and present simulation studies show that the potential of bulk water overestimates the average tetrahedral order parameter compared to those calculated from crystal structures. Hydration waters are found to be more ordered at a distance between the first and second solvation shell from the protein surface. The values of the order parameter decrease sharply when the water molecules are located very near or far away from the protein surface. At small water-water distance, the values of order parameter of water are very low. The average order parameter records a maximum value at a distance equivalent to the first solvation layer with respect to the water-water radial distribution and asymptotically approaches a constant value at large distances. Results from present analysis will help to get a better insight into structure of hydration water around proteins. The analysis will also help to improve the accuracy of water models on the protein surface.

Temperature dependence of hydrogen-bond dynamics in acetic acid-water solutions

An inelastic UV scattering experiment has been carried out on acetic acid-water solutions as a function of temperature and concentration. The analysis of experimental data indicates the presence of a crossover temperature (T(c) approximately 325 +/- 10 K). Above T(c), the energy of hydrogen bonds responsible for water-acetic acid and acetic acid-acetic acid interactions is strongly reduced. This leads to a reduction in the average number of water molecule interacting with acetic acid, as well as to a lower number of acetic acid clusters. The latter behavior can be mainly ascribed to a temperature change in the activation energy of carboxylic groups of acetic acid. These results may be also relevant to better understand the folding mechanism in protein-water solutions.

Structure and dynamics of water in native and tanned collagen fibers: Effect of crosslinking

The influence of crosslinking on the hydration structure of collagen has been investigated. Nuclear magnetic resonance, dielectric relaxation and thermoporometry were used to investigate water structure in native and crosslinked collagen fibers on both wet and dried specimen. Measurements reveal the influence of different chemical treatments on the transverse relaxation time and polarization of the collagen fibers. The frequency dependence of dielectric constant of collagen fibers displays an induction behavior on low frequencies. Bound water constrained in collagen fibers seems to provide signatures for changes induced by crosslinking agents on the pore diameter and distribution in collagen fibers. A correlation of transverse relaxation time of water in dry and wet states presented in this study presents an experimental tool for examining the differences in efficacy of crosslinking agents. Changes in the dielectric relaxation, dynamics of water structure and hydroporometric structure of collagen are dependent on the nature of crosslinking material.

Near infrared spectroscopy and aquaphotomics: Novel approach for rapid in vivo diagnosis of virus infected soybean

Near infrared spectroscopy with aquaphotomics as a novel approach was assessed for the diagnosis of soybean plants (Glycine max) infected with soybean mosaic virus (SMV) at latent symptomless stage of the disease. Near infrared (NIR) leaf spectra (in the range of 730-1025nm) acquired from soybean plants with and without the inoculation of SMV were used. Leaf samples from all plants were assayed with enzyme-linked immunosorbent assay (ELISA) to confirm the infection. Previously reported NIR band for water at 970nm and two new bands at 910nm and 936nm in the water specific region of NIR were found to be markedly sensitive to the SMV infection 2weeks prior to the appearance of visual symptoms on infected leaves. The spectral calibration model soft independent modeling of class analogy (SIMCA), predicted the disease with 91.6% sensitivity and 95.8% specificity when the second order derivative of the individual plant averaged spectra were used. The study shows the potential of NIR spectroscopy with its novel approach to elucidate latent biochemical and biophysical information of an infection as it allowed successful discrimination of SMV infected plant from healthy at the early symptomless stage of the disease.

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.

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.

Fingerprints of amorphous icelike behavior in the vibrational density of states of protein hydration water

The low-frequency modes of protein hydration water are investigated by inelastic neutron scattering. Experiments on both protonated and fully deuterated maltose binding protein samples allow us to unambiguously single out the contribution from water. The low-energy vibrational density of states of hydration water at 100 K is similar to the density of states of high- and low-density amorphous ice, and quite different from that of simple forms of crystalline ice. This result can be related to the picture of hydration water mass density depending on the protein surface curvature, which supports its glassy behavior.

Thermal breaking of spanning water networks in the hydration shell of proteins

The presence of a spanning hydrogen-bonded network of water at the surface of biomolecules is important for their conformational stability, dynamics, and function. We have studied by computer simulations the clustering and percolation of water in the hydration shell of a small elastinlike peptide (ELP) and the medium-size protein staphylococcal nuclease (SNase), in aqueous solution. We have found that in both systems a spanning network of hydration water exists at low temperatures and breaks up with increasing temperature via a quasi-two-dimensional percolation transition. The thermal breaking of the spanning water network occurs at biologically relevant temperatures, in the temperature range, which is close to the temperature of the "inverse temperature transition" of ELP and the unfolding temperature of SNase, respectively.

Ultrafast optical Kerr effect spectroscopy of water confined in nanopores of the gelatin gel

We report on the investigation of a short-time collective dynamics of water confined in the pores of the gelatin gel, using the femtosecond optical Kerr effect spectroscopy. The ultrafast responses of water molecules obtained in bulk liquid and in three concentrations of gelatin gels are explained theoretically, both in a long time and in a short time regime, taking into account all molecular motions. We prove that the contribution of molecules involved in tetrahedral, strongly H-bonded structures stabilizing the gel network increases with the gel concentration. On the other hand the long-time relaxation of water molecules is significantly slowed down in the gel pores.

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.

Protein Dynamics Tightly Connected to the Dynamics of Surrounding and Internal Water Molecules

Proteins are key components of biological cells. For example, enzymes catalyze biochemical reactions, membrane transporters are responsible for uptake and release of critical and superfluous components from the cell environment, and structural proteins are responsible for the stability of the cell wall and cytoskeleton. Many of the diverse protein functions involve dynamic transitions ranging from small local atomic displacements up to large allosteric conformational changes. In any conformation, proteins are in contact with the universal solvent medium of cells, water. Water not only surrounds proteins but is often an integral part of proteins and also is involved in key mechanistic steps. This Minireview discusses recent experimental and theoretical results on the role of water for protein dynamics and function.

Water Dynamics in Graphite Oxide Investigated with Neutron Scattering

Graphite oxide is an inorganic multilayer system that preserves the layered structure of graphite but not the conjugated bond structure. In the past few years, detailed studies of the static structure of graphite oxide were carried out. This was mainly done by NMR investigations and led to a new structural model of graphite oxide. The layer distance of graphite oxide increases with increasing humidity level, giving rise to different spacings of the carbon layers in the range from 6 to 12 A. As a consequence, different types of motions of water and functional groups appear. Information about the mobility of the water molecules is not yet complete but is crucial for the understanding of the structure of the carbon layers as well as the intercalation process. In this paper, the hydration- and temperature-dependent dynamic behavior of graphite oxide will be investigated by quasielastic neutron scattering using the time-of-flight spectrometer NEAT at the Hahn-Meitner-Institut Berlin. The character of the embedded water does not change over a wide range of hydration levels. Especially the interlayer water remains tightly bound and does not show any translational motion. In samples with excess water, however, the water is also distributed in noninterlayer voids, leading to the observation of additional motions of bulklike or confined water. The dynamic behavior of hydrated graphite oxide can be described by a consistent model that combines two two-site jump motions for the motions of the water molecules and the motions of OH groups.

Water structure and interactions with protein surfaces

The structure of liquid water and its interaction with biological molecules is a very active area of experimental and theoretical research. The chemically complex surfaces of protein molecules alter the structure of the surrounding layer of hydrating water molecules. The dynamics of hydration water can be detected by a series of experimental techniques, which show that hydration waters typically have slower correlation times than water in bulk. Specific water-mediated interactions in protein complexes have been studied in detail, and these interactions have been incorporated into potential energy functions for protein folding and design. The subtle changes in the structure of hydration water have been investigated by theoretical studies.

Existence of hydration forces in the interaction between apoferritin molecules adsorbed on silica surfaces

The atomic force microscope, together with the colloid probe technique, has become a very useful instrument to measure interaction forces between two surfaces. Its potential has been exploited in this work to study the interaction between protein (apoferritin) layers adsorbed on silica surfaces and to analyze the effect of the medium conditions (pH, salt concentration, salt type) on such interactions. It has been observed that the interaction at low salt concentrations is dominated by electrical double layer (at large distances) and steric forces (at short distances), the latter being due to compression of the protein layers. The DLVO theory fits these experimental data quite well. However, a non-DLVO repulsive interaction, prior to contact of the protein layers, is observed at high salt concentration above the isoelectric point of the protein. This behavior could be explained if the presence of hydration forces in the system is assumed. The inclusion of a hydration term in the DLVO theory (extended DLVO theory) gives rise to a better agreement between the theoretical fits and the experimental results. These results seem to suggest that the hydration forces play a very important role in the stability of the proteins in the physiological media.

Molecular view of water dynamics near model peptides

Incoherent quasi-elastic neutron scattering (QENS) has been used to measure the dynamics of water molecules in solutions of a model protein backbone, N-acetyl-glycine-methylamide (NAGMA), as a function of concentration, for comparison with results for water dynamics in aqueous solutions of the N-acetyl-leucine-methylamide (NALMA) hydrophobic peptide at comparable concentrations. From the analysis of the elastic incoherent structure factor, we find significant fractions of elastic intensity at high and low concentrations for both solutes, which corresponds to a greater population of protons with rotational time scales outside the experimental resolution (>13 ps). The higher-concentration solutions show a component of the elastic fraction that we propose is due to water motions that are strongly coupled to the solute motions, while for low-concentration solutions an additional component is activated due to dynamic coupling between inner and outer hydration layers. An important difference between the solute types at the highest concentration studied is found from stretched exponential fits to their experimental intermediate scattering functions, showing more pronounced anomalous diffusion signatures for NALMA, including a smaller stretched exponent beta and a longer structural relaxation time tau than those found for NAGMA. The more normal water diffusion exhibited near the hydrophilic NAGMA provides experimental support for an explanation of the origin of the anomalous diffusion behavior of NALMA as arising from frustrated interactions between water molecules when a chemical interface is formed upon addition of a hydrophobic side chain, inducing spatial heterogeneity in the hydration dynamics in the two types of regions of the NALMA peptide. We place our QENS measurements on model biological solutes in the context of other spectroscopic techniques and provide both confirming as well as complementary dynamic information that attempts to give a unifying molecular view of hydration dynamics signatures near peptides and proteins.

Enzyme activity and flexibility at very low hydration

Recent measurements have demonstrated enzyme activity at hydrations as low as 3%. This raises the question of whether hydration-induced enzyme flexibility is important for activity. Here, to address this, picosecond dynamic neutron scattering experiments are performed on pig liver esterase powders at 0%, 3%, 12%, and 50% hydration by weight and at temperatures ranging from 120 to 300 K. At all temperatures and hydrations, significant quasielastic scattering intensity is found in the protein, indicating the presence of anharmonic, diffusive motion. As the hydration increases, a temperature-dependent dynamical transition appears and strengthens involving additional diffusive motion. The implication of these results is that, although the additional hydration-induced diffusive motion in the protein detected here may be related to increased activity, it is not required for the enzyme to function.

Hydration dynamics near a model protein surface

The evolution of water dynamics from dilute to very high concentration solutions of a prototypical hydrophobic amino acid with its polar backbone, N-acetyl-leucine-methylamide (NALMA), is studied by quasi-elastic neutron scattering (QENS) and molecular dynamics (MD) simulation for both the completely deuterated and completely hydrogenated leucine monomer. The NALMA-water system and the QENS data together provide a unique study for characterizing the dynamics of different hydration layers near a prototypical hydrophobic side chain and the backbone of which it is attached. We observe several unexpected features in the dynamics of these biological solutions under ambient conditions. The NALMA dynamics shows evidence of de Gennes narrowing, an indication of coherent long timescale structural relaxation dynamics. The translational and rotational water dynamics at the highest solute concentrations are found to be highly suppressed as characterized by long residential time and slow diffusion coefficients. The analysis of the more dilute concentration solutions models the first hydration shell with the 2.0 M spectra. We find that for outer layer hydration dynamics that the translational diffusion dynamics is still suppressed, although the rotational relaxation time and residential time are converged to bulk-water values. Molecular dynamics analysis of the first hydration shell water dynamics shows spatially heterogeneous water dynamics, with fast water motions near the hydrophobic side chain, and much slower water motions near the hydrophilic backbone. We discuss the hydration dynamics results of this model protein system in the context of protein function and protein-protein recognition.

Analysis of CDK2 active-site hydration: a method to design new inhibitors

The interactions between the protein and the solvent were analyzed, and protein regions with a high density of water molecules, as well as structural water molecules, were determined by using molecular dynamics (MD) simulations. A number of water molecules that were in contact with the protein for the whole trajectory were determined. Their interaction energies and hydrogen bonds with protein residues were analyzed. Altogether, 39, 27, 49, and 32 water molecules bound to the protein were found for trajectories of the free CDK2, CDK2/ATP, CDK2/roscovitine, and CDK2/isopentenyladenine complexes, respectively. Positions of observed water molecules were compared with X-ray crystallography data. Special attention was paid to water molecules in the active site of the enzyme, and especially to the deep pocket, where the N9 roscovitine side-chain is buried. Exchange of active-site water molecules with bulk water through the tunnel from the pocket was observed. In the CDK2/isopentenyladenine complex simulation, two water molecules that arrange interaction between the inhibitor and the enzyme via an H-bond were observed. Two stable water molecules in the trajectory of the free CDK2 were found that occupy the same position as the nitrogens N3 and N9 of the isopentenyladenine or N1 and N6 nitrogens of the adenosine triphosphate (ATP). The positions of structural water molecules were compared with the positions of substrate polar groups and crystallographic water molecules found in the Brookhaven Protein Data Bank for various CDK2 complexes. It was concluded that tracing tightly bound water molecules may substantially help in designing new inhibitors.

Structural effects of drying and rehydration for enzymes in soils: a kinetics-FTIR analysis of α-chymotrypsin adsorbed on montmorillonite

The effects of desiccation and rehydration cycles encountered by extracellular enzymes in soils are studied on -chymotrypsin adsorbed on montmorillonite. The controlled hygrometric FTIR cell used in this study enables to monitor drying and rehydration steps undergone by the -chymotrypsin-montmorillonite suspension or by the enzyme alone. Relative humidity (RH) determines the amount of deuterated water in the FTIR cell atmosphere. The molar water/protein ratio (W/P) as well as the conformational and solvation states of the enzyme have been determined using H/D exchange monitored by FTIR-transmission spectroscopy. When the W/P ratio decreases from 3500 to approximately 400, unfolding of beta-secondary structure in three different domains involves about 8% of the polypeptide backbone with respect to the most solvated states. Desiccation induces beta-unfolding, which opens channels allowing free vapor water molecules to diffuse into the enzyme at 15% RH. On drying to 0% RH, displacements of internal water (H2O) in the enzyme are demonstrated by reverse peptide isotopic exchanges (COND ==> CONH). Specific beta-structures, only formed in highly solvated states, sequester around 20 internal H2O molecules. Indeed, most of the unfolded secondary structures during the drying step are refolded at W/P approximately 1000 during rehydration. However, self-association hinders the recovery of the initial closed tertiary structure. The pD-dependent structural changes controlling inward and outward water diffusion are suppressed, whether the protein is initially in an adsorbed state or in solution. Changes in secondary structures encountered during desiccation/rehydration cycle are similar for the protein either free or in the adsorbed state. Thus domains that are unfolded by adsorption are not concerned by the desiccation/rehydration cycle.

Protein-water interactions in a dynamic world

Protein-water interactions are key to biological function. They have an underlying dynamic component that pervades the functional roles associated both with particular systems and with the properties of proteins in general. This article focuses on the specific ways in which the dynamics of water are important to protein structure, motion and adaptability to changes in the protein environment.

Freeze-substitution: origins and applications

Freeze-substitution is a physicochemical process in which biological specimens are immobilized and stabilized for microscopy. Water frozen within cells is replaced by organic solvents at subzero temperatures. Freeze-substitution is widely used for ultrastructural and immunocytochemical analyses of cells by transmission and scanning electron microscopy. Less well recognized is its superiority over conventional chemical fixation in preserving labile and rare tissue antigens for immunocytochemistry by light microscopy. In the postgenome era, the focus of molecular genetics will shift from analyzing DNA sequence structure to elucidating the function of gene networks, the intercellular effects of polygenetic diseases, and the conformational rearrangements of proteins in situ. Novel strategies will be needed to integrate knowledge of chemical structures of normal and abnormal macromolecules with the physiology and developmental biology of cells and tissues from whole organisms. This review summarizes the progress and future prospects of freeze-substitution for such explorations.