Solution structure of a nitrous acid induced DNA interstrand cross-link

Nitrous acid is a mutagenic agent. It can induce interstrand cross-links in duplex DNA, preferentially at d(CpG) steps: two guanines on opposite strands are linked via a single shared exocyclic imino group. Recent synthetic advances have led to the production of large quantities of such structurally homogenous cross-linked duplex DNA. Here we present the high resolution solution structure of the cross-linked dodecamer [d(GCATCCGGATGC)]2 (the cross-linked guanines are underlined), determined by 2D NMR spectroscopy, distance geometry, restrained molecular dynamics and iterative NOE refinement. The cross-linked guanines form a nearly planar covalently linked 'G:G base pair' with only minor propeller twisting, while the cytidine bases of their normal base pairing partners have been flipped out of the helix and adopt well defined extrahelical positions in the minor groove. On the 5'-side of the cross-link, the minor groove is widened to accommodate these extrahelical bases, and the major groove becomes quite narrow at the cross-link. The cross-linked 'G:G base pair' is well stacked on the spatially adjacent C:G base pairs, particularly on the 3'-side guanines. In addition to providing the first structure of a nitrous acid cross-link in DNA, these studies could be of major importance to the understanding of the mechanisms of nitrous acid cross-linking and mutagenicity, as well as the mechanisms responsible for its repair in intracellular environments. It is also the shortest DNA cross-link structure to be described.     

NMR investigation of Hoogsteen base pairing in quinoxaline antibiotic--DNA complexes: comparison of 2:1 echinomycin, triostin A and [N-MeCys3,N-MeCys7] TANDEM complexes with DNA oligonucleotides.

Hoogsteen base pairs have been demonstrated to occur in base pairs adjacent to the CpG binding sites in complexes of triostin A and echinomycin with a variety of DNA oligonucleotides. To understand the relationship of these unusual base pairs to the sequence specificity of these quinoxaline antibiotics, the conformation of the base pairs flanking the YpR binding sites of the 2:1 drug-DNA complexes of triostin A with [d(ACGTACGT)]2 and of the TpA specific [N-MeCys3, N-MeCys7] TANDEM with [d(ATACGTAT)]2 have been studied by 1H NMR spectroscopy. In both the 2:1 triostin A-DNA complex and the 2:1 [N-MeCys3, N-MeCys7] TANDEM-DNA complex, the terminal A.T base pairs are Hoogsteen base paired with the 5' adenine in the syn conformation. This indicates that both TpA specific and CpG specific quinoxaline antibiotics are capable of inducing Hoogsteen base pairs in DNA. However, in both 2:1 complexes, Hoogsteen base pairing is limited to the terminal base pairs. In the 2:1 triostin A complex, the internal adenines are anti and in the 2:1 [N-MeCys3, N-MeCys7] TANDEM-DNA complex, the internal guanines are anti regardless of pH, which indicates that the central base pairs of both complexes form Watson-Crick base pairs. This indicates that the sequence dependent nature of Hoogsteen base pairing is the same in TpA specific and CpG specific quinoxaline antibiotic-DNA complexes. We have calculated a low resolution three-dimensional structure of the 2triostin A-[d(ACGTACGT)]2 complex and compared it with other CpG specific quinoxaline antibiotic-DNA complexes. The role of stacking in the formation of Hoogsteen base pairs in these complexes is discussed.     

Solution structure of duplex DNA containing the mutagenic lesion 1,N(2)-etheno-2'-deoxyguanine

We have used high-resolution NMR spectroscopy and molecular dynamics simulations to determine the solution structure of DNA containing the genotoxic lesion 1, N (2)-etheno-2'-deoxyguanosine (epsilonG), paired to dC. The NMR data suggest the presence of a major, minimally perturbed structure at neutral pH. NOESY spectra indicate the presence of a right-handed helix with all nucleotides in anti, 2'-deoxyribose conformations within the C2'-endo/C1'-exo range and proper Watson-Crick base pair alignments outside the lesion site. The epsilonG residue remains deeply embedded inside the helix and stacks between the flanking base pairs. The lesion partner dC is extrahelical and is located in the minor groove of the duplex, where it is highly exposed to solvent. Upon acidification of the sample, a second conformation at the lesion site of the duplex emerges, with protonation of the lesion partner dC and possible formation of a Hoogsteen base pair. Restrained molecular dynamics simulations of the neutral-pH structure generated a set of three-dimensional models that show epsilonG inside the helix, where the lesion is stabilized by stacking interactions with flanking bases but without participating in hydrogen bonding. The lesion counterbase dC is displaced in the minor groove of the duplex where it can form a hydrogen bond with the sugar O4' atom of a residue 2 bp away.     

