Thermodynamic characterization of naturally occurring RNA tetraloops

Although tetraloops are one of the most frequently occurring secondary structure motifs in RNA, less than one-third of the 30 most frequently occurring RNA tetraloops have been thermodynamically characterized. Therefore, 24 stem–loop sequences containing common tetraloops were optically melted, and the thermodynamic parameters ΔH°, ΔS°, ΔG°_37, and T_M for each stem–loop were determined. These new experimental values, on average, are 0.7 kcal/mol different from the values predicted for these tetraloops using the model proposed by Vecenie CJ, Morrow CV, Zyra A, Serra MJ. 2006. Biochemistry 45: 1400–1407. The data for the 24 tetraloops reported here were then combined with the data for 28 tetraloops that were published previously. A new model, independent of terminal mismatch data, was derived to predict the free energy contribution of previously unmeasured tetraloops. The average absolute difference between the measured values and the values predicted using this proposed model is 0.4 kcal/mol. This new experimental data and updated predictive model allow for more accurate calculations of the free energy of RNA stem–loops containing tetraloops and, furthermore, should allow for improved prediction of secondary structure from sequence. It was also shown that tetraloops within the sequence 5′-GCCNNNNGGC-3′ are, on average, 0.6 kcal/mol more stable than the same tetraloop within the sequence 5′-GGCNNNNGCC-3′. More systemic studies are required to determine the full extent of non-nearest-neighbor effects on tetraloop stability.     

Unusual nucleotide conformations in GNRA and UNCG type tetraloop hairpins: evidence from Raman markers assignments.

High resolution NMR data on UNCG and GNRA tetraloops (where N is any of the four nucleotides and R is a purine) have shown that they contain ribonucleosides with unusual 2'-endo/anti and 3'-endo/syn conformations, in addition to the 3'-endo/anti ones which are regularly encountered in RNA chains. In the current study, Raman spectroscopy has been used to probe these nucleoside conformations and follow the order (hairpin) to disorder (random chain) structural transitions in aqueous phase in the 5-80 degreesC temperature range. Spectral evolution of GCAA and GAAA tetraloops, as formed in very short hairpins with only three G.C base pairs in their stems (T m >60 degreesC), are reported and compared with those previously published on UUCG and UACG tetraloops, for which the syn orientation of the terminal guanine as well as the 2'-endo/anti conformation of the third rC residue have been confirmed by means of vibrational marker bands. Raman data obtained as a function of temperature show that the first uracil in the UUCG tetraloop is stacked and the two middle residues (rU and rC) are in the 2'-endo/anti conformation, in agreement with the previously published NMR results. As far as the new data concerning the GNRA type tetraloops are concerned, they lead us to conclude that: (i) in both cases (GCAA and GAAA tetraloops) the adenine bases are stacked; (ii) the second rC residue in the GCAA tetraloop has a 3'-endo/anti conformation; (iii) the sugar pucker associated with the third rA residue in both tetraloops possibly undergoes a 3'-endo/2'-endo interconversion as predicted by NMR results; (iv) the stem adopts a regular A-form structure; (v) all other nucleosides of these two GNRA tetraloops possess the usual 3'-endo/anti conformation.     

Thermal stability, structural features, and B-to-Z transition in DNA tetraloop hairpins as determined by optical spectroscopy in d(CG)(3)T(4)(CG)(3) and d(CG)(3)A(4)(CG)(3) oligodeoxynucleotides

