Engineering

Learn how our team incorporated the design engineering cycle!

Introduction

Our team demonstrated engineering success through an iterative approach to metabolic engineering, resulting in the design of composite synthetic biology devices essential for achieving our goal of upcycling PET plastic. This undertaking was informed by our Integrated Human Practices work, where we consulted with Dr. Chitong Rao, Chief Scientist at Bluepha Co., Ltd.—a global leader in using synthetic biology for degradable bio-polymer production—and Dr. Joanne Sadler from the University of Edinburgh, whose work on producing vanillin from the microbial breakdown of PET served as initial inspiration for our project.

From these meetings, we learned that metabolic engineering to increase yields of useful products would be essential for making our concepts both economically viable and environmentally sustainable. This insight was further confirmed through our team's Life Cycle Analysis of our desired products and during the development of our business plan.

Our goal was to increase the production of poly(3-hydroxybutyrate) (PHB), a type of polyhydroxyalkanoate (PHA) that can be downstream depolymerized into beta-hydroxybutyrate (BHB). We chose this target for several reasons:

  1. It is a type of PHA where Dr. Rao's expertise and the information he provided during our collaboration could be leveraged.
  2. PHB can be depolymerized by an enzyme known to the team, as described by the 2024 Concordia iGEM team, thus enabling for the production of BHB as a high-value product from increased PHB levels.
  3. We had already designed genetic circuits in both E. coli and P. putida to produce PHB, allowing our research to focus on enhancing this biosynthesis pathway familiar to the team.

Through literature research, we discovered that increasing the production of NADPH and Acetyl-CoA would result in increased PHB production via the PHB biosynthesis pathway (Shi et al., 1999). We utilized Flux Balance Analysis (FBA), a constraint-based metabolic modeling method, using several target genes, enzymes, and reactions identified from the literature (Lim et a., 2002; Zhang et al., 2020, Zhang et al., 2014).

The first iteration of the Design-Build-Test-Learn (DBTL) cycle tested whether incorporating certain genetic targets (identified from literature) would have the desired effect of increasing NADPH production by increasing metabolic flux. After verifying plausible targets, the second iteration determined whether our designed devices would elicit an increase in NADPH when modeled with expected enzyme concentrations and reaction fluxes.

Our engineering design cycle addressed the following questions:

  1. Which potential gene/protein could be engineered into our devices that will increase the amount of NADPH available to the PHB biosynthesis pathway?
  2. Will the engineering of new genes/proteins into these devices increase the amount of NADPH available to the PHB biosynthesis pathway when the solved enzyme concentrations and reaction rates are incorporated?

Iteration 1

Figure 1: DBTL cycle for Iteration 1.

Design

Suggested by the literature on engineering E.coli to improve poly(3-hydroxybutyrate) production, the team identified genes, enzymes, and reaction targets of which we hypothesized would be effective in increasing the production of Acetyl-CoA and NADPH. Subsequently, we hypothesized an increase in the production of PHB, a high-value product. (Zhang et al., 2014) In the metabolic engineering upregulation process, the reaction targets were G6PD and PDH, as they were proven to be responsible for the production of Acetyl-CoA and NADPH. (Lim et a., 2002; Zhang et al., 2020) The identified gene and enzyme targets included increasing zwf, increasing serA, and more.

Build

The team identified two models: iJN1463 for Pseudomonas putida KT2440 and iEC1356_Bl21DE3 for E. coli BL21(DE3). Using the model of Pseudomonas putida KT2440, we utilized Flux Balance Analysis, a constraints-based modeling, to test the effect of targets. To increase a certain gene or enzyme, the upper and lower bounds of the associated reaction need to be increased. Therefore, we identified the reactions that the selected target genes code for. The upper and lower bounds during the first iteration were theoretical values, as it was impossible to know the exact possible increase in the selected targets.

Test

We implemented a Python package, COBRApy, to conduct a flux balance analysis with the gene, enzyme, and reaction targets.

All influx and outflux values are expressed in mmol/gDW/h.

1. Increase serA, serB, and serC

serA BIGG ID: PP_5155
Associated reaction: Phosphoglycerate dehydrogenase (PGCD)

serB BIGG ID: PP_4909
Associated reaction: Phosphoserine phosphatase (L-serine) (PSP_L) Associated reaction: Phosphoglycerate dehydrogenase (PGCD)

serC BIGG ID: PP_1768
Associated reactions: O-Phospho-4-hydroxy-L-threonine (OHPBAT)

