The solution structure of a d[C(TTCG)G] DNA hairpin and comparison to the unusually stable RNA analogue.

The solution structure of the DNA analogue of the unusually stable r[C(UUCG)G] RNA hairpin, 5'-d[GGA-C(TTCG)GTCC]-3', has been determined by NMR spectroscopy, and its structure has been compared to that of the RNA molecule. The RNA molecule is compact and rigid with a highly structured loop. However, the DNA molecule is much less structured. The DNA hairpin contains a B-form stem of four base pairs. The terminal base pair frays, and the 3'-terminal nucleotides, C11 and C12, are in equilibrium between 2'-endo and 3'-endo conformations. Unlike the RNA loop, the DNA loop contains no syn nucleotides, and there is no evidence for base-base or base-phosphate hydrogen bonding in the loop. The loop is flexible, and reveals no specific internucleotide interactions.     

Experimental and computational studies of the G[UUCG]C RNA tetraloop

In prokaryotic ribosomal RNAs, most UUCG tetraloops are closed by a C-G base-pair. However, this preference is greatly reduced in eukaryotic rRNA species where many UUCG tetraloops are closed by G-C base-pairs. Here, biophysical properties of the C[UUCG]G and G[UUCG]C tetraloops are compared, using experimental and computational methods. Thermal denaturation experiments are used to derive thermodynamic parameters for the wild-type G[UUCG]C tetraloop and variants containing single deoxy substitutions in the loop. A comparison with analogous experiments on the C[UUCG]G motif shows that the two RNA species exhibit similar patterns in response to the substitutions, suggesting that their loop structures are similar. This conclusion is supported by NMR data that suggest that the essential UUCG loop structure is maintained in both tetraloops. However, NMR results show that the G[UUCG]C loop structure is disrupted prior to melting of the stem; this behavior is in contrast to the two-state behavior of the C[UUCG]G molecule. Stochastic dynamics simulations using the GB/SA continuum solvation model, run as a function of temperature, show rare conformational transitions in several G[UUCG]C simulations. These results lead to the conclusion that substitution of a G-C for a C-G closing base-pair increases the intrinsic flexibility of the UUCG loop.     

Extra stable 2',5'-linked RNA loops

We prepared hairpins that differ in the connectivity of phosphodiester linkages in the loop (RNA vs 2', 5'-RNA). We find that the stability of the extra stable RNA hairpin 5'-rGGAC(UUCG)GUCC-3' is the same as that observed for the hairpin containing a 2',5'RNA loop, i.e. 5'-rGGAC(UUCG)GUCC-3' (where UUCG = U2'p5'U2'p5' C2'p5'G2'p5'). Also significant is the finding that when the stem is duplex DNA, duplex 2',5'-RNA, or DNA:2',5'-RNA, hairpins with the UUCG loop are more stable than those with UUCG loop.     

Structure of 2',5'-linked tetraribonucleotide loops: a novel RNA motif

We report on the three dimensional structure of an RNA hairpin containing a 2',5'-linked tetraribonucleotide loop, namely, 5'-rGGAC(UUCG)GUCC-3' (where UUCG = U(2'p5')U(2'p5')C(2'p5')G(2'p5')). We show that the 2',5'-linked RNA loop adopts a conformation that is quite different from that previously observed for the native 3',5'-linked RNA loop. The 2',5'-RNA loop is stabilized by (a) U:G wobble base pairing, with both bases in the anti conformation, (b) extensive base stacking, and (c) sugar-base contacts, all of which contribute to the extra stability of this hairpin structure.     

Use of ultra stable UNCG tetraloop hairpins to fold RNA structures: thermodynamic and spectroscopic applications.

RNA molecules of > 20 nucleotides have been the focus of numerous recent NMR structural studies. Several investigators have used the UNCG family of hairpins to ensure proper folding. We show that th UUCG hairpin has a minimum requirement of a two base-pair stem. Hairpins with a CG loop closing base pair and an initial 5'CG or 5'GC base pair have a melting temperature approximately 55 degrees C in 10 mM sodium phosphate. The high stability of even such small hairpins suggests that the hairpin can serve as a nucleation site for folding. For high resolution NMR work, the UNCG loop family (UACG in particular) provides excellent spectroscopic markers in one-dimensional exchangeable spectra, in two-dimensional COSY spectra and in NOESY spectra that clearly define it as forming a hairpin. This allows straightforward initiation of chemical shift assignments.     

