Structure of DNA hydration shells studied by Raman spectroscopy

We have used Raman scattering to study the water O-H stretching modes at approximately 3450 and approximately 3220 cm-1 in DNA films as a function of relative humidity (r.h.). The intensity of the 3220-cm-1 band vanishes as the r.h. is decreased from 98% to around 80%, which indicates that the hydrogen-bond network of water is disrupted in the primary hydration shell (which therefore cannot have an "ice-like" structure). The number of water molecules in the primary hydration shell was determined from the intensity of the approximately 3200-cm-1 band as about 30 water molecules per nucleotide pair. The approximately 3400-cm-1 O-H stretch band was used for determining the total water content, and this band persists at 0% r.h., implying that 5-6 tightly bound water molecules per nucleotide pair remain. The frequency of the approximately 3400-cm-1 O-H stretch mode is lower by 30 to 45 cm-1 in the primary hydration shell compared to free water. The water content as a function of r.h. obtained from these experiments agrees with gravimetric measurements. The disappearance of the approximately 3200-cm-1 band and the shift of the approximately 3400-cm-1 O-H stretch band provide a reliable way of measuring the hydration number of DNA.     

Hydration of B-DNA: comparison between the water network around poly(dG).poly(dC) and poly(dG-dC).poly(dG-dC) on the basis of Monte Carlo computations

A computational method is elaborated for studying the water environment around regular polynucleotide duplexes; it allows rigorous structural information on the hydration shell of DNA to be obtained. The crucial aspect of this Monte Carlo simulation is the use of periodical boundary conditions. The output data consists of local maxima of water density in the space near the DNA molecule and the properties of one- and two-membered water bridges as function of pairs of polar groups of DNA. In the present paper the results for poly(dG).poly(dC) and poly(dG-dC).poly(dG-dC) are presented. The differences in their hydration shells are of a purely structural nature and are caused by the symmetry of the polar groups of the polymers under study, the symmetry being reflected by the hydration shell. The homopolymer duplex hydration shell mirrors the mononucleotide repeat. The water molecules contacting the polynucleotide in the minor groove are located nearly in the plane midway between the planes of successive base pairs. One water molecule per base pair forms a water bridge facing two polar groups of bases from adjacent base pairs and on different strands making a "spine"-like structure. In contrast, the major groove hydration is stabilized exclusively by two-membered water bridges; the water molecules deepest in the groove are concentrated near the plane of the corresponding base pair. The alternating polymer is characterized by a marked dyad symmetry of the hydration shell corresponding to the axis between two successive base pairs. The minor groove hydration of the dCpdG step resembles the characteristic features of the homopolymer, but the bridge between the O2 oxygens of the other base-stacking type is formed by two water molecules. The major groove hydration is characterized by high probability of one-membered water bridges and by localization of a water molecule on the dyad axis of the dGpdC step. The found structural elements are discussed as reasonable invariants of a dynamic hydration shell.     

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.     

Hydration of counterions interacting with DNA double helix: a molecular dynamics study

In the present work, molecular dynamics simulations have been carried out to study the dependence of counterion distribution around the DNA double helix on the character of ion hydration. The simulated systems consisted of DNA fragment d(CGCGAATTCGCG) in water solution with the counterions Na⁺, K⁺, Cs⁺ or Mg²⁺. The characteristic binding sites of the counterions with DNA and the changes in their hydration shell have been determined. The results show that due to the interaction with DNA at least two hydration shells of the counterions undergo changes. The first hydration shell of Na⁺, K⁺, Cs⁺, and Mg²⁺ counterions in the bulk consists of six, seven, ten, and six water molecules, respectively, while the second one has several times higher values. The Mg²⁺ and Na⁺ counterions, constraining water molecules of the first hydration shell, mostly form with DNA water-mediated contacts. In this case the coordination numbers of the first hydration shell do not change, while the coordination numbers of the second one decrease about twofold. The Cs⁺ and K⁺ counterions that do not constrain surrounding water molecules may be easily dehydrated, and when interacting with DNA their first hydration shell may be decreased by three and five water molecules, respectively. Due to the dehydration effect, these counterions can squeeze through the hydration shell of DNA to the bottom of the double helix grooves. The character of ion hydration establishes the correlation between the coordination numbers of the first and the second hydration shells.

