RNA Structural Motifs Research Articles

De novo discovery of structural motifs in RNA 3D structures through clustering.

As functional components in three-dimensional (3D) conformation of an RNA, the RNA structural motifs provide an easy way to associate the molecular architectures with their biological mechanisms. In the past years, many computational tools have been developed to search motif instances by using the existing knowledge of well-studied families. Recently, with the rapidly increasing number of resolved RNA 3D structures, there is an urgent need to discover novel motifs with the newly presented information. In this work, we classify all the loops in non-redundant RNA 3D structures to detect plausible RNA structural motif families by using a clustering pipeline. Compared with other clustering approaches, our method has two benefits: first, the underlying alignment algorithm is tolerant to the variations in 3D structures. Second, sophisticated downstream analysis has been performed to ensure the clusters are valid and easily applied to further research. The final clustering results contain many interesting new variants of known motif families, such as GNAA tetraloop, kink-turn, sarcin-ricin and T-loop. We have also discovered potential novel functional motifs conserved in ribosomal RNA, sgRNA, SRP RNA, riboswitch and ribozyme.

PATTERNA: transcriptome-wide search for functional RNA elements via structural data signatures

Establishing a link between RNA structure and function remains a great challenge in RNA biology. The emergence of high-throughput structure profiling experiments is revolutionizing our ability to decipher structure, yet principled approaches for extracting information on structural elements directly from these data sets are lacking. We present patteRNA, an unsupervised pattern recognition algorithm that rapidly mines RNA structure motifs from profiling data. We demonstrate that patteRNA detects motifs with an accuracy comparable to commonly used thermodynamic models and highlight its utility in automating data-directed structure modeling from large data sets. patteRNA is versatile and compatible with diverse profiling techniques and experimental conditions.

RNA 3-dimensional structural motifs as a critical constraint of viroid RNA evolution.

Viroids are circular noncoding RNAs infecting plants [1, 2]. During infection, viroids, like RNA viruses, generate swarms of sequence variants called quasispecies [3, 4]. Viroids in Avsunviroidae family replicate in chloroplasts and display the highest mutation rates among all living entities [5]. Viroids in Pospiviroidae family replicate in the nucleus with a relatively lower mutation rate resembling some RNA viruses [6]. Those sequence variants generated during replication are described by the concept of sequence space, which harnesses a geometric representation to illustrate genetic similarities via physical distances. Given the high mutation rate and fast propagation, viroid RNAs have a potentially large sequence space for the evolution of new variants. However, in reality, they use only a small fraction of this space. Constraints of viral sequence space may include genome size, replication fidelity, error thresholds, host or tissue tropism, etc. These factors have been nicely reviewed elsewhere [3, 7, 8] and are not the focus of this Pearl. In addition, RNA secondary structures have been considered, though not adequately, as a constraint factor [8]. Viroids, in contrast to viruses, entirely rely on their RNA structural motifs for function due to their noncoding nature, which offers insights into their capacity to explore regions of sequence space influenced by RNA structures. Here, we describe that 3-dimensional (3D) structural motifs formed by non–Watson-Crick (non-WC) base pairs in viroid RNAs act as a critical constraint for the sequence space of viroid genome evolution. This constraint operates because RNA 3D motifs can play crucial roles by mediating (1) RNA–RNA interactions for the folding of a part or a whole of RNA into a distinct tertiary conformation and (2) RNA–protein interactions. Therefore, mutations in a 3D motif that do not disrupt the structure and function will be retained in the population, whereas mutations that disrupt the 3D structures of motifs, and consequently the function, will be lost.

Insights into Structural and Mechanistic Features of Viral IRES Elements.

Internal ribosome entry site (IRES) elements are cis-acting RNA regions that promote internal initiation of protein synthesis using cap-independent mechanisms. However, distinct types of IRES elements present in the genome of various RNA viruses perform the same function despite lacking conservation of sequence and secondary RNA structure. Likewise, IRES elements differ in host factor requirement to recruit the ribosomal subunits. In spite of this diversity, evolutionarily conserved motifs in each family of RNA viruses preserve sequences impacting on RNA structure and RNA–protein interactions important for IRES activity. Indeed, IRES elements adopting remarkable different structural organizations contain RNA structural motifs that play an essential role in recruiting ribosomes, initiation factors and/or RNA-binding proteins using different mechanisms. Therefore, given that a universal IRES motif remains elusive, it is critical to understand how diverse structural motifs deliver functions relevant for IRES activity. This will be useful for understanding the molecular mechanisms beyond cap-independent translation, as well as the evolutionary history of these regulatory elements. Moreover, it could improve the accuracy to predict IRES-like motifs hidden in genome sequences. This review summarizes recent advances on the diversity and biological relevance of RNA structural motifs for viral IRES elements.

