Computationally reconstructing cotranscriptional RNA folding from experimental data reveals rearrangement of non-native folding intermediates

Non-averaged single-molecule tertiary structures reveal RNA self-folding through individual-particle cryo-electron tomography

Large-scale and continuous conformational changes in the RNA self-folding process present significant challenges for structural studies, often requiring trade-offs between resolution and observational scope. Here, we utilize individual-particle cryo-electron tomography (IPET) to examine the post-transcriptional self-folding process of designed RNA origami 6-helix bundle with a clasp helix (6HBC). By avoiding selection, classification, averaging, or chemical fixation and optimizing cryo-ET data acquisition parameters, we reconstruct 120 three-dimensional (3D) density maps from 120 individual particles at an electron dose of no more than 168 e-Å-2, achieving averaged resolutions ranging from 23 to 35 Å, as estimated by Fourier shell correlation (FSC) at 0.5. Each map allows us to identify distinct RNA helices and determine a unique tertiary structure. Statistical analysis of these 120 structures confirms two reported conformations and reveals a range of kinetically trapped, intermediate, and highly compacted states, demonstrating a maturation folding landscape likely driven by helix-helix compaction interactions.

Sequential structure probing of cotranscriptional RNA folding intermediates

Cotranscriptional RNA folding pathways typically involve the sequential formation of folding intermediates. Existing methods for cotranscriptional RNA structure probing map the structure of nascent RNA in the context of a terminally arrested transcription elongation complex. Consequently, the rearrangement of RNA structures as nucleotides are added to the transcript can be inferred but is not assessed directly. To address this limitation, we have developed linked- m ultipoint T ranscription E longation C omplex RNA structure prob ing (TECprobe-LM), which assesses the cotranscriptional rearrangement of RNA structures by sequentially positioning E. coli RNAP at two or more points within a DNA template so that nascent RNA can be chemically probed. We validated TECprobe-LM by measuring known folding events that occur within the E. coli signal recognition particle RNA, Clostridium beijerinckii pfl ZTP riboswitch, and Bacillus cereus crcB fluoride riboswitch folding pathways. Our findings establish TECprobe-LM as a strategy for detecting cotranscriptional RNA folding events directly using chemical probing.

The role of RNA structure in 3' end processing in eukaryotes

Maturation of pre-mRNA into fully functional mRNA involves a series of highly coordinated steps that are essential for eukaryotic gene expression. RNA structure has been found to play regulatory roles in many of these steps, including cleavage, polyadenylation, and termination. Recent advances in structure probing techniques have been instrumental in revealing how nascent transcript conformation contributes to these dynamic, co-transcriptional processes. In this review, we present examples where RNA structure affects accessibility and/or function of key processing enzymes, thereby influencing the efficiency and precision of 3' end processing machinery. We also discuss emerging technologies that could further enhance our understanding of RNA structure mediated regulation of 3' end processing.

Optimized periphery-core interface increases fitness of the Bacillus subtilisglmS ribozyme

Like other functional RNAs, ribozymes encode a conserved catalytic center supported by peripheral domains that vary among ribozyme sub-families. To understand how core-periphery interactions contribute to ribozyme fitness, we compared the cleavage kinetics of all single base substitutions at 152 sites across the Bacillus subtilis glmS ribozyme by high-throughput sequencing (k-seq). The in vitro activity map mirrored phylogenetic sequence conservation in glmS ribozymes, indicating that biological fitness reports all biochemically important positions. The k-seq results and folding assays showed that most deleterious mutations lower activity by impairing ribozyme self-assembly. All-atom molecular dynamics simulations of the complete ribozyme revealed how individual mutations in the core or the IL4 peripheral loop introduce a non-native tertiary interface that rewires the catalytic center, eliminating activity. We conclude that the need to avoid non-native helix packing powerfully constrains the evolution of tertiary structure motifs in RNA.

Single-molecule force spectroscopy of toehold-mediated strand displacement

Toehold-mediated strand displacement (TMSD) is extensively utilized in dynamic DNA nanotechnology and for a wide range of DNA or RNA-based reaction circuits. Investigation of TMSD kinetics typically relies on bulk fluorescence measurements providing effective, bulk-averaged reaction rates. Information on individual molecules or even base pairs is scarce. In this work, we explore the dynamics of strand displacement processes at the single-molecule level using single-molecule force spectroscopy with a microfluidics-enhanced optical trap supported by state-of-the-art coarse-grained simulations. By applying force, we can trigger and observe TMSD in real-time with microsecond and nanometer resolution. We find TMSD proceeds very rapidly under load with single step times of 1 µs. Tuning invasion efficiency by introducing mismatches allows studying thousands of forward/backward invasion events on a single molecule and analyze the kinetics of the invasion process. Extrapolation to zero force reveals single step times for DNA invading DNA four times faster than for RNA invading RNA. We also study the kinetics of DNA invading RNA, a process that in the absence of force would rarely occur. Our results reveal the importance of sequence effects for the TMSD process and have relevance for a wide range of applications in nucleic acid nanotechnology and synthetic biology.

Cotranscriptional Folding of a 5' Stem-loop in the Escherichia coli tbpA Riboswitch at Single-nucleotide Resolution

Transcription elongation is one of the most important processes in the cell. During RNA polymerase elongation, the folding of nascent transcripts plays crucial roles in the genetic decision. Bacterial riboswitches are prime examples of RNA regulators that control gene expression by altering their structure upon metabolite sensing. It was previously revealed that the thiamin pyrophosphate-sensing tbpA riboswitch in Escherichia coli cotranscriptionally adopts three main structures leading to metabolite sensing. Here, using single-molecule FRET, we characterize the transition in which the first nascent structure, a 5' stem-loop, is unfolded during transcription elongation to form the ligand-binding competent structure. Our results suggest that the structural transition occurs in a relatively abrupt manner, i.e., within a 1-2 nucleotide window. Furthermore, a highly dynamic structural exchange is observed, indicating that riboswitch transcripts perform rapid sampling of nascent co-occurring structures. We also observe that the presence of the RNAP stabilizes the 5' stem-loop along the elongation process, consistent with RNAP interacting with the 5' stem-loop. Our study emphasizes the role of early folding stem-loop structures in the cotranscriptional formation of complex RNA molecules involved in genetic regulation.

