A common motif organizes the structure of multi-helix loops in 16 S and 23 S ribosomal RNAs

Archaeological approaches to RNA virus evolution

Studies with RNA enzymes (ribozymes) and protein enzymes have identified certain structural elements that are present in some cellular mRNAs and viral RNAs. These elements do not share a primary structure and, thus, are not phylogenetically related. However, they have common (secondary/tertiary) structural folds that, according to some lines of evidence, may have an ancient and common origin. The term 'mRNA archaeology' has been coined to refer to the search for such structural/functional relics that may be informative of early evolutionary developments in the cellular and viral worlds and have lasted to the present day. Such identified RNA elements may have developed as biological signals with structural and functional relevance (as if they were buried objects with archaeological value), and coexist with the standard linear information of nucleic acid molecules that is translated into proteins. However, there is a key difference between the methods that extract information from either the primary structure of mRNA or the signals provided by secondary and tertiary structures. The former (sequence comparison and phylogenetic analysis) requires strict continuity of the material vehicle of information during evolution, whereas the archaeological method does not require such continuity. The tools of RNA archaeology (including the use of ribozymes and enzymes to investigate the reactivity of the RNA elements) establish links between the concepts of communication and language theories that have not been incorporated into knowledge of virology, as well as experimental studies on the search for functionally relevant RNA structures.

Graduate Student Literature Review: Hemotropic mycoplasmas in cattle

The objective of this narrative literature review is to better understand bovine hemoplasmosis, an emerging disease that threatens dairy animal health. Several species of hemotropic mycoplasma are known to infect both animals and humans, and Mycoplasma wenyonii and Candidatus Mycoplasma haemobos are the species that infect red blood cells of cattle. These microorganisms are associated with clinical signs in dairy cattle, but the effects of infection on health and productivity of dairy cows are poorly understood. In this paper, we review information about the epidemiology of bovine hemoplasmosis in different countries, including clinical signs associated with hemoplasmosis in cattle, methods of diagnosis, treatment, possible routes of transmission, risk factors for infection, and disease progression. Although hemoplasmas have been reported to infect cattle in many countries, and methods used to detect these organisms have improved, numerous gaps in knowledge were identified. The pathogenesis of the disease and potential effect on animal health and productivity remain unclear. With this review, we seek to contribute to the understanding of hemoplasmosis in cattle and provide insights for further research to improve disease management strategies and overall animal health in the dairy industry.

Genetic disruption of the bacterial raiA motif noncoding RNA causes defects in sporulation and aggregation

Several structured noncoding RNAs in bacteria are essential contributors to fundamental cellular processes. Thus, discoveries of additional ncRNA classes provide opportunities to uncover and explore biochemical mechanisms relevant to other major and potentially ancient processes. A candidate structured ncRNA named the "raiA motif" has been found via bioinformatic analyses in over 2,500 bacterial species. The gene coding for the RNA typically resides between the raiA and comFC genes of many species of Bacillota and Actinomycetota. Structural probing of the raiA motif RNA from the Gram-positive anaerobe Clostridium acetobutylicum confirms key features of its sophisticated secondary structure model. Expression analysis of raiA motif RNA reveals that the RNA is constitutively produced but reaches peak abundance during the transition from exponential growth to stationary phase. The raiA motif RNA becomes the fourth most abundant RNA in C. acetobutylicum, excluding ribosomal RNAs and transfer RNAs. Genetic disruption of the raiA motif RNA causes cells to exhibit substantially decreased spore formation and diminished ability to aggregate. Restoration of normal cellular function in this knock-out strain is achieved by expression of a raiA motif gene from a plasmid. These results demonstrate that raiA motif RNAs normally participate in major cell differentiation processes by operating as a trans-acting factor.

When will RNA get its AlphaFold moment?

The protein structure prediction problem has been solved for many types of proteins by AlphaFold. Recently, there has been considerable excitement to build off the success of AlphaFold and predict the 3D structures of RNAs. RNA prediction methods use a variety of techniques, from physics-based to machine learning approaches. We believe that there are challenges preventing the successful development of deep learning-based methods like AlphaFold for RNA in the short term. Broadly speaking, the challenges are the limited number of structures and alignments making data-hungry deep learning methods unlikely to succeed. Additionally, there are several issues with the existing structure and sequence data, as they are often of insufficient quality, highly biased and missing key information. Here, we discuss these challenges in detail and suggest some steps to remedy the situation. We believe that it is possible to create an accurate RNA structure prediction method, but it will require solving several data quality and volume issues, usage of data beyond simple sequence alignments, or the development of new less data-hungry machine learning methods.

Sarcin/Ricin Domain RNA Retains Its Structure Better Than A-RNA in Adaptively Biased Molecular Dynamics Simulations

Less than one in thirty of the RNA sequences transcribed in humans are translated into protein. The noncoding RNA (ncRNA) functions in catalysis, structure, regulation, and more. However, for the most part, these functions are poorly characterized. RNA is modular and described by motifs that include helical A-RNA with canonical Watson-Crick base-pairing as well as structures with only noncanonical base pairs. Understanding the structure and dynamics of motifs will aid in deciphering functions of specific ncRNAs. We present computational studies on a standard sarcin/ricin domain (SRD), citrus bark cracking viroid SRD, as well as A-RNA. We have applied enhanced molecular dynamics techniques that construct an inverse free-energy surface (iFES) determined by collective variables that monitor base-pairing and backbone conformation. Each SRD RNA is flanked on each side by A-RNA, allowing comparison of the behavior of these motifs in the same molecule. The RNA iFESs have single peaks, indicating that the combined motifs should denature as a single cohesive unit, rather than by regional melting. Local root-mean-square deviation (RMSD) analysis and communication propensity (CProp, variance in distances between residue pairs) reveal distinct motif properties. Our analysis indicates that the standard SRD is more stable than the viroid SRD, which is more stable than A-RNA. Base pairs at SRD to A-RNA transitions have limited flexibility. Application of CProp reveals extraordinary stiffness of the SRD, allowing residues on opposite sides of the motif to sense each other's motions.

RNAloops: a database of RNA multiloops

MOTIVATION

Knowledge of the 3D structure of RNA supports discovering its functions and is crucial for designing drugs and modern therapeutic solutions. Thus, much attention is devoted to experimental determination and computational prediction targeting the global fold of RNA and its local substructures. The latter include multi-branched loops-functionally significant elements that highly affect the spatial shape of the entire molecule. Unfortunately, their computational modeling constitutes a weak point of structural bioinformatics. A remedy for this is in collecting these motifs and analyzing their features.

