Difference between revisions of "Team:Exeter/Model"

Line 6: Line 6:
 
   <script src="https://ajax.googleapis.com/ajax/libs/jquery/1.12.4/jquery.min.js"></script>
 
   <script src="https://ajax.googleapis.com/ajax/libs/jquery/1.12.4/jquery.min.js"></script>
 
   <script src="http://maxcdn.bootstrapcdn.com/bootstrap/3.3.6/js/bootstrap.min.js"></script>
 
   <script src="http://maxcdn.bootstrapcdn.com/bootstrap/3.3.6/js/bootstrap.min.js"></script>
 +
<script type="text/javascript" async
 +
  src="https://cdn.mathjax.org/mathjax/latest/MathJax.js?config=TeX-MML-AM_CHTML">
 +
  MathJax.Hub.Config({
 +
  tex2jax: {
 +
    inlineMath: [['$','$'], ['\\(','\\)']],
 +
    processEscapes: true
 +
  }
 +
});
 +
</script>
 
    
 
    
 
   
 
   
Line 585: Line 594:
 
}
 
}
 
if (distHeight > boxHeight){
 
if (distHeight > boxHeight){
$("#topnav").css("bottom",manualPaddingbot);
+
$("#topnav").
}
+
})
+
}
+
//Causes the navbar logo to switch between the default and
+
//inverse versions when the <a> with #logo_button is hovered over
+
$('#logo_button').hover(function(){
+
$('#logo').toggle();
+
$('#logo_inverse').toggle();
+
});
+
});
+
</script>
+
<!--Counter and topnav ends-->
+
<!--Content begins-->
+
<div class="container">
+
<div class="div_vl backgroundimage">
+
<h1 id="title">Modelling</h1>
+
<!--Contains links to sections on page-->
+
<div class="div_banner">
+
<!--Right picture (the teal line on right)-->
+
<div class="col-sm-2 subdiv_banner left">
+
<img src="https://static.igem.org/mediawiki/2016/8/8a/T--Exeter--Template_Banner_left.png">
+
</div>
+
<!--Contains all information of banner in the middle-->
+
<!--of the two outer images-->
+
<div class="col-xs-12 col-sm-8 subdiv_banner middle">
+
<!--Use bootstrap to add links using all 12 columns-->
+
<!--For each link please give double the coloumns for when-->
+
<!--the screen is xs compared to when it is sm-->
+
+
<!--Give span class "oneline" or "twoline" depending on how llong the section text is-->
+
<a href="#section_1" class="banner_link col-xs-6 col-sm-4"><span class="oneline">Introduction</span></a>
+
<a href="#section_2" class="banner_link col-xs-6 col-sm-4"><span class="twoline">KillerRed/<br />KillerOrange</span></a>
+
<a href="#section_3" class="banner_link col-xs-6 col-sm-4"><span class="twoline">Enzymatic <br /> Killswitches</span></a>
+
+
</div>
+
<!--Left picture (the teal line on left)-->
+
<div class="col-sm-2 subdiv_banner right">
+
<img src="https://static.igem.org/mediawiki/2016/5/56/T--Exeter--Template_Banner_right.png">
+
</div>
+
</div>
+
<div>
+
<a id="Section_link" href="#section_1" style="display:block;margin:34vh auto 0 auto;width:14px;"><span style="color:#47BCC2;font-size: 25px;" class="glyphicon glyphicon-menu-down" aria-hidden="true"></span></a>
+
</div>
+
</div>
+
<div class="col-xs-12 div_content">
+
<div id="section_1" class="link_fix"></div>
+
<div id="contentTitle">
+
Introduction
+
</div>
+
 
+
                <p id="pp">As part of our project we have developed two intracellular models to describe our systems. Both codes that we developed used the same protein production mechanism which takes transcription, translation, degradation, maturation, cell division and multiple ribosome binding sites into account. These features are built on top of code that keeps track of the amount of E. Coli and its respective mRNA. This code was made by Joel Burton-Lowe and Andy Wild. It gave the same results for a specific set of parameters and assumptions, and could be independently checked and verified. </p>
+
 
+
                <p id="pp">The KillerRed/KillerOrange model attempts to look at the phototoxicity caused by reactive oxygen species which are induced by light at specific wavelengths, together with how long it takes for cell death. </p>
+
 
+
                <p id="pp">The enzymatic model looks at the rate of degradation caused by ‘Lysozyme C’, and can be expanded to any other enzyme with a substrate vital to the survival of the cell, providing a certain set of parameters are known. </p>
+
 
+
<div>
+
<a id="Section_link" href="#section_2" style="display:block;margin:20px auto 0 auto;width:14px;"><span style="color:#47BCC2;font-size: 25px;" class="glyphicon glyphicon-menu-down" aria-hidden="true"></span></a>
+
</div>
+
</div>
+
<div class="col-xs-12 div_content">
+
<div id="section_2" class="link_fix"></div>
+
<div id="contentTitle">
+
KillerRed/KillerOrange
+
</div>
+
 
+
                <h6>Protein production:</h6>
+
 
+
                <p id="pp">At the outset, multiple assumptions were made to simplify the system, enabling suitable definition of  appropriate variables and rates. Mutations were also ignored to eradicate this further layer of complexity. Mindful that the timeframe considered was much larger than the replication time of <i>E. coli</i>,  the amount of mRNA, Protein etc was adjusted accordingly by assuming that it splits evenly and halving the amount every time replication occurs, also affecting the amount of plasmids. However,  as cited in (Nordström and Dasgupta, 2006), the frequency of replication is variable to maintain a constant amount. Therefore, in the case of our plasmid pSB1C3, which had a high copy number ~300, the amount was kept constant even after replication of the <i>E. Coli</i>.</p>
