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.     

X-ray crystallographic analysis of the hydration of A- and B-form DNA at atomic resolution

We have determined single crystal structures of an A-DNA decamer and a B-DNA dodecamer at 0.83 and 0.95 A, respectively. The resolution of the former is the highest reported thus far for any right-handed nucleic acid duplex and the quality of the diffraction data allowed determination of the structure with direct methods. The structures reveal unprecedented details of DNA fine structure and hydration; in particular, we have reexamined the overall hydration of A- and B-form DNA, the distribution of water around phosphate groups, and features of the water structure that may underlie the B to A transition.     

An analysis of the relationship between hydration and protein-DNA interactions.

Eleven protein-DNA crystal structures were analyzed to test the hypothesis that hydration sites predicted in the first hydration shell of DNA mark the positions where protein residues hydrogen-bond to DNA. For nine of those structures, protein atoms, which form hydrogen bonds to DNA bases, were found within 1.5 A of the predicted hydration positions in 86% of the interactions. The correspondence of the predicted hydration sites with the hydrogen-bonded protein side chains was significantly higher for bases inside the conserved DNA recognition sequences than outside those regions. In two CAP-DNA complexes, predicted base hydration sites correctly marked 71% (within 1.5 A) of protein atoms, which form hydrogen bonds to DNA bases. Phosphate hydration was compared to actual protein binding sites in one CAP-DNA complex with 78% marked contacts within 2.0 A. These data suggest that hydration sites mark the binding sites at protein-DNA interfaces.     

Solvent bridging sites in A- and B-DNA helices

Nucleotide hydration is important for the understanding of the stability of and the transitions between the different helical conformations of DNA. We have used energy minimization and geometric criteria in order to look for possible sites for solvent which can bridge more than one polar or charged atomic group on a nucleotide. Such bridging sites between phosphate groups have been seen experimentally and used to explain the A to B transition. We show that these phosphate bridging sites occur at energy minima around A-DNA but do not occur around B-DNA. We also find that there are further low energy bridging sites which depend on sequence and which enable the more economical hydration of the A form.     

Water and ions in a high resolution structure of B-DNA.

A detailed picture of hydration and counterion location in the B-DNA duplex d(GCGAATTCG) is presented. Detailed data have been obtained by single crystal x-ray diffraction at atomic resolution (0.89 A) in the presence of Mg(2+). The latter is the highest resolution ever obtained for a B-DNA oligonucleotide. Minor groove hydration is compared with that found in the Na(+) and Ca(2+) crystal forms of the related dodecamer d(CGCGAATTCGCG). High resolution data (1.45 A) of the Ca(2+) form obtained in our laboratory are used for that purpose. The central GAATTC has a very stable hydration spine identical in all cases, independent of duplex length and crystallization conditions (counterions, space group). However, the organization of the water molecules (tertiary and quaternary layers) associated with the central spine vary in each case.     

Sodium and chlorine ions as part of the DNA solvation shell.

The distribution of sodium and chlorine ions around DNA is presented from two molecular dynamics simulations of the DNA fragment d(C(5)T(5)). (A(5)G(5)) in explicit solvent with 0.8 M additional NaCl salt. One simulation was carried out for 10 ns with the CHARMM force field that keeps the DNA structure close to A-DNA, the other for 12 ns with the AMBER force field that preferentially stabilizes B-DNA conformations (, Biophys. J. 75:134-149). From radial distributions of sodium and chlorine ions a primary ion shell is defined. The ion counts and residence times of ions within this shell are compared between conformations and with experiment. Ordered sodium ion sites were found in minor and major grooves around both A and B-DNA conformations. Changes in the surrounding hydration structure are analyzed and implications for the stabilization of A-DNA and B-DNA conformations are discussed.     

A molecular simulation picture of DNA hydration around A- and B-DNA

Recent results from molecular dynamics (MD) simulations on hydration of DNA with respect to conformation are reviewed and compared with experimental data. MD simulations of explicit solvent around DNA can now give a detailed model of DNA that not only matches well with the experimental data but provides additional insight beyond current experimental limitations. Such simulation results are analyzed with a focus on differential hydration properties between A- and B-DNA and between C/G and A/T base pairs. The extent of hydration is determined from the number of waters in the primary shell and compared to experimental numbers from different measurements. High-resolution hydration patterns around the whole DNA are shown and correlated with the conformations. The role of ions associating with DNA is discussed with respect to changes in the hydration structure correlating with DNA conformation.     

