Conformational Analysis of a Hybrid DNA-RNA Double Helical Oligonucleotide in Aqueous Solution: d(CG)r(CG)d(CG) Studied by 1D- and 2D-1H NMR Spectroscopy

Abstract

The double helical structure of the self-complementary DNA-RNA-DNA hybrid d(CG)r(CG) d(CG) was studied in solution by 500 MHz ¹H-NMR spectroscopy. The non-exchangeable base protons and the (deoxy)ribose H1′, H2′ and H2″ protons were unambiguously assigned using 2D-J-correlated (COSY) and 2D-NOE (NOESY) spectroscopy techniques. A general strategy for the sequential assignment of ¹H-NMR spectra of (double) helical DNA and RNA fragments by means of 2D-NMR methods is presented.

Conformational analysis of the sugar rings of d(CG)r(CG)d(CG) at 300 K shows that the central ribonucleotide part of the helix adopts an A-type double helical conformation. The 5′- and 3′-terminal deoxyribose base pairs, however, take up the normal DNA-type conformation. The A-to-B transition in this molecule involves only one (deoxyribose) base pair. It is shown that this A-to-B conformational transition can only be accomodated by two specific sugar pucker combinations for the junction base pair, i.e. N·S (C3′-endo-C2′-endo, 60%, where the pucker given first is that assigned to the junction nucleotide residue of the strand running 5′ → 3′ from A-RNA to B-DNA) and S·S (C2′-endo-C2′-endo, 40%).     

The B-A transition in superhelical DNA.

Relaxation of a DNA superhelical stress due to the B to A transition induced by trifluoroethanol has been studied by assessing the change of DNA orientation in a flow gradient. Using DNAs of different superhelical densities, a decrease in the winding angle during the B----A shift of DNA was found to be 1.5 degrees per base pair in solution. Accepting the winding angle for B-DNA in solution to be 34.1 degrees, that for A-DNA must have a value of 32.6 degrees which agrees with the X-ray data for A-DNA in the condensed state. The date obtained within the B-A transition interval make it possible to conclude that there is an increase in winding at each B/A junction, which is about 5 degrees per one junction.     

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.     

Conformational microheterogeneity in a DNA double helix: structure of restriction endonuclease Bam H1 recognition site

Structural studies using 500 MHz 1H NMR spectroscopy on Bam H1 recognition site d(GGATCC)2 in solution at 19 degrees is reported. The resonances from the sugar ring and base protons have been assigned from the 2D-COSY and NOESY spectra. Analyses of the NOESY cross-peaks between the base protons H8/H6 and sugar protons H2'/H2", H3' reveal that the nucleotide units G2, A3 and C6 adopt (C3'-endo, chi = 200 degrees-220 degrees) conformation while G1, T4 and C5 exhibit (C2'-endo, chi = 240 degrees-260 degrees) conformation. NMR data clearly suggest that the two strands of d(GGATCC)2 are conformationally equivalent and there is a structural two-fold between the two A-T pairs. The above information and the NOESY data are used to generate a structural model of d(GGATCC)2. The important features are: (i) G1-G2 stack, the site of cleavage, shows an alternation in sugar pucker i.e. C2'-endo, C3'-endo as in a B-A junction, (ii) G2-A3 stack adopts a mini A-DNA, both the sugars being C3'-endo, (iii) A3-T4 stack, the site of two-fold, displays an A-B junction with alternation in sugar pucker as C3'-endo, C2'-endo, (iv) T4-C5 stack adopts a mini B-DNA both the sugars being C2'-endo and (v) C5-C6 stack exhibits a B-A junction with C2'-endo, C3'-endo sugar puckers. Thus, our studies demonstrate that conformational microheterogeneity with a structural two fold, is present in the Bam H1 recognition site.     

Structure of the 2-aminopurine-cytosine base pair formed in the polymerase active site of the RB69 Y567A-DNA polymerase

