Mutual backbone phosphate group interactions promote DNA double helix bending at high salt concentrations in solution

Results of free energy calculations connected with the backbone phosphate group interactions upon local bending and helical twist modifications of A-, B- or Z-DNA at high salt concentrations have been reported recently (Jursa and Kypr 1990). Here we calculate energies necessary for DNA bending, using three models based on experimentally determined persistence length values. A comparison of energies following from the two quite different approaches suggests that high salt concentrations induce A- and mainly B-DNA bending into the double helix minor groove at least up to 10 degrees.     

The extended and eccentric E-DNA structure induced by cytosine methylation or bromination

Cytosine methylation or bromination of the DNA sequence d(GGCGCC)2 is shown here to induce a novel extended and eccentric double helix, which we call E-DNA. Like B-DNA, E-DNA has a long helical rise and bases perpendicular to the helix axis. However, the 3'-endo sugar conformation gives the characteristic deep major groove and shallow minor groove of A-DNA. Also, if allowed to crystallize for a period of time longer than that yielding E-DNA, the methylated sequence forms standard A-DNA, suggesting that E-DNA is a kinetically trapped intermediate in the transition to A-DNA. Thus, the structures presented here chart a crystallographic pathway from B-DNA to A-DNA through the E-DNA intermediate in a single sequence. The E-DNA surface is highly accessible to solvent, with waters in the major groove sitting on exposed faces of the stacked nucleotides. We suggest that the geometry of the waters and the stacked base pairs would promote the spontaneous deamination of 5-methylcytosine in the transition mutation of dm5C-dG to dT-dA base pairs.     

Displacements of backbone vibrational modes of A-DNA and B-DNA.

We display the displacement vectors or eigenvectors of calculations of the A- and B-DNA backbones. These calculations are based on a refinement scheme that simultaneously fit several backbone modes of A-DNA, B-DNA, and A-RNA. We discuss the role of symmetry operations in mode calculations and the relevance of these displacement vectors to the interpretation of linear dichroism measurements performed on the A- and B-DNA helix.     

Molecular structure of an A-DNA decamer d(ACCGGCCGGT)

The molecular structure of the DNA decamer d(ACCGGCCGGT) has been solved and refined by single-crystal X-ray-diffraction analysis at 0.20 nm to a final R-factor of 18.0%. The decamer crystallizes as an A-DNA double helical fragment with unit-cell dimensions of a = b = 3.923 nm and c = 7.80 nm in the space group P6(1)22. The overall conformation of this A-DNA decamer is very similar to that of the fiber model for A-DNA which has a large average base-pair tilt and hence a wide and shallow minor groove. This structure is in contrast to that of several A-DNA octamers in which the molecules all have low base-pair-tilt angles (8-12 degrees) resulting in an appearance intermediate between B-DNA and A-DNA. The average helical parameters of this decamer are typical of A-DNA with 10.9 base pairs/turn of helix, an average helical twist angle of 33.1 degrees, and a base-pair-tilt angle of 18.2 degrees. However, the CpG step in this molecule has a low local-twist angle of 24.5 degrees, similar to that seen in other A-DNA oligomers, and therefore appears to be an intrinsic stacking pattern for this step. The molecules pack in the crystal using a recurring binding motif, namely, the terminal base pair of one helix abuts the surface of the shallow minor groove of another helix. In addition, the GC base pairs have large propeller-twist angles, unlike those found most other A-DNA structures.     

Anisotropic flexibility of DNA depends on the base sequence. Conformation calculations of double-stranded tetranucleotides AAAA:TTTT, (AATT)2, (TTAA)2, GGGG:CCCC, (GGCC)2, (CCGG)2

