Systematizing the different notions of growth-coupled product synthesis and a single framework for computing corresponding strain designs

Minimal cut sets in metabolic networks: from conceptual foundations to applications in metabolic engineering and biomedicine

Minimal cut sets (MCSs) have emerged as an important branch of constraint-based metabolic modeling, offering a versatile framework for analyzing and engineering metabolic networks. Over the past two decades, MCSs have evolved from a theoretical concept into a powerful tool for identifying tailored metabolic intervention strategies and studying robustness and failure modes of metabolic networks. Successful (experimental) applications range from designing highly efficient microbial cell factories to targeting cancer cell metabolism. This review highlights key conceptual and algorithmic advancements that have transformed MCSs into a flexible methodology applicable to metabolic models of any size. It also provides a comprehensive overview of their applications and concludes with a perspective on future research directions. The review aims to equip both newcomers and experts with the knowledge needed to effectively leverage MCSs for metabolic network analysis and design, therapeutic targeting, and beyond.

Addressing genome scale design tradeoffs in Pseudomonas putida for bioconversion of an aromatic carbon source

Genome-scale metabolic models (GSMM) are commonly used to identify gene deletion sets that result in growth coupling and pairing product formation with substrate utilization and can improve strain performance beyond levels typically accessible using traditional strain engineering approaches. However, sustainable feedstocks pose a challenge due to incomplete high-resolution metabolic data for non-canonical carbon sources required to curate GSMM and identify implementable designs. Here we address a four-gene deletion design in the Pseudomonas putida KT2440 strain for the lignin-derived non-sugar carbon source, p-coumarate (p-CA), that proved challenging to implement. We examine the performance of the fully implemented design for p-coumarate to glutamine, a useful biomanufacturing intermediate. In this study glutamine is then converted to indigoidine, an alternative sustainable pigment and a model heterologous product that is commonly used to colorimetrically quantify glutamine concentration. Through proteomics, promoter-variation, and growth characterization of a fully implemented gene deletion design, we provide evidence that aromatic catabolism in the completed design is rate-limited by fumarase hydratase (FUM) enzyme activity in the citrate cycle and requires careful optimization of another fumarate hydratase protein (PP_0897) expression to achieve growth and production. A double sensitivity analysis also confirmed a strict requirement for fumarate hydratase activity in the strain where all genes in the growth coupling design have been implemented. Metabolic cross-feeding experiments were used to examine the impact of complete removal of the fumarase hydratase reaction and revealed an unanticipated nutrient requirement, suggesting additional functions for this enzyme. While a complete implementation of the design was achieved, this study highlights the challenge of completely inactivating metabolic reactions encoded by under-characterized proteins, especially in the context of multi-gene edits.

Seven critical challenges in synthetic one-carbon assimilation and their potential solutions

Synthetic C1 assimilation holds the promise of facilitating carbon capture while mitigating greenhouse gas emissions, yet practical implementation in microbial hosts remains relatively limited. Despite substantial progress in pathway design and prototyping, most efforts stay at the proof-of-concept stage, with frequent failures observed even under in vitro conditions. This review identifies seven major barriers constraining the deployment of synthetic C1 metabolism in microorganisms and proposes targeted strategies for overcoming these issues. A primary limitation is the low catalytic activity of carbon-fixing enzymes, particularly carboxylases, which restricts the overall pathway performance. In parallel, challenges in expressing multiple heterologous genes-especially those encoding metal-dependent or oxygen-sensitive enzymes-further hinder pathway functionality. At the systems level, synthetic C1 pathways often exhibit poor flux distribution, limited integration with the host metabolism, accumulation of toxic intermediates, and disruptions in redox and energy balance. These factors collectively reduce biomass formation and compromise product yields in biotechnological setups. Overcoming these interconnected challenges is essential for moving synthetic C1 assimilation beyond conceptual stages and enabling its application in scalable, efficient bioprocesses towards a circular bioeconomy.

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.

FastKnock: an efficient next-generation approach to identify all knockout strategies for strain optimization

