Analyzing and Tuning Ribozyme Activity by Deep Sequencing To Modulate Gene Expression Level in Mammalian Cells

Deep generative design of RNA family sequences

RNA engineering has immense potential to drive innovation in biotechnology and medicine. Despite its importance, a versatile platform for the automated design of functional RNA is still lacking. Here, we propose RNA family sequence generator (RfamGen), a deep generative model that designs RNA family sequences in a data-efficient manner by explicitly incorporating alignment and consensus secondary structure information. RfamGen can generate novel and functional RNA family sequences by sampling points from a semantically rich and continuous representation. We have experimentally demonstrated the versatility of RfamGen using diverse RNA families. Furthermore, we confirmed the high success rate of RfamGen in designing functional ribozymes through a quantitative massively parallel assay. Notably, RfamGen successfully generates artificial sequences with higher activity than natural sequences. Overall, RfamGen significantly improves our ability to design functional RNA and opens up new potential for generative RNA engineering in synthetic biology.

Coupled in vitro transcription/translation based on wheat germ extract for efficient expression from PCR-generated templates in short-time batch reactions

We successfully constructed a coupled in vitro transcription/translation (cIVTT) system based on wheat germ extract (WGE) for efficient expression from PCR-generated DNA templates in short-time (∼3-h) batch reactions. The productivity of this system under optimized conditions was 85 μg (2.8 nmol) per 1 mL of reaction solution (corresponding to 425 μg per 1 mL of WGE), which was about 9-fold higher than that by the conventional batch method using mRNA as a template. The DNA template concentration required for efficient cIVTT was as low as 2.5 nM, which is much lower than those required for other eukaryotic cIVTT systems to maximize their productivity (30-50 nM). The productivity of the present system with a 2.5 nM template was 80-fold and 4-fold higher than that of a commercially available WGE-based cIVTT system with a 2.5 nM and a 40 nM template, respectively. In addition, the present system functioned well in a liposome (i.e., in an artificial cell) without a loss of productivity. Given that WGE-based systems have the advantage of being suitable for the expression of a broad range of proteins, the present cIVTT system is expected to be widely used in future cell-free synthetic biology.

Recent advances in tuning the expression and regulation of genes for constructing microbial cell factories

To overcome environmental problems caused by the use of fossil resources, microbial cell factories have become a promising technique for the sustainable and eco-friendly development of valuable products from renewable resources. Constructing microbial cell factories with high titers, yields, and productivity requires a balance between growth and production; to this end, tuning gene expression and regulation is necessary to optimise and precisely control complicated metabolic fluxes. In this article, we review the current trends and advances in tuning gene expression and regulation and consider their engineering at each of the three stages of gene regulation: genomic, mRNA, and protein. In particular, the technological approaches utilised in a diverse range of genetic-engineering-based tools for the construction of microbial cell factories are reviewed and representative applications of these strategies are presented. Finally, the prospects for strategies and systems for tuning gene expression and regulation are discussed.

High-Throughput Analysis and Engineering of Ribozymes and Deoxyribozymes by Sequencing

Ribozymes and deoxyribozymes are catalytic RNA and DNA, respectively, that catalyze chemical reactions such as self-cleavage or ligation reactions. While some ribozymes are found in nature, a larger variety of ribozymes and deoxyribozymes have been discovered by in vitro selection from random sequences. These catalytic nucleic acids, especially ribozymes, are of fundamental interest because they are crucial for the RNA world hypothesis, which suggests that RNA played a central role in both the propagation of genetic information and catalyzing metabolic reactions in primordial life prior to the emergence of proteins and DNA. On the practical side, catalytic nucleic acids have been extensively engineered for various applications, such as biosensors and genetic devices for synthetic biology. Therefore, it is important to gain a deeper understanding of the sequence-function relationships of ribozymes and deoxyribozymes.Mutational analysis, or measurements of activities of catalytic nucleic acid mutants, is one of the most fundamental approaches for that purpose. Mutations that abolish, reduce, retain, or even increase activity provide useful information about nucleic acid catalysts for engineering and other purposes. However, methods for mutational analysis of ribozymes and deoxyribozymes have not evolved much for decades, requiring tedious and low-throughput assays (e.g., gel electrophoresis) of individually prepared mutants. This has prevented researchers from performing quantitative mutational analysis of ribozymes and deoxyribozymes on a large scale.To address this limitation, we developed a massively parallel ribozyme and deoxyribozyme assay strategy that allows >10⁴ assays using high-throughput sequencing (HTS). We used HTS to literally count the number of cleaved (or ligated) and uncleaved (or unligated) ribozyme (or deoxyribozyme) sequences and calculated the activities of each mutant in a reaction mixture. This simple yet powerful strategy was applied to analyze the mutational effects of various natural and synthetic ribozymes and deoxyribozymes at scales impossible for conventional mutational analysis. These large-scale sequence-function data sets were used to better understand the functional consequences of mutations and to engineer ribozymes for practical applications. Furthermore, these newly available data are motivating researchers to employ more rigorous computational methods to extract additional insights such as structural information and nonlinear effects of multiple mutations. The new HTS-based assay strategy is distinct from and complementary to a related strategy that uses HTS to analyze ribozyme and deoxyribozyme populations subjected to in vitro selection. Postselection sequencing can cover a larger sequence space, although it does not directly quantify the activities of ribozyme and deoxyribozyme mutants. With further advances in DNA sequencing technologies and computational methods, there should be more opportunities to harness the power of HTS to deepen our understanding of catalytic nucleic acids and enhance our ability to engineer them for even more applications.