Open image in new window

Graphical Abstract Hydration of counterions interacting with DNA double helix

     

First-principles study of ammonium ions and their hydration in montmorillonites

Density functional theory calculations were performed to investigate the adsorption and hydration of an ammonium ion (NH₄ ⁺) confined in the interlayer space of montmorillonites (MMT). NH₄ ⁺ is trapped in the six-oxygen-ring on the internal surface and forms a strong binding with the surface O atoms. The hydration of NH₄ ⁺ is affected significantly by the surface. Water molecules prefer the surface sites, and do not bind with the NH₄ ⁺ unless enough water molecules are supplied. Moreover, the water molecules involved in NH₄ ⁺ hydration tend to bind with the surface simultaneously. The hydration energy increases with the intercalated water molecules, in contrast to that in gas phase. In addition, the hydration leads to the extension of MMT basal spacing.

The self-diffusion of water in Artemia cysts

The self-diffusion coefficient of water molecules has been measured by nuclear magnetic resonance in cysts of Artemia over a wide range of hydration. Compared to the value for bulk water, the diffusion coefficient is reduced by a factor of 7 at the highest hydration and by approximately 100 at the lowest hydration. The results are used to evaluate various models that have been proposed to account for the reduction of water self-diffusion coefficients in complex systems.     

Large-scale networks of hydration water molecules around bovine beta-trypsin revealed by cryogenic X-ray crystal structure analysis.

The hydration structure of bovine beta-trypsin was investigated in cryogenic X-ray diffraction experiments. Three crystal forms of the enzyme inhibited by benzamidine with different molecular packing were selected to deduce the hydration structure for the entire surface of the enzyme. The crystal structures in all three of the crystal forms were refined at the resolution of 1.8 A at 100 K and 293 K. The number of hydration water molecules around the enzyme at 100 K was 1.5 to two times larger than that at 293 K, indicating that the motion of hydration water was quenched by cooling. In particular, the increase in the number of hydration water molecules was prominent on flat and electrostatically neutral surface areas. The water-to-protein mass ratio and the radius of gyration of a structural model of hydrated trypsin at 100 K was consistent with the results obtained by other experimental techniques for proteins in solution. Hydration water molecules formed aggregates of various shapes and dimensions, and some of the aggregates even covered hydrophobic residues by forming oligomeric arrangements. In addition, the aggregates brought about large-scale networks of hydrogen bonds. The networks covered a large proportion of the surface of trypsin like a patchwork, and mechanically linked several secondary structures of the enzyme. By merging the hydration structures of the three crystal forms at 100 K, a distribution function of hydration water molecules was introduced to approximate the static hydration structure of trypsin in solution. The function showed that the negatively charged active site of trypsin tended to be easily exposed to bulk solvent. This result is of interest with respect to the solvent shielding effect and the recognition of a positively charged substrate by trypsin.     

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.     

Effect of hydration on the structure of non aqueous ethyl cellulose/propylene glycol dicaprylate gels

