Free energy of imperfect nucleic acid helices: II. Small Hairpin Loops

Abstract

Physical studies of enzymically synthesized oligonucleotides of defined sequence are used to evaluate quantitatively the stability of small RNA hairpin loops and helices. The series (Ap) 4G(pC) N (pU) 4, N = 4, 5 or 6, exists as monomolecular hairpin helices when N ≥ 5, and as imperfect dimer helices when N ≤ 4. In this size range, hairpin loops become more favorable (less destabilizing thermodynamically) as they increase in size from 3 to 4 to 5 unbonded nucleotides. Very small hairpin loops are particularly destabilizing; molecules whose base sequence would imply a hairpin loop of three nucleotides will generally exist with a loop of five, including a broken terminal base pair. Thermodynamic parameters for base pair and loop formation are calculated by a method which makes unnecessary the use of measured enthalpies of polynucleotide melting. Literature data on oligonucleotide double helices yield estimates of the free energy contribution from each of the six types of stacking interactions between three possible neighboring base pairs. The advantage of this approach is that the properties of oligonucleotides are used in predicting the stability of small RNA helices, avoiding the long extrapolation from the properties of high polymers. We provide Tables of temperature-dependent free energies that allow one to predict the stability and thermal transition temperature of many simple RNA secondary structures (applicable to ~1 m-Na + concentration). As an example, we apply the rules to an isolated fragment of tRNA Ser (yeast) (Coutts, 1971), whose properties were not used in calculating the free-energy parameters. The experimental melting temperature of 88 °C is predicted with an error margin of 5 deg. C.