Team:HUST-China/Modelling/MCOS

From 2013.igem.org

(Difference between revisions)
Line 69: Line 69:
<!------------------------------------------Left side navigator start---------------------------------------------------->
<!------------------------------------------Left side navigator start---------------------------------------------------->
          <li><a href="https://2013.igem.org/Team:HUST-China/Modelling">Overview</a></li>
          <li><a href="https://2013.igem.org/Team:HUST-China/Modelling">Overview</a></li>
-
          <li><a href="https://2013.igem.org/Team:HUST-China/Modelling/DDE_Model"></i>Delay Difference Equations</a></li>
+
          <li class="active"><a href="javascript:void;"></i>Delay Difference Equations</a></li>
-
                                         <li class="active"><a href="javascript:void;"></i>Multi Cells Oscillation Simulation</a></li>
+
                                         <li><a href="https://2013.igem.org/Team:HUST-China/Modelling/MCOS"></i>Multi Cells Oscillation Simulation</a></li>
<!------------------------------------------Left side navigator end---------------------------------------------------->
<!------------------------------------------Left side navigator end---------------------------------------------------->
            </ul>
            </ul>
Line 77: Line 77:
              <div id="scrollspy">
              <div id="scrollspy">
<!----------------------------------------------Main content start------------------------------------------------------>
<!----------------------------------------------Main content start------------------------------------------------------>
-
<div id="Part-1">
+
<div id="Part-2">
-
<h1 class="page-header"><strong>Overview</strong></h1>
+
<h1><strong>Multi Cells Oscillation Simulation</strong></h1>
</div>
</div>
-
<p>
+
<h2><strong>Goal</strong></h2>
-
The major goal of modeling this year is stimulating the working condition of the oscillation. The genetic oscillation we used is described as below.
+
In order to see if group would oscillate just like a single cell does.
-
<div style="float:right;">
+
<h2><strong>Methods</strong></h2>
-
<img src="https://static.igem.org/mediawiki/2013/a/ab/Pathway_HUST.png" width="350px;" />
+
1.establish the grow to death program according to our hypothesis;<br />
-
<p class="small">Fig 1.The pathway of genetic oscillator used in this paper</p
+
2.Investigate reasonable parameter sets from previous researches;<br />
-
></div>
+
3.simulation<br />
-
The enzymes that can induce the expression of propanoic acid is replaced with mRFP. Throughout the whole paper, we choose to study AraC instead of mRFP since they are in the same plasmid and expression rate of both protein is similar. By doing this, we can reduce the number of equations. We put our focus mainly on two aspect:
+
 
-
(1)Oscillator of a single cell (DDEs).
+
<h2><strong>Results</strong></h2>
-
(2)Oscillator of a group of cells (MCOS).
+
<img src="https://static.igem.org/mediawiki/2013/e/e6/HUST_Cells.png">
-
To be more specifically, we want to find out if this genetic oscillator can oscillate in the first place, if it is stable against environment change, if its period can be adjusted, whether a large group of these oscillators can oscillate as well.
+
<p class="small">Fig 1.Population of bacteria against time. Average = 105180, largest difference Δ = 5667</p><br />
 +
The population of bacteria is fluctuating within a small range (5667/105180 = 0.0539) and generally steady, showing that the logistic model is feasible and this model is successfully simulating the population within microencapsulation.
 +
 
 +
<img src="https://static.igem.org/mediawiki/2013/b/b0/HUST_AraC_multicell.png">
 +
<p class="small">Fig 2.AraC concentration of simulated multi oscillating cells within microencapsulation from 100000 minutes to 100800 minutes since the simulation started.</p><br />
 +
The multi cells oscillation simulation suggests that even with that amount of cells, the oscillation will still exist just like a single one does. Such result is because of the synchronous of all cells' oscillations throughout the whole process.
 +
 
 +
<h2><strong>Background</strong></h2>
 +
For convenience, we assume that all bacteria are in rapid growth. All bacteria would synthesize protein regardless of cell fission. In terms of death of cells, according to logistic model of cell population, death rate is linear to population itself. Based on this, we set that a certain amount of cells, which is proportion to squares of population, would be killed due to limited food. Those cells are picked randomly. We were also interested in lifespan of bacteria. While in rapid grow, the average lifespan of a bacterium is 20 minutes. Without further information, we assumed that the lifespan of bacterium obey a gaussian distribution with μ= 20 min, σ= 1 min.
 +