Structure of d(CGCGAATTCGCG) in the presence of Ca(2+) ions.

The dodecamer d(CGCGAATTCGCG) was the first oligonucleotide to be crystallized as a B-DNA duplex. Its structure was analyzed in detail in the early 1980s. Here we show that, in the presence of Ca(2+), it crystallizes in a different way (R3 space group). The dodecamers form parallel columns of straight duplexes with ten base pairs in the B form. The terminal cytosines in each molecule are disordered, whereas the terminal guanines are placed in the minor groove of neighbor duplexes. The central GAATTC region is practically identical to that found in the classic structure of the same dodecamer crystallized in the P2(1)2(1)2(1) space group in the presence of Mg(2+) and spermine. Its structure is thus independent of the crystallization conditions which have been used.     

Conformational transitions in thymidine bulge-containing deoxytridecanucleotide duplexes. Role of flanking sequence and temperature in modulating the equilibrium between looped out and stacked thymidine bulge states

Structural features at extra thymidine bulge sites in DNA duplexes have been elucidated from a two-dimensional NMR analysis of through-bond and through-space connectivities in the otherwise self-complementary d(C-C-G-T-G-A-A-T-T-C-C-G-G) (GTG 13-mer) and d(C-C-G-G-A-A-T-T-C-T-C-G-G) (CTC 13-mer) duplexes in aqueous solution. These studies establish that the extra thymidine flanked by guanosines in the GTG 13-mer duplex is in a conformational equilibrium between looped out and stacked states. The looped-out state is favored at low temperature (0 degrees C), whereas the equilibrium shifts in favor of the stacked state at elevated temperatures (35 degrees C) prior to the onset of the duplex-strand transition. By contrast, the extra thymidine flanked by cytidines in the CTC 13-mer duplex is looped out independent of temperature in the duplex state. Our results demonstrate that temperature and flanking sequence modulate the equilibrium between looped-out and stacked conformations of single base thymidine bulges in DNA oligomer duplexes.     

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.     

Extrahelical adenosine stacks into right-handed DNA: solution conformation of the d(C-G-C-A-G-A-G-C-T-C-G-C-G) duplex deduced from distance geometry analysis of nuclear Overhauser effect spectra

This paper reports on features of the three-dimensional structure of the d(C-G-C-A-G-A-G-C-T-C-G-C-G) self-complementary duplex (designated adenosine 13-mer), which contains symmetrical extrahelical adenosines in the interior of the helix. The majority of the protons have been assigned from two-dimensional nuclear Overhauser effect (NOESY) spectra of the adenosine 13-mer in H2O and D2O solution. The measurement of NOESY cross-peak volume integrals as a function of mixing time has yielded a set of 96 short (less than 4.5-A) proton-proton distances defined by lower and upper bounds, which have served as input parameters for a distance geometry analysis of one symmetric half of the adenosine 13-mer duplex. We demonstrate that the extrahelical adenosine stacks into the duplex for all refined structures without disruption of base pairing on either side of the modification site. The distance geometry refinement yields two classes of conformations consistent with distance measurements but which differ in orientation of the stacked extrahelical adenosine at the modification site.     

Dynamics of mismatched base pairs in DNA

