Assessing the Quality of Cotranscriptional Folding Simulations.

Structural changes in RNAs are an important contributor to controlling gene expression not only at the posttranscriptional stage but also during transcription. A subclass of riboswitches and RNA thermometers located in the 5' region of the primary transcript regulates the downstream functional unit - usually an ORF - through premature termination of transcription. Not only such elements occur naturally, but they are also attractive devices in synthetic biology. The possibility to design such riboswitches or RNA thermometers is thus of considerable practical interest. Since these functional RNA elements act already during transcription, it is important to model and understand the dynamics of folding and, in particular, the formation of intermediate structures concurrently with transcription. Cotranscriptional folding simulations are therefore an important step to verify the functionality of design constructs before conducting expensive and labor-intensive wet lab experiments. For RNAs, full-fledged molecular dynamics simulations are far beyond practical reach because of both the size of the molecules and the timescales of interest. Even at the simplified level of secondary structures, further approximations are necessary. The BarMap approach is based on representing the secondary structure landscape for each individual transcription step by a coarse-grained representation that only retains a small set of low-energy local minima and the energy barriers between them. The folding dynamics between two transcriptional elongation steps is modeled as a Markov process on this representation. Maps between pairs of consecutive coarse-grained landscapes make it possible to follow the folding process as it changes in response to transcription elongation. In its original implementation, the BarMap software provides a general framework to investigate RNA folding dynamics on temporally changing landscapes. It is, however, difficult to use in particular for specific scenarios such as cotranscriptional folding. To overcome this limitation, we developed the user-friendly BarMap-QA pipeline described in detail in this contribution. It is illustrated here by an elaborate example that emphasizes the careful monitoring of several quality measures. Using an iterative workflow, a reliable and complete kinetics simulation of a synthetic, transcription-regulating riboswitch is obtained using minimal computational resources. All programs and scripts used in this contribution are free software and available for download as a source distribution for Linux® or as a platform-independent Docker® image including support for Apple macOS® and Microsoft Windows®.

A Guide to Computational Cotranscriptional Folding Featuring the SRP RNA.

Although RNA molecules are synthesized via transcription, little is known about the general impact of cotranscriptional folding in vivo. We present different computational approaches for the simulation of changing structure ensembles during transcription, including interpretations with respect to experimental data from literature. Specifically, we analyze different mutations of the E. coli SRP RNA, which has been studied comparatively well in previous literature, yet the details of which specific metastable structures form as well as when they form are still under debate. Here, we combine thermodynamic and kinetic, deterministic, and stochastic models with automated and visual inspection of those systems to derive the most likely scenario of which substructures form at which point during transcription. The simulations do not only provide explanations for present experimental observations but also suggest previously unnoticed conformations that may be verified through future experimental studies.

Cotranscriptional Assembly and Native Purification of Large RNA-RNA Complexes for Structural Analyses.

Recent technological developments such as cryogenic electron microscopy (Cryo-EM) and X-ray free electron lasers (XFEL) have significantly expanded the available toolkit to visualize large, complex noncoding RNAs and their complexes. Consequently, the quality of the RNA sample, as measured by its chemical monodispersity and conformational homogeneity, has become the bottleneck that frequently precludes effective structural analyses. Here we describe a general RNA sample preparation protocol that combines cotranscriptional RNA folding and RNA-RNA complex assembly, followed by native purification of stoichiometric complexes. We illustrate and discuss the utility of this versatile method in overcoming RNA misfolding and enabling the structural and mechanistic elucidations of the T-box riboswitch-tRNA complexes.

Effect of pausing on the cotranscriptional folding kinetics of RNAs.

The pausing event in RNA molecule folding that occurs during the cotranscription process plays a critical role in cellular RNA function. Based on Conformational Resampling through Kinetic Relaxation (CRKR), this paper investigates a method to reveal the specifics of pausing for RNA cotranscriptional folding of 117-nucleotide E. coli SRP RNA and 73-nucleotide HIV-1 TAR RNA. It can be inferred from the results that pausing events generate valid cotranscriptional conformational rearrangement to protect the function structures and influence the folding pathway, which is remarkably consistent with the experimental results. Additionally, different transcription speeds result in different levels of protection capability. The folding pathway or conformational rearrangement can also be affected by a change in transcription speed after pausing site. These findings illuminate that RNAs with functional structures undergo complex rearrangement folding in pausing events, and different RNAs prefer different speeds for specific folding pathways. Moreover, pausing in particular circumstances can be employed to regulate the population of final functional structures. In conclusion, this paper offers fresh perspectives on the pausing event in the cotranscriptional folding of RNAs.

Dynamic competition between a ligand and transcription factor NusA governs riboswitch-mediated transcription regulation.

