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

Localized Joule heating produced by ion current focusing through micron-size holes

We provide an experimental demonstration that the focusing of ionic currents in a micron size hole connecting two chambers can produce local temperature increases of up to $100^\circ$ C with gradients as large as $1^\circ$ K$\mu m^{-1}$. We find a good agreement between the measured temperature profiles and a finite elements-based numerical calculation. We show how the thermal gradients can be used to measure the full melting profile of DNA duplexes within a region of 40 $\mu$m. The possibility to produce even larger gradients using sub-micron pores is discussed.Comment: 3 pages, accepted to Appl. Phys. Lett

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

Hierarchy and Feedback in the Evolution of the E. coli Transcription Network

The E.coli transcription network has an essentially feedforward structure, with, however, abundant feedback at the level of self-regulations. Here, we investigate how these properties emerged during evolution. An assessment of the role of gene duplication based on protein domain architecture shows that (i) transcriptional autoregulators have mostly arisen through duplication, while (ii) the expected feedback loops stemming from their initial cross-regulation are strongly selected against. This requires a divergent coevolution of the transcription factor DNA-binding sites and their respective DNA cis-regulatory regions. Moreover, we find that the network tends to grow by expansion of the existing hierarchical layers of computation, rather than by addition of new layers. We also argue that rewiring of regulatory links due to mutation/selection of novel transcription factor/DNA binding interactions appears not to significantly affect the network global hierarchy, and that horizontally transferred genes are mainly added at the bottom, as new target nodes. These findings highlight the important evolutionary roles of both duplication and selective deletion of crosstalks between autoregulators in the emergence of the hierarchical transcription network of E.coli.Comment: to appear in PNA

Polymer Induced Bundling of F-actin and the Depletion Force

The inert polymer polyethylene glycol (PEG) induces a "bundling" phenomenon in F-actin solutions when its concentration exceeds a critical onset value C_o. Over a limited range of PEG molecular weight and ionic strength, C_o can be expressed as a function of these two variables. The process is reversible, but hysteresis is also observed in the dissolution of the bundles, with ionic strength having a large influence. Additional actin filaments are able to join previously formed bundles. Little, if any, polymer is associated with the bundle structure. Continuum estimates of the Asakura-Oosawa depletion force, Coulomb repulsion, and van der Waals potential are combined for a partial explanation of the bundling effect and hysteresis. Conjectures are presented concerning the apparent limit in bundle size

Anomalous fluctuations of active polar filaments

Using a simple model, we study the fluctuating dynamics of inextensible, semiflexible polar filaments interacting with active and directed force generating centres such as molecular motors. Taking into account the fact that the activity occurs on time-scales comparable to the filament relaxation time, we obtain some unexpected differences between both the steady-state and dynamical behaviour of active as compared to passive filaments. For the statics, the filaments have a {novel} length-scale dependent rigidity. Dynamically, we find strongly enhanced anomalous diffusion.Comment: 5 pages, 3 figure

Statistical mechanics of secondary structures formed by random RNA sequences

The formation of secondary structures by a random RNA sequence is studied as a model system for the sequence-structure problem omnipresent in biopolymers. Several toy energy models are introduced to allow detailed analytical and numerical studies. First, a two-replica calculation is performed. By mapping the two-replica problem to the denaturation of a single homogeneous RNA in 6-dimensional embedding space, we show that sequence disorder is perturbatively irrelevant, i.e., an RNA molecule with weak sequence disorder is in a molten phase where many secondary structures with comparable total energy coexist. A numerical study of various models at high temperature reproduces behaviors characteristic of the molten phase. On the other hand, a scaling argument based on the extremal statistics of rare regions can be constructed to show that the low temperature phase is unstable to sequence disorder. We performed a detailed numerical study of the low temperature phase using the droplet theory as a guide, and characterized the statistics of large-scale, low-energy excitations of the secondary structures from the ground state structure. We find the excitation energy to grow very slowly (i.e., logarithmically) with the length scale of the excitation, suggesting the existence of a marginal glass phase. The transition between the low temperature glass phase and the high temperature molten phase is also characterized numerically. It is revealed by a change in the coefficient of the logarithmic excitation energy, from being disorder dominated to entropy dominated.Comment: 24 pages, 16 figure