Team:KAIST Korea/Project Modeling

From 2012.igem.org

(Difference between revisions)
 
Line 171: Line 171:
</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 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>
<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 have equal probability to react with invertase, we assume plasmids follow uniform distribution. Final assumption is that each plasmids are mutually independent, that is, each  plasmid cannot affect invertase reaction of other plasmid.</br></br></span>
+
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>
Line 194: Line 194:
</ul>
</ul>
</br>
</br>
-
<span id="little">We also do parameter sensitivity analysis to find what kind of parameters are critically impact on our system. We define sensitivity coefficient and calculate as paramters vary with some ratio. </span></br></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>
<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>
+
<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/1/1d/KAIST_Sensitivityresult.png'/></div></br>
<div align="center" ><img src='https://static.igem.org/mediawiki/2012/c/c8/Sensitivityresult2.png'/></div>
<div align="center" ><img src='https://static.igem.org/mediawiki/2012/c/c8/Sensitivityresult2.png'/></div>
Line 213: Line 213:
<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>
<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>
<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 above reaction, we construct our mathematical model like this.</span></br></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>
<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>
<span id='little'>Using MATLAB we can solve these set of differential equations and get solution curve like below.</span></br></br>
Line 227: Line 227:
<ul><li style="list-style-type:square;font-size:14px;font-weight:bold;">Parameter Sensitivity Analysis</li></ul>
<ul><li style="list-style-type:square;font-size:14px;font-weight:bold;">Parameter Sensitivity Analysis</li></ul>
</br>
</br>
-
<span id='little'>We also do parameter sensitivity analysis and result is like below.</span>
+
<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 align="center" ><img src='https://static.igem.org/mediawiki/2012/6/60/KAIST_sensitivityresult3.png'/></div></br>
</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