Thermodynamics of triple helix formation: spectrophotometric studies on the d(A)10.2d(T)10 and d(C+3T4C+3).d(G3A4G3).d(C3T4C3) triple helices.

We have stabilized the d(A)10.2d(T)10 and d(C+LT4C+3).d(G3A4G3).d(C3T4C3) triple helices with either NaCl or MgCl2 at pH 5.5. UV mixing curves demonstrate a 1:2 stoichiometry of purine to pyrimidine strands under the appropriate conditions of pH and ionic strength. Circular dichroic titrations suggest a possible sequence-independent spectral signature for triplex formation. Thermal denaturation profiles indicate the initial loss of the third strand followed by dissociation of the underlying duplex with increasing temperature. Depending on the base sequence and ionic conditions, the binding affinity of the third strand for the duplex at 25 degrees C is two to five orders of magnitude lower than that of the two strands forming the duplex. Thermodynamic parameters for triplex formation were determined for both sequences in the presence of 50 mM MgCl2 and/or 2.0 M NaCl. Hoogsteen base pairs are 0.22-0.64 kcal/mole less stable than Watson-Crick base pairs, depending on ionic conditions and base composition. C+.G and T.A Hoogsteen base pairs appear to have similar stability in the presence of Mg2+ ions at low pH.     

The DNA sequence at echinomycin binding sites determines the structural changes induced by drug binding: NMR studies of echinomycin binding to [d(ACGTACGT)]2 and [d(TCGATCGA)]2

The complexes formed between the cyclic octadepsipeptide antibiotic echinomycin and the two DNA octamers [d(ACGTACGT)]2 and [d(TCGATCGA)]2 have been investigated by using one- and two-dimensional proton NMR spectroscopy techniques. The results obtained for the two complexes are compared to each other, to the crystal structures of related DNA-echinomycin complexes, and to enzymatic and chemical footprinting results. In the saturated complexes, two echinomycin molecules bind to each octamer by bisintercalation of the quinoxaline moieties on either side of each CpG step. Binding of echinomycin to the octamer [d(ACGTACGT)]2 is cooperative so that only the two-drug complex is observed at lower drug-DNA ratios, but binding to [d(TCGATCGA)]2 is not cooperative. At low temperatures, both the internal and terminal A.T base pairs adjacent to the binding site in the [d(ACGTACGT)]2-2 echinomycin complex are Hoogsteen base paired (Gilbert et al., 1989) as observed in related crystal structures. However, as the temperature is raised, the internal A.T Hoogsteen base pairs are destabilized and are observed to be exchanging between the Hoogsteen base-paired and an open (or Watson-Crick base-paired) state. In contrast, in the [d(TCGATCGA)]2-2 echinomycin complex, no A.T Hoogsteen base pairs are observed, the internal A.T base pairs appear to be stabilized by drug binding, and the structure of the complex does not change significantly from 0 to 45 degrees C. Thus, the structure and stability of the DNA in echinomycin-DNA complexes depends on the sequence at and adjacent to the binding site. While we conclude that no single structural change in the DNA can explain all of the footprinting results, unwinding of the DNA helix in the drug-DNA complexes appears to be an important factor while Hoogsteen base pair formation does not.     

1H NMR study of the binding of bis(acridines) to d(AT)5.d(AT)5. 2. Dynamic aspects

Measurements of the 1H NMR spectra and relaxation rates were used to study the dynamic properties of 9-aminoacridine (9AA) and four bis(acridine) complexes with d(AT)5.d(AT)5. The behavior of the 9AA (monointercalator) and that of C8 (bisintercalator containing an eight-carbon atom linker chain) are entirely similar. For both compounds, the lifetime of the drug in a particular binding site is 2-3 ms at approximately 20 degrees C, and neither affects the A.T base pair opening rates. The complex with C10 (bisintercalator containing a 10-carbon atom linker chain) is slightly more stable than the C8 complex since its estimated binding site lifetime is 5-10 ms at 29 degrees C. Base pairs adjacent to the bound C10 are destabilized, relative to free d(AT)5.d(AT)5, but other base pairs in the C10 complex are little affected. Bis(acridine) pyrazole (BAPY) and bis(acridine) spermine (BAS) considerably stabilize those base pairs that are sandwiched between the two acridine chromophores, but in the BAS complex proton exchange from the two flanking base pairs appears to be accelerated, relative to free d(AT)5.d(AT)5. The lifetime of these drugs in specific binding sites is too long (>10 ms) to be manifested in increased line widths, at least up to 41 degrees C. An important conclusion from this study is that certain bisintercalators rapidly migrate along DNA, despite having large binding constants (K>10(6) M-1). For C8 and C10 complexes, migration rates are little different from those deduced for 9AA. The rigid linker chain in BAPY and the charge interactions in BAS retard migration of these two bisintercalators. These results provide new parameters that are useful in understanding the biochemical and biological properties of these and other bisintercalating drugs.     

