Team:Peking/Modeling/Phototaxis

From 2012.igem.org

(Difference between revisions)
m
Line 1: Line 1:
<html></p></html>{{Template:Peking2012_Color_Prologue}}{{Template:Peking2012_Color_Modeling}}<html>
<html></p></html>{{Template:Peking2012_Color_Prologue}}{{Template:Peking2012_Color_Modeling}}<html>
 +
<script type="text/javascript">
 +
sublists_Now = 2;
 +
var subsubitem=subfirst.getElementsByTagName('ul')[sublists_Now].getElementsByTagName('a')[0];
 +
subsubitem.style.color='#60b0f0';
 +
listTrigger(sublists_Now);
 +
subsubitem.innerHTML=subsubitem.innerHTML+' >>>';
 +
</script>
<div class="PKU_context floatR first">
<div class="PKU_context floatR first">
  <h3 id="title1">Summary</h3>
  <h3 id="title1">Summary</h3>
-
  <div>
+
  <p>
-
  <p>
+
Phototaxis is a light-control bio-system, whose input is the space-distribution of light. By comparison with chemotaxis system, phototaxis has much advantages for application.(see <a href="/Team:Peking/Project/Phototaxis">Phototaxis Page in Project</a>) We have constructed a simple phototaxis system coupling our <i>Luminesensor</i> with the expression level of cheZ protein. In order to confirm the macro light-control to our phototaxis system, we used the Mean-field PDE model. Later we managed to confirm these phenomena by tracing each cells in a micro way. We then did the related experiments to prove the effect of light in this simple system.
-
We have constructed a simple phototaxis system coupling our <i>Luminesensor</i> with the expression level of cheZ protein. After gathering the principles and parameters of the chemotaxis system, we then simulated our phototaxis system in a stochastic way. Based on our simulations, we predicted the outcome of the two demonstrations of phototaxis and presented the mechanism of phototaxis in a quantificational way.
+
</p>
-
  </p>
+
-
</div>
+
</div>
</div>
<div class="PKU_context floatR">
<div class="PKU_context floatR">
  <h3 id="title2">Phototaxis System</h3>
  <h3 id="title2">Phototaxis System</h3>
-
  <div>
+
  <p>
-
  <p>
+
Our phototaxis system functions as <i>Stopping on Light and Running in Dark</i>. As the sketch of this phototaxis system shows (Fig. 1), Light activates the <i>Luminesensor</i> which represses the expression of the CheZ protein. CheZ inactivates CheY<sub>P</sub>, which changes the rotation direction of the flagellum by protein-protein interaction and makes the bacteria tumbling, and reduces the tumbling frequency therefore. Bacteria moves slow with high tumbling frequency and <i>vice versa</i>.
Our phototaxis system functions as <i>Stopping on Light and Running in Dark</i>. As the sketch of this phototaxis system shows (Fig. 1), Light activates the <i>Luminesensor</i> which represses the expression of the CheZ protein. CheZ inactivates CheY<sub>P</sub>, which changes the rotation direction of the flagellum by protein-protein interaction and makes the bacteria tumbling, and reduces the tumbling frequency therefore. Bacteria moves slow with high tumbling frequency and <i>vice versa</i>.
-
  </p>
+
</p>
-
  <div class="floatC">
+
<div class="floatC">
[fig 1: Phototaxis Circuit]
[fig 1: Phototaxis Circuit]
-
  <p class="description">Fig 1. Phototaxis Circuit</p>
+
  <p class="description">Fig 1. Phototaxis Circuit</p>
-
  </div>
+
</div>
-
  <p>
+
<p>
To simplify the calculation, we assume the CheZ component responses immediately. When light reaches the bacteria, the concentration of CheZ behaves as Hill Function:
To simplify the calculation, we assume the CheZ component responses immediately. When light reaches the bacteria, the concentration of CheZ behaves as Hill Function:
-
  </p>
+
</p>
-
  <div class="floatC">
+
<div class="floatC">
[fig: CheZ Equation]
[fig: CheZ Equation]
[CheZ](I) = [CheZ]<sub>0</sub> * I<sub>0</sub> / (I + I<sub>0</sub>)
[CheZ](I) = [CheZ]<sub>0</sub> * I<sub>0</sub> / (I + I<sub>0</sub>)
-
  </div>
+
</div>
-
  <p>where</p><ul><li>