+
 
+
                <p id="pp">We then isolated one cell and listed the major factors in the production of reactive oxygen species (ROS) starting from the transcription of mRNA.</p>
+
                <br>
+
                <br>
+
 
+
                <div class="col-xs-12" style="margin:0;padding:0;">
+
                                <img src="https://static.igem.org/mediawiki/2016/f/f7/T--Exeter--Modelling_Flow1.png" style="max-width:60%;margin:auto;display:block;">
+
                <br>
+
                <br>
+
                </div>
+
<style>
+
.container ul{
+
padding-top:20px;
+
padding-right:40px;
+
padding-left:80px;
+
font-size:150%;
+
}
+
.container ul li {
+
    /* Text color */
+
    color: #333;
+
    list-style-type: none;
+
    list-style-image:none;
+
}
+
 
+
.container ul li:before {
+
    /* Unicode bullet symbol */
+
    content: '\25AA';
+
    /* Bullet color */
+
    color: #47BCC2;
+
    padding-right: 0.5em;
+
}
+
ol{
+
padding-top:20px;
+
padding-right:40px;
+
padding-left:80px;
+
font-size:150%;
+
}
+
ol li {
+
    /* Text color */
+
    color: #333;
+
    list-style-type: none;
+
}
+
ol li {
+
    list-style-type: none;
+
    counter-increment: list;
+
    position: relative;
+
}
+
 
+
ol li:before {
+
    content: counter(list) ".";
+
    position: absolute;
+
    left: -2.5em;
+
    width: 2em;
+
    text-align: right;
+
    color: #47BCC2;
+
}
+
.caption{
+
margin:auto 6vw;
+
display:block;
+
text-align:center;
+
}
+
</style>
+
 
+
                <ul>
+
                    <li>k1 is the rate of transcription.</li>
+
                    <li>k2 is the rate of mRNA degradation</li>
+
                    <li>k3 is the rate of protein production</li>
+
                    <li>k4 is the rate of degradation of the protein</li>
+
                    <li>k6 is the rate of ROS production</li>
+
                </ul>
+
 
+
                <p id="pp">From this, we isolated the Transcription/Translation mechanism (k1 & k3) along with the degradation (k2 & k4) to obtain a protein production time. </p>
+
 
+
                <p id="pp">K1 was calculated to be 22.175s, taking the total number of base pairs of the plasmid (887) and dividing through by the transcription rate of 40 base pairs per second for the T7 promoter (García and Molineux, 1995).</p>
+
 
+
                <p id="pp">K3 was calculated in a similar way to be 28.690s. As the value found for the rate of 8.4 amino acids per seconds (Siwiak and Zielenkiewicz, 2013) is given in amino acids and not base pairs, this result was arrived at by taking the total number of base pairs of the protein coding region (723) and dividing by 3, given that a conversion of 3 base pairs = 1 amino acid is viable.</p>
+
 
+
                <p id="pp">K2 was found to be 3-8 minutes for 80% of mRNA (Bernstein et al., 2002). 5 minutes was used as an approximation of the median time as the actual number fluctuates.</p>
+
 
+
                <p id="pp">K4 was initially estimated to be greater than 10 hours using protparam on expasy.org. 10 hours was used as the degradation time as this exceed the observed time taken for death of the cell.</p>
+
 
+
                <p id="pp">Simbiology, a modelling app for Matlab, was used to model the flow diagram with the above rate parameters. Running the simulation, a rate for mRNA and protein production could be found. Two issues were identified, one, that the degradation times affected the transcription/translation rate directly and two, that as soon as mRNA production started, protein production was initialised. This suggested that when the mRNA was undergoing transcription, translation had already started. This could not be the case, as transcription could not begin with a non integer amount of mRNA.
+
</p>
+
 
+
                <p id="pp">To correct this, the mRNA was altered to be a step function, limiting the amount to whole integers such that production of the protein was restricted until the first mRNA was produced. However, this only served to add a delay to the overall protein production and multiply it by the mRNA amount as seen in Figure 1. This highlights another problem. As shown in in Figure 2 below. It takes ~28 s to produce a protein (K3) after the mRNA has been made, suggesting that it should occur at  ~54s however, the second mRNA is produced in advance of this, increasing the rate at which the proteins are produced, initiating at  ~53s. This correctly describes the overall protein production rate of the system, but is incorrect in finding the time at which they are made, as each mRNA is independent of any other. This means that each mRNA required  individual consideration, based on each having a separate protein production mechanism that contributes to a total protein quantity. </p>
+
<br>
+
<br>
+
 
+
<div class="col-xs-12" style="margin:0;padding:0;">
+
        <img src="https://static.igem.org/mediawiki/2016/1/17/T--Exeter--Modelling_GraphKRone.png" style="max-width:70%;margin:auto;display:block;">
+
<span class="caption">Figure 1: A graph showing the mRNA amount in blue over 100 seconds and the overall rate of Protein production in red over 100 seconds.</span>
+
<br>
+
<br>
+
</div>
+
 
+
<div class="col-xs-12" style="margin:0;padding:0;">
+
        <img src="https://static.igem.org/mediawiki/2016/0/03/T--Exeter--Modelling_GraphKRtwo.png" style="max-width:70%;margin:auto;display:block;">
+
<span class="caption">Figure 2: A plot of the overall protein production in blue and the rate at which the first mRNA independently produces a protein in red. Creating proteins at ~53s and ~54s respectively.</span>
+
 
+
<br>
+
<br>
+
</div>
+
 