Molecular Dynamics Simulations Reveal an Interplay between SHAPE Reagent Binding and RNA Flexibility.

The function of RNA molecules usually depends on their overall fold and on the presence of specific structural motifs. Chemical probing methods are routinely used in combination with nearest-neighbor models to determine RNA secondary structure. Among the available methods, SHAPE is relevant due to its capability to probe all RNA nucleotides and the possibility to be used in vivo. However, the structural determinants for SHAPE reactivity and its mechanism of reaction are still unclear. Here molecular dynamics simulations and enhanced sampling techniques are used to predict the accessibility of nucleotide analogs and larger RNA structural motifs to SHAPE reagents. We show that local RNA reconformations are crucial in allowing reagents to reach the 2′-OH group of a particular nucleotide and that sugar pucker is a major structural factor influencing SHAPE reactivity.

DotAligner: identification and clustering of RNA structure motifs

The diversity of processed transcripts in eukaryotic genomes poses a challenge for the classification of their biological functions. Sparse sequence conservation in non-coding sequences and the unreliable nature of RNA structure predictions further exacerbate this conundrum. Here, we describe a computational method, DotAligner, for the unsupervised discovery and classification of homologous RNA structure motifs from a set of sequences of interest. Our approach outperforms comparable algorithms at clustering known RNA structure families, both in speed and accuracy. It identifies clusters of known and novel structure motifs from ENCODE immunoprecipitation data for 44 RNA-binding proteins.

Viral internal ribosomal entry sites: four classes for one goal.

To ensure efficient propagation, viruses need to rapidly produce viral proteins after cell entrance. Since viral genomes do not encode any components of the protein biosynthesis machinery, viral proteins must be produced by the host cell. To hi-jack the host cellular translation, viruses use a great variety of distinct strategies. Many single-stranded positive-sensed RNA viruses contain so-called internal ribosome entry sites (IRESs). IRESs are structural RNA motifs that have evolved to specific folds that recruit the host ribosomes on the viral coding sequences in order to synthesize viral proteins. In host canonical translation, recruitment of the translation machinery components is essentially guided by the 5' cap (m7 G) of mRNA. In contrast, IRESs are able to promote efficient ribosome assembly internally and in cap-independent manner. IRESs have been categorized into four classes, based on their length, nucleotide sequence, secondary and tertiary structures, as well as their mode of action. Classes I and II require the assistance of cellular auxiliary factors, the eukaryotic intiation factors (eIF), for efficient ribosome assembly. Class III IRESs require only a subset of eIFs whereas Class IV, which are the more compact, can promote translation without any eIFs. Extensive functional and structural investigations of IRESs over the past decades have allowed a better understanding of their mode of action for viral translation. Because viral translation has a pivotal role in the infectious program, IRESs are therefore attractive targets for therapeutic purposes. WIREs RNA 2018, 9:e1458. doi: 10.1002/wrna.1458 This article is categorized under: Translation > Ribosome Structure/Function Translation > Translation Mechanisms RNA Interactions with Proteins and Other Molecules > RNA-Protein Complexes.

Single-Molecule Fluorescence Reveals Commonalities and Distinctions among Natural and in Vitro-Selected RNA Tertiary Motifs in a Multistep Folding Pathway.

Decades of study of the RNA folding problem have revealed that diverse and complex structured RNAs are built from a common set of recurring structural motifs, leading to the perspective that a generalizable model of RNA folding may be developed from understanding of the folding properties of individual structural motifs. We used single-molecule fluorescence to dissect the kinetic and thermodynamic properties of a set of variants of a common tertiary structural motif, the tetraloop/tetraloop-receptor (TL/TLR). Our results revealed a multistep TL/TLR folding pathway in which preorganization of the ubiquitous AA-platform submotif precedes the formation of the docking transition state and tertiary A-minor hydrogen bond interactions form after the docking transition state. Differences in ion dependences between TL/TLR variants indicated the occurrence of sequence-dependent conformational rearrangements prior to and after the formation of the docking transition state. Nevertheless, varying the junction connecting the TL/TLR produced a common kinetic and ionic effect for all variants, suggesting that the global conformational search and compaction electrostatics are energetically independent from the formation of the tertiary motif contacts. We also found that in vitro-selected variants, despite their similar stability at high Mg2+ concentrations, are considerably less stable than natural variants under near-physiological ionic conditions, and the occurrence of the TL/TLR sequence variants in Nature correlates with their thermodynamic stability in isolation. Overall, our findings are consistent with modular but complex energetic properties of RNA structural motifs and will aid in the eventual quantitative description of RNA folding from its secondary and tertiary structural elements.