Weight 0.5 (G6PD) : 0.5 (PDH) Before ser increase After ser increase
NADPH influx 209.080 183.387
NADPH outflux 197.080 524.320
Acetyl-CoA influx 3.898e-14 5.826e-14
Acetyl-CoA outflux 298.160 144
Table 1: Increase serA, serB, and serC (0.5:0.5 weighing)

Weight 0.7 (G6PD) : 0.3 (PDH) Before ser increase After ser increase
NADPH influx 75.387 75.387
NADPH outflux 596.320 596.320
Acetyl-CoA influx 0 0
Acetyl-CoA outflux 1.312e-14 0
Table 2: Increase serA, serB, and serC (0:7:0.3 weighing)

Weight 0.3 (G6PD) : 0.7 (PDH) Before ser increase After ser increase
NADPH influx 197.080 197.080
NADPH outflux 233.080 66.920
Acetyl-CoA influx 3.111e-14 6.539e-16
Acetyl-CoA outflux 310.160 310.160
Table 3: Increase serA, serB, and serC (0.3:0.7 weighing)


2. Increase serA

serA BIGG ID: PP_5155
Associated reaction: Phosphoglycerate dehydrogenase (PGCD)

Weight 0.5 (G6PD) : 0.5 (PDH) Before ser increase After ser increase
NADPH influx 209.080 183.387
NADPH outflux 197.080 524.320
Acetyl-CoA influx 3.898e-14 5.826e-14
Acetyl-CoA outflux 298.160 144
Table 4: Increase serA (0.5:0.5 weighing)

Weight 0.7 (G6PD) : 0.3 (PDH) Before ser increase After ser increase
NADPH influx 75.387 75.387
NADPH outflux 596.320 596.320
Acetyl-CoA influx 0 0
Acetyl-CoA outflux 1.312e-14 2.010e-14
Table 5: Increase serA (0.7:0.3 weighing)

Weight 0.3 (G6PD) : 0.7 (PDH) Before ser increase After ser increase
NADPH influx 197.080 197.080
NADPH outflux 233.080 66.920
Acetyl-CoA influx 3.111e-14 0
Acetyl-CoA outflux 310.160 310.160
Table 6: Increase serA (0.3:0.7 weighing)


3. Increase zwf

zwf BIGG ID: PP_5351
Associated reaction: Beta-D-Glucose-6-phosphate NADP+ 1-oxoreductase (G6PBDH)

Weight 0.5 (G6PD) : 0.5 (PDH) Before zwf increase (G6PBDH) After zwf increase
NADPH influx 209.080 10209.080
NADPH outflux 197.080 10197.080
Acetyl-CoA influx 3.898e-14 3.898e-14
Acetyl-CoA outflux 298.160 298.160
Table 7: Increase zwf (0.5:0.5 weighing)

Weight 0.7 (G6PD) : 0.3 (PDH) Before zwf increase (G6PBDH) After zwf increase
NADPH influx 75.387 10075.387
NADPH outflux 596.320 10596.320
Acetyl-CoA influx 0 0
Acetyl-CoA outflux 1.312e-14 0
Table 8: Increase zwf (0.7:0.3 weighing)

Weight 0.3 (G6PD) : 0.7 (PDH) Before zwf increase (G6PBDH) After zwf increase
NADPH influx 197.080 10197.080
NADPH outflux 233.080 10233.080
Acetyl-CoA influx 3.111e-14 3.181e-12
Acetyl-CoA outflux 310.160 310.160
Table 9: Increase zwf (0.3:0.7 weighing)


4. Increase zwf and serA


Weight 0.5 (G6PD) : 0.5 (PDH) Before zwf and serA increase After zwf and serA increase
NADPH influx 209.080 10183.387
NADPH outflux 197.080 10524.32
Acetyl-CoA influx 3.898e-14 9.702e-14
Acetyl-CoA outflux 298.160 144
Table 10: Increase zwf and serA (0.5:0.5 weighing)

Weight 0.7 (G6PD) : 0.3 (PDH) Before zwf and serA increase After zwf and serA increase
NADPH influx 75.387 10075.387
NADPH outflux 596.320 10596.320
Acetyl-CoA influx 0 0
Acetyl-CoA outflux 1.312e-14 0
Table 11: Increase zwf and serA (0.7:0.3 weighing)