+
                <p id="pp">To enable this we moved away from Simbiology and attempted to use Simulink. However, after careful consideration, we decided instead to write the process in C. This allowed us to handle each mRNA and its creation time separately and have it produce proteins up to the time in which the protein was induced. An array was used with each element representing a second; every time a protein was created the corresponding time element would increase by one. The total of all elements were then taken to find an overall amount of protein.
+
</p>
+
 
+
                <p id="pp">The next step was to include degradation. This was achieved by limiting the protein production of an individual mRNA to its respective degradation time. For protein degradation, the amount of time left from the time of creation until the overall run time (in our case the point at which we started shining light on the cultures) was compared to the degradation time; if the difference was greater, the protein degrades and is subtracted from the total.
+
</p>
+
 
+
                <p id="pp">Further, we took into account multiple ribosomes on one mRNA (polyribosome) and maturation. Polyribosomes were included in the code by restricting the first protein created by each mRNA to the the calculated translation rate found above (~28s). Then for the following proteins we used a modified rate where we divided the rate by the amount of ribosomes, which was found to be 3.46 per 100 codons (Siwiak and Zielenkiewicz, 2013). In this case, it was calculated to be 8 if we rounded down to the nearest integer. As the degradation and the ability to produce ROS applies to only folded proteins, we wanted to know how many mature proteins there would be. Calculating of a new run time was achieved by  taking away both the maturation time in addition to the time taken to make one mRNA and protein,  suggesting that the proteins made in the new run time were created and matured and ready to produce ROS if they had not degraded.  </p>
+
 
+
                <p id="pp">Examining the degradation time for the protein, we adjusted the previously estimated value to be that of green fluorescent protein (GFP), as KillerRed is a homologue to provide a more accurate degradation time. However, able to only find a half-life of the protein, we used equation (1) to calculate the degradation.
+
</p>
+
 
+
<math xmlns="http://www.w3.org/1998/Math/MathML" display="block">
+
  <mi>T</mi>
+
  <mo>=</mo>
+
  <mfrac>
+
    <msub>
+
      <mi>T</mi>
+
      <mfrac>
+
      <mn>1</mn>
+
      <mn>2</mn>
+
      </mfrac>
+
    </msub>
+
    <mrow>
+
      <mi>ln</mi>
+
      <mo>&InvisibleTimes;</mo>
+
      <mn>2</mn>
+
    </mrow>
+
  </mfrac>
+
  <mspace width="2em" />
+
  <mtext>(1)</mtext>
+
</math> 
+
 
+
                <p id="pp">Running the initial model provided an overall protein quantity of under 3 million and mRNA at a maximum of 4200 (Thermofisher.com, 2016), which fits within the expected amount. This suggests  that the rates and process are a good, basic approximation of protein production.
+
</p>
+
 
+
                <h6>ROS production:</h6>
+
 
+
                <p id="pp">The next step was to look at the rate of ROS production (k5). Using the equations provided by (Song et al, 1995).</p>
+
<br>
+
<br>
+
 
+
              <math xmlns="http://www.w3.org/1998/Math/MathML" display="block">
+
  <mtable columnalign="right left right left right left right left right left right left" rowspacing="3pt" columnspacing="0em 2em 0em 2em 0em 2em 0em 2em 0em 2em 0em" displaystyle="true">
+
    <mtr>
+
      <mtd>
+
        <mi>S</mi>
+
        <mo>+</mo>
+
        <mi>h</mi>
+
        <mi>&#x03BD;<!-- ? --></mi>
+
        <mo stretchy="false">&#x2192;<!-- ? --></mo>
+
        <msup>
+
          <mi>S</mi>
+
          <mo>&#x2217;<!-- * --></mo>
+
        </msup>
+
        <mspace width="2em" />
+
        <mspace width="2em" />
+
      </mtd>
+
      <mtd>
+
        <msub>
+
          <mi>k</mi>
+
          <mi>a</mi>
+
        </msub>
+
      </mtd>
+
      <mtd>
+
        <mspace width="2em" />
+
        <mtext>Activation</mtext>
+
      </mtd>
+
    </mtr>
+
    <mtr>
+
      <mtd>
+
        <msup>
+
          <mi>S</mi>
+
          <mo>&#x2217;<!-- * --></mo>
+
        </msup>
+
        <mo stretchy="false">&#x2192;<!-- ? --></mo>
+
        <mi>S</mi>
+
        <mo>+</mo>
+
        <mi>h</mi>
+
        <mi>&#x03BD;<!-- ? --></mi>
+
        <mspace width="2em" />
+
        <mspace width="2em" />
+
      </mtd>
+
      <mtd>
+
        <msub>
+
          <mi>k</mi>
+
          <mi>f</mi>
+
        </msub>
+
      </mtd>
+
      <mtd>
+
        <mspace width="2em" />
+
        <mtext>Fluorescence emission</mtext>
+
      </mtd>
+
    </mtr>
+
    <mtr>
+
      <mtd>
+
        <msup>
+
          <mi>S</mi>
+
          <mo>&#x2217;<!-- * --></mo>
+
        </msup>
+
        <mo stretchy="false">&#x2192;<!-- ? --></mo>
+
        <mi>S</mi>
+
        <mspace width="2em" />
+
        <mspace width="2em" />
+
      </mtd>
+
      <mtd>
+
        <msub>
+
          <mi>k</mi>
+
          <mrow class="MJX-TeXAtom-ORD">
+
            <mi>i</mi>
+
            <mi>c</mi>
+
          </mrow>
+
        </msub>
+
        <mo>=</mo>
+
        <mn>6.9</mn>
+
        <mo>&#x00D7;<!-- × --></mo>
+
        <msup>
