RNA structural motifs: building blocks of a modular biomolecule

A structural map of oncomiR-1 at single-nucleotide resolution

The miR-17-92a cluster, also known as 'oncomiR-1', is an RNA transcript that plays a pivotal regulatory role in cellular processes, including the cell cycle, proliferation and apoptosis. Its dysregulation underlies the development of several cancers. Oncomir-1 comprises six constituent miRNAs, each processed with different efficiencies as a function of both developmental time and tissue type. The structural mechanisms that regulate such differential processing are unknown, and this has impeded our understanding of the dysregulation of oncomiR-1 in pathophysiology. By probing the sensitivity of each nucleotide in oncomiR-1 to reactive small molecules, we present a secondary structural map of this RNA at single-nucleotide resolution. The secondary structure and solvent accessible regions of oncomiR-1 reveal that most of its primary microRNA domains are suboptimal substrates for Drosha-DGCR8, and therefore resistant to microprocessing. The structure indicates that the binding of trans-acting factors is required to remodel the tertiary organization and unmask cryptic primary microRNA domains to facilitate their processing into pre-microRNAs.

Dawn of the in vivo RNA structurome and interactome

RNA is one of the most fascinating biomolecules in living systems given its structural versatility to fold into elaborate architectures for important biological functions such as gene regulation, catalysis, and information storage. Knowledge of RNA structures and interactions can provide deep insights into their functional roles in vivo For decades, RNA structural studies have been conducted on a transcript-by-transcript basis. The advent of next-generation sequencing (NGS) has enabled the development of transcriptome-wide structural probing methods to profile the global landscape of RNA structures and interactions, also known as the RNA structurome and interactome, which transformed our understanding of the RNA structure-function relationship on a transcriptomic scale. In this review, molecular tools and NGS methods used for RNA structure probing are presented, novel insights uncovered by RNA structurome and interactome studies are highlighted, and perspectives on current challenges and potential future directions are discussed. A more complete understanding of the RNA structures and interactions in vivo will help illuminate the novel roles of RNA in gene regulation, development, and diseases.

The super-n-motifs model: a novel alignment-free approach for representing and comparing RNA secondary structures

Motivation

Comparing ribonucleic acid (RNA) secondary structures of arbitrary size uncovers structural patterns that can provide a better understanding of RNA functions. However, performing fast and accurate secondary structure comparisons is challenging when we take into account the RNA configuration (i.e. linear or circular), the presence of pseudoknot and G-quadruplex (G4) motifs and the increasing number of secondary structures generated by high-throughput probing techniques. To address this challenge, we propose the super-n-motifs model based on a latent analysis of enhanced motifs comprising not only basic motifs but also adjacency relations. The super-n-motifs model computes a vector representation of secondary structures as linear combinations of these motifs.

Results

We demonstrate the accuracy of our model for comparison of secondary structures from linear and circular RNA while also considering pseudoknot and G4 motifs. We show that the super-n-motifs representation effectively captures the most important structural features of secondary structures, as compared to other representations such as ordered tree, arc-annotated and string representations. Finally, we demonstrate the time efficiency of our model, which is alignment free and capable of performing large-scale comparisons of 10 000 secondary structures with an efficiency up to 4 orders of magnitude faster than existing approaches.

Availability and Implementation

The super-n-motifs model was implemented in C ++. Source code and Linux binary are freely available at http://jpsglouzon.github.io/supernmotifs/ .

Contact

Shengrui.Wang@Usherbrooke.ca.

Supplementary information

Supplementary data are available at Bioinformatics o nline.

Revisiting GNRA and UNCG folds: U-turns versus Z-turns in RNA hairpin loops

When thinking about RNA three-dimensional structures, coming across GNRA and UNCG tetraloops is perceived as a boon since their folds have been extensively described. Nevertheless, analyzing loop conformations within RNA and RNP structures led us to uncover several instances of GNRA and UNCG loops that do not fold as expected. We noticed that when a GNRA does not assume its "natural" fold, it adopts the one we typically associate with a UNCG sequence. The same folding interconversion may occur for loops with UNCG sequences, for instance within tRNA anticodon loops. Hence, we show that some structured tetranucleotide sequences starting with G or U can adopt either of these folds. The underlying structural basis that defines these two fold types is the mutually exclusive stacking of a backbone oxygen on either the first (in GNRA) or the last nucleobase (in UNCG), generating an oxygen-π contact. We thereby propose to refrain from using sequences to distinguish between loop conformations. Instead, we suggest using descriptors such as U-turn (for "GNRA-type" folds) and a newly described Z-turn (for "UNCG-type" folds). Because tetraloops adopt for the largest part only two (inter)convertible turns, we are better able to interpret from a structural perspective loop interchangeability occurring in ribosomes and viral RNA. In this respect, we propose a general view on the inclination for a given sequence to adopt (or not) a specific fold. We also suggest how long-noncoding RNAs may adopt discrete but transient structures, which are therefore hard to predict.

Aptamer Technology: Adjunct Therapy for Malaria

