Team:Consort Alberta/project

From 2014hs.igem.org

(Difference between revisions)
m
 
(28 intermediate revisions not shown)
Line 164: Line 164:
<p class="auto-style18" style="height: 140px; mso-bidi-font-size: 10.0pt">&nbsp;&nbsp;&nbsp;&nbsp;  
<p class="auto-style18" style="height: 140px; mso-bidi-font-size: 10.0pt">&nbsp;&nbsp;&nbsp;&nbsp;  
-
ECOS (<u>E</u>nvironmental <u>CO</u>ntaminant <u>S</u>ensor) is a biological sensor theorized to detect the presence of aromatic hydrocarbons by producing an indicator protein when exposed to xylene. We chose xylene as our trigger because it is very well correlated with the presence of other more dangerous and carcinogenic compounds such as benzene and its derivatives. The danger of this class of compounds is in the fact that they can intercalate into DNA and cause a variety of mutations. Our biobrick will, in effect, create different amounts of fluorescent protein in the presence of xylene bound to XylR- which the first portion of our plasmid will constantly produce. As the economy in our rural community is based largely on agriculture and oil and gas industries, oil spills have detrimental effects on the environment, economy, and general health. To be able to detect the presence of these aromatic hydrocarbons on site would great benefit all members of our community. This year we will be testing out three different indicators to allow options concerning scale and intensity of colour change. In future years, our project will allow early identification of contamination will facilitate rapid clean-up and minimize health risks to members of our community and to the consumers who rely on the food we produce.</p>
+
ECOS (<u>E</u>nvironmental <u>CO</u>ntaminant <u>S</u>ensor) is a biological sensor theorized to detect the presence of aromatic hydrocarbons by producing an indicator protein when exposed to xylene. We chose xylene as our trigger because it is very well correlated with the presence of other more dangerous and carcinogenic compounds such as benzene and its derivatives. The danger of this class of compounds is in the fact that they can intercalate into DNA and cause a variety of mutations. Our biobrick will, in effect, create different amounts of fluorescent protein in the presence of xylene bound to XylR- which the first portion of our plasmid will constantly produce. As the economy in our rural community is based largely on agriculture and oil and gas industries, oil spills have detrimental effects on the environment, economy, and general health. To be able to detect the presence of these aromatic hydrocarbons on site would greatly benefit all members of our community. This year we will be testing out three different indicators to allow options concerning scale and intensity of colour change. In future years, our project will allow early identification of contamination will facilitate rapid clean-up and minimize health risks to members of our community and to the consumers who rely on the food we produce.</p>
Line 178: Line 178:
<p class="auto-style42" style="height: 40px"><strong><em><a name="The Science...">The Science...</a></em></strong></p>
<p class="auto-style42" style="height: 40px"><strong><em><a name="The Science...">The Science...</a></em></strong></p>
<p class="auto-style18" style="height: 120px; mso-bidi-font-size: 10.0pt">
<p class="auto-style18" style="height: 120px; mso-bidi-font-size: 10.0pt">
-
&nbsp;&nbsp;&nbsp;&nbsp;Our resulting biobrick consists of two main portions. The first portion of our plasmid will be responsible for producing the protein XylR.  This is essential because the bonding of m-xylene undergoes a conformational change when bonded to XylR which then allows it to bind to, and positively regulate the Pu promoter. At a 4:2 XylR to Xylene ratio the bonding starts the production of our reporter protein. The second portion of our plasmid expresses the indicator protein when in the presence of bonded XylR and m-xylene and essentially allow us to know if m-xylene is present in our soil sample or not. The structure of this biobrick is shown below.</p>
+
&nbsp;&nbsp;&nbsp;&nbsp;Our resulting biobrick consists of two main portions. The first portion of our plasmid will be responsible for producing the protein XylR.  This is essential because the bonding of m-xylene undergoes a conformational change when bonded to XylR which then allows it to bind to, and positively regulate, the Pu promoter. At a 4:2 XylR to Xylene ratio the bonding starts the production of our reporter protein. The second portion of our plasmid expresses the indicator protein when in the presence of bonded XylR and m-xylene and essentially allow us to know if m-xylene is present in our soil sample or not. The structure of this biobrick is shown below.</p>
<p class="auto-style18" style="height: 0px; mso-bidi-font-size: 10.0pt">&nbsp;&nbsp;&nbsp;&nbsp;  
<p class="auto-style18" style="height: 0px; mso-bidi-font-size: 10.0pt">&nbsp;&nbsp;&nbsp;&nbsp;  
Line 189: Line 189:
         I723020 - This is the Pu promoter.</li><li>
         I723020 - This is the Pu promoter.</li><li>
         B0030 - RBS.1 strong.</li><li>
         B0030 - RBS.1 strong.</li><li>
