Computational discovery of folded RNA domains in genomes and in vitro selected libraries

Identification of Hammerhead-variant ribozyme sequences in SARS-CoV-2

The SARS-CoV-2 RNA virus and variants, responsible for the COVID-19 pandemic has become endemic, raised a need for further understanding of the viral genome and biology. Despite vast research on SARS-CoV-2, no ribozymes have been found in the virus genome. Here we report the identification of 39 Hammerhead-variant ribozyme sequences (CoV-HHRz) in SARS-CoV-2. These sequences are highly conserved within SARS-CoV-2 variants but show large diversity among other coronaviruses. In vitro CoV-HHRz sequences possess the characteristics of typical ribozymes; cleavage is pH and ion dependent, although their activity is relatively low and Mn2+ is required for cleavage. The cleavage sites of four CoV-HHRz coincide with the breakpoint of expressed subgenomic RNA (sgRNAs) in SARS-CoV-2 transcriptome data suggesting in vivo activity. The CoV-HHRz are involved in processing sgRNAs for ORF7b, ORF 10 and ORF1ab nsp13 which are essential for viral packaging and life cycle.

An efficient and universal In silico screening strategy for acquisition of high-affinity Aptamer and its application in analytical utility

Numerous aptamers against various targets have been identified through the technology of systematic evolution of ligands by exponential enrichment (SELEX), but the affinity of these aptamers are often insufficient due to the limitations of SELEX. Therefore, a more rational in silico screening strategy (ISS) was developed for efficient screening of high affinity aptamers, which took shape complementarity and thermodynamic stability into consideration. Neuron specific enolase (NSE), a tumor marker, was selected as the target molecule. In the screening process, three aptamer candidates with good shape complementarity, lower ΔG values, and higher ZDOCK scores were produced. The dissociation constant (Kd) of these candidates to NSE was determined to be 10.13 nM, 14.82 nM, and 2.76 nM, respectively. Each of them exhibited higher affinity to NSE than the parent aptamer (Kd = 23.83 nM). Finally, an antibody-free fluorescence aptasensor assay, based on the aptamer with the highest affinity, P-5C8G, was conducted, resulting in a limit of detection (LOD) value of 1.8 nM, which was much lower than the parental aptamer (P, LOD = 12.6 nM). The proposed ISS approach provided an efficient and universal strategy to improve the aptamer to have a high affinity and good analytical utility.

Identification of HDV-like theta ribozymes involved in tRNA-based recoding of gut bacteriophages

Trillions of microorganisms, collectively known as the microbiome, inhabit our bodies with the gut microbiome being of particular interest in biomedical research. Bacteriophages, the dominant virome constituents, can utilize suppressor tRNAs to switch to alternative genetic codes (e.g., the UAG stop-codon is reassigned to glutamine) while infecting hosts with the standard bacterial code. However, what triggers this switch and how the bacteriophage manipulates its host is poorly understood. Here, we report the discovery of a subgroup of minimal hepatitis delta virus (HDV)-like ribozymes - theta ribozymes - potentially involved in the code switch leading to the expression of recoded lysis and structural phage genes. We demonstrate their HDV-like self-scission behavior in vitro and find them in an unreported context often located with their cleavage site adjacent to tRNAs, indicating a role in viral tRNA maturation and/or regulation. Every fifth associated tRNA is a suppressor tRNA, further strengthening our hypothesis. The vast abundance of tRNA-associated theta ribozymes - we provide 1753 unique examples - highlights the importance of small ribozymes as an alternative to large enzymes that usually process tRNA 3'-ends. Our discovery expands the short list of biological functions of small HDV-like ribozymes and introduces a previously unknown player likely involved in the code switch of certain recoded gut bacteriophages.

Molecular modeling and simulation studies of SELEX-derived high-affinity DNA aptamers to the Ebola virus nucleoprotein