Overexpression of Ebola virus envelope GP1 protein

Ebola virus uses its envelope GP1 and GP2 for viral attachment and entry into host cells. Due to technical difficulty expressing full-length envelope, many structural and functional studies of Ebola envelope protein have been carried out primarily using GP1 lacking its mucin-like domain. As a result, the viral invasion mechanisms involving the mucin-like domain are not fully understood. To elucidate the role of the mucin-like domain of GP1 in Ebola-host attachment and infection and to facilitate vaccine development, we constructed a GP1 expression vector containing the entire attachment region (1–496). Cysteine 53 of GP1, which forms a disulfide bond with GP2, was mutated to serine to avoid potential disulfide bond mispairing. Stable expression clones using codon optimized open reading frame were developed in human 293-H cells with yields reaching ∼25 mg of GP1 protein per liter of spent medium. Purified GP1 was functional and bound to Ebola attachment receptors, DC-SIGN and DC-SIGNR. The over-expression and easy purification characteristic of this system has implications in Ebola research and vaccine development. To further understand the differential expression yields between the codon optimized and native GP1, we analyzed the presence of RNA structural motifs in the first 100 nucleotides of translational initiation AUG site. RNA structural prediction showed the codon optimization removed two potential RNA pseudoknot structures. This methodology is also applicable to the expression of other difficult virus envelope proteins.

Recurring RNA structural motifs underlie the mechanics of L1 stalk movement

The L1 stalk of the large ribosomal subunit undergoes large-scale movements coupled to the translocation of deacylated tRNA during protein synthesis. We use quantitative comparative structural analysis to localize the origins of L1 stalk movement and to understand its dynamic interactions with tRNA and other structural elements of the ribosome. Besides its stacking interactions with the tRNA elbow, stalk movement is directly linked to intersubunit rotation, rotation of the 30S head domain and contact of the acceptor arm of deacylated tRNA with helix 68 of 23S rRNA. Movement originates from pivoting at stacked non-canonical base pairs in a Family A three-way junction and bending in an internal G-U-rich zone. Use of these same motifs as hinge points to enable such dynamic events as rotation of the 30S subunit head domain and in flexing of the anticodon arm of tRNA suggests that they represent general strategies for movement of functional RNAs.

Predicting RNA Structure with Vfold.

In order to carry out biological functions, RNA molecules must fold into specific three-dimensional (3D) structures. Current experimental methods to determine RNA 3D structures are expensive and time consuming. With the recent advances in computational biology, RNA structure prediction is becoming increasingly reliable. This chapter describes a recently developed RNA structure prediction software, Vfold, a virtual bond-based RNA folding model. The main features of Vfold are the physics-based loop free energy calculations for various RNA structure motifs and a template-based assembly method for RNA 3D structure prediction. For illustration, we use the yybP-ykoY Orphan riboswitch as an example to show the implementation of the Vfold model in RNA structure prediction from the sequence.

Exploring Consensus RNA Substructural Patterns Using Subgraph Mining.

Frequently recurring RNA structural motifs play important roles in RNA folding process and interaction with other molecules. Traditional index-based and shape-based schemas are useful in modeling RNA secondary structures but ignore the structural discrepancy of individual RNA family member. Further, the in-depth analysis of underlying substructure pattern is insufficient due to varied and unnormalized substructure data. This prevents us from understanding RNAs functions and their inherent synergistic regulation networks. This article thus proposes a novel labeled graph-based algorithm RnaGraph to uncover frequently RNA substructure patterns. Attribute data and graph data are combined to characterize diverse substructures and their correlations, respectively. Further, a top-k graph pattern mining algorithm is developed to extract interesting substructure motifs by integrating frequency and similarity. The experimental results show that our methods assist in not only modelling complex RNA secondary structures but also identifying hidden but interesting RNA substructure patterns.

