Team:ETH Zurich/modeling/diffmodel

From 2014.igem.org

(Difference between revisions)
(Reactions)
(Pattern developing)
 
(33 intermediate revisions not shown)
Line 1: Line 1:
{{:Team:ETH Zurich/tpl/head|Diffusion model}}
{{:Team:ETH Zurich/tpl/head|Diffusion model}}
 +
{{:Team:ETH Zurich/tpl/fortables}}
{{:Team:ETH Zurich/tpl/read-more}}
{{:Team:ETH Zurich/tpl/read-more}}
<center>
<center>
Line 14: Line 15:
=== Geometry and boundary conditions ===
=== Geometry and boundary conditions ===
-
[[File:ETH Zurich Compartments Diffusion.png|float|500px|thumb|Millifluidic compartments used for the diffusion model]]
+
[[File:ETH Zurich Compartments Diffusion.png|float|500px|thumb|'''Figure 1''' Millifluidic compartments used for the diffusion model]]
Line 29: Line 30:
-
[[File:ETHZ comsol geometry.png|center|610px|thumb|Geometry used for simulation on Comsol Multiphysics]]
+
[[File:ETHZ comsol geometry.png|center|610px|thumb|'''Figure 2''' Geometry used for simulation on Comsol Multiphysics]]
<br>
<br>
-
In order to fulfill these conditions in Comsol Multiphysics, the cells are not drawn but instead, beads are considered as compartments with a certain density of cells. Rates of diffusion through the membrane are added as reaction rates and intracellular species are contained in beads via a Newmann's boundary condition applied to the surface of the beads, that is :
+
In order to fulfill these conditions in Comsol Multiphysics, the cells are not drawn but instead, beads are considered as compartments with a certain density of cells. Rates of diffusion through the membrane are added as reaction rates and intracellular species are contained in beads via a Neumann boundary condition applied to the surface of the beads, that is :
-
$$ \nabla C . \textbf{n} =0$$
+
$$ \nabla C \cdot \textbf{n} =0$$
$$\text{for all species except extracellular AHL}$$
$$\text{for all species except extracellular AHL}$$
$$\text{for every vector n normal to the surface of a bead}$$
$$\text{for every vector n normal to the surface of a bead}$$
Line 72: Line 73:
=== Reactions ===
=== Reactions ===
-
These reactions happen in every well of the pattern on the grid.
+
These reactions happen in every bead of the pattern on the grid.
$$ \begin{align}
$$ \begin{align}
Line 80: Line 81:
RLux+RLux &\leftrightarrow DRLux\\
RLux+RLux &\leftrightarrow DRLux\\
DRLux+P_{luxOFF} & \leftrightarrow P_{luxON}\\
DRLux+P_{luxOFF} & \leftrightarrow P_{luxON}\\
-
P_{luxON}&\rightarrow P_{luxON}+mRNA_{GFP}\\
+
P_{luxON}&\rightarrow P_{luxON}+mRNA_{GFP}+mRNA_{LuxI}\\
mRNA_{GFP}&\rightarrow GFP \\
mRNA_{GFP}&\rightarrow GFP \\
-
LuxAHL &\rightarrow \\
+
mRNA_{LuxI} &\rightarrow LuxI\\
 +
AHLint &\rightarrow \\
 +