Robust gene expression control in human cells with a novel universal TetR aptamer splicing module

Fine-tuning of gene expression is desirable for a wide range of applications in synthetic biology. In this context, RNA regulatory devices provide a powerful and highly functional tool. We developed a versatile, robust and reversible device to control gene expression by splicing regulation in human cells using an aptamer that is recognized by the Tet repressor TetR. Upon insertion in proximity to the 5′ splice site, intron retention can be controlled via the binding of TetR to the aptamer. Although we were able to demonstrate regulation for different introns, the genomic context had a major impact on regulation. In consequence, we advanced the aptamer to develop a splice device. Our novel device contains the aptamer integrated into a context of exonic and intronic sequences that create and maintain an environment allowing a reliable and robust splicing event. The exon-born, additional amino acids will then be cleaved off by a self-cleaving peptide. This design allows portability of the splicing device, which we confirmed by demonstrating its functionality in different gene contexts. Intriguingly, our splicing device shows a high dynamic range and low basal activity, i.e. desirable features that often prove a major challenge when implementing synthetic biology in mammalian cell lines.

Comprehensive sequence-to-function mapping of cofactor-dependent RNA catalysis in the glmS ribozyme

Massively parallel, quantitative measurements of biomolecular activity across sequence space can greatly expand our understanding of RNA sequence-function relationships. We report the development of an RNA-array assay to perform such measurements and its application to a model RNA: the core glmS ribozyme riboswitch, which performs a ligand-dependent self-cleavage reaction. We measure the cleavage rates for all possible single and double mutants of this ribozyme across a series of ligand concentrations, determining kcat and KM values for active variants. These systematic measurements suggest that evolutionary conservation in the consensus sequence is driven by maintenance of the cleavage rate. Analysis of double-mutant rates and associated mutational interactions produces a structural and functional mapping of the ribozyme sequence, revealing the catalytic consequences of specific tertiary interactions, and allowing us to infer structural rearrangements that permit certain sequence variants to maintain activity.

High-throughput identification of synthetic riboswitches by barcode-free amplicon-sequencing in human cells

Synthetic riboswitches mediating ligand-dependent RNA cleavage or splicing-modulation represent elegant tools to control gene expression in various applications, including next-generation gene therapy. However, due to the limited understanding of context-dependent structure-function relationships, the identification of functional riboswitches requires large-scale-screening of aptamer-effector-domain designs, which is hampered by the lack of suitable cellular high-throughput methods. Here we describe a fast and broadly applicable method to functionally screen complex riboswitch libraries (~1.8 × 104 constructs) by cDNA-amplicon-sequencing in transiently transfected and stimulated human cells. The self-barcoding nature of each construct enables quantification of differential mRNA levels without additional pre-selection or cDNA-manipulation steps. We apply this method to engineer tetracycline- and guanine-responsive ON- and OFF-switches based on hammerhead, hepatitis-delta-virus and Twister ribozymes as well as U1-snRNP polyadenylation-dependent RNA devices. In summary, our method enables fast and efficient high-throughput riboswitch identification, thereby overcoming a major hurdle in the development cascade for therapeutically applicable gene switches.

The L-platform/L-scaffold framework: a blueprint for RNA-cleaving nucleic acid enzyme design

We develop an L-platform/L-scaffold framework we hypothesize may serve as a blueprint to facilitate site-specific RNA-cleaving nucleic acid enzyme design. Building on the L-platform motif originally described by Suslov and coworkers, we identify new critical scaffolding elements required to anchor a conserved general base guanine ("L-anchor") and bind functionally important metal ions at the active site ("L-pocket"). Molecular simulations, together with a broad range of experimental structural and functional data, connect the L-platform/L-scaffold elements to necessary and sufficient conditions for catalytic activity. We demonstrate that the L-platform/L-scaffold framework is common to five of the nine currently known naturally occurring ribozyme classes (Twr, HPr, VSr, HHr, Psr), and intriguingly from a design perspective, the framework also appears in an artificially engineered DNAzyme (8-17dz). The flexibility of the L-platform/L-scaffold framework is illustrated on these systems, highlighting modularity and trends in the variety of known general acid moieties that are supported. These trends give rise to two distinct catalytic paradigms, building on the classifications proposed by Wilson and coworkers and named for the implicated general base and acid. The "G + A" paradigm (Twr, HPr, VSr) exclusively utilizes nucleobase residues for chemistry, and the "G + M + " paradigm (HHr, 8-17dz, Psr) involves structuring of the "L-pocket" metal ion binding site for recruitment of a divalent metal ion that plays an active role in the chemical steps of the reaction. Finally, the modularity of the L-platform/L-scaffold framework is illustrated in the VS ribozyme where the "L-pocket" assumes the functional role of the "L-anchor" element, highlighting a distinct mechanism, but one that is functionally linked with the hammerhead ribozyme.