The bending flexibility of six tetramers was studied in an assumption that they were extended in the both directions by regular double helices. The bends of B-DNA in different directions were considered. The stiffness of the B-DNA double helix when bent into the both grooves proved to be less pronounced than in the perpendicular direction by the order of magnitude. Such an anisotropy is a feature of the sugar-phosphate backbone structure. The calculated fluctuations of the DNA bending along the dyad axis, 5-7 degrees, are in agreement with the experimental value of DNA persistence length. Anisotropy of the double helix is sequence-dependent: most easily bent into the minor groove are the tetramers with purine-pyrimidine dimer (RY) in the middle. In contrast, YR dinucleotides prefer bending into the major groove, moreover, they have an equilibrium bend of 6-12 degrees into this groove. The above inequality is caused by the stacking interaction of the bases. The bend in the central dimers is distributed to some extent between the adjacent links, though the main fraction of the bend remains within the central link. Variation of the sugar-phosphate geometry in the bent helix is unessential, so that DNA remains within the limits of the B-family of forms: namely, when the helical axis is bent by 20 degrees the backbone dihedral angles vary by no more than 15 degrees. The obtained results are in accord with the X-ray structure of B-DNA dodecamer; they further substantiate our earlier model of DNA wrapping in the nucleosome by means of "mini-kinks" separated by a half-pitch of the double helix, i.e. by 5-6 b. p. Sequence-dependent anisotropy of DNA presumably dictates the three-dimensional structure of DNA in solution as well. We have found that nonrandom allocation of YR dimers leads to the systematic bends in the equilibrium structure of certain DNA fragments. To the four "Calladine rules" two more can be added: the minor-groove steric clash of purines in the YR sequences are avoided by: (1) bending of the helix into the major groove; (2) increasing the distance between the base pairs (stretching the double helix).     

Comparison of the solution and crystal conformations of (G + C)-rich fragments of DNA.

DNA fragments crystallize in an unpredictable manner, and relationships between their crystal and solution conformations still are not known. We have studied, using circular dichroism spectroscopy, solution conformations of (G + C)-rich DNA fragments, the crystal structures of which were solved in the laboratory of one of the present authors. In aqueous trifluorethanol (TFE) solutions, all of the examined oligonucleotides adopted the same type of double helix as in the crystal. Specifically, the dodecamer d(CCCCCGCGGGGG) crystalized as A-DNA and isomerized into A-DNA at high TFE concentrations. On the other hand, the hexamer d(CCGCGG) crystallized in Z-form containing tilted base pairs, and high TFE concentrations cooperatively transformed it into the same Z-form as adopted by the RNA hexamer r(CGCGCG), although d(CCGCGG) could isomerize into Z-DNA in the NaCl + NiCl2) aqueous solution. The fragments crystallizing as B-DNA remained B-DNA, regardless of the solution conditions, unless they denatured or aggregated. Effects on the oligonucleotide conformation of 2-methyl-2,4-pentanediol and other crystallization agents were also studied. 2-Methyl-2,4-pentanediol induced the same conformational transitions as TFE but, in addition, caused an oligonucleotide condensation that was also promoted by the other crystallization agents. The present results indicate that the crystal double helices of DNA are stable in aqueous TFE rather than aqueous solution.     

Transitions of poly(dI-dC), poly(dI-methyl5dC) and poly(dI-bromo5dC) among and within the B-, Z-, A- and X-DNA families of conformations.

It is shown, using circular dichroism spectroscopy, that poly(dI-dC) is capable to isomerize into both Z-DNA and A-DNA in concentrated NaCl + NiCl2 and trifluoroethanol solutions, respectively. This polynucleotide also undergoes a cooperative, two-state transition in ethanol into a structure which most probably is a canonical B-DNA. This implies that the conformation of poly(dI-dC) is unusual in low-salt aqueous solution. The canonical B-DNA is also adopted by poly(dI-methyl5dC) in trifluoroethanol while this polynucleotide adopts Z-DNA not only in NaCl + NiCl2 but also in the presence of MgCl2. Poly(dI-methyl5dC) partially adopts X-DNA in concentrated CsF and mainly ethanolic solutions. Poly(dI-bromo5dC) isomerizes into Z-DNA not only in concentrated NaCl even in the absence of NiCl2 but also in concentrated MgCl2. This polynucleotide transforms between two distinct variants of Z-DNA in ethanol or trifluoroethanol solutions.     

Neomycin, spermine and hexaamminecobalt (III) share common structural motifs in converting B- to A-DNA.