RESULTS

RNAloops is a self-updating database that stores multi-branched loops identified in the PDB-deposited RNA structures. A description of each loop includes angular data-planar and Euler angles computed between pairs of adjacent helices to allow studying their mutual arrangement in space. The system enables search and analysis of multiloops, presents their structure details numerically and visually, and computes data statistics.

AVAILABILITY AND IMPLEMENTATION

RNAloops is freely accessible at https://rnaloops.cs.put.poznan.pl.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Sub-3-Å cryo-EM structure of RNA enabled by engineered homomeric self-assembly

High-resolution structural studies are essential for understanding the folding and function of diverse RNAs. Herein, we present a nanoarchitectural engineering strategy for efficient structural determination of RNA-only structures using single-particle cryogenic electron microscopy (cryo-EM). This strategy-ROCK (RNA oligomerization-enabled cryo-EM via installing kissing loops)-involves installing kissing-loop sequences onto the functionally nonessential stems of RNAs for homomeric self-assembly into closed rings with multiplied molecular weights and mitigated structural flexibility. ROCK enables cryo-EM reconstruction of the Tetrahymena group I intron at 2.98-Å resolution overall (2.85 Å for the core), allowing de novo model building of the complete RNA, including the previously unknown peripheral domains. ROCK is further applied to two smaller RNAs-the Azoarcus group I intron and the FMN riboswitch, revealing the conformational change of the former and the bound ligand in the latter. ROCK holds promise to greatly facilitate the use of cryo-EM in RNA structural studies.

A minimal yet flexible likelihood framework to assess correlated evolution

An evolutionary process is reflected in the sequence of changes of any trait (e.g. morphological or molecular) through time. Yet, a better understanding of evolution would be procured by characterizing correlated evolution, or when two or more evolutionary processes interact. Previously developed parametric methods often require significant computing time as they rely on the estimation of many parameters. Here we propose a minimal likelihood framework modelling the joint evolution of two traits on a known phylogenetic tree. The type and strength of correlated evolution is characterized by a few parameters tuning mutation rates of each trait and interdependencies between these rates. The framework can be applied to study any discrete trait or character ranging from nucleotide substitution to gain or loss of a biological function. More specifically, it can be used to 1) test for independence between two evolutionary processes, 2) identify the type of interaction between them and 3) estimate parameter values of the most likely model of interaction. In the current implementation, the method takes as input a phylogenetic tree with discrete evolutionary events mapped on its branches. The method then maximizes the likelihood for one or several chosen scenarios. The strengths and limits of the method, as well as its relative power compared to a few other methods, are assessed using both simulations and data from 16S rRNA sequences in a sample of 54 γ-enterobacteria. We show that, even with datasets of fewer than 100 species, the method performs well in parameter estimation and in evolutionary model selection.

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

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

R2DT is a framework for predicting and visualising RNA secondary structure using templates

Non-coding RNAs (ncRNA) are essential for all life, and their functions often depend on their secondary (2D) and tertiary structure. Despite the abundance of software for the visualisation of ncRNAs, few automatically generate consistent and recognisable 2D layouts, which makes it challenging for users to construct, compare and analyse structures. Here, we present R2DT, a method for predicting and visualising a wide range of RNA structures in standardised layouts. R2DT is based on a library of 3,647 templates representing the majority of known structured RNAs. R2DT has been applied to ncRNA sequences from the RNAcentral database and produced >13 million diagrams, creating the world's largest RNA 2D structure dataset. The software is amenable to community expansion, and is freely available at https://github.com/rnacentral/R2DT and a web server is found at https://rnacentral.org/r2dt .

On the Nature of Nucleobase Stacking in RNA: A Comprehensive Survey of Its Structural Variability and a Systematic Classification of Associated Interactions

The astonishing diversity in folding patterns of RNA three-dimensional (3D) structures is crafted by myriads of noncovalent contacts, of which base pairing and stacking are the most prominent. A systematic and comprehensive classification and annotation of these interactions is necessary for a molecular-level understanding of their roles. However, unlike in the case of base pairing, where a widely accepted nomenclature and classification scheme exists in the public domain, currently available classification schemes for base-base stacking need major enhancements to comprehensively capture the necessary features underlying the rich stacking diversity in RNA. Here, we extend the previous stacking classification based on nucleobase interacting faces by introducing a structurally intuitive geometry-cum topology-based scheme. Specifically, a stack is first classified in terms of the geometry described by the relative orientation of the glycosidic bonds, which generates eight basic stacking geometric families for heterodimeric stacks and six of those for homodimeric stacks. Further annotation in terms of the identity of the bases and the region of involvement of purines (five-membered, six-membered, or both rings) leads to the enumeration of 384 distinct RNA base stacks. Based on our classification scheme, we present an algorithm for automated identification of stacks in RNA crystal structures and analyze the stacking context in selected RNA structures. Overall, the work described here is expected to greatly facilitate the structure-based RNA research.

Neocles B. Leontis (1955 - 2020)

None.

Unboxing the T-box riboswitches-A glimpse into multivalent and multimodal RNA-RNA interactions

The T-box riboswitches are widespread bacterial noncoding RNAs that directly bind specific tRNAs, sense aminoacylation on bound tRNAs, and switch conformations to control amino-acid metabolism and to maintain nutritional homeostasis. The core mechanisms of tRNA recognition, amino acid sensing, and conformational switching by the T-boxes have been recently elucidated, providing a wealth of new insights into multivalent and multimodal RNA-RNA interactions. This review dissects the structures and tRNA-recognition mechanisms by the Stem I, Stem II, and Discriminator domains, which collectively compose the T-box riboswitches. It further compares and contrasts the two classes of T-boxes that regulate transcription and translation, respectively, and integrates recent findings to derive general themes, trends, and insights into complex RNA-RNA interactions. Specifically, the T-box paradigm reveals that noncoding RNAs can interact with each other through multiple coordinated contacts, concatenation of stacked helices, and mutually induced fit. Numerous tertiary contacts, especially those emanating from strings of single-stranded purines, act in concert to reinforce long-range base-pairing and stacking interactions. These coordinated, mixed-mode contacts allow the T-box RNA to sterically sense aminoacylation on the tRNA using a bipartite steric sieve, and to couple this readout to a conformational switch mediated by tRNA-T-box stacking. Together, the insights gleaned from the T-box riboswitches inform investigations into other complex RNA structures and assemblies, development of T-box-targeted antimicrobials, and may inspire design and engineering of novel RNA sensors, regulators, and interfaces. This article is categorized under: RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs Regulatory RNAs/RNAi/Riboswitches > Riboswitches.