Small Molecule-Based Pattern Recognition To Classify RNA Structure.

Three-dimensional RNA structures are notoriously difficult to determine, and the link between secondary structure and RNA conformation is only beginning to be understood. These challenges have hindered the identification of guiding principles for small molecule:RNA recognition. We herein demonstrate that the strong and differential binding ability of aminoglycosides to RNA structures can be used to classify five canonical RNA secondary structure motifs through principal component analysis (PCA). In these analyses, the aminoglycosides act as receptors, while RNA structures labeled with a benzofuranyluridine fluorophore act as analytes. Complete (100%) predictive ability for this RNA training set was achieved by incorporating two exhaustively guanidinylated aminoglycosides into the receptor library. The PCA was then externally validated using biologically relevant RNA constructs. In bulge-stem-loop constructs of HIV-1 transactivation response element (TAR) RNA, we achieved nucleotide-specific classification of two independent secondary structure motifs. Furthermore, examination of cheminformatic parameters and PCA loading factors revealed trends in aminoglycoside:RNA recognition, including the importance of shape-based discrimination, and suggested the potential for size and sequence discrimination within RNA structural motifs. These studies present a new approach to classifying RNA structure and provide direct evidence that RNA topology, in addition to sequence, is critical for the molecular recognition of RNA.

Structural and Functional Annotation of Long Noncoding RNAs.

Protein-coding RNAs represent only a small fraction of the transcriptional output in higher eukaryotes. The remaining RNA species encompass a broad range of molecular functions and regulatory roles, a consequence of the structural polyvalence of RNA polymers. Albeit several classes of small noncoding RNAs are relatively well characterized, the accessibility of affordable high-throughput sequencing is generating a wealth of novel, unannotated transcripts, especially long noncoding RNAs (lncRNAs) that are derived from genomic regions that are antisense, intronic, intergenic, and overlapping protein-coding loci. Parsing and characterizing the functions of noncoding RNAs-lncRNAs in particular-is one of the great challenges of modern genome biology. Here we discuss concepts and computational methods for the identification of structural domains in lncRNAs from genomic and transcriptomic data. In the first part, we briefly review how to identify RNA structural motifs in individual lncRNAs. In the second part, we describe how to leverage the evolutionary dynamics of structured RNAs in a computationally efficient screen to detect putative functional lncRNA motifs using comparative genomics.

A short loop at the heart of a long RNA

A specific structural motif in a long noncoding RNA mediates the interaction with a transcription factor to drive cardiac cell differentiation.

SHAPE reveals transcript-wide interactions, complex structural domains, and protein interactions across the Xist lncRNA in living cells

The 18-kb Xist long noncoding RNA (lncRNA) is essential for X-chromosome inactivation during female eutherian mammalian development. Global structural architecture, cell-induced conformational changes, and protein-RNA interactions within Xist are poorly understood. We used selective 2'-hydroxyl acylation analyzed by primer extension and mutational profiling (SHAPE-MaP) to examine these features of Xist at single-nucleotide resolution both in living cells and ex vivo. The Xist RNA forms complex well-defined secondary structure domains and the cellular environment strongly modulates the RNA structure, via motifs spanning one-half of all Xist nucleotides. The Xist RNA structure modulates protein interactions in cells via multiple mechanisms. For example, repeat-containing elements adopt accessible and dynamic structures that function as landing pads for protein cofactors. Structured RNA motifs create interaction domains for specific proteins and also sequester other motifs, such that only a subset of potential binding sites forms stable interactions. This work creates a broad quantitative framework for understanding structure-function interrelationships for Xist and other lncRNAs in cells.

Quadruplex-Flanking Stem Structures Modulate the Stability and Metal Ion Preferences of RNA Mimics of GFP.

The spinach family of RNA aptamers are RNA mimics of green fluorescent protein (GFP) that have previously been designed to address the challenges of imaging RNA inside living cells. However, relatively low levels of free intracellular magnesium limited the practical use of these aptamers. Recent cell-based selections identified the broccoli RNA aptamer, which requires less magnesium for fluorescence, but the basis for magnesium preference remained unclear. Here, we find that the broccoli RNA structure is very similar to that of baby spinach, a truncated version of the spinach aptamer. Differences in stability and metal ion preferences between these two aptamers, and among broccoli mutants, are primarily associated with the sequence and structure of predicted quadruplex-flanking stem structures. Mutation of purine-purine pairs in broccoli at the terminal stem-quadruplex transition caused reversion of broccoli to a higher magnesium dependence. Unique duplex-to-quadruplex transitions in GFP-mimic RNAs likely explain their sensitivity to magnesium for stability and fluorescence. Thus, optimizations designed to improve aptamers should take into consideration the role of metal ions in stabilizing the transitions and interactions between independently folding RNA structural motifs.

