ORIGINAL RESEARCH article

Front. Bioeng. Biotechnol., 07 January 2022

Sec. Synthetic Biology

Volume 9 - 2021 | https://doi.org/10.3389/fbioe.2021.779405

A Computational Framework to Identify Metabolic Engineering Strategies for the Co-Production of Metabolites

  • 1. Department of Biotechnology, Bhupat and Jyoti Mehta School of Biosciences, Indian Institute of Technology (IIT) Madras, Chennai, India

  • 2. Centre for Integrative Biology and Systems mEdicine (IBSE), IIT Madras, Chennai, India

  • 3. Robert Bosch Centre for Data Science and Artificial Intelligence (RBCDSAI), IIT Madras, Chennai, India

Article metrics

View details

9

Citations

5,1k

Views

1,1k

Downloads

Abstract

Microbial production of chemicals is a more sustainable alternative to traditional chemical processes. However, the shift to bioprocess is usually accompanied by a drop in economic feasibility. Co-production of more than one chemical can improve the economy of bioprocesses, enhance carbon utilization and also ensure better exploitation of resources. While a number of tools exist for in silico metabolic engineering, there is a dearth of computational tools that can co-optimize the production of multiple metabolites. In this work, we propose co-FSEOF (co-production using Flux Scanning based on Enforced Objective Flux), an algorithm designed to identify intervention strategies to co-optimize the production of a set of metabolites. Co-FSEOF can be used to identify all pairs of products that can be co-optimized with ease using a single intervention. Beyond this, it can also identify higher-order intervention strategies for a given set of metabolites. We have employed this tool on the genome-scale metabolic models of Escherichia coli and Saccharomyces cerevisiae, and identified intervention targets that can co-optimize the production of pairs of metabolites under both aerobic and anaerobic conditions. Anaerobic conditions were found to support the co-production of a higher number of metabolites when compared to aerobic conditions in both organisms. The proposed computational framework will enhance the ease of study of metabolite co-production and thereby aid the design of better bioprocesses.

1 Introduction

Recent years have seen several advances in the usage of bioprocessing to produce a wide range of chemicals (Erickson et al., 2012). Microorganisms can produce diverse and complex products from simple carbon sources. Nevertheless, there are many challenges in designing economically feasible bioprocesses. The advancements in synthetic biology have enabled the metabolic engineering of organisms to improve yield and productivity (Yadav et al., 2019). Various computational strain design algorithms have been developed to identify the genetic manipulations required to over-produce a single product (Burgard et al., 2003; Rocha et al., 2008; Yang et al., 2011). Despite the increase in yield achieved through such rational strain design, the bioprocesses are unable to compete with the traditional chemical processes in many cases (Cai and Bennett, 2011). This is due to two main reasons: 1) the cost of raw materials and 2) the maximum yield achievable for a given product in a given organism and environment is limited by the number of genetic manipulations that can be successfully implemented in a single strain (Silva et al., 2012). The former issue can be reduced by using agricultural waste as feedstock instead of a synthetic nutrient medium. The latter can be overcome by co-producing multiple products in the same bioprocess (da Silva et al., 2014).

Co-production equips us to exploit the system in a better fashion and produce more valuable products from the same raw materials. A high-value, low-volume chemical can be co-produced with a low-value, high-volume product in order to increase the economic feasibility, as in the case of riboflavin and butanol, respectively (Cai and Bennett, 2011). Co-production is essential when a cocktail of metabolites need to be produced together, rather than a single metabolite, as in the case of biofuels and fatty acids (Xin et al., 2018). A mixture of different alcohols or fatty acids of varying chain length need to be co-optimized in such cases. It can also balance carbon metabolism, as in the case of uridine, and acetoin (Fan et al., 2018). High carbon inflow towards uridine causes excess production of acetate, which hampers the growth of the organism. Conversion of acetate to acetoin prevents over-acidification of the nutrient medium and thereby improves growth and uridine production. There are many studies that have successfully achieved co-production of a variety of products with/without genetic manipulation of the organisms. Polyhydroxyalkanoates are a common class of metabolites that are co-produced with other metabolites (Li et al., 2017; Kumar and Kim, 2018; Yadav et al., 2021). Butanol and hydrogen have been co-produced in Clostridium beijerinckii (Zhang et al., 2021), and ethanol and xylitol have been co-produced in Candida tropicalis (Raj and Krishnan, 2020; de Souza Queiroz et al., 2021). The carbon source, nutrient medium, pH etc., are optimized in such cases to improve the yield of metabolites. Metabolic engineering can further expand the number of products that are co-produced and also improve their yield significantly. Multiple metabolites like ethanol, isopropanol, butanol and 2,3- butanediol have been co-produced by optimizing the acetone-butanol-ethanol (ABE) fermentation pathway in Clostridium acetobutylicum (Collas et al., 2012). Nisin and 3-phenyllactic acid, two antimicrobial agents, have been co-produced in Lactococcus lactis through genetic manipulation (Julien-Laferrière et al., 2016). Non-native metabolites can also be co-produced with other metabolites, as in the case of butanol and riboflavin, by engineering the heterologous pathway in C. acetobutylicum (Cai and Bennett, 2011).

Although many strain design algorithms have been successfully employed for metabolically engineering organisms to optimize a single product (Pharkya et al., 2003; Kumelj et al., 2019), few studies have applied it for co-production. The studies listed above only use existing literature and readily apparent deletion targets to achieve co-production. This limits the robustness of the bioprocesses that are designed. There is a lack of algorithms that can be easily applied to study co-production. In this study, we have extended the Flux Scanning based on Enforced Objective Flux (FSEOF) (Choi et al., 2010) algorithm to study co-production. Further, while deletion targets can be obtained for metabolites independently using existing algorithms like OptKnock (Burgard et al., 2003), OptGene (Rocha et al., 2008), there are very few algorithms that can identify amplification targets (Ranganathan et al., 2010). In order to identify amplification targets in addition to knock-out targets for co-optimizing a set of metabolites, we propose a new methodology, co-FSEOF, adapting the FSEOF algorithm. Co-FSEOF has a simple computational framework that can be easily modified, and it also provides the entire set of potential intervention strategies in a single run while many algorithms are sequential, returning one intervention target per run. The utility of the potential intervention strategies obtained was further assessed using Flux Variability Analysis (FVA). We applied co-FSEOF to evaluate all possible pairs of secretory metabolites in Escherichia coli and Saccharomyces cerevisiae. The different pairs of metabolites that can be co-produced through a single reaction deletion or amplification were obtained. This analysis helps us choose favorable pairs of metabolites for which higher-order intervention strategies can be obtained. We have demonstrated this by identifying the amplification targets, knock-out targets and mixed intervention strategies of size up to three to co-optimize the production of isobutanol and succinic acid in S. cerevisiae. Higher-order intervention strategies were able to achieve better yield with very little reduction in growth rate. Overall, our analyses provide an overall picture of the biosynthetic capabilities of an organism, particularly highlighting key interdependencies in metabolism.

2 Methods

2.1 Flux Balance Analysis