Weight 0.3 (G6PD) : 0.7 (PDH) Before zwf and serA increase After zwf and serA increase
NADPH influx 197.080 10197.080
NADPH outflux 233.080 10066.920
Acetyl-CoA influx 3.111e-14 0
Acetyl-CoA outflux 310.160 310.160
Table 12: Increase zwf and serA (0.3:0.7 weighing)


5. Increase sdaA

Genes: tdcG-I, tdcG-II, tdcG-III
BIGG ID: PP_0297, PP_0987, PP_3144
Associated reaction: L-serine deaminase (SERD_L)

Weight 0.5 (G6PD) : 0.5 (PDH) Before sdaA increase After sdaA increase
NADPH influx 209.080 183.387
NADPH outflux 197.080 524.320
Acetyl-CoA influx 3.898e-14 0
Acetyl-CoA outflux 298.160 144
Table 13: Increase sdaA (0.5:0.5 weighing)

Weight 0.7 (G6PD) : 0.3 (PDH) Before sdaA increase After sdaA increase
NADPH influx 75.387 75.387
NADPH outflux 596.320 596.320
Acetyl-CoA influx 0 4.704e-15
Acetyl-CoA outflux 1.312e-14 7.116e-13
Table 14: Increase sdaA (0.7:0.3 weighing)

Weight 0.3 (G6PD) : 0.7 (PDH) Before sdaA increase After sdaA increase
NADPH influx 197.080 131.080
NADPH outflux 233.080 167.080
Acetyl-CoA influx 3.111e-14 7.474e-15
Acetyl-CoA outflux 310.160 310.160
Table 15: Increase sdaA (0.3:0.7 weighing)


6. Increase fbaA

fbaA BIGG ID: PP_4960
Associated reaction: Fructose-bisphosphate aldolase (FBA)

Weight 0.5 (G6PD) : 0.5 (PDH) Before fbaA increase After fbaA increase
NADPH influx 209.080 209.080
NADPH outflux 197.080 197.080
Acetyl-CoA influx 3.898e-14 3.898e-14
Acetyl-CoA outflux 298.160 298.160
Table 16: Increase fbaA (0.5:0.5 weighing)

Weight 0.7 (G6PD) : 0.3 (PDH) Before fbaA increase After fbaA increase
NADPH influx 75.387 75.387
NADPH outflux 596.320 596.320
Acetyl-CoA influx 0 0
Acetyl-CoA outflux 1.312e-14 1.312e-14
Table 17: Increase fbaA (0.7:0.3 weighing)

Weight 0.3 (G6PD) : 0.7 (PDH) Before fbaA increase After fbaA increase
NADPH influx 197.080 197.080
NADPH outflux 233.080 233.080
Acetyl-CoA influx 3.111e-14 3.111e-14
Acetyl-CoA outflux 310.160 310.160
Table 18: Increase fbaA (0.3:0.7 weighing)


7. Knock out pgi

Pgi BIGG IDs: PP_4701, PP_1808

Weight 0.5 (G6PD) : 0.5 (PDH) Before knockout After knockout
NADPH influx 209.080 197.080
NADPH outflux 197.080 185.080
Acetyl-CoA influx 3.898e-14 0
Acetyl-CoA outflux 298.160 310.160
Table 19: Knock out pgi (0.5:0.5 weighing)

Weight 0.7 (G6PD) : 0.3 (PDH) Before knockout After knockout
NADPH influx 75.387 197.080
NADPH outflux 596.320 185.080
Acetyl-CoA influx 0 5.227e-47
Acetyl-CoA outflux 1.312e-14 310.160
Table 20: Knock out pgi (0.7:0.3 weighing)

Weight 0.3 (G6PD) : 0.7 (PDH) Before knockout After knockout
NADPH influx 197.080 197.080
NADPH outflux 233.080 185.080
Acetyl-CoA influx 3.111e-14 0
Acetyl-CoA outflux 310.160 310.160
Table 21: Knock out pgi (0.3:0.7 weighing)


8. Coexpression of sdaA, serA, and pgk

Pgk BIGG ID: PP_4963
Associated reaction: Phosphoglycerate kinase (PP_4963)

Weight 0.5 (G6PD) : 0.5 (PDH) Before coexpression After coexpression
NADPH influx 209.080 183.387
NADPH outflux 197.080 524.320
Acetyl-CoA influx 3.898e-14 0
Acetyl-CoA outflux 298.160 144
Table 22: Coexpression of sdaA, serA, and pgk (0.5:0.5 weighing)

