Team:Evry/ODE model

From 2012.igem.org

(Difference between revisions)
 
(69 intermediate revisions not shown)
Line 2: Line 2:
<html>
<html>
-
<h1>Model using Ordinary Differential Equations(ODE)</h1>
+
<h1>Model using Ordinary Differential Equations (ODE)</h1>
<h2>Overview</h2>
<h2>Overview</h2>
-
 
+
<ul>
-
The first model we developed represents the tadpole as a three compartment system:
+
This global model represents the tadpole as a three compartment system:
<ol>
<ol>
-
<li>The skin that produces (or receives) auxins;</li>
+
<li>The skin that produces (or receives) auxins, <strong>denoted S</strong> in the equations;</li>
-
<li>The blood that transport auxins to the organs;</li>
+
<li>The blood that transport auxins to the organs, <strong>denoted B</strong>;</li>
-
<li>The organs (called receptors) that interacts with auxin molecules.</li>
+
<li>The organs (called receptors) that interacts with auxin molecules, <strong>denoted R</strong>.</li>
</ol>
</ol>
<br/>
<br/>
-
The parallel with (electrical) engineering is made easy: the skin represents a generator that will add a quantity to the system;
 
-
The blood represents wires, that convey this quantity throughout the system;
 
-
Finally the organs are the sinks that use the quantity to work.<br/><br/>
 
-
 
-
This very idealized view of the tadpoles allows to make some interesting simplifications: The processes happening in the system can be approximated using Ordinary Differential Equations (ODE), one of the simplest form of differential equations; Plus, the organs repartition and shape are not taken into account.<br/><br/>
 
-
This over-simplication of the problem causes the model to give very imprecise quantitative results but its strength is in allowing us to make some qualitative predictions about the success or failure of some experiments.
+
Using this very idealized view of the tadpoles makes it possible to model it simply : The processes happening in the system can be approximated using Ordinary Differential Equations (ODE), one of the simplest form of differential equations; Plus, the organs repartition and shape are not taken into account. As showed in <a href="https://2012.igem.org/Team:Evry/auxin_pde">this page</a>, these equations can be seen as approximations of more realistic PDEs where space considerations have been removed.<br/><br/>
 +
Although our model is very general and simple, it allows to have a quantitative view. As many parameters would require much experiment to be estimated correctly. At this stage this quantitative model only allows us to make some prediction about the global behaviour of the system and gives us hints about the reasons for the success or failure of experiments.
 +
</p>
<br/>
<br/>
<center>
<center>
Line 27: Line 24:
</center>
</center>
-
<h2>Hypothesis</h2>
+
<h2>Assumptions</h2>
-
There are the different hypothesis we were constrained to make in order to model the system:
+
There are the different assumptions we were constrained to make in order to model the system:
<ol>
<ol>
Line 37: Line 34:
     <li><a name="hyp2"></a>
     <li><a name="hyp2"></a>
-
<strong>No auxins can go from the skin directly to the organs.</strong></li>
+
<strong>All the auxins that enter the receptor come from blood.</strong></li>
-
We have chosen to neglect the exchanges between the skin and the other organs.
+
Considering blood speed and the big density of capillaries, we assumed that the auxin coming from other tissues can be neglected. This assumption will remain while the receptor isn't close to skin.
     <li><a name="hyp3"></a>
     <li><a name="hyp3"></a>
<strong>The auxin flow follows the concentration gradient between compartments.</strong></li>
<strong>The auxin flow follows the concentration gradient between compartments.</strong></li>
-
This hypothesis is based on the fact that a small capillary is nothing more than a hole in a cell layer. So the exchanges between the surrounding cells and the capillary are mutual.
+
This assumption is based on the fact that a small capillary is nothing more than a hole in a cell layer. So the exchanges between the surrounding cells and the capillary are mutual.
 +
  <li><a name="hyp4"></a>
 +
