Linking RNA Sequence, Structure, and Function on Massively Parallel High-Throughput Sequencers

Rapid discovery of high-affinity antibodies via massively parallel sequencing, ribosome display and affinity screening

Developing therapeutic antibodies is laborious and costly. Here we report a method for antibody discovery that leverages the Illumina HiSeq platform to, within 3 days, screen in the order of 108 antibody-antigen interactions. The method, which we named 'deep screening', involves the clustering and sequencing of antibody libraries, the conversion of the DNA clusters into complementary RNA clusters covalently linked to the instrument's flow-cell surface on the same location, the in situ translation of the clusters into antibodies tethered via ribosome display, and their screening via fluorescently labelled antigens. By using deep screening, we discovered low-nanomolar nanobodies to a model antigen using 4 × 106 unique variants from yeast-display-enriched libraries, and high-picomolar single-chain antibody fragment leads for human interleukin-7 directly from unselected synthetic repertoires. We also leveraged deep screening of a library of 2.4 × 105 sequences of the third complementarity-determining region of the heavy chain of an anti-human epidermal growth factor receptor 2 (HER2) antibody as input for a large language model that generated new single-chain antibody fragment sequences with higher affinity for HER2 than those in the original library.

Systematic analysis of cotranscriptional RNA folding using transcription elongation complex display

RNA can fold into structures that mediate diverse cellular functions. Understanding how RNA primary sequence directs the formation of functional structures requires methods that can comprehensively assess how changes in an RNA sequence affect its structure and function. Here we have developed a platform for performing high-throughput cotranscriptional RNA biochemical assays, called Transcription Elongation Complex display (TECdisplay). TECdisplay measures RNA function by fractionating a TEC library based on the activity of cotranscriptionally displayed nascent RNA. In this way, RNA function is measured as the distribution of template DNA molecules between fractions of the transcription reaction. This approach circumvents typical RNA sequencing library preparation steps that can cause technical bias. We used TECdisplay to characterize the transcription antitermination activity of 32,768 variants of the Clostridium beijerinckii pfl ZTP riboswitch designed to perturb steps within its cotranscriptional folding pathway. Our findings establish TECdisplay as an accessible platform for high-throughput RNA biochemical assays.

High-throughput biochemistry in RNA sequence space: predicting structure and function

RNAs are central to fundamental biological processes in all known organisms. The set of possible intramolecular interactions of RNA nucleotides defines the range of alternative structural conformations of a specific RNA that can coexist, and these structures enable functional catalytic properties of RNAs and/or their productive intermolecular interactions with other RNAs or proteins. However, the immense combinatorial space of potential RNA sequences has precluded predictive mapping between RNA sequence and molecular structure and function. Recent advances in high-throughput approaches in vitro have enabled quantitative thermodynamic and kinetic measurements of RNA-RNA and RNA-protein interactions, across hundreds of thousands of sequence variations. In this Review, we explore these techniques, how they can be used to understand RNA function and how they might form the foundations of an accurate model to predict the structure and function of an RNA directly from its nucleotide sequence. The experimental techniques and modelling frameworks discussed here are also highly relevant for the sampling of sequence-structure-function space of DNAs and proteins.

Exploring molecular biology in sequence space: The road to next-generation single-molecule biophysics

Next-generation sequencing techniques have led to a new quantitative dimension in the biological sciences. In particular, integrating sequencing techniques with biophysical tools allows sequence-dependent mechanistic studies. Using the millions of DNA clusters that are generated during sequencing to perform high-throughput binding affinity and kinetics measurements enabled the construction of energy landscapes in sequence space, uncovering relationships between sequence, structure, and function. Here, we review the approaches to perform ensemble fluorescence experiments on next-generation sequencing chips for variations of DNA, RNA, and protein sequences. As the next step, we anticipate that these fluorescence experiments will be pushed to the single-molecule level, which can directly uncover kinetics and molecular heterogeneity in an unprecedented high-throughput fashion. Molecular biophysics in sequence space, both at the ensemble and single-molecule level, leads to new mechanistic insights. The wide spectrum of applications in biology and medicine ranges from the fundamental understanding of evolutionary pathways to the development of new therapeutics.

Purification of synchronized E. coli transcription elongation complexes by reversible immobilization on magnetic beads