Malaria is a life-threatening parasitic infection occurring in the endemic areas, primarily in children under the age of five, pregnant women, and patients with human immunodeficiency virus and acquired immunodeficiency syndrome (HIV)/(AIDS) as well as non-immune individuals. The cytoadherence of infected erythrocytes (IEs) to the host endothelial surface receptor is a known factor that contributes to the increased prevalence of severe malaria cases due to the accumulation of IEs, mainly in the brain and other vital organs. Therefore, further study is needed to discover a new potential anti-adhesive drug to treat severe malaria thus reducing its mortality rate. In this review, we discuss how the aptamer technology could be applied in the development of a new adjunct therapy for current malaria treatment.

Stacking interactions involving non-Watson–Crick basepairs: dispersion corrected density functional theory studies

Stacking interactions between a non Watson–Crick G:A S:HT basepair and C:G basepair is predicted in terms of roll, twist and slide basepair step parameters using DFT-D augmented with coarse-grain energy penalty for sugar–phosphate backbone.

Structural Basis for Ligand Binding to the Guanidine-I Riboswitch

The guanidine-I riboswitch is a conserved RNA element with approximately 2,000 known examples across four phyla of bacteria. It exists upstream of nitrogen metabolism and multidrug resistance transporter genes and alters expression through the specific recognition of a free guanidinium cation. Here we report the structure of a guanidine riboswitch aptamer from Sulfobacillus acidophilus at 2.7 Å resolution. Helices P1, P1a, P1b, and P2 form a coaxial stack that acts as a scaffold for ligand binding. A previously unidentified P3 helix docks into P1a to form the guanidinium binding pocket, which is completely enclosed. Every functional group of the ligand is recognized through hydrogen bonding to guanine bases and phosphate oxygens. Guanidinium binding is further stabilized through cation-π interactions with guanine bases. This allows the riboswitch to recognize guanidinium while excluding other bacterial metabolites with a guanidino group, including the amino acid arginine.

Crystal structure of a DNA aptamer bound to PvLDH elucidates novel single-stranded DNA structural elements for folding and recognition

Structural elements are key elements for understanding single-stranded nucleic acid folding. Although various RNA structural elements have been documented, structural elements of single-stranded DNA (ssDNA) have rarely been reported. Herein, we determined a crystal structure of PvLDH in complex with a DNA aptamer called pL1. This aptamer folds into a hairpin-bulge contact by adopting three novel structural elements, viz, DNA T-loop-like motif, base-phosphate zipper, and DNA G·G metal ion zipper. Moreover, the pL1:PvLDH complex shows unique properties compared with other protein:nucleic acid complexes. Generally, extensive intermolecular hydrogen bonds occur between unpaired nucleotides and proteins for specific recognitions. Although most protein-interacting nucleotides of pL1 are unpaired nucleotides, pL1 recognizes PvLDH by predominant shape complementarity with many bridging water molecules owing to the combination of three novel structural elements making protein-binding unpaired nucleotides stable. Moreover, the additional set of Plasmodium LDH residues which were shown to form extensive hydrogen bonds with unpaired nucleotides of 2008s does not participate in the recognition of pL1. Superimposition of the pL1:PvLDH complex with hLDH reveals steric clashes between pL1 and hLDH in contrast with no steric clashes between 2008s and hLDH. Therefore, specific protein recognition mode of pL1 is totally different from that of 2008s.

Secondary structure model of the naked segment 7 influenza A virus genomic RNA

The influenza A virus (IAV) genome comprises eight negative-sense viral (v)RNA segments. The seventh segment of the genome encodes two essential viral proteins and is specifically packaged alongside the other seven vRNAs. To gain insights into the possible roles of RNA structure both within and without virions, a secondary structure model of a naked (protein-free) segment 7 vRNA (vRNA7) has been determined using chemical mapping and thermodynamic energy minimization. The proposed structure model was validated using microarray mapping, RNase H cleavage and comparative sequence analysis. Additionally, the detailed structures of three vRNA7 fragment constructs - comprising independently folded subdomains - were determined. Much of the proposed vRNA7 structure is preserved between IAV strains, suggesting their importance in the influenza replication cycle. Possible structure rearrangements, which allow or preclude long-range RNA interactions, are also proposed.

Ring Catalog: A resource for designing self-assembling RNA nanostructures

Designing self-assembling RNA ring structures based on known 3D structural elements connected via linker helices is a challenging task due to the immense number of motif combinations, many of which do not lead to ring-closure. We describe an in silico solution to this design problem by combinatorial assembly of RNA 3-way junctions, bulges, and kissing loops, and tabulating the cases that lead to ring formation. The solutions found are made available in the form of a web-accessible Ring Catalog. As an example of a potential use of this resource, we chose a predicted RNA square structure consisting of five RNA strands and demonstrate experimentally that the self-assembly of those five strands leads to the formation of a square-like complex. This is a demonstration of a novel "design by catalog" approach to RNA nano-structure generation. The URL https://rnajunction.ncifcrf.gov/ringdb can be used to access the resource.

Improving NMR Structures of RNA

Here, we show that modern solution nuclear magnetic resonance (NMR) structures of RNA exhibit more steric clashes and conformational ambiguities than their crystallographic X-ray counterparts. To tackle these issues, we developed RNA-ff1, a new force field for structure calculation with Xplor-NIH. Using seven published NMR datasets, RNA-ff1 improves covalent geometry and MolProbity validation criteria for clashes and backbone conformation in most cases, relative to both the previous Xplor-NIH force field and the original structures associated with the experimental data. In addition, with smaller base-pair step rises in helical stems, RNA-ff1 structures enjoy more favorable base stacking. Finally, structural accuracy improves in the majority of cases, as supported by complete residual dipolar coupling cross-validation. Thus, the reported advances show great promise in bridging the quality gap that separates NMR and X-ray structures of RNA.