Preference for ribose over deoxyribose in loop-closing base pairs of extra stable nucleic acid hairpins

We have investigated the effect of switching ribose to deoxyribose at the closing base-pair of an extra-stable RNA hairpin. Specifically, we studied the sequence 5'-GGAC(UUCG)GUCC, a dodecanucleotide that folds into a well-characterized, "extra stable" RNA hairpin structure. Recently, we showed that hairpins containing a 2',5'-linked (UUCG) loop instead of the native 3',5'-linked loop also exhibit extra-stability (Hannoush and Damha, J. Am. Chem. Soc., 2001, 123, 12368-12374). In this article, we show that the ribose units located at the loop-closing positions (i.e., rC4 and rG9) contribute significantly to the stabilization of RNA hairpins, particularly those containing the 3',5'-UUCG loop. Interestingly, the requirement of rC4 and rG9 is more relaxed for DNA hairpins containing the 2',5'-UUCC loop and, in fact, they may be replaced altogether (ribose--> deoxyribose) without affecting stability. The results broaden our understanding of the behavior of highly stable (UUCG) hairpin loops and how they respond to structural perturbation of the loop-closing base pairs.     

NMR study of a heterochiral DNA hairpin:impact of L-enantiomery in the loop.

We carried out a structural study of the DNA heterochiral strand d (AGCTTATCAT(L)CGATAAGCT), -AT(L)C-, where T(L) (L thymine ) replaces T (natural D-thymine). -AT(L)C- is a structural analog of -ATC- that belongs to a strong topoisomerase II DNA cleavage site and which has been shown to resolve into a hairpin structure with a stem formed by eight Waston-Crick base-pairs and a single residue loop closed by an A.C sheared base-pair. Although - AT(L)C-, like its parent -ATC-, folds into a hairpin structure at low and high DNA concentrations it displays a lower stability (Tm of 56 degrees C versus 58.5 degrees C). Several NMR features in -AT(L)C- account for the disruption of the A.C pairing in the loop and a weakening of the C.G base-pair stability at the stem-loop junction. For instance, the exchange of the loop imino protons with solvent is accelerated compared with the natural oligonucleotide -ATC-. The higher flexibility of the heterochiral loop is confirmed by the results of NMR restrained molecular dynamics. In the calculated final structures of -AT(L)C-, the T10(L) residue moves the A9 and C11 residues away, thus preventing the loop closure through a C.A sheared base-pair and the achievement of a good base-base or sugar-base stacking. Actually, most of the stabilizing interactions present in -ATC- are lost in the heterochiral - AT(L)C- explaining its weaker stability.     

Understanding the thermodynamic stability of an RNA hairpin and its mutant.

The thermodynamic stability of RNA hairpin loops has been a subject of considerable interest in the recent past (Wimberly et al., 1991). There have been experimental reports indicating that the hairpins with a C(UUCG)G loop sequence are thermodynamically very stable (Wimberly et al., 1991). We used the solution structure of GGAC(UUCG)GUCC (Cheong et al., 1990; Varani et al., 1991) as the starting conformation in our attempt to understand its thermodynamic stability. We carried out molecular dynamics/free energy simulations to understand the basis for the destabilization of the C(UUCG)G loop by mutating cytosine (C7)-->uracil. Because of the limited length of simulation and the presence of kinetic barriers (solvent intervention) to the uracil-->cytosine mutation, all of our computed free energy differences are based on multiple forward simulations. Based on these calculations we find that the cytosine-->uracil mutation in the loop destabilizes it by approximately 1.5kcal/mol relative to that of the reference state, an A-form RNA but with cytosine (C7) looped out. This is the same sign and magnitude as that observed in the thermodynamic studies carried out by Varani et al.(1991). We have carried out free energy component analysis to understand the effect of mutating the cytosine residue to uracil on the thermodynamic stability of the C(UUCG)G hairpin loops. Our calculations show that the most significant contribution to the stability is from the phosphate group linking U5 and U6, which favors the cytosine residue over uracil by about 6.0 kcal/mol. The residues U5, U6, and G8 in the loop region also contribute significantly to the stability. The contributions from the salt and solvent compensate each other, indicating the dynamic nature of interactions of the environment with the nucleic acid system and the coupling between these two components.     

A thermodynamic study of unusually stable RNA and DNA hairpins.