Molecular dynamics simulation study of oriented polyamine- and Na-DNA: sequence specific interactions and effects on DNA structure

Molecular dynamics (MD) computer simulations have been carried out on four systems that correspond to an infinite array of parallel ordered B-DNA, mimicking the state in oriented DNA fibers and also being relevant for crystals of B-DNA oligonucleotides. The systems were all comprised of a periodical hexagonal cell with three identical DNA decamers, 15 water molecules per nucleotide, and counterions balancing the DNA charges. The sequence of the double helical DNA decamer was d(5'-ATGCAGTCAG)xd(5'-TGACTGCATC). The counterions were the two natural polyamines spermidine(3+) (Spd(3+)) and putrescine(2+) (Put(2+)), the synthetic polyamine diaminopropane(2+) (DAP(2+)), and the simple monovalent cation Na(+). This work compares the specific structures of the polyamine- and Na-DNA systems and how they are affected by counterion interactions. It also describes sequence-specific hydration and interaction of the cations with DNA. The local DNA structure is dependent on the nature of the counterion. Even the very similar polyamines, Put(2+) and DAP(2+), show clear differences in binding to DNA and in effect on hydration and local structure. Generally, the polyamines disorder the hydration of the DNA around their binding sites whereas Na(+) being bound to DNA attracts and organizes water in its vicinity. Cation binding at the selected sites in the minor and in the major groove is compared for the different polyamines and Na(+). We conclude that the synthetic polyamine (DAP(2+)) binds specifically to several structural and sequence-specific motifs on B-DNA, unlike the natural polyamines, Spd(3+) and Put(2+). This specificity of DAP(2+) compared to the more dynamic behavior of Spd(3+) and Put(2+) may explain why the latter polyamines are naturally occurring in cells.     

B-DNA structure is intrinsically polymorphic: even at the level of base pair positions

Increasingly exact measurement of single crystal X-ray diffraction data offers detailed characterization of DNA conformation, hydration and electrostatics. However, instead of providing a more clear and unambiguous image of DNA, highly accurate diffraction data reveal polymorphism of the DNA atomic positions and conformation and hydration. Here we describe an accurate X-ray structure of B-DNA, painstakingly fit to a multistate model that contains multiple competing positions of most of the backbone and of entire base pairs. Two of ten base-pairs of CCAGGCCTGG are in multiple states distinguished primarily by differences in slide. Similarly, all the surrounding ions are seen to fractionally occupy discrete competing and overlapping sites. And finally, the vast majority of water molecules show strong evidence of multiple competing sites. Conventional resolution appears to give a false sense of homogeneity in conformation and interactions of DNA. In addition, conventional resolution yields an average structure that is not accurate, in that it is different from any of the multiple discrete structures observed at high resolution. Because base pair positional heterogeneity has not always been incorporated into model-building, even some high and ultrahigh-resolution structures of DNA do not indicate the full extent of conformational polymorphism.     

Hydration of the phosphate group in double-helical DNA.

Water distributions around phosphate groups in 59 B-, A-, and Z-DNA crystal structures were analyzed. It is shown that the waters are concentrated in six hydration sites per phosphate and that the positions and occupancies of these sites are dependent on the conformation and type of nucleotide. The patterns of hydration that are characteristic of the backbone of the three DNA helical types can be attributed in part to the interactions of these hydration sites.     

Crystallographic studies of metal ion-DNA interactions: different binding modes of cobalt(II), copper(II) and barium(II) to N7 of guanines in Z-DNA and a drug-DNA complex.