+
<p>where</p><ul><li>
-
[CheZ] denotes the concentration of CheZ</li><li>
+
[CheZ] : the concentration of CheZ</li><li>
-
[CheZ]<sub>0</sub> denotes the superior limit of CheZ concentration</li><li>
+
[CheZ]<sub>0</sub> : the superior limit of CheZ concentration</li><li>
-
I<sub>0</sub> denotes the critical illuminance</li><li>
+
I<sub>0</sub> : the critical illuminance</li><li>
-
I denotes the current illuminance</li></ul>
+
I : the current illuminance</li></ul>
-
  <p>
+
<p>
-
Then CheZ dephosphorylates CheY<sub>P</sub> into CheY while CheA phosphorylates CheY back. The typical time of dephosphorylation by CheZ is around 0.5 second and the typical time of phosphorylation by CheA (independent from light) is around 0.05 second.<sup><a href="#ref1" title="Binding of the Escherichia coli response regulator CheY to its target measured in vivo by fluorescence resonance energy transfer, Howard C. Berg, etc. PNAS">[1]</a></sup> By listing ODE equations, we can derive the equilibrium state of CheY<sub>P</sub> concentration as: <a href="/Team:Peking/Modeling/Phototaxis_Detail_1">(detail here)</a>
+
Then CheZ dephosphorylates CheY<sub>P</sub> into CheY while CheA phosphorylates CheY back. The typical time of dephosphorylation by CheZ is around 0.5 second and the typical time of phosphorylation by CheA (independent from light) is around 0.05 second.<sup><a href="#ref1" title="Binding of the Escherichia coli response regulator CheY to its target measured in vivo by fluorescence resonance energy transfer, Howard C. Berg, etc. PNAS">[1]</a></sup> By listing ODE equations, we can derive the equilibrium state of CheY<sub>P</sub> concentration as: <a href="/Team:Peking/Modeling/Appendix/Phototaxis">(detail here)</a>
<!--
<!--
CheY + CheAp -kY-> CheYp + CheA
CheY + CheAp -kY-> CheYp + CheA
Line 38: Line 42:
[CheY] = [CheY]t - [CheYp]
[CheY] = [CheY]t - [CheYp]
-->
-->
-
  </p>
+
</p>
-
  <div class="floatC">
+
<div class="floatC">
-
[fig: CheYp Equation]
+
  <img src="/wiki/images/d/d2/Peking2012_Formula001.png" alt="" />
-
[CheYp(CheZ) = kY * CheAp * CheYt / ( kY * CheAp + kZ * CheZ + gamY )]
+
</div>
-
  </div>
+
<p>where</p><ul><li>
-
  <p>where</p><ul><li>
+
[CheY<sub>P</sub>] denotes the concentration of phosphorylated CheY</li><li>
[CheY<sub>P</sub>] denotes the concentration of phosphorylated CheY</li><li>
[CheA<sub>P</sub>] denotes the steady concentration of active CheA</li><li>
[CheA<sub>P</sub>] denotes the steady concentration of active CheA</li><li>
Line 50: Line 53:
k<sub>Z</sub> denotes the rate constant of CheY dephosphorylation</li><li>
k<sub>Z</sub> denotes the rate constant of CheY dephosphorylation</li><li>
gamma-Y denotes the decay rate constant of CheY<sub>P</sub></li></ul>
gamma-Y denotes the decay rate constant of CheY<sub>P</sub></li></ul>
-
  <p>
+
<p>
CheY<sub>P</sub> can interact the flagellar motor to induce CW (clockwise) rotation. When flagellar motors rotate CCW (counterclockwise), they form a bundle to generate a force similar to a worm wheel. However, if some of the flagellar motors rotate CW (clockwise), the bundle breaks and the cell keeps tumbling. After in CW state for about 0.43s,<sup><a href="#ref2" title="Dependence of Bacterial Chemotaxis on Gradient Shape and Adaptation Rate, Nikita Vladimirov, etc. PLoS Computational Biology">[2]</a></sup> the flagellar motors return to CCW state and reconstruct the bundle to make the cell run. Since the CW state is triggered by CheY<sub>P</sub> molecule stochastically and is independent from its state history, this event is a typical <a href="/Team:Peking/Modeling/PoissonProcess">Possion Process</a> whose average frequency is determined by the concentration of CheY<sub>P</sub> with a Hill Function:<sup><a href="#ref3" title="An Ultrasensitive Bacterial Motor Revealed by Monitoring Signaling Proteins in Single Cells">[3]</a></sup>
CheY<sub>P</sub> can interact the flagellar motor to induce CW (clockwise) rotation. When flagellar motors rotate CCW (counterclockwise), they form a bundle to generate a force similar to a worm wheel. However, if some of the flagellar motors rotate CW (clockwise), the bundle breaks and the cell keeps tumbling. After in CW state for about 0.43s,<sup><a href="#ref2" title="Dependence of Bacterial Chemotaxis on Gradient Shape and Adaptation Rate, Nikita Vladimirov, etc. PLoS Computational Biology">[2]</a></sup> the flagellar motors return to CCW state and reconstruct the bundle to make the cell run. Since the CW state is triggered by CheY<sub>P</sub> molecule stochastically and is independent from its state history, this event is a typical <a href="/Team:Peking/Modeling/PoissonProcess">Possion Process</a> whose average frequency is determined by the concentration of CheY<sub>P</sub> with a Hill Function:<sup><a href="#ref3" title="An Ultrasensitive Bacterial Motor Revealed by Monitoring Signaling Proteins in Single Cells">[3]</a></sup>
-
  </p>
