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.     

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.     

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.     

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.     

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.     

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.     

The structure of the L3 loop from the hepatitis delta virus ribozyme: a syn cytidine.

The structure of the L3 central hairpin loop isolated from the antigenomic sequence of the hepatitis delta virus ribozyme with the P2 and P3 stems from the ribozyme stacked on top of the loop has been determined by NMR spectroscopy. The 26 nt stem-loop structure contains nine base pairs and a 7 nt loop (5'-UCCUCGC-3'). This hairpin loop is critical for efficient catalysis in the intact ribozyme. The structure was determined using homonuclear and heteronuclear NMR techniques on non-labeled and15N-labeled RNA oligonucleotides. The overall root mean square deviation for the structure was 1.15 A (+/- 0.28 A) for the loop and the closing C.G base pair and 0.90 A (+/- 0.18 A) for the loop and the closing C.G base pair but without the lone purine in the loop, which is not well defined in the structure. The structure indicates a U.C base pair between the nucleotides on the 5'- and 3'-ends of the loop. This base pair is formed with a single hydrogen bond involving the cytosine exocyclic amino proton and the carbonyl O4 of the uracil. The most unexpected finding in the loop is a syn cytidine. While not unprecedented, syn pyrimidines are highly unusual. This one can be confidently established by intranucleotide distances between the ribose and the base determined by NMR spectroscopy. A similar study of the structure of this loop showed a somewhat different three-dimensional structure. A discussion of differences in the two structures, as well as possible sites of interaction with the cleavage site, will be presented.     

Effects of osmolytes on stable UUCG tetraloops and their preference for a CG closing base pair

Osmolytes have the potential to affect the stability of secondary structure motifs and alter preferences for conserved nucleic acid sequences in the cell. To contribute to the understanding of the in vivo function of RNA we observed the effects of different classes of osmolytes on the UNCG tetraloop motif. UNCG tetraloops are the most common and stable of the RNA tetraloops and are nucleation sites for RNA folding. They also have a significant thermodynamic preference for a CG closing base pair. The thermal denaturation of model hairpins containing UUCG loops was monitored using UV-Vis spectroscopy in the presence of osmolytes with different chemical properties. Interestingly, all of the osmolytes tested destabilized the hairpins, but all had little effect on the thermodynamic preference for a CG base pair, except for polyethylene glycol (PEG) 200. PEG 200 destabilized the loop with the CG closing base pair relative to the loop with a GC closing base pair. The destabilization was linear with increasing concentrations of PEG 200, and the slope of this relationship was not perturbed by changes in the hairpin stem outside of the closing pair. This result suggests that in the presence of PEG 200, the UUCG loop with a GC closing base pair may retain some preferential interactions with the cosolute that are lost in the presence of the CG closing base pair. These results reveal that relatively small structural changes may influence how osmolytes tune the stability, and thus the function of a secondary structure motif in vivo.     

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.     

Conformation and dynamics of an RNA internal loop

The conformation and the dynamics of an RNA oligonucleotide (26 nucleotides) which is a model for loop E in eukaryotic 5S RNA have been investigated by one- and two-dimensional NMR. The central portion of the oligonucleotide contains two G A oppositions, a common feature of ribosomal RNAs. The exchangeable proton spectrum indicates that an internal loop separates two stems of four and five base pairs. This observation is not consistent with structures for loop E containing mismatched G.A base pairs proposed from chemical and enzymatic studies on Xenopus laevis 5S RNA. The nonexchangeable proton spectrum has been assigned by two-dimensional NMR. Scalar couplings from correlated experiments and interproton distances from NOESY experiments at short mixing times have been used to determine glycosidic angles, sugar puckers, and other conformational features. The conformation of the stems is very close to standard A-form RNA, and extensive base stacking continues into the internal loop. This result provides a structural basis for the large favorable enthalpy of duplex formation determined in thermodynamic studies. Unusual structural and dynamic features are localized in the nucleotides connecting the loop to the stems.     

Structural characterization of d(CAACCCGTTG) and d(CAACGGGTTG) mini-hairpin loops by heteronuclear NMR: the effects of purines versus pyrimidines in DNA hairpins.