NMR and CD data have previously shown the formation of the T(4) tetraloop hairpin in aqueous solutions, as well as the possibility of the B-to-Z transition in its stem in high salt concentration conditions. It has been shown that the stem B-to-Z transition in T(4) hairpins leads to S (south)- to N (north)-type conformational changes in the loop sugars, as well as anti to syn orientations in the loop bases. In this article, we have compared by means of UV absorption, CD, Raman, and Fourier transform infrared (FTIR), the thermodynamic and structural properties of the T(4) and A(4) tetraloop hairpins formed in 5'-d(CGCGCG-TTTT-CGCGCG)-3' and 5'-d(CGCGCG-AAAA-CGCGCG)-3', respectively. In presence of 5M NaClO(4), a complete B-to-Z transition of the stems is first proved by CD spectra. UV melting profiles are consistent with a higher thermal stability of the T(4) hairpin compared to the A(4) hairpin. Order-to-disorder transition of both hairpins has also been analyzed by means of Raman spectra recorded as a function of temperature. A clear Z-to-B transition of the stem has been confirmed in the T(4) hairpin, and not in the A(4) hairpin. With a right-handed stem, Raman and FTIR spectra have confirmed the C2'-endo/anti conformation for all the T(4) loop nucleosides. With a left-handed stem, a part of the T(4) loop sugars adopt a N-type (C3'-endo) conformation, and the C3'-endo/syn conformation seems to be the preferred one for the dA residues involved in the A(4) tetraloop.     

Common Structural Features of UUCG and UACG Tetraloops in Very Short Hairpins Determined by UV Absorption,Raman, IR and NMR Spectroscopies

Abstract

Thermodynamic and structural properties of two UNCG tetraloops in very short hairpin octamers, 5′-r (GCUUCGGC)-3′ and 5′-r (GCUACGGC)-3′. have been studied by means of various physical techniques. Melting profiles of both octamers, obtained from UV absorption spectra taken as a function of temperature, are consistent with a monophasic, progressive and completely reversible order-to-disorder transition and confirm their unusual structural stability (T_m>51^° C). The ¹H, ¹³C and ³¹P NMR chemical shifts and coupling constants of the UACG loop nucleotides are comparable with those reported previously for UUCG loops, i.e. 2′-endo/anti conformation of the second and third nucleotide of the loop as well as the syn orientation of the ultimate guanine base and the A-type double helical conformation of the hairpin stem. Simulation of quantitative NOESY volumes shows that the UACG octamer adopts a very rigid compact structure which is well represented by an average order parameter of 0.9. Three base-pairs and four additional strong hydrogen bonds are undoubtedly responsible for such limited flexibility. Raman and infrared spectra as a function of temperature reflect the order-to-disorder transition, as well. Vibrational conformational markers in low temperature spectra of both octamers indicate the hairpin structure as the major conformer in aqueous phase. These spectra further support the structural features of most of the nucleotides involved in the tetraloops and clearly demonstrate the structural similarities of the phosphodiester backbone in both hairpins. Consequently, on the basis of all present results, one can deduce that the conformational features of the UUCG and UACG tetraloops seem to be inherent to the UNCG type tetraloops, regardless of either the nature of the tetraloop second base or the stem length.     

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.     

Increase in hybridization rates with oligodeoxyribonucleotides containing locked nucleic acids

Locked nucleic acids (LNAs) incorporated into either stable single stranded oligonucleotides containing tetraloops or their complements have been found to increase second order hybridization rate constants by an order of magnitude compared to the all-DNA hybridization rate constants. Model sequences composed of 20 bases in length that can form hairpins due to a stable GAAA tetraloop were used where LNAs were substituted for the nucleotides in the loop, stem, or end regions of the strand and in the complementary strand. Substitution of the LNAs to the loop predictably raised the melting temperatures of the duplex however, the hybridization rates between the tetraloop and the complementary sequence also increased. In contrast, when LNAs were substituted in the stem, the hybridization rate decreased implying the formation of a more stable hairpin. Substitution of LNAs into the end region of the sequence had little effect on the hybridization rate constants although melting temperatures still showed a predictable increase. Rates also increased when LNAs were substituted into complementary strands of DNA tetraloops. The increase in hybridization rate constant is being attributed to changes in the structure of the stable single strands.     

Comprehensive features of natural and in vitro selected GNRA tetraloop-binding receptors

