Sheared Aanti.Aanti base pairs in a destabilizing 2 x 2 internal loop: the NMR structure of 5'(rGGCAAGCCU)2

The 5'(rGGCAAGCCU)(2) duplex contains tandem A.A pairs. The three-dimensional structure of the 5'(rGGCAAGCCU)(2) duplex was modeled by molecular dynamics and energy minimization with NMR-derived distance and dihedral angle restraints. Although the 5'(rCAAG)(2) loop is thermodynamically destabilizing by 1.1 kcal/mol, the tandem A.A pairs adopt a predominant conformation: a sheared anti-anti (A.A trans Hoogsteen/Sugar-edge) alignment similar to that observed in the crystal structure of the P4-P6 domain of the Tetrahymena thermophila intron [Cate, J. H., Gooding, A. R., Podell, E., Zhou, K., Golden, B. L., Kundrot, C. E., Cech, T. R., and Doudna, J. A. (1996) Science 273, 1678-1685]. The NMR-derived structure of the 5'(rGGCAAGCCU)(2) duplex exhibits cross-strand hydrogen bonds from N3 of A4 to an amino hydrogen of A5 and from the 2' oxygen of the A4 sugar to the other amino hydrogen of A5. An intrastrand hydrogen bond is formed from the 2' OH hydrogen of A4 to O5' of A5. The cross-strand A5 bases are stacked. The Watson-Crick G-C regions are essentially A-form. The sheared anti-anti (A.A trans Hoogsteen/Sugar-edge) alignment provides potential contact sites for tertiary interactions and, therefore, is a possible target site for therapeutics. Thus, thermodynamically destabilizing internal loops can be preorganized for tertiary interactions or ligand binding.     

Solution structure of an RNA internal loop with three consecutive sheared GA pairs

Internal loops in RNA are important for folding and function. Many folding motifs are internal loops containing GA base pairs, which are usually thermodynamically stabilizing, i.e., contribute favorable free energy to folding. Understanding the sequence dependence of folding stability and structure in terms of molecular interactions, such as hydrogen bonding and base stacking, will provide a foundation for predicting stability and structure. Here, we report the NMR structure of the oligonucleotide duplex, 5'GGUGGAGGCU3'/3'PCCGAAGCCG5' (P = purine), containing an unusually stable and relatively abundant internal loop, 5'GGA3'/3'AAG5'. This loop contains three consecutive sheared GA pairs (trans Hoogsteen/Sugar edge AG) with separate stacks of three G's and three A's in a row. The thermodynamic consequences of various nucleotide substitutions are also reported. Significant destabilization of approximately 2 kcal/mol at 37 degrees C is found for substitution of the middle GA with AA to form 5'GAA3'/3'AAG5'. This destabilization correlates with a unique base stacking and hydrogen-bonding network within the 5'GGA3'/3'AAG5' loop. Interestingly, the motifs, 5'UG3'/3'GA5' and 5'UG3'/3'AA5', have stability similar to 5'CG3'/3'GA5' even though UG and UA pairs are usually less stable than CG pairs. Consecutive sheared GA pairs in the 5'GGA3'/3'AAG5' loop are preorganized for potential tertiary interactions and ligand binding.     

Proton exchange and base pair opening in a DNA triple helix

Nuclear magnetic resonance spectroscopy has been used to characterize opening reactions and stabilities of individual base pairs in two related DNA structures. The first is the triplex structure formed by the DNA 31-mer 5'-AGAGAGAACCCCTTCTCTCTTTTTCTCTCTT-3'. The structure belongs to the YRY (or parallel) family of triple helices. The second structure is the hairpin double helix formed by the DNA 20-mer 5'-AGAGAGAACCCCTTCTCTCT-3' and corresponds to the duplex part of the YRY triplex. The rates of exchange of imino protons with solvent in the two structures have been measured by magnetization transfer from water and by real-time exchange at 10 degrees C in 100 mM NaCl and 5 mM MgCl2 at pH 5.5 and in the presence of two exchange catalysts. The results indicate that the exchange of imino protons in protonated cytosines is most likely limited by the opening of Hoogsteen C+G base pairs. The base pair opening parameters estimated from imino proton exchange rates suggest that the stability of individual Hoogsteen base pairs in the DNA triplex is comparable to that of Watson-Crick base pairs in double-helical DNA. In the triplex structure, the exchange rates of imino protons in Watson-Crick base pairs are up to 5000-fold lower than those in double-helical DNA. This result suggests that formation of the triplex structure enhances the stability of Watson-Crick base pairs by up to 5 kcal/mol. This stabilization depends on the specific location of each triad in the triplex structure.     

