Dynamic competition between alternative structures in viroid RNAs simulated by an RNA folding algorithm

The Secondary Structure of Potato Spindle Tuber Viroid Determines Its Infectivity in Nicotiana benthamiana

The function of RNAs is determined by their structure. However, studying the relationship between RNA structure and function often requires altering RNA sequences to modify the structures, which leads to the neglect of the importance of RNA sequences themselves. In our research, we utilized potato spindle tuber viroid (PSTVd), a circular-form non-coding infectious RNA, as a model with which to investigate the role of a specific rod-like structure in RNA function. By generating linear RNA transcripts with different start sites, we established 12 PSTVd forms with different secondary structures while maintaining the same sequence. The RNA secondary structures were predicted using the mfold tool and validated through native PAGE gel electrophoresis after in vitro RNA folding. Analysis using plant infection assays revealed that the formation of a correct rod-like structure is crucial for the successful infection of PSTVd. Interestingly, the inability of PSTVd forms with non-rod-like structures to infect plants could be partially compensated by increasing the amount of linear viroid RNA transcripts, suggesting the existence of additional RNA secondary structures, such as the correct rod-like structure, alongside the dominant structure in the RNA inoculum of these forms. Our study demonstrates the critical role of RNA secondary structures in determining the function of infectious RNAs.

Design of nucleic acid strands with long low-barrier folding pathways

A major goal of natural computing is to design biomolecules, such as nucleic acid sequences, that can be used to perform computations. We design sequences of nucleic acids that are "guaranteed" to have long folding pathways relative to their length. This particular sequences with high probability follow low-barrier folding pathways that visit a large number of distinct structures. Long folding pathways are interesting, because they demonstrate that natural computing can potentially support long and complex computations. Formally, we provide the first scalable designs of molecules whose low-barrier folding pathways, with respect to a simple, stacked pair energy model, grow superlinearly with the molecule length, but for which all significantly shorter alternative folding pathways have an energy barrier that is [Formula: see text] times that of the low-barrier pathway for any [Formula: see text] and a sufficiently long sequence.

Cotranscriptional folding kinetics of ribonucleic acid secondary structures

We develop a systematic helix-based computational method to predict RNA folding kinetics during transcription. In our method, the transcription is modeled as stepwise process, where each step is the transcription of a nucleotide. For each step, the kinetics algorithm predicts the population kinetics, transition pathways, folding intermediates, and the transcriptional folding products. The folding pathways, rate constants, and the conformational populations for cotranscription folding show contrastingly different features than the refolding kinetics for a fully transcribed chain. The competition between the transcription speed and rate constants for the transitions between the different nascent structures determines the RNA folding pathway and the end product of folding. For example, fast transcription favors the formation of branch-like structures than rod-like structures and chain elongation in the folding process may reduce the probability of the formation of misfolded structures. Furthermore, good theory-experiment agreements suggest that our method may provide a reliable tool for quantitative prediction for cotranscriptional RNA folding, including the kinetics for the population distribution for the whole conformational ensemble.

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.

Finding a common motif of RNA sequences using genetic programming: the GeRNAMo system

