Design of multistable RNA molecules

We show that the problem of designing RNA sequences that can fold into multiple stable secondary structures can be transformed into a combinatorial optimization problem that can be solved by means of simple heuristics. Hence it is feasible to design RNA switches with prescribed structural alternatives. We discuss the theoretical background and present an efficient tool that allows the design of various types of switches. We argue that both the general properties of the sequence structure map of RNA secondary structures and the ease with which our design tool finds bistable RNAs strongly indicates that RNA switches are easily accessible in evolution. Thus conformational switches are yet another function for which RNA can be employed.     

Genome-Wide Double-Stranded RNA Sequencing Reveals the Functional Significance of Base-Paired RNAs in Arabidopsis

The functional structure of all biologically active molecules is dependent on intra- and inter-molecular interactions. This is especially evident for RNA molecules whose functionality, maturation, and regulation require formation of correct secondary structure through encoded base-pairing interactions. Unfortunately, intra- and inter-molecular base-pairing information is lacking for most RNAs. Here, we marry classical nuclease-based structure mapping techniques with high-throughput sequencing technology to interrogate all base-paired RNA in Arabidopsis thaliana and identify ∼200 new small (sm)RNA–producing substrates of RNA–DEPENDENT RNA POLYMERASE6. Our comprehensive analysis of paired RNAs reveals conserved functionality within introns and both 5′ and 3′ untranslated regions (UTRs) of mRNAs, as well as a novel population of functional RNAs, many of which are the precursors of smRNAs. Finally, we identify intra-molecular base-pairing interactions to produce a genome-wide collection of RNA secondary structure models. Although our methodology reveals the pairing status of RNA molecules in the absence of cellular proteins, previous studies have demonstrated that structural information obtained for RNAs in solution accurately reflects their structure in ribonucleoprotein complexes. Furthermore, our identification of RNA–DEPENDENT RNA POLYMERASE6 substrates and conserved functional RNA domains within introns and both 5′ and 3′ untranslated regions (UTRs) of mRNAs using this approach strongly suggests that RNA molecules are correctly folded into their secondary structure in solution. Overall, our findings highlight the importance of base-paired RNAs in eukaryotes and present an approach that should be widely applicable for the analysis of this key structural feature of RNA.     

Strategies for RNA folding and assembly

RNA is structurally very flexible, which provides the basis for its functional diversity. An RNA molecule can often adopt different conformations, which enables the regulation of its function through folding. Proteins help RNAs reach their functionally active conformation by increasing their structural stability or by chaperoning the folding process. Large, dynamic RNA-protein complexes, such as the ribosome or the spliceosome, require numerous proteins that coordinate conformational switches of the RNA components during assembly and during their respective activities.     

A theophylline responsive riboswitch based on helix slipping controls gene expression in vivo

Riboswitches are newly discovered regulatory elements which control a wide set of basic metabolic pathways. They consist solely of RNA, sense their ligand in a preformed binding pocket and perform a conformational switch in response to ligand binding resulting in altered gene expression. We have utilized the enormous potential of RNA for molecular sensing and conformational changes to develop novel molecular switches with predetermined structural transitions in response to the binding of a small molecule. To validate these in vivo, we exploit the distance-dependent inhibitory potential of secondary structure elements placed close to the bacterial ribosome binding site. We created a translational control element by combining the theophylline aptamer with a helical communication module for which a ligand-dependent one-nucleotide slipping mechanism had been proposed. This structural element was inserted at a position just interfering with translation in the non ligand-bound form. Addition of the ligand then shifts the inhibitory element to a distance which permits efficient translation. We present here a novel regulatory mechanism in the first rationally designed, in vivo active RNA switch. Its use of a slippage mechanism to control gene expression makes it different from natural riboswitches which are based on sequestration or antitermination.     

A novel method for the identification of conserved structural patterns in RNA: From small scale to high-throughput applications