The adenine base analogue 2-aminopurine (2AP) is a potent base substitution mutagen in prokaryotes because of its enhanceed ability to form a mutagenic base pair with an incoming dCTP. Despite more than 50 years of research, the structure of the 2AP-C base pair remains unclear. We report the structure of the 2AP-dCTP base pair formed within the polymerase active site of the RB69 Y567A-DNA polymerase. A modified wobble 2AP-C base pair was detected with one H-bond between N1 of 2AP and a proton from the C4 amino group of cytosine and an apparent bifurcated H-bond between a proton on the 2-amino group of 2-aminopurine and the ring N3 and O2 atoms of cytosine. Interestingly, a primer-terminal region rich in AT base pairs, compared to GC base pairs, facilitated dCTP binding opposite template 2AP. We propose that the increased flexibility of the nucleotide binding pocket formed in the Y567A-DNA polymerase and increased "breathing" at the primer-terminal junction of A+T-rich DNA facilitate dCTP binding opposite template 2AP. Thus, interactions between DNA polymerase residues with a dynamic primer-terminal junction play a role in determining base selectivity within the polymerase active site of RB69 DNA polymerase.     

Solution structure of the dodecamer d-(CATGGGCC-CATG)2 is B-DNA. Experimental and molecular dynamics study

The DNA duplex d-(CATGGGCCCATG)2 has been studied in solution by FTIR, NMR and CD. The experimental approaches have been complemented by series of large-scale unrestrained molecular dynamics simulation with explicit inclusion of solvent and counterions. Typical proton-proton distances extracted from the NMR spectra and the CD spectra are completely in agreement with slightly modified B-DNA. By molecular dynamics simulation, starting from A-type sugar pucker, a spontaneous repuckering to B-type sugar pucker was observed. Both experimental and theoretical approaches suggest for the dodecamer d-(CATGGGCCCATG)2 under solution conditions puckering of all 2'-deoxyribose residues in the south conformation (mostly C2'-endo) and can exclude significant population of sugars in the north conformation (C3'-endo). NMR, FTIR and CD data are in agreement with a B-form of the dodecamer in solution. Furthermore, the duplex shows a cooperative B-A transition in solution induced by addition of trifluorethanol. This contrasts a recently published crystal structure of the same oligonucleotide found as an intermediate between B- and A-DNA where 23 out of 24 sugar residues were reported to adopt the north (N-type) conformation (C3'-endo) like in A-DNA (Ng, H. L., Kopka, M. L. and Dickerson, R. E., Proc. Natl. Acad. Sci. U S A 97, 2035-2039 (2000)). The simulated structures resemble standard B-DNA. They nevertheless show a moderate shift towards A-type stacking similar to that seen in the crystal, despite the striking difference in sugar puckers between the MD and X-ray structures. This is in agreement with preceding MD reports noticing special stacking features of G-tracts exhibiting a tendency towards the A-type stacking supported by the CD spectra also reflecting the G-tract stacking. MD simulations reveal several noticeable local conformational variations, such as redistribution of helical twist and base pair roll between the central GpC steps and the adjacent G-tract segments, as well as a substantial helical twist variability in the CpA(TpG) steps combined with a large positive base pair roll. These local variations are rather different from those seen in the crystal.     

NMR studies for identification of dI:dG mismatch base-pairing structure in DNA.

One- and two-dimensional NMR experiments have been undertaken to investigate deoxyinosine:deoxyguanosine (dI:dG) base pairing in a self-complementary dodecadeoxyribonucleotide, d(C1-G2-C3-I4-A5-A6-T7-T8-G9-G10-G11-G12) (designated IG-12), duplex. The NMR data indicate formation of a dI(syn):dG(anti) base pair in a B-DNA helix. This unusual base pairing results in altered NOE patterns between the base protons (H8 and H2) of the I4 residue and the sugar protons of its own and the 5'-flanking C3 residues. The dI(syn):dG(anti) base pair is accommodated in the B-DNA duplex with only a subtle distortion of the local conformation. Identification of the dI:dG base pairing in this study confirms that a hypoxanthine base can form hydrogen-bonded base pairs with all of the four normal bases, C, A, T, and G, in DNA.     

Crystal structure of the highly distorted chimeric decamer r(C)d(CGGCGCCG)r(G).spermine complex--spermine binding to phosphate only and minor groove tertiary base-pairing.

