Team:Aberdeen Scotland/Probability

From 2010.igem.org

(Difference between revisions)
 
(24 intermediate revisions not shown)
Line 8: Line 8:
<br>
<br>
<center>
<center>
-
<img src="https://static.igem.org/mediawiki/2010/8/8c/Equations.png">
+
<img src="https://static.igem.org/mediawiki/2010/e/e4/Equations3.png">
</center>
</center>
<br>
<br>
<p>We wrote a program which takes the equations for our system, calculates the nullclines and records whether or not bistability is achieved. Whether or not bistability is achieved depends strongly on the parameters  λ<sub style="font-size:10px">i</sub>, K<sub style="font-size:10px">i</sub>, μ<sub style="font-size:10px">i</sub> and n<sub style="font-size:10px">i</sub>.</p>  
<p>We wrote a program which takes the equations for our system, calculates the nullclines and records whether or not bistability is achieved. Whether or not bistability is achieved depends strongly on the parameters  λ<sub style="font-size:10px">i</sub>, K<sub style="font-size:10px">i</sub>, μ<sub style="font-size:10px">i</sub> and n<sub style="font-size:10px">i</sub>.</p>  
<br>
<br>
-
<p>The values for the parameters λ<sub style="font-size:10px">i</sub>, K<sub style="font-size:10px">i</sub> and μ<sub style="font-size:10px">i</sub>  were taken from literature<a href="#ref1"><sup style="font-size:10px">[1]-[13]</sup></a> but have a large uncertainty attached to them. To take this uncertainty into account we selected each parameter value randomly from a range of two orders of magnitude around the literature value for each parameter. We knew that n<sub style="font-size:10px">1</sub> was between 1 and 3 and n<sub style="font-size:10px">2</sub>=1. This is due to the number of stem loops on the two DNA strands (see Fig 1 in <a href="https://2010.igem.org/Team:Aberdeen_Scotland/Equations">Equations</a>). The strand associated with GFP has two stem loops. A protein can bind to one stem loop and this can then hinder the binding of a protein to the second stem loop (thus making n<sub style="font-size:10px">1</sub> less than 2 but more than 1) or encourage the binding of a protein to the second stem loop (thus making n<sub style="font-size:10px">1</sub> more than 2 but less than 3). The strand associated with the CFP has only one stem loop and thus n<sub style="font-size:10px">2</sub>=1. Please see <a href="https://2010.igem.org/Team:Aberdeen_Scotland/Curve_Fitting">Curve Fitting</a> for an explanation on how we narrowed down the value for n<sub style="font-size:10px">1</sub>.
+
<p>The values for the parameters λ<sub style="font-size:10px">i</sub>, K<sub style="font-size:10px">i</sub> and μ<sub style="font-size:10px">i</sub>  were taken from literature<a href="#ref1"><sup style="font-size:10px">[1]-[13]</sup></a> but have a large uncertainty attached to them. To take this uncertainty into account we selected each parameter value randomly from a range of two orders of magnitude around the literature value for each parameter. We knew that n<sub style="font-size:10px">2</sub> was between 1 and 3 and n<sub style="font-size:10px">4</sub>=1. This is due to the number of stem loops on the two DNA strands (see Fig 1 in <a href="https://2010.igem.org/Team:Aberdeen_Scotland/Equations">Equations</a>). The strand associated with GFP has two stem loops. A protein can bind to one stem loop and this can then hinder the binding of a protein to the second stem loop (thus making n<sub style="font-size:10px">2</sub> less than 2 but more than 1) or encourage the binding of a protein to the second stem loop (thus making n<sub style="font-size:10px">2</sub> more than 2 but less than 3). The strand associated with the CFP has only one stem loop and thus n<sub style="font-size:10px">4</sub>=1. Please see <a href="https://2010.igem.org/Team:Aberdeen_Scotland/Curve_Fitting">Determination of the Hill coefficient n<sub style="font-size:10px">2</sub></a> for an explanation on how we narrowed down the value for n<sub style="font-size:10px">2</sub>.
<br>
<br>
<br>
<br>
Line 25: Line 25:
<br>
<br>
We then decided to run the program for various different parameter ranges. We cannot change individual parameters specifically, but in the lab it is possible to vary the order of magnitude of a parameter. We wanted to know, if we varied our parameters and Hill coefficients, what is the best percentage we can possibly get? These results are shown in scenarios 2-4.</p>
We then decided to run the program for various different parameter ranges. We cannot change individual parameters specifically, but in the lab it is possible to vary the order of magnitude of a parameter. We wanted to know, if we varied our parameters and Hill coefficients, what is the best percentage we can possibly get? These results are shown in scenarios 2-4.</p>
-
 
