Team:Evry/Model3

From 2013.igem.org

(Difference between revisions)
 
(21 intermediate revisions not shown)
Line 10: Line 10:
<p>
<p>
-
This model ams to simulate the entire pathway of the enterobactin production in our Iron-coli. It includes the synthetic sensing system, and the chemical reactions leading to the enterobactin.<br/>
+
This model aims to simulate the entire pathway of the enterobactin production in our Iron Coli. It includes the synthetic sensing system, and the chemical reactions leading to the enterobactin.<br/>
This model has two goals:
This model has two goals:
<ul>
<ul>
-
<li>Have an idea of how long it would take to produce siderophores, from the very begining: the iron sensing</li>
+
<li>To the delay needed to produce siderophores starting from the iron-sensing</li>
-
<li>Determine the ideal plasmid configuration</li>
+
<li>To determine the ideal plasmid configuration</li>
</ul>
</ul>
The enterobactin production regulation is based on two consecutives inhibitions, which, in the end, is an activator with a certain delay. The model will follow this principle.
The enterobactin production regulation is based on two consecutives inhibitions, which, in the end, is an activator with a certain delay. The model will follow this principle.
</p>
</p>
 +
<h2>Assumptions</h2>
<h2>Assumptions</h2>
<ul>
<ul>
-
<li>The only limiting reagant in the pathway is the chorismic acid, which we consider being in excess</li>
+
<li>The only limiting reagant in the pathway is the chorismic acid, which we consider being in excess amount</li>
-
<li>The FUR is also considered in exccess</li>
+
<li>The FUR is also considered in excess amount</li>
-
<li>All the enzymes (and the mRNA involved) needed are assimilted to a single variable, since they all play equivalent roles</li>
+
<li>All the enzymes (and the mRNA involved) needed are assimilated to a single variable, since they all play equivalent roles</li>
</ul>
</ul>
Line 53: Line 54:
<img src="https://static.igem.org/mediawiki/2013/f/fe/Reg3.png"/><br/>
<img src="https://static.igem.org/mediawiki/2013/f/fe/Reg3.png"/><br/>
We chose to model the iron input in the bacteria using a linear function of the external iron concentration <i>Ferext</i>, the factor <i>p</i> being the cell-wall permeability for iron.<br/>
We chose to model the iron input in the bacteria using a linear function of the external iron concentration <i>Ferext</i>, the factor <i>p</i> being the cell-wall permeability for iron.<br/>
-
The FUR on the other hand, is produced by the bacteria. It's evolution can also be considered linerar, using a mean production rate <i>Fur0</i>.<br/>
+
The FUR on the other hand is produced by the bacteria. Its evolution can also be considered as linear, using a mean production rate <i>Fur0</i>.<br/>
<img src="https://static.igem.org/mediawiki/2013/5/51/Regfer.png"/><br/>
<img src="https://static.igem.org/mediawiki/2013/5/51/Regfer.png"/><br/>
<img src="https://static.igem.org/mediawiki/2013/6/68/Regfur.png"/><br/>
<img src="https://static.igem.org/mediawiki/2013/6/68/Regfur.png"/><br/>
-
In this model, we only track the free Fe-FUR and not those which are attached to a FUR Binding Site. As <i>LacI</i> is the number of inhibited LacI, we can use this number to express how much Fe-FUR does bind to a FBS per unit of time.<br/>
+
In this model, we only track the free Fe-FUR complex and not those attached to a FUR Binding Site. As <i>LacI</i> is the number of inhibited LacI, we can use this number to express how much Fe-FUR does bind to a FBS per unit of time.<br/>
<img src="https://static.igem.org/mediawiki/2013/6/6d/Regfefur.png"/>
<img src="https://static.igem.org/mediawiki/2013/6/6d/Regfefur.png"/>
</p>
</p>
<p>
<p>
<u>LacI-LacO system:</u><br/>
<u>LacI-LacO system:</u><br/>
-
To simulate the inhibition phenomenon, we chose to use our <a href="https://2013.igem.org/Team:Evry/LogisticFunctions">logistic function</a> under its differential form. Since it is the FeFUR that represses it, thhe LacI can be expressed as a logistic fuction of the FeFUR:<br/>
+
To simulate the inhibition phenomenon, we chose to use our <a href="https://2013.igem.org/Team:Evry/LogisticFunctions">logistic function</a> under its differential form. Since it is the Fe-FUR that represses it, the LacI can be expressed as a logistic fuction of the Fe-FUR:<br/>
<img src="https://static.igem.org/mediawiki/2013/4/49/Dlaci.png"/><br/>
<img src="https://static.igem.org/mediawiki/2013/4/49/Dlaci.png"/><br/>
-
Ki1 is a non-dimensional parameter that repesents the inhibition power, and Kf is the fixation rate of the FeFUR on the FBS. Finally, Nbrpla1 is the number of pasmids containing the LacI.<br/>
+
Ki1 is a non-dimensional parameter that repesents the inhibition power, and Kf is the fixation rate of the Fe-FUR on the FBS. Finally, Nbrpla1 is the number of pasmids containing the LacI.<br/>
In the same way, LacO is modeled with a logistic funtion. LacO is repressed by LacI:<br/>
In the same way, LacO is modeled with a logistic funtion. LacO is repressed by LacI:<br/>
<img src="https://static.igem.org/mediawiki/2013/f/f2/Dlaco.png"/><br/>
<img src="https://static.igem.org/mediawiki/2013/f/f2/Dlaco.png"/><br/>
Line 88: Line 89:
<th>Unite</th>
<th>Unite</th>
<th>Description</th>
<th>Description</th>
-
<th>Reference</th>
+
 