The structural dynamics of mismatched base pairs in duplex DNA have been studied by time-resolved fluorescence anisotropy decay measurements on a series of duplex oligodeoxynucleotides of the general type d[CGG(AP)GGC].d[GCCXCCG], where AP is the fluorescent adenine analogue 2-aminopurine and X = T, A, G, or C. The anisotropy decay is caused by internal rotations of AP within the duplex, which occur on the picosecond time scale, and by overall rotational diffusion of the duplex. The correlation time and angular range of internal rotation of AP vary among the series of AP.X mismatches, showing that the native DNA bases differ in their ability to influence the motion of AP. These differences are correlated with the strength of base-pairing interactions in the various AP.X mismatches. The interactions are strongest with X = T or C. The ability to discern differences in the strength of base-pairing interactions at a specific site in DNA by observing their effect on the dynamics of base motion is a novel aspect of the present study. The extent of AP stacking within the duplex is also determined in this study since it influences the excited-state quenching of AP. AP is thus shown to be extrahelical in the AP.G mismatch. The association state of the AP-containing oligodeoxynucleotide strand is determined from the temperature-dependent tumbling correlation time. An oligodeoxynucleotide triplex is formed with a particular base sequence in a pH-dependent manner.     

The hydrogen-bonding structure in parallel-stranded duplex DNA is reverse Watson-Crick

Raman spectra of the parallel-stranded duplex formed from the deoxyoligonucleotides 5'-d-[(A)10TAATTTTAAATATTT]-3' (D1) and 5'-d[(T)10ATTAAAATTTATAAA]-3' (D2) in H2O and D2O have been acquired. The spectra of the parallel-stranded DNA are then compared to the spectra of the antiparallel double helix formed from the deoxyoligonucleotides D1 and 5'-d(AAATATTTAAAATTA-(T)10]-3' (D3). The Raman spectra of the antiparallel-stranded (aps) duplex are reminiscent of the spectra of poly[d(A)].poly[d(T)] and a B-form structure similar to that adopted by the homopolymer duplex is assigned to the antiparallel double helix. The spectra of the parallel-stranded (ps) and antiparallel-stranded duplexes differ significantly due to changes in helical organization, i.e., base pairing, base stacking, and backbone conformation. Large changes observed in the carbonyl stretching region (1600-1700 cm-1) implicate the involvement of the C(2) carbonyl of thymine in base pairing. The interaction of adenine with the C(2) carbonyl of thymine is consistent wtih formation of reverse Watson-Crick base pairing in parallel-stranded DNA. Phosphate-furanose vibrations similar to those observed for B-form DNA of heterogenous sequence and high A,T content are observed at 843 and 1092 cm-1 in the spectra of the parallel-stranded duplex. The 843-cm-1 band is due to the presence of a sizable population of furanose rings in the C2'-endo conformation. Significant changes observed in the regions from 1150 to 1250 cm-1 and from 1340 to 1400 cm-1 in the spectra of the parallel-stranded duplex are attributed to variations in backbone torsional and glycosidic angles and base stacking.     

Thermodynamics of the various forms of the dodecamer d(ATTACCGGTAAT) and of its constituent hexamers from proton nmr chemical shifts and UV melting curves: three-state and four-state thermodynamic models

Chemical shifts of base and H1' protons of the single-stranded hexamers d(ATTACC) and d(GGTAAT), of the 1:1 mixtures of these complementary hexamers, and of the self-complementary dodecamer d(ATTACCGGTAAT) were measured at various temperatures in aqueous solution. Four different sample concentrations were used in the case of the dodecamer and of the mixture of the complementary hexamers; the individual hexamers were measured at two different DNA concentrations. Absorbance temperature profiles at five different NaCl concentrations were measured for the dodecamer in order to quantify the effect of the ionic strength on the duplex formation. Under suitable conditions of nucleotide concentration, temperature, and ionic strength, the dodecamer adopts either a B-DNA duplex or a hairpin-loop structure. Chemical shift vs temperature profiles, constructed for all samples, were used to obtain thermodynamic parameters either for the various stacking interactions in the single strands or for the duplex or the hairpin-loop formation. In the analysis of the duplex formation of the hexamers, a two-state approach appeared too simple, because systematic deviations were revealed. Therefore, a new three-state model (DUPSTAK) was developed. In order to investigate the magnitude of error arising from the use of the two-state approach in cases where the DUPSTAK model appears more appropriate, a series of test calculations was made. The magnitude of error in the enthalpy and in the entropy of duplex melting is found to depend linearly upon the actual melting temperature and not upon the individual delta Hd degrees and delta Sd degrees values. Thermodynamic analysis of the chemical shift vs temperature profiles in D2O solution (no added salt) yields an average Tmd value of 341 K (1M DNA) and delta Hd degrees of - 121 kJ.mol-1 for the dimer/random-coil transition of the hexamer duplex d(ATTACC).d(GGTAAT). For the duplex in equilibrium random-coil transition of the 12-mer d(ATTACCGGTAAT) an average Tmd value of 336 K (1M DNA) and delta Hd degrees of -372 kJ.mol-1 are found. The hairpin/random-coil transition of d(ATTACCGGTAAT) is characterized by a rather large delta Hh degrees value, -130 kJ.mol-1, and an average Tmh value of 304 K.     