+
          <mn>10</mn>
+
          <mrow class="MJX-TeXAtom-ORD">
+
            <mn>10</mn>
+
          </mrow>
+
        </msup>
+
        <msup>
+
          <mi>s</mi>
+
          <mrow class="MJX-TeXAtom-ORD">
+
            <mo>&#x2212;<!-- - --></mo>
+
            <mn>1</mn>
+
          </mrow>
+
        </msup>
+
        <mspace width="2em" />
+
      </mtd>
+
      <mtd>
+
        <mtext>Internal conversion</mtext>
+
      </mtd>
+
    </mtr>
+
    <mtr>
+
      <mtd>
+
        <msup>
+
          <mi>S</mi>
+
          <mo>&#x2217;<!-- * --></mo>
+
        </msup>
+
        <mo stretchy="false">&#x2192;<!-- ? --></mo>
+
        <msup>
+
          <mi>T</mi>
+
          <mo>&#x2217;<!-- * --></mo>
+
        </msup>
+
        <mspace width="2em" />
+
        <mspace width="2em" />
+
      </mtd>
+
      <mtd>
+
        <msub>
+
          <mi>k</mi>
+
          <mrow class="MJX-TeXAtom-ORD">
+
            <mi>i</mi>
+
            <mi>s</mi>
+
            <mi>c</mi>
+
          </mrow>
+
        </msub>
+
        <mo>=</mo>
+
        <mn>6.6</mn>
+
        <mo>&#x00D7;<!-- × --></mo>
+
        <msup>
+
          <mn>10</mn>
+
          <mn>6</mn>
+
        </msup>
+
        <msup>
+
          <mi>s</mi>
+
          <mrow class="MJX-TeXAtom-ORD">
+
            <mo>&#x2212;<!-- - --></mo>
+
            <mn>1</mn>
+
          </mrow>
+
        </msup>
+
        <mspace width="2em" />
+
      </mtd>
+
      <mtd>
+
        <mtext>Intersystem crossing</mtext>
+
      </mtd>
+
    </mtr>
+
    <mtr>
+
      <mtd>
+
        <msup>
+
          <mi>T</mi>
+
          <mo>&#x2217;<!-- * --></mo>
+
        </msup>
+
        <mo stretchy="false">&#x2192;<!-- ? --></mo>
+
        <msup>
+
          <mi>S</mi>
+
        </msup>
+
        <mspace width="2em" />
+
        <mspace width="2em" />
+
      </mtd>
+
      <mtd>
+
        <msub>
+
          <mi>k</mi>
+
          <mrow class="MJX-TeXAtom-ORD">
+
            <mi>r</mi>
+
            <mi>d</mi>
+
          </mrow>
+
        </msub>
+
        <mo>=</mo>
+
        <mn>50</mn>
+
        <msup>
+
          <mi>s</mi>
+
          <mrow class="MJX-TeXAtom-ORD">
+
            <mo>&#x2212;<!-- - --></mo>
+
            <mn>1</mn>
+
          </mrow>
+
        </msup>
+
        </mrow>
+
        <mspace width="2em" />
+
      </mtd>
+
      <mtd>
+
        <mtext>Radiationless deactivation</mtext>
+
      </mtd>
+
    </mtr>
+
    <mtr>
+
      <mtd>
+
        <msup>
+
          <mi>T</mi>
+
          <mo>&#x2217;<!-- * --></mo>
+
        </msup>
+
        <mo>+</mo>
+
        <msub>
+
          <mi>O</mi>
+
          <mn>2</mn>
+
        </msub>
+
        <mo stretchy="false">&#x2192;<!-- ? --></mo>
+
        <mi>S</mi>
+
        <mo>+</mo>
+
        <msubsup>
+
          <mi>O</mi>
+
          <mn>2</mn>
+
          <mrow class="MJX-TeXAtom-ORD">
+
            <mo>&#x2219;<!-- · --></mo>
+
            <mo>&#x2212;<!-- - --></mo>
+
          </mrow>
+
        </msubsup>
+
        <mspace width="2em" />
+
        <mspace width="2em" />
+
      </mtd>
+
      <mtd>
+
        <msub>
+
          <mi>k</mi>
+
          <mrow class="MJX-TeXAtom-ORD">
+
            <mi>p</mi>
+
            <mi>q</mi>
+
          </mrow>
+
        </msub>
+
        <mo>=</mo>
+
        <mn>1.56</mn>
+
        <mo>&#x00D7;<!-- × --></mo>
+
        <msup>
+
          <mn>10</mn>
+
          <mn>9</mn>
+
        </msup>
+
        <msup>
+
          <mi>M</mi>
+
          <mrow class="MJX-TeXAtom-ORD">
+
            <mo>&#x2212;<!-- - --></mo>
+
            <mn>1</mn>
+
          </mrow>
+
        </msup>
+
        <msup>
+
          <mi>s</mi>
+
          <mrow class="MJX-TeXAtom-ORD">
+
            <mo>&#x2212;<!-- - --></mo>
+
            <mn>1</mn>
+
          </mrow>
+
        </msup>
+
        <mspace width="2em" />
+
      </mtd>
+
      <mtd>
+
        <mtext>Physical quenching</mtext>
+
      </mtd>
+
    </mtr>
+
    <mtr>
+
      <mtd>
+
        <msup>
+
          <mi>T</mi>
+
          <mo>&#x2217;<!-- * --></mo>
+
        </msup>
+
        <mo>+</mo>
+
        <mn>2</mn>
+
        <msub>
+
          <mi>O</mi>
+
          <mn>2</mn>
+
        </msub>
+
        <mo stretchy="false">&#x2192;<!-- ? --></mo>