Two-dimensional NMR studies of d(GGTTAATGCGGT).d(ACCGCATTAACC) complexed with the minor groove binding drug SN-6999.

The dodecadeoxynucleotide duplex d(GGTTAATGCGGT).d(ACCGCATTAACC) and its 1:1 complex with the minor groove binding drug SN-6999 have been prepared and studied by two-dimensional 1H nuclear magnetic resonance spectroscopy. Complete sequence-specific assignments have been obtained for the free duplex by standard methods. The line widths of the resonances in the complex are greater than those observed for the free duplex, which complicates the assignment process. Extensive use of two-quantum spectroscopy was required to determine the scalar correlations for identifying all of the base proton and most of the 1'H-2'H-2'H spin subsystems for the complex. This permitted unambiguous sequence-specific resonance assignments for the complex, which provides the necessary background for a detailed comparison of the structure of the duplex, with and without bound drug. A series of intermolecular NOEs between drug and DNA were identified, providing sufficient structural constraints to position the drug in the minor groove of the duplex. However, the combination of NOEs observed can only be rationalized by a model wherein the drug binds in the minor groove of the DNA in both orientations relative to the long helix axis and exchanges rapidly between the two orientations. The drug binds primarily in the segment of five consecutive dA-dT base pairs d(T3T4A5A6T7).d(A18T19T20A21A22), but surprisingly strong interactions are found to extend one residue in the 3' direction along each strand to G8 and C23. The observation of intermolecular contacts to residues neighboring the AT-rich region demonstrates that the stabilization of the bis(quaternary ammonium) heterocycle family of AT-specific, minor groove binding drugs is not based exclusively on interactions with dA-dT base pairs.     

Sequence-dependent recognition of DNA duplexes. Netropsin complexation to the AATT site of the d(G-G-A-A-T-T-C-C) duplex in aqueous solution

We have investigated intermolecular interactions and conformational features of the netropsin X d(G-G-A-A-T-T-C-C) complex by one- and two-dimensional NMR studies in aqueous solution. Netropsin removes the 2-fold symmetry of the d(G-G-A-A-T-T-C-C) duplex at the AATT binding site and to a lesser extent at adjacent dG X dC base pairs resulting in doubling of resonances for specific positions in the spectrum of the complex at 25 degrees C. We have assigned the amide, pyrrole, and CH2 protons of netropsin, and the base and sugar H1' protons of the nucleic acid from an analysis of the nuclear Overhauser effect (NOESY) and correlated (COSY) spectra of the complex at 25 degrees C. We observe intermolecular nuclear Overhauser effects (NOE) between all three amide and both pyrrole protons on the concave face of the antibiotic and the minor groove adenosine H2 proton of the two central A4 X T5 base pairs of the d(G1-G2-A3-A4-T5-T6-C7-C8) duplex. Weaker intermolecular NOEs are also observed between the pyrrole concave face protons and the sugar H1' protons of residues T5 and T6 in the AATT minor groove of the duplex. We also detect intermolecular NOEs between the guanidino CH2 protons at one end of netropsin and adenosine H2 proton of the two flanking A3 X T6 base pairs of the octanucleotide duplex. These studies establish a set of intermolecular contacts between the concave face of the antibiotic and the minor groove AATT segment of the d(G-G-A-A-T-T-C-C) duplex in solution. The magnitude of the NOEs require that there be no intervening water molecules sandwiched between the antibiotic and the DNA so that release of the minor groove spine of hydration is a prerequisite for netropsin complex formation.     

Interaction between the Z-type DNA duplex and 1,3-propanediamine: crystal structure of d(CACGTG)2 at 1.2 A resolution

