Molecular dynamics simulations of the two disaccharides of hyaluronan in aqueous solution

Hyaluronan is an unusually stiff polymer when in aqueous solution,which has important consequences for its biological function.Molecular dynamics simulations of hyaluronan disaccharides havebeen performed, with explicit inclusion of water, to determinethe molecular basis of this stiffness, and to investigate thedynamics of the glycosidic linkages. Our simulations revealthat stable sets of hydrogen bonds frequently connect the neighboringresidues of hyaluronan. Water caging around the glycosidic linkagewas observed to increase the connectivity between sugars, andfurther constrain them. This, we propose, explains the unusualstiffness of polymeric hyaluronan. It would allow the polysaccharideto maintain local secondary structure, and occupy large solutiondomains consistent with the visco-elastic nature of hyaluronan.Simulations in water showed no significant changes on inclusionof the exo-anomeric effect. This, we deduced, was due to hyaluronandisaccharides ordering first shell water molecules. In somecases these waters were observed to transiently induce con-formationalchange, by breaking intramolecular hydrogen bonds. conformation hyaluronan hydrogen bonds molecular dynamics water     

Long-range DNA-water interactions

DNA functions only in aqueous environments and adopts different conformations depending on the hydration level. The dynamics of hydration water and hydrated DNA leads to rotating and oscillating dipoles that, in turn, give rise to a strong megahertz to terahertz absorption. Investigating the impact of hydration on DNA dynamics and the spectral features of water molecules influenced by DNA, however, is extremely challenging because of the strong absorption of water in the megahertz to terahertz frequency range. In response, we have employed a high-precision megahertz to terahertz dielectric spectrometer, assisted by molecular dynamics simulations, to investigate the dynamics of water molecules within the hydration shells of DNA as well as the collective vibrational motions of hydrated DNA, which are vital to DNA conformation and functionality. Our results reveal that the dynamics of water molecules in a DNA solution is heterogeneous, exhibiting a hierarchy of four distinct relaxation times ranging from ∼8 ps to 1 ns, and the hydration structure of a DNA chain can extend to as far as ∼18 Å from its surface. The low-frequency collective vibrational modes of hydrated DNA have been identified and found to be sensitive to environmental conditions including temperature and hydration level. The results reveal critical information on hydrated DNA dynamics and DNA-water interfaces, which impact the biochemical functions and reactivity of DNA.     

Hydration of enzyme in nonaqueous media is consistent with solvent dependence of its activity

Water plays an important role in enzyme structure and function in aqueous media. That role becomes even more important when one focuses on enzymes in low water media. Here we present results from molecular dynamics simulations of surfactant-solubilized subtilisin BPN' in three organic solvents (octane, tetrahydrofuran, and acetonitrile) and in pure water. Trajectories from simulations are analyzed with a focus on enzyme structure, flexibility, and the details of enzyme hydration. The overall enzyme and backbone structures, as well as individual residue flexibility, do not show significant differences between water and the three organic solvents over a timescale of several nanoseconds currently accessible to large-scale molecular dynamics simulations. The key factor that distinguishes molecular-level details in different media is the partitioning of hydration water between the enzyme and the bulk solvent. The enzyme surface and the active site region are well hydrated in aqueous medium, whereas with increasing polarity of the organic solvent (octane --> tetrahydrofuran --> acetonitrile) the hydration water is stripped from the enzyme surface. Water stripping is accompanied by the penetration of tetrahydrofuran and acetonitrile molecules into crevices on the enzyme surface and especially into the active site. More polar organic solvents (tetrahydrofuran and acetonitrile) replace mobile and weakly bound water molecules in the active site and leave primarily the tightly bound water in that region. In contrast, the lack of water stripping in octane allows efficient hydration of the active site uniformly by mobile and weakly bound water and some structural water similar to that in aqueous solution. These differences in the active site hydration are consistent with the inverse dependence of enzymatic activity on organic solvent polarity and indicate that the behavior of hydration water on the enzyme surface and in the active site is an important determinant of biological function especially in low water media.     

Hydration of polysaccharide hyaluronan observed by IR spectrometry. I. Preliminary experiments and band assignments