+
        <mi>X</mi>
+
        <mo>+</mo>
+
        <msubsup>
+
          <mi>O</mi>
+
          <mn>2</mn>
+
          <mrow class="MJX-TeXAtom-ORD">
+
            <mo>&#x2219;<!-- · --></mo>
+
            <mo>&#x2212;<!-- - --></mo>
+
          </mrow>
+
        </msubsup>
+
        <mspace width="2em" />
+
        <mspace width="2em" />
+
      </mtd>
+
      <mtd>
+
        <msub>
+
          <mi>k</mi>
+
          <mrow class="MJX-TeXAtom-ORD">
+
            <mi>c</mi>
+
            <mi>q</mi>
+
          </mrow>
+
        </msub>
+
        <mo>=</mo>
+
        <mn>1.4</mn>
+
        <mo>&#x00D7;<!-- × --></mo>
+
        <msup>
+
          <mn>10</mn>
+
          <mn>8</mn>
+
        </msup>
+
        <msup>
+
          <mi>M</mi>
+
          <mrow class="MJX-TeXAtom-ORD">
+
            <mo>&#x2212;<!-- - --></mo>
+
            <mn>1</mn>
+
          </mrow>
+
        </msup>
+
        <msup>
+
          <mi>s</mi>
+
          <mrow class="MJX-TeXAtom-ORD">
+
            <mo>&#x2212;<!-- - --></mo>
+
            <mn>1</mn>
+
          </mrow>
+
        </msup>
+
        <mspace width="2em" />
+
      </mtd>
+
      <mtd>
+
        <mtext>Chemical quenching</mtext>
+
      </mtd>
+
    </mtr>
+
  </mtable>
+
</math>
+
 
+
<br>
+
<br>
+
 
+
                <p id="pp">Assumptions were once again made to simplify the system. The first of these being that all singlets are assumed to be in the excited states as there is an excess of photons with an energy exceeding that of activation. As well as this, we assumed that one chromophore produces one ROS before ‘resetting’ to produce the next. Lastly, the damage done by the ROs is negated when the cell splits.</p>
+
+
                <p id="pp">Combining all assumptions in Simbiology gave the value of k5 to be 8.24x10$^{-6}\text{s}$</p>
+
 
+
                <p id="pp">From this we can determine the estimated time of cell death by reading the time off at the first occurrence of the internal ROS threshold. However, a well-documented threshold could not be found meaning that a time could not be estimated.</p>
+
 
+
                <h6>Future</h6>
+
 
+
                <p id="pp">Future improvements of this code would include:</p>
+
 
+
                <ul>
+
                    <li>Implementing a photobleaching rate degrading mature proteins, this will depend on the amount of ROS.</li>
+
                    <li>Doing further research on the internal limit of ROS to determine a threshold. This would allow a time of cell death to be estimated.</li>
+
                    <li>Improving the process to reduce the amount of assumptions, or constrain the model to limits where these assumptions are minimised.</li>
+
                </ul>
+
 
