Team:Manchester/Modelling





Flux Balance Analysis


  • We define Flux Balance Analysis and its importance and value added to our project.
  • We performed a first analysis to determine the importance of tyrosine in the media.
  • We analysed the distribution of growth at a fixed substrate production flux, which can be used to define optimal growth.
  • We predicted the involvement of glycolysis in growth of our engineered strain and related it to the addition of tyrosine (a product precursor) in the medium.

What is Flux Balance Analysis?

The main aim of our project is to overproduce Hipposudoric Acid (HA) in recombinant Escherichia coli. For that we need to achieve optimal production of the HA’s precursor homogentisic acid (HGA). To predict HGA production and explore the metabolic constraints affecting flux towards our product we used Flux Balance Analysis (FBA).

FBA is a method for computing steady-state metabolic flux distributions via genome-scale metabolic models, which are built using known intracellular metabolic reactions and gene-protein-reaction associations (1).

The recently published E. coli genome-scale reconstruction iJO1366 was used in all simulations. (2, 3, 4) FBA was performed using the COBRA Toolbox v.3.0. in MATLAB.

First step: introducing the HGA synthesis pathway

As wild-type E. coli does not produce HGA natively, the first step with FBA is adding the HGA synthesis pathway (simplified to a one-step enzymatic reaction) to the wild-type model:

FBA Reaction

Where 34hpp 3-(4-hydroxyphenyl)pyruvate, a metabolite natively present in E. coli, and HGA is homogentisic acid.

Media Comparison

We first examined and compared the effect of different nutrients (i.e. glucose minimal media only and glucose minimal media with tyrosine supplemented) on the growth of wild type E. coli and HGA-producing E. coli. HGA’s precursor is a tyrosine derivative, naturally we would expect an increased growth of the HGA-producing E. coli if tyrosine was added to the growth medium. To test this assumption we included a tyrosine uptake flux to the model. The tyrosine uptake rate was set to various levels, i.e. zero (glucose only), 20% and 100% of the glucose uptake rate, to investigate the impact on the growth of HGA-production strain. (Figure 1). FBA predicts that the tyrosine supplementation improves the growth of HGA-producing E. coli more significantly than that of the wild-type strain, which validates our assumption.


Figure 1a
Figure 1b
Figure 1c
Figure 1. Growth with fluxes at variable levels of Tyrosine in the media. Minimal media (MM), Tyrosine (tyr), Glucose (glc), Glucose uptake rate (GUR) (a) FBA prediction of maximum growth rate of wild-type and hga-producing E. coli grown on glucose minimal media with no tyrosine supplemented. (b) FBA prediction of maximum growth rate of wild-type and HGA-producing E. coli grown on glucose minimal media with medium-level tyrosine supplemented, the tyrosine uptake rate was set to 20% of the glucose uptake rate. (c) Prediction of maximum growth rate of wild-type and HGA-producing E. coli grown on glucose minimal media with high-level tyrosine supplemented, the tyrosine uptake rate was set to 100% of glucose uptake rate.

E. coli growth at fixed HGA production rate

The second test was designed to investigate the change in maximum growth rate at different glucose and oxygen uptake rate given a fixed HGA synthesis rate.This analysis predicts the optimal growth for every uptake value of glucose and oxygen, from 0 mmol/gdw/h to 20 mmol/gdw/h), at a fixed HGA production rate of 1 mmol/gdw/h (Figure 2).

Figure 2

Figure 2. FBA prediction of the maximum growth rate at various glucose uptake rate (GUR) and oxygen uptake rate (OUR). HGA production rate was set to 1 mmol/gdw/h for all simulations.

Figure 2 can be useful for designing optimal growth conditions. For example, given a fixed GUR of 10 mmol/gdw/h, the optimal OUR that gives the highest growth rate might be around 16 mmol/gdw/h, or alternatively, given a fixed OUR we can infer that the optimal GUR might be 20 mmol/gdw/h or higher (Figure 3).


FBA 5
FBA 6
Figure 3. Projections of the growth distribution from Figure 2. In orange, growth rate against GUR for Constant OUR at 15 mmol/gdw/h. In blue, growth rate against growth rate against OUR for Constant GUR at 10 mmol/gdw/h.

Effect of tyrosine on Glycolysis and growth

The final test performed is the analysis of the metabolic flux carried by enolase (ENO), a key enzyme of glycolysis that catalyses the conversion of 2-phosphoglycerate (2pg) to phosphoenolpyruvate (pep) (5). We choose enolase not only because of its central role in the lower part of glycolysis, but also because the product of ENO, i.e. pep, is an important substrate of the Shikimate pathway, from which HGA derives.


Figure 4a
Figure 4b
Figure 4. Effect of tyrosine addition on glycolysis and growth. (a) Comparison of ENO activity and its effect on growth on the recombinant strain in a glucose minimal medium (blue) , and in a tyrosine enriched minimal medium(orang) at fixed tyrosine uptake rate of 1 mmol/gdw/h.(b) Comparison of growth rate of the recombinant strain in a glucose minimal medium (blue) , and in a tyrosine enriched minimal medium (orange) at fixed tyrosine uptake rate of 1 mmol/gdw/h.

From Figure 4 we can see that the maximum growth rate of the HGA-producing E. coli grown on a glc-tyr mixed medium is higher that of the same strain grown on a glc only minimal medium. Furthermore, as the externally supplemented tyrosine reduces the need of de novo tyrosine synthesis (which requires pep as a precursor), the overall demand on the pep pool declines compared to the case without tyrosine addition. Therefore, at the same growth rate and HGA production rate, the ENO flux of the strain grown on a glc-tyr MM medium is predicted to be lower than that of the strain grown on glc MM.

Conclusion

The three tests performed gave us insight into the possible optimal conditions to produce HGA. The main takeaway is that the addition of tyrosine to the culture medium improves both cell growth and HGA production, which can be beneficial for obtaining a higher yield of hipposudoric acid in a bioreactor.

References

Literature

Orth J, Thiele I, and Palsson, B. (2010) What is flux balance analysis? Nat Biotechnol 28, 245–248
Heirendt L, Arreckx S, Pfau T, et al. (2019) Creation and analysis of biochemical constraint-based models using the COBRA Toolbox v.3.0. Nat Protoc 14, 639–702
Vlassis N, Pacheco MP, Sauter T. (2014) Fast Reconstruction of Compact Context-Specific Metabolic Network Models. PLoS Comput Biol 10: e1003424
Orth, Jeffrey D, Conrad, Tom M, Na, Jessica, Lerman, Joshua A, Nam, Hojung, Feist, Adam M, Palsson, Bernhard Ø, (2011)Orth JD, Conrad TM, Na J, Lerman JA, Nam H, Feist AM, and Palsson BØ, (2011) A comprehensive genome-scale reconstruction of Escherichia coli metabolism. Molecular Systems Biology 7, 535
Zhang E, Brewer JM, Minor W, Carreira LA, Lebioda L. (1997) Mechanism of enolase: The crystal structure of asymmetric dimer enolase−2-phospho-d-glycerate/enolase−phosphoenolpyruvate at 2.0 Å resolution. Biochemistry, 6,12526-34.
Logo 3 white


Logo 4


Logo 5
Logo 6 Logo 7


igem2020manchester@gmail.com


Logo 8 png Logo 9 png
Logo 1


Logo 2