Specific recognitions of GNRA tetraloops by small helical receptors are among the most widespread long-range packing interactions in large ribozymes. However, in contrast to GYRA and GAAA tetraloops, very few GNRA/receptor interactions have yet been identified to involve GGAA tetraloops in nature. A novel in vitro selection scheme based on a rigid self-assembling tectoRNA scaffold designed for isolation of intermolecular interactions with A-minor motifs has yielded new GGAA tetraloop-binding receptors with affinity in the nanomolar range. One of the selected receptors is a novel 12 nt RNA motif, (CCUGUG ... AUCUGG), that recognizes GGAA tetraloop hairpin with a remarkable specificity and affinity. Its physical and chemical characteristics are comparable to those of the well-studied '11nt' GAAA tetraloop receptor motif. A second less specific motif (CCCAGCCC ... GAUAGGG) binds GGRA tetraloops and appears to be related to group IC3 tetraloop receptors. Mutational, thermodynamic and comparative structural analysis suggests that natural and in vitro selected GNRA receptors can essentially be grouped in two major classes of GNRA binders. New insights about the evolution, recognition and structural modularity of GNRA and A-minor RNA-RNA interactions are proposed.     

Isolation and characterization of a family of stable RNA tetraloops with the motif YNMG that participate in tertiary interactions

RNA is known to fold into a variety of structural elements, many of which have sufficient sequence complexity to make the thermodynamic study of each possible variant impractical. We previously reported a method for isolating stable and unstable RNA sequences from combinatorial libraries using temperature gradient gel electrophoresis (TGGE). This method was used herein to analyze a six-nucleotide RNA hairpin loop library. Three rounds of in vitro selection were performed using TGGE, and unusually stable RNAs were identified by cloning and sequencing. Known stable tetraloops were found, including sequences belonging to the UNCG motif closed by a CG base pair, and the CUUG motif closed by a GC base pair. In addition, unknown tetraloops were found that were nearly as stable as cUNCGg, including sequences related through substitution of the U with a C (Y), the C with an A (M), or both. These substitutions allow hydrogen bonding and stacking interactions in the UNCG loop to be maintained. Thermodynamic analysis of YNMG and variant loops confirmed optimal stability with Y at position 1 and M at position 3. Similarity in structure and stability among YNMG loops was further supported by deoxyribose substitution, CD, and NMR experiments. A conserved tertiary interaction in 16S rRNA exists between a YAMG loop at position 343 and two adenines in the loop at position 159 (Escherichia coli numbering). NMR and functional group substitution experiments suggest that YNAG loops in particular have enhanced flexibility, which allows the tertiary interaction to be maintained with diverse loop sequences at position 159. Taken together, these results support the existence of an extended family of UNCG-like tetraloops with the motif cYNMGg that are thermodynamically stable and structurally similar and can engage in tertiary interactions in large RNA molecules.     

Electrical detection of the temperature induced melting transition of a DNA hairpin covalently attached to gold interdigitated microelectrodes

The temperature induced melting transition of a self-complementary DNA strand covalently attached at the 5′ end to the surface of a gold interdigitated microelectrode (GIME) was monitored in a novel, label-free, manner. The structural state of the hairpin was assessed by measuring four different electronic properties of the GIME (capacitance, impedance, dissipation factor and phase angle) as a function of temperature from 25°C to 80°C. Consistent changes in all four electronic properties of the GIME were observed over this temperature range, and attributed to the transition of the attached single-stranded DNA (ssDNA) from an intramolecular, folded hairpin structure to a melted ssDNA. The melting curve of the self-complementary single strand was also measured in solution using differential scanning calorimetry (DSC) and UV absorbance spectroscopy. Temperature dependent electronic measurements on the surface and absorbance versus temperature values measured in solution experiments were analyzed assuming a two-state process. The model analysis provided estimates of the thermodynamic transition parameters of the hairpin on the surface. Two-state analyses of optical melting data and DSC measurements provided evaluations of the thermodynamic transition parameters of the hairpin in solution. Comparison of surface and solution measurements provided quantitative evaluation of the effect of the surface on the thermodynamics of the melting transition of the DNA hairpin.     

