Structural Plasticity and Rapid Evolution in a Viral RNA Revealed by In Vivo Genetic Selection

Satellite RNAs usually lack substantial homology with their helper viruses. The 356-nucleotide satC of Turnip crinkle virus (TCV) is unusual in that its 3′-half shares high sequence similarity with the TCV 3′ end. Computer modeling, structure probing, and/or compensatory mutagenesis identified four hairpins and three pseudoknots in this TCV region that participate in replication and/or translation. Two hairpins and two pseudoknots have been confirmed as important for satC replication. One portion of the related 3′ end of satC that remains poorly characterized corresponds to juxtaposed TCV hairpins H4a and H4b and pseudoknot ψ₃, which are required for the TCV-specific requirement of translation (V. A. Stupina et al., RNA 14:2379-2393, 2008). Replacement of satC H4a with randomized sequence and scoring for fitness in plants by in vivo genetic selection (SELEX) resulted in winning sequences that contain an H4a-like stem-loop, which can have additional upstream sequence composing a portion of the stem. SELEX of the combined H4a and H4b region in satC generated three distinct groups of winning sequences. One group models into two stem-loops similar to H4a and H4b of TCV. However, the selected sequences in the other two groups model into single hairpins. Evolution of these single-hairpin SELEX winners in plants resulted in satC that can accumulate to wild-type (wt) levels in protoplasts but remain less fit in planta when competed against wt satC. These data indicate that two highly distinct RNA conformations in the H4a and H4b region can mediate satC fitness in protoplasts.     

A multifunctional turnip crinkle virus replication enhancer revealed by in vivo functional SELEX

The motif1-hairpin (M1H), located on (-)-strands of Turnip Crinkle Virus (TCV)-associated satellite RNA C (satC), is a replication enhancer and recombination hotspot. Results of in vivo genetic selection (SELEX: systematic evolution of ligands by exponential enrichment), where 28 bases of the M1H were randomized and then subjected to selection in plants, revealed that most winners contained one to three short motifs, many of which in their (-)-sense orientation are found in TCV and satC (-)-strand promoter elements. Ability to replicate in protoplasts correlated with fitness to accumulate in plants with one significant exception. Winner UC, containing only a seven-base replacement sequence, was the second most fit winner, yet replicated no better than a 28-base random replacement sequence. Fitness of satC containing different M1H replacement sequences could be due to enhanced satC replication or enhanced ability to affect TCV movement, since satC interferes with TCV virion accumulation, which is correlated with enhanced movement to younger tissue. Cells inoculated with TCV and UC accumulated fewer virions when compared to other winners that replicated better in protoplasts but were less fit in plants. UC, and other first and second round winners, contained structures that were on average 33% more stable in their (+)-strand orientation, and most formed hairpins with a A-rich sequence at the base. These results suggest that M1H replacement sequences contribute to the fitness of satC by either containing (-)-strand elements that enhance satRNA replication and/or a (+)-strand hairpin flanked with single-stranded sequence that enhances TCV movement.     

Fitness of a turnip crinkle virus satellite RNA correlates with a sequence-nonspecific hairpin and flanking sequences that enhance replication and repress the accumulation of virions

satC, a satellite RNA associated with Turnip crinkle virus (TCV), enhances the ability of the virus to colonize plants by interfering with stable virion accumulation (F. Zhang and A. E. Simon, unpublished data). Previous results suggested that the motif1-hairpin (M1H), a replication enhancer on minus strands, forms a plus-strand hairpin flanked by CA-rich sequence that may be involved in enhancing systemic infection (G. Zhang and A. E. Simon, J. Mol. Biol. 326:35-48, 2003). In this study, sequence and structural requirements of the M1H were further assayed by replacing the 28-base M1H with 10 random bases and then subjecting the pool of satellite RNA to functional selection in plants. Unlike previous results with 28-base replacement sequences (G. Zhang and A. E. Simon, J. Mol. Biol. 326:35-48, 2003), only a few of the 10-base SELEX (systematic evolution of ligands by exponential enrichment) assay winners contained short motifs in their minus-sense orientation that were similar to TCV replication elements. However, all second- and third-round winning replacement sequences folded into hairpins flanked by CA-rich sequence predicted to be more stable on plus strands than minus strands. Plus strands of several of the most fit satellite RNAs contained insertions of CA-rich sequence at the base of their hairpins whose presence correlated with enhanced replication and reduced detection of virions. Deletion of the M1H resulted in no detectable virions despite very low satellite accumulation. These results support the hypothesis that a sequence-nonspecific plus-strand hairpin brings together flanking CA-rich sequences in the M1H region that confers fitness to satC by reducing the accumulation of stable virions.     

