The substitution U<Subscript>475</Subscript> → C with Sabin3-like mutation within the IRES attenuate Coxsackievirus B<Subscript>3</Subscript> cardiovirulence

The Sabin3 mutation in the viral RNA plays an important role in directing attenuation phenotype of Sabin vaccine strain of poliovirus type 1 (PV1). We previously described that Sabin3-like mutation introduced in Coxsackievirus B3 (CVB3) genome led to a defective mutant. However, this mutation do not led to destruction of secondary structure motif C within the stem-loop V of CVB3 RNA because of the presence of one nucleotide difference (C → U) in the region encompassing the Sabin3 mutation at nucleotides 471 of PV1 and 475 of CVB3 RNA. In order to reproduce the same sequence of PV1 sabin3 vaccine strain, we introduce in this study an additional mutation (U₄₇₅ → C) to CVB3 Sabin3-like mutant. Our results demonstrated that Sabin3-like+C mutant displayed a decreased translation initiation defects when translated in cell-free system. This translation initiation defect was correlated with reduced yields of infectious virus particles in HeLa cells in comparison with Sabin3-like mutant and wild-type CVB3 viruses. Inoculation of Swiss mice with mutant viruses resulted in no inflammatory heart disease when compared to heart of mice infected with wild-type. Theses findings indicate that the double mutant could be exploited for the development of a live attenuated vaccine against CVB3.     

Genomic determinants of cardiovirulence in coxsackievirus B3 clinical isolates: localization to the 5' nontranslated region

Coxsackievirus B3 (CVB3) infections can cause myocarditis in humans and are implicated in the pathogenesis of dilated cardiomyopathy. The natural genetic determinants of cardiovirulence for CVB3 have not been identified, although using strains engineered in the laboratory, cardiovirulence determinants have been identified in the CVB3 5' nontranslated region (5'NTR) and capsid. The myocarditic phenotypes of two CVB3 clinical isolates were determined using an established murine model of inflammatory heart disease. The 5'NTRs and capsid proteins of the noncardiovirulent CVB3/CO strain and cardiovirulent CVB3/AS strain were examined to determine their influence on the cardiovirulence phenotype. Six intratypic chimeric viruses were constructed in which 5'NTR and capsid sequences of the infectious cDNA copy of the cardiovirulent CVB3/20 genome were replaced by homologous sequences from CVB3/CO or CVB3/AS. Chimeric strains were tested for cardiovirulence by inoculation of C3H/HeJ mice. Sections of hearts removed at 10 days postinoculation were examined for evidence of myocarditis by light microscopy and assayed for the presence of virus. Replacement of the CVB3/20 capsid coding region by that from the homologous region of CVB3/CO resulted in no change in the cardiovirulent CVB3/20 phenotype, with virus recoverable from the heart at 10 days postinoculation. However, recombinant virus containing the CVB3/CO 5'NTR alone or the 5'NTR and capsid sequences together were not myocarditic, and infectious virus was not recovered from the myocardium. Chimeric viruses containing the CVB3/AS 5'NTR alone, capsid sequence alone, or both together preserved the myocarditic phenotype. These data support the 5'NTR as the primary site in the determination of the natural cardiovirulence phenotype of CVB3.     

The Pro/Hel region is indispensable for packaging non-replicating turnip yellow mosaic virus RNA, but not replicating viral RNA

Turnip yellow mosaic virus (TYMV) is a spherical plant virus that has a single 6.3 kb positive strand RNA. The genomic RNA has a tRNA-like structure (TLS) at the 3′-end. The 3′-TLS and hairpins in the 5′-untranslated region supposedly serve as packaging signals; however, recent studies have shown that they do not play a role in TYMV RNA packaging. In this study, we focused on packaging signals by examining a series of deletion mutants of TYMV. Analysis of encapsidated viral RNA after agroinfiltration of the deletion constructs into Nicotiana benthamiana showed that the mutant RNA lacking the protease (Pro)/helicase (Hel) region was not encapsidated by the coat proteins provided in trans, implicating that a packaging signal lies in the Pro/Hel region. Examination of two Pro⁻Hel⁻ mutants showed that protein activity from the Pro/Hel domains was dispensable for the packaging of the non-replicating TYMV RNA. In contrast, the mutant TYMV RNA lacking the Pro/Hel region was efficiently encapsidated when the mutant TYMV was co-introduced with a wild-type TYMV, suggesting that packaging mechanisms might differ depending on whether the virus is replicating or not.     