+
<br>
<h3>Scenario 1 – unmodified parameters </h3>
<h3>Scenario 1 – unmodified parameters </h3>
<br>
<br>
Line 35: Line 35:
<br>
<br>
<br>
<br>
-
We did not know what concentration of galactose and methionine was optimal, so we ran various tests and found the value for which the percentages were maximised. We found that the optimal percentages occurred if we have 100000 times more galactose than methionine. See figure 1.</p>
+
We did not know what concentration of galactose and methionine was optimal, so we ran various tests and found the value for which the percentages were maximised. We found that the optimal percentages occurred if we have 100000 times more galactose than methionine. See figure 1a and 1b.</p>
<br>
<br>
<br>
<br>
<center>
<center>
-
<img src="https://static.igem.org/mediawiki/2010/9/9f/Original_parameters_colorplot.png">
+
a)<img src="https://static.igem.org/mediawiki/2010/d/d2/Table1.png">
-
<p style="font-size:10px">Figure 1. The percentage of times our system will work for various combinations of Hill coefficients for the original parameters. n<sub style="font-size:5px">1</sub> is the Hill coefficient for the GFP rate-reaction equation and n<sub style="font-size:5px">2</sub> is the Hill coefficient for the CFP rate-reaction equation.</p>
+
<br>
 +
<br>
 +
b)<img src="https://static.igem.org/mediawiki/2010/9/9f/Original_parameters_colorplot.png">
 +
<p style="font-size:10px">Figure 1. a) The percentage of parameter combinations tested that will give us bistability for various combinations of Hill coefficients. n<sub style="font-size:5px">2</sub> is the Hill coefficient for the GFP rate-reaction equation and n<sub style="font-size:5px">4</sub> is the Hill coefficient for the CFP rate-reaction equation. b) Graphical representation of 1a. </p>
</center>
</center>
<br>
<br>
-
<p>Here we can see that if n<sub style="font-size:10px">2</sub>=1, then the best scenario we can hope for is that n<sub style="font-size:10px">1</sub>=5. In this situation, 3.29% of the parameter combinations tested will give bistability and the possibility of switching. However, with n<sub style="font-size:10px">1</sub> estimated to be between 1 and 3, we can say that at most 2.03% of the parameter combinations tested gave bistability. </p>  
+
<br>
-
 
+
<p>Here we can see that if n<sub style="font-size:10px">4</sub>=1, then the best scenario we can hope for is that n<sub style="font-size:10px">2</sub>=5. In this situation, 3.29% of the parameter combinations tested will give bistability and the possibility of switching. However, with n<sub style="font-size:10px">2</sub> estimated to be between 1 and 3, we can say that at most 2.03% of the parameter combinations tested gave bistability. </p>  
 +
<br>
<h3>Scenario 2 – modified parameters 1</h3>
<h3>Scenario 2 – modified parameters 1</h3>
<br>
<br>
Line 52: Line 56:
</center>
</center>
<br>
<br>
-
<p>Percentages are maximised if we have 100000 times more galactose than methionine. See figure 2.</p>
+
<br>
 +
<p>Percentages are maximised if we have 100000 times more galactose than methionine. See figure 2a and 2b.</p>
 +
<br>
<br>
<br>
<center>
<center>
-
<img src="https://static.igem.org/mediawiki/2010/d/dd/Table2b.png">
+
a) <img src="https://static.igem.org/mediawiki/2010/d/dd/Table2b.png">
-
<p style="font-size:10px">Figure 2. The percentage of times our system will work for various combinations of Hill coefficients for the modified parameters.</p>
+
<br>
 +
<br>
 +
b) <img src="https://static.igem.org/mediawiki/2010/3/3b/Mod1.png">
 +
