Structural insights into the effect of hydration and ions on A-tract DNA: a molecular dynamics study

DNA structure is known to be sensitive to hydration and ionic environment. To explore the dynamics, hydration, and ion binding features of A-tract sequences, a 7-ns Molecular dynamics (MD) study has been performed on the dodecamer d(CGCAAATTTGCG)(2). The results suggest that the intrusion of Na(+) ion into the minor groove is a rare event and the structure of this dodecamer is not very sensitive to the location of the sodium ions. The prolonged MD simulation successfully leads to the formation of sequence dependent hydration patterns in the minor groove, often called spine of hydration near the A-rich region and ribbon of hydration near the GC regions. Such sequence dependent differences in the hydration patterns have been seen earlier in the high resolution crystal structure of the Drew-Dickerson sequence, but not reported for the medium resolution structures (2.0 approximately 3.0 A). Several water molecules are also seen in the major groove of the MD simulated structure, though they are not highly ordered over the extended MD. The characteristic narrowing of the minor groove in the A-tract region is seen to precede the formation of the spine of hydration. Finally, the occurrence of cross-strand C2-H2.O2 hydrogen bonds in the minor groove of A-tract sequences is confirmed. These are found to occur even before the narrowing of the minor groove, indicating that such interactions are an intrinsic feature of A-tract sequences.     

An ab initio molecular orbital study on the sequence-dependency of DNA conformation: an evaluation of intra- and inter-strand stacking interaction energy

The various nearest neighbor stacking interaction energies of stacked base pairs in the DNA double helix are calculated for both A- and B-type conformations using an ab initio molecular orbital method. It is demonstrated that the sequence-dependent conformational preference for A- or B-type results from the stacking interaction. In particular, the base sequence showing the highest preference for an A-type conformation is revealed as GC/GC, and the one with the next highest preference, AT/AT; for a B-type conformation, the respective sequences are CG/CG and CA/TG. The overall conformation of a DNA fragment is not determined by these particular sequences only but is influenced by all base pair steps. An intrinsically favorable conformation is predicted from the constituent stacking interaction.     

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.     

The hydration of nucleic acid duplexes as assessed by a combination of volumetric and structural techniques

Using high precision densimetric and ultrasonic measurements, we have determined, at 25°C, the apparent molar volumes ΦV and the apparent molar compressibilities ΦK_S of four nucleic acid duplexes—namely, the DNA duplex, poly(dIdC)poly(dIdC); the RNA duplex, poly(rA)poly(rU); and the two DNA/RNA hybrid duplexes, poly(rA)poly(dT) and poly(dA)poly(rU). Using available fiber diffraction data on these duplexes, we have calculated the molecular volumes as well as the solvent‐accessible surface areas of the constituent charged, polar, and nonpolar atomic groups. We found that the hydration properties of these nucleic acid duplexes do not correlate with the extent and the chemical nature of the solvent‐exposed surfaces, thereby suggesting a more specific set of duplex–water interactions beyond general solvation effects. A comparative analysis of our volumetric data on the four duplexes, in conjunction with available structural information, suggests the following features of duplex hydration: (a) The four duplexes exhibit different degrees of hydration, in the order poly(dIdC)poly(dIdC) > poly(dGdC)poly(dGdC) > poly(dAdT)poly(dAdT) ≈ poly(dA)poly(dT). (b) Repetitive AT and IC sequences within a duplex are solvated beyond general effects by a spine of hydration in the minor groove, with this sequence‐specific water network involving about 8 additional water molecules from the second and, perhaps, even the third hydration layers. (c) Repetitive GC and IC sequences within a duplex are solvated beyond general effects by a “patch of hydration” in the major groove, with this water network involving about 13 additional water molecules from the second and, perhaps, even the third hydration layers. (d) Random sequence, polymeric DNA duplexes, which statistically lack extended regions of repetitive AT, GC, or IC sequences, do not experience such specific enhancements of hydration. Consequently, consistent with our previous observations (T. V. Chalikian, A. P. Sarvazyan, G. E. Plum, and K. J. Breslauer, Biochemistry, 1994, Vol. 33, pp. 2394–2401), duplexes with approximately 50% AT content exhibit the weakest hydration, while an increase or decrease from this AT content causes enhancement of hydration, either due to stronger hydration of the minor groove (an increase in AT content) or due to stronger hydration of the major groove (an increase in GC content). (e) In dilute aqueous solutions, a B‐DNA duplex is more hydrated than an A‐DNA duplex, a volumetric‐based conclusion that is in agreement with previous results obtained on crystals, fibers, and DNA solutions in organic solvent–water mixtures. (f) the A‐like, RNA duplex poly(rA)poly(rU) and the structurally similar A‐like, hybrid duplex poly(rA)poly(dT), exhibit similar hydration properties, while the structurally distinct A‐like, hybrid duplex poly(rA)poly(dT) and non‐A‐like, hybrid duplex poly(dA)poly(rU) exhibit differential hydration properties, consistent with structural features dictating hydration characteristics. We discuss how volumetric characterizations, in conjunction with structural studies, can be used to describe, define, and resolve the general and sequence/conformation‐specific hydration properties of nucleic acid duplexes. © 1999 John Wiley & Sons, Inc. Biopoly 50: 459–471, 1999     