Overproduction of desired native or nonnative biochemical(s) in (micro)organisms can be achieved through metabolic engineering. Appropriate rewiring of cell metabolism is performed by making rational changes such as insertion, up-/down-regulation and knockout of genes and consequently metabolic reactions. Finding appropriate targets (including proper sets of reactions to be knocked out) for metabolic engineering to design optimal production strains has been the goal of a number of computational algorithms. We developed FastKnock, an efficient next-generation algorithm for identifying all possible knockout strategies (with a predefined maximum number of reaction deletions) for the growth-coupled overproduction of biochemical(s) of interest. We achieve this by developing a special depth-first traversal algorithm that allows us to prune the search space significantly. This leads to a drastic reduction in execution time. We evaluate the performance of the FastKnock algorithm using various Escherichia coli genome-scale metabolic models in different conditions (minimal and rich mediums) for the overproduction of a number of desired metabolites. FastKnock efficiently prunes the search space to less than 0.2% for quadruple- and 0.02% for quintuple-reaction knockouts. Compared to the classic approaches such as OptKnock and the state-of-the-art techniques such as MCSEnumerator methods, FastKnock found many more beneficial and important practical solutions. The availability of all the solutions provides the opportunity to further characterize, rank and select the most appropriate intervention strategy based on any desired evaluation index. Our implementation of the FastKnock method in Python is publicly available at https://github.com/leilahsn/FastKnock .

A Fermentation State Marker Rule Design Task in Metabolic Engineering

There are several ways in which mathematical modeling is used in fermentation control, but mechanistic mathematical genome-scale models of metabolism within the cell have not been applied or implemented so far. As part of the metabolic engineering task setting, we propose that metabolite fluxes and/or biomass growth rate be used to search for a fermentation steady state marker rule. During fermentation, the bioreactor control system can automatically detect the desired steady state using a logical marker rule. The marker rule identification can be also integrated with the production growth coupling approach, as presented in this study. A design of strain with marker rule is demonstrated on genome scale metabolic model iML1515 of Escherichia coli MG1655 proposing two gene deletions enabling a measurable marker rule for succinate production using glucose as a substrate. The marker rule example at glucose consumption 10.0 is: IF (specific growth rate μ is above 0.060 h-1, AND CO2 production under 1.0, AND ethanol production above 5.5), THEN succinate production is within the range 8.2-10, where all metabolic fluxes units are mmol ∗ gDW-1 ∗ h-1. An objective function for application in metabolic engineering, including productivity features and rule detecting sensor set characterizing parameters, is proposed. Two-phase approach to implementing marker rules in the cultivation control system is presented to avoid the need for a modeler during production.

OptEnvelope: A target point guided method for growth-coupled production using knockouts

Finding the best knockout strategy for coupling biomass growth and production of a target metabolite using a mathematic model of metabolism is a challenge in biotechnology. In this research, a three-step method named OptEnvelope is presented based on finding minimal set of active reactions for a target point in the feasible solution space (envelope) using a mixed-integer linear programming formula. The method initially finds the reduced desirable solution space envelope in the product versus biomass plot by removing all inactive reactions. Then, with reinsertion of the deleted reactions, OptEnvelope attempts to reduce the number of knockouts so that the desirable production envelope is preserved. Additionally, OptEnvelope searches for envelopes with higher minimum production rates or fewer knockouts by evaluating different target points within the desired solution space. It is possible to limit the maximal number of knockouts. The method was implemented on metabolic models of E. coli and S. cerevisiae to test the method benchmarking the capability of these industrial microbes for overproduction of acetate and glycerol under aerobic conditions and succinate and ethanol under anaerobic conditions. The results illustrate that OptEnvelope is capable to find multiple strong coupled envelopes located in the desired solution space because of its novel target point oriented strategy of envelope search. The results indicate that E. coli is more appropriate to produce acetate and succinate while S. cerevisiae is a better host for glycerol production. Gene deletions for some of the proposed reaction knockouts have been previously reported to increase the production of these metabolites in experiments. Both organisms are suitable for ethanol production, however, more knockouts for the adaptation of E. coli are required. OptEnvelope is available at https://github.com/lv-csbg/optEnvelope.

FastKnock: An efficient next-generation approach to identify all knockout strategies for strain optimization