The first example of a Hoogsteen base-paired DNA duplex in dynamic equilibrium with a Watson-Crick base-paired duplex--a structural (NMR), kinetic and thermodynamic study.

A single-point substitution of the O4' oxygen by a CH2 group at the sugar residue of A6 (i.e. 2'-deoxyaristeromycin moiety) in a self-complementary DNA duplex, 5'-d(C1G2C3G4A5A6T7T8C9G10C11G12)2(-3), has been shown to steer the fully Watson-Crick basepaired DNA duplex (1A), akin to the native counterpart, to a doubly A6:T7 Hoogsteen basepaired (1B) B-type DNA duplex, resulting in a dynamic equilibrium of (1A)<==>(1B): Keq = k1/k(-1) = 0.56+/-0.08. The dynamic conversion of the fully Watson-Crick basepaired (1A) to the partly Hoogsteen basepaired (1B) structure is marginally kinetically and thermodynamically disfavoured [k1 (298K) = 3.9 0.8 sec(-1); deltaHdegrees++ = 164+/-14 kJ/mol; -TdeltaS degrees++ (298K) = -92 kJ/mol giving a deltaG degrees++ 298 of 72 kJ/mol. Ea (k1) = 167 14 kJ/mol] compared to the reverse conversion of the Hoogsteen (1B) to the Watson-Crick (1A) structure [k-1 (298K) = 7.0 0.6 sec-1, deltaH degrees++ = 153 13 kJ/mol; -TdeltaSdegrees++ (298K) = -82 kJ/mol giving a deltaGdegrees++(298) of 71 kJ/mol. Ea (k-1) = 155 13 kJ/mol]. Acomparison of deltaGdegrees++(298) of the forward (k1) and backward (k-1) conversions, (1A)<==>(1B), shows that there is ca 1 kJ/mol preference for the Watson-Crick (1A) over the double Hoogsteen basepaired (1B) DNA duplex, thus giving an equilibrium ratio of almost 2:1 in favour of the fully Watson-Crick basepaired duplex. The chemical environments of the two interconverting DNA duplexes are very different as evident from their widely separated sets of chemical shifts connected by temperature-dependent exchange peaks in the NOESY and ROESY spectra. The fully Watson-Crick basepaired structure (1A) is based on a total of 127 intra, 97 inter and 17 cross-strand distance constraints per strand, whereas the double A6:T7 Hoogsteen basepaired (1B) structure is based on 114 intra, 92 inter and 15 cross-strand distance constraints, giving an average of 22 and 20 NOE distance constraints per residue and strand, respectively. In addition, 55 NMR-derived backbone dihedral constraints per strand were used for both structures. The main effect of the Hoogsteen basepairs in (1B) on the overall structure is a narrowing of the minor groove and a corresponding widening of the major groove. The Hoogsteen basepairing at the central A6:T7 basepairs in (1B) has enforced a syn conformation on the glycosyl torsion of the 2'-deoxyaristeromycin moiety, A6, as a result of substitution of the endocyclic 4'-oxygen in the natural sugar with a methylene group in A6. A comparison of the Watson-Crick basepaired duplex (1A) to the Hoogsteen basepaired duplex (1B) shows that only a few changes, mainly in alpha, sigma and gamma torsions, in the sugar-phosphate backbone seem to be necessary to accommodate the Hoogsteen basepair.     

An alternating sheared AA pair and elements of stability for a single sheared purine-purine pair flanked by sheared GA pairs in RNA

