Engineering of Pseudomonas putida for accelerated co-utilization of glucose and cellobiose yields aerobic overproduction of pyruvate explained by an upgraded metabolic model

A versatile microbial platform as a tunable whole-cell chemical sensor

Biosensors are used to detect and quantify chemicals produced in industrial microbiology with high specificity, sensitivity, and portability. Most biosensors, however, are limited by the need for transcription factors engineered to recognize specific molecules. In this study, we overcome the limitations typically associated with traditional biosensors by engineering Pseudomonas putida for whole-cell sensing of a variety of chemicals. Our approach integrates fluorescent reporters with synthetic auxotrophies within central metabolism that can be complemented by target analytes in growth-coupled setups. This platform enables the detection of a wide array of structurally diverse chemicals under various conditions, including co-cultures of producer cell factories and sensor strains. We also demonstrate the applicability of this versatile biosensor platform for monitoring complex biochemical processes, including plastic degradation by either purified hydrolytic enzymes or engineered bacteria. This microbial system provides a rapid, sensitive, and readily adaptable tool for monitoring cell factory performance and for environmental analyzes.

Systems biology of electrogenic Pseudomonas putida - multi-omics insights and metabolic engineering for enhanced 2-ketogluconate production

BACKGROUND

Pseudomonas putida KT2440 has emerged as a promising host for industrial bioproduction. However, its strictly aerobic nature limits the scope of applications. Remarkably, this microbe exhibits high bioconversion efficiency when cultured in an anoxic bio-electrochemical system (BES), where the anode serves as the terminal electron acceptor instead of oxygen. This environment facilitates the synthesis of commercially attractive chemicals, including 2-ketogluconate (2KG). To better understand this interesting electrogenic phenotype, we studied the BES-cultured strain on a systems level through multi-omics analysis. Inspired by our findings, we constructed novel mutants aimed at improving 2KG production.

RESULTS

When incubated on glucose, P. putida KT2440 did not grow but produced significant amounts of 2KG, along with minor amounts of gluconate, acetate, pyruvate, succinate, and lactate. 13C tracer studies demonstrated that these products are partially derived from biomass carbon, involving proteins and lipids. Over time, the cells exhibited global changes on both the transcriptomic and proteomic levels, including the shutdown of translation and cell motility, likely to conserve energy. These adaptations enabled the cells to maintain significant metabolic activity for several weeks. Acetate formation was shown to contribute to energy supply. Mutants deficient in acetate production demonstrated superior 2KG production in terms of titer, yield, and productivity. The ∆aldBI ∆aldBII double deletion mutant performed best, accumulating 2KG at twice the rate of the wild type and with an increased yield (0.96 mol/mol).

CONCLUSIONS

By integrating transcriptomic, proteomic, and metabolomic analyses, this work provides the first systems biology insight into the electrogenic phenotype of P. putida KT2440. Adaptation to anoxic-electrogenic conditions involved coordinated changes in energy metabolism, enabling cells to sustain metabolic activity for extended periods. The metabolically engineered mutants are promising for enhanced 2KG production under these conditions. The attenuation of acetate synthesis represents the first systems biology-informed metabolic engineering strategy for enhanced 2KG production in P. putida. This non-growth anoxic-electrogenic mode expands our understanding of the interplay between growth, glucose phosphorylation, and glucose oxidation into gluconate and 2KG in P. putida.

Synthetically-primed adaptation of Pseudomonas putida to a non-native substrate D-xylose

To broaden the substrate scope of microbial cell factories towards renewable substrates, rational genetic interventions are often combined with adaptive laboratory evolution (ALE). However, comprehensive studies enabling a holistic understanding of adaptation processes primed by rational metabolic engineering remain scarce. The industrial workhorse Pseudomonas putida was engineered to utilize the non-native sugar D-xylose, but its assimilation into the bacterial biochemical network via the exogenous xylose isomerase pathway remained unresolved. Here, we elucidate the xylose metabolism and establish a foundation for further engineering followed by ALE. First, native glycolysis is derepressed by deleting the local transcriptional regulator gene hexR. We then enhance the pentose phosphate pathway by implanting exogenous transketolase and transaldolase into two lag-shortened strains and allow ALE to finetune the rewired metabolism. Subsequent multilevel analysis and reverse engineering provide detailed insights into the parallel paths of bacterial adaptation to the non-native carbon source, highlighting the enhanced expression of transaldolase and xylose isomerase along with derepressed glycolysis as key events during the process.

Machine learning analysis of RB-TnSeq fitness data predicts functional gene modules in Pseudomonas putida KT2440

