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.     

In-gel probing of individual RNA conformers within a mixed population reveals a dimerization structural switch in the HIV-1 leader

Definitive secondary structural mapping of RNAs in vitro can be complicated by the presence of more than one structural conformer or multimerization of some of the molecules. Until now, probing a single structure of conformationally flexible RNA molecules has typically relied on introducing stabilizing mutations or adjusting buffer conditions or RNA concentration. Here, we present an in-gel SHAPE (selective 2′OH acylation analysed by primer extension) approach, where a mixed structural population of RNA molecules is separated by non-denaturing gel electrophoresis and the conformers are individually probed within the gel matrix. Validation of the technique using a well-characterized RNA stem-loop structure, the HIV-1 trans-activation response element, showed that authentic structure was maintained and that the method was accurate and highly reproducible. To further demonstrate the utility of in-gel SHAPE, we separated and examined monomeric and dimeric species of the HIV-1 packaging signal RNA. Extensive differences in acylation sensitivity were seen between monomer and dimer. The results support a recently proposed structural switch model of RNA genomic dimerization and packaging, and demonstrate the discriminatory power of in-gel SHAPE.     

Ab initio RNA folding by discrete molecular dynamics: from structure prediction to folding mechanisms

RNA molecules with novel functions have revived interest in the accurate prediction of RNA three-dimensional (3D) structure and folding dynamics. However, existing methods are inefficient in automated 3D structure prediction. Here, we report a robust computational approach for rapid folding of RNA molecules. We develop a simplified RNA model for discrete molecular dynamics (DMD) simulations, incorporating base-pairing and base-stacking interactions. We demonstrate correct folding of 150 structurally diverse RNA sequences. The majority of DMD-predicted 3D structures have <4 A deviations from experimental structures. The secondary structures corresponding to the predicted 3D structures consist of 94% native base-pair interactions. Folding thermodynamics and kinetics of tRNA(Phe), pseudoknots, and mRNA fragments in DMD simulations are in agreement with previous experimental findings. Folding of RNA molecules features transient, non-native conformations, suggesting non-hierarchical RNA folding. Our method allows rapid conformational sampling of RNA folding, with computational time increasing linearly with RNA length. We envision this approach as a promising tool for RNA structural and functional analyses.     

Electron microscopy of ribosomal RNA mounted by adsorption in the absence of protein.

A method is described for visualizing ribosomal RNA in the electron microscope by an adsorption technique. The visual configuration of the molecules is sensitive to the composition of the medium from which they are adsorbed, and is reproducible under a given set of conditions. Regular structural features occur in the larger rRNA species similar to those observed with protein monolayer techniques. The dimensions of adsorbed molecules, however, are more consistent with current models of ribosomal RNA as it exists in solution.     

CapR: revealing structural specificities of RNA-binding protein target recognition using CLIP-seq data

RNA-binding proteins (RBPs) bind to their target RNA molecules by recognizing specific RNA sequences and structural contexts. The development of CLIP-seq and related protocols has made it possible to exhaustively identify RNA fragments that bind to RBPs. However, no efficient bioinformatics method exists to reveal the structural specificities of RBP–RNA interactions using these data. We present CapR, an efficient algorithm that calculates the probability that each RNA base position is located within each secondary structural context. Using CapR, we demonstrate that several RBPs bind to their target RNA molecules under specific structural contexts. CapR is available at https://sites.google.com/site/fukunagatsu/software/capr.     

QRNAstruct: a method for extracting secondary structural features of RNA via regression with biological activity

Recent technological advances have enabled the generation of large amounts of data consisting of RNA sequences and their functional activity. Here, we propose a method for extracting secondary structure features that affect the functional activity of RNA from sequence–activity data. Given pairs of RNA sequences and their corresponding bioactivity values, our method calculates position-specific structural features of the input RNA sequences, considering every possible secondary structure of each RNA. A Ridge regression model is trained using the structural features as feature vectors and the bioactivity values as response variables. Optimized model parameters indicate how secondary structure features affect bioactivity. We used our method to extract intramolecular structural features of bacterial translation initiation sites and self-cleaving ribozymes, and the intermolecular features between rRNAs and Shine–Dalgarno sequences and between U1 RNAs and splicing sites. We not only identified known structural features but also revealed more detailed insights into structure–activity relationships than previously reported. Importantly, the datasets we analyzed here were obtained from different experimental systems and differed in size, sequence length and similarity, and number of RNA molecules involved, demonstrating that our method is applicable to various types of data consisting of RNA sequences and bioactivity values.     

