Computational generation and screening of RNA motifs in large nucleotide sequence pools

Although identification of active motifs in large random sequence pools is central to RNA in vitro selection, no systematic computational equivalent of this process has yet been developed. We develop a computational approach that combines target pool generation, motif scanning and motif screening using secondary structure analysis for applications to 10¹²–10¹⁴-sequence pools; large pool sizes are made possible using program redesign and supercomputing resources. We use the new protocol to search for aptamer and ribozyme motifs in pools up to experimental pool size (10¹⁴ sequences). We show that motif scanning, structure matching and flanking sequence analysis, respectively, reduce the initial sequence pool by 6–8, 1–2 and 1 orders of magnitude, consistent with the rare occurrence of active motifs in random pools. The final yields match the theoretical yields from probability theory for simple motifs and overestimate experimental yields, which constitute lower bounds, for aptamers because screening analyses beyond secondary structure information are not considered systematically. We also show that designed pools using our nucleotide transition probability matrices can produce higher yields for RNA ligase motifs than random pools. Our methods for generating, analyzing and designing large pools can help improve RNA design via simulation of aspects of in vitro selection.     

Predicting candidate genomic sequences that correspond to synthetic functional RNA motifs

Riboswitches and RNA interference are important emerging mechanisms found in many organisms to control gene expression. To enhance our understanding of such RNA roles, finding small regulatory motifs in genomes presents a challenge on a wide scale. Many simple functional RNA motifs have been found by in vitro selection experiments, which produce synthetic target-binding aptamers as well as catalytic RNAs, including the hammerhead ribozyme. Motivated by the prediction of Piganeau and Schroeder [(2003) Chem. Biol., 10, 103–104] that synthetic RNAs may have natural counterparts, we develop and apply an efficient computational protocol for identifying aptamer-like motifs in genomes. We define motifs from the sequence and structural information of synthetic aptamers, search for sequences in genomes that will produce motif matches, and then evaluate the structural stability and statistical significance of the potential hits. Our application to aptamers for streptomycin, chloramphenicol, neomycin B and ATP identifies 37 candidate sequences (in coding and non-coding regions) that fold to the target aptamer structures in bacterial and archaeal genomes. Further energetic screening reveals that several candidates exhibit energetic properties and sequence conservation patterns that are characteristic of functional motifs. Besides providing candidates for experimental testing, our computational protocol offers an avenue for expanding natural RNA's functional repertoire.     

Finding specific RNA motifs: function in a zeptomole world?

We have developed a new method for estimating the abundance of any modular (piecewise) RNA motif within a longer random region. We have used this method to estimate the size of the active motifs available to modern SELEX experiments (picomoles of unique sequences) and to a plausible RNA World (zeptomoles of unique sequences: 1 zmole = 602 sequences). Unexpectedly, activities such as specific isoleucine binding are almost certainly present in zeptomoles of molecules, and even ribozymes such as self-cleavage motifs may appear (depending on assumptions about the minimal structures). The number of specified nucleotides is not the only important determinant of a motif's rarity: The number of modules into which it is divided, and the details of this division, are also crucial. We propose three maxims for easily isolated motifs: the Maxim of Minimization, the Maxim of Multiplicity, and the Maxim of the Median. These maxims together state that selected motifs should be small and composed of as many separate, equally sized modules as possible. For evenly divided motifs with four modules, the largest accessible activity in picomole scale (1-1000 pmole) pools of length 100 is about 34 nucleotides; while for zeptomole scale (1-1000 zmole) pools it is about 20 specific nucleotides (50% probability of occurrence). This latter figure includes some ribozymes and aptamers. Consequently, an RNA metabolism apparently could have begun with only zeptomoles of RNA molecules.     

Characterisation of aptamer–target interactions by branched selection and high-throughput sequencing of SELEX pools

Nucleic acid aptamer selection by systematic evolution of ligands by exponential enrichment (SELEX) has shown great promise for use in the development of research tools, therapeutics and diagnostics. Typically, aptamers are identified from libraries containing up to 10¹⁶ different RNA or DNA sequences by 5–10 rounds of affinity selection towards a target of interest. Such library screenings can result in complex pools of many target-binding aptamers. New high-throughput sequencing techniques may potentially revolutionise aptamer selection by allowing quantitative assessment of the dynamic changes in the pool composition during the SELEX process and by facilitating large-scale post-SELEX characterisation. In the present study, we demonstrate how high-throughput sequencing of SELEX pools, before and after a single round of branched selection for binding to different target variants, can provide detailed information about aptamer binding sites, preferences for specific target conformations, and functional effects of the aptamers. The procedure was applied on a diverse pool of 2′-fluoropyrimidine-modified RNA enriched for aptamers specific for the serpin plasminogen activator inhibitor-1 (PAI-1) through five rounds of standard selection. The results demonstrate that it is possible to perform large-scale detailed characterisation of aptamer sequences directly in the complex pools obtained from library selection methods, thus without the need to produce individual aptamers.     

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.     

