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

Exploration of the Structure and Recognition of a G-quadruplex in the her2 Proto-oncogene Promoter and Its Transcriptional Regulation

G-quadruplexes in oncogene promoters provide putative targets for transcriptional regulation. The structure of a putative G-quadruplex sequence (S1: GGAGAAGGAGGAGGTGGAGGAGGAGGG) in potassium solution in the her2 promoter has been resolved mainly through nuclear magnetic resonance (NMR) spectroscopy. By application of various NMR spectra, we proved the formation of a four-layer G-quadruplex composing of two G-tetrads and two G/A-mixed planes with a four-residues loop (A3-G4-A5-A6). Further evidence from a luciferase reporter assay, Q-RT-PCR and Western blotting indicates that S1 G-quadruplex formation can repress her2 promoter activity, and a selected G-quadruplex ligand cβ can enhance the repression by down regulating her2 transcription and expression. These findings provide a G-quadruplex target and perspective implications in her2 transcriptional regulation.

Sequence-dependent structural changes in a self-assembling DNA oligonucleotide

DNA has proved to be a remarkable molecule for the construction of sophisticated two-dimensional and three-dimensional architectures because of its programmability and structural predictability provided by complementary Watson-Crick base pairing. DNA oligonucleotides can, however, exhibit a great deal of local structural diversity. DNA conformation is strongly linked to both environmental conditions and the nucleobase identities inherent in the oligonucleotide sequence, but the exact relationship between sequence and local structure is not completely understood. This study examines how a single-nucleotide addition to a class of self-assembling DNA 13-mers leads to a significantly different overall structure under identical crystallization conditions. The DNA 13-mers self-assemble in the presence of Mg(2+) through a combination of Watson-Crick and noncanonical base-pairing interactions. The crystal structures described here show that all of the predicted Watson-Crick base pairs are present, with the major difference being a significant rearrangement of noncanonical base pairs. This includes the formation of a sheared A-G base pair, a junction of strands formed from base-triple interactions, and tertiary interactions that generate structural features similar to tandem sheared G-A base pairs. The adoption of this alternate noncanonical structure is dependent in part on the sequence in the Watson-Crick duplex region. These results provide important new insights into the sequence-structure relationship of short DNA oligonucleotides and demonstrate a unique interplay between Watson-Crick and noncanonical base pairs that is responsible for crystallization fate.

Automatic workflow for the classification of local DNA conformations

Background

A growing number of crystal and NMR structures reveals a considerable structural polymorphism of DNA architecture going well beyond the usual image of a double helical molecule. DNA is highly variable with dinucleotide steps exhibiting a substantial flexibility in a sequence-dependent manner. An analysis of the conformational space of the DNA backbone and the enhancement of our understanding of the conformational dependencies in DNA are therefore important for full comprehension of DNA structural polymorphism.

Results

A detailed classification of local DNA conformations based on the technique of Fourier averaging was published in our previous work. However, this procedure requires a considerable amount of manual work. To overcome this limitation we developed an automatic classification method consisting of the combination of supervised and unsupervised approaches. A proposed workflow is composed of k-NN method followed by a non-hierarchical single-pass clustering algorithm. We applied this workflow to analyze 816 X-ray and 664 NMR DNA structures released till February 2013. We identified and annotated six new conformers, and we assigned four of these conformers to two structurally important DNA families: guanine quadruplexes and Holliday (four-way) junctions. We also compared populations of the assigned conformers in the dataset of X-ray and NMR structures.

Conclusions

In the present work we developed a machine learning workflow for the automatic classification of dinucleotide conformations. Dinucleotides with unassigned conformations can be either classified into one of already known 24 classes or they can be flagged as unclassifiable. The proposed machine learning workflow permits identification of new classes among so far unclassifiable data, and we identified and annotated six new conformations in the X-ray structures released since our previous analysis. The results illustrate the utility of machine learning approaches in the classification of local DNA conformations.

Moving beyond Watson-Crick models of coarse grained DNA dynamics

DNA produces a wide range of structures in addition to the canonical B-form of double-stranded DNA. Some of these structures are stabilized by Hoogsteen bonds. We developed an experimentally parameterized, coarse-grained model that incorporates such bonds. The model reproduces many of the microscopic features of double-stranded DNA and captures the experimental melting curves for a number of short DNA hairpins, even when the open state forms complicated secondary structures. We demonstrate the utility of the model by simulating the folding of a thrombin aptamer, which contains G-quartets, and strand invasion during triplex formation. Our results highlight the importance of including Hoogsteen bonding in coarse-grained models of DNA.