Embracing proteins: structural themes in aptamer-protein complexes

Understanding the structural rules that govern specific, high-affinity binding characteristic of aptamer-protein interactions is important in view of the increasing use of aptamers across many applications. From the modest number of 16 aptamer-protein structures currently available, trends are emerging. The flexible phosphodiester backbone allows folding into precise three-dimensional structures using known nucleic acid motifs as scaffolds that orient specific functional groups for target recognition. Still, completely novel motifs essential for structure and function are found in modified aptamers with diversity-enhancing side chains. Aptamers and antibodies, two classes of macromolecules used as affinity reagents with entirely different backbones and composition, recognize protein epitopes of similar size and with comparably high shape complementarity.

RNA 3D Modules in Genome-Wide Predictions of RNA 2D Structure

Recent experimental and computational progress has revealed a large potential for RNA structure in the genome. This has been driven by computational strategies that exploit multiple genomes of related organisms to identify common sequences and secondary structures. However, these computational approaches have two main challenges: they are computationally expensive and they have a relatively high false discovery rate (FDR). Simultaneously, RNA 3D structure analysis has revealed modules composed of non-canonical base pairs which occur in non-homologous positions, apparently by independent evolution. These modules can, for example, occur inside structural elements which in RNA 2D predictions appear as internal loops. Hence one question is if the use of such RNA 3D information can improve the prediction accuracy of RNA secondary structure at a genome-wide level. Here, we use RNAz in combination with 3D module prediction tools and apply them on a 13-way vertebrate sequence-based alignment. We find that RNA 3D modules predicted by metaRNAmodules and JAR3D are significantly enriched in the screened windows compared to their shuffled counterparts. The initially estimated FDR of 47.0% is lowered to below 25% when certain 3D module predictions are present in the window of the 2D prediction. We discuss the implications and prospects for further development of computational strategies for detection of RNA 2D structure in genomic sequence.

New in silico approach to assessing RNA secondary structures with non-canonical base pairs

Background

The function of RNA is strongly dependent on its structure, so an appropriate recognition of this structure, on every level of organization, is of great importance. One particular concern is the assessment of base-base interactions, described as the secondary structure, the knowledge of which greatly facilitates an interpretation of RNA function and allows for structure analysis on the tertiary level. The RNA secondary structure can be predicted from a sequence using in silico methods often adjusted with experimental data, or assessed from 3D structure atom coordinates. Computational approaches typically consider only canonical, Watson-Crick and wobble base pairs. Handling of non-canonical interactions, important for a full description of RNA structure, is still very difficult.

Results

We introduce our novel approach to assessing an extended RNA secondary structure, which characterizes both canonical and non-canonical base pairs, along with their type classification. It is based on predicting the RNA 3D structure from a user-provided sequence or a secondary structure that only describes canonical base pairs, and then deriving the extended secondary structure from atom coordinates. In our example implementation, this was achieved by integrating the functionality of two fully automated, high fidelity methods in a computational pipeline: RNAComposer for the 3D RNA structure prediction and RNApdbee for base-pair annotation.

Conclusions

The presented methodology ties together existing applications for RNA 3D structure prediction and base-pair annotation. The example performance, applying RNAComposer and RNApdbee, reveals better accuracy in non-canonical base pair assessment than the compared methods that directly predict RNA secondary structure.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0718-6) contains supplementary material, which is available to authorized users.

RAG-3D: a search tool for RNA 3D substructures

To address many challenges in RNA structure/function prediction, the characterization of RNA's modular architectural units is required. Using the RNA-As-Graphs (RAG) database, we have previously explored the existence of secondary structure (2D) submotifs within larger RNA structures. Here we present RAG-3D—a dataset of RNA tertiary (3D) structures and substructures plus a web-based search tool—designed to exploit graph representations of RNAs for the goal of searching for similar 3D structural fragments. The objects in RAG-3D consist of 3D structures translated into 3D graphs, cataloged based on the connectivity between their secondary structure elements. Each graph is additionally described in terms of its subgraph building blocks. The RAG-3D search tool then compares a query RNA 3D structure to those in the database to obtain structurally similar structures and substructures. This comparison reveals conserved 3D RNA features and thus may suggest functional connections. Though RNA search programs based on similarity in sequence, 2D, and/or 3D structural elements are available, our graph-based search tool may be advantageous for illuminating similarities that are not obvious; using motifs rather than sequence space also reduces search times considerably. Ultimately, such substructuring could be useful for RNA 3D structure prediction, structure/function inference and inverse folding.

MultiSETTER: web server for multiple RNA structure comparison

Background

Understanding the architecture and function of RNA molecules requires methods for comparing and analyzing their tertiary and quaternary structures. While structural superposition of short RNAs is achievable in a reasonable time, large structures represent much bigger challenge. Therefore, we have developed a fast and accurate algorithm for RNA pairwise structure superposition called SETTER and implemented it in the SETTER web server. However, though biological relationships can be inferred by a pairwise structure alignment, key features preserved by evolution can be identified only from a multiple structure alignment. Thus, we extended the SETTER algorithm to the alignment of multiple RNA structures and developed the MultiSETTER algorithm.

