Team:Aberdeen Scotland/Evolution
From 2010.igem.org
(19 intermediate revisions not shown) | |||
Line 3: | Line 3: | ||
<html> | <html> | ||
<h1>What is directed evolution?</h1> | <h1>What is directed evolution?</h1> | ||
- | <p>In protein engineering, the method of <a href="http://en.wikipedia.org/wiki/Directed_evolution">directed evolution</a> is used to harness the power of natural selection. | + | <p>In protein engineering, the method of <a href="http://en.wikipedia.org/wiki/Directed_evolution">directed evolution</a> is used to harness the power of natural selection. It is used to evolve certain proteins or RNA with desirable properties not found in nature. Rather than swapping out different components of the structure, directed evolution can be used to manipulate the behaviour of these components. Then, by screening or selection, we can obtain the desired phenotype.</p> |
<br> | <br> | ||
<br> | <br> | ||
Line 19: | Line 19: | ||
<p>In the experiments, we were able to use the FACS (fluorescence activated cells sorting) machine to separate the proteins with the brightest colours and thus use the strongest ones. Once we had this, we could put the proteins through the system again as initial conditions and then we would have a more distinctive switch.</p> | <p>In the experiments, we were able to use the FACS (fluorescence activated cells sorting) machine to separate the proteins with the brightest colours and thus use the strongest ones. Once we had this, we could put the proteins through the system again as initial conditions and then we would have a more distinctive switch.</p> | ||
<br> | <br> | ||
- | + | <div align="center"><img src="https://static.igem.org/mediawiki/2010/0/09/FACS.jpg"/></div> | |
<h3>Pros</h3> | <h3>Pros</h3> | ||
Line 29: | Line 29: | ||
<p>There are three primary limitations that mathematical modelling can overcome:</p> | <p>There are three primary limitations that mathematical modelling can overcome:</p> | ||
<p><b>1.</b> The evolutionary space for a genetic circuit is too large to explore efficiently.</p> | <p><b>1.</b> The evolutionary space for a genetic circuit is too large to explore efficiently.</p> | ||
- | <p><b>2.</b> | + | <p><b>2.</b> Refining our parameters is much easier than improving (or altering) the system experimentally.</p> |
<p><b>3.</b> Although selecting for independent properties is possible, it usually requires setting up multiple rounds of screening or selection.</p> | <p><b>3.</b> Although selecting for independent properties is possible, it usually requires setting up multiple rounds of screening or selection.</p> | ||
<br> | <br> | ||
<h3>Procedure</h3> | <h3>Procedure</h3> | ||
- | <p><b>1.</b> | + | <p><b>1.</b> We first took the basic program where we integrated the system, plugging in an estimation of the parameters. |
<br> | <br> | ||
- | <p><b>2.</b> | + | <p><b>2.</b> For each parameter, it yielded the value for which we obtain the largest difference between GFP and CFP. |
<br> | <br> | ||
- | <p><b>3.</b> | + | <p><b>3.</b> We then fixed our chosen parameter and continue the program with the next parameter. |
<br> | <br> | ||
- | <p><b>4.</b> | + | <p><b>4.</b> This can be done over and over for each parameter of the model, eventually producing a clearly defined switch that will somewhat resemble our original system. |
<br> | <br> | ||
- | <p><b>5.</b> Repeat these steps for each parameter in the system to | + | <p><b>5.</b> Repeat these steps for each parameter in the system to eventually produce a clearly optimized switch that will yield the maximum difference between GFP and CFP when switching from one state to the other. |
<br> | <br> | ||
<br> | <br> | ||
Line 49: | Line 49: | ||
<p>Fig 1. shows the behaviour of the original system over a set period of time.</p> | <p>Fig 1. shows the behaviour of the original system over a set period of time.</p> | ||
<br> | <br> | ||
- | <div align="center"><img src="https://static.igem.org/mediawiki/2010/ | + | <div align="center"><img src="https://static.igem.org/mediawiki/2010/3/38/Initial.jpg"/></div> |
<div align="center"><p><b>Fig 1.</b></p></div> | <div align="center"><p><b>Fig 1.</b></p></div> | ||
<br> | <br> | ||
<p>As shown, the system has been integrated over a long enough timespan such that it has reached equilibrium.</p> | <p>As shown, the system has been integrated over a long enough timespan such that it has reached equilibrium.</p> | ||
+ | <br> | ||
+ | <p>Once we had this, we could then do this again but change one parameter and instead insert a range of values that can be put through the system.</p> | ||
+ | <br> | ||
+ | <p>Then we found the value for that parameter that produced the largest difference between the two proteins. | ||
+ | We did this for each parameter until we had an “ideal” value for each one. Finally, we entered these values into the system and integrated it again.</p> | ||
+ | <br> | ||
+ | <div align="center"><img src="https://static.igem.org/mediawiki/2010/6/6e/Final.jpg"/></div> | ||
+ | <div align="center"><p><b>Fig 2.</b></p></div> | ||
+ | <div align="center"><img src="http://i252.photobucket.com/albums/hh25/lisa21a/galcopper.jpg"/></div> | ||
+ | <div align="center"><p><b>Fig 3.</b></p></div> | ||
+ | <br> | ||
+ | <p>Fig 2. shows our final result of the mRNAs and the two proteins. We can see at first that CFP was the brightest protein initially, then we added more galactose and, at around time = 50, we can clearly see the decline of CFP and the increase of GFP. Fig 3. shows how copper and galactose are acting in the system. Thus, a clear switch in our system. The results for each parameter were as follows;</p> | ||
+ | <br> | ||
+ | <div align="center"> | ||
+ | <table> | ||
+ | <tr> | ||
+ | <td> | ||
+ | <div align="right"><b><p>Parameter</p></b></div> | ||
+ | </td> | ||
+ | <td> | ||
+ | <p><b>Optimal Valued</b></p> | ||
+ | </td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td> | ||
+ | <div align="right"><p>K_1</p></div> | ||
+ | </td> | ||
+ | <td> | ||
+ | <p>5000</p> | ||
+ | </td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td> | ||
+ | <div align="right"><p>K_2</p></div> | ||
+ | </td> | ||
+ | <td> | ||
+ | <p>10000</p> | ||
+ | </td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td> | ||
+ | <div align="right"><p>K_3</p></div> | ||
+ | </td> | ||
+ | <td> | ||
+ | <p>5000</p> | ||
+ | </td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td> | ||
+ | <div align="right"><p>K_4</p></div> | ||
+ | </td> | ||
+ | <td> | ||
+ | <p>100000</p> | ||
+ | </td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td> | ||
+ | <div align="right"><p>lambda_1</p></div> | ||
+ | </td> | ||
+ | <td> | ||
+ | <p>133.9940</p> | ||
+ | </td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td> | ||
+ | <div align="right"><p>lambda_2</p></div> | ||
+ | </td> | ||
+ | <td> | ||
+ | <p>415.9916</p> | ||
+ | </td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td> | ||
+ | <div align="right"><p>lambda_3</p></div> | ||
+ | </td> | ||
+ | <td> | ||
+ | <p>583.9984</p> | ||
+ | </td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td> | ||
+ | <div align="right"><p>lambda_4</p></div> | ||
+ | </td> | ||
+ | <td> | ||
+ | <p>417.9918</p> | ||
+ | </td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td> | ||
+ | <div align="right"><p>mu_1</p></div> | ||
+ | </td> | ||
+ | <td> | ||
+ | <p>0.00092593</p> | ||
+ | </td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td> | ||
+ | <div align="right"><p>mu_2</p></div> | ||
+ | </td> | ||
+ | <td> | ||
+ | <p>0.00024995</p> | ||
+ | </td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td> | ||
+ | <div align="right"><p>mu_3</p></div> | ||
+ | </td> | ||
+ | <td> | ||
+ | <p>0.0004495</p> | ||
+ | </td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td> | ||
+ | <div align="right"><p>mu_4</p></div> | ||
+ | </td> | ||
+ | <td> | ||
+ | <p>0.00025</p> | ||
+ | </td> | ||
+ | </tr> | ||
+ | </table> | ||
+ | </div> | ||
+ | <h3>References</h3> | ||
+ | <br> | ||
+ | <a name="ref1"></a> | ||
+ | <p><sup style="font-size:10px">[1]</sup> HEric L. Haseltine and Frances H. Arnold (2007) Synthetic Gene Circuits: Design with Directed Evolution. Annual Review Biophys Biomol Struct. 2007;36:1-19. </p> | ||
+ | <br> | ||
+ | <p><sup style="font-size:10px">[2]</sup> Andrianantoandro, E., Basu, S., Karig, D.K., Weiss, R. Synthetic biology: New engineering rules for an emerging discipline (2006) Molecular Systems Biology, 2, art. no. msb4100073, pp. 2006.0028 Biol. Chem., Vol. 266, No.1, pp.71-75 </p> | ||
+ | <br> | ||
+ | <p><sup style="font-size:10px">[3]</sup> Alon, U. (2006) An Introduction to Systems Biology: Design Principles of Biological Circuits. </p> | ||
+ | <br> | ||
+ | <p><sup style="font-size:10px">[4]</sup> Chen, B.-S., Chang, C.-H., Lee, H.-C. Robust synthetic biology design: Stochastic game theory approach (2009) Bioinformatics, 25 (14), pp. 1822-1830. </p> | ||
+ | <br> | ||
+ | <p><sup style="font-size:10px">[5]</sup> Guntas, G; Mansell, TJ, Kim, JR, Ostermeier, M. Directed evolution of protein switches and their application to the creation of ligan-binding proteins (2005). Procedings of the National Academy of Sciences of the United States of America, 102 (32): 11224-11229 </p> | ||
+ | <br> | ||
+ | <p><sup style="font-size:10px">[6]</sup> Marchisio, MA (Marchisio, M. A.); Stelling, J (Stelling, J. Computational design of synthetic gene circuits with composable parts (2008). IOINFORMATICS, 24 (17): 1903-1910 </p> | ||
+ | <br> | ||
+ | <p><sup style="font-size:10px">[7]</sup> Johannes, TW; Zhao, HM. Johannes, TW; Zhao, HM (2006). CURRENT OPINION IN MICROBIOLOGY, 9 (3): 261-267 </p> | ||
+ | <br> | ||
+ | <p><sup style="font-size:10px">[8]</sup> Nowak MA (2006). Evolutionary Dynamics: Exploring the Equations of Life.Harvard University Press. (Excerpt, Nature review, Science review,American Scientist review, Quarterly Review of Biology article, Science Daily, R.R. Hawkins Award, 2007 Book of Distinction) </p> | ||
+ | <br> | ||
+ | <p><sup style="font-size:10px">[9]</sup> <a href="http://www.bio.davidson.edu/courses/genomics/method/FACS.html">Biology @ Davidson, Davidson College Academics.</a> </p> | ||
+ | <br><br> | ||
+ | <hr> | ||
+ | <table class="nav"> | ||
+ | <tr> | ||
+ | <td> | ||
+ | <a href="https://2010.igem.org/Team:Aberdeen_Scotland/Probability"><img src="https://static.igem.org/mediawiki/2010/8/8e/Left_arrow.png"> Return to Determination of the the Hill Coefficient n<sub style=font-size:10px">2</sub></a> | ||
+ | </td> | ||
+ | <td align="right"> | ||
+ | <a href="https://2010.igem.org/Team:Aberdeen_Scotland/Switch_Characterisation">Continue to the Switch Characterisation Lab Diary <img src="https://static.igem.org/mediawiki/2010/3/36/Right_arrow.png"></a> | ||
+ | </td> | ||
+ | </tr> | ||
+ | </table> | ||
</html> | </html> | ||
+ | {{:Team:Aberdeen_Scotland/Footer}} |
Latest revision as of 22:23, 27 October 2010
University of Aberdeen - ayeSwitch
What is directed evolution?
In protein engineering, the method of directed evolution is used to harness the power of natural selection. It is used to evolve certain proteins or RNA with desirable properties not found in nature. Rather than swapping out different components of the structure, directed evolution can be used to manipulate the behaviour of these components. Then, by screening or selection, we can obtain the desired phenotype.
There are three typical steps involved that include:
1. Diversification: The gene that encodes the desired protein is mutated and recombined at random to create a large library of gene mutants. Techniques commonly used in this step involve PCR and DNA shuffling.
2. Selection: The library is tested for the existence of mutants possessing the property of interest using screening or selection.
a. Selections automatically remove all mutants that do not function.
b. Screening enables the researcher to manually detect and isolate the high-performing mutants.
3. Amplifications: The mutations identified int he selection or screening process are replicated several times so it allows the researchers to sequence their DNA and recognize what mutations have occurred.
For optimizing efficiency, altering more than one parameter would be preferable in order to obtain the desired protein. Therefore, we accululate mutations over multiple generations.
In order to examine how changes in the parameters were affecting the stability of the steady states, we could do uniform bifurcation analysis e.g. determining th regions of stability.
In the experiments, we were able to use the FACS (fluorescence activated cells sorting) machine to separate the proteins with the brightest colours and thus use the strongest ones. Once we had this, we could put the proteins through the system again as initial conditions and then we would have a more distinctive switch.
Pros
Directed evolution can systematically perturb and distinguish components then, using these perturbed components, we can gain further insight into the functionality of the normal circuit.
If some of the pathways enhance the desired performance, then iterative rounds of random mutation of the circuit and screening for the desired properties may provide steady improvement. This only works if the initial design is close to a better design.
Why use maths to model directed evolution?
There are three primary limitations that mathematical modelling can overcome:
1. The evolutionary space for a genetic circuit is too large to explore efficiently.
2. Refining our parameters is much easier than improving (or altering) the system experimentally.
3. Although selecting for independent properties is possible, it usually requires setting up multiple rounds of screening or selection.
Procedure
1. We first took the basic program where we integrated the system, plugging in an estimation of the parameters.
2. For each parameter, it yielded the value for which we obtain the largest difference between GFP and CFP.
3. We then fixed our chosen parameter and continue the program with the next parameter.
4. This can be done over and over for each parameter of the model, eventually producing a clearly defined switch that will somewhat resemble our original system.
5. Repeat these steps for each parameter in the system to eventually produce a clearly optimized switch that will yield the maximum difference between GFP and CFP when switching from one state to the other.
Results
Fig 1. shows the behaviour of the original system over a set period of time.
Fig 1.
As shown, the system has been integrated over a long enough timespan such that it has reached equilibrium.
Once we had this, we could then do this again but change one parameter and instead insert a range of values that can be put through the system.
Then we found the value for that parameter that produced the largest difference between the two proteins. We did this for each parameter until we had an “ideal” value for each one. Finally, we entered these values into the system and integrated it again.
Fig 2.
Fig 3.
Fig 2. shows our final result of the mRNAs and the two proteins. We can see at first that CFP was the brightest protein initially, then we added more galactose and, at around time = 50, we can clearly see the decline of CFP and the increase of GFP. Fig 3. shows how copper and galactose are acting in the system. Thus, a clear switch in our system. The results for each parameter were as follows;
Parameter |
Optimal Valued |
K_1 |
5000 |
K_2 |
10000 |
K_3 |
5000 |
K_4 |
100000 |
lambda_1 |
133.9940 |
lambda_2 |
415.9916 |
lambda_3 |
583.9984 |
lambda_4 |
417.9918 |
mu_1 |
0.00092593 |
mu_2 |
0.00024995 |
mu_3 |
0.0004495 |
mu_4 |
0.00025 |
References
[1] HEric L. Haseltine and Frances H. Arnold (2007) Synthetic Gene Circuits: Design with Directed Evolution. Annual Review Biophys Biomol Struct. 2007;36:1-19.
[2] Andrianantoandro, E., Basu, S., Karig, D.K., Weiss, R. Synthetic biology: New engineering rules for an emerging discipline (2006) Molecular Systems Biology, 2, art. no. msb4100073, pp. 2006.0028 Biol. Chem., Vol. 266, No.1, pp.71-75
[3] Alon, U. (2006) An Introduction to Systems Biology: Design Principles of Biological Circuits.
[4] Chen, B.-S., Chang, C.-H., Lee, H.-C. Robust synthetic biology design: Stochastic game theory approach (2009) Bioinformatics, 25 (14), pp. 1822-1830.
[5] Guntas, G; Mansell, TJ, Kim, JR, Ostermeier, M. Directed evolution of protein switches and their application to the creation of ligan-binding proteins (2005). Procedings of the National Academy of Sciences of the United States of America, 102 (32): 11224-11229
[6] Marchisio, MA (Marchisio, M. A.); Stelling, J (Stelling, J. Computational design of synthetic gene circuits with composable parts (2008). IOINFORMATICS, 24 (17): 1903-1910
[7] Johannes, TW; Zhao, HM. Johannes, TW; Zhao, HM (2006). CURRENT OPINION IN MICROBIOLOGY, 9 (3): 261-267
[8] Nowak MA (2006). Evolutionary Dynamics: Exploring the Equations of Life.Harvard University Press. (Excerpt, Nature review, Science review,American Scientist review, Quarterly Review of Biology article, Science Daily, R.R. Hawkins Award, 2007 Book of Distinction)
[9] Biology @ Davidson, Davidson College Academics.
Return to Determination of the the Hill Coefficient n2 | Continue to the Switch Characterisation Lab Diary |