Systematic minimization of RNA ligase ribozyme through large-scale design-synthesis-sequence cycles

Template-directed RNA ligation catalyzed by an RNA enzyme (ribozyme) is a plausible and important reaction that could have been involved in transferring genetic information during prebiotic evolution. Laboratory evolution experiments have yielded several classes of ligase ribozymes, but their minimal sequence requirements remain largely unexplored. Because selection experiments strongly favor highly active sequences, less active but smaller catalytic motifs may have been overlooked in these experiments. We used large-scale DNA synthesis and high-throughput ribozyme assay enabled by deep sequencing to systematically minimize a previously laboratory-evolved ligase ribozyme. After designing and evaluating >10 000 sequences, we identified catalytic cores as small as 18 contiguous bases that catalyze template-directed regiospecific RNA ligation. The fact that such a short sequence can catalyze this critical reaction suggests that similarly simple or even simpler motifs may populate the RNA sequence space which could have been accessible to the prebiotic ribozymes.

Massively parallel RNA device engineering in mammalian cells with RNA-Seq

Synthetic RNA-based genetic devices dynamically control a wide range of gene-regulatory processes across diverse cell types. However, the limited throughput of quantitative assays in mammalian cells has hindered fast iteration and interrogation of sequence space needed to identify new RNA devices. Here we report developing a quantitative, rapid and high-throughput mammalian cell-based RNA-Seq assay to efficiently engineer RNA devices. We identify new ribozyme-based RNA devices that respond to theophylline, hypoxanthine, cyclic-di-GMP, and folinic acid from libraries of ~22,700 sequences in total. The small molecule responsive devices exhibit low basal expression and high activation ratios, significantly expanding our toolset of highly functional ribozyme switches. The large datasets obtained further provide conserved sequence and structure motifs that may be used for rationally guided design. The RNA-Seq approach offers a generally applicable strategy for developing broad classes of RNA devices, thereby advancing the engineering of genetic devices for mammalian systems.

Riboswitch identification using Ligase-Assisted Selection for the Enrichment of Responsive Ribozymes (LigASERR)

In vitro selection of ligand-responsive ribozymes can identify rare, functional sequences from large libraries. While powerful, key caveats of this approach include lengthy and demanding experimental workflows; unpredictable experimental outcomes and unknown functionality of enriched sequences in vivo. To address the first of these limitations, we developed Ligase-Assisted Selection for the Enrichment of Responsive Ribozymes (LigASERR). LigASERR is scalable, amenable to automation and requires less time to implement compared to alternative methods. To improve the predictability of experiments, we modeled the underlying selection process, predicting experimental outcomes based on sequence and population parameters. We applied this new methodology and model to the enrichment of a known, in vitro-selected sequence from a bespoke library. Prior to implementing selection, conditions were optimized and target sequence dynamics accurately predicted for the majority of the experiment. In addition to enriching the target sequence, we identified two new, theophylline-activated ribozymes. Notably, all three sequences yielded riboswitches functional in Escherichia coli, suggesting LigASERR and similar in vitro selection methods can be utilized for generating functional riboswitches in this organism.

Aptamer-based and aptazyme-based riboswitches in mammalian cells

Molecular recognition by RNA aptamers has been exploited to control gene expression in response to small molecules in mammalian cells. These mammalian synthetic riboswitches offer attractive features such as small genetic size and lower risk of immunological complications compared to protein-based transcriptional gene switches. The diversity of gene regulatory mechanisms that involve RNA has also inspired the development of mammalian riboswitches that harness various regulatory mechanisms. In this report, recent advances in synthetic riboswitches that function in mammalian cells are reviewed focusing on the regulatory mechanisms they exploit such as mRNA degradation, microRNA processing, and programmed ribosomal frameshifting.

Applications of high-throughput sequencing to analyze and engineer ribozymes

A large number of catalytic RNAs, or ribozymes, have been identified in the genomes of various organisms and viruses. Ribozymes are involved in biological processes such as regulation of gene expression and viral replication, but biological roles of many ribozymes still remain unknown. Ribozymes have also inspired researchers to engineer synthetic ribozymes that function as sensors or gene switches. To gain deeper understanding of the sequence-function relationship of ribozymes and to efficiently engineer synthetic ribozymes, a large number of ribozyme variants need to be examined which was limited to hundreds of sequences by Sanger sequencing. The advent of high-throughput sequencing technologies, however, has allowed us to sequence millions of ribozyme sequences at low cost. This review focuses on the recent applications of high-throughput sequencing to both characterize and engineer ribozymes, to highlight how the large-scale sequence data can advance ribozyme research and engineering.