Therapeutic applications of cell engineering using mRNA technology

Engineering and reprogramming cells has significant therapeutic potential to treat a wide range of diseases, by replacing missing or defective proteins, to provide transcription factors (TFs) to reprogram cell phenotypes, or to provide enzymes such as RNA-guided Cas9 derivatives for CRISPR-based cell engineering. While viral vector-mediated gene transfer has played an important role in this field, the use of mRNA avoids safety concerns associated with the integration of DNA into the host cell genome, making mRNA particularly attractive for in vivo applications. Widespread application of mRNA for cell engineering is limited by its instability in the biological environment and challenges involved in the delivery of mRNA to its target site. In this review, we examine challenges that must be overcome to develop effective therapeutics.

Co-transcriptional gene regulation in eukaryotes and prokaryotes

Many steps of RNA processing occur during transcription by RNA polymerases. Co-transcriptional activities are deemed commonplace in prokaryotes, in which the lack of membrane barriers allows mixing of all gene expression steps, from transcription to translation. In the past decade, an extraordinary level of coordination between transcription and RNA processing has emerged in eukaryotes. In this Review, we discuss recent developments in our understanding of co-transcriptional gene regulation in both eukaryotes and prokaryotes, comparing methodologies and mechanisms, and highlight striking parallels in how RNA polymerases interact with the machineries that act on nascent RNA. The development of RNA sequencing and imaging techniques that detect transient transcription and RNA processing intermediates has facilitated discoveries of transcription coordination with splicing, 3'-end cleavage and dynamic RNA folding and revealed physical contacts between processing machineries and RNA polymerases. Such studies indicate that intron retention in a given nascent transcript can prevent 3'-end cleavage and cause transcriptional readthrough, which is a hallmark of eukaryotic cellular stress responses. We also discuss how coordination between nascent RNA biogenesis and transcription drives fundamental aspects of gene expression in both prokaryotes and eukaryotes.

Kinetic pathway of HIV-1 TAR cotranscriptional folding

The Trans-Activator Receptor (TAR) RNA, located at the 5'-end untranslated region (5' UTR) of the human immunodeficiency virus type 1 (HIV-1), is pivotal in the virus's life cycle. As the initial functional domain, it folds during the transcription of viral mRNA. Although TAR's role in recruiting the Tat protein for trans-activation is established, the detailed kinetic mechanisms at play during early transcription, especially at points of temporary transcriptional pausing, remain elusive. Moreover, the precise physical processes of transcriptional pause and subsequent escape are not fully elucidated. This study focuses on the folding kinetics of TAR and the biological implications by integrating computer simulations of RNA folding during transcription with nuclear magnetic resonance (NMR) spectroscopy data. The findings reveal insights into the folding mechanism of a non-native intermediate that triggers transcriptional pause, along with different folding pathways leading to transcriptional pause and readthrough. The profiling of the cotranscriptional folding pathway and identification of kinetic structural intermediates reveal a novel mechanism for viral transcriptional regulation, which could pave the way for new antiviral drug designs targeting kinetic cotranscriptional folding pathways in viral RNAs.

The effect of pseudoknot base pairing on cotranscriptional structural switching of the fluoride riboswitch

A central question in biology is how RNA sequence changes influence dynamic conformational changes during cotranscriptional folding. Here we investigated this question through the study of transcriptional fluoride riboswitches, non-coding RNAs that sense the fluoride anion through the coordinated folding and rearrangement of a pseudoknotted aptamer domain and a downstream intrinsic terminator expression platform. Using a combination of Escherichia coli RNA polymerase in vitro transcription and cellular gene expression assays, we characterized the function of mesophilic and thermophilic fluoride riboswitch variants. We showed that only variants containing the mesophilic pseudoknot function at 37°C. We next systematically varied the pseudoknot sequence and found that a single wobble base pair is critical for function. Characterizing thermophilic variants at 65°C through Thermus aquaticus RNA polymerase in vitro transcription showed the importance of this wobble pair for function even at elevated temperatures. Finally, we performed all-atom molecular dynamics simulations which supported the experimental findings, visualized the RNA structure switching process, and provided insight into the important role of magnesium ions. Together these studies provide deeper insights into the role of riboswitch sequence in influencing folding and function that will be important for understanding of RNA-based gene regulation and for synthetic biology applications.

Flipping the script: Understanding riboswitches from an alternative perspective

Riboswitches are broadly distributed regulatory elements most frequently found in the 5'-leader sequence of bacterial mRNAs that regulate gene expression in response to the binding of a small molecule effector. The occupancy status of the ligand-binding aptamer domain manipulates downstream information in the message that instructs the expression machinery. Currently, there are over 55 validated riboswitch classes, where each class is defined based on the identity of the ligand it binds and/or sequence and structure conservation patterns within the aptamer domain. This classification reflects an "aptamer-centric" perspective that dominates our understanding of riboswitches. In this review, we propose a conceptual framework that groups riboswitches based on the mechanism by which RNA manipulates information directly instructing the expression machinery. This scheme does not replace the established aptamer domain-based classification of riboswitches but rather serves to facilitate hypothesis-driven investigation of riboswitch regulatory mechanisms. Based on current bioinformatic, structural, and biochemical studies of a broad spectrum of riboswitches, we propose three major mechanistic groups: (1) "direct occlusion", (2) "interdomain docking", and (3) "strand exchange". We discuss the defining features of each group, present representative examples of riboswitches from each group, and illustrate how these RNAs couple small molecule binding to gene regulation. While mechanistic studies of the occlusion and docking groups have yielded compelling models for how these riboswitches function, much less is known about strand exchange processes. To conclude, we outline the limitations of our mechanism-based conceptual framework and discuss how critical information within riboswitch expression platforms can inform gene regulation.