The DNA decamers, d(CAACCCGTTG) and d(CAACGGGTTG) were studied in solution by proton and heteronuclear NMR. Under appropriate conditions of pH, temperature, salt concentration and DNA concentration, both decamers form hairpin conformations with similar stabilities [Avizonis and Kearns (1995) Biopolymers, 35, 187-200]. Both decamers adopt mini-hairpin loops, where the first and last four nucleotides are involved in Watson-Crick hydrogen bonding and the central two nucleotides, CC or GG respectively, form the loop. Through the use of proton-proton, proton-phosphorus and natural abundance proton-carbon NMR experiments, backbone torsion angles (beta, gamma and epsilon), sugar puckers and interproton distances were measured. The nucleotides forming the loops of these decamers were found to stack upon one another in an L1 type of loop conformation. Both show gamma tr and unusual beta torsion angles in the loop-closing nucleotide G7, as expected for mini-hairpin loop formation. Our results indicate that the beta and epsilon torsion angles of the fifth and sixth nucleotides that form the loop and the loop-closing nucleotide G7 are not in the standard trans conformation as found in B-DNA. Although the loop structures calculated from NMR-derived constraints are not well defined, the stacking of the bases in the two different hairpins is different. This difference in the base stacking of the loop may provide an explanation as to why the cytosine-containing hairpin is thermodynamically more stable than the guanine-containing hairpin.     

Solution structure of a GAAG tetraloop in helix 6 of SRP RNA from Pyrococcus furiosus

The NMR structure of a 12-mer RNA derived from the helix 6 of SRP RNA from Pyrococcus furiosus, whose loop-closing base pair is U.G, was determined, and the structural and thermodynamic properties of the RNA were compared with those of a mutant RNA with the C:G closing base pair. Although the structures of the two RNAs are similar to each other and adopt the GNRR motif the conformational stabilities are significantly different to each other It was suggested that weaker stacking interaction of the GAAG loop with the U:G closing base pair in 12-mer RNA causes the lower conformational stability.     

Structure and function of the conserved 690 hairpin in Escherichia coli 16 S ribosomal RNA: analysis of the stem nucleotides

Nucleotides 680 to 710 of Escherichia coli 16 S rRNA form a distinct structural domain required for ribosome function. The goal of this study was to determine the functional significance of pairing interactions in the 690 region. Two different secondary structures were proposed for this hairpin, based on phylogenetic and chemical modification studies. To study the effect of pairing interactions in the 690 hairpin on ribosome function and to determine which of the proposed secondary structures is biologically significant, we performed an instant-evolution experiment in which the nine nucleotides that form the proposed base-pairs and dangling ends of the 690 stem were randomly mutated, and functional mutant combinations were selected. A total of 96 unique functional mutants were isolated, assayed in vivo, and sequenced. Analysis of these data revealed extensive base-pairing and stacking interactions among the mutated nucleotides. Formation of either a Watson-Crick base-pair or G.U pair between positions 688 and 699 is absolutely required for ribosome function. We also performed NMR studies of a 31-nucleotide RNA which indicate the formation of a functionally important base-pair between nucleotides 688 and 699. Formation of a second base-pair between positions 689 and 698, however, is not essential for ribosome function, but the level of ribosome function correlates with the predicted thermodynamic stability of the nucleotide pairs in these positions. The universally conserved positions G690 and U697 are generally portrayed as forming a G.U mismatch. Our data show co-variation between these positions, but do not support the hypothesis that the G690:U697 pair forms a wobble structure. NMR studies of model 14-nt and 31-nt RNAs support these findings and show that G690 and U697 are involved in unusual stacking interactions but do not form a wobble pair. Preliminary NMR structural analysis reveals that the loop portion of the 690 hairpin folds into a highly structured and novel conformation.     

Studies on the structure and stabilizing factor of the CUUCGG hairpin RNA using chemically synthesized oligonucleotides.

A tridecaribonucleotide, r-UGAGCUUCGGCUC, and two analogues r(UGAGC)d(UUCG)r(GCUC) and r-UGAGCUUCIGCUC, which form a hairpin structure with a four-base-paired stem and a UUCG loop, were synthesized by the solid-phase phosphoramidite method. Properties of these three oligomers and d-TGAGCTTCGGCTC, the DNA analogue, were studied by UV, CD and NMR spectroscopy. The melting temperature (Tm) data suggest that the 2'-hydroxy1 groups and the 2-amino group of guanosine in the loop (9G) stabilize the CUUCGG hairpin which is known to have an unusually high Tm. NMR studies show that this 9G takes a syn conformation and the phosphodiester backbone has a turn at 9G-10G which is a junction of the stem and loop.     

The three-dimensional structure of a DNA hairpin in solution two-dimensional NMR studies and structural analysis of d(ATCCTATTTATAGGAT)