<strong>Permeabilities are symetric.</strong></li>
 +
We have considered that for 2 issues X and Y the permeability from X to Y is equal to the one from Y to X.
 +
 
</ol>
</ol>
Line 48: Line 49:
<h3>Equations</h3>  
<h3>Equations</h3>  
-
Each compartment is modeled by a differential equation representing the evolution of the auxin quantity as a function of the time.
+
<p>Each compartment is modeled by a differential equation representing the evolution of the auxin quantity as a function of the time.
-
Each equations are composed of two kinds of terms: creation and degradation. The creation term can represent either a creation of auxin in the compartment or an arriving quantity of auxin from another one. In the same way, the degradation term can either represent a natural degradation of molecules or a quantity leaving the compartment.<br/><br/>
+
Each equation is composed of two kinds of terms: synthesis or arrival (called <i>born</i>) and degradation (called <i>die</i>). The creation term can represent either a creation of auxin in the compartment or an arriving quantity of auxin from another one. The following equations represent the case where auxin flows in from water. In the case where we consider auxin being synthesized directly in skin, this term would be the mean rate of synthesis. This rate being computed integrating data from the <a href="https://2012.igem.org/Team:Evry/plasmid_splitting">plasmid repartition model</a> and the <a href="https://2012.igem.org/Team:Evry/auxin_production">auxin production model</a>, as described in <a href="https://2012.igem.org/Team:Evry/auxin_pde">our derivation from PDEs</a>. In the same way, the degradation term can either represent a natural degradation of molecules or a quantity leaving the compartment.<br/><br/>
 +
</p>
<center>
<center>
Line 57: Line 59:
<br/>
<br/>
-
In this system, the J terms represent the fluxes between the different compartments. We made them depend on the concentrations of both the in and out compartments as explained in <a href="#hyp2">hypothesis 2</a>. Their mathematical formulation is the following:<br/><br/>
+
<p>
 +
In this system,D<sub>die</sub> is a degradation rate,  the other D terms represent the flows between the different compartments. We made them depend on the concentrations of both the in and out compartments as explained in <a href="#hyp2">assumption 2</a>. Their mathematical formulation is the following:<br/><br/>
 +
</p>
<center>
<center>
-
<img src="https://static.igem.org/mediawiki/2012/7/7e/Fluxes.png" alt="mathematical expression of fluxes" />
+
<img src="https://static.igem.org/mediawiki/2012/7/7e/Fluxes.png" alt="mathematical expression of debits" />
</center>
</center>
<br/>
<br/>
Line 73: Line 77:
</ul>
</ul>
<br/>
<br/>
-
 
+
<p>
These flow equations are based on Newton's law of cooling where the difference between the concentrations of the two compartments gives the direction and magnitude of the flow. This allows us to model in a single equations the two opposite flows between the compartments.
These flow equations are based on Newton's law of cooling where the difference between the concentrations of the two compartments gives the direction and magnitude of the flow. This allows us to model in a single equations the two opposite flows between the compartments.
 +
</p>
 +
<br/>
<h3>Parameters</h3>
<h3>Parameters</h3>
-
The parameters of the model are hard to estimate because they are an aggregation of different physical values. This is due to the high level view provided by the model. The different parameters are:
+
<br/>
 +
<p>
 +
The parameters of the model are hard to estimate because direct experiences concerning auxin's behaviour in the xenopus where never made before, as it is a plant hormone.Moreover,the parameters are an aggregation of different physical values. This is due to the high level view provided by the model. The different parameters are:
 +
</p>
 +
<br/>
 +
<br/>
<center>
<center>
Line 87: Line 98:
     <th>Type</th>
     <th>Type</th>
     <th>Values</th>
     <th>Values</th>
 +
    <th>Experiment Proposal</th>
   </tr>
   </tr>
Line 92: Line 104:
     <td>Permeabilities</td>
     <td>Permeabilities</td>
     <td>P</td>
     <td>P</td>
