Team:Tsinghua-A/Modeling

From 2011.igem.org

(Difference between revisions)
Line 115: Line 115:
the full text </U></I></FONT></A></SPAN></FONT>
the full text </U></I></FONT></A></SPAN></FONT>
</P>
</P>
-
<P ALIGN=LEFT>  <FONT FACE="Arial Unicode MS, sans-serif"><SPAN LANG="en-US">In
+
<P ALIGN=LEFT>  <FONT FACE="Arial Unicode MS, sans-serif"><SPAN LANG="en-US">&nbsp&nbspIn
our project, we designed a quorum-sensing oscillator which consists
our project, we designed a quorum-sensing oscillator which consists
of two types of cells. The expression of the reporter genes (GFP of
of two types of cells. The expression of the reporter genes (GFP of
Line 154: Line 154:
</B></SPAN></FONT></FONT></FONT><A HREF="https://2011.igem.org/Team:Tsinghua-A/Modeling/P2A"><FONT COLOR="#ff950e"><FONT FACE="Times New Roman, serif"><FONT SIZE=3><SPAN LANG="en-US"><I><U><SPAN STYLE="font-weight: normal">Read
</B></SPAN></FONT></FONT></FONT><A HREF="https://2011.igem.org/Team:Tsinghua-A/Modeling/P2A"><FONT COLOR="#ff950e"><FONT FACE="Times New Roman, serif"><FONT SIZE=3><SPAN LANG="en-US"><I><U><SPAN STYLE="font-weight: normal">Read
more</SPAN></U></I></SPAN></FONT></FONT></FONT></A></SPAN></FONT></P>
more</SPAN></U></I></SPAN></FONT></FONT></FONT></A></SPAN></FONT></P>
-
<P ALIGN=LEFT> <FONT FACE="Arial Unicode MS, sans-serif"><SPAN LANG="en-US">The
+
<P ALIGN=LEFT><FONT COLOR="#000000"><FONT FACE="Arial Unicode MS, sans-serif"><SPAN LANG="en-US">The
original model contains too many factors for analyzing the general
original model contains too many factors for analyzing the general
property of system. To understand the essential characters of the
property of system. To understand the essential characters of the
-
oscillator, we <FONT COLOR="#ff0000"><B>simplify the original model</B></FONT>
+
oscillator, we <FONT COLOR="#ff0000"><B>simplify the original model</B></FONT><FONT COLOR="#000000">
according to certain appropriate assumptions, like Quasi-equilibrium
according to certain appropriate assumptions, like Quasi-equilibrium
for fast reactions.</SPAN></FONT></P>
for fast reactions.</SPAN></FONT></P>
-
<P ALIGN=LEFT> <FONT FACE="Arial Unicode MS, sans-serif"><SPAN LANG="en-US">After
+
<P ALIGN=LEFT><FONT COLOR="#000000"><FONT FACE="Arial Unicode MS, sans-serif"><SPAN LANG="en-US">After
series of derivation based on those assumptions (see the attached pdf
series of derivation based on those assumptions (see the attached pdf
file), we came up with the following set of DDEs (Delay Differential
file), we came up with the following set of DDEs (Delay Differential
Equations)</SPAN></FONT></P>
Equations)</SPAN></FONT></P>
<P ALIGN=CENTER><IMG SRC="https://static.igem.org/mediawiki/2011/b/b0/003.png" NAME="图形4" ALT="说明: C:\Users\Harry\AppData\Roaming\Tencent\Users\806995810\QQ\WinTemp\RichOle\RH%W~L%6JQ)ZX06RH76VYV3.jpg" ALIGN=BOTTOM WIDTH=531 HEIGHT=312 BORDER=0></P>
<P ALIGN=CENTER><IMG SRC="https://static.igem.org/mediawiki/2011/b/b0/003.png" NAME="图形4" ALT="说明: C:\Users\Harry\AppData\Roaming\Tencent\Users\806995810\QQ\WinTemp\RichOle\RH%W~L%6JQ)ZX06RH76VYV3.jpg" ALIGN=BOTTOM WIDTH=531 HEIGHT=312 BORDER=0></P>
-
<P ALIGN=LEFT> <FONT FACE="Times New Roman, serif"><SPAN LANG="en-US"><FONT FACE="Arial Unicode MS, sans-serif"><SPAN LANG="en-US">We
+
<P ALIGN=LEFT><FONT COLOR="#000000"><FONT FACE="Times New Roman, serif"><SPAN LANG="en-US"><FONT FACE="Arial Unicode MS, sans-serif"><SPAN LANG="en-US">We
coded the system by DDE description in MATLAB and did simulation
coded the system by DDE description in MATLAB and did simulation
analysis accordingly. The result showed that the system could
analysis accordingly. The result showed that the system could
oscillate under certain parameters .</SPAN></FONT></SPAN></FONT></P>
oscillate under certain parameters .</SPAN></FONT></SPAN></FONT></P>
<P ALIGN=CENTER><IMG SRC="https://static.igem.org/mediawiki/2011/d/d5/004.png" NAME="图形5" ALT="说明: D:\b1.bmp" ALIGN=BOTTOM WIDTH=395 HEIGHT=296 BORDER=0></P>
<P ALIGN=CENTER><IMG SRC="https://static.igem.org/mediawiki/2011/d/d5/004.png" NAME="图形5" ALT="说明: D:\b1.bmp" ALIGN=BOTTOM WIDTH=395 HEIGHT=296 BORDER=0></P>
-
<P ALIGN=LEFT> <FONT FACE="Arial Unicode MS, sans-serif"><SPAN LANG="en-US">To
+
<P ALIGN=LEFT><FONT COLOR="#000000"><FONT FACE="Arial Unicode MS, sans-serif"><SPAN LANG="en-US">To
further understand what parameters could make the system oscillate,
further understand what parameters could make the system oscillate,
we did bifurcation analysis on the Hill parameters. What we had to do
we did bifurcation analysis on the Hill parameters. What we had to do
Line 176: Line 176:
but a little disruption may lead to a steady state like that:</SPAN></FONT></P>
but a little disruption may lead to a steady state like that:</SPAN></FONT></P>
<P ALIGN=CENTER><IMG SRC="https://static.igem.org/mediawiki/2011/2/20/005.png" NAME="图形6" ALT="说明: C:\Users\Harry\AppData\Roaming\Tencent\Users\806995810\QQ\WinTemp\RichOle\)D6IYP[O_UD]}LB_9)SO(Q5.jpg" ALIGN=BOTTOM WIDTH=557 HEIGHT=260 BORDER=0></P>
<P ALIGN=CENTER><IMG SRC="https://static.igem.org/mediawiki/2011/2/20/005.png" NAME="图形6" ALT="说明: C:\Users\Harry\AppData\Roaming\Tencent\Users\806995810\QQ\WinTemp\RichOle\)D6IYP[O_UD]}LB_9)SO(Q5.jpg" ALIGN=BOTTOM WIDTH=557 HEIGHT=260 BORDER=0></P>
-
<P ALIGN=LEFT> <FONT FACE="Arial Unicode MS, sans-serif"><SPAN LANG="en-US">Depicting
+
<P ALIGN=LEFT><FONT COLOR="#000000"><FONT FACE="Arial Unicode MS, sans-serif"><SPAN LANG="en-US">Depicting
all those critical points, as shown in the figure, the system could
all those critical points, as shown in the figure, the system could
oscillate when cellB&rsquo;s Hiill parameters were located in the
oscillate when cellB&rsquo;s Hiill parameters were located in the
<FONT COLOR="#00b0f0"><B>area</B></FONT> named <FONT COLOR="#00b0f0">&lsquo;</FONT><FONT COLOR="#00b0f0"><I><B>Bistable</B></I></FONT><FONT COLOR="#00b0f0"><B>&rsquo;</B></FONT>.</SPAN></FONT></P>
<FONT COLOR="#00b0f0"><B>area</B></FONT> named <FONT COLOR="#00b0f0">&lsquo;</FONT><FONT COLOR="#00b0f0"><I><B>Bistable</B></I></FONT><FONT COLOR="#00b0f0"><B>&rsquo;</B></FONT>.</SPAN></FONT></P>
<P ALIGN=CENTER><IMG SRC="https://static.igem.org/mediawiki/2011/1/18/006.png" NAME="图形7" ALT="说明: D:\b2.bmp" ALIGN=BOTTOM WIDTH=438 HEIGHT=329 BORDER=0></P>
<P ALIGN=CENTER><IMG SRC="https://static.igem.org/mediawiki/2011/1/18/006.png" NAME="图形7" ALT="说明: D:\b2.bmp" ALIGN=BOTTOM WIDTH=438 HEIGHT=329 BORDER=0></P>
-
<P ALIGN=LEFT> <FONT FACE="Arial Unicode MS, sans-serif"><SPAN LANG="en-US">By
+
<P ALIGN=LEFT><FONT COLOR="#000000"><FONT FACE="Arial Unicode MS, sans-serif"><SPAN LANG="en-US">By
adjusting certain parameters, we saw that the oscillation&rsquo;s
adjusting certain parameters, we saw that the oscillation&rsquo;s
period and phase could be controlled properly, which is the most
period and phase could be controlled properly, which is the most

Revision as of 16:58, 4 October 2011

JavaScript Slideshow - TinySlideshow

 

Introduction to Model Download the full text

&nbsp&nbspIn our project, we designed a quorum-sensing oscillator which consists of two types of cells. The expression of the reporter genes (GFP of one cell type and GFP of another) of the cells of the same type can fluctuate synchronously and certain designs were made to adjust the phase and the period of oscillation.

To understand the property of our system, we built a mathematical model based on ODEs (Ordinary Differential Equations) and DDEs (Delayed Differential Equations) to model and characterize this system. The simulation results helped us to deepen into further characteristics of the system.



Original Full Model Read more

Firstly we wanted to describe the system thoroughly without leaving out any seemingly unimportant actions and factors. As a result, the description of the system contains every possible mass actions as well as some hill kinetics, Henri-Michaelis-Menten kinetics, and the parameters were got from literature. The model was represented and simulated in the Matlab toolbox SIMBIOLOGY. We listed all 19 ODEs in the attached pdf file, you can see more details there.



Simplified DDE Model Read more

The original model contains too many factors for analyzing the general property of system. To understand the essential characters of the oscillator, we simplify the original model according to certain appropriate assumptions, like Quasi-equilibrium for fast reactions.

After series of derivation based on those assumptions (see the attached pdf file), we came up with the following set of DDEs (Delay Differential Equations)

说明: C:\Users\Harry\AppData\Roaming\Tencent\Users\806995810\QQ\WinTemp\RichOle\RH%W~L%6JQ)ZX06RH76VYV3.jpg

We coded the system by DDE description in MATLAB and did simulation analysis accordingly. The result showed that the system could oscillate under certain parameters .

说明: D:\b1.bmp

To further understand what parameters could make the system oscillate, we did bifurcation analysis on the Hill parameters. What we had to do was find the critical points where the system can nearly oscillate but a little disruption may lead to a steady state like that:

说明: C:\Users\Harry\AppData\Roaming\Tencent\Users\806995810\QQ\WinTemp\RichOle\)D6IYP[O_UD]}LB_9)SO(Q5.jpg

Depicting all those critical points, as shown in the figure, the system could oscillate when cellB’s Hiill parameters were located in the area named Bistable.

说明: D:\b2.bmp

By adjusting certain parameters, we saw that the oscillation’s period and phase could be controlled properly, which is the most impressive character of our system. Here we present a figure that the oscillation phase was adjusted by adding araC, which could induce the pBad promoter, in cell type B. After adding araC to our system at certain time, the oscillation was interrupted at beginning, but could gradually recover and finally, the phase was changed.



Dimensionless model Read more

In order to make a further analysis on stability of the system, and sensitivity of parameters, we further simplified the model to make them dimensionless. In addition, we tried to introduce feedback to our system and made a brief analysis on different types of feedback we introduced.



Quorum Sensing Effect Read more

What we have analyzed so far is focused on two-cell oscillation. Quorum-sensing oscillator is not simply a matter of expansion in magnitude, but a matter of robustness in allowing difference of each individual cell. Moreover, we test the adjustment of phase and period of oscillation in this part.

As we all know, no two things in this world are exactly the same, so do cells. The major differences between individual cells that we take into consideration include:

Each cell’s activity of promoter is varied, so each cell has different rate to generate AHL.

The initial amount of AHL may be disproportionally distributed among cells.

The rate of generating AHL is closely related to parameter m and n. Therefore, we introduce randomness to both parameters by letting them obey normal distribution, that is:

m(i) = μ1+N(0,σ1);

n(i) = μ2+ N(0,σ2);

μ1 and μ2 are the average ability of generating 30C6HSL and 3012CHSL, and normal distribution-- N(0,σ)--describes the fluctuations of AHL generating rate in individual cell. We then expanded our equations from 2 cells to a population of cells. Each cell share a mutual environment in which we assume that AHL in environment is proportionally distributed.

说明: D:\b1.bmp

The figures indicate that our system can oscillate synchronically being able to tolerate differences at certain range among a population of cells.

We also tested whether the oscillation is dependent on initial distribution of AHL by changing the initial amount drastically by letting them follow uniform distribution. That is:

Initial(i) = U(0,20);

Based on this distribution restraining the initial AHL concentration in each cell, we simulated out a figure as follows.

说明: D:\b1.bmp

The results demonstratively give evidence proving that our system can start to oscillate synchronically given variant initial starting status.



References

[1] Uri Alon, (2007). Network motifs: theory and experimental approaches. Nature.

[2] Chunbo Lou, Xili Liu, Ming Ni, et al. (2010). Synthesizing a novel genetic sequential logic circuit: a push-on push-off switch. Molecular Systems Biology.

[3] Tal Danino, Octavio Mondragon-Palomino, Lev Tsimring & Jeff Hasty (2010). A synchronized quorum of genetic clocks. Nature.

[4] Marcel Tigges, Tatiana T. Marquez-Lago, Jorg Stelling & Martin Fussenegger (2009). A tunable synthetic mammalian oscillator. Nature.

[5] Sergi Regot, Javio Macia el al. (2010). Distributed biological computation with multicellular engineered networks. Nature.

[6] Martin Fussenegger, (2010). Synchronized bacterial clocks. Nature.

[7] Andrew H Babiskin and Christina D Smolke, (2011). A synthetic library of RNA control modules for predictable tuning of gene expression in yeast. Molecular Systems Biology.

[8] Santhosh Palani and Casim A Sarkar, (2011). Synthetic conversion of a graded receptor signal into a tunable, reversible switch. Molecular Systems Biology.

[9] Nancy Kopell, (2002). Synchronizing genetic relaxation oscillation by intercell signaling. PNS