The thermal stability of RNA duplexes differing by a single base pair (bp) substitution or mismatch were investigated by temperature gradient gel electrophoresis (TGGE). All base pair substitutions and mismatches were examined at six sites, and limited changes were investigated at three other sites. DNA templates for in vitro transcription were generated by the polymerase chain reaction (PCR), Transcribed forward and reverse single stranded RNAs were annealed to form 345 bp duper RNA, Solution melting curves of selected RNAs were in good agreement with the predicted three step transitions. Parallel TGGE was used to determine the relative stabilities of the RNAs, and perpendicular TGGE was employed to obtain mobility transitions and midpoint transition temperatures (T-mu) of the RNAs' first melting domain, The gel solvent included formamide and urea, The T-mu values of the first melting domain were influenced by the identity of the base pair substitution or mismatch as well as by the site's neighboring base pairs, The difference In the transition temperatures (delta T-mu) between pairs of RNA ranged from 0 to 5 degrees C. delta T-mu values were used to determine free energy differences (delta Delta G). For RNA pairs distinguished by a base pair substitution, the delta Delta G values were closely correlated with free energy differences calculated from stacking free energies determined from melting studies in 1 M Na+ [Serra, M. J., and Turner, D. H. (1995) Methods Enzymol. 259, 242-261.] An algorithm was developed using the free energies of terminal mismatches [Serra, M, J., and Turner, D. H. (1995) Methods Enzymol. 259, 242-261] that provided very good agreement with experimental free energies for the single internal mismatches.