Discovery of 20 novel ribosomal leader candidates in bacteria and archaea

BACKGROUND

RNAs perform many functions in addition to supplying coding templates, such as binding proteins. RNA-protein interactions are important in multiple processes in all domains of life, and the discovery of additional protein-binding RNAs expands the scope for studying such interactions. To find such RNAs, we exploited a form of ribosomal regulation. Ribosome biosynthesis must be tightly regulated to ensure that concentrations of rRNAs and ribosomal proteins (r-proteins) match. One regulatory mechanism is a ribosomal leader (r-leader), which is a domain in the 5' UTR of an mRNA whose genes encode r-proteins. When the concentration of one of these r-proteins is high, the protein binds the r-leader in its own mRNA, reducing gene expression and thus protein concentrations. To date, 35 types of r-leaders have been validated or predicted.

RESULTS

By analyzing additional conserved RNA structures on a multi-genome scale, we identified 20 novel r-leader structures. Surprisingly, these included new r-leaders in the highly studied organisms Escherichia coli and Bacillus subtilis. Our results reveal several cases where multiple unrelated RNA structures likely bind the same r-protein ligand, and uncover previously unknown r-protein ligands. Each r-leader consistently occurs upstream of r-protein genes, suggesting a regulatory function. That the predicted r-leaders function as RNAs is supported by evolutionary correlations in the nucleotide sequences that are characteristic of a conserved RNA secondary structure. The r-leader predictions are also consistent with the locations of experimentally determined transcription start sites.

CONCLUSIONS

This work increases the number of known or predicted r-leader structures by more than 50%, providing additional opportunities to study structural and evolutionary aspects of RNA-protein interactions. These results provide a starting point for detailed experimental studies.

Root of the Tree: The Significance, Evolution, and Origins of the Ribosome

The ribosome is an ancient molecular fossil that provides a telescope to the origins of life. Made from RNA and protein, the ribosome translates mRNA to coded protein in all living systems. Universality, economy, centrality and antiquity are ingrained in translation. The translation machinery dominates the set of genes that are shared as orthologues across the tree of life. The lineage of the translation system defines the universal tree of life. The function of a ribosome is to build ribosomes; to accomplish this task, ribosomes make ribosomal proteins, polymerases, enzymes, and signaling proteins. Every coded protein ever produced by life on Earth has passed through the exit tunnel, which is the birth canal of biology. During the root phase of the tree of life, before the last common ancestor of life (LUCA), exit tunnel evolution is dominant and unremitting. Protein folding coevolved with evolution of the exit tunnel. The ribosome shows that protein folding initiated with intrinsic disorder, supported through a short, primitive exit tunnel. Folding progressed to thermodynamically stable β-structures and then to kinetically trapped α-structures. The latter were enabled by a long, mature exit tunnel that partially offset the general thermodynamic tendency of all polypeptides to form β-sheets. RNA chaperoned the evolution of protein folding from the very beginning. The universal common core of the ribosome, with a mass of nearly 2 million Daltons, was finalized by LUCA. The ribosome entered stasis after LUCA and remained in that state for billions of years. Bacterial ribosomes never left stasis. Archaeal ribosomes have remained near stasis, except for the superphylum Asgard, which has accreted rRNA post LUCA. Eukaryotic ribosomes in some lineages appear to be logarithmically accreting rRNA over the last billion years. Ribosomal expansion in Asgard and Eukarya has been incremental and iterative, without substantial remodeling of pre-existing basal structures. The ribosome preserves information on its history.

High-affinity recognition of specific tRNAs by an mRNA anticodon-binding groove

T-box riboswitches are modular bacterial noncoding RNAs that sense and regulate amino acid availability through direct interactions with tRNAs. Between the 5' anticodon-binding stem I domain and the 3' amino acid sensing domains of most T-boxes lies the stem II domain of unknown structure and function. Here, we report a 2.8-Å cocrystal structure of the Nocardia farcinica ileS T-box in complex with its cognate tRNA^Ile. The structure reveals a perpendicularly arranged ultrashort stem I containing a K-turn and an elongated stem II bearing an S-turn. Both stems rest against a compact pseudoknot, dock via an extended ribose zipper and jointly create a binding groove specific to the anticodon of its cognate tRNA. Contrary to proposed distal contacts to the tRNA elbow region, stem II locally reinforces the codon-anticodon interactions between stem I and tRNA, achieving low-nanomolar affinity. This study illustrates how mRNA junctions can create specific binding sites for interacting RNAs of prescribed sequence and structure.

A functional genetic screen reveals sequence preferences within a key tertiary interaction in cobalamin riboswitches required for ligand selectivity

Riboswitches are structured mRNA sequences that regulate gene expression by directly binding intracellular metabolites. Generating the appropriate regulatory response requires the RNA rapidly and stably acquire higher-order structure to form the binding pocket, bind the appropriate effector molecule and undergo a structural transition to inform the expression machinery. These requirements place riboswitches under strong kinetic constraints, likely restricting the sequence space accessible by recurrent structural modules such as the kink turn and the T-loop. Class-II cobalamin riboswitches contain two T-loop modules: one directing global folding of the RNA and another buttressing the ligand binding pocket. While the T-loop module directing folding is highly conserved, the T-loop associated with binding is substantially less so, with no clear consensus sequence. To further understand the functional role of the binding-associated module, a functional genetic screen of a library of riboswitches with the T-loop and its interacting nucleotides was used to build an experimental phylogeny comprised of sequences that possess a wide range of cobalamin-dependent regulatory activity. Our results reveal conservation patterns of the T-loop and its interaction with the binding core that allow for rapid tertiary structure formation and demonstrate its importance for generating strong ligand-dependent repression of mRNA expression.

The archaeology of coding RNA