The crystal structure of a hexamer duplex d(CACGTG)(2) has been determined and refined to an R-factor of 18.3% using X-ray data up to 1.2 A resolution. The sequence crystallizes as a left-handed Z-form double helix with Watson-Crick base pairing. There is one hexamer duplex, a spermine molecule, 71 water molecules, and an unexpected diamine (Z-5, 1,3-propanediamine, C(3)H(10)N(2)) in the asymmetric unit. This is the high-resolution non-disordered structure of a Z-DNA hexamer containing two AT base pairs in the interior of a duplex with no modifications such as bromination or methylation on cytosine bases. This structure does not possess multivalent cations such as cobalt hexaammine that are known to stabilize Z-DNA. The overall duplex structure and its crystal interactions are similar to those of the pure-spermine form of the d(CGCGCG)(2) structure. The spine of hydration in the minor groove is intact except in the vicinity of the T5A8 base pair. The binding of the Z-5 molecule in the minor grove of the d(CACGTG)(2) duplex appears to have a profound effect in conferring stability to a Z-DNA conformation via electrostatic complementarity and hydrogen bonding interactions. The successive base stacking geometry in d(CACGTG)(2) is similar to the corresponding steps in d(CG)(3). These results suggest that specific polyamines such as Z-5 could serve as powerful inducers of Z-type conformation in unmodified DNA sequences with AT base pairs. This structure provides a molecular basis for stabilizing AT base pairs incorporated into an alternating d(CG) sequence.     

2D NMR study of the DNA duplex d(CTCTC*A*ACTTCC).d(GGAAGTTGAGAG) cross-linked by the antitumor-active dirhodium(II,II) unit at the cytosine-adenine step

The 2D NMR analysis in solution of the DNA duplex d(CTCTC*A*ACTTCC).d(GGAAGTTGAGAG) binding to the dirhodium unit cis-[Rh2(mu-O2CCH3)2(eta1-O2CCH3)]+ showed that an unprecedented intrastrand adduct, dsII, is formed with the dirhodium unit cross-linking in the major groove residues C5 and A6 (indicated with asterisks), also corroborated by enzyme digestion studies. Formation of the dirhodium complex dsII destabilizes significantly the duplex as indicated by the substantial decrease in its melting temperature (DeltaTm = -22.9 degrees C). The reduced thermal stability of dsII is attributed to the decreased stacking of the bases and the complete disruption and/or weakening of the hydrogen bonds within the base pairs in the immediate vicinity of the metalation site (C5.G20 and A6.T19), but the effects due to the metal binding are more severe for the base pairs in the 5' direction to the lesion site. The NMR spectroscopic data indicate that Watson-Crick hydrogen bonding is completely disrupted for the C5.G20 site and considerably weakened for A6.T19. In dsII, the bases C5 and A6 bind to eq positions of the dirhodium unit cis-[Rh2(mu-O2CCH3)2(eta1-O2CCH3)]+, which retains one monodentate and two bridging acetate groups, presumably due to steric reasons. Binding of A6 takes place via N7, whereas binding of the C5 base takes place via the exocyclic N4 site, resulting in the anti-cytosine rotamer with respect to site N3 in its metal-stabilized rare iminooxo form.     

Crystal structure of a B-DNA dodecamer containing inosine, d(CGCIAATTCGCG), at 2.4 A resolution and its comparison with other B-DNA dodecamers.

The crystal structure of the dodecamer, d(CGCIAATTCGCG), has been determined at 2.4 A resolution by molecular replacement, and refined to an R-factor of 0.174. The structure is isomorphous with that of the B-DNA dodecamer, d(CGCGAATTCGCG), in space group P2(1)2(1)2(1) with cell dimensions of a = 24.9, b = 40.4, and c = 66.4 A. The initial difference Fourier maps clearly indicated the presence of inosine instead of guanine. The structure was refined with 44 water molecules, and compared to the parent dodecamer. Overall the two structures are very similar, and the I:C forms Watson-Crick base pairs with similar hydrogen bond geometry to the G:C base pairs. The propeller twist angle is low for I4:C21 and relatively high for the I16:C9 base pair (-3.2 degrees compared to -23.0 degrees), and the buckle angles alter, probably due to differences in the contacts with symmetry related molecules in the crystal lattice. The central base pairs of d(CGCIAATTCGCG) show the large propeller twist angles, and the narrow minor groove that characterize A-tract DNA, although I:C base pairs cannot form the major groove bifurcated hydrogen bonds that are possible for A:T base pairs.     

Proton NMR study of the [d(ACGTATACGT)]2-2echinomycin complex: conformational changes between echinomycin binding sites.

The interactions of echinomycin and the DNA decamer [d(ACGTATACGT)]2 were studied by proton NMR. Echinomycin binds cooperatively as a bisintercalator at the CpG steps. The terminal A.T base pairs are Hoogsteen base paired, but none of the four central A.T base pairs are Hoogsteen base paired. However, binding of the drug induces unwinding of the DNA which is propagated to the central ApT step. All four central A.T base pairs are destabilized relative to those in the free DNA. Furthermore, based on these and other results from our laboratory, we conclude that the formation of stable Hoogsteen base pairs may not be the relevant structural change in vivo. The structural changes propagated between adjacent ACGT binding sites are the unwinding of the duplex and destabilization of the base pairing between binding sites.     

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

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