Water and ion binding around RNA and DNA (C,G) oligomers

The dynamics, hydration, and ion-binding features of two duplexes, the A(r(CG)(12)) and the B(d(CG)(12)), in a neutralizing aqueous environment with 0.25 M added KCl have been investigated by molecular dynamics (MD) simulations. The regular repeats of the same C=G base-pair motif have been exploited as a statistical alternative to long MD simulations in order to extend the sampling of the conformational space. The trajectories demonstrate the larger flexibility of DNA compared to RNA helices. This flexibility results in less well defined hydration patterns around the DNA than around the RNA backbone atoms. Yet, 22 hydration sites are clearly characterized around both nucleic acid structures. With additional results from MD simulations, the following hydration scale for C=G pairs can be deduced: A-DNA相似文献   

Kinetic and equilibrium analysis of a threading intercalation mode: DNA sequence and ion effects

The interaction of a symmetric naphthalene diimide with alkylamino substituents at each imide position was investigated with the alternating sequence polymers, poly[d(A-T)]2 and poly[d(G-C)]2. Spectrophotometric binding studies indicate strong binding of the diimide to both sequences although the GC binding constant is 20-25 times larger than the AT binding constant. Analysis of the effects of salt concentration on the binding equilibria shows that the diimide forms two ion pairs in its complex with both polymers as expected for a simple dication. Stopped-flow kinetics experiments demonstrate that the diimide both associates and dissociates from DNA more slowly than classical intercalators with similar binding constants. Analysis of salt concentration effects on dissociation kinetics rate constants (kd) reveals that slopes in log kd versus log [Na+] plots are only approximately half the value obtained for classical dicationic intercalators that have both charged groups in the same groove. These kinetics results support a threading intercalation model, with one charged diimide substituent in each of the DNA grooves rather than with both side chains in the same groove, for the diimide complex with DNA. In the rate-determining step of the mechanism for dissociation of a threading complex only one ion pair is broken; the free side chain can then slide between base pairs to put both diimide side chains in the same groove, and this is followed by rapid full dissociation of the diimide. This sequential release of ion pairs makes the dissociation slope for dicationic threading intercalators more similar to the slope for classical monocationic intercalating ligands.(ABSTRACT TRUNCATED AT 250 WORDS)     

Flanking sequence effects within the pyrimidine triple-helix motif characterized by affinity cleaving.

