On the feasibility of growth-coupled product synthesis in microbial strains

A hierarchy of metabolite exchanges in metabolic models of microbial species and communities

The metabolic network of an organism can be analyzed as a constraint-based model. This analysis can be biased, optimizing an objective such as growth rate, or unbiased, aiming to describe the full feasible space of metabolic fluxes through pathway analysis or random flux sampling. In particular, pathway analysis can decompose the flux space into fundamental and formally defined metabolic pathways. Unbiased methods scale poorly with network size due to combinatorial explosion, but a promising approach to improve scalability is to focus on metabolic subnetworks, e.g., cells' metabolite exchanges with each other and the environment, rather than the full metabolic networks. Here, we applied pathway enumeration and flux sampling to metabolite exchanges in microbial species and a microbial community, using models ranging from central carbon metabolism to genome-scale and focusing on pathway definitions that allow direct targeting of subnetworks such as metabolite exchanges (elementary conversion modes, elementary flux patterns, and minimal pathways). Enumerating growth-supporting metabolite exchanges, we found that metabolite exchanges from different pathway definitions were related through a hierarchy, and we show that this hierarchical relationship between pathways holds for metabolic networks and subnetworks more generally. Metabolite exchange frequencies, defined as the fraction of pathways in which each metabolite was exchanged, were similar across pathway definitions, with a few specific exchanges explaining large differences in pathway counts. This indicates that biological interpretation of predicted metabolite exchanges is robust to the choice of pathway definition, and it suggests strategies for more scalable pathway analysis. Our results also signal wider biological implications, facilitating detailed and interpretable analysis of metabolite exchanges and other subnetworks in fields such as metabolic engineering and synthetic biology.

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 .

Eliminating genes for a two-component system increases PHB productivity in Cupriavidus basilensis 4G11 under PHB suppressing, nonstress conditions

Species of bacteria from the genus Cupriavidus are known, in part, for their ability to produce high amounts of poly-hydroxybutyrate (PHB) making them attractive candidates for bioplastic production. The native synthesis of PHB occurs during periods of metabolic stress, and the process regulating the initiation of PHB accumulation in these organisms is not fully understood. Screening an RB-TnSeq transposon library of Cupriavidus basilensis 4G11 allowed us to identify two genes of an apparent, uncharacterized two-component system, which when omitted from the genome enable increased PHB productivity in balanced, nonstress growth conditions. We observe average increases in PHB productivity of 56% and 41% relative to the wildtype parent strain upon deleting each gene individually from the genome. The increased PHB phenotype disappears, however, in nitrogen-free unbalanced growth conditions suggesting the phenotype is specific to fast-growing, replete, nonstress growth. Bioproduction modeling suggests this phenotype could be due to a decreased reliance on metabolic stress induced by nitrogen limitation to initiate PHB production in the mutant strains. Due to uncertainty in the two-component system's input signal and regulon, the mechanism by which these genes impart this phenotype remains unclear. Such strains may allow for the use of single-stage, continuous bioreactor systems, which are far simpler than many PHB bioproduction schemes used previously, given a similar product yield to batch systems in such a configuration. Bioproductivity modeling suggests that omitting this regulation in the cells may increase PHB productivity up to 24% relative to the wildtype organism when using single-stage continuous systems. This work expands our understanding of the regulation of PHB accumulation in Cupriavidus, in particular the initiation of this process upon transition into unbalanced growth regimes.

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.

Tuning a high performing multiplexed-CRISPRi Pseudomonas putida strain to further enhance indigoidine production

In this study, a 14-gene edited Pseudomonas putida KT2440 strain for heterologous indigoidine production was examined using three distinct omic datasets. Transcriptomic data indicated that CRISPR/dCpf1-interference (CRISPRi) mediated multiplex repression caused global gene expression changes, implying potential undesirable changes in metabolic flux. ¹³C-metabolic flux analysis (¹³C-MFA) revealed that the core P. putida flux network after CRISPRi repression was conserved, with moderate reduction of TCA cycle and pyruvate shunt activity along with glyoxylate shunt activation during glucose catabolism. Metabolomic results identified a change in intracellular TCA metabolites and extracellular metabolite secretion profiles (sugars and succinate overflow) in the engineered strains. These omic analyses guided further strain engineering, with a random mutagenesis screen first identifying an optimal ribosome binding site (RBS) for Cpf1 that enabled stronger product-substrate pairing (1.6-fold increase). Then, deletion strains were constructed with excision of the PHA operon (ΔphaAZC-IID) resulting in a 2.2-fold increase in indigoidine titer over the optimized Cpf1-RBS construct at the end of the growth phase (∼6 h). The maximum indigoidine titer (at 72 h) in the ΔphaAZC-IID strain had a 1.5-fold and 1.8-fold increase compared to the optimized Cpf1-RBS construct and the original strain, respectively. Overall, this study demonstrated that integration of omic data types is essential for understanding responses to complex metabolic engineering designs and directly quantified the effect of such modifications on central metabolism.

Metabolic pathway design for growth-associated phenylalanine production using synthetically designed mutualism

Combination of growth-associated pathway engineering based on flux balance analysis (FBA) and adaptive laboratory evolution (ALE) is a powerful approach to enhance the production of useful compounds. However, the feasibility of such growth-associated pathway designs depends on the type of target compound. In the present study, FBA predicted a set of gene deletions (pykA, pykF, ppc, zwf, and adhE) that leads to growth-associated phenylalanine production in Escherichia coli. The knockout strain is theoretically enforced to produce phenylalanine only at high growth yields, and could not be applied to the ALE experiment because of a severe growth defect. To overcome this challenge, we propose a novel approach for ALE based on mutualistic co-culture for coupling growth and production, regardless of the growth rate. We designed a synthetic mutualism of a phenylalanine-producing leucine-auxotrophic strain (KF strain) and a leucine-producing phenylalanine-auxotrophic strain (KL strain) and performed an ALE experiment for approximately 160 generations. The evolved KF strain (KF-E strain) grew in a synthetic medium (with glucose as the main carbon source) supplemented with leucine, while severe growth defects were observed in the parental KF strain. The phenylalanine yield of the KF-E strain was 2.3 times higher than that of the KF strain.

