Team:Aberdeen Scotland/Fixed Points

From 2010.igem.org

(Difference between revisions)
 
(10 intermediate revisions not shown)
Line 3: Line 3:
<html>
<html>
<h1>Fixed Points</h1>
<h1>Fixed Points</h1>
-
<p>Fixed points are the points where and equation's rate of change, or slope, is zero. there are three main types of equilibrium points: stable, unstable, and saddle-node points. A stable equilibrium is a value towards which the function converges, whereas an unstable equilibrium is a value away from which the function will diverge. A saddle-node is here the function both converges towards and diverges away from (Wikipedia).</p>
+
<p>Fixed points are the points where an equation’s rate of change, or its slope, is zero. There are three main types of equilibrium points: stable, unstable, and saddle-node points. A stable equilibrium is a value which the function converges towards, whereas an unstable equilibrium is a value the function will diverge from. A saddle-node is where the function both converges to and diverges from, the equilibrium, depending on the direction from which we approach the point (S. Strogatz, Nonlinear dynamics and Chaos).</p>
<br>
<br>
<h3>Why are fixed points important?</h3>
<h3>Why are fixed points important?</h3>
-
<p>fixed points are determined by setting all differential equations in a system equal to zero and solving for the variable being analysed. The numerical method for directly calculating fixed points is to find the roots of the system of equations using built-in root finding function such as <i>fzero</i> in MATLAB or similar functions in Maple and C. The indirect method is to plot and find the intersections of these nullclines, which represent fixed points.</p>
+
<p>Finding fixed points is the first step in analyzing the stability of a system. In particular, we are interested in the nature of our system’s bistability and how it changes with a variation of the parameters, i.e., its bifurcation analysis. This analysis can be conveyed to the biologists to minimize experimental guessing.
 +
</p>
<br>
<br>
 +
<h3>How are fixed points calculated?</h3>
 +
<p>Fixed points are determined by setting all differential equations in a system equal to zero and solving for the variable being analyzed. The direct method of calculating equilibrium points is to find the roots of the system of equations using built-in root-finding functions such as fzero in MATLAB, and similar functions in Maple and C. The indirect method of calculating the equilibrium points is to plot and find the intersections of these nullclines, which represent equilibrium points.</p>
 +
<h3>Results</h3>
<h3>Results</h3>
-
<p>We were able to find the equilibrium points of the system analytically (for Hill coefficients) and computationally. These points were used for bifurcation analysis, and for analysing the probability that our system would exhibit istable behaviour.</p>
+
<p>We were able to find the fixed points of the system analytically (for small Hill coefficients) and computationally. These points were used for bifurcation analysis, and for analyzing the probability that our system would exhibit bistable behavior.</p>
<br>
<br>
<br>
<br>
Line 18: Line 22:
<br>
<br>
<h3>Why are nullclines important?</h3>
<h3>Why are nullclines important?</h3>
-
<p>the intersections of the cullclines gives the equilibrium points of the system of differential equations. From graphs of the sullclines, it is obvious whether or not a system will be bistable. If the nullclines only intersect in one place the sestem is not viable. If there are more than two interesctions, the middle equilibrium point is often an unstable saddle point.</p>
+
<p>The intersections of the nullclines give the equilibrium points of the system of differential equations. From graphs of the nullclines, it is possible to infer whether or not a system will be bistable. If the nullclines only intersect in one place the system is not bistable, since there is one single equilibrium point. If there are more than two intersections, the middle equilibrium point is often an unstable saddle point.</p>
<br>
<br>
<h3>How are the nullclines calculated?</h3>
<h3>How are the nullclines calculated?</h3>
-
<p>Just as in calculating the equilibrium points, we set the governing differential equations of the system equal to zero and plot the curves generated.</p>
+
<p>Just as in calculating the fixed points, we set the governing differential equations of the system equal to zero and plot the curves generated.</p>
 +
<center>
 +
<img src="https://static.igem.org/mediawiki/2010/e/e4/Equations3.png">
 +
