Team:Groningen/Notebook/Modeling 6June2012

From 2012.igem.org

Revision as of 23:59, 17 September 2012 by Jparrish (Talk | contribs)






Converted SBML xml file to GDX for use in GAMS, then adapted a GAMS script to run flux balance analysis (FBA) with a CO2 objective. Unfortunately the result was a null solution. So, now the task is to reduce the involved reactions until there is a correct result, then keeping adding reactions until the solution changes. In this way I should be able to identify the problem reactions.
Interestingly enough, the SBML model did not include subsystems for the reactions. The subsystems had to be extracted from the xls file included in the article's supplementary information.
Christopher S Henry, Jenifer F Zinner, Matthew P Cohoon, and Rick L Stevens, "iBsu1103: a new genome-scale metabolic model of Bacillus subtilis based on SEED annotations," Genome Biology 2009, vol. 10(6) doi:10.1186/gb-2009-10-6-r69 http://genomebiology.com/2009/10/6/R69

Well, that was easier than expected. Two problems:
1. Before, there were only a couple subsystems. The exchange reactions as well as membrane transport reactions were classed as 'transport', so when I was restricting system input by setting the flux of all reactions in the 'Transport' subsystem to 0 I was actually preventing all membrane transport as well.
2. I was allowing influx of carbon other than glucose.
End result: After adding the subsystems, changing the script to restrict system input based on 'Exchange' rather than 'Transport', and restricting the carbon input to 10 glucose, the system correctly spit out 60 co2.

Back to notebook