Binding of actinomycin D to single-stranded DNA of sequence motifs d(TGTCT(n)G) and d(TGT(n)GTCT)

Our recent binding studies with oligomers derived from base replacements on d(CGTCGTCG) had led to the finding that actinomycin D (ACTD) binds strongly to d(TGTCATTG) of apparent single-stranded conformation without GpC sequence. A fold-back binding model was speculated in which the planar phenoxazone inserts at the GTC site with a loop-out T base whereas the G base at the 3'-terminus folds back to form a basepair with the internal C and stacks on the opposite face of the chromophore. To provide a more concrete support for such a model, ACTD equilibrium binding studies were carried out and the results are reported herein on oligomers of sequence motifs d(TGTCT(n)G) and d(TGT(n)GTC). These oligomers are not expected to form dimeric duplexes and contain no canonical GpC sequences. It was found that ACTD binds strongly to d(TGTCTTTTG), d(TGTTTTGTC), and d(TGTTTTTGTC), all exhibiting 1:1 drug/strand binding stoichiometry. The fold-back binding model with displaced T base is further supported by the finding that appending TC and TCA at the 3'-terminus of d(TGTCTTTTG) results in oligomers that exhibit enhanced ACTD affinities, consequence of the added basepairing to facilitate the hairpin formation of d(TGTCTTTTGTC) and d(TGTCTTTTGTCA) in stabilizing the GTC/GTC binding site for juxtaposing the two G bases for easy stacking on both faces of the phenoxazone chromophore. Further support comes from the observation of considerable reduction in ACTD affinity when GTC is replaced by GTTC in an oligomer, in line with the reasoning that displacing two T bases to form a bulge for ACTD binding is more difficult than displacing a single base. Based on the elucidated binding principle of phenoxazone ring requiring its opposite faces to be stacked by the 3'-sides of two G bases for tight ACTD binding, several oligonucleotide sequences have been designed and found to bind well.     

NMR studies of pH-dependent conformational polymorphism of alternating (C-T)n sequences.

Alternating (C-T)n sequences are involved in the H-DNA structure associated with (GA)n.(CT)n sequences. Low pH values facilitate H-DNA formation. We have undertaken a detailed analysis of the structural consequences of the (C-T)n sequence as a function of pH. The structures of three DNA oligonucleotides, d(CT)4, d(TC)4 and d(TC)15, have been studied by NMR. We found that their conformations are polymorphic and pH dependent. There are at least three major conformational species: an antiparallel-stranded (APS) duplex with entirely C:T base pairs at pH 7, an antiparallel-stranded (APS) duplex with entirely C+:T base pairs at pH 3, and a possible parallel-stranded (PS) duplex with C+:C and T:T base pairs near pH 5. In the intermediate pH range, the APS duplex may have varying numbers of C+:T and C:T base pairs, and there may be a fast exchange going on between APS duplex species involving these two kinds of base pairs. However, the transition between the APS and PS duplexes is slow. Structural refinement of the two octamers, d(TC)4 and d(CT)4, at pH = 6.9 and pH = 3 using 2D-NOE data suggests that the molecules are likely in the duplex form at 5 degrees C. We lack evidence that the structure at pH 3 is a PS structure with T nucleotides residing in the exterior of the helix. Titration of the longer oligonucleotide, d(TC)15, showed a prominent pKa of approximately 6, approaching the value of 7.0 obtained from the titration of poly-(dC).     

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.     

Modelling basic features of specificity in the binding of a dicationic steroid diamine to double-stranded oligonucleotides.

An investigation of the intrinsically preferred binding modes of a steroid diamine, dipyrandium, to the double-stranded hexanucleotides d(TATATA)2, d(ATATAT)2, and d(CGCGCG)2 is carried out by the energy minimization procedure JUMNA. Several alternative binding modes are compared: groove binding in which the conformation of the oligonucleotide remains close to that of B-DNA, intercalation between base-pairs and interaction with variously kinked structures in which base pairs of dinucleoside steps open towards the groove in which the binding occurs. The favored binding configuration occurs at the d(TpA) step of the AT kinked nucleotides in which the kink opens the base pairs towards the minor groove. Thus, for the d(T1A2T3A4T5A6)2 sequences the preferred complexation involves the kink at the T3A4 step facing the cyclohexane rings A, B, and C of the ligand. For the d(A1T2A3T4A5T6)2 sequence, the kink occurs at the T2A3 step facing the cationic pyrrolidine ring linked to ring A. The binding of dipyrandium to d(CGCGCG)2 is found to be considerably less favourable than for either of the two (AT) sequences.     

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

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