</tr>
</tr>
<tr width="100%">
<tr width="100%">
Line 95: Line 96:
<td>min^-1</td>
<td>min^-1</td>
<td>Permeability of cell wall</td>
<td>Permeability of cell wall</td>
-
<td>-</td>
+
 
</tr>
</tr>
<tr width="100%">
<tr width="100%">
Line 102: Line 103:
<td>M^-1.s^-1</td>
<td>M^-1.s^-1</td>
<td>Formation constant of FeFur complex</td>
<td>Formation constant of FeFur complex</td>
-
<td>-</td>
+
 
</tr>
</tr>
<tr width="100%">
<tr width="100%">
Line 109: Line 110:
<td>min^-1</td>
<td>min^-1</td>
<td>-</td>
<td>-</td>
-
<td>-</td>
+
 
</tr>
</tr>
<tr width="100%">
<tr width="100%">
Line 116: Line 117:
<td>min^-1</td>
<td>min^-1</td>
<td>translation rate</td>
<td>translation rate</td>
-
<td>-</td>
+
 
<tr width="100%">
<tr width="100%">
<td>KT</td>
<td>KT</td>
Line 122: Line 123:
<td>M.min^-1</td>
<td>M.min^-1</td>
<td>transcription rate</td>
<td>transcription rate</td>
-
<td>-</td>
+
 
</tr>
</tr>
<tr width="100%">
<tr width="100%">
Line 129: Line 130:
<td>mol^-1</td>
<td>mol^-1</td>
<td>Avogadro's constant</td>
<td>Avogadro's constant</td>
-
<td>-</td>
+
 
</tr>
</tr>
<tr width="100%">
<tr width="100%">
Line 136: Line 137:
<td>m^3</td>
<td>m^3</td>
<td>Volume of a call</td>
<td>Volume of a call</td>
-
<td>-</td>
+
 
</tr>
</tr>
<tr width="100%">
<tr width="100%">
Line 143: Line 144:
<td>min^-1</td>
<td>min^-1</td>
<td>mRNA degradation rate</td>
<td>mRNA degradation rate</td>
-
<td>-</td>
+
 
</tr>
</tr>
<tr width="100%">
<tr width="100%">
Line 150: Line 151:
<td>min^-1</td>
<td>min^-1</td>
<td>Enzyme degradation rate</td>
<td>Enzyme degradation rate</td>
-
<td>-</td>
+
 
</tr>
</tr>
<tr width="100%">
<tr width="100%">
Line 157: Line 158:
<td>min^-1</td>
<td>min^-1</td>
<td>fixation rate of FeFUR</td>
<td>fixation rate of FeFUR</td>
-
<td>-</td>
+
 
</tr>
</tr>
<tr width="100%">
<tr width="100%">
Line 164: Line 165:
<td>mM.min^⁻1</td>
<td>mM.min^⁻1</td>
<td>Fur Production</td>
<td>Fur Production</td>
-
<td>-</td>
+
 
</tr>
</tr>
<tr width="100%">
<tr width="100%">
Line 171: Line 172:
<td>mM.min^⁻1</td>
<td>mM.min^⁻1</td>
<td>Chorismate production</td>
<td>Chorismate production</td>
-
<td>-</td>
 
</tr>
</tr>
</table>
</table>
Line 181: Line 181:
<td>5550</td>
<td>5550</td>
<td>min^-1</td>
<td>min^-1</td>
-
<td>-</td>
+
<td>PMID:2144454</td>
-
<td>-</td>
+
</tr>
</tr>
<tr width="100%">
<tr width="100%">
Line 188: Line 187:
<td>600</td>
<td>600</td>
<td>min^-1</td>
<td>min^-1</td>
-
<td>-</td>
+
<td>PMID:2139796</td>
-
<td>-</td>
+
 
