Electronic structure of DNA by DV-X alpha cluster calculations: II. d(GG).d(CC), d(CG)2, d(GC)2 A and B conformations. (Part 2) Sugars and bases

The electronic structure of d(GG).d(CC), d(CG)2, d(GC)2 which are stacked base pairs in the DNA double helix, are elucidated for both A and B conformations in detail by DV-X alpha cluster calculations. These three DNA double helix fragments are contracted from the same bases, G and C, but the electronic structures of the fragments for both A and B conformations are different from each other characteristically. There are some delicate differences in the admixture of the orbital components and the overlap populations of intra- and inter- strand stacked bases among the stacking isomers. On the other hand, the electronic states of sugars differ in the 5'-3' direction, but are not almost dependent on stacked base pairs.     

A novel end-to-end binding of two netropsins to the DNA decamers d(CCCCCIIIII)2, d(CCCBr5CCIIIII)2and d(CBr5CCCCIIIII)2.

Netropsin is bound to the DNA decamer d(CCCCCIIIII)2, the C-4 bromo derivative d(CCCBr5CCIIIII)2and the C-2 bromo derivative d(CBr5CCCCIIIII)2in a novel 2:1 mode. Complexes of the native decamer and the C-4 bromo derivative are isomorphous, space group P1, unit cell dimensions a = 32.56 A (32.66), b = 32.59 A (32.77), c = 37.64 A (37.71), alpha = 86.30 degrees (86.01 degrees), beta = 84.50 degrees (84.37 degrees), gamma = 68.58 degrees (68.90 degrees) with two independent molecules (A and B) in the asymmetric unit (values in parentheses are for the derivative). The C-2 bromo derivative is hexagonal P61, unit cell dimensions a = b = 32.13 A, c = 143.92, gamma = 120 degrees with one molecule in the asymmetric unit. The structures were solved by the molecular replacement method. The novelty of the structures is that there are two netropsins bound end-to-end in the minor groove of each B-DNA decamer which has nearly a complete turn. The netropsins are held by hydrogen bonding interactions to the base atoms and by sandwiching van der Waal's interactions from the sugar-phosphate backbones of the double helix similar to every other drug.DNA complex. Each netropsin molecule spans approximately 5 bp. The netropsins refined with their guanidinium heads facing each other at the center, although an orientational disorder for the netropsins cannot be excluded. The amidinium ends stretch out toward the junctions and bind to the adjacent duplexes in the columns of stacked symmetry-related complexes. Both cationic ends of netropsin are bridged by water molecules in one of the independent molecules (molecule A) of the triclinic structures and also the hexagonal structure to form pseudo-continuous drug.decamer helices.     

Structure of d(T)n.d(A)n.d(T)n: the DNA triple helix has B-form geometry with C2'-endo sugar pucker.

The polynucleotide helix d(T)n.d(A)n.d(T)n is the only deoxypolynucleotide triple helix for which a structure has been published, and it is generally assumed as the structural basis for studies of DNA triplexes. The helix has been assigned to an A-form conformation with C3'-endo sugar pucker by Arnott and Selsing [1974; cf. Arnott et al. (1976)]. We show here by infrared spectroscopy in D2O solution that the helix is instead B-form and that the sugar pucker is in the C2'-endo region. Distamycin A, which binds only to B-form and not to A-form helices, binds to the triple helix without displacement of the third strand, as demonstrated by CD spectroscopy and gel electrophoresis. Molecular modeling shows that a stereochemically satisfactory structure can be build using C2'-endo sugars and a displacement of the Watson-Crick base-pair center from the helix axis of 2.5 A. Helical constraints of rise per residue (h = 3.26 A) and residues per turn (n = 12) were taken from fiber diffraction experiments of Arnott and Selsing (1974). The conformational torsion angles are in the standard B-form range, and there are no short contacts. In contrast, we were unable to construct a stereochemically allowed model with A-form geometry and C3'-endo sugars. Arnott et al. (1976) observed that their model had short contacts (e.g., 2.3 A between the phosphate-dependent oxygen on the A strand and O2 in the Hoogsteen-paired thymine strand) which are generally known to be outside the allowed range.(ABSTRACT TRUNCATED AT 250 WORDS)     

