Team:NTNU Trondheim/Project/Modelling

From 2014.igem.org

(Difference between revisions)
 
(10 intermediate revisions not shown)
Line 208: Line 208:
<a href="https://2014.igem.org/Team:NTNU_Trondheim/Protocols">Protocols</a>
<a href="https://2014.igem.org/Team:NTNU_Trondheim/Protocols">Protocols</a>
</li>
</li>
-
<li class="pg-outreach pg-outreach_presentations pg-outreach_panels pg-outreach_awareness">
+
<li class="has-flyout pg-outreach pg-outreach_presentations pg-outreach_panels pg-outreach_awareness">
-
<a href="https://2014.igem.org/Team:NTNU_Trondheim/Sponsors">Acknowledgements</a>
+
<a href="https://2014.igem.org/Team:NTNU_Trondheim/Acknowledgements/Sponsors">Acknowledgements</a>
<a href="https://2014.igem.org/Team:NTNU_Trondheim/Project" class="flyout-toggle"><span> </span></a>
<a href="https://2014.igem.org/Team:NTNU_Trondheim/Project" class="flyout-toggle"><span> </span></a>
<ul class="flyout">
<ul class="flyout">
Line 252: Line 252:
<div>
<div>
<div class="twelve columns">
<div class="twelve columns">
-
<h3>Flux balance analysis of the impact of increasing CO2-uptake and introducing Glucose Oxidase into Synechocystis sp. PCC6803 </h3>
+
<h3>Flux balance analysis of the impact of increasing CO<sub>2</sub>-uptake and introducing Glucose Oxidase into Synechocystis sp. PCC6803 </h3>
<p>  
<p>  
-
To aid in the identification of possible target genes for increasing CO2 uptake in Synechocystis, <a href="http://www.ploscompbiol.org/article/info%3Adoi%2F10.1371%2Fjournal.pcbi.1003081">the most recent model of the metabolic network in Synechocystis</a> was acquired and analyzed with  
+
To aid in the identification of possible target genes for increasing CO<sub>2</sub>uptake in Synechocystis, <a href="http://www.ploscompbiol.org/article/info%3Adoi%2F10.1371%2Fjournal.pcbi.1003081">the most recent model of the metabolic network in Synechocystis</a> was acquired and analyzed with  
-
<a href="http://www.nature.com/nbt/journal/v28/n3/abs/nbt.1614.html">flux balance analysis</a>. The approach initially chosen was simply to fix the CO2 uptake flux above that achieved in the optimal FBA solution, and examining the effect on the growth rate (figure 1).  
+
<a href="http://www.nature.com/nbt/journal/v28/n3/abs/nbt.1614.html">flux balance analysis</a>. Computation was performed with the <a href="http://opencobra.sourceforge.net/openCOBRA/Welcome.html">COBRA toolbox</a> in Matlab 2014a.The approach initially chosen was simply to fix the CO<sub>2</sub> uptake flux above that achieved in the optimal FBA solution, and examining the effect on the growth rate (figure 1).  
-
<div class="col4"><a href="https://static.igem.org/mediawiki/2014/0/0b/Growth.jpg"> </a><img src="https://static.igem.org/mediawiki/2014/0/0b/Growth.jpg" width="500">
+
<div class="col4"><a href="https://static.igem.org/mediawiki/2014/0/0b/Growth.jpg"> </a><img src="https://static.igem.org/mediawiki/2014/0/0b/Growth.jpg" width="1000">
-
  <p style="text-align:center; color:black; "><b> Figure:</b> Synechocystis growth rate as a function of CO2 uptake flux in an optimal FBA solution.</p> </div>
+
  <p style="text-align:center; color:black; "><b> Figure 1:</b> Synechocystis growth rate as a function of CO<sub>2</sub> uptake flux in an optimal FBA solution.</p> </div>
  <p>
  <p>