Results

In this paper, we present the updated version of the SETTER web server that implements a user friendly interface to the MultiSETTER algorithm. The server accepts RNA structures either as the list of PDB IDs or as user-defined PDB files. After the superposition is computed, structures are visualized in 3D and several reports and statistics are generated.

Conclusion

To the best of our knowledge, the MultiSETTER web server is the first publicly available tool for a multiple RNA structure alignment. The MultiSETTER server offers the visual inspection of an alignment in 3D space which may reveal structural and functional relationships not captured by other multiple alignment methods based either on a sequence or on secondary structure motifs.

Multiple 3D RNA Structure Superposition Using Neighbor Joining

Recent advances in RNA research and the steady growth of available RNA structures call for bioinformatics methods for handling and analyzing RNA structural data. Recently, we introduced SETTER-a fast and accurate method for RNA pairwise structure alignment. In this paper, we describe MultiSETTER, SETTER extension for multiple RNA structure alignment. MultiSETTER combines SETTER's decomposition of RNA structures into non-overlapping structural subunits with the multiple sequence alignment algorithm ClustalW adapted for the structure alignment. The accuracy of MultiSETTER was assessed by the automatic classification of RNA structures and its comparison to SCOR annotations. In addition, MultiSETTER classification was also compared to multiple sequence alignment-based and secondary structure alignment-based classifications provided by LocARNA and RNADistance tools, respectively. MultiSETTER precompiled Windows libraries, as well as the C++ source code, are freely available from http://siret.cz/multisetter.

Structural fidelity and NMR relaxation analysis in a prototype RNA hairpin

RNA hairpins are widespread and very stable motifs that contribute decisively to RNA folding and biological function. The GTP1G2C3A4C5U6U7C8G9G10U11G12C13C14 construct (with a central UUCG tetraloop) has been extensively studied by solution NMR, and offers and excellent opportunity to evaluate the structure and dynamical description afforded by molecular dynamics (MD) simulations. Here, we compare average structural parameters and NMR relaxation rates estimated from a series of multiple independent explicit solvent MD simulations using the two most recent RNA AMBER force fields (ff99 and ff10). Predicted overall tumbling times are ∼20% faster than those inferred from analysis of NMR data and follow the same trend when temperature and ionic strength is varied. The Watson-Crick stem and the "canonical" UUCG loop structure are maintained in most simulations including the characteristic syn conformation along the glycosidic bond of G9, although some key hydrogen bonds in the loop are partially disrupted. Our analysis pinpoints G9-G10 backbone conformations as a locus of discrepancies between experiment and simulation. In general the results for the more recent force-field parameters (ff10) are closer to experiment than those for the older ones (ff99). This work provides a comprehensive and detailed comparison of state of the art MD simulations against a wide variety of solution NMR measurements.

RNAMotifScanX: a graph alignment approach for RNA structural motif identification

RNA structural motifs are recurrent three-dimensional (3D) components found in the RNA architecture. These RNA structural motifs play important structural or functional roles and usually exhibit highly conserved 3D geometries and base-interaction patterns. Analysis of the RNA 3D structures and elucidation of their molecular functions heavily rely on efficient and accurate identification of these motifs. However, efficient RNA structural motif search tools are lacking due to the high complexity of these motifs. In this work, we present RNAMotifScanX, a motif search tool based on a base-interaction graph alignment algorithm. This novel algorithm enables automatic identification of both partially and fully matched motif instances. RNAMotifScanX considers noncanonical base-pairing interactions, base-stacking interactions, and sequence conservation of the motifs, which leads to significantly improved sensitivity and specificity as compared with other state-of-the-art search tools. RNAMotifScanX also adopts a carefully designed branch-and-bound technique, which enables ultra-fast search of large kink-turn motifs against a 23S rRNA. The software package RNAMotifScanX is implemented using GNU C++, and is freely available from http://genome.ucf.edu/RNAMotifScanX.

Interconversion between parallel and antiparallel conformations of a 4H RNA junction in domain 3 of foot-and-mouth disease virus IRES captured by dynamics simulations

RNA junctions are common secondary structural elements present in a wide range of RNA species. They play crucial roles in directing the overall folding of RNA molecules as well as in a variety of biological functions. In particular, there has been great interest in the dynamics of RNA junctions, including conformational pathways of fully base-paired 4-way (4H) RNA junctions. In such constructs, all nucleotides participate in one of the four double-stranded stem regions, with no connecting loops. Dynamical aspects of these 4H RNAs are interesting because frequent interchanges between parallel and antiparallel conformations are thought to occur without binding of other factors. Gel electrophoresis and single-molecule fluorescence resonance energy transfer experiments have suggested two possible pathways: one involves a helical rearrangement via disruption of coaxial stacking, and the other occurs by a rotation between the helical axes of coaxially stacked conformers. Employing molecular dynamics simulations, we explore this conformational variability in a 4H junction derived from domain 3 of the foot-and-mouth disease virus internal ribosome entry site (IRES); this junction contains highly conserved motifs for RNA-RNA and RNA-protein interactions, important for IRES activity. Our simulations capture transitions of the 4H junction between parallel and antiparallel conformations. The interconversion is virtually barrier-free and occurs via a rotation between the axes of coaxially stacked helices with a transient perpendicular intermediate. We characterize this transition, with various interhelical orientations, by pseudodihedral angle and interhelical distance measures. The high flexibility of the junction, as also demonstrated experimentally, is suitable for IRES activity. Because foot-and-mouth disease virus IRES structure depends on long-range interactions involving domain 3, the perpendicular intermediate, which maintains coaxial stacking of helices and thereby consensus primary and secondary structure information, may be beneficial for guiding the overall organization of the RNA system in domain 3.

