RNA folding during transcription: protocols and studies

Co-transcriptional folding of the glmS ribozyme enables a rapid response to metabolite

The glmS ribozyme riboswitch, located in the 5' untranslated region of the Bacillus subtilis glmS messenger RNA (mRNA), regulates cell wall biosynthesis through ligand-induced self-cleavage and decay of the glmS mRNA. Although self-cleavage of the refolded glmS ribozyme has been studied extensively, it is not known how early the ribozyme folds and self-cleaves during transcription. Here, we combine single-molecule fluorescence with kinetic modeling to show that self-cleavage can occur during transcription before the ribozyme is fully synthesized. Moreover, co-transcriptional folding of the RNA at a physiological elongation rate allows the ribozyme catalytic core to react without the downstream peripheral stability domain. Dimethyl sulfate footprinting further revealed how slow sequential folding favors formation of the native core structure through fraying of misfolded helices and nucleation of a native pseudoknot. Ribozyme self-cleavage at an early stage of transcription may benefit glmS regulation in B. subtilis, as it exposes the mRNA to exoribonuclease before translation of the open reading frame can begin. Our results emphasize the importance of co-transcriptional folding of RNA tertiary structure for cis-regulation of mRNA stability.

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.

Regulation of Gene Expression Through Effector-dependent Conformational Switching by Cobalamin Riboswitches

Riboswitches are an outstanding example of genetic regulation mediated by RNA conformational switching. In these non-coding RNA elements, the occupancy status of a ligand-binding domain governs the mRNA's decision to form one of two mutually exclusive structures in the downstream expression platform. Temporal constraints upon the function of many riboswitches, requiring folding of complex architectures and conformational switching in a limited co-transcriptional timeframe, make them ideal model systems for studying these processes. In this review, we focus on the mechanism of ligand-directed conformational changes in one of the most widely distributed riboswitches in bacteria: the cobalamin family. We describe the architectural features of cobalamin riboswitches whose structures have been determined by x-ray crystallography, which suggest a direct physical role of cobalamin in effecting the regulatory switch. Next, we discuss a series of experimental approaches applied to several model cobalamin riboswitches that interrogate these structural models. As folding is central to riboswitch function, we consider the differences in folding landscapes experienced by RNAs that are produced in vitro and those that are allowed to fold co-transcriptionally. Finally, we highlight a set of studies that reveal the difficulties of studying cobalamin riboswitches outside the context of transcription and that co-transcriptional approaches are essential for developing a more accurate picture of their structure-function relationships in these switches. This understanding will be essential for future advancements in the use of small-molecule guided RNA switches in a range of applications such as biosensors, RNA imaging tools, and nucleic acid-based therapies.

Regulation of the ThiM riboswitch is facilitated by the trapped structure formed during transcription of the wild-type sequence

The ThiM riboswitch from Escherichia coli is a typical mRNA device that modulates downstream gene expression by sensing TPP. The helix-based RNA folding theory is used to investigate its detailed regulatory behaviors in cells. This RNA molecule is transcriptionally trapped in a state with the unstructured SD sequence in the absence of TPP, which induces downstream gene expression. As a key step to turn on gene expression, formation of this trapped state (the genetic ON state) highly depends on the co-transcriptional folding of its wild-type sequence. Instead of stabilities of the genetic ON and OFF states, the transcription rate, pause, and ligand levels are combined to affect the ThiM riboswitch-mediated gene regulation, which is consistent with a kinetic control model.

Role of a hairpin-stabilized pause in the Escherichia coli thiC riboswitch function

Transcriptional pauses have been reported in bacterial riboswitches and, in some cases, their specific positioning has been shown to be important for gene regulation. Here, we show that a hairpin structure in the Escherichia coli thiamin pyrophosphate (TPP) thiC riboswitch is involved in transcriptional pausing and ligand sensitivity. Using in vitro transcription kinetic experiments, we show that all three major transcriptional pauses in the thiC riboswitch are affected by NusA, a transcriptional factor known to stimulate hairpin-stabilized pauses. Using a truncated region of the riboswitch, we isolated the hairpin structure responsible for stabilization of the most upstream pause. Destabilization of this structure led to a weaker pause and a decreased NusA effect. In the context of the full-length riboswitch, this same mutation also led to a weaker pause, as well as a decreased TPP binding affinity. Our work suggests that RNA structures involved in transcriptional pausing in riboswitches are important for ligand sensitivity, most likely by increasing the time allowed to the ligand for binding to the riboswitch.

Viroid research and its significance for RNA technology and basic biochemistry

