Structural Characterization of the Cotranscriptional Folding of the Thiamin Pyrophosphate Sensing thiC Riboswitch in Escherichia coli.

Riboswitches are RNA-regulating elements that mostly rely on structural changes to modulate gene expression at various levels. Recent studies have revealed that riboswitches may control several regulatory mechanisms cotranscriptionally, i.e., during the transcription elongation of the riboswitch or early in the coding region of the regulated gene. Here, we study the structure of the nascent thiamin pyrophosphate (TPP)-sensing thiC riboswitch in Escherichia coli by using biochemical and enzymatic conventional probing approaches. Our chemical (in-line and lead probing) and enzymatic (nucleases S1, A, T1, and RNase H) probing data provide a comprehensive model of how TPP binding modulates the structure of the thiC riboswitch. Furthermore, by using transcriptional roadblocks along the riboswitch sequence, we find that a certain portion of nascent RNA is needed to sense TPP that coincides with the formation of the P5 stem loop. Together, our data suggest that conventional techniques may readily be used to study cotranscriptional folding of nascent RNAs.

A nascent riboswitch helix orchestrates robust transcriptional regulation through signal integration.

Widespread manganese-sensing transcriptional riboswitches effect the dependable gene regulation needed for bacterial manganese homeostasis in changing environments. Riboswitches - like most structured RNAs - are believed to fold co-transcriptionally, subject to both ligand binding and transcription events; yet how these processes are orchestrated for robust regulation is poorly understood. Through a combination of single-molecule and bulk approaches, we discover how a single Mn2+ ion and the transcribing RNA polymerase (RNAP), paused immediately downstream by a DNA template sequence, are coordinated by the bridging switch helix P1.1 in the representative Lactococcus lactis riboswitch. This coordination achieves a heretofore-overlooked semi-docked global conformation of the nascent RNA, P1.1 base pair stabilization, transcription factor NusA ejection, and RNAP pause extension, thereby enforcing transcription readthrough. Our work demonstrates how a central, adaptable RNA helix functions analogous to a molecular fulcrum of a first-class lever system to integrate disparate signals for finely balanced gene expression control.

Assessing the Quality of Cotranscriptional Folding Simulations.

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

Developing Complex RNA Design Applications in the Infrared Framework.

Applications in biotechnology and bio-medical research call for effective strategies to design novel RNAs with very specific properties. Such advanced design tasks require support by computational tools but at the same time put high demands on their flexibility and expressivity to model the application-specific requirements. To address such demands, we present the computational framework Infrared. It supports developing advanced customized design tools, which generate RNA sequences with specific properties, often in a few lines of Python code. This text guides the reader in tutorial format through the development of complex design applications. Thanks to the declarative, compositional approach of Infrared, we can describe this development as a step-by-step extension of an elementary design task. Thus, we start with generating sequences that are compatible with a single RNA structure and go all the way to RNA design targeting complex positive and negative design objectives with respect to single or even multiple target structures. Finally, we present a "real-world" application of computational design to create an RNA device for biotechnology: we use Infrared to generate design candidates of an artificial "AND" riboswitch, which activates gene expression in the simultaneous presence of two different small metabolites. In these applications, we exploit that the system can generate, in an efficient (fixed-parameter tractable) way, multiple diverse designs that satisfy a number of constraints and have high quality w.r.t. to an objective (by sampling from a Boltzmann distribution).

Fluorescent riboswitch-controlled biosensors for the genome scale analysis of metabolic pathways.

Fluorescent detection in cells has been tremendously developed over the years and now benefits from a large array of reporters that can provide sensitive and specific detection in real time. However, the intracellular monitoring of metabolite levels still poses great challenges due to the often complex nature of detected metabolites. Here, we provide a systematic analysis of thiamin pyrophosphate (TPP) metabolism in Escherichia coli by using a TPP-sensing riboswitch that controls the expression of the fluorescent gfp reporter. By comparing different combinations of reporter fusions and TPP-sensing riboswitches, we determine key elements that are associated with strong TPP-dependent sensing. Furthermore, by using the Keio collection as a proxy for growth conditions differing in TPP levels, we perform a high-throughput screen analysis using high-density solid agar plates. Our study reveals several genes whose deletion leads to increased or decreased TPP levels. The approach developed here could be applicable to other riboswitches and reporter genes, thus representing a framework onto which further development could lead to highly sophisticated detection platforms allowing metabolic screens and identification of orphan riboswitches.