+
  </p>
-
  <div class="floatC">
+
-
[fig: CW Triggering Frequency Equation]
+
-
[FreqCW(CheYp) = pow(CheYp/CheYpc,N)/TUMBLE_TIME]
+
-
  </div>
+
-
  <p>where</p><ul><li>
+
-
FreqCW denotes the average frequency of CW (clockwise) rotation inducing</li><li>
+
-
[CheY<sub>Pc</sub>] denotes the critical concentration of phosphorylated CheY in this Hill Function</li><li>
+
-
N denotes the exponential rate of this Hill Function</li><li>
+
-
TUMBLE_TIME denotes the average relaxing time in a tumbling inducing</li></ul>
+
-
  </div>
+
-
</div>
+
-
<div class="PKU_context floatR">
+
-
<h3 id="title3">Phototaxis Simulation</h3>
+
-
<div>
+
-
  <p>
+
-
With the principles above, we construct our simulation system as following:
+
-
  </p>
+
-
  <ul><li>
+
-
(1) There are several bacteria cells in a room.
+
-
  </li><li>
+
-
(2) Cells can not run through the border of room.
+
-
  </li><li>
+
-
(3) The cells can divide in a random cell cycle in uniform distribution between 15min to 30min.
+
-
  </li><li>
+
-
(4) There are only two states of the cells --- running and tumbling.
+
-
  </li><li>
+
-
(5) Cells trigger tumbling as a Poisson Process, the average frequency is set by [CheY<sub>P</sub>] with the equation above.
+
-
  </li><li>
+
-
(6) Cells return running state after tumbling for a fixed time --- TUMBLE_TIME.
+
-
  </li><li>
+
-
(7) Cells run at a fixed speed --- v<sub>0</sub>.
+
-
  </li><li>
+
-
(8) In SPECS model, the running direction after tumbling is independent from previous direction;
+
-
while in RapidCell model, the new running direction performs as:<sup><a href="#ref2" title="Dependence of Bacterial Chemotaxis on Gradient Shape and Adaptation Rate, Nikita Vladimirov, etc. PLoS Computational Biology">[2]</a></sup>
+
-
  </li></ul>
+
-
  <div class="floatC">
+
-
[fig: Tumbling angle distribution in new running]<br/>
+
-
[rho(theta) = (1+cos(theta))*sin(theta)/2]
+
-
[rho(theta,phi) = (1+cos(theta))/4pi]
+
-
  </div>
+
-
  <p>where</p><ul><li>
+
-
theta denotes the tumbling angle (angle from origin direction to new direction)</li><li>
+
-
rho(theta) denotes the probability density of tumbling angle in value</li><li>
+
-
rho(theta,phi) denotes probability density of tumbling angle in the 3D space</li></ul>
+
-
  </p><p>
+
-
Parameters are shown as following:
+
-
  <div class="floatC">
+
-
  <table><tr>
+
-
    <td>Parameter</td><td>Value</td><td>Unit</td><td>Description</td><td>Source</td>
+
-
  </tr><tr>
+
-
    <td>v<sub>0</sub></td><td>20</td><td>um/s</td><td>running speed</td><td><sup><a href="#ref4" title="Chemotaxis in Escherichia coli analysed by Three-dimensional Tracking, Howard C.Berg, Douglas A.Brown, NATURE">[4]</a></sup></td>
+
-
  </tr><tr>
+
-
    <td>TUMBLE_TIME</td><td>0.43</td><td>s</td><td>time during a tumbling</td><td><sup><a href="#ref5" title="Real-Time Imaging of Fluorescent Flagellar Filaments, Linda Turner, etc. JOURNAL OF BACTERIOLOGY">[5]</a></sup></td>
