Team:Tsinghua-A/Modeling

From 2011.igem.org

(Difference between revisions)
(Prototype team page)
(Undo revision 249996 by Xiepei (talk))
 
(85 intermediate revisions not shown)
Line 1: Line 1:
-
No Page Text
+
{{:Team:Tsinghua-A/template/killbanner}}
 +
<html xmlns="http://www.w3.org/1999/xhtml">
 +
<head>
 +
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1" />
 +
<title>JavaScript Slideshow - TinySlideshow</title>
 +
<link rel="stylesheet" href="http://bioinfo.au.tsinghua.edu.cn/member/IGEM/thu_style.css" />
 +
<style type="text/css">
 +
 
 +
body {
 +
margin: 0;
 +
padding: 0;
 +
background: #F7F7F7 url("https://static.igem.org/mediawiki/2011/2/27/Thua-back1.jpg");
 +
background-repeat: no-repeat;
 +
background-attachment: fixed;
 +
font-family: Arial, Helvetica, sans-serif;
 +
font-size: 12px;
 +
color: #787878;
 +
}
 +
div.slider{
 +
float:right;
 +
}
 +
div.nav
 +
{
 +
width:180px;
 +
height:188px;
 +
float:left;
 +
border-style: none;
 +
border-width: 50px;
 +
margin : 0px 0px 0px 0px;
 +
}
 +
div.back{
 +
width:750px;
 +
height:175px;
 +
float:right;
 +
border-style: none;
 +
border-width: 50px;
 +
margin : 0px 0px 0px 0px;
 +
}
 +
div.imgbox{
 +
width:500px;
 +
height:320px;
 +
float:right;
 +
border-style: none;
 +
border-width: 50px;
 +
margin : 10px 10px 15px 10px;
 +
text-align: center;
 +
}
 +
div.imgbox2{
 +
width:450px;
 +
height:320px;
 +
float:left;
 +
border-style: none;
 +
border-width: 50px;
 +
margin : 10px 25px 15px 0px;
 +
text-align: center;
 +
 +
}
 +
div.imgbox3{
 +
width:540px;
 +
height:320px;
 +
float:left;
 +
border-style: none;
 +
border-width: 50px;
 +
margin : 10px 10px 10px 10px;
 +
text-align: center;
 +
}
 +
div.imgbox4{
 +
width:370px;
 +
height:320px;
 +
float:left;
 +
border-style: none;
 +
border-width: 50px;
 +
margin : 10px 10px 15px 10px;
 +
text-align: center;
 +
 +
 
 +
}
 +
div.imgbox5{
 +
width:400px;
 +
height:320px;
 +
float:left;
 +
border-style: none;
 +
border-width: 50px;
 +
margin : 10px 10px 15px 10px;
 +
text-align: center;
 +
position:relative;
 +
top:40px;
 +
}
 +
div.temp{
 +
position:relative;
 +
top:60px
 +
}
 +
h1 {
 +
font-family:"Times New Roman",Georgia,Serif;
 +
color:#548dd4;
 +
font-weight:bold;
 +
font-size:250%;
 +
padding:0;
 +
border:none;
 +
margin:0;
 +
}
 +
p{
 +
text-indent:2em;
 +
line-height:2;
 +
font-family:"Times New Roman",Georgia,Serif;
 +
color:#000000;
 +
font-size:110%
 +
}
 +
p.cite{
 +
line-height:2;
 +
font-family:"Times New Roman",Georgia,Serif;
 +
font-weight:bold;
 +
color:#000000;
 +
font-size:90%
 +
}
 +
h2{
 +
line-height:2;
 +
font-family:"Times New Roman",Georgia,Serif;
 +
color:#005555;
 +
font-size:100%;
 +
}
 +
head{
 +
margin:0;
 +
padding:0;
 +
//background: #F7F7F7 url("https://static.igem.org/mediawiki/2011/e/e2/Thua-back.jpg") repeat-x left top;
 +
}
 +
 
 +
/* Menu */
 +
 
 +
#menu {
 +
width: 940px;
 +
height: 49px;
 +
margin: 0 auto;
 +
padding: 0;
 +
background: #1c1c1c;
 +
}
 +
 
 +
#menu ul {
 +
margin: 0;
 +