alpha-DNA-III. Characterization by high field 1H-NMR, anti-parallel self-recognition and conformation of the unnatural hexadeoxyribonucleotides alpha-[d(CpApTpGpCpG)] and alpha-[d(CpGpCpApTpG)]. Alpha-oligodeoxynucleotides as potential cellular probes for gene control.

High field 2-D-1H-NMR techniques permitted the assignment of all non-exchangeable protons of the unnatural deoxyribonucleotides alpha-[d(CpApTpGpCpG)] and alpha-[d(CpGpCpApTpG)]. 1-D and 2-D NOESY experiments show strong H6H8-H4' dipolar interactions for all nucleotides in both sequences. These data, together with COSY and J-resolved spectra, indicate that these two alpha-oligomers adopt 3'-exo conformations of the sugar moieties in solution with anti conformations of the glycosyl linkages. Both 1H-NMR data, and hypochromocity comparison of alpha-CATGCG and beta-CATGCG, demonstrate a higher degree of base stacking in the case of the alpha-sequence. The UV hyperchromicity at 260 nm, and symmetry considerations in the imino proton NMR experiments reveal antiparallel self-recognition and duplex annealing at positions 1-4 for alpha-[d(CATGCG)] and positions 3-6 for alpha-[d(CGCATG)]. The temperature variation of the imino proton NMR signals suggests that the hydrogen bonding in self-recognition is comparable in strength with that in a beta-DNA duplex, and NOE data are in accord with Watson-Crick rather than Hoogsteen base pairing.     

Solution structure of [d(ATGAGCGAATA)]2. Adjacent G:A mismatches stabilized by cross-strand base-stacking and BII phosphate groups.

The solution structure of a rather unusual B-form duplex [d(ATGAGCGAATA)]2 has been determined using two-dimensional nuclear magnetic resonance (2D-NMR) and distance geometry methods. This sequence forms a stable ten base-pair B-form duplex with 3' overhangs and two pairs of adjacent G:A mismatches paired via a sheared hydrogen-bonding scheme. All non-exchangeable protons, including the stereo-specific H-5'S/H-5'R of the 3G and 7G residues, were assigned by 2D-NMR. The phosphorus spectrum was assigned using heteronuclear correlation with H-3' and H-4' reasonances. The complete assignments reveal several unusual nuclear Overhauser enhancements (NOEs) and unusual chemical shifts for the neighboring G:A mismatch pairs and their adjacent nucleotides. Inter-proton distances were derived from time-dependent NOEs and used to generate initial structures, which were further refined by iterative back-calculation of the two-dimensional nuclear Overhauser enhancement spectra; 22 final structures were calculated from the refined distance bounds. All these final structures exhibit fully wound helical structures with small penalty values against the refined distance bounds and small pair-wise root-mean-square deviation values (typically 0.5 A to 0.9 A). The two helical strands exchange base stacking at both of the two G:A mismatch sites, resulting in base stacking down each side rather than down each strand of the twisted duplex. Very large twist angles (77 degrees) were found at the G:A mismatch steps. All the final structures were found to have BII phosphate conformations at the adjacent G:A mismatch sites, consistent with observed downfield 31P chemical shifts and Monte-Carlo conformational search results. Our results support the hypothesis that 31P chemical shifts are related to backbone torsion angles. These BII phosphate conformations in the adjacent G:A mismatch step suggest that hydrogen bonding of the G:A pair G-NH2 to a nearby phosphate oxygen atom is unlikely. The unusual structure of the duplex may be stabilized by strong interstrand base stacking as well as intrastrand stacking, as indicated by excellent base overlap within the mismatch stacks.     

