Exploring alternative RNA structure sets using MC-flashfold and db2cm

Abstract

We created an accelerated version of MC-Fold called MC-Flashfold that allows us to compute large numbers of competing secondary structures including noncanonical base pairs. We visualize the base pairs in these sets using high quality intuitive dot plots and arc plots. Our new tools allow us to explore RNA dynamics by visualizing the competing structures in free energy bands. Here we describe how to use these tools to generate dot plots that reveal the postulated anti-terminator stem in the E. coli trp operon leader sequence. These plots show the anti-terminator hairpin loop during transcription and as a minor population of the full-length leader sequence. This is a case of switching RNA structure that had been originally postulated based on short dyad inverted repeats. Other switching RNA sequences can be analyzed by using our method.