</tr>
</tr>
<tr width="100%">
<tr width="100%">
Line 195: Line 194:
<td>173</td>
<td>173</td>
<td>min^-1</td>
<td>min^-1</td>
-
<td>-</td>
+
<td></td>
-
<td>-</td>
+
 
</tr>
</tr>
<tr width="100%">
<tr width="100%">
<td>KmA</td>
<td>KmA</td>
<td>300</td>
<td>300</td>
-
<td>M</td>
+
<td>μM</td>
-
<td>-</td>
+
<td>-</td>
<td>-</td>
 +
 +
</tr>
</tr>
<tr width="100%">
<tr width="100%">
<td>KmB</td>
<td>KmB</td>
<td>14.7</td>
<td>14.7</td>
-
<td>M</td>
+
<td>μM</td>
-
<td>-</td>
+
<td>PMID:2139796</td>
-
<td>-</td>
+
 
</tr>
</tr>
<tr width="100%">
<tr width="100%">
<td>KmC</td>
<td>KmC</td>
<td>14</td>
<td>14</td>
-
<td>M</td>
+
<td>μM</td>
-
<td>-</td>
+
<td>-</td>
<td>-</td>
 +
</tr>
</tr>
</table><br/>
</table><br/>
Line 229: Line 229:
<p>
<p>
-
<u>High copy?</u>
+
<u>High copy?</u><br/>
-
Here are 9 gaphs pointing out the relevance of different plasmid configurations. The goal of this simulation is to try answerig a biological question: "should we copy the plLacO plasmids or the LacI plasmids?"<br/>The horizontal axes are the time (minutes) and the iron concentration (logarithmic scale). The vertical axis is the Enzyme production (mmol/L).<br/>
+
Here are 9 gaphs pointing out the relevance of different plasmid configurations. The goal of this simulation is to try answering a biological question: "should we increase the number of plLacO plasmids or the LacI plasmids?"<br/>The horizontal axes are the time (minutes) and the iron concentration (logarithmic scale). The vertical axis is the Enzyme production (mmol/L).<br/>
<div align="center"><img src="https://static.igem.org/mediawiki/2013/d/d5/MetabolismResult1.png"/></div>
<div align="center"><img src="https://static.igem.org/mediawiki/2013/d/d5/MetabolismResult1.png"/></div>
-
The graph is clear: the most important plasmid is the one containing the LacI.
+
The graph is clear: the most important plasmid is the one containing the LacI (the direct iron sensor).
</p>
</p>
<p>
<p>
-
<u>About the KI1/KI2:</u>
+
<u>About the KI1/KI2:</u><br/>
Here are 9 other graphs with a set of different values for KI1 and KI2.
Here are 9 other graphs with a set of different values for KI1 and KI2.
<div align="center"><img src="https://static.igem.org/mediawiki/2013/f/fd/MetaTabregki.png"/></div>
<div align="center"><img src="https://static.igem.org/mediawiki/2013/f/fd/MetaTabregki.png"/></div>
 +
</p><br/><br/>
-
<a id="Conclusion"></a>
 
-
<h2>Conclusion</h2>
 
<p>
<p>
-
</p>
+
<u>Temporal answers</u><br/>
 +
 
 +
 
 +
This curve integrates the chemical pathway developped <a href"https://2013.igem.org/Team:Evry/ChemicalTools">here</a>.
 +
<div><img src="https://static.igem.org/mediawiki/2013/3/34/Systemboss.png"/></div>
 +
<div align="center"><img src="https://static.igem.org/mediawiki/2013/5/5e/MetaEntenz.png" width="900px"/></div>
 +
There is a chemical delay in the enterobactin production, which is about the same length as the bacterial delay.
 +
</p><br/>
 +
 
 +
 
 +
 
 +
<br>
 +
<br>
 +
 
 +
<p id="model_galenic">
 +
<u>Impact on the galenic formulation</u><br/>
 +
The delay of siderophore production is predicted to be way to long if you choose for the flush strategy. As a consequence, these result directly impacted our chelating strategy. Thus, our final decision is to try to implant the bacteria in the intestins (here in the jejunum) to statically enhance its growth and let it enough time to recover for a dessicated environment.
 +
</p><br/>
-
<div id="citation_box">
+
<a href="https://static.igem.org/mediawiki/2013/b/be/Metabolicmodel.zip">Download the model</a>
-
<p id="references">References:</p>
+
-
<ol>
+
-
</ol>
 
-
</div>
 

Latest revision as of 15:55, 22 October 2013

Iron coli project

Metabolic model

Overview

This model aims to simulate the entire pathway of the enterobactin production in our Iron Coli. It includes the synthetic sensing system, and the chemical reactions leading to the enterobactin.
This model has two goals:

  • To the delay needed to produce siderophores starting from the iron-sensing
  • To determine the ideal plasmid configuration