Crystal structure of d(GCGCGCG) with 5''-overhang G residues.

The crystal structure of the DNA heptamer d(GCGCGCG) has been solved at 1.65 A resolution by the molecular replacement method and refined to an R-value of 0.184 for 3598 reflections. The heptamer forms a Z-DNA d(CGCGCG)2 with 5'-overhang G residues instead of an A-DNA d(GCGCGC)2 with 3'-overhang G residues. The overhang G residues from parallel strands of two adjacent duplexes form a trans reverse Hoogsteen G x G basepair that stacks on the six Z-DNA basepairs to produce a pseudocontinuous helix. The reverse Hoogsteen G x G basepair is unusual in that the displacement of one G base relative to the other allows them to participate in a bifurcated (G1)N2 . . . N7(G8) and an enhanced (G8)C8-H . . . O6(G1) hydrogen bond, in addition to the two usual hydrogen bonds. The 5'-overhang G residues are anti and C2'-endo while the 3'-terminal G residues are syn and C2'-endo. The conformations of both G residues are different from the syn/C3'-endo for the guanosine in a standard Z-DNA. The two cobalt hexammine ions bind to the phosphate groups in both GpC and CpG steps in Z(I) and Z(II) conformations. The water structure motif is similar to the other Z-DNA structures.     

Structural studies of DNA fragments: the G.T wobble base pair in A, B and Z DNA; the G.A base pair in B-DNA

The crystal structures of five double helical DNA fragments containing non-Watson-Crick complementary base pairs are reviewed. They comprise four fragments containing G.T base pairs: two deoxyoctamers d(GGGGCTCC) and d(GGGGTCCC) which crystallise as A type helices; a deoxydodecamer d(CGCGAATTTGCG) which crystallises in the B-DNA conformation; and the deoxyhexamer d(TGCGCG), which crystallises as a Z-DNA helix. In all four duplexes the G and T bases form wobble base pairs, with bases in the major tautomer forms and hydrogen bonds linking N1 of G with O2 of T and O6 of G with N3 of T. The X-ray analyses establish that the G.T wobble base pair can be accommodated in the A, B or Z double helix with minimal distortion of the global conformation. There are, however, changes in base stacking in the neighbourhood of the mismatched bases. The fifth structure, d(CGCGAATTAGCG), contains the purine purine mismatch G.A where G is in the anti and A in the syn conformation. The results represent the first direct structure determinations of base pair mismatches in DNA fragments and are discussed in relation to the fidelity of replication and mismatch recognition.     

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.     

Interactions between an anthracycline antibiotic and DNA: molecular structure of daunomycin complexed to d(CpGpTpApCpG) at 1.2-A resolution

The crystal structure of a daunomycin-d(CGTACG) complex has been solved by X-ray diffraction analysis and refined to a final R factor of 0.175 at 1.2-A resolution. The crystals are in a tetragonal crystal system with space group P4(1)2(1)2 and cell dimensions of a = b = 27.86 A and c = 52.72 A. The self-complementary DNA forms a six base pair right-handed double helix with two daunomycin molecules intercalated in the d(CpG) sequences at either end of the helix. Daunomycin in the complex has a conformation different from that of daunomycin alone. The daunomycin aglycon chromophore is oriented at right angles to the long dimension of the DNA base pairs, and the cyclohexene ring A rests in the minor groove of the double helix. Substituents on this ring have hydrogen-bonding interactions to the base pairs above and below the intercalation site. O9 hydroxyl group of the daunomycin forms two hydrogen bonds with N3 and N2 of an adjacent guanine base. Two bridging water molecules between the drug and DNA stabilize the complex in the minor groove. In the major groove, a hydrated sodium ion is coordinated to N7 of the terminal guanine and the O4 and O5 of daunomycin with a distorted octahedral geometry. The amino sugar lies in the minor groove without bonding to the DNA. The DNA double helix is distorted with an asymmetrical rearrangement of the backbone conformation surrounding the intercalator drug. The sugar puckers are C1,C2'-endo, G2,C1'-endo, C11,C1'-endo, and G12,C3'-exo. Only the C1 residue has a normal anti-glycosyl torsion angle (chi = -154 degrees), while the other three residues are all in the high anti range (average chi = -86 degrees). This structure allows us to identify three principal functional components of anthracycline antibiotics: the intercalator (rings B-D), the anchoring functions associated with ring A, and the amino sugar. The structure-function relationships of daunomycin binding to DNA as well as other related anticancer drugs are discussed.     

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.     