RNA-based controllers for engineering gene and cell therapies

Engineered RNA-based genetic controllers provide compact, tunable, post-transcriptional gene regulation. As RNA devices are generally small, these devices are portable to DNA and RNA viral vectors. RNA tools have recently expanded to allow reading and editing of endogenous RNAs for profiling and programming of transcriptional states. With their expanded capabilities and highly compact, modular, and programmable nature, RNA-based controllers will support greater safety, efficacy, and performance in gene and cell-based therapies. In this review, we highlight RNA-based controllers and their potential as user-guided and autonomous systems for control of gene and cell-based therapies.

Regulation of bacterial gene expression by non-coding RNA: It is all about time!

Commensal and pathogenic bacteria continuously evolve to survive in diverse ecological niches by efficiently coordinating gene expression levels in their ever-changing environments. Regulation through the RNA transcript itself offers a faster and more cost-effective way to adapt than protein-based mechanisms and can be leveraged for diagnostic or antimicrobial purposes. However, RNA can fold into numerous intricate, not always functional structures that both expand and obscure the plethora of roles that regulatory RNAs serve within the cell. Here, we review the current knowledge of bacterial non-coding RNAs in relation to their folding pathways and interactions. We posit that co-transcriptional folding of these transcripts ultimately dictates their downstream functions. Elucidating the spatiotemporal folding of non-coding RNAs during transcription therefore provides invaluable insights into bacterial pathogeneses and predictive disease diagnostics. Finally, we discuss the implications of co-transcriptional folding andapplications of RNAs for therapeutics and drug targets.

Tailor made: the art of therapeutic mRNA design

mRNA medicine is a new and rapidly developing field in which the delivery of genetic information in the form of mRNA is used to direct therapeutic protein production in humans. This approach, which allows for the quick and efficient identification and optimization of drug candidates for both large populations and individual patients, has the potential to revolutionize the way we prevent and treat disease. A key feature of mRNA medicines is their high degree of designability, although the design choices involved are complex. Maximizing the production of therapeutic proteins from mRNA medicines requires a thorough understanding of how nucleotide sequence, nucleotide modification and RNA structure interplay to affect translational efficiency and mRNA stability. In this Review, we describe the principles that underlie the physical stability and biological activity of mRNA and emphasize their relevance to the myriad considerations that factor into therapeutic mRNA design.

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.

Programmable modulation of ribosomal frameshifting by mRNA targeting CRISPR-Cas12a system

Minus 1 programmed ribosomal frameshifting (-1 PRF) is a conserved translational regulation event essential for critical biological processes, including the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) replication. Efficient trans-modulation of the structured RNA element crucial to -1 PRF will endow the therapeutic application. Here, we demonstrate that CRISPR RNA can stimulate efficient -1 PRF. Assembled CRISPR-Cas12a, but not CRISPR-Cas9, complex further enhances -1 PRF efficiency through its higher capacity to stall translating ribosomes. We additionally perform CRISPR-Cas12a targeting to impair the SARS-CoV-2 frameshifting pseudoknot structure via a focused screening. We demonstrate that targeting CRISPR-Cas12a results in more than 70% suppression of -1 PRF in vitro and about 50% suppression in mammalian cells. Our results show the expanded function of the CRISPR-Cas12 system in modulating -1 PRF efficiency through stalling ribosomes and deforming frameshifting stimulatory signals, which could serve as a new strategy for future coronavirus pandemics.

The Effect of Pseudoknot Base Pairing on Cotranscriptional Structural Switching of the Fluoride Riboswitch

A central question in biology is how RNA sequence changes influence dynamic conformational changes during cotranscriptional folding. Here we investigated this question through the study of transcriptional fluoride riboswitches, non-coding RNAs that sense the fluoride anion through the coordinated folding and rearrangement of a pseudoknotted aptamer domain and a downstream intrinsic terminator expression platform. Using a combination of E. coli RNA polymerase in vitro transcription and cellular gene expression assays, we characterized the function of mesophilic and thermophilic fluoride riboswitch variants. We showed that only variants containing the mesophilic pseudoknot function at 37 °C. We next systematically varied the pseudoknot sequence and found that a single wobble base pair is critical for function. Characterizing thermophilic variants at 65 °C through Thermus aquaticus RNA polymerase in vitro transcription showed the importance of this wobble pair for function even at elevated temperatures. Finally, we performed all-atom molecular dynamics simulations which supported the experimental findings, visualized the RNA structure switching process, and provided insight into the important role of magnesium ions. Together these studies provide deeper insights into the role of riboswitch sequence in influencing folding and function that will be important for understanding of RNA-based gene regulation and for synthetic biology applications.

Observation of coordinated RNA folding events by systematic cotranscriptional RNA structure probing

RNA begins to fold as it is transcribed by an RNA polymerase. Consequently, RNA folding is constrained by the direction and rate of transcription. Understanding how RNA folds into secondary and tertiary structures therefore requires methods for determining the structure of cotranscriptional folding intermediates. Cotranscriptional RNA chemical probing methods accomplish this by systematically probing the structure of nascent RNA that is displayed from an RNA polymerase. Here, we describe a concise, high-resolution cotranscriptional RNA chemical probing procedure called variable length Transcription Elongation Complex RNA structure probing (TECprobe-VL). We demonstrate the accuracy and resolution of TECprobe-VL by replicating and extending previous analyses of ZTP and fluoride riboswitch folding and mapping the folding pathway of a ppGpp-sensing riboswitch. In each system, we show that TECprobe-VL identifies coordinated cotranscriptional folding events that mediate transcription antitermination. Our findings establish TECprobe-VL as an accessible method for mapping cotranscriptional RNA folding pathways.