Evolutionary Origin and Conserved Structural Building Blocks of Riboswitches and Ribosomal RNAs: Riboswitches as Probable Target Sites for Aminoglycosides Interaction

PURPOSE

Riboswitches, as noncoding RNA sequences, control gene expression through direct ligand binding. Sporadic reports on the structural relation of riboswitches with ribosomal RNAs (rRNA), raises an interest in possible similarity between riboswitches and rRNAs evolutionary origins. Since aminoglycoside antibiotics affect microbial cells through binding to functional sites of the bacterial rRNA, finding any conformational and functional relation between riboswitches/rRNAs is utmost important in both of medicinal and basic research.

METHODS

Analysis of the riboswitches structures were carried out using bioinformatics and computational tools. The possible functional similarity of riboswitches with rRNAs was evaluated based on the affinity of paromomycin antibiotic (targeting "A site" of 16S rRNA) to riboswitches via docking method.

RESULTS

There was high structural similarity between riboswitches and rRNAs, but not any particular sequence based similarity between them was found. The building blocks including "hairpin loop containing UUU", "peptidyl transferase center conserved hairpin A loop"," helix 45" and "S2 (G8) hairpin" as high identical rRNA motifs were detected in all kinds of riboswitches. Surprisingly, binding energies of paromomycin with different riboswitches are considerably better than the binding energy of paromomycin with "16S rRNA A site". Therefore the high affinity of paromomycin to bind riboswitches in comparison with rRNA "A site" suggests a new insight about riboswitches as possible targets for aminoglycoside antibiotics.

CONCLUSION

These findings are considered as a possible supporting evidence for evolutionary origin of riboswitches/rRNAs and also their role in the exertion of antibiotics effects to design new drugs based on the concomitant effects via rRNA/riboswitches.

Structural determinants for geometry and information decoding of tRNA by T box leader RNA

T box riboswitches are cis-acting RNA elements that bind to tRNA and sense its aminoacylation state to influence gene expression. Here, we present the 3.2 Å resolution X-ray crystal structures of the T box Stem I-tRNA complex and tRNA, in isolation. T box Stem I forms an arched conformation with extensive intermolecular contacts to two key points of tRNA, the anticodon and D/T-loops. Free and complexed tRNA exist in significantly different conformations, with the contacts stabilizing flexible D/T-loops and a rearrangement of the D-loop. Using a designed T box RNA/tRNA system, we demonstrate that the T box riboswitch monitors the length and orientation of two essential contacts. Length or orientation mismatches engineered into the T box riboswitch and tRNA disrupt the complex, whereas simultaneous insertion of full helical turns realigns the interfaces and restores interaction between artificially elongated T box riboswitch and tRNA molecules.

Automated identification of RNA 3D modules with discriminative power in RNA structural alignments

Recent progress in predicting RNA structure is moving towards filling the ‘gap’ in 2D RNA structure prediction where, for example, predicted internal loops often form non-canonical base pairs. This is increasingly recognized with the steady increase of known RNA 3D modules. There is a general interest in matching structural modules known from one molecule to other molecules for which the 3D structure is not known yet. We have created a pipeline, metaRNAmodules, which completely automates extracting putative modules from the FR3D database and mapping of such modules to Rfam alignments to obtain comparative evidence. Subsequently, the modules, initially represented by a graph, are turned into models for the RMDetect program, which allows to test their discriminative power using real and randomized Rfam alignments. An initial extraction of 22 495 3D modules in all PDB files results in 977 internal loop and 17 hairpin modules with clear discriminatory power. Many of these modules describe only minor variants of each other. Indeed, mapping of the modules onto Rfam families results in 35 unique locations in 11 different families. The metaRNAmodules pipeline source for the internal loop modules is available at http://rth.dk/resources/mrm.

RNA structural motif recognition based on least-squares distance

RNA structural motifs are recurrent structural elements occurring in RNA molecules. RNA structural motif recognition aims to find RNA substructures that are similar to a query motif, and it is important for RNA structure analysis and RNA function prediction. In view of this, we propose a new method known as RNA Structural Motif Recognition based on Least-Squares distance (LS-RSMR) to effectively recognize RNA structural motifs. A test set consisting of five types of RNA structural motifs occurring in Escherichia coli ribosomal RNA is compiled by us. Experiments are conducted for recognizing these five types of motifs. The experimental results fully reveal the superiority of the proposed LS-RSMR compared with four other state-of-the-art methods.

An RNA folding motif: GNRA tetraloop-receptor interactions