-
     <td>calculated</td>
+
     <td>similar molecule in litterature</td>
-
     <td><a href="https://2012.igem.org/Team:Evry/parameters">here</a></td>
+
     <td><a href="https://2012.igem.org/Team:Evry/parameters#permeabilities">here</a></td>
 +
    <td><a href="https://2012.igem.org/Team:Evry/Experimentals_Parameters#2">link</a></td>
   </tr>
   </tr>
Line 100: Line 113:
     <td>S</td>
     <td>S</td>
     <td>calculated</td>
     <td>calculated</td>
-
     <td><a href="https://2012.igem.org/Team:Evry/parameters">here</a></td>
+
     <td><a href="https://2012.igem.org/Team:Evry/parameters#surfaces">here</a></td>
 +
    <td> - </td>
 +
   
   </tr>
   </tr>
Line 107: Line 122:
     <td>V</td>
     <td>V</td>
     <td>calculated</td>
     <td>calculated</td>
-
     <td><a href="https://2012.igem.org/Team:Evry/parameters">here</a></td>
+
     <td><a href="https://2012.igem.org/Team:Evry/parameters#volumes">here</a></td>
 +
    <td> - </td>
   </tr>
   </tr>
   <tr>
   <tr>
     <td>Degradation rate</td>
     <td>Degradation rate</td>
-
     <td>D<sub>Die</sub>
+
     <td>D<sub>die</sub>
-
     <td>estimated</td>
+
     <td>unknown</td>
-
     <td>?</td>
+
     <td>unknown</td>
 +
    <td><a href="https://2012.igem.org/Team:Evry/Experimentals_Parameters#4">link</a></td>
   </tr>
   </tr>
   <tr>
   <tr>
     <td>Creation rate</td>
     <td>Creation rate</td>
-
     <td>D<sub>Born</td>
+
     <td>D<sub>born</td>
     <td>computed</td>
     <td>computed</td>
-
     <td><a href="https://2012.igem.org/Team:Evry/plasmid_splitting">plasmid repartition model</a></td>
+
     <td><a href="https://2012.igem.org/Team:Evry/auxin_pde#Dborn">our derivation from PDEs</a></td>
 +
    <td> - </td>
   </tr>
   </tr>
</table>
</table>
</center>
</center>
 +
<h2> Sensitivity Analysis</h2>
 +
<p>There are 2 kinds of parameters in this model; those concerning the geometry of the tadpole and those concerning the behaviour of auxin in the tadpole's tissues. As auxin is a plant hormone, we had to evaluate these parameters by taking those of similar molecules. That's why the uncertainties on these measures are very big. For this category, we observed variations from reference value (ref) from ref*10<sup>-3</sup> to ref*10<sup>3</sup>.
 +
You can see the effect of these variations:
 +
<p>
 +
<ul>
 +
  <li>For permeabilities:
 +
<p>Playing with skin permeability influences highly the auxin quantity in each compartment. This is expected since the skin is the emitter in our system. As for receptor permeability, the higher it is the less we get auxin in skin, leeding us to a conservative behaviour in auxin quantity between skin and receptor compartments. Tuning the blood permeability seems to be less significative for auxin quantity in the compartments giving that the blood plays a transport role.<p>
 +
  <center>
 +
  <table>
 +
    <tr>
 +
      <td><img src="https://static.igem.org/mediawiki/2012/2/2d/P1.gif" width="450" alt="" /></td>
 +
      <td><img src="https://static.igem.org/mediawiki/2012/0/02/P2.gif" width="450" alt="" /></td>
 +
    </tr>
 +
    <tr>
 +
      <td><img src="https://static.igem.org/mediawiki/2012/2/25/P3.gif" width="450" alt="" /></td>
 +
    </tr>
 +
  </table>
 +
  </center>
 +
 +
  </li>
 +
  <li>For die rate:
 +