Pronounced instability of tandem IU base pairs in RNA

Optical melting was used to determine the stabilities of three series of RNA oligomers containing tandem XU base pairs, GGCXUGCC (5′XU3′), GGCUXGCC (5′UX3′) and GGCXXGGC/CCGUUCCG (5′XX3′), where X is either A, G or I (inosine). The helices containing tandem AU base pairs were the most stable in the first two series (5′XU3′ and 5′UX3′), with an average melting temperature ~11°C higher than the helices with tandem 5′GU3′ base pairs and 25°C higher than the helices with tandem 5′IU3′ base pairs. For the third series (5′XX3′), the helix containing tandem GG is the most stable, with an average melting temperature ~2°C higher than the helix with tandem AA base pairs and ~24°C higher than the helix with tandem II base pairs. The thermodynamic stability of the oligomers with tandem IU base pairs was also investigated as a function of magnesium ion concentration. As with normal A–U or G–U tandem duplexes, the data could best be interpreted as non-specific binding of magnesium ions to the inosine-containing RNA oligonucleotides.     

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.     

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

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

The GANC tetraloop: a novel motif in the group IIC intron structure

Tetraloops are a common building block for RNA tertiary structure, and most tetraloops fall into one of three well-characterized classes: GNRA, UNCG, and CUYG. Here, we present the sequence and structure of a fourth highly conserved class of tetraloop that occurs only within the ζ-ζ′ interaction of group IIC introns. This GANC tetraloop was identified, along with an unusual cognate receptor, in the crystal structure of the group IIC intron and through phylogenetic analysis of intron RNA sequence alignments. Unlike conventional tetraloop-receptor interactions, which are stabilized by extensive hydrogen-bonding interactions, the GANC-receptor interaction is limited to a single base stack between the conserved adenosine of the tetraloop and a single purine of the receptor, which consists of a one- to three-nucleotide bulge and does not contain an A-platform. Unlike GNRA tetraloops, the GANC tetraloop forms a sharp angle relative to the adjacent helix, bending by approximately 45° toward the major groove side of the helix. These structural attributes allow GANC tetraloops to fit precisely within the group IIC intron core, thereby demonstrating that structural motifs can adapt to function in a specific niche.     

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.     

RNA tetraloops as minimalist substrates for aminoacylation.

Previous work established that seven-base-pair hairpin microhelices with sequences based on the acceptor stems of alanine, glycine, methionine, and histidine tRNAs can be aminoacylated specifically with their cognate amino acids. To obtain "minimalist" substrates with fewer base pairs, we took advantage of the high thermodynamic stability of RNA tetraloop motifs that are found in ribosomal RNAs. We show here that rationally designed RNA tetraloops with as few as four base pairs are substrates for aminoacylation. Major nucleotide determinants for recognition by the class II synthetases were incorporated into each of the respective tetraloop substrates, resulting in specific aminoacylation by the alanine, glycine, and histidine tRNA synthetases. An analysis of the kinetics of aminoacylation shows that, for the alanine system, the majority of the transition-state stabilization provided by the synthetase-tRNA interaction is reproduced by the interaction of the synthetase with nucleotides in its minimalist tetraloop substrate. In an extension of this work, we also observed specific aminoacylation with the class I methionine tRNA synthetase of RNA tetraloops based on sequences in the acceptor stem of methionine tRNA. Thus, the results demonstrate four different examples where specific aminoacylation is directed by sequences/structures contained in less than half of a turn of an RNA helix.     

稳定的RNA发夹结构及其生物学功能

以UNCG、GNRA、CUUG(N=A、U、C或G,R=G或A)为端环能够形成稳定的、保守的发夹结构。高分辨率的溶液结构、晶体结构和计算机模拟等方法从原子水平上解析了这些发夹特殊的结构特征。在体内,它们发挥着重要的生物学功能:在折叠过程中作为折叠的起始位置帮助组织RNA分子正确折叠;与核酸受体结合参与三级相互作用;与蛋白质发生相互作用;阻止逆转录酶的延伸等等。另外,由于C(UUCG)G发夹极其稳定的特征,在体外RNA分子的实验测定中它还是稳定核酸结构的理想工具。这些稳定的发夹广泛分布于体内rRNA、催化RNA和非编码mRNA中。但在对人类编码区mRNA结构特征的研究当中,却未发现C(UUCG)G发夹。     