Functional RNA regions are often related to recurrent secondary structure patterns (or motifs), which can exert their role in several different ways, particularly in dictating the interaction with RNA-binding proteins, and acting in the regulation of a large number of cellular processes. Among the available motif-finding tools, the majority focuses on sequence patterns, sometimes including secondary structure as additional constraints to improve their performance. Nonetheless, secondary structures motifs may be concurrent to their sequence counterparts or even encode a stronger functional signal. Current methods for searching structural motifs generally require long pipelines and/or high computational efforts or previously aligned sequences. Here, we present BEAM (BEAr Motif finder), a novel method for structural motif discovery from a set of unaligned RNAs, taking advantage of a recently developed encoding for RNA secondary structure named BEAR (Brand nEw Alphabet for RNAs) and of evolutionary substitution rates of secondary structure elements. Tested in a varied set of scenarios, from small- to large-scale, BEAM is successful in retrieving structural motifs even in highly noisy data sets, such as those that can arise in CLIP-Seq or other high-throughput experiments.     

Identification and characterization of yeast mutants that overcome an experimentally introduced block to splicing at the 3' splice site.

An experimentally introduced secondary structure in exon 2 adjacent to the 3' splice site of a yeast ACT-Escherichia coli lacZ fusion gene abolishes splicing in vivo and inhibits beta-galactosidase production. We have devised a genetic screen to isolate both cis and trans-acting mutants that restore beta-galactosidase activity. Two cis-acting mutants potentially destabilize the stem in the region close to the 3' splice site. One trans-acting mutant, designated rss1-1, partially restores beta-galactosidase activity by both increasing the splicing efficiency and stabilizing the precursor and lariat intermediate. The trans-acting suppression activity of rss1-1 is specific for a particular structure because another artificially introduced secondary structure, which also blocks splicing, is not suppressed by this mutant allele. We have cloned the gene encoding the trans-acting mutant protein. The RSS1 gene is located on Saccharomyces cerevisiae chromosome V and is a single copy, essential gene. The predicted RSS1 protein has marked similarity to members of the putative ATP-dependent RNA helicase family. At the nonpermissive temperature, the rss1-1 mutant allele decreases the steady-state levels of several endogenous messenger RNAs and increases the ratio of pre-mRNA to mRNA of specific messages. RSS1 is likely to play an interesting role in RNA processing.     

Hydrostatic and osmotic pressure study of the hairpin ribozyme

The recent discovery of numerous catalytically active RNAs in various living species as well as the in vitro selection of a large series of RNA aptamers able to bind specifically various molecules such as metabolites and co-factors, emphasize the adaptability of RNAs through the plasticity of their secondary structure. Furthermore, all these observations give support to the "RNA world" hypothesis as a step in the primitive development of life on Earth. On this background, we used high pressure to study the mechanism of action of a model hairpin ribozyme which exhibits self-cleavage and ligation. The activation volume (DeltaV( not equal)) of the cleavage reaction (34+/-4 ml/mol) indicates that an important compaction of the RNA molecule occurs during the reaction and must be accompanied by a significant movement of water molecules . Indeed, such a release of 78+/-4 water molecules per RNA molecule could be measured by complementary osmotic shock experiments. These results are consistent with the information provided by the structural studies which indicate that two loops of the RNA molecule should come into contact for the reaction to occur .The high pressure study of a modified form of the ribozyme whose activity is strictly dependent on the presence of adenine as a co-factor should bring some information about the structural significance of this important DeltaV( not equal) of activation.     

Magnesium dependence of the amplified conformational switch in the trans-acting hepatitis delta virus ribozyme

The ability of divalent metal ions to participate in both structure formation and catalytic chemistry of RNA enzymes (ribozymes) has made it difficult to separate their cause and effect in ribozyme function. For example, the recently solved crystal structures of precursor and product forms of the cis-cleaving genomic hepatitis delta virus (HDV) ribozyme show a divalent metal ion bound in the active site that is released upon catalysis due to an RNA conformational change. This conformational switch is associated with a repositioning of the catalytically involved base C75 in the active-site cleft, thus controlling catalysis. These findings confirm previous data from fluorescence resonance energy transfer (FRET) on a trans-acting form of the HDV ribozyme that found a global conformational change to accompany catalysis. Here, we further test the conformational switch model by measuring the Mg(2+) dependence of the global conformational change of the trans-acting HDV ribozyme, using circular dichroism and time-resolved FRET as complementary probes of secondary and tertiary structure formation, respectively. We observe significant differences in both structure and Mg(2+) affinity of the precursor and product forms, in the presence and absence of 300 mM Na(+) background. The precursor shortens while the product extends with increasing Mg(2+) concentration, essentially amplifying the structural differences observed in the crystal structures. In addition, the precursor has an approximately 2-fold and approximately 13-fold lower Mg(2+) affinity than the product in secondary and tertiary structure formation, respectively. We also have compared the C75 wild-type with the catalytically inactive C75U mutant and find significant differences in global structure and Mg(2+) affinity for both their precursor and product forms. Significantly, the Mg(2+) affinity of the C75 wild-type is 1.7-2.1-fold lower than that of the C75U mutant, in accord with the notion that C75 is essential for a catalytic conformational change that leads to a decrease in the local divalent metal ion affinity and release of a catalytic metal. Thus, a consistent picture emerges in which divalent metal ions and RNA functional groups are intimately intertwined in affecting structural dynamics and catalysis in the HDV ribozyme.     