Moreover, we assumed that the quantity of AraC are linear with plasmid copies. Given a specific environment, the number of initial plasmid copies are a constant. The replication of plasmids can be thought to be completed instantly (within 0.05min). When cells started to divide, the plasmids are allocated into two filial cells evenly. During the lifespan of a single cell, the number of replicated plasmids obey Poisson distribution.<br/>
 +
Lastly, we assumed that when started simulation, all the cell are in the oscillation phase and remain exactly the same oscillation rhythm with each other.<br/>
 +
Based on these information, we set several rules for our computer model:<br/>
 +
(1)Each round of simulation represent 0.1min. <br/>
 +
(2)The lifespan of cells $x\sim (20, 1^2).$<br/>
 +
(3)When reached their lifespan, cells division. Division is completed immediately. <br/>
 +
(4)The expression of AraC would not be affected by division, which means the phase of AraC does not change.<br/>
 +
(5)Certain amount of cells, which is proportion to square of population, is 'sentenced' to dead in every round. They are picked randomly.<br/>
 +
(6)AraC's concentration is proportion to plasmid copies.<br/>
 +
(7)When cells are dividing, the plasmid copies would increase by y, and y~Pois(50), then split evenly into two filial cells.<br/>
 +
 
 +
<h2><strong>Colon-Specific Drug Delivery System </strong></h2>
 +
To simulate real-life situation, we researched for how our genetic oscillator will enter our bodies. As propionate is absorbed in human colon, so the best way to delivery those genetic engineered bacteria into colon. However, considered the possible safety issued brought by bacteria, we sought to solve it by adopting Oral Colon-Specific Drug Delivery System (OCDDS). Such system encapsulate bacteria with semipermeable polymer membrane. Patient take such microencapsulation orally and it will stay in colon for certain days. Due to the semipermeable property, small molecules such as propionate will penetrate the membrane through pores while bacteria are too large to pass through. Bacteria simply stays in microencapsulation and eventually excreted from human body. Such system can prevent bacteria from directly contacting human body and enhance the security.<br/>
 +
Modified microencapsulation can stay in colon for 70 days. Bacteria concentration in microencapsulation can reach 70 days and diameter of microencapsulation can reach 433 67μm. For convenience, we assumed that each monocolony originated from a single bacteria. Based on that, the numbers of bacteria are approximately 81713~208333. Since numbers of bacteria remain relatively constant within microencapsulation, the feasibility and stability of supposed therapy are guaranteed.
 +
 
 +
 
-
</p>
 
<!----------------------------------------------Main content end------------------------------------------------------>
<!----------------------------------------------Main content end------------------------------------------------------>
Line 125: Line 148:
     var scroll = $(document).scrollTop();
     var scroll = $(document).scrollTop();
     var navli = $("#navbar .nav li");
     var navli = $("#navbar .nav li");
-
+
/* var part1= $("#Part-1").outerHeight();
 +
                                var part2= $("#Part-2").outerHeight();
 +
                                var part3= $("#Part-3").outerHeight();
 +
var futurework= $("#Future-Work").outerHeight();
 +
                                var judge= $("#Judging-Criteria").outerHeight();*/
if(scroll > 454){
if(scroll > 454){
                                         $("#up").slideDown("ease-out");
                                         $("#up").slideDown("ease-out");
                                         $("#mainNavBar").addClass("navbar-fixed-top");
                                         $("#mainNavBar").addClass("navbar-fixed-top");
$(".bs-docs-sidenav").removeClass("affix-top").addClass("affix");
$(".bs-docs-sidenav").removeClass("affix-top").addClass("affix");
-
+
/* if (scroll < 454 + part1){
 +
navli.removeClass("active").eq(0).addClass("active");
 +
}else if (!(scroll < 454 + part1|| scroll > 454 + part1+ part2)){
 +
navli.removeClass("active").eq(1).addClass("active");
 +
}else if (!(scroll < 454 + part1 + part2 || scroll > 454 + part1+ part2 + part3)){
 +
navli.removeClass("active").eq(2).addClass("active");
 +
}else if (!(scroll < 454 + part1 + part2 + part3 || scroll > 454 + part1+ part2 + part3 + futurework)){
 +
navli.removeClass("active").eq(3).addClass("active");
 +
}else if (!(scroll < 454 + part1 + part2 + part3  + futurework || scroll > 454 + part1+ part2 + part3 + judge)){
 +
navli.removeClass("active").eq(4).addClass("active");
 +
}*/
}else{
}else{
                                         $("#up").slideUp("ease-in");
                                         $("#up").slideUp("ease-in");
                                         $("#mainNavBar").removeClass("navbar-fixed-top");
                                         $("#mainNavBar").removeClass("navbar-fixed-top");
$(".bs-docs-sidenav").removeClass("affix").addClass("affix-top");
$(".bs-docs-sidenav").removeClass("affix").addClass("affix-top");
-
+
//navli.removeClass("active").eq(0).addClass("active");
}
}
});
});