Weight 0.7 (G6PD) : 0.3 (PDH) Before coexpression After coexpression
NADPH influx 75.387 75.387
NADPH outflux 596.320 596.320
Acetyl-CoA influx 0 4.704e-15
Acetyl-CoA outflux 1.312e-14 7.116e-13
Table 23: Coexpression of sdaA, serA, and pgk (0.7:0.3 weighing)

Weight 0.3 (G6PD) : 0.7 (PDH) Before coexpression After coexpression
NADPH influx 197.080 131.080
NADPH outflux 233.080 167.080
Acetyl-CoA influx 3.111e-14 7.474e-15
Acetyl-CoA outflux 310.160 310.160
Table 24: Coexpression of sdaA, serA, and pgk (0.3:0.7 weighing)


Learn

serA and zwf were shown to be successful in increasing the NADPH yield, though they had minimal effect on Acetyl-CoA yield. However, increasing NADPH yield can have the desired effect of downstream upregulation of key metabolites in PHB synthesis. Other targets showed insignificant changes in yields. Hence, serA and zwf were chosen as effective targets, specifically zwf being the most successful one.


Iteration 2

Figure 2: DBTL cycle for Iteration 2.

Design

From Iteration 1, the team designed genetic circuits using targets and standard parts, including inducible promoter systems T7 (E.coli) and Xysl (P.putida). The genetic circuits included serA and zwf as they were shown to be effective.

Build

To find yields of key enzymes, we aimed to find exact upper reaction bounds for our Constraint-Based Modeling. This was done through ODEs and enzyme kinetics for our designed circuits. We modelled PHGDH production by serA to obtain the maximum enzyme concentration. A similar work was done with the second target, zwf. Then, utilizing enzyme kinetics for our designed circuits, that value was converted to an upper bound that could be used in the FBA. Full working can be found here: Model

Test

To test how our system will behave after modifying the model, we used the upper bounds for our model from previous steps in the engineering cycle and conducted a second flux balance analysis for zwf.

Weight 0.5 (G6PD) : 0.5 (PDH) Before zwf increase (G6PBDH) After zwf increase
NADPH influx 209.080 955.940
NADPH outflux 197.080 983.949
Acetyl-CoA influx 3.898e-14 3.287e-15
Acetyl-CoA outflux 298.160 298.160
Table 25: Increase zwf (0.5:0.5 weighing)

Weight 0.7 (G6PD) : 0.3 (PDH) Before zwf increase (G6PBDH) After zwf increase
NADPH influx 75.387 862.247
NADPH outflux 596.320 1383.180
Acetyl-CoA influx 0 0
Acetyl-CoA outflux 1.312e-14 9.714e-14
Table 26: Increase zwf (0.7:0.3 weighing)

Weight 0.3 (G6PD) : 0.7 (PDH) Before zwf increase (G6PBDH) After zwf increase
NADPH influx 197.080 983.940
NADPH outflux 233.080 1019.940
Acetyl-CoA influx 3.111e-14 3.111e-14
Acetyl-CoA outflux 310.160 310.160
Table 27: Increase zwf (0.3:0.7 weighing)

Learn

The team learned that zwf increase is effective in increasing NADPH yield in our designed inducible systems. We confirmed our decision to utilize zwf in the experimentations, since it was also most effective in increasing NADPH and following PHB yield, a high-value product produced from PET plastic. We then proceeded to the final design based on data from FBA models.


Final Product (Designed Composite Parts)

Figure 3: ZWF gene (pink) which codes for Glucose-6-Phosphate-Dehydrogenase under the influence of T7 promoter. PhaC/A/B code for the biosynthesis of PHB. Bba_K2260002 is Phasin.

Figure 4: Zwf in P. putida under the Xysl/Pm expression system. SerA gene (pink) which codes for D-3-phosphoglycerate dehydrogenase. PhaC/A/B code for the biosyntehsis of PHB.

Through constraint-based metabolic modeling, our team demonstrated engineering success by showing that the increase in zwf and serA increased the amount of NADPH available to the PHB biosynthesis pathway. Moreover, our designed devices elicited an increase in NADPH when modeled with expected enzyme concentrations and reaction fluxes, thus confirming that, theoretically, our gene engineering result in increased PHB production.

References