This article is the first one in a series dedicated to the study of hyaluronan as observed by IR spectrometry. The goal is to determine its hydration mechanism and the structural changes this mechanism implies. Hyaluronan is a natural polysaccharide that is widely used in biomedical applications and cosmetics. Its macroscopic properties are significantly dependent on its degree of hydration. In this article we record the IR spectrum of a several micron thick dried film and deduce that four or five residual H(2)O molecules remain around each disaccharide repeat unit in the dried film. We then compare the spectra of sodium hyaluronan and its acid form to assign vibrational bands linked to the carboxylate group. We proceed with a qualitative analysis of the spectral changes induced by changes of temperature and hygroscopicity, two independent parameters that act by modifying the hydrogen bond network of the sample. This enables us to assign most of the vibrational bands of the hydrophilic groups and to distinguish the bands that are due to these hydrophilic groups when they are or are not hydrogen bonded. It constitutes a prerequisite for the quantitative analysis of hydration spectra that will be described in the following articles of this series.     

The dynamics of protein hydration water: a quantitative comparison of molecular dynamics simulations and neutron-scattering experiments

We present results from an extensive molecular dynamics simulation study of water hydrating the protein Ribonuclease A, at a series of temperatures in cluster, crystal, and powder environments. The dynamics of protein hydration water appear to be very similar in crystal and powder environments at moderate to high hydration levels. Thus, we contend that experiments performed on powder samples are appropriate for discussing hydration water dynamics in native protein environments. Our analysis reveals that simulations performed on cluster models consisting of proteins surrounded by a finite water shell with free boundaries are not appropriate for the study of the solvent dynamics. Detailed comparison to available x-ray diffraction and inelastic neutron-scattering data shows that current generation force fields are capable of accurately reproducing the structural and dynamical observables. On the time scale of tens of picoseconds, at room temperature and high hydration, significant water translational diffusion and rotational motion occur. At low hydration, the water molecules are translationally confined but display appreciable rotational motion. Below the protein dynamical transition temperature, both translational and rotational motions of the water molecules are essentially arrested. Taken together, these results suggest that water translational motion is necessary for the structural relaxation that permits anharmonic and diffusive motions in proteins. Furthermore, it appears that the exchange of protein-water hydrogen bonds by water rotational/librational motion is not sufficient to permit protein structural relaxation. Rather, the complete exchange of protein-bound water molecules by translational displacement seems to be required.     

Contribution of solvent water to the solution X-ray scattering profile of proteins

A theoretical framework is presented to analyze how solvent water contributes to the X-ray scattering profile of protein solution. Molecular dynamics simulations were carried out on pure water and an aqueous solution of myoglobin to determine the spatial distribution of water molecules in each of them. Their solution X-ray scattering (SXS) profiles were numerically evaluated with obtained atomic-coordinate data. It is shown that two kinds of contributions from solvent water must be considered to predict the SXS profile of a solution accurately. One is the excluded solvent scattering originating in exclusion of water molecules from the space occupied by solutes. The other is the hydration effect resulting from formation of a specific distribution of water around solutes. Explicit consideration of only two molecular layers of water is practically enough to incorporate the hydration effect. Care should be given to using an approximation in which an averaged electron density distribution is assumed for the structure factor because it may predict profiles considerably deviating from the correct profile at large K.     

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.     

NMR observation of individual molecules of hydration water bound to DNA duplexes: direct evidence for a spine of hydration water present in aqueous solution.

The residence times of individual hydration water molecules in the major and minor grooves of DNA were measured by nuclear magnetic resonance (NMR) spectroscopy in aqueous solutions of d-(CGCGAATTCGCG)2 and d-(AAAAATTTTT)2. The experimental observations were nuclear Overhauser effects (NOE) between water protons and the protons of the DNA. The positive sign of NOEs with the thymine methyl groups shows that the residence times of the hydration water molecules near these protons in the major groove of the DNA must be shorter than about 500 ps, which coincides with the behavior of surface hydration water in peptides and proteins. Negative NOEs were observed with the hydrogen atoms in position 2 of adenine in both duplexes studied. This indicates that a 'spine of hydration' in the minor groove, as observed by X-ray diffraction in DNA crystals, is present also in solution, with residence times significantly longer than 1 ns. Such residence times are reminiscent of 'interior' hydration water molecules in globular proteins, which are an integral part of the molecular architecture both in solution and in crystals.     

Large-scale networks of hydration water molecules around proteins investigated by cryogenic X-ray crystallography.