AHLext &\rightarrow \\
LuxR &\rightarrow  \\
LuxR &\rightarrow  \\
RLux &\rightarrow\\
RLux &\rightarrow\\
Line 142: Line 145:
\begin{align*}
\begin{align*}
\frac{d[AHLext]}{dt} &= \frac{N_0 N_m e^{rt}}{N_m+N_0(e^{rt}-1)} \alpha\  D_m \ (AHL_{int}-AHL_{ext}) -d_{AHLext}[AHL_{ext}]\\
\frac{d[AHLext]}{dt} &= \frac{N_0 N_m e^{rt}}{N_m+N_0(e^{rt}-1)} \alpha\  D_m \ (AHL_{int}-AHL_{ext}) -d_{AHLext}[AHL_{ext}]\\
-
\frac{d[AHLint]}{dt} &= Dm (AHL_{ext}-AHL_{int}) + k_{-RLux}[R_{Lux}]-k_{RLux}[LuxAHL_{int}][LuxR]-d_{AHLint}[AHL_{int}]\\
+
\frac{d[AHLint]}{dt} &= Dm (AHL_{ext}-AHL_{int}) + k_{-RLux}[R_{Lux}]-k_{RLux}[AHL_{int}][LuxR]-d_{AHLint}[AHL_{int}]\\
\frac{d[LuxR]}{dt} &= \alpha_{LuxR} -k_{RLux}[AHL_{int}][LuxR] + k_{-RLux}[RLux] - d_{LuxR}[LuxR] \\
\frac{d[LuxR]}{dt} &= \alpha_{LuxR} -k_{RLux}[AHL_{int}][LuxR] + k_{-RLux}[RLux] - d_{LuxR}[LuxR] \\
\frac{d[RLux]}{dt} &=  k_{RLux}[AHL_{int}][LuxR] - k_{-RLux}[RLux]  - d_{RLux} [RLux] \\  
\frac{d[RLux]}{dt} &=  k_{RLux}[AHL_{int}][LuxR] - k_{-RLux}[RLux]  - d_{RLux} [RLux] \\  
\frac{d[mRNA_{GFP}]}{dt} &=  \frac{k_{mRNA_{GFP}}[RLux]^2}{K_{mLux}^2 + [RLux]^2}- d_{mRNA_{GFP}} [mRNA_{GFP}]\\
\frac{d[mRNA_{GFP}]}{dt} &=  \frac{k_{mRNA_{GFP}}[RLux]^2}{K_{mLux}^2 + [RLux]^2}- d_{mRNA_{GFP}} [mRNA_{GFP}]\\
-
\frac{d[GFP]}{dt} &=  k_{GFP} [mRNA_{GFP}] - d_{GFP}[GFP]\\
+
\frac{d[GFP]}{dt} &=  k_{GFP} [mRNA_{GFP}] - d_{GFP}[GFP]
\end{align*}
\end{align*}
$$
$$
Line 163: Line 166:
$$ [LuxR]_i=\frac{a_{LuxR}}{d_{LuxR}} $$
$$ [LuxR]_i=\frac{a_{LuxR}}{d_{LuxR}} $$
The cells of the first row are induced by 10 &mu; M AHL just before they are encapsulated in alginate. Then, bead formation takes half-an-hour before these cells are added in the first well. During this half-an-hour AHL can diffuse out. Therefore initial conditions in the first cell completely depend on these previous steps. We have also simulated this bead preparation in order to get these initial conditions.
The cells of the first row are induced by 10 &mu; M AHL just before they are encapsulated in alginate. Then, bead formation takes half-an-hour before these cells are added in the first well. During this half-an-hour AHL can diffuse out. Therefore initial conditions in the first cell completely depend on these previous steps. We have also simulated this bead preparation in order to get these initial conditions.
 +
 +
[[File:ETHZ Beads in storage solution.jpg|center|300px|thumb|'''Figure 3''' Beads stored in calcium chloride, 10mM]]
The geometry is a simple bead in a bigger compartment filled with calcium chloride. We look at the concentrations at the middle of the bead. Initial concentration of LuxR before bead diffusion starts is a classical steady state a<sub>LuxR</sub>/d<sub>LuxR</sub>
The geometry is a simple bead in a bigger compartment filled with calcium chloride. We look at the concentrations at the middle of the bead. Initial concentration of LuxR before bead diffusion starts is a classical steady state a<sub>LuxR</sub>/d<sub>LuxR</sub>
-
[[File:ETHZ Signaling molecules.png|600px|center]]
+
[[File:ETHZ Signaling molecules.png|600px|center|thumb|'''Figure 4''' Signaling molecules]]
-
[[File:ETHZ Regulators and GFP.png|600px|center]]
+
[[File:ETHZ Regulators and GFP.png|600px|center|thumb|'''Figure 5''' Regulators and GFP]]
Extracellular AHL diffuses very fast through the membrane, which makes intracellular AHL increase, until both reach the same value and diffusion through the membrane reaches an equilibrium. From this point on AHL starts to bind to LuxR to form RLux, so AHL decreases, until RLux triggers enough production of LuxI which catalyses again AHL production, and makes AHL concentration become stable.
Extracellular AHL diffuses very fast through the membrane, which makes intracellular AHL increase, until both reach the same value and diffusion through the membrane reaches an equilibrium. From this point on AHL starts to bind to LuxR to form RLux, so AHL decreases, until RLux triggers enough production of LuxI which catalyses again AHL production, and makes AHL concentration become stable.
-
After 30 minutes, almost all LuxR is bound to AHL and converted to RLux, so that we have  
+
After 30 minutes, all LuxR is bound to AHL and converted to RLux, so that we have  
-
$$ [LuxR]_i= 10nM$$
+
$$ [LuxR]_i= 0nM$$
-
$$ [RLux]_i= 160 nM$$
+
$$ [RLux]_i= 200 nM$$
 +