The 3' proximal translational enhancer of Turnip crinkle virus binds to 60S ribosomal subunits

During cap-dependent translation of eukaryotic mRNAs, initiation factors interact with the 5′ cap to attract ribosomes. When animal viruses translate in a cap-independent fashion, ribosomes assemble upstream of initiation codons at internal ribosome entry sites (IRES). In contrast, many plant viral genomes do not contain 5′ ends with substantial IRES activity but instead have 3′ translational enhancers that function by an unknown mechanism. A 393-nucleotide (nt) region that includes the entire 3′ UTR of the Turnip crinkle virus (TCV) synergistically enhances translation of a reporter gene when associated with the TCV 5′ UTR. The major enhancer activity was mapped to an internal region of ~140 nt that partially overlaps with a 100-nt structural domain previously predicted to adopt a form with some resemblance to a tRNA, according to a recent study by J.C. McCormack and colleagues. The T-shaped structure binds to 80S ribosomes and 60S ribosomal subunits, and binding is more efficient in the absence of surrounding sequences and in the presence of a pseudoknot that mimics the tRNA-acceptor stem. Untranslated TCV satellite RNA satC, which contains the TCV 3′ end and 6-nt differences in the region corresponding to the T-shaped element, does not detectably bind to 80S ribosomes and is not predicted to form a comparable structure. Binding of the TCV T-shaped element by 80S ribosomes was unaffected by salt-washing, reduced in the presence of AcPhe-tRNA, which binds to the P-site, and enhanced binding of Phe-tRNA to the ribosome A site. Mutations that reduced translation in vivo had similar effects on ribosome binding in vitro. This strong correlation suggests that ribosome entry in the 3′ UTR is a key function of the 3′ translational enhancer of TCV and that the T-shaped element contains some tRNA-like properties.     

Biphasic folding kinetics of RNA pseudoknots and telomerase RNA activity

Using a combined master equation and kinetic cluster approach, we investigate RNA pseudoknot folding and unfolding kinetics. The energetic parameters are computed from a recently developed Vfold model for RNA secondary structure and pseudoknot folding thermodynamics. The folding kinetics theory is based on the complete conformational ensemble, including all the native-like and non-native states. The predicted folding and unfolding pathways, activation barriers, Arrhenius plots, and rate-limiting steps lead to several findings. First, for the PK5 pseudoknot, a misfolded 5' hairpin emerges as a stable kinetic trap in the folding process, and the detrapping from this misfolded state is the rate-limiting step for the overall folding process. The calculated rate constant and activation barrier agree well with the experimental data. Second, as an application of the model, we investigate the kinetic folding pathways for human telomerase RNA (hTR) pseudoknot. The predicted folding and unfolding pathways not only support the proposed role of conformational switch between hairpin and pseudoknot in hTR activity, but also reveal molecular mechanism for the conformational switch. Furthermore, for an experimentally studied hTR mutation, whose hairpin intermediate is destabilized, the model predicts a long-lived transient hairpin structure, and the switch between the transient hairpin intermediate and the native pseudoknot may be responsible for the observed hTR activity. Such finding would help resolve the apparent contradiction between the observed hTR activity and the absence of a stable hairpin.     