Ebola viral disease (EVD) is a highly infectious and potentially fatal illness with a case fatality rate ranging from 25% to 90%. To effectively control its spread, there is a need for rapid, reliable and lowcost point-of-care (P OC) diagnostic tests. While various EVD diagnostic tests exist, few are P OC tests, and many are not cost-effective. The use of antibodies in these tests has limitations, prompting the exploration of aptamers as potential alternatives. Various proteins from the Ebola virus (EBOV) proteome, including EBOV nucleoprotein (NP), are considered viable targets for diagnostic assays. A previous study identified three aptamers (Apt1. Apt2 and Apt3) with high affinity for EBOV NP using systemic evolution of ligands by exponential enrichment (SELEX). This study aimed to employ in silico methods, such as Phyre2, RNAfold, RNAComposer, HADDOCK and GROMACS, to model the structures of EBOV NP and the aptamers, and to investigate their binding. The in silico analysis revealed successful binding of all the three aptamers to EBOV NP, with a suggested ranking of Apt1 > Apt2 > Apt3 based on binding affinity. Microscale thermophoresis (MST) analysis confirmed the binding, providing dissociation constants of 25 ± 2.84, 56 ± 2.76 and 140 ±3.69 nM for Apt1, Apt2 and Apt3, respectively. The study shows that the findings of the in silico analysis was in agreement with the MST analysis. Inclusion of these in silico approaches in diagnostic assay development can expedite the selection of candidate aptamers, potentially overcoming challenges associated with aptamer application in diagnostics.Communicated by Ramaswamy H. Sarma.

Aptamer-Based Technologies for Parasite Detection

Centuries of scientific breakthroughs have brought us closer to understanding and managing the spread of parasitic diseases. Despite ongoing technological advancements in the detection, treatment, and control of parasitic illnesses, their effects on animal and human health remain a major concern worldwide. Aptamers are single-stranded oligonucleotides whose unique three-dimensional structures enable them to interact with high specificity and affinity to a wide range of targets. In recent decades, aptamers have emerged as attractive alternatives to antibodies as therapeutic and diagnostic agents. Due to their superior stability, reusability, and modifiability, aptamers have proven to be effective bioreceptors for the detection of toxins, contaminants, biomarkers, whole cells, pathogens, and others. As such, they have been integrated into a variety of electrochemical, fluorescence, and optical biosensors to effectively detect whole parasites and their proteins. This review offers a summary of the various types of parasite-specific aptamer-based biosensors, their general mechanisms and their performance.

Review on the Selection of Aptamers and Application in Paper-Based Sensors

An aptamer is a synthetic oligonucleotide, referring to a single-stranded deoxyribonucleic acid or ribonucleic acid ligand produced by synthesis from outside the body using systematic evolution of ligands by exponential enrichment (SELEX) technology. Owing to their special screening process and adjustable tertiary structures, aptamers can bind to multiple targets (small molecules, proteins, and even whole cells) with high specificity and affinity. Moreover, due to their simple preparation and stable modification, they have been widely used to construct biosensors for target detection. The paper-based sensor is a product with a low price, short detection time, simple operation, and other superior characteristics, and is widely used as a rapid detection method. This review mainly focuses on the screening methods of aptamers, paper-based devices, and applicable sensing strategies. Furthermore, the design of the aptamer-based lateral flow assay (LFA), which underlies the most promising devices for commercialization, is emphasized. In addition, the development prospects and potential applications of paper-based biosensors using aptamers as recognition molecules are also discussed.

A Novel ssDNA Aptamer Targeting Carcinoembryonic Antigen: Selection and Characterization

