RNA folding in living cells

RNA G-quadruplexes and stress: emerging mechanisms and functions

RNA G-quadruplexes (rG4s) are noncanonical secondary structures formed by guanine-rich sequences that are found in different regions of RNA molecules. These structures have been implicated in diverse biological processes, including translation, splicing, and RNA stability. Recent studies have suggested that rG4s play a role in the cellular response to stress. This review summarizes the current knowledge on rG4s under stress, focusing on their formation, regulation, and potential functions in stress response pathways. We discuss the molecular mechanisms that regulate the formation of rG4 under different stress conditions and the impact of these structures on RNA metabolism, gene expression, and cell survival. Finally, we highlight the potential therapeutic implications of targeting rG4s for the treatment of stress-related diseases through modulating cell survival.

Rock, scissors, paper: How RNA structure informs function

RNA can fold back on itself to adopt a wide range of structures. These range from relatively simple hairpins to intricate 3D folds and can be accompanied by regulatory interactions with both metabolites and macromolecules. The last 50 yr have witnessed elucidation of an astonishing array of RNA structures including transfer RNAs, ribozymes, riboswitches, the ribosome, the spliceosome, and most recently entire RNA structuromes. These advances in RNA structural biology have deepened insight into fundamental biological processes including gene editing, transcription, translation, and structure-based detection and response to temperature and other environmental signals. These discoveries reveal that RNA can be relatively static, like a rock; that it can have catalytic functions of cutting bonds, like scissors; and that it can adopt myriad functional shapes, like paper. We relate these extraordinary discoveries in the biology of RNA structure to the plant way of life. We trace plant-specific discovery of ribozymes and riboswitches, alternative splicing, organellar ribosomes, thermometers, whole-transcriptome structuromes and pan-structuromes, and conclude that plants have a special set of RNA structures that confer unique types of gene regulation. We finish with a consideration of future directions for the RNA structure-function field.

RNAloops: a database of RNA multiloops

MOTIVATION

Knowledge of the 3D structure of RNA supports discovering its functions and is crucial for designing drugs and modern therapeutic solutions. Thus, much attention is devoted to experimental determination and computational prediction targeting the global fold of RNA and its local substructures. The latter include multi-branched loops-functionally significant elements that highly affect the spatial shape of the entire molecule. Unfortunately, their computational modeling constitutes a weak point of structural bioinformatics. A remedy for this is in collecting these motifs and analyzing their features.

RESULTS

RNAloops is a self-updating database that stores multi-branched loops identified in the PDB-deposited RNA structures. A description of each loop includes angular data-planar and Euler angles computed between pairs of adjacent helices to allow studying their mutual arrangement in space. The system enables search and analysis of multiloops, presents their structure details numerically and visually, and computes data statistics.

AVAILABILITY AND IMPLEMENTATION

RNAloops is freely accessible at https://rnaloops.cs.put.poznan.pl.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Predicting RNA Secondary Structure Using In Vitro and In Vivo Data

The new flow of high-throughput RNA secondary structure data coming from different techniques allowed the further development of machine learning approaches. We developed CROSS and CROSSalive, two algorithms trained on experimental data able to predict the RNA secondary structure propensity both in vitro and in vivo. Since the in vivo folding of RNA molecules depends on multiple factors due to the cellular crowded environment, prediction is a complex problem that needs additional calculations for the interaction with proteins and other molecules. In the following chapter, we will describe the differences in predicting RNA secondary structure propensity using experimental data as input for an Artificial Neural Network (ANN) in vitro and in vivo.

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

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

Action and function of helicases on RNA G-quadruplexes

•

Methodological progresses and piling evidence prove the rG4 biology in vivo.

•

rG4s step in virtually every aspect of RNA biology.

•

Helicases unwinding of rG4s is a fine regulatory layer to the downstream processes and general cell homeostasis.

•

The current knowledge is however limited to a few cell lines.

•

The regulation of helicases themselves is delineating as a important question.

•

Non-helicase rG4-processing proteins likely play a role.