3'-End stem-loops of the subviral RNAs associated with turnip crinkle virus are involved in symptom modulation and coat protein binding

Many plant RNA viruses are associated with one or more subviral RNAs. Two subviral RNAs, satellite RNA C (satC) and defective interfering RNA G (diG) intensify the symptoms of their helper, turnip crinkle virus (TCV). However, when the coat protein (CP) of TCV was replaced with that of the related Cardamine chlorotic fleck virus (CCFV), both subviral RNAs attenuated symptoms of the hybrid virus TCV-CP(CCFV). In contrast, when the translation initiation codon of the TCV CP was altered to ACG and reduced levels of CP were synthesized, satC attenuated symptoms while diG neither intensified nor attenuated symptoms. The determinants for this differential symptom modulation were previously localized to the 3'-terminal 100 bases of the subviral RNAs, which contain six positional differences (Q. Kong, J.-W. Oh, C. D. Carpenter, and A. E. Simon, Virology 238:478-485, 1997). In the current study, we have determined that certain sequences within the 3'-terminal stem-loop structures of satC and diG, which also serve as promoters for complementary strand synthesis, are critical for symptom modulation. Furthermore, the ability to attenuate symptoms was correlated with weakened binding of TCV CP to the hairpin structure.     

Three-dimensional structure of a flavivirus dumbbell RNA reveals molecular details of an RNA regulator of replication

Mosquito-borne flaviviruses (MBFVs) including dengue, West Nile, yellow fever, and Zika viruses have an RNA genome encoding one open reading frame flanked by 5′ and 3′ untranslated regions (UTRs). The 3′ UTRs of MBFVs contain regions of high sequence conservation in structured RNA elements known as dumbbells (DBs). DBs regulate translation and replication of the viral RNA genome, functions proposed to depend on the formation of an RNA pseudoknot. To understand how DB structure provides this function, we solved the x-ray crystal structure of the Donggang virus DB to 2.1Å resolution and used structural modeling to reveal the details of its three-dimensional fold. The structure confirmed the predicted pseudoknot and molecular modeling revealed how conserved sequences form a four-way junction that appears to stabilize the pseudoknot. Single-molecule FRET suggests that the DB pseudoknot is a stable element that can regulate the switch between translation and replication during the viral lifecycle by modulating long-range RNA conformational changes.     

Structural domains within the 3' untranslated region of Turnip crinkle virus

The genomes of positive-strand RNA viruses undergo conformational shifts that complicate efforts to equate structures with function. We have initiated a detailed analysis of secondary and tertiary elements within the 3′ end of Turnip crinkle virus (TCV) that are required for viral accumulation in vivo. MPGAfold, a massively parallel genetic algorithm, suggested the presence of five hairpins (H4a, H4b, and previously identified hairpins H4, H5, and Pr) and one H-type pseudoknot (Ψ₃) within the 3′-terminal 194 nucleotides (nt). In vivo compensatory mutagenesis analyses confirmed the existence of H4a, H4b, Ψ₃ and a second pseudoknot (Ψ₂) previously identified in a TCV satellite RNA. In-line structure probing of the 194-nt fragment supported the coexistence of H4, H4a, H4b, Ψ₃ and a pseudoknot that connects H5 and the 3′ end (Ψ₁). Stepwise replacements of TCV elements with the comparable elements from Cardamine chlorotic fleck virus indicated that the complete 142-nt 3′ end, and subsets containing Ψ₃, H4a, and H4b or Ψ₃, H4a, H4b, H5, and Ψ₂, form functional domains for virus accumulation in vivo. A new 3-D molecular modeling protocol (RNA2D3D) predicted that H4a, H4b, H5, Ψ₃, and Ψ₂ are capable of simultaneous existence and bears some resemblance to a tRNA. The related Japanese iris necrotic ring virus does not have comparable domains. These results provide a framework for determining how interconnected elements participate in processes that require 3′ untranslated region sequences such as translation and replication.     

