Structural differences within the loop E motif imply alternative mechanisms of viroid processing

Viroids replicate via a rolling circle mechanism, and cleavage/ligation requires extensive rearrangement of the highly base-paired native structure. For Potato spindle tuber viroid (PSTVd), the switch from cleavage to ligation is driven by the change from a multibranched tetraloop structure to a loop E conformation. Here we present evidence that processing of Citrus viroid III (CVd-III), a member of a related group of viroids that also replicate in the nucleus, may proceed via a distinct pathway. Chemical probing of PSTVd and CVd-III miniRNAs with DMS and CMCT revealed that the loop E motifs of these two viroids have quite different tertiary structures. As shown by temperature gradient gel electrophoresis, the presence of two likely Watson-Crick GC pairs results in a significant overall stabilization of the CVd-III loop E-like motif. Unlike PSTVd, the upper strand of the CVd-III loop E-like motif cannot fold into a GNRA tetraloop, and comparison of suboptimal structures indicates that the initial cleavage event could occur on the 5' side of the only GU wobble pair in a helix involving a nearby pair of inverted repeats. According to our model, rearrangement of 3' sequences into a hairpin stem containing an identical arrangement of GC, GU, and CG base pairs and a second cleavage event is followed by formation of loop E, which serves to align the 5' and 3' termini of the CVd-III monomer prior to ligation. Because ligation would occur within loop E itself, stabilization of this motif may be needed to hold the 5' and 3' termini of CVd-III in position for the host ligase.     

Automated classification of RNA 3D motifs and the RNA 3D Motif Atlas