padding: 0;
 +
list-style: none;
 +
line-height: normal;
 +
}
 +
 
 +
#menu li {
 +
float: left;
 +
border-left: 1px solid #FFFFFF;
 +
}
 +
 
 +
#menu a {
 +
display: block;
 +
width: 133px;
 +
height: 33px;
 +
padding-top: 16px;
 +
text-decoration: none;
 +
text-align: center;
 +
font-family: Arial, Helvetica, sans-serif;
 +
font-size: 13px;
 +
font-weight: bold;
 +
color: #f7f7f7;
 +
border: none;
 +
}
 +
 
 +
#menu a:hover, #menu .current_page_item a {
 +
background: #727272;
 +
text-decoration: none;
 +
}
 +
 
 +
#menu .current_page_item a {
 +
background: #727272;
 +
padding-left: 0;
 +
}
 +
#menu .normal_item a {
 +
background: #1c1c1c;
 +
padding-left: 0;
 +
}
 +
</style>
 +
 
 +
 
 +
</head>
 +
<body>
 +
<div style="position:fixed;right:0px;bottom: 0px;display:inline"><a href="#back-top">
 +
<img src="https://static.igem.org/mediawiki/2011/5/52/Back-top.png" width="60px" height="60px"></a></div>
 +
<div style="position:fixed;right:80px;bottom: 0px;display:inline"><a href="https://2011.igem.org/Team:Tsinghua-A/Homepage">
 +
<img src="https://static.igem.org/mediawiki/2011/8/82/Back-home.png" width="60px" height="60px"></a></div>
 +
 
 +
<p id="back-top" STYLE="text-indent:0em" align="CENTER"><img src="https://static.igem.org/mediawiki/2011/f/f9/THU_banner_mod.jpg" alt="" width="960"/></p>
 +
<div id="menu">
 +
<ul>
 +
<li ><a href="https://2011.igem.org/Team:Tsinghua-A/Homepage">Home</a></li>
 +
<li ><a href="https://2011.igem.org/Team:Tsinghua-A/Team">Team</a></li>
 +
<li ><a href="https://2011.igem.org/Team:Tsinghua-A/Project">Project</a></li>
 +
<li class="current_page_item"><a href="https://2011.igem.org/Team:Tsinghua-A/Modeling">Modeling</a></li>
 +
<li ><a href="https://2011.igem.org/Team:Tsinghua-A/Wetlab">Wetlab</a></li>
 +
<li><a href="https://2011.igem.org/Team:Tsinghua-A/Safety">Safety</a></li>
 +
<li ><a href="https://2011.igem.org/Team:Tsinghua-A/Notebook">Notebook</a></li>
 +
</ul>
 +
</div>
 +
<br>
 +
<!-- end #menu -->
 +
 
 +
 
 +
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
 +
<HTML>
 +
<HEAD>
 +
<META HTTP-EQUIV="CONTENT-TYPE" CONTENT="text/html; charset=gb2312">
 +
<TITLE></TITLE>
 +
<STYLE TYPE="text/css">
 +
<!--
 +
@page { margin: 2cm }
 +
P { margin-bottom: 0.21cm }
 +
A:link { so-language: zxx }
 +
-->
 +
</STYLE>
 +
</HEAD>
 +
<BODY LANG="zh-CN" DIR="LTR">
 +
 
 +
<div style="position:relative;left:30px;display:inline"><A HREF="https://2011.igem.org/Team:Tsinghua-A/Modeling"><img src="https://static.igem.org/mediawiki/2011/b/b7/ThuA_A2.png" width="150px" height="100px"></A></div>
 +
 
 +
<div style="position:relative;left:50px;display:inline"><A HREF="https://2011.igem.org/Team:Tsinghua-A/Modeling/P1A"><img src="https://static.igem.org/mediawiki/2011/0/0d/ThuA_B1.png" width="110px" height="100px"></A></div>
 +
 
 +
<div style="position:relative;left:70px;display:inline"><A HREF="https://2011.igem.org/Team:Tsinghua-A/Modeling/P2A"><img src="https://static.igem.org/mediawiki/2011/b/b6/ThuA_C1.png" width="110px" height="60px"></A></div>
 +
 
 +