The crystal structure of the self-complementary chimeric decamer duplex r(C)d(CGGCGCCG)r(G), with RNA base pairs at both termini, has been solved at 1.9 A resolution by the molecular replacement method and refined to an R value of 0.145 for 2,314 reflections. The C3'-endo sugar puckers of the terminal riboses apparently drive the entire chimeric duplex into an A-DNA conformation, in contrast to the B-DNA conformation adopted by the all-deoxy decamer of the same sequence. Five symmetry related duplexes encapsulate a spermine molecule which interacts with ten phosphate groups, both directly and through water molecules to form multiple ionic and hydrogen bonding interactions. The spermine interaction severely bends the duplexes by 31 degrees into the major groove at the fourth base pair G(4).C(17), jolts it and slides the 'base plate' into the minor groove. This base pair, together with the adjacent base pair in the top half and the corresponding pseudo two-fold related base pairs in the bottom half, form four minor groove base-paired multiples with the terminal base pairs of two neighboring duplexes.     

Base Dependence of B-DNA Sugar Conformation in Solution and in the Solid State

Abstract

Analysis of 1H-NOESY solution data for eight short DNA duplexes has revealed pronounced differences between the sugar conformations of purine and pyrimidine nucleotides. It was found that the H1′-H4′ interproton distance is less than ca. 3.0 A. in pyrimidine sugars, while in purine sugars it is more than ca. 3.OA. This difference has been analyzed by comparison with the sugar conformations of highly resolved B-DNA crystal structures and model sugar conformations. The conclusion can be drawn that the deoxyribose conformation is of the general C2′-endo type but pyrimidine sugars are characterized by smaller phase angles of pseudorotation P (90°<P<150°), while purine sugars have larger P values that are greater than ca. 140° (140°<P<180°). There is no such clear base dependence of sugar conformation in highly resolved B-DNA crystal structures; however the similar trend can be seen as in the solution studies. Based on B-type DNA crystal structures, I-coupling constants have been calculated, and the applicability of experimental coupling measurements to the determination of sugar conformation is discussed.     

Characterization of DNA structures by Raman spectroscopy: high-salt and low-salt forms of double helical poly(dG-dC) in H2O and D2O solutions and application to B, Z and A-DNA.

Raman spectra of poly(dG-dC) . poly(dG-dC) in D2O solutions of high (4.0M NaCl) and low-salt (0.1M NaCl) exhibit differences due to different nucleotide conformations and secondary structures of Z and B-DNA. Characteristic carbonyl modes in the 1600-1700 cm-1 region also reflect differences in base pair hydrogen bonding of the respective GC complexes. Comparison with A-DNA confirms the uniqueness of C = O stretching frequencies in each of the three DNA secondary structures. Most useful for qualitative identification of B, Z and A-DNA structures are the intense Raman lines of the phosphodiester backbone in the 750-850 cm-1 region. A conformation-sensitive guanine mode, which yields Raman lines near 682, 668, or 625 cm-1 in B (C2'-endo, anti), A (C3'-endo, anti) or Z (C3'-endo, syn) structures, respectively, is the most useful for quantitative analysis. In D2O, the guanine line of Z-DNA is shifted to 615 cm-1, permitting its detection even in the presence of proteins.     

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.     

Conformational sub-states in B-DNA.

Theoretical studies of the sequence-dependent conformation of B-DNA have been carried out using Jumna, a helicoidal co-ordinate minimization algorithm. The results obtained for a series of six oligomers with repetitive sequences show that, with the exception of the homopolymers (dA)n.(dT)n and (dG)n.(dC)n, all sequences can adopt a variety of conformations characterized by considerable changes in helicoidal parameters and also in sugar puckers which adopt C(2')-endo (falling into 2 classes) or, in the case of pyrimidine nucleotides, O(1')-endo forms. These studies lead to an improved understanding of the role of base sequence on DNA conformation and point to a number of interesting correlations between the various structural parameters describing the double helix.     

Invariant and variable base stacking geometries in B-DNA and A-DNA

We calculated the interatomic distances between all couples of non-hydrogen atoms belonging to the neighboring Watson-Crick base pairs in the available crystal structures of DNA. Their standard deviations revealed remarkably large differences in the variability of the base stacking geometries of the particular steps. In line with experimental studies in solution, (CpA)-(TpG) and (TpA).(TpA) were identified as the most variable or flexible steps in the crystal structures of B-DNA. On the other hand, base stacking geometries of the (ApT).(ApT) steps were the most invariant, which was very surprising because all three steps composed only of C and G were much more flexible. This finding suggests that conformational stability of DNA and the rigidity have different origins. Furthermore, the nucleotide sequence dependence of the flexibility was almost reversed in A-DNA because the most flexible steps in B-DNA were the least flexible in A-DNA. The most invariant steps of B-DNA were variable in A-DNA. The (ApT).(ApT) step was a notable exception to this rule because it belonged to the most rigid steps in both B-DNA and A-DNA. The present results are fully consistent with the properties that poly(dA-dT).poly(dA-dT), poly(dA).poly(dT), poly(dAdC).poly(dG-dT) and poly(dA-dG).poly(dC-dT) exhibit in solution.     

A reduced set of coordinates for modeling DNA structures: (I). A B-to-A transition pathway driven by pseudorotational angle.

The A-DNA and the B-DNA are two well characterized polymorphous forms of DNA duplex. By using Metropolis Monte Carlo Simulations in a reduced coordinate space, we have shown that the B in equilibrium with A transitions can be induced by forcing pseudorotational angle (W) to change between C3'-endo and C2'-endo puckerings. The energy barrier for the transition pathway is less than 10 Kcal.mol-1. Base-pair parameters x-displacement (Dx) and roll (rho), which have the largest differences between the two forms of structures, cannot drive the transition. Our results support the view that the bistable states of the DNA duplex are due to the bistable structures of the sugar ring.     

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.     

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)     

