Elucidating sequence-function relationships in a template-independent polymerase to enable novel DNA recording applications.

Harnessing DNA as a high-density storage medium for information storage and molecular recording of signals has been of increasing interest in the biotechnology field. Recently, progress in enzymatic DNA synthesis, DNA digital data storage, and DNA-based molecular recording has been made by leveraging the activity of the template-independent DNA polymerase, terminal deoxynucleotidyl transferase (TdT). TdT adds deoxyribonucleotides to the 3' end of single-stranded DNA, generating random sequences of single-stranded DNA. TdT can use several divalent cations for its enzymatic activity and exhibits shifts in deoxyribonucleotide incorporation frequencies in response to changes in its reaction environment. However, there is limited understanding of sequence-structure-function relationships regarding these properties, which in turn limits our ability to modulate TdT to further advance TdT-based tools. Most TdT literature to-date explores the activity of murine, bovine or human TdTs; studies probing TdT sequence and structure largely focus on strictly conserved residues that are functionally critical to TdT activity. Here, we explore non-conserved TdT sequence space by surveying the natural diversity of TdT. We characterize a diverse set of TdT homologs from different organisms and identify several TdT residues/regions that confer differences in TdT behavior between homologs. The observations in this study can design rules for targeted TdT libraries, in tandem with a screening assay, to modulate TdT properties. Moreover, the data can be useful in guiding further studies of potential residues of interest. Overall, we characterize TdTs that have not been previously studied in the literature, and we provide new insights into TdT sequence-function relationships.

Engineering the Activity of a Template-Independent DNA Polymerase.

Enzymatic DNA writing technologies based on the template-independent DNA polymerase terminal deoxynucleotidyl transferase (TdT) have the potential to advance DNA information storage. TdT is unique in its ability to synthesize single-stranded DNA de novo but has limitations, including catalytic inhibition by ribonucleotide presence and slower incorporation rates compared to replicative polymerases. We anticipate that protein engineering can improve, modulate, and tailor the enzyme's properties, but there is limited information on TdT sequence-structure-function relationships to facilitate rational approaches. Therefore, we developed an easily modifiable screening assay that can measure the TdT activity in high-throughput to evaluate large TdT mutant libraries. We demonstrated the assay's capabilities by engineering TdT mutants that exhibit both improved catalytic efficiency and improved activity in the presence of an inhibitor. We screened for and identified TdT variants with greater catalytic efficiency in both selectively incorporating deoxyribonucleotides and in the presence of deoxyribonucleotide/ribonucleotide mixes. Using this information from the screening assay, we rationally engineered other TdT homologues with the same properties. The emulsion-based assay we developed is, to the best of our knowledge, the first high-throughput screening assay that can measure TdT activity quantitatively and without the need for protein purification.

A massively parallel in vivo assay of TdT mutants yields variants with altered nucleotide insertion biases.

Terminal deoxynucleotidyl transferase (TdT) is a unique DNA polymerase capable of template-independent extension of DNA with random nucleotides. TdT's de novo DNA synthesis ability has found utility in DNA recording, DNA data storage, oligonucleotide synthesis, and nucleic acid labeling, but TdT's intrinsic nucleotide biases limit its versatility in such applications. Here, we describe a multiplexed assay for profiling and engineering the bias and overall activity of TdT variants in high throughput. In our assay, a library of TdTs is encoded next to a CRISPR-Cas9 target site in HEK293T cells. Upon transfection of Cas9 and sgRNA, the target site is cut, allowing TdT to intercept the double strand break and add nucleotides. Each resulting insertion is sequenced alongside the identity of the TdT variant that generated it. Using this assay, 25,623 unique TdT variants, constructed by site-saturation mutagenesis at strategic positions, were profiled. This resulted in the isolation of several altered-bias TdTs that expanded the capabilities of our TdT-based DNA recording system, Cell History Recording by Ordered Insertion (CHYRON), by increasing the information density of recording through an unbiased TdT and achieving dual-channel recording of two distinct inducers (hypoxia and Wnt) through two differently biased TdTs. Select TdT variants were also tested in vitro , revealing concordance between each variant's in vitro bias and the in vivo bias determined from the multiplexed high throughput assay. Overall, our work, and the multiplex assay it features, should support the continued development of TdT-based DNA recorders, in vitro applications of TdT, and further study of the biology of TdT.

Post-translation digital data encoding into the genomes of mammalian cell populations.

