Encoding folding paths of RNA switches

RNA co-transcriptional folding has long been suspected to play an active role in helping proper native folding of ribozymes and structured regulatory motifs in mRNA untranslated regions. Yet, the underlying mechanisms and coding requirements for efficient co-transcriptional folding remain unclear. Traditional approaches have intrinsic limitations to dissect RNA folding paths, as they rely on sequence mutations or circular permutations that typically perturb both RNA folding paths and equilibrium structures. Here, we show that exploiting sequence symmetries instead of mutations can circumvent this problem by essentially decoupling folding paths from equilibrium structures of designed RNA sequences. Using bistable RNA switches with symmetrical helices conserved under sequence reversal, we demonstrate experimentally that native and transiently formed helices can guide efficient co-transcriptional folding into either long-lived structure of these RNA switches. Their folding path is controlled by the order of helix nucleations and subsequent exchanges during transcription, and may also be redirected by transient antisense interactions. Hence, transient intra- and intermolecular base pair interactions can effectively regulate the folding of nascent RNA molecules into different native structures, provided limited coding requirements, as discussed from an information theory perspective. This constitutive coupling between RNA synthesis and RNA folding regulation may have enabled the early emergence of autonomous RNA-based regulation networks.Comment: 9 pages, 6 figure

Kinefold web server for RNA/DNA folding path and structure prediction including pseudoknots and knots

The Kinefold web server provides a web interface for stochastic folding simulations of nucleic acids on second to minute molecular time scales. Renaturation or co-transcriptional folding paths are simulated at the level of helix formation and dissociation in agreement with the seminal experimental results. Pseudoknots and topologically ‘entangled’ helices (i.e. knots) are efficiently predicted taking into account simple geometrical and topological constraints. To encourage interactivity, simulations launched as immediate jobs are automatically stopped after a few seconds and return adapted recommendations. Users can then choose to continue incomplete simulations using the batch queuing system or go back and modify suggested options in their initial query. Detailed output provide (i) a series of low free energy structures, (ii) an online animated folding path and (iii) a programmable trajectory plot focusing on a few helices of interest to each user. The service can be accessed at

Probing complex RNA structures by mechanical force

RNA secondary structures of increasing complexity are probed combining single molecule stretching experiments and stochastic unfolding/refolding simulations. We find that force-induced unfolding pathways cannot usually be interpretated by solely invoking successive openings of native helices. Indeed, typical force-extension responses of complex RNA molecules are largely shaped by stretching-induced, long-lived intermediates including non-native helices. This is first shown for a set of generic structural motifs found in larger RNA structures, and then for Escherichia coli's 1540-base long 16S ribosomal RNA, which exhibits a surprisingly well-structured and reproducible unfolding pathway under mechanical stretching. Using out-of-equilibrium stochastic simulations, we demonstrate that these experimental results reflect the slow relaxation of RNA structural rearrangements. Hence, micromanipulations of single RNA molecules probe both their native structures and long-lived intermediates, so-called "kinetic traps", thereby capturing -at the single molecular level- the hallmark of RNA folding/unfolding dynamics.Comment: 9 pages, 9 figure

Prediction and statistics of pseudoknots in RNA structures using exactly clustered stochastic simulations

Ab initio RNA secondary structure predictions have long dismissed helices interior to loops, so-called pseudoknots, despite their structural importance. Here, we report that many pseudoknots can be predicted through long time scales RNA folding simulations, which follow the stochastic closing and opening of individual RNA helices. The numerical efficacy of these stochastic simulations relies on an O(n^2) clustering algorithm which computes time averages over a continously updated set of n reference structures. Applying this exact stochastic clustering approach, we typically obtain a 5- to 100-fold simulation speed-up for RNA sequences up to 400 bases, while the effective acceleration can be as high as 100,000-fold for short multistable molecules (<150 bases). We performed extensive folding statistics on random and natural RNA sequences, and found that pseudoknots are unevenly distributed amongst RNAstructures and account for up to 30% of base pairs in G+C rich RNA sequences (Online RNA folding kinetics server including pseudoknots : http://kinefold.u-strasbg.fr/ ).Comment: 6 pages, 5 figure

Transat—A Method for Detecting the Conserved Helices of Functional RNA Structures, Including Transient, Pseudo-Knotted and Alternative Structures

The prediction of functional RNA structures has attracted increased interest, as it allows us to study the potential functional roles of many genes. RNA structure prediction methods, however, assume that there is a unique functional RNA structure and also do not predict functional features required for in vivo folding. In order to understand how functional RNA structures form in vivo, we require sophisticated experiments or reliable prediction methods. So far, there exist only a few, experimentally validated transient RNA structures. On the computational side, there exist several computer programs which aim to predict the co-transcriptional folding pathway in vivo, but these make a range of simplifying assumptions and do not capture all features known to influence RNA folding in vivo. We want to investigate if evolutionarily related RNA genes fold in a similar way in vivo. To this end, we have developed a new computational method, Transat, which detects conserved helices of high statistical significance. We introduce the method, present a comprehensive performance evaluation and show that Transat is able to predict the structural features of known reference structures including pseudo-knotted ones as well as those of known alternative structural configurations. Transat can also identify unstructured sub-sequences bound by other molecules and provides evidence for new helices which may define folding pathways, supporting the notion that homologous RNA sequence not only assume a similar reference RNA structure, but also fold similarly. Finally, we show that the structural features predicted by Transat differ from those assuming thermodynamic equilibrium. Unlike the existing methods for predicting folding pathways, our method works in a comparative way. This has the disadvantage of not being able to predict features as function of time, but has the considerable advantage of highlighting conserved features and of not requiring a detailed knowledge of the cellular environment

Attenuation of loop-receptor interactions with pseudoknot formation

RNA tetraloops can recognize receptors to mediate long-range interactions in stable natural RNAs. In vitro selected GNRA tetraloop/receptor interactions are usually more ‘G/C-rich’ than their ‘A/U-rich’ natural counterparts. They are not as widespread in nature despite comparable biophysical and chemical properties. Moreover, while AA, AC and GU dinucleotide platforms occur in natural GAAA/11 nt receptors, the AA platform is somewhat preferred to the others. The apparent preference for ‘A/U-rich’ GNRA/receptor interactions in nature might stem from an evolutionary adaptation to avoid folding traps at the level of the larger molecular context. To provide evidences in favor of this hypothesis, several riboswitches based on natural and artificial GNRA receptors were investigated in vitro for their ability to prevent inter-molecular GNRA/receptor interactions by trapping the receptor sequence into an alternative intra-molecular pseudoknot. Extent of attenuation determined by native gel-shift assays and co-transcriptional assembly is correlated to the G/C content of the GNRA receptor. Our results shed light on the structural evolution of natural long-range interactions and provide design principles for RNA-based attenuator devices to be used in synthetic biology and RNA nanobiotechnology