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.     

Improved deoxyribozymes for synthesis of covalently branched DNA and RNA

A covalently branched nucleic acid can be synthesized by joining the 2′-hydroxyl of the branch-site ribonucleotide of a DNA or RNA strand to the activated 5′-phosphorus of a separate DNA or RNA strand. We have previously used deoxyribozymes to synthesize several types of branched nucleic acids for experiments in biotechnology and biochemistry. Here, we report in vitro selection experiments to identify improved deoxyribozymes for synthesis of branched DNA and RNA. Each of the new deoxyribozymes requires Mn²⁺ as a cofactor, rather than Mg²⁺ as used by our previous branch-forming deoxyribozymes, and each has an initially random region of 40 rather than 22 or fewer combined nucleotides. The deoxyribozymes all function by forming a three-helix-junction (3HJ) complex with their two oligonucleotide substrates. For synthesis of branched DNA, the best new deoxyribozyme, 8LV13, has k_obs on the order of 0.1 min⁻¹, which is about two orders of magnitude faster than our previously identified 15HA9 deoxyribozyme. 8LV13 also functions at closer-to-neutral pH than does 15HA9 (pH 7.5 versus 9.0) and has useful tolerance for many DNA substrate sequences. For synthesis of branched RNA, two new deoxyribozymes, 8LX1 and 8LX6, were identified with broad sequence tolerances and substantial activity at pH 7.5, versus pH 9.0 for many of our previous deoxyribozymes that form branched RNA. These experiments provide new, and in key aspects improved, practical catalysts for preparation of synthetic branched DNA and RNA.     

In vitro selection of small RNA-cleaving deoxyribozymes that cleave pyrimidine-pyrimidine junctions

Herein, we sought new or improved endoribonucleases based on catalytic DNA molecules known as deoxyribozymes. The current repertoire of RNA-cleaving deoxyribozymes can cleave nearly all of the 16 possible dinucleotide junctions with rates of at least 0.1/min, with the exception of pyrimidine–pyrimidine (pyr–pyr) junctions, which are cleaved 1–3 orders of magnitude slower. We conducted four separate in vitro selection experiments to target each pyr–pyr dinucleotide combination (i.e. CC, UC, CT and UT) within a chimeric RNA/DNA substrate. We used a library of DNA molecules containing only 20 random-sequence nucleotides, so that all possible sequence permutations could be sampled in each experiment. From a total of 245 clones, we identified 22 different sequence families, of which 21 represented novel deoxyribozyme motifs. The fastest deoxyribozymes exhibited k_obs values (single-turnover, intermolecular format) of 0.12/min, 0.04/min, 0.13/min and 0.15/min against CC, UC, CT and UT junctions, respectively. These values represent a 6- to 8-fold improvement for CC and UC junctions, and a 1000- to 1600-fold improvement for CT and UT junctions, compared to the best rates reported previously under identical reaction conditions. The same deoxyribozymes exhibited ~1000-fold lower activity against all RNA substrates, but could potentially be improved through further in vitro evolution and engineering.     

Characterization of non-8-17 sequences uncovers structurally diverse RNA-cleaving deoxyribozymes

RNA-cleaving deoxyribozymes (DNAzymes) can be isolated from random-sequence DNA pools via the process of in vitro selection. However, small and simple catalytic motifs, such as the 8-17 DNAzyme, are commonly observed in sequence space, presenting a challenge in discovering large and complex DNAzymes. In an effort to investigate underrepresented molecular species derived from in vitro selection, in this study we sought to characterize non-8-17 sequences obtained from a previous in vitro selection experiment wherein the 8-17 deoxyribozyme was the dominant motif. We examined 9 sequence families from 21 motifs by characterizing their structural and functional features. We discovered 9 novel deoxyribozyme classes with large catalytic domains (>40 nucleotides) utilizing three-way or four-way junction structural frameworks. Kinetic studies revealed that these deoxyribozymes exhibit moderate to excellent catalytic rates (k(obs) from 0.003 to 1 min(-1)), compared to other known RNA-cleaving DNAzymes. Although chemical probing experiments, site-directed mutational analyses, and metal cofactor dependency tests suggest unique catalytic cores for each deoxyribozyme, common dinucleotide junction selectivity was observed between DNAzymes with similar secondary structural features. Together, our findings indicate that larger, structurally more complex, and diverse catalytic motifs are able to survive the process of in vitro selection despite a sequence space dominated by smaller and structurally simpler catalysts.     

脱氧核酶的研究进展