<div style="position:relative;left:90px;display:inline"><A HREF="https://2011.igem.org/Team:Tsinghua-A/Modeling/P3A"><img src="https://static.igem.org/mediawiki/2011/6/6e/ThuA_D1.png" width="110px" height="60px"></A></div>
 +
 
 +
<div style="position:relative;left:110px;display:inline"><A HREF="https://2011.igem.org/Team:Tsinghua-A/Modeling/P4"><img src="https://static.igem.org/mediawiki/2011/6/6a/ThuA_E1.png" width="110px" height="120px"></A></div>
 +
 
 +
<div class="slider" style="padding-buttom:20px">
 +
<p><IMG SRC="https://static.igem.org/mediawiki/2011/1/13/ModelBack.jpg" ALIGN=bottom WIDTH=750 HEIGHT=175 BORDER=0 ISMAP></div>
 +
<div class="nav">
 +
<table id="toc" class="toc">
 +
<tr>
 +
<td>
 +
<div id="toctitle">
 +
<h2>Contents</h2>
 +
</div>
 +
<ul>
 +
<li class="toclevel-1 tocsection-1"><a href="#Introduction"><span class="tocnumber">1</span> <span class="toctext">Introduction to Model</span></a></li>
 +
<li class="toclevel-1 tocsection-2"><a href="#Original Full Model"><span class="tocnumber">2</span> <span class="toctext">Original Full Model</span></a></li>
 +
<li class="toclevel-1 tocsection-3"><a href="#Simplified DDE Model"><span class="tocnumber">3</span> <span class="toctext">Simplified DDE Model</span></a></li>
 +
<li class="toclevel-1 tocsection-4"><a href="#Dimensionless Model"><span class="tocnumber">4</span> <span class="toctext">Dimensionless Model</span></a></li>
 +
<li class="toclevel-1 tocsection-5"><a href="#Quorum Sensing"><span class="tocnumber">5</span> <span class="toctext">Quorum Sensing Effect</span></a></li>
 +
<li class="toclevel-1 tocsection-6"><a href="#References"><span class="tocnumber">6</span> <span class="toctext">References</span></a></li>
 +
</ul>
 +
</td>
 +
</tr>
 +
</table>
 +
</div>
 +
<br>
 +
<div class="slider">
 +
<p><IMG SRC="https://static.igem.org/mediawiki/2011/2/2b/000.png" NAME="graph1" ALIGN=bottom WIDTH=20 HEIGHT=20 BORDER=0 ISMAP><A HREF="https://static.igem.org/mediawiki/2011/9/9a/Modeling_Wiki.pdf"><U><I>Download the full text </I></U></A><IMG SRC="https://static.igem.org/mediawiki/2011/0/08/Thu_matlab.png" NAME="graph2" ALIGN=BOTTOM WIDTH=20 HEIGHT=20 BORDER=0 ISMAP><A HREF="https://static.igem.org/mediawiki/2011/c/c6/Thu_A_Matlab_Code.zip"><U><I>Download the source code package(Matlab)</I></U></A></P></div>
 +
<h1 id="Introduction">Introduction to Model</h1><hr width="100%" size=2 color=gray>
 +
 
 +
<div class="imgbox">
 +
<img class="border" style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/2011/0/03/001.png" width = "430px" height="148px"/>
 +
<img class="border" style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/2011/b/b0/002.png" width = "430px" height="148px" />
 +
<p class="cite">Designed gene circuit in cell A and cell B</p>
 +
</div>
 +
 
 +
<p>In 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.
 +
</p>
 +
 
 +
<br><br>
 +
<div class="slider">
 +
<p><A HREF="https://2011.igem.org/Team:Tsinghua-A/Modeling/P1A"><U><I>Read
 +
more</I></U></A></P>
 +
</div>
 +
<P id="Original Full Model"><h1>Original Full Model</h1></P>
 +
<hr width="100%" size=2 color=gray>
 +
 
 +
<P>Firstly we
 +
described 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, but too many parameters make
 +
it difficult to do further analyses, So here we only listed all 19 ODEs
 +
and a reletive parameter table( see <A HREF="https://static.igem.org/mediawiki/2011/9/9a/Modeling_Wiki.pdf">attached pdf file</A>).
 +
</P>
 +
 
 +
<br>
 +