A transient intermediate RNA structure underlies the regulatory function of the E. coli thiB TPP translational riboswitch

Riboswitches are cis-regulatory RNA elements that regulate gene expression in response to ligand binding through the coordinated action of a ligand-binding aptamer domain (AD) and a downstream expression platform (EP). Previous studies of transcriptional riboswitches have uncovered diverse examples that utilize structural intermediates that compete with the AD and EP folds to mediate the switching mechanism on the timescale of transcription. Here we investigate whether similar intermediates are important for riboswitches that control translation by studying the Escherichia coli thiB thiamin pyrophosphate (TPP) riboswitch. Using cellular gene expression assays, we first confirmed that the riboswitch acts at the level of translational regulation. Deletion mutagenesis showed the importance of the AD-EP linker sequence for riboswitch function. Sequence complementarity between the linker region and the AD P1 stem suggested the possibility of an intermediate nascent RNA structure called the antisequestering stem that could mediate the thiB switching mechanism. Experimentally informed secondary structure models of the thiB folding pathway generated from chemical probing of nascent thiB structures in stalled transcription elongation complexes confirmed the presence of the antisequestering stem, and showed it may form cotranscriptionally. Additional mutational analysis showed that mutations to the antisequestering stem break or bias thiB function according to whether the antisequestering stem or P1 is favored. This work provides an important example of intermediate structures that compete with AD and EP folds to implement riboswitch mechanisms.

Linking folding dynamics and function of SAM/SAH riboswitches at the single molecule level

Riboswitches are regulatory elements found in bacterial mRNAs that control downstream gene expression through ligand-induced conformational changes. Here, we used single-molecule FRET to map the conformational landscape of the translational SAM/SAH riboswitch and probe how co-transcriptional ligand-induced conformational changes affect its translation regulation function. Riboswitch folding is highly heterogeneous, suggesting a rugged conformational landscape that allows for sampling of the ligand-bound conformation even in the absence of ligand. The addition of ligand shifts the landscape, favoring the ligand-bound conformation. Mutation studies identified a key structural element, the pseudoknot helix, that is crucial for determining ligand-free conformations and their ligand responsiveness. We also investigated ribosomal binding site accessibility under two scenarios: pre-folding and co-transcriptional folding. The regulatory function of the SAM/SAH riboswitch involves kinetically favoring ligand binding, but co-transcriptional folding reduces this preference with a less compact initial conformation that exposes the Shine-Dalgarno sequence and takes min to redistribute to more compact conformations of the pre-folded riboswitch. Such slow equilibration decreases the effective ligand affinity. Overall, our study provides a deeper understanding of the complex folding process and how the riboswitch adapts its folding pattern in response to ligand, modulates ribosome accessibility and the role of co-transcriptional folding in these processes.

DrForna: visualization of cotranscriptional folding

MOTIVATION

Understanding RNA folding at the level of secondary structures can give important insights concerning the function of a molecule. We are interested to learn how secondary structures change dynamically during transcription, as well as whether particular secondary structures form already during or only after transcription. While different approaches exist to simulate cotranscriptional folding, the current strategies for visualization are lagging behind. New, more suitable approaches are necessary to help with exploring the generated data from cotranscriptional folding simulations.

RESULTS

We present DrForna, an interactive visualization app for viewing the time course of a cotranscriptional RNA folding simulation. Specifically, users can scroll along the time axis and see the population of structures that are present at any particular time point.

AVAILABILITY AND IMPLEMENTATION

DrForna is a JavaScript project available on Github at https://github.com/ViennaRNA/drforna and deployed at https://viennarna.github.io/drforna.

Pre-mRNA splicing and its cotranscriptional connections

Transcription of eukaryotic genes by RNA polymerase II (Pol II) yields RNA precursors containing introns that must be spliced out and the flanking exons ligated together. Splicing is catalyzed by a dynamic ribonucleoprotein complex called the spliceosome. Recent evidence has shown that a large fraction of splicing occurs cotranscriptionally as the RNA chain is extruded from Pol II at speeds of up to 5 kb/minute. Splicing is more efficient when it is tethered to the transcription elongation complex, and this linkage permits functional coupling of splicing with transcription. We discuss recent progress that has uncovered a network of connections that link splicing to transcript elongation and other cotranscriptional RNA processing events.

RNA folding pathways from all-atom simulations with a variationally improved history-dependent bias

Atomically detailed simulations of RNA folding have proven very challenging in view of the difficulties of developing realistic force fields and the intrinsic computational complexity of sampling rare conformational transitions. As a step forward in tackling these issues, we extend to RNA an enhanced path-sampling method previously successfully applied to proteins. In this scheme, the information about the RNA's native structure is harnessed by a soft history-dependent biasing force promoting the generation of productive folding trajectories in an all-atom force field with explicit solvent. A rigorous variational principle is then applied to minimize the effect of the bias. Here, we report on an application of this method to RNA molecules from 20 to 47 nucleotides long and increasing topological complexity. By comparison with analog simulations performed on small proteins with similar size and architecture, we show that the RNA folding landscape is significantly more frustrated, even for relatively small chains with a simple topology. The predicted RNA folding mechanisms are found to be consistent with the available experiments and some of the existing coarse-grained models. Due to its computational performance, this scheme provides a promising platform to efficiently gather atomistic RNA folding trajectories, thus retain the information about the chemical composition of the sequence.

Tertiary structure of single-instant RNA molecule reveals folding landscape

