Team:Tsinghua-A/Modeling/P1A

From 2011.igem.org

(Difference between revisions)
Line 26: Line 26:
}
}
div.imgbox{
div.imgbox{
-
width:1000px;
+
width:500px;
-
height:800px;
+
height:360px;
-
border-style: solid;  
+
float:right;
-
border-width: 1px;  
+
border-style: none;  
 +
border-width: 50px;
 +
margin : 30px 0px 10px 10px;
text-align: center;
text-align: center;
-
padding:100px,100px,300px,100px;
 
}
}
div.imgbox2{
div.imgbox2{
-
width:450px;
+
width:710px;
-
height:320px;
+
height:1200px;
-
float:left;
+
border-style: none;  
border-style: none;  
-
border-width: 50px;
+
margin : 0px 0px 0px 0px;
-
margin : 10px 25px 15px 0px;
+
text-align: center;
text-align: center;
}
}
div.imgbox3{
div.imgbox3{
-
width:540px;
+
width:750px;
-
height:320px;
+
height:1800px;
-
float:left;
+
border-style: none;  
border-style: none;  
-
border-width: 50px;
 
-
margin : 10px 10px 10px 10px;
 
text-align: center;
text-align: center;
}
}
Line 100: Line 96:
font-size:100%;
font-size:100%;
}
}
 +
 +
/* Menu */
/* Menu */
Line 161: Line 159:
<div id="wrapper">
<div id="wrapper">
   <div id="header">
   <div id="header">
-
 
+
<p align="left"><img src="https://static.igem.org/mediawiki/2011/f/f9/THU_banner_mod.jpg" alt="" width="960"/></p></div>
-
 
+
-
<p 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>
+
<!-- end #header -->
<!-- end #header -->
<div id="menu">
<div id="menu">
Line 177: Line 171:
</ul>
</ul>
</div>
</div>
-
<p>&nbsp;</p>
+
<p>&nbsp;</p>
-
<!-- end #menu -->
+
<!-- end #menu --><!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
-
 