High resolution cellular signal encoding is critical for better understanding of complex biological phenomena. DNA-based biosignal encoders alter genomic or plasmid DNA in a signal dependent manner. Current approaches involve the signal of interest affecting a DNA edit by interacting with a signal specific promoter which then results in expression of the effector molecule (DNA altering enzyme). Here, we present the proof of concept of a biosignal encoding system where the enzyme terminal deoxynucleotidyl transferase (TdT) acts as the effector molecule upon directly interacting with the signal of interest. A template independent DNA polymerase (DNAp), TdT incorporates nucleotides at the 3' OH ends of DNA substrate in a signal dependent manner. By employing CRISPR-Cas9 to create double stranded breaks in genomic DNA, we make 3'OH ends available to act as substrate for TdT. We show that this system can successfully resolve and encode different concentrations of various biosignals into the genomic DNA of HEK-293T cells. Finally, we develop a simple encoding scheme associated with the tested biosignals and encode the message "HELLO WORLD" into the genomic DNA of HEK-293T cells at a population level with 91% accuracy. This work demonstrates a simple and engineerable system that can reliably store local biosignal information into the genomes of mammalian cell populations.

Engineering Ca2+-Dependent DNA Polymerase Activity.

Advancements in synthetic biology have provided new opportunities in biosensing, with applications ranging from genetic programming to diagnostics. Next generation biosensors aim to expand the number of accessible environments for measurements, increase the number of measurable phenomena, and improve the quality of the measurement. To this end, an emerging area in the field has been the integration of DNA as an information storage medium within biosensor outputs, leveraging nucleic acids to record the biosensor state over time. However, slow signal transduction steps, due to the time scales of transcription and translation, bottleneck many sensing-DNA recording approaches. DNA polymerases (DNAPs) have been proposed as a solution to the signal transduction problem by operating as both the sensor and responder, but there is presently a lack of DNAPs with functional sensitivity to many desirable target ligands. Here, we engineer components of the Pol δ replicative polymerase complex of Saccharomyces cerevisiae to sense and respond to Ca2+, a metal cofactor relevant to numerous biological phenomena. Through domain insertion and binding site grafting to Pol δ subunits, we demonstrate functional allosteric sensitivity to Ca2+. Together, this work provides an important foundation for future efforts in the development of DNAP-based biosensors.

Aromatic natural products synthesis from aromatic lignin monomers using Acinetobacter baylyi ADP1.

Achieving sustainable chemical synthesis and a circular economy will require process innovation to minimize or recover existing waste streams. Valorization of lignin biomass has the ability to advance this goal. While lignin has proved a recalcitrant feedstock for upgrading, biological approaches can leverage native microbial metabolism to simplify complex and heterogeneous feedstocks to tractable starting points for biochemical upgrading. Recently, we demonstrated that one microbe with lignin relevant metabolism, Acinetobacter baylyi ADP1, is both highly engineerable and capable of undergoing rapid design-build-test-learn cycles, making it an ideal candidate for these applications. Here, we utilize these genetic traits and ADP1's native ß-ketoadipate metabolism to convert mock alkali pretreated liquor lignin (APL) to two valuable natural products, vanillin-glucoside and resveratrol. En route, we create strains with up to 22 genetic modifications, including up to 8 heterologously expressed enzymes. Our approach takes advantage of preexisting aromatic species in APL (vanillate, ferulate, and p-coumarate) to create shortened biochemical routes to end products. Together, this work demonstrates ADP1's potential as a platform for upgrading lignin waste streams and highlights the potential for biosynthetic methods to maximize the existing chemical potential of lignin aromatic monomers.

Pickaxe: a Python library for the prediction of novel metabolic reactions.