具有RNA裂解活性的DNA分子称为脱氧核酶。它是经过体外选择技术经多次筛选获得的。脱氧核酶在切割与其互补的RNA底物分子时具有极高的特异性和切割效率,有望成为新的RNA灭活工具。     

Characterization of long RNA-cleaving deoxyribozymes with short catalytic cores: the effect of excess sequence elements on the outcome of in vitro selection

We previously conducted an in vitro selection experiment for RNA-cleaving deoxyribozymes, using a combinatorial DNA library containing 80 random nucleotides. Ultimately, 110 different sequence classes were isolated, but the vast majority contained a short14-15 nt catalytic DNA motif commonly known as 8-17. Herein, we report extensive truncation experiments conducted on multiple sequence classes to confirm the suspected catalytic role played by 8-17 and to determine the effect of excess sequence elements on the activity of this motif and the outcome of selection. Although we observed beneficial, detrimental and neutral consequences for activity, the magnitude of the effect rarely exceeded 2-fold. These deoxyribozymes appear to have survived increasing selection pressure despite the presence of additional sequence elements, rather than because of them. A new deoxyribozyme with comparable activity, called G15-30, was approximately 2.5-fold larger and experienced a approximately 4-fold greater inhibitory effect from excess sequence elements than the average 8-17 motif. Our results suggest that 8-17 may be less susceptible to the potential inhibitory effects of excess arbitrary sequence than larger motifs, which represents a previously unappreciated selective advantage that may contribute to its widespread recurrence.     

Secondary-structure characterization of two proficient kinase deoxyribozymes

Dk1 and Dk2 are two catalytically proficient, manganese-dependent, guanine-rich deoxyribozymes previously isolated for DNA phosphorylation. In this study, we carried out a series of experiments that aimed to understand the structural properties of Dk1 and Dk2 and compare the structural similarities or differences of these two distinct deoxyribozymes that carry out similar catalytic functions. First, we performed reselections from two partially randomized DNA libraries on the basis of the original Dk1 and Dk2 sequences to isolate catalytically active sequence variants and identify nucleotides that are invariable, well-conserved, or highly mutagenized. Sequence analysis of these variants assisted by secondary-structure predictions led to the identification of possible Watson-Crick base-pairing regions within each deoxyribozyme. Sequence truncation and base-pair partner exchange experiments were conducted to confirm, or rule out, the existence of the predicted secondary-structure elements. Finally, methylation interference experiments were applied to identify nucleotides that are potentially important for the tertiary structure folding of the deoxyribozymes. Our data suggest that Dk1 and Dk2, despite the differences in their primary sequences and NTP requirements, use an analogous stem-loop element to anchor a structural domain of substantial tertiary interactions to execute their catalytic functions.     

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).     

Monitoring Genomic Sequences during SELEX Using High-Throughput Sequencing: Neutral SELEX

Background

SELEX is a well established in vitro selection tool to analyze the structure of ligand-binding nucleic acid sequences called aptamers. Genomic SELEX transforms SELEX into a tool to discover novel, genomically encoded RNA or DNA sequences binding a ligand of interest, called genomic aptamers. Concerns have been raised regarding requirements imposed on RNA sequences undergoing SELEX selection.

Methodology/Principal Findings

To evaluate SELEX and assess the extent of these effects, we designed and performed a Neutral SELEX experiment omitting the selection step, such that the sequences are under the sole selective pressure of SELEX''s amplification steps. Using high-throughput sequencing, we obtained thousands of full-length sequences from the initial genomic library and the pools after each of the 10 rounds of Neutral SELEX. We compared these to sequences obtained from a Genomic SELEX experiment deriving from the same initial library, but screening for RNAs binding with high affinity to the E. coli regulator protein Hfq. With each round of Neutral SELEX, sequences became less stable and changed in nucleotide content, but no sequences were enriched. In contrast, we detected substantial enrichment in the Hfq-selected set with enriched sequences having structural stability similar to the neutral sequences but with significantly different nucleotide selection.

Conclusions/Significance

Our data indicate that positive selection in SELEX acts independently of the neutral selective requirements imposed on the sequences. We conclude that Genomic SELEX, when combined with high-throughput sequencing of positively and neutrally selected pools, as well as the gnomic library, is a powerful method to identify genomic aptamers.     

An in vitro-selected RNA receptor for the GAAC loop: modular receptor for non–GNRA-type tetraloop