The nucleic acid structure called G-quadruplex (G4) is currently discussed to function in nucleic acid-based mechanisms that influence several cellular processes. They can modulate the cellular machinery either positively or negatively, both at the DNA and RNA level. The majority of what we know about G4 biology comes from DNA G4 (dG4) research. RNA G4s (rG4), on the other hand, are gaining interest as researchers become more aware of their role in several aspects of cellular homeostasis. In either case, the correct regulation of G4 structures within cells is essential and demands specialized proteins able to resolve them. Small changes in the formation and unfolding of G4 structures can have severe consequences for the cells that could even stimulate genome instability, apoptosis or proliferation. Helicases are the most relevant negative G4 regulators, which prevent and unfold G4 formation within cells during different pathways. Yet, and despite their importance only a handful of rG4 unwinding helicases have been identified and characterized thus far. This review addresses the current knowledge on rG4s-processing helicases with a focus on methodological approaches. An example of a non-helicase rG4s-unwinding protein is also briefly described.

The Role of RNA Secondary Structure in Regulation of Gene Expression in Bacteria

Due to the high exposition to changing environmental conditions, bacteria have developed many mechanisms enabling immediate adjustments of gene expression. In many cases, the required speed and plasticity of the response are provided by RNA-dependent regulatory mechanisms. This is possible due to the very high dynamics and flexibility of an RNA structure, which provide the necessary sensitivity and specificity for efficient sensing and transduction of environmental signals. In this review, we will discuss the current knowledge about known bacterial regulatory mechanisms which rely on RNA structure. To better understand the structure-driven modulation of gene expression, we describe the basic theory on RNA structure folding and dynamics. Next, we present examples of multiple mechanisms employed by RNA regulators in the control of bacterial transcription and translation.

Regulation of polyphosphate glucokinase gene expression through co-transcriptional processing in Mycobacterium tuberculosis H37Rv

Transcription is a molecular process that involves the synthesis of RNA chain into the 5'-3' direction, and simultaneously nascent RNA chain tends to form geometric structures, known as co-transcriptional folding. This folding determines the functional properties of RNA molecules and possibly has a critical role during the synthesis. This functioning includes the characterized properties of riboswitches and ribozymes, which are significant when the transcription rate is comparable to the cellular environment. This study reports a novel non-coding region important in the genetic expression of polyphosphate glucokinase (ppgk) in Mycobacterium tuberculosis. This non-coding element of ppgk gene undergoes cleavage activity during the transcriptional process in Mycobacterium tuberculosis. We revealed that cleavage occurs within the nascent RNA, and the resultant cleaved 3'RNA fragment carries the Shine- Dalgarno (SD) sequence and expression platform. We concluded co-transcriptional processing at the non-coding region as the required mechanism for ppgk expression that remains constitutive within the bacterial environment. This study defines the molecular mechanism dependent on the transient but highly active structural features of the nascent RNA.

RNA Remodeling by RNA Chaperones Monitored by RNA Structure Probing

RNA structure probing enables the characterization of RNA secondary structures by established procedures such as the enzyme- or chemical-based detection of single- or double-stranded regions. A specific type of application involves the detection of changes of RNA structures and conformations that are induced by proteins with RNA chaperone activity. This chapter outlines a protocol to analyze RNA structures in vitro in the presence of an RNA-binding protein with RNA chaperone activity. For this purpose, we make use of the methylating agents dimethyl sulfate (DMS) and 1-cyclohexyl-3-(2-morpholinoethyl) carbodiimide metho-p-toluenesulfonate (CMCT). DMS and CMCT specifically modify nucleotides that are not involved in base-pairing or tertiary structure hydrogen bonding and that are not protected by a ligand such as a protein. Modified bases are identified by primer extension. As an example, we describe how the RNA chaperone activity of an isoform of the RNA-binding protein AUF1 induces the flaviviral RNA switch required for viral genome cyclization and viral replication.This chapter includes comprehensive protocols for in vitro synthesis of RNA, ³²P-5'-end labeling of DNA primers, primer extension, as well as the preparation and running of analytical gels. The described methodology should be applicable to any other RNA and protein of interest to identify protein-directed RNA remodeling.

Transcription Regulation Through Nascent RNA Folding

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

Base-pair ambiguity and the kinetics of RNA folding

Background