BACKGROUND: Biochemical reaction prediction tools leverage enzymatic promiscuity rules to generate reaction networks containing novel compounds and reactions. The resulting reaction networks can be used for multiple applications such as designing novel biosynthetic pathways and annotating untargeted metabolomics data. It is vital for these tools to provide a robust, user-friendly method to generate networks for a given application. However, existing tools lack the flexibility to easily generate networks that are tailor-fit for a user's application due to lack of exhaustive reaction rules, restriction to pre-computed networks, and difficulty in using the software due to lack of documentation. RESULTS: Here we present Pickaxe, an open-source, flexible software that provides a user-friendly method to generate novel reaction networks. This software iteratively applies reaction rules to a set of metabolites to generate novel reactions. Users can select rules from the prepackaged JN1224min ruleset, derived from MetaCyc, or define their own custom rules. Additionally, filters are provided which allow for the pruning of a network on-the-fly based on compound and reaction properties. The filters include chemical similarity to target molecules, metabolomics, thermodynamics, and reaction feasibility filters. Example applications are given to highlight the capabilities of Pickaxe: the expansion of common biological databases with novel reactions, the generation of industrially useful chemicals from a yeast metabolome database, and the annotation of untargeted metabolomics peaks from an E. coli dataset. CONCLUSION: Pickaxe predicts novel metabolic reactions and compounds, which can be used for a variety of applications. This software is open-source and available as part of the MINE Database python package ( https://pypi.org/project/minedatabase/ ) or on GitHub ( https://github.com/tyo-nu/MINE-Database ). Documentation and examples can be found on Read the Docs ( https://mine-database.readthedocs.io/en/latest/ ). Through its documentation, pre-packaged features, and customizable nature, Pickaxe allows users to generate novel reaction networks tailored to their application.

A dynamic kinetic model captures cell-free metabolism for improved butanol production.

Cell-free systems are useful tools for prototyping metabolic pathways and optimizing the production of various bioproducts. Mechanistically-based kinetic models are uniquely suited to analyze dynamic experimental data collected from cell-free systems and provide vital qualitative insight. However, to date, dynamic kinetic models have not been applied with rigorous biological constraints or trained on adequate experimental data to the degree that they would give high confidence in predictions and broadly demonstrate the potential for widespread use of such kinetic models. In this work, we construct a large-scale dynamic model of cell-free metabolism with the goal of understanding and optimizing butanol production in a cell-free system. Using a combination of parameterization methods, the resultant model captures experimental metabolite measurements across two experimental conditions for nine metabolites at timepoints between 0 and 24 h. We present analysis of the model predictions, provide recommendations for butanol optimization, and identify the aldehyde/alcohol dehydrogenase as the primary bottleneck in butanol production. Sensitivity analysis further reveals the extent to which various parameters are constrained, and our approach for probing valid parameter ranges can be applied to other modeling efforts.

MINE 2.0: enhanced biochemical coverage for peak identification in untargeted metabolomics.

SUMMARY: Although advances in untargeted metabolomics have made it possible to gather data on thousands of cellular metabolites in parallel, identification of novel metabolites from these datasets remains challenging. To address this need, Metabolic in silico Network Expansions (MINEs) were developed. A MINE is an expansion of known biochemistry which can be used as a list of potential structures for unannotated metabolomics peaks. Here, we present MINE 2.0, which utilizes a new set of biochemical transformation rules that covers 93% of MetaCyc reactions (compared to 25% in MINE 1.0). This results in a 17-fold increase in database size and a 40% increase in MINE database compounds matching unannotated peaks from an untargeted metabolomics dataset. MINE 2.0 is thus a significant improvement to this community resource. AVAILABILITY AND IMPLEMENTATION: The MINE 2.0 website can be accessed at https://minedatabase.ci.northwestern.edu. The MINE 2.0 web API documentation can be accessed at https://mine-api.readthedocs.io/en/latest/. The data and code underlying this article are available in the MINE-2.0-Paper repository at https://github.com/tyo-nu/MINE-2.0-Paper. MINE 2.0 source code can be accessed at https://github.com/tyo-nu/MINE-Database (MINE construction), https://github.com/tyo-nu/MINE-Server (backend web API) and https://github.com/tyo-nu/MINE-app (web app). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Enabling commercial success of industrial biotechnology.

Commercial-scale research translation has been muted.

Recording Temporal Signals with Minutes Resolution Using Enzymatic DNA Synthesis.

Employing DNA as a high-density data storage medium has paved the way for next-generation digital storage and biosensing technologies. However, the multipart architecture of current DNA-based recording techniques renders them inherently slow and incapable of recording fluctuating signals with subhour frequencies. To address this limitation, we developed a simplified system employing a single enzyme, terminal deoxynucleotidyl transferase (TdT), to transduce environmental signals into DNA. TdT adds nucleotides to the 3'-ends of single-stranded DNA (ssDNA) in a template-independent manner, selecting bases according to inherent preferences and environmental conditions. By characterizing TdT nucleotide selectivity under different conditions, we show that TdT can encode various physiologically relevant signals such as Co2+, Ca2+, and Zn2+ concentrations and temperature changes in vitro. Further, by considering the average rate of nucleotide incorporation, we show that the resulting ssDNA functions as a molecular ticker tape. With this method we accurately encode a temporal record of fluctuations in Co2+ concentration to within 1 min over a 60 min period. Finally, we engineer TdT to allosterically turn off in the presence of a physiologically relevant concentration of calcium. We use this engineered TdT in concert with a reference TdT to develop a two-polymerase system capable of recording a single-step change in the Ca2+ signal to within 1 min over a 60 min period. This work expands the repertoire of DNA-based recording techniques by developing a novel DNA synthesis-based system that can record temporal environmental signals into DNA with a resolution of minutes.

