Monte Carlo Simulation of Hydration of the Nucleic Acid Fragments

Abstract

Systems containing a base or a base pair and 25 water molecules, as well as a helical stack and 30 water molecules per base pair, have been simulated. Changes in the base hydration shell structure, after the bases have been included into the pair and then into the base pair stack are discussed. Hydration shells of several configurations of the base pair stacks are discussed. Probabilities of formation of the hydrogen-bonded bridges of 1, 2 and 3 water molecules between hydrophilic centres have been estimated. The hydration shell structure was shown to depend on the nature of the base pair and on the stack configuration, while dependence of the global hydration shell characteristics on the stack configuration has been proved to be rather slight. The most typical structural elements of hydration shells, in the glycosidic (minor in B-like conformation) and non-glycosidic (major) grooves, for different configurations of AU and GC stacks, have been found and discussed. The number of hydrogen bonds between water molecules and bases per water molecule was shown to change upon transformation of the stack from A to B configuration. This result is discussed in connection with the reasons for B to A conformational transition and the concept of “water economy”. Hydration shell patterns of NH₂-groups of AU and GC helical stacks differ significantly.     

Monte-Carlo simulation of hydration of nucleic acid fragments

Monte-Carlo simulation of the systems containing a stack of 6 complementary base pairs and 180 water molecules has been performed. Characteristic of the hydration shell structure in major and minor grooves has been found for the stacks of repeating A : U and G : C base pairs as well as alternating (A : U, U : A) and (G : C, C : G) ones. Probabilities of the formation of bridges, formed by 1, 2 and 3 water molecules, between hydrophilic centres of the bases have been estimated. One water molecule forms an H-bonded bridge between two adjacent hydrophilic centres with high probability if N...N, N...O or O...O distance between these centres is close to 4.3 A. Hydration shell structure was found to depend significantly on the stack sequence and configuration, while global hydration characteristics (average energy, the number of water-water and water-base H-bonds) are only slightly dependent on the stack sequence and configuration. For the stacks in A conformation the number of water molecules forming more than one H-bonds with the bases is greater in comparison with the stacks in B-like conformation. This result is discussed in connection with the concept of hydration economy during B to A transition.     

AT base pairs are less stable than GC base pairs in Z-DNA: The crystal structure of d(m5CGTAm5CG)

Two hexanucleoside pentaphosphates, 5-methyl and 5-bromo cytosine derivatives of d(CpGpTp-ApCpG) have been synthesized, crystallized, and their three-dimensional structure solved. They both form left-handed Z-DNA and the methylated derivative has been refined to 1.2 Å resolution. These are the first crystal Z-DNA structures that contain AT base pairs. The overall form of the molecule is very similar to that of the unmethylated or the fully methylated (dC-dG)₃ hexamer although there are slight changes in base stacking. However, significant differences are found in the hydration of the helical groove. When GC base pairs are present, the helical groove is systematically filled with two water molecules per base pair hydrogen bonded to the bases. Both of these water molecules are not seen in the electron density map in the segments of the helix containing AT base pairs, probably because of solvent disorder. This could be one of the features that makes AT base pairs form Z-DNA less readily than GC base pairs.     

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

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

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.     

Structure of the hydration shells of oligo(dA-dT).oligo(dA-dT) and oligo(dA).oligo(dT) tracts in B-type conformation on the basis of Monte Carlo calculations