Viroids were described 47 years ago as the smallest RNA molecules capable of infecting plants and autonomously self-replicating without an encoded protein. Work on viroids initiated the development of a number of innovative methods. Novel chromatographic and gelelectrophoretic methods were developed for the purification and characterization of viroids; these methods were later used in molecular biology, gene technology and in prion research. Theoretical and experimental studies of RNA folding demonstrated the general biological importance of metastable structures, and nuclear magnetic resonance spectroscopy of viroid RNA showed the partially covalent nature of hydrogen bonds in biological macromolecules. RNA biochemistry and molecular biology profited from viroid research, such as in the detection of RNA as template of DNA-dependent polymerases and in mechanisms of gene silencing. Viroids, the first circular RNA detected in nature, are important for studies on the much wider spectrum of circular RNAs and other non-coding RNAs.

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.

Computational Methods for Modeling Aptamers and Designing Riboswitches

Riboswitches, which are located within certain noncoding RNA region perform functions as genetic "switches", regulating when and where genes are expressed in response to certain ligands. Understanding the numerous functions of riboswitches requires computation models to predict structures and structural changes of the aptamer domains. Although aptamers often form a complex structure, computational approaches, such as RNAComposer and Rosetta, have already been applied to model the tertiary (three-dimensional (3D)) structure for several aptamers. As structural changes in aptamers must be achieved within the certain time window for effective regulation, kinetics is another key point for understanding aptamer function in riboswitch-mediated gene regulation. The coarse-grained self-organized polymer (SOP) model using Langevin dynamics simulation has been successfully developed to investigate folding kinetics of aptamers, while their co-transcriptional folding kinetics can be modeled by the helix-based computational method and BarMap approach. Based on the known aptamers, the web server Riboswitch Calculator and other theoretical methods provide a new tool to design synthetic riboswitches. This review will represent an overview of these computational methods for modeling structure and kinetics of riboswitch aptamers and for designing riboswitches.

Distributed biotin-streptavidin transcription roadblocks for mapping cotranscriptional RNA folding

RNA folding during transcription directs an order of folding that can determine RNA structure and function. However, the experimental study of cotranscriptional RNA folding has been limited by the lack of easily approachable methods that can interrogate nascent RNA structure at nucleotide resolution. To address this, we previously developed cotranscriptional selective 2΄-hydroxyl acylation analyzed by primer extension sequencing (SHAPE-Seq) to simultaneously probe all intermediate RNA transcripts during transcription by stalling elongation complexes at catalytically dead EcoRIE111Q roadblocks. While effective, the distribution of elongation complexes using EcoRIE111Q requires laborious PCR using many different oligonucleotides for each sequence analyzed. Here, we improve the broad applicability of cotranscriptional SHAPE-Seq by developing a sequence-independent biotin-streptavidin (SAv) roadblocking strategy that simplifies the preparation of roadblocking DNA templates. We first determine the properties of biotin-SAv roadblocks. We then show that randomly distributed biotin-SAv roadblocks can be used in cotranscriptional SHAPE-Seq experiments to identify the same RNA structural transitions related to a riboswitch decision-making process that we previously identified using EcoRIE111Q. Lastly, we find that EcoRIE111Q maps nascent RNA structure to specific transcript lengths more precisely than biotin-SAv and propose guidelines to leverage the complementary strengths of each transcription roadblock in cotranscriptional SHAPE-Seq.

The determinants of alternative RNA splicing in human cells

Alternative splicing represents an important level of the regulation of gene function in eukaryotic organisms. It plays a critical role in virtually every biological process within an organism, including regulation of cell division and cell death, differentiation of tissues in the embryo and the adult organism, as well as in cellular response to diverse environmental factors. In turn, studies of the last decade have shown that alternative splicing itself is controlled by different mechanisms. Unfortunately, there is no clear understanding of how these diverse mechanisms, or determinants, regulate and constrain the set of alternative RNA species produced from any particular gene in every cell of the human body. Here, we provide a consolidated overview of alternative splicing determinants including RNA-protein interactions, epigenetic regulation via chromatin remodeling, coupling of transcription-to-alternative splicing, effect of secondary structures in pre-RNA, and function of the RNA quality control systems. We also extensively and critically discuss some mechanistic insights on coordinated inclusion/exclusion of exons during the formation of mature RNA molecules. We conclude that the final structure of RNA is pre-determined by a complex interplay between cis- and trans-acting factors. Altogether, currently available empirical data significantly expand our understanding of the functioning of the alternative splicing machinery of cells in normal and pathological conditions. On the other hand, there are still many blind spots that require further deep investigations.