A pairings of nucleotide sequences. Given this forbidding free-energy landscape, mechanisms have evolved that contribute to a directed and efficient folding process, including catalytic proteins and error-detecting chaperones. Among structural RNA molecules we make a distinction between “bound” molecules, which are active as part of ribonucleoprotein (RNP) complexes, and “unbound,” with physiological functions performed without necessarily being bound in RNP complexes. We hypothesized that unbound molecules, lacking the partnering structure of a protein, would be more vulnerable than bound molecules to kinetic traps that compete with native stem structures. We defined an “ambiguity index”—a normalized function of the primary and secondary structure of an individual molecule that measures the number of kinetic traps available to nucleotide sequences that are paired in the native structure, presuming that unbound molecules would have lower indexes. The ambiguity index depends on the purported secondary structure, and was computed under both the comparative (“gold standard”) and an equilibrium-based prediction which approximates the minimum free energy (MFE) structure. Arguing that kinetically accessible metastable structures might be more biologically relevant than thermodynamic equilibrium structures, we also hypothesized that MFE-derived ambiguities would be less effective in separating bound and unbound molecules.

Results

We have introduced an intuitive and easily computed function of primary and secondary structures that measures the availability of complementary sequences that could disrupt the formation of native stems on a given molecule—an ambiguity index. Using comparative secondary structures, the ambiguity index is systematically smaller among unbound than bound molecules, as expected. Furthermore, the effect is lost when the presumably more accurate comparative structure is replaced instead by the MFE structure.

Conclusions

A statistical analysis of the relationship between the primary and secondary structures of non-coding RNA molecules suggests that stem-disrupting kinetic traps are substantially less prevalent in molecules not participating in RNP complexes. In that this distinction is apparent under the comparative but not the MFE secondary structure, the results highlight a possible deficiency in structure predictions when based upon assumptions of thermodynamic equilibrium.

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.

Allele-specific SHAPE-MaP assessment of the effects of somatic variation and protein binding on mRNA structure

The impact of inherited and somatic mutations on messenger RNA (mRNA) structure remains poorly understood. Recent technological advances that leverage next-generation sequencing to obtain experimental structure data, such as SHAPE-MaP, can reveal structural effects of mutations, especially when these data are incorporated into structure modeling. Here, we analyze the ability of SHAPE-MaP to detect the relatively subtle structural changes caused by single-nucleotide mutations. We find that allele-specific sorting greatly improved our detection ability. Thus, we used SHAPE-MaP with a novel combination of clone-free robotic mutagenesis and allele-specific sorting to perform a rapid, comprehensive survey of noncoding somatic and inherited riboSNitches in two cancer-associated mRNAs, TPT1 and LCP1 Using rigorous thermodynamic modeling of the Boltzmann suboptimal ensemble, we identified a subset of mutations that change TPT1 and LCP1 RNA structure, with approximately 14% of all variants identified as riboSNitches. To confirm that these in vitro structures were biologically relevant, we tested how dependent TPT1 and LCP1 mRNA structures were on their environments. We performed SHAPE-MaP on TPT1 and LCP1 mRNAs in the presence or absence of cellular proteins and found that both mRNAs have similar overall folds in all conditions. RiboSNitches identified within these mRNAs in vitro likely exist under biological conditions. Overall, these data reveal a robust mRNA structural landscape where differences in environmental conditions and most sequence variants do not significantly alter RNA structural ensembles. Finally, predicting riboSNitches in mRNAs from sequence alone remains particularly challenging; these data will provide the community with benchmarks for further algorithmic development.

Computational Methods for Modeling Aptamers and Designing Riboswitches

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

Predicting RNA secondary structures from sequence and probing data

RNA secondary structures have proven essential for understanding the regulatory functions performed by RNA such as microRNAs, bacterial small RNAs, or riboswitches. This success is in part due to the availability of efficient computational methods for predicting RNA secondary structures. Recent advances focus on dealing with the inherent uncertainty of prediction by considering the ensemble of possible structures rather than the single most stable one. Moreover, the advent of high-throughput structural probing has spurred the development of computational methods that incorporate such experimental data as auxiliary information.

Regulatory effects of cotranscriptional RNA structure formation and transitions

RNAs, which play significant roles in many fundamental biological processes of life, fold into sophisticated and precise structures. RNA folding is a dynamic and intricate process, which conformation transition of coding and noncoding RNAs form the primary elements of genetic regulation. The cellular environment contains various intrinsic and extrinsic factors that potentially affect RNA folding in vivo, and experimental and theoretical evidence increasingly indicates that the highly flexible features of the RNA structure are affected by these factors, which include the flanking sequence context, physiochemical conditions, cis RNA-RNA interactions, and RNA interactions with other molecules. Furthermore, distinct RNA structures have been identified that govern almost all steps of biological processes in cells, including transcriptional activation and termination, transcriptional mutagenesis, 5'-capping, splicing, 3'-polyadenylation, mRNA export and localization, and translation. Here, we briefly summarize the dynamic and complex features of RNA folding along with a wide variety of intrinsic and extrinsic factors that affect RNA folding. We then provide several examples to elaborate RNA structure-mediated regulation at the transcriptional and posttranscriptional levels. Finally, we illustrate the regulatory roles of RNA structure and discuss advances pertaining to RNA structure in plants. WIREs RNA 2016, 7:562-574. doi: 10.1002/wrna.1350 For further resources related to this article, please visit the WIREs website.

