Kinetics and thermodynamics make different contributions to RNA folding in vitro and in yeast

Multiplexed assessment of engineered bacterial constructs for intracellular β-galactosidase expression by redox amplification on catechol-chitosan modified nanoporous gold

Synthetic biology approaches for rewiring of bacterial constructs to express particular intracellular factors upon induction with the target analyte are emerging as sensing paradigms for applications in environmental and in vivo monitoring. To aid in the design and optimization of bacterial constructs for sensing analytes, there is a need for lysis-free intracellular detection modalities that monitor the signal level and kinetics of expressed factors within different modified bacteria in a multiplexed manner, without requiring cumbersome surface immobilization. Herein, an electrochemical detection system on nanoporous gold that is electrofabricated with a biomaterial redox capacitor is presented for quantifying β-galactosidase expressed inside modified Escherichia coli constructs upon induction with dopamine. This nanostructure-mediated redox amplification approach on a microfluidic platform allows for multiplexed assessment of the expressed intracellular factors from different bacterial constructs suspended in distinct microchannels, with no need for cell lysis or immobilization. Since redox mediators present over the entire depth of the microchannel can interact with the electrode and with the E. coli construct in each channel, the platform exhibits high sensitivity and enables multiplexing. We envision its application in assessing synthetic biology-based approaches for comparing specificity, sensitivity, and signal response time upon induction with target analytes of interest.

Probing RNA Conformational Equilibria within the Functional Cellular Context

Low-abundance short-lived non-native conformations referred to as excited states (ESs) are increasingly observed in vitro and implicated in the folding and biological activities of regulatory RNAs. We developed an approach for assessing the relative abundance of RNA ESs within the functional cellular context. Nuclear magnetic resonance (NMR) spectroscopy was used to estimate the degree to which substitution mutations bias conformational equilibria toward the inactive ES in vitro. The cellular activity of the ES-stabilizing mutants was used as an indirect measure of the conformational equilibria within the functional cellular context. Compensatory mutations that restore the ground-state conformation were used to control for changes in sequence. Using this approach, we show that the ESs of two regulatory RNAs from HIV-1, the transactivation response element (TAR) and the Rev response element (RRE), likely form in cells with abundances comparable to those measured in vitro, and their targeted stabilization may provide an avenue for developing anti-HIV therapeutics.

Tying the knot in the tetrahydrofolate (THF) riboswitch: A molecular basis for gene regulation

Effective gene regulation by the tetrahydrofolate riboswitch depends not only on ligand affinity but also on the kinetics of ligand association, which involves two cooperative binding sites. We have determined a 1.9-Å resolution crystal structure of the ligand-free THF riboswitch aptamer. The pseudoknot binding site 'unwinds' in the absence of ligand, whereby the adjacent helical domains (P1, P2, and P3) become disjointed, resulting in rotation and misalignment of the gene-regulatory P1 helix with respect to P3. In contrast, the second binding site at the three-way junction, which is the first to fold, is structurally conserved between apo and holo forms. This suggests a kinetic role for this site, in which binding of the first ligand molecule to the stably folded three-way junction promotes formation of the regulatory pseudoknot site and subsequent binding of the second molecule. As such, these findings provide a molecular basis for both conformational switching and kinetic control.

RNA structures in alternative splicing and back-splicing

Alternative splicing greatly expands the transcriptomic and proteomic diversities related to physiological and developmental processes in higher eukaryotes. Splicing of long noncoding RNAs, and back- and trans- splicing further expanded the regulatory repertoire of alternative splicing. RNA structures were shown to play an important role in regulating alternative splicing and back-splicing. Application of novel sequencing technologies made it possible to identify genome-wide RNA structures and interaction networks, which might provide new insights into RNA splicing regulation in vitro to in vivo. The emerging transcription-folding-splicing paradigm is changing our understanding of RNA alternative splicing regulation. Here, we review the insights into the roles and mechanisms of RNA structures in alternative splicing and back-splicing, as well as how disruption of these structures affects alternative splicing and then leads to human diseases. This article is categorized under: RNA Processing > Splicing Regulation/Alternative Splicing RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems.

Engineering a Circular Riboregulator in Escherichia coli

RNAs of different shapes and sizes, natural or synthetic, can regulate gene expression in prokaryotes and eukaryotes. Circular RNAs have recently appeared to be more widespread than previously thought, but their role in prokaryotes remains elusive. Here, by inserting a riboregulatory sequence within a group I permuted intron-exon ribozyme, we created a small noncoding RNA that self-splices to produce a circular riboregulator in Escherichia coli. We showed that the resulting riboregulator can trans-activate gene expression by interacting with a cis-repressed messenger RNA. We characterized the system with a fluorescent reporter and with an antibiotic resistance marker, and we modeled this novel posttranscriptional mechanism. This first reported example of a circular RNA regulating gene expression in E. coli adds to an increasing repertoire of RNA synthetic biology parts, and it highlights that topological molecules can play a role in the case of prokaryotic regulation.

Urea-aromatic interactions in biology

Noncovalent interactions are key determinants in both chemical and biological processes. Among such processes, the hydrophobic interactions play an eminent role in folding of proteins, nucleic acids, formation of membranes, protein-ligand recognition, etc.. Though this interaction is mediated through the aqueous solvent, the stability of the above biomolecules can be highly sensitive to any small external perturbations, such as temperature, pressure, pH, or even cosolvent additives, like, urea-a highly soluble small organic molecule utilized by various living organisms to regulate osmotic pressure. A plethora of detailed studies exist covering both experimental and theoretical regimes, to understand how urea modulates the stability of biological macromolecules. While experimentalists have been primarily focusing on the thermodynamic and kinetic aspects, theoretical modeling predominantly involves mechanistic information at the molecular level, calculating atomistic details applying the force field approach to the high level electronic details using the quantum mechanical methods. The review focuses mainly on examples with biological relevance, such as (1) urea-assisted protein unfolding, (2) urea-assisted RNA unfolding, (3) urea lesion interaction within damaged DNA, (4) urea conduction through membrane proteins, and (5) protein-ligand interactions those explicitly address the vitality of hydrophobic interactions involving exclusively the urea-aromatic moiety.