Nearly two decades after Westhof and Michel first proposed that RNA tetraloops may interact with distal helices, tetraloop–receptor interactions have been recognized as ubiquitous elements of RNA tertiary structure. The unique architecture of GNRA tetraloops (N=any nucleotide, R=purine) enables interaction with a variety of receptors, e.g., helical minor grooves and asymmetric internal loops. The most common example of the latter is the GAAA tetraloop–11 nt tetraloop receptor motif. Biophysical characterization of this motif provided evidence for the modularity of RNA structure, with applications spanning improved crystallization methods to RNA tectonics. In this review, we identify and compare types of GNRA tetraloop–receptor interactions. Then we explore the abundance of structural, kinetic, and thermodynamic information on the frequently occurring and most widely studied GAAA tetraloop–11 nt receptor motif. Studies of this interaction have revealed powerful paradigms for structural assembly of RNA, as well as providing new insights into the roles of cations, transition states and protein chaperones in RNA folding pathways. However, further research will clearly be necessary to characterize other tetraloop–receptor and long-range tertiary binding interactions in detail – an important milestone in the quantitative prediction of free energy landscapes for RNA folding.

RNA structure and dynamics: a base pairing perspective

RNA is now known to possess various structural, regulatory and enzymatic functions for survival of cellular organisms. Functional RNA structures are generally created by three-dimensional organization of small structural motifs, formed by base pairing between self-complementary sequences from different parts of the RNA chain. In addition to the canonical Watson-Crick or wobble base pairs, several non-canonical base pairs are found to be crucial to the structural organization of RNA molecules. They appear within different structural motifs and are found to stabilize the molecule through long-range intra-molecular interactions between basic structural motifs like double helices and loops. These base pairs also impart functional variation to the minor groove of A-form RNA helices, thus forming anchoring site for metabolites and ligands. Non-canonical base pairs are formed by edge-to-edge hydrogen bonding interactions between the bases. A large number of theoretical studies have been done to detect and analyze these non-canonical base pairs within crystal or NMR derived structures of different functional RNA. Theoretical studies of these isolated base pairs using ab initio quantum chemical methods as well as molecular dynamics simulations of larger fragments have also established that many of these non-canonical base pairs are as stable as the canonical Watson-Crick base pairs. This review focuses on the various structural aspects of non-canonical base pairs in the organization of RNA molecules and the possible applications of these base pairs in predicting RNA structures with more accuracy.

RNA tectonics (tectoRNA) for RNA nanostructure design and its application in synthetic biology

RNA molecules are versatile biomaterials that act not only as DNA-like genetic materials but also have diverse functions in regulation of cellular biosystems. RNA is capable of regulating gene expression by sequence-specific hybridization. This feature allows the design of RNA-based artificial gene regulators (riboregulators). RNA can also build complex two-dimensional (2D) and 3D nanostructures, which afford protein-like functions and make RNA an attractive material for nanobiotechnology. RNA tectonics is a methodology in RNA nanobiotechnology for the design and construction of RNA nanostructures/nanoobjects through controlled self-assembly of modular RNA units (tectoRNAs). RNA nanostructures designed according to the concept of RNA tectonics are also attractive as tools in synthetic biology, but in vivo RNA tectonics is still in the early stages. This review presents a summary of the achievements of RNA tectonics and its related researches in vitro, and also introduces recent developments that facilitated the use of RNA nanostructures in bacterial cells.

Backbone-base interactions critical to quantum stabilization of transfer RNA anticodon structure

Transfer RNA (tRNA) anticodons adopt a highly ordered 3'-stack without significant base overlap. Density functional theory at the M06-2X/6-31+G(d,p) level in combination with natural bond orbital analysis was utilized to calculate the intramolecular interactions within the tRNA anticodon that are responsible for stabilizing the stair-stepped conformation. Ten tRNA X-ray crystal structures were obtained from the PDB databank and were trimmed to include only the anticodon bases. Hydrogenic positions were added and optimized for the structures in the stair-stepped conformation. The sugar-phosphate backbone has been retained for these calculations, revealing the role it plays in RNA structural stability. It was found that electrostatic interactions between the sugar-phosphate backbone and the base provide the most stability, rather than the traditionally studied interbase stacking. Base-stacking interactions, though present, were weak and inconsistent. Aqueous solvation was found to have little effect on the intramolecular interactions.

Structure and function of the T-loop structural motif in noncoding RNAs

The T-loop is a frequently occurring five-nucleotide motif found in the structure of noncoding RNAs where it is commonly assumed to play an important role in stabilizing the tertiary RNA structure by facilitating long-range interactions between different regions of the molecule. T-loops were first identified in tRNA(Phe) and a formal consensus sequence for this motif was formulated and later revised based on analyses of the crystal structures of prokaryotic ribosomal RNAs and RNase P and the corresponding primary sequence of their orthologues. In the past decade, several new structures of large RNA molecules have been added to the RCSB Protein Data Bank, including the eukaryotic ribosome, a self-splicing group II intron, numerous synthetases in complex with their cognate transfer RNAs (tRNAs), transfer-messenger RNA (tmRNA) in complex with SmpB, several riboswitches, and a complex of bacterial RNase P bound to its tRNA substrate. In this review, the search for T-loops is extended to these new RNA molecules based on the previously established structure-based criteria. The review highlights and discusses the function and additional roles of T-loops in four broad categories of RNA molecules, namely tRNAs, ribosomal RNAs (rRNAs), P RNAs, and RNA genetic elements. Additionally, the potential application for T-loops as interaction modules is also discussed.

The GA-minor submotif as a case study of RNA modularity, prediction, and design