-
The growth rate increases from 0 as the CO2 uptake flux is increased from zero to its optimal value (1.19), and decreases after this point. The decrease in growth rate after this optimal CO2 uptake flux is linear with a slope of -0.0083. As the growth rate decreases with increasing CO2 fixaton after this point,
+
The growth rate increases from 0 as the CO<sub>2</sub> uptake flux is increased from zero to its optimal value (1.2 mmol/gDW/h), and decreases after this point. The decrease in growth rate after this optimal CO<sub>2</sub> uptake flux is linear with a slope of -0.0083. As the growth rate decreases with increasing CO<sub>2</sub> fixaton after this point,
-
this means that any metabolic reaction added to the model that increases the CO2 fixation without also increasing the growth rate, will reduce the growth rate by at least -0.0089/h per unit of flux the CO2 uptake is increased by.  
+
this means that any metabolic reaction added to the model that increases the CO<sub>2</sub> fixation without also increasing the growth rate, will reduce the growth rate by at least -0.0089/h per unit of flux the CO<sub>2</sub> uptake is increased by.  
</p>
</p>
<p>The glucose oxidase gene was added to the model in order to examine its effect on cellular growth. This was done by implementing the reaction:</p>
<p>The glucose oxidase gene was added to the model in order to examine its effect on cellular growth. This was done by implementing the reaction:</p>
<p> <img src="https://static.igem.org/mediawiki/2014/6/62/CodeCogsEqn_%282%29.png"></p>
<p> <img src="https://static.igem.org/mediawiki/2014/6/62/CodeCogsEqn_%282%29.png"></p>
-
<p> Since the metabolite D-glucono-1,5-lactone was not present in the original metabolic model, this issue was addressed by the addition of a passive export reaction. The flux through this reaction was then increased from 0, and the CO2 uptake predicted by the FBA solution recorded (figure 2)</p>
+
<p> Since the metabolite D-glucono-1,5-lactone was not present in the original metabolic model, this issue was addressed by the addition of a passive export reaction. The flux through this reaction was then increased from 0, and the CO<sub>2</sub> uptake predicted by the FBA solution recorded (figure 2)</p>
-
<div class="col4"><a href="https://static.igem.org/mediawiki/2014/6/61/C02_uptake.jpg"> </a><img src="https://static.igem.org/mediawiki/2014/6/61/C02_uptake.jpg" width="500">
+
<div class="col4"><a href="https://static.igem.org/mediawiki/2014/6/61/C02_uptake.jpg"> </a><img src="https://static.igem.org/mediawiki/2014/6/61/C02_uptake.jpg" width="1000">
-
  <p style="text-align:center; color:black; "><b> Figure 2:</b> CO2 uptake as a function of glucose oxidase flux in an optimal FBA solution.</p> </div>
+
  <p style="text-align:center; color:black; "><b> Figure 2:</b> CO<sub>2</sub> uptake as a function of glucose oxidase flux in an optimal FBA solution.</p> </div>
   
   
-
  <p> The FBA solution then predicts that the CO2 uptake increases as the glucose oxidase flux increases. In accordance with the previous observation that the reduction in growth rate would be at least -0.0089/h per per unit of the flux the CO2 uptake was increased by, the reduction in growth rate was found to be -0.0533 per unit of CO2 uptake flux.</p>
+
  <p> The FBA solution then predicts that the CO<sub>2</sub> uptake increases as the glucose oxidase flux increases. In accordance with the previous observation that the reduction in growth rate would be at least -0.0089/h per per unit of the flux the CO<sub>2</sub> uptake was increased by, the reduction in growth rate was found to be -0.0533 per unit of CO<sub>2</sub> uptake flux.</p>