<p>Tuning the degradation rate has leed us to a signficant value 0.00001 mol after which, we get a different behaviour of the auxin quantity evolution in receptor. As it represents the destination compartment of auxin emitted, it may explain why it was the most affected.<p>
 +
 +
  <center>
 +
  <table>
 +
    <tr>
 +
    <td><img src="https://static.igem.org/mediawiki/2012/archive/f/f7/20121026235650!D.gif" width="450" alt="" /></td>
 +
    </tr>
 +
  </table>
 +
  </center>
 +
 +
  </li>
 +
</ul>
 +
 +
For the geometrical parameters, we observed variations going from ref*0.3 to ref*2.We've obtained:
 +
<ul>
 +
  <li>For Surfaces:
 +
<p> Simulations with different surfaces values do not really show big differences in the main landscape which keep the emitter with the higher quantity of auxin. Tuning the skin surface gives a proportional evolution of auxin quantity in receptor and in emitter, in the other hand tuning the receptor surface produces an antiproportional evolution of auxin quantity in the two compartments. <p>
 +
 +
  <center>
 +
  <table>
 +
    <tr>
 +
      <td><img src="https://static.igem.org/mediawiki/2012/6/6e/S1.gif" width="450" alt="" /></td>
 +
      <td><img src="https://static.igem.org/mediawiki/2012/2/26/S2.gif" width="450" alt="" /></td>
 +
    </tr>
 +
    <tr>
 +
      <td><img src="https://static.igem.org/mediawiki/2012/0/04/S3.gif" width="450" alt="" /></td>
 +
    </tr>
 +
  </table>
 +
  </center>
 +
 +
  </li>
 +
  <li>For Volumes:
 +
<p> The auxin quantity in each compartment seems to be sensitive to the volume of the latter, except for the receptor. The system is not sensitive to small variations of the receptor volume.<p>     
 +
 +
  <center>
 +
  <table>
 +
    <tr>
 +
      <td><img src="https://static.igem.org/mediawiki/2012/a/a1/V1.gif" width="450" alt="" /></td>
 +
      <td><img src="https://static.igem.org/mediawiki/2012/6/62/V2.gif" width="450" alt="" /></td>
 +
    </tr>
 +
    <tr>
 +
      <td><img src="https://static.igem.org/mediawiki/2012/a/af/V3.gif" width="450" alt="" /></td>
 +
    </tr>
 +
  </table>
 +
  </center>
 +
 +
  </li>
 +
</ul>
-
<h2>Calibration</h2>
 
<h2>Results</h2>
<h2>Results</h2>
 +
By modifying the initial conditions and c<sub>w</sub> we have observed globaly 3 different behaviours of the system:
 +
<ul>
 +
  <li>
 +
    <p>This first graph is obtained considering no auxin in the tadpole and c<sub>w</sub>=5.10<sup>-3</sup>mol.L<sup>-1</sup> which is the concentration we wanted to use for experiences. The result shows that the quantity of auxin in skin and in the receptor grow until a limit value.
 +
    </p>
 +
    <img src="https://static.igem.org/mediawiki/2012/4/41/Refodegraph.png" width="450px" >
 +
  </li>
 +
    <li><p>If there has been injected auxin in the receptor and in skin: n<sub>s</sub>(t=0)=10 mmol n<sub>r</sub>(t=0)=10 mmol, we can observe a decrease of auxin in receptor before an establishment of the equilibrium.  </p>
 +
<img src="https://static.igem.org/mediawiki/2012/1/19/Init10inskinandreceptor.png" width="450px">
 +
  </li>
 +
  <li><p> Here , the concentration in water has been devided by ten c<sub>w</sub>=0.5.10<sup>-3</sup>mol.L<sup>-1</sup>, while the initial value of skin is still n<sub>s</sub>(t=0)=10 mmol. We observe first a big increase in the receptor due to the diffusion of the auxin intially injected in skin then a decrease because of the small water concentration in auxin.
 +
      </p>
 +
  <img src="https://static.igem.org/mediawiki/2012/4/41/OdegraphInit10Cwater500.png" width="450px" >
 +
  </li>
 +
 +