Massively parallel RNA device engineering in mammalian cells with RNA-Seq

Synthetic RNA-based genetic devices dynamically control a wide range of gene-regulatory processes across diverse cell types. However, the limited throughput of quantitative assays in mammalian cells has hindered fast iteration and interrogation of sequence space needed to identify new RNA devices. Here we report developing a quantitative, rapid and high-throughput mammalian cell-based RNA-Seq assay to efficiently engineer RNA devices. We identify new ribozyme-based RNA devices that respond to theophylline, hypoxanthine, cyclic-di-GMP, and folinic acid from libraries of ~22,700 sequences in total. The small molecule responsive devices exhibit low basal expression and high activation ratios, significantly expanding our toolset of highly functional ribozyme switches. The large datasets obtained further provide conserved sequence and structure motifs that may be used for rationally guided design. The RNA-Seq approach offers a generally applicable strategy for developing broad classes of RNA devices, thereby advancing the engineering of genetic devices for mammalian systems.

A highly ordered, nonprotective MALAT1 ENE structure is adopted prior to triplex formation

The 3' end of the ∼7 kb lncRNA MALAT1 contains an evolutionarily and structurally conserved element for nuclear expression (ENE) which confers protection from cellular degradation pathways. Formation of an ENE triple helix is required to support transcript accumulation, leading to persistent oncogenic activity of MALAT1 in multiple cancer types. Though the specific mechanism of triplex-mediated protection remains unknown, the MALAT1 ENE triplex has been identified as a promising target for therapeutic intervention. Interestingly, a maturation step of the nascent lncRNA 3' end is required prior to triplex formation. We hypothesize that disruption of the maturation or folding process may be a viable mechanism of inhibition. To assess putative cotranscriptional ENE conformations prior to triplex formation, we perform microsecond MD simulations of a partially folded ENE conformation and the ENE triplex. We identify a highly ordered ENE structure prior to triplex formation. Extensive formation of U•U base pairs within the large U-rich internal loops produces a global rod-like architecture. We present a three-dimensional structure of the isolated ENE motif, the global features of which are consistent with small angle X-ray scattering (SAXS) experiments. Our structural model represents a nonprotective conformation of the MALAT1 ENE, providing a molecular description useful for future mechanistic and inhibition studies. We anticipate that targeting stretches of U•U pairs within the ENE motif will prove advantageous for the design of therapeutics targeting this oncogenic lncRNA.

Conformational Dynamics of the RNA G-Quadruplex and its Effect on Translation Efficiency

During translation, intracellular mRNA folds co-transcriptionally and must refold following the passage of ribosome. The mRNAs can be entrapped in metastable structures during these folding events. In the present study, we evaluated the conformational dynamics of the kinetically favored, metastable, and hairpin-like structure, which disturbs the thermodynamically favored G-quadruplex structure, and its effect on co-transcriptional translation in prokaryotic cells. We found that nascent mRNA forms a metastable hairpin-like structure during co-transcriptional folding instead of the G-quadruplex structure. When the translation progressed co-transcriptionally before the metastable hairpin-like structure transition to the G-quadruplex, function of the G-quadruplex as a roadblock of the ribosome was sequestered. This suggested that kinetically formed RNA structures had a dominant effect on gene expression in prokaryotes. The results of this study indicate that it is critical to consider the conformational dynamics of RNA-folding to understand the contributions of the mRNA structures in controlling gene expression.

Challenges and approaches to predicting RNA with multiple functional structures

The revolution in sequencing technology demands new tools to interpret the genetic code. As in vivo transcriptome-wide chemical probing techniques advance, new challenges emerge in the RNA folding problem. The emphasis on one sequence folding into a single minimum free energy structure is fading as a new focus develops on generating RNA structural ensembles and identifying functional structural features in ensembles. This review describes an efficient combinatorially complete method and three free energy minimization approaches to predicting RNA structures with more than one functional fold, as well as two methods for analysis of a thermodynamics-based Boltzmann ensemble of structures. The review then highlights two examples of viral RNA 3'-UTR regions that fold into more than one conformation and have been characterized by single molecule fluorescence energy resonance transfer or NMR spectroscopy. These examples highlight the different approaches and challenges in predicting structure and function from sequence for RNA with multiple biological roles and folds. More well-defined examples and new metrics for measuring differences in RNA structures will guide future improvements in prediction of RNA structure and function from sequence.

mRNAs and lncRNAs intrinsically form secondary structures with short end-to-end distances

The 5' and 3' termini of RNA play important roles in many cellular processes. Using Förster resonance energy transfer (FRET), we show that mRNAs and lncRNAs have an intrinsic propensity to fold in the absence of proteins into structures in which the 5' end and 3' end are ≤7 nm apart irrespective of mRNA length. Computational estimates suggest that the inherent proximity of the ends is a universal property of most mRNA and lncRNA sequences. Only guanosine-depleted RNA sequences with low sequence complexity are unstructured and exhibit end-to-end distances expected for the random coil conformation of RNA. While the biological implications remain to be explored, short end-to-end distances could facilitate the binding of protein factors that regulate translation initiation by bridging mRNA 5' and 3' ends. Furthermore, our studies provide the basis for measuring, computing and manipulating end-to-end distances and secondary structure in RNA in research and biotechnology.

RNA Regulations and Functions Decoded by Transcriptome-wide RNA Structure Probing

RNA folds into intricate structures that are crucial for its functions and regulations. To date, a multitude of approaches for probing structures of the whole transcriptome, i.e., RNA structuromes, have been developed. Applications of these approaches to different cell lines and tissues have generated a rich resource for the study of RNA structure–function relationships at a systems biology level. In this review, we first introduce the designs of these methods and their applications to study different RNA structuromes. We emphasize their technological differences especially their unique advantages and caveats. We then summarize the structural insights in RNA functions and regulations obtained from the studies of RNA structuromes. And finally, we propose potential directions for future improvements and studies.