Complex natural RNAs such as the ribosome, group I and group II introns, and RNase P exemplify the fact that three-dimensional (3D) RNA structures are highly modular and hierarchical in nature. Tertiary RNA folding typically takes advantage of a rather limited set of recurrent structural motifs that are responsible for controlling bends or stacks between adjacent helices. Herein, the GA minor and related structural motifs are presented as a case study to highlight several structural and folding principles, to gain further insight into the structural evolution of naturally occurring RNAs, as well as to assist the rational design of artificial RNAs.

Cleavage mediated by the catalytic domain of bacterial RNase P RNA

Like other RNA molecules, RNase P RNA (RPR) is composed of domains, and these have different functions. Here, we provide data demonstrating that the catalytic (C) domain of Escherichia coli (Eco) RPR when separated from the specificity (S) domain mediates cleavage using various model RNA hairpin loop substrates. Compared to full-length Eco RPR, the rate constant, k(obs), of cleavage for the truncated RPR (CP RPR) was reduced 30- to 13,000-fold depending on substrate. Specifically, the structural architecture of the -1/+73 played a significant role where a C(-1)/G(+73) pair had the most dramatic effect on k(obs). Substitution of A(248) (E. coli numbering), positioned near the cleavage site in the RNase P-substrate complex, with G in the CP RPR resulted in 30-fold improvement in rate. In contrast, strengthening the interaction between the RPR and the 3' end of the substrate only had a modest effect. Interestingly, although deleting the S-domain gave a reduction in the rate, it resulted in a less erroneous RPR with respect to cleavage site selection. These data support and extend our understanding of the coupling between the distal interaction between the S-domain and events at the active site. Our findings will also be discussed with respect to the structure of RPR derived from different organisms.

Mapping targetable sites on human telomerase RNA pseudoknot/template domain using 2'-OMe RNA-interacting polynucleotide (RIPtide) microarrays

Most cellular RNAs engage in intrastrand base-pairing that gives rise to complex three-dimensional folds. This self-pairing presents an impediment toward binding of the RNA by nucleic acid-based ligands. An important step in the discovery of RNA-targeting ligands is therefore to identify those regions in a folded RNA that are accessible toward the nucleic acid-based ligand. Because the folding of RNA targets can involve interactions between nonadjacent regions and employ both Watson-Crick and non-Watson-Crick base-pairing, screening of candidate binder ensembles is typically necessary. Microarray-based screening approaches have shown great promise in this regard and have suggested that achieving complete sequence coverage would be a valuable attribute of a next generation system. Here, we report a custom microarray displaying a library of RNA-interacting polynucleotides comprising all possible 2'-OMe RNA sequences from 4- to 8-nucleotides in length. We demonstrate the utility of this array in identifying RNA-interacting polynucleotides that bind tightly and specifically to the highly conserved, functionally essential template/pseudoknot domain of human telomerase RNA and that inhibit telomerase function in vitro.

Cleavage mediated by the P15 domain of bacterial RNase P RNA

Independently folded domains in RNAs frequently adopt identical tertiary structures regardless of whether they are in isolation or are part of larger RNA molecules. This is exemplified by the P15 domain in the RNA subunit (RPR) of the universally conserved endoribonuclease P, which is involved in the processing of tRNA precursors. One of its domains, encompassing the P15 loop, binds to the 3'-end of tRNA precursors resulting in the formation of the RCCA-RNase P RNA interaction (interacting residues underlined) in the bacterial RPR-substrate complex. The function of this interaction was hypothesized to anchor the substrate, expose the cleavage site and result in re-coordination of Mg(2+) at the cleavage site. Here we show that small model-RNA molecules (~30 nt) carrying the P15-loop mediated cleavage at the canonical RNase P cleavage site with significantly reduced rates compared to cleavage with full-size RPR. These data provide further experimental evidence for our model that the P15 domain contributes to both substrate binding and catalysis. Our data raises intriguing evolutionary possibilities for 'RNA-mediated' cleavage of RNA.

Hydration sites of unpaired RNA bases: a statistical analysis of the PDB structures

BACKGROUND

Hydration is crucial for RNA structure and function. X-ray crystallography is the most commonly used method to determine RNA structures and hydration and, therefore, statistical surveys are based on crystallographic results, the number of which is quickly increasing.

RESULTS

A statistical analysis of the water molecule distribution in high-resolution X-ray structures of unpaired RNA nucleotides showed that: different bases have the same penchant to be surrounded by water molecules; clusters of water molecules indicate possible hydration sites, which, in some cases, match those of the major and minor grooves of RNA and DNA double helices; complex hydrogen bond networks characterize the solvation of the nucleotides, resulting in a significant rigidity of the base and its surrounding water molecules. Interestingly, the hydration sites around unpaired RNA bases do not match, in general, the positions that are occupied by the second nucleotide when the base-pair is formed.

CONCLUSIONS

The hydration sites around unpaired RNA bases were found. They do not replicate the atom positions of complementary bases in the Watson-Crick pairs.

Structural, mechanical, and thermodynamic properties of a coarse-grained DNA model