Riboswitch and small RNAs modulate btuB translation initiation in Escherichia coli and trigger distinct mRNA regulatory mechanisms.

Small RNAs (sRNAs) and riboswitches represent distinct classes of RNA regulators that control gene expression upon sensing metabolic or environmental variations. While sRNAs and riboswitches regulate gene expression by affecting mRNA and protein levels, existing studies have been limited to the characterization of each regulatory system in isolation, suggesting that sRNAs and riboswitches target distinct mRNA populations. We report that the expression of btuB in Escherichia coli, which is regulated by an adenosylcobalamin (AdoCbl) riboswitch, is also controlled by the small RNAs OmrA and, to a lesser extent, OmrB. Strikingly, we find that the riboswitch and sRNAs reduce mRNA levels through distinct pathways. Our data show that while the riboswitch triggers Rho-dependent transcription termination, sRNAs rely on the degradosome to modulate mRNA levels. Importantly, OmrA pairs with the btuB mRNA through its central region, which is not conserved in OmrB, indicating that these two sRNAs may have specific targets in addition to their common regulon. In contrast to canonical sRNA regulation, we find that OmrA repression of btuB is lost using an mRNA binding-deficient Hfq variant. Together, our study demonstrates that riboswitch and sRNAs modulate btuB expression, providing an example of cis- and trans-acting RNA-based regulatory systems maintaining cellular homeostasis.

Characterization of the dual regulation by a c-di-GMP riboswitch Bc1 with a long expression platform from Bacillus thuringiensis.

A riboswitch generally regulates the expression of its downstream genes through conformational change in its expression platform (EP) upon ligand binding. The cyclic diguanosine monophosphate (c-di-GMP) class I riboswitch Bc1 is widespread and conserved among Bacillus cereus group species. In this study, we revealed that Bc1 has a long EP with two typical ρ-independent terminator sequences 28 bp apart. The upstream terminator T1 is dominant in vitro, while downstream terminator T2 is more efficient in vivo. Through mutation analysis, we elucidated that Bc1 exerts a rare and incoherent "transcription-translation" dual regulation with T2 playing a crucial role. However, we found that Bc1 did not respond to c-di-GMP under in vitro transcription conditions, and the expressions of downstream genes did not change with fluctuation in intracellular c-di-GMP concentration. To explore this puzzle, we conducted SHAPE-MaP and confirmed the interaction of Bc1 with c-di-GMP. This shows that as c-di-GMP concentration increases, T1 unfolds but T2 remains almost intact and functional. The presence of T2 masks the effect of T1 unwinding, resulting in no response of Bc1 to c-di-GMP. The high Shannon entropy values of EP region imply the potential alternative structures of Bc1. We also found that zinc uptake regulator can specifically bind to the dual terminator coding sequence and slightly trigger the response of Bc1 to c-di-GMP. This work will shed light on the dual-regulation riboswitch and enrich our understanding of the RNA world.IMPORTANCEIn nature, riboswitches are involved in a variety of metabolic regulation, most of which preferentially regulate transcription termination or translation initiation of downstream genes in specific ways. Alternatively, the same or different riboswitches can exist in tandem to enhance regulatory effects or respond to multiple ligands. However, many putative conserved riboswitches have not yet been experimentally validated. Here, we found that the c-di-GMP riboswitch Bc1 with a long EP could form a dual terminator and exhibit non-canonical and incoherent "transcription-translation" dual regulation. Besides, zinc uptake regulator specifically bound to the coding sequence of the Bc1 EP and slightly mediated the action of Bc1. The application of SHAPE-MaP to the dual regulation mechanism of Bc1 may establish the foundation for future studies of such complex untranslated regions in other bacterial genomes.

Development and validation of a generic methyltransferase enzymatic assay based on an SAH riboswitch.