Crystal and molecular structure of a DNA fragment: d(CGTGAATTCACG)

The crystal structure of the dodecanucleotide d(CGTGAATTCACG) has been determined to a resolution of 2.7 A and refined to an R factor of 17.0% for 1532 reflections. The sequence crystallizes as a B-form double helix, with Watson-Crick base pairing. This sequence contains the EcoRI restriction endonuclease recognition site, GAATTC, and is flanked by CGT on the 5'-end and ACG on the 3'-end, in contrast to the CGC on the 5'-end and GCG on the 3'-end in the parent dodecamer d(CGCGAATTCGCG). A comparison with the isomorphous parent compound shows that any changes in the structure induced by the change in the sequence in the flanking region are highly localized. The global conformation of the duplex is conserved. The overall bend in the helix is 10 degrees. The average helical twist values for the present and the parent structures are 36.5 degrees and 36.4 degrees, respectively, corresponding to 10 base pairs per turn. The buckle at the substituted sites are significantly different from those seen at the corresponding positions in the parent dodecamer. Step 2 (GpT) is underwound with respect to the parent structure (27 degrees vs 36 degrees) and step 3 (TpG) is overwound (34 degrees vs 27 degrees). There is a spine of hydration in the narrow minor groove. The N3 atom of adenine on the substituted A10 and A22 bases are involved in the formation of hydrogen bonds with other duplexes or with water; the N3 atom of guanine on G10 and G22 bases in the parent structure does not form hydrogen bonds.     

Effects of bulge composition and flanking sequence on the kinking of DNA by bulged bases

We recently showed that bulged bases kink duplex DNA, with the degree of kinking increasing in roughly equal increments as the number of bases in the bulge increases from one to four [Hsieh, C.-H., & Griffith, J.D. (1989) Proc. Natl. Acad. Sci. U.S.A. 86, 4833-4837]. Here we have examined the kinking of DNA by single A, C, G, or T bulges with different neighboring base pairs. Synthetic 30 base pair (bp) duplex DNAs containing 2 single-base bulges spaced by 10 bp were ligated head to tail, and their electrophoretic behavior in highly cross-linked gels was examined. All bulge-containing DNAs showed marked electrophoretic retardations as compared to non-bulge-containing DNA. Regardless of the sequence of the flanking base pairs, purine bulges produced greater retardations than pyrimidine bulges. Furthermore, C and T bulges produced the same retardations as did G and A bulges. Bulged DNA containing different flanking base pairs showed marked differences in electrophoretic mobility. For C-bulged DNA, the greatest retardations were observed with G.C neighbors, the least with T.A neighbors, and an intermediate amount with a mixture of neighboring base pairs. For A-bulged DNA, the retardations were greatest with G.C neighbors, less with T.A neighbors, even less with a mixture of neighboring base pairs, and finally least with C.G neighbors. Thus flanking base pairs affect C-bulged DNA and A-bulged DNA differently, and G.C and C.G flanking base pairs were seen to have very different effects. These results imply an important role of base stacking in determining how neighboring base pairs influence the kinking of DNA by a single-base bulge.     

Effect of DNA flanking sequence on charge transport in short DNA duplexes