Although artificial RNA motifs that can functionally replace the GNRA/receptor interaction, a class of RNA–RNA interacting motifs, were isolated from RNA libraries and used to generate designer RNA structures, receptors for non-GNRA tetraloops have not been found in nature or selected from RNA libraries. In this study, we report successful isolation of a receptor motif interacting with GAAC, a non-GNRA tetraloop, from randomized sequences embedded in a catalytic RNA. Biochemical characterization of the GAAC/receptor interacting motif within three structural contexts showed its binding affinity, selectivity and structural autonomy. The motif has binding affinity comparable with that of a GNRA/receptor, selectivity orthogonal to GNRA/receptors and structural autonomy even in a large RNA context. These features would be advantageous for usage of the motif as a building block for designer RNAs. The isolated motif can also be used as a query sequence to search for unidentified naturally occurring GANC receptor motifs.     

Prospective identification of parasitic sequences in phage display screens

Phage display empowered the development of proteins with new function and ligands for clinically relevant targets. In this report, we use next-generation sequencing to analyze phage-displayed libraries and uncover a strong bias induced by amplification preferences of phage in bacteria. This bias favors fast-growing sequences that collectively constitute <0.01% of the available diversity. Specifically, a library of 10⁹ random 7-mer peptides (Ph.D.-7) includes a few thousand sequences that grow quickly (the ‘parasites’), which are the sequences that are typically identified in phage display screens published to date. A similar collapse was observed in other libraries. Using Illumina and Ion Torrent sequencing and multiple biological replicates of amplification of Ph.D.-7 library, we identified a focused population of 770 ‘parasites’. In all, 197 sequences from this population have been identified in literature reports that used Ph.D.-7 library. Many of these enriched sequences have confirmed function (e.g. target binding capacity). The bias in the literature, thus, can be viewed as a selection with two different selection pressures: (i) target-binding selection, and (ii) amplification-induced selection. Enrichment of parasitic sequences could be minimized if amplification bias is removed. Here, we demonstrate that emulsion amplification in libraries of ∼10⁶ diverse clones prevents the biased selection of parasitic clones.     

In vitro selection, characterization, and application of deoxyribozymes that cleave RNA

Over the last decade, many catalytically active DNA molecules (deoxyribozymes; DNA enzymes) have been identified by in vitro selection from random-sequence DNA pools. This article focuses on deoxyribozymes that cleave RNA substrates. The first DNA enzyme was reported in 1994 and cleaves an RNA linkage. Since that time, many other RNA-cleaving deoxyribozymes have been identified. Most but not all of these deoxyribozymes require a divalent metal ion cofactor such as Mg²⁺ to catalyze attack by a specific RNA 2′-hydroxyl group on the adjacent phosphodiester linkage, forming a 2′,3′-cyclic phosphate and a 5′-hydroxyl group. Several deoxyribozymes that cleave RNA have utility for in vitro RNA biochemistry. Some DNA enzymes have been applied in vivo to degrade mRNAs, and others have been engineered into sensors. The practical impact of RNA-cleaving deoxyribozymes should continue to increase as additional applications are developed.     

Directing the outcome of deoxyribozyme selections to favor native 3'-5' RNA ligation

Previous experiments have identified numerous RNA ligase deoxyribozymes, each of which can synthesize either 2',5'-branched RNA, linear 2'-5'-linked RNA, or linear 3'-5'-linked RNA. These products may be formed by reaction of a 2'-hydroxyl or 3'-hydroxyl of one RNA substrate with the 5'-triphosphate of a second RNA substrate. Here the inherent propensities for nucleophilic reactivity of specific hydroxyl groups were assessed using RNA substrates related to the natural sequences of spliceosome substrates and group II introns. With the spliceosome substrates, nearly half of the selected deoxyribozymes mediate a ligation reaction involving the natural branch-point adenosine as the nucleophile. In contrast, mostly linear RNA is obtained with the group II intron substrates. Because the two sets of substrates differ at only three nucleotides, we conclude that the location of the newly created ligation junction in DNA-catalyzed branch formation depends sensitively on the RNA substrate sequences. During the experiment that led primarily to branched RNA, we abruptly altered the selection strategy to demand that the deoxyribozymes create linear 3'-5' linkages by introducing an additional selection step involving the 3'-5'-selective 8-17 deoxyribozyme. Although no 3'-5' linkages (相似文献   

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.     

New DNA-hydrolyzing DNAs isolated from an ssDNA library carrying a terminal hybridization stem