<div class="slider">
 +
<p><A HREF="https://2011.igem.org/Team:Tsinghua-A/Modeling/P2A"><U><I>Read
 +
more</I></U></A></P>
 +
</div>
 +
<P id="Simplified DDE Model"><h1>Simplified DDE</h1></P>
 +
<hr width="100%" size=2 color=gray>
 +
 
 +
<div class="imgbox">
 +
<img class="border" style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/2011/b/b0/003.png" width = "440px" height="310px"/>
 +
<p class="cite">Simplified DDEs</p>
 +
</div>
 +
 
 +
<P>
 +
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.</P>
 +
<P>After series of derivation based on those assumptions, we came
 +
up with the following set of DDEs (Delay Differential Equations)
 +
which contains only 6 equations, see the figure right. And it
 +
would be much more convenient for us to do some neccessary analyses
 +
and research into the mathematical essence of our model.
 +
</P>
 +
 
 +
<div class="imgbox2">
 +
<img class="border" style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/2011/d/d5/004.png" width = "440px" height="310px"/>
 +
<p class="cite">Figure shows all variables are oscillating</p>
 +
</div>
 +
 
 +
<P>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.</P>
 +
<P>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.</P>
 +
<P>Depicting
 +
all those critical points, as shown in the figure, the system could
 +
oscillate when cellB's Hill parameters were located in the
 +
area named <FONT COLOR="#00b0f0"></FONT><FONT COLOR="#00b0f0"><I><B>Bistable</B></I></FONT><FONT COLOR="#00b0f0"></B></I>.</P>
 +
 
 +
<div class="imgbox3">
 +
<img class="border" style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/2011/2/20/005.png" width = "540px" height="300px"/>
 +
<p class="cite">The figure right is the phase trajectory of two signal molecules in environment,
 +
the more it looks like a circle, the more steadily our system will oscillate.</p>
 +
</div>
 +
<div class="imgbox4">
 +
<img class="border" style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/2011/1/18/006.png" width = "370px" height="300px"/>
 +
<p class="cite">Our system oscillates when parameters were selected in the area named 'Bistable'</p>
 +
</div>
 +
 
 +
<div class="imgbox5">
 +
<img class="border" style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/2011/e/e8/007.png" width = "370px" height="280px"/>
 +
</div>
 +
<br><br>
 +
<div class="temp"><P>By
 +
adjusting certain parameters, we saw that the oscillation&rsquo;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.</P></div>
 +
<br><br><br><br><br><br>
 +
 
 +
<div class="slider">
 +
<p><A HREF="https://2011.igem.org/Team:Tsinghua-A/Modeling/P3A"><U><I>Read
 +
more</I></U></A></P>
 +
</div>
 +
 
 +
<P id="Dimensionless Model"><h1>Dimensionless Model</h1></P>
 +
<hr width="100%" size=2 color=gray>
 +
 
 +
<P>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. Some analyses were similar to the simplified
 +
DDE model, and you can see more details by clicking
 +
<A HREF="https://2011.igem.org/Team:Tsinghua-A/Modeling/P4">read
 +
more</A></P>
 +
 
 +
<br>
 +
<div class="slider">
 +
<p><A HREF="https://2011.igem.org/Team:Tsinghua-A/Modeling/P4"><U><I>Read
 +
more</I></U></A></P>
 +
</div>
 +
<P id="Quorum Sensing"><h1>Quorum Sensing Effect</h1></P>
 +
<hr width="100%" size=2 color=gray>
 +
 
 +
<P>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.</P>
 +
<P>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:</P>
 +
<P><B>●Each cell's activity of promoter is varied, so each cell has
 +
different rate to generate AHL.</B></P>
 +
<P><B>●The initial amount of AHL may be disproportionally distributed among
 +
cells.</B></P>
 +
<P>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:
 +
</P>
 +
<P ALIGN=CENTER style="text-intent:0em">m(i)= &mu;1+<I>N</I>(0,&sigma;1);</P>
 +
<P ALIGN=CENTER style="text-intent:0em">n(i)= &mu;2+<I>N</I>(0,&sigma;2);</P>
 +
<P ALIGN=LEFT><A NAME="OLE_LINK63"></A><A NAME="OLE_LINK62"></A>&mu;<SPAN LANG="en-US">1
 +
