Difference between revisions of "Team:RHIT/GeneticsModel"

 
(60 intermediate revisions by 3 users not shown)
Line 2: Line 2:
 
<html>
 
<html>
  
<h1>Genetics Model </h1>
+
<script>
<p>The DNA coding for the 6 enzymes required for breakdown and assimilation of PET was too long to fit on one plasmid. To rectify this and to be able to test smaller subsystems, the PETase and MHETase genes were placed on Backbone 1, Plasmid 1. The Glycolaldehyde Reductase, Glycolaldehyde Dehydrogenase, and Glycolate Oxidase, and Malate Synthase were placed in sequence on Backbone 2, Plasmid 2.</p>
+
 
 +
$(document).ready(function(){
 +
    $("#gen").click(function(){
 +
        $("#genpic").slideToggle("slow");
 +
    });
 +
});
 +
</script>
 +
 
 +
<div class = "column third_size"></div>
 +
<div class = "column third_size">
 +
<center>
 +
<h4 style="cursor:pointer" id="gen">Genetics Model </h4>
 +
</center>
 +
</div>
 +
<div class="column third_size">
 +
<div class="popup" onclick="PopUp()"> <img src="https://static.igem.org/mediawiki/2018/8/8e/T--RHIT--PebbleHelp.jpg" style="width:40%">
 +
  <span class="popuptext" id="myPopup" style="position:absolute; top:20px; left:300px; width:500px">  Click the "Genetics Model" title to toggle our visual
 +
</span>
 +
 
 +
</div>
 +
</div>
 +
<div class = "column full_size">
 +
<center>
 +
<img style="display:none" id="genpic" src = "https://static.igem.org/mediawiki/2018/0/0c/T--RHIT--GenModPic.png">
 +
</center>
 +
 
 +
<p>The DNA coding for the 6 enzymes required for breakdown and assimilation of PET was too long to fit on one plasmid. To rectify this and to be able to test smaller subsystems, the PETase and MHETase genes were placed on Backbone 1, Plasmid 1. The Glycolaldehyde Reductase, Glycolaldehyde Dehydrogenase, Glycolate Oxidase, and Malate Synthase were placed in sequence on Backbone 2, Plasmid 2.</p>
 +
</div>
 
<div class = "clear extra_space"></div>
 
<div class = "clear extra_space"></div>
  
 +
<div class = "column full_size">
 
<h1>Plasmid 1</h1>
 
<h1>Plasmid 1</h1>
 +
</div>
 
<div class="column half_size">
 
<div class="column half_size">
 
<div class="container">
 
<div class="container">
         <img src="https://static.igem.org/mediawiki/2018/thumb/e/ea/T--RHIT--Plasmid1off.png/800px-T--RHIT--Plasmid1off.png"  
+
         <img src="https://static.igem.org/mediawiki/2018/thumb/d/d4/T--RHIT--Plasmid1off2.png/800px-T--RHIT--Plasmid1off2.png"  
 
         class="image">
 
         class="image">
 
         <div class="overlay">
 
         <div class="overlay">
Line 18: Line 47:
 
<div class="column half_size">
 
<div class="column half_size">
 
<div class="container">
 
<div class="container">
           <img src="https://static.igem.org/mediawiki/2018/d/d9/T--RHIT--Plasmid1on1.png" class="image">
+
           <img src="https://static.igem.org/mediawiki/2018/thumb/f/f2/T--RHIT--Plasmid1on2.png/800px-T--RHIT--Plasmid1on2.png" class="image">
 
           <div class="overlay">
 
           <div class="overlay">
           <div class="text"> Sir Richard Henhathel's Activated System </div>
+
           <div class="text"> Activated System </div>
 
           </div></div>
 
           </div></div>
<br /><br /><br /><br /><br />
+
<br /><br /><br /><br />
 
</div>
 
</div>
  
