Team:Evry/Modeltr2

From 2013.igem.org

(Difference between revisions)
 
(18 intermediate revisions not shown)
Line 14: Line 14:
<h2>Observations</h2>
<h2>Observations</h2>
<p>
<p>
-
Once our genetically modified bacteria are released in the duodenum, they produce siderophores to chelate the solved iron, thus making it unavailable for intestinal absorption. Then, <b>they eventually flush out of the duodenum</b>. The main hypothesis in this model is that the bacteria don't colonize the duodenum : they only flow through. <br/>
+
Once our genetically modified bacteria are released in the duodenum lumen, they produce siderophores to chelate the solved iron, thus making it unavailable for intestinal absorption. Then, <b>they eventually flush out of the duodenum</b>. The main hypothesis in this model is that the bacteria do not colonize the duodenum : they only flow through. They do not even have time to grow, for the time required to flush through is close to 40 seconds.<br/>
-
The goal of this model is to measure how efficient could this form of treatment be. Because too much parameters remain unknown, it is a theoretical simulation which will <b>not</b> give any numerical results.  
+
Because some of the parameters remain unknown, it is a theoretical simulation which will <b>not</b> give any numerical results, only qualitative results.
</p>
</p>
-
<h2>Goals</h2>
+
<h2>Goal</h2>
<p>
<p>
This model was made to check if our first strategy was viable. It aims to answer the following question: <br/>
This model was made to check if our first strategy was viable. It aims to answer the following question: <br/>
Line 30: Line 30:
<li>Our bacteria don't settle in the duodenum</li>
<li>Our bacteria don't settle in the duodenum</li>
<li>No regulation in the patient's iron absorption</li>
<li>No regulation in the patient's iron absorption</li>
-
<li>Constant iron flow</li>
+
<li>Constant iron flow in the duodenum lumen</li>
<li>Homogeneous fluid</li>
<li>Homogeneous fluid</li>
<li>The bacterial quantity is constant</li>  
<li>The bacterial quantity is constant</li>  
Line 41: Line 41:
<b><span style="color:#0000FF;">A</span></b> : Total quantity of iron absorbed by the duodenum (mol)
<b><span style="color:#0000FF;">A</span></b> : Total quantity of iron absorbed by the duodenum (mol)
<br/>
<br/>
-
<b><span style="color:#FF0000;">S</span></b> : Quantity of solubilized iron (mol)
+
<b><span style="color:#FF0000;">S</span></b> : Quantity of solved iron (mol)
<br/>
<br/>
<b>P</b> : Total quantity of enterobactin produced by our population of bacteria (mol)
<b>P</b> : Total quantity of enterobactin produced by our population of bacteria (mol)
Line 49: Line 49:
<b>N</b> : Number of bacteria
<b>N</b> : Number of bacteria
</p>
</p>
-
<img src="https://static.igem.org/mediawiki/2013/9/93/EquadiffModele1.png" width=300px /><br/>
+
<img src="https://static.igem.org/mediawiki/2013/1/1b/Duoeqcoul.png" width=300px /><br/>
<p>
<p>
 +
 +
<br/>
 +
<table border="1" style="border-collapse:collapse;">
 +
<tr>
 +
<td><b>S<sub>p</sub></b></td>
 +
<td>mol.s<sup>-1</sup></td>
 +
<td>Iron pulse</td>
 +
</tr>
 +
<tr>
 +
<td><b>v</b></td>
 +
<td>m.s<sup>-1</sup></td>
 +
<td>Chyme's flow average speed</td>
 +
</tr>
 +
<tr>
 +
<td><b>L</b></td>
 +
<td>m</td>
 +
<td>Duodenum length</td>
 +
</tr>
 +
<tr>
 +
<td><b>α</b></td>
 +
<td>s<sup>-1</sup></td>
 +
<td>Duodenum absorption rate</td>
 +
</tr>
 +
<tr>
 +
<td><b>K</b></td>
 +
<td>mol/s</td>
 +
<td>Activator Magnitude</td>
 +
</tr>
 +
<tr>
 +
<td><b>p</b></td>
 +
<td>mol.s<sup>-1</sup></td>
 +
<td>Value at zero of the activator</td>
 +