Changes in the structural properties of ethyl cellulose/propylene glycol dicaprylate systems (EC/PGD), intended for topical drug delivery, upon addition of water were investigated. Although designed to be a non-aqueous vehicle for moisture sensitive drugs, these systems are expected to experience an aqueous environment during production, storage and application on the skin. Hence, the interaction of water molecules with the non aqueous gel system and their distribution within the gel network is of interest and critical to its application. Experimental techniques of this study were small-deformation dynamic oscillation in shear, modulated differential scanning calorimetry (MDSC), (2)H NMR spectroscopy, ATR-infrared spectroscopy, wide-angle X-ray diffraction patterns and light microscopy. Rheological profiles of the gels containing moisture from 0.1 to 40.0% (w/w) deviated considerably from that of the non aqueous system at levels of water above 10.0% in preparations. Gradual replacement of the EC/PGD dipole interactions with stronger hydrogen bonding between ethyl cellulose chains, as the level of hydration increased, contributed to these observations. Formation of clusters of ethyl cellulose, observed under a light microscope, was thus ensued. X-ray diffraction patterns showed that the rearrangement of the polymer chains led to the loss of liquid crystal structures found in the anhydrous gel. MDSC and (2)H NMR were used to further shed light on the thermodynamic state of added water molecules in the gels. Plots of enthalpy obtained calorimetrically and a good correlation between MDSC and (2)H NMR data indicate that gels with less than two percent hydration contain water in a non-freezable bound state, whereas freezable moieties are obtained at levels of hydration above five percent in composite (EC/PGD/water) gels.     

Role of flexibility and polarity as determinants of the hydration of internal cavities and pockets in proteins

Molecular dynamics simulations of Staphylococcal nuclease and of 10 variants with internal polar or ionizable groups were performed to investigate systematically the molecular determinants of hydration of internal cavities and pockets in proteins. In contrast to apolar cavities in rigid carbon structures, such as nanotubes or buckeyballs, internal cavities in proteins that are large enough to house a few water molecules will most likely be dehydrated unless they contain a source of polarity. The water content in the protein interior can be modulated by the flexibility of protein elements that interact with water, which can impart positional disorder to water molecules, or bias the pattern of internal hydration that is stabilized. This might explain differences in the patterns of hydration observed in crystal structures obtained at cryogenic and room temperature conditions. The ability of molecular dynamics simulations to determine the most likely sites of water binding in internal pockets and cavities depends on its efficiency in sampling the hydration of internal sites and alternative protein and water conformations. This can be enhanced significantly by performing multiple molecular dynamics simulations as well as simulations started from different initial hydration states.     

State of water and its diffusion across lipid bilayers: effect of hydration degree

The state of adsorbed water (estimated from the dependence of the shape of the 1H NMR spectrum on the angle between the normal to the bilayers and the direction of the magnetic field) and the diffusion of water molecules in the direction of the normal to the bilayers (estimated by 1H NMR spectroscopy with the impulse gradient of magnetic field) in microscopically oriented dioleoylphosphatidylcholine bilayers have been studied depending on hydration. The dependences of the shape of the NMR spectrum on angle differ qualitatively only at concentrations of water greater and less than the concentration that is achieved upon hydration from saturated vapors chi(eq) (about 23 weight %). At concentrations below chi(eq), all water present in samples enters the hydrate shells of polar "heads" of lipids or is in the state of "rapid exchange" with the water of hydrate shells, with the result that the signal of spin echo for water is observed only in a narrow range of angles close to the "magic angle", 54 degrees C. At concentrations above xhi(eq), the signal of spin echo for water is retained at all orientations, indicating probably that part of water between the bilayers ("quasi-free water") is in the state of a "slow exchange" with water "bound" to polar "heads". It was found that the coefficient of self-diffusion of water across the system of bilayers inversely depends on the degree of hydration, which is described in the Tanner model with consideration of the self-diffusion of water molecules in the hydrophobic moiety of the bilayer. The permeability of the bilayer, the coefficient of distribution of molecules between the water and lipid phases, and the coefficient of self-diffusion of water in the hydrophobic moiety of the bilayer were estimated.     

Hydration of nucleic acid bases studied using novel atom-atom potential functions