DNA-hydrolyzing DNAs represent an attractive type of DNA-processing catalysts distinctive from the protein-based restriction enzymes. The innate DNA property has enabled them to readily join DNA-based manipulations to promote the development of DNA biotechnology. A major in vitro selection strategy to identify these DNA catalysts relies tightly on the isolation of linear DNAs processed from a circular single-stranded (ss) DNA sequence library by self-hydrolysis. Herein, we report that by programming a terminal hybridization stem in the library, other than the previously reported classes (I & II) of deoxyribozymes, two new classes (III & IV) were identified with the old selection strategy to site-specifically hydrolyze DNA in the presence of Zn²⁺. Their representatives own a catalytic core consisting of ∼20 conserved nucleotides and a half-life of ∼15 min at neutral pH. In a bimolecular construct, class III exhibits unique broad generality on the enzyme strand, which can be potentially harnessed to engineer DNA-responsive DNA hydrolyzers for detection of any target ssDNA sequence. Besides the new findings, this work should also provide an improved approach to select for DNA-hydrolyzing deoxyribozymes that use various molecules and ions as cofactors.     

In silico selection of RNA aptamers

In vitro selection of RNA aptamers that bind to a specific ligand usually begins with a random pool of RNA sequences. We propose a computational approach for designing a starting pool of RNA sequences for the selection of RNA aptamers for specific analyte binding. Our approach consists of three steps: (i) selection of RNA sequences based on their secondary structure, (ii) generating a library of three-dimensional (3D) structures of RNA molecules and (iii) high-throughput virtual screening of this library to select aptamers with binding affinity to a desired small molecule. We developed a set of criteria that allows one to select a sequence with potential binding affinity from a pool of random sequences and developed a protocol for RNA 3D structure prediction. As verification, we tested the performance of in silico selection on a set of six known aptamer–ligand complexes. The structures of the native sequences for the ligands in the testing set were among the top 5% of the selected structures. The proposed approach reduces the RNA sequences search space by four to five orders of magnitude—significantly accelerating the experimental screening and selection of high-affinity aptamers.     

Enhanced crossover SCRATCHY: construction and high-throughput screening of a combinatorial library containing multiple non-homologous crossovers

SCRATCHY is a methodology for the construction of libraries of chimeras between genes that display low sequence homology. We have developed a strategy for library creation termed enhanced crossover SCRATCHY, that significantly increases the number of clones containing multiple crossovers. Complementary chimeric gene libraries generated by incremental truncation (ITCHY) of two distinct parental sequences are created, and are then divided into arbitrarily defined sections. The respective sections are amplified by skewed sets of primers (i.e. a combination of gene A specific forward primer and gene B specific reverse primer, etc.) allowing DNA fragments containing non-homologous crossover points to be amplified. The amplified chimeric sections are then subjected to a DNA shuffling process generating an enhanced crossover SCRATCHY library. We have constructed such a library using the rat theta 2 glutathione transferase (rGSTT2) and the human theta 1 glutathione transferase (hGSTT1) genes (63% DNA sequence identity). DNA sequencing analysis of unselected library members revealed a greater diversity than that obtained by canonical family shuffling or with conventional SCRATCHY. Expression and high-throughput flow cytometric screening of the chimeric GST library identified several chimeric progeny that retained rat-like parental substrate specificity.     

Revitalization of six abandoned catalytic DNA species reveals a common three-way junction framework and diverse catalytic cores

A library containing as many as 10(16) nucleic acid candidates is typically used to isolate artificial ribozymes and deoxyribozymes (DNAzymes) in an in vitro selection experiment, with only a handful of sequences surviving many rounds of stringent selection steps. These winning species are generally the focus of interest whereas the less competitive contenders are usually not examined. Nevertheless, molecular species abandoned during the selection process might still represent a rich pool of catalytic motifs that are useful for the examination of DNA's inherent catalytic ability, and for the design of molecular tools for practical applications. Here we report a study of six RNA-cleaving, fluorescence-signaling deoxyribozymes that appeared in the early generations of a previous in vitro selection experiment, using the combined approaches of reselection, rational structural analysis, and reaction condition optimization. All six deoxyribozymes were found to use a three-way junction as a common structural framework for catalysis. However, disparities observed in the conserved nucleotide allocations, methylation interference patterns and metal-ion selectivities, pointed to distinct catalytic cores. The rate constants of the optimized deoxyribozymes fell in the range of approximately 0.2 to 1.6 min(-1), which are comparable to those of similar ribozymes. Our findings indicate that deoxyribozymes eliminated by harsh selection criteria are structurally simple molecules that can be tailored into efficient catalysts.     

一种结合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结合蛋白识别底物的序列特异性,并进一步了解其在生物体内的调控机制有重要意义。