Nearest neighbor interactions affect the stabilities of triple-helical complexes. Within a pyrimidine triple-helical motif, the relative stabilities of natural base triplets T.AT, C + GC, and G.TA, as well as triplets, D3.TA and D3.CG, containing the nonnatural deoxyribonucleoside 1-(2-deoxy-beta-D-ribofuranosyl)-4-(3-benzamido)phenylimidazole (D3) were characterized by the affinity cleaving method in the context of different flanking triplets (T.AT, T.AT: T.AT, C + GC: C + GC, T.AT: G + GC, C + GC). The to be insensitive to substitutions in either the 3' or 5' directions, while the relative stabilities of triple helices containing C + GC triplets decreased as the number of adjacent C + GC triplets increased. Triple helices incorporating a G.TA interaction were most stable when this triplet was flanked by two T.AT triplets and were adversely affected when a C + GC triplet was placed in the adjacent 5' direction. Similarly, complexes containing D3.TA or D3.CG triplets were most stable when the triplet was flanked by two T.AT triplets but were destabilized when the adjacent 3' neighbor position was occupied with a C + GC triplet. This information regarding sequence composition effects in triple-helix formation establishes a set of guidelines for targeting sequences of double-helical DNA by the pyrimidine triple-helix motif.     

Site-resolved energetics in DNA triple helices containing G*TA and T*CG triads

Recognition of specific sites in double-helical DNA by triplex-forming oligonucleotides has been limited until recently to sites containing homopurine-homopyrimidine sequences. G*TA and T*CG triads, in which TA and CG base pairs are specifically recognized by guanine or by thymine, have now extended this recognition code to DNA target sites of mixed base sequences. In the present work, we have obtained a characterization of the stabilities of G*TA and T*CG triads, and of the effects of these triads upon canonical triads, in triple-helical DNA. The three DNA triplexes investigated are formed by the folding of the 31-mers d(GAAXAGGT(5)CCTYTTCT(5)CTTZTCC) with X = G, T, or C, Y = C, A, or G, and Z = C, G, or T. We have measured the exchange rates of imino protons in each triad of the three triplexes using nuclear magnetic resonance spectroscopy. The exchange rates are used to map the local free energy of structural stabilization in each triplex. The results indicate that the stability of Watson-Crick base pairs in the G*TA and T*CG triads is comparable to that of Watson-Crick base pairs in canonical triads. The presence of G*TA and T*CG triads, however, destabilizes neighboring canonical triads, two or three positions removed from the G*TA/T*CG site. Moreover, the long-range destabilizing effects induced by the T*CG triad are larger than those induced by the G*TA triad. These findings reveal the molecular basis for the lower overall stability of G*TA- and T*CG-containing triplexes.     

Mutagenic target for hydroxyl radicals generated in Escherichia coli mutant deficient in Mn- and Fe- superoxide dismutases and Fur,a repressor for iron-uptake systems

We previously reported that mutations in Mn- and Fe-superoxide dismutases and Fur, a repressor for iron uptake systems, simulated generation of hydroxyl radicals, and caused hypermutability in Escherichia coli. The predominant type of spontaneous mutation was GC --> TA, followed by AT --> CG, suggesting the involvement of 7,8-dihydro-8-oxoguanine (8-oxoG) and 1,2-dihydro-2-oxoadenine (2-oxoA) in DNA as well as 7,8-dihydro-8-oxodeoxyguanosine triphosphate (8-oxodGTP) and 1,2-dihydro-2-oxodeoxyadenosine triphosphate (2-oxodATP) in the nucleotide pool. To determine the targets contributing to oxidative mutagenesis, DNA or nucleotides, we characterized spontaneous mutations and compared the distribution to those in mutMY and mutT strains, in which GC --> TA and AT --> CG were predominantly induced, respectively. The hotspots and sequence contexts where AT --> CG occurred frequently in sodAB fur strain were almost identical to those in mutT strain,whereas, those where GC --> TA occurred frequently in sodAB fur strain were quite different from those in mutMY strain. These observations suggested that AT --> CG is due to 8-oxodGTP, while GC --> TA is produced by some other lesion(s). The 2-oxodATP is also a major oxidative lesion in nucleotides, and strongly induces GC --> TA. The expression of cDNA for MTH1, which can hydrolyze 2-oxodATP as well as 8-oxodGTP, partially but significantly, suppressed the GC --> TA mutator phenotype of the sodAB fur strain, whereas, it did not for the mutMY strain. Additionally, the sequence contextby 2-oxodATP in E. coli was similar to that in sodAB fur strain. These results suggested that the targets contributing to oxidative mutagenesis in sodAB fur strain are nucleotides such as dGTP and dATP, rather than 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 nonlinear dynamic model of DNA with a sequence-dependent stacking term

