Nearest Neighbor Parameters for Estimating the Folding Stability of RNA Including Pseudouridine

Nearest neighbor parameters are widely used in software for estimating the conformational stability of an RNA sequence folding into a specific structure. Folding stability for RNA with canonical nucleotides A, C, G, and U has been widely studied, but the same is not true for most modified nucleotides. In this work, we present a comprehensive set of nearest neighbor parameters for estimating the folding stability of RNA including pseudouridine in helical or loop contexts. These parameters are derived from 210 optical melting experiments involving helices with pseudouridine-A and pseudouridine-G pairs and with pseudouridine in loop motifs. The experiments include sequences with pseudouridine and U in the same strand, including U-A and U-G pairs, allowing us to consider the folding stability of sequences with both U and pseudouridine. On average, pseudouridine stabilizes RNA folding compared to U in an analogous motif, although this effect is sequence-context dependent. These parameters improve the modeling of folding stability for RNA secondary structures containing pseudouridine. We demonstrate that these parameters successfully model the secondary structure change for Saccharomyces cerevisiae U2 snRNA when two additional inducible pseudouridines are present. These parameters are freely available and incorporated into the RNAstructure software package.