<img src="https://static.igem.org/mediawiki/2018/5/53/T--RHIT--Plasmid1Parameters1.png">
 
  
<p>Plasmid 1 uses an AraC and pBAD promoter to regulate expression of PETase and MHETase. The transcription factor made from AraC usually binds and represses the pBAD promoter, halting transcription of the plasmid. The inducer, Arabinose, can be added to the media, and this molecule binds to the AraC transcription factor on the DNA strand and changes it conformation so that transcription can occur [bmcsys]. The reaction scheme on the left explains a more complete mechanism of the transcription/translation of these proteins. The creation of AraC protein is related to a constitutive promoter which we assume enters the system as a constant rate, K. This method was also used by the UC Davis team in 2012. We assumed fast-equilibrium hypothesis on the formation of the dimer and that there is essentially a constant pool of arabinose in the environment. We can also streamline the binding of the two arabinose to the AraC dimer into one reaction determined by the rate parameters k3+ and k3-. Since the amount of RNA polymerase does not change relative to these molecules and the frequent assumption used literature to group the transcription and translation rate into one overall rate of protein production, we simplify the system further. From these assumptions, we can simplify the system down into the system shown on the right.</p>
+
<div class = "column half_size">
 +
<img src="https://static.igem.org/mediawiki/2018/5/53/T--RHIT--GeneticsParameters.png">
 +
</div>
 +
 
 +
<div class = "column full_size">
 +
<p>Plasmid 1 uses an AraC and pBAD promoter to regulate the expression of PETase and MHETase. The transcription factor made from AraC usually binds and represses the pBAD promoter, halting transcription of the plasmid. The inducer, arabinose, can be added to the media, and this molecule binds to the AraC transcription factor on the DNA strand and changes its conformation so that transcription can occur [1]. The reaction scheme on the left explains a more complete mechanism of the transcription/translation of these proteins. The creation of AraC protein is related to a constitutive promoter which we assume enters the system as a constant rate, K. This method was also used by the UC Davis iGEM team in 2012. We assumed a fast-equilibrium hypothesis on the formation of the dimer and an essentially constant pool of arabinose in the environment. We can also streamline the binding of the two arabinose inducers to the AraC dimer into one reaction determined by the rate parameters k3+ and k3-. Since the amount of RNA polymerase does not change relative to these molecules, and since the frequent assumption used in literature is to group the transcription and translation rate into one overall rate of protein production, we simplify the system further. From these assumptions, we can simplify the system down into the system shown on the right.</p>
  
 
<center>
 
<center>
<img src="https://static.igem.org/mediawiki/2018/9/96/T--RHIT--Plasmid1Eq.png">
+
<img src="https://static.igem.org/mediawiki/2018/9/96/T--RHIT--Plasmid1Eq.png" style="width:75%">
 
</center>
 
</center>
  
Line 36: Line 69:
  
 
<center>
 
<center>
<img src="https://static.igem.org/mediawiki/2018/d/d8/T--RHIT--Plasmid1Eqnew.png">
+
<img src="https://static.igem.org/mediawiki/2018/a/a3/T--RHIT--EQGEN.png" style="width:50%">
 
</center>
 
</center>
<p> <a href = "https://bmcsystbiol.biomedcentral.com/articles/10.1186/1752-0509-5-111">
 
https://bmcsystbiol.biomedcentral.com/articles/10.1186/1752-0509-5-111 </a> useful description and numbers for AraC promoter</p> </a>
 
  
<div class = "clear extra_space"></div>
 
<h1>Plasmid 2</h1>
 
<div class="column full_size">
 
<div class="container">
 
          <img src="https://static.igem.org/mediawiki/2018/thumb/0/09/T--RHIT--Plasmid2off.png/800px-T--RHIT--Plasmid2off.png"
 
            class = "image">
 
          <div class="overlay">
 
          <div class="text"> Repressed System </div>
 
          </div></div>
 
</div>
 
 
<div class="column full_size">
 