Cotranscriptional RNA folding is widely assumed to influence the timely control of gene expression, but our understanding remains limited. In bacteria, the fluoride (F-)-sensing riboswitch is a transcriptional control element essential to defend against toxic F- levels. Using this model riboswitch, we find that its ligand F- and essential bacterial transcription factor NusA compete to bind the cotranscriptionally folding RNA, opposing each other's modulation of downstream pausing and termination by RNA polymerase. Single-molecule fluorescence assays probing active transcription elongation complexes discover that NusA unexpectedly binds highly reversibly, frequently interrogating the complex for emerging, cotranscriptionally folding RNA duplexes. NusA thus fine-tunes the transcription rate in dependence of the ligand-responsive higher-order structure of the riboswitch. At the high NusA concentrations found intracellularly, this dynamic modulation is expected to lead to adaptive bacterial transcription regulation with fast response times.

A roadmap for rRNA folding and assembly during transcription.

Ribonucleoprotein (RNP) assembly typically begins during transcription when folding of the newly synthesized RNA is coupled with the recruitment of RNA-binding proteins (RBPs). Upon binding, the proteins induce structural rearrangements in the RNA that are crucial for the next steps of assembly. Focusing primarily on bacterial ribosome assembly, we discuss recent work showing that early RNA-protein interactions are more dynamic than previously supposed, and remain so, until sufficient proteins are recruited to each transcript to consolidate an entire domain of the RNP. We also review studies showing that stable assembly of an RNP competes against modification and processing of the RNA. Finally, we discuss how transcription sets the timeline for competing and cooperative RNA-RBP interactions that determine the fate of the nascent RNA. How this dance is coordinated is the focus of this review.

Cotranscriptional Kinetic Folding of RNA Secondary Structures Including Pseudoknots.

Computational prediction of ribonucleic acid (RNA) structures is an important problem in computational structural biology. Studies of RNA structure formation often assume that the process starts from a fully synthesized sequence. Experimental evidence, however, has shown that RNA folds concurrently with its elongation. We investigate RNA secondary structure formation, including pseudoknots, that takes into account the cotranscriptional effects. We propose a single-nucleotide resolution kinetic model of the folding process of RNA molecules, where the polymerase-driven elongation of an RNA strand by a new nucleotide is included as a primitive operation, together with a stochastic simulation method that implements this folding concurrently with the transcriptional synthesis. Numerical case studies show that our cotranscriptional RNA folding model can predict the formation of conformations that are favored in actual biological systems. Our new computational tool can thus provide quantitative predictions and offer useful insights into the kinetics of RNA folding.

Cotranscriptional Folding of a Bio-orthogonal Fluorescent Scaffolded RNA Origami.

The scaffolded origami technique is an attractive tool for engineering nucleic acid nanostructures. This paper demonstrates scaffolded RNA origami folding in vitro in which, for the first time, all components are transcribed simultaneously in a single-pot reaction. Double-stranded DNA sequences are transcribed by T7 RNA polymerase into scaffold and staple strands able to correctly fold in a high synthesis yield into the nanoribbon. Synthesis is successfully confirmed by atomic force microscopy, and the unpurified transcription reaction mixture is analyzed by an in gel-imaging assay where the transcribed RNA nanoribbons are able to capture the specific dye through the reconstituted split Broccoli aptamer showing a clear green fluorescent band. Finally, we simulate the RNA origami in silico using the nucleotide-level coarse-grained model oxRNA to investigate the thermodynamic stability of the assembled nanostructure in isothermal conditions over a period of time. Our work suggests that the scaffolded origami technique is a viable, and potentially more powerful, assembly alternative to the single-stranded origami technique for future in vivo applications.

Single-Molecule FRET Assay for Studying Cotranscriptional RNA Folding.

Cotranscriptional RNA folding plays important roles in gene regulation steps such as splicing, transcription termination, and translation initiation. Progression of our understanding of cotranscriptional RNA folding mechanisms is still retarded by the lacking of experimental tools to study the kinetics of cotranscriptional RNA folding properly. In this chapter, we describe fluorescence resonance energy transfer (FRET) assay that enables the study of RNA cotranscriptional folding at the single-molecule level.

Regulation of the thiamine pyrophosphate (TPP)-sensing riboswitch in NMT1 mRNA from Neurospora crassa.

The expression of Neurospora crassa NMT1 involved in thiamine pyrophosphate (TPP) metabolism is regulated at the level of mRNA splicing by a TPP-sensing riboswitch within the precursor NMT1 mRNA. Here, using the systematic helix-based computational method, we investigated the regulation of this riboswitch. We find that the function of the riboswitch does not depend on the transcription process. Whether TPP is present or not, the riboswitch predominately folds into the ON state, while the OFF state aptamer structure does not appear during transcription. Since the transition from the ON state to the aptamer structure is extremely slow, TPP may interact with the RNA before full formation of the aptamer structure, promoting the switch flipping. The potential to fully form helix P0 of the ON state is necessary to restore ligand-dependent gene control by the riboswitch.

Ruleset Optimization on Isomorphic Oritatami Systems.