EFMlrs: a Python package for elementary flux mode enumeration via lexicographic reverse search

Background

Elementary flux mode (EFM) analysis is a well-established, yet computationally challenging approach to characterize metabolic networks. Standard algorithms require huge amounts of memory and lack scalability which limits their application to single servers and consequently limits a comprehensive analysis to medium-scale networks. Recently, Avis et al. developed mplrs—a parallel version of the lexicographic reverse search (lrs) algorithm, which, in principle, enables an EFM analysis on high-performance computing environments (Avis and Jordan. mplrs: a scalable parallel vertex/facet enumeration code. arXiv:1511.06487, 2017). Here we test its applicability for EFM enumeration.

Results

We developed EFMlrs, a Python package that gives users access to the enumeration capabilities of mplrs. EFMlrs uses COBRApy to process metabolic models from sbml files, performs loss-free compressions of the stoichiometric matrix, and generates suitable inputs for mplrs as well as efmtool, providing support not only for our proposed new method for EFM enumeration but also for already established tools. By leveraging COBRApy, EFMlrs also allows the application of additional reaction boundaries and seamlessly integrates into existing workflows.

Conclusion

We show that due to mplrs’s properties, the algorithm is perfectly suited for high-performance computing (HPC) and thus offers new possibilities for the unbiased analysis of substantially larger metabolic models via EFM analyses. EFMlrs is an open-source program that comes together with a designated workflow and can be easily installed via pip.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04417-9.

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

A widely used design principle for metabolic engineering of microorganisms aims to introduce interventions that enforce growth-coupled product synthesis such that the product of interest becomes a (mandatory) by-product of growth. However, different variants and partially contradicting notions of growth-coupled production (GCP) exist. Herein, we propose an ontology for the different degrees of GCP and clarify their relationships. Ordered by coupling degree, we distinguish four major classes: potentially, weakly, and directionally growth-coupled production (pGCP, wGCP, dGCP) as well as substrate-uptake coupled production (SUCP). We then extend the framework of Minimal Cut Sets (MCS), previously used to compute dGCP and SUCP strain designs, to allow inclusion of implicit optimality constraints, a feature required to compute pGCP and wGCP designs. This extension closes the gap between MCS-based and bilevel-based strain design approaches and enables computation (and comparison) of designs for all GCP classes within a single framework. By computing GCP strain designs for a range of products, we illustrate the hierarchical relationships between the different coupling degrees. We find that feasibility of coupling is not affected by the chosen GCP degree and that strongest coupling (SUCP) requires often only one or two more interventions than wGCP and dGCP. Finally, we show that the principle of coupling can be generalized to couple product synthesis with other cellular functions than growth, for example, with net ATP formation. This work provides important theoretical results and algorithmic developments and a unified terminology for computational strain design based on GCP.

Overproduction of D-pantothenic acid via fermentation conditions optimization and isoleucine feeding from recombinant Escherichia coli W3110

D-pantothenic acid (D-PA), as a crucial vitamin, is widely used in food, animal feed, cosmetics, and pharmaceutical industries. In our previous work, recombinant Escherichia coli W3110 for production of D-PA was constructed through metabolic pathway modification. In this study, to enhance D-PA production, statistical optimization techniques including Plackett-Burman (PB) design and Box-Behnken design (BBD) first were adopted to optimize the culture condition. The results showed that the glucose, β-alanine and (NH₄)₂SO₄ have the most significant effects on D-PA biosynthesis. The response surface model based on BBD predicted that the optimal concentration is glucose 56.0 g/L, β-alanine 2.25 g/L and (NH4)₂SO₄ 11.8 g/L, the D-PA titer increases from 3.2 g/L to 6.73 g/L shake flask fermentation. For the fed-batch fermentation in 5 L fermenter, the isoleucine feeding strategy greatly increased the titer and productivity of D-PA. As a result, titer (31.6 g/L) and productivity (13.2 g/L·d) of D-PA were achieved, they increased by 4.66 times and 2.65 times, respectively, compared with batch culture. At the same time, the accumulation of acetate reduced from 29.79 g/L to 8.55 g/L in the fed-batch fermentation. These results demonstrated that the optimization of medium composition and the cell growth rate are important to increase the concentration of D-PA for microbial fermentation. This work laid the foundation for further research on the application of D-PA microbial synthesis.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02773-0.

Advanced strategies and tools to facilitate and streamline microbial adaptive laboratory evolution

Adaptive laboratory evolution (ALE) has served as a historic microbial engineering method that mimics natural selection to obtain desired microbes. The past decade has witnessed improvements in all aspects of ALE workflow, in terms of growth coupling, genotypic diversification, phenotypic selection, and genotype-phenotype mapping. The developing growth-coupling strategies facilitate ALE to a wider range of appealing traits. In vivo mutagenesis methods and multiplexed automated culture platforms open new gates to streamline its execution. Meanwhile, the application of multi-omics analyses and multiplexed genetic engineering promote efficient knowledge mining. This article provides a comprehensive and updated review of these advances, highlights newest significant applications, and discusses future improvements, aiming to provide a practical guide for implementation of novel, effective, and efficient ALE experiments.

Understanding FBA Solutions under Multiple Nutrient Limitations

Genome-scale stoichiometric modeling methods, in particular Flux Balance Analysis (FBA) and variations thereof, are widely used to investigate cell metabolism and to optimize biotechnological processes. Given (1) a metabolic network, which can be reconstructed from an organism’s genome sequence, and (2) constraints on reaction rates, which may be based on measured nutrient uptake rates, FBA predicts which reactions maximize an objective flux, usually the production of cell components. Although FBA solutions may accurately predict the metabolic behavior of a cell, the actual flux predictions are often hard to interpret. This is especially the case for conditions with many constraints, such as for organisms growing in rich nutrient environments: it remains unclear why a certain solution was optimal. Here, we rationalize FBA solutions by explaining for which properties the optimal combination of metabolic strategies is selected. We provide a graphical formalism in which the selection of solutions can be visualized; we illustrate how this perspective provides a glimpse of the logic that underlies genome-scale modeling by applying our formalism to models of various sizes.