Dynamic Control of Gene Expression with Riboregulated Switchable Feedback Promoters.

One major challenge in synthetic biology is the deleterious impacts of cellular stress caused by expression of heterologous pathways, sensors, and circuits. Feedback control and dynamic regulation are broadly proposed strategies to mitigate this cellular stress by optimizing gene expression levels temporally and in response to biological cues. While a variety of approaches for feedback implementation exist, they are often complex and cannot be easily manipulated. Here, we report a strategy that uses RNA transcriptional regulators to integrate additional layers of control over the output of natural and engineered feedback responsive circuits. Called riboregulated switchable feedback promoters (rSFPs), these gene expression cassettes can be modularly activated using multiple mechanisms, from manual induction to autonomous quorum sensing, allowing control over the timing, magnitude, and autonomy of expression. We develop rSFPs in Escherichia coli to regulate multiple feedback networks and apply them to control the output of two metabolic pathways. We envision that rSFPs will become a valuable tool for flexible and dynamic control of gene expression in metabolic engineering, biological therapeutic production, and many other applications.

Curating a comprehensive set of enzymatic reaction rules for efficient novel biosynthetic pathway design.

Enzyme substrate promiscuity has significant implications for metabolic engineering. The ability to predict the space of possible enzymatic side reactions is crucial for elucidating underground metabolic networks in microorganisms, as well as harnessing novel biosynthetic capabilities of enzymes to produce desired chemicals. Reaction rule-based cheminformatics platforms have been implemented to computationally enumerate possible promiscuous reactions, relying on existing knowledge of enzymatic transformations to inform novel reactions. However, past versions of curated reaction rules have been limited by a lack of comprehensiveness in representing all possible transformations, as well as the need to prune rules to enhance computational efficiency in pathway expansion. To this end, we curated a set of 1224 most generalized reaction rules, automatically abstracted from atom-mapped MetaCyc reactions and verified to uniquely cover all common enzymatic transformations. We developed a framework to systematically identify and correct redundancies and errors in the curation process, resulting in a minimal, yet comprehensive, rule set. These reaction rules were capable of reproducing more than 85% of all reactions in the KEGG and BRENDA databases, for which a large fraction of reactions is not present in MetaCyc. Our rules exceed all previously published rule sets for which reproduction was possible in this coverage analysis, which allows for the exploration of a larger space of known enzymatic transformations. By leveraging the entire knowledge of possible metabolic reactions through generalized enzymatic reaction rules, we are able to better utilize underground metabolic pathways and accelerate novel biosynthetic pathway design to enable bioproduction towards a wider range of new molecules.

Development of a genetic toolset for the highly engineerable and metabolically versatile Acinetobacter baylyi ADP1.

One primary objective of synthetic biology is to improve the sustainability of chemical manufacturing. Naturally occurring biological systems can utilize a variety of carbon sources, including waste streams that pose challenges to traditional chemical processing, such as lignin biomass, providing opportunity for remediation and valorization of these materials. Success, however, depends on identifying micro-organisms that are both metabolically versatile and engineerable. Identifying organisms with this combination of traits has been a historic hindrance. Here, we leverage the facile genetics of the metabolically versatile bacterium Acinetobacter baylyi ADP1 to create easy and rapid molecular cloning workflows, including a Cas9-based single-step marker-less and scar-less genomic integration method. In addition, we create a promoter library, ribosomal binding site (RBS) variants and test an unprecedented number of rationally integrated bacterial chromosomal protein expression sites and variants. At last, we demonstrate the utility of these tools by examining ADP1's catabolic repression regulation, creating a strain with improved potential for lignin bioprocessing. Taken together, this work highlights ADP1 as an ideal host for a variety of sustainability and synthetic biology applications.

Bayesian inference of metabolic kinetics from genome-scale multiomics data.

