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.     

Structural equilibrium of DNA represented with different force fields.

We have recently indicated preliminary evidence of different equilibrium average structures with the CHARMM and AMBER force fields in explicit solvent molecular dynamics simulations on the DNA duplex d(C5T5) . d(A5G5) (Feig, M. and B.M. Pettitt, 1997, Experiment vs. Force Fields: DNA conformation from molecular dynamics simulations. J. Phys. Chem. B. (101:7361-7363). This paper presents a detailed comparison of DNA structure and dynamics for both force fields from extended simulation times of 10 ns each. Average structures display an A-DNA base geometry with the CHARMM force field and a base geometry that is intermediate between A- and B-DNA with the AMBER force field. The backbone assumes B form on both strands with the AMBER force field, while the CHARMM force field produces heterogeneous structures with the purine strand in A form and the pyrimidine strand in dynamical equilibrium between A and B conformations. The results compare well with experimental data for the cytosine/guanine part but fail to fully reproduce an overall B conformation in the thymine/adenine tract expected from crystallographic data, particularly with the CHARMM force field. Fluctuations between A and B conformations are observed on the nanosecond time scale in both simulations, particularly with the AMBER force field. Different dynamical behavior during the first 4 ns indicates that convergence times of several nanoseconds are necessary to fully establish a dynamical equilibrium in all structural quantities on the time scale of the simulations presented here.     

Molecular dynamic simulations of environment and sequence dependent DNA conformations: the development of the BMS nucleic acid force field and comparison with experimental results

Molecular dynamic (MD) simulations using the BMS nucleic acid force field produce environment and sequence dependent DNA conformations that closely mimic experimentally derived structures. The parameters were initially developed to reproduce the potential energy surface, as defined by quantum mechanics, for a set of small molecules that can be used as the building blocks for nucleic acid macromolecules (dimethyl phosphate, cyclopentane, tetrahydrofuran, etc.). Then the dihedral parameters were fine tuned using a series of condensed phase MD simulations of DNA and RNA (in zero added salt, 4M NaCl, and 75% ethanol solutions). In the tuning process the free energy surface for each dihedral was derived from the MD ensemble and fitted to the conformational distributions and populations observed in 87 A- and B-DNA x-ray and 17 B-DNA NMR structures. Over 41 nanoseconds of MD simulations are presented which demonstrate that the force field is capable of producing stable trajectories, in the correct environments, of A-DNA, double stranded A-form RNA, B-DNA, Z-DNA, and a netropsin-DNA complex that closely reproduce the experimentally determined and/or canonical DNA conformations. Frequently the MD averaged structure is closer to the experimentally determined structure than to the canonical DNA conformation. MD simulations of A- to B- and B- to A-DNA transitions are also shown. A-DNA simulations in a low salt environment cleanly convert into the B-DNA conformation and converge into the RMS space sampled by a low salt simulation of the same sequence starting from B-DNA. In MD simulations using the BMS force field the B-form of d(GGGCCC)2 in a 75% ethanol solution converts into the A-form. Using the same methodology, parameters, and conditions the A-form of d(AAATTT)2 correctly converts into the B-DNA conformation. These studies demonstrate that the force field is capable of reproducing both environment and sequence dependent DNA structures. The 41 nanoseconds (nsec) of MD simulations presented in this paper paint a global picture which suggests that the DNA structures observed in low salt solutions are largely due to the favorable internal energy brought about by the nearly uniform screening of the DNA electrostatics. While the conformations sampled in high salt or mixed solvent environments occur from selective and asymmetric screening of the phosphate groups and DNA grooves, respectively, brought about by sequence induced ion and solvent packing.     

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.     

Characterization of the base stacking interactions in DNA by means of Lennard-Jones empirical potentials

Three empirical potentials of the Lennard-Jones type taken from literature were used to calculate van der Waals contributions to the base-pair couples stacking energies in B-DNA and A-DNA type double helical conformations. The information obtained can be summarized as follows: (1) Purine-pyrimidine and purine-purine (pyrimidine-pyrimidine in the complementary strand) sequences preferred right-handed helical arrangement, whereas pyrimidine-purine sequences favoured left-handed (C-G) or unwound (T-A) stacking geometry; in the latter case this only held for B- but not A-DNA (the C-G sequence was not studied in A-DNA owing to difficulties (see below) with the G amino group in B-DNA); (2) Positive propeller twist of base-pairs was stable in both B- and A-DNA; the thymine methyl group promoted the propeller and this effect was strongest in the A-T step; (3) Tilt of base pairs occurred around zero in B-DNA and between 15-20 degrees C in A-DNA, in agreement with the experimental observations; (4) Vertical separation of base pairs was optimal within 0.33-0.34 nm for B-DNA and around 0.29 nm for A-DNA using the 9-6 potential. The 12-6 potential gave similar results with B-DNA as the 9-6 potential if, however, base pairs were separated by 0.35-0.36 nm; (5) The calculated effect of the guanine amino group was substantially stronger than expected on the basis of data derived from X-ray diffraction studies of oligonucleotide single crystals; (6) In comparison with the 9-6 potential, the 12-6 potential provided more strict energy minima. In summary, the empirical potentials reproduce, at least semiquantitatively, many but not all DNA properties; this should be taken into account whenever the potentials are used for prediction purposes.     

A note on sequence-dependence of DNA structure

A circular dichroism study was conducted on the solution structure of several different oligonucleotides, whose X-ray structures have been solved. It is suggested that in aqueous solution the oligonucleotides can form structures that maintain geometrical elements which are typical of B-DNA, A-DNA, and their intermediate forms. It is shown that 5'GGATGGGAG:5'CTCCCATCC, which forms an A-DNA helix in the crystal state (McCall et al. 1986), in aqueous solution maintains an A-DNA like structure at temperatures below 10 degrees C. At temperatures between 10 degrees C and 25 degrees C it shows a tendency to form an intermediate structure between A-DNA and B-DNA. Also, it is shown that TFE does not cause a transition from B-DNA to A-DNA helix in short DNA fragments, but instead disrupts the helix.     

The potentially Z-DNA-forming sequence d(GTGTACAC) crystallizes as A-DNA

(GT)n/(CA)n sequences have stimulated much interest because of their frequent occurrence in eukaryotic DNA and their potential for forming the left-handed Z-DNA structure. We here report the X-ray crystal structure of a self-complementary octadeoxynucleotide, d(GTGTACAC), at 2.5 A resolution. The molecule adopts a right-handed double-helical conformation belonging to the A-DNA family. In this alternating purine-pyrimidine DNA minihelix the roll and twist angles show alternations qualitatively consistent with Calladine's rules. The average tilt angle of 9.3 degrees is between the values found in A-DNA (19 degrees) and B-DNA (-6 degrees) fibers. It is envisaged that such intermediate conformations may render diversity to genomic DNA. The base-pair tilt angles and the base-pair displacements from the helix axis are found to be correlated for the known A-DNA double-helical fragments.     

The B- to A-DNA transition and the reorganization of solvent at the DNA surface

DNA geometry depends on relative humidity. Using the CHARMM22 force field to push B-DNA to A-DNA, a molecular dynamics simulation of a mixed-sequence 24-basepair DNA double-stranded oligomer, starting from B-DNA, was carried out to explore both the mechanism of the transition and the evolution of hydration patterns on the surface of DNA. Over the 11-ns trajectory, the transition recapitulates the slide-first, roll-later mechanism, is opposed by DNA electrostatics, and is favored by an increasing amount of condensed sodium ions. Hydration was characterized by counting the hydrogen bonds between water and DNA, and by the number of water bridges linking two DNA atoms. The number of hydrogen bonds between water and DNA remains constant during the transition, but there is a 40% increase in the number of water bridges, in agreement with the principle of economy of hydration. Water bridges emerge as delicate sensors of both structure and dynamics of DNA. Both local flexibility and the frustration of the water network on the surface of DNA probably account for the low populations and short residence times of the bridges, and for the lubricant role of water in ligand-DNA interactions.     

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.     

DNA flexibility as a function of allomorphic conformation and of base sequence.

Systematic theoretical modeling of symmetric DNA oligomers, carried out earlier for the B conformation, is now extended to A-DNA. In contrast to the previous results, it is found that A-DNA shows no multiplicity of low-energy substate conformations. The possibilities of the Jumna algorithm are subsequently applied to studying deformations of the oligomers. Controlled winding and stretching deformations are used to study how the two allomorphs and different base sequences absorb such external stress. The results help explain the internal mechanics of the DNA double helix and the extent to which fine structure influences this behavior. The results point to some differences between the A and B double helices, but also to many similarities. Sequence effects on flexibility are relatively limited compared to their impact on optimal energy conformations. It is also shown that the conformational substates detected for B-DNA oligomers are preserved under deformation, but have little influence on its energetics.     

Crystal structures of A-DNA duplexes

All crystal structures of A-DNA duplexes exhibit a typical crystal packing, with the termini of one molecule abutting the shallow grooves of symmetry related neighbors, while all other forms (B, Z, and RNA) tend to form infinitely stacked helices. The A-DNA arrangement leads to the formation of shallow groove base multiples that have implications for the structure of DNA in compacted states. The characteristic packing leaves big solvent channels, which can be sometimes occupied by B-DNA duplexes. Comparisons of the structures of the same oligomer crystallizing in two different space groups and of different sequences crystallizing in the same space group show that the lattice forces dominate the A-DNA conformation in the crystals, complicating the effort to elucidate the influence of the base sequence on the structures. Nevertheless, in both alternating and nonalternating fragments some sequence effects can still be uncovered. Furthermore, several studies have started to define the minimal sequence changes or chemical modifications that can interconvert the oligomers between different double-helical conformers (A-, B-, and Z-form). Overall, it is seen that the rigid nucleotide principle applies to the oligomeric fragments. Besides the structures of the naked DNAs, their interactions with water, polyamines, and metal ions have attracted considerable attention. There are conserved patterns in the hydration, involving both the grooves and the backbone, which are different from those of B-DNA or Z-DNA. Overall, A-DNA seems to be more economically hydrated than B-DNA, particularly around the sugar-phosphate backbone. Spermine was found to be able to bind exclusively to either of the grooves or to the phosphate groups of the backbone, or exhibit a mixed binding mode. The located metal cations prefer binding to guanine bases and phosphate groups. The only mispairs investigated in A-DNA are the wobble pairs, yielding structural insight into their effects on helix stabilities and hydration. G · T wobble pairs have been determined in various sequence contexts, where they differentially affect the conformations and stableness of the duplexes. The structure of a G · ^m5C base pair, which surprisingly also adopted the wobble conformation, suggests that a similar geometry may transiently exist for G · C pairs. These results from the crystalline state will be compared to the solution state and discussed in relation to their relevance in biology. © 1997 John Wiley & Sons, Inc. Biopoly 44: 45–63, 1997     

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.     

Crystal and solution structures of the B-DNA dodecamer d(CGCAAATTTGCG) probed by Raman spectroscopy: heterogeneity in the crystal structure does not persist in the solution structure

The self-complementary dodecamer d(CGCAAATTTGCG) crystallizes as a double helix of the B form and manifests a Raman spectrum with features not observed in Raman spectra of either DNA solutions or wet DNA fibers. A number of Raman bands are assigned to specific nucleoside sugar and phosphodiester conformations associated with this model B-DNA crystal structure. The Raman bands proposed as markers of the crystalline B-DNA structure are compared and contrasted with previously proposed markers of Z-DNA and A-DNA crystals. The results indicate that the three canonical forms of DNA can be readily distinguished by Raman spectroscopy. However, unlike Z-DNA and A-DNA, which retain their characteristic Raman fingerprints in aqueous solution, the B-DNA Raman spectrum is not completely conserved between crystal and solution states. The Raman spectra reveal greater heterogeneity of nucleoside conformations (sugar puckers) in the DNA molecules of the crystal structure than in those of the solution structure. The results are consistent with conversion of one-third of the dG residues from the C2'-endo/anti conformation in the solution structure to another conformation, deduced to be C1'-exo/anti, in the crystal. The dodecamer crystal also exhibits unusually broad Raman bands at 790 and 820 cm-1, associated with the geometry of the phosphodiester backbone and indicating a wider range of (alpha, zeta) backbone torsion angles in the crystal than in the solution structure. The results suggest that backbone torsion angles in the CGC and GCG sequences, which flank the central AAATTT sequence, are significantly different for crystal and solution structures, the former containing the greater diversity.(ABSTRACT TRUNCATED AT 250 WORDS)     

Investigation of fibrous structures. II. General method for the elucidation of the conformations of the complementary nucleic acids

The problem of deducing the DNA structure with correct base pairing and specific symmetry is formulated in the form of algebraic equations. In this way the number of independent variables determining double-helix conformation can be reduced from six to four. The effectiveness of the method is illustrated by the computations for A-DNA and B-DNA. The method allows one to predict all possible conformations of the complementary nucleic acids.     

The molecular electrostatic potential and steric accessibility of A-DNA.

The molecular electrostatic potential and steric accessibility of A-DNA are computed for base sequences (dG.dC)n and (dA.dT)n. An interpretation of the results in terms of the structure of A-DNA is provided and differences with respect to other forms of DNA, namely B-DNA and Z-DNA, are discussed.     

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.     

Influence of divalent magnesium ion on DNA: molecular dynamics simulation studies

A large amount of experimental evidence is available on the effect of magnesium ions on the structure and stability of DNA double helix. Less is known, however, on how these ions affect the stability and dynamics of the molecule. The static time average pictures from X-ray structures or the quantum chemical energy minimized structures lack understanding of the dynamic DNA–ion interaction. The present work addresses these questions by molecular dynamics simulation studies on two DNA duplexes and their interaction with magnesium ions. Results show typical B-DNA character with occasional excursions to deviated states. We detected expected stability of the duplexes in terms of backbone conformations and base pair parameter by the CHARMM-27 force field. Ion environment analysis shows that Mg²⁺ retains the coordination sphere throughout the simulation with a preference for major groove over minor. An extensive analysis of the influence of the Mg²⁺ ion shows no evidence of the popular predictions of groove width narrowing by dipositive metal ion. The major groove atoms show higher occupancy and residence time compared to minor groove for magnesium, where no such distinction is found for the charge neutralizing Na⁺ ions. The determining factor of Mg²⁺ ion’s choice in DNA binding site evolves as the steric hindrance faced by the bulky hexahydrated cation where wider major groove gets the preference. We have shown that in case of binding of Mg²⁺ to DNA non electrostatic contributions play a major role.

An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:5