One of the major causes of a drastically shorter life expectancy and one of the most prevalent diseases in the world today is cancer. Given the data on the rise in cancer cases throughout the world, it is obvious that, despite the diagnostic techniques currently being used, there is a pressing need to develop precise and sensitive techniques for early diagnosis of the disease. A high degree of affinity and specificity towards particular targets is maintained by the short nucleic acid molecules known as aptamers. Aptamers outperform antibodies due to their unique benefits, such as their simplicity in synthesis and modification, lack of toxicity, and long-term stability. Utilizing an accurate recognition element and a robust signal transduction mechanism, molecular diagnostics can be extremely sensitive and specific. In this study, development of new single-stranded DNA aptamers against CEA for use in cancer diagnostics was accomplished using SELEX and NGS methods. As a result of 12 iterative SELEX rounds, nine aptamer candidates against CEA were developed. NGS comparative analysis revealed that round twelve had an enriched number of aptamers that were specifically bound, as opposed to round eight. Among the selected nine sequences characterized by bioinformatics analysis and ELONA, an aptamer sequence with the highest specificity and affinity for the target protein was identified and further examined. Aptamer sequence (6) was screened in a concentration-dependent assay, specificity analysis was performed, and its potential secondary and tertiary structures were predicted, which enabled us to test one of the possible putative interactions with CEA. Finally, aptamer sequence (6) labelled with a Cy5 fluorescent tag was used in confocal microscopy to observe its binding towards the CEA expressed in HT-29 human colon adenocarcinoma cell line.

The function of twister ribozyme variants in non-LTR retrotransposition in Schistosoma mansoni

The twister ribozyme is widely distributed over numerous organisms and is especially abundant in Schistosoma mansoni, but has no confirmed biological function. Of the 17 non-LTR retrotransposons known in S. mansoni, none have thus far been associated with ribozymes. Here we report the identification of novel twister variant (T-variant) ribozymes and their function in S. mansoni non-LTR retrotransposition. We show that T-variant ribozymes are located at the 5′ end of Perere-3 non-LTR retrotransposons in the S. mansoni genome. T-variant ribozymes were demonstrated to be catalytically active in vitro. In reporter constructs, T-variants were shown to cleave in vivo, and cleavage of T-variants was sufficient for the translation of downstream reporter genes. Our analysis shows that the T-variants and Perere-3 are transcribed together. Target site duplications (TSDs); markers of target-primed reverse transcription (TPRT) and footmarks of retrotransposition, are located adjacent to the T-variant cleavage site and suggest that T-variant cleavage has taken place inS. mansoni. Sequence heterogeneity in the TSDs indicates that Perere-3 retrotransposition is not site-specific. The TSD sequences contribute to the 5′ end of the terminal ribozyme helix (P1 stem). Based on these results we conclude that T-variants have a functional role in Perere-3 retrotransposition.

A database of flavivirus RNA structures with a search algorithm for pseudoknots and triple base interactions

Motivation

The Flavivirus genus includes several important pathogens, such as Zika, dengue and yellow fever virus. Flavivirus RNA genomes contain a number of functionally important structures in their 3′ untranslated regions (3′UTRs). Due to the diversity of sequences and topologies of these structures, their identification is often difficult. In contrast, predictions of such structures are important for understanding of flavivirus replication cycles and development of antiviral strategies.

Results

We have developed an algorithm for structured pattern search in RNA sequences, including secondary structures, pseudoknots and triple base interactions. Using the data on known conserved flavivirus 3′UTR structures, we constructed structural descriptors which covered the diversity of patterns in these motifs. The descriptors and the search algorithm were used for the construction of a database of flavivirus 3′UTR structures. Validating this approach, we identified a number of domains matching a general pattern of exoribonuclease Xrn1-resistant RNAs in the growing group of insect-specific flaviviruses.

Availability and implementation

The Leiden Flavivirus RNA Structure Database is available at https://rna.liacs.nl. The search algorithm is available at https://github.com/LeidenRNA/SRHS.

Supplementary information

Supplementary data are available at Bioinformatics online.

In silico molecular docking in DNA aptamer development

Aptamers are single-stranded DNA or RNA oligonucleotides generated by SELEX that exhibit binding affinity and specificity against a wide variety of target molecules. Compared to RNA aptamers, DNA aptamers are much more stable and therefore are widely adopted in a number of applications especially in diagnostics. The tediousness and rigor associated with certain steps of the SELEX intensify the efforts to adopt in silico molecular docking approaches together with in vitro SELEX procedures in developing DNA aptamers. Inspired by these endeavors, we carry out an overview of the in silico molecular docking approaches in DNA aptamer generation, by detailing the stepwise procedures as well as shedding some light on the various softwares used. The in silico maturation strategy and the limitations of the in silico approaches are also underscored.