The folding of RNA and protein molecules during their synthesis is a crucial self-assembly process that nature employs to convert genetic information into the complex molecular machinery that supports life. Misfolding events are the cause of several diseases, and the folding pathway of central biomolecules, such as the ribosome, is strictly regulated by programmed maturation processes and folding chaperones. However, the dynamic folding processes are challenging to study because current structure determination methods heavily rely on averaging, and existing computational methods do not efficiently simulate non-equilibrium dynamics. Here we utilize individual-particle cryo-electron tomography (IPET) to investigate the folding landscape of a rationally designed RNA origami 6-helix bundle that undergoes slow maturation from a "young" to "mature" conformation. By optimizing the IPET imaging and electron dose conditions, we obtain 3D reconstructions of 120 individual particles at resolutions ranging from 23-35 Å, enabling us first-time to observe individual RNA helices and tertiary structures without averaging. Statistical analysis of 120 tertiary structures confirms the two main conformations and suggests a possible folding pathway driven by helix-helix compaction. Studies of the full conformational landscape reveal both trapped states, misfolded states, intermediate states, and fully compacted states. The study provides novel insight into RNA folding pathways and paves the way for future studies of the energy landscape of molecular machines and self-assembly processes.

Design Approaches to Expand the Toolkit for Building Cotranscriptionally Encoded RNA Strand Displacement Circuits

Cotranscriptionally encoded RNA strand displacement (ctRSD) circuits are an emerging tool for programmable molecular computation, with potential applications spanning in vitro diagnostics to continuous computation inside living cells. In ctRSD circuits, RNA strand displacement components are continuously produced together via transcription. These RNA components can be rationally programmed through base pairing interactions to execute logic and signaling cascades. However, the small number of ctRSD components characterized to date limits circuit size and capabilities. Here, we characterize over 200 ctRSD gate sequences, exploring different input, output, and toehold sequences and changes to other design parameters, including domain lengths, ribozyme sequences, and the order in which gate strands are transcribed. This characterization provides a library of sequence domains for engineering ctRSD components, i.e., a toolkit, enabling circuits with up to 4-fold more inputs than previously possible. We also identify specific failure modes and systematically develop design approaches that reduce the likelihood of failure across different gate sequences. Lastly, we show the ctRSD gate design is robust to changes in transcriptional encoding, opening a broad design space for applications in more complex environments. Together, these results deliver an expanded toolkit and design approaches for building ctRSD circuits that will dramatically extend capabilities and potential applications.

How does precursor RNA structure influence RNA processing and gene expression?

RNA is a fundamental biomolecule that has many purposes within cells. Due to its single-stranded and flexible nature, RNA naturally folds into complex and dynamic structures. Recent technological and computational advances have produced an explosion of RNA structural data. Many RNA structures have regulatory and functional properties. Studying the structure of nascent RNAs is particularly challenging due to their low abundance and long length, but their structures are important because they can influence RNA processing. Precursor RNA processing is a nexus of pathways that determines mature isoform composition and that controls gene expression. In this review, we examine what is known about human nascent RNA structure and the influence of RNA structure on processing of precursor RNAs. These known structures provide examples of how other nascent RNAs may be structured and show how novel RNA structures may influence RNA processing including splicing and polyadenylation. RNA structures can be targeted therapeutically to treat disease.

A trailing ribosome speeds up RNA polymerase at the expense of transcript fidelity via force and allostery

In prokaryotes, translation can occur on mRNA that is being transcribed in a process called coupling. How the ribosome affects the RNA polymerase (RNAP) during coupling is not well understood. Here, we reconstituted the E. coli coupling system and demonstrated that the ribosome can prevent pausing and termination of RNAP and double the overall transcription rate at the expense of fidelity. Moreover, we monitored single RNAPs coupled to ribosomes and show that coupling increases the pause-free velocity of the polymerase and that a mechanical assisting force is sufficient to explain the majority of the effects of coupling. Also, by cryo-EM, we observed that RNAPs with a terminal mismatch adopt a backtracked conformation, while a coupled ribosome allosterically induces these polymerases toward a catalytically active anti-swiveled state. Finally, we demonstrate that prolonged RNAP pausing is detrimental to cell viability, which could be prevented by polymerase reactivation through a coupled ribosome.

Observation of coordinated cotranscriptional RNA folding events

RNA begins to fold as it is transcribed by an RNA polymerase. Consequently, RNA folding is constrained by the direction and rate of transcription. Understanding how RNA folds into secondary and tertiary structures therefore requires methods for determining the structure of cotranscriptional folding intermediates. Cotranscriptional RNA chemical probing methods accomplish this by systematically probing the structure of nascent RNA that is displayed from RNA polymerase. Here, we have developed a concise, high-resolution cotranscriptional RNA chemical probing procedure called Transcription Elongation Complex RNA structure probing-Multilength (TECprobe-ML). We validated TECprobe-ML by replicating and extending previous analyses of ZTP and fluoride riboswitch folding, and mapped the folding pathway of a ppGpp-sensing riboswitch. In each system, TECprobe-ML identified coordinated cotranscriptional folding events that mediate transcription antitermination. Our findings establish TECprobe-ML as an accessible method for mapping cotranscriptional RNA folding pathways.

Integrating experimental data with molecular simulations to investigate RNA structural dynamics

Conformational dynamics is crucial for ribonucleic acid (RNA) function. Techniques such as nuclear magnetic resonance, cryo-electron microscopy, small- and wide-angle X-ray scattering, chemical probing, single-molecule Förster resonance energy transfer, or even thermal or mechanical denaturation experiments probe RNA dynamics at different time and space resolutions. Their combination with accurate atomistic molecular dynamics (MD) simulations paves the way for quantitative and detailed studies of RNA dynamics. First, experiments provide a quantitative validation tool for MD simulations. Second, available data can be used to refine simulated structural ensembles to match experiments. Finally, comparison with experiments allows for improving MD force fields that are transferable to new systems for which data is not available. Here we review the recent literature and provide our perspective on this field.

Translational regulation by Hfq-Crc assemblies emerges from polymorphic ribonucleoprotein folding