No simple model exists that accurately describes the melting behavior and breathing dynamics of double-stranded DNA as a function of nucleotide sequence. This is especially true for homogenous and periodic DNA sequences, which exhibit large deviations in melting temperature from predictions made by additive thermodynamic contributions. Currently, no method exists for analysis of the DNA breathing dynamics of repeats and of highly G/C- or A/T-rich regions, even though such sequences are widespread in vertebrate genomes. Here, we extend the nonlinear Peyrard–Bishop–Dauxois (PBD) model of DNA to include a sequence-dependent stacking term, resulting in a model that can accurately describe the melting behavior of homogenous and periodic sequences. We collect melting data for several DNA oligos, and apply Monte Carlo simulations to establish force constants for the 10 dinucleotide steps (CG, CA, GC, AT, AG, AA, AC, TA, GG, TC). The experiments and numerical simulations confirm that the GG/CC dinucleotide stacking is remarkably unstable, compared with the stacking in GC/CG and CG/GC dinucleotide steps. The extended PBD model will facilitate thermodynamic and dynamic simulations of important genomic regions such as CpG islands and disease-related repeats.     

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.     

Effects of base substituents on the hydration of B- and Z-DNA: correlations to the B- to Z-DNA transition.

We present a study of how substituent groups of naturally occurring and modified nucleotide bases affect the degree of hydration of right-handed B-DNA and left-handed Z-DNA. A comparison of poly(dG-dC) and poly(dG-dm5C) titrations with the lipotropic salts of the Hofmeister series infers that the methyl stabilization of cytosines as Z-DNA is primarily a hydrophobic effect. The hydration free energies of various alternating pyrimidine-purine sequences in the two DNA conformations were calculated as solvent free energies from solvent accessible surfaces. Our analysis focused on the N2 amino group of purine bases that sits in the minor groove of the double helix. Removing this amino group from guanine to form inosine (I) destabilizes Z-DNA, while adding this group to adenines to form 2-aminoadenine (A') stabilizes Z-DNA. These predictions were tested by comparing the salt concentrations required to crystallize hexanucleotide sequences that incorporate d(CG), d(CI), d(TA) and d(TA') base pairs as Z-DNA. Combining the current results with our previous analysis of major groove substituents, we derived a thermodynamic cycle that relates the systematic addition, deletion, or substitution of each base substituent to the B- to Z-DNA transition free energy.     

Effects of an abasic site on triple helix formation characterized by affinity cleaving.

The stability of triple helical complexes of pyrimidine oligodeoxyribonucleotides containing one abasic 1,2-dideoxy-D-ribose (phi) residue was examined by affinity cleaving. Within a pyrimidine third strand, the triplets phi.AT, phi.GC, phi.TA and phi.CG are significantly less stable than the triplets, T.AT, C+GC and G.TA. The decrease in binding produced by an abasic residue is similar to that observed with imperfectly matched natural base triplets, with phi.AT and phi.GC being less stable than phi.TA and phi.CG triplets for the sequences studied.     

Molecular dynamics simulations of A. T-rich oligomers: sequence-specific binding of Na+ in the minor groove of B-DNA

Molecular dynamics simulations have been employed to probe the sequence-specific binding of sodium ions to the minor groove of B-DNA of three A. T-rich oligomers having identical compositions but different orders of the base pairs: C(AT)(4)G, CA(4)T(4)G, and CT(4)A(4)G. Recent experimental investigations, either in crystals or in solution, have shown that monovalent cations bind to DNA in a sequence-specific mode, preferentially in the narrow minor groove regions of uninterrupted sequences of four or more adenines (A-tracts), replacing a water molecule of the ordered hydration structure, the hydration spine. Following this evidence, it has been hypothesized that in A-tracts these events may be responsible for structural peculiarities such as a narrow minor groove and a curvature of the helix axis. The present simulations confirm a sequence specificity of the binding of sodium ions: Na(+) intrusions in the first layer of hydration of the minor groove, with long residence times, up to approximately 3 ns, are observed only in the minor groove of A-tracts but not in the alternating sequence. The effects of these intrusions on the structure of DNA depend on the ion coordination: when the ion replaces a water molecule of the spine, the minor groove becomes narrower. Ion intrusions may also disrupt the hydration spine modifying the oligomer structure to a large extent. However, in no case intrusions were observed to locally bend the axis toward the minor groove. The simulations also show that ions may reside for long time periods in the second layer of hydration, particularly in the wider regions of the groove, often leading to an opening of the groove.     