In silico methods for co-transcriptional RNA secondary structure prediction and for investigating alternative RNA structure expression

RNA transcripts are the primary products of active genes in any living organism, including many viruses. Their cellular destiny not only depends on primary sequence signals, but can also be determined by RNA structure. Recent experimental evidence shows that many transcripts can be assigned more than a single functional RNA structure throughout their cellular life and that structure formation happens co-transcriptionally, i.e. as the transcript is synthesised in the cell. Moreover, functional RNA structures are not limited to non-coding transcripts, but can also feature in coding transcripts. The picture that now emerges is that RNA structures constitute an additional layer of information that can be encoded in any RNA transcript (and on top of other layers of information such as protein-context) in order to exert a wide range of functional roles. Moreover, different encoded RNA structures can be expressed at different stages of a transcript's life in order to alter the transcript's behaviour depending on its actual cellular context. Similar to the concept of alternative splicing for protein-coding genes, where a single transcript can yield different proteins depending on cellular context, it is thus appropriate to propose the notion of alternative RNA structure expression for any given transcript. This review introduces several computational strategies that my group developed to detect different aspects of RNA structure expression in vivo. Two aspects are of particular interest to us: (1) RNA secondary structure features that emerge during co-transcriptional folding and (2) functional RNA structure features that are expressed at different times of a transcript's life and potentially mutually exclusive.

Bridging the gap between in vitro and in vivo RNA folding

Deciphering the folding pathways and predicting the structures of complex three-dimensional biomolecules is central to elucidating biological function. RNA is single-stranded, which gives it the freedom to fold into complex secondary and tertiary structures. These structures endow RNA with the ability to perform complex chemistries and functions ranging from enzymatic activity to gene regulation. Given that RNA is involved in many essential cellular processes, it is critical to understand how it folds and functions in vivo. Within the last few years, methods have been developed to probe RNA structures in vivo and genome-wide. These studies reveal that RNA often adopts very different structures in vivo and in vitro, and provide profound insights into RNA biology. Nonetheless, both in vitro and in vivo approaches have limitations: studies in the complex and uncontrolled cellular environment make it difficult to obtain insight into RNA folding pathways and thermodynamics, and studies in vitro often lack direct cellular relevance, leaving a gap in our knowledge of RNA folding in vivo. This gap is being bridged by biophysical and mechanistic studies of RNA structure and function under conditions that mimic the cellular environment. To date, most artificial cytoplasms have used various polymers as molecular crowding agents and a series of small molecules as cosolutes. Studies under such in vivo-like conditions are yielding fresh insights, such as cooperative folding of functional RNAs and increased activity of ribozymes. These observations are accounted for in part by molecular crowding effects and interactions with other molecules. In this review, we report milestones in RNA folding in vitro and in vivo and discuss ongoing experimental and computational efforts to bridge the gap between these two conditions in order to understand how RNA folds in the cell.

A DEAD-box RNA helicase promotes thermodynamic equilibration of kinetically trapped RNA structures in vivo

RNAs must assemble into specific structures in order to carry out their biological functions, but in vitro RNA folding reactions produce multiple misfolded structures that fail to exchange with functional structures on biological time scales. We used carefully designed self-cleaving mRNAs that assemble through well-defined folding pathways to identify factors that differentiate intracellular and in vitro folding reactions. Our previous work showed that simple base-paired RNA helices form and dissociate with the same rate and equilibrium constants in vivo and in vitro. However, exchange between adjacent secondary structures occurs much faster in vivo, enabling RNAs to quickly adopt structures with the lowest free energy. We have now used this approach to probe the effects of an extensively characterized DEAD-box RNA helicase, Mss116p, on a series of well-defined RNA folding steps in yeast. Mss116p overexpression had no detectable effect on helix formation or dissociation kinetics or on the stability of interdomain tertiary interactions, consistent with previous evidence that intracellular factors do not affect these folding parameters. However, Mss116p overexpression did accelerate exchange between adjacent helices. The nonprocessive nature of RNA duplex unwinding by DEAD-box RNA helicases is consistent with a branch migration mechanism in which Mss116p lowers barriers to exchange between otherwise stable helices by the melting and annealing of one or two base pairs at interhelical junctions. These results suggest that the helicase activity of DEAD-box proteins like Mss116p distinguish intracellular RNA folding pathways from nonproductive RNA folding reactions in vitro and allow RNA structures to overcome kinetic barriers to thermodynamic equilibration in vivo.

Secondary Structure Prediction of Single Sequences Using RNAstructure

RNA secondary structure is often predicted using folding thermodynamics. RNAstructure is a software package that includes structure prediction by free energy minimization, prediction of base pairing probabilities, prediction of structures composed of highly probably base pairs, and prediction of structures with pseudoknots. A user-friendly graphical user interface is provided, and this interface works on Windows, Apple OS X, and Linux. This chapter provides protocols for using RNAstructure for structure prediction.

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.

Challenge of mimicking the influences of the cellular environment on RNA structure by PEG-induced macromolecular crowding

There are large differences between the cellular environment and the conditions widely used to study RNA in vitro. SHAPE RNA structure probing in Escherichia coli cells has shown that the cellular environment stabilizes both long-range and local tertiary interactions in the adenine riboswitch aptamer domain. Synthetic crowding agents are widely used to understand the forces that stabilize RNA structure and in efforts to recapitulate the cellular environment under simplified experimental conditions. Here, we studied the structure and ligand binding ability of the adenine riboswitch in the presence of the macromolecular crowding agent, polyethylene glycol (PEG). Ethylene glycol and low-molecular mass PEGs destabilized RNA structure and caused the riboswitch to sample secondary structures different from those observed in simple buffered solutions or in cells. In the presence of larger PEGs, longer-range loop-loop interactions were more similar to those in cells than in buffer alone, consistent with prior work showing that larger PEGs stabilize compact RNA states. Ligand affinity was weakened by low-molecular mass PEGs but increased with high-molecular mass PEGs, indicating that PEG cosolvents exert complex chemical and steric effects on RNA structure. Regardless of polymer size, however, nucleotide-resolution structural characteristics observed in cells were not recapitulated in PEG solutions. Our results reveal that the cellular environment is difficult to recapitulate in vitro; mimicking the cellular state will likely require a combination of crowding agents and other chemical species.