Metal ion coordination to nucleic acids is not only required for charge neutralization, it is also essential for the biological function of nucleic acids. The structural impact of different metal ion coordinations of DNA helices is an open question. We carried out X-ray diffraction analyses of the interactions of the two transition metal ions Co(II) and Cu(II) and an alkaline earth metal ion Ba(II), with DNA of different conformations. In crystals, Co(II) ion binds exclusively at the N7 position of guanine bases by direct coordination. The coordination geometry around Co(II) is octahedral, although some sites have an incomplete hydration shell. The averaged Co-N7 bond distance is 2.3 A. The averaged Co-N7-C8 angle is 121 degrees, significantly smaller than the value of 128 degrees if the Co-N7 vector were to bisect the C5-N7-C8 bond angle. Model building of Co(II) binding to guanine N7 in B-DNA indicates that the coordinated waters in the axial positions would have a van der Waals clash with the neighboring base on the 5' side. In contrast, the major groove of A-DNA does not have enough room to accommodate the entire hydration shell. This suggests that Co(II) binding to either B-DNA or A-DNA may induce significant conformational changes. The Z-DNA structure of Cu(II)-soaked CGCGTG crystal revealed that the Cu(II) ion is bis-coordinated to N7 position of G10 and #G12 (# denotes a symmetry-related position) bases with a trigonal bipyramid geometry, suggesting a possible N7-Cu-N7 crosslinking mechanism. A similar bis-coordination to two guanines has also been seen in the interaction of Cu(II) in m5CGUAm5CG Z-DNA crystal and of Ba(II) with two other Z-DNA crystals.     

Hydration of nucleic acid fragments: comparison of theory and experiment for high-resolution crystal structures of RNA, DNA, and DNA-drug complexes.

A computationally efficient method to describe the organization of water around solvated biomolecules is presented. It is based on a statistical mechanical expression for the water-density distribution in terms of particle correlation functions. The method is applied to analyze the hydration of small nucleic acid molecules in the crystal environment, for which high-resolution x-ray crystal structures have been reported. Results for RNA [r(ApU).r(ApU)] and DNA [d(CpG).d(CpG) in Z form and with parallel strand orientation] and for DNA-drug complexes [d(CpG).d(CpG) with the drug proflavine intercalated] are described. A detailed comparison of theoretical and experimental data shows positional agreement for the experimentally observed water sites. The presented method can be used for refinement of the water structure in x-ray crystallography, hydration analysis of nuclear magnetic resonance structures, and theoretical modeling of biological macromolecules such as molecular docking studies. The speed of the computations allows hydration analyses of molecules of almost arbitrary size (tRNA, protein-nucleic acid complexes, etc.) in the crystal environment and in aqueous solution.     

Helix geometry and hydration in an A-DNA tetramer: IC-C-G-G

The DNA oligomer of sequence IC-C-G-G has been synthesized, and its X-ray crystal structure solved at a resolution of 2.0 A, using anomalous scattering from iodines in phase analysis: 48 cycles of Jack-Levitt restrained least-squares refinement resulted in a residual error of 19.9% over all data, or 16.5% for two-sigma data. Two double-helical tetramers stack in the crystal to form a continuous octamer, except for the two missing phosphate connections across the center. The octamer has a mean helix rotation of 33.7 degrees (10.7 base-pairs per turn), rise of 2.87 A, mean inclination angle of base-pairs of 14 degrees, and mean base-pair propeller twist of +16.3 degrees. Local variations in both helix rotation and base plane roll angles, including those across the center of the octamer, are as predicted from base sequence by sum functions sigma 1 and sigma 2. The three known DNA octamers: IC-C-G-G/IC-C-G-G, G-G-T-A-T-A-C-C and G-G-C-C-G-G-C-C, make up a graded series in this order, with monotonically changing structural parameters. An exhaustive comparison of torsion angle correlations among the known A helices confirms some structural expectations and reveals some new features. 86 water molecules have been located per double-helical IC-C-G-G tetramer (the asymmetric unit), of which 451/2 per tetramer lie within a first hydrogen-bonded shell of hydration. No ordered water structure is observed comparable to the minor groove spine of hydration in B-DNA.     

Mediation of the A/B-DNA helix transition by G-tracts in the crystal structure of duplex CATGGGCCCATG