Speeding up the core algorithm for the dual calculation of minimal cut sets in large metabolic networks

Background

The concept of minimal cut sets (MCS) has become an important mathematical framework for analyzing and (re)designing metabolic networks. However, the calculation of MCS in genome-scale metabolic models is a complex computational problem. The development of duality-based algorithms in the last years allowed the enumeration of thousands of MCS in genome-scale networks by solving mixed-integer linear problems (MILP). A recent advancement in this field was the introduction of the MCS² approach. In contrast to the Farkas-lemma-based dual system used in earlier studies, the MCS² approach employs a more condensed representation of the dual system based on the nullspace of the stoichiometric matrix, which, due to its reduced dimension, holds promise to further enhance MCS computations.

Results

In this work, we introduce several new variants and modifications of duality-based MCS algorithms and benchmark their effects on the overall performance. As one major result, we generalize the original MCS² approach (which was limited to blocking the operation of certain target reactions) to the most general case of MCS computations with arbitrary target and desired regions. Building upon these developments, we introduce a new MILP variant which allows maximal flexibility in the formulation of MCS problems and fully leverages the reduced size of the nullspace-based dual system. With a comprehensive set of benchmarks, we show that the MILP with the nullspace-based dual system outperforms the MILP with the Farkas-lemma-based dual system speeding up MCS computation with an averaged factor of approximately 2.5. We furthermore present several simplifications in the formulation of constraints, mainly related to binary variables, which further enhance the performance of MCS-related MILP. However, the benchmarks also reveal that some highly condensed formulations of constraints, especially on reversible reactions, may lead to worse behavior when compared to variants with a larger number of (more explicit) constraints and involved variables.

Conclusions

Our results further enhance the algorithmic toolbox for MCS calculations and are of general importance for theoretical developments as well as for practical applications of the MCS framework.

Genome-scale metabolic rewiring improves titers rates and yields of the non-native product indigoidine at scale

High titer, rate, yield (TRY), and scalability are challenging metrics to achieve due to trade-offs between carbon use for growth and production. To achieve these metrics, we take the minimal cut set (MCS) approach that predicts metabolic reactions for elimination to couple metabolite production strongly with growth. We compute MCS solution-sets for a non-native product indigoidine, a sustainable pigment, in Pseudomonas putida KT2440, an emerging industrial microbe. From the 63 solution-sets, our omics guided process identifies one experimentally feasible solution requiring 14 simultaneous reaction interventions. We implement a total of 14 genes knockdowns using multiplex-CRISPRi. MCS-based solution shifts production from stationary to exponential phase. We achieve 25.6 g/L, 0.22 g/l/h, and ~50% maximum theoretical yield (0.33 g indigoidine/g glucose). These phenotypes are maintained from batch to fed-batch mode, and across scales (100-ml shake flasks, 250-ml ambr®, and 2-L bioreactors).

Control strategies to manage trade-offs during microbial production

When engineering microbes to overproduce a target molecule, engineers face multiple layers of trade-offs to allocate limited cellular resources between the target pathway and native cellular systems. These trade-offs arise from limited free ribosomes during translation, competition for metabolic precursors, as well as the negative relationship between production and growth rate. To achieve high production performance, microbes need to spontaneously make decisions in the dynamic and heterogeneous fermentation environment. In this review, we discuss recent advances in microbial control strategies that are used to manage these trade-offs and to improve microbial production. This review focuses on design principles and compares different implementations, with the hope to provide guidelines to future microbial engineering.

Harnessing Natural Modularity of Metabolism with Goal Attainment Optimization to Design a Modular Chassis Cell for Production of Diverse Chemicals

Modular design is key to achieve efficient and robust systems across engineering disciplines. Modular design potentially offers advantages to engineer microbial systems for biocatalysis, bioremediation, and biosensing, overcoming the slow and costly design-build-test-learn cycles in the conventional cell engineering approach. These systems consist of a modular (chassis) cell compatible with exchangeable modules that enable programmed functions such as overproduction of a desirable chemical. We previously proposed a multiobjective optimization framework coupled with metabolic flux models to design modular cells and solved it using multiobjective evolutionary algorithms. However, such approach might not achieve solution optimality and hence limits design options for experimental implementation. In this study, we developed the goal attainment formulation compatible with optimization algorithms that guarantee solution optimality. We applied goal attainment to design an Escherichia coli modular cell capable of synthesizing all molecules in a biochemically diverse library at high yields and rates with only a few genetic manipulations. To elucidate modular organization of the designed cells, we developed a flux variance clustering (FVC) method by identifying reactions with high flux variance and clustering them to identify metabolic modules. Using FVC, we identified reaction usage patterns for different modules in the modular cell, revealing that its broad pathway compatibility is enabled by the natural modularity and flexible flux capacity of endogenous core metabolism. Overall, this study not only sheds light on modularity in metabolic networks from their topology and metabolic functions but also presents a useful synthetic biology toolbox to design modular cells with broad applications.

An extended and generalized framework for the calculation of metabolic intervention strategies based on minimal cut sets