RNA structure: merging chemistry and genomics for a holistic perspective

The advent of deep sequencing technology has unexpectedly advanced our structural understanding of molecules composed of nucleic acids. A significant amount of progress has been made recently extrapolating the chemical methods to probe RNA structure into sequencing methods. Herein we review some of the canonical methods to analyze RNA structure, and then we outline how these have been used to probe the structure of many RNAs in parallel. The key is the transformation of structural biology problems into sequencing problems, whereby sequencing power can be interpreted to understand nucleic acid proximity, nucleic acid conformation, or nucleic acid-protein interactions. Utilizing such technologies in this way has the promise to provide novel structural insights into the mechanisms that control normal cellular physiology and provide insight into how structure could be perturbed in disease.

Structure-guided mutational analysis of gene regulation by the Bacillus subtilis pbuE adenine-responsive riboswitch in a cellular context

Riboswitches are a broadly distributed form of RNA-based gene regulation in Bacteria and, more rarely, Archaea and Eukarya. Most often found in the 5'-leader sequence of bacterial mRNAs, they are generally composed of two functional domains: a receptor (aptamer) domain that binds an effector molecule and a regulatory domain (or expression platform) that instructs the expression machinery. One of the most studied riboswitches is the Bacillus subtilis adenine-responsive pbuE riboswitch, which regulates gene expression at the transcriptional level, up-regulating expression in response to increased intracellular effector concentrations. In this work, we analyzed sequence and structural elements that contribute to efficient ligand-dependent regulatory activity in a co-transcriptional and cellular context. Unexpectedly, we found that the P1 helix, which acts as the antitermination element of the switch in this RNA, supported ligand-dependent activation of a reporter gene over a broad spectrum of lengths from 3 to 10 bp. This same trend was also observed using a minimal in vitro single-turnover transcription assay, revealing that this behavior is intrinsic to the RNA sequence. We also found that the sequences at the distal tip of the terminator not directly involved in alternative secondary structure formation are highly important for efficient regulation. These data strongly support a model in which the switch is highly localized to the P1 helix adjacent to the ligand-binding pocket that likely presents a local kinetic block to invasion of the aptamer by the terminator.

Probing viral genomic structure: alternative viewpoints and alternative structures for satellite tobacco mosaic virus RNA

Viral RNA structure prediction is a valuable tool for development of drugs against viral disease. This work discusses different approaches to predicting encapsidated viral RNA and highlights satellite tobacco mosaic virus (STMV) RNA as a model system with excellent crystallography data. Fundamentally important issues for debate include thermodynamic versus kinetic control of virus assembly and the possible consequences of quasi-species in the primary structure on RNA secondary structure prediction of a single structure or an ensemble of structures. Multiple computational tools and chemical reagents are now available for improved viral RNA structure prediction. Two different predicted structures for encapsidated STMV RNA result from differences in three main areas: a different approach and philosophy to studying encapsidated viral RNA, an emphasis on different RNA motifs, and technical differences in computational methods and chemical reagents. The experiments with traditional chemical probing and SHAPE reagents are compared in terms of chemistry, results, and interpretation for STMV RNA as well as other RNA protein assemblies, such as the 5'UTR of HIV and the ribosome. This discussion of the challenges of viral RNA structure prediction will lead to new experiments and improved future predictions for viral RNA.

RNA structural analysis by evolving SHAPE chemistry

RNA is central to the flow of biological information. From transcription to splicing, RNA localization, translation, and decay, RNA is intimately involved in regulating every step of the gene expression program, and is thus essential for health and understanding disease. RNA has the unique ability to base-pair with itself and other nucleic acids to form complex structures. Hence the information content in RNA is not simply its linear sequence of bases, but is also encoded in complex folding of RNA molecules. A general chemical functionality that all RNAs have is a 2'-hydroxyl group in the ribose ring, and the reactivity of the 2'-hydroxyl in RNA is gated by local nucleotide flexibility. In other words, the 2'-hydroxyl is reactive at single-stranded and conformationally flexible positions but is unreactive at nucleotides constrained by base-pairing. Recent efforts have been focused on developing reagents that modify RNA as a function of RNA 2' hydroxyl group reactivity. Such RNA structure probing techniques can be read out by primer extension in experiments termed RNA SHAPE (selective 2'- hydroxyl acylation and primer extension). Herein, we describe the efforts devoted to the design and utilization of SHAPE probes for characterizing RNA structure. We also describe current technological advances that are being applied to utilize SHAPE chemistry with deep sequencing to probe many RNAs in parallel. The merging of chemistry with genomics is sure to open the door to genome-wide exploration of RNA structure and function.

On the importance of cotranscriptional RNA structure formation

The expression of genes, both coding and noncoding, can be significantly influenced by RNA structural features of their corresponding transcripts. There is by now mounting experimental and some theoretical evidence that structure formation in vivo starts during transcription and that this cotranscriptional folding determines the functional RNA structural features that are being formed. Several decades of research in bioinformatics have resulted in a wide range of computational methods for predicting RNA secondary structures. Almost all state-of-the-art methods in terms of prediction accuracy, however, completely ignore the process of structure formation and focus exclusively on the final RNA structure. This review hopes to bridge this gap. We summarize the existing evidence for cotranscriptional folding and then review the different, currently used strategies for RNA secondary-structure prediction. Finally, we propose a range of ideas on how state-of-the-art methods could be potentially improved by explicitly capturing the process of cotranscriptional structure formation.

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.

Transient RNA structure features are evolutionarily conserved and can be computationally predicted

Functional RNA structures tend to be conserved during evolution. This finding is, for example, exploited by comparative methods for RNA secondary structure prediction that currently provide the state-of-art in terms of prediction accuracy. We here provide strong evidence that homologous RNA genes not only fold into similar final RNA structures, but that their folding pathways also share common transient structural features that have been evolutionarily conserved. For this, we compile and investigate a non-redundant data set of 32 sequences with known transient and final RNA secondary structures and devise a dedicated computational analysis pipeline.

COFOLD: an RNA secondary structure prediction method that takes co-transcriptional folding into account

Existing state-of-the-art methods that take a single RNA sequence and predict the corresponding RNA secondary structure are thermodynamic methods. These aim to predict the most stable RNA structure. There exists by now ample experimental and theoretical evidence that the process of structure formation matters and that sequences in vivo fold while they are being transcribed. None of the thermodynamic methods, however, consider the process of structure formation. Here, we present a conceptually new method for predicting RNA secondary structure, called CoFold, that takes effects of co-transcriptional folding explicitly into account. Our method significantly improves the state-of-art in terms of prediction accuracy, especially for long sequences of >1000 nt in length.

RNA SHAPE analysis in living cells

RNA structure has important roles in practically every facet of gene regulation, but the paucity of in vivo structural probes limits current understanding. Here we design, synthesize and demonstrate two new chemical probes that enable selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) in living cells. RNA structures in human, mouse, fly, yeast and bacterial cells are read out at single-nucleotide resolution, revealing tertiary contacts and RNA-protein interactions.

