The energetics of small internal loops in RNA

The energetics of small internal loops are important for prediction of RNA secondary and tertiary structure, selection of drug target sites, and understanding RNA structure-function relationships. Hydrogen bonding, base stacking, electrostatic interactions, backbone distortion, and base-pair size compatibility all contribute to the energetics of small internal loops. Thus, the sequence dependence of these energetics are idiosyncratic. Current approximations for predicting the free energies of internal loops consider size, asymmetry, closing base pairs, and the potential to form GA, GG, or UU pairs. The database of known three-dimensional structures allows for comparison with the models used for predicting stability from sequence.     

Thermodynamic stability and structural features of the J4/5 loop in a Pneumocystis carinii group I intron

The J4/5 loop of the group I intron in the mouse-derived fungal pathogen Pneumocystis carinii is the docking site for the first step of the RNA-catalyzed self-splicing reaction and thus is a model of a potential drug target. This purine-rich asymmetric internal loop, 5'GGAAG/3'UAGU, is also thermodynamically more stable than other internal loops with two GU closing pairs and three nucleotides opposite two nucleotides. The results from optical melting, nuclear magnetic resonance spectroscopy, and functional group substitution experiments suggest that the GU closing pairs form and that sheared GA pairs form in the internal loop. The NMR spectra show evidence of conformational dynamics, and several GA pairings are possible. Thus, this dynamic loop presents several possible structures for potential binding of drugs that target group I self-splicing introns. The results also contribute to understanding the structural and dynamic basis for the function and thermodynamic stability of this loop.     

Thermodynamic stabilities of internal loops with GU closing pairs in RNA

Many internal loops that form tertiary contacts in natural RNAs have GU closing pairs; examples include the tetraloop receptor and P1 helix docking site in group I introns. Thus, thermodynamic parameters of internal loops with GU closing pairs can contribute to the prediction of both secondary and tertiary structure. Oligoribonucleotide duplexes containing small internal loops with GU closing pairs were studied by optical melting, one-dimensional imino proton NMR, and one-dimensional phosphorus NMR. The thermodynamic stabilities of asymmetric internal loops with GU closing pairs relative to those of loops with GC closing pairs may be explained by hydrogen bonds. In contrast, the free energy increments for symmetric internal loops of two noncanonical pairs with GU closing pairs relative to loops with GC closing pairs show much more sequence dependence. Imino proton and phosphorus NMR spectra suggest that some GA pairs adjacent to GU closing pairs may form an overall thermodynamically stable but non-A-form conformation.     

Factors affecting thermodynamic stabilities of RNA 3 x 3 internal loops

Internal loops in RNA are important for folding and function. The 3 x 3 nucleotide internal loops are the smallest size symmetric loops with a potential noncanonical base pair (middle pair) flanked on both sides by a noncanonical base pair (loop-terminal pair). Thermodynamic and structural insights acquired for 3 x 3 loops should improve approximations for stabilities of 3 x 3 and larger internal loops. Most natural 3 x 3 internal loops are purine rich, which is also true of other internal loops. A series of oligoribonucleotides containing different 3 x 3 internal loops were studied by UV melting and imino proton NMR. Both loop-terminal and middle pairs contribute to the thermodynamic stabilities of 3 x 3 loops. Extra stabilization of -1.2 kcal/mol was found for a GA middle pair when flanked by at least one non-pyrimidine-pyrimidine loop-terminal pair. A penalty of approximately 1 kcal/mol was found for loops with a single loop-terminal GA pair that has a U 3' to the G of the GA pair. A revised model for predicting stabilities of 3 x 3 loops is derived by multiple linear regression.     

Stabilities of consecutive A.C, C.C, G.G, U.C, and U.U mismatches in RNA internal loops: Evidence for stable hydrogen-bonded U.U and C.C.+ pairs.