</tr>
 +
<tr>
 +
<td><b>h</b></td>
 +
<td>-</td>
 +
<td>Activator efficiency</td>
 +
</tr>
 +
<tr>
 +
<td><b>d</b></td>
 +
<td>mol</td>
 +
<td>Activator threshold</td>
 +
</tr>
 +
<tr>
 +
<td><b>δ</b></td>
 +
<td>mol<sup>-1</sup></td>
 +
<td>Chelation rate</td>
 +
</tr>
 +
</table>
 +
<br/><br/>
The graph on the right explains the reasoning: for instance, the arrow with a + between N and P means that the variation of P has a positive linear term in N.
The graph on the right explains the reasoning: for instance, the arrow with a + between N and P means that the variation of P has a positive linear term in N.
</p>
</p>
-
Where <img src="https://static.igem.org/mediawiki/2013/0/01/LogistiqueDuodenum.png" /> is our <a href="https://2013.igem.org/Team:Evry/LogisticFunctions">logistic function</a> (which we will here abusively call activator).<br/>
+
Where <img src="https://static.igem.org/mediawiki/2013/0/01/LogistiqueDuodenum.png" /> is our <a href="https://2013.igem.org/Team:Evry/LogisticFunctions">Logistic function</a> (which we will here abusively call activator).<br/>
<h2>Results</h2>
<h2>Results</h2>
<p>
<p>
 +
 +
<div align="center">
 +
<a id="Fig1"></a>
<div class="captionedPicture" style="width:100%;">
<div class="captionedPicture" style="width:100%;">
   <a title="Absorption" href="https://static.igem.org/mediawiki/2013/4/44/Model_absorption_treatment.png">
   <a title="Absorption" href="https://static.igem.org/mediawiki/2013/4/44/Model_absorption_treatment.png">
Line 65: Line 118:
   </div>
   </div>
  </div>
  </div>
 +
</div>
 +
<br/>
<br/>
-
<p>Figure 1 represent the iron absorbed by the duodenum during a meal. We can see the reduction of iron absorption in the duodenum</p>
+
<p>The <a href="#Fig1">Figure 1</a> represents the total iron absorbed by the duodenum during an average meal (with and without bacterial treatment), whereas the <a href="#Fig2">Figure 2</a> represents the instant ambient iron in the duodenum lumen(with and without bacterial treatment).
 +
</p>
<br/>
<br/>
 +
 +
<div align="center">
 +
<a id="Fig2"></a>
  <div class="captionedPicture" style="width:100%;">
  <div class="captionedPicture" style="width:100%;">
   <a title="Quantity" href="https://static.igem.org/mediawiki/2013/0/07/Model_iron_treatment.png">
   <a title="Quantity" href="https://static.igem.org/mediawiki/2013/0/07/Model_iron_treatment.png">
Line 76: Line 135:
   </div>
   </div>
  </div>
  </div>
-
<br/>
+
</div>
-
<p>Figure 2 represent iron disolved in the chyme. This graph point out the influence of bacteria on the middle. The blue curve decrease significantly under the bacteria flux. That underline the fact that the treatment don't affect the absorption directly. After 150 secondes, chyme is no longer in the duodenum, thus the iron absorption is reduce in both cases.<em>Baptiste: En fait je suis pas sûr que ça soit exactement ça, mais faudrait l'expliquer en tout cas même si c'est un artefact.</em></p>
+
-
 
+
 +
<br/>
 +
<p>
 +
The <a href="#Fig2">Figure 2</a> points out the influence of our bacteria on the middle. The blue curve decreases significantly under the bacterial flow. That underlines the fact that the treatment doesn't directly affect the absorption. After 150 seconds, the chyme is considered to progressively leave the duodenum (entering the jejunum), thus reducing the iron absorption in both cases, which finally leads to a steady state.
 +