-
         Our indicators which we will be testing so that we can see which indicator for our project will be the best:</li><ol class="auto-style18" style="margin-bottom: 5px"><li>
+
         Reporter:</li><ol class="auto-style18" style="margin-bottom: 5px"><li>
               E0040 - GFP wild-type, will glow when exposed to UV light.</li><li>
               E0040 - GFP wild-type, will glow when exposed to UV light.</li><li>
-
               K592009 - amilCFP, This chromoprotein from the coral Acropora millepora, naturally exhibits strong color when expressed.</li><li>
+
               K592009 - amilCP, This chromoprotein from the coral Acropora millepora, naturally exhibits strong color when expressed.</li><li>
               I732005 - lacZ, Beta-galactosidase cleaves X-gal and ONPG into colorful products.</li></ol>
               I732005 - lacZ, Beta-galactosidase cleaves X-gal and ONPG into colorful products.</li></ol>
         <li>B0015 - lastly another double stop codon.</li><li>
         <li>B0015 - lastly another double stop codon.</li><li>
         pSB1C3 - Our plasmid backbone which includes the resistance towards the antibiotic, chloramphenicol.</li></ol></p>
         pSB1C3 - Our plasmid backbone which includes the resistance towards the antibiotic, chloramphenicol.</li></ol></p>
 +
 +
<!--<img alt="" class="auto-style24" height="140" src="https://static.igem.org/mediawiki/2014hs/5/5d/Sequence.JPG" width="700" />-->
 +
 +
<p align="middle"><img alt=""  class="auto-style24" height="135" src="https://static.igem.org/mediawiki/2014hs/thumb/5/5d/Sequence.JPG/800px-Sequence.JPG" width="900"></p>
 +
Line 208: Line 213:
<p class="auto-style42" style="height: 40px"><em><strong><a name="Mathematical_Modelling...">Mathematical Modelling...</a></em></strong></p>
<p class="auto-style42" style="height: 40px"><em><strong><a name="Mathematical_Modelling...">Mathematical Modelling...</a></em></strong></p>
 +
<p class="auto-style18" style="height: 90px; mso-bidi-font-size: 10.0pt">&nbsp;&nbsp;&nbsp;&nbsp; Making a sensor involves more than just the creation of the basic biological circuit, as there could be many factors that might improve the performance of the system. An example of these could be the strength of the RBS, the spacing between the promoter and the reporter gene, or the copy number of the plasmid used to house the circuit. While we couldn't nearly address all of these points in our project, we did want to look into the optimization of our system. To do this, we turned to mathematical modelling to create a visual representation of our system.</p>
 +
<p class="auto-style18" style="height: 220px; mso-bidi-font-size: 10.0pt">&nbsp;&nbsp;&nbsp;&nbsp; We worked off of a basic framework published by Koutinas, M., <em>et al.</em>,<a class="nav" href="#References"><sup>[1,2]</sup></a> which was modelling the expression of the Ps promoter with XylR induction. The Ps promoter is the natural promoter found with the <em>xylR</em> gene, and the XylR protein can interact with both it and the Pu promoter we used in our project. Due to similarities between the Ps and Pu promoters, we assumed the deactivation rate of the two components were alike. For remaining values, we simply replaced the values of the Ps promoter with known values of the Pu promoter, keeping constants reported for XylR the same. In the end, we created five equations to represent the action of our system. Our model shows our biobrick producing the protein XylR, binding with xylene and the relationship to the output of our respective proteins. Each of our formulas are constructed as a rate in which a concentration or output is given of a specific substance in respect to time- otherwise known as a derivative. We will now expand upon the equations we modified and used.</p>
-
<p class="auto-style18" style="height: 190px; mso-bidi-font-size: 10.0pt">&nbsp;&nbsp;&nbsp;&nbsp;
 