Synchronized transcription elongation complexes (TECs) are a fundamental tool for in vitro studies of transcription and RNA folding. Transcription elongation can be synchronized by omitting one or more nucleoside triphosphates (NTPs) from an in vitro transcription reaction so that RNA polymerase can only transcribe to the first occurrence of the omitted nucleotide(s) in the coding DNA strand. This approach was developed over four decades ago and has been applied extensively in biochemical investigations of RNA polymerase enzymes, but has not been optimized for RNA-centric assays. In this work, we describe the development of a system for isolating synchronized TECs from an in vitro transcription reaction. Our approach uses a custom 5' leader sequence, called C3-SC1, to reversibly capture synchronized TECs on magnetic beads. We first show using electrophoretic mobility shift and high-resolution in vitro transcription assays that complexes isolated by this procedure, called ^C3-SC1TECs, are >95% pure, >98% active, highly synchronous (94% of complexes chase in <15s upon addition of saturating NTPs), and compatible with solid-phase transcription; the yield of this purification is ∼8%. We then show that ^C3-SC1TECs perturb, but do not interfere with, the function of ZTP (5-aminoimidazole-4-carboxamide riboside 5'-triphosphate)-sensing and ppGpp (guanosine-3',5'-bisdiphosphate)-sensing transcriptional riboswitches. For both riboswitches, transcription using ^C3-SC1TECs improved the efficiency of transcription termination in the absence of ligand but did not inhibit ligand-induced transcription antitermination. Given these properties, ^C3-SC1TECs will likely be useful for developing biochemical and biophysical RNA assays that require high-performance, quantitative bacterial in vitro transcription.

High-Throughput Selection and Characterisation of Aptamers on Optical Next-Generation Sequencers

Aptamers feature a number of advantages, compared to antibodies. However, their application has been limited so far, mainly because of the complex selection process. 'High-throughput sequencing fluorescent ligand interaction profiling' (HiTS-FLIP) significantly increases the selection efficiency and is consequently a very powerful and versatile technology for the selection of high-performance aptamers. It is the first experiment to allow the direct and quantitative measurement of the affinity and specificity of millions of aptamers simultaneously by harnessing the potential of optical next-generation sequencing platforms to perform fluorescence-based binding assays on the clusters displayed on the flow cells and determining their sequence and position in regular high-throughput sequencing. Many variants of the experiment have been developed that allow automation and in situ conversion of DNA clusters into base-modified DNA, RNA, peptides, and even proteins. In addition, the information from mutational assays, performed with HiTS-FLIP, provides deep insights into the relationship between the sequence, structure, and function of aptamers. This enables a detailed understanding of the sequence-specific rules that determine affinity, and thus, supports the evolution of aptamers. Current variants of the HiTS-FLIP experiment and its application in the field of aptamer selection, characterisation, and optimisation are presented in this review.

High-throughput dissection of the thermodynamic and conformational properties of a ubiquitous class of RNA tertiary contact motifs

Despite RNA's diverse secondary and tertiary structures and its complex conformational changes, nature utilizes a limited set of structural "motifs"-helices, junctions, and tertiary contact modules-to build diverse functional RNAs. Thus, in-depth descriptions of a relatively small universe of RNA motifs may lead to predictive models of RNA tertiary conformational landscapes. Motifs may have different properties depending on sequence and secondary structure, giving rise to subclasses that expand the universe of RNA building blocks. Yet we know very little about motif subclasses, given the challenges in mapping conformational properties in high throughput. Previously, we used "RNA on a massively parallel array" (RNA-MaP), a quantitative, high-throughput technique, to study thousands of helices and two-way junctions. Here, we adapt RNA-MaP to study the thermodynamic and conformational properties of tetraloop/tetraloop receptor (TL/TLR) tertiary contact motifs, analyzing 1,493 TLR sequences from different classes. Clustering analyses revealed variability in TL specificity, stability, and conformational behavior. Nevertheless, natural GAAA/11ntR TL/TLRs, while varying in tertiary stability by ∼2.5 kcal/mol, exhibited conserved TL specificity and conformational properties. Thus, RNAs may tune stability without altering the overall structure of these TL/TLRs. Furthermore, their stability correlated with natural frequency, suggesting thermodynamics as the dominant selection pressure. In contrast, other TL/TLRs displayed heterogenous conformational behavior and appear to not be under strong thermodynamic selection. Our results build toward a generalizable model of RNA-folding thermodynamics based on the properties of isolated motifs, and our characterized TL/TLR library can be used to engineer RNAs with predictable thermodynamic and conformational behavior.

µIVC-Useq: a microfluidic-assisted high-throughput functionnal screening in tandem with next generation sequencing and artificial neural network to rapidly characterize RNA molecules

The function of an RNA is intimately linked to its three-dimensional structure. X-ray crystallography or NMR allow the fine structural characterization of small RNA (e.g., aptamers) with a precision down to atomic resolution. Yet, these technics are time consuming, laborious and do not inform on mutational robustness and the extent to which a sequence can be modified without altering RNA function, an important set of information to assist RNA engineering. On another hand, thought powerful, in silico predictions still lack the required accuracy. These limitations can be overcome by using high-throughput microfluidic-assisted functional screening technologies, as they allow exploring large mutant libraries in a rapid and cost-effective manner. Among them, we recently introduced the microfluidic-assisted In Vitro Compartmentalization (µIVC), an efficient screening strategy in which reactions are performed in picoliter droplets at rates of several thousand per second. We later improved µIVC efficiency by using in tandem with high throughput sequencing, thought a laborious bioinformatic step was still required at the end of the process. In the present work, we strongly increased the automation level of the pipeline by implementing an artificial neural network enabling unsupervised bioinformatic analysis. We demonstrate the efficiency of this "µIVC-Useq" technology by rapidly identifying a set of sequences readily accepted by a key domain of the light-up RNA aptamer SRB-2. This work not only shed some new light on the way this aptamer can be engineered, but it also allowed us to easily identify new variants with an up-to 10-fold improved performance.