</p>
<br/><br/>
<br/><br/>
-
<u>Parameters:</u><br/><br/>
+
<h2>Parameters:</h2>
-
The simulations was made with following parameters:<br/>
+
The simulations was made with following parameters:<br/><br/>
<table width="100%" border="1" style="border-collapse:collapse;">
<table width="100%" border="1" style="border-collapse:collapse;">
<tr width="100%">
<tr width="100%">
Line 94: Line 155:
<tr width="100%">
<tr width="100%">
<td>Sp</td>
<td>Sp</td>
-
<td>4.5*10^-9</td>
+
<td>4.5*10<sup>-9</sup></td>
<td>mol.s<sup>-1</sup></td>
<td>mol.s<sup>-1</sup></td>
<td>Iron pulse</td>
<td>Iron pulse</td>
Line 119: Line 180:
<td>Duodenum absorption rate</td>
<td>Duodenum absorption rate</td>
<td><a href="https://2013.igem.org/Team:Evry/Modeltr1">Disease model</a></td>
<td><a href="https://2013.igem.org/Team:Evry/Modeltr1">Disease model</a></td>
-
</tr>
 
-
<tr width="100%">
 
-
<td>δ</td>
 
-
<td>2.65*10^-8</td>
 
-
<td>mol<sup>-1</sup></td>
 
-
<td>Dimensional parameter</td>
 
-
<td><b>REF</b></td>
 
</tr>
</tr>
<tr width="100%">
<tr width="100%">
Line 143: Line 197:
<tr width="100%">
<tr width="100%">
<td>h</td>
<td>h</td>
-
<td>10^-5</td>
+
<td>10<sup>-5</sup></td>
<td>-</td>
<td>-</td>
<td>Activator efficiency</td>
<td>Activator efficiency</td>
Line 150: Line 204:
<tr width="100%">
<tr width="100%">
<td>d</td>
<td>d</td>
-
<td>10^-9</td>
+
<td>10<sup>-9</sup></td>
<td>mol</td>
<td>mol</td>
<td>Activator threshold</td>
<td>Activator threshold</td>
 +
<td><a href="#Ref">[4]</a></td>
 +
</tr>
 +
<tr width="100%">
 +
<td>δ</td>
 +
<td>2.65*10<sup>-8</sup></td>
 +
<td>mol<sup>-1</sup></td>
 +
<td>Chelation rate</td>
<td><a href="#Ref">[4]</a></td>
<td><a href="#Ref">[4]</a></td>
</tr>
</tr>
Line 158: Line 219:
</p>
</p>
<br/>
<br/>
-
<p>There is a slight difference between the two cases, which would possibly allow us to fit the treatment to the patient's characteristics.</p>
 
<h2>Conclusion</h2>
<h2>Conclusion</h2>
<p>
<p>
-
It is possible to significantly reduce intestinal iron intake if the patient takes one pill before or during each meal. This means that the patient would endure a lighter treatment : less bloodletting for people suffering from hemochromatosis, and less iron chelator's side effects for the thalassemia.
+
The combination of the two curves clearly shows that bacteria would be able to <b>neutralize a significant quantity of iron</b> before it is absorbed by the duodenum.<br/>
 +
It would thus be possible to approximately <b>divide the intestinal iron intake by 2</b> if the patient takes one pill during each meal. This means that the patient would endure a lighter treatment : less bloodletting for people suffering from hemochromatosis, and less iron chelator's side effects for the thalassemia.
</p>
</p>
<h2>Models and scripts</h2>
<h2>Models and scripts</h2>
<p>
<p>
-
 
+
This model was made using the Python language. You can <a href="https://static.igem.org/mediawiki/2013/a/ae/Duodenummodel.zip"> download the python script here</a>.
</p>
</p>
Line 177: Line 238:
   <li>Calculated from : Computational Modeling and Simulation of the Human Duodenum - B. Hari, S. Bakalis, P. Fryer - 2012</li>
   <li>Calculated from : Computational Modeling and Simulation of the Human Duodenum - B. Hari, S. Bakalis, P. Fryer - 2012</li>
   <li>Wikipedia</li>
   <li>Wikipedia</li>
-
   <li>Calculated from : http://onlinelibrary.wiley.com/doi/10.1016/S0168-6445%2803%2900055-X/pdf p.217 </li>
+
   <li>Calculated from : Bacterial iron homeostasis - Simon C. Andrews, Andrea K. Robinson, Francisco Rodriguez-Quinones</li>
  </ol>
  </ol>
</div>
</div>

Latest revision as of 03:36, 29 October 2013

Iron coli project