Basis for ligand discrimination between ON and OFF state riboswitch conformations: the case of the SAM-I riboswitch

Riboswitches are RNA elements that bind to effector ligands and control gene expression. Most consist of two domains. S-Adenosyl Methionine (SAM) binds the aptamer domain of the SAM-I riboswitch and induces conformational changes in the expression domain to form an intrinsic terminator (transcription OFF state). Without SAM the riboswitch forms the transcription ON state, allowing read-through transcription. The mechanistic link between the SAM/aptamer recognition event and subsequent secondary structure rearrangement by the riboswitch is unclear. We probed for those structural features of the Bacillus subtilis yitJ SAM-I riboswitch responsible for discrimination between the ON and OFF states by SAM. We designed SAM-I riboswitch RNA segments forming "hybrid" structures of the ON and OFF states. The choice of segment constrains the formation of a partial P1 helix, characteristic of the OFF state, together with a partial antiterminator (AT) helix, characteristic of the ON state. For most choices of P1 vs. AT helix lengths, SAM binds with micromolar affinity according to equilibrium dialysis. Mutational analysis and in-line probing confirm that the mode of SAM binding by hybrid structures is similar to that of the aptamer. Altogether, binding measurements and in-line probing are consistent with the hypothesis that when SAM is present, stacking interactions with the AT helix stabilize a partially formed P1 helix in the hybrids. Molecular modeling indicates that continuous stacking between the P1 and the AT helices is plausible with SAM bound. Our findings raise the possibility that conformational intermediates may play a role in ligand-induced aptamer folding.

An in vivo selection method to optimize trans-splicing ribozymes

Group I intron ribozymes can repair mutated mRNAs by replacing the 3'-terminal portion of the mRNA with their own 3'-exon. This trans-splicing reaction has the potential to treat genetic disorders and to selectively kill cancer cells or virus-infected cells. However, these ribozymes have not yet been used in therapy, partially due to a low in vivo trans-splicing efficiency. Previous strategies to improve the trans-splicing efficiencies focused on designing and testing individual ribozyme constructs. Here we describe a method that selects the most efficient ribozymes from millions of ribozyme variants. This method uses an in vivo rescue assay where the mRNA of an inactivated antibiotic resistance gene is repaired by trans-splicing group I intron ribozymes. Bacterial cells that express efficient trans-splicing ribozymes are able to grow on medium containing the antibiotic chloramphenicol. We randomized a 5'-terminal sequence of the Tetrahymena thermophila group I intron and screened a library with 9 × 10⁶ ribozyme variants for the best trans-splicing activity. The resulting ribozymes showed increased trans-splicing efficiency and help the design of efficient trans-splicing ribozymes for different sequence contexts. This in vivo selection method can now be used to optimize any sequence in trans-splicing ribozymes.

Sequence requirements for RNA binding by HuR and AUF1

The stability of RNAs bearing AU-rich elements in their 3'-UTRs, and thus the level of expression of their protein products, is regulated by interactions with cytoplasmic RNA-binding proteins. Binding by HuR generally leads to mRNA stabilization and increased protein production, whereas binding by AUF1 isoforms generally lead to rapid degradation of the mRNA and reduced protein production. The exact nature of the interplay between these and other RNA-binding proteins remains unclear, although recent studies have shown close interactions between them and even suggested competition between the two for binding to their cognate recognition sequences. Other recent reports have suggested that the sequences recognized by the two proteins are different. We therefore performed a detailed in vitro analysis of the binding site(s) for HuR and AUF1 present in androgen receptor mRNA to define their exact target sequences, and show that the same sequence is contacted by both proteins. Furthermore, we analysed a proposed HuR target within the 3'-UTR of MTA1 mRNA, and show that the contacted bases lie outside of the postulated motif and are a better match to a classical ARE than the postulated motif. The defining features of these HuR binding sites are their U-richness and single strandedness.

Attenuation of loop-receptor interactions with pseudoknot formation