The (dG)n.(dC)n-containing 34mer DNA duplex [d(A2G15C15T2)]2 can be effectively converted from the B-DNA to the A-DNA conformation by neomycin, spermine and Co(NH3)6(3+). Conversion is demonstrated by a characteristic red shift in the circular dichroism spectra and dramatic NMR spectral changes in chemical shifts. Additional support comes from the substantially stronger CH6/GH8-H3'NOE intensities of the ligand-DNA complexes than those from the native DNA duplex. Such changes are consistent with a deoxyribose pucker transition from the predominate C2'-endo (S-type) to the C3'-endo (N-type). The changes for all three ligand-DNA complexes are identical, suggesting that those three complex cations share common structural motifs for the B- to A-DNA conversion. The A-DNA structure of the 4:1 complex of Co(NH3)6(3+)/d(ACCCGCGGGT) has been analyzed by NOE-restrained refinement. The structural basis of the transition may be related to the closeness of the two negatively charged sugar-phosphate backbones along the major groove in A-DNA, which can be effectively neutralized by the multivalent positively charged amine functions of these ligands. In addition, ligands like spermine or Co(NH3)6(3+) can adhere to guanine bases in the deep major groove of the double helix, as is evident from the significant direct NOE cross-peaks from the protons of Co(NH3)6(3+) to GH8, GH1 (imino) and CH4 (amino) protons. Our results point to future directions in preparing more potent derivatives of Co(NH3)6(3+) for RNA binding or the induction of A-DNA.     

Helix bending as a factor in protein/DNA recognition

Normal vectors perpendicular to individual base pairs are a powerful tool for studying the bending behavior of B-DNA, both in the form of normal vector plots and in matrices that list angles between vectors for all possible base pair combinations. A new analysis program, FREEHELIX, has been written for this purpose, and applied to 86 examples of sequence-specific protein/DNA complexes whose coordinates are on deposit in the Nucleic Acid Data Base. Bends in this sample of 86 structures almost invariably follow from roll angles between adjacent base pairs; tilt makes no net contribution. Roll in a direction compressing the broad major groove is much more common than that which compresses the minor groove. Three distinct types of B-DNA bending are observed, each with a different molecular origin: (1) Localized kinking is produced by large roll at single steps or at two steps separated by one turn of helix. (2) Smooth, planar curvature is produced by positive and negative roll angles spaced a half-turn apart, with random side-to-side zigzag roll at intermediate points, rather than a tilt contribution that might have been expected theoretically. (3) Three-dimensional writhe results from significant roll angles at a continuous series of steps. Writhe need not change the overall direction of helix axis, if it is continued indefinitely or for an integral number of helical turns. A-DNA itself can be formally considered as possessing uniform, continuous writhe that yields no net helix bending. Smooth curvature is the most intricate deformation of the three, and is least common. Writhe is the simplest deformation and is most common; indeed, a low level of continuous writhe is the normal condition of an otherwise unbent B-DNA helix of general sequence. With one exception, every example of major kinking in this sample of 86 structures involves a pyrimidine–purine step: C–A/T–G, T–A, or C–G. Purine–purine steps, especially A–A, show the least tendency toward roll deformations. © 1998 John Wiley & Sons, Inc. Biopoly 44: 361–403, 1997     

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.     

DNA structure from A to Z--biological implications of structural diversity of DNA

Deoxyribonucleic acid (DNA) is a biopolymer of nucleotides, usually adopting a double-stranded helical form in cells, with complementary base pairing holding the two strands together. The most stable is B-DNA conformation, although numerous other double helical structures can occur under specific conditions (A-DNA, Z-DNA, P-DNA). The existence of multiple-stranded (triplex, tetraplex) forms in vivo and their biological function in cells are subject of intensive studies.     

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.     

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.     

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)     

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.     

Sequence-dependent anisotropic flexibility of B-DNA. A conformational study