-
In order to create a visual representation of our output, we turned to mathematical modelling. We created a modification of a model done by Koutinas, M., <em>et al.</em><a class="nav" href="#References"><sup>[1]</sup></a>  which was modelling the expression of the Ps promoter. The Ps promoter is the natural promoter found with the <em>xylR</em> gene. Due to similarities between the Ps and Pu promoters, we assumed the deactivation rate of the two components were alike. For remaining values, we simply replaced the values of the Ps promoter with known values of the Pu promoter. <a class="nav" href="#References"><sup>[2,3]</sup></a> In the end, we created five equations to represent the action of our biobrick. Our model shows our biobrick producing the protein XylR, binding with xylene and the relationship to the output of our respective proteins. Each of our formulas are constructed as a rate in which a concentration or output is given of a specific substance in respect to time- otherwise known as a derivative. We will now expand upon the equations we modified and used. </p>
 
-
<img alt="" class="auto-style24" height="100" src="https://static.igem.org/mediawiki/2014hs/2/27/ModellingEquation1-1.PNG" width="530" /></p>
+
<p align="middle"><img alt="" class="auto-style24" height="100" src="https://static.igem.org/mediawiki/2014hs/2/27/ModellingEquation1-1.PNG" width="530" /></p>
<p class="auto-style18" style="height: 80px; mso-bidi-font-size: 10.0pt">
<p class="auto-style18" style="height: 80px; mso-bidi-font-size: 10.0pt">
-
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; The first equation above represents the translation of RNA into the protein XylR. RNA is the representation of the rate of DNA being translated into the protein XylR. P represents the number of plasmids inside a cell. J23100<sub>TC</sub> represents the activity of our promoter. t<sub>a</sub> is the transcription rate of the <em>xylR</em> into RNA, and K<sub>RNAdeg</sub> represents the degradation rate in the E.Coli.</p>
+
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; The first equation above represents the transcription of DNA into the XylR RNA transcript. RNA is the representation of the transcript, with P representing the number of plasmids inside a cell. J23100<sub>TC</sub> represents the activity of our promoter controlling XylR. ta is the transcription rate of the xylR into RNA, based on the transcription rate of <i>E. Coli</i> and the size of XylR, and K<sub>RNAdeg</sub> represents the degradation rate in the <i>E. Coli</i>.</p>
-
<img alt="" class="auto-style24" height="140" src="https://static.igem.org/mediawiki/2014hs/f/f4/ModellingEquation2-1%263-1.PNG" width="600" /></p>
+
<p align="middle"><img alt="" class="auto-style24" height="140" src="https://static.igem.org/mediawiki/2014hs/f/f4/ModellingEquation2-1%263-1.PNG" width="600" /></p>
-
<p class="auto-style18" style="height: 140px; mso-bidi-font-size: 10.0pt">
+
<p class="auto-style18" style="height: 80px; mso-bidi-font-size: 10.0pt">
-
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; The second and third equations represent the inactive and active forms of XylR- or the concentration of XylR that has not bonded to xylene and the concentration of XylR that has bonded to xylene. XylR<sub>i</sub> and XylR<sub>a</sub> refer to the concentrations of the inactive and active forms of XylR. RNA is the substitution of our first equation. t<sub>r</sub> is the translation rate of our RNA into XylR. r<sub>XylR</sub> is the oligomerization constant of XylR, r<sub>R,XylR</sub> is the dissociation constant of inactive XylR and xyl in the total concentration of xylene. &alpha;<sub>XylRi</sub> accounts for degradation and dilution because of cellular volume increase. Xylene represents the concentration of xylene present.</p>
+
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; The second and third equations represent the inactive and active forms of XylR- or the concentration of XylR that has not bonded to xylene (XylR<sub>i</sub>) and the concentration of XylR that has bonded to xylene (XylR<sub>a</sub>). RNA is the substitution of our first equation. t<sub>r</sub> is the translation rate of our RNA into XylR. r<sub>XylR</sub> is the oligomerization constant of XylR, r<sub>R,XylR</sub> is the dissociation constant of active XylR, and xyl in the total concentration of xylene. &alpha;<sub>XylRi</sub> accounts for degradation of XylR.</p>
-
<img alt="" class="auto-style24" height="70" src="https://static.igem.org/mediawiki/2014hs/1/11/ModellingEquation4-1.PNG" width="700" /></p>
+
<p align="middle"><img alt="" class="auto-style24" height="70" src="https://static.igem.org/mediawiki/2014hs/1/11/ModellingEquation4-1.PNG" width="700" /></p>
-
<p class="auto-style18" style="height: 110px; mso-bidi-font-size: 10.0pt">
+
<p class="auto-style18" style="height: 93px; mso-bidi-font-size: 10.0pt">
-
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; The fourth equation represents the concentration of xylene present. XylR<sub>i</sub> and XylR<sub>a</sub> refer to the concentrations of the inactive and active forms of XylR. r<sub>XylR</sub> is the oligomerization constant of XylR, r<sub>R,XylR</sub> is the dissociation constant of inactive XylR and xyl in the total concentration of xylene. &alpha;<sub>XylRi</sub> accounts for degradation and dilution because of cellular volume increase. Xylene represents the concentration of xylene present. We multiplied by seconds per hour to get a domain more realistic to our needs.</p>
+
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; The fourth equation represents the concentration of xylene present. XylR<sub>i</sub> and XylR<sub>a</sub> refer to the concentrations of the inactive and active forms of XylR. r<sub>XylR</sub> is the oligomerization constant of XylR, r<sub>R,XylR</sub> is the dissociation constant of inactive XylR and xyl is the total concentration of xylene. α<sub>XylRi</sub> accounts for the degradation of XylR<sub>a</sub>, which would release the xylene it held. We multiplied by seconds per hour to get a domain more realistic to our needs.</p>
-
<img alt="" class="auto-style24" height="88" style="top-margin: 5px" src="https://static.igem.org/mediawiki/2014hs/8/86/ModellingEquation5.JPG" width="600" /></p>
+
<p align="middle"><img alt="" class="auto-style24" height="88" style="top-margin: 5px" src="https://static.igem.org/mediawiki/2014hs/8/86/ModellingEquation5.JPG" width="600" /></p>
<p class="auto-style18" style="height: 80px; mso-bidi-font-size: 10.0pt">
<p class="auto-style18" style="height: 80px; mso-bidi-font-size: 10.0pt">
-
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; The fifth equation represents the output of the pu promoter. Pu<sub>TC</sub> represents the activity of the Pu promoter or the output of protein. &alpha;<sub>pu</sub> is the deactivation rate of the Pu promoter. K<sub>XylRa</sub> is the repression coefficients of the Pu due to binding and n<sub>pr,i</sub> and n<sub>pr,a</sub> are the hill coefficients of the Pu promoter. &beta;<sub>0</sub> and &beta;<sub>PS</sub> are the basal and maximal expression respectively of the Pu promoter.</p>
+
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; The fifth equation represents the output of the pu promoter. Pu<sub>TC</sub> represents the activity of the Pu promoter or the output of protein. &alpha;<sub>pu</sub> is the deactivation rate of the Pu promoter. K<sub>XylRa</sub> is the activation coefficients of the Pu due to binding and n<sub>ps,a</sub> is the hill coefficient for the interaction. &beta;<sub>0</sub> and &beta;<sub>PS</sub> are the basal and maximal expression, respectively, of the Pu promoter.</p>
-
<p class="auto-style18" style="height: 80px; mso-bidi-font-size: 10.0pt">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; In order to create an output, we used the program <a class="nav" href="http://www.scilab.org/">Scilab</a>. By inputting these equations into code we were able to have the program calculate our estimates and graph the results. In the end our output was around 5.5 mPoPS. This is a logical output and tells us our biobrick will successfully over express the protein. Below is a copy of our code that we had used in our program. The resulting graph is the output of our Pu promoter.</p>
+
<p class="auto-style18" style="height: 150px; mso-bidi-font-size: 10.0pt">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;   In order to create an output, we used the program <a class="nav" href="http://www.scilab.org/" target="_blank">Scilab</a>. By inputting these equations into code we were able to have the program calculate our estimates and graph the results. In the end our output was around 5.5 mPoPS. This is a logical output and tells us our biobrick will successfully over express the protein. We set the output of our model in PoPS due to our system's flexibility when it comes to our reporter protein. This allows us to quickly add in an extra equation to turn our current output into fluorescence, LacZ output, or pigment production without having to worry about editing other parts of our model. Below is a copy of our code that we had used in our program. The resulting graph is the output of our Pu promoter.</p>
 +