Perspective on the Future Role of Aptamers in Analytical Chemistry

It has been almost 30 years since the invention of Systematic Evolution of Ligands by Exponential Enrichment (SELEX) methodology and the description of the first aptamers. In retrospect over the past 30 years, advances in aptamer development and application have demonstrated that aptamers are potentially useful reagents that can be employed in diverse areas within analytical chemistry, biotechnology, biomedicine, and molecular biology. While often touted as artificial antibodies with an ability to be selected for any target, aptamer development, unfortunately, lags behind development of analytical methodologies that employ aptamers, hindering deeper integration into the application of analytical tool development. This perspective covers recent advances in SELEX methodology for improving efficiency of the SELEX procedure and enhancing affinity and specificity of the selected aptamers, what we view as a critical barrier in the future role of aptamers in analytical chemistry. We discuss postselection modifications that can be used for enhancing performance of the selected aptamers in an analytical device by including understanding intermolecular interaction forces in the binding domain. While highlighting promising properties of aptamers that enable several analytical advances, we provide discussion on the challenges of penetration of aptamers in the analytical field.

Investigations on the interface of nucleic acid aptamers and binding targets

Nucleic acid aptamers are single-stranded DNA or RNA of 20-100 nucleotides in length that have attracted substantial scientific interest due to their ability to specifically bind to target molecules via the formation of three-dimensional structures. Compared to traditional protein antibodies, aptamers have several advantages, such as their small size, high binding affinity, specificity, flexible structure, being chemical synthesizable and modifiable, good biocompatibility, high stability and low immunogenicity, which all contribute to their widely applications in the biomedical field. To date, much progress has been made in the study and applications of aptamers, however, detailed information on how aptamers bind to their targets is still scarce. Over the past few decades, many methods have been introduced to investigate the aptamer-target binding process, such as measuring the main kinetic or thermodynamic parameters, detecting the structural changes of the binding complexes, etc. Apart from traditional physicochemical methods, various types of molecular docking programs have been applied to simulate the aptamer-target interactions, while these simulations also have limitations. To facilitate the further research on the interactions, herein, we provide a brief review to illustrate the recent advances in the study of aptamer-target interactions. We summarize the binding targets of aptamers, such as small molecules, macromolecules, and even cells. Their binding constants (KD) are also summarized. Methods to probe the aptamer-target binding process, such as surface plasmon resonance (SPR), circular dichroism spectroscopy (CD), isothermal titration calorimetry (ITC), footprinting assay, truncation and mutation assay, nuclear magnetic resonance spectroscopy (NMR), X-ray crystallography and molecular docking simulation are indicated. The binding forces mediating the aptamer-target interactions, such as hydrogen bonding, electrostatic interaction, the hydrophobic effect, π-π stacking and van der Waals forces are summarized. The challenges and future perspectives are also discussed.

Kinetic Parameters of trans Scission by Extended HDV-like Ribozymes and the Prospect for the Discovery of Genomic trans-Cleaving RNAs

Hepatitis delta virus (HDV)-like ribozymes are self-cleaving catalytic RNAs with a widespread distribution in nature and biological roles ranging from self-scission during rolling-circle replication in viroids to co-transcriptional processing of eukaryotic retrotransposons, among others. The ribozymes fold into a double pseudoknot with a common catalytic core motif and highly variable peripheral domains. Like other self-cleaving ribozymes, HDV-like ribozymes can be converted into trans-acting catalytic RNAs by bisecting the self-cleaving variants at non-essential loops. Here we explore the trans-cleaving activity of ribozymes derived from the largest examples of the ribozymes (drz-Agam-2 family), which contain an extended domain between the substrate strand and the rest of the RNA. When this peripheral domain is bisected at its distal end, the substrate strand is recognized through two helices, rather than just one 7 bp helix common among the HDV ribozymes, resulting in stronger binding and increased sequence specificity. Kinetic characterization of the extended trans-cleaving ribozyme revealed an efficient trans-cleaving system with a surprisingly high K_M', supporting a model that includes a recently proposed activation barrier related to the assembly of the catalytically competent ribozyme. The ribozymes also exhibit a very long k_off for the products (∼2 weeks), resulting in a trade-off between sequence specificity and turnover. Finally, structure-based searches for the catalytic cores of these ribozymes in the genome of the mosquito Anopheles gambiae, combined with sequence searches for their putative substrates, revealed two potential ribozyme-substrate pairs that may represent examples of natural trans-cleaving ribozymes.