Different theories concerning the origin of RNA (and, in particular, mRNA) point to the concatenation and expansion of proto-tRNA-like structures. Different biochemical and biophysical tools have been used to search for ancient-like RNA elements with a specific structure in genomic viral RNAs, including that of the hepatitis C virus, as well as in cellular mRNA populations, in particular those of human hepatocytes. We define this method as "archaeological," and it has been designed to discover evolutionary patterns through a nonphylogenetic and nonrepresentational strategy. tRNA-like elements were found in structurally or functionally relevant positions both in viral RNA and in one of the liver mRNAs examined, the antagonist interferon-alpha subtype 5 (IFNA5) mRNA. Additionally, tRNA-like elements are highly represented within the hepatic mRNA population, which suggests that they could have participated in the formation of coding RNAs in the distant past. Expanding on this finding, we have observed a recurring dsRNA-like motif next to the tRNA-like elements in both viral RNAs and IFNA5 mRNA. This suggested that the concatenation of these RNA motifs was an activity present in the RNA pools that might have been relevant in the RNA world. The extensive alteration of sequences that likely triggered the transition from the predecessors of coding RNAs to the first fully functional mRNAs (which was not the case in the stepwise construction of noncoding rRNAs) hinders the phylogeny-based identification of RNA elements (both sequences and structures) that might have been active before the advent of protein synthesis. Therefore, our RNA archaeological method is presented as a way to better understand the structural/functional versatility of a variety of RNA elements, which might represent "the losers" in the process of RNA evolution as they had to adapt to the selective pressures favoring the coding capacity of the progressively longer mRNAs.

rRNA Mimicry in RNA Regulation of Gene Expression

The rRNA is the largest and most abundant RNA in bacterial and archaeal cells. It is also one of the best-characterized RNAs in terms of its structural motifs and sequence variation. Production of ribosome components including >50 ribosomal proteins (r-proteins) consumes significant cellular resources. Thus, RNA cis-regulatory structures that interact with r-proteins to repress further r-protein synthesis play an important role in maintaining appropriate stoichiometry between r-proteins and rRNA. Classically, such mRNA structures were thought to directly mimic the rRNA. However, more than 30 years of research has demonstrated that a variety of different recognition and regulatory paradigms are present. This review will demonstrate how structural mimicry between the rRNA and mRNA cis-regulatory structures may take many different forms. The collection of mRNA structures that interact with r-proteins to regulate r-protein operons are best characterized in Escherichia coli, but are increasingly found within species from nearly all phyla of bacteria and several archaea. Furthermore, they represent a unique opportunity to assess the plasticity of RNA structure in the context of RNA-protein interactions. The binding determinants imposed by r-proteins to allow regulation can be fulfilled in many ways. Some r-protein-interacting mRNAs are immediately obvious as rRNA mimics from primary sequence similarity, others are identifiable only after secondary or tertiary structure determination, and some show no obvious similarity. In addition, across different bacterial species a host of different mechanisms of action have been characterized, showing that there is no simple one-size-fits-all solution.

Structural basis for activation of fluorogenic dyes by an RNA aptamer lacking a G-quadruplex motif

The DIR2s RNA aptamer, a second-generation, in-vitro selected binder to dimethylindole red (DIR), activates the fluorescence of cyanine dyes, DIR and oxazole thiazole blue (OTB), allowing detection of two well-resolved emission colors. Using Fab BL3-6 and its cognate hairpin as a crystallization module, we solved the crystal structures of both the apo and OTB-SO₃ bound forms of DIR2s at 2.0 Å and 1.8 Å resolution, respectively. DIR2s adopts a compact, tuning fork-like architecture comprised of a helix and two short stem-loops oriented in parallel to create the ligand binding site through tertiary interactions. The OTB-SO₃ fluorophore binds in a planar conformation to a claw-like structure formed by a purine base-triple, which provides a stacking platform for OTB-SO_3, and an unpaired nucleotide, which partially caps the binding site from the top. The absence of a G-quartet or base tetrad makes the DIR2s aptamer unique among fluorogenic RNAs with known 3D structure.

The DIR2s RNA aptamer activates the fluorescence of cyanine dyes allowing detection of two well-resolved emission colors. Here authors solve the crystal structures of the apo and OTB-SO3 fluorophore-bound DIR2s and show how the fluorophore ligand is bound.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

New tRNA contacts facilitate ligand binding in a Mycobacterium smegmatis T box riboswitch

T box riboswitches are RNA regulatory elements widely used by organisms in the phyla Firmicutes and Actinobacteria to regulate expression of amino acid-related genes. Expression of T box family genes is down-regulated by transcription attenuation or inhibition of translation initiation in response to increased charging of the cognate tRNA. Three direct contacts with tRNA have been described; however, one of these contacts is absent in a subclass of T box RNAs and the roles of several structural domains conserved in most T box RNAs are unknown. In this study, structural elements of a Mycobacterium smegmatis ileS T box riboswitch variant with an Ultrashort (US) Stem I were sequentially deleted, which resulted in a progressive decrease in binding affinity for the tRNA^Ile ligand. Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) revealed structural changes in conserved riboswitch domains upon interaction with the tRNA ligand. Cross-linking and mutational analyses identified two interaction sites, one between the S-turn element in Stem II and the T arm of tRNA^Ile and the other between the Stem IIA/B pseudoknot and the D loop of tRNA^Ile These newly identified RNA contacts add information about tRNA recognition by the T box riboswitch and demonstrate a role for the S-turn and pseudoknot elements, which resemble structural elements that are common in many cellular RNAs.

Structural insights into species-specific features of the ribosome from the human pathogen Mycobacterium tuberculosis

Ribosomes from Mycobacterium tuberculosis (Mtb) possess species-specific ribosomal RNA (rRNA) expansion segments and ribosomal proteins (rProtein). Here, we present the near-atomic structures of the Mtb 50S ribosomal subunit and the complete Mtb 70S ribosome, solved by cryo-electron microscopy. Upon joining of the large and small ribosomal subunits, a 100-nt long expansion segment of the Mtb 23S rRNA, named H54a or the ‘handle’, switches interactions from with rRNA helix H68 and rProtein uL2 to with rProtein bS6, forming a new intersubunit bridge ‘B9’. In Mtb 70S, bridge B9 is mostly maintained, leading to correlated motions among the handle, the L1 stalk and the small subunit in the rotated and non-rotated states. Two new protein densities were discovered near the decoding center and the peptidyl transferase center, respectively. These results provide a structural basis for studying translation in Mtb as well as developing new tuberculosis drugs.

RNA Structure: Advances and Assessment of 3D Structure Prediction

Biological functions of RNA molecules are dependent upon sustained specific three-dimensional (3D) structures of RNA, with or without the help of proteins. Understanding of RNA structure is frequently based on 2D structures, which describe only the Watson-Crick (WC) base pairs. Here, we hierarchically review the structural elements of RNA and how they contribute to RNA 3D structure. We focus our analysis on the non-WC base pairs and on RNA modules. Several computer programs have now been designed to predict RNA modules. We describe the RNA-Puzzles initiative, which is a community-wide, blind assessment of RNA 3D structure prediction programs to determine the capabilities and bottlenecks of current predictions. The assessment metrics used in RNA-Puzzles are briefly described. The detection of RNA 3D modules from sequence data and their automatic implementation belong to the current challenges in RNA 3D structure prediction.