<p align="middle"><img alt=""  class="auto-style24" height="380" src="https://static.igem.org/mediawiki/2014hs/3/3d/E5_Output_Graph.PNG" width="450">
 +
<img alt=""  class="auto-style24" height="380" src="https://static.igem.org/mediawiki/2014hs/4/4b/All_Code.PNG" width="450"></p>
 +
 +
<p class="auto-style18" style="height: 110px; mso-bidi-font-size: 10.0pt">&nbsp;&nbsp;&nbsp;&nbsp; While we have yet to apply our modelling framework to the optimization of our sensor, we have created the foundation of a model that can be easily tweaked to test the effect of different variables on sensor performance. Through this we aim to do a sensitivity analysis of the various factors involved in our model in the future so that we can predict what to change for version 2 of our sensor. From there we can return to our model, feeding new data into it to create many future versions of our system, improving performance each time.</p>
<!--PROTOTYPING-->
<!--PROTOTYPING-->
-
<p class="auto-style42" style="height: 0px"><strong><em><a name="Protoyping...">Prototyping...</a></em></strong></p>
+
<p class="auto-style42" style="height: 50px"><strong><em><a name="Protoyping...">Prototyping...</a></em></strong></p>
 +
 
 +
<p class="auto-style18" style="height: 90px; mso-bidi-font-size: 10.0pt">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ECOS was meant to be more than just a theory; more than just a project to play out at school- we wanted ECOS to have real world applications by the time we finished making our construct. We had one major problem: How were we going to test our product in an immediate and effective manner? To answer this we designed two prototypes to design and test. Unfortunately we were unable to finish the wet-lab portion of this project, but the designing is complete thanks to the Maya program via Autodesk. These prototypes are beneficial because they 1) are cheap to build and 2) relatively quick to test. Here are how they turned out:</p>
<p class="auto-style42" style="height: 19px">&nbsp;</p>
<p class="auto-style42" style="height: 19px">&nbsp;</p>
-
<img alt="" src="https://static.igem.org/mediawiki/2014hs/e/e9/Photo_1.GIF" /></p>
+
<img alt="" style="float: left; margin-right: 20px" src="https://static.igem.org/mediawiki/2014hs/e/e9/Photo_1.GIF" /></p>
-
<p class="auto-style18" style="height: 100px; mso-bidi-font-size: 10.0pt"> &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
+
<p class="auto-style18" style="height: 180px; mso-bidi-font-size: 10.0pt"> &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
-
PROTOTYPE #1: We designed our first prototype so that it has a positive pressure source which creates a current for our xylene to travel from our sample, which is in a heated container so that the xylene is vapourized, to our ECOS container in which we bubble it through the E. Coli culture to get results. This was a plausible idea, as it would require fairly non-expensive materials and could be maintained by the average business person. However, we thought it could be more efficient in it's design, and so we kept looking for ideas. Our biggest problem with this was the fact that it is a relatively complex prototype. So we challenged ourselves to try and create a simpler prototype.</p>
+
PROTOTYPE #1: We designed our first prototype so that it has a positive pressure source which creates a current for our xylene to travel from our sample, which is in a heated container so that the xylene is vapourized, to our ECOS container in which we bubble it through the <i>E. Coli</i> culture to get results. This was a plausible idea, as it would require fairly non-expensive materials and could be maintained by the average business person. However, we thought it could be more efficient in it's design, and so we kept looking for ideas. Our biggest problem with this was the fact that it is a relatively complex prototype. So we challenged ourselves to try and create a simpler prototype.</p>
<p class="auto-style42" style="height: 19px">&nbsp;</p>
<p class="auto-style42" style="height: 19px">&nbsp;</p>
-
<img alt="" src="https://static.igem.org/mediawiki/2014hs/d/dc/Photo.GIF" /></p>
+
<img style="float: right; margin-left: 5px; margin-bottom: 5px" alt="" src="https://static.igem.org/mediawiki/2014hs/d/dc/Photo.GIF" /></p>
-
<FONT COLOR=white><p class="auto-style18" style="height: 180px; mso-bidi-font-size: 10.0pt"> &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
+
<FONT COLOR=white><p class="auto-style18" style="height: 240px; mso-bidi-font-size: 10.0pt"> &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
PROTOTYPE #2: Another solution we theorized was to have our ECOS impregnated into alginate and then made into small spheres so that we can simply add them into a heated sample to get our results. In a container, we combine our soil sample and our alginate beads, mix them around, and then sift out the beads. The beads produce our indicator protein from the contamination levels of the soil. This design has been inspired by the <a class="nav" <a href="https://2013.igem.org/Team:Peking/Project/Devices#Milestone2" target="_blank"> Peking 2013 Team</a> who used alginate encapsulation beads. Their reasoning and justification for using alginate beads stood out to us, as they stated that it was stable and inexpensive and easy to shape and manipulate. Upon further research into the use of these beads, we discovered that this prototype could be more efficient to use than our first prototype that required a few different pieces to the system. This revised prototype would also be easier for the average oil worker to use compared to the larger prototype.</p>
PROTOTYPE #2: Another solution we theorized was to have our ECOS impregnated into alginate and then made into small spheres so that we can simply add them into a heated sample to get our results. In a container, we combine our soil sample and our alginate beads, mix them around, and then sift out the beads. The beads produce our indicator protein from the contamination levels of the soil. This design has been inspired by the <a class="nav" <a href="https://2013.igem.org/Team:Peking/Project/Devices#Milestone2" target="_blank"> Peking 2013 Team</a> who used alginate encapsulation beads. Their reasoning and justification for using alginate beads stood out to us, as they stated that it was stable and inexpensive and easy to shape and manipulate. Upon further research into the use of these beads, we discovered that this prototype could be more efficient to use than our first prototype that required a few different pieces to the system. This revised prototype would also be easier for the average oil worker to use compared to the larger prototype.</p>
<p class="auto-style42" style="height: 30px">&nbsp;</p>
<p class="auto-style42" style="height: 30px">&nbsp;</p>
 +