The concept of minimal cut sets (MCS) provides a flexible framework for analyzing properties of metabolic networks and for computing metabolic intervention strategies. In particular, it has been used to support the targeted design of microbial strains for bio-based production processes. Herein we present a number of major extensions that generalize the existing MCS approach and broaden its scope for applications in metabolic engineering. We first introduce a modified approach to integrate gene-protein-reaction associations (GPR) in the metabolic network structure for the computation of gene-based intervention strategies. In particular, we present a set of novel compression rules for GPR associations, which effectively speedup the computation of gene-based MCS by a factor of up to one order of magnitude. These rules are not specific for MCS and as well applicable to other computational strain design methods. Second, we enhance the MCS framework by allowing the definition of multiple target (undesired) and multiple protected (desired) regions. This enables precise tailoring of the metabolic solution space of the designed strain with unlimited flexibility. Together with further generalizations such as individual cost factors for each intervention, direct combinations of reaction/gene deletions and additions as well as the possibility to search for substrate co-feeding strategies, the scope of the MCS framework could be broadly extended. We demonstrate the applicability and performance benefits of the described developments by computing (gene-based) Escherichia coli strain designs for the bio-based production of 2,3-butanediol, a chemical, that has recently received much attention in the field of metabolic engineering. With our extended framework, we could identify promising strain designs that were formerly unpredictable, including those based on substrate co-feeding.

Characterizing and ranking computed metabolic engineering strategies

MOTIVATION

The computer-aided design of metabolic intervention strategies has become a key component of an integrated metabolic engineering approach and a broad range of methods and algorithms has been developed for this task. Many of these algorithms enforce coupling of growth with product synthesis and may return thousands of possible intervention strategies from which the most suitable strategy must then be selected.

RESULTS

This work focuses on how to evaluate and rank, in a meaningful way, a given pool of computed metabolic engineering strategies for growth-coupled product synthesis. Apart from straightforward criteria, such as a preferably small number of necessary interventions, a reasonable growth rate and a high product yield, we present several new criteria useful to pick the most suitable intervention strategy. Among others, we investigate the robustness of the intervention strategies by searching for metabolites that may disrupt growth coupling when accumulated or secreted and by checking whether the interventions interrupt pathways at their origin (preferable) or at downstream steps. We also assess thermodynamic properties of the pathway(s) favored by the intervention strategy. Furthermore, strategies that have a significant overlap with alternative solutions are ranked higher because they provide flexibility in implementation. We also introduce the notion of equivalence classes for grouping intervention strategies with identical solution spaces. Our ranking procedure involves in total ten criteria and we demonstrate its applicability by assessing knockout-based intervention strategies computed in a genome-scale model of E.coli for the growth-coupled synthesis of l-methionine and of the heterologous product 1,4-butanediol.

AVAILABILITY AND IMPLEMENTATION

The MATLAB scripts that were used to characterize and rank the example intervention strategies are available at http://www2.mpi-magdeburg.mpg.de/projects/cna/etcdownloads.html.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A Systems-Based Approach for Cyanide Overproduction by Bacillus megaterium for Gold Bioleaching Enhancement

With the constant accumulation of electronic waste, extracting precious metals contained therein is becoming a major challenge for sustainable development. Bacillus megaterium is currently one of the microbes used for the production of cyanide, which is the main leaching agent for gold recovery. The present study aimed to propose a strategy for metabolic engineering of B. megaterium to overproduce cyanide, and thus ameliorate the bioleaching process. For this, we employed constraint-based modeling, running in silico simulations on iJA1121, the genome-scale metabolic model of B. megaterium DSM319. Flux balance analysis (FBA) was initially used to identify amino acids to be added to the culture medium. Considering cyanide as the desired product, we used growth-coupled methods, constrained minimal cut sets (cMCSs) and OptKnock to identify gene inactivation targets. To identify gene overexpression targets, flux scanning based on enforced objective flux (FSEOF) was performed. Further analysis was carried out on the identified targets to determine compounds with beneficial regulatory effects. We have proposed a chemical-defined medium for accelerating cyanide production on the basis of microplate assays to evaluate the components with the greatest improving effects. Accordingly, the cultivation of B. megaterium DSM319 in a chemically-defined medium with 5.56 mM glucose as the carbon source, and supplemented with 413 μM cysteine, led to the production of considerably increased amounts of cyanide. Bioleaching experiments were successfully performed in this medium to recover gold and copper from telecommunication printed circuit boards. The results of inductively coupled plasma (ICP) analysis confirmed that gold recovery peaked out at around 55% after 4 days, whereas copper recovery continued to increase for several more days, peaking out at around 85%. To further validate the bioleaching results, FESEM, XRD, FTIR, and EDAX mapping analyses were performed. We concluded that the proposed strategy represents a viable route for improving the performance of the bioleaching processes.

Key Challenges in Designing CHO Chassis Platforms

Following the success of and the high demand for recombinant protein-based therapeutics during the last 25 years, the pharmaceutical industry has invested significantly in the development of novel treatments based on biologics. Mammalian cells are the major production systems for these complex biopharmaceuticals, with Chinese hamster ovary (CHO) cell lines as the most important players. Over the years, various engineering strategies and modeling approaches have been used to improve microbial production platforms, such as bacteria and yeasts, as well as to create pre-optimized chassis host strains. However, the complexity of mammalian cells curtailed the optimization of these host cells by metabolic engineering. Most of the improvements of titer and productivity were achieved by media optimization and large-scale screening of producer clones. The advances made in recent years now open the door to again consider the potential application of systems biology approaches and metabolic engineering also to CHO. The availability of a reference genome sequence, genome-scale metabolic models and the growing number of various “omics” datasets can help overcome the complexity of CHO cells and support design strategies to boost their production performance. Modular design approaches applied to engineer industrially relevant cell lines have evolved to reduce the time and effort needed for the generation of new producer cells and to allow the achievement of desired product titers and quality. Nevertheless, important steps to enable the design of a chassis platform similar to those in use in the microbial world are still missing. In this review, we highlight the importance of mammalian cellular platforms for the production of biopharmaceuticals and compare them to microbial platforms, with an emphasis on describing novel approaches and discussing still open questions that need to be resolved to reach the objective of designing enhanced modular chassis CHO cell lines.

Automatic construction of metabolic models with enzyme constraints

Background