Methylation of proteins and nucleic acids plays a fundamental role in epigenetic regulation, and discovery of methyltransferase (MT) inhibitors is an area of intense activity. Because of the diversity of MTs and their products, assay methods that detect S-adenosylhomocysteine (SAH) - the invariant product of S-adenosylmethionine (SAM)-dependent methylation reactions - offer some advantages over methods that detect specific methylation events. However, direct, homogenous detection of SAH requires a reagent capable of discriminating between SAH and SAM, which differ by a single methyl group. Moreover, MTs are slow enzymes and many have submicromolar affinities for SAM; these properties translate to a need for detection of SAH at low nanomolar concentrations in the presence of excess SAM. To meet these needs, we leveraged the exquisite molecular recognition properties of a naturally occurring SAH-sensing RNA aptamer, or riboswitch. By splitting the riboswitch into two fragments, such that SAH binding induces assembly of a trimeric complex, we engineered sensors that transduce binding of SAH into positive fluorescence polarization (FP) and time resolved Förster resonance energy transfer (TR-FRET) signals. The split riboswitch configuration, called the AptaFluor™ SAH Methyltransferase Assay, allows robust detection of SAH (Z' > 0.7) at concentrations below 10 nM, with overnight signal stability in the presence of typical MT assay components. The AptaFluor assay tolerates diverse MT substrates, including histones, nucleosomes, DNA and RNA, and we demonstrated its utility as a robust, enzymatic assay method for several methyltransferases with SAM Km values < 1 µM. The assay was validated for HTS by performing a pilot screen of 1,280 compounds against the SARS-CoV-2 RNA capping enzyme, nsp14. By enabling direct, homogenous detection of SAH at low nanomolar concentrations, the AptaFluor assay provides a universal platform for screening and profiling MTs at physiologically relevant SAM concentrations.

Molecular dynamics in multidimensional space explains how mutations affect the association path of neomycin to a riboswitch.

Riboswitches are messenger RNA (mRNA) fragments binding specific small molecules to regulate gene expression. A synthetic N1 riboswitch, inserted into yeast mRNA controls the translation of a reporter gene in response to neomycin. However, its regulatory activity is sensitive to single-point RNA mutations, even those distant from the neomycin binding site. While the association paths of neomycin to N1 and its variants remain unknown, recent fluorescence kinetic experiments indicate a two-step process driven by conformational selection. This raises the question of which step is affected by mutations. To address this, we performed all-atom two-dimensional replica-exchange molecular dynamics simulations for N1 and U14C, U14C[Formula: see text], U15A, and A17G mutants, ensuring extensive conformational sampling of both RNA and neomycin. The obtained neomycin association and binding paths, along with multidimensional free-energy profiles, revealed a two-step binding mechanism, consisting of conformational selection and induced fit. Neomycin binds to a preformed N1 conformation upon identifying a stable upper stem and U-turn motif in the riboswitch hairpin. However, the positioning of neomycin in the binding site occurs at different RNA-neomycin distances for each mutant, which may explain their different regulatory activities. The subsequent induced fit arises from the interactions of the neomycin's N3 amino group with RNA, causing the G9 backbone to rearrange. In the A17G mutant, the critical C6-A17/G17 stacking forms at a closer RNA-neomycin distance compared to N1. These findings together with estimated binding free energies coincide with experiments and elucidate why the A17G mutation decreases and U15A enhances N1 activity in response to neomycin.

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.

The discovery of novel noncoding RNAs in 50 bacterial genomes.

Structured noncoding RNAs (ncRNAs) contribute to many important cellular processes involving chemical catalysis, molecular recognition and gene regulation. Few ncRNA classes are broadly distributed among organisms from all three domains of life, but the list of rarer classes that exhibit surprisingly diverse functions is growing. We previously developed a computational pipeline that enables the near-comprehensive identification of structured ncRNAs expressed from individual bacterial genomes. The regions between protein coding genes are first sorted based on length and the fraction of guanosine and cytidine nucleotides. Long, GC-rich intergenic regions are then examined for sequence and structural similarity to other bacterial genomes. Herein, we describe the implementation of this pipeline on 50 bacterial genomes from varied phyla. More than 4700 candidate intergenic regions with the desired characteristics were identified, which yielded 44 novel riboswitch candidates and numerous other putative ncRNA motifs. Although experimental validation studies have yet to be conducted, this rate of riboswitch candidate discovery is consistent with predictions that many hundreds of novel riboswitch classes remain to be discovered among the bacterial species whose genomes have already been sequenced. Thus, many thousands of additional novel ncRNA classes likely remain to be discovered in the bacterial domain of life.

Towards parsimonious generative modeling of RNA families.

Generative probabilistic models emerge as a new paradigm in data-driven, evolution-informed design of biomolecular sequences. This paper introduces a novel approach, called Edge Activation Direct Coupling Analysis (eaDCA), tailored to the characteristics of RNA sequences, with a strong emphasis on simplicity, efficiency, and interpretability. eaDCA explicitly constructs sparse coevolutionary models for RNA families, achieving performance levels comparable to more complex methods while utilizing a significantly lower number of parameters. Our approach demonstrates efficiency in generating artificial RNA sequences that closely resemble their natural counterparts in both statistical analyses and SHAPE-MaP experiments, and in predicting the effect of mutations. Notably, eaDCA provides a unique feature: estimating the number of potential functional sequences within a given RNA family. For example, in the case of cyclic di-AMP riboswitches (RF00379), our analysis suggests the existence of approximately 1039 functional nucleotide sequences. While huge compared to the known <4000 natural sequences, this number represents only a tiny fraction of the vast pool of nearly 1082 possible nucleotide sequences of the same length (136 nucleotides). These results underscore the promise of sparse and interpretable generative models, such as eaDCA, in enhancing our understanding of the expansive RNA sequence space.

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.

T7 phage-assisted evolution of riboswitches using error-prone replication and dual selection.

Leveraging riboswitches, non-coding mRNA fragments pivotal to gene regulation, poses a challenge in effectively selecting and enriching these functional genetic sensors, which can toggle between ON and OFF states in response to their cognate inducers. Here, we show our engineered phage T7, enabling the evolution of a theophylline riboswitch. We have replaced T7's DNA polymerase with a transcription factor controlled by a theophylline riboswitch and have created two types of host environments to propagate the engineered phage. Both types host an error-prone T7 DNA polymerase regulated by a T7 promoter along with another critical gene-either cmk or pifA, depending on the host type. The cmk gene is necessary for T7 replication and is used in the first host type for selection in the riboswitch's ON state. Conversely, the second host type incorporates the pifA gene, leading to abortive T7 infections and used for selection in the riboswitch's OFF state. This dual-selection system, termed T7AE, was then applied to a library of 65,536 engineered T7 phages, each carrying randomized riboswitch variants. Through successive passage in both host types with and without theophylline, we observed an enrichment of phages encoding functional riboswitches that conferred a fitness advantage to the phage in both hosts. The T7AE technique thereby opens new pathways for the evolution and advancement of gene switches, including non-coding RNA-based switches, setting the stage for significant strides in synthetic biology.

Direct and indirect control of Rho-dependent transcription termination by the Escherichia coli lysC riboswitch.

Bacterial riboswitches are molecular structures that play a crucial role in controlling gene expression to maintain cellular balance. The Escherichia coli lysC riboswitch has been previously shown to regulate gene expression through translation initiation and mRNA decay. Recent research suggests that lysC gene expression is also influenced by Rho-dependent transcription termination. Through a series of in silico, in vitro, and in vivo experiments, we provide experimental evidence that the lysC riboswitch directly and indirectly modulates Rho transcription termination. Our study demonstrates that Rho-dependent transcription termination plays a significant role in the cotranscriptional regulation of lysC expression. Together with previous studies, our work suggests that lysC expression is governed by a lysine-sensing riboswitch that regulates translation initiation, transcription termination, and mRNA degradation. Notably, both Rho and RNase E target the same region of the RNA molecule, implying that RNase E may degrade Rho-terminated transcripts, providing a means to selectively eliminate these incomplete messenger RNAs. Overall, this study sheds light on the complex regulatory mechanisms used by bacterial riboswitches, emphasizing the role of transcription termination in the control of gene expression and mRNA stability.

A spermidine riboswitch class in bacteria exploits a close variant of an aptamer for the enzyme cofactor S-adenosylmethionine.

Natural polyamines such as spermidine and spermine cations have characteristics that make them highly likely to be sensed by riboswitches, such as their general affinity to polyanionic RNA and their broad contributions to cell physiology. Despite previous claims that polyamine riboswitches exist, evidence of their biological functions has remained unconvincing. Here, we report that rare variants of bacterial S-adenosylmethionine-I (SAM-I) riboswitches reject SAM and have adapted to selectively sense spermidine. These spermidine-sensing riboswitch variants are associated with genes whose protein products are directly involved in the production of spermidine and other polyamines. Biochemical and genetic assays demonstrate that representatives of this riboswitch class robustly function as genetic "off" switches, wherein spermidine binding causes premature transcription termination to suppress the expression of polyamine biosynthetic genes. These findings confirm the existence of natural spermidine-sensing riboswitches in bacteria and expand the list of variant riboswitch classes that have adapted to bind different ligands.

Creation of Architecturally Minimal Transcriptionally Activating Riboswitches Responsive to Theophylline Reveals an Unconventional Design Strategy.

Riboswitches are noncoding RNA switches that are largely utilized in bacteria and play a significant role in synthetic biology. Nonetheless, their natural counterparts possess lengthy sequences and intricate structures, posing challenges for their modular integration into complex gene circuits. Consequently, it is imperative to develop simplified synthetic riboswitches that can be effortlessly incorporated into gene circuits. The conventional approach to generate synthetic riboswitches entails tedious library construction and extensive screening, which frequently yields suboptimal performance. To overcome this obstacle, alternative methods are urgently needed. In this study, we created a novel approach to designing a diverse set of transcription-activating riboswitches that exhibit high performance and broad compatibility. The strategy involved starting with a synthetic theophylline RNA aptamer and designing an expression platform that forms a transcriptional terminator in its inactive state but switches to an antiterminator when it is activated. Several sequences were designed, constructed, and subjected to virtual screening, resulting in the identification of two transcription-activating riboswitches. These riboswitches were then engineered to reduce the basal leakage and increase the activation level through extending the hairpin region using a screened random sequence. These architecturally minimal synthetic riboswitches were highly adapted to different constitutive promoters in a modular manner, generating a differentially responsive output to theophylline. As a proof-of-principle, the synthetic riboswitches were applied to rewire a synthetic quorum-sensing circuit (QSC). The reprogrammed QSC successfully modulated the temporal responsive profile against the activation. This strategy is expected to expand the variety of high-performance riboswitches that are responsive to different ligands, thereby further facilitating the design of complex genetic circuits.

Engineering U1-Based Tetracycline-Inducible Riboswitches to Control Gene Expression in Mammals.

Synthetic riboswitches are promising regulatory devices due to their small size, lack of immunogenicity, and ability to fine-tune gene expression in the absence of exogenous trans-acting factors. Based on a gene inhibitory system developed at our lab, termed U1snRNP interference (U1i), we developed tetracycline (TC)-inducible riboswitches that modulate mRNA polyadenylation through selective U1 snRNP recruitment. First, we engineered different TC-U1i riboswitches, which repress gene expression unless TC is added, leading to inductions of gene expression of 3-to-4-fold. Second, we developed a technique called Systematic Evolution of Riboswitches by Exponential Enrichment (SEREX), to isolate riboswitches with enhanced U1 snRNP binding capacity and activity, achieving inducibilities of up to 8-fold. Interestingly, by multiplexing riboswitches we increased inductions up to 37-fold. Finally, we demonstrated that U1i-based riboswitches are dose-dependent and reversible and can regulate the expression of reporter and endogenous genes in culture cells and mouse models, resulting in attractive systems for gene therapy applications. Our work probes SEREX as a much-needed technology for the in vitro identification of riboswitches capable of regulating gene expression in vivo.

Structural and dynamic mechanisms for coupled folding and tRNA recognition of a translational T-box riboswitch.

T-box riboswitches are unique riboregulators where gene regulation is mediated through interactions between two highly structured RNAs. Despite extensive structural insights, how RNA-RNA interactions drive the folding and structural transitions of T-box to achieve functional conformations remains unclear. Here, by combining SAXS, single-molecule FRET and computational modeling, we elaborate the folding energy landscape of a translational T-box aptamer consisting of stems I, II and IIA/B, which Mg2+-induced global folding and tRNA binding are cooperatively coupled. smFRET measurements reveal that high Mg2+ stabilizes IIA/B and its stacking on II, which drives the pre-docking of I and II into a competent conformation, subsequent tRNA binding promotes docking of I and II to form a high-affinity tRNA binding groove, of which the essentiality of IIA/B and S-turn in II is substantiated with mutational analysis. We highlight a delicate balance among Mg2+, the intra- and intermolecular RNA-RNA interactions in modulating RNA folding and function.

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.