Structural consensus of RNA templates replicated by Q beta replicase.

Q beta replicase replicates a variety of enzyme-specific small RNAs in addition to the phage genomic RNA. The sequence analysis has revealed that all these RNAs are potentially capable of forming a consensus secondary structure element. It represents a stalk which is formed by the 5'-GGG ... and ... CCCA-3' complementary stretches at the termini of the replicating RNA molecules and adjacent 5'- and 3'-hairpins, which may form a stacking with the stalk. The structure found is rather similar to the analogous structure in the tRNA molecule. The genomic RNA of the Q beta phage and other related phages can also form a similar structural element.     

Biologically-supported structural model for a viral satellite RNA

Satellite RNAs (satRNAs) are a class of small parasitic RNA replicon that associate with different viruses, including plus-strand RNA viruses. Because satRNAs do not encode a polymerase or capsid subunit, they rely on a companion virus to provide these proteins for their RNA replication and packaging. SatRNAs recruit these and other required factors via their RNA sequences and structures. Here, through a combination of chemical probing analysis of RNA structure, phylogenetic structural comparisons, and viability assays of satRNA mutants in infected cells, the biological importance of a deduced higher-order structure for a 619 nt long tombusvirus satRNA was assessed. Functionally-relevant secondary and tertiary RNA structures were identified throughout the length of the satRNA. Notably, a 3′-terminal segment was found to adopt two mutually-exclusive RNA secondary structures, both of which were required for efficient satRNA accumulation. Accordingly, these alternative conformations likely function as a type of RNA switch. The RNA switch was also found to engage in a required long-range kissing-loop interaction with an upstream sequence. Collectively, these results establish a high level of conformational complexity within this small parasitic RNA and provide a valuable structural framework for detailed mechanistic studies.     

A structured viroid RNA serves as a substrate for dicer-like cleavage to produce biologically active small RNAs but is resistant to RNA-induced silencing complex-mediated degradation

RNA silencing is a potent means of antiviral defense in plants and animals. A hallmark of this defense response is the production of 21- to 24-nucleotide viral small RNAs via mechanisms that remain to be fully understood. Many viruses encode suppressors of RNA silencing, and some viral RNAs function directly as silencing suppressors as counterdefense. The occurrence of viroid-specific small RNAs in infected plants suggests that viroids can trigger RNA silencing in a host, raising the question of how these noncoding and unencapsidated RNAs survive cellular RNA-silencing systems. We address this question by characterizing the production of small RNAs of Potato spindle tuber viroid (srPSTVds) and investigating how PSTVd responds to RNA silencing. Our molecular and biochemical studies provide evidence that srPSTVds were derived mostly from the secondary structure of viroid RNAs. Replication of PSTVd was resistant to RNA silencing, although the srPSTVds were biologically active in guiding RNA-induced silencing complex (RISC)-mediated cleavage, as shown with a sensor system. Further analyses showed that without possessing or triggering silencing suppressor activities, the PSTVd secondary structure played a critical role in resistance to RISC-mediated cleavage. These findings support the hypothesis that some infectious RNAs may have evolved specific secondary structures as an effective means to evade RNA silencing in addition to encoding silencing suppressor activities. Our results should have important implications in further studies on RNA-based mechanisms of host-pathogen interactions and the biological constraints that shape the evolution of infectious RNA structures.     