RNA tetraloops can recognize receptors to mediate long-range interactions in stable natural RNAs. In vitro selected GNRA tetraloop/receptor interactions are usually more ‘G/C-rich’ than their ‘A/U-rich’ natural counterparts. They are not as widespread in nature despite comparable biophysical and chemical properties. Moreover, while AA, AC and GU dinucleotide platforms occur in natural GAAA/11 nt receptors, the AA platform is somewhat preferred to the others. The apparent preference for ‘A/U-rich’ GNRA/receptor interactions in nature might stem from an evolutionary adaptation to avoid folding traps at the level of the larger molecular context. To provide evidences in favor of this hypothesis, several riboswitches based on natural and artificial GNRA receptors were investigated in vitro for their ability to prevent inter-molecular GNRA/receptor interactions by trapping the receptor sequence into an alternative intra-molecular pseudoknot. Extent of attenuation determined by native gel-shift assays and co-transcriptional assembly is correlated to the G/C content of the GNRA receptor. Our results shed light on the structural evolution of natural long-range interactions and provide design principles for RNA-based attenuator devices to be used in synthetic biology and RNA nanobiotechnology.

The shape-shifting quasispecies of RNA: one sequence, many functional folds

E Unus pluribum, or "Of One, Many", may be at the root of decoding the RNA sequence-structure-function relationship. RNAs embody the large majority of genes in higher eukaryotes and fold in a sequence-directed fashion into three-dimensional structures that perform functions conserved across all cellular life forms, ranging from regulating to executing gene expression. While it is the most important determinant of the RNA structure, the nucleotide sequence is generally not sufficient to specify a unique set of secondary and tertiary interactions due to the highly frustrated nature of RNA folding. This frustration results in folding heterogeneity, a common phenomenon wherein a chemically homogeneous population of RNA molecules folds into multiple stable structures. Often, these alternative conformations constitute misfolds, lacking the biological activity of the natively folded RNA. Intriguingly, a number of RNAs have recently been described as capable of adopting multiple distinct conformations that all perform, or contribute to, the same function. Characteristically, these conformations interconvert slowly on the experimental timescale, suggesting that they should be regarded as distinct native states. We discuss how rugged folding free energy landscapes give rise to multiple native states in the Tetrahymena Group I intron ribozyme, hairpin ribozyme, sarcin-ricin loop, ribosome, and an in vitro selected aptamer. We further describe the varying degrees to which folding heterogeneity impacts function in these RNAs, and compare and contrast this impact with that of heterogeneities found in protein folding. Embracing that one sequence can give rise to multiple native folds, we hypothesize that this phenomenon imparts adaptive advantages on any functionally evolving RNA quasispecies.

Transient RNA-protein interactions in RNA folding

The RNA folding trajectory features numerous off-pathway folding traps, which represent conformations that are often equally as stable as the native functional ones. Therefore, the conversion between these off-pathway structures and the native correctly folded ones is the critical step in RNA folding. This process, referred to as RNA refolding, is slow, and is represented by a transition state that has a characteristic high free energy. Because this kinetically limiting process occurs in vivo, proteins (called RNA chaperones) have evolved that facilitate the (re)folding of RNA molecules. Here, we present an overview of how proteins interact with RNA molecules in order to achieve properly folded states. In this respect, the discrimination between static and transient interactions is crucial, as different proteins have evolved a multitude of mechanisms for RNA remodeling. For RNA chaperones that act in a sequence-unspecific manner and without the use of external sources of energy, such as ATP, transient RNA–protein interactions represent the basis of the mode of action. By presenting stretches of positively charged amino acids that are positioned in defined spatial configurations, RNA chaperones enable the RNA backbone, via transient electrostatic interactions, to sample a wider conformational space that opens the route for efficient refolding reactions.

Determination of intracellular RNA folding rates using self-cleaving RNAs

We have developed a system that relies on RNA self-cleavage to report quantitatively on assembly of RNA structures in vivo. Self-cleaving RNA sequences are inserted into mRNAs or snoRNAs and expressed in yeast under the control of a regulated promoter. Chimeric RNAs that contain self-cleaving ribozymes turn over faster than chimeric RNAs that contain a mutationally inactivated ribozyme by an amount that reflects the rate at which the ribozyme folds and self-cleaves. A key feature of this system is the choice of assay conditions that selectively monitor intracellular assembly and self-cleavage by suppressing further ribozyme activity during the analysis.

RNA folding in living cells

RNA folding is the most essential process underlying RNA function. While significant progress has been made in understanding the forces driving RNA folding in vitro, exploring the rules governing intracellular RNA structure formation is still in its infancy. The cellular environment hosts a great diversity of factors that potentially influence RNA folding in vivo. For example, the nature of transcription and translation is known to shape the folding landscape of RNA molecules. Trans-acting factors such as proteins, RNAs and metabolites, among others, are also able to modulate the structure and thus the fate of an RNA. Here we summarize the ongoing efforts to uncover how RNA folds in living cells.

DEAD-box protein facilitated RNA folding in vivo

In yeast mitochondria the DEAD-box helicase Mss116p is essential for respiratory growth by acting as group I and group II intron splicing factor. Here we provide the first structure-based insights into how Mss116p assists RNA folding in vivo. Employing an in vivo chemical probing technique, we mapped the structure of the ai5γ group II intron in different genetic backgrounds to characterize its intracellular fold. While the intron adopts the native conformation in the wt yeast strain, we found that the intron is able to form most of its secondary structure, but lacks its tertiary fold in the absence of Mss116p. This suggests that ai5γ is largely unfolded in the mss116-knockout strain and requires the protein at an early step of folding. Notably, in this unfolded state misfolded substructures have not been observed. As most of the protein-induced conformational changes are located within domain D1, Mss116p appears to facilitate the formation of this largest domain, which is the scaffold for docking of other intron domains. These findings suggest that Mss116p assists the ordered assembly of the ai5γ intron in vivo.

Taming free energy landscapes with RNA chaperones