In order to improve the accuracy of constraint-based metabolic models, several approaches have been developed which intend to integrate additional biological information. Two of these methods, MOMENT and GECKO, incorporate enzymatic (kcat) parameters and enzyme mass constraints to further constrain the space of feasible metabolic flux distributions. While both methods have been proven to deliver useful extensions of metabolic models, they may considerably increase size and complexity of the models and there is currently no tool available to fully automate generation and calibration of such enzyme-constrained models from given stoichiometric models.

Results

In this work we present three major developments. We first conceived short MOMENT (sMOMENT), a simplified version of the MOMENT approach, which yields the same predictions as MOMENT but requires significantly fewer variables and enables direct inclusion of the relevant enzyme constraints in the standard representation of a constraint-based model. When measurements of enzyme concentrations are available, these can be included as well leading in the extreme case, where all enzyme concentrations are known, to a model representation that is analogous to the GECKO approach. Second, we developed the AutoPACMEN toolbox which allows an almost fully automated creation of sMOMENT-enhanced stoichiometric metabolic models. In particular, this includes the automatic read-out and processing of relevant enzymatic data from different databases and the reconfiguration of the stoichiometric model with embedded enzymatic constraints. Additionally, tools have been developed to adjust (kcat and enzyme pool) parameters of sMOMENT models based on given flux data. We finally applied the new sMOMENT approach and the AutoPACMEN toolbox to generate an enzyme-constrained version of the E. coli genome-scale model iJO1366 and analyze its key properties and differences with the standard model. In particular, we show that the enzyme constraints improve flux predictions (e.g., explaining overflow metabolism and other metabolic switches) and demonstrate, for the first time, that these constraints can markedly change the spectrum of metabolic engineering strategies for different target products.

Conclusions

The methodological and tool developments presented herein pave the way for a simplified and routine construction and analysis of enzyme-constrained metabolic models.

Exploiting Day- and Night-Time Metabolism of Synechocystis sp. PCC 6803 for Fitness-Coupled Fumarate Production around the Clock

Cyanobacterial cell factories are widely researched for the sustainable production of compounds directly from CO2. Their application, however, has been limited for two reasons. First, traditional approaches have been shown to lead to unstable cell factories that lose their production capability when scaled to industrial levels. Second, the alternative approaches developed so far are mostly limited to growing conditions, which are not always the case in industry, where nongrowth periods tend to occur (e.g., darkness). We tackled both by generalizing the concept of growth-coupled production to fitness coupling. The feasibility of this new approach is demonstrated for the production of fumarate by constructing the first stable dual-strategy cell factory. We exploited circadian metabolism using both systems and synthetic biology tools, resulting in the obligatorily coupling of fumarate to either biomass or energy production. Resorting to laboratory evolution experiments, we show that this engineering approach is more stable than conventional methods.

Determination of growth-coupling strategies and their underlying principles

BACKGROUND

Metabolic coupling of product synthesis and microbial growth is a prominent approach for maximizing production performance. Growth-coupling (GC) also helps stabilizing target production and allows the selection of superior production strains by adaptive laboratory evolution. To support the implementation of growth-coupling strain designs, we seek to identify biologically relevant, metabolic principles that enforce strong growth-coupling on the basis of reaction knockouts.

RESULTS

We adapted an established bilevel programming framework to maximize the minimally guaranteed production rate at a fixed, medium growth rate. Using this revised formulation, we identified various GC intervention strategies for metabolites of the central carbon metabolism, which were examined for GC generating principles under diverse conditions. Curtailing the metabolism to render product formation an essential carbon drain was identified as one major strategy generating strong coupling of metabolic activity and target synthesis. Impeding the balancing of cofactors and protons in the absence of target production was the underlying principle of all other strategies and further increased the GC strength of the aforementioned strategies.

CONCLUSION

Maximizing the minimally guaranteed production rate at a medium growth rate is an attractive principle for the identification of strain designs that couple growth to target metabolite production. Moreover, it allows for controlling the inevitable compromise between growth coupling strength and the retaining of microbial viability. With regard to the corresponding metabolic principles, generating a dependency between the supply of global metabolic cofactors and product synthesis appears to be advantageous in enforcing strong GC for any metabolite. Deriving such strategies manually, is a hard task, due to which we suggest incorporating computational metabolic network analyses in metabolic engineering projects seeking to determine GC strain designs.

The emergence of adaptive laboratory evolution as an efficient tool for biological discovery and industrial biotechnology

Harnessing the process of natural selection to obtain and understand new microbial phenotypes has become increasingly possible due to advances in culturing techniques, DNA sequencing, bioinformatics, and genetic engineering. Accordingly, Adaptive Laboratory Evolution (ALE) experiments represent a powerful approach both to investigate the evolutionary forces influencing strain phenotypes, performance, and stability, and to acquire production strains that contain beneficial mutations. In this review, we summarize and categorize the applications of ALE to various aspects of microbial physiology pertinent to industrial bioproduction by collecting case studies that highlight the multitude of ways in which evolution can facilitate the strain construction process. Further, we discuss principles that inform experimental design, complementary approaches such as computational modeling that help maximize utility, and the future of ALE as an efficient strain design and build tool driven by growing adoption and improvements in automation.

Developing a pyruvate-driven metabolic scenario for growth-coupled microbial production

Microbial-based chemical synthesis serves as a promising approach for sustainable production of industrially important products. However, limited production performance caused by metabolic burden or genetic variations poses one of the major challenges in achieving an economically viable biomanufacturing process. To address this issue, one superior strategy is to couple the product synthesis with cellular growth, which renders production obligatory for cell survival. Here we create a pyruvate-driven metabolic scenario in engineered Escherichia coli for growth-coupled bioproduction, with which we demonstrate its application in boosting production of anthranilate and its derivatives. Deletion of a minimal set of endogenous pyruvate-releasing pathways engenders anthranilate synthesis as the salvage route for pyruvate generation to support cell growth, concomitant with simultaneous anthranilate production. Further introduction of native and non-native downstream pathways affords production enhancement of two anthranilate-derived high-value products including L-tryptophan and cis, cis-muconic acid from different carbon sources. The work reported here presents a new growth-coupled strategy with demonstrated feasibility for promoting microbial production.