+
-
  </tr><tr>
+
-
    <td>CELL_PERIOD</td><td>15~30</td><td>min</td><td>period of a cell cycle</td><td></td>
+
-
  </tr><tr>
+
-
    <td>[CheA]<sub>T</sub></td><td>5.3</td><td>u mol/L</td><td>total concentration of CheA</td><td><sup><a href="#ref6" title="Relationship between cellular response and behavioral variability in bacterial chemotaxis, Thierry Emonet, Philippe Cluzel. PNAS">[6]</a></sup></td>
+
-
  </tr><tr>
+
-
    <td>[CheZ]<sub>c</sub></td><td>1.1</td><td>u mol/L</td><td>typical concentration of CheZ</td><td><sup><a href="#ref1" title="Binding of the Escherichia coli response regulator CheY to its target measured in vivo by fluorescence resonance energy transfer, Victor Sourjik and Howard C. Berg, PNAS">[1]</a></sup></td>
+
-
  </tr><tr>
+
-
    <td>[CheY]<sub>T</sub></td><td>9.7</td><td>u mol/L</td><td>total concentration of CheY</td><td><sup><a href="#ref6" title="Relationship between cellular response and behavioral variability in bacterial chemotaxis, Thierry Emonet, Philippe Cluzel. PNAS">[6]</a></sup></td>
+
-
  </tr><tr>
+
-
    <td>k<sub>Y</sub></td><td>100</td><td>(u mol/L)<sup>-1</sup> s<sup>-1</sup></td><td>phosphorylation rate constant of CheY</td><td><sup><a href="#ref6" title="Relationship between cellular response and behavioral variability in bacterial chemotaxis, Thierry Emonet, Philippe Cluzel. PNAS">[6]</a></sup></td>
+
-
  </tr><tr>
+
-
    <td>k<sub>Z</sub></td><td>30/[CheZ]<sub>c</sub></td><td>(u mol/L)<sup>-1</sup> s<sup>-1</sup></td><td>dephosphorylation rate constant of CheY</td><td><sup><a href="#ref6" title="Relationship between cellular response and behavioral variability in bacterial chemotaxis, Thierry Emonet, Philippe Cluzel. PNAS">[6]</a></sup></td>
+
-
  </tr><tr>
+
-
    <td>gamma<sub>Y</sub></td><td>0.1</td><td>s<sup>-1</sup></td><td>decay rate constant of CheY<sub>P</sub></td><td><sup><a href="#ref2" title="Dependence of Bacterial Chemotaxis on Gradient Shape and Adaptation Rate, Nikita Vladimirov, etc. PLoS Computational Biology">[2]</a></sup></td>
+
-
  </tr><tr>
+
-
    <td>N</td><td>10.3</td><td>1</td><td>the exponential rate of Hill Function of CW (clockwise) bias</td><td><sup><a href="#ref7" title="An Ultrasensitive Bacterial Motor Revealed by Monitoring Signaling Proteins in Single Cells, Philippe Cluzel, etc. Science">[7]</a></sup></td>
+
-
  </tr><tr>
+
-
    <td>[CheY]<sub>Pc</sub></td><td>3.1</td><td>u mol/L</td><td>the critical concentration of phosphorylated CheY of Hill Function of CW (clockwise) bias</td><td><sup><a href="#ref7" title="An Ultrasensitive Bacterial Motor Revealed by Monitoring Signaling Proteins in Single Cells, Philippe Cluzel, etc. Science">[7]</a></sup></td>
+
-
  </tr><tr>
+
-
    <td>r<sub>A</sub></td><td>1/3</td><td>1</td><td>phosphorylation rate of CheA</td><td></td>
+
-
  </tr></table>
+
-
[table 1: Simulation Parameters]
+
-
  <p class="description">Tab 1. </p>
+
-
  </div>
+
-
</div>
+
-
</div>
+
-
<div class="PKU_context floatR">
+
-
<h3 id="title4">Result 1: Half-light-half-dark Room</h3>
+
-
<div>
+
-
  <p>
+
-
Our first Demonstration is in a Half-light-half-dark plate, and we would like to see how cells behave differently in such a high contrast environment.
+
-
The lighting of light room is set to 0.8 unit while the dark is set to 0.1 unit with I<sub>0</sub> = 0.5. Here goes the results:
+
-
  </p>
+
-
  <div>
+
-
[fig 2: diffusion from center]
+
-
[fig 3: initial uniform distribution]
+
-
  <p class="description">Fig 2,3 </p>
+
-
  </div>
+
-
  <p>