<div class="container">
 
          <img src="https://static.igem.org/mediawiki/2018/thumb/5/59/T--RHIT--Plasmid2on.png/800px-T--RHIT--Plasmid2on.png"
 
            class="image">
 
          <div class="overlay">
 
          <div class="text"> Activated System </div>
 
          </div></div>
 
</div>
 
 
<img src="https://static.igem.org/mediawiki/2018/d/d0/T--RHIT--Plasmid2Parameters1.png">
 
 
<p>Plasmid 2 uses a LacI and pTrc promoter is its backbone, which is requires a slightly different model structure from the other promoter. The transcription factor constitutively expressed by the LacI sequence creates a homotetramer and binds to the pTrc sequence repressing transcription of the Reductase, Dehydrogenase, Oxidase, and Synthase genes. We used a similar expression equation for this protein to the repressor protein in plasmid 1. However when IPTG, a lactose analog, is added to the cells, two of these molecule binds to the LacI protein causing it to change conformation and fall off the promoter. This then allows RNA polymerase to bind and start transcription. We used the same assumption as in plasmid 1, grouping the transcription and translation rate into one rate parameter, which simplified to the system on the right. We also grouped the binding of 2 IPTG to the LacI tetramer into one reaction instead of two separate events. Finally, since we are uncertain about the ordering of binding or the rate parameters of the formation of the LacI tetramer, we kept it with separate rate parameter for each subunit binding. This was important to include as the tetramer is the molecule that actually represses the genes. The last thing to note is that for there to be a functional Glycolaldehyde Reductase protein it needs two copies of itself to dimerize. This explains the factor multiplying its transcription rate.</p>
 
  
 +
<h3> Simulations: </h3>
 
<center>
 
<center>
<img src="https://static.igem.org/mediawiki/2018/1/13/T--RHIT--Plasmid2Rxn.png">
+
Table 1. List of values used in simulation parameters<br />
 +
<img src = "https://static.igem.org/mediawiki/2018/4/49/T--RHIT--GenModParam.png" style="width:60%">
 
</center>
 
</center>
 +
<p>Figure 1 maps out the expression level of the free repressor protein AraC and Plasmid 1 genes, PETase and MHETase over 12 hours. The population of cells in the simulation was 1e+08. While free AraC protein, which would normally suppress expression, is relatively low there is a high expression of the objective proteins. The MHETase and PETase plots both follow the same trend, although, MHETase does decay slower than PETase and also reaches a higher maximum concentration. </p>
 +
<center>
 +
<img src = "https://static.igem.org/mediawiki/2018/c/c5/T--RHIT--GenModGraphs.png" style="width:70%"><br />
 +
Figure 1. Expression of AraC promoter in presence of Arabinose in Plasmid 1<br /><br /><br />
 +
</center>
 +
<h3> Sensitivity Analysis: </h3>
 +
<p>In Table 1, there are three values that say they are estimated in the calculated column. This is due to the lack of comprehensive material on them, and that these values create realistic behavior of the system. There is a small range in each value where the system remains stable and biologically relevant, so these values were pulled from that range. The behavior of the system throughout the range is consistent with only slight variation on amounts. </p>
  
<h3>Model Equations:</h3>
+
<h3>References:</h3>
 
+
<ul>
<center>
+
<li>[1] D. Madar, E. Dekel, A. Bren and U. Alon. “Negative auto-regulation increases the input dynamic-range of the arabinose system of Escherichia coli,” BMC Systems Biology. 2011. [Online]. https://bmcsystbiol.biomedcentral.com/articles/10.1186/1752-0509-5-111 </li>
<img src="https://static.igem.org/mediawiki/2018/b/bb/T--RHIT--Plasmid2EqNew1.png">
+
<li>[2] L. Wang, Y. J. Zhou, D. Ji, and Z.K. Zhao, “An accurate method for estimation of the intracellular aqueous volume of Escherichia coli cells,” Journal of Microbiological Methods, 2013, p. 8. [Online].http://bionumbers.hms.harvard.edu/bionumber.aspx?id=108813&ver=3&trm=e%20coli%20cell%20volume&org</li>
</center>
+
<li>[3] H. Bremer and P. P. Dennis, “Modulation of chemical composition and other parameters of the cell by growth rate.” Escherichia coli and Salmonella typhimurium: Cellular and Molecular Biology, 2nd ed. [Online]. Neidhardt, et al. eds. chapter 97, 1991, p. 1559. [Online]. http://bionumbers.hms.harvard.edu/bionumber.aspx?&id=100059&ver=39</li>
 +