</center>
<br>
<br>
<h3>Simulations involved</h3>
<h3>Simulations involved</h3>
 +
 +
 +
<table>
<table>
   <tr>
   <tr>
     <td>
     <td>
-
     nullclines1.png
+
     <p><img src="https://static.igem.org/mediawiki/2010/1/18/Nullclines_diagram1.png"/></p>
     </td>
     </td>
     <td>
     <td>
-
     <p><b>Figure 1.</b> This figure is a plot of how much galactose (in green) and methionine (in blue) must be added to produce GFP and CFP. In this plot, both Hill coefficients are equal to two. The nullclines cross over at three equilibrium points, where the middle is a saddle-node fixed point. This is an 'ideal' bistability plot.</p>
+
     <p><b>Figure 1.</b>This figure is a plot of the nullclines of the differential equations for CFP and GFP, where we solve them for CFP as a function of GFP. In this plot, both Hill coefficients are two, ie. n2 and n4 are both 2. The nullclines cross over at three fixed points, where the middle is a saddle-node fixed point. This is an ‘ideal’ bistability plot.
 +
</p>
     </td>
     </td>
   </tr>
   </tr>
   <tr>
   <tr>
     <td>
     <td>
-
     nullclines2.png
+
     <p><img src="https://static.igem.org/mediawiki/2010/5/53/Nullclines_diagram2.png"></p>
     </td>
     </td>
     <td>
     <td>
-
     <p><b>Figure 2.</b> This figure again shows the required amounts of galactose and methionine to produce GFP and CFP respectively. In this plot, both Hill cefficients are equal to one. Note that the nullclines only cross once, resulting in only one fixed point and hence no possibility of bistability. This result allowed us to tell the biologists that the toggle switch would definitely <i>not</i> work if both of the Hill coefficients were equal to one.</p>
+
     <p><b>Figure 2.</b>This figure is a plot, again, of the nullclines of the differential equations for CFP and GFP. In this plot, both Hill coefficients are one. Note that the nullclines only cross once, resulting in only one fixed point and hence no possibility of bistability. This result allowed us to tell the biologists that the toggle switch would definitely not work if both of the Hill coefficients were one.
 +
</p>
     </td>
     </td>
   </tr>
   </tr>
   <tr>
   <tr>
     <td>
     <td>
-
     <p><b>Figure 3.</b> This figure again shows the concentrations of galactose and methionine to produce GFP and CFP. In this graph, we experimented with setting the initial values of galactose and methionine and seeing how the intersection points are offset. The blue line represents the GFP = CFP line.  The further from this line the equilibrium points are, the harder it will be to switch between stable states.</p>
+
     <p><img src="https://static.igem.org/mediawiki/2010/d/da/Nullclines_diagram3.png"/></p>
     </td>
     </td>
     <td>
     <td>
-
     nullclines3.png
+
     <p><b>Figure 3.</b> Figure 3: This figure is of CFP as a function of GFP for the two nullclines when the initial galactose and copper concentrations are changed. From this we can see how the fixed points change with different initial galactose and methionine concentrations. The blue line represents the GFP = CFP line. The further from this line equilibrium points are, the harder it will be to switch between stable states.</p>
     </td>
     </td>
   </tr>
   </tr>
Line 52: Line 64:
<br>
<br>
<h3>Results</h3>
<h3>Results</h3>
-
<p>We focused our efforts on plotting the mullclines of the system for a range of Hill coefficient combinations in order to get a general idea of which combination would most likely produce stable bistability. We cound that all combinations gave bistability, except when both Hill coefficents were equal to one. the optimal Hill coefficient combination occurred when both Hill coefficients were equal to two. We informed the biology team of this, and that is we wanted our system to successfully switch, we could not have both Hill coefficents equal to one.</p>
+
<p>We focused our efforts on plotting the nullclines of the system for a range of Hill coefficient combinations in order to get a general idea of which combinations would most likely produce robust bistability. We found that all combinations gave bistability, except when both Hill coefficients were one. The optimal Hill coefficient combination occurred when both Hill coefficients were two. We passed this information onto the biology team, letting them know that if we wanted the system to successfully switch, we could not have both Hill coefficients with a value of one.
 +