Real-Time Monitoring of G-Quadruplex Formation during Transcription

Cotranscriptional folding of an RNA transcript enables formation of metastable RNA structures. Thermodynamic and kinetic properties of RNA G-quadruplex formation have previously been investigated using purified guanine-rich oligonucleotides. Here, we describe a method for analysis of cotranscriptional dynamics of the G-quadruplex formation based on real-time monitoring of the fluorescence of G-quadruplex ligands. For RNA sequences with the potential to form mutually exclusive hairpin or G-quadruplex structures, the efficiency of G-quadruplex formation during transcription depended on position of the hairpin forming sequence. The real-time monitoring enabled evaluation of environmental effects on RNA dynamics, as we demonstrated facilitation of post-transcriptional G-quadruplex formation under molecular crowding conditions. The strategy demonstrated here provides folding insights into the G-quadruplex during transcription that should be involved in gene regulation.

Understanding the kinetic mechanism of RNA single base pair formation

RNA functions are intrinsically tied to folding kinetics. The most elementary step in RNA folding is the closing and opening of a base pair. Understanding this elementary rate process is the basis for RNA folding kinetics studies. Previous studies mostly focused on the unfolding of base pairs. Here, based on a hybrid approach, we investigate the folding process at level of single base pairing/stacking. The study, which integrates molecular dynamics simulation, kinetic Monte Carlo simulation, and master equation methods, uncovers two alternative dominant pathways: Starting from the unfolded state, the nucleotide backbone first folds to the native conformation, followed by subsequent adjustment of the base conformation. During the base conformational rearrangement, the backbone either retains the native conformation or switches to nonnative conformations in order to lower the kinetic barrier for base rearrangement. The method enables quantification of kinetic partitioning among the different pathways. Moreover, the simulation reveals several intriguing ion binding/dissociation signatures for the conformational changes. Our approach may be useful for developing a base pair opening/closing rate model.

Traditional Chemical Mapping of RNA Structure In Vitro and In Vivo

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

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

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

Genome-wide profiling of mouse RNA secondary structures reveals key features of the mammalian transcriptome

Background

The understanding of RNA structure is a key feature toward the comprehension of RNA functions and mechanisms of action. In particular, non-coding RNAs are thought to exert their functions by specific secondary structures, but an efficient annotation on a large scale of these structures is still missing.

Results

By using a novel high-throughput method, named chemical inference of RNA structures, CIRS-seq, that uses dimethyl sulfate, and N-cyclohexyl- N'-(2-morpholinoethyl)carbodiimide metho-p-toluenesulfonate to modify RNA residues in single-stranded conformation within native deproteinized RNA secondary structures, we investigate the structural features of mouse embryonic stem cell transcripts. Our analysis reveals an unexpected higher structuring of the 5′ and 3′ untranslated regions compared to the coding regions, a reduced structuring at the Kozak sequence and stop codon, and a three-nucleotide periodicity across the coding region of messenger RNAs. We also observe that ncRNAs exhibit a higher degree of structuring with respect to protein coding transcripts. Moreover, we find that the Lin28a binding protein binds selectively to RNA motifs with a strong preference toward a single stranded conformation.

Conclusions

This work defines for the first time the complete RNA structurome of mouse embryonic stem cells, revealing an extremely distinct RNA structural landscape. These results demonstrate that CIRS-seq constitutes an important tool for the identification of native deproteinized RNA structures.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0491-2) contains supplementary material, which is available to authorized users.

Ribonucleic acid (RNA) biosynthesis in human cancer