Transcription factor clusters regulate genes in eukaryotic cells

Transcription is regulated through binding factors to gene promoters to activate or repress expression, however, the mechanisms by which factors find targets remain unclear. Using single-molecule fluorescence microscopy, we determined in vivo stoichiometry and spatiotemporal dynamics of a GFP tagged repressor, Mig1, from a paradigm signaling pathway of Saccharomyces cerevisiae. We find the repressor operates in clusters, which upon extracellular signal detection, translocate from the cytoplasm, bind to nuclear targets and turnover. Simulations of Mig1 configuration within a 3D yeast genome model combined with a promoter-specific, fluorescent translation reporter confirmed clusters are the functional unit of gene regulation. In vitro and structural analysis on reconstituted Mig1 suggests that clusters are stabilized by depletion forces between intrinsically disordered sequences. We observed similar clusters of a co-regulatory activator from a different pathway, supporting a generalized cluster model for transcription factors that reduces promoter search times through intersegment transfer while stabilizing gene expression.

Improving the odds: Influence of starting pools on in vitro selection outcomes

As with any outcome of an evolutionary process, the success of in vitro selection experiments depends critically on the starting population. In vitro selections isolate functional nucleic acids that fold into specific structures and form unique binding and catalytic sites. The selection outcomes therefore depend on the molecular and structural diversity of the initial pools. In addition, the experiments are strongly influenced by the length of the starting pool. Longer randomized regions support the formation of more complex structures and presumably allow formation of more intricate tertiary interactions, but they also tend to misfold and aggregate, whereas shorter pools are sufficient to yield simpler motifs. Furthermore, introducing a sequence bias that promotes secondary structure formation appears to prejudice the population towards more functional macromolecules. We review the literature on the influence of the starting pools on the predicted and actual outcomes of laboratory evolution experiments.

Hammerhead Ribozymes in Archaeal Genomes: A Computational Hunt

Hammerhead ribozymes (HHRs) are small self-cleaving RNAs, first discovered in viroids and satellite RNAs of plant viruses. They are composed of a catalytic core of conserved nucleotides flanked by three helices. More recently, hammerhead-encoding sequences have been identified in the genomes of many eukaryotes, prokaryotes and other non-viral species regulating various functions. In this study we have explored the Archaeal domain to identify HHRs using three different bioinformatics approach. Our study reveals four putative hits of HHRs type I and type II in the group Thaumarchaeota and Euryarchaeota in the Archaeal domain, one of which is the instance of HHR 1 in C. symbiosum A, already identified in a previous study. These HHRs are very similar to those previously described in terms of the conservation of their catalytic core. Based on 3-D structure analysis and free energy, these instances were concluded as putative HHRs. Our findings reveal that the catalytic core contains the conserved motifs that are essential for cleavage activity, but there are some instances in which compensatory core variations are present. However, no instances of HHRs have been found in Crenarchaeota. This study reveals a very scarce presence of HHRs in Archaea which suggests the involvement of other ncRNA elements in gene regulatory system like RNase P which are abundantly found in the Archaeal domain.

Chemistry and Biology of Self-Cleaving Ribozymes

Self-cleaving ribozymes were discovered 30 years ago, but their biological distribution and catalytic mechanisms are only beginning to be defined. Each ribozyme family is defined by a distinct structure, with unique active sites accelerating the same transesterification reaction across the families. Biochemical studies show that general acid-base catalysis is the most common mechanism of self-cleavage, but metal ions and metabolites can be used as cofactors. Ribozymes have been discovered in highly diverse genomic contexts throughout nature, from viroids to vertebrates. Their biological roles include self-scission during rolling-circle replication of RNA genomes, co-transcriptional processing of retrotransposons, and metabolite-dependent gene expression regulation in bacteria. Other examples, including highly conserved mammalian ribozymes, suggest that many new biological roles are yet to be discovered.

