Helix-Based RNA Landscape Partition and Alternative Secondary Structure Determination

cRNAsp12 Web Server for the Prediction of Circular RNA Secondary Structures and Stabilities

Circular RNAs (circRNAs) are a novel class of non-coding RNA that, unlike linear RNAs, form a covalently closed loop without the 5' and 3' ends. Growing evidence shows that circular RNAs play important roles in life processes and have great potential implications in clinical and research fields. The accurate modeling of circRNAs structure and stability has far-reaching impact on our understanding of their functions and our ability to develop RNA-based therapeutics. The cRNAsp12 server offers a user-friendly web interface to predict circular RNA secondary structures and folding stabilities from the sequence. Through the helix-based landscape partitioning strategy, the server generates distinct ensembles of structures and predicts the minimal free energy structures for each ensemble with the recursive partition function calculation and backtracking algorithms. For structure predictions in the limited structural ensemble, the server also provides users with the option to set the structural constraints of forcing the base pairs and/or forcing the unpaired bases, such that only structures that meet the criteria are enumerated recursively.

Landscape Zooming toward the Prediction of RNA Cotranscriptional Folding

RNA molecules fold as they are transcribed. Cotranscriptional folding of RNA plays a critical role in RNA functions in vivo. Present computational strategies focus on simulations where large structural changes may not be completely sampled. Here, we describe an alternative approach to predicting cotranscriptional RNA folding by zooming in and out of the RNA folding energy landscape. By classifying the RNA structural ensemble into "partitions" based on long, stable helices, we zoom out of the landscape and predict the overall slow folding kinetics from the interpartition kinetic network, and for each interpartition transition, we zoom in on the landscape to simulate the kinetics. Applications of the model to the 117-nucleotide E. coli SRP RNA and the 59-nucleotide HIV-1 TAR RNA show agreements with the experimental data and new structural and kinetic insights into biologically significant conformational switches and pathways for these important systems. This approach, by zooming in/out of an RNA folding landscape at different resolutions, might allow us to treat large RNAs in vivo with transcriptional pause, transcription speed, and other in vivo effects.