<p style="font-size:10px">Figure 2. a) The percentage of parameter combinations tested that will give us bistability for various combinations of Hill coefficients. b) Graphical representation of 2a.</p>
</center>
</center>
<br>
<br>
-
<p>The best scenario is if we have n<sub style="font-size:10px">1</sub>=5 and n<sub style="font-size:10px">2</sub>=5.</p>
+
<br>
-
 
+
<p>The best scenario is if we have n<sub style="font-size:10px">2</sub>=5 and n<sub style="font-size:10px">4</sub>=5.</p>
 +
<br>
<h3>Scenario 3 – modified parameters 2</h3>
<h3>Scenario 3 – modified parameters 2</h3>
<br>
<br>
Line 67: Line 77:
</center>
</center>
<br>
<br>
-
<p>Percentages are maximised if there is 100000 times more galactose than methionine. See table 3.</p>
+
<br>
 +
<p>Percentages are maximised if there is 100000 times more galactose than methionine. See figure 3a and 3b.</p>
 +
<br>
<br>
<br>
<center>
<center>
-
<img src="https://static.igem.org/mediawiki/2010/2/29/Table3b.png">
+
a) <img src="https://static.igem.org/mediawiki/2010/2/29/Table3b.png">
-
<p style="font-size:10px">Table 3. The percentage of times our system will work for various combinations of Hill coefficients for the modified parameters.</p>
+
<br>
 +
<br>
 +
b) <img src="https://static.igem.org/mediawiki/2010/d/da/Mod2.png">
 +
<p style="font-size:10px">Figure 3. a) The percentage of parameter combinations tested that will give us bistability for various combinations of Hill coefficients. b) Graphical representation of 3a.</p>
</center>
</center>
<br>
<br>
-
<p>Here we find that the optimal scenario is n<sub style="font-size:10px">1</sub>=3 and n<sub style="font-size:10px">2</sub>=5. The trend which has developed in this table will be discussed at the end of scenario 4. </p>
+
<br>
-
 
+
<p>Here we find that the optimal scenario is n<sub style="font-size:10px">2</sub>=3 and n<sub style="font-size:10px">4</sub>=5. The trend which has developed in this table will be discussed at the end of scenario 4. </p>
 +
<br>
<h3>Scenario 4 – modified parameters 3 </h3>
<h3>Scenario 4 – modified parameters 3 </h3>
<br>
<br>
Line 82: Line 98:
</center>
</center>
<br>
<br>
-
<p>Percentages are maximised if we have 10000000 times more galactose than methionine. See table 4.</p>
+
<br>
 +
<p>Percentages are maximised if we have 10000000 times more galactose than methionine. See Figure 4a and 4b.</p>
 +
<br>
<br>
<br>
<center>
<center>
-
<img src="https://static.igem.org/mediawiki/2010/2/22/Table4b.png">
+
a) <img src="https://static.igem.org/mediawiki/2010/2/22/Table4b.png">
-
<p style="font-size:10px">Table 4. The percentage of times our system will work for various combinations of Hill coefficients for the modified parameters.</p>
+
<br>
 +
<br>
 +
b) <img src="https://static.igem.org/mediawiki/2010/1/1b/Modified_parameters_colorplot.png">
 +
<p style="font-size:10px">Figure 4. a) The percentage of parameter combinations tested that will give us bistability for various combinations of Hill coefficients. b) Graphical representation of 4a.</p>
</center>
</center>
<br>
<br>
-
<p>An interesting trend is seen in scenarios 3 and 4. The percentages increase to a maximum at n<sub style="font-size:10px">1</sub>=2 and n<sub style="font-size:10px">2</sub>=5 (n<sub style="font-size:10px">1</sub>=3 and n<sub style="font-size:10px">2</sub>=5 in scenario 3) and then decrease, quite rapidly, towards Hill coefficients of n<sub style="font-size:10px">1</sub>=5 and n<sub style="font-size:10px">2</sub>=5. In the parameter space, there is a region where bistability is possible. For the smaller Hill coefficients, the randomly chosen parameters (as described above) would allow the nullclines to have a shape such that they would cross three or more times and bistability is achieved. The higher Hill coefficients, however, for the same range of parameters, would have nullclines which were shaped so that they did not cross three or more times – in parameter space, they fall outside the region of bistability. </p>
+
<br>
-
 