<p class="auto-style18" style="height: 240px; mso-bidi-font-size: 10.0pt"> &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;As time progresses with our project we plan to test and find what concentrations of m-Xylene trigger what kind of report. To do this we have purchased a couple containers of Xylene, each with a specified concentration, which we will apply to our prototypes. In-lab trials, in combination with our mathemantical model, will allow us to predict what [Xylene] is around and therefore what level of Benzene and Benzene derivatives are around as well. Once we manage to find this correlation between reporter and [Xylene] then we will be moving the project out to attempt some real-world application. Depending on our results here, ECOS may move away from rural Consort Alberta and out accross Canada and potentially the world. Of course there will always be things to fix- expecially when it comes to synthetic biology, so just like a computer programer we would try to fine tune and bring out newer version of this construct. The changes, of course, would come from our observations in-lab and customer suggestions as time progressed.</p>
 +
 +
<!--REFERENCES-->
<p class="auto-style41" style="height: 0px; mso-bidi-font-size: 10.0pt;">
<p class="auto-style41" style="height: 0px; mso-bidi-font-size: 10.0pt;">

Latest revision as of 03:52, 21 June 2014

Our Project

 Home | Our Project | Our Team | Human Practices | Our Notebook | Safety | Attributions | Sponsors| Contact Us