The crystal structure of the DNA dodecamer duplex CATGGGCCCATG lies on a structural continuum along the transition between A- and B-DNA. The dodecamer possesses the normal vector plot and inclination values typical of B-DNA, but has the crystal packing, helical twist, groove width, sugar pucker, slide and x-displacement values typical of A-DNA. The structure shows highly ordered water structures, such as a double spine of water molecules against each side of the major groove, stabilizing the GC base pairs in an A-like conformation. The different hydration of GC and AT base pairs provides a physical basis for solvent-dependent facilitation of the A↔B helix transition by GC base pairs. Crystal structures of CATGGGCCCATG and other A/B-DNA intermediates support a ‘slide first, roll later’ mechanism for the B→A helix transition. In the distribution of helical parameters in protein–DNA crystal structures, GpG base steps show A-like properties, reflecting their innate predisposition for the A conformation.     

A systematic study of patterns of hydration in nucleic acids:(I) guanine and cytosine

The hydration sites of guanine and cytosine are defined by examination of the crystal structures of bases, nucleosides, nucleotides, and three dinucleoside phosphate salts. The patterns of hydration for two guanine and cytosine containing oligonucleotides are then predicted. The relationship between these structural motifs and thermodynamic parameters is discussed.     

The role of a minor groove spine of hydration in stabilizing poly(dA).poly(dT) against fluctuational interbase H-bond disruption in the premelting temperature regime.

Experimental estimates of the premelting Adenine-Thymine base pair opening probability for some B-DNA sequences are two orders of magnitude smaller than those of other B-DNA sequences. The AT pairs in the sequence with smaller open probability seem to be those that have a well defined spine of hydration in the minor groove. We show that this spine of hydration can significantly enhance the thermal stability of the base pairs to which they are attached. The effect of this spine of hydration coupled with the possible stabilization effect contributed from neighboring GC pairs can explain the differences in the observed AT pair opening probability for different AT containing B-DNA sequences.     

Crystal structure of the B-DNA hexamer d(CTCGAG): model for an A-to-B transition.

The crystal structure of the B-DNA hexamer d(CTCGAG) has been solved at 1.9 A resolution by iterative single isomorphous replacement, using the brominated derivative d(CG5BrCGAG), and refined to an R-factor of 18.6% for 120 nonhydrogen nucleic acid atoms and 32 water molecules. Although the central four base pairs form a typical B-form helix, several parameters suggest a transition to an A-like conformation at the termini. Based on this observation, a B-to-A transition was modeled, maintaining efficient base stacking across the junction. The wide minor groove (approximately 6.9 A) is reminiscent of that in the side-by-side double drug-DNA complexes and hosts a double spine of hydration. The global helix axes of the pseudo-continuous helices are at an acute angle of 60 degrees. The pseudocontinuous stacking is reinforced by the minor groove water structure extending between the two duplexes. The crossover point of two pairs of stacked duplexes is at the stacking junction, unlike that observed in the B-DNA decamers and dodecamers. This arrangement may have implications for the structure of a four-way DNA junction. The duplexes are arranged around a large (approximately 20 A diameter) channel centered on a 6(2) screw axis.     

A Systematic Study of Patterns of Hydration in Nucleic Acids:(I) Guanine and Cytosine

Abstract

The hydration sites of guanine and cytosine are defined by examination of the crystal structures of bases, nucleosides, nucleotides, and three dinucleoside phosphate salts. The patterns of hydration for two guanine and cytosine containing oligonucleotides are then predicted. The relationship between these structural motifs and thermodynamic parameters is discussed.     

Theoretical study of hydration of RNA.

The hydration phenomena of A-RNA double helix and the anticodon loop of transfer RNA have been theoretically investigated using the empirical potential energy functions. The hydration schemes of a model compound of A-RNA and a polynucleotide which has the structure of the anticodon loop of yeast tRNAPhe have been determined, and their stabilization energies produced by the introduction of water in the first hydration shell were calculated by considering the hydrated ones as supermolecules. The results indicate that hydration scheme of A-RNA considerably differs from that of B-DNA and stabilization energy due to hydration of A-RNA is not so great as B-DNA. In the anticodon loop structure, however, stabilizing effect of the bound water molecules upon the structure is significant. From the results, the reason why the structure of RNA remains unchanged with the change of hydration degree while that of DNA is altered was studied.