Team:UPO-Sevilla/Modeling/Signaling
From 2010.igem.org
(7 intermediate revisions not shown) | |||
Line 16: | Line 16: | ||
<p> | <p> | ||
- | The signaling circuit 3 described in the <a href="https://2010.igem.org/Team:UPO-Sevilla/Biobricks/Circuits">Circuit Section</a> has been modeled using Matlab Simbiology desktop. The following diagram shows the different parts of the model we have simulated: | + | The signaling circuit 3 described in the <a href="https://2010.igem.org/Team:UPO-Sevilla/Biobricks/Circuits">Circuit Section</a> has been modeled using <a href="http://www.mathworks.com/products/simbiology/" target="_blank">Matlab Simbiology</a> desktop. The following diagram shows the different parts of the model we have simulated: |
</p> | </p> | ||
+ | |||
<center> | <center> | ||
<a href="https://2010.igem.org/Image:UPOModelv2.png"> | <a href="https://2010.igem.org/Image:UPOModelv2.png"> | ||
Line 65: | Line 66: | ||
<h1>Simulations</h1> | <h1>Simulations</h1> | ||
- | <p>The following figure shows the typical evolution of the output of the system (the generated chemoattractant) againts the inputs (the wall cells ligand and the FecA-PrhA components on the outer membrane)</p> | + | <p>The following figure shows the typical evolution of the output of the system (the generated chemoattractant <i>medium.L_aspartate</i>) againts the inputs (the wall cells ligand and the FecA-PrhA components on the outer membrane)</p> |
Line 98: | Line 99: | ||
<p> | <p> | ||
- | This analysis is referred to the steady-state of the system. Some parameters do not affect the final steady-state number of molecules, but on the other hand | + | This analysis is referred to the steady-state of the system. Some parameters do not affect the final steady-state number of molecules, but on the other hand they affect the velocity of the system in the transitory. This can be seen in the following paragraphs. |
</p> | </p> | ||
Line 104: | Line 105: | ||
<p> | <p> | ||
- | If we perform a scan over several values of some parameters it can be seen the influence of these parameters on the output. For instance, in the following figure it can be seen how the parameter <i>kCellBinding</i> (the "force" of the binding with the cell wall) affects the final output of the system (<i>medium.L_aspartate</i>, the amount of chemoattractant) | + | If we perform a scan over several values of some parameters, it can be seen the influence of these parameters on the output. For instance, in the following figure it can be seen how the parameter <i>kCellBinding</i> (the "force" of the binding with the cell wall) affects the final output of the system (<i>medium.L_aspartate</i>, the amount of chemoattractant). The parameter has been scanned through several orders of magnitude. |
</p> | </p> | ||
Line 114: | Line 115: | ||
<p> | <p> | ||
- | + | It can be seen that this parameter affects the velocity of the system, but not in great deal. In this case, the parameter mainly influences the speed of the signaling circuit on the cell wall. But once this signal is activated, the rest of the system (in this case, the final production of aspartate) is not affected very much. | |
+ | </p> | ||
+ | |||
+ | <p> | ||
+ | In the following, the same scan of parameters is done for the four main steps in the generation of the chemoattractant: | ||
</p> | </p> | ||
Line 131: | Line 136: | ||
</center> | </center> | ||
+ | <p> | ||
+ | Again, it can be seen that the evolution of the output is quite unsensitive to the changes of parameters. In order to illustrate this, a simple mathematical analysis is performed in the next paragraphs. | ||
+ | </p> | ||
+ | |||
+ | <h2></h2> | ||
+ | |||
+ | <p> | ||
+ | The system can be seen as a set of smaller subsystems (listed above) that are coupled by the enzymatic promoters. This subsystems are qualitatively quite similar. They are usually modeled as mass action driven kynetic equations. | ||
+ | </p> | ||
+ | |||
+ | <p> | ||
+ | If we analyze the FecR activation reaction, induced by the binding between the FecA-PrhA and the cell wall ligand: | ||
+ | </p> | ||
+ | |||
+ | |||
+ | <center> | ||
+ | ecoli.FecR + ecoli.[ligand:FecA-PrhA] <-> ecoli.FecR_a + ecoli.[ligand:FecA-PrhA] | ||
+ | </center> | ||
+ | |||
+ | <p> | ||
+ | The differential equation governing the reaction is: | ||
+ | </p> | ||
+ | |||
+ | <center> | ||
+ | <img src="https://static.igem.org/mediawiki/2010/e/e9/UPOModel-React1.png" width="500" alt="Simbiology model"/> | ||
+ | </center> | ||
+ | <br> | ||
+ | <p> | ||
+ | were we have changed the names of some quantities for the sake of clarity: | ||
+ | </p> | ||
+ | |||
+ | <table xmlns="" id="modelContentsSpeciesTable" width="100%" border="5" cellpadding="5" cellspacing="0" class="dataTable"><tbody> | ||
+ | <tr><th>Name</th><th>Variable Name</th></tr> | ||
+ | |||
+ | <tr class="oddRow"><td>ligand:FecA-PrhA</td><td>bind</td></tr> | ||
+ | <tr class="oddRow"><td>kFecRActivation</td><td>k<sub>9</sub></td></tr> | ||
+ | <tr class="oddRow"><td>kFecRDeactivation</td><td>k<sub>10</sub></td></tr> | ||
+ | </table> | ||
+ | |||
+ | |||
+ | <p> | ||
+ | Knowing that <i>FecR</i>+<i>FecRa</i>=<i>FecR</i><sub>0</sub> (<i>FecR</i><sub>0</sub> being the initial amount of <i>FecR</i>), and some simple calculations, then: | ||
+ | </p> | ||
+ | |||
+ | |||
+ | |||
+ | <center> | ||
+ | <img src="https://static.igem.org/mediawiki/2010/a/a0/UPOModel-React2.png" width="500" alt="Simbiology model"/> | ||
+ | </center> | ||
+ | <br> | ||
+ | |||
+ | If we assume that the inducer <i>bind</i> is constant, then the equation is just a first order differential equation, and the evolution of the quantity of <i>FecRa</i> is given by: | ||
+ | |||
+ | <center> | ||
+ | <img src="https://static.igem.org/mediawiki/2010/8/8c/UPOModel-React3.png" width="500" alt="Simbiology model"/> | ||
+ | </center> | ||
+ | <br> | ||
+ | |||
+ | <p> | ||
+ | The steady-state solution of this equation is k<sub>9</sub><i>FecR</i><sub>0</sub>/(k<sub>9</sub>+k<sub>10</sub>), and it does not depend on <i>bind</i>. | ||
+ | </p> | ||
+ | |||
+ | <p> | ||
+ | The <b>time constant</b> of the system (its speed) is governed by (k<sub>9</sub>+k<sub>10</sub>) and <i>bind</i>. | ||
+ | </p> | ||
+ | |||
+ | <p> | ||
+ | In the actual system, <i>bind</i> evolves with time. If its evolution is much quicker than this system, reaching its steady-state, it can be considered constant and the same analysis can be applied. It acts as an activator of the system; moreover, the higher the quantity, the quicker the system. | ||
+ | </p> | ||
+ | |||
+ | <p> | ||
+ | If <i>bind</i> evolves with time the solution of the equation is bit more involved. A simplified analysis can be done assuming that <i>bind</i> is timewise constant for small amounts of time, and applying the same idea. The main effect of <i>bind</i> is to accelerate or decelerate the evolution. | ||
+ | </p> | ||
+ | |||
+ | |||
+ | <p> | ||
+ | The simulation of the FecR system, with no approximation, is very similar to the first order system described: | ||
+ | </p> | ||
+ | |||
+ | |||
+ | <center> | ||
+ | <a href="https://2010.igem.org/Image:UPOModel-FecR.png"> | ||
+ | <img src="https://static.igem.org/mediawiki/2010/1/1f/UPOModel-FecR.png" width="800" alt="Simbiology model"/> | ||
+ | </a> | ||
+ | </center> | ||
Latest revision as of 15:29, 27 October 2010
The signaling circuit
The signaling circuit 3 described in the Circuit Section has been modeled using Matlab Simbiology desktop. The following diagram shows the different parts of the model we have simulated:
![Simbiology model](https://static.igem.org/mediawiki/2010/5/5b/UPOModelv2.png)
For the level of detail considered, the main parts simulated are the following (the number correspond to the equations listed in the table in the next section):
- Generation of L_aspartate induced by AAL
- Diffusion of L_aspartate through the cell wall
- Transcription of the aspA, promoted by FecI_a (active)
- Translation of aspA
- Activation of FecI, induced by the activation of FecR
- Activation of FecR induced by FecA-PrhA
- Plant cell wall lingand, FecA-PrhA binding
Reactions
The reaction equations for the previous parts, and the reactions rates associated, are summarized in the following table:
# | Reaction | ReactionRate | Active |
---|---|---|---|
1 | ecoli.ammonia + ecoli.fumarate + ecoli.AAL <-> ecoli.L_aspartate + ecoli.AAL | k1*ecoli.ammonia*ecoli.fumarate*ecoli.AAL - k2*ecoli.L_aspartate*ecoli.AAL | true |
2 | ecoli.L_aspartate <-> medium.L_aspartate | kWallDiffusion*ecoli.L_aspartate - kWallDiffusionBack*medium.L_aspartate | true |
3 | ecoli.DNAaspA + ecoli.FecI_a -> ecoli.ARNm_aspA + ecoli.DNAaspA + ecoli.FecI_a | kTranscript*ecoli.DNAaspA*ecoli.FecI_a | true |
4 | ecoli.ARNm_aspA -> ecoli.AAL + ecoli.ARNm_aspA | kTranslation*ecoli.ARNm_aspA | true |
5 | ecoli.FecR_a + ecoli.FecI <-> ecoli.FecI_a + ecoli.FecR_a | kFecIActivation*ecoli.FecR_a*ecoli.FecI - kFecIDeactivation*ecoli.FecI_a*ecoli.FecR_a | true |
6 | ecoli.FecR + ecoli.[ligand:FecA-PrhA] <-> ecoli.FecR_a + ecoli.[ligand:FecA-PrhA] | kFecRActivation*ecoli.FecR*ecoli.[ligand:FecA-PrhA] - kFecRDeactivation*ecoli.FecR_a*ecoli.[ligand:FecA-PrhA] | true |
7 | plant_cell_wall.ligand + ecoli.[FecA-PrhA] <-> ecoli.[ligand:FecA-PrhA] | kCellBinding*plant_cell_wall.ligand*ecoli.[FecA-PrhA] - kCellUnbinding*ecoli.[ligand:FecA-PrhA] | true |
Simulations
The following figure shows the typical evolution of the output of the system (the generated chemoattractant medium.L_aspartate) againts the inputs (the wall cells ligand and the FecA-PrhA components on the outer membrane)
![Simbiology model](https://static.igem.org/mediawiki/2010/9/9f/UPOSimulation.png)
Analysis
Sensibility
Simbiology allows to compute the sensibility of the system against the different parameters.
The following figure shows the sensibility of all state variables (molecules of the different species considered) with respect to all the parameters.
![Simbiology model](https://static.igem.org/mediawiki/2010/2/21/UPOSensibilityAll.png)
What the analysis reveal is that the system is quite insensitive to changes in the parameters. This is due mainly to the nature of the transduction signals. The promoters act as a kind of "switch". This means that, provided these promoters reach certain levels, the other parts of the circuits are activated, even if the levels are not equal.
This analysis is referred to the steady-state of the system. Some parameters do not affect the final steady-state number of molecules, but on the other hand they affect the velocity of the system in the transitory. This can be seen in the following paragraphs.
Scanning of Parameters
If we perform a scan over several values of some parameters, it can be seen the influence of these parameters on the output. For instance, in the following figure it can be seen how the parameter kCellBinding (the "force" of the binding with the cell wall) affects the final output of the system (medium.L_aspartate, the amount of chemoattractant). The parameter has been scanned through several orders of magnitude.
![Simbiology model](https://static.igem.org/mediawiki/2010/8/84/UPOScanBinding.png)
It can be seen that this parameter affects the velocity of the system, but not in great deal. In this case, the parameter mainly influences the speed of the signaling circuit on the cell wall. But once this signal is activated, the rest of the system (in this case, the final production of aspartate) is not affected very much.
In the following, the same scan of parameters is done for the four main steps in the generation of the chemoattractant:
- Transcription of the aspA, promoted by FecI_a (active)
- Translation of aspA
- Activation of FecI, induced by the activation of FecR
- Activation of FecR induced by FecA-PrhA
![Simbiology model](https://static.igem.org/mediawiki/2010/1/19/UPOScanBindFecIFecRTranscriptTranslation.png)
Again, it can be seen that the evolution of the output is quite unsensitive to the changes of parameters. In order to illustrate this, a simple mathematical analysis is performed in the next paragraphs.
The system can be seen as a set of smaller subsystems (listed above) that are coupled by the enzymatic promoters. This subsystems are qualitatively quite similar. They are usually modeled as mass action driven kynetic equations.
If we analyze the FecR activation reaction, induced by the binding between the FecA-PrhA and the cell wall ligand:
The differential equation governing the reaction is:
![Simbiology model](https://static.igem.org/mediawiki/2010/e/e9/UPOModel-React1.png)
were we have changed the names of some quantities for the sake of clarity:
Name | Variable Name |
---|---|
ligand:FecA-PrhA | bind |
kFecRActivation | k9 |
kFecRDeactivation | k10 |
Knowing that FecR+FecRa=FecR0 (FecR0 being the initial amount of FecR), and some simple calculations, then:
![Simbiology model](https://static.igem.org/mediawiki/2010/a/a0/UPOModel-React2.png)
If we assume that the inducer bind is constant, then the equation is just a first order differential equation, and the evolution of the quantity of FecRa is given by:
![Simbiology model](https://static.igem.org/mediawiki/2010/8/8c/UPOModel-React3.png)
The steady-state solution of this equation is k9FecR0/(k9+k10), and it does not depend on bind.
The time constant of the system (its speed) is governed by (k9+k10) and bind.
In the actual system, bind evolves with time. If its evolution is much quicker than this system, reaching its steady-state, it can be considered constant and the same analysis can be applied. It acts as an activator of the system; moreover, the higher the quantity, the quicker the system.
If bind evolves with time the solution of the equation is bit more involved. A simplified analysis can be done assuming that bind is timewise constant for small amounts of time, and applying the same idea. The main effect of bind is to accelerate or decelerate the evolution.
The simulation of the FecR system, with no approximation, is very similar to the first order system described:
![Simbiology model](https://static.igem.org/mediawiki/2010/1/1f/UPOModel-FecR.png)