About 70% of the RNA tetra-loop sequences identified in ribosomal RNAs from different organisms fall into either (UNCG) or (GNRA) families (where N = A, C, G, or U; and R = A or G). RNA hairpins with these loop sequences form unusually stable tetra-loop structures. We have studied the RNA hairpin GGAC(UUCG)GUCC and several sequence variants to determine the effect of changing the loop sequence and the loop-closing base pair on the thermodynamic stability of (UNCG) tetra-loops. The hairpin GGAG(CUUG)CUCC with the conserved loop G(CUUG)C was also unusually stable. We have determined melting temperatures (Tm), and obtained thermodynamic parameters for DNA hairpins with sequences analogous to stable RNA hairpins with (UNCG), C(GNRA)G, C(GAUA)G, and G(CUUG)C loops. DNA hairpins with (TTCG), (dUdUCG), and related sequences in the loop, unlike their RNA counterparts, did not form unusually stable hairpins. However, DNA hairpins with the consensus loop sequence C(GNRA)G were very stable compared to hairpins with C(TTTT)G or C(AAAA)G loops. The C(GATA)G and G(CTTG)C loops were also extra stable. The relative stabilities of the unusually stable DNA hairpins are similar to those observed for their RNA analogs.     

Preference for Ribose Over Deoxyribose in Loop-Closing Base Pairs of Extra Stable Nucleic Acid Hairpins

We have investigated the effect of switching ribose to deoxyribose at the closing base-pair of an extra-stable RNA hairpin. Specifically, we studied the sequence 5′-GGAC(UUCG)GUCC, a dodecanucleotide that folds into a well-characterized, “extra stable” RNA hairpin structure. Recently, we showed that hairpins containing a 2′,5′-linked (UUCG) loop instead of the native 3′,5′-linked loop also exhibit extra-stability (Hannoush and Damha, J. Am. Chem. Soc., 2001, 123, 12368–12374). In this article, we show that the ribose units located at the loop-closing positions (i.e., rC ₄ and rG ₉ ) contribute significantly to the stabilization of RNA hairpins, particularly those containing the 3′,5′-UUCG loop. Interestingly, the requirement of rC4 and rG9 is more relaxed for DNA hairpins containing the 2′,5′-UUCG loop and, in fact, they may be replaced altogether (ribose → deoxyribose) without affecting stability. The results broaden our understanding of the behavior of highly stable (UUCG) hairpin loops and how they respond to structural perturbation of the loop-closing base pairs.     

Structural characterization of a six-nucleotide RNA hairpin loop found in Escherichia coli, r(UUAAGU)

The binding region of the Escherichia coli S2 ribosomal protein contains a conserved UUAAGU hairpin loop. The structure of the hairpin formed by the oligomer r(GCGU4U5A6A7G8U9CGCA), which has an r(UUAAGU) hairpin loop, was determined by NMR and molecular modeling techniques as part of a study aimed at characterizing the structure and thermodynamics of RNA hairpin loops. Thermodynamic data obtained from melting curves for this RNA oligomer show that it forms a hairpin in solution with the following parameters: DeltaH degrees = -42.8 +/- 2.2 kcal/mol, DeltaS degrees = -127.6 +/- 6.5 eu, and DeltaG degrees (37) = -3.3 +/- 0.2 kcal/mol. Two-dimensional NOESY WATERGATE spectra show an NOE between U imino protons, which suggests that U4 and U9 form a hydrogen bonded U.U pair. The U5(H2') proton shows NOEs to both the A6(H8) proton and the A7(H8) proton, which is consistent with formation of a "U" turn between nucleotides U5 and A6. An NOE between the A7(H2) proton and the U9(H4') proton shows the proximity of the A7 base to the U9 sugar, which is consistent with the structure determined for the six-nucleotide loop. In addition to having a hydrogen-bonded U.U pair as the first mismatch and a U turn, the r(UUAAGU) loop has the G8 base protruding into the solvent. The solution structure of the r(UUAAGU) loop is essentially identical to the structure of an identical loop found in the crystal structure of the 30S ribosomal subunit where the guanine in the loop is involved in tertiary interactions with RNA bases from adjacent regions [Wimberly, B. T., Brodersen, D. E., Clemons, W. M., Morgan-Warren, R. J., Carter, A. P., Vonrhein, C., Hartsch, T., and Ramakrishnan, V. (2000) Nature 407, 327-339]. The similarity of the solution and solid-state structures of this hairpin loop suggests that formation of this hairpin may facilitate folding of 16S RNA.     

