Team:Penn State/Modelling
From 2011.igem.org
(13 intermediate revisions not shown) | |||
Line 39: | Line 39: | ||
- | p, * {font-family: times, Times New Roman, times-roman, georgia, serif; | + | p, li, * {font-family: times, Times New Roman, times-roman, georgia, serif; |
- | color: # | + | color: #333; |
- | font-size: | + | font-size: 14px;} |
- | + | ||
a:link, a:visited, a:hover {font-family: "Lucida Grande", Tahoma; | a:link, a:visited, a:hover {font-family: "Lucida Grande", Tahoma; | ||
font-size: 12px; | font-size: 12px; | ||
Line 91: | Line 90: | ||
Initially, to model the OR operating region a statistical thermodynamic model developed | Initially, to model the OR operating region a statistical thermodynamic model developed | ||
by Gary Ackers, Alexander Johnosn, and Madeline Shea was used. This model appeared | by Gary Ackers, Alexander Johnosn, and Madeline Shea was used. This model appeared | ||
- | in a 1981 paper entitled Quantitative | + | in a 1981 paper entitled <u>Quantitative Model For Gene Regulation by Lambda Phage |
- | + | Repressor</u>. The model consisted of the following equation, which related the probability | |
- | of the OR operator region in configuration s | + | of the OR operator region in configuration 's' to the Gibbs free energy of binding of |
- | repressor to that site and the concentration of unbound repressor dimer | + | repressor to that site and the concentration of unbound repressor dimer. |
</br> | </br> | ||
<center><img src="https://static.igem.org/mediawiki/2011/3/36/PSUequ1.JPG"></center> | <center><img src="https://static.igem.org/mediawiki/2011/3/36/PSUequ1.JPG"></center> | ||
Line 108: | Line 107: | ||
graph below. | graph below. | ||
</br> | </br> | ||
- | <center><img src="https://static.igem.org/mediawiki/2011/ | + | <center><img src="https://static.igem.org/mediawiki/2011/d/db/Fraction_repressed_plot.png"></center> |
</br> | </br> | ||
This model however was too simplistic in that it only took in to account repressor | This model however was too simplistic in that it only took in to account repressor | ||
Line 117: | Line 116: | ||
</br> | </br> | ||
This new model, also developed by Madeline Shea and Gary Ackers appeared in | This new model, also developed by Madeline Shea and Gary Ackers appeared in | ||
- | the paper The OR Control System of Bacteriophage Lambda A Physical-Chemical Model | + | the paper <u>The OR Control System of Bacteriophage Lambda A Physical-Chemical Model |
- | for Gene Regulation published in 1983. This model had 40 configurations of operator | + | for Gene Regulation</u> published in 1983. This model had 40 configurations of operator |
region in it while the previous paper only had 8, which makes for a more accurate | region in it while the previous paper only had 8, which makes for a more accurate | ||
representation of real life. Once again it relied on a statistical thermodynamic model of | representation of real life. Once again it relied on a statistical thermodynamic model of | ||
Line 176: | Line 175: | ||
repressor approaches zero due to the inhibitory effects of the CRO protein on OR3. | repressor approaches zero due to the inhibitory effects of the CRO protein on OR3. | ||
</p> | </p> | ||
+ | <hr> | ||
+ | <h3> | ||
+ | Modeling Done Since the Regional Jamboree | ||
+ | </h3> | ||
+ | A paper found on the effect of ionizing radiation on D. Radiodurans called <u>A Model of Interactions Between Radiation-induced Oxidative Stress, Protein and DNA Damage in Deinococcus Radiodurans</u> by Igor Shuryak and David J. Brenner was used to model the relationship between double strand breaks and intensity of ionizing radiation. A value of 200 Gy was chosen for the radiation intensity as that seems to be a slightly less than lethal dose for E. Coli. Modifications were made to the parameters in the model due to the increased sensitivity to ionizing radiation of Escherichia Coli. | ||
+ | </br> | ||
+ | <center><img src="https://static.igem.org/mediawiki/2011/7/7d/DSBs_over_time.png" width = 700></center> | ||
+ | </br> | ||
+ | <p> | ||
+ | At a constant dose rate the equilibrium number of double strand breaks can be calculated. This can then be applied to the number of single strand breaks that occur and then to the rate of activation of RecA, which is yet to be done. | ||
+ | </p> | ||
+ | </br> | ||
+ | <center><img src="https://static.igem.org/mediawiki/2011/6/6e/DSBeq.png"></center> | ||
+ | </br> | ||
+ | <h2> | ||
+ | Modeling yet to be done | ||
+ | </h2> | ||
+ | <p> | ||
+ | The overall goal of modeling the system would comprise of the following: a fully developed kinetic model of the system starting from initiation of radiative exposure to the cleavage of catechol or x-gal and the subsequent expression of color. With a fully developed model and known data from testing, future devices could be made with varying sensitivities, which could be accomplished through varying the RBS strengths of the Pr promoter. | ||
+ | </p> | ||
+ | <hr> | ||
+ | The following files consist of the Matlab code used to model the OR operator region based off of the statistical thermodynamic model by Shea and Ackers. | ||
+ | <a href="https://2011.igem.org/File:NewIGEMCode.zip"></br>Matlab code for the masses!!</a> | ||
+ | |||
+ | <hr> | ||
+ | <h3> | ||
+ | Bacterial Dosimeter Cost Analysis | ||
+ | </h3> | ||
+ | <p>The purpose of this report is to emphasis the monetary cost per device that it would cost to produce such a device for sale and compare against major competitors already in the market of selling dosimeters.</p> | ||
+ | <p>When buying our bacterial dosimeter, there are two kinds of cost associated with the device: fixed cost and variable costs. The fixed cost occurs when buying the tray because the tray can be used over and over again with a new tray put in after every use. The variable costs are the cost associated with buying the individual trays per use. Additionally, the cost of catechol must be factored into the equation because it is required to show the results of the dosimeter.</p> | ||
+ | </br> | ||
+ | <h4> | ||
+ | <center>Equation 1: Cost on User for Using our Bacterial Dosimeter to Detect Radiation | ||
+ | </h4> | ||
+ | </br> | ||
+ | <center>Cdosimeter = Ccase + ( N*Ctray) + Ccatechol</br></center> | ||
+ | Cdosimeter = Cost per Bacterial Dosimeter </br> | ||
+ | Ctray = Cost per Tray</br> | ||
+ | Ccase = Cost per Case </br> | ||
+ | N = Number of Trays</br> | ||
+ | Ccatechol= cost of catechol per use</br></br> | ||
+ | <h4> | ||
+ | <center>Equation 2: Cost of building each individual Case in the Dosimeter</br> | ||
+ | </h4> | ||
+ | <center>Ccase = V*Cplastic</br></center> | ||
+ | Cplastic= Cost per cubic meter of plastic </br> | ||
+ | V = Volume of Plastic</br> | ||
+ | <h4> | ||
+ | <center>Equation 3: Cost of building each individual tray in the dosimeter.</br></br> | ||
+ | </h4> | ||
+ | <center>Ctray = Cecoli + Cgellan + Cplastic</br></center> | ||
+ | Cecoli = Cost of Growing the Bacteria </br> | ||
+ | Cgellan= Cost of Phytagel™</br> | ||
+ | <p> | ||
+ | Taking the price of the gellan gum (PHYTAGEL™) in bulk, the average cost would be that of $1,003.20 USD for 5kg of powder. At this price, it would cost approximately Considering that the dimensions for one sheet of gellan gum in our model has the dimensions of 0.035m x 0.02 m x 0.00075 m it has a volume of 5.25e-7 m3. At this volume it would cost approximately $0.20 to produce the gellan gum per dosimeter. | ||
+ | </p> | ||
+ | </br> | ||
+ | <h4> | ||
+ | <center>Equation 4: Final Cost Analysis of Building a Bacterial Dosimeter</br> | ||
+ | </h4> | ||
+ | <center>Cdosimeter = (V*Cplastic) + [n*( Cecoli + Cgellan + Cplastic ) + Ccatechol </br> | ||
+ | </center> | ||
+ | |||
<a href="https://2011.igem.org/Team:Penn_State/sitemap"> <img src="https://static.igem.org/mediawiki/2011/5/51/Sitemap.jpg" /> </a> | <a href="https://2011.igem.org/Team:Penn_State/sitemap"> <img src="https://static.igem.org/mediawiki/2011/5/51/Sitemap.jpg" /> </a> | ||
+ | |||
</div> | </div> | ||
</html> | </html> |
Latest revision as of 03:19, 29 October 2011
Initially, to model the OR operating region a statistical thermodynamic model developed by Gary Ackers, Alexander Johnosn, and Madeline Shea was used. This model appeared in a 1981 paper entitled Quantitative Model For Gene Regulation by Lambda Phage Repressor. The model consisted of the following equation, which related the probability of the OR operator region in configuration 's' to the Gibbs free energy of binding of repressor to that site and the concentration of unbound repressor dimer.
![](https://static.igem.org/mediawiki/2011/9/98/Psumodelgraph1.png)
![](https://static.igem.org/mediawiki/2011/d/db/Fraction_repressed_plot.png)
![](https://static.igem.org/mediawiki/2011/c/c2/PSUEqu2.png)
![](https://static.igem.org/mediawiki/2011/f/fc/PSUmodelgraph3.png)
![](https://static.igem.org/mediawiki/2011/8/8f/PSUmodelgraph4.png)
![](https://static.igem.org/mediawiki/2011/7/7e/PSUequ3.png)
![](https://static.igem.org/mediawiki/2011/d/db/PSUequ4.png)
![](https://static.igem.org/mediawiki/2011/e/e9/PSUmodelgraph5.png)
Modeling Done Since the Regional Jamboree
A paper found on the effect of ionizing radiation on D. Radiodurans called A Model of Interactions Between Radiation-induced Oxidative Stress, Protein and DNA Damage in Deinococcus Radiodurans by Igor Shuryak and David J. Brenner was used to model the relationship between double strand breaks and intensity of ionizing radiation. A value of 200 Gy was chosen for the radiation intensity as that seems to be a slightly less than lethal dose for E. Coli. Modifications were made to the parameters in the model due to the increased sensitivity to ionizing radiation of Escherichia Coli.![](https://static.igem.org/mediawiki/2011/7/7d/DSBs_over_time.png)
At a constant dose rate the equilibrium number of double strand breaks can be calculated. This can then be applied to the number of single strand breaks that occur and then to the rate of activation of RecA, which is yet to be done.
![](https://static.igem.org/mediawiki/2011/6/6e/DSBeq.png)
Modeling yet to be done
The overall goal of modeling the system would comprise of the following: a fully developed kinetic model of the system starting from initiation of radiative exposure to the cleavage of catechol or x-gal and the subsequent expression of color. With a fully developed model and known data from testing, future devices could be made with varying sensitivities, which could be accomplished through varying the RBS strengths of the Pr promoter.
The following files consist of the Matlab code used to model the OR operator region based off of the statistical thermodynamic model by Shea and Ackers. Matlab code for the masses!!
Bacterial Dosimeter Cost Analysis
The purpose of this report is to emphasis the monetary cost per device that it would cost to produce such a device for sale and compare against major competitors already in the market of selling dosimeters.
When buying our bacterial dosimeter, there are two kinds of cost associated with the device: fixed cost and variable costs. The fixed cost occurs when buying the tray because the tray can be used over and over again with a new tray put in after every use. The variable costs are the cost associated with buying the individual trays per use. Additionally, the cost of catechol must be factored into the equation because it is required to show the results of the dosimeter.
Equation 1: Cost on User for Using our Bacterial Dosimeter to Detect Radiation
Equation 2: Cost of building each individual Case in the Dosimeter
Equation 3: Cost of building each individual tray in the dosimeter.
Taking the price of the gellan gum (PHYTAGEL™) in bulk, the average cost would be that of $1,003.20 USD for 5kg of powder. At this price, it would cost approximately Considering that the dimensions for one sheet of gellan gum in our model has the dimensions of 0.035m x 0.02 m x 0.00075 m it has a volume of 5.25e-7 m3. At this volume it would cost approximately $0.20 to produce the gellan gum per dosimeter.
Equation 4: Final Cost Analysis of Building a Bacterial Dosimeter
![](https://static.igem.org/mediawiki/2011/5/51/Sitemap.jpg)