The widely occurring bacterial RNA chaperone Hfq is a key factor in the post-transcriptional control of hundreds of genes in Pseudomonas aeruginosa. How this broadly acting protein can contribute to the regulatory requirements of many different genes remains puzzling. Here, we describe cryo-EM structures of higher order assemblies formed by Hfq and its partner protein Crc on control regions of different P. aeruginosa target mRNAs. Our results show that these assemblies have mRNA-specific quaternary architectures resulting from the combination of multivalent protein-protein interfaces and recognition of patterns in the RNA sequence. The structural polymorphism of these ribonucleoprotein assemblies enables selective translational repression of many different target mRNAs. This system elucidates how highly complex regulatory pathways can evolve with a minimal economy of proteinogenic components in combination with RNA sequence and fold.

DrTransformer: heuristic cotranscriptional RNA folding using the nearest neighbor energy model

MOTIVATION

Folding during transcription can have an important influence on the structure and function of RNA molecules, as regions closer to the 5' end can fold into metastable structures before potentially stronger interactions with the 3' end become available. Thermodynamic RNA folding models are not suitable to predict structures that result from cotranscriptional folding, as they can only calculate properties of the equilibrium distribution. Other software packages that simulate the kinetic process of RNA folding during transcription exist, but they are mostly applicable for short sequences.

RESULTS

We present a new algorithm that tracks changes to the RNA secondary structure ensemble during transcription. At every transcription step, new representative local minima are identified, a neighborhood relation is defined and transition rates are estimated for kinetic simulations. After every simulation, a part of the ensemble is removed and the remainder is used to search for new representative structures. The presented algorithm is deterministic (up to numeric instabilities of simulations), fast (in comparison with existing methods), and it is capable of folding RNAs much longer than 200 nucleotides.

AVAILABILITY AND IMPLEMENTATION

This software is open-source and available at https://github.com/ViennaRNA/drtransformer.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.

Discovering functional motifs in long noncoding RNAs

Long noncoding RNAs (lncRNAs) are products of pervasive transcription that closely resemble messenger RNAs on the molecular level, yet function through largely unknown modes of action. The current model is that the function of lncRNAs often relies on specific, typically short, conserved elements, connected by linkers in which specific sequences and/or structures are less important. This notion has fueled the development of both computational and experimental methods focused on the discovery of functional elements within lncRNA genes, based on diverse signals such as evolutionary conservation, predicted structural elements, or the ability to rescue loss-of-function phenotypes. In this review, we outline the main challenges that the different methods need to overcome, describe the recently developed approaches, and discuss their respective limitations. This article is categorized under: RNA Evolution and Genomics > Computational Analyses of RNA RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs.

The dynamic, motile and deformative properties of RNA nanoparticles facilitate the third milestone of drug development

Besides mRNA, rRNA, and tRNA, cells contain many other noncoding RNA that display critical roles in the regulation of cellular functions. Human genome sequencing revealed that the majority of non-protein-coding DNA actually codes for non-coding RNAs. The dynamic nature of RNA results in its motile and deformative behavior. These conformational transitions such as the change of base-pairing, breathing within complemented strands, and pseudoknot formation at the 2D level as well as the induced-fit and conformational capture at the 3D level are important for their biological functions including regulation, translation, and catalysis. The dynamic, motile and catalytic activity has led to a belief that RNA is the origin of life. We have recently reported that the deformative property of RNA nanoparticles enhances their penetration through the leaky blood vessel of cancers which leads to highly efficient tumor accumulation. This special deformative property also enables RNA nanoparticles to pass the glomerulus, overcoming the filtration size limit, resulting in fast renal excretion and rapid body clearance, thus low or no toxicity. The biodistribution of RNA nanoparticles can be further improved by the incorporation of ligands for cancer targeting. In addition to the favorable biodistribution profiles, RNA nanoparticles possess other properties including self-assembly, negative charge, programmability, and multivalency; making it a great material for pharmaceutical applications. The intrinsic negative charge of RNA nanoparticles decreases the toxicity of drugs by preventing nonspecific binding to the negative charged cell membrane and enhancing the solubility of hydrophobic drugs. The polyvalent property of RNA nanoparticles allows the multi-functionalization which can apply to overcome drug resistance. This review focuses on the summary of these unique properties of RNA nanoparticles, which describes the mechanism of RNA dynamic, motile and deformative properties, and elucidates and prepares to welcome the RNA therapeutics as the third milestone in pharmaceutical drug development.

How does RNA fold dynamically?

Recent advances in interrogating RNA folding dynamics have shown the classical model of RNA folding to be incomplete. Here, we pose three prominent questions for the field that are at the forefront of our understanding of the importance of RNA folding dynamics for RNA function. The first centers on the most appropriate biophysical framework to describe changes to the RNA folding energy landscape that a growing RNA chain encounters during transcriptional elongation. The second focuses on the potential ubiquity of strand displacement - a process by which RNA can rapidly change conformations - and how this process may be generally present in broad classes of seemingly different RNAs. The third raises questions about the potential importance and roles of cellular protein factors in RNA conformational switching. Answers to these questions will greatly improve our fundamental knowledge of RNA folding and function, drive biotechnological advances that utilize engineered RNAs, and potentially point to new areas of biology yet to be discovered.

Computational approaches for RNA structure ensemble deconvolution from structure probing data

RNA structure probing experiments have emerged over the last decade as a straightforward way to determine the structure of RNA molecules in a number of different contexts. Although powerful, the ability of RNA to dynamically interconvert between, and to simultaneously populate, alternative structural configurations, poses a nontrivial challenge to the interpretation of data derived from these experiments. Recent efforts aimed at developing computational methods for the reconstruction of coexisting alternative RNA conformations from structure probing data are paving the way to the study of RNA structure ensembles, even in the context of living cells. In this review, we critically discuss these methods, their limitations and possible future improvements.

Global mapping of RNA homodimers in living cells