A previous NMR structure of the duplex 5'GGU GGA GGCU/PCCG AAG CCG5' revealed an unusually stable RNA internal loop with three consecutive sheared GA pairs. Here, we report NMR studies of two duplexes, 5'GGU GGA GGCU/PCCA AAG CCG5' (replacing the UG pair with a UA closing pair) and 5'GGU GAA GGCU/PCCG AAG CCG5' (replacing the middle GA pair with an AA pair). An unusually stable loop with three consecutive sheared GA pairs forms in the duplex 5'GGU GGA GGCU/PCCA AAG CCG5'. The structure contrasts with that reported for this loop in the crystal structure of the large ribosomal subunit of Deinococcus radiodurans [Harms, J., Schluenzen, F., Zarivach, R., Bashan, A., Gat, S., Agmon, I., Bartels, H., Franceschi, F., and Yonath, A. (2001) Cell 107, 679-688]. The middle AA pair in the duplex 5'GGU GAA GGCU/PCCG AAG CCG5' rapidly exchanges orientations, resulting in alternative base stacking and pseudosymmetry with exclusively sheared pairs. The U GAA G/G AAG C internal loop is 2.1 kcal/mol less stable than the U GGA G/G AAG C internal loop at 37 degrees C. Structural, energetic, and dynamic consequences upon functional group substitutions within related 3 x 3 and 3 x 6 internal loops are also reported.     

Nearest neighbor parameters for inosine x uridine pairs in RNA duplexes

An enzyme family known as adenosine deaminases that act on RNA (ADARs) catalyzes adenosine deamination in RNA. ADARs act on RNA that is largely double-stranded and convert adenosine to inosine, resulting, in many cases, in an I x U pair. Thermodynamic parameters derived from optical melting studies are reported for a series of 14 oligoribonucleotides containing single I x U pairs adjacent to Watson-Crick pairs. In order to determine unique linearly independent nearest neighbor parameters for I x U pairs, four duplexes containing 3'-terminal I x U pairs and four duplexes containing 5'-terminal I x U pairs have also been thermodynamically characterized. This data was combined with previously published data of seven duplexes containing internal, terminal, or tandem I x U pairs from Strobel et al. [Strobel, S. A., Cech, T. R., Usman, N., and Beigelman, L. (1994) Biochemistry 33, 13824-13838] and Serra et al. [Serra, M. J., Smolter, P. E., and Westhof, E. (2004) Nucleic Acids Res. 32, 1824-1828]. On average, a duplex with an internal I x U pair is 2.3 kcal/mol less stable than the same duplex with an A-U pair, however, a duplex with a terminal I x U pair is 0.8 kcal/mol more stable than the same duplex with an A-U pair. Although isosteric with a G-U pair, on average, a duplex with an internal I x U pair is 1.9 kcal/mol less stable than the same duplex with a G-U pair, however, a duplex with a terminal I x U pair is 0.9 kcal/mol more stable than the same duplex with a G-U pair. Duplexes with tandem I x U pairs are on average 5.9 and 3.8 kcal/mol less stable than the same duplex with tandem A-U or tandem G-U pairs, respectively. Using the combined thermodynamic data and a complete linear least-squares fitting routine, nearest neighbor parameters for all nearest neighbor combinations of I x U pairs and an additional parameter for terminal I x U pairs have been derived.     

Non-nearest-neighbor dependence of the stability for RNA bulge loops based on the complete set of group I single-nucleotide bulge loops

