Search

Sustainable aviation fuel (SAF) will significantly impact global warming in the aviation sector, and important SAF targets are emerging. Isoprenol is a precursor for a promising SAF compound DMCO (1,4-dimethylcyclooctane) and has been produced in several engineered microorganisms. Recently, Pseudomonas putida has gained interest as a future host for isoprenol bioproduction as it can utilize carbon sources from inexpensive plant biomass. Here, we engineer metabolically versatile host P. putida for isoprenol production. We employ two computational modeling approaches (Bilevel optimization and Constrained Minimal Cut Sets) to predict gene knockout targets and optimize the "IPP-bypass" pathway in P. putida to maximize isoprenol production. Altogether, the highest isoprenol production titer from P. putida was achieved at 3.5 g/L under fed-batch conditions. This combination of computational modeling and strain engineering on P. putida for an advanced biofuels production has vital significance in enabling a bioproduction process that can use renewable carbon streams.

Metabolic fluxes, the number of metabolites traversing each biochemical reaction in a cell per unit time, are crucial for assessing and understanding cell function. 13C Metabolic Flux Analysis (13C MFA) is considered to be the gold standard for measuring metabolic fluxes. 13C MFA typically works by leveraging extracellular exchange fluxes as well as data from 13C labeling experiments to calculate the flux profile which best fit the data for a small, central carbon, metabolic model. However, the nonlinear nature of the 13C MFA fitting procedure means that several flux profiles fit the experimental data within the experimental error, and traditional optimization methods offer only a partial or skewed picture, especially in "non-gaussian" situations where multiple very distinct flux regions fit the data equally well. Here, we present a method for flux space sampling through Bayesian inference (BayFlux), that identifies the full distribution of fluxes compatible with experimental data for a comprehensive genome-scale model. This Bayesian approach allows us to accurately quantify uncertainty in calculated fluxes. We also find that, surprisingly, the genome-scale model of metabolism produces narrower flux distributions (reduced uncertainty) than the small core metabolic models traditionally used in 13C MFA. The different results for some reactions when using genome-scale models vs core metabolic models advise caution in assuming strong inferences from 13C MFA since the results may depend significantly on the completeness of the model used. Based on BayFlux, we developed and evaluated novel methods (P-13C MOMA and P-13C ROOM) to predict the biological results of a gene knockout, that improve on the traditional MOMA and ROOM methods by quantifying prediction uncertainty.

Biology has changed radically in the past two decades, growing from a purely descriptive science into also a design science. The availability of tools that enable the precise modification of cells, as well as the ability to collect large amounts of multimodal data, open the possibility of sophisticated bioengineering to produce fuels, specialty and commodity chemicals, materials, and other renewable bioproducts. However, despite new tools and exponentially increasing data volumes, synthetic biology cannot yet fulfill its true potential due to our inability to predict the behavior of biological systems. Here, we showcase a set of computational tools that, combined, provide the ability to store, visualize, and leverage multiomics data to predict the outcome of bioengineering efforts. We show how to upload, visualize, and output multiomics data, as well as strain information, into online repositories for several isoprenol-producing strain designs. We then use these data to train machine learning algorithms that recommend new strain designs that are correctly predicted to improve isoprenol production by 23%. This demonstration is done by using synthetic data, as provided by a novel library, that can produce credible multiomics data for testing algorithms and computational tools. In short, this paper provides a step-by-step tutorial to leverage these computational tools to improve production in bioengineered strains.

Caprolactam is an important polymer precursor to nylon traditionally derived from petroleum and produced on a scale of 5 million tons per year. Current biological pathways for the production of caprolactam are inefficient with titers not exceeding 2 mg/L, necessitating novel pathways for its production. As development of novel metabolic routes often require thousands of designs and result in low product titers, a highly sensitive biosensor for the final product has the potential to rapidly speed up development times. Here we report a highly sensitive biosensor for valerolactam and caprolactam from Pseudomonas putida KT2440 which is >1000× more sensitive to an exogenous ligand than previously reported sensors. Manipulating the expression of the sensor oplR (PP_3516) substantially altered the sensing parameters, with various vectors showing Kd values ranging from 700 nM (79.1 μg/L) to 1.2 mM (135.6 mg/L). Our most sensitive construct was able to detect in vivo production of caprolactam above background at ∼6 μg/L. The high sensitivity and range of OplR is a powerful tool toward the development of novel routes to the biological synthesis of caprolactam.