Prediction of Rate-Limiting Reactions for Growth-Associated Production Using a Constraint-Based Approach

Identification of a rate-limiting step in pathways is a key challenge in metabolic engineering. Although the prediction of rate-limiting steps using a kinetic model is a powerful approach, there are several technical hurdles for developing a kinetic model. In this study, an in silico screening algorithm of key enzyme for metabolic engineering is developed to identify the possible rate-limiting reactions for the growth-coupled target production using a stoichiometric model without any experimental data and kinetic parameters. In this method, for each reaction, an upper-bound flux constraint is imposed and the target production is predicted by linear programming. When the constraint decreases the target production at the optimal growth state, the reaction is thought to be a possible rate-limiting step. For validation, this method is applied to the production of succinate or 1,4-butanediol (1,4-BDO) in Escherichia coli, in which the experimental engineering for eliminating rate-limiting steps has been previously reported. In succinate production from glycerol, nine reactions including phosphoenolpyruvate carboxylase are predicted as the rate-limiting steps. In 1,4-BDO production from glucose, eight reactions including pyruvate dehydrogenase are predicted as the rate-limiting steps. These predictions include experimentally identified rate-limiting steps, which would contribute to metabolic engineering as a practical tool for screening candidates of rate-limiting reactions.

Comparison of pathway analysis and constraint-based methods for cell factory design

BACKGROUND

Computational strain optimisation methods (CSOMs) have been successfully used to exploit genome-scale metabolic models, yielding strategies useful for allowing compound overproduction in metabolic cell factories. Minimal cut sets are particularly interesting since their definition allows searching for intervention strategies that impose strong growth-coupling phenotypes, and are not subject to optimality bias when compared with simulation-based CSOMs. However, since both types of methods have different underlying principles, they also imply different ways to formulate metabolic engineering problems, posing an obstacle when comparing their outputs.

RESULTS

In this work, we perform an in-depth analysis of potential strategies that can be obtained with both methods, providing a critical comparison of performance, robustness, predicted phenotypes as well as strategy structure and size. To this end, we devised a pipeline including enumeration of strategies from evolutionary algorithms (EA) and minimal cut sets (MCS), filtering and flux analysis of predicted mutants to optimize the production of succinic acid in Saccharomyces cerevisiae. We additionally attempt to generalize problem formulations for MCS enumeration within the context of growth-coupled product synthesis. Strategies from evolutionary algorithms show the best compromise between acceptable growth rates and compound overproduction. However, constrained MCSs lead to a larger variety of phenotypes with several degrees of growth-coupling with production flux. The latter have proven useful in revealing the importance, in silico, of the gamma-aminobutyric acid shunt and manipulation of cofactor pools in growth-coupled designs for succinate production, mechanisms which have also been touted as potentially useful for metabolic engineering.

CONCLUSIONS

The two main groups of CSOMs are valuable for finding growth-coupled mutants. Despite the limitations in maximum growth rates and large strategy sizes, MCSs help uncover novel mechanisms for compound overproduction and thus, analyzing outputs from both methods provides a richer overview on strategies that can be potentially carried over in vivo.

A critical comparison of cellular and cell-free bioproduction systems

Conversion of biological feedstocks into value-added chemicals is mostly performed via microbial fermentation. An emerging alternative approach is the use of cell-free systems, consisting of purified enzymes and cofactors. Unfortunately, the in vivo and in vitro research communities rarely interact, which leads to oversimplifications and exaggerations that do not permit fair comparison of the two strategies and impede synergistic interactions. Here, we provide a comprehensive account for the advantages and drawbacks associated with each strategy, and further discuss recent research efforts that aim to breach the limits of cellular and cell-free production. We also explore emerging hybrid solutions that integrate the benefits of both worlds and could expand the boundaries of biosynthesis.

Modular design: Implementing proven engineering principles in biotechnology

Modular design is at the foundation of contemporary engineering, enabling rapid, efficient, and reproducible construction and maintenance of complex systems across applications. Remarkably, modularity has recently been discovered as a governing principle in natural biological systems from genes to proteins to complex networks within a cell and organism communities. The convergent knowledge of natural and engineered modular systems provides a key to drive modern biotechnology to address emergent challenges associated with health, food, energy, and the environment. Here, we first present the theory and application of modular design in traditional engineering fields. We then discuss the significance and impact of modular architectures on systems biology and biotechnology. Next, we focus on the very recent theoretical and experimental advances in modular cell engineering that seeks to enable rapid and systematic development of microbial catalysts capable of efficiently synthesizing a large space of useful chemicals. We conclude with an outlook towards theoretical and practical opportunities for a more systematic and effective application of modular cell engineering in biotechnology.

OptCouple: Joint simulation of gene knockouts, insertions and medium modifications for prediction of growth-coupled strain designs

Biological production of chemicals is an attractive alternative to petrochemical-based production, due to advantages in environmental impact and the spectrum of feasible targets. However, engineering microbial strains to overproduce a compound of interest can be a long, costly and painstaking process. If production can be coupled to cell growth it is possible to use adaptive laboratory evolution to increase the production rate. Strategies for coupling production to growth, however, are often not trivial to find. Here we present OptCouple, a constraint-based modeling algorithm to simultaneously identify combinations of gene knockouts, insertions and medium supplements that lead to growth-coupled production of a target compound. We validated the algorithm by showing that it can find novel strategies that are growth-coupled in silico for a compound that has not been coupled to growth previously, as well as reproduce known growth-coupled strain designs for two different target compounds. Furthermore, we used OptCouple to construct an alternative design with potential for higher production. We provide an efficient and easy-to-use implementation of the OptCouple algorithm in the cameo Python package for computational strain design.

On the use of oxygenic photosynthesis for the sustainable production of commodity chemicals