Comparison of Turnip Crinkle Virus RNA-Dependent RNA Polymerase Preparations Expressed in Escherichia coli or Derived from Infected Plants

Turnip crinkle virus (TCV) is a small, plus-sense, single-stranded RNA virus of plants. A virus-coded protein, p88, which is required for replication has been expressed and purified from Escherichia coli. In vitro assays revealed that the recombinant p88 has an RNA-dependent RNA polymerase (RdRp) activity and can also bind to RNA. Deletion of the N-terminal region in p88 resulted in a more active RdRp, while further deletions abolished RdRp activity. Comparison of the E. coli-expressed p88, the N-terminal deletion mutant of p88, and a TCV RdRp preparation obtained from infected plants revealed that these preparations show remarkable similarities in RNA template recognition and usage. Both the recombinant and the plant TCV RdRp preparations are capable of de novo initiation on both plus- and minus-strand satC and satD templates, which are small parasitic RNAs associated with TCV infections. In addition, these RdRp preparations can efficiently recognize the related Tomato bushy stunt virus promoter sequences, including the minus- and plus-strand initiation promoters. Heterologous viral and artificial promoters are recognized poorly by the recombinant and the plant TCV RdRps. Further comparison of the single-component recombinant TCV RdRp and the multicomponent plant TCV RdRp will help dissect the functions of various components of the TCV replicase.     

Evaluation of the conformational switch model for alfalfa mosaic virus RNA replication

Key elements of the conformational switch model describing regulation of alfalfa mosaic virus (AMV) replication (R. C. Olsthoorn, S. Mertens, F. T. Brederode, and J. F. Bol, EMBO J. 18:4856-4864, 1999) have been tested using biochemical assays and functional studies in nontransgenic protoplasts. Although comparative sequence analysis suggests that the 3' untranslated regions of AMV and ilarvirus RNAs have the potential to fold into pseudoknots, we were unable to confirm that a proposed pseudoknot forms or has a functional role in regulating coat protein-RNA binding or viral RNA replication. Published work has suggested that the pseudoknot is part of a tRNA-like structure (TLS); however, we argue that the canonical sequence and functional features that define the TLS are absent. We suggest here that the absence of the TLS correlates directly with the distinctive requirement for coat protein to activate replication in these viruses. Experimental data are evidence that elevated magnesium concentrations proposed to stabilize the pseudoknot structure do not block coat protein binding. Additionally, covarying nucleotide changes proposed to reestablish pseudoknot pairings do not rescue replication. Furthermore, as described in the accompanying paper (L. M. Guogas, S. M. Laforest, and L. Gehrke, J. Virol. 79:5752-5761, 2005), coat protein is not, by definition, inhibitory to minus-strand RNA synthesis. Rather, the activation of viral RNA replication by coat protein is shown to be concentration dependent. We describe the 3' organization model as an alternate model of AMV replication that offers an improved fit to the available data.     

The 3′ end of Turnip crinkle virus contains a highly interactive structure including a translational enhancer that is disrupted by binding to the RNA-dependent RNA polymerase

Precise temporal control is needed for RNA viral genomes to translate sufficient replication-required products before clearing ribosomes and initiating replication. A 3′ translational enhancer in Turnip crinkle virus (TCV) overlaps an internal T-shaped structure (TSS) that binds to 60S ribosomal subunits. The higher-order structure in the region was examined through alteration of critical sequences revealing novel interactions between an H-type pseudoknot and upstream residues, and between the TSS and internal and terminal loops of an upstream hairpin. Our results suggest that the TSS forms a stable scaffold that allows for simultaneous interactions with external sequences through base pairings on both sides of its large internal symmetrical loop. Binding of TCV RNA-dependent RNA polymerase (RdRp) to the region potentiates a widespread conformational shift with substantial rearrangement of the TSS region, including the element required for efficient ribosome binding. Degrading the RdRp caused the RNA to resume its original conformation, suggesting that the initial conformation is thermodynamically favored. These results suggest that the 3′ end of TCV folds into a compact, highly interactive structure allowing RdRp access to multiple elements including the 3′ end, which causes structural changes that potentiate the shift between translation and replication.     