+
-
Since the frequency of tumbling in light area is much higher than in dark area, the diffusion of population in light area is much smaller. If we initialize the room with cells in uniform distribution, a high population band will emerge at the border in light area. Our experiments show:
+
-
  </p>
+
-
  <div>
+
-
[fig 4: diffusion from center (experiment)]
+
-
[fig 5: initial uniform distribution (experiment)]
+
-
  <p class="description">Fig 4,5 </p>
+
-
  </div>
+
-
  <p>
+
-
which fit our predictions by modeling.
+
-
  </p>
+
-
</div>
+
-
</div>
+
-
<div class="PKU_context floatR">
+
-
<h3 id="title5">Result 2: Light Gradient Room</h3>
+
-
<div>
+
-
  <p>
+
-
Phototaxis is designed to move cells in a given direction. Just like diffusion (SPECS model in a large population can derive the diffusion equation<sup><a href="#ref8" title="A Pathway-based Mean-field Model for Escherichia coli Chemotaxis, Tailin Wu, etc. PACS">[8]</a></sup>), the movement order requires a gradient lighting field in the room. We set the lighting from 0 to 1 unit in 1 mm, then discovered the directed movement bias towards light area in this simulation.
+
-
  </p>
+
-
  <div>
+
-
[fig 6: Gradient Lighting]
+
-
  <p class="description">Fig 6 </p>
+
-
  </div>
+
-
  <p>
+
-
Then we do this movement experiment in a much larger scale, and the bacteria successfully response with their motion.
+
-
  </p>
+
-
</div>
+
</div>
</div>
</html>{{Template:Peking2012_Color_Epilogue}}
</html>{{Template:Peking2012_Color_Epilogue}}

Revision as of 07:50, 24 September 2012

Summary

Phototaxis is a light-control bio-system, whose input is the space-distribution of light. By comparison with chemotaxis system, phototaxis has much advantages for application.(see Phototaxis Page in Project) We have constructed a simple phototaxis system coupling our Luminesensor with the expression level of cheZ protein. In order to confirm the macro light-control to our phototaxis system, we used the Mean-field PDE model. Later we managed to confirm these phenomena by tracing each cells in a micro way. We then did the related experiments to prove the effect of light in this simple system.

Phototaxis System

Our phototaxis system functions as Stopping on Light and Running in Dark. As the sketch of this phototaxis system shows (Fig. 1), Light activates the Luminesensor which represses the expression of the CheZ protein. CheZ inactivates CheYP, which changes the rotation direction of the flagellum by protein-protein interaction and makes the bacteria tumbling, and reduces the tumbling frequency therefore. Bacteria moves slow with high tumbling frequency and vice versa.

[fig 1: Phototaxis Circuit]

Fig 1. Phototaxis Circuit

To simplify the calculation, we assume the CheZ component responses immediately. When light reaches the bacteria, the concentration of CheZ behaves as Hill Function:

[fig: CheZ Equation] [CheZ](I) = [CheZ]0 * I0 / (I + I0)

where

  • [CheZ] : the concentration of CheZ
  • [CheZ]0 : the superior limit of CheZ concentration
  • I0 : the critical illuminance
  • I : the current illuminance

Then CheZ dephosphorylates CheYP into CheY while CheA phosphorylates CheY back. The typical time of dephosphorylation by CheZ is around 0.5 second and the typical time of phosphorylation by CheA (independent from light) is around 0.05 second.[1] By listing ODE equations, we can derive the equilibrium state of CheYP concentration as: (detail here)

where

  • [CheYP] denotes the concentration of phosphorylated CheY
  • [CheAP] denotes the steady concentration of active CheA
  • [CheYT] denotes the total concentration of CheY
  • kY denotes the rate constant of CheY phosphorylation
  • kZ denotes the rate constant of CheY dephosphorylation
  • gamma-Y denotes the decay rate constant of CheYP

CheYP can interact the flagellar motor to induce CW (clockwise) rotation. When flagellar motors rotate CCW (counterclockwise), they form a bundle to generate a force similar to a worm wheel. However, if some of the flagellar motors rotate CW (clockwise), the bundle breaks and the cell keeps tumbling. After in CW state for about 0.43s,[2] the flagellar motors return to CCW state and reconstruct the bundle to make the cell run. Since the CW state is triggered by CheYP molecule stochastically and is independent from its state history, this event is a typical Possion Process whose average frequency is determined by the concentration of CheYP with a Hill Function:[3]

  • Totop Totop

Retrieved from "http://2012.igem.org/Team:Peking/Modeling/Phototaxis"