The structures at protein-water interface, i.e. the hydration structure of proteins, have been investigated by cryogenic X-ray crystal structure analyses. Hydration structures appeared far clearer at cryogenic temperature than at ambient temperature, presumably because the motions of hydration water molecules were quenched by cooling. Based on the structural models obtained, the hydration structures were systematically analyzed with respect to the amount of water molecules, the interaction modes between water molecules and proteins, the local and the global distribution of them on the surface of proteins. The standard tetrahedral interaction geometry of water in bulk retained at the interface and enabled the three-dimensional chain connection of hydrogen bonds between hydration water molecules and polar protein atoms. Large-scale networks of hydrogen bonds covering the entire surface of proteins were quite flexible to accommodate to the large-scale conformational changes of proteins and seemed to have great influences on the dynamics and function of proteins. The present observation may provide a new concept for discussing the dynamics of proteins in aqueous solution.     

Dynamics of hydration in hen egg white lysozyme

We investigate the hydration dynamics of a small globular protein, hen egg-white lysozyme. Extensive simulations (two trajectories of 9 ns each) were carried out to identify the time-scales and mechanism of water attachment to this protein. The location of the surface and integral water molecules in lysozyme was also investigated. Three peculiar temporal scales of the hydration dynamics can be discerned: two among these, with sub-nanosecond mean residence time, tau(w), are characteristic of surface hydration water; the slower time-scale (tau(w) approximately 2/3 ns) is associated with buried water molecules in hydrophilic pores and in superficial clefts. The computed tau(w) values in the two independent runs fall in a similar range and are consistent with each other, thus adding extra weight to our result. The tau(w) of surface water obtained from the two independent trajectories is 20 and 24 ps. In both simulations only three water molecules are bound to lysozyme for the entire length of the trajectories, in agreement with nuclear magnetic relaxation dispersion estimates. Locations other than those identified in the protein crystal are found to be possible for these long-residing water molecules. The dynamics of the hydration water molecules observed in our simulations implies that each water molecule visits a multitude of residues during the lifetime of its bound with the protein. The number of residues seen by a single water molecule increases with the time-scale of its residence time and, on average, is equal to one only for the water molecules with shorter residence time. Thus, tau(w) values obtained from inelastic neutron scattering and based on jump-diffusion models are likely not to account for the contribution of water molecules with longer residence time.     

MC-PHS: A Monte Carlo Implementation of the Primary Hydration Shell for Protein Folding and Design

A primary hydration shell (PHS) approach is developed for Monte Carlo simulations of conformationally rich macromolecular systems in an environment that efficiently captures principal solvation effects. It has been previously demonstrated that molecular dynamics using PHS is an efficient method to study peptide structure and dynamics in aqueous solution. Here, we extend the PHS approach to Monte Carlo simulations, whereby a stable shell of water molecules is maintained with a flexible, nonspherical, half-harmonic potential, tuned to maintain a constant restraining energy, with the difference between the restraint and shell energies used to dynamically adjust the shell radius. Examination of the shell and system size dependence of the restraining potential reveals its robustness. Moreover, its suitability for biomolecular simulations is evaluated using small spheres of water, hydration properties of small biological molecules, and configurational sampling of β-hairpin pentapeptide YPGDV. This method, termed MC-PHS, appears to provide efficient representation of dominant solvation effects and should prove useful in the study of protein folding and design.     

Sequence Dependencies of DNA Deformability and Hydration in the Minor Groove

DNA deformability and hydration are both sequence-dependent and are essential in specific DNA sequence recognition by proteins. However, the relationship between the two is not well understood. Here, systematic molecular dynamics simulations of 136 DNA sequences that differ from each other in their central tetramer revealed that sequence dependence of hydration is clearly correlated with that of deformability. We show that this correlation can be illustrated by four typical cases. Most rigid basepair steps are highly likely to form an ordered hydration pattern composed of one water molecule forming a bridge between the bases of distinct strands, but a few exceptions favor another ordered hydration composed of two water molecules forming such a bridge. Steps with medium deformability can display both of these hydration patterns with frequent transition. Highly flexible steps do not have any stable hydration pattern. A detailed picture of this correlation demonstrates that motions of hydration water molecules and DNA bases are tightly coupled with each other at the atomic level. These results contribute to our understanding of the entropic contribution from water molecules in protein or drug binding and could be applied for the purpose of predicting binding sites.     