<li>[4] X. Zhang, T. Reeder, and R. Schleif. Transcription Activation Parameters at ara pBAD.” Journal of Molecular Biology. Vol 258, 1996, p. 14-24. [Online]. http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.417.6037&rep=rep1&type=pdf</li>
 +
<li>[5] D. Kolodrubetz and R. Schleif. “Identifical of araC protein on two-dimensional gels, its in vivo instability and normal level.” Journal of Molecular Biology. Vol 149, issue 1, pp.133-139, 1981. [Online].  https://www.sciencedirect.com/science/article/pii/0022283681902655 </li>
 +
<li>[6] J.A. Mergerle, F. Georg, et al. “Timing and Dynamics of Single Cell Gene Expression in the Arabinose Utilization System.” Biophysical Journal. Vol. 95, Issue 4, pp.2103-2115, Aug 2008. [Online]. https://www.sciencedirect.com/science/article/pii/S0006349508701681 </li>
 +
<li>[7] R. Schleif. “AraC protein, regulation of the L-arabinose operon in Escherichia coli, and the light switch mechanism of AraC action.” FEMS Microbiology Review. vol. 34, issue 5, September 2010, pp. 779-796. [Online]. https://academic.oup.com/femsre/article/34/5/779/797770</li>
 +
<li>[8] K. Martin, L. Huo, and R.Schleif. “The DNA Loop model for ara repression: AraC Protein occupies the proposed loop sites in vivo and repression-negative mutations lie in these same sites.” In USA Proceedings National Academy of Science. vol. 82, June 1986, pp. 3654-3658. [Available: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC323581/pdf/pnas00315-0095.pdf]</li>
 +
<li>[9] “Part: pSB1C3.”iGEM Registry of Standard Biological Parts. Sept 2008. [Online]. http://parts.igem.org/Part:pSB1C3</li>
 +
</ul>
 +
</div>
 +
</html>

Latest revision as of 20:22, 2 August 2018




Genetics Model

The DNA coding for the 6 enzymes required for breakdown and assimilation of PET was too long to fit on one plasmid. To rectify this and to be able to test smaller subsystems, the PETase and MHETase genes were placed on Backbone 1, Plasmid 1. The Glycolaldehyde Reductase, Glycolaldehyde Dehydrogenase, Glycolate Oxidase, and Malate Synthase were placed in sequence on Backbone 2, Plasmid 2.

Plasmid 1

Repressed System
Activated System




Plasmid 1 uses an AraC and pBAD promoter to regulate the expression of PETase and MHETase. The transcription factor made from AraC usually binds and represses the pBAD promoter, halting transcription of the plasmid. The inducer, arabinose, can be added to the media, and this molecule binds to the AraC transcription factor on the DNA strand and changes its conformation so that transcription can occur [1]. The reaction scheme on the left explains a more complete mechanism of the transcription/translation of these proteins. The creation of AraC protein is related to a constitutive promoter which we assume enters the system as a constant rate, K. This method was also used by the UC Davis iGEM team in 2012. We assumed a fast-equilibrium hypothesis on the formation of the dimer and an essentially constant pool of arabinose in the environment. We can also streamline the binding of the two arabinose inducers to the AraC dimer into one reaction determined by the rate parameters k3+ and k3-. Since the amount of RNA polymerase does not change relative to these molecules, and since the frequent assumption used in literature is to group the transcription and translation rate into one overall rate of protein production, we simplify the system further. From these assumptions, we can simplify the system down into the system shown on the right.