High throughput approaches to study RNA-protein interactions in vitro

To understand the regulation of gene expression it is critical to determine how proteins interact with and discriminate between different RNAs. In this review, we discuss experimental techniques that utilize high throughput approaches to characterize the interactions of proteins with large numbers of RNAs in vitro. We describe the underlying principles for the main methods, briefly discuss their scope and limitations, and outline how insight from the techniques contributes to our understanding of specificity for RNA-protein interactions.

Application of magnetic nanoparticles in nucleic acid detection

Nucleic acid is the main material for storing, copying, and transmitting genetic information. Gene sequencing is of great significance in DNA damage research, gene therapy, mutation analysis, bacterial infection, drug development, and clinical diagnosis. Gene detection has a wide range of applications, such as environmental, biomedical, pharmaceutical, agriculture and forensic medicine to name a few. Compared with Sanger sequencing, high-throughput sequencing technology has the advantages of larger output, high resolution, and low cost which greatly promotes the application of sequencing technology in life science research. Magnetic nanoparticles, as an important part of nanomaterials, have been widely used in various applications because of their good dispersion, high surface area, low cost, easy separation in buffer systems and signal detection. Based on the above, the application of magnetic nanoparticles in nucleic acid detection was reviewed.

Zipcode Binding Protein 1 (ZBP1; IGF2BP1): A Model for Sequence-Specific RNA Regulation

The fate of an RNA, from its localization, translation, and ultimate decay, is dictated by interactions with RNA binding proteins (RBPs). β-actin mRNA has functioned as the classic example of RNA localization in eukaryotic cells. Studies of β-actin mRNA over the past three decades have allowed understanding of how RBPs, such as ZBP1 (IGF2BP1), can control both RNA localization and translational status. Here, we summarize studies of β-actin mRNA and focus on how ZBP1 serves as a model for understanding interactions between RNA and their binding protein(s). Central to the study of RNA and RBPs were technological developments that occurred along the way. We conclude with a future outlook highlighting new technologies that may be used to address still unanswered questions about RBP-mediated regulation of mRNA during its life cycle, within the cell.

Approaches for measuring the dynamics of RNA-protein interactions

RNA-protein interactions are pivotal for the regulation of gene expression from bacteria to human. RNA-protein interactions are dynamic; they change over biologically relevant timescales. Understanding the regulation of gene expression at the RNA level therefore requires knowledge of the dynamics of RNA-protein interactions. Here, we discuss the main experimental approaches to measure dynamic aspects of RNA-protein interactions. We cover techniques that assess dynamics of cellular RNA-protein interactions that accompany biological processes over timescales of hours or longer and techniques measuring the kinetic dynamics of RNA-protein interactions in vitro. This article is categorized under: RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications RNA Interactions with Proteins and Other Molecules > Protein-RNA Recognition RNA Interactions with Proteins and Other Molecules > RNA-Protein Complexes RNA Evolution and Genomics > Ribonomics.

Automated Design of Diverse Stand-Alone Riboswitches

Riboswitches that couple binding of ligands to conformational changes offer sensors and control elements for RNA synthetic biology and medical biotechnology. However, design of these riboswitches has required expert intuition or software specialized to transcription or translation outputs; design has been particularly challenging for applications in which the riboswitch output cannot be amplified by other molecular machinery. We present a fully automated design method called RiboLogic for such “stand-alone” riboswitches and test it via high-throughput experiments on 2875 molecules using RNA-MaP (RNA on a massively parallel array) technology. These molecules consistently modulate their affinity to the MS2 bacteriophage coat protein upon binding of flavin mononucleotide, tryptophan, theophylline, and microRNA miR-208a, achieving activation ratios of up to 20 and significantly better performance than control designs. By encompassing a wide diversity of stand-alone switches and highly quantitative data, the resulting ribologic-solves experimental data set provides a rich resource for further improvement of riboswitch models and design methods.

The Story of RNA Folding, as Told in Epochs

The past decades have witnessed tremendous developments in our understanding of RNA biology. At the core of these advances have been studies aimed at discerning RNA structure and at understanding the forces that influence the RNA folding process. It is easy to take the present state of understanding for granted, but there is much to be learned by considering the path to our current understanding, which has been tortuous, with the birth and death of models, the adaptation of experimental tools originally developed for characterization of protein structure and catalysis, and the development of novel tools for probing RNA. In this review we tour the stages of RNA folding studies, considering them as "epochs" that can be generalized across scientific disciplines. These epochs span from the discovery of catalytic RNA, through biophysical insights into the putative primordial RNA World, to characterization of structured RNAs, the building and testing of models, and, finally, to the development of models with the potential to yield generalizable predictive and quantitative models for RNA conformational, thermodynamic, and kinetic behavior. We hope that this accounting will aid others as they navigate the many fascinating questions about RNA and its roles in biology, in the past, present, and future.