In many respects, the most remarkable chemical substances within the genome of eukaryotic cells are remarkable proteins which are the critical structural and functional units of living cells. The specifications for everything that goes in the cell are natural digital-to-digital decoding process in an archive sequence by deoxyribonucleic acid (DNA) and an articulate construction by ribonucleic acid (RNA). The products of DNA transcription are long polymers of ribonucleotides rather than deoxyribonucleotides and are termed ribonucleic acids. Certain deoxyribonucleotide sequences, or genes, give rise to transfer RNA (tRNA) and other ribosomal RNA (rRNA) when transcribed. The ribonucleotide sequences fold extensively and rRNA is associated with specific proteins to yield the essential cell components, ribosomes. Transcription of other special sequences yields messenger RNAs (mRNAs) that contain ribonucleotide sequences that will be ultimately translated into new types of amino acid sequences of functional cellular protein molecules. This switch to a different variety of cellular molecular sequences is complex, but each sequence of the three ribonucleotides specifies the insertion of one particular amino acid into the polypeptide chain under production. Whilst mRNA is considered the vehicle by which genetic information is transmitted from the genome and allocated in the appropriate cytoplasmic sites for translation into protein via cap-dependent mechanism, the actual translation depends also on the presence of other so-called household and luxury protein molecules. Recent evidence suggests RNA species are required at initiation, because treatment of cells with antibiotics or drugs that inhibit RNA synthesis cause a decrease in protein synthesis. The rRNA is necessary as a structural constituent of the ribosomes upon which translation takes place, whereas tRNA is necessary as an adaptor in amino acid activation and elongation protein chains to ribosomes. In this article, we review malignant tumor, with stem like properties, and recent technical advances into the phenomenon of micro-particles and micro-vesicles containing cell-free nucleic acids that circulate plasma. New areas of research have been opened into screening tumor telomerase progression, prognosis of aptamers targeting cell surface, monitoring the efficacy of anticancer therapies, oncogenic transformation of host cell, and RNA polymerases role in the cell cycle progression and differentiation.

Detecting riboSNitches with RNA folding algorithms: a genome-wide benchmark

Ribonucleic acid (RNA) secondary structure prediction continues to be a significant challenge, in particular when attempting to model sequences with less rigidly defined structures, such as messenger and non-coding RNAs. Crucial to interpreting RNA structures as they pertain to individual phenotypes is the ability to detect RNAs with large structural disparities caused by a single nucleotide variant (SNV) or riboSNitches. A recently published human genome-wide parallel analysis of RNA structure (PARS) study identified a large number of riboSNitches as well as non-riboSNitches, providing an unprecedented set of RNA sequences against which to benchmark structure prediction algorithms. Here we evaluate 11 different RNA folding algorithms’ riboSNitch prediction performance on these data. We find that recent algorithms designed specifically to predict the effects of SNVs on RNA structure, in particular remuRNA, RNAsnp and SNPfold, perform best on the most rigorously validated subsets of the benchmark data. In addition, our benchmark indicates that general structure prediction algorithms (e.g. RNAfold and RNAstructure) have overall better performance if base pairing probabilities are considered rather than minimum free energy calculations. Although overall aggregate algorithmic performance on the full set of riboSNitches is relatively low, significant improvement is possible if the highest confidence predictions are evaluated independently.

Unzippers, resolvers and sensors: a structural and functional biochemistry tale of RNA helicases

The centrality of RNA within the biological world is an irrefutable fact that currently attracts increasing attention from the scientific community. The panoply of functional RNAs requires the existence of specific biological caretakers, RNA helicases, devoted to maintain the proper folding of those molecules, resolving unstable structures. However, evolution has taken advantage of the specific position and characteristics of RNA helicases to develop new functions for these proteins, which are at the interface of the basic processes for transference of information from DNA to proteins. RNA helicases are involved in many biologically relevant processes, not only as RNA chaperones, but also as signal transducers, scaffolds of molecular complexes, and regulatory elements. Structural biology studies during the last decade, founded in X-ray crystallography, have characterized in detail several RNA-helicases. This comprehensive review summarizes the structural knowledge accumulated in the last two decades within this family of proteins, with special emphasis on the structure-function relationships of the most widely-studied families of RNA helicases: the DEAD-box, RIG-I-like and viral NS3 classes.

Exploiting post-transcriptional regulation to probe RNA structures in vivo via fluorescence