Discovering human RNA aptamers by structure-based bioinformatics and genome-based in vitro selection

In vitro selection and structure-based searches have emerged as useful techniques for the discoveries of structurally complex RNAs with high affinity and specificity toward metabolites. Here, we focus on the design of a human genomic library that serves as the DNA template for in vitro selection of RNA aptamers. In addition, the structural solutions obtained from the in vitro selection can be used for structure-based searches for discovery of analogous aptamers in various genomic databases.

HDV family of self-cleaving ribozymes

The hepatitis delta virus (HDV) ribozymes are catalytic RNAs capable of cleaving their own sugar-phosphate backbone. The HDV virus possesses the ribozymes in both sense and antisense genomic transcripts, where they are essential for processing during replication. These ribozymes have been the subject of intense biochemical scrutiny and have yielded a wealth of mechanistic insights. In recent years, many HDV-like ribozymes have been identified in nearly all branches of life. The ribozymes are implicated in a variety of biological events, including episodic memory in mammals and retrotransposition in many eukaryotes. Detailed analysis of additional HDV-like ribozyme isolates will likely reveal many more biological functions and provide information about the evolution of this unique RNA.

Identification of minimal HDV-like ribozymes with unique divalent metal ion dependence in the human microbiome

HDV-like self-cleaving ribozymes have been found in a wide variety of organisms, implicated in diverse biological processes, and their activity typically shows a strong divalent metal dependence, but little metal specificity. Recent studies suggested that very short variants of these ribozymes exist in nature, but their distribution and biochemical properties have not been established. To map out the distribution of small HDV-like ribozymes, the drz-Spur-3 sequence was minimized to yield a core construct for structure-based bioinformatic searches. These searches revealed several microbial ribozymes, particularly in the human microbiome. Kinetic profile of the smallest ribozyme revealed two distinct metal binding sites, only one of which promotes fast catalysis. Furthermore, this ribozyme showed markedly reduced activity in Ca(2+), even in the presence of physiological Mg(2+) concentrations. Our study substantially expands the number of microbial HDV-like ribozymes and provides an example of cleavage regulation by divalent metals.

Next-generation sequencing reveals how RNA catalysts evolve from random space

Catalytic RNAs are attractive objects for studying molecular evolution. To understand how RNA libraries can evolve from randomness toward highly active catalysts, we analyze the original samples that led to the discovery of Diels-Alderase ribozymes by next-generation sequencing. Known structure-activity relationships are used to correlate abundance with catalytic performance. We find that efficient catalysts arose not just from selection for reactivity among the members of the starting library, but from improvement of less potent precursors by mutations. We observe changes in the ribozyme population in response to increasing selection pressure. Surprisingly, even after many rounds of enrichment, the libraries are highly diverse, suggesting that potential catalysts are more abundant in random space than generally thought. To highlight the use of next-generation sequencing as a tool for in vitro selections, we also apply this technique to a recent, less characterized ribozyme selection. Making use of the correlation between sequence evolution and catalytic activity, we predict mutations that improve ribozyme activity and validate them biochemically. Our study reveals principles underlying ribozyme in vitro selections and provides guidelines to render future selections more efficient, as well as to predict the conservation of key structural elements, allowing the rational improvement of catalysts.

Detecting and comparing non-coding RNAs in the high-throughput era

In recent years there has been a growing interest in the field of non-coding RNA. This surge is a direct consequence of the discovery of a huge number of new non-coding genes and of the finding that many of these transcripts are involved in key cellular functions. In this context, accurately detecting and comparing RNA sequences has become important. Aligning nucleotide sequences is a key requisite when searching for homologous genes. Accurate alignments reveal evolutionary relationships, conserved regions and more generally any biologically relevant pattern. Comparing RNA molecules is, however, a challenging task. The nucleotide alphabet is simpler and therefore less informative than that of amino-acids. Moreover for many non-coding RNAs, evolution is likely to be mostly constrained at the structural level and not at the sequence level. This results in very poor sequence conservation impeding comparison of these molecules. These difficulties define a context where new methods are urgently needed in order to exploit experimental results to their full potential. This review focuses on the comparative genomics of non-coding RNAs in the context of new sequencing technologies and especially dealing with two extremely important and timely research aspects: the development of new methods to align RNAs and the analysis of high-throughput data.

