Team:KAIST Korea/Project Modeling

From 2012.igem.org

(Difference between revisions)
 
(6 intermediate revisions not shown)
Line 139: Line 139:
<div id="kaistcontent">
<div id="kaistcontent">
<div>
<div>
-
<a href="#1st"><span id="tab"><b>1] Cell Growth Curve</b> </span></a>
+
<a href="#1st"><span id="tab"><b>Cell Growth Curve</b> </span></a>
-
<a href="#2nd"><span id="tab"><b>2] title??</b></span></a>
+
<a href="#2nd"><span id="tab"><b>Proof of concept</b></span></a>
-
<a href="#3rd"><span id="tab"><b>3] title??</b></span></a>
+
<a href="#3rd"><span id="tab"><b>Auto Regulation</b></span></a>
</br></br>
</br></br>
<img id="starter-grad" style="height:80px" src="https://static.igem.org/mediawiki/2012/9/95/Starter_gradient_kaist.png"></img>
<img id="starter-grad" style="height:80px" src="https://static.igem.org/mediawiki/2012/9/95/Starter_gradient_kaist.png"></img>
Line 151: Line 151:
<section id="1st">
<section id="1st">
<div>
<div>
-
<span id ="sub-sub-title">1] Cell Growth Curve</span></br></br>
+
<span id ="sub-sub-title">Cell Growth Curve</span></br></br>
-
<span id="little">Cell growth can be modeled using Logistic differential equation. The logistic equation is shown below.</span></br></br>
+
<span id="little">Cell growth can be modeled using Logistic differential equation as shown below.</span></br></br>
<div align="center" ><img src='https://static.igem.org/mediawiki/2012/a/a2/KAIST_Logistic_equation.png'/></div>
<div align="center" ><img src='https://static.igem.org/mediawiki/2012/a/a2/KAIST_Logistic_equation.png'/></div>
</br>
</br>
-
<span id="little">When we solve above equation with appropriate parameters(using MATLAB), we can get solution curve like below. This curve matches with our knowledge about cell growth.</span></br></br>
+
<span id="little">When we solve this equation with appropriate parameters(using MATLAB), we can get solution curve as shown below. This curve matches with our knowledge about cell growth.</span></br></br>
-
<div align="center" ><img src='#'/></div>
+
<div align="center" ><img src='https://static.igem.org/mediawiki/2012/4/4a/KAIST_Cell_growth_curve.png'/></div>
-
</br>
+
-
+
</div>
</div>
</section>
</section>
</br></br>
</br></br>
-
<div align='right'><a href="#top">▲Back to the top</a></div></br>
+
<div align='right'><a href="#top">▲ Back to the top</a></div></br>
</br>
</br>
<section id="2nd">
<section id="2nd">
<div>
<div>
-
<span id="sub-sub-title">2] Title???</span></br></br>
+
<span id="sub-sub-title">Proof of concept</span></br></br>
-
<span id="little">With GFP and RFP, we want to check whether our system really works or not. Since we cannot consider all of complicated biological phenomena, we assume our system simply follows reaction rate theory and mass balance equations. Below reactions are simplified version of our system, and we CONSIDER the production and degradation of mRNAs and Genes, but we IGNORE the polysome phenomena and any gene regulatory system that is occur in real biological system. Also, we CONSIDER that there is many copies of plasmids in <i>E.coli</i> and each plasmid can react with invertase to invert their gene sequence.</br>
+
<ul>
-
Because every plasmid have equal probability that react with invertase, we assume plasmids follow uniform distribution. Also assume each plasmids are mutually independent, that is each  plasmid cannot affect invertase reaction of other plasmid.</br></br>
+
<li style="list-style-type:square;font-size:14px;font-weight:bold;">Mathematical model</li>
 +
</ul>
 +
</br>
 +