The stability and structure of RNA duplexes with consecutive A.C, C.A, C.C, G.G, U.C, C.U, and U.U mismatches were studied by UV melting, CD, and NMR. The results are compared to previous results for GA and AA internal loops [SantaLucia, J., Kierzek, R., & Turner, D. H. (1990) Biochemistry 29, 8813-8819; Peritz, A., Kierzek, R., & Turner, D.H. (1991) Biochemistry 30, 6428-6436)]. The observed order for stability increments of internal loop formation at pH 7 is AG = GA approximately UU greater than GG greater than or equal to CA greater than or equal to AA = CU = UC greater than or equal to CC greater than or equal to AC. The results suggest two classes for internal loops with consecutive mismatches: (1) loops that stabilize duplexes and have strong hydrogen bonding and (2) loops that destabilize duplexes and may not have strong hydrogen bonding. Surprisingly, rCGCUUGCG forms a very stable duplex at pH 7 in 1 M NaCl with a TM of 44.8 degrees C at 1 x 10(-4) M and a delta G degrees 37 of -7.2 kcal/mol. NOE studies of the imino protons indicate hydrogen bonding within the U.U mismatches in a wobble-type structure. Resonances corresponding to the hydrogen-bonded uridines are located at 11.3 and 10.4 ppm. At neutral pH, rCGCCCGCG is one of the least stable duplexes with a TM of 33.2 degrees C and delta G degrees 37 of -5.1 kcal/mol. Upon lowering the pH to 5.5, however, the TM increases by 12 degrees C, and delta G degrees 37 becomes more favorable by 2.5 kcal/mol. The pH dependence of rCGCCCGCG may be due to protonation of the internal loop C's, since no changes in thermodynamic parameters are observed for rCGCUUGCG between pH 7 and 5.5. Furthermore, two broad imino proton resonances are observed at 10.85 and 10.05 ppm for rCGCCCGCG at pH 5.3, but not at pH 6.5. This is also consistent with C.C+ base pairs forming at pH 5.5. rCGCCAGCG and rGGCACGCC have a small pH dependence, with TM increases of 5 and 3 degrees C, respectively, upon lowering the pH from 7 to 5.5. rCGCCUGCG and rCGCUCGCG also show little pH dependence, with TM increases of 0.8 and 1.4 degrees C, respectively, upon lowering the pH to 5.5.(ABSTRACT TRUNCATED AT 400 WORDS)     

NMR structures of (rGCUGAGGCU)2 and (rGCGGAUGCU)2: probing the structural features that shape the thermodynamic stability of GA pairs

The NMR structures of [see text] and [see text] are reported. The internal loop, [see text], is about 2 kcal/mol more stable than [see text] at 37 degrees C. The duplexes assemble into similar global folds characterized by the formation of tandem sheared GA pairs. The different stabilities of the loops are accompanied by differences in the local structure of the closing GU pairs. In the [see text] internal loop, the GU pairs form canonical wobble configurations with two hydrogen bonds, whereas in [see text], the GU pairs form a single hydrogen bond involving the amino group, GH22, and the carbonyl group, UO4. This pairing is similar to the GU closing pair of the 690 hairpin loop found in E. coli 16S rRNA. The [see text] and [see text] structures reveal how the subtle interplay between stacking and hydrogen bonding determines sequence dependent conformation and thermodynamic stability. Thus, this work provides structural and thermodynamic benchmarks for theoreticians in the ongoing effort to understand the sequence dependence of RNA physicochemical properties.     

Molecular recognition of 6'-N-5-hexynoate kanamycin A and RNA 1x1 internal loops containing CA mismatches

In our previous study to identify the RNA internal loops that bind an aminoglycoside derivative, we determined that 6'-N-5-hexynoate kanamycin A prefers to bind 1x1 nucleotide internal loops containing C·A mismatches. In this present study, the molecular recognition between a variety of RNAs that are mutated around the C·A loop and the ligand was investigated. Studies show that both loop nucleotides and loop closing pairs affect binding affinity. Most interestingly, it was shown that there is a correlation between the thermodynamic stability of the C·A internal loops and ligand affinity. Specifically, C·A loops that had relatively high or low stability bound the ligand most weakly whereas loops with intermediate stability bound the ligand most tightly. In contrast, there is no correlation between the likelihood that a loop forms a C-A(+) pair at lower pH and ligand affinity. It was also found that a 1x1 nucleotide C·A loop that bound to the ligand with the highest affinity is identical to the consensus site in RNAs that are edited by adenosine deaminases acting on RNA type 2 (ADAR2). These studies provide a detailed investigation of factors affecting small molecule recognition of internal loops containing C·A mismatches, which are present in a variety of RNAs that cause disease.     

Stability of RNA hairpin loops closed by AU base pairs