RNA templating of molecular assembly and covalent modification patterning in early molecular evolution and modern biosystems

The Direct RNA Template (DRT) hypothesis proposes that an early stage of genetic code evolution involved RNA molecules acting as stereochemical recognition templates for assembly of specific amino acids in sequence-ordered arrays, providing a framework for directed covalent peptide bond formation. It is hypothesized here that modern biological precedents may exist for RNA-based structural templating with functional analogies to hypothetical DRT systems. Beyond covalent molecular assembly, an extension of the DRT concept can include RNA molecules acting as dynamic structural template guides for the specific non-covalent assembly of multi-subunit complexes, equivalent to structural assembly chaperones. However, despite numerous precedents for RNA molecules acting as scaffolds for protein complexes, true RNA-mediated assembly chaperoning appears to be absent in modern biosystems. Another level of function with parallels to a DRT system is possible if RNA structural motifs dynamically guided specific patterns of catalytic modifications within multiple target sites in a pre-formed polymer or macromolecular complex. It is suggested that this type of structural RNA templating could logically play a functional role in certain areas of biology, one of which is the glycome of complex organisms. If any such RNA templating processes are shown to exist, they would share no necessary evolutionary relationships with events during early molecular evolution, but may promote understanding of the practical limits of biological RNA functions now and in the ancient RNA World. Awareness of these formal possibilities may also assist in the current search for functions of extensive non-coding RNAs in complex organisms, or for efforts towards artificial rendering of DRT systems.     

Pure multiple RNA secondary structure alignments: a progressive profile approach

In functional, noncoding RNA, structure is often essential to function. While the full 3D structure is very difficult to determine, the 2D structure of an RNA molecule gives good clues to its 3D structure, and for molecules of moderate length, it can be predicted with good reliability. Structure comparison is, in analogy to sequence comparison, the essential technique to infer related function. We provide a method for computing multiple alignments of RNA secondary structures under the tree alignment model, which is suitable to cluster RNA molecules purely on the structural level, i.e., sequence similarity is not required. We give a systematic generalization of the profile alignment method from strings to trees and forests. We introduce a tree profile representation of RNA secondary structure alignments which allows reasonable scoring in structure comparison. Besides the technical aspects, an RNA profile is a useful data structure to represent multiple structures of RNA sequences. Moreover, we propose a visualization of RNA consensus structures that is enriched by the full sequence information.     

Identifying constraints on the higher-order structure of RNA: continued development and application of comparative sequence analysis methods.

Comparative sequence analysis addresses the problem of RNA folding and RNA structural diversity, and is responsible for determining the folding of many RNA molecules, including 5S, 16S, and 23S rRNAs, tRNA, RNAse P RNA, and Group I and II introns. Initially this method was utilized to fold these sequences into their secondary structures. More recently, this method has revealed numerous tertiary correlations, elucidating novel RNA structural motifs, several of which have been experimentally tested and verified, substantiating the general application of this approach. As successful as the comparative methods have been in elucidating higher-order structure, it is clear that additional structure constraints remain to be found. Deciphering such constraints requires more sensitive and rigorous protocols, in addition to RNA sequence datasets that contain additional phylogenetic diversity and an overall increase in the number of sequences. Various RNA databases, including the tRNA and rRNA sequence datasets, continue to grow in number as well as diversity. Described herein is the development of more rigorous comparative analysis protocols. Our initial development and applications on different RNA datasets have been very encouraging. Such analyses on tRNA, 16S and 23S rRNA are substantiating previously proposed associations and are now beginning to reveal additional constraints on these molecules. A subset of these involve several positions that correlate simultaneously with one another, implying units larger than a basepair can be under a phylogenetic constraint.     

Molecular mimicry: quantitative methods to study structural similarity between protein and RNA

