Understanding the thermodynamic stability of an RNA hairpin and its mutant

The thermodynamic stability of RNA hairpin loops has been a subject of considerable interest in the recent past (Wimberly et al., 1991). There have been experimental reports indicating that the hairpins with a C(UUCG)G loop sequence are thermodynamically very stable (Wimberly et al., 1991). We used the solution structure of GGAC(UUCG)GUCC (Cheong et al., 1990; Varani et al., 1991) as the starting conformation in our attempt to understand its thermodynamic stability. We carried out molecular dynamics/free energy simulations to understand the basis for the destabilization of the C(UUCG)G loop by mutating cytosine (C7)-->uracil. Because of the limited length of simulation and the presence of kinetic barriers (solvent intervention) to the uracil-->cytosine mutation, all of our computed free energy differences are based on multiple forward simulations. Based on these calculations we find that the cytosine-->uracil mutation in the loop destabilizes it by approximately 1.5kcal/mol relative to that of the reference state, an A-form RNA but with cytosine (C7) looped out. This is the same sign and magnitude as that observed in the thermodynamic studies carried out by Varani et al.(1991). We have carried out free energy component analysis to understand the effect of mutating the cytosine residue to uracil on the thermodynamic stability of the C(UUCG)G hairpin loops. Our calculations show that the most significant contribution to the stability is from the phosphate group linking U5 and U6, which favors the cytosine residue over uracil by about 6.0 kcal/mol. The residues U5, U6, and G8 in the loop region also contribute significantly to the stability. The contributions from the salt and solvent compensate each other, indicating the dynamic nature of interactions of the environment with the nucleic acid system and the coupling between these two components.