Bending flexibility of the six tetrameric duplexes was investigated d(AAAA):d(TTTT), d(AATT)2, d(TTAA)2, d(GGGG):d(CCCC), d(GGCC)2 and d(CCGG)2,. The tetramers were extended in the both directions by regular double helices. The stiffness of the B-DNA double helix when bent into the both grooves proved to be less than that in the perpendicular direction by an order of magnitude. Such an anisotropy is a property of the sugar-phosphate backbone structure. The calculated fluctuations of the DNA bending along the dyad axis, 5-7 degree, are in agreement with experimental value of the DNA persistence length. Anisotropy of the double helix is sequence-dependent: most easily bent into the minor groove are the tetramers with purine-pyrimidine dimer (RY) in the middle. In contrast, YR dinucleotides prefer bending into the major groove. Moreover, they have an equilibrium bend of 6-12 degree into this groove. The above inequality is caused by stacking interaction of the bases. The bend in the central dimer is distributed to some extent between the adjacent links, though the main fraction of the bend remains within the central link. Variation of the sugar-phosphate geometry in the bent helix is inessential, so that DNA remains within the B-family of forms: namely, when the helical axis is bent by 20 degree. the backbone dihedral angles vary by no more than 15 degree. The obtained results are in accord with x-ray structure of the B-DNA dodecamer; they further substantiate our early model of DNA wrapping in the nucleosome by means of "mini-kinks" separated by a half-pitch of the double helix, i.e. by 5-6 b.p. Sequence-dependent anisotropy of DNA presumably dictates the three-dimensional structure of DNA in solution as well. We have found that nonrandom allocation of YR dimers leads to the systematic bends in equilibrium structure of certain DNA fragments.     

The electrostatic potential and steric accessibility of reactive sites within Z-DNA.

A theoretical study of indices potentially useful for investigation of the reactivity of the recently discovered Z-DNA double helix is presented. The electrostatic potential minima and the steric accessibility of reactive sites are calculated. The effect of screening the phosphate groups by metal cations is investigated. The results are compared with those for the B-DNA double helix.     

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.     

Spatial translational motions of base pairs in DNA molecules: Application of the extended matrix generator method

We have used the elementary generator matrices outlined in the preceding paper to examine the conformational plasticity of the nucleic acid double helix. Here we investigate kinked DNA structures made up of alternating B- and A-type helices and intrinsically curved duplexes perturbed by the intercalation of ligands. We model the B-to-A transition by the lateral translation of adjacent base pairs, and the intercalation of ligands by the vertical displacement of neighboring residues. We report a complete set of average configuration-dependent parameters, ranging from scalars (i.e., persistence lengths) to first- and second-order tensor parameters (i.e., average second moments of inertia), as well as approximations of the associated spatial distributions of the DNA and their angular correlations. The average structures of short chains (of lengths less than 100 base pairs) with local kinks or intrinsically curved sequences are essentially rigid rods. At the smallest chain lengths (10 base pairs), the kinked and curved chains exhibit similar average properties, although they are structurally perturbed compared to the standard B-DNA duplex. In contrast, at lengths of 200 base pairs, the curved and kinked chains are more compact on average and are located in a different space from the standard B- or A-DNA helix. While A-DNA is shorter and thicker than B-DNA in x-ray models, the long flexible A-DNA helix is thinner and more extended on average than its B-DNA counterpart because of more limited fluctuations in local structure. Curved polymers of 50 base pairs or longer also show significantly greater asymmetry than other DNAs (in terms of the distribution of base pairs with respect to the center of gravity of the chain). The intercalation of drugs in the curved DNA straightens and extends the smoothly deformed template. The dimensions of the average ellipsoidal boundaries defining the configurations of the intercalated polymers are roughly double those of the intrinsically curved chain. The altered proportions and orientations of these density functions reflect the changing shape and flexibility of the double helix. The calculations shed new light on the possible structural role of short A-DNA fragments in long B-type duplexes and also offer a model for understanding how GC-specific intercalative ligands can straighten naturally curved DNA. The mechanism is not immediately obvious from current models of DNA curvature, which attribute the bending of the chain to a perturbed structure in repeating tracts of A · T base pairs. © 1994 John Wiley & Sons, Inc.     

Vacuum ultraviolet circular dichroism as an indicator of helical handedness in nucleic acids.

Calculated circular dichroism spectra are presented for double-stranded polynucleotides of regular sequences in A-RNA, A-DNA, B-DNA, and Z-DNA conformations. Quantum mechanical matrix method calculations were carried out in the near and vacuum ultraviolet regions. In the near UV, the calculated spectra agreed qualitatively with the measured spectra. However in the far and vacuum UV, the calculated CD compared nearly quantitatively with the experimental spectra. The calculations show that the sign of the CD in the vacuum UV, in contrast to that in the near UV, can be correlated with the handedness of the helix.