In Silico Prediction of RNA Secondary Structure

The secondary structure of an RNA molecule represents the base-pairing interactions within the molecule and fundamentally determines its overall structure. In this chapter, we overview the main approaches and existing tools for predicting RNA secondary structures, as well as methods for identifying noncoding RNAs from genomic sequences or RNA sequencing data. We then focus on the identification of a well-known class of small noncoding RNAs, namely microRNAs, which play very important roles in many biological processes through regulating post-transcriptionally the expression of genes and which dysregulation has been shown to be involved in several human diseases.

Accumulation of Stable Full-Length Circular Group I Intron RNAs during Heat-Shock

Group I introns in nuclear ribosomal RNA of eukaryotic microorganisms are processed by splicing or circularization. The latter results in formation of full-length circular introns without ligation of the exons and has been proposed to be active in intron mobility. We applied qRT-PCR to estimate the copy number of circular intron RNA from the myxomycete Didymium iridis. In exponentially growing amoebae, the circular introns are nuclear and found in 70 copies per cell. During heat-shock, the circular form is up-regulated to more than 500 copies per cell. The intron harbours two ribozymes that have the potential to linearize the circle. To understand the structural features that maintain circle integrity, we performed chemical and enzymatic probing of the splicing ribozyme combined with molecular modeling to arrive at models of the inactive circular form and its active linear counterpart. We show that the two forms have the same overall structure but differ in key parts, including the catalytic core element P7 and the junctions at which reactions take place. These differences explain the relative stability of the circular species, demonstrate how it is prone to react with a target molecule for circle integration and thus supports the notion that the circular form is a biologically significant molecule possibly with a role in intron mobility.

The role of RNA conformation in RNA-protein recognition

Interactions between protein and RNA play a key role in many biological processes in the gene expression pathway. Those interactions are mediated through a variety of RNA-binding protein domains, among them the highly abundant RNA recognition motif (RRM). Here we studied protein-RNA complexes from different RNA binding domain families solved by NMR and x-ray crystallography. Characterizing the structural properties of the RNA at the binding interfaces revealed an unexpected number of nucleotides with unusual RNA conformations, specifically found in RNA-RRM complexes. Moreover, we observed that the RNA nucleotides that are directly involved in interactions with the RRM domains, via hydrogen bonds and hydrophobic contacts, are significantly enriched with unique RNA conformations. Further examination of the sequences binding the RRM domain showed a preference for G nucleotides in syn conformation to precede or to follow U nucleotides in the anti-conformation, and U nucleotides in C2' endo conformation to precede U and G nucleotides possessing the more common C3' endo conformation. These findings imply a possible mode of RNA recognition by the RRM domains which enables the recognition of a wide variety of different RNA sequences and shapes. Overall, this study suggests an additional way by which the RRM domain recognizes its RNA target, involving a conformational readout.

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.

It's a loop world - single strands in RNA as structural and functional elements

Unpaired regions in RNA molecules - loops - are centrally involved in defining the characteristic three-dimensional (3D) architecture of RNAs and are of high interest in RNA engineering and design. Loops adopt diverse, but specific conformations stabilised by complex tertiary structural interactions that provide structural flexibility to RNA structures that would otherwise not be possible if they only consisted of the rigid A-helical shapes usually formed by canonical base pairing. By participating in sequence-non-local contacts, they furthermore contribute to stabilising the overall fold of RNA molecules. Interactions between RNAs and other nucleic acids, proteins, or small molecules are also generally mediated by RNA loop structures. Therefore, the function of an RNA molecule is generally dependent on its loops. Examples include intermolecular interactions between RNAs as part of the microRNA processing pathways, ribozymatic activity, or riboswitch-ligand interactions. Bioinformatics approaches have been successfully applied to the identification of novel RNA structural motifs including loops, local and global RNA 3D structure prediction, and structural and conformational analysis of RNAs and have contributed to a better understanding of the sequence-structure-function relationships in RNA loops.

An introduction to recurrent nucleotide interactions in RNA

RNA secondary structure diagrams familiar to molecular biologists summarize at a glance the folding of RNA chains to form Watson–Crick paired double helices. However, they can be misleading: First of all, they imply that the nucleotides in loops and linker segments, which can amount to 35% to 50% of a structured RNA, do not significantly interact with other nucleotides. Secondly, they give the impression that RNA molecules are loosely organized in three-dimensional (3D) space. In fact, structured RNAs are compactly folded as a result of numerous long-range, sequence-specific interactions, many of which involve loop or linker nucleotides. Here, we provide an introduction for students and researchers of RNA on the types, prevalence, and sequence variations of inter-nucleotide interactions that structure and stabilize RNA 3D motifs and architectures, using Escherichia coli (E. coli) 16S ribosomal RNA as a concrete example. The picture that emerges is that almost all nucleotides in structured RNA molecules, including those in nominally single-stranded loop or linker regions, form specific interactions that stabilize functional structures or mediate interactions with other molecules. The small number of noninteracting, ‘looped-out’ nucleotides make it possible for the RNA chain to form sharp turns. Base-pairing is the most specific interaction in RNA as it involves edge-to-edge hydrogen bonding (H-bonding) of the bases. Non-Watson–Crick base pairs are a significant fraction (30% or more) of base pairs in structured RNAs.

The ancient heart of the ribosomal large subunit: a response to Caetano-Anolles

Our recent Accretion Model of ribosomal evolution uses insertion fingerprints and a "trunk-branch" formalism to recapitulate the building up of common core rRNA of the Large Ribosomal Subunit. The Accretion Model is a conservative and natural extension of a method developed by Bokov and Steinberg (Nature 457:977-80, 2009), which confirms the correctness of lower resolution models by Fox and others. In each of these models, the LSU originates with the peptidyl transferase center (PTC), consistent with expectations that the ribosome is the source of defined-sequence functional proteins. In an adjacent note, Caetano-Anolles (J Mol Evol 80:162-165, 2015) disparages the Accretion Model, because it controverts the 'Growth Inferred by Genothermal Ordering' (GIGO) model. GIGO analyzes secondary structures, assigns the origin of the ribosome to a region outside of the PTC, and assumes or deduces that (i) large protein enzymes of defined amino acid sequence predate ribosomal synthesis of proteins, (ii) proteins directly replicate by non-ribosomal mechanisms, (iii) rRNA unfailingly increases in thermodynamic stability over time, and (iv) the Woese and Fox canonical tree of life is mis-rooted. Much of the specific GIGO critique of the Accretion Model is based on confusion about the three-dimensional nature of RNA and trunk-branch polymorphism; the Accretion Model incorporates several types of trunk-branch relationships.