While RNA structures have been extensively characterized in vitro, very few techniques exist to probe RNA structures inside cells. Here, we have exploited mechanisms of post-transcriptional regulation to synthesize fluorescence-based probes that assay RNA structures in vivo. Our probing system involves the co-expression of two constructs: (i) a target RNA and (ii) a reporter containing a probe complementary to a region in the target RNA attached to an RBS-sequestering hairpin and fused to a sequence encoding the green fluorescent protein (GFP). When a region of the target RNA is accessible, the area can interact with its complementary probe, resulting in fluorescence. By using this system, we observed varied patterns of structural accessibility along the length of the Tetrahymena group I intron. We performed in vivo DMS footprinting which, along with previous footprinting studies, helped to explain our probing results. Additionally, this novel approach represents a valuable tool to differentiate between RNA variants and to detect structural changes caused by subtle mutations. Our results capture some differences from traditional footprinting assays that could suggest that probing in vivo via oligonucleotide hybridization facilitates the detection of folding intermediates. Importantly, our data indicate that intracellular oligonucleotide probing can be a powerful complement to existing RNA structural probing methods.

How to find RNA thermometers

Temperature is one of the decisive signals that a mammalian pathogen has entered its warm-blooded host. Among the many ways to register temperature changes, bacteria often use temperature-modulated structures in the untranslated region of mRNAs. In this article, we describe how such RNA thermometers (RNATs) have been discovered one by one upstream of heat shock and virulence genes in the past, and how next-generation sequencing approaches are able to reveal novel temperature-responsive RNA structures on a global scale.

Noncanonical structures and their thermodynamics of DNA and RNA under molecular crowding: beyond the Watson-Crick double helix

How does molecular crowding affect the stability of nucleic acid structures inside cells? Water is the major solvent component in living cells, and the properties of water in the highly crowded media inside cells differ from that in buffered solution. As it is difficult to measure the thermodynamic behavior of nucleic acids in cells directly and quantitatively, we recently developed a cell-mimicking system using cosolutes as crowding reagents. The influences of molecular crowding on the structures and thermodynamics of various nucleic acid sequences have been reported. In this chapter, we discuss how the structures and thermodynamic properties of nucleic acids differ under various conditions such as highly crowded environments, compartment environments, and in the presence of ionic liquids, and the major determinants of the crowding effects on nucleic acids are discussed. The effects of molecular crowding on the activities of ribozymes and riboswitches on noncanonical structures of DNA- and RNA-like quadruplexes that play important roles in transcription and translation are also described.

Chemical probing of RNA in living cells

RNAs need to adopt a specific architecture to exert their task in cells. While significant progress has been made in describing RNA folding landscapes in vitro, understanding intracellular RNA structure formation is still in its infancy. This is in part due to the complex nature of the cellular environment but also to the limited availability of suitable methodologies. To assess the intracellular structure of large RNAs, we recently applied a chemical probing technique and a metal-induced cleavage assay in vivo. These methods are based on the fact that small molecules, like dimethyl sulfate (DMS), or metal ions, such as Pb(2+), penetrate and spread throughout the cell very fast. Hence, these chemicals are able to modify accessible RNA residues or to induce cleavage of the RNA strand in the vicinity of a metal ion in living cells. Mapping of these incidents allows inferring information on the intracellular conformation, metal ion binding sites or ligand-induced structural changes of the respective RNA molecule. Importantly, in vivo chemical probing can be easily adapted to study RNAs in different cell types.

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.

Perspectives on the mechanism of transcriptional regulation by long non-coding RNAs

Long non-coding RNAs (lncRNAs) are increasingly being recognized as epigenetic regulators of gene transcription. The diversity and complexity of lncRNA genes means that they exert their regulatory effects by a variety of mechanisms. Although there is still much to be learned about the mechanism of lncRNA function, general principles are starting to emerge. In particular, the application of high throughput (deep) sequencing methodologies has greatly advanced our understanding of lncRNA gene function. lncRNAs function as adaptors that link specific chromatin loci with chromatin-remodeling complexes and transcription factors. lncRNAs can act in cis or trans to guide epigenetic-modifier complexes to distinct genomic sites, or act as scaffolds which recruit multiple proteins simultaneously, thereby coordinating their activities. In this review we discuss the genomic organization of lncRNAs, the importance of RNA secondary structure to lncRNA functionality, the multitude of ways in which they interact with the genome, and what evolutionary conservation tells us about their function.

The DEAD-box protein Dbp2 functions with the RNA-binding protein Yra1 to promote mRNP assembly