New simple atom-atom potential functions for simulating behavior of nucleic acids and their fragments in aqueous solutions are suggested. These functions contains terms which are inversely proportional to the first (electrostatics), sixth (or tenth for the atoms, forming hydrogen bonds) and twelfth (repulsion of all the atoms) powers of interatomic distance. For the refinement of the potential function parameters calculations of ice lattice energy, potential energy and configuration of small clusters consisting of water and nucleic acid base molecules as well as Monte Carlo simulation of liquid water were performed. Calculations using new potential functions give rise to more linear hydrogen bonds between water and base molecules than using other potentials. Sites of preferential hydration of five nucleic bases - uracil, thymine, cytosine, guanine and adenine as well as of 6,6,9-trimethyladenine were found. In the most energetically favourable sites water molecular interacts with two adjacent hydrophilic centres of the base. Studies of interaction of the bases with several water molecules showed that water-water interactions play an important role in the arrangement of the nearest to the base water molecules. Hydrophilic centres are connected by "bridges" formed by hydrogen bonded water molecules. The results obtained are consistent with crystallographic and mass-spectrometric data.     

Monte Carlo simulation of hydration of the guanine-uracil pairs with guanine in two tautomeric forms: contribution of water bridging to relative stability of mispairs.

Results are presented from Monte Carlo simulation of hydration of guanine-uracil mispairs by 25 and 50 water molecules. The hydration shells of three mispairs formed between "normal" dioxo form of uracil (U) and three forms of guanine ("normal" amino-oxo tautomer G and two rotamers of the "rare" amino-hydroxy tautomer G*) depend on the tautomeric forms of the guanine molecule. The simulation shows the important role of hydration effects on the relative stability of the mispairs.     

Hydration properties of DNA-lysine gels by microwave dielectric measurements as a function of temperature

Dielectric measurements by a cavity perturbation method at 10 GHz in the temperature range from-20°C to +45°C are reported for aqueous gels of herring sperm DNA in the presence of 1 or 3 lysine molecules per nucleotide. Measurements for lysine-water and DNA-water systems are also reported. The experimental results can be accounted for by the presence of interfacial water, with dielectric properties different from those of bulk water, and are analyzed in terms of a three component equation (solute molecules, interfacial water and bulk water) to calculate hydration parameters of the systems. The lysine molecule is found to coordinate a particular number of water molecules, in agreement with the literature. The specific hydration of DNA is reduced by the presence of lysine, indicating a direct interaction between the polyion and the aminoacid: a decrease to about 50% was observed at a ratio of one molecule of lysine per nucleotide. A suggestion is made that the interaction is mainly electrostatic in nature.     

Some ways of looking at compensatory kosmotropes and different water environments

Hydration of macromolecular structures determines biological activity. Stabilizing solutes are kosmotropic (increase order of water) rather than chaotropic (decrease order). Preferential hydration of surfaces is a thermodynamic consequence of the solution behavior of kosmotropic solutes, but inconsistencies imply interactions such as the hydration of specific sites within macromolecules. Thermodynamic measures require bulk pure solutes; here simpler measures of the effects on bulk water, water at surfaces and hydration water of probes have been applied to solutes including natural stabilizers, analogues and example chaotropes. Changes in the near-infrared spectra, water proton NMR chemical shifts and relaxation times measure changes in the bulk liquid; HPLC-column retention of solutes indicate interactions with hydration water at different surfaces, and fluorescence probes detect effects on functional group hydration water. Ab initio calculations and Monte-Carlo simulations of the solutes in water measure the energetics of the solute-water interactions, the dipole moments of these molecules, their charge distributions and the effect of the solute molecules on the structure of water. The rankings of the test solutes by these measures are not consistent. Thus, stabilizing solutes are not interchangeable in biological systems and the intracellular replacement of one by another could affect the integration of cell metabolism.     

Hydration of RNA base pairs