Fifty-nine RNA duplexes containing single-nucleotide bulge loops were optically melted in 1 M NaCl, and the thermodynamic parameters DeltaH degrees, DeltaS degrees, DeltaG 37 degrees, and TM for each sequence were determined. Sequences from this study were combined with sequences from previous studies [Longfellow, C. E., et al. (1990) Biochemistry 29, 278-285; Znosko, B. M., et al. (2002) Biochemistry 41, 10406-10417], thus examining all possible group I single-nucleotide bulge loop and nearest-neighbor sequence combinations. The free energy increments at 37 degrees C for the introduction of a group I single-nucleotide bulge loop range between 1.3 and 5.2 kcal/mol. The combined data were used to develop a model for predicting the free energy of a RNA duplex containing a single-nucleotide bulge. For bulge loops with adjacent Watson-Crick base pairs, neither the identity of the bulge nor the nearest-neighbor base pairs had an effect on the influence of the bulge loop on duplex stability. The proposed model for prediction of the stability of a duplex containing a bulged nucleotide was primarily affected by non-nearest-neighbor interactions. The destabilization of the duplex by the bulge was related to the stability of the stems adjacent to the bulge. Specifically, there was a direct correlation between the destabilization of the duplex and the stability of the less stable duplex stem. The stability of a duplex containing a bulged nucleotide adjacent to a wobble base pair also was primarily affected by non-nearest-neighbor interactions. Again, there was a direct correlation between the destabilization of the duplex and the stability of the less stable duplex stem. However, when one or both of the bulge nearest neighbors was a wobble base pair, the free energy increment for insertion of a bulge loop is dependent upon the position and orientation of the wobble base pair relative the bulged nucleotide. Bulge sequences of the type ((5'UBX)(3'GY)), ((5'GBG)(3'UU)) and ((5'UBU)(3'GG)) are less destabilizing by 0.6 kcal/mol, and bulge sequences of the type ((5'GBX)(3'UY)) and ((5'XBU)(3'YG)) are more destabilizing by 0.4 kcal/mol than bulge loops adjacent to Watson-Crick base pairs.     

The NMR structure of an internal loop from 23S ribosomal RNA differs from its structure in crystals of 50s ribosomal subunits

Internal loops play an important role in structure and folding of RNA and in recognition of RNA by other molecules such as proteins and ligands. An understanding of internal loops with propensities to form a particular structure will help predict RNA structure, recognition, and function. The structures of internal loops 5' 1009CUAAG1013 3'/3' 1168GAAGC1164 5' and 5' 998CUAAG1002 3'/3' 1157GAAGC1153 5' from helix 40 of the large subunit rRNA in Deinococcus radiodurans and Escherichia coli, respectively, are phylogenetically conserved, suggesting functional relevance. The energetics and NMR solution structure of the loop were determined in the duplex 5' 1GGCUAAGAC9 3'/3' 18CCGAAGCUG10 5'. The internal loop forms a different structure in solution and in the crystal structures of the ribosomal subunits. In particular, the crystal structures have a bulged out adenine at the equivalent of position A15 and a reverse Hoogsteen UA pair (trans Watson-Crick/Hoogsteen UA) at the equivalent of U4 and A14, whereas the solution structure has a single hydrogen bond UA pair (cis Watson-Crick/sugar edge A15U4) between U4 and A15 and a sheared AA pair (trans Hoogsteen/sugar edge A14A5) between A5 and A14. There is cross-strand stacking between A6 and A14 (A6/A14/A15 stacking pattern) in the NMR structure. All three structures have a sheared GA pair (trans Hoogsteen/sugar edge A6G13) at the equivalent of A6 and G13. The internal loop has contacts with ribosomal protein L20 and other parts of the RNA in the crystal structures. These contacts presumably provide the free energy to rearrange the base pairing in the loop. Evidently, molecular recognition of this internal loop involves induced fit binding, which could confer several advantages. The predicted thermodynamic stability of the loop agrees with the experimental value, even though the thermodynamic model assumes a Watson-Crick UA pair.     

Base-pair dynamics in an antiparallel DNA triplex measured by catalyzed imino proton exchange monitored via 1H NMR spectroscopy

Using (1)H NMR spectroscopy, the base-pair opening dynamics of an antiparallel foldback DNA triplex and the corresponding duplex has been characterized via catalyzed imino proton exchange. The triplex system was found to be in an equilibrium between a duplex and a triplex form. The exchange rate between the two forms (i.e., the on/off-rate of the third strand) was measured to be 5 s(-1) at 1 degrees C, and the base-pair dynamics of both forms were investigated separately. Both Watson-Crick and reverse Hoogsteen base pairs were found to have base-pair lifetimes in the order of milliseconds. The stability of the Watson-Crick base pairs was, however, substantially increased in the presence of the third strand. In the DNA triplex, the opening dynamics of the reverse Hoogsteen base pairs was significantly faster than the dynamics of the Watson-Crick pairs. We were able to conclude that, for both Watson-Crick and reverse Hoogsteen base pairs, spontaneous and individual opening from within the closed base triplet is the dominating opening pathway.     