Role of RNA Structure Motifs in IRES-Dependent Translation Initiation of the Coxsackievirus B3: New Insights for Developing Live-Attenuated Strains for Vaccines and Gene Therapy

Internal ribosome entry site (IRES) elements are highly structured RNA sequences that function to recruit ribosomes for the initiation of translation. In contrast to the canonical cap-binding, the mechanism of IRES-mediated translation initiation is still poorly understood. Translation initiation of the coxsackievirus B3 (CVB3), a causative agent of viral myocarditis, has been shown to be mediated by a highly ordered structure of the 5′ untranslated region (5′UTR), which harbors an IRES. Taking into account that efficient initiation of mRNA translation depends on temporally and spatially orchestrated sequence of RNA–protein and RNA–RNA interactions, and that, at present, little is known about these interactions, we aimed to describe recent advances in our understanding of molecular structures and biochemical functions of the translation initiation process. Thus, this review will explore the IRES elements as important RNA structures and the significance of these structures in providing an alternative mechanism of translation initiation of the CVB3 RNA. Since translation initiation is the first intracellular step during the CVB3 infection cycle, the IRES region provides an ideal target for antiviral therapies. Interestingly, the 5′ and 3′UTRs represent promising candidates for the study of CVB3 cardiovirulence and provide new insights for developing live-attenuated vaccines.     

Requirements for RNA Replication of a Poliovirus Replicon by Coxsackievirus B3 RNA Polymerase

A chimeric poliovirus type 1 (PV1) genome was constructed in which the 3D RNA polymerase (3D(pol)) coding sequences were replaced with those from coxsackievirus B3 (CVB3). No infectious virus was produced from HeLa cells transfected with the chimeric RNA. Processing of the PV1 capsid protein precursor was incomplete, presumably due to inefficient recognition of the P1 protein substrate by the chimeric 3CD proteinase containing CVB3 3D sequences. The ability of the chimeric RNA to replicate in the absence of capsid formation was measured after replacement of the P1 region with a luciferase reporter gene. No RNA synthesis was detected, despite efficient production of enzymatically active 3D(pol) from the 3D portion of the chimeric 3CD. The chimeric 3CD protein was unable to efficiently bind to the cloverleaf-like structure (CL) at the 5' end of PV1 RNA, which has been demonstrated previously to be required for viral RNA synthesis. The CVB3 3CD protein bound the PV1 CL as well as PV1 3CD. An additional chimeric PV1 RNA that contained CVB3 3CD sequences also failed to produce virus after transfection. Since processing of PV1 capsid protein precursors by the CVB3 3CD was again incomplete, a luciferase-containing replicon was also analyzed for RNA replication. The 3CD chimera replicated at 33 degrees C, but not at 37 degrees C. Replacement of the PV1 5'-terminal CL with that of CVB3 did not rescue the temperature-sensitive phenotype. Thus, there is an essential interaction(s) between 3CD and other viral P2 or P3 protein products required for efficient RNA replication which is not fully achieved between proteins from the two different members of the same virus genus.     

Role of GNRA motif mutations within stem-loop V of internal ribosome entry segment in coxsackievirus B3 molecular attenuation

The lengthy 5' nontranslated region of coxsackievirus B3 (CVB3) forms a highly ordered secondary structure containing an internal ribosome entry segment (IRES), which plays an important role in controlling viral translation and pathogenesis. The stem-loop V (SL-V) of this IRES contains a large lateral bulge loop which encompasses two conserved GNRA motifs. In this study, we analyzed the effects of point mutations within the GNRA motifs of the CVB3 IRES. We characterized in vitro virus production and translation efficiency and we tested in vivo virulence of two CVB3 mutants produced by site-directed mutagenesis. The GNAA1 and GNAA2 RNAs displayed decreased translation initiation efficiency when translated in rabbit reticulocyte lysates. This translation defect was correlated with reduced yields of infectious virus particles in HeLa cells in comparison with the wild type. When inoculated orally into Swiss mice, both mutant viruses were avirulent and caused neither inflammation nor necrosis in hearts. These results highlight the important role of the GNRA motifs within the SL-V of the IRES of CVB3, in directing translation initiation.     

Extending the family of UNCG-like tetraloop motifs: NMR structure of a CACG tetraloop from coxsackievirus B3