Structures of (5'S)-8,5'-Cyclo-2'-deoxyguanosine Mismatched with dA or dT

Diastereomeric 8,5'-cyclopurine 2'-deoxynucleosides, containing a covalent bond between the deoxyribose and the purine base, are induced in DNA by ionizing radiation. They are suspected to play a role in the etiology of neurodegeneration in xeroderma pigmentosum patients. If not repaired, the S-8,5'-cyclo-2'-deoxyguanosine lesion (S-cdG) induces Pol V-dependent mutations at a frequency of 34% in Escherichia coli. Most are S-cdG → A transitions, suggesting mis-incorporation of dTTP opposite the lesion during replication bypass, although low levels of S-cdG → T transversions, arising from mis-incorporation of dATP, are also observed. We report the structures of 5'-d(GTGCXTGTTTGT)-3'·5'-d(ACAAACAYGCAC)-3', where X denotes S-cdG and Y denotes either dA or dT, corresponding to the situation following mis-insertion of either dTTP or dATP opposite the S-cdG lesion. The S-cdG·dT mismatch pair adopts a wobble base pairing. This provides a plausible rationale for the S-cdG → A transitions. The S-cdG·dA mismatch pair differs in conformation from the dG·dA mismatch pair. For the S-cdG·dA mismatch pair, both S-cdG and dA intercalate, but no hydrogen bonding is observed between S-cdG and dA. This is consistent with the lower levels of S-cdG → T transitions in E. coli.

Tandem GA residues on opposite sides of the loop in molecular beacon-like DNA hairpins compact the loop and increase hairpin stability

The free solution electrophoretic mobilities and thermal stabilities of hairpins formed by two complementary 26-nucleotide oligomers have been measured by capillary electrophoresis. The oligomers are predicted to form molecular beacon-like hairpins with 5 bp stems and 16 nucleotides in the loop. One hairpin, called hairpin2 (hp2), migrates with a relatively fast free solution mobility and exhibits melting temperatures that are reasonably well predicted by the popular structure-prediction program Mfold. Its complement, called hairpin1 (hp1), migrates with a slower free solution mobility and forms a stable hairpin only in solutions containing ≥200 mM Na(+). The melting temperatures observed for hp1 are ~18 °C lower than those observed for hp2 and ~20 °C lower than those predicted by Mfold. The greater thermal stability of hp2 is due to the presence of tandem GA residues on opposite sides of the loop. If the corresponding TC residues in the hp1 loop are replaced by tandem GA residues, the melting temperatures of the modified hairpin are close to those observed for hp2. Eliminating the tandem GA residues in the hp2 loop significantly decreases the thermal stability of hp2. If the loops are replaced by a loop of 16 thymine residues, the free solution mobilities and thermal stabilities of the T-loop hairpin are equal to those observed for hp1. Hence, the loop of hp1 appears to be relatively unstructured, with few base-base stacking interactions. Interactions between tandem GA residues on opposite sides of the hp2 loop appear to compact the loop and increase hairpin stability.

Structural properties of DNA oligomers containing (GACX)(n) and (GAXC)(n) tandem repeats

The antisense DNA sequence of mature mouse micro RNA, miR341, includes three repeats of the tetranucleotide (GACC). The -GAC- repeat is known to form a parallel duplex, in acidic environments. The thermal melting profile of miR341 DNA, at pH 4, 5, and 6 indicates the formation of a very stable structure, which loses its stability when pH is increased. Thus, the addition of a cytosine at the 3' end of the (GAC) motif preserves the molecule's potential to fold into an unusual structure at low pH. The effect of modifying the nucleotide composition of the GACC sequence on the secondary structures formed by oligomers containing seven tandem repeats of the altered motifs was examined here. UV melting profiles were determined, as a function of pH, for 28-mers of the two series (GAXC)(7) and (GACX)(7) (X= A/C/T/G)(.) The sequence (GACC)(7) was found to be extremely sensitive to pH variations, with a stable structure formed at pH 5 (T(m) ≥ 60°C). NMR spectroscopy established that the low pH structure is not B-DNA. (GACA)(7) and (GACT)(7) also formed stable structures at low pH but the addition of guanine at the 3'end, as seen in the (GACG) series resulted in the loss of this property. Introducing a break in the 5'-GAC-3' motif, explored in the (GAXC) series, also inhibits formation of stable structures under acidic conditions.

Study on intercalations between double-stranded DNA and pyrene by single-molecule force spectroscopy: toward the detection of mismatch in DNA

Intercalation interactions between planar aromatic molecules and double-stranded DNA (dsDNA) relate to not only the structure of the guest molecules but also the structure of the dsDNA. In this letter, we have comparatively studied the intercalation between pyrene and fully matched or mismatched dsDNA using single-molecule force spectroscopy (SMFS). The significant difference in rupture forces, upon pyrene unbinding from 25-mer dsDNA with or without mismatches, is observed at the single-molecule level, indicating the influence of mismatches on the interaction between pyrene and dsDNA. In the analysis of the dynamic force spectra, two transition barriers are revealed for pyrene unbinding from matched sites in dsDNA and for pyrene unbinding from mismatched sites as well. These results suggest that SMFS is a useful single-molecule method for the detection of mismatches in dsDNA by the intercalation of pyrene.

DNA structures from phosphate chemical shifts

For B-DNA, the strong linear correlation observed by nuclear magnetic resonance (NMR) between the (31)P chemical shifts (deltaP) and three recurrent internucleotide distances demonstrates the tight coupling between phosphate motions and helicoidal parameters. It allows to translate deltaP into distance restraints directly exploitable in structural refinement. It even provides a new method for refining DNA oligomers with restraints exclusively inferred from deltaP. Combined with molecular dynamics in explicit solvent, these restraints lead to a structural and dynamical view of the DNA as detailed as that obtained with conventional and more extensive restraints. Tests with the Jun-Fos oligomer show that this deltaP-based strategy can provide a simple and straightforward method to capture DNA properties in solution, from routine NMR experiments on unlabeled samples.

Structure-specific binding of [Co(phen)(2)(HPIP)](3+) to a DNA duplex containing sheared G:A mismatch base pairs

The binding of a Co(III) complex to the decanucleotide d(CCGAATGAGG)(2) containing two pairs of G:A mismatches was studied by 2D-NMR, UV absorption, and molecular modeling. NMR investigations indicate that racemic [Co(phen)(2)(HPIP)]Cl(3) [HPIP=2-(2-hydroxyphenyl) imidazo [4,5-f][1,10] phenanthroline] binds the decanucleotide by intercalation: the HPIP ligand selectively inserts between the stacked bases from the minor groove at the terminal regions and from the major groove at the sheared region. Further, molecular modeling revealed that the recognition shows strong enantioselectivity: the Lambda-isomer preferentially intercalates into the T(6)G(7):A(5)A(4) region from the DNA major groove, while Delta-isomer favors the terminal C(1)C(2):G(10)G(9) region and intercalates from the minor groove. Detailed energy analysis suggests that the steric interaction, especially the electrostatic effect, is the primary determinants of the recognition event. Melting experiments indicate that binding stabilizes the DNA duplex and increases the melting temperature by 9.5 degrees C. The intrinsic binding constant of the complex to the mismatched duplex was determined to be 3.5x105M(-1).

Molecular dynamics of a kappaB DNA element: base flipping via cross-strand intercalative stacking in a microsecond-scale simulation

The sequence-dependent structural variability and conformational dynamics of DNA play pivotal roles in many biological milieus, such as in the site-specific binding of transcription factors to target regulatory elements. To better understand DNA structure, function, and dynamics in general, and protein···DNA recognition in the ‘κB’ family of genetic regulatory elements in particular, we performed molecular dynamics simulations of a 20-bp DNA encompassing a cognate κB site recognized by the proto-oncogenic ‘c-Rel’ subfamily of NF-κB transcription factors. Simulations of the κB DNA in explicit water were extended to microsecond duration, providing a broad, atomically detailed glimpse into the structural and dynamical behavior of double helical DNA over many timescales. Of particular note, novel (and structurally plausible) conformations of DNA developed only at the long times sampled in this simulation—including a peculiar state arising at ≈0.7 μs and characterized by cross-strand intercalative stacking of nucleotides within a longitudinally sheared base pair, followed (at ≈1 μs) by spontaneous base flipping of a neighboring thymine within the A-rich duplex. Results and predictions from the microsecond-scale simulation include implications for a dynamical NF-κB recognition motif, and are amenable to testing and further exploration via specific experimental approaches that are suggested herein.

C-->G base mutations in the CArG box of c-fos serum response element alter its bending flexibility. Consequences for core-SRF recognition

By binding to the CArG box sequence, the serum response factor (SRF) activates several muscle-specific genes, as well as genes that respond to mitogens. The core domain of the SRF (core-SRF) binds as a dimer to the CArG box C-5C-4A-3T-2A-1T+1T+2A+3G+4G+5 of the c-fos serum response element (SREfos). However, previous studies using 20-mer DNAs have shown that the binding stoichiometry of core-SRF is significantly altered by mutations C-5-->G (SREGfos) and C-5C-4-->GG (SREGGfos) of the CArG box [A Huet, A Parlakian, M-C Arnaud, J-M Glandières, P Valat, S Fermandjian, D Paulin, B Alpert & C Zentz (2005) FEBS J272, 3105-3119]. To understand these effects, we carried out a comparative analysis of the three 20-mer DNAs SREfos, SREGfos and SREGGfos in aqueous solution. Their CD spectra were of the B-DNA type with small differences generated by variations in the mutual arrangement of the base pairs. Analysis by singular value decomposition of a set of Raman spectra recorded as a function of temperature, revealed a premelting transition associated with a conformational shift in the DNA double helices from a bent to a linear form. Time-resolved fluorescence anisotropy shows that the fluorescein reporter linked to the oligonucleotide 5'-ends experiences twisting motions of the double helices related to the interconversion between bent and linear conformers. The three SREs present various bent populations submitted, however, to particular internal dynamics, decisive for the mutual adjustment of binding partners and therefore specific complex formation.

Z-DNA's conformer substates revealed by FT-IR difference spectroscopy of nonoriented left-handed double helical poly(dG-dC)

Nonoriented hydrated films of double helical poly(dG-dC) in the Z-form were studied by Fourier transform infrared (FT-IR) spectroscopy either as equilibrated slow-cooled samples between 290 and 220 K or, after quenching into the glassy state, as nonequilibrated film isothermally at 200, 220, and 240 K. IR spectral changes on isothermal relaxation at 200 and 220 K toward equilibrium, caused by interconversion of two conformer substates (CS) called Z1 and Z2, are revealed by IR difference spectra. Pronounced spectral changes on Z1-to-Z2 interconversion occur between approximately 750-1250 cm(-1) and these are attributed to structural changes of the phosphodiester-sugar backbone caused by changes of torsion angles, and to decreasing hydrogen-bonding of the ionic phosphate group with water molecules. These spectral changes on Z1-to-Z2 transition can be related to structural differences between ZI and ZII CS observed in single crystals. ZI/ZII CS occurs only at (dGpdC) base steps, and similar behavior is assumed for Z1/Z2. The Z1/Z2 population ratio was determined via curve resolution of marker bands for Z1 and Z2 centered at 785 and 779 cm(-1). This ratio is 0.64 at 290 K, corresponding to 39% of the phosphates of the (dGpdC) base steps in Z1 and 61% in Z2, and it increases to 1.24 on cooling to 220 K. For the Z2<=>Z1 equilibrium, an enthalpy change of -4.9 +/- 0.2 kJ mol(dGpdC)(-1) is obtained from the temperature dependence of the equilibrium constant. Z1 interconverts into Z2 at isothermal relaxation at 200 and 220 K, whereas on slow cooling from ambient temperature, Z2 interconverts into Z1. This unexpected reversal of CS interconversion is attributed to slow restructuring of hydration shells of the CS on quenching, in the same manner reported by Pichler et al. for the BI and BII CS of B-DNA (J. Phys. Chem. B 106, 3263-3274 (2002)). IR difference curves demonstrate two time scales on isothermal relaxation of Z1-->Z2 interconversion, a fast one for structural relaxation of the sugar-phosphate backbone, and a slow one for relaxation of the hydration shells. This slowing down of restructuring of CS hydration shells at approximately 220-240 K could be the cause for the suppression of biological functions at low temperatures.

Structural study of DNA duplex containing an N-(2-deoxy-beta-D-erythro-pentofuranosyl) formamide frameshift by NMR and restrained molecular dynamics

The presence of an N-(2-deoxy-beta-D-erythro-pentofuranosyl) formamide (F) residue, a ring fragmentation product of thymine, in a frameshift context in the sequence 5'-d-(AGGACCACG)*d(CGTGGFTCCT) has been studied by 1H and 31P nuclear magnetic resonance (NMR) and molecular dynamics. Two-dimensional NMR studies show that the formamide residue, whether the cis or trans isomer, is rotated out of the helix and that the bases on either side of the formamide residue in the sequence, G14 and T16, are stacked over each other in a way similar to normal B-DNA. The cis and trans isomers were observed in the ratio 3:2 in solution. Information extracted from 31P NMR data reveal a modification of the phosphodiester backbone conformation at the extrahelical site, which is also observed during the molecular dynamics simulations.

Sheared-type G(anti).C(syn) base-pair: a unique d(GXC) loop closure motif

Stable DNA loop structures closed by a novel G.C base-pair have been determined for the single-residue d(GXC) loops (X=A, T, G or C) in low-salt solution by high-resolution nuclear magnetic resonance (NMR) techniques. The closing G.C base-pair in these loops is not of the canonical Watson-Crick type, but adopts instead a unique sheared-type (trans Watson-Crick/sugar-edge) pairing, like those occurring in the sheared mismatched G.A or A.C base-pair, to draw the two opposite strands together. The cytidine residue in the closing base-pair is transformed into the rare syn domain to form two H-bonds with the guanine base and to prevent the steric clash between the G 2NH(2) and the C H-5 protons. Besides, the sugar pucker of the syn cytidine is still located in the regular C2'-endo domain, unlike the C3'-endo domain adopted for the pyrimidines of the out-of-alternation left-handed Z-DNA structure. The facile formation of the compact d(GXC) loops closed by a unique sheared-type G(anti).C(syn) base-pair demonstrates the great potential of the single-stranded d(GXC) triplet repeats to fold into stable hairpins.

Unusual DNA duplex and hairpin motifs

Single-stranded DNA or double-stranded DNA has the potential to adopt a wide variety of unusual duplex and hairpin motifs in the presence (trans) or absence (cis) of ligands. Several principles for the formation of those unusual structures have been established through the observation of a number of recurring structural motifs associated with different sequences. These include: (i) internal loops of consecutive mismatches can occur in a B-DNA duplex when sheared base pairs are adjacent to each other to confer extensive cross- and intra-strand base stacking; (ii) interdigitated (zipper-like) duplex structures form instead when sheared G*A base pairs are separated by one or two pairs of purine*purine mismatches; (iii) stacking is not restricted to base, deoxyribose also exhibits the potential to do so; (iv) canonical G*C or A.T base pairs are flexible enough to exhibit considerable changes from the regular H-bonded conformation. The paired bases become stacked when bracketed by sheared G.A base pairs, or become extruded out and perpendicular to their neighboring bases in the presence of interacting drugs; (v) the purine-rich and pyrimidine-rich loop structures are notably different in nature. The purine-rich loops form compact triloop structures closed by a sheared G*A, A*A, A*C or sheared-like G(anti)*C(syn) base pair that is stacked by a single residue. On the other hand, the pyrimidine-rich loops with a thymidine in the first position exhibit no base pairing but are characterized by the folding of the thymidine residue into the minor groove to form a compact loop structure. Identification of such diverse duplex or hairpin motifs greatly enlarges the repertoire for unusual DNA structural formation.

Cleavage of fragments containing DNA mismatches by enzymic and chemical probes

We prepared synthetic 50-mer DNA duplexes, each containing four mismatched base-pairs in similar positions. We examined their cleavage by DNases I and II, micrococcal nuclease (MNase), methidiumpropyl-EDTA-Fe(II) [MPE-Fe(II)] and hydroxyl radicals. We find that single mismatches only produce subtle changes in the DNase I-cleavage pattern, the most common of which is attenuated cleavage at locations 2-3 bases on the 3'-side of the mismatch. Subtle changes are also observed in most of the DNase II-cleavage patterns, although GT and GG inhibit the cleavage over longer regions and generate patterns that resemble footprints. MNase cleaves the heteroduplexes at the mismatches themselves (except for CC), and in some cases cleaves CpG and CpC steps. None of the mismatches causes any change in the cleavage patterns produced by hydroxyl radicals or MPE-Fe(II). We also examined the cleavage patterns of fragments containing tandem GA mismatches in the sequences RGAY/RGAY and YGAR/YGAR (R, purine; Y, pyrimidine). RGAY causes only subtle changes in the cleavage patterns, which are similar to those seen with single mismatches, except that there are no changes in MNase cleavage. However, YGAR inhibits DNases I and II cleavage over 4-6 bases, and attenuates MPE-Fe(II) and hydroxyl radical cleavage at 2 bases. These changes suggest that this mismatch has a more pronounced effect on the local DNA structure. These changes are discussed in terms of the structural and dynamic effects of each mismatch.

The conformer substates of nonoriented B-type DNA in double helical poly(dG-dC)

A nonoriented hydrated film of poly(dG-dC) with ?20 water molecules per nucleotide (called B* by Loprete and Hartman (Biochem. 32, 4077-4082 (1993)) was studied by Fourier transform infrared (FT-IR) spectroscopy either as equilibrated sample between 290 and 270 K or, after quenching into the glassy state, as nonequilibrated film isothermally at 200 and 220 K. IR spectral changes on isothermal relaxation at 200 and 220 K, caused by interconversion of two conformer substates, are revealed by difference spectra. Comparison with difference curves obtained in the same manner from two classical B-DNA forms, namely the d(CGCGAATTCGCG)(2) dodecamer and polymeric NaDNA from salmon testes, revealed that the spectral changes on B(I)-to-B(II) interconversion in the classical B-DNA forms are very similar to those in the B*-form, and that the spectroscopic differences between the B(I) and B(II) features from classical B-DNA and those from the modified B*-form are minor. Nonexponential kinetics of the B(I)-->B(II) transition in the B*-form of poly(dG-dC) at 200 K showed that the structural relaxation time is about three times of that in the classical B-DNA forms (approximately equal to 30 versus approximately equal to 10 min at 200 K). The unexpected reversal of conformer substates interconversion (that is B(II)-->B(I) transition on cooling from 290 K and B(I)-->B(II) transition on isothermal relaxation at 200 K) observed for classical B-DNA occurs also in the modified B*-form. We therefore conclude that restructuring of hydration shells rules the low-temperature dynamics of the B*-form via its two conformer substates in the same manner reported for classical B-DNA by Pichler et al. (J. Phys. Chem. B 106, 3263-3274 (2002)).

1H NMR determination of base-pair lifetimes in oligonucleotides containing single base mismatches

Proton nuclear magnetic resonance (NMR) spectroscopy is employed to characterize the kinetics of base-pair opening in a series of 9mer duplexes containing different single base mismatches. The imino protons from the different mismatched, as well as fully matched, duplexes are assigned from the imino-imino region in the WATERGATE NOESY spectra. The exchange kinetics of the imino protons are measured from selective longitudinal relaxation times. In the limit of infinite exchange catalyst concentration, the exchange times of the mismatch imino protons extrapolate to much shorter lifetimes than are commonly observed for an isolated GC base pair. Different mismatches exhibit different orders of base-pair lifetimes, e.g. a TT mismatch has a shorter base-pair lifetime than a GG mismatch. The effect of the mismatch was observed up to a distance of two neighboring base pairs. This indicates that disruption in the duplex caused by the mismatch is quite localized. The overall order of base-pair lifetimes in the selected sequence context of the base pair is GC > GG > AA > CC > AT > TT. Interestingly, the fully matched AT base pair has a shorter base-pair lifetime relative to many of the mismatches. Thus, in any given base pair, the exchange lifetime can exhibit a strong dependence on sequence context. These findings may be relevant to the way mismatch recognition is accomplished by proteins and small molecules.

Quadruple intercalated G-6 stack: a possible motif in the fold-back structure of the Drosophila centromeric dodeca-satellite?

The purine-rich strand d(GTACGGGACCGA)(n) of the Drosophila centromeric dodeca-satellite sequence is highly conserved and was found to form stable fold-back structures in which the homopurine 5'-GGGA-3' sequence was determined to play a crucial role. Here, we report the stable formation of the d(GGGA)(2) motif in the stem of a DNA hairpin closed by a single-residue d(ACC) loop. Similar to the zipper-like d(GGA)(2) motif observed in the human centromeric (TGGAA)(n) sequence, the central four guanosine bases in the d(GGGA)(2) motif do not pair, but interdigitate to form an elongated zipper-like quadruple-intercalated G-6 stack bracketed by sheared G.A base-pairs. Comparison between the current d(GGGA)(2) structure and the published crystal d(GAAA)(2) structure implies that the alignment of the unpaired purine bases plays an important role in determining the minor groove width of the purine-rich d(GPuPuA)(2) motif. Similarity between the zipper-like motifs possibly present in the Drosophila centromeric dodeca-satellite sequence and in the human centromeric (TGGAA)(n) sequence led us to propose that these special zipper-like motifs may constitute common cores in organizing eukaryotic centromeres.

An intramolecular quadruplex of (GGA)(4) triplet repeat DNA with a G:G:G:G tetrad and a G(:A):G(:A):G(:A):G heptad, and its dimeric interaction

The structure of d(GGAGGAGGAGGA) containing four tandem repeats of a GGA triplet sequence has been determined under physiological K(+) conditions. d(GGAGGAGGAGGA) folds into an intramolecular quadruplex composed of a G:G:G:G tetrad and a G(:A):G(:A):G(:A):G heptad. Four G-G segments of d(GGAGGAGGAGGA) are aligned parallel with each other due to six successive turns of the main chain at each of the GGA and GAGG segments. Two quadruplexes form a dimer stabilized through a stacking interaction between the heptads of the two quadruplexes. Comparison of the structure of d(GGAGGAGGAGGA) with the reported structure of d(GGAGGAN) (N=G or T) containing two tandem repeats of the GGA triplet revealed that although the two structures resemble each other to some extent, the extension of the repeats of the GGA triplet leads to distinct structural differences: intramolecular quadruplex for 12-mer versus intermolecular quadruplex for 7-mer; heptad versus hexad in the quadruplex; and three sheared G:A base-pairs versus two sheared G:A base-pairs plus one A:A base-pair per quadruplex. It was also suggested that d(GGAGGAGGAGGA) forms a similar quadruplex under low salt concentration conditions. This is in contrast to the case of d(GGAGGAN) (N=G or T), which forms a duplex under low salt concentration conditions. On the basis of these results, the structure of naturally occurring GGA triplet repeat DNA is discussed.

Zipper-like Watson-Crick base-pairs

A series of DNA heptadecamers containing the DNA analogues of RNA E-like 5'-d(GXA)/(AYG)-5' motifs (X/Y is complementary T/A, A/T, C/G, or G/C pair) were studied using nuclear magnetic resonance (NMR) methodology and distance geometry (DG)/molecular dynamics (MD) approaches. Such oligomers reveal excellent resolution in NMR spectra and exhibit many unusual nuclear Overhauser effects (NOEs) that allow for good characterization of an unusual zipper-like conformation with zipper-like Watson-Crick base-pairs; the potential canonical X.Y H-bonding is not present, and the central X/Y pairs are transformed instead into inter-strand stacks that are bracketed by sheared G.A base-pairs. Such phenomenal structural change is brought about mainly through two backbone torsional angle adjustments, i.e. delta from C2'-endo to C3'-endo for the sugar puckers of unpaired residues and gamma from gauche(+) to trans for the following 3'-adenosine residues. Such motifs are analogous to the previously studied (GGA)(2) motif presumably present in the human centromeric (TGGAA)(n) tandem repeat sequence. The novel zipper-like motifs are only 4-7 deg. C less stable than the (GGA)(2) motif, suggesting that inter-strand base stacking plays an important role in stabilizing unusual nucleic acid structures. The discovery that canonical Watson-Crick G.C or A.T hydrogen-bonded pairs can be transformed into stacking pairs greatly increases the repertoire for unusual nucleic acid structural motifs.

Impact of CpG methylation on structure, dynamics and solvation of cAMP DNA responsive element

Methylation of CpG motifs in DNA is involved in the control of gene expression and in several other epigenic effects. It suppresses also the immuno-stimulation properties of bacterial or viral DNAs that contain CPGS: However, effects of methylation on the DNA structure and dynamics are not clear. Here we carried out a 10 ns MD simulation, confronted to an NMR analysis, of a hexadecanucleotide with the cAMP responsive element (CRE) DNA methylated at its center: d(GAGATGAmCGTCATCTC)(2) (CREmet). Methylation does not introduce significant structure modification but reduces the dynamics. Molecular mechanics and generalized Born solvation energy calculations showed that the stiffness of CREmet arises from both a restriction of the conformational space by the bulky methyl groups and a folding of DNA around the hydrophobic methyls. The latter effect is favored when the GpA steps belonging to the TGA binding half-sites adopt the BII conformation. The inability of the methylated DNAs to interact with their protein partners-either transcription factors for gene regulation or a Toll-like receptor for immunostimulation-could result from both the obstacle created by methyls, preventing crucial interactions, and the loss of DNA flexibility, reducing its adaptability. Results are discussed in the light of NMR and crystallographic data.

Stable formation of a pyrimidine-rich loop hairpin in a cruciform promoter

We have determined the solution structure of a TCC-loop hairpin in the cruciform promoter for the bacteriophage N4 virion RNA polymerase (N4 vRNAP). This hairpin and its complementary GGA-loop hairpin are extruded at physiological superhelical density and are required for vRNAP recognition. Contrary to its complementary GGA-loop, the three pyrimidines in the TCC-loop are all unpaired. However, with the help of two juxtaposed stem Watson-Crick G.C base-pairs, each nucleotide in the loop employs a special method to stabilize the hairpin structure. The resulting structures display extensive loop base-stacking rearrangement yet minor backbone distortion, which is largely accomplished through some loop zeta and alpha torsional angle changes. Consistent with the structural studies, UV melting of the GAAGCTCCGCTTC hairpin revealed a higher melting temperature (66 degrees C) than that of the GAACGTCCCGTTC hairpin (58 degrees C) with reversed stem G.C base-pairs, indicating significant contribution from the extra three loop-stem H-bonds. Thermodynamic parameters DeltaG degrees 25of the GAAGCTCCGCTTC hairpin and its complementary GAAGCGGAGCTTC hairpin are -4.1 and -4. 3 kcal/mol respectively, indicating approximately equal contribution of each hairpin to the cruciform formation of the N4 virion RNA polymerase promoter. No significant loop dynamics in the microsecond to millisecond NMR time-scale was observed, and the abundant well-defined exchangeable and non-exchangeable proton NOEs allowed us to efficiently determine a well-converged family for the final structures of the TCC-loop hairpin.

Unexpected BII conformer substate population in unoriented hydrated films of the d(CGCGAATTCGCG)2 dodecamer and of native B-DNA from salmon testes

Conformational substates of B-DNA had been observed so far in synthetic oligonucleotides but not in naturally occurring highly polymeric B-DNA. Our low-temperature experiments show that native B-DNA from salmon testes and the d(CGCGAATTCGCG)2 dodecamer have the same BI and BII substates. Nonequilibrium distribution of conformer population was generated by quenching hydrated unoriented films to 200 K, and isothermal structural relaxation toward equilibrium by interconversion of substates was followed by Fourier transform infrared spectroscopy. BI interconverts into BII on isothermal relaxation at 200 K, whereas on slow cooling from ambient temperature, BII interconverts into BI. Our estimation of the dodecamer's BI-to-BII conformer substate population by curve resolution of the symmetrical stretching vibration of the ionic phosphate is 2.4 +/- 0.5 to 1 at 200 K, and it is 1.3 +/- 0.5 to 1 between 270 and 290 K. Pronounced spectral changes upon BI-to-BII interconversion are consistent with base destacking coupled with migration of water from ionic phosphate toward the phosphodiester and sugar moieties. Nonspecific interaction of proteins with the DNA backbone could become specific by induced-fit-type interactions with either BI or BII backbone conformations. This suggests that the BI-to-BII substate interconversion could be a major contributor to the protein recognition process.

Stable sheared A.C pair in DNA hairpins

Single-residue d(Pu1NPu2) (Pu1.Pu2=G.A, G.G or A.A) hairpin loops can be stably closed by sheared purine.purine pairs. These special motifs have been found in several important biological systems. We now extend these loop-closing base-pairs to a sheared purine. pyrimidine (A.C) pair at a neutral pH condition. High-resolution NMR spectroscopy, distance geometry, and molecular dynamics methods were used to study d(GTACANCGTAC) oligomers. Numerous idiosyncratic nuclear Overhauser enhancements, especially those across the A.C base-pair between C4NH2left and right arrow AH1', C4NH2left and right arrow AH2, and CH5left and right arrow AH2 proton pairs, clearly define the novel sheared nature of the closing A.C base-pair. This novel base-pair is possibly present in several biological systems and in two single-stranded DNA aptamers selected from oligonucleotide libraries.

Solution-state structure of a DNA dodecamer duplex containing a Cis-syn thymine cyclobutane dimer, the major UV photoproduct of DNA

The solution structures of a duplex DNA dodecamer containing a cis-syn cyclobutane thymine dimer d(GCACGAAT[cs]TAAG).d(CTTAATTCG TGC) and its native parent sequence were determined using NMR data collected at 750 MHz. The dodecamer sequence corresponds to the section of a site-specific cis-syn dimer containing 49-mer that was found to be the binding site for the dimer-specific T4 denV endonuclease V repair enzyme by chemical and enzymatic footprinting experiments. Structures of both sequences were derived from NOE restrained molecular dynamics/simulated annealing calculations using a fixed outer layer of water and an inner dynamic layer of water with sodium counterions. The resulting structures reveal a subtle distortion to the phosphodiester backbone in the dimer-containing sequence which includes a BII phosphate at the T9pA10 junction immediately 3' to the dimer. The BII phosphate, established experimentally by analysis of the 31P chemical shifts and interpretation of the 3JP-H3' values using an optimized Karplus relationship, enables the DNA helix to accommodate the dimer by destacking the base 3' to the dimer. Furthermore, the structures provide explanations for the unusually shifted T8-N3H imino, A16-H2 and T8-Me proton resonances and T9pA10 (31)P NMR resonance and are consistent with bending, unwinding, and thermodynamic data. The implications of the structural data for the mechanism by which cis-syn dimers are recognized by repair enzymes and bypassed by DNA polymerases are also discussed.

Solution structure of a non-palindromic 16 base-pair DNA related to the HIV-1 kappa B site: evidence for BI-BII equilibrium inducing a global dynamic curvature of the duplex

1H and 31P NMR spectroscopy have been used together with molecular modelling to determine the fine structure of a non-palindromic 16 bp DNA containing the NF-kappa B binding site. Much emphasis has been placed upon NMR optimization of both two-dimensional 31P NMR techniques to extract structural information defining the phosphodiester backbone conformation and selective homonuclear 2D COSY experiments to determine sugar conformations. NMR data show evidence for a dynamic behaviour of steps flanking the ten base-pairs of the NF-kappa B binding site. A BI-BII equilibrium at these steps is demonstrated and two models for each extreme conformation are proposed in agreement with NMR data. In the refined BII structures, the NF-kappa B binding site exhibits an intrinsic curvature towards the major groove that is magnified by the four flanking steps in the BII conformation. Furthermore, the base-pairs are translated into the major groove. Thus, we present a novel mode of dynamic intrinsic curvature compatible with the DNA curvature observed in the X-ray structure of the p50-DNA complex.

Centromeric pyrimidine strands fold into an intercalated motif by forming a double hairpin with a novel T:G:G:T tetrad: solution structure of the d(TCCCGTTTCCA) dimer

The solution structures of the oligodeoxynucleotides d(CCCGTTTCC) and d(TCCCGTTTCCA) have been determined by two-dimensional NMR spectroscopy. These oligomers are part of a DNA box in human centromeric alpha satellite targeted by the centromere protein B (CENP-B). Both CENP-B and its recognition box in alphoid DNA are conserved in mammals, suggesting an important biological role. At acidic pH, d(CCCGTTTCC), d(TCCCGTTTCCA) and the full d(TCCCGTTTCCAACGAAG) CENP-B box strand all fold and dimerize in solution forming a stable bimolecular structure containing two GTTT hairpin loops that interact through a novel T : G : G : T tetrad. The stem region of the dimer is a four-stranded intercalated motif in which the hairpin monomers are parallel and held together by C : C+ hydrogen-bonding and intercalation. The loops are at the same end of the dimer and lie across the narrow grooves of the tetraplex. They are remarkably structured and stabilized by base-base cross-stacking, sugar-base stacking, and parallel G:G and antiparallel G:T pairing. In the d(TCCCGTTTCCA)2 structure, the intercalated motif is continued at the other end of the dimer with unpaired but stacked adenine and thymine bases. The possible biological implications of these structures are discussed.

Bioactive and nuclease-resistant L-DNA ligand of vasopressin

In vitro selection experiments have produced nucleic acid ligands (aptamers) that bind tightly and specifically to a great variety of target biomolecules. The utility of aptamers is often limited by their vulnerability to nucleases present in biological materials. One way to circumvent this problem is to select an aptamer that binds the enantiomer of the target, then synthesize the enantiomer of the aptamer as a nuclease-insensitive ligand of the normal target. We have so identified a mirror-image single-stranded DNA that binds the peptide hormone vasopressin and have demonstrated its stability to nucleases and its bioactivity as a vasopressin antagonist in cell culture.

The detailed structure of tandem G.A mismatched base-pair motifs in RNA duplexes is context dependent

The solution structure of the RNA duplex (rGGGCUGAAGCCCU), containing tandem G.A mismatches has been determined by NMR spectroscopy and restrained molecular dynamics. A homonuclear 3D TOCSY-NOESY was used to derive 18 to 30 distance restraints per nucleotide, as well as all gamma torsion angles and sugar puckers for the central UGAA part of the molecule. Using these constraints, together with cross-strand distances, involving exchangeable imino protons, and essentially all other torsion angles that can accurately be determined (i.e. beta, epsilon) otherwise, the structure of the UGAA domain could be determined with high precision (r.m.s.d. 0.62 A), without the aid of isotopically enriched RNA. The G.A base-pairs are of the sheared pairing type, with both nucleotides in the anti conformation, and hydrogen bonds between the guanine 2-amino and the adenine N7 and between the guanine N3 and the adenine 6-amino. Surprisingly the sugar of the guanosine of the G.A. mismatch adopts a 2'-endo sugar pucker conformation. Comparison with other RNA structures, in which two such G.A base-pairs are formed reveals that this detailed structure depends on the identity of the base 5' to the guanosine in the tandem G.A base-pairs. A geometrical model for the incorporation of sheared tandem G.A base-pairs in A-form helices is formulated, which explains the distinct different stacking properties and helical parameters in sequences containing tandem, sheared G.A base-pairs.

Homopurine and homopyrimidine strands complementary in parallel orientation form an antiparallel duplex at neutral pH with A-C, G-T, and T-C mismatched base pairs

DNA sequences d-TGAGGAAAGAAGGT (a 14-mer) and d-CTCCTTTCTTCC (a 12-mer) are complementary in parallel orientation forming either Donohue (reverse Watson-Crick) base pairing at neutral pH or Hoogsteen base pairing at slightly acidic pH. The structure of the complex formed by dissolving the two strands in equimolar ratio in water has been investigated by nmr. At neutral pH, the system forms an ordered antiparallel duplex with five A : T and four G : C Watson-Crick base pairs and three mismatches, namely G-T, A-C, and T-C. The nuclear Overhauser effect cross-peak pattern suggests an overall B-DNA conformation with major structural perturbations near the mismatches. The duplex has a low melting point and dissociates directly into single strands with a broad melting profile. The hydrogen-bonding schemes in the mismatched base pairs have been investigated. It has been shown earlier that in acidic pH, the system prefers a triple-stranded structure with two pyrimidine strands and one purine strand. One of the pyrimidine strands has protonated cytosines, forms Hoogsteen base pairing, and is aligned parallel to the purine strand; the other has nonprotonated cytosines and has base-pairing scheme similar to the one discussed in this paper. The parallel duplex is therefore less stable than either the antiparallel duplex or the triplex, in spite of its perfect complementarity.

Context dependent RNA-RNA recognition in a three-dimensional model of the 16S rRNA core

A 3-D model of the core of the 16S rRNA of Escherichia coli containing 328 residues has been built in the protein map derived from neutron scattering data with the help of all the available phylogenetic, biochemical, and cross-linking data. The three pseudoknots of the 16S-core cluster, through the arrangement of complex three-, four- and five-way junctions, around the neck and at the subunit interface. The roles in assembly, initiation or elongation of the three pseudoknots in ribosomal dynamics are emphasized. The 530-loop, localized on the periphery of the 30S particle, could be built with and without a pseudoknot independently of the state of the particle. The pseudoknot of the central domain controls the dynamics of an helix connected to the subunit interface which could trigger some mechanism during translation. The process of the model construction is compatible with a folding scenario in which the 5'-terminal pseudoknot controls the assembly of the central junction and the subsequent folding of the 3'-major domain. The modelling, together with the phylogenetic analysis and the experimental data, point to several potential RNA-RNA contacts which depend on the structural and sequence context in which they occur.

Sheared purine x purine pairing in biology

The Watson-Crick G x C and A x T base-paired DNA duplex has been the single most important milestone in modem molecular biology. However, it is possible that other types of stable DNA structures besides the double helix might exist, since only about 5% of the human chromosome is transcribed and expressed. Stable, four-stranded G-tetraplex DNA structures occur in the extensive tandem repeated sequences at the telomeres of chromosome. Formation of stable triplexes of the Py x Pu x Py or Pu x Pu x Py type have been implicated at the control regions of certain human genes. We review and discuss the various types of DNA duplex structures containing stable sheared base-pairs and compare their structural characteristics with that of B-DNA. Pu x Pu structural motifs are found in the highly conserved sequences at the replication origins of several single-stranded DNA viruses and in the peri-centromeric regions of human chromosomes, and may be involved in important biological functions, such as viral DNA replication and centromere formation.

Structural features of the DNA hairpin d(ATCCTA-GTTA-TAGGAT): formation of a G-A base pair in the loop

The three-dimensional structure of the hairpin formed by d(ATCCTA-GTTA-TAGGAT) has been determined by means of two-dimensional NMR studies, distance geometry and molecular dynamics calculations. The first and the last residues of the tetraloop of this hairpin form a sheared G-A base pair on top of the six Watson-Crick base pairs in the stem. The glycosidic torsion angles of the guanine and adenine residues in the G-A base pair reside in the anti and high- anti domain ( approximately -60 degrees ) respectively. Several dihedral angles in the loop adopt non-standard values to accommodate this base pair. The first and second residue in the loop are stacked in a more or less normal helical fashion; the fourth loop residue also stacks upon the stem, while the third residue is directed away from the loop region. The loop structure can be classified as a so-called type-I loop, in which the bases at the 5'-end of the loop stack in a continuous fashion. In this situation, loop stability is unlikely to depend heavily on the nature of the unpaired bases in the loop. Moreover, the present study indicates that the influence of the polarity of a closing A.T pair is much less significant than that of a closing C.G base pair.

The dependence of DNase I activity on the conformation of oligodeoxynucleotides

We have developed a sensitive continuous assay for nucleases using proton release. The assay has been applied to the determination of the kinetics of DNase I acting on short, defined deoxyoligonucleotides. The dependence of Kcat/K(m) on sequence and structure of short oligonucleotide substrates has been measured: increasing lengths of AnTn sequences decrease the rate of cleavage. G.A mismatches in which the bases pair using imino protons are cleaved quite effectively by DNase I. In contrast, tandem G.A mismatches which use amino pairing and have BII phosphodiesters, are refractory to DNase I. Also, the DNA strands of DNA.RNA hybrid duplexes are not cleaved by DNase I. These results show that the global conformation of a duplex and the details of its minor groove affect the cleavage efficiency by DNase I. The assay has also been used to measure the inhibition constant of the minor-groove-binding ligand propamidine. A value of 3 microM has been determined for binding to the sequence d(CGCGAATTCGCG)2, showing that dissociation constants can be determined even when there are no convenient optical signals for titrations.

Hairpin loops consisting of single adenine residues closed by sheared A.A and G.G pairs formed by the DNA triplets AAA and GAG: solution structure of the d(GTACAAAGTAC) hairpin

The DNA undecamers GTACAAAGTAC (AAA 11-mer) and GTACGAGGTAC (GAG 11-mer) have been studied in solution by high-resolution NMR spectroscopy. Both duplexes form stable hairpins containing single deoxyadenosine loops and stems containing five base-pairs that are closed at the loop end by sheared AxA and GxC pairs, respectively. These molecules thus contain new AAA and GAG loop turn motifs. All protons, including the chiral H5'/H5" protons of the loop residues, were assigned using NOESY, DQF-COSY and heteronuclear 1H-31P COSY experiments. The backbone torsion angles were constrained using experimental data from NOE crosspeaks, three-bond 1H-1H coupling constants and four-bond 1H-31P coupling constants and four-bond 1H-31P coupling constants. The AAA and GAG 11-mers form similar structures in solution. The detailed structure of the AAA 11-mer was determined by the combined use of NMR, distance geometry and energy minimization methods. This structure exhibits good stacking of the loop adenosine base on the closing 5Ax7A sheared pair, with the 6A base stacking on the 5A base and the 6A deoxyribose stacking with the 7A base. All sugars in the AAA 11-mer hairpin adopt the typical DNA C2'-endo conformation and a sharp backbone turn occurs between residues 6A and 7A. This loop turn is brought about mainly by a change in the backbone phosphate torsion angles from zeta(g-) alpha(g-) to zeta(g+) alphat(g+) at the turn. The gamma torsion angle of residue 7A in the closing sheared pair also changes from gauche+ to trans. In Pu1NPu2 loop turns of the GCA, AAA and GAG types, the chemical shift of the H4' proton of the loop deoxyribose depends on the nature of Pu2; this reflects the stacking of the loop sugar on the Pu2 base and the different ring current effects of A or G in this position.

DNA structural forms

In the years that have passed since the publication of Wolfram Saenger's classic book on nucleic acid structure (Saenger, 1984), a considerable amount of new data has been accumulated on the range of conformations which can be adopted by DNA. Many unusual species have joined the DNA zoo, including new varieties of two, three and four stranded helices. Much has been learnt about intrinsic DNA curvature, dynamics and conformational transitions and many types of damaged or deformed DNA have been investigated. In this article, we will try to summarise this progress, pointing out the scope of the various experimental techniques used to study DNA structure, and, where possible, trying to discern the rules which govern the behaviour of this subtle macromolecule. The article is divided into six major sections which begin with a general discussion of DNA structure and then present successively, B-DNA, DNA deformations, A-DNA, Z-DNA and DNARNA hybrids. An extensive set of references is included and should serve the reader who wishes to delve into greater detai.

A single G-to-C change causes human centromere TGGAA repeats to fold back into hairpins

Recently, we established that satellite III (TGGAA)n tandem repeats, which occur at the centromeres of human chromosomes, pair with themselves to form an unusual "self-complementary" antiparallel duplex containing (GGA)2 motifs in which two unpaired guanines from opposite strands intercalate between sheared G.A base pairs. In separate studies, we have also established that the GCA triplet does not form bimolecular (GCA)2 motifs but instead promotes the formation of hairpins containing a GCA-turn motif in which the loop contains a single cytidine closed by a sheared G.A pair. Since TGCAA is the most frequent variant of TGGAA found in satellite III repeats, we reasoned that the potential of this variant to form GCA-turn miniloop fold-back structures might be an important factor in modulating the local structure in natural (TGGAA)n repeats. We report here the NMR-derived solution structure of the heptadecadeoxynucleotide (G)TGGAATGCAATGGAA(C) in which a central TGCAA pentamer is flanked by two TGGAA pentamers. This 17-mer forms a rather unusual and very stable hairpin structure containing eight base pairs in the stem, only four of which are Watson-Crick pairs, and a loop consisting of a single cytidine residue. The stem contains a (GGA)2 motif with intercalative 14G/4G stacking between two sheared G.A base pairs; the loop end of the stem consists of a sheared 8G.10A closing pair with the cytosine base of the 9C loop stacked on 8G. The remarkable stability of this unusual hairpin structure (Tm = 63 degrees C) suggests that it probably plays an important role in modulating the folding of satellite III (TGGAA)n repeats at the centromere.

Nonplanar DNA base pairs

Three-dimensional structures of a representative set of more than 30 hydrogen-bonded nucleic acids pairs have been studied by reliable ab initio quantum mechanical methods. We show that many hydrogen-bonded nucleic acid base pairs are intrinsically nonplanar, mainly due to the partial sp3 hybridization of nitrogen atoms of their amino groups and secondary electrostatic interactions. This finding extends the variability of intermolecular interactions of DNA bases in that i) flexibility of the base pairs is larger than has been assumed before, and ii) attractive proton-proton acceptor interactions oriented out of the base pair plane are allowed. For example, all four G-A mismatch base pairs are propeller twisted, and the energy preferences for the nonplanar structures range from less than 0.1 kcal/mol to 1.8 kcal/mol. We predict that nonplanarity of the amino group of guanine in the G(anti)...A(anti) pair of the ApG step of the d(CCAAGATTGG)2 crystal structure is an important stabilizing factor that improves the energy of this structure by almost 3 kcal/mol. Currently used empirical potentials are not accurate enough to properly cover the interactions associated with amino-group and base-pair nonplanarity.

