Monte Carlo simulation of hydration of the nucleic acid fragments

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 NH2-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.     

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.     

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

     

Effect of the sulphur atom on geometry and spectra of the biomolecule 2-thiouracil and in the WC base pair 2-thiouridine-adenosine. Influence of water in the first hydration shell

The effect of the sulphur atom on 2-thiouracil (2TU) and 2-thiouridine molecules, as compared with uracil and uridine molecules, respectively, was carried out in several environments. The predicted IR spectrum of 2TU in the isolated state was compared with that obtained for uracil molecule and with those reported experimentally in matrix isolation. Its crystal unit cell in the solid state was simulated through a tetramer form using DFT methods for the first time. The calculated Raman spectrum was compared to the experimental ones in the solid state. A linear scaling procedure was used for this task. The first hydration shell was simulated by explicit number of water molecules surrounding 2TU up to 30 and was compared with that obtained in uracil molecule. Water molecules ‘distributed’ around 2TU was preferred over that ‘clustering’, because it can better reproduce the hydration and their effects on different parameters of the molecular structure of 2TU and uracil. The total atomic charges and several calculated thermodynamic parameters were discussed. The effect of the sulphur atom on the Watson-Crick (WC) and reverse WC base pair uridine-adenosine was estimated, and the CP corrected interaction energies were calculated. 2-thiouridine has a weaker WC pair than that with uridine, although its slight higher dipole moment (μ) facilitates the interaction with the water molecules. Several helical parameters were determined.     

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.     

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.     

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.     

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.     

Hydration of Cytidine, 2′-Deoxycytidine and their Phosphate Salts in the Aqueous Solutions. A Molecular Dynamics Computer Simulation Study

Abstract

To compare the hydration pattern of the cytidine (Cyd) and 2′-deoxycytidine (dCyd) in the aqueous solutions at the level of microscopic interactions, Molecular Dynamics (MD) computer simulations have been undertaken. The results indicate that the hydration of the heterocyclic base moiety in cytidine and 2′-deoxycytidine has a hydrophobic character. None of the three potential Watson—;Crick base pair centres hydrogen bonds with the water molecules and the formation of something akin to a clathrate cage structure of water around base moieties of nucleosides in the aqueous solution is suggested. In contrast, the hydration of Cyd and dCyd sugar moieties shows a hydrophilic character and the three-dimensional networks of H-bonds involving all hydrophilic centres are formed differently around the ribose and 2′-deoxyribose. The sugar hydroxyl groups participate in the hydrogen bonding with water both as H-donor and as H-acceptor. Their donor-acceptor abilities have been evaluated and compared. The coordination numbers, the geometrical data of the first hydration shell, and the number of hydrogen bonds have been calculated. The changes in the pattern of hydration with the increased concentration of nucleosides and upon nucleoside protonation are discussed. The analysis of the pairwise interaction energies are also presented.     

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.     

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.     

Solvent-accessible surfaces of nucleic acids

Static solvent-accessible surface areas were calculated for DNA and RNA double helices of varied conformation, composition and sequence, for the single helix of poly(rC), and for a transfer RNA. The results show that for DNA and RNA double helices, two thirds of the water-accessible surface area become buried on double helix formation; phosphate oxygens retain near maximal exposure while the bases are 80% buried. Transfer RNA exposes slightly less surface per residue than does double-helical RNA, despite the presence of several additional “modified” groups, all of which are exposed significantly.When a probe corresponding to a single water molecule is used, both the total and atom type exposures are very similar for A-DNA and B-DNA, although marked differences appear in the major and minor groove exposures between the two conformations. For a given base-pair, the accessible surface area buried upon double-helical stacking is nearly constant (within 5%) for different sequences of neighboring base-pairs.For probes larger than single water molecules, there exist considerable differences in the total and atom type exposures of A-DNA and B-DNA. Conformational transitions between the A-DNA and B-DNA helical forms can thus be related to differences in the accessible areas for “structured” water, or a secondary hydration shell, rather than to interactions with individual water molecules of the primary hydration shell. The base-composition dependence of DNA helical conformation can be explained in terms of the opposing effects of thymine methyl groups of A · T base-pairs and the amino groups of G · C base-pairs upon the solvent within the grooves.The area calculations show that primarily the major groove of B-DNA and the minor groove of A-DNA have sufficient accessible surface area to be recognized by a probe size corresponding to the side-chains of amino acids.     

Nature of the Stacking of Nucleic Acid Bases in Water: a Monte Carlo Simulation

Abstract