RNA aptamers to the adenosine moiety of S-adenosyl methionine: structural inferences from variations on a theme and the reproducibility of SELEX.

We used in vitro selection (SELEX) to isolate RNA 'aptamers' to S-adenosyl methionine (SAM). Individual aptamer sequences conform to the structural element noted previously for adenosine binding in selections for aptamers to ATP and NAD+. When we compare the patterns of sequence conservation among 65 adenosine-binding sequences to the published structure of the adenosine aptamer, we find that the most highly conserved nucleotides contact the bound adenosine directly, and that one conserved nucleotide outside the binding pocket is in position to stabilize nucleotides within the binding pocket. The aptamer's ability to bind diverse adenosine-containing cofactors is easily understood in terms of its mode of binding, which leaves the 5'position exposed to solvent. We propose that aptamers that bind their targets away from the reactive moiety may be particularly well suited for catalysis. Finally, we estimate that one sequence in 10(11) may be able to form this structural motif, and that there may be many other adenosine-binding motifs that have escaped detection because of their lower representation in the starting random pools.     

Size, constant sequences, and optimal selection

Because the abundance of functional molecules in RNA sequence space has many unexplored aspects, we compared the outcome of 11 independent selections, performed using the same affinity selection protocol and contiguous randomized regions of 16, 22, 26, 50, 70, and 90 nucleotides. All affinity selections targeted the simplest isoleucine aptamer, an asymmetric internal loop. This loop should be abundant in all selections, so that it can be compared across all experiments. In some cases, two primer sets intended to favor selection of different structures have also been compared. The simplest isoleucine aptamer dominates all selections except with the shortest tract, 16 contiguous randomized nucleotides. Here the isoleucine aptamer cannot be accommodated and no other motif can be selected. Our results suggest an optimum length for selection; surprisingly, both the shortest and the longest randomized tracts make it more difficult to recover the motif. Estimated apparent initial abundances suggest that the simplest isoleucine motif was 20- to 40-fold more frequent in selection with 50- or 70-nucleotide randomized regions than with any other length. Considering primer sets, a pre-formed stable stem within fixed flanking sequences had a five-to 10-fold negative effect on apparent motif abundance at all lengths. Differing random tract lengths also determined the probable motif permutation and the most abundant helix lengths. These data support a significant but lesser role for primer sequences in the outcome of selections.     

An RNA aptamer that discriminates bovine factor IX from human factor IX

An RNA aptamer has been selected by SELEX against bovine factor IX using an RNA pool containing 74-nucleotides randomized region. Selected RNA aptamer (Clone 5) could discriminate bovine factor IX effectively from human factor IX. Interestingly, the nucleotide regions 73-78 and 80-83 of the selected aptamer were determined to be important for bovine factor IX-binding using phosphate interference. Based on phosphate interference and binding studies the minimal motif for aptamer with discriminating ability is found with the nucleotide regions from 65 to 106. The discriminating ability of this mini aptamer is calculated as more than 1,000 fold. The equilibrium dissociation constant (K(d)) for the above complex was 10 nM as determined by surface plasmon resonance. Based on the available structural informations, probable binding site of aptamer on the target was predicted.     

Computational approaches toward the design of pools for the in vitro selection of complex aptamers

It is well known that using random RNA/DNA sequences for SELEX experiments will generally yield low-complexity structures. Early experimental results suggest that having a structurally diverse library, which, for instance, includes high-order junctions, may prove useful in finding new functional motifs. Here, we develop two computational methods to generate sequences that exhibit higher structural complexity and can be used to increase the overall structural diversity of initial pools for in vitro selection experiments. Random Filtering selectively increases the number of five-way junctions in RNA/DNA pools, and Genetic Filtering designs RNA/DNA pools to a specified structure distribution, whether uniform or otherwise. We show that using our computationally designed DNA pool greatly improves access to highly complex sequence structures for SELEX experiments (without losing our ability to select for common one-way and two-way junction sequences).     

Enzymatic and antisense effects of a specific anti-Ki-ras ribozyme in vitro and in cell culture.