Cotranscriptional folding of a riboswitch at nucleotide resolution

RNAs can begin to fold immediately after emerging from RNA polymerase during transcription. Interactions between nascent RNAs and ligands during cotranscriptional folding can direct the formation of alternative RNA structures, a feature exploited by non-coding RNAs called riboswitches to make gene regulatory decisions. Despite their importance, cotranscriptional folding pathways have yet to be uncovered with sufficient resolution to reveal how cotranscriptional folding governs RNA structure and function. To access cotranscriptional folding at nucleotide resolution, we extend selective 2’-hydroxyl acylation analyzed by primer extension sequencing (SHAPE-Seq) to measure structural information of nascent RNAs during transcription. With cotranscriptional SHAPE-Seq, we determine how the B. cereus crcB fluoride riboswitch cotranscriptional folding pathway undergoes a ligand-dependent bifurcation that delays or promotes terminator formation via a series of coordinated structural transitions. Our results directly link cotranscriptional RNA folding to a genetic decision and establish a framework for cotranscriptional analysis of RNA structure at nucleotide resolution.

Conformational Rearrangements of Individual Nucleotides during RNA-Ligand Binding Are Rate-Differentiated

A pronounced rate differentiation has been found for conformational rearrangements of individual nucleobases that occur during ligand recognition of the preQ1 class-I riboswitch aptamer from Thermoanaerobacter tengcongensis. Rate measurements rely on the 2ApFold approach by analyzing the fluorescence response of riboswitch variants, each with a single, strategically positioned 2-aminopurine nucleobase substitution. Observed rate discrimination between the fastest and the slowest conformational adaption is 22-fold, with the largest rate observed for the rearrangement of a nucleoside directly at the binding site and the smallest rate observed for the 3'-unpaired nucleoside that stacks onto the pseudo-knot-closing Watson-Crick base pair. Our findings provide novel insights into how compact, prefolded RNAs that follow the induced-fit recognition mechanism adapt local structural elements in response to ligand binding on a rather broad time scale and how this process culminates in a structural signal that is responsible for efficient downregulation of ribosomal translation.

Quantitative and predictive model of kinetic regulation by E. coli TPP riboswitches

Riboswitches are non-coding elements upstream or downstream of mRNAs that, upon binding of a specific ligand, regulate transcription and/or translation initiation in bacteria, or alternative splicing in plants and fungi. We have studied thiamine pyrophosphate (TPP) riboswitches regulating translation of thiM operon and transcription and translation of thiC operon in E. coli, and that of THIC in the plant A. thaliana. For all, we ascertained an induced-fit mechanism involving initial binding of the TPP followed by a conformational change leading to a higher-affinity complex. The experimental values obtained for all kinetic and thermodynamic parameters of TPP binding imply that the regulation by A. thaliana riboswitch is governed by mass-action law, whereas it is of kinetic nature for the two bacterial riboswitches. Kinetic regulation requires that the RNA polymerase pauses after synthesis of each riboswitch aptamer to leave time for TPP binding, but only when its concentration is sufficient. A quantitative model of regulation highlighted how the pausing time has to be linked to the kinetic rates of initial TPP binding to obtain an ON/OFF switch in the correct concentration range of TPP. We verified the existence of these pauses and the model prediction on their duration. Our analysis also led to quantitative estimates of the respective efficiency of kinetic and thermodynamic regulations, which shows that kinetically regulated riboswitches react more sharply to concentration variation of their ligand than thermodynamically regulated riboswitches. This rationalizes the interest of kinetic regulation and confirms empirical observations that were obtained by numerical simulations.

A Two-Way Street: Regulatory Interplay between RNA Polymerase and Nascent RNA Structure

The vectorial (5'-to-3' at varying velocity) synthesis of RNA by cellular RNA polymerases (RNAPs) creates a rugged kinetic landscape, demarcated by frequent, sometimes long-lived, pauses. In addition to myriad gene-regulatory roles, these pauses temporally and spatially program the co-transcriptional, hierarchical folding of biologically active RNAs. Conversely, these RNA structures, which form inside or near the RNA exit channel, interact with the polymerase and adjacent protein factors to influence RNA synthesis by modulating pausing, termination, antitermination, and slippage. Here, we review the evolutionary origin, mechanistic underpinnings, and regulatory consequences of this interplay between RNAP and nascent RNA structure. We categorize and rationalize the extensive linkage between the transcriptional machinery and its product, and provide a framework for future studies.

Traditional Chemical Mapping of RNA Structure In Vitro and In Vivo