A local, interactive network of 3' RNA elements supports translation and replication of Turnip crinkle virus

The majority of the 3' untranslated region (UTR) of Turnip crinkle virus (TCV) was previously identified as forming a highly interactive structure with a ribosome-binding tRNA-shaped structure (TSS) acting as a scaffold and undergoing a widespread conformational shift upon binding to RNA-dependent RNA polymerase (RdRp). Tertiary interactions in the region were explored by identifying two highly detrimental mutations within and adjacent to a hairpin H4 upstream of the TSS that reduce translation in vivo and cause identical structural changes in the loop of the 3' terminal hairpin Pr. Second-site changes that compensate for defects in translation/accumulation and reverse the structural differences in the Pr loop were found in the Pr stem, as well as in a specific stem within the TSS and within the capsid protein (CP) coding region, suggesting that the second-site changes were correcting a conformational defect and not restoring specific base pairing. The RdRp-mediated conformational shift extended upstream through this CP open reading frame (ORF) region after bypassing much of an intervening, largely unstructured region, supporting a connection between 3' elements and coding region elements. These data suggest that the Pr loop, TSS, and H4 are central elements in the regulation of translation and replication in TCV and allow for development of an RNA interactome that maps the higher-order structure of a postulated RNA domain within the 3' region of a plus-strand RNA virus.     

Predicting RNA pseudoknot folding thermodynamics

Based on the experimentally determined atomic coordinates for RNA helices and the self-avoiding walks of the P (phosphate) and C₄ (carbon) atoms in the diamond lattice for the polynucleotide loop conformations, we derive a set of conformational entropy parameters for RNA pseudoknots. Based on the entropy parameters, we develop a folding thermodynamics model that enables us to compute the sequence-specific RNA pseudoknot folding free energy landscape and thermodynamics. The model is validated through extensive experimental tests both for the native structures and for the folding thermodynamics. The model predicts strong sequence-dependent helix-loop competitions in the pseudoknot stability and the resultant conformational switches between different hairpin and pseudoknot structures. For instance, for the pseudoknot domain of human telomerase RNA, a native-like and a misfolded hairpin intermediates are found to coexist on the (equilibrium) folding pathways, and the interplay between the stabilities of these intermediates causes the conformational switch that may underlie a human telomerase disease.     

In vivo structural analysis of spliced leader RNAs in Trypanosoma brucei and Leptomonas collosoma: a flexible structure that is independent of cap4 methylations.

The formation of the mRNA 5' end in trypanosomatid protozoa is carried out by trans-splicing, which transfers a spliced leader (SL) sequence and its hypermethylated cap (cap4) from the SL RNA to the pre-mRNA. Previous in vitro studies with synthetic uncapped RNAs have shown that the SL sequence of Leptomonas collosoma can assume two alternate conformations, Form 1 and Form 2, with Form 1 being the dominant one. To gain information about the structure of the SL RNA in vivo, in its protein-rich environment, we have used permeable Trypanosoma brucei and L. collosoma cells for chemical modification experiments. We introduce the use in vivo of the water-soluble reagents CMCT and kethoxal. In contrast to the in vitro results, the Form 2 secondary structure predominates. However, there are chemically accessible regions that suggest conformational flexibility in SL RNPs and a chemically inaccessible region suggestive of protection by protein or involvement in tertiary interactions. Using complementary 2'-O-methyl RNA oligonucleotides, we show that T. brucei SL RNA can be induced to switch conformation in vivo. SL RNA stripped of proteins and probed in vitro does not display the same Form 2 bias, indicating that SL RNA structure is determined, at least in part, by its RNP context. Finally, the methyl groups of the cap4 do not seem to affect the secondary structure of T. brucei SL RNA, as shown by chemical modification of undermethylated SL RNA probed in vivo.