FBA is a widely used steady-state constraint-based modelling approach to predict the metabolic capabilities of a variety of organisms (Varma and Palsson, 1994; Kauffman et al., 2003; Orth et al., 2010). The metabolic network of an organism, which comprises all reactions known to occur in the organism, is represented as a stoichiometric matrix , of size , where is the number of metabolites, and is the number of reactions (Figure 1A). The entries in the th column of represent the stoichiometric coefficients of the metabolites that participate in the th reaction. The minimum and maximum values of flux that any reaction can assume are constrained by the lower and upper bounds, respectively. The flux through a reaction under a given set of conditions, at steady-state, is calculated by solving a linear programming (LP) problem. The LP problem is formulated as:where is a vector of weights denoting the contribution of each of the reactions to the objective function, is the vector of metabolic fluxes, and are vectors representing the lower and upper bounds for the reaction fluxes, respectively (Orth et al., 2010).

FIGURE 1

2.2 Flux Variability Analysis

FVA is used to identify the range of fluxes of each reaction that still satisfy the constraints, where two optimization problems are solved for each flux of interest.where is the vector of metabolic fluxes, and are the maximum and minimum values of fluxes, respectively for each reaction flux (Gudmundsson and Thiele, 2010).

2.3 Flux Scanning Based on Enforced Objective Flux

FSEOF (Choi et al., 2010) is a method used to identify potential reaction deletion and amplification targets in metabolic networks by observing the change in the reaction fluxes when the system moves from the wild-type flux of the target product to the theoretical maximum flux of the product (Figure 1B). The maximum biomass and maximum product fluxes are obtained by performing FBA with the biomass reaction and the exchange reaction of the product as the objective, respectively. The flux of the target reaction, is pinned to The change in the flux of a reaction, , is studied as the product flux, is increased, and it is classified as a potential deletion or amplification target based on the decrease or increase in its flux, respectively. The reactions that undergo no change or oscillations in the fluxes are discarded from the set of potential intervention strategies. The set of potential intervention strategies obtained are assessed by simulating each intervention and performing FVA on the mutant.

2.4 Co-FSEOF: Co-Optimization of Metabolites