The enterobactin production regulation is based on two consecutives inhibitions, which, in the end, is an activator with a certain delay. The model will follow this principle.

Assumptions

  • The only limiting reagant in the pathway is the chorismic acid, which we consider being in excess amount
  • The FUR is also considered in excess amount
  • All the enzymes (and the mRNA involved) needed are assimilated to a single variable, since they all play equivalent roles

Model Description

Variables:

  • [Fe] : Iron concentration inside the bacteria
  • [Fur] : FUR concentration inside the bacteria
  • [FeFur] : Iron-FUR complex concentration inside the bacteria
  • LacI : Number of inhibited LacI
  • LacO : Number of non-inhibited LacO
  • [mRNA]: mRNA (from LacO) concentration
  • [Enz] : Enzyme concentration : EntA,-B,-C,-D,-E,-F
All those concentrations are expressed in mmol/L

Fe, FUR and FeFUR:
The iron-FUR complex is simply formed that way:
 
We reduced this equation to:

Which is not annoying, since we just have to divide our [FeFur] by to to get the real complex concentration.
We can easily write down both the formation (v) and the dissociation (v') speed:

We chose to model the iron input in the bacteria using a linear function of the external iron concentration Ferext, the factor p being the cell-wall permeability for iron.
The FUR on the other hand is produced by the bacteria. Its evolution can also be considered as linear, using a mean production rate Fur0.


In this model, we only track the free Fe-FUR complex and not those attached to a FUR Binding Site. As LacI is the number of inhibited LacI, we can use this number to express how much Fe-FUR does bind to a FBS per unit of time.

LacI-LacO system:
To simulate the inhibition phenomenon, we chose to use our logistic function under its differential form. Since it is the Fe-FUR that represses it, the LacI can be expressed as a logistic fuction of the Fe-FUR:

Ki1 is a non-dimensional parameter that repesents the inhibition power, and Kf is the fixation rate of the Fe-FUR on the FBS. Finally, Nbrpla1 is the number of pasmids containing the LacI.
In the same way, LacO is modeled with a logistic funtion. LacO is repressed by LacI:

Ki2 is the inhibition power and Nbrpla2 is the number of plasmimds containing LacO.
LacI and LacO are both ruled by a normal logistic function. If we were to track the number of expressed LacI/LacO, we would be using two inverted logistic fuctions to model a double inverter. The thing is, since LacI represents the number of repressed genes and LacO the number of expressed genes, the double inverter is still there, but the calculations are easier.

mRNA and Enzymes:
The [mRNA] and [Enz] equations are alike. The prodction rates are Kr for the mRNA and Kp for the enzymmes, and both variables have a negative degadation term:

Finally:

Parameters

Name Value Unite Description
p 0.1 min^-1 Permeability of cell wall
KfeFUR 0.01 M^-1.s^-1 Formation constant of FeFur complex
Dff 0.001 min^-1 -
Kp 0.5 min^-1 translation rate
KT 2.0 M.min^-1 transcription rate
milliNa 6.02.10^20 mol^-1 Avogadro's constant
V 6.5.10^-16 m^3 Volume of a call
Dmrna 0.001 min^-1 mRNA degradation rate
Denz 0.001 min^-1 Enzyme degradation rate
Kf 10^-4 min^-1 fixation rate of FeFUR
Fur0 0.01 mM.min^⁻1 Fur Production
vE1 0.01 mM.min^⁻1 Chorismate production
Enzymatic Parameters:
KcatA 5550 min^-1 PMID:2144454
KcatB 600 min^-1 PMID:2139796
KcatC 173 min^-1
KmA 300 μM -
KmB 14.7 μM PMID:2139796
KmC 14 μM -

Note that we don't have any informations about the KI1 and KI2 parameters, since they are purely artificial. We thus had to run simulations to study the effects of different sets of values on the model.

Results

High copy?
Here are 9 gaphs pointing out the relevance of different plasmid configurations. The goal of this simulation is to try answering a biological question: "should we increase the number of plLacO plasmids or the LacI plasmids?"
The horizontal axes are the time (minutes) and the iron concentration (logarithmic scale). The vertical axis is the Enzyme production (mmol/L).

The graph is clear: the most important plasmid is the one containing the LacI (the direct iron sensor).

About the KI1/KI2:
Here are 9 other graphs with a set of different values for KI1 and KI2.



Temporal answers
This curve integrates the chemical pathway developped here.

There is a chemical delay in the enterobactin production, which is about the same length as the bacterial delay.




Impact on the galenic formulation
The delay of siderophore production is predicted to be way to long if you choose for the flush strategy. As a consequence, these result directly impacted our chelating strategy. Thus, our final decision is to try to implant the bacteria in the intestins (here in the jejunum) to statically enhance its growth and let it enough time to recover for a dessicated environment.


Download the model