</p>
<hr>
<hr>
<table class="nav">
<table class="nav">
Line 60: Line 73:
</td>
</td>
<td align="right">
<td align="right">
-
<a href="https://2010.igem.org/Team:Aberdeen_Scotland/Bifurcation">Continue to Bifurcation&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/Bifurcation">Continue to Bifurcation and Stability&nbsp;&nbsp;<img src="https://static.igem.org/mediawiki/2010/3/36/Right_arrow.png"></a>
</td>
</td>
</tr>
</tr>

Latest revision as of 20:31, 27 October 2010

University of Aberdeen - ayeSwitch - iGEM 2010

Fixed Points

Fixed points are the points where an equation’s rate of change, or its slope, is zero. There are three main types of equilibrium points: stable, unstable, and saddle-node points. A stable equilibrium is a value which the function converges towards, whereas an unstable equilibrium is a value the function will diverge from. A saddle-node is where the function both converges to and diverges from, the equilibrium, depending on the direction from which we approach the point (S. Strogatz, Nonlinear dynamics and Chaos).


Why are fixed points important?

Finding fixed points is the first step in analyzing the stability of a system. In particular, we are interested in the nature of our system’s bistability and how it changes with a variation of the parameters, i.e., its bifurcation analysis. This analysis can be conveyed to the biologists to minimize experimental guessing.


How are fixed points calculated?

Fixed points are determined by setting all differential equations in a system equal to zero and solving for the variable being analyzed. The direct method of calculating equilibrium points is to find the roots of the system of equations using built-in root-finding functions such as fzero in MATLAB, and similar functions in Maple and C. The indirect method of calculating the equilibrium points is to plot and find the intersections of these nullclines, which represent equilibrium points.

Results

We were able to find the fixed points of the system analytically (for small Hill coefficients) and computationally. These points were used for bifurcation analysis, and for analyzing the probability that our system would exhibit bistable behavior.



Nullclines

In a system of differential equations, the nullclines are the solution curves for which all of the differential equations are equal to zero (Wikipedia).


Why are nullclines important?

The intersections of the nullclines give the equilibrium points of the system of differential equations. From graphs of the nullclines, it is possible to infer whether or not a system will be bistable. If the nullclines only intersect in one place the system is not bistable, since there is one single equilibrium point. If there are more than two intersections, the middle equilibrium point is often an unstable saddle point.


How are the nullclines calculated?

Just as in calculating the fixed points, we set the governing differential equations of the system equal to zero and plot the curves generated.


Simulations involved

Figure 1.This figure is a plot of the nullclines of the differential equations for CFP and GFP, where we solve them for CFP as a function of GFP. In this plot, both Hill coefficients are two, ie. n2 and n4 are both 2. The nullclines cross over at three fixed points, where the middle is a saddle-node fixed point. This is an ‘ideal’ bistability plot.

Figure 2.This figure is a plot, again, of the nullclines of the differential equations for CFP and GFP. In this plot, both Hill coefficients are one. Note that the nullclines only cross once, resulting in only one fixed point and hence no possibility of bistability. This result allowed us to tell the biologists that the toggle switch would definitely not work if both of the Hill coefficients were one.

Figure 3. Figure 3: This figure is of CFP as a function of GFP for the two nullclines when the initial galactose and copper concentrations are changed. From this we can see how the fixed points change with different initial galactose and methionine concentrations. The blue line represents the GFP = CFP line. The further from this line equilibrium points are, the harder it will be to switch between stable states.


Results

We focused our efforts on plotting the nullclines of the system for a range of Hill coefficient combinations in order to get a general idea of which combinations would most likely produce robust bistability. We found that all combinations gave bistability, except when both Hill coefficients were one. The optimal Hill coefficient combination occurred when both Hill coefficients were two. We passed this information onto the biology team, letting them know that if we wanted the system to successfully switch, we could not have both Hill coefficients with a value of one.



Back to the Top