$$ [AHLint]_i=[AHLint]_i=2 &mu; M$$
RLux has already started to trigger the production of LuxI and GFP, so that we have  
RLux has already started to trigger the production of LuxI and GFP, so that we have  
-
$$ [mRNAGFP]_i =  24nM$$
+
$$ [mRNAGFP]_i =  40nM$$
-
$$ [GFP]_i =  180nM$$
+
$$ [GFP]_i =  900nM$$
-
$$ [mRNALuxI]_i = 24nM$$
+
$$ [mRNALuxI]_i = 40nM$$
-
$$ [LuxI]_i = 3300nM$$
+
$$ [LuxI]_i = 16000nM$$
<html><div id='Estimation'></html>
<html><div id='Estimation'></html>
 +
=== Estimation of parameters from literature ===
=== Estimation of parameters from literature ===
-
The initial number of beads is 10 million. According to Lars Müller's master thesis<sup>[[Team:ETH_Zurich/project/references|[29]]]</sup>, in picoliter beads, cells doubling time is 30 minutes. Here we are using beads with a volume in the microliter range. Because of bead volume, oxygen and nutrients are much less accessible. Therefore, we multiplied this doubling time by 4. We have a rgowth rate of 0.006 min<sup>-1</sup> which is still above the growth rate in anaerobic conditions (0.004 min<sup>-1</sup> according to [http://bionumbers.hms.harvard.edu/search.aspx?log=y&task=searchbytrmorg&trm=growth+rate+e+coli&org= Bionumbers]) )
+
The initial number of beads is 10 million. According to Lars Müller's master thesis<sup>[[Team:ETH_Zurich/project/references|[29]]]</sup>, in picoliter beads, cells doubling time is 30 minutes. Here we are using beads with a volume in the microliter range. Because of bead volume, oxygen and nutrients are much less accessible. Therefore, we multiplied this doubling time by 4. We have a growth rate of 0.006 min<sup>-1</sup> which is still above the growth rate in anaerobic conditions (0.004 min<sup>-1</sup> according to [http://bionumbers.hms.harvard.edu/search.aspx?log=y&task=searchbytrmorg&trm=growth+rate+e+coli&org= Bionumbers])  
-
Lars Müller also mentions in his master thesis<sup>[[Team:ETH_Zurich/project/references|[29]]]</sup> that the maximum capacity of his 34 pL beads is 3000 cells, which would correspond to  maximum of N<sub>m</sub> = 8. 10<sup>8</sup> cells per bead in our case (10 &mu;L beads).  
+
In the conditions of lack of oxygen and nutrients where our cells are, we consider that they might at maximum double or triple. So we take 8 . 10<sup>7</sup> cells per bead for N<sub>m</sub>.
Line 198: Line 205:
-
<html></article></html>
+
<html></div></html>
Line 212: Line 219:
The background in our experimental setup is very high and this is high we see actual fluorescence appear only after 11 hours. In order to account for this, we also set up a background in the simulated pattern by adjusting the scale.  
The background in our experimental setup is very high and this is high we see actual fluorescence appear only after 11 hours. In order to account for this, we also set up a background in the simulated pattern by adjusting the scale.  
-
------------ Video  
+
 