Stable RNA tetraloop motifs are found frequently in biologically active RNAs. These motifs carry out a wide variety of functions in RNA folding, in RNA-RNA and RNA-protein interactions. A great deal of knowledge about the structures and functions of tetraloop motifs has accumulated largely due to intensive theoretical, biochemical, and biophysical studies on three most frequently occurring families of tetraloop sequences, namely, the cUNCGg, the cGNRAg, and the gCUUGc sequences. Our knowledge surely is not exhaustive, and efforts are still being made to gain a better understanding. Here we report the NMR structure of a uCACGg tetraloop that occurs naturally within the cloverleaf RNA structure of the 5'-UTR of coxsackievirus B3. This tetraloop is the major determinant for interaction between the cloverleaf RNA and viral 3C protease, which is an essential part of a ribonucleoprotein complex that plays a critical role in the regulation of viral translation and replication. Our structure shows that the CACG tetraloop is closed by a wobble U.G base pair. The structure of the CACG tetraloop is stabilized by extensive base stacking and hydrogen bonding interactions strikingly similar to those previously reported for the cUUCGg tetraloop. Identification of these hallmark structural features strongly supports the existence of an extended YNCG tetraloop family. The U.G base pair closing the stem and the A residue in the loop introduce some small structural and themodynamic distinctions from the canonical cUUCGg tetraloop that may be important for recognition by the viral 3C protease.     

Ribosomal RNA identity elements for ricin A-chain recognition and catalysis. Analysis with tetraloop mutants.

Ricin is a cytotoxic protein that inactivates ribosomes by hydrolyzing the N-glycosidic bond between the base and the ribose of the adenosine at position 4324 in eukaryotic 28 S rRNA. Ricin A-chain will also catalyze depurination in naked prokaryotic 16 S rRNA; the adenosine is at position 1014 in a GAGA tetraloop. The rRNA identity elements for recognition by ricin A-chain and for the catalysis of cleavage were examined using synthetic GAGA tetraloop oligoribonucleotides. The RNA designated wild-type, an oligoribonucleotide (19-mer) that approximates the structure of the ricin-sensitive site in 16 S rRNA, and a number of mutants were transcribed in vitro from synthetic DNA templates with phage T7 RNA polymerase. With the wild-type tetraloop oligoribonucleotide the ricin A-chain-catalyzed reaction has a Km of 5.7 microM and a Kcat of 0.01 min-1. The toxin alpha-sarcin, which cleaves the phosphodiester bond on the 3' side of G4325 in 28 S rRNA, does not recognize the tetraloop RNA, although alpha-sarcin does affect a larger synthetic oligoribonucleotide that has a 17-nucleotide loop with a GAGA sequence; thus, there is a clear divergence in the identity elements for the two toxins. Mutants were constructed with all of the possible transitions and transversions of each nucleotide in the GAGA tetraloop; none was recognized by ricin A-chain. Thus, there is an absolute requirement for the integrity of the GAGA sequence in the tetraloop. The helical stem of the tetraloop oligoribonucleotide can be reduced to three base-pairs, indeed, to two base-pairs if the temperature is decreased, without affecting recognition; the nature of these base-pairs does not influence recognition or catalysis by ricin A-chain. If the tetraloop is opened so as to form a GAGA-containing hexaloop, recognition by ricin A-chain is lost. This suggests that during the elongation cycle, a GAGA tetraloop either exists or is formed in the putative 17-member single-stranded region of the ricin domain in 28 S rRNA and this bears on the mechanism of protein synthesis.     

Nucleotide sequence of an attenuated mutant of coxsackievirus B3 compared with the cardiovirulent wildtype: assessment of candidate mutations by analysis of a revertant to cardiovirulence

Background: Coxsackievirus B3 (CVB3) causes myocarditis in the SWR (H^2q) mouse model and persistence of CVB3 in myocardium disposes to the development of dilated cardiomyopathy. An attenuated strain of CVB3 has been isolated, sequenced and several candidate mutations for attenuation identified. Derivation of a revertant to cardiovirulence allows the significance of these mutations to be assessed.Objectives: To ascertain which candidate mutation(s) determine(s) the attenuated phenotype.Study design: A revertant to cardiovirulence was isolated following passage through severe combined immunodeficient disease (SCID) mouse heart. The 5′-non-translated region (NTR) and region coding for capsid proteins were sequenced and compared to the wildtype and attenuant.Results: There are five candidates for attenuation: (1) A–G at base 580 in the 5′-NTR; (2) A–T at base 690 in the 5′-NTR; (3) CG–GC at bases 1401/2 (Thr to Ser at amino acid 151 in VP2); (4) AA–GT at bases 2691/2 (Lys to Ser at amino acid 80 in VP1); (5) A–G at base 2916 (Asp to Gly at amino acid 155 in VP1). It was shown previously that mutations at 580, 690 and 2691/2 are not important in attenuation. Additionally, there are three novel mutations in the coding region of the revertant and one in the 5′-NTR which are unlikely to be relevant for attenuation as they are not present in the attenuant. Of nucleotide changes seen at 1401/2 and 2916 in the attenuant, only 2916 reverts to the wildtype sequence and so is a strong candidate for a determinant of attenuation.Conclusions: The A–G mutation at 2916 (Asp to Gly at amino acid 155 in VP1) is a strong candidate for attenuation. It is located at the top of the receptor binding cleft and mutation of the Asp to a Gly may destabilise the receptor binding site.     