Comparative sequence analysis as a tool for studying the secondary structure of mRNAs.

Analysis of phylogenetically conserved secondary structure has been important in the development of models for the secondary structure of structural RNAs. In this paper, we apply this type of analysis to several families of informational RNAs to evaluate its usefulness in developing secondary structure models for mRNAs and mRNA precursors. We observed many conserved helices in all mRNA groups analyzed. Three criteria were used to identify potential helices which were not conserved solely because of coding sequence constraints, and may therefore be important for the structure and function of the RNA. These results suggest that this approach will be useful in deriving secondary structure models for informational RNAs when used in conjunction with other complementary techniques, and in designing experiments to determine the functional significance of conserved base pairing interactions.     

Analysis of trans-acting response decoy RNA-mediated inhibition of human immunodeficiency virus type 1 transactivation.

Overexpression of trans-acting response element (TAR)-containing sequences (TAR decoys) in CEM SS cells renders cells resistant to human immunodeficiency type 1 (HIV-1) replication. Mutagenesis of TAR was used to investigate the molecular mechanism underlying the observed inhibition. A nucleotide change which disrupts the stem structure of TAR or sequence alterations in the loop abolish the ability of the corresponding TAR decoy RNAs to inhibit HIV replication. A compensatory mutation which restores the stem structure also restores TAR decoy RNA function. Synthesis of viral RNA is drastically reduced in cells expressing a functional TAR decoy RNA, but it is unaffected in cells expressing a mutant form of TAR decoy RNA. It is therefore concluded that overexpression of TAR-containing sequences in CEM SS cells interferes with the process of Tat-mediated transactivation of viral gene expression. However, the phenotype of several mutations suggests that TAR decoy RNA does not inhibit HIV-1 gene expression by simply sequestering Tat but rather does so by sequestering a transactivation protein complex, implying that transactivation requires the cooperative binding of both Tat and a loop-binding cellular factor(s) to TAR. Expression of wild-type or mutant forms of TAR had no discernible effects on cell viability, thus reducing concerns about using TAR decoy RNAs as part of an intracellular immunization protocol for the treatment of AIDS.     

An algorithm for searching RNA motifs in genomic sequences

RNA molecules, which are found in all living cells, fold into characteristic structures that account for their diverse functional activities. Many of these RNA structures consist of a collection of fundamental RNA motifs. The various combinations of RNA basic components form different RNA classes and define their unique structural and functional properties. The availability of many genome sequences makes it possible to search computationally for functional RNAs. Biological experiments indicate that functional RNAs have characteristic RNA structural motifs represented by specific combinations of base pairings and conserved nucleotides in the loop regions. The searching for those well-ordered RNA structures and their homologues in genomic sequences is very helpful for the understanding of RNA-based gene regulation. In this paper, we consider the following problem: given an RNA sequence with a known secondary structure, efficiently determine candidate segments in genomic sequences that can potentially form RNA secondary structures similar to the given RNA secondary structure. Our new bottom-up approach searches all potential stem-loops similar to ones of the given RNA secondary structure first, and then based on located stem-loops, detects potential homologous structural RNAs in genomic sequences.     

In-gel probing of individual RNA conformers within a mixed population reveals a dimerization structural switch in the HIV-1 leader

Definitive secondary structural mapping of RNAs in vitro can be complicated by the presence of more than one structural conformer or multimerization of some of the molecules. Until now, probing a single structure of conformationally flexible RNA molecules has typically relied on introducing stabilizing mutations or adjusting buffer conditions or RNA concentration. Here, we present an in-gel SHAPE (selective 2′OH acylation analysed by primer extension) approach, where a mixed structural population of RNA molecules is separated by non-denaturing gel electrophoresis and the conformers are individually probed within the gel matrix. Validation of the technique using a well-characterized RNA stem-loop structure, the HIV-1 trans-activation response element, showed that authentic structure was maintained and that the method was accurate and highly reproducible. To further demonstrate the utility of in-gel SHAPE, we separated and examined monomeric and dimeric species of the HIV-1 packaging signal RNA. Extensive differences in acylation sensitivity were seen between monomer and dimer. The results support a recently proposed structural switch model of RNA genomic dimerization and packaging, and demonstrate the discriminatory power of in-gel SHAPE.