The hydration patterns around the RNA Watson-Crick and non-Watson-Crick base pairs in crystals are analyzed and described. The results indicate that (i) the base pair hydration is mostly "in-plane"; (ii) eight hydration sites surround the Watson-Crick G-C and A-U base pairs, with five in the deep and three in the shallow groove, an observation which extends the characteristic isostericity of Watson-Crick pairs; (iii) while the hydration around G-C base pairs is well defined, the hydration around A-U base pairs is more diffuse; (iv) the hydration sites close to the phosphate groups are the best defined and the most recurrent ones; (v) a string of water molecules links the two shallow groove 2'-hydroxyl groups, and (vi) the water molecules fit into notches, the size and accessibility of which are almost as important as the number and strength of the hydrophilic groups lining the cavity. Residence times of water molecules at specific hydration sites, inferred from molecular dynamics simulations, are discussed in the light of present data.     

Structured water layers adjacent to biological membranes

Water amid the restricted space of crowded biological macromolecules and at membrane interfaces is essential for cell function, though the structure and function of this "biological water" itself remains poorly defined. The force required to remove strongly bound water is referred to as the hydration force and due to its widespread importance, it has been studied in numerous systems. Here, by using a highly sensitive dynamic atomic force microscope technique in conjunction with a carbon nanotube probe, we reveal a hydration force with an oscillatory profile that reflects the removal of up to five structured water layers from between the probe and biological membrane surface. Further, we find that the hydration force can be modified by changing the membrane fluidity. For 1,2-dipalmitoyl-sn-glycero-3-phosphocholine gel (Lbeta) phase bilayers, each oscillation in the force profile indicates the force required to displace a single layer of water molecules from between the probe and bilayer. In contrast, 1,2-dipalmitoyl-sn-glycero-3-phosphocholine fluid (Lalpha) phase bilayers at 60 degrees C and 1,2-dioleoyl-sn-glycero-3-phosphocholine fluid (Lalpha) phase bilayers at 24 degrees C seriously disrupt the molecular ordering of the water and result predominantly in a monotonic force profile.     

Influence of water clustering on the dynamics of hydration water at the surface of a lysozyme

Dynamics of hydration water at the surface of a lysozyme molecule is studied by computer simulations at various hydration levels in relation with water clustering and percolation transition. Increase of the translational mobility of water molecules at the surface of a rigid lysozyme molecule upon hydration is governed by the water-water interactions. Lysozyme dynamics strongly affect translational motions of water and this dynamic coupling is maximal at hydration levels, corresponding to the formation of a spanning water network. Anomalous diffusion of hydration water does not depend on hydration level up to monolayer coverage and reflects spatial disorder. Rotational dynamics of water molecules show stretched exponential decay at low hydrations. With increasing hydration, we observe appearance of weakly bound water molecules with bulklike rotational dynamics, whose fraction achieves 20-25% at the percolation threshold.     

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.     

Crystallographic study of hydration of an internal cavity in engineered proteins with buried polar or ionizable groups

Although internal water molecules are essential for the structure and function of many proteins, the structural and physical factors that govern internal hydration are poorly understood. We have examined the molecular determinants of internal hydration systematically, by solving the crystal structures of variants of staphylococcal nuclease with Gln-66, Asn-66, and Tyr-66 at cryo (100 K) and room (298 K) temperatures, and comparing them with existing cryo and room temperature structures of variants with Glu-66, Asp-66, Lys-66, Glu-92 or Lys-92 obtained under conditions of pH where the internal ionizable groups are in the neutral state. At cryogenic temperatures the polar moieties of all these internal side chains are hydrated except in the cases of Lys-66 and Lys-92. At room temperature the internal water molecules were observed only in variants with Glu-66 and Tyr-66; water molecules in the other variants are probably present but they are disordered and therefore undetectable crystallographically. Each internal water molecule establishes between 3 and 5 hydrogen bonds with the protein or with other internal water molecules. The strength of interactions between internal polar side chains and water molecules seems to decrease from carboxylic acids to amides to amines. Low temperature, low cavity volume, and the presence of oxygen atoms in the cavity increase the positional stability of internal water molecules. This set of structures and the physical insight they contribute into internal hydration will be useful for the development and benchmarking of computational methods for artificial hydration of pockets, cavities, and active sites in proteins.