Chemical probing is often used to gain knowledge on the secondary and tertiary structures of RNA molecules either free or engaged in complexes with ligands. The method monitors the reactivity of each nucleotide towards chemicals of various specificities reflecting the hydrogen bonding environment of each nucleotide within the RNA molecule. In addition, information can be obtained on the binding site of a ligand (noncoding RNAs, protein, metabolites), and on RNA conformational changes that accompanied ligand binding or perturbation of the environmental cues. The detection of the modifications can be obtained either by using end-labeled RNA molecules or by primer extension using reverse transcriptase. The goal of this chapter is to provide the reader with an experimental guide to probe the structure of RNA in vitro and in vivo with the most suitable chemical probes.

Progress and challenges for chemical probing of RNA structure inside living cells

Proper gene expression is essential for the survival of every cell. Once thought to be a passive transporter of genetic information, RNA has recently emerged as a key player in nearly every pathway in the cell. A full description of its structure is critical to understanding RNA function. Decades of research have focused on utilizing chemical tools to interrogate the structures of RNAs, with recent focus shifting to performing experiments inside living cells. This Review will detail the design and utility of chemical reagents used in RNA structure probing. We also outline how these reagents have been used to gain a deeper understanding of RNA structure in vivo. We review the recent merger of chemical probing with deep sequencing. Finally, we outline some of the hurdles that remain in fully characterizing the structure of RNA inside living cells, and how chemical biology can uniquely tackle such challenges.

Trigger-helix folding pathway and SI3 mediate catalysis and hairpin-stabilized pausing by Escherichia coli RNA polymerase

The conformational dynamics of the polymorphous trigger loop (TL) in RNA polymerase (RNAP) underlie multiple steps in the nucleotide addition cycle and diverse regulatory mechanisms. These mechanisms include nascent RNA hairpin-stabilized pausing, which inhibits TL folding into the trigger helices (TH) required for rapid nucleotide addition. The nascent RNA pause hairpin forms in the RNA exit channel and promotes opening of the RNAP clamp domain, which in turn stabilizes a partially folded, paused TL conformation that disfavors TH formation. We report that inhibiting TH unfolding with a disulfide crosslink slowed multiround nucleotide addition only modestly but eliminated hairpin-stabilized pausing. Conversely, a substitution that disrupts the TH folding pathway and uncouples establishment of key TH–NTP contacts from complete TH formation and clamp movement allowed rapid catalysis and eliminated hairpin-stabilized pausing. We also report that the active-site distal arm of the TH aids TL folding, but that a 188-aa insertion in the Escherichia coli TL (sequence insertion 3; SI3) disfavors TH formation and stimulates pausing. The effect of SI3 depends on the jaw domain, but not on downstream duplex DNA. Our results support the view that both SI3 and the pause hairpin modulate TL folding in a constrained pathway of intermediate states.

RNA polymerase pausing and nascent-RNA structure formation are linked through clamp-domain movement

The rates of RNA synthesis and the folding of nascent RNA into biologically active structures are linked via pausing by RNA polymerase (RNAP). Structures that form within the RNA-exit channel can either increase pausing by interacting with RNAP or decrease pausing by preventing backtracking. Conversely, pausing is required for proper folding of some RNAs. Opening of the RNAP clamp domain has been proposed to mediate some effects of nascent-RNA structures. However, the connections among RNA structure formation and RNAP clamp movement and catalytic activity remain uncertain. Here, we assayed exit-channel structure formation in Escherichia coli RNAP with disulfide cross-links that favor closed- or open-clamp conformations and found that clamp position directly influences RNA structure formation and RNAP catalytic activity. We report that exit-channel RNA structures slow pause escape by favoring clamp opening through interactions with the flap that slow translocation.

A personal perspective on chemistry-driven RNA research

In this mini review, we discuss how our understanding of ribonucleic acid (RNA) properties becomes significantly deepened when a broad range of modern chemical and biophysical methods is applied. We span our perspective from RNA solid-phase synthesis and site-specific labeling to single-molecule fluorescence-resonance-energy-transfer imaging and NMR spectroscopy approaches to explore the dynamics of RNA over a broad timescale. We then move on to Fourier-transform-ion-cyclotron-resonance mass spectrometry (FT-ICR-MS) as a powerful technique for RNA sequencing and modification analysis. The novel methodological developments are discussed for selected biological systems that include the thiamine-pyrophosphate riboswitch, HIV and ribosomal A-site RNA, and transfer RNA.

Changes in transcriptional pausing modify the folding dynamics of the pH-responsive RNA element