There is growing interest in engineering Pseudomonas putida KT2440 as a microbial chassis for the conversion of renewable and waste-based feedstocks, and metabolic engineering of P. putida relies on the understanding of the functional relationships between genes. In this work, independent component analysis (ICA) was applied to a compendium of existing fitness data from randomly barcoded transposon insertion sequencing (RB-TnSeq) of P. putida KT2440 grown in 179 unique experimental conditions. ICA identified 84 independent groups of genes, which we call fModules ("functional modules"), where gene members displayed shared functional influence in a specific cellular process. This machine learning-based approach both successfully recapitulated previously characterized functional relationships and established hitherto unknown associations between genes. Selected gene members from fModules for hydroxycinnamate metabolism and stress resistance, acetyl coenzyme A assimilation, and nitrogen metabolism were validated with engineered mutants of P. putida. Additionally, functional gene clusters from ICA of RB-TnSeq data sets were compared with regulatory gene clusters from prior ICA of RNAseq data sets to draw connections between gene regulation and function. Because ICA profiles the functional role of several distinct gene networks simultaneously, it can reduce the time required to annotate gene function relative to manual curation of RB-TnSeq data sets.

IMPORTANCE

This study demonstrates a rapid, automated approach for elucidating functional modules within complex genetic networks. While Pseudomonas putida randomly barcoded transposon insertion sequencing data were used as a proof of concept, this approach is applicable to any organism with existing functional genomics data sets and may serve as a useful tool for many valuable applications, such as guiding metabolic engineering efforts in other microbes or understanding functional relationships between virulence-associated genes in pathogenic microbes. Furthermore, this work demonstrates that comparison of data obtained from independent component analysis of transcriptomics and gene fitness datasets can elucidate regulatory-functional relationships between genes, which may have utility in a variety of applications, such as metabolic modeling, strain engineering, or identification of antimicrobial drug targets.

Genome-scale and pathway engineering for the sustainable aviation fuel precursor isoprenol production in Pseudomonas putida

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.

Reconstruction, simulation and analysis of enzyme-constrained metabolic models using GECKO Toolbox 3.0

Genome-scale metabolic models (GEMs) are computational representations that enable mathematical exploration of metabolic behaviors within cellular and environmental constraints. Despite their wide usage in biotechnology, biomedicine and fundamental studies, there are many phenotypes that GEMs are unable to correctly predict. GECKO is a method to improve the predictive power of a GEM by incorporating enzymatic constraints using kinetic and omics data. GECKO has enabled reconstruction of enzyme-constrained metabolic models (ecModels) for diverse organisms, which show better predictive performance than conventional GEMs. In this protocol, we describe how to use the latest version GECKO 3.0; the procedure has five stages: (1) expansion from a starting metabolic model to an ecModel structure, (2) integration of enzyme turnover numbers into the ecModel structure, (3) model tuning, (4) integration of proteomics data into the ecModel and (5) simulation and analysis of ecModels. GECKO 3.0 incorporates deep learning-predicted enzyme kinetics, paving the way for improved metabolic models for virtually any organism and cell line in the absence of experimental data. The time of running the whole protocol is organism dependent, e.g., ~5 h for yeast.

Evaluation of enzyme-constrained genome-scale model through metabolic engineering of anaerobic co-production of 2,3-butanediol and glycerol by Saccharomyces cerevisiae

Enzyme-constrained genome-scale models (ecGEMs) have potential to predict phenotypes in a variety of conditions, such as growth rates or carbon sources. This study investigated if ecGEMs can guide metabolic engineering efforts to swap anaerobic redox-neutral ATP-providing pathways in yeast from alcoholic fermentation to equimolar co-production of 2,3-butanediol and glycerol. With proven pathways and low product toxicity, the ecGEM solution space aligned well with observed phenotypes. Since this catabolic pathway provides only one-third of the ATP of alcoholic fermentation (2/3 versus 2 ATP per glucose), the ecGEM predicted a growth decrease from 0.36 h-1 in the reference to 0.175 h-1 in the engineered strain. However, this <3-fold decrease would require the specific glucose consumption rate to increase. Surprisingly, after the pathway swap the engineered strain immediately grew at 0.15 h-1 with a glucose consumption rate of 29 mmol (g CDW)-1 h-1, which was indeed higher than reference (23 mmol (g CDW)-1 h-1) and one of the highest reported for S. cerevisiae. The accompanying 2,3-butanediol- (15.8 mmol (g CDW)-1 h-1) and glycerol (19.6 mmol (g CDW)-1 h-1) production rates were close to predicted values. Proteomics confirmed that this increased consumption rate was facilitated by enzyme reallocation from especially ribosomes (from 25.5 to 18.5 %) towards glycolysis (from 28.7 to 43.5 %). Subsequently, 200 generations of sequential transfer did not improve growth of the engineered strain, showing the use of ecGEMs in predicting opportunity space for laboratory evolution. The observations in this study illustrate both the current potential, as well as future improvements, of ecGEMs as a tool for both metabolic engineering and laboratory evolution.