Imaging and tracking of single hyaluronan molecules diffusing in solution

Hyaluronan is an important soluble component of the extracellular matrix of many tissues with well known space-filling, lubricating and signaling functions. As such, hyaluronan can regulate cell adhesion, migration, differentiation and proliferation. Ultrastructural studies showed the existence of fibers and networks of hyaluronan molecules at surfaces, while bulk studies of hyaluronan in solution indicated that the polymer forms random coils. Here, we show that single hyaluronan molecules can be visualized and tracked in three-dimensional samples at room temperature in aqueous buffer. Using a wide-field fluorescence microscope equipped with laser excitation and an sensitive and fast EMCCD camera for fluorescence detection, single FITC-labeled hyaluronan molecules from rooster comb were detected in aqueous solutions. Freely moving hyaluronan-FITC could be tracked over up to 20 images acquired at a frame rate of 98 Hz. Analysis of the trajectories revealed Brownian motion of hyaluronan in tris-buffered saline with an average diffusion coefficient D = 3.0 ± 0.2 μm²/s. These observations confirm the concept that hyaluronan molecules form random coils in solution. The possibility of following the tracks of single hyaluronan molecules in solution facilitates the analysis of processes that lead to the formation of more organized forms of hyaluronan and its interactions with cells with very high spatial and temporal accuracy. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

The role of hydration on the mechanism of allosteric regulation: in situ measurements of the oxygen-linked kinetics of water binding to hemoglobin

We report here the first direct measurements of changes in protein hydration triggered by a functional binding. This task is achieved by weighing hemoglobin (Hb) and myoglobin films exposed to an atmosphere of 98% relative humidity during oxygenation. The binding of the first oxygen molecules to Hb tetramer triggers a change in protein conformation, which increases binding affinity to the remaining empty sites giving rise to the appearance of cooperative phenomena. Although crystallographic data have evidenced that this structural change increases the protein water-accessible surface area, isobaric osmotic stress experiments in aqueous cosolutions have shown that water binding is linked to Hb oxygenation. Now we show that the differential hydration between fully oxygenated and fully deoxygenated states of these proteins, determined by weighing protein films with a quartz crystal microbalance, agree with the ones determined by osmotic stress in aqueous cosolutions, from the linkage between protein oxygen affinity and water activity. The agreements prove that the changes in water activity brought about by adding osmolytes to the buffer solution shift biochemical equilibrium in proportion to the number of water molecules associated with the reaction. The concomitant kinetics of oxygen and of water binding to Hb have been also determined. The data show that the binding of water molecules to the extra protein surface exposed on the transition from the low-affinity T to the high-affinity R conformations of hemoglobin is the rate-limiting step of Hb cooperative reaction. This evidences that water binding is a crucial step on the allosteric mechanism regulating cooperative interactions, and suggests the possibility that environmental water activity might be engaged in the kinetic control of some important reactions in vivo.     

A molecular dynamics study of branched α-cyclodextrin

A branched α-cyclodextrin is a derivative of an α-cyclodextrin with a branch consisting of an extra glucose unit. Its water solubility is considerably higher than that of the unbranched one. We have studied the high solubility of the molecule in aqueous solution by molecular dynamics simulations. Trajectories of the molecule at 293 K were calculated using GROMOS programs in three different environments, i.e., in vacuo, in the crystalline state, and in aqueous solution. A simulation in vacuo was carried out to explore stable conformations of the molecule in the isolated system. The quality of the simulations were examined by comparing the X-ray and the simulated crystal structures.The results of the simulations show three remarkable structural features of the molecule: self-inclusion with its branched portion, twist-boat conformation of a glucose ring, and wobbling of its macrocycle. Among these, the last feature is closely related to the water solubility of the molecule. The solubility of cyclodextrin appears to be mainly governed by its intramolecular interglucose hydrogen bonds, which inhibit hydration by solvent water molecules. The results of our simulations indicate that the capability to form hydrogen bonds in branched α-cyclodextrin decreased as the macrocycle of the molecule lost its regular circular shape. Such wobbling of the macrocycle was observed on a relatively short time scale (several picoseconds). An extra glucose unit introduced to α-cyclodextrin may cause the improved water solubility of the molecule through the greater wobbling motion of its macrocycle.