Temperature-Sensitive Mutants of Murine Leukemia Virus IV. Further Physiological Characterization and Evidence for Genetic Recombination

Several temperature-sensitive mutants of the Rauscher strain of murine leukemia virus representing three distinct physiological groups have been further characterized. Genetic recombination between different pairs of these mutants has been demonstrated. Several representative genetic recombinants were isolated and shown to replicate equally well at the permissive (31 C) and nonpermissive (38 C) temperatures and to have serological characteristics of the wild-type parental virus. Alternative models for the mechanisms involved in recombination between type C RNA viruses are discussed.     

The nonstructural protein 2C of Coxsackie B virus has RNA helicase and chaperoning activities

RNA-remodeling proteins, including RNA helicases and chaperones, play vital roles in the remodeling of structured RNAs. During viral replication, viruses require RNA-remodeling proteins to facilitate proper folding and/or re-folding the viral RNA elements. Coxsackieviruses B3 (CVB3) and Coxsackieviruses B5 (CVB5), belonging to the genus Enterovirus in the family Picornaviridae, have been reported to cause various infectious diseases such as hand-foot-and-mouth disease, aseptic meningitis, and viral myocarditis. However, little is known about whether CVB3 and CVB5 encode any RNA remodeling proteins. In this study, we showed that 2C proteins of CVB3 and CVB5 contained the conserved SF3 helicase A, B, and C motifs, and functioned not only as RNA helicase that unwound RNA helix bidirectionally in an NTP-dependent manner, but also as RNA chaperone that remodeled structured RNAs and facilitated RNA strand annealing independently of NTP. In addition, we determined that the NTPase activity and RNA helicase activity of 2C proteins of CVB3 and CVB5 were dependent on the presence of divalent metallic ions. Our findings demonstrate that 2C proteins of CVBs possess RNA-remodeling activity and underline the functional importance of 2C protein in the life cycle of CVBs.     

In Vitro Molecular Characterization of RNA–Proteins Interactions During Initiation of Translation of a Wild-Type and a Mutant Coxsackievirus B3 RNAs

Translation initiation of Coxsackievirus B3 (CVB3) RNA is directed by an internal ribosome entry site (IRES) within the 5′ untranslated region. Host cell factors involved in this process include some canonical translation factors and additional RNA-binding proteins. We have, previously, described that the Sabin3-like mutation (U475 → C) introduced in CVB3 genome led to a defective mutant with a serious reduction in translation efficiency. With the aim to identify proteins interacting with CVB3 wild-type and Sabin3-like IRESes and to study interactions between HeLa cell or BHK-21 protein extracts and CVB3 RNAs, UV-cross-linking assays were performed. We have observed a number of proteins that specifically interact with both RNAs. In particular, molecular weights of five of these proteins resemble to those of the eukaryotic translation initiation factors 4G, 3b, 4B, and PTB. According to cross-linking patterns obtained, we have demonstrated a better affinity of CVB3 RNA binding to BHK-21 proteins and a reduced interaction of the mutant RNA with almost cellular polypeptides compared to the wild-type IRES. On the basis of phylogeny of some initiation factors and on the knowledge of the initiation of translation process, we focused on the interaction of both IRESes with eIF3, p100 (eIF4G), and 40S ribosomal subunit by filter-binding assays. We have demonstrated a better affinity of binding to the wild-type CVB3 IRES. Thus, the reduction efficiency of the mutant RNA to bind to cellular proteins involved in the translation initiation could be the reason behind inefficient IRES function.     

5'-Terminal deletions occur in coxsackievirus B3 during replication in murine hearts and cardiac myocyte cultures and correlate with encapsidation of negative-strand viral RNA

