Some experimental results for the thermodynamics of RNA folding cannot be explained by simple pairwise hydrogen-bonding models. Such effects include the stabilities of isoguanosineisocytidine (iG-iC) base pairs and of various 2 × 2 nucleotide internal loops. Presumably, these results can be explained by base stacking effects, which can be partitioned into Coulombic and overlap effects. We review experimental measurements that provide benchmarks for testing the approximations and theories used for modeling nucleic acids. Quantitative agreement between experiment and theory will indicate understanding of the interactions determining RNA stability and structure.An understanding of the physical-chemical interactions underlying RNA folding would allow predictions of structure and perhaps function from sequence. The large number of atoms in an RNA molecule necessitates the use of approximate methods, rather than the rigorous equations of quantum mechanics. Many approaches are possible, including coarsegrained potentials (1), which use residue-centered force fields; molecular mechanics, which uses atom-centered force fields such as AMBER (2), CHARMM (3, 4), and GROMOS (5, 6); approximate quantum mechanics, which considers interactions of electrons and nuclei (7,8); and QM/MM, which combines quantum mechanics with molecular mechanics (9-16). Depending on the property to be predicted, the size of the RNA, and the domain of interest, different approaches will provide acceptable approximations.The secondary structure of RNA can sometimes be deduced by sequence comparison, which relies on the Watson-Crick rules for base pairing and the assumption that secondary structures are more conserved than sequences for function (17). Often, however, there are not enough sequences to determine a definitive secondary structure. For these cases, freeenergy minimization with a nearest neighbor model is the most popular method for predicting secondary structures (18-38). In this method, possible secondary structure motifs are assigned free-energy parameters, and these values are added to predict the total free energy of forming a secondary structure. The structure with the lowest free energy is assumed to dominate in solution. Rigorously, however, the concentrations of the various possible structures are predicted to be weighted by a Boltzmann factor, exp(−ΔG°/RT), where ΔG° is the free energy change for folding, R is the gas constant, 1.987 cal K −1 mol −1 , and T is the temperature in kelvins. Understanding intermolecular interactions such as hydrogen bonding, stacking, and so forth would allow accurate prediction of ΔG° and therefore secondary structure for an RNA. † This work was supported by NIH Grant GM22939 (D.H.T. Watson-Crick base pairs are the most common and most extensively studied motif in RNA structures. Configurations with the same base compositions but different permutations of base pairs generally have different free energies. For example, the duplexes (5′CGCG3′) 2 and (5′GGCC3′) 2 both have four GC base pairs, bu...