1 5.2: Enzyme Parameters. (2019, July 15). Chemistry LibreTexts.
2 7.2: Derivation of Michaelis-Menten equation. (2021, August 23). Biology LibreTexts.
3 Batianis, C., Van Rosmalen, R. P., Major, M., Van Ee, C., Kasiotakis, A., Weusthuis, R. A., & Martins Dos Santos, V. A. P. (2023). A tunable metabolic valve for precise growth control and increased product formation in Pseudomonas putida. Metabolic Engineering, 75, 47–57.
4 Bernstein, J. A., Lin, P.-H., Cohen, S. N., & Lin-Chao, S. (2004). Global analysis of Escherichia coli RNA degradosome function using DNA microarrays. Proceedings of the National Academy of Sciences, 101(9), 2758–2763.
5 Bintu, L., Buchler, N. E., Garcia, H. G., Gerland, U., Hwa, T., Kondev, J., & Phillips, R. (2005). Transcriptional regulation by the numbers: models. Current opinion in genetics & development, 15(2), 116–124.
6 Chang, A., Jeske, L., Ulbrich, S., Hofmann, J., Koblitz, J., Schomburg, I., Neumann-Schaal, M., Jahn, D., & Schomburg, D. (2021). BRENDA, the ELIXIR core data resource in 2021: New developments and updates. Nucleic Acids Research, 49(D1), D498–D508.
7 Daun, S., Rubin, J., Vodovotz, Y., & Clermont, G. (2008). Equation-based models of dynamic biological systems. Journal of critical care, 23(4), 585–594.
8 Ebrahim, A., Lerman, J. A., Palsson, B. O., & Hyduke, D. R. (2013). COBRApy: COnstraints-Based Reconstruction and Analysis for Python. BMC Systems Biology, 7(1), 74.
9 Gibson, B., Wilson, D. J., Feil, E., & Eyre-Walker, A. (2018). The distribution of bacterial doubling times in the wild. Proceedings. Biological sciences, 285(1880), 20180789.
10 Giordano, T. J., Deuschle, U., Bujard, H., & McAllister, W. T. (1989). Regulation of coliphage T3 and T7 RNA polymerases by the lac represser-operator system. Gene, 84(2), 209–219.
11 Hussey, B. J., & McMillen, D. R. (2018). Programmable T7-based synthetic transcription factors. Nucleic acids research, 46(18), 9842–9854.
12 Jeong, H., Kim, H. J., & Lee, S. J. (2015). Complete Genome Sequence of Escherichia coli Strain BL21. Genome announcements, 3(2), e00134-15.
13 King, Z. A., Lu, J., Dräger, A., Miller, P., Federowicz, S., Lerman, J. A., Ebrahim, A., Palsson, B. O., & Lewis, N. E. (2016). BiGG Models: A platform for integrating, standardizing and sharing genome-scale models. Nucleic Acids Research, 44(D1), D515–D522.
14 Kuzmine, I., Gottlieb, P. A., & Martin, C. T. (2003). Binding of the priming nucleotide in the initiation of transcription by T7 RNA polymerase. The Journal of biological chemistry, 278(5), 2819–2823.
15 Lim, S.-J., Jung, Y.-M., Shin, H.-D., & Lee, Y.-H. (2002). Amplification of the NADPH-related genes zwf and gnd for the oddball biosynthesis of PHB in an E. coli transformant harboring a cloned phbCAB operon. Journal of Bioscience and Bioengineering, 93(6), 543–549.
16 Lutz, R., & Bujard, H. (1997). Independent and tight regulation of transcriptional units in Escherichia coli via the LacR/O, the TetR/O and AraC/I1-I2 regulatory elements. Nucleic Acids Research, 25(6), 1203–1210.
17 Maurizi M. R. (1992). Proteases and protein degradation in Escherichia coli. Experientia, 48(2), 178–201.
18 McClure W. R. (1980). Rate-limiting steps in RNA chain initiation. Proceedings of the National Academy of Sciences of the United States of America, 77(10), 5634–5638.
19 Michaelis-Menten Equation—An overview | ScienceDirect Topics. (n.d.).
20 Michaelis-Menten Kinetics. (2013, October 2). Chemistry LibreTexts.
21 Miyanaga, K., & Unno, H. (2011). Reaction Kinetics and Stoichiometry. In Comprehensive Biotechnology (pp. 33–46). Elsevier.
22 Molpeceres-García, F. J., Sanz-Mata, D., García-Miro, A., Prieto, A., & Barriuso, J. (2025). Towards polyethylene terephthalate valorisation into PHB using an engineered Comamonas testosteroni strain. New biotechnology, 85, 75–83.
23 Nogales, J., Mueller, J., Gudmundsson, S., Canalejo, F. J., Duque, E., Monk, J., Feist, A. M., Ramos, J. L., Niu, W., & Palsson, B. O. (2020). High-quality genome-scale metabolic modelling of Pseudomonas putida highlights its broad metabolic capabilities. Environmental Microbiology, 22(1), 255–269.
24 Olavarria, K., Marone, M. P., da Costa Oliveira, H., Roncallo, J. C., da Costa Vasconcelos, F. N., da Silva, L. F., & Gomez, J. G. C. (2015). Quantifying NAD(P)H production in the upper Entner–Doudoroff pathway from Pseudomonas putida KT2440. FEBS Open Bio, 5(1), 908–915.
25 Olsson, Sanna & Penacho, Vanessa & Puente-Sánchez, Fernando & Díaz, Silvia & Gonzalez-Pastor, Eduardo & Aguilera, Angeles. (2017). Horizontal Gene Transfer of Phytochelatin Synthases from Bacteria to Extremophilic Green Algae. Microbial Ecology, 73.
26 Prinz H. (2010). Hill coefficients, dose-response curves and allosteric mechanisms. Journal of chemical biology, 3(1), 37–44.
27 Roberts, E., Magis, A., Ortiz, J. O., Baumeister, W., & Luthey-Schulten, Z. (2011). Noise contributions in an inducible genetic switch: A whole-cell simulation study. PLoS Computational Biology, 7(3).
28 Rowe, E., Palsson, B. O., & King, Z. A. (2018). Escher-FBA: A web application for interactive flux balance analysis. BMC Systems Biology, 12(1), 84.
29 Ruiz, N., & Silhavy, T. J. (2022). How Escherichia coli Became the Flagship Bacterium of Molecular Biology. Journal of bacteriology, 204(9), e0023022.
30 Shah, B. A., Kasarlawar, S. T., & Phale, P. S. (2022). Glucose-6-Phosphate Dehydrogenase, ZwfA, a Dual Cofactor-Specific Isozyme Is Predominantly Involved in the Glucose Metabolism of Pseudomonas bharatica CSV86T. Microbiology Spectrum, 10(6), e03818-22.
31 Shi, H., Nikawa, J., & Shimizu, K. (1999). Effect of modifying metabolic network on poly-3-hydroxybutyrate biosynthesis in recombinant Escherichia coli. Journal of Bioscience and Bioengineering, 87(5), 666–677.
32 Sites.duke.edu/bossbackup/files/2013/02/nonlime. (n.d.).
33 Srinivasan, B. (2022). A guide to the Michaelis–Menten equation: Steady state and beyond. The FEBS Journal, 289(20), 6086–6098.
34 Tabor S. (2001). Expression using the T7 RNA polymerase/promoter system. Current protocols in molecular biology, Chapter 16, Unit16.2.
35 Tokic, M., Hatzimanikatis, V., & Miskovic, L. (2020). Large-scale kinetic metabolic models of Pseudomonas putida KT2440 for consistent design of metabolic engineering strategies. Biotechnology for Biofuels, 13(1), 33.
36 Verhoef, S., Ballerstedt, H., Volkers, R. J. M., de Winde, J. H., & Ruijssenaars, H. J. (2010). Comparative transcriptomics and proteomics of p-hydroxybenzoate producing Pseudomonas putida S12: Novel responses and implications for strain improvement. Applied Microbiology and Biotechnology, 87(2), 679–690.
37 Zhang, Y., Lin, Z., Liu, Q., Li, Y., Wang, Z., Ma, H., Chen, T., & Zhao, X. (2014). Engineering of Serine-Deamination pathway, Entner-Doudoroff pathway and pyruvate dehydrogenase complex to improve poly(3-hydroxybutyrate) production in Escherichia coli. Microbial Cell Factories, 13(1), 172
38 Zhang, Y., Su, M., Qin, N., Nielsen, J., & Liu, Z. (2020). Expressing a cytosolic pyruvate dehydrogenase complex to increase free fatty acid production in Saccharomyces cerevisiae. Microbial Cell Factories, 19(1), 226.
39 Ziegler, A. L., Filbig, M., Parschau, J., Tiso, T., Blank, L. M., & Mitsos, A. (2023). Combined metabolic modeling and experimental data for enhanced biotechnological production. In A. C. Kokossis, M. C. Georgiadis, & E. Pistikopoulos (Eds.), 33rd european symposium on computer aided process engineering (Vol. 52, pp. 2685–2690). Elsevier.