+
<p>An interesting trend is seen in scenarios 3 and 4. The percentages increase to a maximum at n<sub style="font-size:10px">2</sub>=2 and n<sub style="font-size:10px">4</sub>=5 (n<sub style="font-size:10px">2</sub>=3 and n<sub style="font-size:10px">4</sub>=5 in scenario 3) and then decrease, quite rapidly, towards Hill coefficients of n<sub style="font-size:10px">2</sub>=5 and n<sub style="font-size:10px">4</sub>=5. In the parameter space, there is a region where bistability is possible. For the smaller Hill coefficients, the randomly chosen parameters (as described above) would allow the nullclines to have a shape such that they would cross three or more times and bistability is achieved. The higher Hill coefficients, however, for the same range of parameters, would have nullclines which were shaped so that they did not cross three or more times – in parameter space, they fall outside the region of bistability. </p>
 +
<br>
<h3>Conclusions</h3>
<h3>Conclusions</h3>
-
<p>The ideal scenario is to change the parameters to that in scenario 4 and have Hill coefficients of n<sub style="font-size:10px">1</sub>=2 and n<sub style="font-size:10px">2</sub>=5. Here 98.17% of the parameter combinations tested gave bistability. Contrary to our initial assumptions, depending on the range of parameters, higher Hill coefficients are not always better. The parameters will dictate the shape and position of the nullclines and whether or not bistability is possible.  
+
<p>The ideal scenario is to change the parameters to that in scenario 4 and have Hill coefficients of n<sub style="font-size:10px">2</sub>=2 and n<sub style="font-size:10px">4</sub>=5. Here 98.17% of the parameter combinations tested gave bistability. Contrary to our initial assumptions, depending on the range of parameters, higher Hill coefficients are not always better. The parameters will dictate the shape and position of the nullclines and whether or not bistability is possible.  
<br>
<br>
<br>
<br>
It is possible to determine the optimal values of the parameters for an ideal working of the switch. However, to change the parameters in the experimental construct so much would require much more time than we had available. </p>
It is possible to determine the optimal values of the parameters for an ideal working of the switch. However, to change the parameters in the experimental construct so much would require much more time than we had available. </p>
-
 
