A model for the stabilities of RNA hairpins based on a study of the sequence dependence of stability for hairpins of six nucleotides.

Thermodynamic parameters are reported for hairpin formation in 1 M NaCl by RNA sequences of the type GGCXUAAUYGCC, where XY is the set of 10 possible mismatch base pairs. A nearest neighbor analysis of the data indicates the free energy for loop formation at 37 C varies from 2.9 to 4.5 kcal/mol. Thermodynamic parameters are also reported for hairpin formation by RNA sequences of the type GGXGUAAUAYCC (where XY are CG, GC, AU, UA, GU, and UG), with the common naturally occurring GA first mismatch (45% of small and large subunit rRNA loops of six). These results allow the development of a model to predict the stability of RNA hairpin loops. The model includes the size of the loop, the identity of the closing base pair, the free energy increment (delta G zero 37MM) for interaction of the closing base pair with the first mismatch, and an additional stabilization term for GA and UU first mismatches. delta G zero 37L(n) = delta G zero 37i(n) + delta G zero 37MM + 0.4 (if closed by AU or UA) -0.7 (if first mismatch is GA or UU). Here delta G zero 37i(n) is the free energy for initiating a loop of n nucleotides. delta G zero 37i(n) for n = 4-9 is 4.9, 4.4, 5.0, 5.0, 5.1, and 5.2 kcal/mol, respectively. The delta G zero 37MM is derived from measurements of model duplexes with terminal mismatches. The model gives good agreement when tested against four naturally occurring hairpin sequences.

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.