Model Equations:

Simulations:

Table 1. List of values used in simulation parameters

Figure 1 maps out the expression level of the free repressor protein AraC and Plasmid 1 genes, PETase and MHETase over 12 hours. The population of cells in the simulation was 1e+08. While free AraC protein, which would normally suppress expression, is relatively low there is a high expression of the objective proteins. The MHETase and PETase plots both follow the same trend, although, MHETase does decay slower than PETase and also reaches a higher maximum concentration.


Figure 1. Expression of AraC promoter in presence of Arabinose in Plasmid 1


Sensitivity Analysis:

In Table 1, there are three values that say they are estimated in the calculated column. This is due to the lack of comprehensive material on them, and that these values create realistic behavior of the system. There is a small range in each value where the system remains stable and biologically relevant, so these values were pulled from that range. The behavior of the system throughout the range is consistent with only slight variation on amounts.

References:

  • [1] D. Madar, E. Dekel, A. Bren and U. Alon. “Negative auto-regulation increases the input dynamic-range of the arabinose system of Escherichia coli,” BMC Systems Biology. 2011. [Online]. https://bmcsystbiol.biomedcentral.com/articles/10.1186/1752-0509-5-111
  • [2] L. Wang, Y. J. Zhou, D. Ji, and Z.K. Zhao, “An accurate method for estimation of the intracellular aqueous volume of Escherichia coli cells,” Journal of Microbiological Methods, 2013, p. 8. [Online].http://bionumbers.hms.harvard.edu/bionumber.aspx?id=108813&ver=3&trm=e%20coli%20cell%20volume&org
  • [3] H. Bremer and P. P. Dennis, “Modulation of chemical composition and other parameters of the cell by growth rate.” Escherichia coli and Salmonella typhimurium: Cellular and Molecular Biology, 2nd ed. [Online]. Neidhardt, et al. eds. chapter 97, 1991, p. 1559. [Online]. http://bionumbers.hms.harvard.edu/bionumber.aspx?&id=100059&ver=39
  • [4] X. Zhang, T. Reeder, and R. Schleif. Transcription Activation Parameters at ara pBAD.” Journal of Molecular Biology. Vol 258, 1996, p. 14-24. [Online]. http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.417.6037&rep=rep1&type=pdf
  • [5] D. Kolodrubetz and R. Schleif. “Identifical of araC protein on two-dimensional gels, its in vivo instability and normal level.” Journal of Molecular Biology. Vol 149, issue 1, pp.133-139, 1981. [Online]. https://www.sciencedirect.com/science/article/pii/0022283681902655
  • [6] J.A. Mergerle, F. Georg, et al. “Timing and Dynamics of Single Cell Gene Expression in the Arabinose Utilization System.” Biophysical Journal. Vol. 95, Issue 4, pp.2103-2115, Aug 2008. [Online]. https://www.sciencedirect.com/science/article/pii/S0006349508701681
  • [7] R. Schleif. “AraC protein, regulation of the L-arabinose operon in Escherichia coli, and the light switch mechanism of AraC action.” FEMS Microbiology Review. vol. 34, issue 5, September 2010, pp. 779-796. [Online]. https://academic.oup.com/femsre/article/34/5/779/797770
  • [8] K. Martin, L. Huo, and R.Schleif. “The DNA Loop model for ara repression: AraC Protein occupies the proposed loop sites in vivo and repression-negative mutations lie in these same sites.” In USA Proceedings National Academy of Science. vol. 82, June 1986, pp. 3654-3658. [Available: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC323581/pdf/pnas00315-0095.pdf]
  • [9] “Part: pSB1C3.”iGEM Registry of Standard Biological Parts. Sept 2008. [Online]. http://parts.igem.org/Part:pSB1C3