Ligand-Induced Stabilization of a Duplex-like Architecture Is Crucial for the Switching Mechanism of the SAM-III Riboswitch.

Riboswitches are structured RNA motifs that control gene expression by sensing the concentrations of specific metabolites and make up a promising new class of antibiotic targets. S-Adenosylmethionine (SAM)-III riboswitch, mainly found in lactic acid bacteria, is involved in regulating methionine and SAM biosynthetic pathways. SAM-III riboswitch regulates the gene expression by switching the translation process on and off with respect to the absence and presence of the SAM ligand, respectively. In this study, an attempt is made to understand the key conformational transitions involved in ligand binding using atomistic molecular dynamics (MD) simulations performed in an explicit solvent environment. G26 is found to recognize the SAM ligand by forming hydrogen bonds, whereas the absence of the ligand leads to opening of the binding pocket. Consistent with experimental results, the absence of the SAM ligand weakens the base pairing interactions between the nucleobases that are part of the Shine-Dalgarno (SD) and anti-Shine-Dalgarno (aSD) sequences, which in turn facilitates recognition of the SD sequence by ribosomes. Detailed analysis reveals that a duplex-like structure formed by nucleotides from different parts of the RNA and the adenine base of the ligand is crucial for the stability of the completely folded state in the presence of the ligand. Previous experimental studies have shown that the SAM-III riboswitch exists in equilibrium between the unfolded and partially folded states in the absence of the ligand, which completely folds upon binding of the ligand. Comparison of the results presented here to the available experimental data indicates the structures obtained using the MD simulations resemble the partially folded state. Thus, this study provides a detailed understanding of the fully and partially folded structures of the SAM-III riboswitch in the presence and absence of the ligand, respectively. This study hypothesizes a dual role for the SAM ligand, which facilitates conformational switching between partially and fully folded states by forming a stable duplex-like structure and strengthening the interactions between SD and aSD nucleotides.

Abstract IA25: Regulatory control of lincRNA function through formation of complex RNA structural motifs

Abstract While the regulatory role of lincRNAs is becoming increasingly appreciated, the functional role of RNA structures within these molecules is not well understood. We have embarked on a study to understand the folding of lincRNA molecules, and to identify functional structural modules within large RNA transcripts. Using approaches to promote the native folding of large RNA molecules, we have been purifying lincRNAs, mapping their structures with chemical probes, and then using this information to construct structural maps of lincRNA architecture. This work is complemented by phylogenetic and genetic approaches for exploring the functional relevance of RNA structures identified in our study. In this way, we have mapped the structure of lincRNAs such as HOTAIR and lincRNA p21, identifying novel structural motifs, enabling evolutionary comparisons and facilitating targeted mutational studies of lincRNA function. Citation Format: Isabel Chillon, Srinivas Somarowthu, Fei Liu, Anna Marie Pyle. Regulatory control of lincRNA function through formation of complex RNA structural motifs. [abstract]. In: Proceedings of the AACR Special Conference on Noncoding RNAs and Cancer: Mechanisms to Medicines ; 2015 Dec 4-7; Boston, MA. Philadelphia (PA): AACR; Cancer Res 2016;76(6 Suppl):Abstract nr IA25.

The Coding Region of the HCV Genome Contains a Network of Regulatory RNA Structures

RNA is a versatile macromolecule that accommodates functional information in primary sequence and secondary and tertiary structure. We use a combination of chemical probing, RNA structure modeling, comparative sequence analysis, and functional assays to examine the role of RNA structure in the hepatitis C virus (HCV) genome. We describe a set of conserved but functionally diverse structural RNA motifs that occur in multiple coding regions of the HCV genome, and we demonstrate that conformational changes in these motifs influence specific stages in the virus' life cycle. Our study shows that these types of structures can pervade a genome, where they play specific mechanistic and regulatory roles, constituting a "code within the code" for controlling biological processes.