<h2>Conclusion</h2>
<h2>Conclusion</h2>
 +
<p>This model being simplified to the maximum, has let us, thanks to an estimation of very few parameters, to see the global behaviour of the system. It helps us to determine expected behaviours and helps us to establish possible reasons of dysfunctions. Nevertheless, quantitatively, this model is too imprecise because it doesn't take tadpole's geomoetry into account. This is the aim of our most elaborated model which will use partial equations.
 +
</p>
 +
 +
You can find the model <a href="https://2012.igem.org/File:General_ode.zip"> here</a>.
<h2>References</h2>
<h2>References</h2>
Line 137: Line 243:
<p id="references">References:</p>
<p id="references">References:</p>
<ol>
<ol>
-
<li></li>
+
<li><a href="http://www.uni-leipzig.de/diffusion/journal/pdf/volume2/diff_fund_2%282005%291.pdf"> Jean Philibert, One and a Half Century of Diffusion: Fick, Einstein, before and beyond, Diffusion Fundamentals 2, 2005</a></li>
</ol>
</ol>
</div>
</div>

Latest revision as of 03:26, 27 October 2012

Model using Ordinary Differential Equations (ODE)

Overview

    This global model represents the tadpole as a three compartment system:
    1. The skin that produces (or receives) auxins, denoted S in the equations;
    2. The blood that transport auxins to the organs, denoted B;
    3. The organs (called receptors) that interacts with auxin molecules, denoted R.

    Using this very idealized view of the tadpoles makes it possible to model it simply : The processes happening in the system can be approximated using Ordinary Differential Equations (ODE), one of the simplest form of differential equations; Plus, the organs repartition and shape are not taken into account. As showed in this page, these equations can be seen as approximations of more realistic PDEs where space considerations have been removed.

    Although our model is very general and simple, it allows to have a quantitative view. As many parameters would require much experiment to be estimated correctly. At this stage this quantitative model only allows us to make some prediction about the global behaviour of the system and gives us hints about the reasons for the success or failure of experiments.


    tadpole compartments

    Assumptions

    There are the different assumptions we were constrained to make in order to model the system:
    1. The auxin concentration inside a compartment is homogeneous
    2. This condition is inherent to this kind of modeling.
    3. All the auxins that enter the receptor come from blood.
    4. Considering blood speed and the big density of capillaries, we assumed that the auxin coming from other tissues can be neglected. This assumption will remain while the receptor isn't close to skin.
    5. The auxin flow follows the concentration gradient between compartments.
    6. This assumption is based on the fact that a small capillary is nothing more than a hole in a cell layer. So the exchanges between the surrounding cells and the capillary are mutual.
    7. Permeabilities are symetric.
    8. We have considered that for 2 issues X and Y the permeability from X to Y is equal to the one from Y to X.

    Model description

    Equations

    Each compartment is modeled by a differential equation representing the evolution of the auxin quantity as a function of the time. Each equation is composed of two kinds of terms: synthesis or arrival (called born) and degradation (called die). The creation term can represent either a creation of auxin in the compartment or an arriving quantity of auxin from another one. The following equations represent the case where auxin flows in from water. In the case where we consider auxin being synthesized directly in skin, this term would be the mean rate of synthesis. This rate being computed integrating data from the plasmid repartition model and the auxin production model, as described in our derivation from PDEs. In the same way, the degradation term can either represent a natural degradation of molecules or a quantity leaving the compartment.

    ODE system

    In this system,Ddie is a degradation rate, the other D terms represent the flows between the different compartments. We made them depend on the concentrations of both the in and out compartments as explained in assumption 2. Their mathematical formulation is the following:

    mathematical expression of debits

    Where:
    • S in m2, represents the area of the exchange surface between the two compartments.
    • P in m2, represents the permeability of the membrane between the specified compartments.
    • C in [quantity] / m3, represents the concentration of auxin in the specified compartment

    These flow equations are based on Newton's law of cooling where the difference between the concentrations of the two compartments gives the direction and magnitude of the flow. This allows us to model in a single equations the two opposite flows between the compartments.


    Parameters


    The parameters of the model are hard to estimate because direct experiences concerning auxin's behaviour in the xenopus where never made before, as it is a plant hormone.Moreover,the parameters are an aggregation of different physical values. This is due to the high level view provided by the model. The different parameters are:



    Description Symbol Type Values Experiment Proposal
    Permeabilities P similar molecule in litterature here link
    Contact surfaces S calculated here -
    Volumes V calculated here -
    Degradation rate Ddie unknown unknown link
    Creation rate Dborn computed our derivation from PDEs -

    Sensitivity Analysis

    There are 2 kinds of parameters in this model; those concerning the geometry of the tadpole and those concerning the behaviour of auxin in the tadpole's tissues. As auxin is a plant hormone, we had to evaluate these parameters by taking those of similar molecules. That's why the uncertainties on these measures are very big. For this category, we observed variations from reference value (ref) from ref*10-3 to ref*103. You can see the effect of these variations:

    • For permeabilities:

      Playing with skin permeability influences highly the auxin quantity in each compartment. This is expected since the skin is the emitter in our system. As for receptor permeability, the higher it is the less we get auxin in skin, leeding us to a conservative behaviour in auxin quantity between skin and receptor compartments. Tuning the blood permeability seems to be less significative for auxin quantity in the compartments giving that the blood plays a transport role.

    • For die rate:

      Tuning the degradation rate has leed us to a signficant value 0.00001 mol after which, we get a different behaviour of the auxin quantity evolution in receptor. As it represents the destination compartment of auxin emitted, it may explain why it was the most affected.

    For the geometrical parameters, we observed variations going from ref*0.3 to ref*2.We've obtained:
    • For Surfaces:

      Simulations with different surfaces values do not really show big differences in the main landscape which keep the emitter with the higher quantity of auxin. Tuning the skin surface gives a proportional evolution of auxin quantity in receptor and in emitter, in the other hand tuning the receptor surface produces an antiproportional evolution of auxin quantity in the two compartments.

    • For Volumes:

      The auxin quantity in each compartment seems to be sensitive to the volume of the latter, except for the receptor. The system is not sensitive to small variations of the receptor volume.

    Results

    By modifying the initial conditions and cw we have observed globaly 3 different behaviours of the system:
    • This first graph is obtained considering no auxin in the tadpole and cw=5.10-3mol.L-1 which is the concentration we wanted to use for experiences. The result shows that the quantity of auxin in skin and in the receptor grow until a limit value.

    • If there has been injected auxin in the receptor and in skin: ns(t=0)=10 mmol nr(t=0)=10 mmol, we can observe a decrease of auxin in receptor before an establishment of the equilibrium.

    • Here , the concentration in water has been devided by ten cw=0.5.10-3mol.L-1, while the initial value of skin is still ns(t=0)=10 mmol. We observe first a big increase in the receptor due to the diffusion of the auxin intially injected in skin then a decrease because of the small water concentration in auxin.

    • Conclusion

      This model being simplified to the maximum, has let us, thanks to an estimation of very few parameters, to see the global behaviour of the system. It helps us to determine expected behaviours and helps us to establish possible reasons of dysfunctions. Nevertheless, quantitatively, this model is too imprecise because it doesn't take tadpole's geomoetry into account. This is the aim of our most elaborated model which will use partial equations.

      You can find the model here.

      References

      Other possible topologies

      With auxin in the external medium:
      tadpole + external compartments

      With a specific receptor organ:
      tadpole + other compartments