<span id="little">With GFP and RFP, we want to check whether our system really works or not. Since we cannot consider all the complicated biological phenomena, we assume our system simply follows reaction rate theory and mass balance equations. Reactions below are simplified version of our system. We CONSIDERED the production and degradation of mRNAs and Genes, while we IGNORED the polysome phenomena and any gene regulatory system that occurs in real biological system. Also, we CONSIDERED that there are many copies of plasmids in <i>E.coli</i> and each plasmid can react with invertase to invert their gene sequence.</br>
 +
Because every plasmid has equal probability to react with invertase, we assume plasmids follow uniform distribution. Final assumption is that each plasmid is mutually independent, that is, each  plasmid cannot affect invertase reaction of other plasmid.</br></br></span>
<div align="center" ><img src='https://static.igem.org/mediawiki/2012/d/d7/KAIST_Gfprfp.png'/></div>
<div align="center" ><img src='https://static.igem.org/mediawiki/2012/d/d7/KAIST_Gfprfp.png'/></div>
</br>
</br>
-
<span id='little'>Using these reaction, we construct mathematical model of our system that is shown below. Pg(gene probability), in our model, represents the number of plasmid which is inverted. And the rate of producing inverted gene is reduced as the remaining non-inverted gene is reduced.</span></br></br>
+
<span id='little'>Using these reaction, we constructed mathematical model of our system as shown below. Pg(gene probability), in our model, represents the number of plasmid which is inverted. And the rate of producing inverted gene is reduced as the remaining non-inverted gene is reduced.</span></br></br>
<div align="center" ><img src='https://static.igem.org/mediawiki/2012/b/b0/KAIST_Gfprfp2.png'/></div></br>
<div align="center" ><img src='https://static.igem.org/mediawiki/2012/b/b0/KAIST_Gfprfp2.png'/></div></br>
 +
 +
<span id='little'>
 +
Using MATLAB we can solve these set of differential equations and get solution curve like below.
 +
</span></br></br>
 +
<!--graph-->
 +
<div align="center" ><img src='https://static.igem.org/mediawiki/2012/1/13/Result_graph_gfprfp1.png'/></div>
 +
</div>
 +
</br>
 +
<span id='little'>We also consider concentration of GFP/RFP of cell colony. The solution curve shown below represents that result.</span></br>
 +
<div align="center" ><img src='https://static.igem.org/mediawiki/2012/b/b8/KAIST_Result_graph_gfprfp2.png'/></div>
 +
</section>
 +
</br>
 +
<div align='right'><a href="#top">▲ Back to the top</a></div></br></br>
 +
<div>
 +
<ul>
 +
<li style="list-style-type:square;font-size:14px;font-weight:bold;">Parameter Sensitivity Analysis</li>
 +
</ul>
 +
</br>
 +
<span id="little">We also did parameter sensitivity analysis to find what kind of parameters critically impact on our system. We defined sensitivity coefficient and calculated as paramters vary with some ratio. </span></br></br>
 +
<div align="center" ><img src='https://static.igem.org/mediawiki/2012/d/d3/KAIST_Sensitivity_coeff.png'/></div></br>
 +
<span id="little">And we plot the result using 3D bar graph. As graph represents, some parameters are critical to change the output of our system and some are not.</span></br></br>
 +
<div align="center" ><img src='https://static.igem.org/mediawiki/2012/1/1d/KAIST_Sensitivityresult.png'/></div></br>
 +
<div align="center" ><img src='https://static.igem.org/mediawiki/2012/c/c8/Sensitivityresult2.png'/></div>
-
</span>
+
 
 +
</div>
 +
</br></br>
 +
<div align='right'><a href="#top">▲ Back to the top</a></div></br></br>
 +
<section id="3rd">
 +
<div>
 +
<span id="sub-sub-title">Auto Regulation</span></br></br>
 +
<ul>
 +
<li style="list-style-type:square;font-size:14px;font-weight:bold;">Mathematical model</li>
 +
</ul>
 +
</br>
 +
<span id="little">Similar to 'Proof of concept, mathematical model', we can write simplified equation of our system, that is shown below.</span></br></br>
 +