Many non-coding RNAs fold into complex three-dimensional structures, yet the self-assembly of RNA structure is hampered by mispairing, weak tertiary interactions, electrostatic barriers, and the frequent requirement that the 5' and 3' ends of the transcript interact. This rugged free energy landscape for RNA folding means that some RNA molecules in a population rapidly form their native structure, while many others become kinetically trapped in misfolded conformations. Transient binding of RNA chaperone proteins destabilize misfolded intermediates and lower the transition states between conformations, producing a smoother landscape that increases the rate of folding and the probability that a molecule will find the native structure. DEAD-box proteins couple the chemical potential of ATP hydrolysis with repetitive cycles of RNA binding and release, expanding the range of conditions under which they can refold RNA structures.

mRNA secondary structures fold sequentially but exchange rapidly in vivo

RNAs adopt defined structures to perform biological activities, and conformational transitions among alternative structures are critical to virtually all RNA-mediated processes ranging from metabolite-activation of bacterial riboswitches to pre-mRNA splicing and viral replication in eukaryotes. Mechanistic analysis of an RNA folding reaction in a biological context is challenging because many steps usually intervene between assembly of a functional RNA structure and execution of a biological function. We developed a system to probe mechanisms of secondary structure folding and exchange directly in vivo using self-cleavage to monitor competition between mutually exclusive structures that promote or inhibit ribozyme assembly. In previous work, upstream structures were more effective than downstream structures in blocking ribozyme assembly during transcription in vitro, consistent with a sequential folding mechanism. However, upstream and downstream structures blocked ribozyme assembly equally well in vivo, suggesting that intracellular folding outcomes reflect thermodynamic equilibration or that annealing of contiguous sequences is favored kinetically. We have extended these studies to learn when, if ever, thermodynamic stability becomes an impediment to rapid equilibration among alternative RNA structures in vivo. We find that a narrow thermodynamic threshold determines whether kinetics or thermodynamics govern RNA folding outcomes in vivo. mRNA secondary structures fold sequentially in vivo, but exchange between adjacent secondary structures is much faster in vivo than it is in vitro. Previous work showed that simple base-paired RNA helices dissociate at similar rates in vivo and in vitro so exchange between adjacent structures must occur through a different mechanism, one that likely involves facilitation of branch migration by proteins associated with nascent transcripts.

Catalytic activity as a probe of native RNA folding

As RNAs fold to functional structures, they traverse complex energy landscapes that include many partially folded and misfolded intermediates. For structured RNAs that possess catalytic activity, this activity can provide a powerful means of monitoring folding that is complementary to biophysical approaches. RNA catalysis can be used to track accumulation of the native RNA specifically and quantitatively, readily distinguishing the native structure from intermediates that resemble it and may not be differentiated by other approaches. Here, we outline how to design and interpret experiments using catalytic activity to monitor RNA folding, and we summarize adaptations of the method that have been used to probe aspects of folding well beyond determination of the folding rates.

Probing RNA structure within living cells

RNA folding is the most fundamental process underlying RNA function. RNA structure and associated folding paradigms have been intensively studied in vitro. However, in vivo RNA structure formation has only been explored to a limited extent. To determine the influence of the cellular environment, which differs significantly from the in vitro refolding conditions, on RNA architecture, we have applied a chemical probing technique to assess the structure of catalytic RNAs in living cells. This method is based on the fact that chemicals like dimethyl sulfate readily penetrate cells and modify specific atoms of RNA bases (N1-A, N3-C), provided that these positions are solvent accessible. By mapping the modified residues, one gains substantial information on the architecture of the target RNA on the secondary and tertiary structure level. This method also allows exploration of interactions of the target RNA with ligands such as proteins, metabolites, or other RNA molecules and associated conformational changes. In brief, in vivo chemical probing is a powerful tool to investigate RNA structure in its natural environment and can be easily adapted to study RNAs in different cell types.

The jerky and knotty dynamics of RNA

RNA is known to exhibit a jerky dynamics, as intramolecular thermal motion, on <0.1 micros time scales, is punctuated by infrequent structural rearrangements on much longer time scales, i.e. from >10 micros up to a few minutes or even hours. These rare stochastic events correspond to the formation or dissociation of entire stems through cooperative base pairing/unpairing transitions. Such a clear separation of time scales in RNA dynamics has made it possible to implement coarse grained RNA simulations, which predict RNA folding and unfolding pathways including kinetically trapped structures on biologically relevant time scales of seconds to minutes. RNA folding simulations also enable to predict the formation of pseudoknots, that is, helices interior to loops, which mechanically restrain the relative orientations of other non-nested helices. But beyond static structural constraints, pseudoknots can also strongly affect the folding and unfolding dynamics of RNA, as the order by which successive helices are formed and dissociated can lead to topologically blocked transition intermediates. The resulting knotty dynamics can enhance the stability of RNA switches, improve the efficacy of co-transcriptional folding pathways and lead to unusual self-assembly properties of RNA.

Modulation of alternative splicing by long-range RNA structures in Drosophila

Accurate and efficient recognition of splice sites during pre-mRNA splicing is essential for proper transcriptome expression. Splice site usage can be modulated by secondary structures, but it is unclear if this type of modulation is commonly used or occurs to a significant degree with secondary structures forming over long distances. Using phlyogenetic comparisons of intronic sequences among 12 Drosophila genomes, we elucidated a group of 202 highly conserved pairs of sequences, each at least nine nucleotides long, capable of forming stable stem structures. This set was highly enriched in alternatively spliced introns and introns with weak acceptor sites and long introns, and most occurred over long distances (>150 nucleotides). Experimentally, we analyzed the splicing of several of these introns using mini-genes in Drosophila S2 cells. Wild-type splicing patterns were changed by mutations that opened the stem structure, and restored by compensatory mutations that re-established the base-pairing potential, demonstrating that these secondary structures were indeed implicated in the splice site choice. Mechanistically, the RNA structures masked splice sites, brought together distant splice sites and/or looped out introns. Thus, base-pairing interactions within introns, even those occurring over long distances, are more frequent modulators of alternative splicing than is currently assumed.

RNA folding during transcription: protocols and studies