Modern biological tools generate a wealth of data on metabolite and protein concentrations that can be used to help inform new strain designs. However, learning from these data to predict how a cell will respond to genetic changes, a key need for engineering, remains challenging. A promising technique for leveraging omics measurements in metabolic modeling involves the construction of kinetic descriptions of the enzymatic reactions that occur within a cell. Parameterizing these models from biological data can be computationally difficult, since methods must also quantify the uncertainty in model parameters resulting from the observed data. While the field of Bayesian inference offers a wide range of methods for efficiently estimating distributions in parameter uncertainty, such techniques are poorly suited to traditional kinetic models due to their complex rate laws and resulting nonlinear dynamics. In this paper, we employ linear-logarithmic kinetics to simplify the calculation of steady-state flux distributions and enable efficient sampling and inference methods. We demonstrate that detailed information on the posterior distribution of parameters can be obtained efficiently at a variety of problem scales, including nearly genome-scale kinetic models trained on multiomics datasets. These results allow modern Bayesian machine learning tools to be leveraged in understanding biological data and in developing new, efficient strain designs.

Model-guided mechanism discovery and parameter selection for directed evolution.

Directed evolution is frequently applied to identify genetic variants with improvements in a single or multiple properties. When used to improve multiple properties simultaneously, a common strategy is to apply alternating rounds of selection criteria to enrich for variants with each desirable trait. In particular, counterselection, or selection against undesired traits rather than for desired ones, has been successfully employed in many studies. Although the sequence and stringency of alternating selective pressures for different traits are known to be highly consequential for the outcome of the screen, the effects of these parameters have not been systematically evaluated. We developed a method for producing a statistical modeling framework to elucidate these effects. The model uses single-cell fluorescence intensity distributions to estimate the proportions of phenotypic populations within a library and then predicts the changes in these proportions depending on specified positive selective or counterselective pressures. We validated the approach using recently described systems for metabolite-responsive bacterial transcription factors and yeast G-protein-coupled receptors. Finally, we applied the model to identify biological sources that exert undesirable selective pressure on libraries during sorting. Notably, these pressures produce substantial artifacts that, if unaddressed, can lead to failure of the screen. This method for model generation can be applied to FACS-based directed evolution experiments to create a quantitative framework that identifies subtle population effects. Such models can guide the choice of experimental design parameters to better enrich for true positive genetic variants and improve the chance of successful directed evolution.

Metabolic In Silico Network Expansions to Predict and Exploit Enzyme Promiscuity.

There is a growing consensus that enzymes are capable of catalyzing not just one canonical reaction but entire families of related reactions. These capacities often go unnoticed in the enzyme's native context but can become apparent in engineered metabolism when the enzyme is exposed to novel substrates or high concentrations of pathway intermediates. This chapter describes how to use metabolic in silico network expansion (MINE) databases to predict novel biotransformations and their resulting metabolites. In particular, searching MINEs by structural similarity or with metabolomics data allows scientists to detect, exploit, or avoid these predicted transformations.

Chemically Inducible Chromosomal Evolution (CIChE) for Multicopy Metabolic Pathway Engineering.

Chemically inducible chromosomal evolution (CIChE) was developed for stable multicopy chromosomal integration of heterologous genes. In this technique, flanking an antibiotic selection marker and a gene of interest with identical regions of homology permits gene duplication via recA mediated homologous recombination. A strong selective pressure for gene duplication can be applied by increasing antibiotic concentration, and in a week's time one can create a set of strains with a wide range of cassette copy numbers (upward of 20×), which can be made stable by deletion of recA. Herein, we describe a generalized workflow for this methodology.

High-resolution mapping of DNA polymerase fidelity using nucleotide imbalances and next-generation sequencing.

DNA polymerase fidelity is affected by both intrinsic properties and environmental conditions. Current strategies for measuring DNA polymerase error rate in vitro are constrained by low error subtype sensitivity, poor scalability, and lack of flexibility in types of sequence contexts that can be tested. We have developed the Magnification via Nucleotide Imbalance Fidelity (MagNIFi) assay, a scalable next-generation sequencing assay that uses a biased deoxynucleotide pool to quantitatively shift error rates into a range where errors are frequent and hence measurement is robust, while still allowing for accurate mapping to error rates under typical conditions. This assay is compatible with a wide range of fidelity-modulating conditions, and enables high-throughput analysis of sequence context effects on base substitution and single nucleotide deletion fidelity using a built-in template library. We validate this assay by comparing to previously established fidelity metrics, and use it to investigate neighboring sequence-mediated effects on fidelity for several DNA polymerases. Through these demonstrations, we establish the MagNIFi assay for robust, high-throughput analysis of DNA polymerase fidelity.