</p>
</p>
</div>
</div>
Line 283: Line 283:
</p>
</p>
<p>
<p>
-
The system was simulated deterministically using <a href="http://www.nanodrop.com/library/nd-1000-v3.7-users-manual-8.5x11.pdf">Cain</a> with the parameters in table 1. Parameters were taken from  <a href="http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2996131/?tool=pubmed">[1]</a>.
+
The system was simulated deterministically using <a href="http://cain.sourceforge.net/">Cain</a> with the parameters in table 1. Parameters were taken from  <a href="http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2996131/?tool=pubmed">[1]</a>.
</p>
</p>
Line 297: Line 297:
<tr>
<tr>
<td>k<sub>1</sub>
<td>k<sub>1</sub>
-
</td><td>0.23 nM minute <sup>-1</sup>
+
</td><td>0.23 nM min <sup>-1</sup>
</td><td>Transcription of lacI mRNA
</td><td>Transcription of lacI mRNA
</td></tr>
</td></tr>
<tr>
<tr>
<td>k<sub>2</sub>
<td>k<sub>2</sub>
-
</td><td>15 minute <sup>-1</sup>
+
</td><td>15 min <sup>-1</sup>
</td><td>Translation of LacI
</td><td>Translation of LacI
</td></tr>
</td></tr>
Line 312: Line 312:
<tr>
<tr>
<td>k<sub>4</sub>
<td>k<sub>4</sub>
-
</td><td>30 minute <sup>-1</sup>
+
</td><td>30 min <sup>-1</sup>
</td><td>Translation of GOx
</td><td>Translation of GOx
</td></tr>
</td></tr>
<tr>
<tr>
<td>k<sub>5</sub>
<td>k<sub>5</sub>
-
</td><td>0.92 minute <sup>-1</sup>
+
</td><td>0.92 min <sup>-1</sup>
</td><td>Passive IPTG influx
</td><td>Passive IPTG influx
</td></tr>
</td></tr>
<tr>
<tr>
<td>k<sub>-5</sub>
<td>k<sub>-5</sub>
-
</td><td>0.92 minute <sup>-1</sup>
+
</td><td>0.92 min <sup>-1</sup>
</td><td>Passive IPTG efflux
</td><td>Passive IPTG efflux
</td></tr>
</td></tr>
<tr>
<tr>
<td>k<sub>6</sub>
<td>k<sub>6</sub>
-
</td><td>50 nM <sup>-1</sup> minute <sup>-1</sup>
+
</td><td>50 nM <sup>-1</sup> min <sup>-1</sup>
</td><td>LacI dimerization
</td><td>LacI dimerization
</td></tr>
</td></tr>
Line 333: Line 333:
<tr>
<tr>
<td>k<sub>-6</sub>
<td>k<sub>-6</sub>
-
</td><td>1e-3  minute <sup>-1</sup>
+
</td><td>1e-3  min <sup>-1</sup>
</td><td>LacI dimer dissociation  
</td><td>LacI dimer dissociation  
</td></tr>
</td></tr>
Line 340: Line 340:
<tr>
<tr>
<td>k<sub>7</sub>
<td>k<sub>7</sub>
-
</td><td>3e-7  nM<sup>-2</sup> minute <sup>-1</sup>
+
</td><td>3e-7  nM<sup>-2</sup> min <sup>-1</sup>
</td><td>LacI dimer IPTG binding
</td><td>LacI dimer IPTG binding
</td></tr>
</td></tr>
Line 347: Line 347:
<tr>
<tr>
<td>k<sub>-7</sub>
<td>k<sub>-7</sub>
-
</td><td>12 minute <sup>-1</sup>
+
</td><td>12 min <sup>-1</sup>
</td><td>LacI dimer IPTG unbinding
</td><td>LacI dimer IPTG unbinding
</td></tr>
</td></tr>
Line 353: Line 353:
<tr>
<tr>
<td>k<sub>8</sub>
<td>k<sub>8</sub>
-
</td><td>960 nM<sup>-1</sup> minute <sup>-1</sup>
+
</td><td>960 nM<sup>-1</sup> min <sup>-1</sup>
</td><td>Repression of GOx
</td><td>Repression of GOx
</td></tr>
</td></tr>
Line 360: Line 360:
<tr>
<tr>
<td>k<sub>-8</sub>
<td>k<sub>-8</sub>
-
</td><td>2.4  minute <sup>-1</sup>
+
</td><td>2.4  min <sup>-1</sup>
</td><td>Derepression of GOx
</td><td>Derepression of GOx
</td></tr>
</td></tr>
Line 367: Line 367:
<tr>
<tr>
<td>d<sub>1</sub>
<td>d<sub>1</sub>
-
</td><td> 0.2 minute <sup>-1</sup>
+
</td><td> 0.2 min <sup>-1</sup>
</td><td>Degradation of GOx
</td><td>Degradation of GOx
</td></tr>
</td></tr>
Line 373: Line 373:
<tr>
<tr>
<td>d<sub>2</sub>
<td>d<sub>2</sub>
-
</td><td> 0.2 minute <sup>-1</sup>
+
</td><td> 0.2 min <sup>-1</sup>
</td><td>Degradation of LacI
</td><td>Degradation of LacI
</td></tr>
</td></tr>
Line 379: Line 379:
<tr>
<tr>
<td>d<sub>3</sub>
<td>d<sub>3</sub>
-
</td><td> 0.462 minute <sup>-1</sup>
+
</td><td> 0.462 min <sup>-1</sup>
</td><td>Degradation of GOx mRNA
</td><td>Degradation of GOx mRNA
</td></tr>
</td></tr>
Line 385: Line 385:
<tr>
<tr>
<td>d<sub>4</sub>
<td>d<sub>4</sub>
-
</td><td> 0.462 minute <sup>-1</sup>
+
</td><td> 0.462 min <sup>-1</sup>
</td><td>Degradation of LacI mRNA
</td><td>Degradation of LacI mRNA
</td></tr>
</td></tr>
Line 391: Line 391:
<tr>
<tr>
<td>d<sub>5</sub>
<td>d<sub>5</sub>
-
</td><td> 0.2 minute <sup>-1</sup>
+
</td><td> 0.2 min <sup>-1</sup>
</td><td>Degradation of LacI dimer
</td><td>Degradation of LacI dimer
</td></tr>
</td></tr>
Line 398: Line 398:
<tr>
<tr>
<td>d<sub>6</sub>
<td>d<sub>6</sub>
-
</td><td> 0.2 minute <sup>-1</sup>
+
</td><td> 0.2 min <sup>-1</sup>
</td><td>Degradation of LacI dimer-IPTG complex
</td><td>Degradation of LacI dimer-IPTG complex
</td></tr>
</td></tr>
-
 