The analysis of atomic-resolution RNA three-dimensional (3D) structures reveals that many internal and hairpin loops are modular, recurrent, and structured by conserved non-Watson–Crick base pairs. Structurally similar loops define RNA 3D motifs that are conserved in homologous RNA molecules, but can also occur at nonhomologous sites in diverse RNAs, and which often vary in sequence. To further our understanding of RNA motif structure and sequence variability and to provide a useful resource for structure modeling and prediction, we present a new method for automated classification of internal and hairpin loop RNA 3D motifs and a new online database called the RNA 3D Motif Atlas. To classify the motif instances, a representative set of internal and hairpin loops is automatically extracted from a nonredundant list of RNA-containing PDB files. Their structures are compared geometrically, all-against-all, using the FR3D program suite. The loops are clustered into motif groups, taking into account geometric similarity and structural annotations and making allowance for a variable number of bulged bases. The automated procedure that we have implemented identifies all hairpin and internal loop motifs previously described in the literature. All motif instances and motif groups are assigned unique and stable identifiers and are made available in the RNA 3D Motif Atlas (http://rna.bgsu.edu/motifs), which is automatically updated every four weeks. The RNA 3D Motif Atlas provides an interactive user interface for exploring motif diversity and tools for programmatic data access.     

Characterization of the RNA motif responsible for the specific interaction of potato spindle tuber viroid RNA (PSTVd) and the tomato protein Virp1

Viroids are small non-coding parasitic RNAs that are able to infect their host plants systemically. This circular naked RNA makes use of host proteins to accomplish its proliferation. Here we analyze the specific binding of the tomato protein Virp1 to the terminal right domain of potato spindle tuber viroid RNA (PSTVd). We find that two asymmetric internal loops within the PSTVd (+) RNA, each composed of the sequence elements 5′-ACAGG and CUCUUCC-5′, are responsible for the specific RNA–protein interaction. In view of the nucleotide composition we call this structural element an ‘RY motif’. The RY motif located close to the terminal right hairpin loop of the PSTVd secondary structure has an ~5-fold stronger binding affinity than the more centrally located RY motif. Simultaneous sequence alterations in both RY motifs abolished the specific binding to Virp1. Mutations in any of the two RY motifs resulted in non-infectious viroid RNA, with the exception of one case, where reversion to sequence wild type took place. In contrast, the simultaneous exchange of two nucleotides within the terminal right hairpin loop of PSTVd had only moderate influence on the binding to Virp1. This variant was infectious and sequence changes were maintained in the progeny. The relevance of the phylogenetic conservation of the RY motif, and sequence elements therein, amongst various genera of the family Pospiviroidae is discussed.     

Direct role of a viroid RNA motif in mediating directional RNA trafficking across a specific cellular boundary

The plasmodesmata and phloem form a symplasmic network that mediates direct cell-cell communication and transport throughout a plant. Selected endogenous RNAs, viral RNAs, and viroids traffic between specific cells or organs via this network. Whether an RNA itself has structural motifs to potentiate trafficking is not well understood. We have used mutational analysis to identify a motif that the noncoding Potato spindle tuber viroid RNA evolved to potentiate its efficient trafficking from the bundle sheath into mesophyll that is vital to establishing systemic infection in tobacco (Nicotiana tabacum). Surprisingly, this motif is not necessary for trafficking in the reverse direction (i.e., from the mesophyll to bundle sheath). It is not required for trafficking between other cell types either. We also found that the requirement for this motif to mediate bundle sheath-to-mesophyll trafficking is dependent on leaf developmental stages. Our results provide genetic evidence that (1) RNA structural motifs can play a direct role in mediating trafficking across a cellular boundary in a defined direction, (2) the bundle sheath-mesophyll boundary serves as a novel regulatory point for RNA trafficking between the phloem and nonvascular tissues, and (3) the symplasmic network remodels its capacity to traffic RNAs during plant development. These findings may help further studies to elucidate the interactions between RNA motifs and cellular factors that potentiate directional trafficking across specific cellular boundaries.     

Tertiary structure and function of an RNA motif required for plant vascular entry to initiate systemic trafficking

Vascular entry is a decisive step for the initiation of long-distance movement of infectious and endogenous RNAs, silencing signals and developmental/defense signals in plants. However, the mechanisms remain poorly understood. We used Potato spindle tuber viroid (PSTVd) as a model to investigate the direct role of the RNA itself in vascular entry. We report here the identification of an RNA motif that is required for PSTVd to traffic from nonvascular into the vascular tissue phloem to initiate systemic infection. This motif consists of nucleotides U/C that form a water-inserted cis Watson-Crick/Watson-Crick base pair flanked by short helices that comprise canonical Watson-Crick/Watson-Crick base pairs. This tertiary structural model was inferred by comparison with X-ray crystal structures of similar motifs in rRNAs and is supported by combined mutagenesis and covariation analyses. Hydration pattern analysis suggests that water insertion induces a widened minor groove conducive to protein and/or RNA interactions. Our model and approaches have broad implications to investigate the RNA structural motifs in other RNAs for vascular entry and to study the basic principles of RNA structure-function relationships.     

Secondary structure model of the Mason-Pfizer monkey virus 5' leader sequence: identification of a structural motif common to a variety of retroviruses.

A stable secondary structure model is presented for the region 3' of the primer-binding site to 130 bases into the gag sequence of the prototype type D retrovirus Mason-Pfizer monkey virus. Using biochemical probing of RNA from this region in association with free energy minimization, we have identified a stem-loop structure in the region, which from other studies has been shown to be important for genomic RNA encapsidation. The structure involves a highly stable stem of five G-C pairs terminating in a heptaloop. Comparison of the Mason-Pfizer monkey virus structure with one predicted for squirrel monkey retrovirus demonstrates an identical stem and a common ACC motif in the loop. Free energy studies of the secondary structure of the 5' regions of eight other retroviruses predict stem loops which have similar GAYC motifs. We believe this may represent a common structural and sequence motif which among other functions may be involved in genomic RNA packaging in these viruses.     

Identifying novel sequence variants of RNA 3D motifs

Predicting RNA 3D structure from sequence is a major challenge in biophysics. An important sub-goal is accurately identifying recurrent 3D motifs from RNA internal and hairpin loop sequences extracted from secondary structure (2D) diagrams. We have developed and validated new probabilistic models for 3D motif sequences based on hybrid Stochastic Context-Free Grammars and Markov Random Fields (SCFG/MRF). The SCFG/MRF models are constructed using atomic-resolution RNA 3D structures. To parameterize each model, we use all instances of each motif found in the RNA 3D Motif Atlas and annotations of pairwise nucleotide interactions generated by the FR3D software. Isostericity relations between non-Watson–Crick basepairs are used in scoring sequence variants. SCFG techniques model nested pairs and insertions, while MRF ideas handle crossing interactions and base triples. We use test sets of randomly-generated sequences to set acceptance and rejection thresholds for each motif group and thus control the false positive rate. Validation was carried out by comparing results for four motif groups to RMDetect. The software developed for sequence scoring (JAR3D) is structured to automatically incorporate new motifs as they accumulate in the RNA 3D Motif Atlas when new structures are solved and is available free for download.     

Nuclear magnetic resonance spectroscopy and molecular modeling reveal that different hydrogen bonding patterns are possible for G.U pairs: one hydrogen bond for each G.U pair in r(GGCGUGCC)(2) and two for each G.U pair in r(GAGUGCUC)(2)

G.U pairs occur frequently and have many important biological functions. The stability of symmetric tandem G.U motifs depends both on the adjacent Watson-Crick base pairs, e.g., 5'G > 5'C, and the sequence of the G.U pairs, i.e., 5'-UG-3' > 5'-GU-3', where an underline represents a nucleotide in a G.U pair [Wu, M., McDowell, J. A., and Turner, D. H. (1995) Biochemistry 34, 3204-3211]. In particular, at 37 degrees C, the motif 5'-CGUG-3' is less stable by approximately 3 kcal/mol compared with other symmetric tandem G.U motifs with G-C as adjacent pairs: 5'-GGUC-3', 5'-GUGC-3', and 5'-CUGG-3'. The solution structures of r(GAGUGCUC)(2) and r(GGCGUGCC)(2) duplexes have been determined by NMR and restrained simulated annealing. The global geometry of both duplexes is close to A-form, with some distortions localized in the tandem G.U pair region. The striking discovery is that in r(GGCGUGCC)(2) each G.U pair apparently has only one hydrogen bond instead of the two expected for a canonical wobble pair. In the one-hydrogen-bond model, the distance between GO6 and UH3 is too far to form a hydrogen bond. In addition, the temperature dependence of the imino proton resonances is also consistent with the different number of hydrogen bonds in the G.U pair. To test the NMR models, U or G in various G.U pairs were individually replaced by N3-methyluridine or isoguanosine, respectively, thus eliminating the possibility of hydrogen bonding between GO6 and UH3. The results of thermal melting studies on duplexes with these substitutions support the NMR models.     

Predicting U-turns in ribosomal RNA with comparative sequence analysis

The U-turn is a well-known RNA motif characterized by a sharp reversal of the RNA backbone following a single-stranded uridine base. In experimentally determined U-turn motifs, the nucleotides 3' to the turn are frequently involved in tertiary interactions, rendering this motif particularly attractive in RNA modeling and functional studies. The U-turn signature is composed of an UNR sequence pattern flanked by a Y:Y, Y:A (Y=pyrimidine) or G:A base juxtaposition. We have identified 33 potential UNR-type U-turns and 25 related GNRA-type U-turns in a large set of aligned 16 S and 23 S rRNA sequences. U-turn candidates occur in hairpin loops (34 times) as well as in internal and multi-stem loops (24 times). These are classified into ten families based on loop type, sequence pattern (UNR or GNRA) and the nature of the closing base juxtaposition. In 13 cases, the bases on the 3' side of the turn, or on the immediate 5' side, are involved in tertiary covariations, making these sites strong candidates for tertiary interactions.     

A novel RNA structural motif in the selenocysteine insertion element of eukaryotic selenoprotein mRNAs.

In eukaryotes, co-translational insertion of selenocysteine into selenoproteins necessitates the participation of the selenocysteine insertion sequence (SECIS), an element lying in the 3'-untranslated region of selenoprotein mRNAs. We report a detailed experimental study of the secondary structures of the SECIS elements of three selenoprotein mRNAs, the rat and human type I iodothyronine deiodinase (5'DI) and rat glutathione peroxidase (GPx). Based on RNase and chemical probing, a new secondary structure model is established. It is characterized by a stem-loop structure, comprising two helices (I and II) separated by an internal loop, with an apical loop surmounting helix II. Sequence comparisons of 20 SECIS elements, arising from 2 5'DI, 13 GPx, 2 selenoprotein P, and 1 selenoprotein W mRNAs, confirm the secondary structure model. The most striking finding of the experimental study concerns a set of conserved sequences in helix II that interact to form a novel RNA structural motif consisting of a quartet composed of non-Watson-Crick base pairs 5'UGAY3': 5'UGAU3'. The potential for forming the quartet is preserved in 15 SECIS elements, but three consecutive non-Watson-Crick base pairs can nevertheless form in the other five SECIS, the central G.A tandem being invariant in all cases. A 3D model, derived by computer modeling with the use of the solution data, suggests that the base pairing interactions in the G.A tandem are of the type found in GNRA loops. The 3D model displays the quartet lying in an accessible position at the foot of helix II, which is bent at the internal loop, suggesting that the non-Watson-Crick base pair arrangement provides an unusual pattern of chemical groups for putative ligand interaction.     