-
--------------------
+
{|class="wikitable" style="background-color: white; text-align:center; width:auto; margin: auto; font-size:10pt;"
 +
|colspan="2" style='font-size:10pt';text-align:left|{{:Team:ETH_Zurich/Templates/Video|width=1080px|id=video3|ratio=1920/720|srcMP4=<html>https://static.igem.org/mediawiki/2014/b/b1/ETH_Zurich_2014_signal_propagation_with_simulation.mp4</html>|poster=<html>https://static.igem.org/mediawiki/2014/6/69/ETH_Zurich_2014_signal_propagation_with_simulation_preview.png</html>}}
 +
|-
 +
|colspan="2" style='font-size:10pt';text-align:left|'''Video 1''' '''Row wise, self-propagating [https://2014.igem.org/Team:ETH_Zurich/project/background/biotools#Quorum_Sensing cell-to-cell communication] of ''E. coli'' cells confined in [https://2014.igem.org/Team:ETH_Zurich/lab/bead alginate beads] (d=3 mm, initially 10<sup>7</sup> cells/bead) on a [https://2014.igem.org/Team:ETH_Zurich/lab/chip custom-made millifluidic PDMS chip].'''
 +
|-
 +
|style="width:50%"|All fluorescent cells contained [https://2014.igem.org/Team:ETH_Zurich/expresults/rr#Riboregulators riboregulated] sfGFP followed by [http://parts.igem.org/Part:BBa_C0161 LuxI (BBa_C0161)] together under the control of the [http://parts.igem.org/Part:BBa_R0062 pLux promoter (BBa_R0062)], and [http://parts.igem.org/Part:BBa_J23100 constitutively (BBa_J23100)] expressed [http://parts.igem.org/Part:BBa_C0062 LuxR (BBa_C0062)]. LuxI catalyzes the production of the autoinducer 3OC6-HSL, which is then diffusing from cell to cell. For initialization, the cells in one bead of the top row were induced with 3OC6-HSL before encapsulation. Imaging was implemented with a [https://2014.igem.org/Team:ETH_Zurich/lab/protocols#Biostep_Dark-Hood_DH-50.E2.84.A2__and_the_Argus-X1.E2.84.A2_software Biostep Dark-Hood DH-50 (Argus X1 software)] fitted with a Canon EOS 500D DSLR camera and a fluorescence filter (545 nm filter). Pictures were usually taken every 2 min at an excitation wavelength of 470 nm with the standard Canon EOS Utility software. Time-lapse movies were created with Adobe After Effects CC software. 1950x faster than real-time, the video shown starts 10 h after the initiation of the experiment (however the time scale shown corresponds to minutes after loading of the chip). For precise experimental setup, check the [https://2014.igem.org/Team:ETH_Zurich/expresults#Diffusion Results] page.
 +
||Simulation of the propagation of the pattern in the millifluidic chip. [http://www.comsol.com/comsol-multiphysics Comsol Multiphysics Simulation software] was used in order to simulate a detailed diffusion model including quorum sensing steps in colonies and cell growth. Overall GFP concentration in beads has been scaled in order to account for the high background of the experimental setup. Green Fluorescence Protein is produced earlier in the wells, but can be seen only above a certain threshold.<br>Accurate prediction of experimental data by the model has been achieved, with parameters from our own fittings or from the literature. Experimental observation combined with simulation enables to show that a pattern is able to develop in the millifluidic chip in a reasonable time scale.
 +
|}
=== Dynamics by row ===
=== Dynamics by row ===
Line 221: Line 235:
<br>
<br>
-
<!--    Update graphs with new parameters
+
[[File:ETHZ AHL-regulators-pattern.png|center|900px|thumb|'''Figure 6''']]
-
 
+
-
[[File:ETHZ AHL-regulators-pattern.png|center]]
+
-
-->
+
<br>
<br>
Line 232: Line 243:
<br>
<br>
-
<!--    Update graphs with new parameters
+
[[File:ETHZ GFP in beads.png|600px|center|thumb|'''Figure 7''' GFP in beads ]]
-
[[File:ETHZ GFP in beads.png|600px|center]]
+
[[File:ETHZ GFP in single cell.png|600px|center|thumb|'''Figure 8''' GFP in single cell]]
-
[[File:ETHZ GFP in single cell.png|600px|center]]
+
 
-
-->
+
{{:Team:ETH_Zurich/tpl/topbutton|green}}
{{:Team:ETH_Zurich/tpl/topbutton|green}}

Latest revision as of 01:13, 18 October 2014

iGEM ETH Zurich 2014