A sustainable society will have to largely refrain from the use of fossil carbon deposits. In such a regime, renewable electricity can be harvested as a primary source of energy. However, as for the synthesis of carbon-based materials from bulk chemicals, an alternative is required. A sustainable approach towards this is the synthesis of commodity chemicals from CO₂ , water and sunlight. Multiple paths to achieve this have been designed and tested in the domains of chemistry and biology. In the latter, the use of both chemotrophic and phototrophic organisms has been advocated. 'Direct conversion' of CO₂ and H₂ O, catalyzed by an oxyphototroph, has excellent prospects to become the most economically competitive of these transformations, because of the relative ease of scale-up of this process. Significantly, for a wide range of energy and commodity products, a proof of principle via engineering of the corresponding production organism has been provided. In the optimization of a cyanobacterial production organism, a wide range of aspects has to be addressed. Of these, here we will put our focus on: (1) optimizing the (carbon) flux to the desired product; (2) increasing the genetic stability of the producing organism and (3) maximizing its energy conversion efficiency. Significant advances have been made on all these three aspects during the past 2 years and these will be discussed: (1) increasing the carbon partitioning to >50%; (2) aligning product formation with the growth of the cells and (3) expanding the photosynthetically active radiation region for oxygenic photosynthesis.

Yeast Genome-Scale Metabolic Models for Simulating Genotype-Phenotype Relations

Understanding genotype-phenotype dependency is a universal aim for all life sciences. While the complete genotype-phenotype relations remain challenging to resolve, metabolic phenotypes are moving within the reach through genome-scale metabolic model simulations. Genome-scale metabolic models are available for commonly investigated yeasts, such as model eukaryote and domesticated fermentation species Saccharomyces cerevisiae, and automatic reconstruction methods facilitate obtaining models for any sequenced species. The models allow for investigating genotype-phenotype relations through simulations simultaneously considering the effects of nutrient availability, and redox and energy homeostasis in cells. Genome-scale models also offer frameworks for omics data integration to help to uncover how the translation of genotypes to the apparent phenotypes is regulated at different levels. In this chapter, we provide an overview of the yeast genome-scale metabolic models and the simulation approaches for using these models to interrogate genotype-phenotype relations. We review the methodological approaches according to the underlying biological reasoning in order to inspire formulating novel questions and applications that the genome-scale metabolic models could contribute to. Finally, we discuss current challenges and opportunities in the genome-scale metabolic model simulations.

An Optimized Bistable Metabolic Switch To Decouple Phenotypic States during Anaerobic Fermentation

Metabolic engineers aim to genetically modify microorganisms to improve their ability to produce valuable compounds. Despite the prevalence of growth-coupled production processes, these strategies can significantly limit production rates. Instead, rates can be improved by decoupling and optimizing growth and production independently, and operating with a growth stage followed by a production stage. Here, we implement a bistable transcriptional controller to decouple and switch between these two states. We optimize the controller in anaerobic conditions, typical of industrial fermentations, to ensure stability and tight expression control, while improving switching dynamics. The stability of this controller can be maintained through a simulated seed train scale-up from 5 mL to 500 000 L, indicating industrial feasibility. Finally, we demonstrate a two-stage production process using our optimal construct to improve the instantaneous rate of lactate production by over 50%, motivating the use of these systems in broad metabolic engineering applications.

MoVE identifies metabolic valves to switch between phenotypic states

Metabolism is highly regulated, allowing for robust and complex behavior. This behavior can often be achieved by controlling a small number of important metabolic reactions, or metabolic valves. Here, we present a method to identify the location of such valves: the metabolic valve enumerator (MoVE). MoVE uses a metabolic model to identify genetic intervention strategies which decouple two desired phenotypes. We apply this method to identify valves which can decouple growth and production to systematically improve the rate and yield of biochemical production processes. We apply this algorithm to the production of diverse compounds and obtained solutions for over 70% of our targets, identifying a small number of highly represented valves to achieve near maximal growth and production. MoVE offers a systematic approach to identify metabolic valves using metabolic models, providing insight into the architecture of metabolic networks and accelerating the widespread implementation of dynamic flux redirection in diverse systems.

Integrated constraints based analysis of an engineered violacein pathway in Escherichia coli

Strategies towards optimal violacein biosynthesis, a potential drug molecule, need systems level coordination of enzymatic activities of individual genes in a multigene operon vioABCDE. Constraints-based flux balance analysis of an extended iAF1260 model (iAF1260vio) with a reconstructed violacein module predicted growth and violacein yields in Escherichia coli accurately. Shadow price (SP) analysis identified tryptophan metabolism and NADPH as limiting. Increased tryptophan levels in Δpgi & ΔpheA were validated using in silico gene deletion analysis. Phenotypic phase plane (PhPP) analysis highlighted sensitivity between tryptophan and NADPH for violacein synthesis at molar growth yields. A synthetic VioABCDE operon (SYNO) sequence was designed to maximize Codon Adaptive Index (CAI: 0.9) and tune translation initiation rates (TIR: 2-50 fold higher) in E. coli. All pSYN E. coli transformants produced higher violacein, with a maximum six-fold increase in yields. The rational design E. coli: ΔpheA SYN: gave the highest violacein titers (33.8 mg/l). Such integrated approaches targeting multiple molecular hierarchies in the cell can be extended further to increase violacein yields.

Selecting the Best: Evolutionary Engineering of Chemical Production in Microbes

Microbial cell factories have proven to be an economical means of production for many bulk, specialty, and fine chemical products. However, we still lack both a holistic understanding of organism physiology and the ability to predictively tune enzyme activities in vivo, thus slowing down rational engineering of industrially relevant strains. An alternative concept to rational engineering is to use evolution as the driving force to select for desired changes, an approach often described as evolutionary engineering. In evolutionary engineering, in vivo selections for a desired phenotype are combined with either generation of spontaneous mutations or some form of targeted or random mutagenesis. Evolutionary engineering has been used to successfully engineer easily selectable phenotypes, such as utilization of a suboptimal nutrient source or tolerance to inhibitory substrates or products. In this review, we focus primarily on a more challenging problem-the use of evolutionary engineering for improving the production of chemicals in microbes directly. We describe recent developments in evolutionary engineering strategies, in general, and discuss, in detail, case studies where production of a chemical has been successfully achieved through evolutionary engineering by coupling production to cellular growth.