Thermodynamic parameters are reported for hairpin formation in 1 M NaCl by RNA sequence of the type GCAXUAAUYUGC, where XY is the set of 10 possible mismatch base pairs. A nearest-neighbor analysis of the data indicates that the free energy of loop formation at 37 degrees C varies from 3.2 to 5.0 kcal/mol. These results combined with the model previously developed [Dale et al. (2000) RNA 6, 608] allow improvements in the model to predict the stability of RNA hairpin loops: DeltaG degrees (37L(n) = DeltaG degrees (37i(n)) + DeltaG degrees (37MM) - 0.8 (if first mismatch is GA or UU) - 0.8 (if first mismatch is GG and loop is closed on 5' side by a purine). Here, DeltaG degrees (37i(n) is the free energy for initiating a loop of n nucleotides, and DeltaG degrees (37MM) is the free energy for the interaction of the first mismatch with the closing base pair. Hairpins with GG first mismatches were found to vary in stability depending upon the orientation of the closing base pair (5' or 3' purine relative to the loop). The model gives good agreement when tested against four naturally occurring hairpin sequences.     

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.     

A test of the model to predict unusually stable RNA hairpin loop stability

To investigate the accuracy of a model [Giese et al., 1998, Biochemistry37:1094-1100 and Mathews et al., 1999, JMol Biol 288:911-940] that predicts the stability of RNA hairpin loops, optical melting studies were conducted on sets of hairpins previously determined to have unusually stable thermodynamic parameters. Included were the tetraloops GNRA and UNCG (where N is any nucleotide and R is a purine), hexaloops with UU first mismatches, and the hairpin loop of the iron responsive element, CAGUGC. The experimental values for the GNRA loops are in excellent agreement (deltaG degrees 37 within 0.2 kcal/mol and melting temperature (TM) within 4 degrees C) with the values predicted by the model. When the UNCG hairpin loops are treated as tetraloops, and a bonus of 0.8 kcal/mol included in the prediction to account for the extra stable first mismatch (UG), the measured and predicted values are also in good agreement (deltaG degrees 37 within 0.7 kcal/mol and TM within 3 degrees C). Six hairpins with unusually stable UU first mismatches also gave good agreement with the predictions (deltaG degrees 37 within 0.5 kcal/mol and TM within 8 degrees C), except for hairpins closed by wobble base pairs. For these hairpins, exclusion of the additional stabilization term for UU first mismatches improved the prediction (AG degrees 37 within 0.1 kcal/mol and TM within 3 degrees C). Hairpins with the iron-responsive element loop were not predicted well by the model, as measured deltaG degrees 37 values were at least 1 kcal/mol greater than predicted.     

Sequence dependence of the stability of RNA hairpin molecules with six nucleotide loops

Thermodynamic parameters are reported for hairpin formation in 1 M NaCl by RNA sequence of the types GCGXUAAUYCGC and GGUXUAAUYACC with Watson-Crick loop closure, where XY is the set of 10 possible mismatch base pairs. A nearest-neighbor analysis of the data indicates the free energy of loop formation at 37 degrees C varies from 3.1 to 5.1 kcal/mol. These results agree with the model previously developed [Vecenie, C. J., and Serra, M. J. (2004) Biochemistry 43, 11813] to predict the stability of RNA hairpin loops: DeltaG degrees (37L(n) = DeltaG degrees (37i(n) + DeltaG degrees (37MM) - 0.8 (if first mismatch is GA or UU) - 0.8 (if first mismatch is GG and loop is closed on the 5' side by a purine). Here, DeltaG degrees (37i(n) is the free energy for initiating a loop of n nucleotides, and DeltaG degrees (37MM) is the free energy for the interaction of the first mismatch with the closing base pair. Thermodynamic parameters are also reported for hairpin formation in 1 M NaCl by RNA sequence of the types GACGXUAAUYUGUC and GGUXUAAUYGCC with GU base pair closure, where XY is the set of 10 possible mismatch base pairs. A nearest-neighbor analysis of the data indicates the free energy of loop formation at 37 degrees C varies from 3.6 to 5.3 kcal/mol. These results allow the development of a model for predicting the stability of hairpin loops closed by GU base pairs. DeltaG degrees (37L(n) (kcal/mol) = DeltaG degrees (37i(n) - 0.8 (if the first mismatch is GA) - 0.8 (if the first mismatch is GG and the loop is closed on the 5' side by a purine). Note that for these hairpins, the stability of the loops does not depend on DeltaG degrees (37MM). For hairpin loops closed by GU base pairs, the DeltaG degrees (37i(n) values, when n = 4, 5, 6, 7, and 8, are 4.9, 5.0, 4.6, 5.0, and 4.8 kcal/mol, respectively. The model gives good agreement when tested against six naturally occurring hairpin sequences. Thermodynamic values for terminal mismatches adjacent to GC, GU, and UG base pairs are also reported.     

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.     

Thermodynamics of single mismatches in RNA duplexes

The thermodynamic properties and structures of single mismatches in short RNA duplexes were studied in optical melting and imino proton NMR experiments. The free energy increments at 37 degrees C measured for non-GU single mismatches range from -2.6 to 1.7 kcal/mol. These increments depend on the identity of the mismatch, adjacent base pairs, and the position in the helix. UU and AA mismatches are more stable close to a helix end, but GG mismatch stability is essentially unaffected by the position in the helix. Approximations are suggested for predicting stabilities of single mismatches in short RNA duplexes.     

NMR reveals the absence of hydrogen bonding in adjacent UU and AG mismatches in an isolated internal loop from ribosomal RNA

NMR studies provide insights into structural features of internal loops. These insights can be combined with thermodynamic studies to generate models for predicting structure and energetics. The tandem mismatch internal loop, 5'GUGG3'(3'CUAC5'), has been studied by NMR. The NMR structure reveals an internal loop with no hydrogen bonding between the loop bases and with the G in the AG mismatch flipped out of the helix. The sequence of this internal loop is highly conserved in rRNA. The loop is located in the large ribosomal subunit and is part of a conserved 58-nt fragment that is the binding domain of ribosomal protein L11. Structural comparisons between variants of this internal loop in crystal structures of the 58-nt domain complexed with L11 protein and of the large ribosomal subunit (LSU) suggest that this thermodynamically destabilizing internal loop is partially preorganized for tertiary interactions and for binding L11. A model for predicting the base pairing and free energy of 2 x 2 nucleotide internal loops with a purine-purine mismatch next to a pyrimidine-pyrimidine mismatch is proposed on the basis of the present NMR structure and previously reported thermodynamics.     

Consecutive GA pairs stabilize medium-size RNA internal loops

Internal loops in RNA are important for folding and function. Consecutive noncanonical pairs can form in internal loops having at least two nucleotides on each side. Thermodynamic and structural insights into such internal loops should improve approximations for their stabilities and predictions of secondary and three-dimensional structures. Most natural internal loops are purine rich. A series of oligoribonucleotides that form purine-rich internal loops of 5-10 nucleotides, including kink-turn loops, were studied by UV melting, exchangeable proton and phosphorus NMR. Three consecutive GA pairs with the motif 5' Y GGA/3' R AAG or GGA R 3'/AAG Y 5' (i.e., 5' GGA 3'/3' AAG 5' closed on at least one side with a CG, UA, or UG pair with Y representing C or U and R representing A or G) stabilize internal loops having 6-10 nucleotides. Certain motifs with two consecutive GA pairs are also stabilizing. In internal loops with three or more nucleotides on each side, the motif 5' U G/3' G A has stability similar to 5' C G/3' G A. A revised model for predicting stabilities of internal loops with 6-10 nucleotides is derived by multiple linear regression. Loops with 2 x 3 nucleotides are predicted well by a previous thermodynamic model.     

Thermodynamic parameters for loop formation in RNA and DNA hairpin tetraloops.

We determined the melting temperatures (Tm) and thermodynamic parameters of 15 RNA and 19 DNA hairpins at 1 M NaCl, 0.01 M sodium phosphate, 0.1 mM EDTA, at pH 7. All these hairpins have loops of four bases, the most common loop size in 16S and 23S ribosomal RNAs. The RNA hairpins varied in loop sequence, loop-closing base pair (A.U, C.G, or G.C), base sequence of the stem, and stem size (four or five base pairs). The DNA hairpins varied in loop sequence, loop-closing base pair (C.G, or G.C), and base sequence of the four base-pair stem. Thermodynamic properties of a hairpin may be represented by nearest-neighbor interactions of the stem plus contributions from the loop. Thus, we obtained thermodynamic parameters for the formation of RNA and DNA tetraloops. For the tetraloops we studied, a free energy of loop formation (at 37 degrees C) of about +3 kcal/mol is most common for either RNA or DNA. There are extra stable loops with delta G degrees 37 near +1 kcal/mol, but the sequences are not necessarily the same for RNA and DNA. The closing base pair is also important; changing from C.G to G.C lowered the stability of several tetraloops in both RNA and DNA. These values will be useful in predicting RNA and DNA secondary structures.     

Thermodynamics of RNA duplexes with tandem mismatches containing a uracil-uracil pair flanked by C.G/G.C or G.C/A.U closing base pairs

The thermodynamics governing the denaturation of RNA duplexes containing 8 bp and a central tandem mismatch or 10 bp were evaluated using UV absorbance melting curves. Each of the eight tandem mismatches that were examined had one U-U pair adjacent to another noncanonical base pair. They were examined in two different RNA duplex environments, one with the tandem mismatch closed by G.C base pairs and the other with G.C and A.U closing base pairs. The free energy increments (Delta Gdegrees(loop)) of the 2 x 2 loops were positive, and showed relatively small differences between the two closing base pair environments. Assuming temperature-independent enthalpy changes for the transitions, (Delta Gdegrees(loop)) for the 2 x 2 loops varied from 0.9 to 1.9 kcal/mol in 1 M Na(+) at 37 degrees C. Most values were within 0.8 kcal/mol of previously estimated values; however, a few sequences differed by 1.2-2.0 kcal/mol. Single strands employed to form the RNA duplexes exhibited small noncooperative absorbance increases with temperature or transitions indicative of partial self-complementary duplexes. One strand formed a partial self-complementary duplex that was more stable than the tandem mismatch duplexes it formed. Transitions of the RNA duplexes were analyzed using equations that included the coupled equilibrium of self-complementary duplex and non-self-complementary duplex denaturation. The average heat capacity change (DeltaC(p)) associated with the transitions of two RNA duplexes was estimated by plotting DeltaH degrees and DeltaS degrees evaluated at different strand concentrations as a function of T(m) and ln T(m), respectively. The average DeltaC(p) was 70 +/- 5 cal K(-)(1) (mol of base pairs)(-)(1). Consideration of this heat capacity change reduced the free energy of formation at 37 degrees C of the 10 bp control RNA duplexes by 0.3-0.6 kcal/mol, which may increase Delta Gdegrees(loop) values by similar amounts.     

Conformational determinants of tandem GU mismatches in RNA: insights from molecular dynamics simulations and quantum mechanical calculations

Structure and energetic properties of base pair mismatches in duplex RNA have been the focus of numerous investigations due to their role in many important biological functions. Such efforts have contributed to the development of models for secondary structure prediction of RNA, including the nearest-neighbor model. In RNA duplexes containing GU mismatches, 5'-GU-3' tandem mismatches have a different thermodynamic stability than 5'-UG-3' mismatches. In addition, 5'-GU-3' mismatches in some sequence contexts do not follow the nearest-neighbor model for stability. To characterize the underlying atomic forces that determine the structural and thermodynamic properties of GU tandem mismatches, molecular dynamics (MD) simulations were performed on a series of 5'-GU-3' and 5'-UG-3' duplexes in different sequence contexts. Overall, the MD-derived structural models agree well with experimental data, including local deviations in base step helicoidal parameters in the region of the GU mismatches and the model where duplex stability is associated with the pattern of GU hydrogen bonding. Further analysis of the simulations, validated by data from quantum mechanical calculations, suggests that the experimentally observed differences in thermodynamic stability are dominated by GG interstrand followed by GU intrastrand base stacking interactions that dictate the one versus two hydrogen bonding scenarios for the GU pairs. In addition, the inability of 5'-GU-3' mismatches in different sequence contexts to all fit into the nearest-neighbor model is indicated to be associated with interactions of the central four base pairs with the surrounding base pairs. The results emphasize the role of GG and GU stacking interactions on the structure and thermodynamics of GU mismatches in RNA.     

Isolation and characterization of thermodynamically stable and unstable RNA hairpins from a triloop combinatorial library

Hairpins are the most common elements of RNA secondary structure, playing important roles in RNA tertiary architecture and forming protein binding sites.Triloops are common in a variety of naturally occurring RNA hairpins, but little is known about their thermodynamic stability. Reported here are the sequences and thermodynamic parameters for a variety of stable and unstable triloop hairpins. Temperature gradient gel electrophoresis (TGGE) can be used to separate a simple RNA combinatorial library based on thermal stability [Bevilacqua, J. M., and Bevilacqua, P. C. (1998) Biochemistry 45, 15877-15884]. Here we introduce the application of TGGE to separating and analyzing a complex RNA combinatorial library based on thermal stability, using an RNA triloop library. Several rounds of in vitro selection of an RNA triloop library were carried out using TGGE, and preferences for exceptionally stable and unstable closing base pairs and loop sequences were identified. For stable hairpins, the most common closing base pair is CG, and U-rich loop sequences are preferred. Closing base pairs of GC and UA result in moderately stable hairpins when combined with a stable loop sequence. For unstable hairpins, the most common closing base pairs are AU and UG, and U-rich loop sequences are no longer preferred. In general, the contributions of the closing base pair and loop sequence to overall hairpin stability appear to be additive. Thermodynamic parameters for individual hairpins determined by UV melting are generally consistent with outcomes from selection experiments, with hairpins containing a CG closing base pair having a DeltaDeltaG degrees (37) 2.1-2.5 kcal/mol more favorable than hairpins with other closing base pairs. Sequences and thermodynamic rules for triloop hairpins should aid in RNA structure prediction and determination of whether naturally occurring triloop hairpins are thermodynamically stable.