Pseudomonas putida as a synthetic biology chassis and a metabolic engineering platform

The soil bacterium Pseudomonas putida, especially the KT2440 strain, is increasingly being utilized as a host for biotransformations of both industrial and environmental interest. The foundations of such performance include its robust redox metabolism, ability to tolerate a wide range of physicochemical stresses, rapid growth, versatile metabolism, nonpathogenic nature, and the availability of molecular tools for advanced genetic programming. These attributes have been leveraged for hosting engineered pathways for production of valuable chemicals or degradation/valorization of environmental pollutants. This has in turn pushed the boundaries of conventional enzymology toward previously unexplored reactions in nature. Furthermore, modifications to the physical properties of the cells have been made to enhance their catalytic performance. These advancements establish P. putida as bona fide chassis for synthetic biology, on par with more traditional metabolic engineering platforms.

Metabolic impact of heterologous protein production in Pseudomonas putida: Insights into carbon and energy flux control

For engineered microorganisms, the production of heterologous proteins that are often useless to host cells represents a burden on resources, which have to be shared with normal cellular processes. Within a certain metabolic leeway, this competitive process has no impact on growth. However, once this leeway, or free capacity, is fully utilized, the extra load becomes a metabolic burden that inhibits cellular processes and triggers a broad cellular response, reducing cell growth and often hindering the production of heterologous proteins. In this study, we sought to characterize the metabolic rearrangements occurring in the central metabolism of Pseudomonas putida at different levels of metabolic load. To this end, we constructed a P. putida KT2440 strain that expressed two genes encoding fluorescent proteins, one in the genome under constitutive expression to monitor the free capacity, and the other on an inducible plasmid to probe heterologous protein production. We found that metabolic fluxes are considerably reshuffled, especially at the level of periplasmic pathways, as soon as the metabolic load exceeds the free capacity. Heterologous protein production leads to the decoupling of anabolism and catabolism, resulting in large excess energy production relative to the requirements of protein biosynthesis. Finally, heterologous protein production was found to exert a stronger control on carbon fluxes than on energy fluxes, indicating that the flexible nature of P. putida's central metabolic network is solicited to sustain energy production.

Improving microbial bioproduction under low-oxygen conditions

Microbial bioconversion provides access to a wide range of sustainably produced chemicals and commodities. However, industrial-scale bioproduction process operations are preferred to be anaerobic due to the cost associated with oxygen transfer. Anaerobic bioconversion generally offers limited substrate utilization profiles, lower product yields, and reduced final product diversity compared with aerobic processes. Bioproduction under conditions of reduced oxygen can overcome the limitations of fully aerobic and anaerobic bioprocesses, but many microbial hosts are not developed for low-oxygen bioproduction. Here, we describe advances in microbial strain engineering involving the use of redox cofactor engineering, genome-scale metabolic modeling, and functional genomics to enable improved bioproduction processes under low oxygen and provide a viable path for scaling these bioproduction systems to industrial scales.

Time-resolved, deuterium-based fluxomics uncovers the hierarchy and dynamics of sugar processing by Pseudomonas putida

Pseudomonas putida, a microbial host widely adopted for metabolic engineering, processes glucose through convergent peripheral pathways that ultimately yield 6-phosphogluconate. The periplasmic gluconate shunt (PGS), composed by glucose and gluconate dehydrogenases, sequentially transforms glucose into gluconate and 2-ketogluconate. Although the secretion of these organic acids by P. putida has been extensively recognized, the mechanism and spatiotemporal regulation of the PGS remained elusive thus far. To address this challenge, we adopted a dynamic 13C- and 2H-metabolic flux analysis strategy, termed D-fluxomics. D-fluxomics demonstrated that the PGS underscores a highly dynamic metabolic architecture in glucose-dependent batch cultures of P. putida, characterized by hierarchical carbon uptake by the PGS throughout the cultivation. Additionally, we show that gluconate and 2-ketogluconate accumulation and consumption can be solely explained as a result of the interplay between growth rate-coupled and decoupled metabolic fluxes. As a consequence, the formation of these acids in the PGS is inversely correlated to the bacterial growth rate-unlike the widely studied overflow metabolism of Escherichia coli and yeast. Our findings, which underline survival strategies of soil bacteria thriving in their natural environments, open new avenues for engineering P. putida towards efficient, sugar-based bioprocesses.