With rapidly increasing availability of three-dimensional structures, one major challenge for the post-genome era is to infer the functions of biological molecules based on their structural similarity. While quantitative studies of structural similarity between the same type of biological molecules (e.g., protein vs. protein) have been carried out intensively, the comparable study of structural similarity between different types of biological molecules (e.g., protein vs. RNA) remains unexplored. Here we have developed a new bioinformatics approach to quantitatively study the structural similarity between two different types of biopolymers--proteins and RNA--based on the spatial distribution of conserved elements. We applied it to two previously proposed tRNA-protein mimicry pairs whose functional relatedness between two molecules has been recently determined experimentally. Our method detected the biologically meaningful signals, which are consistent with experimental evidence.     

Nanomanipulation of Single RNA Molecules by Optical Tweezers

A large portion of the human genome is transcribed but not translated. In this post genomic era, regulatory functions of RNA have been shown to be increasingly important. As RNA function often depends on its ability to adopt alternative structures, it is difficult to predict RNA three-dimensional structures directly from sequence. Single-molecule approaches show potentials to solve the problem of RNA structural polymorphism by monitoring molecular structures one molecule at a time. This work presents a method to precisely manipulate the folding and structure of single RNA molecules using optical tweezers. First, methods to synthesize molecules suitable for single-molecule mechanical work are described. Next, various calibration procedures to ensure the proper operations of the optical tweezers are discussed. Next, various experiments are explained. To demonstrate the utility of the technique, results of mechanically unfolding RNA hairpins and a single RNA kissing complex are used as evidence. In these examples, the nanomanipulation technique was used to study folding of each structural domain, including secondary and tertiary, independently. Lastly, the limitations and future applications of the method are discussed.     

Getting to the end of RNA: Structural analysis of protein recognition of 5′ and 3′ termini

The specific recognition by proteins of the 5′ and 3′ ends of RNA molecules is an important facet of many cellular processes, including RNA maturation, regulation of translation initiation and control of gene expression by degradation and RNA interference. The aim of this review is to survey recent structural analyses of protein binding domains that specifically bind to the extreme 5′ or 3′ termini of RNA. For reasons of space and because their interactions are also governed by catalytic considerations, we have excluded enzymes that modify the 5′ and 3′ extremities of RNA. It is clear that there is enormous structural diversity among the proteins that have evolved to bind to the ends of RNA molecules. Moreover, they commonly exhibit conformational flexibility that appears to be important for binding and regulation of the interaction. This flexibility has sometimes complicated the interpretation of structural results and presents significant challenges for future investigations.     

Site-specific labeling of RNA with fluorophores and other structural probes

Site-specific probes provide a powerful tool for structure and function studies of nucleic acids, especially in elucidating tertiary structures of large ribozymes and other folded RNA molecules. Among many types of extrinsic labels, fluorophores are most attractive because they can provide structural information at millisecond time resolution, thus allowing real-time observation of structural transition during biological function. Methods for introducing fluorophores in RNA molecules are summarized here. These methods are robust and readily applicable to the labeling of other types of probes. However, as each case of RNA modification is unique, fine tuning of the general methodology is beneficial.     

The use of calorimetry in the biophysical characterization of small molecule alkaloids binding to RNA structures

BackgroundRNA has now emerged as a potential target for therapeutic intervention. RNA targeted drug design requires detailed thermodynamic characterization that provides new insights into the interactions and this together with structural data, may be used in rational drug design. The use of calorimetry to characterize small molecule–RNA interactions has emerged as a reliable and sensitive tool after the recent advancements in biocalorimetry.Scope of the reviewThis review summarizes the recent advancements in thermodynamic characterization of small molecules, particularly some natural alkaloids binding to various RNA structures. Thermodynamic characterization provides information that can supplement structural data leading to more effective drug development protocols.Major conclusionsThis review provides a concise report on the use of isothermal titration calorimetry (ITC) and differential scanning calorimetry (DSC) techniques in characterizing small molecules, mostly alkaloids–RNA interactions with particular reference to binding of tRNA, single stranded RNA, double stranded RNA, poly(A), triplex RNA.General significanceIt is now apparent that a combination of structural and thermodynamic data is essential for rational design of specific RNA targeted drugs. Recent advancements in biocalorimetry instrumentation have led to detailed understanding of the thermodynamics of small molecules binding to various RNA structures paving the path for the development of many new natural and synthetic molecules as specific binders to various RNA structures. RNA targeted drug design, that remained unexplored, will immensely benefit from the calorimetric studies leading to the development of effective drugs for many diseases. This article is part of a Special Issue entitled Microcalorimetry in the BioSciences — Principles and Applications, edited by Fadi Bou-Abdallah.     