+
<div>
+
<a id="Section_link" href="#section_3" style="display:block;margin:20px auto 0 auto;width:14px;"><span style="color:#47BCC2;font-size: 25px;" class="glyphicon glyphicon-menu-down" aria-hidden="true"></span></a>
+
</div>
+
</div>
+
<div class="col-xs-12 div_content">
+
<div id="section_3" class="link_fix"></div>
+
<div id="contentTitle">
+
Enzymatic Models
+
</div>
+
<h5>mRNA production in Enzymatic kill switches</h5>
+
<p id="pp">
+
The change from software such as Simbiology and Simulink called for a more fundamental method of modelling cell death. Preliminary research
+
showed kill switches producing the proteins “Lysozyme c” and “DNase 1” both had very similar mechanisms; as both are enzymes. Therefore, it
+
was decided that the two models would use the same code to simulate mRNA and protein production.
+
</p>
+
<p id="pp">
+
Initially, following advice from biologists and biochemists on our team, the first code incorporated two step functions, the first of these
+
modelled mRNA production by a single <i>E. coli</i> cell. To calculate protein made by each mRNA a secondary step function triggered each time an mRNA
+
was produced, the sum of these functions gave the total amount of protein. This program was simple with the only input variables being the production time of
+
mRNA and the respective protein.
+
The initial model was presented to the rest of the team receiving plenty of feedback, the most prevalent point being that the code modelled a single
+
cell system with a single plasmid whilst it should model a single cell system that duplicates and has multiple plasmids. It was suggested the following factors were added to the model:
+
</p>
+
<ul>
+
<li>Plasmid production</li>
+
<li>Degradation of mRNA and protein</li>
+
<li>Maturation of protein</li>
+
<li>Duplication rate of <i>E. coli</i></li>
+
</ul>
+
<p id="pp">
+
The initial feedback called for a re-write of the code as a considerable amount of the suggestions came in stages before mRNA production.
+
A second version of the code was written which incorporated all main steps between the splitting of <i>E. coli</i> to the degradation rate of protein.
+
</p>
+
<h6>Assumptions</h6>
+
<p id="pp">
+
It is important to outline the factors that were overlooked due to research finding their effect on the model would be negligible. Firstly, after
+
researching the production time of plasmids in <i>E. coli</i> it was found that plasmids will reproduce at a variable rate to maintain a constant population
+
determined by their copy number (Nordström and Dasgupta, 2006). To address this, the production rate of plasmids was overlooked, allowing for a
+
constant value of plasmids to be maintained throughout the simulation. In experiments a pSB1C3 strain was used - a high copy number plasmid;
+
therefore the copy number was set to 300.
+
Secondly, both enzymatic models will assume that travel time of protein to the substrate is negligible. Protein diffusing through the cytoplasm of
+
<i>E. coli</i> have a diffusion coefficient on the scale of $10 \mu \text{m}^2 \text{s}^{-1}$ (Elowitz et al., 1998). Considering the surface area of <i>E. coli</i> is approximately
+
$10 \mu \text{m}^2$, the protein will reach the cell wall in a several seconds which is several magnitudes of order smaller than the simulation time. Lastly, the models
+
will assume no mutations occur; the aim of the simulations is to determine whether kill switches are a plausible method of biosafety, if the models show
+
a kill switch is not a reliable way to terminate GMO’s then accounting for mutations will only that enforce statement.
+
</p>
+
<h6>Features</h6>
+
<p id="pp">
+
Research showed that there is an upper limit of mRNA in <i>E. coli</i>,
+
therefore an upper limit of $4 \times 10^3$ mRNA per <i>E. coli</i> cell (Thermofisher.com, 2016) has been included in the model. The lifetime of mRNA can be found
+
from the observed half life of approximately 5 minutes or $300\text{s}$ (Bernstein et al., 2002), resulting in an average lifetime of $430\text{s}$. The production rate of
+
mRNA along with ribosomes per coding region will be worked out for both the lysozyme and DNase models independently. In addition to this, both
+
enzymatic models use the well known duplication time of <i>E. coli</i> - 17 minutes, at which point both the mRNA and protein is assumed to split among the two cells equally,.
+
</p>
+
<p id="pp">
+
<span class="equation">$T_{(mRNA)} = \frac{T_{(mRNA) \frac{1}{2}}}{ln(2)} = \frac{300\text{s}}{ln(2)} = 430\text{s (2sf)}$<span class="equation_ref">(Hyperphysics.phy-astr.gsu.edu, 2016)</span></span>
+
<span class="equation_key">
+
$T_{(mRNA)}$: Lifetime of mRNA [$\text{s}$]<br />
+
$T_{(mRNA) \frac{1}{2}}$: Half life of mRNA [$\text{s}$]
+
</span>
+
</p>
+
<h5>Lysozyme Model</h5>
+
<p id="pp">
+
In the case of the “Lysozyme c” kill switch the mechanisms beyond producing mRNA need to be modelled separately from the DNase model. There are several
+
assumptions that will be made, the first of these is that enzymatic reactions can be modelled by Michaelis-Menten kinetics. Secondly, the temperature
+
of the constants taken imply that this model is running in the range of $37-40^o\text{C}$ at an optimal pH for <i>E. coli</i> growth. The model will assume that when
+
<i>E. coli</i> splits, the contents of the cell and damage of the cell wall is shared equally among the two resulting <i>E. coli</i>. Lastly, it will assume that lysozyme does not degrade
+
throughout the simulation.
+
</p>
+
<p id="pp">
+
The plasmid used to produce lysozyme has a PCR with a length of $507\text{bp}$ with the promoter, RBS and terminator totalling a further $164\text{bp}$. Therefore the
+
production rates of mRNA and lysozyme can be calculated using translation and transcription rates of $V_{translation} = 8.4\text{aas}^{-1}$
+
(Siwiak and Zielenkiewicz, 2013) and $V_{transcription} = 40\text{bps}^{-1}$ (García and Molineux, 1995) respectively. The translation time is for <i>E. coli</i> at $37^o\text{C}$, other values
+
have been taken at $40^o\text{C}$ as this was the closest temperature that could be found. All reaction rates have been rounded to the nearest second as this
+
reduces calculation times.
+
</p>
+
<p id="pp">
+
<span class="equation">$t_{lysozyme} = \frac{L_{protein}}{V_{translation}} = \frac{\frac{507\text{bp}}{3}}{8.4\text{aas}^{-1}} = 20\text{s (2sf)}$</span><br />
+
<span class="equation">$t_{mRNA} = \frac{L_{plasmid}}{V_{transcription}} = \frac{507\text{bp} + 164\text{bp}}{40\text{bps}^{-1}} = 17\text{s (2sf)}$</span>
+
<span class="equation_key">
+
$t_{lysozyme}$: Time to produce one lysozyme protein [$\text{s}$]<br />
+
$t_{mRNA}$: Time to produce one mRNA [$\text{s}$]<br />
+
$L_{plasmid}\text{, }L_{protein}$: Length of plasmid and protein coding region [$\text{bp}$]
+
</span>
+
</p>
+
<p id="pp">
+
To supplement the production of lysozyme, the program will implement the affects of multiple ribosomes on the coding site of lysozyme. For <i>E. coli</i>
+
it has been found there are 3.46 codons per $100\text{bp}$ (Siwiak and Zielenkiewicz, 2013). The PCR used for lysozyme production has a length of $507\text{bp}$, meaning
+
there are approximately 5.8 codons which will be rounded down to 5 as not to overproduce lysozyme.
+
</p>
+
<p id="pp">
+
“Lysozyme c” is an enzyme that hydrolyses bonds holding together peptidoglycan in the cell wall, this is explained in detail on the <a href="https://2016.igem.org/Team:Exeter/Project">project page</a>. The
+
next task of the lysozyme model is to connect the amount of protein at each time to the degradation of the <i>E. coli</i> cell wall, to do this
+
Michaelis-Menten kinetics were applied, which gives the reaction rate of one enzyme (Berg et al., 2002).
+
</p>
+
<p id="pp">
+
<span class="equation">$k_{(cat)} = \frac{[S]k_{(cat)max}}{[S] + K_M}$</span>
+
<span class="equation_key">
+
$k_{(cat)}$: Reaction rate of one lysozyme [$\text{s}^{-1}$]<br />
+
$k_{(cat)max}$: Maximum reaction rate of one lysozyme [$\text{s}^{-1}$]<br />
+
$[S]$: Substrate concentration [$\text{M}$]<br />
+
$K_M$: Michaelis constant [$\text{M}$]
+
</span>
+
</p>
+
<p id="pp">
+
To use this model two constants are required, the Michaelis constant ($K_M$), the concentration of the substrate when the reaction rate is exactly one half of the maximum reaction rate.
+
The average reaction rate assumed to be $k_{(cat)avg} = (k_{(cat)max}/2$). The logarithms of both of these
+
values has been calculated at $40^oC$ to be $-log(K_M) = 5.18 \pm 0.3$M and $-log(k_{cat}^{obs}) = 0.15 \pm 0.005$s$^{-1}$ (Banerjee et al., 1975). Giving values of:
+
</p>
+
<p id="pp">
+
<span class="equation">$K_M = 5.6\text{mM}$ (2sf)</span><br />
+
<span class="equation">$k_{(cat)avg} = \frac{k_{(cat)max}}{2} \approx \frac{k_{(cat)}^{obs}}{2} = \frac{0.86\text{s}^{-1}}{2} = 0.43\text{s}^{-1}$ (2sf)</span>
+
</p>
+
<p id="pp">
+
Lastly, the initial concentration of the substrate peptidoglycan is calculated. An assumption is made that determines that all the peptidoglycan in
+
<i>E. coli</i> is spread out over the entire volume of the cell. Peptidoglycan or murein amount has been calculated to be approximately $3.5\text{x}10^6$ molecules
+
per cell in a strain of <i>E. coli</i> (Vollmer and Höltje, 2004). Using an approximate volume of <i>E. coli</i> of $0.7 \mu \text{m}^3$, the concentration of peptidoglycan is:
+
</p>
+
<p id="pp">
+
<span class="equation">$[Pep]_{int} = \frac{\frac{N_{pep}}{N_A}}{V_{E. coli}} = \frac{\frac{3.5\text{x}10^6}{6.02\text{x}10^{23}}}{0.7 \mu \text{m}^3} = 8.3\text{mM}$ (2sf)</span>
+
<span class="equation_key">
+
$[Pep]_{int}$: Initial concentration of peptidoglycan [$\text{mM}$]<br />
+
$N_{pep}$: Amount of peptidoglycan [molecules]<br />
+
$N_A$: Avogadro's constant [molecules/mole]<br />
+
$V_{\textit{E. coli}}$: Volume of <i>E. coli</i> [$\mu\text{m}^3$]
+
</span>
+
</p>
+
<p id="pp">
+
This calculation gives a value in the same order of magnitude as the Michaelis constant, which represents the concentration of substrate when the
+
reaction rate is at half of its maximum value.
+
</p>
+
<h5>Results</h5>
+
<div class="col-xs-12" style="width:100%;position:relative;margin:auto;padding:0;">
+
<div class="graph_box_single col-xs-12">
+
<img src="https://static.igem.org/mediawiki/2016/2/24/T--Exeter--Modelling_Enzyme_Graph3.png">
+
<span>Fig. 1. Using Michaelis-Menten kinetics the reaction rate of each lysozyme enzyme has
+
been plotted for each peptidoglycan substrate concentration. The
+
average reaction rate of $0.43\text{s}^{-1}$ occurs when the concentration is equal to $K_M = 0.0056\text{M}$.
+
The maximum or initial concentration $[Pep]_{int} = 0.0083\text{M}$ of the substrate causes a reaction rate of $0.51\text{s}^{-1}$.</span>
+
</div>
+
</div>
+
<p id="pp">
+
Enzyme reaction rates have been modelled by the Michaelis-Menten kinetics model in Fig. 1, therefore the reaction
+
rate decreases as the substrate concentration decreases. This graph shows that the reaction
+
rate will be greatest at the beginning of the simulation and approach zero when the cell wall is most damaged.
+
</p>
+
<div class="col-xs-12" style="width:100%;position:relative;margin:auto;padding:0;">
+
<div class="graph_box col-xs-12">
+
<img src="https://static.igem.org/mediawiki/2016/8/84/T--Exeter--Modelling_Enzyme_Graph1.png">
+
<span>Fig. 2. The percentage of peptidoglycan compared to the original concentration plotted against time.</span>
+
</div>
+
<div class="graph_box col-xs-12">
+
<img src="https://static.igem.org/mediawiki/2016/a/a1/T--Exeter--Modelling_Enzyme_Graph2.png">
+
<span>Fig. 3. Plots a smaller range of times as Fig. 2. To show the rapid decrease in peptidoglycan concentration.</span>
+
</div>
+
</div>
+
<p id="pp">
+
The model predicted the complete degradation of the cell wall to be within the first generation
+
of <i>E. coli</i>, Fig. 2. The reaction rate is slow at first due to the cell having no initial lysozyme,
+
this slowly increases, until the low concentration of the substrate casues
+
the reaction rate of lysozyme to slow considerably. The peptidoglycan concentration in the cell is
+
negligible until 17 minutes at which point the <i>E. coli</i> splits sharing the cell wall damage equally
+
between the two daughter cells, hence cell damage drops from almost 100% to 50%. The immediate concern
+
is that in this model cell death would occur far before it is able to duplicate, meaning that
+
assuming no mutations the cell would terminate before 17 minutes.
+
</p>
+
<p id="pp">
+
The cell death threshold of peptidoglycan concentration in the cell wall is not well defined from
+
research. Fig. 3 demonstrates that any threshold that is chosen is likely to fall in between 2
+
and 5 minutes of the simulation which is well before the reproduction rate of <i>E. coli</i> of 17
+
minutes. Therefore it is a reasonable assumption that given no mutations were to occur that the
+
cell would be terminated before the <i>E. coli</i> could reproduce.
+
</p>
+
<h5>References</h5>
+
<ol style="font-size:100%;">
+
<li>Nordström, K. and Dasgupta, S. (2006). Copy-number control of the <i>Escherichia coli</i> chromosome: a plasmidologist's view. EMBO Rep, [online] 7(5), pp.484-489. Available at: http://www.ncbi.nlm.nih.gov/pmc/articles/PMC1479556/ [Accessed 6 Sep. 2016].</li>
+
<li>Elowitz, M., Surette, M., Wolf, P., Stock, J. and Leibler, S. (1998) ‘Protein mobility in the cytoplasm of Escherichia coli’, Journal of bacteriology., 181(1), pp. 197–203 [Accessed 6 Sep. 2016].</li>
+
<li>Thermofisher.com. (2016). Macromolecular Components of E. coli and HeLa Cells | Thermo Fisher Scientific. [online] Available at: https://www.thermofisher.com/uk/en/home/references/ambion-tech-support/rna-tools-and-calculators/macromolecular-components-of-e.html# [Accessed 6 Sep. 2016].</li>
+
<li>Bernstein, J., Khodursky, A., Lin, P., Lin-Chao, S. and Cohen, S. (2002). Global analysis of mRNA decay and abundance in Escherichia coli at single-gene resolution using two-color fluorescent DNA microarrays. Proceedings of the National Academy of Sciences, [online] 99(15), pp.9697-9702. Available at: http://www.pnas.org/content/99/15/9697.long [Accessed 6 Sep. 2016].</li>
+
<li>Berg, J., Tymoczko, J., Stryer, L. and Stryer, L. (2002). Biochemistry. New York: W.H. Freeman, pp.Section 8.4 - Equation (23).</li>
+
<li>Banerjee, S., Holler, E., Hess, G. and Rupley, J. (1975). Reaction of N-acetylglucosamine oligosaccharides with lysozyme. Temperature, pH, and solvent deuterium isotope effects; equilbrium, steady state, and pre-steady state measurements*. Journal of Biological Chemistry, [online] 250(11), pp.4357, Figure 2 and 4359, Table I. Available at: http://www.jbc.org/content/250/11/4355.long [Accessed 7 Sep. 2016].</li>
+
<li>Vollmer, W. and Höltje, J. (2004). The Architecture of the Murein (Peptidoglycan) in Gram-Negative Bacteria: Vertical Scaffold or Horizontal Layer(s)?. Journal of Bacteriology, [online] 186(18), p.5980. Available at: http://jb.asm.org/content/186/18/5978 [Accessed 7 Sep. 2016].</li>
+
<li>Siwiak, M. and Zielenkiewicz, P. (2013) ‘Transimulation - protein Biosynthesis web service’, PLoS ONE, 8(9), p. e73943. doi: 10.1371/journal.pone.0073943. [Accessed 7 Sep. 2016]</li>
+
<li>García, L. and Molineux, I. (1995) ‘Rate of translocation of bacteriophage T7 DNA across the membranes of Escherichia coli’, Journal of bacteriology., 177(14), pp. 4066–76.[Accessed 7 Sep. 2016]</li>
+
<li>Siwiak, M. and Zielenkiewicz, P. (2013). Transimulation - Protein Biosynthesis Web Service. PLoS ONE, [online] 8(9), p.3, left column, second paragraph. Available at: http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3764131/ [Accessed 7 Sep. 2016].</li>
+
<li>Hyperphysics.phy-astr.gsu.edu. (2016). Mean Lifetime for Particle Decay. [online] Available at: http://hyperphysics.phy-astr.gsu.edu/hbase/nuclear/meanlif.html [Accessed 10 Sep. 2016].</li>
+
                    <li>Song L, Hennink EJ, Young IT and Tanke HJ. Photobleaching Kinetics of Fluorescein in Quantitative Fluorescence Microscopy. Biophys. J. 1995 68: 2588-2600. [Accessed 10 Sep. 2016].</li>