A mathematical framework for yield (vs. rate) optimization in constraint-based modeling and applications in metabolic engineering

BACKGROUND

The optimization of metabolic rates (as linear objective functions) represents the methodical core of flux-balance analysis techniques which have become a standard tool for the study of genome-scale metabolic models. Besides (growth and synthesis) rates, metabolic yields are key parameters for the characterization of biochemical transformation processes, especially in the context of biotechnological applications. However, yields are ratios of rates, and hence the optimization of yields (as nonlinear objective functions) under arbitrary linear constraints is not possible with current flux-balance analysis techniques. Despite the fundamental importance of yields in constraint-based modeling, a comprehensive mathematical framework for yield optimization is still missing.

RESULTS

We present a mathematical theory that allows one to systematically compute and analyze yield-optimal solutions of metabolic models under arbitrary linear constraints. In particular, we formulate yield optimization as a linear-fractional program. For practical computations, we transform the linear-fractional yield optimization problem to a (higher-dimensional) linear problem. Its solutions determine the solutions of the original problem and can be used to predict yield-optimal flux distributions in genome-scale metabolic models. For the theoretical analysis, we consider the linear-fractional problem directly. Most importantly, we show that the yield-optimal solution set (like the rate-optimal solution set) is determined by (yield-optimal) elementary flux vectors of the underlying metabolic model. However, yield- and rate-optimal solutions may differ from each other, and hence optimal (biomass or product) yields are not necessarily obtained at solutions with optimal (growth or synthesis) rates. Moreover, we discuss phase planes/production envelopes and yield spaces, in particular, we prove that yield spaces are convex and provide algorithms for their computation. We illustrate our findings by a small example and demonstrate their relevance for metabolic engineering with realistic models of E. coli.

CONCLUSIONS

We develop a comprehensive mathematical framework for yield optimization in metabolic models. Our theory is particularly useful for the study and rational modification of cell factories designed under given yield and/or rate requirements.

Alignment of microbial fitness with engineered product formation: obligatory coupling between acetate production and photoautotrophic growth

BACKGROUND

Microbial bioengineering has the potential to become a key contributor to the future development of human society by providing sustainable, novel, and cost-effective production pipelines. However, the sustained productivity of genetically engineered strains is often a challenge, as spontaneous non-producing mutants tend to grow faster and take over the population. Novel strategies to prevent this issue of strain instability are urgently needed.

RESULTS

In this study, we propose a novel strategy applicable to all microbial production systems for which a genome-scale metabolic model is available that aligns the production of native metabolites to the formation of biomass. Based on well-established constraint-based analysis techniques such as OptKnock and FVA, we developed an in silico pipeline-FRUITS-that specifically 'Finds Reactions Usable in Tapping Side-products'. It analyses a metabolic network to identify compounds produced in anabolism that are suitable to be coupled to growth by deletion of their re-utilization pathway(s), and computes their respective biomass and product formation rates. When applied to Synechocystis sp. PCC6803, a model cyanobacterium explored for sustainable bioproduction, a total of nine target metabolites were identified. We tested our approach for one of these compounds, acetate, which is used in a wide range of industrial applications. The model-guided engineered strain shows an obligatory coupling between acetate production and photoautotrophic growth as predicted. Furthermore, the stability of acetate productivity in this strain was confirmed by performing prolonged turbidostat cultivations.

CONCLUSIONS

This work demonstrates a novel approach to stabilize the production of target compounds in cyanobacteria that culminated in the first report of a photoautotrophic growth-coupled cell factory. The method developed is generic and can easily be extended to any other modeled microbial production system.

Challenges in the Application of Synthetic Biology Toward Synthesis of Commodity Products by Cyanobacteria via "Direct Conversion"

Cyanobacterial direct conversion of CO₂ to several commodity chemicals has been recognized as a potential contributor to support the much-needed sustainable development of human societies. However, the feasibility of this "green conversion" hinders on our ability to overcome the hurdles presented by the natural evolvability of microbes. The latter may result in the genetic instability of engineered cyanobacterial strains leading to impaired productivity. This challenge is general to any "cell factory" approach in which the cells grow for multiple generations, and based on several studies carried out in different microbial hosts, we could identify that three distinct strategies have been proposed to tackle it. These are (1) to reduce microbial evolvability by decreasing the native mutation rate, (2) to align product formation with cell growth/fitness, and, paradoxically, (3) to efficiently reallocate cellular resources to product formation by uncoupling it from growth. The implementation of either of these strategies requires an advanced synthetic biology toolkit. Here, we review the existing methods available for cyanobacteria and identify areas of focus in which specific developments are still needed. Furthermore, we discuss how potentially stabilizing strategies may be used in combination leading to further increases of productivity while ensuring the stability of the cyanobacterial-based direct conversion process.

Growth-coupled overproduction is feasible for almost all metabolites in five major production organisms

Computational modelling of metabolic networks has become an established procedure in the metabolic engineering of production strains. One key principle that is frequently used to guide the rational design of microbial cell factories is the stoichiometric coupling of growth and product synthesis, which makes production of the desired compound obligatory for growth. Here we show that the coupling of growth and production is feasible under appropriate conditions for almost all metabolites in genome-scale metabolic models of five major production organisms. These organisms comprise eukaryotes and prokaryotes as well as heterotrophic and photoautotrophic organisms, which shows that growth coupling as a strain design principle has a wide applicability. The feasibility of coupling is proven by calculating appropriate reaction knockouts, which enforce the coupling behaviour. The study presented here is the most comprehensive computational investigation of growth-coupled production so far and its results are of fundamental importance for rational metabolic engineering.