EasyDIVER: A Pipeline for Assembling and Counting High-Throughput Sequencing Data from In Vitro Evolution of Nucleic Acids or Peptides

Discovery of Red-Shifting Mutations in Firefly Luciferase Using High-Throughput Biochemistry

Photinus pyralis luciferase (FLuc) has proven a valuable tool for bioluminescence imaging, but much of the light emitted from the native enzyme is absorbed by endogenous biomolecules. Thus, luciferases displaying red-shifted emission enable higher resolution during deep-tissue imaging. A robust model of how protein structure determines emission color would greatly aid the engineering of red-shifted mutants, but no consensus has been reached to date. In this work, we applied deep mutational scanning to systematically assess 20 functionally important amino acid positions on FLuc for red-shifting mutations, predicting that an unbiased approach would enable novel contributions to this debate. We report dozens of red-shifting mutations as a result, a large majority of which have not been previously identified. Further characterization revealed that mutations N229T and T352M, in particular, bring about unimodal emission with the majority of photons being >600 nm. The red-shifting mutations identified by this high-throughput approach provide strong biochemical evidence for the multiple-emitter mechanism of color determination and point to the importance of a water network in the enzyme binding pocket for altering the emitter ratio. This work provides a broadly applicable mutational data set tying FLuc structure to emission color that contributes to our mechanistic understanding of emission color determination and should facilitate further engineering of improved probes for deep-tissue imaging.

Advanced Methods for Natural Products Discovery: Bioactivity Screening, Dereplication, Metabolomics Profiling, Genomic Sequencing, Databases and Informatic Tools, and Structure Elucidation

Natural Products (NP) are essential for the discovery of novel drugs and products for numerous biotechnological applications. The NP discovery process is expensive and time-consuming, having as major hurdles dereplication (early identification of known compounds) and structure elucidation, particularly the determination of the absolute configuration of metabolites with stereogenic centers. This review comprehensively focuses on recent technological and instrumental advances, highlighting the development of methods that alleviate these obstacles, paving the way for accelerating NP discovery towards biotechnological applications. Herein, we emphasize the most innovative high-throughput tools and methods for advancing bioactivity screening, NP chemical analysis, dereplication, metabolite profiling, metabolomics, genome sequencing and/or genomics approaches, databases, bioinformatics, chemoinformatics, and three-dimensional NP structure elucidation.

Emergent properties as by-products of prebiotic evolution of aminoacylation ribozymes

Systems of catalytic RNAs presumably gave rise to important evolutionary innovations, such as the genetic code. Such systems may exhibit particular tolerance to errors (error minimization) as well as coding specificity. While often assumed to result from natural selection, error minimization may instead be an emergent by-product. In an RNA world, a system of self-aminoacylating ribozymes could enforce the mapping of amino acids to anticodons. We measured the activity of thousands of ribozyme mutants on alternative substrates (activated analogs for tryptophan, phenylalanine, leucine, isoleucine, valine, and methionine). Related ribozymes exhibited shared preferences for substrates, indicating that adoption of additional amino acids by existing ribozymes would itself lead to error minimization. Furthermore, ribozyme activity was positively correlated with specificity, indicating that selection for increased activity would also lead to increased specificity. These results demonstrate that by-products of ribozyme evolution could lead to adaptive value in specificity and error tolerance.

REVERSE: a user-friendly web server for analyzing next-generation sequencing data from in vitro selection/evolution experiments