Several factors can influence charge transport (CT)-mediated DNA, such as sequence, distance, base stacking, base pair mismatch, conformation, tether length, etc. However, the DNA context effect or how flanking sequences influence redox active drugs in the DNA CT reaction and later in DNA enzymatic repair and synthesis is still not well understood. The set of seven DNA molecules in this study have been characterized well for the study of flanking sequence effects. These DNA duplexes are formed from self-complementary strands and contain the common central four-base sequence 5'-A-G-C-T-3', flanked on both sides by either (AT)(n) or (AA)(n) (n = 2, 3, or 4) or AA(AT)(2). UV-vis, fluorescence, UV melting, circular dichroism, and cyclic voltammetry experiments were used to study the flanking sequence effect on CT-mediated DNA by using daunomycin or adriamycin cross-linked with these seven DNA molecules. Our results showed that charge transport was related to the flanking sequence, DNA melting free energy, and ionic strength. For (AA)(n) or (AT)(n) species of the same length, (AA)(n) series were more stable and more efficient CT was observed through the (AA)(n) series. The same trend was observed for (AA)(n)() and (AT)(n) series at different ionic strengths, further supporting the idea that flanking sequence can result in different base stacking and modulate charge transport through these seven DNA molecules.     

A conserved pseudouridine modification in eukaryotic U2 snRNA induces a change in branch-site architecture.

The removal of noncoding sequences (introns) from eukaryotic precursor mRNA is catalyzed by the spliceosome, a dynamic assembly involving specific and sequential RNA-RNA and RNA-protein interactions. An essential RNA-RNA pairing between the U2 small nuclear (sn)RNA and a complementary consensus sequence of the intron, called the branch site, results in positioning of the 2'OH of an unpaired intron adenosine residue to initiate nucleophilic attack in the first step of splicing. To understand the structural features that facilitate recognition and chemical activity of the branch site, duplexes representing the paired U2 snRNA and intron sequences from Saccharomyces cerevisiae were examined by solution NMR spectroscopy. Oligomers were synthesized with pseudouridine (psi) at a conserved site on the U2 snRNA strand (opposite an A-A dinucleotide on the intron strand, one of which forms the branch site) and with uridine, the unmodified analog. Data from NMR spectra of nonexchangeable protons demonstrated A-form helical backbone geometry and continuous base stacking throughout the unmodified molecule. Incorporation of psi at the conserved position, however, was accompanied by marked deviation from helical parameters and an extrahelical orientation for the unpaired adenosine. Incorporation of psi also stabilized the branch-site interaction, contributing -0.7 kcal/mol to duplex deltaG degrees 37. These findings suggest that the presence of this conserved U2 snRNA pseudouridine induces a change in the structure and stability of the branch-site sequence, and imply that the extrahelical orientation of the branch-site adenosine may facilitate recognition of this base during splicing.     

Extrahelical bases in duplex DNA

The two non-complementary synthetic DNAs, poly[d(G-G-A)] and poly[d(T-C)] can form a duplex structure at moderate ionic strengths in which every other T is extrahelical:

Only this structure is consistent with the stoichiometry of formation and with the production of C? photodimers on ultraviolet light irradiation. Further characterization with respect to melting temperature and interaction with nuclease or damage-recognizing proteins is described.     

Deuteriation of sugar protons simplify NMR assignments and structure determination of large oligonucleotide by the 1H-NMR window approach.

The concept of the 1H-NMR window has been developed and examined through a comparative study of NOESY spectra of a self-complementary Dickerson's dodecamer (I) [5'd(5C6G7C8G9A10A11T12T13C-14G15C16G)2(3')], a self-complementary 20-mer (II) [(5'd(1C2G3C4G5C6G7C8G9A10A11T12T13C14G15C16G17C18G19C20G)2(3')] in which the core part consists of the same Dickerson's dodecamer sequence with the flanking CGCG residues at both 3' and 5'-ends, and the partly-deuteriated (shown by underlined CGCG residues at both 3' and 5'-ends) analogous duplex (III) [5'd(1C2G3C4G5C6G7C8G9A10A11T12T13C14G15C16G17C18G19C20G)2(3')] in which the core 5C to 16G part (i.e. 1H-NMR window) consists of the natural Dickerson's dodecamer sequence. A comparison of their NOESY spectra clearly demonstrates that the severe overlap of proton resonances in the larger DNA duplex (II) has been successfully reduced in the partly-deuterated duplex (III) as a result of specific incorporations of the sugar-deuteriated nucleotide residues in the latter [stereospecific > 97 atom % 2H enrichment at H2', H2' and H3' sites, approximately 85 atom % 2H enrichment at H4' and approximately 20 atom % 2H enrichment at H1' (see refs. 10 and 11) in the 20-mer duplex (III)]. These simplifications of the resonance overlap by the deuteriation approach have enabled unequivocal chemical shift assignments and extraction of the quantitative NOE data in the 1H-NMR window part of duplex (III). A comparison of the 12-nucleotide long 1H-NMR window in (III) with that of the 12-mer duplex (I) also shows the scope of studying the changes in conformation and dynamics of the core 12-mer region in (III) which result from the increase of molecular weight due to the DNA chain extension. It is noteworthy that such a study is clearly impossible for the natural 20-mer (II) because of the inherent problem of the overlap of resonances.     