The crystal structure of UUCG tetraloop

All large structured RNAs contain hairpin motifs made of a stem closed by several looped nucleotides. The most frequent loop motif is the UUCG one. This motif belongs to the tetraloop family and has the peculiarity of being highly thermodynamically stable. Here, we report the first crystal structure of two UUCG tetraloops embedded in a larger RNA-protein complex solved at 2.8 A resolution. The two loops present in the asymmetric unit are in a different crystal packing environment but, nevertheless, have an identical conformation. The observed structure is globally close to that obtained in solution by nuclear magnetic resonance. However, subtle differences point to a more detailed picture of the role played by 2'-hydroxyl groups in stabilising this tetraloop.     

Common and distinctive features of GNRA tetraloops based on a GUAA tetraloop structure at 1.4 A resolution

GNRA tetraloops (N is A, C, G, or U; R is A or G) are basic building blocks of RNA structure that often interact with proteins or other RNA structural elements. Understanding sequence-dependent structural variation among different GNRA tetraloops is an important step toward elucidating the molecular basis of specific GNRA tetraloop recognition by proteins and RNAs. Details of the geometry and hydration of this motif have been based on high-resolution crystallographic structures of the GRRA subset of tetraloops; less is known about the GYRA subset (Y is C or U). We report here the structure of a GUAA tetraloop determined to 1.4 A resolution to better define these details and any distinctive features of GYRA tetraloops. The tetraloop is part of a 27-nt structure that mimics the universal sarcin/ricin loop from Escherichia coli 23S ribosomal RNA in which a GUAA tetraloop replaces the conserved GAGA tetraloop. The adenosines of the GUAA tetraloop form an intermolecular contact that is a commonplace RNA tertiary interaction called an A-minor motif. This is the first structure to reveal in great detail the geometry and hydration of a GUAA tetraloop and an A-minor motif. Comparison of tetraloop structures shows a common backbone geometry for each of the eight possible tetraloop sequences and suggests a common hydration. After backbone atom superposition, equivalent bases from different tetraloops unexpectedly depart from coplanarity by as much as 48 degrees. This variation displaces the functional groups of tetraloops implicated in protein and RNA binding, providing a recognition feature.     

A comparison of thermodynamic foldings with comparatively derived structures of 16S and 16S-like rRNAs.

To increase our understanding of the dynamics and complexities of the RNA folding process, and therewith to improve our ability to predict RNA secondary structure by computational means, we have examined the foldings of a large number of phylogenetically and structurally diverse 16S and 16S-like rRNAs and compared these results with their comparatively derived secondary structures. Our initial goals are to establish the range of prediction success for this class of rRNAs, and to begin comparing and contrasting the foldings of these RNAs. We focus here on structural features that are predicted with confidence as well as those that are poorly predicted. Whereas the large set of Archaeal and (eu)Bacterial 16S rRNAs all fold well (69% and 55% respectively), some as high as 80%, many Eucarya and mitochondrial 16S rRNAs are poorly predicted (approximately 30%), with a few of these predicted as low as 10-20%. In general, base pairs interacting over a short distance and, in particular, those closing hairpin loops, are predicted significantly better than long-range base pairs and those closing multistem loops and bulges. The prediction success of hairpin loops varies, however, with their size and context. Analysis of some of the RNAs that do not fold well suggests that the composition of some hairpin loops (e.g., tetraloops) and the higher frequency of noncanonical pairs in their comparatively derived structures might contribute to these lower success rates. Eucarya and mitochondrial rRNAs reveal further novel tetraloop motifs, URRG/A and CRRG, that interchange with known stable tetraloop in the procaryotes.     

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.