We study an optimization problem of a computational folding model, proving its hardness and proposing heuristic algorithms. RNA cotranscriptional folding refers to the phenomenon in which an RNA transcript folds upon itself while being synthesized out of a gene. An oritatami model (OM) is a computational model of this phenomenon that lets its sequence of beads (abstract molecules) fold cotranscriptionally by the interactions between beads, according to its ruleset. We study the problem of reducing the ruleset size, while keeping the terminal conformations geometrically the same. We first prove the hardness of finding the smallest ruleset, and then suggest two approaches that reduce the ruleset size efficiently.

An emergent understanding of strand displacement in RNA biology.

DNA and RNA are generally regarded as central molecules in molecular biology. Recent advancements in the field of DNA/RNA nanotechnology successfully used DNA/RNA as programmable molecules to construct molecular machines and nanostructures with predefined shapes and functions. The key mechanism for dynamic control of the conformations of these DNA/RNA nanodevices is a reaction called strand displacement, in which one strand in a formed duplex is replaced by a third invading strand. While DNA/RNA strand displacement has mainly been used to de novo design molecular devices, we argue in this review that this reaction is also likely to play a key role in multiple cellular events such as gene recombination, CRISPR-based genome editing, and RNA cotranscriptional folding. We introduce the general mechanism of strand displacement reaction, give examples of its use in the construction of molecular machines, and finally review natural processes having characteristic which suggest that strand displacement is occurring.

Out-of-Plane Aptamer Functionalization of RNA Three-Helix Tiles.

Co-transcriptionally folding RNA nanostructures have great potential as biomolecular scaffolds, which can be used to organize small molecules or proteins into spatially ordered assemblies. Here, we develop an RNA tile composed of three parallel RNA double helices, which can associate into small hexagonal assemblies via kissing loop interactions between its two outer helices. The inner RNA helix is modified with an RNA motif found in the internal ribosome entry site (IRES) of the hepatitis C virus (HCV), which provides a 90° bend. This modification is used to functionalize the RNA structures with aptamers pointing perpendicularly away from the tile plane. We demonstrate modifications with the fluorogenic malachite green and Spinach aptamers as well with the protein-binding PP7 and streptavidin aptamers. The modified structures retain the ability to associate into larger assemblies, representing a step towards RNA hybrid nanostructures extending in three dimensions.

Single-molecule FRET studies on the cotranscriptional folding of a thiamine pyrophosphate riboswitch.

Because RNAs fold as they are being synthesized, their transcription rate can affect their folding. Here, we report the results of single-molecule fluorescence studies that characterize the ligand-dependent cotranscriptional folding of the Escherichia coli thiM riboswitch that regulates translation. We found that the riboswitch aptamer folds into the "off" conformation independent of its ligand, but switches to the "on" conformation during transcriptional pausing near the translational start codon. Ligand binding maintains the riboswitch in the off conformation during transcriptional pauses. We expect our assay will permit the controlled study of the two main physical mechanisms that regulate cotranscriptional folding: transcriptional pausing and transcriptional speed.

Kinetic folding design of aptazyme-regulated expression devices as riboswitches for metabolic engineering.

Recent developments in the fields of synthetic biology and metabolic engineering have opened the doors for the microbial production of biofuels and other valuable organic compounds. There remain, however, significant metabolic hurdles to the production of these compounds in cost-effective quantities. This is due, in part, to mismatches between the metabolic engineer's desire for high yields and the microbe's desire to survive. Many valuable compounds, or the intermediates necessary for their biosynthesis, prove deleterious at the desired production concentrations. One potential solution to these toxicity-related issues is the implementation of nonnative dynamic genetic control mechanisms that sense excessively high concentrations of metabolic intermediates and respond accordingly to alleviate their impact. One potential class of dynamic regulator is the riboswitch: cis-acting RNA elements that regulate the expression of downstream genes based on the presence of an effector molecule. Here, we present combined methods for constructing aptazyme-regulated expression devices (aREDs) through computational cotranscriptional kinetic folding design and experimental validation. These approaches can be used to engineer aREDs within novel genetic contexts for the predictable, dynamic regulation of gene expression in vivo.

Cotranscriptionally folded RNA nanostructures pave the way to intracellular nanofabrication.

A very attractive goal in nanotechnology is to manufacture smart nanodevices that integrate multiple biological/biomedical functions and autonomously function in vivo in a predefined and well-controlled manner. For decades, researchers have been developing many different ways toward this target, using bottom-up assembly of types of nanomaterials or top-down fabrication of devices with nanometer-scale precision. However, the practical application of these nanodevices remains challenging. One possible barrier lies in the spatiotemporal separation between fabrication and use, which poses a great challenge for the non-invasive delivery of fully functional nanodevice into live cells. Indeed, cells themselves are highly complex natural machines with membrane barriers and finely regulated pathways for intracellular delivery. However, there is plenty of evidence that nanomaterials or nanodevices are easily aggregated or trapped inside of the cells.