Convergent evolution of adenosine aptamers spanning bacterial, human, and random sequences revealed by structure-based bioinformatics and genomic SELEX

Aptamers are structured macromolecules in vitro evolved to bind molecular targets, whereas in nature they form the ligand-binding domains of riboswitches. Adenosine aptamers of a single structural family were isolated several times from random pools, but they have not been identified in genomic sequences. We used two unbiased methods, structure-based bioinformatics and human genome-based in vitro selection, to identify aptamers that form the same adenosine-binding structure in a bacterium, and several vertebrates, including humans. Two of the human aptamers map to introns of RAB3C and FGD3 genes. The RAB3C aptamer binds ATP with dissociation constants about 10 times lower than physiological ATP concentration, while the minimal FGD3 aptamer binds ATP only cotranscriptionally.

Computational biology of RNA interactions

The biodiversity of the RNA world has been underestimated for decades. RNA molecules are key building blocks, sensors, and regulators of modern cells. The biological function of RNA molecules cannot be separated from their ability to bind to and interact with a wide space of chemical species, including small molecules, nucleic acids, and proteins. Computational chemists, physicists, and biologists have developed a rich tool set for modeling and predicting RNA interactions. These interactions are to some extent determined by the binding conformation of the RNA molecule. RNA binding conformations are approximated with often acceptable accuracy by sequence and secondary structure motifs. Secondary structure ensembles of a given RNA molecule can be efficiently computed in many relevant situations by employing a standard energy model for base pair interactions and dynamic programming techniques. The case of bi-molecular RNA-RNA interactions can be seen as an extension of this approach. However, unbiased transcriptome-wide scans for local RNA-RNA interactions are computationally challenging yet become efficient if the binding motif/mode is known and other external information can be used to confine the search space. Computational methods are less developed for proteins and small molecules, which bind to RNA with very high specificity. Binding descriptors of proteins are usually determined by in vitro high-throughput assays (e.g., microarrays or sequencing). Intriguingly, recent experimental advances, which are mostly based on light-induced cross-linking of binding partners, render in vivo binding patterns accessible yet require new computational methods for careful data interpretation. The grand challenge is to model the in vivo situation where a complex interplay of RNA binders competes for the same target RNA molecule. Evidently, bioinformaticians are just catching up with the impressive pace of these developments.

The ubiquitous hammerhead ribozyme

The hammerhead ribozyme is a small catalytic RNA motif capable of endonucleolytic (self-) cleavage. It is composed of a catalytic core of conserved nucleotides flanked by three helices, two of which form essential tertiary interactions for fast self-scission under physiological conditions. Originally discovered in subviral plant pathogens, its presence in several eukaryotic genomes has been reported since. More recently, this catalytic RNA motif has been shown to reside in a large number of genomes. We review the different approaches in discovering these new hammerhead ribozyme sequences and discuss possible biological functions of the genomic motifs.

Structure-based search and in vitro analysis of self-cleaving ribozymes

Detecting functional RNAs is increasingly accomplished through structure-based searches for patterns of conserved secondary structure. With large amounts of new sequencing data becoming available, there is a greater demand for efficient methods of identifying new RNAs. Here we present a method of identifying self-cleaving ribozymes and characterizing the in vitro activity.

Discovery of RNA motifs using a computational pipeline that allows insertions in paired regions and filtering of candidate sequences

The enormous impact of noncoding RNAs on biology and biotechnology has motivated the development of systematic approaches to their discovery and characterization. Here we present a methodology for reliable detection of genomic ribozymes that centers on pipelined structure-based searches, utilizing two versatile algorithms for structure prediction. RNArobo is a prototype structure-based search package that enables a single search to return all sequences matching a designated motif descriptor, taking into account the possibility of single nucleotide insertions within base-paired regions. These outputs are then filtered through a structure prediction algorithm based on free energy minimization in order to maximize the proportion of catalytically active RNA motifs. This pipeline provides a fast approach to uncovering new catalytic RNAs with known secondary structures and verifying their activity in vitro.