Molecular dynamics simulations of the 136 unique tetranucleotide sequences of DNA oligonucleotides. II: sequence context effects on the dynamical structures of the 10 unique dinucleotide steps

Molecular dynamics (MD) simulations including water and counterions on B-DNA oligomers containing all 136 unique tetranucleotide basepair steps are reported. The objective is to obtain the calculated dynamical structure for at least two copies of each case, use the results to examine issues with regard to convergence and dynamical stability of MD on DNA, and determine the significance of sequence context effects on all unique dinucleotide steps. This information is essential to understand sequence effects on DNA structure and has implications on diverse problems in the structural biology of DNA. Calculations were carried out on the 136 cases embedded in 39 DNA oligomers with repeating tetranucleotide sequences, capped on both ends by GC pairs and each having a total length of 15 nucleotide pairs. All simulations were carried out using a well-defined state-of-the-art MD protocol, the AMBER suite of programs, and the parm94 force field. In a previous article (Beveridge et al. 2004. Biophysical Journal. 87:3799-3813), the research design, details of the simulation protocol, and informatics issues were described. Preliminary results from 15 ns MD trajectories were presented for the d(CpG) step in all 10 unique sequence contexts. The results indicated the sequence context effects to be small for this step, but revealed that MD on DNA at this length of trajectory is subject to surprisingly persistent cooperative transitions of the sugar-phosphate backbone torsion angles alpha and gamma. In this article, we report detailed analysis of the entire trajectory database and occurrence of various conformational substates and its impact on studies of context effects. The analysis reveals a possible direct correspondence between the sequence-dependent dynamical tendencies of DNA structure and the tendency to undergo transitions that "trap" them in nonstandard conformational substates. The difference in mean of the observed basepair step helicoidal parameter distribution with different flanking sequence sometimes differs by as much as one standard deviation, indicating that the extent of sequence effects could be significant. The observations reveal that the impact of a flexible dinucleotide such as CpG could extend beyond the immediate basepair neighbors. The results in general provide new insight into MD on DNA and the sequence-dependent dynamical structural characteristics of DNA.     

Structural change in a B-DNA helix with hydrostatic pressure

Study of the effects of pressure on macromolecular structure improves our understanding of the forces governing structure, provides details on the relevance of cavities and packing in structure, increases our understanding of hydration and provides a basis to understand the biology of high-pressure organisms. A study of DNA, in particular, helps us to understand how pressure can affect gene activity. Here we present the first high-resolution experimental study of B-DNA structure at high pressure, using NMR data acquired at pressures up to 200 MPa (2 kbar). The structure of DNA compresses very little, but is distorted so as to widen the minor groove, and to compress hydrogen bonds, with AT pairs compressing more than GC pairs. The minor groove changes are suggested to lead to a compression of the hydration water in the minor groove.     

The influence of single base triplet changes on the stability of a pur.pur.pyr triple helix determined by affinity cleaving.

The influence of sixteen base triplet changes at a single position within a pur.pur.pyr triple helix was examined by affinity cleaving. For the 15 base pair target site studied here, G.GC, A.AT and T.AT triplets stabilize a triple helix to a greater extent than the other 13 natural triplets (pH = 7.4, 25 degrees C). Weaker interactions were detected for the C.AT, A.GC and T.CG triplets. The absence of specific, highly stabilizing interactions between third strand bases and the CG or TA base pairs demonstrates a current sequence limitation to formation of this structure. Models for the two dimensional base triplet interactions for all possible 16 natural triplets are presented.