+
<br>
<h3>References</h3>
<h3>References</h3>
<br>
<br>
Line 103: Line 125:
<p><sup style="font-size:10px">[1]</sup> Hershey, J., Sonenberg, N., and Mathews, M. (2007), Origins and Principles of Translational Control. In Translational Control (pp.1-29). Cold Spring Harbor, NY: Cold Spring Harbor Press. </p>
<p><sup style="font-size:10px">[1]</sup> Hershey, J., Sonenberg, N., and Mathews, M. (2007), Origins and Principles of Translational Control. In Translational Control (pp.1-29). Cold Spring Harbor, NY: Cold Spring Harbor Press. </p>
<br>
<br>
-
<p><sup style="font-size:10px">[2]</sup> Edwards, A.M., Kane, C.M., Young, R.A., and Kornberg, R.D. (1991), Two Dissociable Subunits of Yeast RNA Polymerase II Stimulate the Initiation of Transcription at a Promoter in Vitro, J. Biol. Chem., Vol. 266, No.1, pp.71-75 </p>
+
<p><sup style="font-size:10px">[2]</sup> Edwards, A.M., Kane, C.M., Young, R.A., and Kornberg, R.D. (1991), Two Dissociable Subunits of Yeast RNA Polymerase II Stimulate the Initiation of Transcription at a Promoter in Vitro, <i>J. Biol. Chem.</i>, Vol. 266, No.1, pp.71-75 </p>
<br>
<br>
-
<p><sup style="font-size:10px">[3]</sup> Fürst, P., and Hamer, D. (1989), Cooperative activation of a eukaryotic transcription factor: Interaction between Cu(I) and yeast ACE1 protein, Proc. Natl. Acad. Sci. USA, Vol. 86, pp. 5267-5271 </p>
+
<p><sup style="font-size:10px">[3]</sup> Fürst, P., and Hamer, D. (1989), Cooperative activation of a eukaryotic transcription factor: Interaction between Cu(I) and yeast ACE1 protein, <i>Proc. Natl. Acad. Sci. USA</i>, Vol. 86, pp. 5267-5271 </p>
<br>
<br>
-
<p><sup style="font-size:10px">[4]</sup> Jensen, L.T., Posewitz, M.C., Srinivasan, C., and Winge, D.R. (1998), Mapping of the DNA binding domain of the copper-responsive transcription factor Mac1 from Saccharomyces cerevisiae, J. Biol. Chem., Vol. 273, No. 37, pp.23805-23811 </p>
+
<p><sup style="font-size:10px">[4]</sup> Jensen, L.T., Posewitz, M.C., Srinivasan, C., and Winge, D.R. (1998), Mapping of the DNA binding domain of the copper-responsive transcription factor Mac1 from Saccharomyces cerevisiae, <i>J. Biol. Chem.</i>, Vol. 273, No. 37, pp.23805-23811 </p>
<br>
<br>
-
<p><sup style="font-size:10px">[5]</sup> Rodgers, K.K., and Coleman, J.E. (1994), DNA binding and bending by the transcription factors GAL4(62*) and GAL4(149*), Protein Science, Vol.3, No. 4, pp. 608-619 </p>
+
<p><sup style="font-size:10px">[5]</sup> Rodgers, K.K., and Coleman, J.E. (1994), DNA binding and bending by the transcription factors GAL4(62*) and GAL4(149*), <i>Protein Science</i>, Vol.3, No. 4, pp. 608-619 </p>
<br>
<br>
-
<p><sup style="font-size:10px">[6]</sup> Garcia-Martinez, J., Aranda, A., and Perez-Ortin, J.E. (2004), Genomic run-on evaluates transcription rates for all yeast genes and identifies gene regulatory mechanisms, Mol. Cell., Vol. 15, No. 2, pp.303-313 </p>
+
<p><sup style="font-size:10px">[6]</sup> Garcia-Martinez, J., Aranda, A., and Perez-Ortin, J.E. (2004), Genomic run-on evaluates transcription rates for all yeast genes and identifies gene regulatory mechanisms, <i>Mol. Cell.</i>, Vol. 15, No. 2, pp.303-313 </p>
<br>
<br>
-
<p><sup style="font-size:10px">[7]</sup> Wang, Y., Liu, C.L., Storey, J.D., Tibshirani, R.J., Herschlag, D., and Brown, P.O. (2002), Precision and functional specificity in mRNA decay, Proc. Natl. Acad. Sci. USA, Vol. 99, No.9, pp.5860-5865 </p>
+
<p><sup style="font-size:10px">[7]</sup> Wang, Y., Liu, C.L., Storey, J.D., Tibshirani, R.J., Herschlag, D., and Brown, P.O. (2002), Precision and functional specificity in mRNA decay, <i>Proc. Natl. Acad. Sci. USA</i>, Vol. 99, No.9, pp.5860-5865 </p>
<br>
<br>
-
<p><sup style="font-size:10px">[8]</sup> Beyer, A. et al. (2004), Post-transcriptional sxpression regulation in the Yeast Saccharomyces Cerevisiae on a genomic scale, Molecular and Cellular Proteomics, Vol. 3, No. 11, pp.1083-1092 </p>
+
<p><sup style="font-size:10px">[8]</sup> Beyer, A. et al. (2004), Post-transcriptional sxpression regulation in the Yeast Saccharomyces Cerevisiae on a genomic scale, <i>Molecular and Cellular Proteomics</i>, Vol. 3, No. 11, pp.1083-1092 </p>
<br>
<br>
-
<p><sup style="font-size:10px">[9]</sup> Mateus, C. and Avery, S.V. (2000), Destabilized green fluorescent protein for monitoring dynamic changes in yeast gene expression with flow cytometry, Yeast, Vol.16, No.14, pp.1313-1323 </p>
+
<p><sup style="font-size:10px">[9]</sup> Mateus, C. and Avery, S.V. (2000), Destabilized green fluorescent protein for monitoring dynamic changes in yeast gene expression with flow cytometry, <i>Yeast</i>, Vol.16, No.14, pp.1313-1323 </p>
<br>
<br>
-
<p><sup style="font-size:10px">[10]</sup> Carey, J., Cameron, V., de Haseth, P.L., and Uhlenbeck, O.C. (1983), Sequence-specific interaction of R17 coat protein with its ribonucleic acid binding site, Biochemistry, Vol.22, No.11, pp.2601-2610 </p>
+
<p><sup style="font-size:10px">[10]</sup> Carey, J., Cameron, V., de Haseth, P.L., and Uhlenbeck, O.C. (1983), Sequence-specific interaction of R17 coat protein with its ribonucleic acid binding site, <i>Biochemistry</i>, Vol.22, No.11, pp.2601-2610 </p>
<br>
<br>
-
<p><sup style="font-size:10px">[11]</sup> Talbot, S.J., Goodman, S., Bates, S.R., Fishwick, C.W., and Stockley, P.G. (1990), Use of synthetic oligoribonucleotides to probe RNA-protein interactions in the MS2 translational operator complex, Nucleic Acids Res., Vol.18, No.12, pp.3521-3528 </p>
+
<p><sup style="font-size:10px">[11]</sup> Talbot, S.J., Goodman, S., Bates, S.R., Fishwick, C.W., and Stockley, P.G. (1990), Use of synthetic oligoribonucleotides to probe RNA-protein interactions in the MS2 translational operator complex, <i>Nucleic Acids Res.</i>, Vol.18, No.12, pp.3521-3528 </p>
<br>
<br>
-
<p><sup style="font-size:10px">[12]</sup> Arava, Y., Wang, Y., Storey, J.D., Liu, C.L., Brown, P.O., and Herschlag, D. (2003), Genome-wide analysis of mRNA translation profiles in Saccharomyces cerevisiae, Proc. Natl. Acad. Sci. USA, Vol. 100, No.7, pp.3889-3894 </p>
+
<p><sup style="font-size:10px">[12]</sup> Arava, Y., Wang, Y., Storey, J.D., Liu, C.L., Brown, P.O., and Herschlag, D. (2003), Genome-wide analysis of mRNA translation profiles in Saccharomyces cerevisiae, <i>Proc. Natl. Acad. Sci. USA</i>, Vol. 100, No.7, pp.3889-3894 </p>
<br>
<br>
-
<p><sup style="font-size:10px">[13]</sup> Chattopadhyay, S., Garcia-Mena, J., DeVito, J., Wolska, K., and Das, A. (1995), Bipartite function of a small RNA hairpin in transcription antitermination in bacteriophage lambda, Proc. Natl. Acad. Sci. USA, Vol.92, No.9, pp.4061-4065 </p>
+
<p><sup style="font-size:10px">[13]</sup> Chattopadhyay, S., Garcia-Mena, J., DeVito, J., Wolska, K., and Das, A. (1995), Bipartite function of a small RNA hairpin in transcription antitermination in bacteriophage lambda, <i>Proc. Natl. Acad. Sci. USA</i>, Vol.92, No.9, pp.4061-4065 </p>
<br><br>
<br><br>
Line 132: Line 154:
<tr>
<tr>
<td>
<td>
-
<a href="https://2010.igem.org/Team:Aberdeen_Scotland/Stability"><img src="https://static.igem.org/mediawiki/2010/8/8e/Left_arrow.png">&nbsp;&nbsp;Return to Stability</a>
+
<a href="https://2010.igem.org/Team:Aberdeen_Scotland/Stochastic_Model"><img src="https://static.igem.org/mediawiki/2010/8/8e/Left_arrow.png">&nbsp;&nbsp;Return to the Stochastic Model</a>
</td>
</td>
<td align="right">
<td align="right">
-
<a href="https://2010.igem.org/Team:Aberdeen_Scotland/Curve_Fitting">Continue to Determination of the the Hill Coefficient n<sub style=font-size:10px">1</sub>&nbsp;&nbsp;<img src="https://static.igem.org/mediawiki/2010/3/36/Right_arrow.png"></a>
+
<a href="https://2010.igem.org/Team:Aberdeen_Scotland/Curve_Fitting">Continue to Determination of the Hill Coefficient n<sub style=font-size:10px">2</sub>&nbsp;&nbsp;<img src="https://static.igem.org/mediawiki/2010/3/36/Right_arrow.png"></a>
</td>
</td>
</tr>
</tr>
Line 142: Line 164:
</html>
</html>
 +
{{:Team:Aberdeen_Scotland/Footer}}

