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:
- It is a type of PHA where Dr. Rao's expertise and the information he provided during our collaboration could be leveraged.
- 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.
- 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:
- Which potential gene/protein could be engineered into our devices that will increase the amount of NADPH available to the PHB biosynthesis pathway?
- 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

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_5155Associated 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 |
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 |
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 |
2. Increase serA
serA BIGG ID: PP_5155Associated 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 |
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 |
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 |
3. Increase zwf
zwf BIGG ID: PP_5351Associated 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 |
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 |
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 |
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 |
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 |
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 |
5. Increase sdaA
Genes: tdcG-I, tdcG-II, tdcG-IIIBIGG 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 |
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 |
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 |
6. Increase fbaA
fbaA BIGG ID: PP_4960Associated 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 |
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 |
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 |
7. Knock out pgi
Pgi BIGG IDs: PP_4701, PP_1808Weight 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 |
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 |
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 |
8. Coexpression of sdaA, serA, and pgk
Pgk BIGG ID: PP_4963Associated 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 |
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 |
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 |
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

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 |
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 |
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 |
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)


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.