Eukaryotic gene expression involves numerous biochemical steps that are dependent on RNA structure and ribonucleoprotein (RNP) complex formation. The DEAD-box class of RNA helicases plays fundamental roles in formation of RNA and RNP structure in every aspect of RNA metabolism. In an effort to explore the diversity of biological roles for DEAD-box proteins, our laboratory recently demonstrated that the DEAD-box protein Dbp2 associates with actively transcribing genes and is required for normal gene expression in Saccharomyces cerevisiae. We now provide evidence that Dbp2 interacts genetically and physically with the mRNA export factor Yra1. In addition, we find that Dbp2 is required for in vivo assembly of mRNA-binding proteins Yra1, Nab2, and Mex67 onto poly(A)+ RNA. Strikingly, we also show that Dbp2 is an efficient RNA helicase in vitro and that Yra1 decreases the efficiency of ATP-dependent duplex unwinding. We provide a model whereby messenger ribonucleoprotein (mRNP) assembly requires Dbp2 unwinding activity and once the mRNP is properly assembled, inhibition by Yra1 prevents further rearrangements. Both Yra1 and Dbp2 are conserved in multicellular eukaryotes, suggesting that this constitutes a broadly conserved mechanism for stepwise assembly of mature mRNPs in the nucleus.

The long-range P3 helix of the Tetrahymena ribozyme is disrupted during folding between the native and misfolded conformations

RNAs are prone to misfolding, but how misfolded structures are formed and resolved remains incompletely understood. The Tetrahymena group I intron ribozyme folds in vitro to a long-lived misfolded conformation (M) that includes extensive native structure but is proposed to differ in topology from the native state (N). A leading model predicts that exchange of the topologies requires unwinding of the long-range, core helix P3, despite the presence of P3 in both conformations. To test this model, we constructed 16 mutations to strengthen or weaken P3. Catalytic activity and in-line probing showed that nearly all of the mutants form the M state before folding to N. The P3-weakening mutations accelerated refolding from M (3- to 30-fold) and the P3-strengthening mutations slowed refolding (6- to 1400-fold), suggesting that P3 indeed unwinds transiently. Upon depletion of Mg(2+), the mutations had analogous effects on unfolding from N to intermediates that subsequently fold to M. The magnitudes for the P3-weakening mutations were larger than in refolding from M, and small-angle X-ray scattering showed that the ribozyme expands rapidly to intermediates from which P3 is disrupted subsequently. These results are consistent with previous results indicating unfolding of native peripheral structure during refolding from M, which probably permits rearrangement of the core. Together, our results demonstrate that exchange of the native and misfolded conformations requires loss of a core helix in addition to peripheral structure. Further, the results strongly suggest that misfolding arises from a topological error within the ribozyme core, and a specific topology is proposed.

Mss116p: a DEAD-box protein facilitates RNA folding

RNA folding is an essential aspect underlying RNA-mediated cellular processes. Many RNAs, including large, multi-domain ribozymes, are capable of folding to the native, functional state without assistance of a protein cofactor in vitro. In the cell, trans-acting factors, such as proteins, are however known to modulate the structure and thus the fate of an RNA. DEAD-box proteins, including Mss116p, were recently found to assist folding of group I and group II introns in vitro and in vivo. The underlying mechanism(s) have been studied extensively to explore the contribution of ATP hydrolysis and duplex unwinding in helicase-stimulated intron splicing. Here we summarize the ongoing efforts to understand the novel role of DEAD-box proteins in RNA folding.

Structure and mechanism of purine-binding riboswitches

A riboswitch is a non-protein coding sequence capable of directly binding a small molecule effector without the assistance of accessory proteins to regulate expression of the mRNA in which it is embedded. Currently, over 20 different classes of riboswitches have been validated in bacteria with the promise of many more to come, making them an important means of regulating the genome in the bacterial kingdom. Strikingly, half of the known riboswitches recognize effector compounds that contain a purine or related moiety. In the last decade, significant progress has been made to determine how riboswitches specifically recognize these compounds against the background of many other similar cellular metabolites and transduce this signal into a regulatory response. Of the known riboswitches, the purine family containing guanine, adenine and 2'-deoxyguanosine-binding classes are the most extensively studied, serving as a simple and useful paradigm for understanding how these regulatory RNAs function. This review provides a comprehensive summary of the current state of knowledge regarding the structure and mechanism of these riboswitches, as well as insights into how they might be exploited as therapeutic targets and novel biosensors.

BioMe: biologically relevant metals