Our Project

   Introduction    Science   Comparison 

  

Mathematical Modelling    Prototyping

 

 

 

Introduction...

     ECOS (Environmental COntaminant Sensor) is a biological sensor theorized to detect the presence of aromatic hydrocarbons by producing an indicator protein when exposed to xylene. We chose xylene as our trigger because it is very well correlated with the presence of other more dangerous and carcinogenic compounds such as benzene and its derivatives. The danger of this class of compounds is in the fact that they can intercalate into DNA and cause a variety of mutations. Our biobrick will, in effect, create different amounts of fluorescent protein in the presence of xylene bound to XylR- which the first portion of our plasmid will constantly produce. As the economy in our rural community is based largely on agriculture and oil and gas industries, oil spills have detrimental effects on the environment, economy, and general health. To be able to detect the presence of these aromatic hydrocarbons on site would greatly benefit all members of our community. This year we will be testing out three different indicators to allow options concerning scale and intensity of colour change. In future years, our project will allow early identification of contamination will facilitate rapid clean-up and minimize health risks to members of our community and to the consumers who rely on the food we produce.

The Science...

    Our resulting biobrick consists of two main portions. The first portion of our plasmid will be responsible for producing the protein XylR. This is essential because the bonding of m-xylene undergoes a conformational change when bonded to XylR which then allows it to bind to, and positively regulate, the Pu promoter. At a 4:2 XylR to Xylene ratio the bonding starts the production of our reporter protein. The second portion of our plasmid expresses the indicator protein when in the presence of bonded XylR and m-xylene and essentially allow us to know if m-xylene is present in our soil sample or not. The structure of this biobrick is shown below.

     Parts:

  1. J23100 - which is a constituent promoter.
  2. B0034 - the RBS for our XylR gene.
  3. I723017 - the XylR coding region which encodes for the transcriptional regulator XylR protein.
  4. B0015 - the double stop codon for this sequence.
  5. I723020 - This is the Pu promoter.
  6. B0030 - RBS.1 strong.
  7. Reporter:
    1. E0040 - GFP wild-type, will glow when exposed to UV light.
    2. K592009 - amilCP, This chromoprotein from the coral Acropora millepora, naturally exhibits strong color when expressed.
    3. I732005 - lacZ, Beta-galactosidase cleaves X-gal and ONPG into colorful products.
  8. B0015 - lastly another double stop codon.
  9. pSB1C3 - Our plasmid backbone which includes the resistance towards the antibiotic, chloramphenicol.