3' Cap-independent translation enhancers of positive-strand RNA plant viruses

Positive-strand RNA plant viruses that are neither 5'-capped nor 3'-polyadenylated use nontraditional mechanisms to recruit ribosomes to the 5'-end of their viral genomes. One strategy employed by some of these viruses involves a type of RNA element, termed the 3' cap-independent translation enhancer (3'CITE), located in or near the 3'-untranslated region of viral RNA genomes. 3'CITEs function to mediate efficient translation of 5'-proximally encoded viral proteins and function by recruiting either translation initiation factors or the 60S ribosomal subunit to the viral RNA. Recent mechanistic and structural studies have revealed important new insights and details of how 3'CITEs are able to facilitate viral translation and allow these viruses to compete efficiently against cellular mRNAs for the host translational machinery.

Processing and translation initiation of non-long terminal repeat retrotransposons by hepatitis delta virus (HDV)-like self-cleaving ribozymes

Many non-long terminal repeat (non-LTR) retrotransposons lack internal promoters and are co-transcribed with their host genes. These transcripts need to be liberated before inserting into new loci. Using structure-based bioinformatics, we show that several classes of retrotransposons in phyla-spanning arthropods, nematodes, and chordates utilize self-cleaving ribozymes of the hepatitis delta virus (HDV) family for processing their 5' termini. Ribozyme-terminated retrotransposons include rDNA-specific R2, R4, and R6, telomere-specific SART, and Baggins and RTE. The self-scission of the R2 ribozyme is strongly modulated by the insertion site sequence in the rDNA, with the most common insertion sequences promoting faster processing. The ribozymes also promote translation initiation of downstream open reading frames in vitro and in vivo. In some organisms HDV-like and hammerhead ribozymes appear to be dedicated to processing long and short interspersed elements, respectively. HDV-like ribozymes serve several distinct functions in non-LTR retrotransposition, including 5' processing, translation initiation, and potentially trans-templating.

HDV-like self-cleaving ribozymes

HDV ribozymes catalyze their own scission from the transcript during rolling circle replication of the hepatitis delta virus. In vitro selection of self-cleaving ribozymes from a human genomic library revealed an HDV-like ribozyme in the second intron of the human CPEB3 gene and recent results suggest that this RNA affects episodic memory performance. Bioinformatic searches based on the secondary structure of the HDV/CPEB3 fold yielded numerous functional ribozymes in a wide variety of organisms. Genomic mapping of these RNAs suggested several biological roles, one of which is the 5' processing of non-LTR retrotransposons. The family of HDV-like ribozymes thus continues to grow in numbers and biological importance.

Structure-based search reveals hammerhead ribozymes in the human microbiome

Deep sequencing of viral or bacterial nucleic acids monitors the presence and diversity of microbes in select populations and locations. Metagenomic study of mammalian viromes can help trace paths of viral transmissions within or between species. High throughput sequencing of patient and untreated sewage microbiomes showed many sequences with no similarity to genomic sequences of known function or origin. To estimate the distribution of functional RNAs in these microbiomes, we used the hammerhead ribozyme (HHR) motif to search for sequences capable of assuming its three-way junction fold. Although only two of the three possible natural HHR topologies had been known, our analysis revealed highly active ribozymes that terminated in any of the three stems. The most abundant of these are type II HHRs, one of which is the fastest natural cis-acting HHR yet discovered. Altogether, 13 ribozymes were confirmed in vitro, but only one showed sequence similarity to previously described HHRs. Sequences surrounding the ribozymes do not generally show similarity to known genes, except in one case, where a ribozyme is immediately preceded by a bacterial RadC gene. We demonstrate that a structure-based search for a known functional RNA is a powerful tool for analysis of metagenomic datasets, complementing sequence alignments.