Through advanced mechanistic modeling and the generation of large high-quality datasets, machine learning is becoming an integral part of understanding and engineering living systems. Here we show that mechanistic and machine learning models can be combined to enable accurate genotype-to-phenotype predictions. We use a genome-scale model to pinpoint engineering targets, efficient library construction of metabolic pathway designs, and high-throughput biosensor-enabled screening for training diverse machine learning algorithms. From a single data-generation cycle, this enables successful forward engineering of complex aromatic amino acid metabolism in yeast, with the best machine learning-guided design recommendations improving tryptophan titer and productivity by up to 74 and 43%, respectively, compared to the best designs used for algorithm training. Thus, this study highlights the power of combining mechanistic and machine learning models to effectively direct metabolic engineering efforts.

The Design-Build-Test-Learn (DBTL) cycle, facilitated by exponentially improving capabilities in synthetic biology, is an increasingly adopted metabolic engineering framework that represents a more systematic and efficient approach to strain development than historical efforts in biofuels and biobased products. Here, we report on implementation of two DBTL cycles to optimize 1-dodecanol production from glucose using 60 engineered Escherichia coli MG1655 strains. The first DBTL cycle employed a simple strategy to learn efficiently from a relatively small number of strains (36), wherein only the choice of ribosome-binding sites and an acyl-ACP/acyl-CoA reductase were modulated in a single pathway operon including genes encoding a thioesterase (UcFatB1), an acyl-ACP/acyl-CoA reductase (Maqu_2507, Maqu_2220, or Acr1), and an acyl-CoA synthetase (FadD). Measured variables included concentrations of dodecanol and all proteins in the engineered pathway. We used the data produced in the first DBTL cycle to train several machine-learning algorithms and to suggest protein profiles for the second DBTL cycle that would increase production. These strategies resulted in a 21% increase in dodecanol titer in Cycle 2 (up to 0.83 g/L, which is more than 6-fold greater than previously reported batch values for minimal medium). Beyond specific lessons learned about optimizing dodecanol titer in E. coli, this study had findings of broader relevance across synthetic biology applications, such as the importance of sequencing checks on plasmids in production strains as well as in cloning strains, and the critical need for more accurate protein expression predictive tools.

Metabolic models can estimate intrinsic product yields for microbial factories, but such frameworks struggle to predict cell performance (including product titer or rate) under suboptimal metabolism and complex bioprocess conditions. On the other hand, machine learning, complementary to metabolic modeling necessitates large amounts of data. Building such a database for metabolic engineering designs requires significant manpower and is prone to human errors and bias. We propose an approach to integrate data-driven methods with genome scale metabolic model for assessment of microbial bio-production (yield, titer and rate). Using engineered E. coli as an example, we manually extracted and curated a data set comprising about 1200 experimentally realized cell factories from ~100 papers. We furthermore augmented the key design features (e.g., genetic modifications and bioprocess variables) extracted from literature with additional features derived from running the genome-scale metabolic model iML1515 simulations with constraints that match the experimental data. Then, data augmentation and ensemble learning (e.g., support vector machines, gradient boosted trees, and neural networks in a stacked regressor model) are employed to alleviate the challenges of sparse, non-standardized, and incomplete data sets, while multiple correspondence analysis/principal component analysis are used to rank influential factors on bio-production. The hybrid framework demonstrates a reasonably high cross-validation accuracy for prediction of E.coli factory performance metrics under presumed bioprocess and pathway conditions (Pearson correlation coefficients between 0.8 and 0.93 on new data not seen by the model).