Final flush treatment model

Introduction

This model comes right after the Disease model. We now want to model the flush treatment by simulating a flush of iron-chelating bacteria.

Observations

Once our genetically modified bacteria are released in the duodenum lumen, they produce siderophores to chelate the solved iron, thus making it unavailable for intestinal absorption. Then, they eventually flush out of the duodenum. The main hypothesis in this model is that the bacteria do not colonize the duodenum : they only flow through. They do not even have time to grow, for the time required to flush through is close to 40 seconds.
Because some of the parameters remain unknown, it is a theoretical simulation which will not give any numerical results, only qualitative results.

Goal

This model was made to check if our first strategy was viable. It aims to answer the following question:
"Is it possible to chelate a significant amount of iron with a flush strategy?"

Assumptions

The same assumptions as in the previous model apply:

  • Our bacteria don't settle in the duodenum
  • No regulation in the patient's iron absorption
  • Constant iron flow in the duodenum lumen
  • Homogeneous fluid
  • The bacterial quantity is constant
  • The bacterial natural absorption is insignificant compared to the chelation
  • The patient ingests 20mg of iron per day (Guideline Daily Amounts)

Materials and methods

A : Total quantity of iron absorbed by the duodenum (mol)
S : Quantity of solved iron (mol)
P : Total quantity of enterobactin produced by our population of bacteria (mol)
Q : Total quantity of chelated iron (mol)
N : Number of bacteria



Sp mol.s-1 Iron pulse
v m.s-1 Chyme's flow average speed
L m Duodenum length
α s-1 Duodenum absorption rate
K mol/s Activator Magnitude
p mol.s-1 Value at zero of the activator
h - Activator efficiency
d mol Activator threshold
δ mol-1 Chelation rate


The graph on the right explains the reasoning: for instance, the arrow with a + between N and P means that the variation of P has a positive linear term in N.

Where is our Logistic function (which we will here abusively call activator).

Results

Absorption
Figure 1 : Iron Absorption in the duodenum with and without treatment.

The Figure 1 represents the total iron absorbed by the duodenum during an average meal (with and without bacterial treatment), whereas the Figure 2 represents the instant ambient iron in the duodenum lumen(with and without bacterial treatment).


Quantity
Figure 2 : Iron in the duodenum during the simulation. The grey strip represent the presence of the bacteria in the duodenum.

The Figure 2 points out the influence of our bacteria on the middle. The blue curve decreases significantly under the bacterial flow. That underlines the fact that the treatment doesn't directly affect the absorption. After 150 seconds, the chyme is considered to progressively leave the duodenum (entering the jejunum), thus reducing the iron absorption in both cases, which finally leads to a steady state.



Parameters:

The simulations was made with following parameters:

Name Value Unit Description Reference
Sp 4.5*10-9 mol.s-1 Iron pulse [1][2]
v 0.007 m.s-1 Chyme's flow average speed [2]
L 0.3 m Duodenum length [3]
α 0.03 s-1 Duodenum absorption rate Disease model
K 100 mol/s Activator Magnitude set
p 0.005 mol.s-1 Value at zero of the activator set
h 10-5 - Activator efficiency set
d 10-9 mol Activator threshold [4]
δ 2.65*10-8 mol-1 Chelation rate [4]


Conclusion

The combination of the two curves clearly shows that bacteria would be able to neutralize a significant quantity of iron before it is absorbed by the duodenum.
It would thus be possible to approximately divide the intestinal iron intake by 2 if the patient takes one pill during each meal. This means that the patient would endure a lighter treatment : less bloodletting for people suffering from hemochromatosis, and less iron chelator's side effects for the thalassemia.

Models and scripts

This model was made using the Python language. You can download the python script here.

References:

  1. Physiol Rev 93: 1721–1741, 2013 doi:10.1152/physrev.00008.2013 - Tomas Ganz "SYSTEMIC IRON HOMEOSTASIS"
  2. Calculated from : Computational Modeling and Simulation of the Human Duodenum - B. Hari, S. Bakalis, P. Fryer - 2012
  3. Wikipedia
  4. Calculated from : Bacterial iron homeostasis - Simon C. Andrews, Andrea K. Robinson, Francisco Rodriguez-Quinones