RNA folds during transcription in the cell. Compared to most in vitro studies where the focus is generally on Mg(2+)-initiated refolding of fully synthesized transcripts, cotranscriptional RNA folding studies better replicate how RNA folds in a cellular environment. Unique aspects of cotranscriptional folding include the 5'- to 3'-polarity of RNA, the transcriptional speed, pausing properties of the RNA polymerase, the effect of the transcriptional complex and associated factors, and the effect of RNA-binding proteins. Identifying strategic pause sites can reveal insights on the folding pathway of the nascent transcript. Structural mapping of the paused transcription complexes identifies important folding intermediates along these pathways. Oligohybridization assays and the appearance of the catalytic activity of a ribozyme either in trans or in cis can be used to monitor cotranscriptional folding under a wide range of conditions. In our laboratory, these methodologies have been applied to study the folding of three highly conserved RNAs (RNase P, SRP, and tmRNA), several circularly permuted forms of a bacterial RNase P RNA, a riboswitch (thiM), and an aptamer-activated ribozyme (glmS).

Predicting the sizes of large RNA molecules

We present a theory of the dependence on sequence of the three-dimensional size of large single-stranded (ss) RNA molecules. The work is motivated by the fact that the genomes of many viruses are large ssRNA molecules-often several thousand nucleotides long-and that these RNAs are spontaneously packaged into small rigid protein shells. We argue that there has been evolutionary pressure for the genome to have overall spatial properties-including an appropriate radius of gyration, R(g)-that facilitate this assembly process. For an arbitrary RNA sequence, we introduce the (thermal) average maximum ladder distance (MLD) and use it as a measure of the "extendedness" of the RNA secondary structure. The MLD values of viral ssRNAs that package into capsids of fixed size are shown to be consistently smaller than those for randomly permuted sequences of the same length and base composition, and also smaller than those of natural ssRNAs that are not under evolutionary pressure to have a compact native form. By mapping these secondary structures onto a linear polymer model and by using MLD as a measure of effective contour length, we predict the R(g) values of viral ssRNAs are smaller than those of nonviral sequences. More generally, we predict the average MLD values of large nonviral ssRNAs scale as N(0.67+/-0.01), where N is the number of nucleotides, and that their R(g) values vary as MLD(0.5) in an ideal solvent, and hence as N(0.34). An alternative analysis, which explicitly includes all branches, is introduced and shown to yield consistent results.

Kinetic redistribution of native and misfolded RNAs by a DEAD-box chaperone

DExD/H-box proteins are ubiquitously involved in RNA-mediated processes and use ATP to accelerate conformational changes in RNA. However, their mechanisms of action, and what determines which RNA species are targeted, are not well understood. Here we show that the DExD/H-box protein CYT-19, a general RNA chaperone, mediates ATP-dependent unfolding of both the native conformation and a long-lived misfolded conformation of a group I catalytic RNA with efficiencies that depend on the stabilities of the RNA species but not on specific structural features. CYT-19 then allows the RNA to refold, changing the distribution from equilibrium to kinetic control. Because misfolding is favoured kinetically, conditions that allow unfolding of the native RNA yield large increases in the population of misfolded species. Our results suggest that DExD/H-box proteins act with sufficient breadth and efficiency to allow structured RNAs to populate a wider range of conformations than would be present at equilibrium. Thus, RNAs may face selective pressure to stabilize their active conformations relative to inactive ones to avoid significant redistribution by DExD/H-box proteins. Conversely, RNAs whose functions depend on forming multiple conformations may rely on DExD/H-box proteins to increase the populations of less stable conformations, thereby increasing their overall efficiencies.

Folding of noncoding RNAs during transcription facilitated by pausing-induced nonnative structures

RNA folding in the cell occurs during transcription. Expedient RNA folding must avoid the formation of undesirable structures as the nascent RNA emerges from the RNA polymerase. We show that efficient folding during transcription of three conserved noncoding RNAs from Escherichia coli, RNase P RNA, signal-recognition particle RNA, and tmRNA is facilitated by their cognate polymerase pausing at specific locations. These pause sites are located between the upstream and downstream portions of all of the native long-range helices in these noncoding RNAs. In the paused complexes, the nascent RNAs form labile structures that sequester these upstream portions in a manner to possibly guide folding. Both the pause sites and the secondary structure of the nonnative portions of the paused complexes are phylogenetically conserved among gamma-proteobacteria. We propose that specific pausing-induced structural formation is a general strategy to facilitate the folding of long-range helices. This polymerase-based mechanism may result in portions of noncoding RNA sequences being evolutionarily conserved for efficient folding during transcription.

Communication between RNA folding domains revealed by folding of circularly permuted ribozymes

To study the role of sequence and topology in RNA folding, we determined the kinetic folding pathways of two circularly permuted variants of the Tetrahymena group I ribozyme, using time-resolved hydroxyl radical footprinting. Circular permutation changes the distance between interacting residues in the primary sequence, without changing the native structure of the RNA. In the natural ribozyme, tertiary interactions in the P4-P6 domain form in 1 s, while interactions in the P3-P9 form in 1-3 min at 42 degrees C. Permutation of the 5' end to G111 in the P4 helix allowed the stable P4-P6 domain to fold in 200 ms at 30 degrees C, five times faster than in the wild-type RNA, while the other domains folded five times more slowly (5-8 min). By contrast, circular permutation of the 5' end to G303 in J8/7 decreased the folding rate of the P4-P6 domain. In this permuted RNA, regions joining P2, P3 and P4 were protected in 500 ms, while the P3-P9 domain was 60-80% folded within 30 s. RNase T(1) digestion and FMN photocleavage showed that circular permutation of the RNA sequence alters the initial ensemble of secondary structures, thereby changing the tertiary folding pathways. Our results show that the natural 5'-to-3' order of the structural domains in group I ribozymes optimizes structural communication between tertiary domains and promotes self-assembly of the catalytic center.