Robust fermentation of biomass-derived sugars into bioproducts demands the reliable microbial expression of metabolic pathways. Plasmid-based expression systems may suffer from instability and result in highly variable titers, rates, and yields. An established mitigation approach, chemical induced chromosomal expansion (CIChE), expands a singly integrated pathway to plasmid-like copy numbers while maintaining stability in the absence of antibiotic selection pressure. Here, we report parallel integration and chromosomal expansion (PIACE), extensions to CIChE that enable independent expansions of pathway components across multiple loci, use suicide vectors to achieve high-efficiency site-specific integration of sequence-validated multigene components, and introduce a heat-curable plasmid to obviate recA deletion post pathway expansion. We applied PIACE to stabilize an isopentenol pathway across three loci in E. coli DH1 and then generate libraries of pathway component copy number variants to screen for improved titers. Polynomial regressor statistical modeling of the production screening data suggests that increasing copy numbers of all isopentenol pathway components would further improve titers.

Genome scale modeling (GSM) predicts the performance of microbial workhorses and helps identify beneficial gene targets. GSM integrated with intracellular flux dynamics, omics, and thermodynamics have shown remarkable progress in both elucidating complex cellular phenomena and computational strain design (CSD). Nonetheless, these models still show high uncertainty due to a poor understanding of innate pathway regulations, metabolic burdens, and other factors (such as stress tolerance and metabolite channeling). Besides, the engineered hosts may have genetic mutations or non-genetic variations in bioreactor conditions and thus CSD rarely foresees fermentation rate and titer. Metabolic models play important role in design-build-test-learn cycles for strain improvement, and machine learning (ML) may provide a viable complementary approach for driving strain design and deciphering cellular processes. In order to develop quality ML models, knowledge engineering leverages and standardizes the wealth of information in literature (e.g., genomic/phenomic data, synthetic biology strategies, and bioprocess variables). Data driven frameworks can offer new constraints for mechanistic models to describe cellular regulations, to design pathways, to search gene targets, and to estimate fermentation titer/rate/yield under specified growth conditions (e.g., mixing, nutrients, and O2). This review highlights the scope of information collections, database constructions, and machine learning techniques (such as deep learning and transfer learning), which may facilitate "Learn and Design" for strain development.

Flowers of the hop plant provide both bitterness and "hoppy" flavor to beer. Hops are, however, both a water and energy intensive crop and vary considerably in essential oil content, making it challenging to achieve a consistent hoppy taste in beer. Here, we report that brewer's yeast can be engineered to biosynthesize aromatic monoterpene molecules that impart hoppy flavor to beer by incorporating recombinant DNA derived from yeast, mint, and basil. Whereas metabolic engineering of biosynthetic pathways is commonly enlisted to maximize product titers, tuning expression of pathway enzymes to affect target production levels of multiple commercially important metabolites without major collateral metabolic changes represents a unique challenge. By applying state-of-the-art engineering techniques and a framework to guide iterative improvement, strains are generated with target performance characteristics. Beers produced using these strains are perceived as hoppier than traditionally hopped beers by a sensory panel in a double-blind tasting.

BackgroundGlycolysis breakdowns glucose into essential building blocks and ATP/NAD(P)H for the cell, occupying a central role in its growth and bio-production. Among glycolytic pathways, the Entner Doudoroff pathway (EDP) is a more thermodynamically favorable pathway with fewer enzymatic steps than either the Embden-Meyerhof-Parnas pathway (EMPP) or the oxidative pentose phosphate pathway (OPPP). However, Escherichia coli do not use their native EDP for glucose metabolism.ResultsOverexpression of edd and eda in E. coli to enhance EDP activity resulted in only a small shift in the flux directed through the EDP (~20 % of glycolysis flux). Disrupting the EMPP by phosphofructokinase I (pfkA) knockout increased flux through OPPP (~60 % of glycolysis flux) and the native EDP (~14 % of glycolysis flux), while overexpressing edd and eda in this ΔpfkA mutant directed ~70 % of glycolytic flux through the EDP. The downregulation of EMPP via the pfkA deletion significantly decreased the growth rate, while EDP overexpression in the ΔpfkA mutant failed to improve its growth rates due to metabolic burden. However, the reorganization of E. coli glycolytic strategies did reduce glucose catabolite repression. The ΔpfkA mutant in glucose medium was able to cometabolize acetate via the citric acid cycle and gluconeogenesis, while EDP overexpression in the ΔpfkA mutant repressed acetate flux toward gluconeogenesis. Moreover, 13C-pulse experiments in the ΔpfkA mutants showed unsequential labeling dynamics in glycolysis intermediates, possibly suggesting metabolite channeling (metabolites in glycolysis are pass from enzyme to enzyme without fully equilibrating within the cytosol medium).ConclusionsWe engineered E. coli to redistribute its native glycolytic flux. The replacement of EMPP by EDP did not improve E. coli glucose utilization or biomass growth, but alleviated catabolite repression. More importantly, our results supported the hypothesis of channeling in the glycolytic pathways, a potentially overlooked mechanism for regulating glucose catabolism and coutilization of other substrates. The presence of channeling in native pathways, if proven true, would affect synthetic biology applications and metabolic modeling.