Comparison...

     Where the comparison portion comes into play is the second part, which will be chosen from any one of three indicators: amilCP which when produced expresses a strong blue color, GFP when under a UV light source fluoresces green or, LacZ which when a substrate is added it cleaves it into strongly colored blue products. We will be looking at which indicator functions with the most efficiency and produces results of the highest clarity.

Mathematical Modelling...

     Making a sensor involves more than just the creation of the basic biological circuit, as there could be many factors that might improve the performance of the system. An example of these could be the strength of the RBS, the spacing between the promoter and the reporter gene, or the copy number of the plasmid used to house the circuit. While we couldn't nearly address all of these points in our project, we did want to look into the optimization of our system. To do this, we turned to mathematical modelling to create a visual representation of our system.

     We worked off of a basic framework published by Koutinas, M., et al.,[1,2] which was modelling the expression of the Ps promoter with XylR induction. The Ps promoter is the natural promoter found with the xylR gene, and the XylR protein can interact with both it and the Pu promoter we used in our project. Due to similarities between the Ps and Pu promoters, we assumed the deactivation rate of the two components were alike. For remaining values, we simply replaced the values of the Ps promoter with known values of the Pu promoter, keeping constants reported for XylR the same. In the end, we created five equations to represent the action of our system. Our model shows our biobrick producing the protein XylR, binding with xylene and the relationship to the output of our respective proteins. Each of our formulas are constructed as a rate in which a concentration or output is given of a specific substance in respect to time- otherwise known as a derivative. We will now expand upon the equations we modified and used.

      The first equation above represents the transcription of DNA into the XylR RNA transcript. RNA is the representation of the transcript, with P representing the number of plasmids inside a cell. J23100TC represents the activity of our promoter controlling XylR. ta is the transcription rate of the xylR into RNA, based on the transcription rate of E. Coli and the size of XylR, and KRNAdeg represents the degradation rate in the E. Coli.

      The second and third equations represent the inactive and active forms of XylR- or the concentration of XylR that has not bonded to xylene (XylRi) and the concentration of XylR that has bonded to xylene (XylRa). RNA is the substitution of our first equation. tr is the translation rate of our RNA into XylR. rXylR is the oligomerization constant of XylR, rR,XylR is the dissociation constant of active XylR, and xyl in the total concentration of xylene. αXylRi accounts for degradation of XylR.

      The fourth equation represents the concentration of xylene present. XylRi and XylRa refer to the concentrations of the inactive and active forms of XylR. rXylR is the oligomerization constant of XylR, rR,XylR is the dissociation constant of inactive XylR and xyl is the total concentration of xylene. αXylRi accounts for the degradation of XylRa, which would release the xylene it held. We multiplied by seconds per hour to get a domain more realistic to our needs.

      The fifth equation represents the output of the pu promoter. PuTC represents the activity of the Pu promoter or the output of protein. αpu is the deactivation rate of the Pu promoter. KXylRa is the activation coefficients of the Pu due to binding and nps,a is the hill coefficient for the interaction. β0 and βPS are the basal and maximal expression, respectively, of the Pu promoter.

      In order to create an output, we used the program Scilab. By inputting these equations into code we were able to have the program calculate our estimates and graph the results. In the end our output was around 5.5 mPoPS. This is a logical output and tells us our biobrick will successfully over express the protein. We set the output of our model in PoPS due to our system's flexibility when it comes to our reporter protein. This allows us to quickly add in an extra equation to turn our current output into fluorescence, LacZ output, or pigment production without having to worry about editing other parts of our model. Below is a copy of our code that we had used in our program. The resulting graph is the output of our Pu promoter.

     While we have yet to apply our modelling framework to the optimization of our sensor, we have created the foundation of a model that can be easily tweaked to test the effect of different variables on sensor performance. Through this we aim to do a sensitivity analysis of the various factors involved in our model in the future so that we can predict what to change for version 2 of our sensor. From there we can return to our model, feeding new data into it to create many future versions of our system, improving performance each time.