Functional conservation despite structural divergence in ligand-responsive RNA switches

An internal ribosome entry site (IRES) initiates protein synthesis in RNA viruses, including the hepatitis C virus (HCV). We have discovered ligand-responsive conformational switches in viral IRES elements. Modular RNA motifs of greatly distinct sequence and local secondary structure have been found to serve as functionally conserved switches involved in viral IRES-driven translation and may be captured by identical cognate ligands. The RNA motifs described here constitute a new paradigm for ligand-captured switches that differ from metabolite-sensing riboswitches with regard to their small size, as well as the intrinsic stability and structural definition of the constitutive conformational states. These viral RNA modules represent the simplest form of ligand-responsive mechanical switches in nucleic acids.

Principles of ion recognition in RNA: insights from the group II intron structures

Metal ions promote both RNA folding and catalysis, thus being essential in stabilizing the structure and determining the function of large RNA molecules, including group II introns. The latter are self-splicing metalloribozymes, containing a heteronuclear four-metal-ion center within the active site. In addition to these catalytic ions, group II introns bind many other structural ions, including delocalized ions that bind the RNA diffusively and well-ordered ions that bind the RNA tightly with high occupancy. The latter ions, which can be studied by biophysical methods, have not yet been analyzed systematically. Here, we compare crystal structures of the group IIC intron from Oceanobacillus iheyensis and classify numerous site-bound ions, which are primarily localized in the intron core and near long-range tertiary contacts. Certain ion-binding sites resemble motifs observed in known RNA structures, while others are idiosyncratic to the group II intron. Particularly interesting are (1) ions proximal to the active site, which may participate in splicing together with the catalytic four-metal-ion center, (2) organic ions that bind regions predicted to interact with intron-encoded proteins, and (3) unusual monovalent ions bound to GU wobble pairs, GA mismatches, the S-turn, the tetraloop-receptor, and the T-loop. Our analysis extends the general principles by which ions participate in RNA structural organization and it will aid in the determination and interpretation of future RNA structures.

Secondary structures of rRNAs from all three domains of life