We focus on finding a consensus motif of a set of homologous or functionally related RNA molecules. Recent approaches to this problem have been limited to simple motifs, require sequence alignment, and make prior assumptions concerning the data set. We use genetic programming to predict RNA consensus motifs based solely on the data set. Our system -- dubbed GeRNAMo (Genetic programming of RNA Motifs) -- predicts the most common motifs without sequence alignment and is capable of dealing with any motif size. Our program only requires the maximum number of stems in the motif, and if prior knowledge is available the user can specify other attributes of the motif (e.g., the range of the motif's minimum and maximum sizes), thereby increasing both sensitivity and speed. We describe several experiments using either ferritin iron response element (IRE); signal recognition particle (SRP); or microRNA sequences, showing that the most common motif is found repeatedly, and that our system offers substantial advantages over previous methods.

Encoding folding paths of RNA switches

RNA co-transcriptional folding has long been suspected to play an active role in helping proper native folding of ribozymes and structured regulatory motifs in mRNA untranslated regions (UTRs). Yet, the underlying mechanisms and coding requirements for efficient co-transcriptional folding remain unclear. Traditional approaches have intrinsic limitations to dissect RNA folding paths, as they rely on sequence mutations or circular permutations that typically perturb both RNA folding paths and equilibrium structures. Here, we show that exploiting sequence symmetries instead of mutations can circumvent this problem by essentially decoupling folding paths from equilibrium structures of designed RNA sequences. Using bistable RNA switches with symmetrical helices conserved under sequence reversal, we demonstrate experimentally that native and transiently formed helices can guide efficient co-transcriptional folding into either long-lived structure of these RNA switches. Their folding path is controlled by the order of helix nucleations and subsequent exchanges during transcription, and may also be redirected by transient antisense interactions. Hence, transient intra- and inter-molecular base pair interactions can effectively regulate the folding of nascent RNA molecules into different native structures, provided limited coding requirements, as discussed from an information theory perspective. This constitutive coupling between RNA synthesis and RNA folding regulation may have enabled the early emergence of autonomous RNA-based regulation networks.

Computer simulation of chaperone effects of Archaeal C/D box sRNA binding on rRNA folding

Archaeal C/D box small RNAs (sRNAs) are homologues of eukaryotic C/D box small nucleolar RNAs (snoRNAs). Their main function is guiding 2'-O-ribose methylation of nucleotides in rRNAs. The methylation requires the pairing of an sRNA antisense element to an rRNA target site with formation of an RNA-RNA duplex. The temporary formation of such a duplex during rRNA maturation is expected to influence rRNA folding in a chaperone-like way, in particular in thermophilic Archaea, where multiple sRNAs with two binding sites are found. Here we investigate possible mechanisms of chaperone function of Archaeoglobus fulgidus and Pyrococcus abyssi C/D box sRNAs using computer simulations of rRNA secondary structure formation by genetic algorithm. The effects of sRNA binding on rRNA structure are introduced as temporary structural constraints during co-transcriptional folding. Comparisons of the final predictions with simulations without sRNA binding and with phylogenetic structures show that sRNAs with two antisense elements may significantly facilitate the correct formation of long-range interactions in rRNAs, in particular at elevated temperatures. The simulations suggest that the main mechanism of this effect is a transient restriction of folding in rRNA domains where the termini are brought together by binding to double-guide sRNAs.

Comparison of P-RnaPredict and mfold--algorithms for RNA secondary structure prediction

MOTIVATION

Ribonucleic acid is vital in numerous stages of protein synthesis; it also possesses important functional and structural roles within the cell. The function of an RNA molecule within a particular organic system is principally determined by its structure. The current physical methods available for structure determination are time-consuming and expensive. Hence, computational methods for structure prediction are sought after. The energies involved by the formation of secondary structure elements are significantly greater than those of tertiary elements. Therefore, RNA structure prediction focuses on secondary structure.

RESULTS

We present P-RnaPredict, a parallel evolutionary algorithm for RNA secondary structure prediction. The speedup provided by parallelization is investigated with five sequences, and a dramatic improvement in speedup is demonstrated, especially with longer sequences. An evaluation of the performance of P-RnaPredict in terms of prediction accuracy is made through comparison with 10 individual known structures from 3 RNA classes (5S rRNA, Group I intron 16S rRNA and 16S rRNA) and the mfold dynamic programming algorithm. P-RnaPredict is able to predict structures with higher true positive base pair counts and lower false positives than mfold on certain sequences.

AVAILABILITY

P-RnaPredict is available for non-commercial usage. Interested parties should contact Kay C. Wiese (wiese@cs.sfu.ca).

Kinefold web server for RNA/DNA folding path and structure prediction including pseudoknots and knots

The Kinefold web server provides a web interface for stochastic folding simulations of nucleic acids on second to minute molecular time scales. Renaturation or co-transcriptional folding paths are simulated at the level of helix formation and dissociation in agreement with the seminal experimental results. Pseudoknots and topologically ‘entangled’ helices (i.e. knots) are efficiently predicted taking into account simple geometrical and topological constraints. To encourage interactivity, simulations launched as immediate jobs are automatically stopped after a few seconds and return adapted recommendations. Users can then choose to continue incomplete simulations using the batch queuing system or go back and modify suggested options in their initial query. Detailed output provide (i) a series of low free energy structures, (ii) an online animated folding path and (iii) a programmable trajectory plot focusing on a few helices of interest to each user. The service can be accessed at .

Cis and trans requirements for rolling circle replication of a satellite RNA

Satellite RNAs usurp the replication machinery of their helper viruses, even though they bear little or no sequence similarity to the helper virus RNA. In Cereal yellow dwarf polerovirus serotype RPV (CYDV-RPV), the 322-nucleotide satellite RNA (satRPV RNA) accumulates to high levels in the presence of the CYDV-RPV helper virus. Rolling circle replication generates multimeric satRPV RNAs that self-cleave via a double-hammerhead ribozyme structure. Alternative folding inhibits formation of a hammerhead in monomeric satRPV RNA. Here we determine helper virus requirements and the effects of mutations and deletions in satRPV RNA on its replication in oat cells. Using in vivo selection of a satRPV RNA pool randomized at specific bases, we found that disruption of the base pairing necessary to form the non-self-cleaving conformation reduced satRPV RNA accumulation. Unlike other satellite RNAs, both the plus and minus strands proved to be equally infectious. Accordingly, very similar essential replication structures were identified in each strand. A different region is required only for encapsidation. The CYDV-RPV RNA-dependent RNA polymerase (open reading frames 1 and 2), when expressed from the nonhelper Barley yellow dwarf luteovirus, was capable of replicating satRPV RNA. Thus, the helper virus's polymerase is the sole determinant of the ability of a virus to replicate a rolling circle satellite RNA. We present a framework for functional domains in satRPV RNA with three types of function: (i) conformational control elements comprising an RNA switch, (ii) self-functional elements (hammerhead ribozymes), and (iii) cis-acting elements that interact with viral proteins.

A mini-RNA containing the tetraloop, wobble-pair and loop E motifs of the central conserved region of potato spindle tuber viroid is processed into a minicircle

A Mini-RNA from potato spindle tuber viroid (PSTVd) was constructed specifically for cleavage and ligation to circles in vitro. It contains the C-domain with the so-called central conserved region (CCR) of PSTVd with a 17 nt duplication in the upper strand and hairpin structures at the left and rights ends of the secondary structure. The CCR was previously shown to be essential for processing of in vitro transcripts. When folded under conditions which favor formation of a kinetically controlled conformation and incubated in a potato nuclear extract, the Mini-RNA is cleaved correctly at the 5'- and the 3'-end and ligated to a circle. Thus, the CCR obviously contains all structural and functional requirements for correct processing and therefore may be regarded as 'processing domain' of PSTVd. Using the Mini-RNA as a model substrate, the structural and functional relevance of its conserved non-canonical motifs GAAA tetraloop, loop E and G:U wobble base pair were studied by mutational analysis. It was found that (i) the conserved GAAA tetraloop is essential for processing by favoring the kinetically controlled conformation, (ii) a G:U wobble base pair at the 5'-cleavage site contributes to its correct recognition and (iii) an unpaired nucleotide in loop E, which is different from the corresponding nucleotide in the conserved loop E motif, is essential for ligation of the 5'- with the 3'-end. Hence all three structural motifs are functional elements for processing in a potato nuclear extract.

Discovery of RNA structural elements using evolutionary computation

RNA molecules fold into characteristic secondary and tertiary structures that account for their diverse functional activities. Many of these RNA structures, or certain structural motifs within them, are thought to recur in multiple genes within a single organism or across the same gene in several organisms and provide a common regulatory mechanism. Search algorithms, such as RNAMotif, can be used to mine nucleotide sequence databases for these repeating motifs. RNAMotif allows users to capture essential features of known structures in detailed descriptors and can be used to identify, with high specificity, other similar motifs within the nucleotide database. However, when the descriptor constraints are relaxed to provide more flexibility, or when there is very little a priori information about hypothesized RNA structures, the number of motif 'hits' may become very large. Exhaustive methods to search for similar RNA structures over these large search spaces are likely to be computationally intractable. Here we describe a powerful new algorithm based on evolutionary computation to solve this problem. A series of experiments using ferritin IRE and SRP RNA stem-loop motifs were used to verify the method. We demonstrate that even when searching extremely large search spaces, of the order of 10(23) potential solutions, we could find the correct solution in a fraction of the time it would have taken for exhaustive comparisons.

Replication of Potato spindle tuber viroid in cultured cells of tobacco and Nicotiana benthamiana: the role of specific nucleotides in determining replication levels for host adaptation

We have developed an electroporation protocol to inoculate cultured cells of tobacco and Nicotiana benthamiana with in vitro transcripts of Potato spindle tuber viroid (PSTVd) to characterize viroid structural features that determine replication efficiency at the cellular level. Both (+)- and (-)-strands of PSTVd were detected by Northern blots as early as 6 h postinoculation (h.p.i.). Accumulation of the (+)-circular PSTVd increased very rapidly starting at 24 h.p.i. and continued beyond 6 days postinoculation. Viroid accumulation in individual cells was visualized by in situ hybridization, which showed that 60-70% of the cells were infected. Previous work showed that C259 --> U substitution converted tomato isolate PSTVd(KF440-2) into a strain that is infectious on tobacco (M. Wassenegger, R. L. Spieker, S. Thalmeir, F.-U. Gast, L. Riedel, and H. L. Sänger, 1996. Virology 226, 191-197). Similarly, C259 --> U or U257 --> A substitution in the Intermediate strain (PSTVd(Int)) conferred infectivity in tobacco (Y. Zhu, Y. Qi, Y. Xun, R. Owens, and B. Ding, 2002. Plant Physiol. 130, 138-146). Our replication assays in tobacco-cultured cells demonstrated that U257 --> A and C259 --> U substitutions each enhanced PSTVd replication by 5- to 10-fold. Replacement of U257 with C, but not with G, also led to enhanced replication in tobacco cells. Replacement of C259 with nucleotide A or G did not enhance replication. Elevated accumulation of the (-)- and (+)-strands of these mutants was in part due to enhanced transcription. Interestingly, all of the nucleotide changes did not alter PSTVd replication levels in N. benthamiana cells. These results provide insights about PSTVd structures that modulate replication efficiency in adapting to a specific host.

Detection and analysis of hairpin II, an essential metastable structural element in viroid replication intermediates

In (-)-stranded replication intermediates of the potato spindle tuber viroid (PSTVd) a thermodynamically metastable structure containing a specific hairpin structure (HP II) has been proposed to be essential for viroid replication. In the present work a method was devised allowing the direct detection of the HP II structure in vitro and in vivo using a biophysical approach. An RNA oligonucleotide was constructed which specifically binds to the HP II loop region in transient (-)-strand intermediates. Analysis of the resulting oligonucleotide/HP II complexes on temperature-gradient gels enabled us to follow the formation of HP II during in vitro transcription by T7 RNA polymerase. Moreover, we were able to demonstrate the formation of HP II during viroid replication in potato (Solanum tuberosum) cells.

A pH-jump approach for investigating secondary structure refolding kinetics in RNA

It has been shown that premature translation of the plasmid-mediated toxin in hok/sok of plasmid R1 and pnd/pndB of plasmid R483 is prevented during transcription of the hok and pnd mRNAs by the formation of metastable hairpins at the 5'-end of the mRNA. Here, an experimental approach is presented, which allows the accurate measurement of the refolding kinetics of the 5'-end RNA fragments in vitro without chemically modifying the RNA. The method is based on acid denaturation followed by a pH-jump to neutral pH as a novel way to trap kinetically favoured RNA secondary structures, allowing the measurement of a wide range of biologically relevant refolding rates, with or without the use of standard stopped-flow equipment. The refolding rates from the metastable to the stable conformation in both the hok74 and pnd58 5'-end RNA fragments were determined by using UV absorbance changes corresponding to the structural rearrangements. The measured energy barriers showed that the refolding path does not need complete unfolding of the metastable structures before the formation of the final structures. Two alternative models of such a pathway are discussed.

Temporal translational control by a metastable RNA structure

Programmed cell death by the hok/sok locus of plasmid R1 relies on a complex translational control mechanism. The highly stable hok mRNA is activated by 3'-end exonucleolytical processing. Removal of the mRNA 3' end releases a 5'-end sequence that triggers refolding of the mRNA. The refolded hok mRNA is translatable but can also bind the inhibitory Sok antisense RNA. Binding of Sok RNA leads to irreversible mRNA inactivation by an RNase III-dependent mechanism. A coherent model predicts that during transcription hok mRNA must be refractory to translation and antisense RNA binding. Here we provide genetic evidence for the existence of a 5' metastable structure in hok mRNA that locks the nascent transcript in an inactive configuration in vivo. Consistently, the metastable structure reduces the rate of Sok RNA binding and completely blocks hok translation in vitro. Structural analyses of native RNAs strongly support that the 5' metastable structure exists in the nascent transcript. Further structural analyses reveal that the mRNA 3' end triggers refolding of the mRNA 5' end into the more stable tac-stem conformation. These results provide a profound understanding of an unusual and intricate post-transcriptional control mechanism.

RNA folding pathway functional intermediates: their prediction and analysis

The massively parallel genetic algorithm (GA) for RNA structure prediction uses the concepts of mutation, recombination, and survival of the fittest to evolve a population of thousands of possible RNA structures toward a solution structure. As described below, the properties of the algorithm are ideally suited to use in the prediction of possible folding pathways and functional intermediates of RNA molecules given their sequences. Utilizing Stem Trace, an interactive visualization tool for RNA structure comparison, analysis of not only the solution ensembles developed by the algorithm, but also the stages of development of each of these solutions, can give strong insight into these folding pathways. The GA allows the incorporation of information from biological experiments, making it possible to test the influence of particular interactions between structural elements on the dynamics of the folding pathway. These methods are used to reveal the folding pathways of the potato spindle tuber viroid (PSTVd) and the host killing mechanism of Escherichia coli plasmid R1, both of which are successfully explored through the combination of the GA and Stem Trace. We also present novel intermediate folds of each molecule, which appear to be phylogenetically supported, as determined by use of the methods described below.

Characterization of the initiation sites of both polarity strands of a viroid RNA reveals a motif conserved in sequence and structure

Viroids replicate through a rolling-circle mechanism in which the infecting circular RNA and its complementary (-) strand are transcribed. The precise site at which transcription starts was investigated for the avocado sunblotch viroid (ASBVd), the type species of the family of viroids with hammerhead ribozymes. Linear ASBVd (+) and (-) RNAs begin with a UAAAA sequence that maps to similar A+U-rich terminal loops in their predicted quasi-rod-like secondary structures. The sequences around the initiation sites of ASBVd, which replicates and accumulates in the chloroplast, are similar to the promoters of a nuclear-encoded chloroplastic RNA polymerase (NEP), supporting the involvement of an NEP-like activity in ASBVd replication. Since RNA folding appears to be kinetically determined, the specific location of both ASBVd initiation sites provides a mechanistic insight into how the nascent ASBVd strands may fold in vivo. The approach used here, in vitro capping and RNase protection assays, may be useful for investigating the initiation sites of other small circular RNA replicons.

Satellite cereal yellow dwarf virus-RPV (satRPV) RNA requires a douXble hammerhead for self-cleavage and an alternative structure for replication

The 110 nt hammerhead ribozyme in the satellite RNA of cereal yellow dwarf virus-RPV (satRPV RNA) folds into an alternative conformation that inhibits self-cleavage. This alternative structure comprises a pseudoknot with base-pairing between loop (L1) and a single-stranded bulge (L2a), which are located in hammerhead stems I and II, respectively. Mutations that disrupt this base-pairing, or otherwise cause the ribozyme to more closely resemble a canonical hammerhead, greatly increase self-cleavage. In a more natural multimeric sequence context containing the full-length satRPV RNA and two copies of the hammerhead, wild-type RNA cleaves much more efficiently than in the 110 nt context. Mutations in the upstream hammerhead, including a knock-out in the catalytic core, affect cleavage at the downstream cleavage site, indicating that multimers of satRPV RNA cleave via a double hammerhead. The double hammerhead includes base-pairing between two copies of the L1 sequence which extends stem I. Disruption of L1-L1 base-pairing slows cleavage of the multimer. L1-L2a base-pairing is required for efficient replication of satRPV RNA in oat protoplasts. Mutations that affect self-cleavage of the multimer do not correlate with replication efficiency, indicating that the ability to self-cleave is not a primary determinant of replication. We present a replication model in which multimeric satRPV RNA folds into alternative conformations that cannot form in the monomer. One potential metastable intermediate conformation involves L1-L2a base-pairing that may facilitate formation of the double hammerhead. However, we conclude that L1-L2a also performs some other essential function in the satRPV RNA replication cycle, because the L1-L2a base-pairing is more important than efficient self-cleavage for replication.