RNAMotifScan: automatic identification of RNA structural motifs using secondary structural alignment

Recent studies have shown that RNA structural motifs play essential roles in RNA folding and interaction with other molecules. Computational identification and analysis of RNA structural motifs remains a challenging task. Existing motif identification methods based on 3D structure may not properly compare motifs with high structural variations. Other structural motif identification methods consider only nested canonical base-pairing structures and cannot be used to identify complex RNA structural motifs that often consist of various non-canonical base pairs due to uncommon hydrogen bond interactions. In this article, we present a novel RNA structural alignment method for RNA structural motif identification, RNAMotifScan, which takes into consideration the isosteric (both canonical and non-canonical) base pairs and multi-pairings in RNA structural motifs. The utility and accuracy of RNAMotifScan is demonstrated by searching for kink-turn, C-loop, sarcin-ricin, reverse kink-turn and E-loop motifs against a 23S rRNA (PDBid: 1S72), which is well characterized for the occurrences of these motifs. Finally, we search these motifs against the RNA structures in the entire Protein Data Bank and the abundances of them are estimated. RNAMotifScan is freely available at our supplementary website (http://genome.ucf.edu/RNAMotifScan).     

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.     

A method for discovering common patterns from two RNA secondary structures and its application to structural repeat detection

We propose an ab initio method, named DiscoverR, for finding common patterns from two RNA secondary structures. The method works by representing RNA secondary structures as ordered labeled trees and performs tree pattern discovery using an efficient dynamic programming algorithm. DiscoverR is able to identify and extract the largest common substructures from two RNA molecules having different sizes without prior knowledge of the locations and topologies of these substructures. We also extend DiscoverR to find repeated regions in an RNA secondary structure, and apply this extended method to detect structural repeats in the 3'-untranslated region of a protein kinase gene. We describe the biological significance of a repeated hairpin found by our method, demonstrating the usefulness of the method. DiscoverR is implemented in Java; a jar file including the source code of the program is available for download at http://bioinformatics.njit.edu/DiscoverR.     

Ribonucleopeptides: functional RNA-peptide complexes

Selection of functional RNAs from randomized pool of RNA molecules successfully affords RNA aptamers that specifically bind to small molecules, and that have catalytic activities. Recent structural analyses of the ribosomal RNA complex suggest that the RNA-protein complex would be a new structural candidate for the design of tailor-made receptors and enzymes. We have designed an ATP binding domain that consists of an RNA subunit and a peptide subunit by means of structure-based design approach and successive in vitro selection method. The RNA subunit is designed to consist of two functional domains; an ATP binding domain with 20 randomized nucleotides and an adjacent stem region that serves as a binding site for the RNA-binding peptide. The randomized nucleotide region was placed next to the HIV-1 Rev response element to enable the formation of "ribonucleopeptide" pools in the presence of the Rev peptide. In vitro selection of RNA oligonucleotides from the randomized pool afforded a ribonucleopeptide receptor specific for ATP. The ATP-binding ribonucleopeptide did not share the known consensus nucleotide sequence for ATP aptamers, and completely lost its ATP-binding ability in the absence of the Rev peptide. The ATP-binding activity of the ribonucleopeptide was increased by a substitution of the N-terminal amino acid of the Rev peptide. These results demonstrate that the peptide stabilizes the functional structure of RNA and suggest that amino acids outside the RNA binding region of the peptide participate in the ATP binding. Our approach would provide a new strategy for the design of tailor-made ribonucleopeptide receptors.     

Scanning tunnelling microscopy of 16S ribosomal RNA in water

The scanning tunnelling microscope has been used to image 16S ribosomal RNA molecules in water electrophoretically deposited on graphite surface. Two kinds of images have been obtained: images showing aggregates of 16S ribosomal RNA molecules similar to those obtained from DNA solutions and others showing individual 16S ribosomal RNA molecules. An interesting characteristic of these images, recorded in constant current mode, is that the 16S ribosomal RNA molecules appear to be located below the graphite surface. The morphology and several structural parameters of the molecules were consistent with the data obtained from electron microscopy.