Prototyping...

      ECOS was meant to be more than just a theory; more than just a project to play out at school- we wanted ECOS to have real world applications by the time we finished making our construct. We had one major problem: How were we going to test our product in an immediate and effective manner? To answer this we designed two prototypes to design and test. Unfortunately we were unable to finish the wet-lab portion of this project, but the designing is complete thanks to the Maya program via Autodesk. These prototypes are beneficial because they 1) are cheap to build and 2) relatively quick to test. Here are how they turned out:

 

      PROTOTYPE #1: We designed our first prototype so that it has a positive pressure source which creates a current for our xylene to travel from our sample, which is in a heated container so that the xylene is vapourized, to our ECOS container in which we bubble it through the E. Coli culture to get results. This was a plausible idea, as it would require fairly non-expensive materials and could be maintained by the average business person. However, we thought it could be more efficient in it's design, and so we kept looking for ideas. Our biggest problem with this was the fact that it is a relatively complex prototype. So we challenged ourselves to try and create a simpler prototype.

 

      PROTOTYPE #2: Another solution we theorized was to have our ECOS impregnated into alginate and then made into small spheres so that we can simply add them into a heated sample to get our results. In a container, we combine our soil sample and our alginate beads, mix them around, and then sift out the beads. The beads produce our indicator protein from the contamination levels of the soil. This design has been inspired by the Peking 2013 Team who used alginate encapsulation beads. Their reasoning and justification for using alginate beads stood out to us, as they stated that it was stable and inexpensive and easy to shape and manipulate. Upon further research into the use of these beads, we discovered that this prototype could be more efficient to use than our first prototype that required a few different pieces to the system. This revised prototype would also be easier for the average oil worker to use compared to the larger prototype.

 

     As time progresses with our project we plan to test and find what concentrations of m-Xylene trigger what kind of report. To do this we have purchased a couple containers of Xylene, each with a specified concentration, which we will apply to our prototypes. In-lab trials, in combination with our mathemantical model, will allow us to predict what [Xylene] is around and therefore what level of Benzene and Benzene derivatives are around as well. Once we manage to find this correlation between reporter and [Xylene] then we will be moving the project out to attempt some real-world application. Depending on our results here, ECOS may move away from rural Consort Alberta and out accross Canada and potentially the world. Of course there will always be things to fix- expecially when it comes to synthetic biology, so just like a computer programer we would try to fine tune and bring out newer version of this construct. The changes, of course, would come from our observations in-lab and customer suggestions as time progressed.

    References:

  1. The regulatory logic of m-xylene biodegradation by Pseudomonas putida mt-2 exposed by dynamic modelling of the principal node Ps/Pr of the TOL plasmid
  2. Calgary 2012 ECHEM Values, 
  3. Wiley Value References