and &mu;2 are the average ability of generating 30C6HSL and 3012CHSL,
 +
and normal distribution-- </SPAN><SPAN LANG="en-US"><I>N</I></SPAN><SPAN LANG="en-US">(0,&sigma;)--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.</SPAN></SPAN></P>
 +
<P ALIGN="center" style="text-intent:0em"><IMG SRC="https://static.igem.org/mediawiki/2011/1/17/008.png" WIDTH=800 HEIGHT=600 BORDER=0></P>
 +
<P>The
 +
figures indicate that our system can oscillate synchronically being
 +
able to tolerate differences at certain range among a population of
 +
cells.</P>
 +
<P>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:</P>
 +
<P ALIGN="center" style="text-intent:0em">Initial(i)= <I>U</I>(0,20);</P>
 +
<P>Based
 +
on this distribution restraining the initial AHL concentration in
 +
each cell, we simulated out a figure as follows.</P>
 +
<P ALIGN="center" style="text-intent:0em"><IMG SRC="https://static.igem.org/mediawiki/2011/d/d2/009.png" WIDTH=800 HEIGHT=600 BORDER=0></P>
 +
<P>The
 +
results demonstratively give evidence proving that our system can
 +
start to oscillate synchronically given variant initial starting
 +
status.</P>
 +
<br>
 +
 
 +
<P id="References"><h1>References</h1></P>
 +
<hr width="100%" size=2 color=gray>
 +
<P>[1]
 +
Uri Alon, (2007). Network motifs: theory and experimental approaches.
 +
Nature.</P>
 +
<P>[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.</P>
 +
<P>[3]
 +
Tal Danino, Octavio Mondragon-Palomino, Lev Tsimring &amp; Jeff Hasty
 +
(2010). A synchronized quorum of genetic clocks. Nature.</P>
 +
<P>[4]
 +
Marcel Tigges, Tatiana T. Marquez-Lago, Jorg Stelling &amp; Martin
 +
Fussenegger (2009). A tunable synthetic mammalian oscillator. Nature.</P>
 +
<P>[5]
 +
Sergi Regot, Javio Macia el al. (2010). Distributed biological
 +
computation with multicellular engineered networks. Nature.</P>
 +
<P>[6]
 +
Martin Fussenegger, (2010). Synchronized bacterial clocks. Nature.</P>
 +
<P>[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.</P>
 +
<P>[8]
 +
Santhosh Palani and Casim A Sarkar, (2011). Synthetic conversion of a
 +
graded receptor signal into a tunable, reversible switch. Molecular
 +
Systems Biology.</P>
 +
<P>[9]
 +
Nancy Kopell, (2002). Synchronizing genetic relaxation oscillation by
 +
intercell signaling. PNS</P>
 +
 
 +
<br><br><br><br>
 +
<p style="text-indent:0em" align="CENTER"><a href="https://2011.igem.org/Team:Tsinghua-A"><img src="https://static.igem.org/mediawiki/2011/9/92/Killbanner_header.jpg" alt="" width="960"/><a href="https://2011.igem.org"><img src="https://static.igem.org/mediawiki/igem.org/2/29/Killbanner_header2.jpg" alt="" width="960"/></p>
 +
 
 +
</BODY>
 +
</HTML>

Latest revision as of 05:20, 28 October 2011

JavaScript Slideshow - TinySlideshow



Introduction to Model


Designed gene circuit in cell A and cell B

In 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


Firstly we described 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, but too many parameters make it difficult to do further analyses, So here we only listed all 19 ODEs and a reletive parameter table( see attached pdf file).


Simplified DDE


Simplified DDEs

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, we came up with the following set of DDEs (Delay Differential Equations) which contains only 6 equations, see the figure right. And it would be much more convenient for us to do some neccessary analyses and research into the mathematical essence of our model.

Figure shows all variables are oscillating

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.

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.

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

The figure right is the phase trajectory of two signal molecules in environment, the more it looks like a circle, the more steadily our system will oscillate.

Our system oscillates when parameters were selected in the area named 'Bistable'



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


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. Some analyses were similar to the simplified DDE model, and you can see more details by clicking read more


Quorum Sensing Effect


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.

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.

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