In this article, we introduce BioMe (biologically relevant metals), a web-based platform for calculation of various statistical properties of metal-binding sites. Users can obtain the following statistical properties: presence of selected ligands in metal coordination sphere, distribution of coordination numbers, percentage of metal ions coordinated by the combination of selected ligands, distribution of monodentate and bidentate metal-carboxyl, bindings for ASP and GLU, percentage of particular binuclear metal centers, distribution of coordination geometry, descriptive statistics for a metal ion–donor distance and percentage of the selected metal ions coordinated by each of the selected ligands. Statistics is presented in numerical and graphical forms. The underlying database contains information about all contacts within the range of 3 Å from a metal ion found in the asymmetric crystal unit. The stored information for each metal ion includes Protein Data Bank code, structure determination method, types of metal-binding chains [protein, ribonucleic acid (RNA), deoxyribonucleic acid (DNA), water and other] and names of the bounded ligands (amino acid residue, RNA nucleotide, DNA nucleotide, water and other) and the coordination number, the coordination geometry and, if applicable, another metal(s). BioMe is on a regular weekly update schedule. It is accessible at http://metals.zesoi.fer.hr.

RNA catalysis as a probe for chaperone activity of DEAD-box helicases

DEAD-box proteins are vitally important to cellular processes and make up the largest class of helicases. Many DEAD-box proteins function as RNA chaperones by accelerating structural transitions of RNA, which can result in the resolution of misfolded conformers or conversion between functional structures. While the biological importance of chaperone proteins is clear, their mechanisms are incompletely understood. Here, we illustrate how the catalytic activity of certain RNAs can be used to measure RNA chaperone activity. By measuring the amount of substrate converted to product, the fraction of catalytically active molecules is measured over time, providing a quantitative measure of the formation or loss of native RNA. The assays are described with references to group I and group II introns and their ribozyme derivatives, and examples are included that illustrate potential complications and indicate how catalytic activity measurements can be combined with physical approaches to gain insights into the mechanisms of DEAD-box proteins as RNA chaperones.

ATP-dependent roles of the DEAD-box protein Mss116p in group II intron splicing in vitro and in vivo

The yeast DEAD-box protein Mss116p functions as a general RNA chaperone in splicing mitochondrial group I and group II introns. For most of its functions, Mss116p is thought to use ATP-dependent RNA unwinding to facilitate RNA structural transitions, but it has been suggested to assist in the folding of one group II intron (aI5γ) primarily by stabilizing a folding intermediate. Here we compare three aI5γ constructs: one with long exons, one with short exons, and a ribozyme construct lacking exons. The long exons result in slower splicing, suggesting that they misfold and/or stabilize nonnative intronic structures. Nevertheless, Mss116p acceleration of all three constructs depends on ATP and is inhibited by mutations that compromise RNA unwinding, suggesting similar mechanisms. Results of splicing assays and a new two-stage assay that separates ribozyme folding and catalysis indicate that maximal folding of all three constructs by Mss116p requires ATP-dependent RNA unwinding. ATP-independent activation is appreciable for only a subpopulation of the minimal ribozyme construct and not for constructs containing exons. As expected for a general RNA chaperone, Mss116p can also disrupt the native ribozyme, which can refold after Mss116p removal. Finally, using yeast strains with mitochondrial DNA containing only the single intron aI5γ, we show that Mss116p mutants promote splicing in vivo to degrees that correlate with their residual ATP-dependent RNA-unwinding activities. Together, our results indicate that, although DEAD-box proteins play multiple roles in RNA folding, the physiological function of Mss116p in aI5γ splicing includes a requirement for ATP-dependent local unfolding, allowing the conversion of nonfunctional RNA structure into functional RNA structure.

Understanding the transcriptome through RNA structure

RNA structure is crucial for gene regulation and function. In the past, transcriptomes have largely been parsed by primary sequences and expression levels, but it is now becoming feasible to annotate and compare transcriptomes based on RNA structure. In addition to computational prediction methods, the recent advent of experimental techniques to probe RNA structure by high-throughput sequencing has enabled genome-wide measurements of RNA structure and has provided the first picture of the structural organization of a eukaryotic transcriptome - the 'RNA structurome'. With additional advances in method refinement and interpretation, structural views of the transcriptome should help to identify and validate regulatory RNA motifs that are involved in diverse cellular processes and thereby increase understanding of RNA function.