Previously, we described a novel pH-responsive RNA element in Escherichia coli that resides in the 5′ untranslated region of the alx gene and controls its translation in a pH-dependent manner. Under normal growth conditions, this RNA region forms a translationally inactive structure, but when transcribed under alkaline conditions, it forms an active structure producing the Alx protein. We identified two distinct transcriptional pause sites and proposed that pausing at these sites interfered with the formation of the inactive structure while facilitating folding of the active one. Alkali increases the longevity of pausing at these sites, thereby promoting folding of the translationally active form of alx RNA. We show here that mutations that modify the extent and/or position of pausing, although silent with regard to structure stability per se, greatly influence the dynamics of folding and thereby translation. Our data illustrate the mechanistic design of alx regulation, relying on precise temporal and spatial characteristics. We propose that this unique design provides an opportunity for environmental signals such as pH to introduce structural changes in the RNA and thereby modulate expression.

Folding and ligand recognition of the TPP riboswitch aptamer at single-molecule resolution

Thiamine pyrophosphate (TPP)-sensitive mRNA domains are the most prevalent riboswitches known. Despite intensive investigation, the complex ligand recognition and concomitant folding processes in the TPP riboswitch that culminate in the regulation of gene expression remain elusive. Here, we used single-molecule fluorescence resonance energy transfer imaging to probe the folding landscape of the TPP aptamer domain in the absence and presence of magnesium and TPP. To do so, distinct labeling patterns were used to sense the dynamics of the switch helix (P1) and the two sensor arms (P2/P3 and P4/P5) of the aptamer domain. The latter structural elements make interdomain tertiary contacts (L5/P3) that span a region immediately adjacent to the ligand-binding site. In each instance, conformational dynamics of the TPP riboswitch were influenced by ligand binding. The P1 switch helix, formed by the 5' and 3' ends of the aptamer domain, adopts a predominantly folded structure in the presence of Mg(2+) alone. However, even at saturating concentrations of Mg(2+) and TPP, the P1 helix, as well as distal regions surrounding the TPP-binding site, exhibit an unexpected degree of residual dynamics and disperse kinetic behaviors. Such plasticity results in a persistent exchange of the P3/P5 forearms between open and closed configurations that is likely to facilitate entry and exit of the TPP ligand. Correspondingly, we posit that such features of the TPP aptamer domain contribute directly to the mechanism of riboswitch-mediated translational regulation.

Transcriptional pausing coordinates folding of the aptamer domain and the expression platform of a riboswitch

Riboswitches are cis-acting elements that regulate gene expression by affecting transcriptional termination or translational initiation in response to binding of a metabolite. A typical riboswitch is made of an upstream aptamer domain and a downstream expression platform. Both domains participate in the folding and structural rearrangement in the absence or presence of its cognate metabolite. RNA polymerase pausing is a fundamental property of transcription that can influence RNA folding. Here we show that pausing plays an important role in the folding and conformational rearrangement of the Escherichia coli btuB riboswitch during transcription by the E. coli RNA polymerase. This riboswitch consists of an approximately 200 nucleotide, coenzyme B12 binding aptamer domain and an approximately 40 nucleotide expression platform that controls the ribosome access for translational initiation. We found that transcriptional pauses at strategic locations facilitate folding and structural rearrangement of the full-length riboswitch, but have minimal effect on the folding of the isolated aptamer domain. Pausing at these regulatory sites blocks the formation of alternate structures and plays a chaperoning role that couples folding of the aptamer domain and the expression platform. Pausing at strategic locations may be a general mechanism for coordinated folding and conformational rearrangements of riboswitch structures that underlie their response to environmental cues.

Folding of a transcriptionally acting preQ1 riboswitch

7-Aminomethyl-7-deazaguanine (preQ(1)) sensitive mRNA domains belong to the smallest riboswitches known to date. Although recent efforts have revealed the three-dimensional architecture of the ligand-aptamer complex less is known about the molecular details of the ligand-induced response mechanism that modulates gene expression. We present an in vitro investigation on the ligand-induced folding process of the preQ(1) responsive RNA element from Fusobacterium nucleatum using biophysical methods, including fluorescence and NMR spectroscopy of site-specifically labeled riboswitch variants. We provide evidence that the full-length riboswitch domain adopts two different coexisting stem-loop structures in the expression platform. Upon addition of preQ(1), the equilibrium of the competing hairpins is significantly shifted. This system therefore, represents a finely tunable antiterminator/terminator interplay that impacts the in vivo cellular response mechanism. A model is presented how a riboswitch that provides no obvious overlap between aptamer and terminator stem-loop solves this communication problem by involving bistable sequence determinants.