Thermodynamic determination of RNA duplex stability in magnesium solutions

Abstract

The prediction of RNA secondary structure and thermodynamics from sequence relies on free energy minimization and nearest neighbor parameters. Currently, algorithms used to make these predictions are based on parameters from optical melting studies performed in 1 M NaCl. However, many physiological and biochemical buffers containing RNA include much lower concentrations of monovalent cations and the presence of divalent cations. In order to improve these algorithms, thermodynamic data was previously collected for RNA duplexes in solutions containing 71, 121, 221, and 621 mM Na⁺. From this data, correction factors for free energy (ΔG^°₃₇) and melting temperature (T_m) were derived. Despite these newly derived correction factors for sodium, the stabilizing effects of magnesium have been ignored. Here, the same RNA duplexes were melted in solutions containing 0.5, 1.5, 3.0, and 10.0 mM Mg²⁺ in the absence of monovalent cations. Correction factors for T_m and ΔG^°₃₇ were derived to scale the current parameters to a range of magnesium concentrations. The T_m correction factor predicts the melting temperature within 1.2°C, and the ΔG^°₃₇ correction factor predicts the free energy within 0.30 kcalmol. These newly derived magnesium correction factors can be incorporated into algorithms that predict RNA secondary structure and stability from sequence.