Comparison of the B- and Z-form hairpin loop structures formed by d(CG)5T4(CG)5

The partially self-complementary synthetic DNA oligonucleotide d(CG)5T4(CG)5 has been studied by using 1H and 31P NMR and circular dichroism. Results show that, under low-salt conditions (120 mM NaCl buffer), an intramolecular hairpin loop exists in which the double-helical stem region is B-form and the thymidine loop residues have predominantly southern (C2'-endo) sugar conformations. The thymidine glycosidic torsion angles are intermediate between syn and anti or exist as an equilibrium mixture of residues in the two extremes. NOESY data indicate that the structure of the loop region is very similar to that found for d(CG)2T4(CG)2 [Hare, D. R., & Reid, B. R. (1986) Biochemistry 25, 5341-5350]. Under high-salt conditions (6 M NaClO4 buffer), the dominant form (approximately equal to 85%) is an intramolecular hairpin structure in which the stem region forms a Z-form double helix. As in the B-form, the loop thymidine residues are intermediate between the syn and anti conformations or exist as an equilibrium mixture of the two, but the thymidine sugar conformations differ in that they are biased toward northern (C3'-endo) conformations.     

Conformation of DNA modified at a d(GG) or a d(AG) site by the antitumor drug cis-diamminedichloroplatinum(II)

The purpose of this work was the comparison of the conformational changes induced in the double helix by the adducts formed at d(GG) and d(AG) sites in the reaction between the antitumor drug cis-diamminedichloroplatinum(II) (cis-DDP) and DNA. Two duplexes (20-mer) containing either a single d(A*G*) or a single d(G*G) adduct were studied by means of gel electrophoresis and artificial nuclease and chemical probes. It is shown that the d(G*G*) and the d(A*G*) adducts bend DNA similarly, but at the nucleotide level they distort differently the double helix. We suggest that the weaker interactions between platinated A residues and the other nucleotides, as compared to the interactions between platinated G residues and the other nucleotides, are largely responsible for the differences in the distortions induced in DNA by the d(A*G) and d(G*G*) adducts. This suggestion is supported by the study of the distortions induced in duplexes by the d(G*G*) adducts, one of the platinated G residues being paired with a T residue.     

X-ray and NMR studies of the DNA oligomer d(ATATAT): Hoogsteen base pairing in duplex DNA

We present and analyze the structure of the oligonucleotide d(ATATAT) found in two different forms by X-ray crystallography and in solution by NMR. We find that in both crystal lattices the oligonucleotide forms an antiparallel double helical duplex in which base pairing is of the Hoogsteen type. The double helix is apparently very similar to the standard B-form of DNA, with about 10 base pairs per turn. However, the adenines in the duplex are flipped over; as a result, the physicochemical features of both grooves of the helix are changed. In particular, the minor groove is narrow and hydrophobic. On the other hand, d(ATATAT) displays a propensity to adopt the B conformation in solution. These results confirm the polymorphism of AT-rich sequences in DNA. Furthermore, we show that extrahelical adenines and thymines can be minor groove binders in Hoogsteen DNA.     

100ps molecular dynamic simulation of d(TATCACC)2.