Solution structure of a purine rich hexaloop hairpin belonging to PGY/MDR1 mRNA and targeted by antisense oligonucleotides

A preferential target of antisense oligonucleotides directed against human PGY/MDR1 mRNA is a hairpin containing a stem with a G*U wobble pair, capped by the purine-rich 5'r(GGGAUG)3' hexaloop. This hairpin is studied by multidimensional NMR and restrained molecular dynamics, with special emphasis on the conformation of south sugars and non-standard phosphate linkages evidenced in both the stem and the loop. The hairpin is found to be highly structured. The G*U wobble pair, a strong counterion binding site, displays structural particularities that are characteristic of this type of mismatch. The upper part of the stem undergoes distortions that optimize its interactions with the beginning of the loop. The loop adopts a new fold in which the single-stranded GGGA purine tract is structured in A-like conformation stacked in continuity of the stem and displays an extensive hydrogen bonding surface for recognition. The remarkable hairpin stability results from classical inter- and intra-strand interactions reinforced by numerous hydrogen bonds involving unusual backbone conformations and ribose 2'-hydroxyl groups. Overall, this work emphasizes numerous features that account for the well-ordered structure of the whole hairpin and highlights the loop properties that facilitate interaction with antisense oligonucleotides.     

Thermodynamics of 2'-ribose substitutions in UUCG tetraloops

The ribose 2'-hydroxyl group confers upon RNA many unique molecular properties. To better appreciate its contribution to structure and stability and to monitor how substitutions of the 2' hydroxyl can alter an RNA molecule, each loop pyrimidine ribonucleotide in the UUCG tetraloop was substituted with a nucleotide containing either a fluorine (2'-F), hydrogen (2'-H), amino (2'-NH2), or methoxy (2'-OCH3) group, in the context of both the C:G and G:C loop-closing base pair. The thermodynamic parameters of these tetraloop variants have been determined and NMR experiments used to monitor the structural changes resulting from the substitutions. The modified riboses are better tolerated in the G[UUCG]C tetraloop, which may be due to its increased loop flexibility relative to the C[UUCG]G loop. Even for these simple substitutions, the free-energy change reflects a complex interplay of hydrogen bonding, solvation effects, and intrinsic pucker preferences of the nucleotides.     

EXTRA STABLE 2′,5′-LINKED RNA LOOPS

We prepared hairpins that differ in the connectivity of phosphodiester linkages in the loop (RNA vs 2′, 5′-RNA). We find that the stability of the extra stable RNA hairpin 5′-rGGAC(UUCG)GUCC-3′ is the same as that observed for the hairpin containing a 2′,5′RNA loop, i.e. 5′-rGGAC(UUCG)GUCC-3′ (where UUCG = U_2′p5′U_2′p5′ C_2′p5′G_2′p5′). Also significant is the finding that when the stem is duplex DNA, duplex 2′,5′-RNA, or DNA:2′,5′-RNA, hairpins with the UUCG loop are more stable than those with UUCG loop.     

Continuum solvent studies of the stability of RNA hairpin loops and helices

We apply continuum solvent models to investigate the relative stability of various conformational forms for two RNA sequences, GGAC(UUCG)GUCC and GGUG(UGAA)CACC. In the first part, we compare alternate hairpin conformations to explore the reliability of these models to discriminate between different local conformations. A second part looks at the hairpin-duplex conversion for the UUCG sequence, identifying major contributors to the thermodynamics of a much large scale transition. Structures were taken as snapshots from multi-nanosecond molecular dynamics simulations computed in a consistent fashion using explicit solvent and with long-range electrostatics accounted for using the Particle-Mesh Ewald procedure. The electrostatic contribution to solvation energies were computed using both a finite-difference Poisson-Boltzmann (PB) model and a pairwise Generalized Born model; non-electrostatic contributions were estimated with a surface-area dependent term. To these solvation free energies were added the mean solute internal energies (determined from a molecular mechanics potential) and estimates of the solute entropy (from a harmonic analysis). Consistent with experiment and with earlier solvated molecular dynamics simulations, the UUCG hairpin was found to prefer conformers close to a recent NMR structure determination in preference to those from an earlier NMR study. Similarly, results for the UGAA hairpin favored an NMR-derived structure over that to be expected for a generic GNRA hairpin loop. Experimental free energies are not known for the hairpin/duplex conversion, but must be close to zero since hairpins are seen in solution and duplexes in crystals; out calculations find a value near zero and illustrate the expected interplay of solvation, salt effects and entropy in affecting this equilibrium.     

