"
Team:Groningen/Notebook/Modeling 6June2012
From 2012.igem.org
m |
m |
||
| Line 13: | Line 13: | ||
1. Before, there were only a couple subsystems. The exchange 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. | 1. Before, there were only a couple subsystems. The exchange 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. | 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 carbon. | 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 carbon. | ||
Revision as of 22:44, 6 June 2012
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 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 carbon.
