Team:Groningen/Notebook/Modeling 6June2012


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.

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