Analysis of the possible helical structures of nucleic acids and polynucleotides. Application of (n-h) plots.

The two helical parameters n and h where n is the number of nucleotide residues per turn and h is the height per nucleotide residue have been evaluated for single stranded helical polynucleotide chains comprising C(3') -endo and C(2') endo class of nucleotides. The helical parameters are found to be especially sensitive to the C(4')-C(3') (sugar pucker) and the C(4')-C(5') torsions. The (n-h) plots display only one important helix forming domain for each class of nucleotides characterized by the sugar pucker and the C(4')-C(5') torsion. A correlation between the (n-h) plots and the known RNA (A,A') and DNA (A,B,C) helical forms has been established. It is found that all forms of helices except the C-DNA possess a favorable combination of P-O torsions. The analysis of the (n-h) plots suggests that C-DNA can have a conformation very similar to B-DNA. Although the (n-h) plots predict the stereochemical possibility of both right-handed and left-handed helices, nucleic acids apparently prefer right-handed conformation because of the energetics associated with the sugar-phosphate backbone and the base.     

The solution structure of [d(CGC)r(aaa)d(TTTGCG)](2): hybrid junctions flanked by DNA duplexes

The solution structure and hydration of the chimeric duplex [d(CGC)r(aaa)d(TTTGCG)]₂, in which the central hybrid segment is flanked by DNA duplexes at both ends, was determined using two-dimensional NMR, simulated annealing and restrained molecular dynamics. The solution structure of this chimeric duplex differs from the previously determined X-ray structure of the analogous B-DNA duplex [d(CGCAAATTTGCG)]₂ as well as NMR structure of the analogous A-RNA duplex [r(cgcaaauuugcg)]₂. Long-lived water molecules with correlation time τ_c longer 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 with the other two thymines (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, while the other DNA residues including 3C in the DNA–RNA junction, adopt C1′-exo or C2′-endo conformations. The exchange rates for RNA C2′-OH were found to be ~5–20 s^–1. This slow exchange rate may be due to the narrow minor groove width of [d(CGC)r(aaa)d(TTTGCG)]₂, which may trap the water molecules and restrict the dynamic motion of hydroxyl protons. The minor groove width of [d(CGC)r(aaa)d(TTTGCG)]₂ is wider than its B-DNA analog but narrower than that of the A-RNA analog. It was further confirmed by its titration with the minor groove binding drug distamycin. A possible 2:1 binding mode was found by the titration experiments, suggesting that this chimeric duplex contains a wider minor groove than its B-DNA analog but still narrow enough to hold two distamycin molecules. These distinct structural features and hydration patterns of this chimeric duplex provide a molecular basis for further understanding the structure and recognition of DNA·RNA hybrid and chimeric duplexes.