Team:Aberdeen Scotland/Stability
From 2010.igem.org
University of Aberdeen - ayeSwitch
Bistability Analysis
The literal meaning of bifurcation is,m "division into two parts, branches, or aspects" (Answers). In mathematics, this 'division' is in reference to when a non-linear equation suddenly switches from having one solution to having multiple solutions when a parameter of the system is varied. The point at which bifurcation is referred to as the 'bifurcation point.' If a system is bistable it can be stable, or at rest in two distinct states.
bistability1.png
figure 1. This is a graphical representation of bistability. Ball 2 can be equally 'happy' in either the 'one' or the 'three' positions. Once it falls into one of these states, ball 2 will not bin inclined to leave this state, hence making it a 'stable state.'
Why are they important?
Analysing the stability of a system is ket to the success of a project. In biological terms, bistability represetns the potential for our system to switch between the expression fo gree and cyan fluarescent proteins. Our goal is to analyse the parameter space in order to determione which paramenters andre needed to allow for bistability in the system.
a bifurcation diagram has branches which define the possbile states a system can rest in. By analysing the bifurcation plots, we can determing how far apart the two 'branches,' or stable states, are. the further apart the branches, the harder it will be to attacin a switch between the two stable states, however the switch also becomes more stable with increased distance between the branchs. Therefore, we want to manipulate the parameters in the system to guarantee stable bistability.
How do you analyse the stability of a system?
In order to analyse the stability of the system, we utilised XPPAuto, C programming and MATLAB. We will attempt to detail the maths behind what the programme does computationally.
The process involved with general stability is as follws. We have taken a simple example system that is easy to analyse in order to demonstrate the general process:
Take the following example system of two differential equations:
bistability2.png
1. Set the differential equations equal to zero:
bistability3.png
2. Find the fixed points of the system of differential equations:
There is a fixed point at x=y=0. By observation, we can see that there are fixed points at x=0, y=2, and x=3, y=0. Finally, by solving this system of equations, we see that there is also a fixed point at x=1 and y=1.
3. Evaluate the Jacobian matrix of the system at these fixed points:
If we label the top fuction f1 and the bottom f2, the Jacobian matrix of the system will be of the firm:
bistability4.png
When we evaluate the Jacobian matrix for our example system of equations, we get:
bistability5.png
4. Using the Jacobian matrix, find the corresponding eigenvalues:
Definition: an eigenvalue is a value λ such that Ax = λx, where x is said to be an eigenvector of A. In order to solve for λ, solve the following equation:
bistability6.png
Where det is the deterinant of the matrix and I is the identity matrix. THis is known as the characteristic equation of a matrix.
In the case of our example system, the eigenvalues are:
5. Take the signs of these eigenvalues in order to analyse the stability:
The eigenvalues of a system give important information about the stability of the system. The stability of each of the fixed points can be analysed based on the following criteria:
Positive real and complex parts | = unstable |
Negative real and complex parts | = stable |
Positive real part | = unstable spiral (spiraling out) |
Negative real part | stable spiral (spiralling in) |
So, for our ixample system, it can be deduced that the stability at the fixed poitns is as follows:
One positive and one negative eigenvalue at (0,0) means that | (0,0) is a saddle node | |
Two negative eigenvalues at (0,2) means that | ||
(0,2) is a stable node | ||
Two negative eigenvalues at (3,0) means that | (3,0) is a stable node | |
One positive and one negative eigenvalue at (1,1) means that | (1,1) is a saddle node |
Simulations involved
Figure 2. This is a bifurcation diagram that we generated using Maple. It ranges over the parameter space of P1 (GFP) and C (a dimensionless group of parameters) and plots the fixed points o the system at those values.
Return to the Stochastic Model | Continue to Parameter Space Analysis |