The thermal stability of DNA fragments with tandem mismatches at a d(CXYG).d(CY'X'G) site

Temperature-Gradient Gel Electrophoresis (TGGE) was employed to determine the thermal stabilities of 28 DNA fragments, 373 bp long, with two adjacent mismatched base pairs, and eight DNAs with Watson-Crick base pairs at the same positions. Heteroduplex DNAs containing two adjacent mismatches were formed by melting and reannealing pairs of homologous 373 bp DNA fragments differing by two adjacent base pairs. Product DNAs were separated based on their thermal stability by parallel and perpendicular TGGE. The polyacrylamide gel contained 3.36 M urea and 19.2 % formamide to lower the DNA melting temperatures. The order of stability was determined in the sequence context d(CXYG).d(CY'X'G) where X.X' and Y.Y" represent the mismatched or Watson-Crick base pairs. The identity of the mismatched bases and their stacking interactions influence DNA stability. Mobility transition melting temperatures (T u) of the DNAs with adjacent mismatches were 1.0-3.6 degrees C (+/-0.2 degree C) lower than the homoduplex DNA with the d(CCAG).d(CTGG) sequence. Two adjacent G.A pairs, d(CGAG).d(CGAG), created a more stable DNA than DNAs with Watson-Crick A.T pairs at the same sites. The d(GA).d(GA) sequence is estimated to be 0.4 (+/-30%) kcal/mol more stable in free energy than d(AA).d(TT) base pairs. This result confirms the unusual stability of the d(GA).d(GA) sequence previously observed in DNA oligomers. All other DNAs with adjacent mismatched base pairs were less stable than Watson-Crick homoduplex DNAs. Their relative stabilities followed an order expected from previous results on single mismatches. Two homoduplex DNAs with identical nearest neighbor sequences but different next-nearest neighbor sequences had a small but reproducible difference in T u value. This result indicates that sequence dependent next neighbor stacking interactions influence DNA stability.