Overproduction of desired native or nonnative biochemical(s) in (micro)organisms can be achieved through metabolic engineering. Appropriate rewiring of cell metabolism is performed making rational changes such as insertion, up-/down-regulation and knockout of genes and consequently metabolic reactions. Finding appropriate targets (including proper sets of reactions to be knocked out) for metabolic engineering to design optimal production strains has been the goal of a number of computational algorithms. We developed FastKnock, an efficient next-generation algorithm for identifying all possible knockout strategies for the growth-coupled overproduction of biochemical(s) of interest. We achieve this by developing a special depth-first traversal algorithm that allows us to prune the search space significantly. This leads to a drastic reduction in execution time. We evaluate the performance of the FastKnock algorithm using three Escherichia coli genome-scale metabolic models in different conditions (minimal and rich mediums) for the overproduction of a number of desired metabolites. FastKnock efficiently prunes the search space to less than 0.2% for quadruple and 0.02% for quintuple-reaction knockouts. Compared to the classic approaches such as OptKnock and the state-of-the-art techniques such as MCSEnumerator methods, FastKnock found many more useful and important practical solutions. The availability of all the solutions provides the opportunity to further characterize and select the most appropriate intervention strategy based on any desired evaluation index. Our implementation of the FastKnock method in Python is publicly available at https://github.com/leilahsn/FastKnock.

Assessing and reducing phenotypic instability in cyanobacteria

Cyanobacteria have promising potential as sustainable cell factories. However, one challenge that is still largely unreported in scaling-up cyanobacteria bioproduction is phenotypic instability, where the emergence and selection of nonproducing cells leading to loss in production has longer evolutionary timescales to take place in industrial-scale bioreactors. Quantifying phenotypic instability early on in strain development allows researchers to make informed decisions on whether to proceed with scalable designs, or if present, devise countermeasures to reduce instability. One particularly effective strategy to mitigate instability is the use of genome-scale metabolic models to design growth-coupled production strains. In silico studies have predicted that creating certain cofactor imbalances or removing recycling reactions in cyanobacteria can be exploited to stably produce a wide variety of metabolites.

StrainDesign: a comprehensive Python package for computational design of metabolic networks

SUMMARY

Various constraint-based optimization approaches have been developed for the computational analysis and design of metabolic networks. Herein, we present StrainDesign, a comprehensive Python package that builds upon the COBRApy toolbox and integrates the most popular metabolic design algorithms, including nested strain optimization methods such as OptKnock, RobustKnock and OptCouple as well as the more general minimal cut sets approach. The optimization approaches are embedded in individual modules, which can also be combined for setting up more elaborate strain design problems. Advanced features, such as the efficient integration of GPR rules and the possibility to consider gene and reaction additions or regulatory interventions, have been generalized and are available for all modules. The package uses state-of-the-art preprocessing methods, supports multiple solvers and provides a number of enhanced tools for analyzing computed intervention strategies including 2D and 3D plots of user-selected metabolic fluxes or yields. Furthermore, a user-friendly interface for the StrainDesign package has been implemented in the GUI-based metabolic modeling software CNApy. StrainDesign provides thus a unique and rich framework for computational strain design in Python, uniting many algorithmic developments in the field and allowing modular extension in the future.

AVAILABILITY AND IMPLEMENTATION

The StrainDesign package can be retrieved from PyPi, Anaconda and GitHub (https://github.com/klamt-lab/straindesign) and is also part of the latest CNApy package.

Examining organic acid production potential and growth-coupled strategies in Issatchenkia orientalis using constraint-based modeling

Growth-coupling product formation can facilitate strain stability by aligning industrial objectives with biological fitness. Organic acids make up many building block chemicals that can be produced from sugars obtainable from renewable biomass. Issatchenkia orientalis is a yeast strain tolerant to acidic conditions and is thus a promising host for industrial production of organic acids. Here, we use constraint-based methods to assess the potential of computationally designing growth-coupled production strains for I. orientalis that produce 22 different organic acids under aerobic or microaerobic conditions. We explore native and engineered pathways using glucose or xylose as the carbon substrates as proxy constituents of hydrolyzed biomass. We identified growth-coupled production strategies for 37 of the substrate-product pairs, with 15 pairs achieving production for any growth rate. We systematically assess the strain design solutions and categorize the underlying principles involved.

Synthetic metabolism for biohalogenation

The pressing need for novel bioproduction approaches faces a limitation in the number and type of molecules accessed through synthetic biology. Halogenation is widely used for tuning physicochemical properties of molecules and polymers, but traditional halogenation chemistry often lacks specificity and generates harmful by-products. Here, we pose that deploying synthetic metabolism tailored for biohalogenation represents an unique opportunity towards economically attractive and environmentally friendly organohalide production. On this background, we discuss growth-coupled selection of functional metabolic modules that harness the rich repertoire of biosynthetic and biodegradation capabilities of environmental bacteria for in vivo biohalogenation. By rationally combining these approaches, the chemical landscape of living cells can accommodate bioproduction of added-value organohalides which, as of today, are obtained by traditional chemistry.