We explore in detail the structural, mechanical, and thermodynamic properties of a coarse-grained model of DNA similar to that recently introduced in a study of DNA nanotweezers [T. E. Ouldridge, A. A. Louis, and J. P. K. Doye, Phys. Rev. Lett. 134, 178101 (2010)]. Effective interactions are used to represent chain connectivity, excluded volume, base stacking, and hydrogen bonding, naturally reproducing a range of DNA behavior. The model incorporates the specificity of Watson-Crick base pairing, but otherwise neglects sequence dependence of interaction strengths, resulting in an "average base" description of DNA. We quantify the relation to experiment of the thermodynamics of single-stranded stacking, duplex hybridization, and hairpin formation, as well as structural properties such as the persistence length of single strands and duplexes, and the elastic torsional and stretching moduli of double helices. We also explore the model's representation of more complex motifs involving dangling ends, bulged bases and internal loops, and the effect of stacking and fraying on the thermodynamics of the duplex formation transition.

Riboswitches: structures and mechanisms

A critical feature of the hypothesized RNA world would have been the ability to control chemical processes in response to environmental cues. Riboswitches present themselves as viable candidates for a sophisticated mechanism of regulatory control in RNA-based life. These regulatory elements in the modern world are most commonly found in the 5'-untranslated regions of bacterial mRNAs, directly interacting with metabolites as a means of regulating expression of the coding region via a secondary structural switch. In this review, we focus on recent insights into how these RNAs fold into complex architectures capable of both recognizing a specific small molecule compound and exerting regulatory control over downstream sequences, with an emphasis on transcriptional regulation.

RAG: an update to the RNA-As-Graphs resource

Background

In 2004, we presented a web resource for stimulating the search for novel RNAs, RNA-As-Graphs (RAG), which classified, catalogued, and predicted RNA secondary structure motifs using clustering and build-up approaches. With the increased availability of secondary structures in recent years, we update the RAG resource and provide various improvements for analyzing RNA structures.

Description

Our RAG update includes a new supervised clustering algorithm that can suggest RNA motifs that may be "RNA-like". We use this utility to describe RNA motifs as three classes: existing, RNA-like, and non-RNA-like. This produces 126 tree and 16,658 dual graphs as candidate RNA-like topologies using the supervised clustering algorithm with existing RNAs serving as the training data. A comparison of this clustering approach to an earlier method shows considerable improvements. Additional RAG features include greatly expanded search capabilities, an interface to better utilize the benefits of relational database, and improvements to several of the utilities such as directed/labeled graphs and a subgraph search program.

Conclusions

The RAG updates presented here augment the database's intended function - stimulating the search for novel RNA functionality - by classifying available motifs, suggesting new motifs for design, and allowing for more specific searches for specific topologies. The updated RAG web resource offers users a graph-based tool for exploring available RNA motifs and suggesting new RNAs for design.

A quantitative measure of chirality inside nucleic acid databank

We show the capability of a chirality index (Pietropaolo et al., Proteins 2008;70:667-677) to investigate nucleic acid structures because of its high sensitivity to helical conformations. By analyzing selected structures of DNA and RNA, we have found that sequences rich in cytosine and guanine have a tendency to left-handed chirality, in contrast to regions rich in adenine or thymine which show strong negative, right-handed, chirality values. We also analyze RNA structures, where specific loops and hairpin motifs are characterized by a well-defined chirality value. We find that in nucleosome the chirality is exalted, whereas in ribosome it is reduced. Our results illustrate the sensitivity of this descriptor for nucleic acid conformations.

An expanded collection and refined consensus model of glmS ribozymes

Self-cleaving glmS ribozymes selectively bind glucosamine-6-phosphate (GlcN6P) and use this metabolite as a cofactor to promote self-cleavage by internal phosphoester transfer. Representatives of the glmS ribozyme class are found in Gram-positive bacteria where they reside in the 5' untranslated regions (UTRs) of glmS messenger RNAs that code for the essential enzyme L-glutamine:D-fructose-6-phosphate aminotransferase. By using comparative sequence analyses, we have expanded the number of glmS ribozyme representatives from 160 to 463. All but two glmS ribozymes are present in glmS mRNAs and most exhibit striking uniformity in sequence and structure, which are features that make representatives attractive targets for antibacterial drug development. However, our discovery of rare variants broadens the consensus sequence and structure model. For example, in the Deinococcus-Thermus phylum, several structural variants exist that carry additional stems within the catalytic core and changes to the architecture of core-supporting substructures. These findings reveal that glmS ribozymes have a broader phylogenetic distribution than previously known and suggest that additional rare structural variants may remain to be discovered.

Engineering cooperative tecto-RNA complexes having programmable stoichiometries

High affinity and specificity RNA-RNA binding interfaces can be constructed by combining pairs of GNRA loop/loop-receptor interaction motifs. These interactions can be fused using flexible four-way junction motifs to create divalent, self-assembling scaffolding units ('tecto-RNA') that have favorable properties for nanomedicine and other applications. We describe the design and directed assembly of tecto-RNA units ranging from closed, cooperatively assembling ring-shaped complexes of programmable stoichiometries (dimers, trimers and tetramers) to open multimeric structures. The novelty of this work is that tuning of the stoichiometries of self-assembled complexes is achieved by precise positioning of the interaction motifs in the monomer units rather than changing their binding specificities. Structure-probing and transmission electron microscopy studies as well as thermodynamic analysis support formation of closed cooperative complexes that are highly resistant to nuclease digestion. The present designs provide two helical arms per RNA monomer for further functionalization aims.