Team:Grenoble/Modeling/Signaling
From 2012.igem.org
(Difference between revisions)
(30 intermediate revisions not shown) | |||
Line 5: | Line 5: | ||
<div id="cadre"> | <div id="cadre"> | ||
<section> | <section> | ||
- | <h1> | + | <h1>Signaling module</h1> |
</br> | </br> | ||
</br> | </br> | ||
Line 14: | Line 14: | ||
</br> | </br> | ||
</br> | </br> | ||
- | Once the dipeptide molecule is fixed to the TAP receptor, it activates the phosphorylation of | + | Once the dipeptide molecule is fixed to the TAP receptor, it activates the phosphorylation of ompR. ompR* (the phosphorylated form of ompR) is the transcription factor that activates the gene expression of <i>cyaA</i>. |
</br> | </br> | ||
- | + | Further details about the receptor design, are available <a href="https://2012.igem.org/Team:Grenoble/Biology/Network#10">here</a>. | |
</br> | </br> | ||
</br> | </br> | ||
- | + | A computer model was developped to answer 2 questions: | |
</br> | </br> | ||
</br> | </br> | ||
- | 1- How sensitive is the signaling module? | + | 1- How sensitive is the signaling module? Wich dipeptide concentration should be present in order to trigger a response? |
</br> | </br> | ||
</br> | </br> | ||
Line 28: | Line 28: | ||
</br> | </br> | ||
</br> | </br> | ||
- | Answering these two questions will help us assess the sensitivity and the | + | Answering these two questions will help us assess the sensitivity and the response time of the whole system (ie the signaling module, the amplification module and the quorum sensing). |
</br> | </br> | ||
</br> | </br> | ||
- | + | First, we define the ordinary differential equations that govern adenylate cyclase (Ca) evolution. Thus we plot the evolution of Ca concentration versus initial dipeptide concentration (sensitivity) as well as the temporal evolution of Ca concentration (time response). | |
- | + | ||
</section> | </section> | ||
<section> | <section> | ||
- | <h1> | + | <h1>Ordinary Differential Equations</h1> |
- | Let’s begin by considering the cya gene activation by the transcription factor | + | Let’s begin by considering the cya gene activation by the transcription factor ompR*. |
</br> | </br> | ||
- | As it is a gene activation, the transcription rate is usually modelized by a | + | As it is a gene activation, the regulated transcription rate is usually modelized by a Hill function: |
</br> | </br> | ||
</br> | </br> | ||
Line 44: | Line 43: | ||
</br> | </br> | ||
</br> | </br> | ||
- | Where Vm is the maximal transcription rate, k is the activation coefficient, p is the basal production coefficient and | + | Where Vm is the maximal transcription rate, k is the activation coefficient, p is the basal production coefficient and α the degradation coefficient. |
</br> | </br> | ||
- | For | + | For ompR phosphorylation, we considered a model that takes into account the enzymatic mechanism of the Histidine Kinase EnvZ as well as the phosphotransfer and the phosphatase.<!--Référence à mettre ici !!!--> |
</br> | </br> | ||
- | The model we use is a phenomenological extension of the | + | The model we use is a phenomenological extension of the Goldbeter-Koshland biochemical switch model.<a href="#ref">[1]</a></br> |
- | The resulting equation governs [ | + | The resulting equation governs the temporal evolution of [ompR*] and highlights the fact that the phosphorylation is activated by a dipeptide. |
</br> | </br> | ||
</br> | </br> | ||
Line 57: | Line 56: | ||
where K and K’ are the dimensionless Michaelis-Menten coefficients. | where K and K’ are the dimensionless Michaelis-Menten coefficients. | ||
</br> | </br> | ||
- | Since the process involved in the production of the new protein | + | Since the process involved in the production of the new protein proceed at much slower timescale than the phosphorylation process that aims at chemically modifying the existing protein OmpR, the time derivative of [ompR*] is null (click here for more explanation) and [ompR] and [ompR*] are linked by a conservation law: |
</br> | </br> | ||
- | <center>[ | + | <center>[ompR*]+[ompR]=[ompR]<SUB>tot</SUB></center> |
</br> | </br> | ||
- | Once we have set the derivative equation of | + | Once we have set the derivative equation of ompR* equals to zero and replaced the value of [ompR] and [ompR*] by their expressions involving the total quantity of ompR, we get a second order polynomial equation of [ompR*]: |
</br> | </br> | ||
</br> | </br> | ||
Line 73: | Line 72: | ||
</br> | </br> | ||
</br> | </br> | ||
- | We notice that the coefficient “a” is always negative because the dephosphorylation rate of | + | We notice that the coefficient “a” is always negative because the dephosphorylation rate of ompR* is lower in value than the phosphorylation rate .As a product of kinetic parameters, the coefficient “c” is positive. |
</br> | </br> | ||
For the reasons given above the determinant is positive: | For the reasons given above the determinant is positive: | ||
Line 81: | Line 80: | ||
</br> | </br> | ||
</br> | </br> | ||
- | We have then two | + | We have then two roots x and y for [ompR*]. Their product is given by x.y = c/a. |
- | As “a” and “c” have opposite signs, x and y have opposite signs: We choose the positive root. | + | As “a” and “c” have opposite signs, x and y have opposite signs: We choose the positive root, because [ompR*] is always positive. |
</section> | </section> | ||
<section> | <section> | ||
Line 94: | Line 93: | ||
</br> | </br> | ||
- | <center><img src="https://static.igem.org/mediawiki/2012/f/f1/Parameters_nadia_grenoble.png" alt="" /></center> | + | <!--<center><img src="https://static.igem.org/mediawiki/2012/f/f1/Parameters_nadia_grenoble.png" alt="" /></center> |
- | <center><img src="https://static.igem.org/mediawiki/2012/7/7a/Parameters_nadia_grenoble2.png" alt="" /></center> | + | <center><img src="https://static.igem.org/mediawiki/2012/7/7a/Parameters_nadia_grenoble2.png" alt="" /></center>--> |
<table id="param_3"> | <table id="param_3"> | ||
Line 107: | Line 106: | ||
<tbody> | <tbody> | ||
<tr> | <tr> | ||
- | <td class="column1">Total quantity [ | + | <td class="column1">Total quantity [ompR]<span class="indice">tot</span></td> |
- | <td class="column2">6.8 | + | <td class="column2">6.8 10<span class="exposant">-8</span> mol.L<span class="exposant">-1</span></td> |
- | <td class="column3">The average number of | + | <td class="column3">The average number of ompR molecules per cell is 80.769 ± 0.719 <a href="#ref">[2]</a>. Knowing the cell volume (<i>v<span class="indice">c</span> = 1.1 10<span class="exposant">-15</span> L</i><a href="#ref">[3]</a>)<br/> and the Avogadro number <i>N<span class="indice">A</span> = 6.02 10<span class="exposant">-23</span> mol.L<span class="exposant">-1<span></i>, we deduce <br/>[ompR]<span class="indice">tot</span> = 80/(N<span class="indice">A</span>*v<span class="indice">c</span>) = 6.8 10<span class="exposant">-8</span> mol.L<span class="exposant">-1</span></td> |
</tr> | </tr> | ||
<tr> | <tr> | ||
<td class="column1">Goldbeter-Koshland model constants</td> | <td class="column1">Goldbeter-Koshland model constants</td> | ||
<td class="column2">v = 80 L<span class="exposant">-1</span>.min<span class="exposant">-1</span><br/> | <td class="column2">v = 80 L<span class="exposant">-1</span>.min<span class="exposant">-1</span><br/> | ||
- | V' = 7 | + | V' = 7 10<span class="exposant">-8</span> mol.L<span class="exposant">-1</span>.min<span class="exposant">-1</span><br/> |
- | K = 7 | + | K = 7 10<span class="exposant">-7</span> mol.L<span class="exposant">-1</span><br/> |
- | K' = 9 | + | K' = 9 10<span class="exposant">-8</span> mol.L<span class="exposant">-1</span></td> |
<td class="column3">We could not find these parameters in literature and we hope we will be able to conduct the necessary experiments to set them. Nevertheless, we could use a simple approach to estimate them : <br/> | <td class="column3">We could not find these parameters in literature and we hope we will be able to conduct the necessary experiments to set them. Nevertheless, we could use a simple approach to estimate them : <br/> | ||
- | The receptor should at least be sensitive to [dipeptide]=10<span class="exposant">-5</span>mol.L<span class="exposant">-1</span> (it represents the maximum concentration expected <a href="#ref">[4]</a>). We consider | + | The receptor should at least be sensitive to [dipeptide]=10<span class="exposant">-5</span> mol.L<span class="exposant">-1</span> (it represents the maximum concentration expected <a href="#ref">[4]</a>). We consider V = v[dipeptide], the equation is given by :<br/> |
<center><img src="https://static.igem.org/mediawiki/2012/1/16/DOmpR.png" alt="" /></center> | <center><img src="https://static.igem.org/mediawiki/2012/1/16/DOmpR.png" alt="" /></center> | ||
- | First of all, the value of K (resp K') should be in the same range of concentration as [ | + | First of all, the value of K (resp K') should be in the same range of concentration as [ompR]<span class="indice">tot</span>. Indeed, if K»[ompR]<span class="indice">tot</span> the phosphorylation term becomes negligible and the curve has not the desired evolution. Else if K«[ompR]<span class="indice">tot</span>, <img src="https://static.igem.org/mediawiki/2012/6/60/OmpR-K.png" alt="" /><br/> |
- | + | and we have a high phosphorylation rate even if almost all ompR has been phosphorylated. We chose K<span class="indice"><i>cyaA</i> </span> = 7 10<span class="exposant">-7</span> mol.L<span class="exposant">-1</span> and K'<span class="indice"><i>cyaA</i></span> = 9 10<span class="exposant">-8</span> mol.L<span class="exposant">-1</span><br/> | |
+ | Given that <img src="https://static.igem.org/mediawiki/2012/0/05/DOmpR-0.png" alt="" />,<br/> if we consider [ompR]~[omprR]<span class="indice">tot</span>≅6.8 10<span class="exposant">-8</span> and [ompR*]≅10<span class="exposant">-11</span>»[omprR]<span class="indice">tot</span> we find that V≥10<span class="exposant">4</span> V'. | ||
</td> | </td> | ||
</tr> | </tr> | ||
<tr> | <tr> | ||
- | <td class="column1">Maximal transcription rate of | + | <td class="column1">Maximal transcription rate of <i>cyaA</i></td> |
- | <td class="column2">V<span class="indice"> | + | <td class="column2">V<span class="indice">m<i>cyaA</i></span> = 2 10<span class="exposant">-9</span> mol.L<span class="exposant">-1</span>.min<span class="exposant">-1</span></td> |
<td class="column3">The value of this constant should be understood in the continuity of the network. For full details, consider the <a href="https://2012.igem.org/Team:Grenoble/Modeling/Amplification/Sensitivity#exp2">amplification section, parameters, explanation2</a></td> | <td class="column3">The value of this constant should be understood in the continuity of the network. For full details, consider the <a href="https://2012.igem.org/Team:Grenoble/Modeling/Amplification/Sensitivity#exp2">amplification section, parameters, explanation2</a></td> | ||
</tr | </tr | ||
<tr> | <tr> | ||
- | <td class="column1">Basal production of | + | <td class="column1">Basal production of Ca</td> |
- | <td class="column2">p<span class="indice"> | + | <td class="column2">p<span class="indice">Ca</span> = 2*10<span class="exposant">-12</span> mol.L<span class="exposant">-1</span>.min<span class="exposant">-1</span></td> |
<td class="column3">The value of this constant should be understood in the continuity of the network. For full details, consider the <a href="https://2012.igem.org/Team:Grenoble/Modeling/Amplification/Sensitivity#exp1">amplification section, parameters, explanation1</a></td> | <td class="column3">The value of this constant should be understood in the continuity of the network. For full details, consider the <a href="https://2012.igem.org/Team:Grenoble/Modeling/Amplification/Sensitivity#exp1">amplification section, parameters, explanation1</a></td> | ||
</tr> | </tr> | ||
<tr> | <tr> | ||
- | <td class="column1">Degradation rate of | + | <td class="column1">Degradation rate of Ca</td> |
- | <td class="column2">α<span class="indice"> | + | <td class="column2">α<span class="indice">Ca</span> = 6 10<span class="exposant">-3</span> min<span class="exposant">-1</span></td> |
<td class="column3">The value of this constant should be understood in the continuity of the network. For full details, consider the <a href="https://2012.igem.org/Team:Grenoble/Modeling/Amplification/Sensitivity#exp4">amplification section, parameters, explanation4</a></td> | <td class="column3">The value of this constant should be understood in the continuity of the network. For full details, consider the <a href="https://2012.igem.org/Team:Grenoble/Modeling/Amplification/Sensitivity#exp4">amplification section, parameters, explanation4</a></td> | ||
</tr> | </tr> | ||
<tr> | <tr> | ||
- | <td class="column1">Activation coefficient of | + | <td class="column1">Activation coefficient of <i>cyaA</i></td> |
- | <td class="column2">K<span class="indice"> | + | <td class="column2">K<span class="indice"><i>cyaA</i></span> = 10<span class="exposant">-7</span> mol.L<span class="exposant">-1</span></td> |
- | <td class="column3">The value K<span class="indice"> | + | <td class="column3">The value K<span class="indice">Ca</span> was set considering the maximum value of [Ca]. Indeed if we consider the steady state and assume that <i>p<span class="indice">Ca</span></i> is negligible compared to the other terms we have : <img src="https://static.igem.org/mediawiki/2012/f/f8/AC.png" alt="" /> where h stands for the Hill function, 0<h<1.<br/> |
- | We have then : .<br>K<span class="indice"> | + | We have then : <img src="https://static.igem.org/mediawiki/2012/a/a4/ACmax.png" alt="" /><br>K<span class="indice">Ca</span> should be in the same range as [Ca]<span class="indice">max</span> not too high otherwise the gene would never be expressed and not too low otherwise the protein is always produced. We chose K<span class="indice"><i>cyaA</i></span> = 10<span class="exposant">-7</span> mol.L<span class="exposant">-1</span></td> |
</tr> | </tr> | ||
<tr> | <tr> | ||
Line 170: | Line 170: | ||
<section> | <section> | ||
<h1>Sensitivity</h1> | <h1>Sensitivity</h1> | ||
- | + | ( [ompR])/( [ompR]<SUB>tot</SUB> ) is plotted as a function of the initial concentration of dipeptide. | |
</br> | </br> | ||
</br> | </br> | ||
- | <center><img src="https://static.igem.org/mediawiki/2012/ | + | <center><img src="https://static.igem.org/mediawiki/2012/4/48/OmpR_wiki.jpg" alt="" /></center> |
</br> | </br> | ||
</br> | </br> | ||
- | We notice that the receptor has a | + | We notice that the receptor has a saturated response at an initial dipeptide concentration of 10<SUP>-6</SUP> mol.L<span class="exposant">-1</span> . This concentration represents the expected concentration of dipeptide out of the cell [4]. Thus the sensitivity of the phosphorylation cascade is sufficient. |
</br> | </br> | ||
</br> | </br> | ||
- | + | In the next figure [Ca] is a represented as a function of initial dipeptide concentration. | |
</br> | </br> | ||
</br> | </br> | ||
- | <center><img src="https://static.igem.org/mediawiki/2012/c/ | + | <center><img src="https://static.igem.org/mediawiki/2012/c/ca/Cya_wiki.jpg" alt="" /></center> |
</br> | </br> | ||
</br> | </br> | ||
- | The sensitivity of the | + | The sensitivity of the sensor is 10<span class="exposant">-8</span> mol.L<span class="exposant">-1</span> of initial dipeptide concentration. However, at this level we are not able to assess the whole system sensitivity. We have got to wait for the amplification modeling results and link the two models to get an answer to the sensitivity question. |
</br> | </br> | ||
</br> | </br> | ||
Line 192: | Line 192: | ||
<section> | <section> | ||
<h1>Time response</h1> | <h1>Time response</h1> | ||
- | In order to | + | In order to determine the time response, we plotted the temporal evolution of Ca for two different initial dipeptide concentrations: 10<SUP>8</SUP> mol.L<span class="exposant">-1</span> and 10<SUP>-6</SUP> mol.L<span class="exposant">-1</span>. |
</br> | </br> | ||
</br> | </br> | ||
- | <center><img src="https://static.igem.org/mediawiki/2012/ | + | <center><img src="https://static.igem.org/mediawiki/2012/2/21/AC_time_8_wiki.jpg" alt="" /></center> |
</br> | </br> | ||
</br> | </br> | ||
</br> | </br> | ||
</br> | </br> | ||
- | <center><img src="https://static.igem.org/mediawiki/2012/ | + | <center><img src="https://static.igem.org/mediawiki/2012/2/21/AC_time_6_wiki.jpg" alt="" /></center> |
</br> | </br> | ||
</br> | </br> | ||
- | As we expected, we notice that we have the expected steady state values (see the graph above): 0.55 10-7 | + | <!--As we expected, we notice that we have the expected steady state values (see the graph above): 0.55 10-7 mol.L<span class="exposant">-1</span> for an initial dipeptide concentration of 10<SUP>-8</SUP> mol/L<span class="exposant">-1</span> and 1.1 10<SUP>-7</SUP> mol.L<span class="exposant"></span> for an initial dipeptide concentration of 10<SUP>-6</SUP> mol.L<span class="exposant">-1</span>.--> |
+ | In both cases, steady states is reached after 10h. | ||
</br> | </br> | ||
</br> | </br> | ||
- | Moreover, we see that to reach an | + | Moreover, we see that to reach an Ca production of 10<SUP>-7</SUP> mol.L<span class="exposan">-1</span> we need approximately 600 min. |
<a href="https://static.igem.org/mediawiki/2012/9/9e/Cya_vs_time.zip">Click here</a> to download the commented matlab code that gave us the temporal evolution. | <a href="https://static.igem.org/mediawiki/2012/9/9e/Cya_vs_time.zip">Click here</a> to download the commented matlab code that gave us the temporal evolution. | ||
Line 215: | Line 216: | ||
</br> | </br> | ||
</br> | </br> | ||
- | + | Since the amplification loop will both increase and accelerate the response, we expect that the sensitivity and response time determinated without the loop are the minimum compared to the whole system. | |
</section> | </section> | ||
</div> | </div> |
Latest revision as of 02:55, 27 September 2012
Signaling module
The design of signaling module is given by the figure below:Ordinary Differential Equations
Let’s begin by considering the cya gene activation by the transcription factor ompR*. As it is a gene activation, the regulated transcription rate is usually modelized by a Hill function:Parameters
Here is the link to the parameters of the amplification module we sometimes refer to.Constants | Value | Derivation |
---|---|---|
Total quantity [ompR]tot | 6.8 10-8 mol.L-1 | The average number of ompR molecules per cell is 80.769 ± 0.719 [2]. Knowing the cell volume (vc = 1.1 10-15 L[3]) and the Avogadro number NA = 6.02 10-23 mol.L-1, we deduce [ompR]tot = 80/(NA*vc) = 6.8 10-8 mol.L-1 |
Goldbeter-Koshland model constants | v = 80 L-1.min-1 V' = 7 10-8 mol.L-1.min-1 K = 7 10-7 mol.L-1 K' = 9 10-8 mol.L-1 |
We could not find these parameters in literature and we hope we will be able to conduct the necessary experiments to set them. Nevertheless, we could use a simple approach to estimate them : The receptor should at least be sensitive to [dipeptide]=10-5 mol.L-1 (it represents the maximum concentration expected [4]). We consider V = v[dipeptide], the equation is given by : and we have a high phosphorylation rate even if almost all ompR has been phosphorylated. We chose KcyaA = 7 10-7 mol.L-1 and K'cyaA = 9 10-8 mol.L-1 Given that , if we consider [ompR]~[omprR]tot≅6.8 10-8 and [ompR*]≅10-11»[omprR]tot we find that V≥104 V'. |
Maximal transcription rate of cyaA | VmcyaA = 2 10-9 mol.L-1.min-1 | The value of this constant should be understood in the continuity of the network. For full details, consider the amplification section, parameters, explanation2 | Basal production of Ca | pCa = 2*10-12 mol.L-1.min-1 | The value of this constant should be understood in the continuity of the network. For full details, consider the amplification section, parameters, explanation1 |
Degradation rate of Ca | αCa = 6 10-3 min-1 | The value of this constant should be understood in the continuity of the network. For full details, consider the amplification section, parameters, explanation4 |
Activation coefficient of cyaA | KcyaA = 10-7 mol.L-1 | The value KCa was set considering the maximum value of [Ca]. Indeed if we consider the steady state and assume that pCa is negligible compared to the other terms we have : where h stands for the Hill function, 0<h<1. We have then : KCa should be in the same range as [Ca]max not too high otherwise the gene would never be expressed and not too low otherwise the protein is always produced. We chose KcyaA = 10-7 mol.L-1 |
Hill Coefficient | n = 2 | We took a number greater than one to indicate positive cooperativity. |
References
- [1] Alejandra C.Ventura, Jacques-A. Sepulchre, Sofia D.Merajver. A Hidden Feedback in Signaling Cascades Is Revealed. PLOS Computational Biology, 2008, 4, 3, e1000041.
- [4] Michael D.Manson, Volker BlanK and Gabriele Brade. Peptide chemotaxis in E.Coli involves the Tap signal transducer and the dipeptide permease.Nature,15 May 1986,321,253-256.
- [5] Edith Gstrein-Reider and Manfred Schweiger, Institut fur Biochemie (nat. Fak.),UniversitAt Innsbruck, A-6020 Innsbruck, Austria. Regulation of adenylate cyclase in E. coli.