Latest revision as of 20:36, 27 October 2010

University of Aberdeen - ayeSwitch - iGEM 2010

Parameter Space Analysis

Recall that the equations of our system are as follows (see Equations page for parameter definitions):



We wrote a program which takes the equations for our system, calculates the nullclines and records whether or not bistability is achieved. Whether or not bistability is achieved depends strongly on the parameters λi, Ki, μi and ni.


The values for the parameters λi, Ki and μi were taken from literature[1]-[13] but have a large uncertainty attached to them. To take this uncertainty into account we selected each parameter value randomly from a range of two orders of magnitude around the literature value for each parameter. We knew that n2 was between 1 and 3 and n4=1. This is due to the number of stem loops on the two DNA strands (see Fig 1 in Equations). The strand associated with GFP has two stem loops. A protein can bind to one stem loop and this can then hinder the binding of a protein to the second stem loop (thus making n2 less than 2 but more than 1) or encourage the binding of a protein to the second stem loop (thus making n2 more than 2 but less than 3). The strand associated with the CFP has only one stem loop and thus n4=1. Please see Determination of the Hill coefficient n2 for an explanation on how we narrowed down the value for n2.

We ran the program for various combinations of Hill coefficients (ni) between 1 and 5. Each time the program was run, the Hill coefficients remained constant but the other parameters varied.