Revision as of 17:04, 27 September 2013

Multi Cells Oscillation Simulation

Goal

In order to see if group would oscillate just like a single cell does.

Methods

1.establish the grow to death program according to our hypothesis;
2.Investigate reasonable parameter sets from previous researches;
3.simulation

Results

Fig 1.Population of bacteria against time. Average = 105180, largest difference Δ = 5667


The population of bacteria is fluctuating within a small range (5667/105180 = 0.0539) and generally steady, showing that the logistic model is feasible and this model is successfully simulating the population within microencapsulation.

Fig 2.AraC concentration of simulated multi oscillating cells within microencapsulation from 100000 minutes to 100800 minutes since the simulation started.


The multi cells oscillation simulation suggests that even with that amount of cells, the oscillation will still exist just like a single one does. Such result is because of the synchronous of all cells' oscillations throughout the whole process.

Background

For convenience, we assume that all bacteria are in rapid growth. All bacteria would synthesize protein regardless of cell fission. In terms of death of cells, according to logistic model of cell population, death rate is linear to population itself. Based on this, we set that a certain amount of cells, which is proportion to squares of population, would be killed due to limited food. Those cells are picked randomly. We were also interested in lifespan of bacteria. While in rapid grow, the average lifespan of a bacterium is 20 minutes. Without further information, we assumed that the lifespan of bacterium obey a gaussian distribution with μ= 20 min, σ= 1 min. Moreover, we assumed that the quantity of AraC are linear with plasmid copies. Given a specific environment, the number of initial plasmid copies are a constant. The replication of plasmids can be thought to be completed instantly (within 0.05min). When cells started to divide, the plasmids are allocated into two filial cells evenly. During the lifespan of a single cell, the number of replicated plasmids obey Poisson distribution.
Lastly, we assumed that when started simulation, all the cell are in the oscillation phase and remain exactly the same oscillation rhythm with each other.
Based on these information, we set several rules for our computer model:
(1)Each round of simulation represent 0.1min.
(2)The lifespan of cells $x\sim (20, 1^2).$
(3)When reached their lifespan, cells division. Division is completed immediately.
(4)The expression of AraC would not be affected by division, which means the phase of AraC does not change.
(5)Certain amount of cells, which is proportion to square of population, is 'sentenced' to dead in every round. They are picked randomly.
(6)AraC's concentration is proportion to plasmid copies.
(7)When cells are dividing, the plasmid copies would increase by y, and y~Pois(50), then split evenly into two filial cells.

Colon-Specific Drug Delivery System

To simulate real-life situation, we researched for how our genetic oscillator will enter our bodies. As propionate is absorbed in human colon, so the best way to delivery those genetic engineered bacteria into colon. However, considered the possible safety issued brought by bacteria, we sought to solve it by adopting Oral Colon-Specific Drug Delivery System (OCDDS). Such system encapsulate bacteria with semipermeable polymer membrane. Patient take such microencapsulation orally and it will stay in colon for certain days. Due to the semipermeable property, small molecules such as propionate will penetrate the membrane through pores while bacteria are too large to pass through. Bacteria simply stays in microencapsulation and eventually excreted from human body. Such system can prevent bacteria from directly contacting human body and enhance the security.
Modified microencapsulation can stay in colon for 70 days. Bacteria concentration in microencapsulation can reach 70 days and diameter of microencapsulation can reach 433 67μm. For convenience, we assumed that each monocolony originated from a single bacteria. Based on that, the numbers of bacteria are approximately 81713~208333. Since numbers of bacteria remain relatively constant within microencapsulation, the feasibility and stability of supposed therapy are guaranteed.