<div align="center" ><img src='https://static.igem.org/mediawiki/2012/thumb/1/1b/KAIST_ReactionofbFMO.png/623px-KAIST_ReactionofbFMO.png'/></div></br>
 +
<span id='little'>Based on the reaction above, we constructed our mathematical model like this.</span></br></br>
 +
<div align="center" ><img src='https://static.igem.org/mediawiki/2012/e/ed/KAIST_ODEbFMO.png'/></div></br>
 +
<span id='little'>Using MATLAB we can solve these set of differential equations and get solution curve like below.</span></br></br>
 +
<div align="center" ><img src='https://static.igem.org/mediawiki/2012/8/8e/KAIST_Result_graph_bFMO1.png'/></div></br>
 +
<span id='little'>When we consider the cell colony instead of one cell, we simply multiply cell growth curve to original result. And derived curve is like below.</span></br></br>
 +
<div align="center" ><img src='https://static.igem.org/mediawiki/2012/7/77/KAIST_Result_graph_bFMO2.png'/></div></br>
</div>
</div>
</section>
</section>
 +
</br>
 +
<div align='right'><a href="#top">▲ Back to the top</a></div></br></br>
 +
<div>
 +
<ul><li style="list-style-type:square;font-size:14px;font-weight:bold;">Parameter Sensitivity Analysis</li></ul>
 +
</br>
 +
<span id='little'>We also did parameter sensitivity analysis and the result is shown below.</span>
 +
<div align="center" ><img src='https://static.igem.org/mediawiki/2012/6/60/KAIST_sensitivityresult3.png'/></div></br>
 +
</div>
 +
</br></br>
 +
<div align='right'><a href="#top">▲ Back to the top</a></div></br></br>
 +
</div>
</div>
</div>
</div>

Latest revision as of 17:39, 26 October 2012

KAIST Korea 2012 iGEM

Project : Modeling
Cell Growth Curve Proof of concept Auto Regulation

Modeling

Computational modeling of our project

Cell Growth Curve

Cell growth can be modeled using Logistic differential equation as shown below.


When we solve this equation with appropriate parameters(using MATLAB), we can get solution curve as shown below. This curve matches with our knowledge about cell growth.





Proof of concept

  • Mathematical model

With GFP and RFP, we want to check whether our system really works or not. Since we cannot consider all the complicated biological phenomena, we assume our system simply follows reaction rate theory and mass balance equations. Reactions below are simplified version of our system. We CONSIDERED the production and degradation of mRNAs and Genes, while we IGNORED the polysome phenomena and any gene regulatory system that occurs in real biological system. Also, we CONSIDERED that there are many copies of plasmids in E.coli and each plasmid can react with invertase to invert their gene sequence.
Because every plasmid has equal probability to react with invertase, we assume plasmids follow uniform distribution. Final assumption is that each plasmid is mutually independent, that is, each plasmid cannot affect invertase reaction of other plasmid.


Using these reaction, we constructed mathematical model of our system as shown below. Pg(gene probability), in our model, represents the number of plasmid which is inverted. And the rate of producing inverted gene is reduced as the remaining non-inverted gene is reduced.


Using MATLAB we can solve these set of differential equations and get solution curve like below.


We also consider concentration of GFP/RFP of cell colony. The solution curve shown below represents that result.



  • Parameter Sensitivity Analysis

We also did parameter sensitivity analysis to find what kind of parameters critically impact on our system. We defined sensitivity coefficient and calculated as paramters vary with some ratio.


And we plot the result using 3D bar graph. As graph represents, some parameters are critical to change the output of our system and some are not.






Auto Regulation

  • Mathematical model

Similar to 'Proof of concept, mathematical model', we can write simplified equation of our system, that is shown below.


Based on the reaction above, we constructed our mathematical model like this.


Using MATLAB we can solve these set of differential equations and get solution curve like below.


When we consider the cell colony instead of one cell, we simply multiply cell growth curve to original result. And derived curve is like below.





  • Parameter Sensitivity Analysis

We also did parameter sensitivity analysis and the result is shown below.





Kaist Footer