Monte Carlo simulations [(N, V, T)-ensemble] were performed for the hydration shell of poly(dA-dT).poly(dA-dT) in canonical B form and for the hydration shell of poly(dA).poly(dT) in canonical B conformation and in a conformation with narrow minor groove, highly inclined bases, but with a nearly zero-inclined base pair plane (B' conformation). We introduced helical periodic boundary conditions with a rather small unit cell and a limited number of water molecules to reduce the dimensionality of the configuration space. The coordinates of local maxima of water density and the properties of one- and two-membered water bridges between polar groups of the DNA were obtained. The AT-alternating duplex hydration mirrors the dyad symmetry of polar group distribution. At the dApdT step, a water bridge between the two carbonyl oxygens O2 of thymines is formed as in the central base-pair step of Dickerson's dodecamer. In the major groove, 5-membered water chains along the tetranucleotide pattern d(TATA).d(TATA) are observed. The hydration geometry of poly(dA).poly(dT) in canonical B conformation is distinguished by autonomous primary hydration of the base-pair edges in both grooves. When this polymer adopts a conformation with highly inclined bases and narrow minor groove, the water density distribution in the minor groove is in excellent agreement with Dickerson's spine model. One local maximum per base pair of the first layer is located near the dyad axis between adjacent base pairs, and one local maximum per base pair in the second shell lies near the dyad axis of the base pair itself. The water bridge between the two strands formed within the first layer was observed with high probability. But the water molecules of the second layer do not have a statistically favored orientation necessary for bridging first layer waters. In the major groove, the hydration geometry of the (A.T) base-pair edge resembles the main features of the AT-pair hydration derived from other sequences for the canonical B form. The preference of the B' conformation for oligo(dA).oligo(dT) tracts may express the tendency to common hydration of base-pair edges of successive base pairs in the grooves of B-type DNA. The mean potential energy of hydration of canonical B-DNA was estimated to be -60 to -80 kJ/mole nucleotides in dependence on the (G.C) contents. Because of the small system size, this estimation is preliminary.     

Interactions of hydrated IIa and IIb group metal cations with thioguanine-cytosine DNA base pair: Ab initio and density functional theory investigation of polarization effects, differences among cations, and flexibility of the cation hydration shell.

The structures and energies of the thioguanine-cytosine Watson-Crick (thioGC WC) base pair interacting with hydrated IIa (Mg2+, Ca2+, Ba2+) and IIb group (Zn2+, Cd2+, Hg2+) cations have been studied using ab initio techniques. Furthermore, complexes between guanine and thioguanine with hydrated cations have been characterized assuming various structures of the hydration shells. The complexes of the thioGC WC base pair with hydrated cations have similar properties as the previously studied GC WC base pair. There is substantial polarization stabilization of the base pairing due to cation binding which amounts to 7 - 11 kcal/mol. Soft Cd2+ and Hg2+ cations have a uniquely strong interaction with the thiogroup and induce substantial nonplanarity of the pairing. The thiogroup tends to reduce the number of water molecules in the first hydration shell of the cation. All complexes were optimized within the Hartree-Fock (HF) approximation while their energetics has been evaluated using the second-order Moller-Plesset perturbational method (MP2). All interaction energy evaluations and a substantial portion of the optimizations of the hydrated cation-(thio)guanine complexes have been repeated using Becke-3LYP Density Functional Theory method. All three approximations used (HF, Becke-3LYP, and MP2) give qualitatively the same results for the present cationic complexes. The results demonstrate specific differences among the cations and provide a set of reference structures and energies for verification and/or parametrization of empirical potentials and other theoretical methods.     

Structural principles of B-DNA grooves hydration in fibers as revealed by Monte Carlo simulations and X-ray diffraction

Being interested in possible effects of sequence-dependent hydration of B-DNA with mixed sequence in fibers, we performed a series of Monte Carlo calculations of hydration of polydeoxyribonucleotides in B form, considering all sequences with dinucleotide repeat. The computational results allow the ten base-stacking types to be classified in accordance with their primary hydration in the minor groove. As a rule, the minor groove is occupied by two water molecules per base pair in the depth of the groove, which are located nearly midway between the planes of successive base pairs and symmetrically according to the dyad there. The primary hydration of the major groove depends on the type of the given base pair. The coordinates of 3 water molecules per base pair in the depth of the major groove are determined by the type of this pair together with its position and orientation in the helix, and are practically independent on the adjacent base pairs. A/T-homopolymer tracts do not fit into this hydration pattern; the base pair edges are hydrated autonomously in both grooves. Analysis of the Li-B-DNA x-ray diffraction intensities reveals those two water positions in the minor groove. In the major groove, no electronic density peaks in sufficient distance from the base edges were found, thus confirming the absence of any helical invariance of primary hydration in this region. With the help of the rules proposed in this paper it is possible to position the water molecules of the first hydration shell in the grooves of canonical B-DNA for any given sequence.     

Structure of a T4 hairpin loop on a Z-DNA stem and comparison with A-RNA and B-DNA loops

The synthetic DNA oligomer C-G-C-G-C-G-T-T-T-T-C-G-C-G-C-G crystallizes as a Z-DNA hexamer, capped at one end by a T4 loop. The crystals are monoclinic, space group C2, with a = 57.18 A, b = 21.63 A, c = 36.40 A, beta = 95.22 degrees, and one hairpin molecule per asymmetric unit. The structure of the z-hexamer stem was determined by molecular replacement, and the T4 loop was positioned by difference map methods. The final R factor at 2.1 A resolution for hairpin plus 70 water molecules is 20% for 2 sigma data, with a root-mean-square error of 0.26 A. The (C-G)3 stem resembles the free Z-DNA hexamer with minor crystal packing effects. The T4 loop differs from that observed on a B-DNA stem in solution, or in longer loops in tRNA, in that it shows intraloop and intermolecular interactions rather than base stacking on the final base-pair of the stem. Bases T7, T8 and T9 stack with one another and with the sugar of T7. Two T10 bases from different molecules stack between the C1-G12 terminal base-pairs of a third and fourth molecule, to simulate a T.T "base-pair". Distances between thymine N and O atoms suggest that the two thymine bases are hydrogen bonded, and a keto-enol tautomer pair is favored over disordered keto-keto wobble pairs. The hairpin molecules pack in the crystal in herringbone columns in a manner that accounts well for the observed relative crystal growth rates in a, b and c directions. Hydration seems to be most extensive around the phosphate groups, with lesser hydration within the grooves.     

Modeling of hydration of incorrect nucleic acid base pairs by the Monte Carlo method

The hydration of water bridged base pairs of nucleic acids have been simulated via the Monte Carlo method. The simulation have shown that water molecules forming H-bonds with both bases preserve this H-bonding with large probability in the water surrounding. This fact supports the supposition about the important role of water molecules in wrong base pair formation and about the role of these base pairs in the structure and functioning of nucleic acids.     

Crystal structure of 2'-O-Me(CGCGCG)2, an RNA duplex at 1.30 A resolution. Hydration pattern of 2'-O-methylated RNA.

The molecular and crystal structure of 2'-O-Me (CGCGCG)2 has been determined using synchrotron radiation at near-atomic resolution (1.30 A), the highest resolution to date in the RNA field. The crystal structure is a half-turn A-type helix with some helical parameters deviating from canonical A-RNA, such as low base pair rise, elevated helical twist and inclination angles. In CG steps, inter-strand guanines are parallel while cytosines are not parallel. In steps GC this motif is reversed. This type of regularity is not seen in other RNA crystal structures. The structure includes 44 water molecules and two hydrated Mg2+ions one of which lies exactly on the crystallographic 2-fold axis. There are distinct patterns of hydration in the major and the minor grooves. The major groove is stabilised by water clusters consisting of fused five- and six-membered rings. Minor groove contains only a single row of water molecules; each water bridges either two self-parallel cytosines or two self-parallel guanines by a pair of hydrogen bonds. The structure provides the first view of the hydration scheme of 2'-O-methylated RNA duplex.     

Dependence of the hydration shell structure in the minor groove of the DNA double helix on the groove width as revealed by Monte Carlo simulation.

The hydration shell of several conformations of the polynucleotides poly(dA).poly(dT), poly(dA).poly(dU), and poly(dA-dI).poly(dT-dC) has been simulated using the Monte Carlo method (Metropolis sampling). Calculations have shown that the structure of the hydration shell of the minor groove greatly depends on its width. In conformations with a narrowed minor groove, the first layer of the hydration shell of this groove has only one molecule per nucleotide pair that forms H bonds with purine N3 of one pair and pyrimidine O2 of the next pair. The second layer of the hydration shell of such conformations contains molecules that form H bonds between two adjacent molecules of the first layer. The probability of formation of hydration spine is about 20% while the bridges of the first layer are formed with a probability of about 70%. In the first layer of the minor groove of the B-DNA conformation with wide minor groove there are approximately two water molecules per base pair that form H bonds with purine N3 or pyrimidine O2 and with the sugar ring oxygen of the adjacent nucleotide. The probability of simultaneous H bonding of a water molecule with N3 (or O2) and O of sugar ring is about 30%. The results of simulation suggest that hydration spine proposed for the narrowed minor groove of oligonucleotide crystals [H. R. Drew, and R. E. Dickerson (1981) Journal of Molecular Biology, Vol. 151, pp. 535-556] can be formed in fibers of poly(dA).poly(dT), poly(dA).poly(dU), and poly(dA-dI).poly(dT-dC) as well as in DNA fragments of these sequences in solution.(ABSTRACT TRUNCATED AT 250 WORDS)     

Intercalation-mediated synthesis and replication: a new approach to the origin of life

We propose that a molecular midwife, a flat molecule approximately 10 Ax10 A with two hydrophobic faces, was essential to the origin of life. This molecule was positively charged, water soluble and did not strongly associate with itself in solution. It may have been a derivative of phthalocyanine that no longer exists on the Earth today, and might have been formed solely from hydrogen cyanide and formaldehyde. The midwife tended to intercalate between side groups (bases, similar to those in RNA) of polymers to form stacks, which incorporated bare bases. The midwife alternated in these stacks with hydrogen-bonded tetrads of bases. Under conditions of low water activity, as in a desert during the day, bare bases in the stacks were joined together by neutral and chemically heterogeneous backbones of no fixed chirality. The components of the backbones were the products of the formose reaction of formaldehyde, and were involved in the reversible formation of N -glycosides and acetals catalysed by divalent metal ions. The final product of this assemblage was a fully intercalated quadruplex of four information-containing polymer strands (four proto -RNA molecules). This process constituted replication of the original polymer that had seeded the formation of the stack. The stack structure ensured that the polymer's base sequence was replicated faithfully despite the lack of both homochirality and chemical homogeneity in the backbone. At night, water from condensing dew would suddenly come in contact with these products, quenching all chemical reactions and releasing midwife molecules and single- or double-stranded proto-RNA. Evaporation of water during the day then gave new stacks containing one or two proto-RNA strands, bare bases, and midwife molecules, which could begin a new replication cycle. Our model also allows for the generation of new stacks and the extension of existing ones, without restricting the base sequence of either, thereby providing a source of genetic information. The proto-RNA replication cycle is driven purely by concentration changes caused by the Sun and the rotation of the Earth. We propose that this system as a whole could have gradually evolved into the RNA World.     

Modeling high-resolution hydration patterns in correlation with DNA sequence and conformation

Hydration around the DNA fragment d(C5T5).(A5G5) is presented from two molecular dynamics simulations of 10 and 12 ns total simulation time. The DNA has been simulated as a flexible molecule with both the CHARMM and AMBER force fields in explicit solvent including counterions and 0.8 M additional NaCl salt. From the previous analysis of the DNA structure B-DNA conformations were found with the AMBER force-field and A-DNA conformations with CHARMM parameters. High-resolution hydration patterns are compared between the two conformations and between C.G and T.A base-pairs from the homopolymeric parts of the simulated sequence. Crystallographic results from a statistical analysis of hydration sites around DNA crystal structures compare very well with the simulation results. Differences between the crystal sites and our data are explained by variations in conformation, sequence, and limitations in the resolution of water sites by crystal diffraction. Hydration layers are defined from radial distribution functions and compared with experimental results. Excellent agreement is found when the measured experimental quantities are compared with the equivalent distribution of water molecules in the first hydration shell. The number of water molecules bound to DNA was found smaller around T.A base-pairs and around A-DNA as compared to B-DNA. This is partially offset by a larger number of water molecules in hydrophobic contact with DNA around T.A base-pairs and around A-DNA. The numbers of water molecules in minor and major grooves have been correlated with helical roll, twist, and inclination angles. The data more fully explain the observed B-->A transition at low humidity.     

Structure and hydration of BamHI DNA recognition site: a molecular dynamics investigation

The results of a 3-ns molecular dynamics simulation of the dodecamer duplex d(TATGGATCCATA)(2) recognized by the BamHI endonuclease are presented here. The DNA has been simulated as a flexible molecule using an AMBER force field and the Ewald summation method, which eliminates the undesired effects of truncation and permits evaluation of the full effects of electrostatic forces. The starting B conformation evolves toward a configuration quite close to that observed through x-ray diffraction in its complex with BamHI. This configuration is fairly stable and the Watson-Crick hydrogen bonds are well maintained over the simulation trajectory. Hydration analysis indicates a preferential hydration for the phosphate rather than for the ester oxygens. Hydration shells in both the major and minor groove were observed. In both grooves the C-G pairs were found to be more hydrated than A-T pairs. The "spine of hydration" in the minor groove was clear. Water residence times are longer in the minor groove than in the major groove, although relatively short in both cases. No special long values are observed for sites where water molecules were observed by x-ray diffraction, indicating that water molecules having a high probability of being located in a specific site are also fast-exchanging.     

Modeling DNA Hydration: Comparison of Calculated and Experimental Hydration Properties of Nuclic Acid Bases

Abstract

Hydration properties of individual nucleic acid bases were calculated and compared with the available experimental data. Three sets of classical potential functions (PF) used in simulations of nucleic acid hydration were juxtaposed: (i) the PF developed by Poltev and Malenkov (PM), (ii) the PF of Weiner and Kollman (WK), which together with Jorgensen's TIP3P water model are widely used in the AMBER program, and (HI) OPLS (optimized potentials for liquid simulations) developed by Jorgensen (J). The global minima of interaction energy of single water molecules with all the natural nucleic acid bases correspond to the formation of two water-base hydrogen bonds (water bridging of two hydrophilic atoms of the base). The energy values of these minima calculated via PM potentials are in somewhat better conformity with mass-spectrometric data than the values calculated via WK PF. OPLS gave much weaker water-base interactions for all compounds considered, thus these PF were not used in further computations. Monte Carlo simulations of the hydration of 9- methyladenine, 1-methyluracil and 1-methylthymine were performed in systems with 400 water molecules and periodic boundary conditions. Results of simulations with PM potentials give better agreement with experimental data on hydration energies than WK PF. Computations with PM PF of the hydration energy of keto and enol tautomers of 9-methyl- guanine can account for the shift in the tautomeric equilibrium of guanine in aqueous media to a dominance of the keto form in spite of nearly equal intrinsic stability of keto and enol tautomers. The results of guanine hydration computations are discussed in relation to mechanisms of base mispairing errors in nucleic acid biosynthesis. The data presented in this paper along with previous results on simulation of hydration shell structures in DNA duplex grooves provide ample evidence for the advantages of PM PF in studies of nucleic-acid hydration.     

The A--B transition: temperature and base composition effects on hydration of DNA

Natural DNAs and some polynucleotides organised in fiber present the A--B form transition at a relative humidity (r.h.) which depends on the temperature. A shift of the midpoint of that helix--helix transition to higher r.h. values is observed when the temperature is risen. It is shown that the average number of water molecules associated to a nucleotide pair is the relevant parameter for the A-B transition and that this parameter can be given a precise value by a combination of different r.h. and temperature values. The minimum number of water molecules necessary to get the B form depends on the base composition of the DNA. It is observed that AT base pairs have a higher affinity toward water molecules than GC base pairs. In the B form there are 27 water molecules per GC nucleotide pair and 44 per AT pair. Moreover, we noted that the fraction of nucleotides in the B form as a function of the average number of water molecules associated per base pair does not depend on the temperature. The A helical form is obtained with about 11 water molecules per nucleotide pair and this number is not very sensitive to the base composition of DNA.     

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

     

Crystal structure of paromomycin docked into the eubacterial ribosomal decoding A site

BACKGROUND: Aminoglycoside antibiotics interfere with translation in both gram-positive and gram-negative bacteria by binding to the tRNA decoding A site of the 16S ribosomal RNA. RESULTS: Crystals of complexes between oligoribonucleotides incorporating the sequence of the ribosomal A site of Escherichia coli and the aminoglycoside paromomycin have been solved at 2.5 A resolution. Each RNA fragment contains two A sites inserted between Watson-Crick pairs. The paromomycin molecules interact in an enlarged deep groove created by two bulging and one unpaired adenines. In both sites, hydroxyl and ammonium side chains of the antibiotic form 13 direct hydrogen bonds to bases and backbone atoms of the A site. In the best-defined site, 8 water molecules mediate 12 other hydrogen bonds between the RNA and the antibiotics. Ring I of paromomycin stacks over base G1491 and forms pseudo-Watson-Crick contacts with A1408. Both the hydroxyl group and one ammonium group of ring II form direct and water-mediated hydrogen bonds to the U1495oU1406 pair. The bulging conformation of the two adenines A1492 and A1493 is stabilized by hydrogen bonds between phosphate oxygens and atoms of rings I and II. The hydrophilic sites of the bulging A1492 and A1493 contact the shallow groove of G=C pairs in a symmetrical complex. CONCLUSIONS: Water molecules participate in the binding specificity by exploiting the antibiotic hydration shell and the typical RNA water hydration patterns. The observed contacts rationalize the protection, mutation, and resistance data. The crystal packing mimics the intermolecular contacts induced by aminoglycoside binding in the ribosome.     

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.