Choosing the parameters in this way meant that each time the program was run we would sometimes achieve a combination of parameters that allowed bistability, and sometimes not.

The program was run 100 times and the number of times bistability was achieved was output to the screen as a percentage.

We then decided to run the program for various different parameter ranges. We cannot change individual parameters specifically, but in the lab it is possible to vary the order of magnitude of a parameter. We wanted to know, if we varied our parameters and Hill coefficients, what is the best percentage we can possibly get? These results are shown in scenarios 2-4.


Scenario 1 – unmodified parameters



These parameters were taken from, or derived from, the information contained within various sources of literature[1]-[13]. The parameters were converted from concentrations to numbers of molecules.

We did not know what concentration of galactose and methionine was optimal, so we ran various tests and found the value for which the percentages were maximised. We found that the optimal percentages occurred if we have 100000 times more galactose than methionine. See figure 1a and 1b.



a)

b)

Figure 1. a) The percentage of parameter combinations tested that will give us bistability for various combinations of Hill coefficients. n2 is the Hill coefficient for the GFP rate-reaction equation and n4 is the Hill coefficient for the CFP rate-reaction equation. b) Graphical representation of 1a.



Here we can see that if n4=1, then the best scenario we can hope for is that n2=5. In this situation, 3.29% of the parameter combinations tested will give bistability and the possibility of switching. However, with n2 estimated to be between 1 and 3, we can say that at most 2.03% of the parameter combinations tested gave bistability.


Scenario 2 – modified parameters 1




Percentages are maximised if we have 100000 times more galactose than methionine. See figure 2a and 2b.



a)

b)

Figure 2. a) The percentage of parameter combinations tested that will give us bistability for various combinations of Hill coefficients. b) Graphical representation of 2a.



The best scenario is if we have n2=5 and n4=5.


Scenario 3 – modified parameters 2




Percentages are maximised if there is 100000 times more galactose than methionine. See figure 3a and 3b.



a)

b)

Figure 3. a) The percentage of parameter combinations tested that will give us bistability for various combinations of Hill coefficients. b) Graphical representation of 3a.



Here we find that the optimal scenario is n2=3 and n4=5. The trend which has developed in this table will be discussed at the end of scenario 4.


Scenario 4 – modified parameters 3




Percentages are maximised if we have 10000000 times more galactose than methionine. See Figure 4a and 4b.



a)

b)

Figure 4. a) The percentage of parameter combinations tested that will give us bistability for various combinations of Hill coefficients. b) Graphical representation of 4a.



An interesting trend is seen in scenarios 3 and 4. The percentages increase to a maximum at n2=2 and n4=5 (n2=3 and n4=5 in scenario 3) and then decrease, quite rapidly, towards Hill coefficients of n2=5 and n4=5. In the parameter space, there is a region where bistability is possible. For the smaller Hill coefficients, the randomly chosen parameters (as described above) would allow the nullclines to have a shape such that they would cross three or more times and bistability is achieved. The higher Hill coefficients, however, for the same range of parameters, would have nullclines which were shaped so that they did not cross three or more times – in parameter space, they fall outside the region of bistability.