Accurate secondary structures are important for understanding ribosomes, which are extremely large and highly complex. Using 3D structures of ribosomes as input, we have revised and corrected traditional secondary (2°) structures of rRNAs. We identify helices by specific geometric and molecular interaction criteria, not by co-variation. The structural approach allows us to incorporate non-canonical base pairs on parity with Watson-Crick base pairs. The resulting rRNA 2° structures are up-to-date and consistent with three-dimensional structures, and are information-rich. These 2° structures are relatively simple to understand and are amenable to reproduction and modification by end-users. The 2° structures made available here broadly sample the phylogenetic tree and are mapped with a variety of data related to molecular interactions and geometry, phylogeny and evolution. We have generated 2° structures for both large subunit (LSU) 23S/28S and small subunit (SSU) 16S/18S rRNAs of Escherichia coli, Thermus thermophilus, Haloarcula marismortui (LSU rRNA only), Saccharomyces cerevisiae, Drosophila melanogaster, and Homo sapiens. We provide high-resolution editable versions of the 2° structures in several file formats. For the SSU rRNA, the 2° structures use an intuitive representation of the central pseudoknot where base triples are presented as pairs of base pairs. Both LSU and SSU secondary maps are available (http://apollo.chemistry.gatech.edu/RibosomeGallery). Mapping of data onto 2° structures was performed on the RiboVision server (http://apollo.chemistry.gatech.edu/RiboVision).

Isosteric and nonisosteric base pairs in RNA motifs: molecular dynamics and bioinformatics study of the sarcin-ricin internal loop

The sarcin-ricin RNA motif (SR motif) is one of the most prominent recurrent RNA building blocks that occurs in many different RNA contexts and folds autonomously, that is, in a context-independent manner. In this study, we combined bioinformatics analysis with explicit-solvent molecular dynamics (MD) simulations to better understand the relation between the RNA sequence and the evolutionary patterns of the SR motif. A SHAPE probing experiment was also performed to confirm the fidelity of the MD simulations. We identified 57 instances of the SR motif in a nonredundant subset of the RNA X-ray structure database and analyzed their base pairing, base-phosphate, and backbone-backbone interactions. We extracted sequences aligned to these instances from large rRNA alignments to determine the frequency of occurrence for different sequence variants. We then used a simple scoring scheme based on isostericity to suggest 10 sequence variants with a highly variable expected degree of compatibility with the SR motif 3D structure. We carried out MD simulations of SR motifs with these base substitutions. Nonisosteric base substitutions led to unstable structures, but so did isosteric substitutions which were unable to make key base-phosphate interactions. The MD technique explains why some potentially isosteric SR motifs are not realized during evolution. We also found that the inability to form stable cWW geometry is an important factor in the case of the first base pair of the flexible region of the SR motif. A comparison of structural, bioinformatics, SHAPE probing, and MD simulation data reveals that explicit solvent MD simulations neatly reflect the viability of different sequence variants of the SR motif. Thus, MD simulations can efficiently complement bioinformatics tools in studies of conservation patterns of RNA motifs and provide atomistic insight into the role of their different signature interactions.

Secondary structure and domain architecture of the 23S and 5S rRNAs

We present a de novo re-determination of the secondary (2°) structure and domain architecture of the 23S and 5S rRNAs, using 3D structures, determined by X-ray diffraction, as input. In the traditional 2° structure, the center of the 23S rRNA is an extended single strand, which in 3D is seen to be compact and double helical. Accurately assigning nucleotides to helices compels a revision of the 23S rRNA 2° structure. Unlike the traditional 2° structure, the revised 2° structure of the 23S rRNA shows architectural similarity with the 16S rRNA. The revised 2° structure also reveals a clear relationship with the 3D structure and is generalizable to rRNAs of other species from all three domains of life. The 2° structure revision required us to reconsider the domain architecture. We partitioned the 23S rRNA into domains through analysis of molecular interactions, calculations of 2D folding propensities and compactness. The best domain model for the 23S rRNA contains seven domains, not six as previously ascribed. Domain 0 forms the core of the 23S rRNA, to which the other six domains are rooted. Editable 2° structures mapped with various data are provided (http://apollo.chemistry.gatech.edu/RibosomeGallery).

Conformational readout of RNA by small ligands

RNA molecules have highly versatile structures that can fold into myriad conformations, providing many potential pockets for binding small molecules. The increasing number of available RNA structures, in complex with proteins, small ligands and in free form, enables the design of new therapeutically useful RNA-binding ligands. Here we studied RNA ligand complexes from 10 RNA groups extracted from the protein data bank (PDB), including adaptive and non-adaptive complexes. We analyzed the chemical, physical, structural and conformational properties of binding pockets around the ligand. Comparing the properties of ligand-binding pockets to the properties of computed pockets extracted from all available RNA structures and RNA-protein interfaces, revealed that ligand-binding pockets, mainly the adaptive pockets, are characterized by unique properties, specifically enriched in rare conformations of the nucleobase and the sugar pucker. Further, we demonstrate that nucleotides possessing the rare conformations are preferentially involved in direct interactions with the ligand. Overall, based on our comprehensive analysis of RNA-ligand complexes, we suggest that the unique conformations adopted by RNA nucleotides play an important role in RNA recognition by small ligands. We term the recognition of a binding site by a ligand via the unique RNA conformations "RNA conformational readout." We propose that "conformational readout" is a general way by which RNA binding pockets are recognized and selected from an ensemble of different RNA states.

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.

RNA secondary structure prediction using soft computing

Prediction of RNA structure is invaluable in creating new drugs and understanding genetic diseases. Several deterministic algorithms and soft computing-based techniques have been developed for more than a decade to determine the structure from a known RNA sequence. Soft computing gained importance with the need to get approximate solutions for RNA sequences by considering the issues related with kinetic effects, cotranscriptional folding, and estimation of certain energy parameters. A brief description of some of the soft computing-based techniques, developed for RNA secondary structure prediction, is presented along with their relevance. The basic concepts of RNA and its different structural elements like helix, bulge, hairpin loop, internal loop, and multiloop are described. These are followed by different methodologies, employing genetic algorithms, artificial neural networks, and fuzzy logic. The role of various metaheuristics, like simulated annealing, particle swarm optimization, ant colony optimization, and tabu search is also discussed. A relative comparison among different techniques, in predicting 12 known RNA secondary structures, is presented, as an example. Future challenging issues are then mentioned.

The DNA and RNA sugar-phosphate backbone emerges as the key player. An overview of quantum-chemical, structural biology and simulation studies

Knowledge of geometrical and physico-chemical properties of the sugar-phosphate backbone substantially contributes to the comprehension of the structural dynamics, function and evolution of nucleic acids. We provide a side by side overview of structural biology/bioinformatics, quantum chemical and molecular mechanical/simulation studies of the nucleic acids backbone. We highlight main features, advantages and limitations of these techniques, with a special emphasis given to their synergy. The present status of the research is then illustrated by selected examples which include classification of DNA and RNA backbone families, benchmark structure-energy quantum chemical calculations, parameterization of the dihedral space of simulation force fields, incorporation of arsenate into DNA, sugar-phosphate backbone self-cleavage in small RNA enzymes, and intricate geometries of the backbone in recurrent RNA building blocks. Although not apparent from the current literature showing limited overlaps between the QM, simulation and bioinformatics studies of the nucleic acids backbone, there in fact should be a major cooperative interaction between these three approaches in studies of the sugar-phosphate backbone.

Structural architecture of an RNA that competitively inhibits RNase L

Activation of RNase L endonuclease activity is part of the mammalian innate immune response to viral infection. The poliovirus RNA genome contains a sequence in its protein-coding region that can act as a competitive inhibitor of RNase L. Mutation, sequence, and functional analysis of this competitive inhibitor RNA (ciRNA) revealed that its activity depends on specific sequences, showed that a loop-loop hairpin interaction forms in the ciRNA, and suggested the presence of a loop E motif. These features lead to the hypothesis that the ciRNA's function is conferred in part by a specific three-dimensional folded RNA architecture. By using a combination of biophysical, mutational, and functional studies, we have mapped features of the three-dimensional architecture of the ciRNA in its unbound form. We show that the loop-loop interaction forms in the free ciRNA and affects the overall structure, perhaps forming long-range tertiary interactions with the loop E motif. Local tight RNA-RNA backbone packing occurs in parts of the structure, but the fold appears to be less stable than many other tightly packed RNAs. This feature may allow the ciRNA to accommodate the translocation of ribosomes and polymerase across this multifunctional region of the viral RNA but also to function as an RNase L inhibitor.

DISCOVERING COMMON STRUCTURAL MOTIFS OF RIBOSOMAL RNA SECONDARY STRUCTURES IN PROKARYOTES

Certain structural motifs, like tetra-loops, in ribosomal RNA are known to functionally implicate in virtually every aspect of protein synthesis. Ribosomal RNA molecules were also widely used as a tool in molecular evolutionary studies because of their ubiquity, size and low evolutionary rate. In this study, we adapt a data mining approach to discover common structural motifs, and then we use a machine learning approach to identify discriminating CSMs from groups of organisms. Finally, we construct phylogeneitc trees to investigate the evolution of ribosomal RNA by serving the CSMs discovered as targets, which are used to estimate the evolutionary relatedness between organisms. The aim of this study is to discover common structural motifs (CSMs), i.e., those single-strain regions shared in ribosomal RNA secondary structures by several organisms, which are related to specific domains or functions. We discover a set of common structural motifs from several data sets of Archaea and Bacteria. Significant CSMs are then induced by a decision tree. Furthermore, phylogenetic trees are constructed based on CSMs and primary sequences of SSU 16 S ribosomal RNA.

Comprehensive survey and geometric classification of base triples in RNA structures

Base triples are recurrent clusters of three RNA nucleobases interacting edge-to-edge by hydrogen bonding. We find that the central base in almost all triples forms base pairs with the other two bases of the triple, providing a natural way to geometrically classify base triples. Given 12 geometric base pair families defined by the Leontis-Westhof nomenclature, combinatoric enumeration predicts 108 potential geometric base triple families. We searched representative atomic-resolution RNA 3D structures and found instances of 68 of the 108 predicted base triple families. Model building suggests that some of the remaining 40 families may be unlikely to form for steric reasons. We developed an on-line resource that provides exemplars of all base triples observed in the structure database and models for unobserved, predicted triples, grouped by triple family, as well as by three-base combination (http://rna.bgsu.edu/Triples). The classification helps to identify recurrent triple motifs that can substitute for each other while conserving RNA 3D structure, with applications in RNA 3D structure prediction and analysis of RNA sequence evolution.

RNA self-cleavage activated by ultraviolet light-induced oxidation

A novel UV-C-light-induced ribozyme activity was discovered within the highly structured 5'-genomic regions of both Hepatitis C Virus (HCV) and the related Classic Swine Fever Virus (CSFV). Cleavage is mediated by exposure to UV-C light but not by exogenous oxygen radicals. It is also very selective, occurring at base positions HCV C(79) and CSFV A(45) in some molecules and at the immediately adjacent 5'-positions HCV U(78) and CSFV U(44) in others. Among other reaction products, the majority of biochemically active products detected contained 3'-phosphate and 5'-phosphate-end groups at the newly generated termini, along with a much lower amount of 3'-hydroxyl end group. While preservation of an E-loop RNA structure in the vicinity of the cleavage site was a requisite for HCV RNA self-cleavage, this was not the case for CSFV RNA. The short size of the reactive domains (~33 nt), which are compatible with primitive RNA motifs, and the lack of sequence homology, indicate that as-yet unidentified UV-activated ribozymes are likely to be found throughout structured RNAs, thereby providing clues to whether early RNA self-cleavage events were mediated by photosensitive RNA structures.

Sequence-based identification of 3D structural modules in RNA with RMDetect

Structural RNA modules, sets of ordered non-Watson-Crick base pairs embedded between Watson-Crick pairs, have central roles as architectural organizers and sites of ligand binding in RNA molecules, and are recurrently observed in RNA families throughout the phylogeny. Here we describe a computational tool, RNA three-dimensional (3D) modules detection, or RMDetect, for identifying known 3D structural modules in single and multiple RNA sequences in the absence of any other information. Currently, four modules can be searched for: G-bulge loop, kink-turn, C-loop and tandem-GA loop. In control test sequences we found all of the known modules with a false discovery rate of 0.23. Scanning through 1,444 publicly available alignments, we identified 21 yet unreported modules and 141 known modules. RMDetect can be used to refine RNA 2D structure, assemble RNA 3D models, and search and annotate structured RNAs in genomic data.

Solution structure of an alternate conformation of helix27 from Escherichia coli16S rRNA

Helix (H)27 of 16S ribosomal (r)RNA from Escherichia coli was dubbed the "switch helix" when mutagenesis suggested that two alternative base pair registers may have distinct functional roles in the bacterial ribosome. Although more recent genetic analyses suggest that H27 conformational switching is not required for translation, previous solution studies demonstrated that the isolated E. coli H27 can dynamically convert between the 885 and 888 conformations. Here, we have solved the nuclear magnetic resonance solution structure of a locked 888 conformation. NOE and residual dipolar coupling restraints reveal an architecture that markedly differs from that of the 885 conformation found in crystal structures of the bacterial ribosome. In place of the loop E motif that characterizes the 885 conformer and that the 888 conformer cannot adopt, we find evidence for an asymmetrical A-rich internal loop stabilized by stacking interactions among the unpaired A's. Comparison of the isolated H27 888 solution structure with the 885 crystal structure within the context of the ribosome suggests a difference in overall length of H27 that presents one plausible reason for the absence of H27 conformational switching within the sterically confining ribosome.

Solution structure of the K-turn and Specifier Loop domains from the Bacillus subtilis tyrS T-box leader RNA

In Gram-positive bacteria, the RNA transcripts of many amino acid biosynthetic and aminoacyl tRNA synthetase genes contain 5' untranslated regions, or leader RNAs, that function as riboswitches. These T-box riboswitches bind cognate tRNA molecules and regulate gene expression by a transcription attenuation mechanism. The Specifier Loop domain of the leader RNA contains nucleotides that pair with nucleotides in the tRNA anticodon loop and is flanked on one side by a kink-turn (K-turn), or GA, sequence motif. We have determined the solution NMR structure of the K-turn sequence element within the context of the Specifier Loop domain. The K-turn sequence motif has several noncanonical base pairs typical of K-turn structures but adopts an extended conformation. The Specifier Loop domain contains a loop E structural motif, and the single-strand Specifier nucleotides stack with their Watson-Crick edges displaced toward the minor groove. Mg(2+) leads to a significant bending of the helix axis at the base of the Specifier Loop domain, but does not alter the K-turn. Isothermal titration calorimetry indicates that the K-turn sequence causes a small enhancement of the interaction between the tRNA anticodon arm and the Specifier Loop domain. One possibility is that the K-turn structure is formed and stabilized when tRNA binds the T-box riboswitch and interacts with Stem I and the antiterminator helix. This motif in turn anchors the orientation of Stem I relative to the 3' half of the leader RNA, further stabilizing the tRNA-T box complex.

Predicting and modeling RNA architecture

A general approach for modeling the architecture of large and structured RNA molecules is described. The method exploits the modularity and the hierarchical folding of RNA architecture that is viewed as the assembly of preformed double-stranded helices defined by Watson-Crick base pairs and RNA modules maintained by non-Watson-Crick base pairs. Despite the extensive molecular neutrality observed in RNA structures, specificity in RNA folding is achieved through global constraints like lengths of helices, coaxiality of helical stacks, and structures adopted at the junctions of helices. The Assemble integrated suite of computer tools allows for sequence and structure analysis as well as interactive modeling by homology or ab initio assembly with possibilities for fitting within electronic density maps. The local key role of non-Watson-Crick pairs guides RNA architecture formation and offers metrics for assessing the accuracy of three-dimensional models in a more useful way than usual root mean square deviation (RMSD) values.

NMR structure of the A730 loop of the Neurospora VS ribozyme: insights into the formation of the active site

The Neurospora VS ribozyme is a small nucleolytic ribozyme with unique primary, secondary and global tertiary structures, which displays mechanistic similarities to the hairpin ribozyme. Here, we determined the high-resolution NMR structure of a stem–loop VI fragment containing the A730 internal loop, which forms part of the active site. In the presence of magnesium ions, the A730 loop adopts a structure that is consistent with existing biochemical data and most likely reflects its conformation in the VS ribozyme prior to docking with the cleavage site internal loop. Interestingly, the A730 loop adopts an S-turn motif that is also present in loop B within the hairpin ribozyme active site. The S-turn appears necessary to expose the Watson–Crick edge of a catalytically important residue (A756) so that it can fulfill its role in catalysis. The A730 loop and the cleavage site loop of the VS ribozyme display structural similarities to internal loops found in the active site of the hairpin ribozyme. These similarities provided a rationale to build a model of the VS ribozyme active site based on the crystal structure of the hairpin ribozyme.