Understanding the complex interactions that occur between heterologous and native biochemical pathways represents a major challenge in metabolic engineering and synthetic biology. We present a workflow that integrates metabolomics, proteomics, and genome-scale models of Escherichia coli metabolism to study the effects of introducing a heterologous pathway into a microbial host. This workflow incorporates complementary approaches from computational systems biology, metabolic engineering, and synthetic biology; provides molecular insight into how the host organism microenvironment changes due to pathway engineering; and demonstrates how biological mechanisms underlying strain variation can be exploited as an engineering strategy to increase product yield. As a proof of concept, we present the analysis of eight engineered strains producing three biofuels: isopentenol, limonene, and bisabolene. Application of this workflow identified the roles of candidate genes, pathways, and biochemical reactions in observed experimental phenomena and facilitated the construction of a mutant strain with improved productivity. The contributed workflow is available as an open-source tool in the form of iPython notebooks.

BackgroundGlycolysis breakdowns glucose into essential building blocks and ATP/NAD(P)H for the cell, occupying a central role in its growth and bio-production. Among glycolytic pathways, the Entner Doudoroff pathway (EDP) is a more thermodynamically favorable pathway with fewer enzymatic steps than either the Embden-Meyerhof-Parnas pathway (EMPP) or the oxidative pentose phosphate pathway (OPPP). However, Escherichia coli do not use their native EDP for glucose metabolism.ResultsOverexpression of edd and eda in E. coli to enhance EDP activity resulted in only a small shift in the flux directed through the EDP (~20 % of glycolysis flux). Disrupting the EMPP by phosphofructokinase I (pfkA) knockout increased flux through OPPP (~60 % of glycolysis flux) and the native EDP (~14 % of glycolysis flux), while overexpressing edd and eda in this ΔpfkA mutant directed ~70 % of glycolytic flux through the EDP. The downregulation of EMPP via the pfkA deletion significantly decreased the growth rate, while EDP overexpression in the ΔpfkA mutant failed to improve its growth rates due to metabolic burden. However, the reorganization of E. coli glycolytic strategies did reduce glucose catabolite repression. The ΔpfkA mutant in glucose medium was able to cometabolize acetate via the citric acid cycle and gluconeogenesis, while EDP overexpression in the ΔpfkA mutant repressed acetate flux toward gluconeogenesis. Moreover, 13C-pulse experiments in the ΔpfkA mutants showed unsequential labeling dynamics in glycolysis intermediates, possibly suggesting metabolite channeling (metabolites in glycolysis are pass from enzyme to enzyme without fully equilibrating within the cytosol medium).ConclusionsWe engineered E. coli to redistribute its native glycolytic flux. The replacement of EMPP by EDP did not improve E. coli glucose utilization or biomass growth, but alleviated catabolite repression. More importantly, our results supported the hypothesis of channeling in the glycolytic pathways, a potentially overlooked mechanism for regulating glucose catabolism and coutilization of other substrates. The presence of channeling in native pathways, if proven true, would affect synthetic biology applications and metabolic modeling.