The hairpin formed by d(ATCCTATTTATAGGAT) was studied by means of two-dimensional NMR spectroscopy and conformational analysis. Almost all 1H resonances of the stem region could be assigned, while the 1H and 31P spectra of the loop region were interpreted completely; this includes the stereospecific assignment of the H5' and H5" resonances. The derivation of the detailed loop structure was carried out in a stepwise fashion including some improved and new methods for structure determination from NMR data. In the first step, the mononucleotide structures were examined. The conformational space available to the mononucleotide was scanned systematically by varying the glycosidic torsion angle and pseudorotational parameters. Each generated conformer was tested against the experimental J coupling constants and NOE parameters. In the following stage, the structures of dinucleotides and longer fragments were derived. Inter-residue distances between protons were calculated by means of a procedure in which the simulated NOEs, obtained via a relaxation-matrix approach, were fitted to the experimental NOEs without the introduction of a molecular model. In addition, the backbone torsion angles beta, gamma and epsilon were deduced from homocoupling and heterocoupling constants. These data served as constraints in the next step, in which the loop sequence was subjected to a multi-conformer generation procedure. The resulting structures were tested against the mentioned constraints and disregarded if these constraints were violated. This yielded a family of structures for the loop region, confined to a relatively narrow conformational space. A representative conformation was subsequently docked on a B-type stem which fulfilled the structural constraints (derived from the NMR experiments for the stem region) to yield the hairpin structure. Results obtained from subsequent restrained-molecular-mechanics as well as free-molecular-mechanics calculations are in accordance with those obtained by means of the analysis described above. The structure of the hairpin loop is a compactly folded conformation and the first base of the central TTTA region forms a Hoogsteen T-A pair with the fourth base. This Hoogsteen base pair is stacked upon the sixth base pair of the B-type double-helical stem. The second base of the loop is folded into the minor groove, whereas the third base of the loop is partly stacked on the first and fourth bases. The phosphate backbone exhibits a sharp turn between the third and fourth nucleotides of the loop. The peculiar structure of this hairpin loop is discussed in relation to loop folding in DNA and RNA hairpins and in relation to a general model for loop folding.     

Solution Structure of a GAAG Tetraloop in Helix 6 of SRP RNA from Pyrococcus furiosus

The NMR structure of a 12-mer RNA derived from the helix 6 of SRP RNA from Pyrococcus furiosus, whose loop-closing base pair is U:G, was determined, and the structural and thermodynamic properties of the RNA were compared with those of a mutant RNA with the C:G closing base pair. Although the structures of the two RNAs are similar to each other and adopt the GNRR motif, the conformational stabilities are significantly different to each other. It was suggested that weaker stacking interaction of the GAAG loop with the U:G closing base pair in 12-mer RNA causes the lower conformational stability.     

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.     

A phenanthrene modified RNA hairpin

The influence of hairpin loop replacement with the phenanthrene moiety in RNA was investigated. The stability of this novel structure was compared to a hairpin with a U(4) loop, an extra stable tetra-loop (UUCG), and an analogous phenanthrene modified DNA hairpin. Thermal denaturation experiments and CD spectra were used to study the structure and stability of the modified hairpin.     

Non-nearest-neighbor dependence of stability for group III RNA single nucleotide bulge loops

Thirty-five RNA duplexes containing single nucleotide bulge loops were optically melted and the thermodynamic parameters for each duplex determined. The bulge loops were of the group III variety, where the bulged nucleotide is either a AG/U or CU/G, leading to ambiguity to the exact position and identity of the bulge. All possible group III bulge loops with Watson–Crick nearest-neighbors were examined. The data were used to develop a model to predict the free energy of an RNA duplex containing a group III single nucleotide bulge loop. The destabilization of the duplex by the group III bulge could be modeled so that the bulge nucleotide leads to the formation of the Watson–Crick base pair rather than the wobble base pair. The destabilization of an RNA duplex caused by the insertion of a group III bulge is primarily dependent upon non-nearest-neighbor interactions and was shown to be dependent upon the stability of second least stable stem of the duplex. In-line structure probing of group III bulge loops embedded in a hairpin indicated that the bulged nucleotide is the one positioned further from the hairpin loop irrespective of whether the resulting stem formed a Watson–Crick or wobble base pair. Fourteen RNA hairpins containing group III bulge loops, either 3′ or 5′ of the hairpin loop, were optically melted and the thermodynamic parameters determined. The model developed to predict the influence of group III bulge loops on the stability of duplex formation was extended to predict the influence of bulge loops on hairpin stability.     

RNA hairpin loop stability depends on closing base pair.

Thermodynamic parameters are reported for hairpin formation in 1 M NaCl by RNA sequences of the type GGXAUAAUAYCC, where X and Y are CG, GC, AU, UA, GU, or UG. A nearest neighbor analysis of the data indicates the free energy change for loop formation at 37 degrees C, delta degrees Gl,37, averages 3.4 kcal/mol for hairpin loops closed with C.G, G.C, and G.U pairs. In contrast, delta G degree l,37 averages 4.6 kcal/mol for loops closed with A.U, U.A, or U.G pairs. Thus the stability of an RNA hairpin depends on the closing base pair. The hairpin with a GA mismatch that is formed by GGCGUAAUAGCC is more stable than the corresponding hairpin with an AA mismatch. Thus hairpin stability also depends on loop sequence. These effects are not included in current algorithms for prediction of RNA structure from sequence.