Thermodynamics of DNA duplexes with adjacent G.A mismatches.

The sequence 5'-d(ATGAGCGAAT) forms a very stable self-complementary duplex with four G.A mismatch base pairs (underlined) out of ten total base pairs [Li et al. (1991) Proc. Natl. Acad. Sci. U.S.A. 88, 26-30]. The conformation is in the general B-family and is stabilized by base-pair hydrogen bonding of an unusual type, by favorable base dipole orientations, and by extensive purine-purine stacking at the mismatched sites. We have synthesized 13 decamers with systematic variations in the sequence above to determine how the flanking sequences, the number of G.A mismatches, and the mismatch sequence order (5'-GA-3' or 5'-AG-3') affect the duplex stability. Changing A.T to G.C base pairs in sequences flanking the mismatches stabilizes the duplexes, but only to the extent observed with B-form DNA. The sequence 5'-pyrimidine-GA-purine-3', however, is considerably more stable than 5'-purine-GA-pyrimidine-3'. The most stable sequences with two pairs of adjacent G.A mismatches have thermodynamic parameters for duplex formation that are comparable to those for fully Watson-Crick base-paired duplexes. Similar sequences with single G.A pairs are much less stable than sequences with adjacent G.A mismatches. Reversing the mismatch order from 5'-GA-3' to 5'-AG-3' results in an oligomer that does not form a duplex. These results agree with predictions from the model derived from NMR and molecular mechanics and indicate that the sequence 5'-pyrimidine-GA-purine-3' forms a stable conformational unit that fits quite well into a B-form double helix.(ABSTRACT TRUNCATED AT 250 WORDS)     

Oxetane locked thymidine in the Dickerson-Drew dodecamer causes local base pairing distortions -- an NMR structure and hydration study

The introduction of a North-type sugar conformation constrained oxetane T block, 1-(1',3'-O-anhydro-beta-D-psicofuranosyl) thymine, at the T(7) position of the self-complementary Dickerson-Drew dodecamer, d[(5'-C(1)G(2)C(3)G(4)A(5)A(6)T(7)T(8)C(9)G(10)C(11)G(12)-3')](2), considerably perturbs the conformation of the four central base pairs, reducing the stability of the structure. UV spectroscopy and 1D NMR display a drop in melting temperature of approximately 10 degrees C per modification for the T(7) oxetane modified duplex, where the T(7) block has been introduced in both strands, compared to the native Dickerson-Drew dodecamer. The three dimensional structure has been determined by NMR spectroscopy and has subsequently been compared with the results of 2.4 ns MD simulations of the native and the T(7) oxetane modified duplexes. The modified T(7) residue is found to maintain its constrained sugar- and the related glycosyl torsion conformations in the duplex, resulting in staggered and stretched T(7).A(6) and A(6).T(7) non-linear base pairs. The stacking is less perturbed, but there is an increased roll between the two central residues compared to the native counterpart, which is compensated by tilts of the neighboring base steps. The one dimensional melting profile of base protons of the T(7) and T(8) residues reveals that the introduction of the North-type sugar constrained thymine destabilizes the core of the modified duplex, promoting melting to start simultaneously from the center as well as from the ends. Temperature dependent hydration studies by NMR demonstrate that the central T(7).A(6)/A(6).T(7) base pairs of the T(7) oxetane modified Dickerson-Drew dodecamer have at least one order of magnitude higher water exchange rates (correlated to the opening rate of the base pair) than the corresponding base pairs in the native duplex.