RNA homodimerization is important for various physiological processes, including the assembly of membraneless organelles, RNA subcellular localization, and packaging of viral genomes. However, understanding RNA dimerization has been hampered by the lack of systematic in vivo detection methods. Here, we show that CLASH, PARIS, and other RNA proximity ligation methods detect RNA homodimers transcriptome-wide as "overlapping" chimeric reads that contain more than one copy of the same sequence. Analyzing published proximity ligation data sets, we show that RNA:RNA homodimers mediated by direct base-pairing are rare across the human transcriptome, but highly enriched in specific transcripts, including U8 snoRNA, U2 snRNA, and a subset of tRNAs. Mutations in the homodimerization domain of U8 snoRNA impede dimerization in vitro and disrupt zebrafish development in vivo, suggesting an evolutionarily conserved role of this domain. Analysis of virus-infected cells reveals homodimerization of SARS-CoV-2 and Zika genomes, mediated by specific palindromic sequences located within protein-coding regions of N gene in SARS-CoV-2 and NS2A gene in Zika. We speculate that regions of viral genomes involved in homodimerization may constitute effective targets for antiviral therapies.

Cotranscriptional RNA strand exchange underlies the gene regulation mechanism in a purine-sensing transcriptional riboswitch

RNA folds cotranscriptionally to traverse out-of-equilibrium intermediate structures that are important for RNA function in the context of gene regulation. To investigate this process, here we study the structure and function of the Bacillus subtilis yxjA purine riboswitch, a transcriptional riboswitch that downregulates a nucleoside transporter in response to binding guanine. Although the aptamer and expression platform domain sequences of the yxjA riboswitch do not completely overlap, we hypothesized that a strand exchange process triggers its structural switching in response to ligand binding. In vivo fluorescence assays, structural chemical probing data and experimentally informed secondary structure modeling suggest the presence of a nascent intermediate central helix. The formation of this central helix in the absence of ligand appears to compete with both the aptamer's P1 helix and the expression platform's transcriptional terminator. All-atom molecular dynamics simulations support the hypothesis that ligand binding stabilizes the aptamer P1 helix against central helix strand invasion, thus allowing the terminator to form. These results present a potential model mechanism to explain how ligand binding can induce downstream conformational changes by influencing local strand displacement processes of intermediate folds that could be at play in multiple riboswitch classes.

Structural basis for impaired 5' processing of a mutant tRNA associated with defects in neuronal homeostasis

SignificanceUnderstanding and treating neurological disorders are global priorities. Some of these diseases are engendered by mutations that cause defects in the cellular synthesis of transfer RNAs (tRNAs), which function as adapter molecules that translate messenger RNAs into proteins. During tRNA biogenesis, ribonuclease P catalyzes removal of the transcribed sequence upstream of the mature tRNA. Here, we focus on a cytoplasmic tRNAArgUCU that is expressed specifically in neurons and, when harboring a particular point mutation, contributes to neurodegeneration in mice. Our results suggest that this mutation favors stable alternative structures that are not cleaved by mouse ribonuclease P and motivate a paradigm that may help to understand the molecular basis for disease-associated mutations in other tRNAs.

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.

Cotranscriptional RNA Chemical Probing

Cotranscriptional folding is a fundamental step in RNA biogenesis and the basis for many RNA-mediated gene regulation systems. Understanding how RNA folds as it is synthesized requires experimental methods that can systematically identify intermediate RNA structures that form during transcription. Cotranscriptional RNA chemical probing experiments achieve this by applying high-throughput RNA structure probing to an in vitro transcribed array of cotranscriptionally folded intermediate transcripts. In this chapter, we present guidelines and procedures for integrating single-round in vitro transcription using E. coli RNA polymerase with high-throughput RNA chemical probing workflows. We provide an overview of key concepts including DNA template design, transcription roadblocking strategies, single-round in vitro transcription with E. coli RNA polymerase, and RNA chemical probing and describe procedures for DNA template preparation, cotranscriptional RNA chemical probing, RNA purification, and 3' adapter ligation. The end result of these procedures is a purified RNA library that can be prepared for Illumina sequencing using established high-throughput RNA structure probing library construction strategies.

Computer-aided comprehensive explorations of RNA structural polymorphism through complementary simulation methods

While RNA folding was originally seen as a simple problem to solve, it has been shown that the promiscuous interactions of the nucleobases result in structural polymorphism, with several competing structures generally observed for non-coding RNA. This inherent complexity limits our understanding of these molecules from experiments alone, and computational methods are commonly used to study RNA. Here, we discuss three advanced sampling schemes, namely Hamiltonian-replica exchange molecular dynamics (MD), ratchet-and-pawl MD and discrete path sampling, as well as the HiRE-RNA coarse-graining scheme, and highlight how these approaches are complementary with reference to recent case studies. While all computational methods have their shortcomings, the plurality of simulation methods leads to a better understanding of experimental findings and can inform and guide experimental work on RNA polymorphism.

Detection of pks Island mRNAs Using Toehold Sensors in Escherichia coli

Synthetic biologists have applied biomolecular engineering approaches toward the goal of novel biological devices and have shown progress in diverse areas of medicine and biotechnology. Especially promising is the application of synthetic biological devices towards a novel class of molecular diagnostics. As an example, a de-novo-designed riboregulator called toehold switch, with its programmability and compatibility with field-deployable devices showed promising in vitro applications for viral RNA detection such as Zika and Corona viruses. However, the in vivo application of high-performance RNA sensors remains challenging due to the secondary structure of long mRNA species. Here, we introduced ‘Helper RNAs’ that can enhance the functionality of toehold switch sensors by mitigating the effect of secondary structures around a target site. By employing the helper RNAs, previously reported mCherry mRNA sensor showed improved fold-changes in vivo. To further generalize the Helper RNA approaches, we employed automatic design pipeline for toehold sensors that target the essential genes within the pks island, an important target of biomedical research in connection with colorectal cancer. The toehold switch sensors showed fold-changes upon the expression of full-length mRNAs that apparently depended sensitively on the identity of the gene as well as the predicted local structure within the target region of the mRNA. Still, the helper RNAs could improve the performance of toehold switch sensors in many instances, with up to 10-fold improvement over no helper cases. These results suggest that the helper RNA approaches can further assist the design of functional RNA devices in vivo with the aid of the streamlined automatic design software developed here. Further, our solutions for screening and stabilizing single-stranded region of mRNA may find use in other in vivo mRNA-sensing applications such as cas13 crRNA design, transcriptome engineering, and trans-cleaving ribozymes.

