Team:York/Modeling

Modeling

(Click Images to Enlarge)

To guide the design of our final organism the team has employed the use of Flux Balance Analysis (FBA) modeling. FBA calculates estimates of the rates at which metabolites in an organism will pass through various reaction pathways. As phosphate is used in nearly every part of the cells metabolome, we felt it was important to determine how the accumulation of large amounts might affect the overall function of our organism. Using the iAF1260 model[1] of E. coli as a basis, we first attempted to look at the reactions controlled by the genes we had found as candidates for manipulation during the team's dry lab period. These include, Polyphosphate Kinase (PPK), Exopolyphosphatase (PPX), and the Phosphate Specific Transporters (Pst). One of the immediate problems we found is that phosphate is used by many different parts of the cell, and as such altering any one part has knock on effects on many others. Attempting to alter phosphate pathways requires a fine hand.

Figure 1. Metabolic pathway map for E. Coli generated from KEGG[2]. Each of the red nodes represents a metabolic reaction in which phosphate is a major component. View Original

The FBA computer model looks for the most efficient, if not the most realistic, pathways to use, meaning it was difficult to force any amount of flux through the pathways we had created to accumulate phosphate. Results from multiple different attempts to achieve theoretical uptake showed that to be able to accumulate phosphate, growth needed to be compromised. By forcing the model to compromise a fraction of it's growth potential we were able to show much higher storage rates. These rates increase in a linear way as further restrictions on growth are put in place, showing direct links between growth and our organisms ability to accumulate phosphate.

Figure 2. Graph showing correlation between theoretical growth restriction and accumulation of polyphosphate.

The results of our model allowed us to better direct the project towards the genes that would give us our desired result, and provide a good starting point for optimising the levels of gene expression. The outcomes of the model also point to potential further developments that are currently out of reach due to the iGEM timeframe, but would greatly develop our final product. Some of these ideas are explored in our business plan.

References

  • [1] Feist AM, Henry CS, Reed JL, et al. A genome-scale metabolic reconstruction for Escherichia coli K-12 MG1655 that accounts for 1260 ORFs and thermodynamic information. Molecular Systems Biology. 2007;3:121. doi:10.1038/msb4100155.
  • [2] Kanehisa M, Goto S (2000). "KEGG: Kyoto Encyclopedia of Genes and Genomes". Nucleic Acids Res 28 (1): 27–30. doi:10.1093/nar/28.1.27