Adult human enteroviral heart disease is often associated with the detection of enteroviral RNA in cardiac muscle tissue in the absence of infectious virus. Passage of coxsackievirus B3 (CVB3) in adult murine cardiomyocytes produced CVB3 that was noncytolytic in HeLa cells. Detectable but noncytopathic CVB3 was also isolated from hearts of mice inoculated with CVB3. Sequence analysis revealed five classes of CVB3 genomes with 5' termini containing 7, 12, 17, 30, and 49 nucleotide deletions. Structural changes (assayed by chemical modification) in cloned, terminally deleted 5'-nontranslated regions were confined to the cloverleaf domain and localized within the region of the deletion, leaving key functional elements of the RNA intact. Transfection of CVB3 cDNA clones with the 5'-terminal deletions into HeLa cells generated noncytolytic virus (CVB3/TD) which was neutralized by anti-CVB3 serum. Encapsidated negative-strand viral RNA was detected using CsCl-purified CVB3/TD virions, although no negative-strand virion RNA was detected in similarly treated parental CVB3 virions. The viral protein VPg was detected on CVB3/TD virion RNA molecules which terminate in 5' CG or 5' AG. Detection of viral RNA in mouse hearts from 1 week to over 5 months postinoculation with CVB3/TD demonstrated that CVB3/TD virus strains replicate and persist in vivo. These studies describe a naturally occurring genomic alteration to an enteroviral genome associated with long-term viral persistence.     

A Predicted Secondary Structural Domain within the Internal Ribosome Entry Site of Echovirus 12 Mediates a Cell-Type-Specific Block to Viral Replication