Zipper-like Watson-Crick base-pairs

A series of DNA heptadecamers containing the DNA analogues of RNA E-like 5'-d(GXA)/(AYG)-5' motifs (X/Y is complementary T/A, A/T, C/G, or G/C pair) were studied using nuclear magnetic resonance (NMR) methodology and distance geometry (DG)/molecular dynamics (MD) approaches. Such oligomers reveal excellent resolution in NMR spectra and exhibit many unusual nuclear Overhauser effects (NOEs) that allow for good characterization of an unusual zipper-like conformation with zipper-like Watson-Crick base-pairs; the potential canonical X.Y H-bonding is not present, and the central X/Y pairs are transformed instead into inter-strand stacks that are bracketed by sheared G.A base-pairs. Such phenomenal structural change is brought about mainly through two backbone torsional angle adjustments, i.e. delta from C2'-endo to C3'-endo for the sugar puckers of unpaired residues and gamma from gauche(+) to trans for the following 3'-adenosine residues. Such motifs are analogous to the previously studied (GGA)(2) motif presumably present in the human centromeric (TGGAA)(n) tandem repeat sequence. The novel zipper-like motifs are only 4-7 deg. C less stable than the (GGA)(2) motif, suggesting that inter-strand base stacking plays an important role in stabilizing unusual nucleic acid structures. The discovery that canonical Watson-Crick G.C or A.T hydrogen-bonded pairs can be transformed into stacking pairs greatly increases the repertoire for unusual nucleic acid structural motifs.     