-
 
</table>
</table>
<p>
<p>
-
A .zip file with the model file can be found a <href="https://static.igem.org/mediawiki/2014/1/1d/Model.zip">here</a>.
+
A .zip file with the model file can be found <a href="https://static.igem.org/mediawiki/2014/1/1d/Model.zip">here</a>. Simulating the system with various levels of IPTG (figure 3) yields increasing steady state levels of GOx.  
</p>
</p>
 +
<div class="col4"><a href="https://static.igem.org/mediawiki/2014/7/71/Iptg.png"> </a><img src="https://static.igem.org/mediawiki/2014/7/71/Iptg.png" width="1000">
 +
<p style="text-align:center; color:black; "><b> Figure 3:</b> Cellular GOx concentration as a function of IPTG concentration</p> </div>
</div>
</div>
</div>
</div>

Latest revision as of 23:53, 17 October 2014

Team:NTNU_Trondheim/Project/Modelling - 2014.igem.org

 

Team:NTNU_Trondheim/Home

From 2014.igem.org

NTNU Genetically Engineered Machines

Modelling

Flux balance analysis of the impact of increasing CO2-uptake and introducing Glucose Oxidase into Synechocystis sp. PCC6803

To aid in the identification of possible target genes for increasing CO2uptake in Synechocystis, the most recent model of the metabolic network in Synechocystis was acquired and analyzed with flux balance analysis. Computation was performed with the COBRA toolbox in Matlab 2014a.The approach initially chosen was simply to fix the CO2 uptake flux above that achieved in the optimal FBA solution, and examining the effect on the growth rate (figure 1).

Figure 1: Synechocystis growth rate as a function of CO2 uptake flux in an optimal FBA solution.

The growth rate increases from 0 as the CO2 uptake flux is increased from zero to its optimal value (1.2 mmol/gDW/h), and decreases after this point. The decrease in growth rate after this optimal CO2 uptake flux is linear with a slope of -0.0083. As the growth rate decreases with increasing CO2 fixaton after this point, this means that any metabolic reaction added to the model that increases the CO2 fixation without also increasing the growth rate, will reduce the growth rate by at least -0.0089/h per unit of flux the CO2 uptake is increased by.

The glucose oxidase gene was added to the model in order to examine its effect on cellular growth. This was done by implementing the reaction:

Since the metabolite D-glucono-1,5-lactone was not present in the original metabolic model, this issue was addressed by the addition of a passive export reaction. The flux through this reaction was then increased from 0, and the CO2 uptake predicted by the FBA solution recorded (figure 2)

Figure 2: CO2 uptake as a function of glucose oxidase flux in an optimal FBA solution.

The FBA solution then predicts that the CO2 uptake increases as the glucose oxidase flux increases. In accordance with the previous observation that the reduction in growth rate would be at least -0.0089/h per per unit of the flux the CO2 uptake was increased by, the reduction in growth rate was found to be -0.0533 per unit of CO2 uptake flux.

Simulation of genetic circuit

The construct chosen to express glucose odxidase (GOx) is a lac promotor regulated by LacI. GOx is therefore inducible by the addition of IPTG. The genetic circuit can be modelled by the following chemical reaction equations:

The system was simulated deterministically using Cain with the parameters in table 1. Parameters were taken from [1].

Table 1: Parameters for the inducible circuit

Reaction rate constant Propensity Comment
k1 0.23 nM min -1 Transcription of lacI mRNA
k2 15 min -1 Translation of LacI
k3 0.5 Transcription of GOx mRNA
k4 30 min -1 Translation of GOx
k5 0.92 min -1 Passive IPTG influx
k-5 0.92 min -1 Passive IPTG efflux
k6 50 nM -1 min -1 LacI dimerization
k-6 1e-3 min -1 LacI dimer dissociation
k7 3e-7 nM-2 min -1 LacI dimer IPTG binding
k-7 12 min -1 LacI dimer IPTG unbinding
k8 960 nM-1 min -1 Repression of GOx
k-8 2.4 min -1 Derepression of GOx
d1 0.2 min -1 Degradation of GOx
d2 0.2 min -1 Degradation of LacI
d3 0.462 min -1 Degradation of GOx mRNA
d4 0.462 min -1 Degradation of LacI mRNA
d5 0.2 min -1 Degradation of LacI dimer
d6 0.2 min -1 Degradation of LacI dimer-IPTG complex

A .zip file with the model file can be found here. Simulating the system with various levels of IPTG (figure 3) yields increasing steady state levels of GOx.

Figure 3: Cellular GOx concentration as a function of IPTG concentration