A heptanucleotide sequence d(TATCACC)2 from OR3 region of bacteriophage lambda is considered sufficient for the recognition of Cro protein. We present here results on molecular dynamic simulations on this sequence for 100 ps in 0.02 ps interval. The simulations are done using computer program GROMOS. The conformational results are averaged over each ps. The IUPAC torsional parameters for 100 conformations are illustrated using a wheal and a dial systems. Several other stereochemical parameters such as H-bonding lengths and angles, sugar puckers, helix twist and roll angles as also distances between opposite strand phosphorus are depicted graphically. We find that there is rupture of terminal H-bonds. The bases are tilted and shifted away from the helix axis giving rise to bifurcated H-bonds. H-bonds are seen even in between different base pairs. The role of these dynamic structural changes in the recognition of OR3 operator by Cro protein is discussed in the paper.     

Crystal structure of the 2:1 complex between d(GAAGCTTC) and the anticancer drug actinomycin D.

The crystal structures of the 2:1 complex of the self-complementary DNA octamer d(GAAGCTTC) with actinomycin D has been determined at 3.0 A resolution. This is the first example of a crystal structure of a DNA-drug complex in which the drug intercalates into the middle of a relatively long DNA segment. The results finally confirmed the DNA-actinomycin intercalation model proposed by Sobell & co-workers in 1971. The DNA molecule adopts a severely distorted and slightly kinked B-DNA-like structure with an actinomycin D molecule intercalated in the middle sequence, GC. The two cyclic depsipeptides, which differ from each other in overall conformation, lie in the minor groove. The complex is further stabilized by forming base-peptide and chromophore-backbone hydrogen bonds. The DNA helix appears to be unwound by rotating one of the base-pairs at the intercalation site. This single base-pair unwinding motion generates a unique asymmetrically wound helix at the binding site of the drug, i.e. the helix is loosened at one end of the intercalation site and tightened at the other end. The large unwinding of the DNA by the drug intercalation is absorbed mostly in a few residues adjacent to the intercalation site. The asymmetrical twist of the DNA helix, the overall conformation of the two cyclic depsipeptides and their interaction mode with DNA are correlated to each other and rationally explained.     

Molecular dynamics simulation of the DNA triplex d(TC)5.d(GA)5.d(C+T)5.

A molecular dynamics simulation of the DNA triple helix d(TC)5.d(GA)5.d(C+T)5 is described (C+ represents a protonated cytosine residue). The simulation has been performed using the program AMBER 3.1 and includes counterions and explicit solvent under periodic boundary conditions. Both the dynamic and time-averaged behaviour of the system has been analysed. Considerable deviations from the fibre-diffraction model for DNA triple helix structure are observed, including the repuckering of the purine strand sugars that has been identified in some nuclear magnetic resonance (n.m.r.) studies. The simulation suggests that this conformational change may be driven by the possibility of improved interactions between the phosphate groups of this strand and both the solvent and counterions. Several examples of a particular conformational transition are observed, involving correlated changes in the backbone angles alpha and gamma. These transitions provide a possible explanation for some unusual n.m.r. data that have been reported. The structure of the triple helix major groove also suggests an explanation for the observed stabilization of DNA triplexes by polyvalent cations, and their ability to interact with drugs that bind in the minor groove of DNA duplexes.     

Conformational properties and thermodynamics of the RNA duplex r(CGCAAAUUUGCG)2: comparison with the DNA analogue d(CGCAAATTTGCG)2.

The thermodynamic stability of nine dodecamers (four DNA and five RNA) of the same base composition has been compared by UV-melting. TheDeltaG of stabilisation were in the order: r(GACUGAUCAGUC)2>r(CGCAAATTTGCG)2 approximately r(CGCAUAUAUGCG)2>d(CGCAAATTTGCG)2 approximately r(CGCAAAUUUGCG)2>d(CGCATATATGCG)2 approximately d(GACTGATCAGTC)2>r(CGCUUUAAAGCG)2 approximately d(CGCTTTAAAGCG)2. Compared with the mixed sequences, both r(AAAUUU) and r(UUUAAA) are greatly destablising in RNA, whereas in DNA, d(TTTAAA) is destabilising but d(AAATTT) is stabilising, which has been attributed to the formation of a special B'structure involving large propeller twists of the A-T base pairs. The solution structure of the RNA dodecamer r(CGCAAAUUUGCG)2has been determined using NMR and restrained molecular dynamics calculations to assess the conformational reasons for its stability in comparison with d(CGCAAATTTGCG)2. The structures refined to a mean pairwise r.m.s.d. of 0.89+/-0.29 A. The nucleotide conformations are typical of the A family of structures. However, although the helix axis displacement is approximately 4.6 A into the major groove, the rise (3.0 A) and base inclination ( approximately 6 degrees ) are different from standard A form RNA. The extensive base-stacking found in the AAATTT tract of the DNA homologue that is largely responsible for the higher thermodynamic stability of the DNA duplex is reduced in the RNA structure, which may account for its low relative stability.     