The enterovirus 5' nontranslated region (NTR) contains an internal ribosome entry site (IRES), which facilitates translation initiation of the viral open reading frame in a 5' (m(7)GpppN) cap-independent manner, and cis-acting signals for positive-strand RNA replication. For several enteroviruses, the 5' NTR has been shown to determine the virulence phenotype. We have constructed a chimera consisting of the putative IRES element from the Travis strain of echovirus 12 (ECV12), a wild-type, relatively nonvirulent human enterovirus, exchanged with the homologous region of a full-length infectious clone of coxsackievirus B3 (CBV3). The resulting chimera, known as ECV12(5'NTR)CBV3, replicates similarly to CBV3 in human and simian cell lines yet, unlike CBV3, is completely restricted for growth on two primary murine cell lines at 37 degrees C. By utilizing a reverse-genetics approach, the growth restriction phenotype was localized to the predicted stem-loop II within the IRES of ECV12. In addition, a revertant of ECV12(5'NTR)CBV3 was isolated which possessed three transition mutations and had restored capability for replication in the utilized murine cell lines. Assays for cardiovirulence indicated that the ECV12 IRES is responsible for a noncardiovirulent phenotype in a murine model for acute myocarditis. The results indicate that the 5' NTRs of ECV12 and CBV3 exhibit variable intracellular requirements for function and serve as secondary determinants of tissue or species tropism.     

Effects of the Sabin-like mutations in domain V of the internal ribosome entry segment on translational efficiency of the Coxsackievirus B3

The domain V within the internal ribosome entry segment (IRES) of poliovirus (PV) is expected to be important in its own neurovirulence because it contains an attenuating mutation in each of the Sabin vaccine strains. In this study, we try to find out if the results observed in the case of Sabin vaccine strains of PV can be extrapolated to another virus belonging to the same genus of enteroviruses but with a different tropism. To test this hypothesis, we used the coxsackievirus B3 (CVB3), known to be the most common causal agent of viral myocarditis. The introduction of the three PV Sabin-like mutations in the equivalent positions (nucleotides 484, 485, and 473) to the domain V of the CVB3 IRES results in significant reduced viral titer of the Sabin3-like mutant (Sab3-like) but not on those of Sab1- and Sab2-like mutants. This low titer was correlated with poor translation efficiency in vitro when all mutants were translated in rabbit reticulocyte lysates. However, elucidation by biochemical probing of the secondary structure of the entire domain V of the IRES of Sabin-like mutants reveals no distinct profiles in comparison with the wild-type counterpart. Prediction of secondary structure by MFOLD program indicates a structural perturbation of the stem containing the Sab3-like mutation, suggesting that specific protein-viral RNA interactions are disrupted, preventing efficient viral translation.     

Echovirus 22 is an atypical enterovirus

Although echovirus 22 (EV22) is classified as an enterovirus in the family Picornaviridae, it is atypical of the enterovirus paradigm, typified by the polioviruses and the coxsackie B viruses. cDNA reverse transcribed from coxsackievirus B3 (CVB3) RNA does not hybridize to genomic RNA of EV22, and conversely, cDNA made to EV22 does not hybridize to CVB3 genomic RNA or to molecular clones of CVB3 or poliovirus type 1. EV22 cDNA does not hybridize to viral RNA of encephalomyocarditis virus or to a molecular clone of Theiler's murine encephalomyelitis virus, members of the cardiovirus genus. The genomic RNA of EV22 cannot be detected by the polymerase chain reaction using generic enteroviral primers. EV22 does not shut off host cell protein synthesis, and the RNA of EV22 is efficiently translated in vitro in rabbit reticulocyte lysates. Murine enterovirus-immune T cells recognize and proliferate against EV22 as an antigen in vitro, demonstrating that EV22 shares an epitope(s) common to enteroviruses but not found among other picornaviruses.     

Exploring the repertoire of RNA secondary motifs using graph theory; implications for RNA design

Understanding the structural repertoire of RNA is crucial for RNA genomics research. Yet current methods for finding novel RNAs are limited to small or known RNA families. To expand known RNA structural motifs, we develop a two-dimensional graphical representation approach for describing and estimating the size of RNA’s secondary structural repertoire, including naturally occurring and other possible RNA motifs. We employ tree graphs to describe RNA tree motifs and more general (dual) graphs to describe both RNA tree and pseudoknot motifs. Our estimates of RNA’s structural space are vastly smaller than the nucleotide sequence space, suggesting a new avenue for finding novel RNAs. Specifically our survey shows that known RNA trees and pseudoknots represent only a small subset of all possible motifs, implying that some of the ‘missing’ motifs may represent novel RNAs. To help pinpoint RNA-like motifs, we show that the motifs of existing functional RNAs are clustered in a narrow range of topological characteristics. We also illustrate the applications of our approach to the design of novel RNAs and automated comparison of RNA structures; we report several occurrences of RNA motifs within larger RNAs. Thus, our graph theory approach to RNA structures has implications for RNA genomics, structure analysis and design.     

Nucleotide sequence of the 5'nontranslated and virion polypeptides regions of coxsackievirus B6.

The nucleotide sequence of coxsackievirus B6 (CVB6) has been determined, and the nucleotides encoding the 5' nontranslated region (5' NTR) and virion polypeptides (VP4, 2, 3 and 1) were compared with other serotype CVBs. An Unweighted Pair-Group Method Analysis (UPGMA) of phylogenetic trees indicated that the 5' NTR of CVB6 locates on an independent branch from the other CVBs. The tree based on the amino acid sequences showed that CVB6 has close correlation with CVB4 in the VP4 and VP2 regions, with CVB1 and CVB5 in the VP3 region, and with CVB5 in the VP1 region. Amino acid sequences of variable regions within the VP2, VP3, and VP1 of CVB6 were unique among CVBs. Thus, by comparison of the nucleotide and amino acid sequences of these variable regions, CVB6 can be easily distinguished from other serotypes. In addition, serine, instead of glycine, was found to locate at the amino-terminus of the VP1 region of CVB6, indicating that CVB6 has a unique cleavage site (i.e., glutamine/serine instead of glutamine/glycine) for proteinase 3C of Picornaviridae.     

Clustering RNA structural motifs in ribosomal RNAs using secondary structural alignment

RNA structural motifs are the building blocks of the complex RNA architecture. Identification of non-coding RNA structural motifs is a critical step towards understanding of their structures and functionalities. In this article, we present a clustering approach for de novo RNA structural motif identification. We applied our approach on a data set containing 5S, 16S and 23S rRNAs and rediscovered many known motifs including GNRA tetraloop, kink-turn, C-loop, sarcin-ricin, reverse kink-turn, hook-turn, E-loop and tandem-sheared motifs, with higher accuracy than the state-of-the-art clustering method. We also identified a number of potential novel instances of GNRA tetraloop, kink-turn, sarcin-ricin and tandem-sheared motifs. More importantly, several novel structural motif families have been revealed by our clustering analysis. We identified a highly asymmetric bulge loop motif that resembles the rope sling. We also found an internal loop motif that can significantly increase the twist of the helix. Finally, we discovered a subfamily of hexaloop motif, which has significantly different geometry comparing to the currently known hexaloop motif. Our discoveries presented in this article have largely increased current knowledge of RNA structural motifs.