The Genome-Scale Metabolic Models (GSMMs) of E. coli iML1515 and S. cerevisiae iMM904 were obtained from the BiGG models database (http://bigg.ucsd.edu/) (Raj and Krishnan, 2020). The simulations were done with the following constraints on uptakes: −10 mmol/gDW/h glucose and −2 mmol/gDW/h oxygen for aerobic conditions and −10 mmol/gDW/h glucose and zero oxygen uptake for anaerobic conditions. Exchange and transport reactions were removed from the search space for FSEOF to increase the relevance of the results and to reduce the computational time. The potential intervention strategies for all secretory metabolites (metabolites that can be secreted into the medium) in the organism were obtained using FSEOF as described in Section 2.3. All possible pairs of secretory metabolites were examined for co-production by identifying their common intervention strategies which were obtained through FSEOF (Figure 1C). These common deletion or amplification targets constitute the potential intervention strategies for a set of metabolites. The implementation steps of the algorithm are summarized in Figure 2.

FIGURE 2

The reliability of the potential intervention targets is analyzed by comparing the flux values of biomass and product in the mutant with those of the wild-type organism. The mutant model is obtained by deleting the reaction in case of knock-out targets or by fixing the flux bounds of the amplification target to its theoretical maximum. FVA is performed on this mutant model with biomass as objective to obtain the range of flux values for the products and biomass. Any reaction with more than a 5% increase in the maximum product flux and less than 75% decrease in biomass flux is considered a promising intervention strategy. FVA was performed using fastFVA to reduce the computational time (Gudmundsson and Thiele, 2010).

To obtain higher-order intervention strategies, the potential targets obtained earlier for a given set of metabolites were combined, and all possible combinations of intervention strategies of a certain size (up to three) were evaluated using FVA. The score for each product (total number of products, ) and the overall score are calculated aswhere and are the mutant and wild-type maximum fluxes of the exchange reaction of the product , and and are the mutant and wild-type fluxes of the biomass reaction. denotes the score for the individual product while denotes the cumulative score for the set of metabolites. All simulations were performed in MATLAB R2018a (MathWorks Inc., United States) using the COBRA Toolbox v3.0 (Heirendt et al., 2019) and IBM ILOG CPLEX 12.8 as the linear programming solver.

3 Results

Metabolic engineering strategies for the co-production of all pairs of secreted metabolites in E. coli and S. cerevisiae were obtained using co-FSEOF as described in Section 2.4. We identified the intervention strategies required to optimize the co-production of metabolites in both aerobic and anaerobic conditions. Anaerobic conditions favor the co-production of more pairs of metabolites when compared to aerobic conditions. The intervention strategy for each pair of metabolites is scored as in Section 2.4. The best intervention strategy can be chosen using the . In cases where one metabolite might be favored over the others, the individual scores, can be used to choose the best intervention strategies. Some of the intervention strategies obtained have been successfully verified through experimental studies in literature. This shows the credibility of the intervention strategies obtained. We discuss a few of the industrially significant metabolites and their intervention strategies, along with supporting literature. We also propose many other intervention strategies, which form a ready short-list for experimental validation. We were able to identify other hitherto unexplored intervention strategies, which may be better alternatives to those in existing literature, further demonstrating the utility of the algorithm.

3.1 Co-Production in Escherichia coli

E. coli is one of the well-studied model organisms and has high-quality GSMMs available. The latest GSMM, iML1515 (Monk et al., 2017), was used in this study, and the co-production of 337 secretory metabolites was studied in both aerobic and anaerobic conditions.

3.1.1 Aerobic Fermentation

Co-production of all pairs of metabolites was studied in E. coli, using co-FSEOF and FVA as described in Section 2.3, 2.4. Out of pairs of secretory metabolites, only 237 could be successfully overproduced through deletion or amplification of a single reaction. The intervention strategies for a few industrially significant pairs of metabolites are listed in Table 1. One of the important pairs of metabolites that can be easily co-produced is L-lysine, a food additive and drug additive and cadaverine, which is essential for polyamide production. co-FSEOF was able to identify several reactions from the diaminopimelate pathway (DAP), which can be over-expressed to co-produce L-lysine and cadaverine. An experimental study by Xu et al. (2019) demonstrates the effect of engineering the DAP pathway in E. coli for the production of L-lysine. This indicates the reliability of the results obtained through our computational approach. Another significant result is the co-production of succinate and ethanol through the amplification of glyceraldehyde-3-phosphate dehydrogenase. Other studies have also successfully co-produced ethanol and succinate by other genetic manipulations (Liang et al., 2019).

TABLE 1

#Product AWT flux AProduct BWT flux BInterventionMutant product flux AMutant product flux BMutant biomass fluxScore AScore BScore A + BKO/Amp
1L-lysine*1,5-Diamino pentane*Diaminopimelate decarboxylase5.815.810.228.848.8417.68Amp
Diaminopimelate epimerase5.815.810.228.848.8417.68Amp
Dihydrodipicolinate reductase5.815.810.228.848.8417.68Amp
Dihydrodipicolinate synthase5.815.810.228.848.8417.68Amp
Succinyl-diaminopimelate desuccinylase5.815.810.228.848.8417.68Amp
Tetrahydrodipicolinate succinylase5.815.810.228.848.8417.68Amp
2Succinate*Ethanol*Glyceraldehyde-3-phosphate dehydrogenase12.8615.840.2219.5624.1043.67Amp
3Spermidine*5-Methylthio-D-ribose*Adenosylmethionine decarboxylase2.012.010.223.063.066.12Amp
Methylthioadenosine nucleosidase2.012.010.223.063.066.12Amp
Spermidine synthase2.012.010.223.063.066.12Amp
4Xanthine*D-Lactate*Glyceraldehyde-3-phosphate dehydrogenase8.9015.840.2213.5524.1037.65Amp
5Glycine*L-Asparagine*Glyceraldehyde-3-phosphate dehydrogenase20.9111.820.2231.8017.9949.79Amp
6Fe-enterobactin*Enterobactin*2,3-Dihydro-2,3-dihydroxybenzoate dehydrogenase1.291.290.221.971.973.93Amp
Isochorismatase1.291.290.221.971.973.93Amp
7Pyruvate*L-Asparagine*Glyceraldehyde-3-phosphate dehydrogenase18.0111.820.2227.3917.9945.38Amp

Intervention strategies for co-production of pairs of metabolites in E. coli under aerobic conditions.

WT, wild type; *, less than 10−5 mmol/gDW/h.

3.1.2 Anaerobic Fermentation

Anaerobic conditions support the co-production of more metabolites when compared to aerobic conditions. More than 1,000 pairs of metabolites can be co-produced, out of which few are listed in Table 2. L-lysine and cadaverine can be co-produced under anaerobic conditions too. But the maximum flux achievable is lower when compared to aerobic conditions. The yield of metabolites like acetate, formate, and hexanoate can be co-optimized by deleting acetaldehyde dehydrogenase or alcohol dehydrogenase. We also found that succinate and lactate can be co-produced by the knock-out of pyruvate formate lyase. The effect of deletion of pflB gene encoding pyruvate formate lyase has been experimentally verified in E. coli for succinate production (Zhang et al., 2009) and lactate production (Utrilla et al., 2009) through separate studies. This shows that there are multiple co-production strategies available in the existing literature that can be easily utilized to design an efficient process.

TABLE 2

#Product AWT flux AProduct BWT flux BInterventionMutant product flux AMutant product flux BMutant biomass fluxScore AScore BScore A + BKO/Amp
1Acetate8.83Formate18.22Acetaldehyde dehydrogenase18.1836.790.12237.95472.31710.26KO
Alcohol dehydrogenase18.1836.790.12237.95472.31710.26KO
2Succinate0.05D-Lactate3.76x10−4Glyceraldehyde-3-phosphate dehydrogenase14.9819.410.04126.78164.92291.71Amp
Acetaldehyde dehydrogenase9.0718.120.12229.50461.04690.55KO
Alcohol dehydrogenase9.0718.120.12229.50461.04690.55KO
Pyruvate formate lyase0.4917.760.1210.57426.65437.22KO
3L-Alanine*Xanthine*Glyceraldehyde-3-phosphate dehydrogenase13.560.970.04115.188.23123.41Amp
4Spermidine*5-Methylthio-D-ribose*Adenosylmethionine decarboxylase0.640.640.045.465.4610.93Amp
Methylthioadenosine nucleosidase0.640.640.045.465.4610.93Amp
Spermidine synthase0.640.640.045.465.4610.93Amp
5L-Aspartate*L-Glutamate*Glyceraldehyde-3-phosphate dehydrogenase6.034.520.0451.1938.3989.59Amp
6L-Lysine*1,5-Diamino pentane*Diaminopimelate decarboxylase2.722.720.0423.0423.0446.07Amp
Diaminopimelate epimerase2.722.720.0423.0423.0446.07Amp
Dihydrodipicolinate synthase2.722.720.0423.0423.0446.07Amp
Succinyl-diaminopimelate desuccinylase2.722.720.0423.0423.0446.07Amp
Tetrahydrodipicolinate succinylase2.722.720.0423.0423.0446.07Amp
7Acetate8.83Hexanoate*Acetaldehyde dehydrogenase18.184.530.12237.95115.26353.21KO
Alcohol dehydrogenase18.184.530.12237.95115.26353.21KO

Intervention strategies for co-production of pairs of metabolites in E. coli under anaerobic conditions.

WT, wild type; *, less than 10−5 mmol/gDW/h.

3.2 Co-Production in Saccharomyces cerevisiae

Another industrially relevant and well-studied model organism is S. cerevisiae. Though heterologous pathways have not been analyzed in this study, one can easily modify the GSMM and apply co-FSEOF to identify co-production strategies for heterologous metabolites. Since S. cerevisiae is a better candidate for recombinant protein production, it is essential to study co-production in yeast (Bill, 2014). It can also produce more complex metabolites when compared to E. coli and is, therefore, a favorable candidate for bio-production. The latest GSMM iMM904 (Mo et al., 2009) was used, and the ability to optimize the co-production of 164 secretory metabolites was studied in both aerobic and anaerobic conditions.

3.2.1 Aerobic Fermentation

We found that many industrially important metabolites like ethanol and L-alanine, and 4-aminobutanoate and L-serine can be co-produced in S. cerevisiae under aerobic conditions. We were also able to co-optimize isobutyl alcohol and 2-methyl propanal, which are long-chain alcohols that are used as biofuels. The deletion of pyruvate dehydrogenase increases the production of pyruvate and acetate, as shown in Table 3. Although the deletion of pyruvate dehydrogenase has not been experimentally verified as yet, a similar study has been carried out in E. coli (Moxley and Eiteman, 2021). In this study, it has been shown that the deletion of the genes encoding pyruvate dehydrogenase improves pyruvate production (Moxley and Eiteman, 2021). In addition to pyruvate dehydrogenase, co-FSEOF was able to identify several other amplification targets, which can also improve the production of pyruvate and acetate.

TABLE 3

#Product AWT flux AProduct BWT flux BInterventionMutant product flux AMutant product flux BMutant biomass fluxScore AScore BScore A + BKO/Amp
1Ethanol15.81L-Alanine1.69x10−4Sedoheptulose 1,7-bisphosphate D-glyceraldehyde-3-phosphate-lyase18.450.150.0712.220.7112.93Amp
Phosphofructokinase (s7p)18.450.150.0712.220.7112.93Amp
2Acetate3.55×10−32,3-Butanediol3.21x10−4Pyruvate dehydrogenase1.480.250.28230.9438.48269.42KO
Enolase10.7110.270.0749.6147.6197.22Amp
Fructose-bisphosphate aldolase10.7010.270.0749.6247.6497.26Amp
Glyceraldehyde-3-phosphate dehydrogenase10.7110.270.0749.6147.6197.22Amp
Triose-phosphate isomerase10.7010.270.0749.6247.6497.26Amp
3Isobutyl alcohol*Succinate*Pyruvate decarboxylase7.635.450.1979.2956.66135.95KO
Enolase9.5612.970.0744.3060.10104.40Amp
Glyceraldehyde-3-phosphate dehydrogenase9.5612.970.0744.3060.13104.43Amp
Triose-phosphate isomerase9.5612.970.0744.3360.17104.50Amp
42-Methylpropanal*Isobutyl alcohol*3-Methyl-2-oxobutanoate decarboxylase6.509.360.0730.1443.3973.54Amp
Acetolactate synthase mitochondrial6.509.360.0730.1443.4073.54Amp
Dihydroxy acid dehydratase 2,3-dihydroxy-3-methylbutanoate mitochondrial6.509.360.0730.1443.4073.54Amp
Enolase8.959.560.0741.4744.3085.77Amp
Glyceraldehyde-3-phosphate dehydrogenase8.959.560.0741.4744.3085.77Amp
Acetohydroxy acid isomeroreductase mitochondrial6.509.360.0730.1443.4073.54Amp
Triose-phosphate isomerase8.949.560.0741.4944.3385.82Amp
5L-Glutamate*2-Oxoglutarate*Citrate synthase2.931.950.0713.589.0522.63Amp
Enolase4.504.040.0720.8518.7339.58Amp
Fructose-bisphosphate aldolase4.504.040.0720.8618.7339.60Amp
Glyceraldehyde-3-phosphate dehydrogenase4.504.040.0720.8518.7339.58Amp
Isocitrate dehydrogenase2.992.000.0713.879.2523.12Amp
Triose-phosphate isomerase4.504.040.0720.8618.7339.60Amp
6Acetate3.55×10−3Pyruvate3.24x10−4Pyruvate dehydrogenase1.480.250.28230.9438.48269.42KO
Aspartate-semialdehyde dehydrogenase8.7911.850.0741.1455.4996.63Amp
Enolase10.7114.420.0749.6166.83116.44Amp
Fructose-bisphosphate aldolase10.7014.410.0749.6266.85116.47Amp
Glyceraldehyde-3-phosphate dehydrogenase10.7114.420.0749.6166.83116.44Amp
Triose-phosphate isomerase10.7014.410.0749.6266.85116.47Amp
74-Aminobutanoate*L-Serine*Glyceraldehyde-3-phosphate dehydrogenase4.776.980.0722.1032.3454.44Amp
Triose-phosphate isomerase4.776.970.0722.1132.3554.45Amp
8L-Alanine1.69x10−4L-Cysteine*Glucose-6-phosphate dehydrogenase10.601.200.0749.115.5654.67Amp
Phosphogluconate dehydrogenase10.601.200.0749.115.5654.67Amp
6-phosphogluconolactonase10.601.200.0749.115.5654.67Amp
Ribulose-5-phosphate-3-epimerase10.601.200.0749.115.5654.67Amp
Transketolase10.601.200.0749.125.5654.68Amp
Ribose-5-phosphate isomerase18.162.210.0885.4410.4095.84Amp
9sn-Glycero-3-phosphocholine*L-Methionine*Methionine synthase0.022.190.070.1110.1710.27Amp
5,10-Methylene-tetrahydrofolate reductase0.022.190.070.1110.1710.27Amp
Ribose-5-phosphate isomerase0.661.710.083.118.0311.14Amp
102-Methylbutyl acetate*2-Methyl-1-butanol*Ribose-5-phosphate isomerase3.004.420.0814.1320.7834.92Amp

Intervention strategies for co-production of pairs of metabolites in S. cerevisiae under aerobic conditions.

WT, wild type; *, less than 10−5 mmol/gDW/h.

3.2.2 Anaerobic Fermentation

As in the case of E. coli, anaerobic fermentation enables the co-production of more pairs of metabolites in S. cerevisiae when compared to aerobic fermentation. 2-methyl-1-butanol, which is an important solvent used in the manufacture of pesticides and paints and isobutyl alcohol, which is a biofuel, can be co-produced by the amplification of malic enzyme, as shown in Table 4. Formate, which is used in dyeing and printing, can be co-produced with spermidine, a metabolite increasingly studied for its anti-ageing properties (Minois, 2014), through the amplification of a number of reactions. These strategies not only include readily apparent reactions that are involved in spermidine synthesis like spermidine synthase and adenosylmethionine decarboxylase but also provide some non-intuitive strategies like the amplification of aspartate transaminase or 2-keto-4-methylthiobutyrate transaminase. We also found that the deletion of pyruvate decarboxylase improves the production of succinate, isobutyl alcohol and pyruvate. The effect of deletion of pyruvate decarboxylase has been studied in S. cerevisiae, and the improvement in the production of pyruvate (van Maris et al., 2004) and succinate (Zahoor et al., 2019) has been verified experimentally in separate studies in literature.

TABLE 4

#Product AWT flux AProduct BWT flux BInterventionMutant product flux AMutant product flux BMutant biomass fluxScore AScore BScore A + BKO/Amp
12-Methyl-1-butanol*Isobutyl alcohol*Malic enzyme NADP mitochondrial0.039.640.050.2061.3061.51Amp
2Isobutyl alcohol*Pyruvate*Pyruvate decarboxylase8.735.660.1188.7457.57146.31KO
Enolase9.689.970.0561.5963.45125.03Amp
Fructose-bisphosphate aldolase9.689.960.0561.6763.47125.14Amp
Glyceraldehyde-3-phosphate dehydrogenase9.689.970.0561.6063.45125.05Amp
Triose-phosphate isomerase9.689.960.0561.6663.47125.13Amp
3Formate6.26x10−4Spermidine*Adenosylmethionine decarboxylase1.221.220.057.767.7615.52Amp
Aspartate transaminase1.211.210.057.757.7615.51Amp
2,3-Diketo-5-methylthio-1-phosphopentane degradation1.221.220.057.767.7615.52Amp
5-Methylthio-5-deoxy-D-ribulose-1-phosphate dehydratase1.221.220.057.767.7615.52Amp
5-Methylthioadenosine phosphorylase1.221.220.057.767.7615.52Amp
5-Methylthioribose-1-phosphate isomerase1.221.220.057.767.7615.52Amp
Spermidine synthase1.221.220.057.767.7615.52Amp
2-Keto-4-methylthiobutyrate transamination1.221.220.057.767.7615.52Amp
44-Amino butanoate*Isobutyl acetate*Enolase2.842.910.0518.0418.5036.53Amp
Fructose-bisphosphate aldolase2.842.900.0518.1118.5036.61Amp
Glyceraldehyde-3-phosphate dehydrogenase2.852.910.0518.1118.5036.61Amp
Triose-phosphate isomerase2.842.900.0518.1118.5036.61Amp
52-Methyl-1-butanol*Glycine*Aspartate kinase0.054.550.050.3228.9529.26Amp
Glucose-6-phosphate dehydrogenase0.084.150.050.5226.3926.91Amp
Phosphogluconate dehydrogenase0.084.150.050.5226.3926.91Amp
Homoserine dehydrogenase NADH irreversible0.054.550.050.3228.9529.26Amp
Homoserine kinase0.054.550.050.3228.9129.22Amp
6-phosphogluconolactonase0.084.150.050.5226.3926.91Amp
Ribulose-5-phosphate-3-epimerase0.084.150.050.5226.3926.91Amp
Ribose-5-phosphate isomerase3.314.010.0621.5226.0447.56Amp
Threonine synthase0.054.550.050.3228.9129.22Amp
Transketolase0.084.150.050.5226.3926.91Amp
Transketolase0.084.150.050.5226.3926.91Amp
6Isobutyl alcohol*Succinate0.67Pyruvate decarboxylase8.736.700.1188.7461.27150.02KO
7L-Glutamate*Xanthine*Fructose-bisphosphate aldolase2.810.950.0517.906.0423.94Amp
Glyceraldehyde-3-phosphate dehydrogenase2.810.950.0517.906.0423.93Amp
Triose-phosphate isomerase2.810.950.0517.906.0423.94Amp
8Sorbitol*L-Methionine*Ribose-5-phosphate isomerase5.611.280.0636.458.2944.74Amp
92,3-Butanediol*L-Serine*Fructose-bisphosphate aldolase9.304.740.0559.2130.1989.40Amp
Glyceraldehyde-3-phosphate dehydrogenase9.304.740.0559.1930.1989.37Amp
Triose-phosphate isomerase9.304.740.0559.2130.1989.40Amp

Intervention strategies for co-production of pairs of metabolites in S. cerevisiae under anaerobic conditions.

WT, wild type; *, less than 10−5 mmol/gDW/h.

3.2.3 Higher-Order Intervention StrategiesCo-Production of Isobutanol and Succinate

Higher-order intervention strategies can increase the maximum yield achievable for any product with a little difference in growth rate when compared to single interventions. But they are more cumbersome to identify, as the problem becomes time-consuming and computationally expensive. Instead of identifying higher-order targets for all metabolites in an organism, we have used the previous analysis to explore the metabolic capabilities of the organism and chose one set of metabolites to demonstrate the power of higher-order intervention strategies.

Isobutanol is a long-chain alcohol that is an attractive biofuel (Nanda et al., 2017). Succinic acid is an important metabolite essential for the production of various other products like biodegradable polymers, fatty acids, butyrolactone and tetrahydrofuran (Akhtar et al., 2014). The co-production of isobutanol and succinate has been proposed as a sustainable and economical process by Xu et al. (2018). They have discussed the development of various strains for the production of isobutanol and succinate separately. They emphasize how the co-production of isobutanol and succinate is not only of economic significance, but the high amount of carbon dioxide released from long-chain alcohol fermentation can be used for succinate production, and is hence also of ecological importance. But the article does not discuss any strategy to co-optimize the production of isobutanol and succinate.

Here we identified the higher-order intervention strategies (size up to three) for co-production of isobutanol and succinate in S. cerevisiae in aerobic conditions. More than 3,700 interventions can improve the yield of both the metabolites when compared to the wild-type strain. Table 5 lists a few examples of each type of intervention strategy obtained, which are also represented in Figure 3. Though most of the amplification and deletion targets are integral components of the target product and by-product synthesis respectively, co-FSEOF is also able to find intervention targets in distant pathways like those in pentose phosphate pathway, shikimate pathway, and nucleotide metabolism. This shows the ability of the algorithm to identify non-intuitive targets. Also, a number of the targets predicted here are also found in experimental studies reported in literature. For example, the deletion of pyruvate decarboxylase has been shown to improve the production of isobutanol by Kondo et al. (2012). Zahoor et al. (2019) have shown that both pyruvate decarboxylase deletion and fumarase deletion can increase the production of succinate. This shows the dependability of the results obtained using the algorithm.

TABLE 5

#Intervention 1Intervention 2Intervention 3A ⊕/K ⊝Mutant flux 1Mutant flux 2Biomass fluxScore AScore BScore A + B
1Glyceraldehyde-3-phosphate dehydrogenase (GAPD)Pyruvate kinase (PYK)NA⊕⊕7.198.600.2081.7397.87179.59
2Enolase (ENO)Pyruvate kinase (PYK)NA⊕⊕7.138.600.2081.0897.87178.95
3Glutamate-5-kinase (GLU5K)Phosphoglycerate dehydrogenase (PGCD)Pyruvate decarboxylase (PYRDC)⊝⊝⊝7.965.800.1979.8158.12137.93
4Fumarase (FUMm)Phosphoserine phosphatase (PSP_L)Pyruvate decarboxylase (PYRDC)⊝⊝⊝7.955.780.1979.7858.05137.84
5Aldehyde dehydrogenase (ALCD23x)Glycerol-3-phosphate dehydrogenase (G3PD1iR)NA⊕⊝2.992.520.2243.7636.7780.54
6Succinate CoA ligase ADP forming (SUCOASm)Pyruvate decarboxylase (PYRDC)NA⊕⊝7.635.450.1979.3156.69136.00
7Acetolactate synthase (ACLSm)Oxoglutarate dehydrogenase lipoamide (AKGDam)Guanylate kinase (GK2)⊕⊕⊝4.726.890.1534.5150.3884.89
8Citrate synthase (CSm)Dihydroxy-acid dehydratase (DHAD1im)Prephenate dehydrogenase (PPND)⊕⊕⊝4.736.980.1432.3147.7180.02
9Aldehyde dehydrogenase (ALCD23x)Ribonucleoside-diphosphate reductase (RNDR1)Ribulose 5-phosphate 3-epimerase (RPE)⊕⊝⊝2.992.480.2243.8436.3280.16
10Succinate CoA ligase ADP forming (SUCOASm)Phosphoserine phosphatase (PSP_L)Pyruvate decarboxylase (PYRDC)⊕⊝⊝7.955.780.1979.7858.05137.84

Higher-order intervention strategies for co-production of isobutanol and succinate in S. cerevisiae under aerobic conditions.

⊕, Amplification of reaction; ⊝, Knock-out of reaction.

FIGURE 3

4 Discussion

Chemical processes based on fossil fuels are cheaper when compared to bioprocesses, which leads to reluctance in the adoption of sustainable bioprocesses in industries. To improve the economic feasibility of a bioprocess, we can optimize the process variables and/or genetically engineer the microbes (Dzurendova et al., 2020). Even then, in some cases, the bioprocess might be less lucrative when compared to their chemical counterparts (Zhang et al., 2021). In such cases, we can co-produce multiple metabolites to improve the economic feasibility and efficiency of a bioprocess. For example, in the case of biofuels and fatty acids, we need to design processes that can support the production of multiple metabolites of similar nature (Xin et al., 2018). Co-production also ensures better utilization of microbial capabilities, and better balance in the carbon metabolism (Fan et al., 2018). While there are multiple computational tools and algorithms to identify intervention strategies for a single product, there is a lack of readily appliable algorithms for co-production. As a result, almost no co-production study in existing literature was found to use computational algorithms to aid rational strain design. All of the studies rely on previous findings or readily apparent strategies to achieve co-production. This limits the intervention strategies designed.

In this study, we present co-FSEOF, by adapting the effective FSEOF algorithm to study the co-optimization of a set of metabolites. FSEOF is a well-established constraint-based modelling algorithm, which has been used to reliably predict metabolic engineering strategies for a variety of systems (Choi et al., 2010; Boghigian et al., 2012; Badri et al., 2019; Srinivasan et al., 2019). It has a simple and efficient framework and can identify both deletion and amplification targets. Flux Coupling Analysis (FCA) is yet another interesting algorithm that can identify which metabolites can be coupled together. But very few metabolites are innately coupled without interventions. Moreover, it excludes all the reactions that do not carry flux under a given set of conditions from the analysis. This affects the applicability of FCA because these reactions, though not coupled, can have an effect on target production in the presence of other interventions. co-FSEOF is able to identify more combinations of products that can be co-produced and also provides a wider range of intervention strategies for a given set of metabolites.

Using co-FSEOF, we examined the co-production of multiple pairs of metabolites, and both deletion and amplification targets were obtained in E. coli and S. cerevisiae under both aerobic and anaerobic conditions. Anaerobic fermentation enabled the co-production of a higher number of metabolites when compared to aerobic fermentation in both organisms. This could be due to the incomplete respiration in the absence of oxygen that leads to the formation of multiple by-products. Also, S. cerevisiae produces more industrially significant metabolites when compared to E. coli. Some of these proposed intervention strategies have been verified experimentally by other studies in literature, as mentioned in Section 3. This shows the efficacy of the algorithm in furnishing reliable targets. In addition to readily apparent intervention strategies, co-FSEOF also provides non-intuitive intervention strategies that are present in auxiliary biochemical pathways (as discussed in Section 3.2.2, 3.2.3).

The co-optimization analysis for all possible pairs of metabolites in the network is intended to be exploratory in order to give a larger picture of the metabolic capabilities of the organism. This analysis showed that around 200 pairs could be co-optimized in E. coli under aerobic conditions, and around 1,000 pairs of metabolites could be co-optimized in the other cases. An important class of metabolites observed in the analysis are alcohols such as ethanol, isobutanol, and 2,3-butanediol, which can be co-produced using various interventions. Co-optimization can thus enable the efficient production of biofuels (Collas et al., 2012). We can explore the metabolic capabilities of the organism to identify all possible pairs of metabolites that can be co-produced. Following this, important and commercially valuable pairs of metabolites can be further studied to obtain higher-order intervention strategies, as shown in Section 3.2.3. Exploring the higher-order strategies can expand the efficiency of the intervention targets obtained. Since the evaluation of higher-order intervention strategies is laborious and computationally expensive, we have limited the size to a maximum of three manipulations at a time. It not only enhances yield, but also provides alternate routes to achieve a similar yield. The advantageous strategies can be chosen based on the ease of manipulation in an experimental setup in such cases.

The evaluation of the results is carried out using FVA, which ensures the robustness of the targets obtained. While FBA provides one optimal solution from the solution space, FVA gives the entire range of values the flux can take up. This is a significant difference that sets co-FSEOF apart from other existing algorithms like OptKnock (Burgard et al., 2003) and OptReg (Pharkya and Maranas, 2006). Also, the algorithm validates and returns all the intervention strategies in a single run, contrary to the existing algorithms, most of which are sequential and require a separate run for each strategy obtained. The set of intervention strategies validated through FVA can be short-listed for experimental verification using the scores. The can be used to compare the effectiveness of different intervention strategies. The score is designed to incorporate both increase in product flux and decrease in biomass flux, so that both biomass and product production are favored in the mutant. Moreover, co-FSEOF uses biomass as the objective throughout the analysis, and the effect of the intervention strategies on product synthesis is studied when the organism optimizes growth. Thereby, the intervention strategies result in the co-optimization of product production and biomass formation. An intervention strategy with better ensures better product synthesis along with good biomass formation. If one product is more favored than the others economically or otherwise, we can use the individual scores to choose the appropriate strategy for the process that is formulated. The products can be chosen based on their economic value, or ease of co-production. One drawback of co-production is the cost associated with downstream processing. But this can be overcome by choosing easily separable products or choosing metabolites such that one is accumulated in the cell and one is secreted out, as in the case of polyhydroxy butyrate and succinate, respectively (Kang et al., 2010). However, this problem does not occur in the case of biofuels where the alcohol mixture is optimized for and therefore does not require extensive separation of the products. Co-FSEOF not only identifies intervention strategies for co-production of a given set of metabolites, but also allows us to explore the different combinations of products that can be co-produced in an easy and efficient manner.

5 Conclusion

Co-production can open new avenues for the sustainable production of chemicals. Designing bioprocesses for co-production using laboratory experiments alone is cumbersome and can result in sub-optimal strategies. co-FSEOF empowers us to explore and exploit microbial systems in a better fashion. It can be used to computationally study and optimize co-production by identifying intervention strategies for multiple metabolites and thereby improve the efficiency of bioprocesses. It should be noted that the co-optimization analysis was limited to pairs of metabolites to reduce the computational time. Nevertheless, co-FSEOF can be easily extended to co-optimize a more extensive set of metabolites. To conclude, this study can be used to identify various genetic manipulations that can co-optimize a set of products, which might be challenging to achieve through pure experimentation. It provides a novel and critical approach to study co-production computationally. We hope this study will aid the design and development of more sustainable bioprocesses.

Statements

Data availability statement

Publicly available datasets were analyzed in this study. This data can be found here: http://bigg.ucsd.edu/. All models used in this work and the codes used for our analysis are available at https://github.com/RamanLab/co-FSEOF.

Author contributions

Conceptualization—LR, KR; Funding—KR; Writing—original draft—LR; Writing—review and editing—LR, KR. Signed by all authors as follows: Lavanya Raajaraam (LR), Karthik Raman (KR)* *Corresponding author.

Acknowledgments

LR acknowledges the HTRA fellowship from the Ministry of Education, Government of India. KR acknowledges support from the Science and Engineering Board (SERB) MATRICS Grant MTR/2020/000490.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

  • 1

    AkhtarJ.IdrisA.Abd AzizR. (2014). Recent Advances in Production of Succinic Acid from Lignocellulosic Biomass. Appl. Microbiol. Biotechnol.98 (3), 9871000. 10.1007/s00253-013-5319-6

  • 2

    BadriA.RamanK.JayaramanG. (2019). Uncovering Novel Pathways for Enhancing Hyaluronan Synthesis in Recombinant Lactococcus Lactis: Genome-Scale Metabolic Modeling and Experimental Validation. Processes7 (6), 343. 10.3390/pr7060343

  • 3

    BillR. M. (2014). Playing Catch-Up with Escherichia coli: Using Yeast to Increase success Rates in Recombinant Protein Production Experiments. Front. Microbiol.5, 85. 10.3389/fmicb.2014.00085

  • 4

    BoghigianB. A.ArmandoJ.SalasD.PfeiferB. A. (2012). Computational Identification of Gene Over-expression Targets for Metabolic Engineering of Taxadiene Production. Appl. Microbiol. Biotechnol.93 (5), 20632073. 10.1007/s00253-011-3725-1

  • 5

    BurgardA. P.PharkyaP.MaranasC. D. (2003). Optknock: A Bilevel Programming Framework for Identifying Gene Knockout Strategies for Microbial Strain Optimization. Biotechnol. Bioeng.84 (6), 647657. 10.1002/bit.10803

  • 6

    CaiX.BennettG. N. (2011). Improving the Clostridium Acetobutylicum Butanol Fermentation by Engineering the Strain for Co-production of Riboflavin. J. Ind. Microbiol. Biotechnol.38 (8), 10131025. 10.1007/s10295-010-0875-6

  • 7

    ChoiH. S.LeeS. Y.KimT. Y.WooH. M. (2010). In Silico Identification of Gene Amplification Targets for Improvement of Lycopene Production. Appl. Environ. Microbiol.76 (10), 30973105. 10.1128/aem.00115-10

  • 8

    CollasF.KuitW.ClémentB.MarchalR.López-ContrerasA. M.MonotF. (2012). Simultaneous Production of Isopropanol, Butanol, Ethanol and 2,3-butanediol by Clostridium Acetobutylicum ATCC 824 Engineered Strains. AMB Expr.2 (1), 45. 10.1186/2191-0855-2-45

  • 9

    da SilvaT. L.GouveiaL.ReisA. (2014). Integrated Microbial Processes for Biofuels and High Value-Added Products: the Way to Improve the Cost Effectiveness of Biofuel Production. Appl. Microbiol. Biotechnol.98 (3), 10431053. 10.1007/s00253-013-5389-5

  • 10

    de Souza QueirozS.JofreF. M.dos SantosH. A.Hernández-PérezA. F.de Almeida FelipeM. D. G. (2021). Xylitol and Ethanol Co-production from Sugarcane Bagasse and Straw Hemicellulosic Hydrolysate Supplemented with Molasses. Biomass Conv. Bioref.110. 10.1007/s13399-021-01493-y

  • 11

    DzurendovaS.ZimmermannB.KohlerA.TafintsevaV.SlanyO.CertikM.et al (2020). Microcultivation and FTIR Spectroscopy-Based Screening Revealed a Nutrient-Induced Co-production of High-Value Metabolites in Oleaginous Mucoromycota Fungi. Plos One15 (6), e0234870. 10.1371/journal.pone.0234870

  • 12

    EricksonB.NelsonWintersP. (2012). Perspective on Opportunities in Industrial Biotechnology in Renewable Chemicals. Biotechnol. J.7 (2), 176185. 10.1002/biot.201100069

  • 13

    FanX.WuH.JiaZ.LiG.LiQ.ChenN.et al (2018). Metabolic Engineering of Bacillus Subtilis for the Co-production of Uridine and Acetoin. Appl. Microbiol. Biotechnol.102 (20), 87538762. 10.1007/s00253-018-9316-7

  • 14

    GudmundssonS.ThieleI. (2010). Computationally Efficient Flux Variability Analysis. BMC Bioinformatics11 (1), 489. 10.1186/1471-2105-11-489

  • 15

    HeirendtL.ArreckxS.PfauT.MendozaS. N.RichelleA.HeinkenA.et al (2019). Creation and Analysis of Biochemical Constraint-Based Models Using the COBRA Toolbox v.3.0. Nat. Protoc.14 (3), 639702. 10.1038/s41596-018-0098-2

  • 16

    Julien-LaferrièreA.BulteauL.ParrotD.Marchetti-SpaccamelaA.StougieL.VingaS.et al (2016). A Combinatorial Algorithm for Microbial Consortia Synthetic Design. Sci. Rep.6 (1), 29182. 10.1038/srep29182

  • 17

    KangZ.GaoC.WangQ.LiuH.QiQ. (2010). A Novel Strategy for Succinate and Polyhydroxybutyrate Co-production in Escherichia coli. Bioresour. Technol.101 (19), 76757678. 10.1016/j.biortech.2010.04.084

  • 18

    KauffmanK. J.PrakashP.EdwardsJ. S. (2003). Advances in Flux Balance Analysis. Curr. Opin. Biotechnol.14 (5), 491496. 10.1016/j.copbio.2003.08.001

  • 19

    KondoT.TezukaH.IshiiJ.MatsudaF.OginoC.KondoA. (2012). Genetic Engineering to Enhance the Ehrlich Pathway and Alter Carbon Flux for Increased Isobutanol Production from Glucose by Saccharomyces cerevisiae. J. Biotechnol.159 (1), 3237. 10.1016/j.jbiotec.2012.01.022

  • 20

    KumarP.KimB. S. (2018). Valorization of Polyhydroxyalkanoates Production Process by Co-synthesis of Value-Added Products. Bioresour. Technol.269, 544556. 10.1016/j.biortech.2018.08.120

  • 21

    KumeljT.SulheimS.WentzelA.AlmaasE. (2019). Predicting Strain Engineering Strategies Using iKS1317: A Genome-Scale Metabolic Model of Streptomyces Coelicolor. Biotechnol. J.14 (4), e1800180. 10.1002/biot.201800180

  • 22

    LiT.ElhadiD.ChenG.-Q. (2017). Co-Production of Microbial Polyhydroxyalkanoates with Other Chemicals. Metab. Eng.43, 2936. 10.1016/j.ymben.2017.07.007

  • 23

    LiangC. Y.XuJ. L.XuH. J.QiW.ZhangY.LuoW.et al (2019). Gene Cloning and Characterization of an Organic Solvent-Stimulated β-glucosidase and its Application for the Co-production of Ethanol and Succinic Acid. Cellulose26 (15), 82378248. 10.1007/s10570-019-02477-y

  • 24

    MinoisN. (2014). Molecular Basis of the 'Anti-Aging' Effect of Spermidine and Other Natural Polyamines - A Mini-Review. Gerontology60 (4), 319326. 10.1159/000356748

  • 25

    MoM. L.PalssonB. Ø.HerrgårdM. J. (2009). Connecting Extracellular Metabolomic Measurements to Intracellular Flux States in Yeast. BMC Syst. Biol.3 (1), 37. 10.1186/1752-0509-3-37

  • 26

    MonkJ. M.LloydC. J.BrunkE.MihN.SastryA.KingZ.et al (2017). iML1515, a Knowledgebase that Computes Escherichia coli Traits. Nat. Biotechnol.35 (10), 904908. 10.1038/nbt.3956

  • 27

    MoxleyW. C.EitemanM. A. (2021). Pyruvate Production by Escherichia coli by Use of Pyruvate Dehydrogenase Variants. Appl. Environ. Microbiol.87 (13), e0048721. 10.1128/aem.00487-21

  • 28

    NandaS.Golemi-KotraD.McDermottJ. C.DalaiA. K.GökalpI.KozinskiJ. A. (2017). Fermentative Production of Butanol: Perspectives on Synthetic Biology. New Biotechnol.37, 210221. 10.1016/j.nbt.2017.02.006

  • 29

    OrthJ. D.ThieleI.PalssonB. Ø. (2010). What Is Flux Balance Analysis?Nat. Biotechnol.28 (3), 245248. 10.1038/nbt.1614

  • 30

    PharkyaP.MaranasC. D. (2006). An Optimization Framework for Identifying Reaction Activation/inhibition or Elimination Candidates for Overproduction in Microbial Systems. Metab. Eng.8 (1), 113. 10.1016/j.ymben.2005.08.003

  • 31

    PharkyaP.BurgardA. P.MaranasC. D. (2003). Exploring the Overproduction of Amino Acids Using the Bilevel Optimization Framework OptKnock. Biotechnol. Bioeng.84 (7), 887899. 10.1002/bit.10857

  • 32

    RajK.KrishnanC. (2020). Improved Co-production of Ethanol and Xylitol from Low-Temperature Aqueous Ammonia Pretreated Sugarcane Bagasse Using Two-Stage High Solids Enzymatic Hydrolysis and Candida tropicalis. Renew. Energ.153, 392403. 10.1016/j.renene.2020.02.042

  • 33

    RanganathanS.SuthersP. F.MaranasC. D. (2010). OptForce: An Optimization Procedure for Identifying All Genetic Manipulations Leading to Targeted Overproductions. Plos Comput. Biol.6 (4), e1000744. 10.1371/journal.pcbi.1000744

  • 34

    RochaI.MaiaP.RochaM.FerreiraE. C. (2008). OptGene – A Framework for in Silico Metabolic Engineering in 10th International Conference on Chemical and Biological Engineering (Portugal: University of Minho), 218219.

  • 35

    SilvaF.QueirozJ. A.DominguesF. C. (2012). Evaluating Metabolic Stress and Plasmid Stability in Plasmid DNA Production by Escherichia coli. Biotechnol. Adv.30 (3), 691708. 10.1016/j.biotechadv.2011.12.005

  • 36

    SrinivasanA.SV.RamanK.SrivastavaS. (2019). Rational Metabolic Engineering for Enhanced Alpha-Tocopherol Production in Helianthus Annuus Cell Culture. Biochem. Eng. J.151, 107256. 10.1016/j.bej.2019.107256

  • 37

    UtrillaJ.GossetG.MartinezA. (2009). ATP Limitation in a Pyruvate Formate Lyase Mutant of Escherichia coli MG1655 Increases Glycolytic Flux to D-Lactate. J. Ind. Microbiol. Biotechnol.36 (8), 10571062. 10.1007/s10295-009-0589-9

  • 38

    van MarisA. J. A.GeertmanJ.-M. A.VermeulenA.GroothuizenM. K.WinklerA. A.PiperM. D. W.et al (2004). Directed Evolution of Pyruvate Decarboxylase-Negative Saccharomyces cerevisiae, Yielding a C 2 -Independent, Glucose-Tolerant, and Pyruvate-Hyperproducing Yeast. Appl. Environ. Microbiol.70 (1), 159166. 10.1128/aem.70.1.159-166.2004

  • 39

    VarmaA.PalssonB. O. (1994). Metabolic Flux Balancing: Basic Concepts, Scientific and Practical Use. Nat. Biotechnol.12 (10), 994998. 10.1038/nbt1094-994

  • 40

    XinF.DongW.JiangY.MaJ.ZhangW.WuH.et al (2018). Recent Advances on Conversion and Co-production of Acetone-Butanol-Ethanol into High Value-Added Bioproducts. Crit. Rev. Biotechnol.38 (4), 529540. 10.1080/07388551.2017.1376309

  • 41

    XuC.ZhangJ.ZhangY.GuoY.XuH.XuJ. (2018). Long Chain Alcohol and Succinic Acid Co-production Process Based on Full Utilization of Lignocellulosic Materials. Curr. Opin. Green Sustain. Chem.14, 19. 10.1016/j.cogsc.2018.04.012

  • 42

    XuJ.-Z.RuanH.-Z.LiuL.-M.WangL.-P.ZhangW.-G. (2019). Overexpression of Thermostable Meso-Diaminopimelate Dehydrogenase to Redirect Diaminopimelate Pathway for Increasing L-Lysine Production in Escherichia coli. Sci. Rep.9 (1), 2423. 10.1038/s41598-018-37974-w

  • 43

    YadavA. N.KourD.RanaK. L.YadavN.SinghB.ChauhanV. S.et al (2019). “Metabolic Engineering to Synthetic Biology of Secondary Metabolites Production,” in New and Future Developments in Microbial Biotechnology and Bioengineering. Editors GuptaVKPandeyA, 279320. 10.1016/b978-0-444-63504-4.00020-7

  • 44

    YadavB.TalanA.TyagiR. D.DroguiP. (2021). Concomitant Production of Value-Added Products with Polyhydroxyalkanoate (PHA) Synthesis: A Review. Bioresour. Technol.337, 125419. 10.1016/j.biortech.2021.125419

  • 45

    YangL.CluettW. R.MahadevanEMILiOR. (2011). EMILiO: A Fast Algorithm for Genome-Scale Strain Design. Metab. Eng.13 (3), 272281. 10.1016/j.ymben.2011.03.002

  • 46

    ZahoorA.KüttnerF. T. F.BlankL. M.EbertB. E. (2019). Evaluation of Pyruvate decarboxylase‐negativeSaccharomyces Cerevisiaestrains for the Production of Succinic Acid. Eng. Life Sci.19 (10), 711720. 10.1002/elsc.201900080

  • 47

    ZhangX.JantamaK.ShanmugamK. T.IngramL. O. (2009). Reengineering Escherichia coli for Succinate Production in Mineral Salts Medium. Appl. Environ. Microbiol.75 (24), 78077813. 10.1128/aem.01758-09

  • 48

    ZhangY.LiJ.MengJ.SunK.YanH. (2021). A Neutral Red Mediated Electro-Fermentation System of Clostridium Beijerinckii for Effective Co-production of Butanol and Hydrogen. Bioresour. Technol.332, 125097. 10.1016/j.biortech.2021.125097

Summary

Keywords

metabolic modelling, genome-scale models, bioproduction, concomitant production, co-synthesis, constraint-based modelling

Citation

Raajaraam L and Raman K (2022) A Computational Framework to Identify Metabolic Engineering Strategies for the Co-Production of Metabolites. Front. Bioeng. Biotechnol. 9:779405. doi: 10.3389/fbioe.2021.779405

Received

18 September 2021

Accepted

02 December 2021

Published

07 January 2022

Volume

9 - 2021

Edited by

Hua Ling, National University of Singapore, Singapore

Reviewed by

Meiyappan Lakshmanan, Bioprocessing Technology Institute (A∗STAR), Singapore

Guoqiang Xu, Jiangnan University, China

Updates

Copyright

*Correspondence: Karthik Raman,

This article was submitted to Synthetic Biology, a section of the journal Frontiers in Bioengineering and Biotechnology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics