Team:RHIT/Modeling

From 2012.igem.org

(Difference between revisions)
 
(23 intermediate revisions not shown)
Line 14: Line 14:
                                         $('#logged_in').show();
                                         $('#logged_in').show();
                                 }   
                                 }   
-
$('#rhit-sheet2').hide();   
+
$('#rhit-sheet2').hide();   
 +
var split = location.search.replace('?', '').split('=');
 +
window[split[1]]();
});
});
function page1() {
function page1() {
Line 42: Line 44:
<div class="rhit-content">
<div class="rhit-content">
<div class="rhit-header">
<div class="rhit-header">
-
<img src="https://static.igem.org/mediawiki/igem.org/8/87/Banner_only2.png" width="100%"/>
+
<a href="https://2012.igem.org/Team:RHIT"><img src="https://static.igem.org/mediawiki/igem.org/8/87/Banner_only2.png" width="100%"/></a>
<div id="logged_out">
<div id="logged_out">
<a href="https://2012.igem.org/Main_Page">iGEM Home</a>
<a href="https://2012.igem.org/Main_Page">iGEM Home</a>
Line 96: Line 98:
<p>In addition to the synthetic biology project, the team also pursued the creation of a mathematical model of the modified biological system. The first step of the modeling process was to gain an understanding of the various chemical and biological agents, how they interacted with each other, and ultimately how these interactions caused the physiological changes in the cell. This was done by studying the mechanism presented in Bradwell's paper, <i>A walk-through of the yeast pheromone response pathway</i>. Once the key components were identified, a simplified explanation of the process of interest was created. Below is the mechanism the team used for the basis of the model:</p><br /><br />
<p>In addition to the synthetic biology project, the team also pursued the creation of a mathematical model of the modified biological system. The first step of the modeling process was to gain an understanding of the various chemical and biological agents, how they interacted with each other, and ultimately how these interactions caused the physiological changes in the cell. This was done by studying the mechanism presented in Bradwell's paper, <i>A walk-through of the yeast pheromone response pathway</i>. Once the key components were identified, a simplified explanation of the process of interest was created. Below is the mechanism the team used for the basis of the model:</p><br /><br />
<img src="https://static.igem.org/mediawiki/igem.org/3/32/Quiagen_Yeast_Mating_Pheromone-Response_Pathway.png" width="100%" />
<img src="https://static.igem.org/mediawiki/igem.org/3/32/Quiagen_Yeast_Mating_Pheromone-Response_Pathway.png" width="100%" />
-
<ul="rhit-pathwayList">
 
-
<li>Mating Factor Binds to Ste2/Ste3 homodimer (Ste2 for mat alpha / Ste3 for mat a)</li>
 
-
<li>Gpa1 subunit exchanges GDP for GTP</li>
 
-
<li>Gpa1 releases Ste4/G-gamma heterodimer</li>
 
-
<li>Ste4/Ste18 (through Ste4), transmits the signal to the following complexes, Ste5/Ste11, Ste20 protein kinase, and a Far1/Cdc24 complex</li>
 
-
<li>Ste4 binds to Cdc24, Far1 brings Cdc24 to Cdc42, Cdc24 facilitates the exchange of GDP for GTP, causing activated Cdc42 to bind to the complex between Ste20 and Bem1, resulting in the activation of Ste20
 
-
* Ste20 is in a low-activity state before being bound by Cdc42</li>
 
-
<li>Ste4 binds to Ste5, Ste5 acts as an adaptor to bring Ste4-beta and Ste11 toward Ste20</li>
 
-
<li>Ste20 phosphorylates Ste11 (aided by Ste50), causing it to activate</li>
 
-
<li>Ste11 then activates Ste7 by phosphorylating it</li>
 
-
<li>Ste7 in turn activates Fus3 and Kss1 by phosphorylating them
 
-
*Ste5 also plays a role in all both of these events</li>
 
-
<li>Fus3 and Kss1 phosphorylate Ste12, Dig1 and Dig2, in their phosphorylated state Dig1 and Dig2 release activated Ste12
 
-
*Kss1 holds the Dig1/Dig2 complex together with Ste12</li>
 
-
<li>Activated Ste12 binds to DNA to facilitate transcription of pheromone response genes</li><br />
 
-
</ul>
 
<p>After laying out this process, the team began making evaluations as to which parts of the process would be important, and what kind of questions could be answered by the model. After careful deliberation and discussion, the team decided that the questions of interest to the project were to determine the sensitivity of the circuit and the time after exposure that fluorescence is detectable. With those questions in mind, the system was broken into four distinct portions: initiation, transduction, production, and regulation.</p>
<p>After laying out this process, the team began making evaluations as to which parts of the process would be important, and what kind of questions could be answered by the model. After careful deliberation and discussion, the team decided that the questions of interest to the project were to determine the sensitivity of the circuit and the time after exposure that fluorescence is detectable. With those questions in mind, the system was broken into four distinct portions: initiation, transduction, production, and regulation.</p>
<p>The initiation portion of the model pertained to the interactions between the pheromone and the receptor. Since one of the goals of the model was to address the sensitivity of the construct, this portion was mostly conserved and incorporated into the model. The transduction portion covered most of the unmodified kinase cascade covered in the mechanism. Because of the relatively short time period of these interactions in comparison to the approximated time frame of the circuit, and the well characterized and unchanged nature of this part of the pathway, these interactions were removed from the model. The third portion of the model covered the transcriptional and translational activities necessary for the production of the construct. The final module of the model was accounting for the regulatory functions contained in the construct.</p>
<p>The initiation portion of the model pertained to the interactions between the pheromone and the receptor. Since one of the goals of the model was to address the sensitivity of the construct, this portion was mostly conserved and incorporated into the model. The transduction portion covered most of the unmodified kinase cascade covered in the mechanism. Because of the relatively short time period of these interactions in comparison to the approximated time frame of the circuit, and the well characterized and unchanged nature of this part of the pathway, these interactions were removed from the model. The third portion of the model covered the transcriptional and translational activities necessary for the production of the construct. The final module of the model was accounting for the regulatory functions contained in the construct.</p>
</div></a>
</div></a>
<div class="rhit-redSheet" id="rhit-redSheet">
<div class="rhit-redSheet" id="rhit-redSheet">
-
<p>Having defined what the model was supposed to answer, and what parts of the process were going to be included, the next step was to design a system of differential equations that would account for the various pieces of the model. The first draft of this system is shown below, along with a brief description of what each equation represents, and what the various terms in each equation account for.</p><br />
+
<p>For the team’s mathematical modeling section, two differential models of the system were created: a case with one or more independent binding sites and a case that exhibits cooperativity with various binding sites. Both of these systems of equations are shown below. The behavior of the steady-state solutions of the system of equations were analyzed using algebraic and graphical methods. This analysis provided insight into the possible results from this kind of system.</p><br />
-
<br />dP/dt eqn img here<br /><br />
+
<div align="center"><img src="https://static.igem.org/mediawiki/igem.org/2/24/DiffRHIT1.png" width="40%"/></div>
-
<p>This equation represented the pheromone concentration, (matA or mata) available to interact with the receptors of the cell. The equation was designed to account for the pheromones secreted by the yeast cells, the binding of the pheromone to the receptor, and the possible unbinding of the pheromone from the receptor. The secretion of the pheromone from one yeast cell was assumed to be independent of the changes exhibited by the other yeast cells. It was also assumed that this pheromone was equally distributed around the whole environment of the cells.</p><br />
+
<div align="center"><h4>Figure 1. Circuit diagram.</h4></div>
-
<br />dB/dt eqn img here<br />
+
<div align="center"><img src="https://static.igem.org/mediawiki/igem.org/e/ed/DiffRHITnew.PNG" width="40%"/></div>
-
<br />dU/dt eqn img here<br /><br />
+
<div align="center"><h4>Figure 2. Independent binding sites.</h4></div>
-
These two equations represented the two possible states of the pheromone receptor, (Ste2 or Ste3) bound to pheromone or unbound. While these proteins like most proteins will decay overtime and must be reproduced, due to the relatively small time scale that these proteins significantly contribute to the behavior of the system, these terms were assumed to be negligible. Furthermore, this model also assumes that all unbound receptors are capable of being bound and that all bound receptors are capable of signal transduction.</p><br />
+
<div align="center"><img src="https://static.igem.org/mediawiki/igem.org/4/44/RHITdiff3.PNG" width="40%"/></div>
-
<br />dS/dt eqn img here<br />
+
<div align="center"><h4>Figure 3. Perfect cooperativity.</h4></div>
-
<br />dA/dt eqn img here<br /><br />
+
<p>During the analysis of these systems, the value of S became an area of interest. In order to determine the nature of this function, the biological system from which the S is derived was studied. In order to study this system, a paper written by Kofahl and Klipp’s titled <i>Modeling the dynamics of the yeast pheromone pathway</i> was analyzed. There were two methods employed to analyze this model: the deterministic model and a stochastic model. The deterministic model detailed in the paper was recreated and the level of activated Ste12 was measured as a function of the amount of mating factor. The stochastic model was created by the NTNU team. Using Mohan, Shahrezaei, et. al.’s paper <i>The scaffold protein Ste5 directly controls a switch-like mating decision in yeast</i>, the Norwegian team made a direct correlation between the activation of Fus3 and the amount of mating factor, and measured activated Ste12 activity as a function of Fus3 activation. The results of these models indicated that for a given amount of mating factor, a constant level of Ste12 is produced. The diagrams of these models are shown below:</p><br />
-
<p>These two equations represented a hypothetical signal that directly stimulated the production of the team's fluorescent hetero-transcription factor. The S equation signified the series of protein kinase interactions, caused by the production of bound pheromone receptor and a decay of the signal as the proteins are turned off. The A equation represented the total sum of the factors contributing to the synthesis of the construct, including the signal from the S equation, the auto-regulation of the construct, and the loss of the signal as the construct is synthesized.</p><br />
+
<div align="center"><img src="https://static.igem.org/mediawiki/igem.org/3/32/DiffRHIT6.png" width="65%"/></div>
-
<br />dX/dt eqn img here<br /><br />
+
<div align="center"><h4>Figure 4. Steady state response diagrams.</h4></div>
-
<p>Finally this equation represented the concentration of the synthetic hetero-transcription factor that is produced by the signal, and is lost over time due to decay.</p><br />
+
<p>The figures below illustrate the three possible solutions predicted by the mathematical model. The x-axis is a measure of the external mating pheromone concentration, and the y-axis represents the amount of protein. The first illustration depicts a system where once the signal passes a particular threshold the circuit is turned on and will continue to produce protein. The second depicts a system where there is a particular threshold of mating pheromone required to bring about a stable level of protein; if the signal falls below this level, the protein production is transient and returns to zero. The third depicts a system in which there is no discontinuity in protein production, and signal is required for all protein production.</p><br />
-
<p>After creating this simple model, the team began doing research on the kinetic rate constants and binding affinities for the various parts of the model. During this research, the team came across a published model from Kofahl and Klipp's paper, <i>Modeling the dynamics of the yeast pheromone pathway</i>.</p><br />
+
<div align="center"><img src="https://static.igem.org/mediawiki/igem.org/d/d9/DiffRHIT7.png" width="80%"/></div>
 +
<div align="center"><h4>Figure 5. Bifurcation diagrams.</h4></div>
 +
<p>While it is not yet possible to verify which, if any, accurately represents the actual system, the created models predict possible scenarios that would allow for the success of the project. Furthermore, the model predicts that the success of this project is dependent solely on the values of the parameters; specifically, the decay terms must be less than the production terms of the system. The results of the model make good biological and intuitive sense. In order for the project to be successful, either the first or second depictions must hold true. The analysis, derivation, future work, and all the work leading to these conclusions are listed below in the named sections.</p>
 +
<p>For an in depth look at the work and analysis that went into developing the mathematical model for this system download <a href="https://static.igem.org/mediawiki/igem.org/8/85/Differential_Model.pdf">this pdf</a>.</p>
</div>
</div>
<div class="rhit-grnSheet" id="rhit-grnSheet">
<div class="rhit-grnSheet" id="rhit-grnSheet">
-
Stochastic text
+
<p>The team had a collaborative effort with the iGEM team from Norwegian University of Science and Technology (NTNU). Their contribution to the project was the development of a stochastic model of the yeast pheromone response pathway presented in the Kofahl and Klipp’s paper. In particular, the NTNU team was looking at the system’s Ste12 activation response from varying amounts of initial mating factor. Using a paper from <i>Nature</i>, the team discovered that there was a directly proportional relationship between the amount of activated Ste12 and the amount of free Fus3. This allowed for a substantial simplification of the system. The result of their stochastic simulations for 250 micromolar mating factor is shown below:</p><br />
 +
<div align="center"><img src="https://static.igem.org/mediawiki/igem.org/e/e5/Stochastic_model_1.png" width=40%/></div><br />
 +
<p>This plot shows that a given amount of mating factor will result in a fairly steady amount of activated Ste12. Taking the average of the various values produced by the simulation allowed for the production of a dose-response curve of the system. This curve is shown below:</p>
 +
<div align="center"><img src="https://static.igem.org/mediawiki/igem.org/1/17/Stochastic_model_2.png" width=40%/></div><br />
 +
<p>This result is significant because it supports the differential model’s conclusion that for a particular amount of mating factor, a relatively constant level of Ste12 activation is observed. This provided the team with more evidence for the assumption that the original input signal can be treated as a gradually decaying constant.</p><br /><br />
</div>
</div>
</div>
</div>

Latest revision as of 03:54, 4 October 2012

Home
Team
Project
BioBricks
Modeling
Notebook
Outreach
Planning Process
Differential Model
Stochastic Model

In addition to the synthetic biology project, the team also pursued the creation of a mathematical model of the modified biological system. The first step of the modeling process was to gain an understanding of the various chemical and biological agents, how they interacted with each other, and ultimately how these interactions caused the physiological changes in the cell. This was done by studying the mechanism presented in Bradwell's paper, A walk-through of the yeast pheromone response pathway. Once the key components were identified, a simplified explanation of the process of interest was created. Below is the mechanism the team used for the basis of the model:



After laying out this process, the team began making evaluations as to which parts of the process would be important, and what kind of questions could be answered by the model. After careful deliberation and discussion, the team decided that the questions of interest to the project were to determine the sensitivity of the circuit and the time after exposure that fluorescence is detectable. With those questions in mind, the system was broken into four distinct portions: initiation, transduction, production, and regulation.

The initiation portion of the model pertained to the interactions between the pheromone and the receptor. Since one of the goals of the model was to address the sensitivity of the construct, this portion was mostly conserved and incorporated into the model. The transduction portion covered most of the unmodified kinase cascade covered in the mechanism. Because of the relatively short time period of these interactions in comparison to the approximated time frame of the circuit, and the well characterized and unchanged nature of this part of the pathway, these interactions were removed from the model. The third portion of the model covered the transcriptional and translational activities necessary for the production of the construct. The final module of the model was accounting for the regulatory functions contained in the construct.

For the team’s mathematical modeling section, two differential models of the system were created: a case with one or more independent binding sites and a case that exhibits cooperativity with various binding sites. Both of these systems of equations are shown below. The behavior of the steady-state solutions of the system of equations were analyzed using algebraic and graphical methods. This analysis provided insight into the possible results from this kind of system.


Figure 1. Circuit diagram.

Figure 2. Independent binding sites.

Figure 3. Perfect cooperativity.

During the analysis of these systems, the value of S became an area of interest. In order to determine the nature of this function, the biological system from which the S is derived was studied. In order to study this system, a paper written by Kofahl and Klipp’s titled Modeling the dynamics of the yeast pheromone pathway was analyzed. There were two methods employed to analyze this model: the deterministic model and a stochastic model. The deterministic model detailed in the paper was recreated and the level of activated Ste12 was measured as a function of the amount of mating factor. The stochastic model was created by the NTNU team. Using Mohan, Shahrezaei, et. al.’s paper The scaffold protein Ste5 directly controls a switch-like mating decision in yeast, the Norwegian team made a direct correlation between the activation of Fus3 and the amount of mating factor, and measured activated Ste12 activity as a function of Fus3 activation. The results of these models indicated that for a given amount of mating factor, a constant level of Ste12 is produced. The diagrams of these models are shown below:


Figure 4. Steady state response diagrams.

The figures below illustrate the three possible solutions predicted by the mathematical model. The x-axis is a measure of the external mating pheromone concentration, and the y-axis represents the amount of protein. The first illustration depicts a system where once the signal passes a particular threshold the circuit is turned on and will continue to produce protein. The second depicts a system where there is a particular threshold of mating pheromone required to bring about a stable level of protein; if the signal falls below this level, the protein production is transient and returns to zero. The third depicts a system in which there is no discontinuity in protein production, and signal is required for all protein production.


Figure 5. Bifurcation diagrams.

While it is not yet possible to verify which, if any, accurately represents the actual system, the created models predict possible scenarios that would allow for the success of the project. Furthermore, the model predicts that the success of this project is dependent solely on the values of the parameters; specifically, the decay terms must be less than the production terms of the system. The results of the model make good biological and intuitive sense. In order for the project to be successful, either the first or second depictions must hold true. The analysis, derivation, future work, and all the work leading to these conclusions are listed below in the named sections.

For an in depth look at the work and analysis that went into developing the mathematical model for this system download this pdf.

The team had a collaborative effort with the iGEM team from Norwegian University of Science and Technology (NTNU). Their contribution to the project was the development of a stochastic model of the yeast pheromone response pathway presented in the Kofahl and Klipp’s paper. In particular, the NTNU team was looking at the system’s Ste12 activation response from varying amounts of initial mating factor. Using a paper from Nature, the team discovered that there was a directly proportional relationship between the amount of activated Ste12 and the amount of free Fus3. This allowed for a substantial simplification of the system. The result of their stochastic simulations for 250 micromolar mating factor is shown below:



This plot shows that a given amount of mating factor will result in a fairly steady amount of activated Ste12. Taking the average of the various values produced by the simulation allowed for the production of a dose-response curve of the system. This curve is shown below:


This result is significant because it supports the differential model’s conclusion that for a particular amount of mating factor, a relatively constant level of Ste12 activation is observed. This provided the team with more evidence for the assumption that the original input signal can be treated as a gradually decaying constant.



Retrieved from "http://2012.igem.org/Team:RHIT/Modeling"