The results of a Monte Carlo simulation of the hydration of uracil and thymine molecules, their stacked dimers and hydrogen-bonded base pairs are presented. Simulations have been performed in a cluster approximation. The semiempirical atom-atom potential functions have been used (cluster consisting of 200 water molecules). It has been shown that the stacking interactions of uracil and thymine molecules in water arise mainly due to the increase in the water-water interaction during the transition from monomers to dimer. It has been found out that stacked base associates are more preferable than base pairs in water. This preference is mainly due to the energetically more favourable structure of water around the stack.     

Monte-Carlo Simulation of DNA Duplex Hydration. B and B′ Conformations of Poly(dA) · Poly(dT) Have Different Hydration Shells

Abstract

Monte-Carlo simulation of poly(dA) · poly(dT) hydration by 30 water molecules per nucleotide pair has been performed. Two B-family conformations, both with a 36° helical twist but with different minor groove widths, were considered. One conformation is Arnott's standard B form, the other one is specific for poly(dA) · poly(dT) B′ form with a narrowed minor groove. The mean energies and the mean numbers of water-water and water-DNA hydrogen bonds are close for the two conformations. Nevertheless, the hydration shell of the B' form differs drastically from that of the standard B form. The water arrangement in the minor groove of the B′ form resembles the spine of hydration in the central part of Dickerson's dodecamer d(CGCGAATTCGCG). No such spine is formed in the hydration shell of the usual B form with a wider minor groove. In this conformation water bridges between adenine N3 or thymine O2 and oxygen of the sugar ring of the neighbouring nucleotide along the chain can be formed (“strings” in Dickerson's decamer d(CCAAGATTGG)).     

Transcription of the his3 gene region in Saccharomyces cerevisiae

The dodecamer d(CpGpCpGpApApTpTpCpGpCpG) or C-G-C-G-A-A-T-T-C-G-C-G crystallizes as slightly more than one full turn of right-handed B-DNA. It is surrounded in the crystal by one bound spermine molecule and 72 ordered water molecules, most of which associate with polar N and O atoms at the exposed edges of base-pairs. Hydration within the major groove is principally confined to a monolayer of water molecules associated with exposed N and O groups on the bases, with most association being monodentate. Waters hydrating backbone phosphate oxygens tend not to be ordered, except where they are immobilized by 5-methyl groups from nearby thymines. In contrast, the minor groove is hydrated in an extensive and regular manner, with a zigzag “spine” of first- and second-shell hydration along the floor of the groove serving as a foundation for less-regular outer shells extending beyond the radius of the phosphate backbone. This spine network bridges purine N-3 and pyrimidine O-2 atoms in adjacent base-pairs. It is particularly regular in the A-A-T-T center, and is disrupted at the C-G-C-G ends, in part by the presence of the N-2 amino groups on guanine residues. The minor groove hydration spine may be responsible for the stability of the B form of polymers containing only A · T and I · C base-pairs, and its disruption may explain the ease of transition to the A form of polymers with G · C pairs.     

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.     

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.     

Ordered water structure around a B-DNA dodecamer: A quantitative study

The crystal structure of the double-helical B-DNA dodecamer of sequence C-G-C-G-A-A-T-T-C-G-C-G has been solved and refined independently in three forms: (1) the parent sequence at room temperature; (2) the same sequence at 16 K; and (3) the 9-bromo variant C-G-C-G-A-A-T-T^BrC-G-C-G at 7 °C in 60% (v/v) 2-methyl-2.4-pentanediol. The latter two structures show extensive hydration along the phosphate backbone, a feature that was invisible in the native structure because of high temperature factors (indicating thermal or static disorder) of the backbone atoms. Sixty-five solvent peaks are associated with the phosphate backbone, or an average of three per phosphate group. Nineteen other molecules form a first shell of hydration to base edge N and O atoms within the major groove, and 36 more are found in upper hydration layers. The latter tend to occur in strings or clusters spanning the major groove from one phosphate group to another. A single spermine molecule also spans the major groove. In the minor groove, the zig-zag spine of hydration that we believe to be principally responsible for stabilizing the B form of DNA is found in all three structures. Upper level hydration in the minor groove is relatively sparse, and consists mainly of strings of water molecules extending across the groove, with few contacts to the spine below. Sugar O-1′ atoms are closely associated with water molecules, but these are chiefly molecules in the spine, so the association may reflect the geometry of the minor groove rather than any intrinsic attraction of O-1′ atoms for hydration. The phosphate O-3′ and O-5′ atoms within the backbone chain are least hydrated of all, although no physical or steric impediment seems to exist that would deny access to these oxygen atoms by water molecules.