Next-generation sequencing (NGS) enables the identification of functional nucleic acid sequences from in vitro selection/evolution experiments and illuminates the evolutionary process at single-nucleotide resolution. However, analyzing the vast output from NGS can be daunting, especially with limited programming skills. We developed REVERSE (Rapid EValuation of Experimental RNA Selection/Evolution) (https://www.reverseserver.org/), a web server that implements an integrated computational pipeline through a graphical user interface, which performs both pre-processing and detailed sequence level analyses within minutes. Raw FASTQ files are quality-filtered, dereplicated, and trimmed before being analyzed by either of two pipelines. The first pipeline counts, sorts, and tracks enrichment of unique sequences and user-defined sequence motifs. It also identifies mutational intermediates present in the sequence data that connect two input sequences. The second pipeline sorts similar sequences into clusters and tracks enrichment of peak sequences. It also performs nucleotide conservation analysis on the cluster of choice and generates a consensus sequence. Both pipelines generate downloadable spreadsheets and high-resolution figures. Collectively, REVERSE is a one-stop solution for the rapid analysis of NGS data obtained from in vitro selection/evolution experiments that obviates the need for computational expertise.

Sequencing the Origins of Life

One goal of origins of life research is to understand how primitive informational and catalytic biopolymers emerged and evolved. Recently, a number of sequencing techniques have been applied to analysis of replicating and evolving primitive biopolymer systems, providing a sequence-specific and high-resolution view of primitive chemical processes. Here, we review application of sequencing techniques to analysis of synthetic and primitive nucleic acids and polypeptides. This includes next-generation sequencing of primitive polymerization and evolution processes, followed by discussion of other novel biochemical techniques that could contribute to sequence analysis of primitive biopolymer driven chemical systems. Further application of sequencing to origins of life research, perhaps as a life detection technology, could provide insight into the origin and evolution of informational and catalytic biopolymers on early Earth or elsewhere.

PacBio sequencing output increased through uniform and directional fivefold concatenation

Advances in sequencing technology have allowed researchers to sequence DNA with greater ease and at decreasing costs. Main developments have focused on either sequencing many short sequences or fewer large sequences. Methods for sequencing mid-sized sequences of 600-5,000 bp are currently less efficient. For example, the PacBio Sequel I system yields ~ 100,000-300,000 reads with an accuracy per base pair of 90-99%. We sought to sequence several DNA populations of ~ 870 bp in length with a sequencing accuracy of 99% and to the greatest depth possible. We optimised a simple, robust method to concatenate genes of ~ 870 bp five times and then sequenced the resulting DNA of ~ 5,000 bp by PacBioSMRT long-read sequencing. Our method improved upon previously published concatenation attempts, leading to a greater sequencing depth, high-quality reads and limited sample preparation at little expense. We applied this efficient concatenation protocol to sequence nine DNA populations from a protein engineering study. The improved method is accompanied by a simple and user-friendly analysis pipeline, DeCatCounter, to sequence medium-length sequences efficiently at one-fifth of the cost.

Directing evolution of novel ligands by mRNA display

mRNA display is a powerful biological display platform for the directed evolution of proteins and peptides. mRNA display libraries covalently link the displayed peptide or protein (phenotype) with the encoding genetic information (genotype) through the biochemical activity of the small molecule puromycin. Selection for peptide/protein function is followed by amplification of the linked genetic material and generation of a library enriched in functional sequences. Iterative selection cycles are then performed until the desired level of function is achieved, at which time the identity of candidate peptides can be obtained by sequencing the genetic material. The purpose of this review is to discuss the development of mRNA display technology since its inception in 1997 and to comprehensively review its use in the selection of novel peptides and proteins. We begin with an overview of the biochemical mechanism of mRNA display and its variants with a particular focus on its advantages and disadvantages relative to other biological display technologies. We then discuss the importance of scaffold choice in mRNA display selections and review the results of selection experiments with biological (e.g., fibronectin) and linear peptide library architectures. We then explore recent progress in the development of "drug-like" peptides by mRNA display through the post-translational covalent macrocyclization and incorporation of non-proteogenic functionalities. We conclude with an examination of enabling technologies that increase the speed of selection experiments, enhance the information obtained in post-selection sequence analysis, and facilitate high-throughput characterization of lead compounds. We hope to provide the reader with a comprehensive view of current state and future trajectory of mRNA display and its broad utility as a peptide and protein design tool.

Kinetic sequencing (k-Seq) as a massively parallel assay for ribozyme kinetics: utility and critical parameters

Characterizing genotype-phenotype relationships of biomolecules (e.g. ribozymes) requires accurate ways to measure activity for a large set of molecules. Kinetic measurement using high-throughput sequencing (e.g. k-Seq) is an emerging assay applicable in various domains that potentially scales up measurement throughput to over 106 unique nucleic acid sequences. However, maximizing the return of such assays requires understanding the technical challenges introduced by sequence heterogeneity and DNA sequencing. We characterized the k-Seq method in terms of model identifiability, effects of sequencing error, accuracy and precision using simulated datasets and experimental data from a variant pool constructed from previously identified ribozymes. Relative abundance, kinetic coefficients, and measurement noise were found to affect the measurement of each sequence. We introduced bootstrapping to robustly quantify the uncertainty in estimating model parameters and proposed interpretable metrics to quantify model identifiability. These efforts enabled the rigorous reporting of data quality for individual sequences in k-Seq experiments. Here we present detailed protocols, define critical experimental factors, and identify general guidelines to maximize the number of sequences and their measurement accuracy from k-Seq data. Analogous practices could be applied to improve the rigor of other sequencing-based assays.

Encapsulation of ribozymes inside model protocells leads to faster evolutionary adaptation

Functional biomolecules, such as RNA, encapsulated inside a protocellular membrane are believed to have comprised a very early, critical stage in the evolution of life, since membrane vesicles allow selective permeability and create a unit of selection enabling cooperative phenotypes. The biophysical environment inside a protocell would differ fundamentally from bulk solution due to the microscopic confinement. However, the effect of the encapsulated environment on ribozyme evolution has not been previously studied experimentally. Here, we examine the effect of encapsulation inside model protocells on the self-aminoacylation activity of tens of thousands of RNA sequences using a high-throughput sequencing assay. We find that encapsulation of these ribozymes generally increases their activity, giving encapsulated sequences an advantage over nonencapsulated sequences in an amphiphile-rich environment. In addition, highly active ribozymes benefit disproportionately more from encapsulation. The asymmetry in fitness gain broadens the distribution of fitness in the system. Consistent with Fisher's fundamental theorem of natural selection, encapsulation therefore leads to faster adaptation when the RNAs are encapsulated inside a protocell during in vitro selection. Thus, protocells would not only provide a compartmentalization function but also promote activity and evolutionary adaptation during the origin of life.

In vitro selected GUAA tetraloop-binding receptors with structural plasticity and evolvability towards natural RNA structural modules

GNRA tetraloop-binding receptor interactions are key components in the macromolecular assembly of a variety of functional RNAs. In nature, there is an apparent bias for GAAA/11nt receptor and GYRA/helix interactions, with the former interaction being thermodynamically more stable than the latter. While past in vitro selections allowed isolation of novel GGAA and GUGA receptors, we report herein an in vitro selection that revealed several novel classes of specific GUAA receptors with binding affinities comparable to those from natural GAAA/11nt interactions. These GUAA receptors have structural homology with double-locked bulge RNA modules naturally occurring in ribosomal RNAs. They display mutational robustness that enables exploration of the sequence/phenotypic space associated to GNRA/receptor interactions through epistasis. Their thermodynamic self-assembly fitness landscape is characterized by a rugged neutral network with possible evolutionary trajectories toward natural GNRA/receptor interactions. High throughput sequencing analysis revealed synergetic mutations located away from the tertiary interactions that positively contribute to assembly fitness. Our study suggests that the repertoire of GNRA/receptor interactions is much larger than initially thought from the analysis of natural stable RNA molecules and also provides clues for their evolution towards natural GNRA/receptors.