IGEM:IMPERIAL/2009/M3/Modelling/

From OpenWetWare
Jump to: navigation, search

<html> <style type="text/css"> /*----------------------------------------------------------------------------*/ /* Structure */ /*----------------------------------------------------------------------------*/

  1. II09_mainWrapper{
   width: 800px;
   align: center;
   border: none;
   margin: 0px auto;
   padding: 0px;

}

/*----------------------------------------------------------------------------*/ /* Content */ /*----------------------------------------------------------------------------*/

  1. II09_content_wrapper {
   width: 550px;
   float: left;
   margin-top: 5px;
   margin-bottom: 5px;
   position:relative;

}

  1. II09_content_boxtop {
   width: 550px;
   height: 10px;
   background:url(http://2009.igem.org/wiki/images/c/c8/II09_content_top.png) no-repeat bottom left;
   position: relative;
   float: left;

}

  1. II09_content_boxbottom {
   width: 550px;
   height: 10px;
   background:url(http://2009.igem.org/wiki/images/b/b4/II09_content_bottom.png) no-repeat top left;
   position: relative;
   float: left;

}

  1. II09_content {
   width: 550px;
   padding: 10px;
   background-color:#ffffff;
   border-left: 1px solid #ccc;
   border-right: 1px solid #ccc;
   position: relative;
   float: left;
   font-size: 9pt;
   text-align: justify;

}

/*----------------------------------------------------------------------------*/ /* Tabs */ /*----------------------------------------------------------------------------*/

.II09_Tabs_Blue {

 float:left;
 width:100%;
 background:#CCCCFF url("http://openwetware.org/images/b/bb/SDbg.gif") repeat-x bottom;
 font-size:93%;
 line-height:normal;
 margin: 0;
 padding: 0;

}

.II09_Tabs_Blue ul {

 margin:0;
 padding:10px 10px 0;
 list-style:none;

} .II09_Tabs_Blue li {

 float:left;
 background:url("http://openwetware.org/images/b/b7/SDleft.gif") no-repeat left top;
 margin:0;
 padding:0 0 0 9px;

}

.II09_Tabs_Blue b {

 float:left;
 display:block;
 background:url("http://openwetware.org/images/f/f8/SDright.gif") no-repeat right top;
 padding:5px 15px 4px 6px;
 text-decoration:none;
 font-weight:bold;
 color:#765;

} </style> </html>

<html> <body> <!--<div id="II09_mainWrapper">--> </html> <html>

<head> <style type="text/css">

  1. II09_Header{
 border-top: 1px solid #7b7b7b;
 border-left: 1px solid #7b7b7b;
 border-right: 1px solid #7b7b7b;
 border-bottom: none;
 background: url(http://2009.igem.org/wiki/images/3/39/II09_Banner.jpg);
 background-position: right -75px;
 width: 783px;
 height: 150px;
 
 margin-left:15px;
 align: right;

} </style> </head>

<body>

<div id="II09_Header"> <center>The E.ncapsulator: the disease killer</center> </div>

</body> </html> <html> <head> <style type="text/css">

  1. II09_Breadcrumb{
  backgrounewd: #7b7b7b;
  background: url(http://2008.igem.org/wiki/images/c/c8/Tudelft_menu_header_bg_1.gif);
  border-top: none;
  border-left: 1px #7b7b7b solid;
  border-right: 1px #7b7b7b solid;
  border-bottom: 1px #7b7b7b solid;
  
  filter:alpha(opacity=100);
  -moz-opacity:1.00;
  opacity:1.0;
  text-color: #000000;
  padding-left: 10px;
  margin-left: 15px;
  width: 773px;

} </style> </head>

<body> <div id="II09_Breadcrumb">

<script type="text/javascript">

var address = document.location.href; var addressSplit = address.split("/");

document.write('<a href="http://2009.igem.org/' +addressSplit[3]+ '">Home (' +addressSplit[3]+ ')</a>');

for (i=4;i<8;i++){

  if (addressSplit[i]!=undefined){
     document.write('<span align="left"><font color=#000000> | </font>');
     document.write('<a href="' +addressSplit[i]+ '">' +addressSplit[i]+ '</a>');
  }

}

  document.write('</span>');

</script> <!--<span align="right">iuvbivb</span>--> </div>

</body> </html> <html> <style type="text/css"> /*----------------------------------------------------------------------------*/ /* Menu */ /*----------------------------------------------------------------------------*/

  1. menu {
   width: 200px;
   margin: 5px 15px 0 15px;
   position: relative;
   float: left;

}

  1. menu ul {
   padding: 0;
   margin: 0 0 10px 0;
   list-style: none;

}

  1. menu li {
   margin:0;

}

  1. menu a {
   display: block;
   width: 185px;
   height: 19px;
   margin-top: 1px;
   margin-left:15px;
   background: url(http://openwetware.org/images/4/46/II09_Menu_ItemHover.png) no-repeat right top;
   padding: 2px 4px;
   text-decoration: none;
   color: #fff;

}

  1. menu a:hover {
   display: block;
   width: 185px;
   height: 19px;
   margin-top: 1px;
   margin-left:15px;
   background: url(http://openwetware.org/images/b/bf/II09_Menu_Item.png) no-repeat right top;
   padding: 2px 4px;
   color: #000;

}

  1. menu .selflink {
   /*font-weight:bold;
   width: 200px;
   border: 1px solid #bbb;
   background: url(http://2008.igem.org/wiki/images/c/c8/Tudelft_menu_header_bg_1.gif) repeat-x bottom;
   color: #111;
   padding: 2px 4px;*/
   display: block;
   width: 200px;
   height: 30px;
   margin-top: 1px;
   border: 1px solid #bbb;
   background: url(http://2008.igem.org/wiki/images/c/c8/Tudelft_menu_header_bg_1.gif) repeat-x left bottom;
   padding: 2px 4px;
   color: #000;

}

  1. menu a.external {
   margin: 0;
   padding: 0;
   width: 156px;

}

.menu_header {

   font-weight:bold;
   width: 200px;
   border: 1px solid #bbb;
   background: url(http://2008.igem.org/wiki/images/c/c8/Tudelft_menu_header_bg_1.gif) repeat-x bottom;
   color: #111;
   padding: 2px 4px;
   text-align:left;

} </style>

<div id="menu"> <ul>

<li><strong class="selflink"><img src="http://2009.igem.org/wiki/images/8/83/II09_Thumb_home.png" width="15%" align="left">Home</strong></li> <li><a href="http://openwetware.org/wiki/IGEM:IMPERIAL/2009/Project"> Project Description</a></li> <li><a href="http://openwetware.org/wiki/IGEM:IMPERIAL/2009/Major_results"> Major results</a></li> <li><a href="http://openwetware.org/wiki/IGEM:IMPERIAL/2009/Biobricks"> Submitted parts</a></li> <li><a href="http://openwetware.org/wiki/IGEM:IMPERIAL/2009/Achievements"> Achievements</a></li> </ul>

<ul> <li class="menu_header"><img src="http://2009.igem.org/wiki/images/c/c1/II09_Thumb_m1.png" width="15%" align="left"> Module 1: Protein production</li> <li><a href="http://openwetware.org/wiki/IGEM:IMPERIAL/2009/M1"> Overview</a></li> <li><a href="http://openwetware.org/wiki/IGEM:IMPERIAL/2009/M1/Genetic"> Genetic circuit</a></li> <li><a href="http://openwetware.org/wiki/IGEM:IMPERIAL/2009/M1/Results"> Results</a></li> <li><a href="http://openwetware.org/wiki/IGEM:IMPERIAL/2009/M1/Modelling"> Modelling</a></li> </ul>

<ul> <li class="menu_header"><img src="http://2009.igem.org/wiki/images/d/df/II09_Thumb_m2.png" width="15%" align="left">Module 2: Encapsulation &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</li> <li><a href="http://openwetware.org/wiki/IGEM:IMPERIAL/2009/M2"> Overview</a></li> <li><a href="http://openwetware.org/wiki/IGEM:IMPERIAL/2009/M2/Genetic"> Genetic circuit</a></li> <li><a href="http://openwetware.org/wiki/IGEM:IMPERIAL/2009/M2/Results"> Results</a></li> <li><a href="http://openwetware.org/wiki/IGEM:IMPERIAL/2009/M2/Modelling"> Modelling</a></li> <li><a href="http://openwetware.org/wiki/IGEM:IMPERIAL/2009/M2/Pictures"> Picture gallery</a></li> </ul>

<ul> <li class="menu_header"><img src="http://2009.igem.org/wiki/images/2/25/II09_Thumb_m3.png" width="15%" align="left">Module 3: Genomic neutralisation</li> <li><a href="http://openwetware.org/wiki/IGEM:IMPERIAL/2009/M3"> Overview</a></li> <li><a href="http://openwetware.org/wiki/IGEM:IMPERIAL/2009/M3/Genetic"> Genetic circuit</a></li> <li><a href="http://openwetware.org/wiki/IGEM:IMPERIAL/2009/M3/Results"> Results</a></li> <li><a href="http://openwetware.org/wiki/IGEM:IMPERIAL/2009/M3/Modelling"> Modelling</a></li> </ul>

<ul> <li class="menu_header"> Resources</li> <li><a href="http://openwetware.org/wiki/IGEM:IMPERIAL/2009/Wetlab"> Wetlab</a></li> <li><a href="http://openwetware.org/wiki/IGEM:IMPERIAL/2009/Drylab"> Drylab</a></li> </ul>

<ul> <li class="menu_header"> Ethics project</li> <li><a href="http://openwetware.org/wiki/IGEM:IMPERIAL/2009/Ethics"> Overview</a></li> <li><a href="http://openwetware.org/wiki/IGEM:IMPERIAL/2009/Ethics/Videos"> Video gallery</a></li> <li><a href="http://openwetware.org/wiki/IGEM:IMPERIAL/2009/Ethics/Contribute"> Contribute!</a></li> </ul>

<ul> <li class="menu_header"> Team</li> <li><a href="http://openwetware.org/wiki/IGEM:IMPERIAL/2009/Contact"> Contact</a></li> <li><a href="http://openwetware.org/wiki/IGEM:IMPERIAL/2009/Team"> Who are we?</a></li> <li><a href="http://openwetware.org/wiki/IGEM:IMPERIAL/2009/Friends"> Our friends</a></li> </ul>

<ul> <li class="menu_header"><a href="http://openwetware.org/wiki/Template:Imperial/09/Menu/OLDLinks">More links (from old menu)</a></li> </ul>

</div> </html> <html>


<div id="II09_content_wrapper"> <img src="http://2008.igem.org/wiki/images/a/a1/Tudelft_content_top.png" width="572"> <div id="II09_content">


<script type="text/javascript"> var gaJsHost = (("https:" == document.location.protocol) ? "https://ssl." : "http://www."); document.write(unescape("%3Cscript src='" + gaJsHost + "google-analytics.com/ga.js' type='text/javascript'%3E%3C/script%3E")); </script> <script type="text/javascript"> try { var pageTracker = _gat._getTracker("UA-11206997-1"); pageTracker._trackPageview(); } catch(err) {}</script>

</html>



Background

Module 3 is on cell killing. An increase in temperature, from 28°C to 42°C will kick off M3. Restriction enzymes DpnII and TaqI will be produced, and these will result in cell death.

We would like to model two things:

1) Thermoinduction with promoter characterisation

The cell death mechanism is triggered by an increase in temperature. An increase in temperature activates the production of restriction enzymes. In our model, we will look at how temperature correlates to restriction enzyme concentration, and seeing how it affects the population of live cells. This model will be linked to data obtained from the staining assay for tracking live and dead cells over time.

From the model, we hope to determine the temperature for optimum cell death.

In promoter characterisation, we would be modelling how the PoPS output of the promoters varies with time. This will help us characterise the promoters for future users' benefit.


2) Killing strategy

In our assays, we have determined quantitatively the enzymatic activity of our restriction enzymes and methylases. The balance between the amount of methylases needed to protect against the restriction enzymes, and the amount of methylase that will stifle restriction enzyme activity, can be modelled.

From our models, we will be able to determine the optimum amount of methylation for killing for different values of promoter leakiness.

Furthermore, from the graph of the percentage of live cells over time time, we can determine the rate of cell death. Hence, we would be able to model the efficiency of the restriction enzymes. We can then predict cell death over a certain a time period.


  • MB's FB:
    • OK but your model does nt include these pieces of info....
    • Drawings + more explanations please

Preliminary Model 1: Thermoinduction: How increase in temperature affects cell growth

Hypothesis

Assumptions

  • The production of restriction enzymes DpnII and TaqI is temperature dependent.
    • Therefore, cells will have a greater probability of dying at higher temperatures (37°C to 42°C).
  • Production is repressed at low temperatures.
  • When temperatures increase beyond threshold, repression on the production of restriction enzymes stops.
  • Hence, restriction enzymes will be produced at higher temperatures.
  • There is a methylase, Dam, to prevent cell death for basal levels of production of restriction enzymes. Dam is produced constitutively.
  • However, as the concentration of restriction enzymes increase at higher temperatures, Dam will no longer be sufficient to protect against cell death.
  • Restriction enzymes DpnII and TaqI will cut the cell DNA. Therefore, the increase in restriction enzyme concentration will increase the rate of cell death.
  • The death term is modelled by an activating hill function of temperature.
    • The hill function is chosen because when temperatures are low, there is no cell death
    • When temperatures increase above threshold, there is suddenly a large increase in cell death, akin to a switch

Predictions

  • If the temperature is below threshold, little restriction enzymes are being produced.
    • Population of live cells will increase exponentially
    • Population of dead cells will also increase exponentially, but should be magnitudes less than live cells.
    • This means that total dead cell population is negligible at low temperatures.
  • If the temperature is above threshold, the population growth of live cells is constrained by the concentration of restriction enzyme.
    • Initially, cells will grow exponentially as the restriction enzymes take time to be produced.
    • When enough restriction enzymes are being produced, cell death will begin.
    • This results in a "bump" in the cell number vs time graph, due to an initial increase in cell number and a subsequent killing of cells.
    • When we increase the temperature, we produce more restriction enzymes at a faster rate. Hence, the higher the temperature, the lower the "bump".


[N.B. We have used arbitrary values in our simulations]



  • MB'sFB:
    • make it clear that the results of promoter characterisation will inform this section and that in the meantime you are using the basic assumption of a Hill function for the effect of temperature
    • Make distinction between assumptions and predictions!!!!
    • Explain what your assumptions are regarding the killing - there are many...

Equations

same remarks as module 1 please explain to your reader what these equations mean

  • they're explained at the bottom. shld we move that up?
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://api.formulasearchengine.com/v1/":): {\displaystyle \begin{alignat}{1} \frac{d[M]}{dt} & = \frac{k_1*{T}^{nt}}{{T}^{nt} + {K_{mt}}^{nt}} - d_m * [M] \qquad (1)\\ \frac{d[E]}{dt} & = k_2 * [M] - d_e * [E] \qquad (2)\\ \frac{d[N]}{dt} & = g * N - d * \frac{k_2*{[E]}^{ne}}{{[E]}^{ne} + {K_{me}}^{ne}} * N \qquad (3)\\ \frac{d[D]}{dt} & = d * \frac{k_2*{[E]}^{ne}}{{[E]}^{ne} + {K_{me}}^{ne}} * N \qquad (4)\\ \end{alignat} }


Legend

Term Variable/Parameter Meaning/description
M Variable mRNA of restriction enzymes DpnII or TaqI
E Variable Restriction enzymes DpnII or TaqI
D Variable Number of dead cells
g Parameter growth constant
d Parameter death constant
k1 Parameter Rate of production of mRNA
k2 Parameter Rate of production of enzyme
T Variable Temperature
nt Variable Hill coefficient for thermal induction
de Parameter Rate of decay of DpnII or TaqI
dm Parameter Rate of decay of mRNA
ne Parameter Hill coefficient for cell death
Kme Parameter Threshold of hill function for cell death
kmt Parameter Threshold of hill function for thermoinduction


Explanation of equations

  • In (1), we have an activating hill function of temperature. This is because an increase in temperature (T) will trigger the production of mRNA and restriction enzymes (DpnII and TaqI) needed for killing. A hill function is chosen because initially, there is very little killing with the increase of temperatures. When the threshold temperature is reached, we start killing cells.
    • At steady state:Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://api.formulasearchengine.com/v1/":): {\displaystyle [M]=\frac{k_{1}*\frac{T^n}{{[K_m]}^n+ T^n}}{d_{m}}} , so when Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://api.formulasearchengine.com/v1/":): {\displaystyle T < K_m} , [M]='off'.
    • When Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://api.formulasearchengine.com/v1/":): {\displaystyle T > K_mt} we produce more [M], which triggers production of restriction enzymes.
  • In (2), the rate of production of enzyme is a balance of rate of transcription (PoPS) and rate of degradation of the enzyme itself.
    • At steady state: Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://api.formulasearchengine.com/v1/":): {\displaystyle [E]=\frac{k_{m}*k_{p}}{d_{m}*d_{p}}*[M]} , so once we produce MRNA, as triggered by the temperature threshold, we will also get restriction enzyme.
  • In (3), the growth term shows that exponential cell growth occurs when there is no increase in temperature. The death term is proportional to the population(N) and is a function of enzyme concentration. When we hit the threshold of enzyme necessary for killing, there will be a increase in the probability of cell death , hence an activating hill function is usedv(effectively, the enzyme hill function acts on increasing the death term so cells are killed over time past the threshold temperature).
    • Consider the death term :Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://api.formulasearchengine.com/v1/":): {\displaystyle d * \frac{k_2*{[E]}^{ne}}{{[E]}^{ne} + {K_{me}}^{ne}} * N } .
    • When restriction enzyme concentrations are subthreshold, the death term is close to 0, so the population grows exponentially (as we see for T=1-4 in simulations)
    • When restriction enzyme concentrations are suprathreshold, this has the effect of increasing the influence of the death term, so now bacterial growth is bounded and saturates (see simulations for T=5-8)
  • In (4)we are describing the time trajectory of dead cells, which is simply the death term in equation (3), but now it is positive.


Simulations

Simulation 1: Effects of increased temperatures on cell population

Very ugly split your results over 2 rows (mRNS and Protein on top) if need be Explain what is special about T=5 and in general how you chose your parameters Explain both cases separately

Temperature below threshold
M3.00.jpg
Temperature above threshold
M3.01.jpg

Explanations

  • The killing term is set to be greater than the growth term in our parameters.

why? explain what this means etc....


  • In the first simulation, we set the arbitrary temperature to be low, varying from 0 to 4. From the simulation, we can tell that when the temperature is low, cell growth will be exponential, ie, there will be no cell death.

explain what T=4 means wrt your activation threshold Give qualitative explanation for system behaviour generalize your results for a temperature dependence that is not Hill-like Comment dead-popultion plot and live-population plot


  • In the second simulation, when we increase the temperature, ie, varying the arbitrary temperature from 5 to 9, cell start to die after an initial exponential cell death. Hence temperature indeed will induce cell death. Our hypothesis is proven right.
  • Also note that at the "bump", this is when cell death remain constant.

what does it mean? is is informative? etc...


N.B. Note that the values we have used for temperature are arbitrary, and they are for the purpose of illustrating the behaviour of the system. In reality, this will be linked to our data, and the parameters will change.

Data Fitting

Matlab files

<html> </div> <img src="http://2008.igem.org/wiki/images/8/89/Tudelft_content_bottom.png" width="572"> </div> <!--</div>--> </body> </html>







































<html> <style type="text/css">

  1. II09_Footer{
 width:763px;
 margin-left: 15px;
 background: url(http://2008.igem.org/wiki/images/c/c8/Tudelft_menu_header_bg_1.gif) repeat-x;
 background-position: bottom;
 padding: 10px;
 border-top: 1px #7b7b7b solid;
 border-left: 1px #7b7b7b solid;
 border-right: 1px #7b7b7b solid;
 border-bottom: 1px #7b7b7b solid;
 position:relative;
 float: left;

} </style>

<body> <div id="II09_Footer"> <a href="http://www.mrgene.com"><img src="http://mrgene.com/portaldata/1/resources/_system/img/logo.gif" alt="Mr. Gene" width="5%"></a> &nbsp; <a href="http://www.geneart.com"><img src="http://www.geneart.com/fileadmin/templates/img/logo.gif" alt="Geneart" width="10%"></a> </div> </body> </html>