+
-
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
+
<HTML>
<HTML>
<HEAD>
<HEAD>
Line 186: Line 178:
</HEAD>
</HEAD>
<BODY LANG="zh-CN" DIR="LTR">
<BODY LANG="zh-CN" DIR="LTR">
-
<P><FONT COLOR="#000000"><FONT FACE="Times New Roman, serif"><SPAN LANG="en-US"><FONT FACE="Arial, sans-serif"><FONT SIZE=5 STYLE="font-size: 20pt"><SPAN LANG="en-US"><A HREF="https://2011.igem.org/Team:Tsinghua-A/Modeling">Modeling</A></SPAN></FONT></FONT></FONT><FONT FACE="Arial, sans-serif"><FONT SIZE=5 STYLE="font-size: 20pt"><SPAN LANG="en-US">::</SPAN></FONT></FONT><FONT COLOR="#666600"><FONT FACE="Arial, sans-serif"><FONT SIZE=5 STYLE="font-size: 20pt"><SPAN LANG="en-US">Quorum-sensing Effect</SPAN></FONT></FONT></FONT></SPAN></FONT></P><br>
+
<P><FONT FACE="Times New Roman, serif"><SPAN LANG="en-US"><FONT FACE="Arial, sans-serif"><FONT SIZE=5 STYLE="font-size: 20pt"><SPAN LANG="en-US"><A HREF="https://2011.igem.org/Team:Tsinghua-A/Modeling">Modeling</A>:: </SPAN></FONT></FONT><FONT COLOR="#355e00"><FONT FACE="Arial, sans-serif"><FONT SIZE=5 STYLE="font-size: 20pt"><SPAN LANG="en-US">Original Full Model</SPAN></FONT></FONT></FONT></SPAN></FONT></P>
-
<p>What we have done insofar 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 amplitude of oscillation in this part.</p>
+
<table id="toc" class="toc">
-
<p>As we all know, no two things in this world are the same, so do cells. The major difference of individual cell that we take into considerations is twofold:</p>
+
<tr>
-
<P><B>●Each cell's activity of promoter is varied, so each cell has
+
<td>
-
different rate to generate AHL.</B></P>
+
<div id="toctitle">
-
<P><B>●The initial amount of AHL may be disproportionally distributed among
+
<h2>Contents</h2>
-
cells.</B></P>
+
</div>
-
<P>The rate of generating AHL is closely related to parameter m and n.
+
<ul>
-
Therefore, we introduce randomness to both parameters by letting them
+
<li class="toclevel-1 tocsection-1"><a href="#Construction"><span class="tocnumber">1</span> <span class="toctext">Construction of ODE equation</span></a></li>
-
obey normal distribution, that is:
+
<li class="toclevel-1 tocsection-2"><a href="#Parameters"><span class="tocnumber">2</span> <span class="toctext">Parameters</span></a></li>
-
</P>
+
<li class="toclevel-1 tocsection-3"><a href="#Results"><span class="tocnumber">3</span> <span class="toctext">Results</span></a></li>
-
<P ALIGN=CENTER>m(i)= &mu;1+<I>N</I>(0,&sigma;1);</P>
+
</ul>
-
<P ALIGN=CENTER>n(i)= &mu;2+<I>N</I>(0,&sigma;2);</P>
+
</td>
-
<P>&mu;1 and &mu;2 are the average ability of generating 30C6HSL and 3012CHSL,
+
</tr>
-
and normal distribution--(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.</P>
+
</table>
 +
<br>
 +
<h1 id="Construction">Construction of ODE equation</h1><hr width="100%" size=2 color=gray>
 +
<p>At our first step, 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. We came up a set of ODEs with 19 equations.</p>
-
<p align="center" style="text-indent:0em"><img src="https://static.igem.org/mediawiki/igem.org/0/00/Part4-1.png" width = "900px" height="675px"></p>
+
<div class="imgbox">
-
<p align="center" style="text-indent:0em"><b>Figure 16 100 Cells Varied in parameter m and n</b></p>
+
<img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/2011/0/03/001.png" width = "430px" height="150px"/>
-
<p>The figures indicate that our system can oscillate synchronically being able to tolerate differences among a population of cells. Furthermore, the figures prove that different ability of generating AHLs of cells have nothing to do with the period and phase of the oscillation. We can also see that the oscillation amplitude of each cell is to a greater extent varied when the Variance of interruption is enlarged.</p>
+
<p class="cite">Figure 1 designed circuit of cell I</p>
-
<p>Moreover, we test 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>
+
<img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/2011/b/b0/002.png" width = "430px" height="150px" />
-
<P ALIGN=CENTER>Initial(i)= <I>U</I>(0,20);</P>
+
<p class="cite">Figure 2 designed circuit of cell II</p>
-
<p>The results would give evidence to prove that our system can start to oscillate synchronically given variant initial starting numbers.</p>
+
</div>
-
<p>Based on this distribution restraining the initial AHL concentration in each cell, we simulated out a figure as follows.</p>
+
<p>Promoter 1 and promoter 2 preceding lasR and luxR genes respectively are constant promoters, which will transcribe and translate into protein PlasR and PluxR. LA1 is the binding association of lasR and 30C12HSL(A2C1) and it can affect the subsequent promoter 2 which can be described by Hill Equation. The same goes to LA2. Gene luxI will be translated into protein PluxI which would generate 30C6HSL(A1C1) through enzymatic reaction. The AHL will diffuse through the membrane to the environment(A1e) and finally enter into Cell 2(A1C2). Protein PtetR which is translated from gene tetR represses promoter 5 which is responsible for transcription of gene lasI. Promoter 6 is constant for translation of protein PlasI. 30C12HSL(A2C2) is generated from Protein PlasI through enzymatic reaction. 30C12HSL in the environment is called A2e which will diffuse to Cell 1. aTc is added manipulatively to change the phase of oscillation by binding the protein PTetR. Therefore, we have these following ODEs:<br><br></p>
-
 
+
<p align="center" style="text-intend:0em"><img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/9/9a/ThuAModel_1_1.png" width = "707px" height="1115px" /><br></p>
-
<p align="center" style="text-indent:0em"><img src="https://static.igem.org/mediawiki/igem.org/2/2d/Part4-2.png" width = "900px" height="675px"></p>
+
<h1 id="Parameters">Parameters</h1><hr width="100%" size=2 color=gray>
-
<p align="center" style="text-indent:0em"><b>Figure 17 100 Cells Varied in initial AHL concentration</b></p>
+
<p>The parameters are inherent factors determining the behaviors, properties of a system. We selected the quantities thoughtfully from previous iGEM teams and some others were found from published papers.</p>
-
<p>The results demonstratively give evidence proving that our system can start to oscillate synchronically given variant initial starting numbers.</p>
+
-
<br><br>
+
 +
<p align="center" style="text-intend:0em"><img style="border-color:#B2B2B2;"src="https://static.igem.org/mediawiki/igem.org/c/c8/ThuAModel_1_2.png" width = "746px" height="1794px" /><br></p>
 +
<p align="center" style="text-intend:0em"><b>Table 1 Parameters of ODEs</b></p>
 +
<br>
 +
<h1 id="Results">Results</h1><hr width="100%" size=2 color=gray>
 +
<p>We simulated this system by SIMBIOLOGY, a toolbox embedded in MATLAB. However, unaware of the key parameters to which the system is sensitive, we felt difficult to control or adjust properly, and the simulation result of the system came into a damped oscillation. We ascribed the inability of our model to the fact that the precise descriptions contain too many equations and parameters and we felt obliged to establish a simplified model in place of the precise one for simulation and further analysis.</p>
 +
<br><br><br>
<p align="CENTER" style="text-indent:0em"><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>
<p align="CENTER" style="text-indent:0em"><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>
</BODY>
</HTML>
</HTML>

Revision as of 18:47, 25 October 2011

Republic by Free CSS Templates

 

Modeling:: Original Full Model

Contents

Construction of ODE equation

At our first step, 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. We came up a set of ODEs with 19 equations.

Figure 1 designed circuit of cell I

Figure 2 designed circuit of cell II

Promoter 1 and promoter 2 preceding lasR and luxR genes respectively are constant promoters, which will transcribe and translate into protein PlasR and PluxR. LA1 is the binding association of lasR and 30C12HSL(A2C1) and it can affect the subsequent promoter 2 which can be described by Hill Equation. The same goes to LA2. Gene luxI will be translated into protein PluxI which would generate 30C6HSL(A1C1) through enzymatic reaction. The AHL will diffuse through the membrane to the environment(A1e) and finally enter into Cell 2(A1C2). Protein PtetR which is translated from gene tetR represses promoter 5 which is responsible for transcription of gene lasI. Promoter 6 is constant for translation of protein PlasI. 30C12HSL(A2C2) is generated from Protein PlasI through enzymatic reaction. 30C12HSL in the environment is called A2e which will diffuse to Cell 1. aTc is added manipulatively to change the phase of oscillation by binding the protein PTetR. Therefore, we have these following ODEs:


Parameters

The parameters are inherent factors determining the behaviors, properties of a system. We selected the quantities thoughtfully from previous iGEM teams and some others were found from published papers.


Table 1 Parameters of ODEs


Results

We simulated this system by SIMBIOLOGY, a toolbox embedded in MATLAB. However, unaware of the key parameters to which the system is sensitive, we felt difficult to control or adjust properly, and the simulation result of the system came into a damped oscillation. We ascribed the inability of our model to the fact that the precise descriptions contain too many equations and parameters and we felt obliged to establish a simplified model in place of the precise one for simulation and further analysis.




Retrieved from "http://2011.igem.org/Team:Tsinghua-A/Modeling/P1A"