Due to their mode of action, ribozymes show antisense effects in addition to their specific cleavage activity. In the present study we investigated whether a hammerhead ribozyme is capable of cleaving mutated Ki-ras mRNA in a pancreatic carcinoma cell line and whether antisense effects contribute to the activity of the ribozyme. A 2[prime]-O-allyl modified hammerhead ribozyme was designed to cleave specifically the mutated form of the Ki- ras mRNA (GUU motif in codon 12). The activity was monitored by RT-PCR on Ki- ras RNA expression by determination of the relative amount of wild type to mutant Ki-ras mRNA, by 5-bromo-2[prime]-deoxy-uridine incorporation on cell proliferation and by colony formation in soft agar on malignancy in the human pancreatic adenocarcinoma cell line CFPAC-1, which is heterozygous for the Ki-ras mutation. A catalytically inactive ribozyme was used as control to differentiate between antisense and cleavage activity and a ribozyme with random guide sequences as negative control. The catalytically active anti-Ki-ras ribozyme was at least 2-fold more potent in decreasing cellular Ki-ras mRNA levels, inhibiting cell proliferation and colony formation in soft agar than the catalytically inactive ribozyme. The catalytically active anti-Ki-ras ribozyme, but not the catalytically inactive or random ribozyme, increased the ratio of wild type to mutated Ki-ras mRNA in CFPAC-1 cells. In conclusion, both cleavage activity and antisense effects contribute to the activity of the catalytically active anti-Ki-ras hammerhead ribozyme. Specific ribozymes might be useful in the treatment of pancreatic carcinomas containing an oncogenic GTT mutation in codon 12 of the Ki-ras gene.     

Metal ion specificities for folding and cleavage activity in the Schistosoma hammerhead ribozyme

The effects of various metal ions on cleavage activity and global folding have been studied in the extended Schistosoma hammerhead ribozyme. Fluorescence resonance energy transfer was used to probe global folding as a function of various monovalent and divalent metal ions in this ribozyme. The divalent metals ions Ca²⁺, Mg²⁺, Mn²⁺, and Sr²⁺ have a relatively small variation (less than sixfold) in their ability to globally fold the hammerhead ribozyme, which contrasts with the very large difference (>10,000-fold) in apparent rate constants for cleavage for these divalent metal ions in single-turnover kinetic experiments. There is still a very large range (>4600-fold) in the apparent rate constants for cleavage for these divalent metal ions measured in high salt (2 M NaCl) conditions where the ribozyme is globally folded. These results demonstrate that the identity of the divalent metal ion has little effect on global folding of the Schistosoma hammerhead ribozyme, whereas it has a very large effect on the cleavage kinetics. Mechanisms by which the identity of the divalent metal ion can have such a large effect on cleavage activity in the Schistosoma hammerhead ribozyme are discussed.     

一种结合SELEX与二代测序技术筛选RNA结合蛋白特异识别序列新方法的建立

RNA结合蛋白通过特异识别RNA底物发挥重要的生物学作用。指数富集的配体系统进化(Systematic evolution of ligands by exponential enrichment,SELEX)技术是一种体外筛选核酸底物的基本方法,SELEX技术通过重复多轮筛选从随机核酸序列库中筛选出特异性与靶物质高度亲和的核酸底物,本研究将利用该技术与二代高通量测序(NGS)相结合,体外合成含有20个随机碱基的RNA文库,将所要研究的蛋白构建到带有可被链亲和酶素磁珠捕获的SBP标记的载体上去,显著提高筛选效率,仅需1轮筛选即可获得所需RNA底物motif。通过该方法获得了人的hn RNP A1的UP1结构域特异识别AGG和AG二种RNA序列,并通过EMSA实验证实其可以与获得的RNA motif结合。这一方法的建立对于研究RNA结合蛋白识别底物的序列特异性,并进一步了解其在生物体内的调控机制有重要意义。     

A computational proposal for designing structured RNA pools for in vitro selection of RNAs

Although in vitro selection technology is a versatile experimental tool for discovering novel synthetic RNA molecules, finding complex RNA molecules is difficult because most RNAs identified from random sequence pools are simple motifs, consistent with recent computational analysis of such sequence pools. Thus, enriching in vitro selection pools with complex structures could increase the probability of discovering novel RNAs. Here we develop an approach for engineering sequence pools that links RNA sequence space regions with corresponding structural distributions via a "mixing matrix" approach combined with a graph theory analysis. We define five classes of mixing matrices motivated by covariance mutations in RNA; these constructs define nucleotide transition rates and are applied to chosen starting sequences to yield specific nonrandom pools. We examine the coverage of sequence space as a function of the mixing matrix and starting sequence via clustering analysis. We show that, in contrast to random sequences, which are associated only with a local region of sequence space, our designed pools, including a structured pool for GTP aptamers, can target specific motifs. It follows that experimental synthesis of designed pools can benefit from using optimized starting sequences, mixing matrices, and pool fractions associated with each of our constructed pools as a guide. Automation of our approach could provide practical tools for pool design applications for in vitro selection of RNAs and related problems.     