DNA interstrand cross-links of the novel antitumor trinuclear platinum complex BBR3464. Conformation,recognition by high mobility group domain proteins,and nucleotide excision repair

The novel phase II antitumor polynuclear platinum drug BBR3464 ([(trans-PtCl(NH(3))(2))(2)(mu-trans-Pt(NH(3))(2)(NH(2)(CH(2))(6)NH(2))(2))](NO(3))(4)) forms intra- and interstrand cross-links (CLs) on DNA (which is the pharmacological target of platinum drugs). We examined first in our recent work how various intrastrand CLs of BBR3464 affect the conformation of DNA and its recognition by cellular components (Zehnulova, J., Kasparkova, J., Farrell, N., and Brabec, V. (2001) J. Biol. Chem. 276, 22191-22199). In the present work, we have extended the studies on the DNA interstrand CLs of this drug. The results have revealed that the interstrand CLs are preferentially formed between guanine residues separated by 2 base pairs in both the 3' --> 3' and 5' --> 5' directions. The major 1,4-interstrand CLs distort DNA, inducing a directional bending of the helix axis and local unwinding of the duplex. Although such distortions represent a potential structural motif for recognition by high mobility group proteins, these proteins do not recognize 1,4-interstrand CLs of BBR3464. On the other hand, in contrast to intrastrand adducts of BBR3464, 1,4-interstrand CLs are not removed from DNA by nucleotide excision repair. It has been suggested that interstrand CLs of BBR3464 could persist considerably longer in cells compared with intrastrand adducts, which would potentiate the toxicity of the interstrand lesions to tumors sensitive to this polynuclear drug.     

Nuclear magnetic resonance spectroscopy and molecular modeling reveal that different hydrogen bonding patterns are possible for G.U pairs: one hydrogen bond for each G.U pair in r(GGCGUGCC)(2) and two for each G.U pair in r(GAGUGCUC)(2)

G.U pairs occur frequently and have many important biological functions. The stability of symmetric tandem G.U motifs depends both on the adjacent Watson-Crick base pairs, e.g., 5'G > 5'C, and the sequence of the G.U pairs, i.e., 5'-UG-3' > 5'-GU-3', where an underline represents a nucleotide in a G.U pair [Wu, M., McDowell, J. A., and Turner, D. H. (1995) Biochemistry 34, 3204-3211]. In particular, at 37 degrees C, the motif 5'-CGUG-3' is less stable by approximately 3 kcal/mol compared with other symmetric tandem G.U motifs with G-C as adjacent pairs: 5'-GGUC-3', 5'-GUGC-3', and 5'-CUGG-3'. The solution structures of r(GAGUGCUC)(2) and r(GGCGUGCC)(2) duplexes have been determined by NMR and restrained simulated annealing. The global geometry of both duplexes is close to A-form, with some distortions localized in the tandem G.U pair region. The striking discovery is that in r(GGCGUGCC)(2) each G.U pair apparently has only one hydrogen bond instead of the two expected for a canonical wobble pair. In the one-hydrogen-bond model, the distance between GO6 and UH3 is too far to form a hydrogen bond. In addition, the temperature dependence of the imino proton resonances is also consistent with the different number of hydrogen bonds in the G.U pair. To test the NMR models, U or G in various G.U pairs were individually replaced by N3-methyluridine or isoguanosine, respectively, thus eliminating the possibility of hydrogen bonding between GO6 and UH3. The results of thermal melting studies on duplexes with these substitutions support the NMR models.     

Theoretical study of the Hoogsteen-Watson-Crick junctions in DNA