Molecular dynamics study of intrinsic stability in six RNA terminal loop motifs

Single stranded RNA molecules can assume a wide range of tertiary structures beyond the canonical A-form double helix. Certain sequences, termed motifs, are more common than a random distribution would suggest. The existence of such motifs can be rationalized in structural terms. In this study, we have investigated the intrinsic structural stability of RNA terminal loop motifs using multiple MD simulations in explicit water. Representative loops were chosen from the major tetraloop motifs, including also the U-turn motif. Not all loops retain their folded starting structure, but lowering the temperature to 277 K, or adding adjacent base pairs from the stem to which the motif is attached, helps stabilizing the folded loop structure.     

The crystal structure at 1.5 angstroms resolution of an RNA octamer duplex containing tandem G.U basepairs

The crystal structure of the RNA octamer, 5'-GGCGUGCC-3' has been determined from x-ray diffraction data to 1.5 angstroms resolution. In the crystal, this oligonucleotide forms five self-complementary double-helices in the asymmetric unit. Tandem 5'GU/3'UG basepairs comprise an internal loop in the middle of each duplex. The NMR structure of this octameric RNA sequence is also known, allowing comparison of the variation among the five crystallographic duplexes and the solution structure. The G.U pairs in the five duplexes of the crystal form two direct hydrogen bonds and are stabilized by water molecules that bridge between the base of guanine (N2) and the sugar (O2') of uracil. This contrasts with the NMR structure in which only one direct hydrogen bond is observed for the G.U pairs. The reduced stability of the r(CGUG)2 motif relative to the r(GGUC)2 motif may be explained by the lack of stacking of the uracil bases between the Watson-Crick and G.U pairs as observed in the crystal structure.     

Context-sensitivity of isosteric substitutions of non-Watson–Crick basepairs in recurrent RNA 3D motifs

Sequence variation in a widespread, recurrent, structured RNA 3D motif, the Sarcin/Ricin (S/R), was studied to address three related questions: First, how do the stabilities of structured RNA 3D motifs, composed of non-Watson–Crick (non-WC) basepairs, compare to WC-paired helices of similar length and sequence? Second, what are the effects on the stabilities of such motifs of isosteric and non-isosteric base substitutions in the non-WC pairs? And third, is there selection for particular base combinations in non-WC basepairs, depending on the temperature regime to which an organism adapts? A survey of large and small subunit rRNAs from organisms adapted to different temperatures revealed the presence of systematic sequence variations at many non-WC paired sites of S/R motifs. UV melting analysis and enzymatic digestion assays of oligonucleotides containing the motif suggest that more stable motifs tend to be more rigid. We further found that the base substitutions at non-Watson–Crick pairing sites can significantly affect the thermodynamic stabilities of S/R motifs and these effects are highly context specific indicating the importance of base-stacking and base-phosphate interactions on motif stability. This study highlights the significance of non-canonical base pairs and their contributions to modulating the stability and flexibility of RNA molecules.     

A potential RNA drug target in the hepatitis C virus internal ribosomal entry site

Subdomain IlId from the hepatitis C virus (HCV) internal ribosome entry site (IRES) has been shown to be essential for cap-independent translation. We have conducted a structural study of a 27-nt fragment, identical in sequence to IlId, to explore the structural features of this subdomain. The proposed secondary structure of IlId is comprised of two 3 bp helical regions separated by an internal loop and closed at one end by a 6-nt terminal loop. NMR and molecular modeling were used interactively to formulate a validated model of the three-dimensional structure of IlId. We found that this fragment contains several noncanonical structural motifs and non-Watson-Crick base pairs, some of which are common to other RNAs. In particular, a motif characteristic of the rRNA alpha-sarcin/ricin loop was located in the internal loop. The terminal loop, 5'-UUGGGU, was found to fold to form a trinucleotide loop closed by a trans-wobble U.G base pair. The sixth nucleotide was bulged out to allow stacking of this U.G pair on the adjacent helical region. In vivo mutational analysis in the context of the full IRES confirmed the importance of each structural motif within IIId for IRES function. These findings may provide clues as to host cellular proteins that play a role in IRES-directed translation and, in particular, the mechanism through which host ribosomes are sequestered for viral function.     

Consecutive GA pairs stabilize medium-size RNA internal loops

Internal loops in RNA are important for folding and function. Consecutive noncanonical pairs can form in internal loops having at least two nucleotides on each side. Thermodynamic and structural insights into such internal loops should improve approximations for their stabilities and predictions of secondary and three-dimensional structures. Most natural internal loops are purine rich. A series of oligoribonucleotides that form purine-rich internal loops of 5-10 nucleotides, including kink-turn loops, were studied by UV melting, exchangeable proton and phosphorus NMR. Three consecutive GA pairs with the motif 5' Y GGA/3' R AAG or GGA R 3'/AAG Y 5' (i.e., 5' GGA 3'/3' AAG 5' closed on at least one side with a CG, UA, or UG pair with Y representing C or U and R representing A or G) stabilize internal loops having 6-10 nucleotides. Certain motifs with two consecutive GA pairs are also stabilizing. In internal loops with three or more nucleotides on each side, the motif 5' U G/3' G A has stability similar to 5' C G/3' G A. A revised model for predicting stabilities of internal loops with 6-10 nucleotides is derived by multiple linear regression. Loops with 2 x 3 nucleotides are predicted well by a previous thermodynamic model.