Conclusions

The ideal scenario is to change the parameters to that in scenario 4 and have Hill coefficients of n2=2 and n4=5. Here 98.17% of the parameter combinations tested gave bistability. Contrary to our initial assumptions, depending on the range of parameters, higher Hill coefficients are not always better. The parameters will dictate the shape and position of the nullclines and whether or not bistability is possible.

It is possible to determine the optimal values of the parameters for an ideal working of the switch. However, to change the parameters in the experimental construct so much would require much more time than we had available.


References


[1] Hershey, J., Sonenberg, N., and Mathews, M. (2007), Origins and Principles of Translational Control. In Translational Control (pp.1-29). Cold Spring Harbor, NY: Cold Spring Harbor Press.


[2] Edwards, A.M., Kane, C.M., Young, R.A., and Kornberg, R.D. (1991), Two Dissociable Subunits of Yeast RNA Polymerase II Stimulate the Initiation of Transcription at a Promoter in Vitro, J. Biol. Chem., Vol. 266, No.1, pp.71-75


[3] Fürst, P., and Hamer, D. (1989), Cooperative activation of a eukaryotic transcription factor: Interaction between Cu(I) and yeast ACE1 protein, Proc. Natl. Acad. Sci. USA, Vol. 86, pp. 5267-5271


[4] Jensen, L.T., Posewitz, M.C., Srinivasan, C., and Winge, D.R. (1998), Mapping of the DNA binding domain of the copper-responsive transcription factor Mac1 from Saccharomyces cerevisiae, J. Biol. Chem., Vol. 273, No. 37, pp.23805-23811


[5] Rodgers, K.K., and Coleman, J.E. (1994), DNA binding and bending by the transcription factors GAL4(62*) and GAL4(149*), Protein Science, Vol.3, No. 4, pp. 608-619


[6] Garcia-Martinez, J., Aranda, A., and Perez-Ortin, J.E. (2004), Genomic run-on evaluates transcription rates for all yeast genes and identifies gene regulatory mechanisms, Mol. Cell., Vol. 15, No. 2, pp.303-313


[7] Wang, Y., Liu, C.L., Storey, J.D., Tibshirani, R.J., Herschlag, D., and Brown, P.O. (2002), Precision and functional specificity in mRNA decay, Proc. Natl. Acad. Sci. USA, Vol. 99, No.9, pp.5860-5865


[8] Beyer, A. et al. (2004), Post-transcriptional sxpression regulation in the Yeast Saccharomyces Cerevisiae on a genomic scale, Molecular and Cellular Proteomics, Vol. 3, No. 11, pp.1083-1092


[9] Mateus, C. and Avery, S.V. (2000), Destabilized green fluorescent protein for monitoring dynamic changes in yeast gene expression with flow cytometry, Yeast, Vol.16, No.14, pp.1313-1323


[10] Carey, J., Cameron, V., de Haseth, P.L., and Uhlenbeck, O.C. (1983), Sequence-specific interaction of R17 coat protein with its ribonucleic acid binding site, Biochemistry, Vol.22, No.11, pp.2601-2610


[11] Talbot, S.J., Goodman, S., Bates, S.R., Fishwick, C.W., and Stockley, P.G. (1990), Use of synthetic oligoribonucleotides to probe RNA-protein interactions in the MS2 translational operator complex, Nucleic Acids Res., Vol.18, No.12, pp.3521-3528


[12] Arava, Y., Wang, Y., Storey, J.D., Liu, C.L., Brown, P.O., and Herschlag, D. (2003), Genome-wide analysis of mRNA translation profiles in Saccharomyces cerevisiae, Proc. Natl. Acad. Sci. USA, Vol. 100, No.7, pp.3889-3894


[13] Chattopadhyay, S., Garcia-Mena, J., DeVito, J., Wolska, K., and Das, A. (1995), Bipartite function of a small RNA hairpin in transcription antitermination in bacteriophage lambda, Proc. Natl. Acad. Sci. USA, Vol.92, No.9, pp.4061-4065





Back to the Top