+
</ol>
+
+
+
+
</div>
+
+
</div>
+
 
+
<!--Content Ends-->
+
<nav id="topnav" class="nav-navbar" style="right:2px;">
+
  <div class="container-fluid">
+
    <div class="navbar-header">
+
      <a class="navbar-brand" id="Top_link" href="#top_page">
+
<span class="glyphicon glyphicon-menu-up" aria-hidden="true" style="max-width:1%;"></span>
+
  </a>
+
    </div>
+
  </div>
+
</nav>
+
<div class="flare"></div>
+
<footer class="footer">
+
<a style="margin-left:3%" href="https://www.youtube.com/channel/UC31qfG4hnm8gRHDCkrBtAiQ">
+
<img  src="https://static.igem.org/mediawiki/2016/2/23/You.png"/>
+
</a>
+
<a href="https://www.facebook.com/ExeteriGEM2016/?ref=aymt_homepage_panel">
+
<img src="https://static.igem.org/mediawiki/2016/5/51/Fb.png"/>
+
</a>
+
<a href="https://twitter.com/exeter_igem">
+
<img src="https://static.igem.org/mediawiki/2016/7/76/Twit.jpg"/>
+
</a>
+
</footer>
+
<!--numbr 1:  i cant get your toplink to show number 2: weird grey bar on toggle menue-->
+
</body>
+
 
+
</html>
+

Revision as of 01:35, 20 October 2016