Crystal structure at 1.5-A resolution of d(CGCICICG), an octanucleotide containing inosine, and its comparison with d(CGCG) and d(CGCGCG) structures.

The octadeoxyribonucleotide d(CGCICICG) has been crystallized in space group P(6)5(22) with unit cell dimensions of a = b = 31.0 A and c = 43.7 A, and X-ray diffraction data have been collected to 1.5-A resolution. Precession photographs and the self-Patterson function indicate that 12 base pairs of Z-conformation DNA stack along the c-axis, and the double helices pack in a hexagonal array similar to that seen in other crystals of Z-DNA. The structure has been solved by both Patterson deconvolution and molecular replacement methods and refined in space group P(6)5 to an R factor of 0.225 using 2503 unique reflections greater than 3.0 sigma (F). Comparison of the molecules within the hexagonal lattice with highly refined crystal structures of other Z-DNA reveals only minor conformational differences, most notably in the pucker of the deoxyribose of the purine residues. The DNA has multiple occupancy of C:I and C:G base pairs, and C:I base pairs adopt a conformation similar to that of C:G base pairs.     

Kinetics for exchange of imino protons in the d(C-G-C-G-A-A-T-T-C-G-C-G) double helix and in two similar helices that contain a G . T base pair, d(C-G-T-G-A-A-T-T-C-G-C-G), and an extra adenine, d(C-G-C-A-G-A-A-T-T-C-G-C-G)

The relaxation lifetimes of imino protons from individual base pairs were measured in (I) a perfect helix, d(C-G-C-G-A-A-T-T-C-G-C-G), (II) this helix with a G . C base pair replaced with a G . T base pair, d(C-G-T-G-A-A-T-T-C-G-C-G), and (III) the perfect helix with an extra adenine base in a mismatch, d(C-G-C-A-G-A-A-T-T-C-G-C-G). The lifetimes were measured by saturation recovery proton nuclear magnetic resonance experiments performed on the imino protons of these duplexes. The measured lifetimes of the imino protons were shown to correspond to chemical exchange lifetimes at higher temperatures and spin-lattice relaxation times at lower temperatures. Comparison of the lifetimes in these duplexes showed that the destabilizing effect of the G . T base pair in II affected the opening rate of only the nearest-neighbor base pairs. For helix III, the extra adenine affected the opening rates of all the base pairs in the helix and thus was a larger perturbation for opening of the base pairs than the G . T base pair. The temperature dependence of the exchange rates of the imino proton in the perfect helix gives values of 14-15 kcal/mol for activation energies of A . T imino protons. These relaxation rates were shown to correspond to exchange involving individual base pair opening in this helix, which means that one base-paired imino proton can exchange independent of the others. For the other two helices that contain perturbations, much larger activation energies for exchange of the imino protons were found, indicating that a cooperative transition involving exchange of at least several base pairs was the exchange mechanism of the imino protons. The effects of a perturbation in a helix on the exchange rates and the mechanisms for exchange of imino protons from oligonucleotide helices are discussed.     

Selection for thermodynamically stable DNA tetraloops using temperature gradient gel electrophoresis reveals four motifs: d(cGNNAg), d(cGNABg),d(cCNNGg), and d(gCNNGc)

Hairpins play important roles in the function of DNA, forming cruciforms and affecting processes such as replication and recombination. Temperature gradient gel electrophoresis (TGGE) and in vitro selection have been used to isolate thermodynamically stable DNA hairpins from a six-nucleotide random library. The TGGE-selection process was optimized such that known stable DNA tetraloops were recovered, and the selection appears to be exhaustive. In the selection, four families of exceptionally stable DNA loops were identified: d(cGNNAg), d(cGNABg), d(cCNNGg), and d(gCNNGc). (Lowercase denotes the closing base pair; N = A, C, G, or T; and B = C, G, or T.) It appears that the known stable d(cGNAg) triloop motif can be embedded into a tetraloop, with the extra nucleotide inserted into either the middle of the loop, d(cGNNAg), or at the 3'-end of the loop, d(cGNABg). For d(cGNNAg) and d(cGNABg), a CG closing base pair was strongly preferred over a GC, with DeltaDeltaG degrees (37) approximately 2 kcal/mol. Members of the two families, d(cCNNGg) and d(gCNNGc), are similar in stability. The loop sequences and closing base pairs identified for exceptionally stable DNA tetraloops show many similarities to those known for exceptionally stable RNA tetraloops. These data provide an expanded set of thermodynamic rules for the formation of tetraloops in DNA.