Conformational and model-building studies of the hairpin form of the mismatched DNA octamer d(m5C-G-m5C-G-T-G-m5C-G)

The hairpin form of the mismatched octamer d(m5C-G-m5C-G-T-G-m5C-G) was studied by means of NMR spectroscopy. In a companion study it is shown that the hairpin form of this DNA fragment consists of a structure with a stem of three Watson-Crick-type base pairs and a loop consisting of only two nucleotides. The non-exchangeable proton resonances were assigned by means of two-dimensional correlation spectroscopy and two-dimensional nuclear Overhauser effect spectroscopy. Proton-proton coupling constants were used for the conformational analysis of the deoxyribose ring and for some of the backbone torsion angles. From the two-dimensional NMR spectra and the coupling-constant analysis it is concluded that: (i) the stem of the hairpin exhibits B-DNA characteristics; (ii) the sugar rings are not conformationally pure, but display a certain amount of conformational flexibility; (iii) the stacking interaction in the stem of the hairpin is elongated from the 3'-side in a more or less regular fashion with the two loop nucleotides; (iv) at the 5'-side of the stem a stacking discontinuity occurs between the stem and the loop; (v) at the 5'-side of the stem the loop is closed by means of a sharp backbone turn which involves unusual gamma' and beta+ torsion angles in residue dG(6). The NMR results led to the construction of a hairpin-loop model which was energy-minimized by means of a molecular-mechanics program. The results clearly show that a DNA hairpin-loop structure in which the loop consists of only two nucleotides bridging the minor groove in a straightforward fashion, (i) causes no undue steric strain, and (ii) involves well-known conformational principles throughout the course of the backbone. The hairpin form of the title compound is compared with the hairpin form of d(A-T-C-C-T-A-T4-T-A-G-G-A-T), in which the central -T4- part forms a loop of four nucleotides. Both models display similarities as far as stacking interactions are concerned.     

Solution conformation of an RNA hairpin loop.

The hairpin conformation adopted by the RNA sequence 5'GCGAUUUCUGACCGCC3' has been studied by one- and two-dimensional NMR spectroscopy. Exchangeable imino spectra in 60 mM Na+ indicate that the hairpin has a stem of six base pairs (indicated by boldface type) and a loop of three nucleotides. NOESY spectra of nonexchangeable protons confirm the formation of the stem region. The duplex has an A-conformation and contains an A.C apposition; a G.U base pair closes the loop region. The stem nucleotides have C3'-endo sugar conformations, as expected of an A-form duplex, whereas the three loop nucleotides adopt C2'-endo sugar puckers. Stacking within the loop, C8 upon the sugar of U7, stabilizes the structure. The pH dependence of both the exchangeable and nonexchangeable NMR spectra is consistent with the formation of an A+.C base pair, protonated at the N1 position of adenine. The stability of the hairpin was probed by using absorbance melting curves. The hairpin structure with the A+.C base pair is about +2 kcal/mol less stable in free energy at 37 degrees C than the hairpin formed with an A.U pair replacing the A+.C pair.     

Multinuclear NMR studies of DNA hairpins. 1. Structure and dynamics of d(CGCGTTGTTCGCG)

The solution structure of the hairpin formed by d(CGCGTTGTTCGCG) has been examined in detail by a wide variety of NMR techniques. The hairpin was characterized by proton NMR to obtain interproton distances and torsion angle information. An energy-minimized model was constructed that is consistent with these data. The hairpin consists of a B-DNA stem of four C-G base pairs and a loop region consisting of five unpaired bases. Three bases in the 5' of the loop are stacked over the 3' end of the stem, and the other two bases in the 3' of the loop are stacked over the 5' end of the stem. The phosphorus NMR spectrum revealed a phosphate in the stem region with an unusual conformation, and two phosphates, P9 and P10, were found to undergo intermediate exchange between conformations. The hairpin was also synthesized with a carbon-13 label in each of the thymidine C6 carbons, and relaxation measurements were performed to determine the extent of internal motions in the loop region. The loop bases are more flexible than the stem bases and exhibit subnanosecond motions with an amplitude corresponding to diffusion in a cone of approximately 30 degrees.