Team:Groningen/Notebook/Modeling 6June2012

From 2012.igem.org

(Difference between revisions)
m
m
 
(5 intermediate revisions not shown)
Line 1: Line 1:
-
{{HeaderGroningen2012}}  
+
{{HeaderGroningen2012}}
 +
<html>
 +
<head>
 +
<style type="text/css">
 +
p.margin{
 +
font-size:12pt;
 +
line-height:12pt;
 +
color:white;
 +
margin-top:0px;
 +
margin-bottom:20px;
 +
margin-left:150px;
 +
margin-right:250px;
 +
}
 +
</style>
 +
</head>
-
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.
+
<body>
 +
<br>
 +
<br>
 +
<br>
 +
<p class="margin">
 +
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.
 +
<br>
 +
<br>
 +
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 <A HREF="http://genomebiology.com/2009/10/6/R69"><FONT COLOR=#ff6700>article's</FONT></A> supplementary information.
 +
<br>
 +
<br>
 +
Well, that was easier than expected. Two problems:
 +
<br>
 +
<br>
 +
<b>1.</b> 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.
 +
<br>
 +
<br>
 +
2. I was allowing influx of carbon other than glucose.
 +
<br>
 +
<br>
 +
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.
 +
<br>
 +
<A HREF="https://2012.igem.org/Team:Groningen/Notebook"><FONT COLOR=#ff6700>Back to notebook</FONT></A>
 +
</p>
 +
<br>
 +
<br>
 +
</body>
 +
</html>
-
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.
+
{{Template:SponsorsGroningen2012}}
-
 
+
-
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.
+

Latest revision as of 14:24, 26 September 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.

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