A series of d (AT)(n) oligonucleotides containing mixtures of normal B-type Watson-Crick and antiparallel Hoogsteen helices have been studied using molecular dynamics simulation techniques to analyze the structural and thermodynamic impact of the junction between Watson-Crick and antiparallel Hoogsteen structures. Analysis of molecular dynamics simulations strongly suggests that for all oligonucleotides studied the antiparallel Hoogsteen appears as a reasonable conformation, only slightly less stable than the canonical B-type Watson-Crick one. The junctions between the Watson-Crick and Hoogsteen structures introduces a priori a sharp discontinuity in the helix, because the properties of each type of conformation are very well preserved in the corresponding fragments. However, and quite counterintuitively, junctions do not largely distort the duplex in structural, dynamics or energetic terms. Our results strongly support the possibility that small fragments of antiparallel Hoogsteen duplex might be embedded into large fragments of B-type Watson-Crick helices, making possible protein-DNA interactions that are specific of the antiparallel Hoogsteen conformation.     

Thermodynamic parameters for DNA sequences with dangling ends

The thermodynamic contributions to duplex formation of all 32 possible single-nucleotide dangling ends on a Watson-Crick pair are reported. In most instances, dangling ends are stabilizing with free energy contributions ranging from +0.48 (GT(A)) to-0.96 kcal/mol (). In comparison, Watson-Crick nearest-neighbor increments range from -0. 58 (TA/AT) to -2.24 (GC/CG) kcal/mol. Hence, in some cases, a dangling end contributes as much to duplex stability as a Watson-Crick A-T base pair. The implications of these results for DNA probe design are discussed. Analysis of the sequence dependence of dangling-end stabilities show that the nature of the closing base pair largely determines the stabilization. For a given closing base pair, however, adenine dangling ends are always more or equally as stable as the other dangling nucleotides. Moreover, 5' dangling ends are more or equally as stabilizing as their 3' counterparts. Comparison of DNA with RNA dangling-end motifs shows that DNA motifs with 5' dangling ends contribute to stability equally or more than their RNA counterparts. Conversely, RNA 3' dangling ends contribute to stability equally or more than their DNA counterparts. This data set has been incorporated into a DNA secondary structure prediction algorithm (DNA MFOLD) (http://mfold2.wustl.edu/mfold/dna/for m1.cgi) as well as a DNA hybridization prediction algorithm (HYTHERtrade mark) (http://jsl1.chem.wayne.edu/Hyther/hythermenu .html).     

Comparative Studies of the Thermodynamic Stabilities Between Sheared A:G and Watson-Crick A:U(T) Base Pairs in RNA and DNA

Abstract

Thermodynamic parameters for duplex formation were determined from CD melting curves for r(GGACGAGUCC)₂ and d(GGACGAGTCC)₂, both of which form two consecutive ‘sheared’ A:G base pairs at the center [Katahira et al. (1993) Nucleic Acids Res. 21, 5418–5424; Katahira et al., (1994) Nucleic Acids Res. 22, 2752–27591. The parameters were determined also for r(GGACUAGUCC)₂ and d(GGACTAGTCC)₂, where the A:G mismatches are replaced by Watson-Crick A:U(T) base pairs. Thermodynamic properties for duplex formation are compared between the sheared and the Watson-Crick base pairs, and between RNA and DNA. Difference in the thermodynamic stability is analyzed and discussed in terms of enthalpy and entropy changes. The characteristic features in CD spectra of RNA and DNA containing the sheared A:G base pairs are also reported.

     

N 2-methylguanosine is iso-energetic with guanosine in RNA duplexes and GNRA tetraloops.

Modified nucleotides are resource-intensive alternatives to the four nucleotides that constitute the bulk of natural RNAs. Yet, even in cases where modifications are highly conserved, their functions are difficult to identify. One possible function might be to modulate the stability of RNA structures. To investigate this possibility for N 2-methylguanosine (m2G), which is present in a wide variety of RNAs, we have determined the thermodynamic consequences of substituting m2G for G in G-C Watson-Crick pairs and G@U wobble pairs within RNA duplexes. The m2G substitution is iso-energetic with G in all cases, except for aninternal m2G@U pair, where it has a modest (0.3 kcal/mol) stabilizing effect. We have also examined theconsequences of replacing G by m2G, and A by N 6, N 6-dimethyladenosine (m26A) in the helix 45 tetraloop of 16S rRNA, which would otherwise be a standard GNRA tetraloop. This loop is a conserved, hypermethylated region of the ribosome where methylation appears to modulate activity. m26A substitution destabilizes the tetraloop, presumably because it prevents the formation of the G@A sheared pair it would otherwise contain. m2G substitution has no effect on tetraloop stability. Together, these results suggest that m2G is equally stable as either the s-cis or s-trans rotamer. The lack of a significant effect on secondary structural stability in these systems suggests that m2G is introduced into naturally occurring RNAs for reasons other than modulation of duplex stability.     

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.     

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 First Example of a Hoogsteen Basepaired DNA Duplex in Dynamic Equilibrium with a Watson-Crick Basepaired Duplex—A Structural (NMR), Kinetic and Thermodynamic Study

Abstract

A single-point substitution of the O4′ oxygen by a CH₂ group at the sugar residue of A ⁶ (i.e. 2′-deoxyaristeromycin moiety) in a self-complementary DNA duplex, 5′- d(C¹G²C³G⁴A⁵A⁶T⁷T⁸C⁹G¹⁰C¹¹G¹²)₂ ^?3, has been shown to steer the fully Watson-Crick basepaired DNA duplex (1A), akin to the native counterpart, to a doubly A ⁶:T⁷ Hoogsteen basepaired (1B) B-type DNA duplex, resulting in a dynamic equilibrium of (1A)→←(1B): K_eq = k₁/k_-1 = 0.56±0.08. The dynamic conversion of the fully Watson-Crick basepaired (1A) to the partly Hoogsteen basepaired (1B) structure is marginally kinetically and thermodynamically disfavoured [k₁ (298K) = 3.9± 0.8 sec^?1; δH°? = 164±14 kJ/mol;-TδS°? (298K) = ?92 kJ/mol giving a δG₂₉₈°? of 72 kJ/mol. Ea (k1) = 167±14 kJ/mol] compared to the reverse conversion of the Hoogsteen (1B) to the Watson-Crick (1A) structure [k-1 (298K) = 7.0±0.6 sec-1, δH°? = 153±13 kJ/mol;-TδS°? (298K) = ?82 kJ/mol giving a δG₂₉₈°? of 71 kJ/mol. Ea (k-1) = 155±13 kJ/mol]. A comparison of δG₂₉₈°? of the forward (k1) and backward (k-1) conversions, (1A)→←(1B), shows that there is ca 1 kJ/mol preference for the Watson-Crick (1A) over the double Hoogsteen basepaired (1B) DNA duplex, thus giving an equilibrium ratio of almost 2:1 in favour of the fully Watson-Crick basepaired duplex. The chemical environments of the two interconverting DNA duplexes are very different as evident from their widely separated sets of chemical shifts connected by temperature-dependent exchange peaks in the NOESY and ROESY spectra. The fully Watson-Crick basepaired structure (1A) is based on a total of 127 intra, 97 inter and 17 cross-strand distance constraints per strand, whereas the double A ⁶:T⁷ Hoogsteen basepaired (1B) structure is based on 114 intra, 92 inter and 15 cross-strand distance constraints, giving an average of 22 and 20 NOE distance constraints per residue and strand, respectively. In addition, 55 NMR-derived backbone dihedral constraints per strand were used for both structures. The main effect of the Hoogsteen basepairs in (1B) on the overall structure is a narrowing of the minor groove and a corresponding widening of the major groove. The Hoogsteen basepairing at the central A ⁶:T⁷ basepairs in (1B) has enforced a syn conformation on the glycosyl torsion of the 2′- deoxyaristeromycin moiety, A ⁶, as a result of substitution of the endocyclic 4′-oxygen in the natural sugar with a methylene group in A ⁶. A comparison of the Watson-Crick basepaired duplex (1A) to the Hoogsteen basepaired duplex (1B) shows that only a few changes, mainly in α, σ and γ torsions, in the sugar-phosphate backbone seem to be necessary to accommodate the Hoogsteen basepair.     

Thermodynamic, spectroscopic, and equilibrium binding studies of DNA sequence context effects in four 40 base pair deoxyoligonucleotides

Effects of different end sequences on melting, circular dichroism spectra (CD), and enzyme binding properties were investigated for four 40 base pair, non-self-complementary duplex DNA oligomers. The center sequences of these oligoduplexes have either of two 22 base pair modules flanked on both sides by sequences differing in AT content. Temperature-induced melting transitions monitored by differential scanning calorimetry (DSC) and ultraviolet absorbance were measured for the six duplexes in buffered 115 mM Na(+) solutions. Values of the melting transition enthalpy, DeltaH(cal), and entropy, DeltaS(cal), were obtained directly from DSC experiments. Melting transition parameters, DeltaH(vH) and DeltaS(vH), were also estimated from a van't Hoff analysis of optical melting curves collected as a function of DNA concentration, assuming that the melting transition is two-state. Melting free energies (20 degrees C) evaluated from DSC melting experiments on the four duplex DNAs ranged from -52.2 to -77.5 kcal/mol. Free energies based on the van't Hoff analysis were -37.9 to -58.8 kcal/mol. Although the values are different, trends in the melting free energies of the four duplex DNAs as a function of sequence were identical in both DSC and optical analyses. Subject to several assumptions, values for the initiation free energy were estimated for each duplex, defined as DeltaG(int) = DeltaG(cal) - DeltaG(pred), where DeltaG(cal) is the experimental free energy at 20 degrees C determined from the experimentially measured values of the transition enthalpy, DeltaH(cal), and entropy, DeltaS(cal). The predicted free energy of the sequence, DeltaG(pred)(20 degrees C), is obtained using published nearest-neighbor sequence stability values. For three of the four duplexes, values of DeltaG(int) are essentially nil. In contrast, the duplex with 81.8% GC has a considerably higher estimate of DeltaG(int) = 7.1 kcal/mol. The CD spectra for the six duplexes collected over the wavelength range from 200 to 320 nm are also sequence-dependent. Factor analysis of the CD spectra by singular value decomposition revealed that the experimental CD spectra could be reconstructed from linear combinations of two minor and one major subspectra. Changes in the coefficients of the major subspectrum for different sequences reflect incremental sequence-dependent variations of the CD spectra. Equilibrium binding by BamHI restriction endonuclease to the 40 base pair DNAs whose central eight base pairs contain the recognition sequence for BamHI restriction enzyme bounded by A.T base pairs, 5'-A-GGATCC-A-3' was investigated. Binding assays were performed by titering BamHI against a constant concentration of each of the duplex DNA substrates, in the absence of Mg(2+), followed by analysis by gel retardation. Under the conditions employed, the enzyme binds but does not cleave the DNAs. Results of the assays revealed two binding modes with retarded gel mobilities. Binding isotherms for the fraction of bound DNA species versus enzyme concentration for each binding mode were constructed and analyzed with a simple two-step equilibrium binding model. This analysis provided semiquantitative estimates on the equilibrium binding constants for BamHI to the four DNAs. Values obtained for the binding constants varied only 7-fold and ranged from 6 x 10(-)(8) to 42 x 10(-)(8) M, with binding free energies from -8.6 to -9.7 (+/- 0.2) kcal/mol depending on the sequence that flanks the enzyme binding site. Unlike what was found earlier in binding studies of the 22 base pair duplexes that constitute the core modules of the present 40-mers [Riccelli, P. V., Vallone, P. M., Kashin, I., Faldasz, B. D., Lane, M. J., and Benight, A. S. (1999) Biochemistry 38, 11197-11208], no obvious relationship between binding and stability was found for these longer DNAs. Apparently, effects of flanking sequence stability on restriction enzyme binding may only be measurable in very short duplex deoxyoligonucl     

Thermodynamics of unpaired terminal nucleotides on short RNA helixes correlates with stacking at helix termini in larger RNAs.

Free energies for stacking of unpaired nucleotides (dangling ends) at the termini of oligoribonucleotide Watson-Crick helixes (DeltaG(0)37,stack) depend on sequence for 3' ends but are always small for 5' ends. Here, these free energies are correlated with stacking at helix termini in a database of 34 RNA structures determined by X-ray crystallography and NMR spectroscopy. Stacking involving GA pairs is considered separately. A base is categorized as stacked by its distance from (相似文献