Exploring the sequence space of a DNA aptamer using microarrays

The relationship between sequence and binding properties of an aptamer for immunoglobulin E (IgE) was investigated using custom DNA microarrays. Single, double and some triple mutations of the aptamer sequence were created to evaluate the importance of specific base composition on aptamer binding. The majority of the positions in the aptamer sequence were found to be immutable, with changes at these positions resulting in more than a 100-fold decrease in binding affinity. Improvements in binding were observed by altering the stem region of the aptamer, suggesting that it plays a significant role in binding. Results obtained for the various mutations were used to estimate the information content and the probability of finding a functional aptamer sequence by selection from a random library. For the IgE-binding aptamer, this probability is on the order of 10⁻¹⁰ to 10⁻⁹. Results obtained for the double and triple mutations also show that there are no compensatory mutations within the space defined by those mutations. Apparently, at least for this particular aptamer, the functional sequence space can be represented as a rugged landscape with sharp peaks defined by highly constrained base compositions. This makes the rational optimization of aptamer sequences using step-wise mutagenesis approaches very challenging.     

Selection of tetracycline inducible self-cleaving ribozymes as synthetic devices for gene regulation in yeast

Synthetic regulatory devices are key components for the development of complex biological systems and the reprogramming of cellular functions and networks. Here we describe the selection of tetracycline inducible hammerhead ribozymes. A tetracycline aptamer was fused to the full-length hammerhead ribozyme via a variable linker region. 11 rounds of in vitro selection were applied to isolate linker sequences that mediate tetracycline dependent hammerhead cleavage. We identified allosteric ribozymes that cleave in the presence of 1 μM tetracycline as fast as the full-length ribozyme whereas cleavage is inhibited up to 333-fold in the absence of tetracycline. Reporter gene assays indicate that the allosteric ribozymes can be employed to control gene expression in yeast.     

Comparative single-turnover kinetic analyses of trans-cleaving hammerhead ribozymes with naturally derived non-conserved sequence motifs

trans-Cleaving hammerhead ribozyme variants were generated with mimicked non-conserved internal loop motifs derived from five structurally diverse natural cis-cleaving ribozymes. Most modified trans-cleaving variants showed enhanced single-turnover cleavage rates relative to minimal counterparts that lack tertiary interactions between internal loop motifs I and II, and relative to controls with sequence changes in loop I. The trans-cleaving ribozyme derived from the positive strand of peach latent mosaic viroid had the highest observed cleavage rate, suggesting a structurally optimized motif that facilitates rapid formation of the ribozyme catalytic center in a trans-reaction.     

Selection of the simplest RNA that binds isoleucine

We have identified the simplest RNA binding site for isoleucine using selection-amplification (SELEX), by shrinking the size of the randomized region until affinity selection is extinguished. Such a protocol can be useful because selection does not necessarily make the simplest active motif most prominent, as is often assumed. We find an isoleucine binding site that behaves exactly as predicted for the site that requires fewest nucleotides. This UAUU motif (16 highly conserved positions; 27 total), is also the most abundant site in successful selections on short random tracts. The UAUU site, now isolated independently at least 63 times, is a small asymmetric internal loop. Conserved loop sequences include isoleucine codon and anticodon triplets, whose nucleotides are required for amino acid binding. This reproducible association between isoleucine and its coding sequences supports the idea that the genetic code is, at least in part, a stereochemical residue of the most easily isolated RNA-amino acid binding structures.     

Design and optimization of effector-activated ribozyme ligases

A selected ribozyme ligase, L1, has been engineered to respond to small organic effectors. Residues important for ribozyme catalysis were mapped to a compact core structure. Aptamers that bound adenosine and theophylline were appended to the core structure, and the resultant aptazymes proved to be responsive to their cognate effectors. Rational sequence substitutions in the joining region between the aptamer and the ribozyme yielded aptazymes whose activities were enhanced from 800–1600-fold in the presence of 1 mM ATP or theophylline, respectively. However, when an anti-flavin aptamer was appended to the core ribozyme structure flavin-responsivity was minimal. The joining region between the aptamer and the ribozyme core was randomized and a series of negative and positive selection steps yielded aptazymes that were activated by up to 260-fold in the presence of 100 µM FMN. The selected joining regions proved to be ‘communication modules’ that could be used to join other aptamers to the ribozyme core to form aptazymes. These results show that ribozyme ligases can be readily engineered to function as allosteric enzymes, and reveal that many of the techniques and principles previously demonstrated during the development of hammerhead aptazymes may be generalizable.