Kinetics of RNA and RNA:DNA Hybrid Strand Displacement

In nucleic acid nanotechnology, strand displacement is a widely used mechanism where one strand from a hybridized duplex is exchanged with an invading strand that binds to a toehold, a single-stranded region on the duplex. It is used to perform logic operations on a molecular level, initiate cascaded reactions, or even for in vivo diagnostics and treatments. While systematic experimental studies have been carried out to probe the kinetics of strand displacement in DNA with different toehold lengths, sequences, and mismatch positions, there has not been a comparable investigation of RNA or RNA-DNA hybrid systems. Here, we experimentally study how toehold length, toehold location (5' or 3' end of the strand), and mismatches influence the strand displacement kinetics. We observe reaction acceleration with increasing toehold length and placement of the toehold at the 5' end of the substrate. We find that mismatches closer to the interface of toehold and duplex slow down the reaction more than remote mismatches. A comparison of RNA and DNA displacement with hybrid displacement (RNA invading DNA or DNA invading RNA) is partly explainable by the thermodynamic stabilities of the respective toehold regions, but also suggests that the rearrangement from B-form to A-form helix in the case of RNA invading DNA might play a role in the kinetics.

Transcription complexes as RNA chaperones

To exert their functions, RNAs adopt diverse structures, ranging from simple secondary to complex tertiary and quaternary folds. In vivo, RNA folding starts with RNA transcription, and a wide variety of processes are coupled to co-transcriptional RNA folding events, including the regulation of fundamental transcription dynamics, gene regulation by mechanisms like attenuation, RNA processing or ribonucleoprotein particle formation. While co-transcriptional RNA folding and associated co-transcriptional processes are by now well accepted as pervasive regulatory principles in all organisms, investigations into the role of the transcription machinery in co-transcriptional folding processes have so far largely focused on effects of the order in which RNA regions are produced and of transcription kinetics. Recent structural and structure-guided functional analyses of bacterial transcription complexes increasingly point to an additional role of RNA polymerase and associated transcription factors in supporting co-transcriptional RNA folding by fostering or preventing strategic contacts to the nascent transcripts. In general, the results support the view that transcription complexes can act as RNA chaperones, a function that has been suggested over 30 years ago. Here, we discuss transcription complexes as RNA chaperones based on recent examples from bacterial transcription.

The many faces of RNA-based RNase P, an RNA-world relic

RNase P is an essential enzyme that catalyzes removal of the 5' leader from precursor transfer RNAs. The ribonucleoprotein (RNP) form of RNase P is present in all domains of life and comprises a single catalytic RNA (ribozyme) and a variable number of protein cofactors. Recent cryo-electron microscopy structures of representative archaeal and eukaryotic (nuclear) RNase P holoenzymes bound to tRNA substrate/product provide high-resolution detail on subunit organization, topology, and substrate recognition in these large, multisubunit catalytic RNPs. These structures point to the challenges in understanding how proteins modulate the RNA functional repertoire and how the structure of an ancient RNA-based catalyst was reshaped during evolution by new macromolecular associations that were likely necessitated by functional/regulatory coupling.

Conformational Ensembles of Noncoding Elements in the SARS-CoV-2 Genome from Molecular Dynamics Simulations

The 5' untranslated region (UTR) of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome is a conserved, functional and structured genomic region consisting of several RNA stem-loop elements. While the secondary structure of such elements has been determined experimentally, their three-dimensional structures are not known yet. Here, we predict structure and dynamics of five RNA stem loops in the 5'-UTR of SARS-CoV-2 by extensive atomistic molecular dynamics simulations, more than 0.5 ms of aggregate simulation time, in combination with enhanced sampling techniques. We compare simulations with available experimental data, describe the resulting conformational ensembles, and identify the presence of specific structural rearrangements in apical and internal loops that may be functionally relevant. Our atomic-detailed structural predictions reveal a rich dynamics in these RNA molecules, could help the experimental characterization of these systems, and provide putative three-dimensional models for structure-based drug design studies.

Transcription Regulation Through Nascent RNA Folding

Folding of RNA into secondary structures through intramolecular base pairing determines an RNA's three-dimensional architecture and associated function. Simple RNA structures like stem loops can provide specialized functions independent of coding capacity, such as protein binding, regulation of RNA processing and stability, stimulation or inhibition of translation. RNA catalysis is dependent on tertiary structures found in the ribosome, tRNAs and group I and II introns. While the extent to which non-coding RNAs contribute to cellular maintenance is generally appreciated, the fact that both non-coding and coding RNA can assume relevant structural states has only recently gained attention. In particular, the co-transcriptional folding of nascent RNA of all classes has the potential to regulate co-transcriptional processing, RNP (ribonucleoprotein particle) formation, and transcription itself. Riboswitches are established examples of co-transcriptionally folded coding RNAs that directly regulate transcription, mainly in prokaryotes. Here we discuss recent studies in both prokaryotes and eukaryotes showing that structure formation may carry a more widespread regulatory logic during RNA synthesis. Local structures forming close to the catalytic center of RNA polymerases have the potential to regulate transcription by reducing backtracking. In addition, stem loops or more complex structures may alter co-transcriptional RNA processing or its efficiency. Several examples of functional structures have been identified to date, and this review provides an overview of physiologically distinct processes where co-transcriptionally folded RNA plays a role. Experimental approaches such as single-molecule FRET and in vivo structural probing to further advance our insight into the significance of co-transcriptional structure formation are discussed.