N.m.r. and c.d. studies of the DNA fragments d(TATATATA) and d(TATATA) in solution

DNA fragments d(TATATATA) and d(TATATA) were studied in low-salt aqueous solutions and found to coexist in more than one conformer. 1H-n.m.r. demonstrates that single-stranded and double-stranded states are involved in the conformational coexistence. Circular dichroism spectroscopy indicates a global B-DNA stacking of bases in the fragments. 31P-n.m.r. resonances of the TpA and ApT phosphodiester bonds are substantially separated in the spectra of both d(TATATATA) and d(TATATA) duplexes to suggest an alternating architecture of their backbones. In fact, the oligonucleotide duplexes are much more alternating than the corresponding polynucleotide under the same solution conditions. The alternating character of the d(TATATATA) double helix is further enhanced in molar caesium fluoride solutions. The oligonucleotide isomerization into X-DNA is, however, accompanied by gel formation, which makes high resolution n.m.r. measurements impossible.     

Hydration of [d(CGC)r(aaa)d(TTTGCG)](2)

We have studied the hydration and dynamics of RNA C2'-OH in a DNA. RNA hybrid chimeric duplex [d(CGC)r(aaa)d(TTTGCG)](2). Long-lived water molecules with correlation time tau(c) larger than 0.3 ns were found close to the RNA adenine H2 and H1' protons in the hybrid segment. A possible long-lived water molecule was also detected close to the methyl group of 7T in the RNA-DNA junction but not to the other two thymine bases (8T and 9T). This result correlates with the structural studies that only DNA residue 7T in the RNA-DNA junction adopts an O4'-endo sugar conformation (intermediate between B-form and A-form), while the other DNA residues including 3C in the DNA-RNA junction, adopt C1'-exo or C2'-endo conformations (in the B-form domain). Based on the NOE cross-peak patterns, we have found that RNA C2'-OH tends to orient toward the O3' direction, forming a possible hydrogen bond with the 3'-phosphate group. The exchange rates for RNA C2'-OH were found to be around 5-20 s(-1), compared to 26.7(+/-13.8) s(-1) reported previously for the other DNA.RNA hybrid duplex. This slow exchange rate may be due to the narrow minor groove width of [d(CGC)r(aaa)d(TTTGCG)](2), which may trap the water molecules and restrict the dynamic motion of hydroxyl protons. The distinct hydration patterns of the RNA adenine H2 and H1' protons and the DNA 7T methyl group in the hybrid segment, as well as the orientation and dynamics of the RNA C2'-OH protons, may provide a molecular basis for further understanding the structure and recognition of DNA.RNA hybrid and chimeric duplexes.     

Poly(d2NH2A-dT): two-dimensional NMR shows a B to A conversion in high salt

Poly(d2NH2A-dT) forms a structure in high salt that is clearly distinct from the B form present in low salt. Two-dimensional nuclear Overhauser effect (2D NOE) NMR spectra establish that the conformation of the high-salt form is not Z. Correlations of observed cross peaks in the 2D NOE spectra and estimated interproton distances of the common DNA conformations are consistent only with an A form or closely related structure. This interpretation is also consistent with the negative circular dichroic band observed in the high-salt form of poly(d2NH2A-dT) and in A-form ribohomopolymer helices containing 2NH2A.     

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)