Stability analysis

The

IN THE CASE OF STE2 DISSIPATION

In figure 3.1, there is only one steady state and no bifurcation. Moreover, the system is stable around the steady state.

Because of mass balance and dying out Ste2a (Ste2), the following pathway will recover to the initial concentration which leads to GFP degrade to zero asymptotically.

When analysing the system by linearization as above, the eigenvalues of full Jacobian A contains both zeros and values of negative real-part, leading to undertermined conclusion for original nonlinear system. However, because of mass balance, 6 out of total 17 species are conserved values including Sst2a (Sst2), Fus3a (Fus3), Ste12 (Ste12a), and GaGTP, C, D in the big cycle of combining G-protein and MAPK cascade. They are linear combinations of particle numbers that do not vary during the time evolution of the system. Each conservation relation can be used to eliminate one variable of the system, leading to a reduced system with a smaller number of variables. These variables are called the independent variables of the system; the dependent variables are defined as linear combinations of independent variables.Technically finding the conservation relation means finding rows in the stoichiometry matrix that can be expressed as linear combinations of other rows.

By using COPASI which uses Householder QR factorization [9] to reduce model, the Jacobian of reduced model can be calculated and its eigenvalues are

IN THE CASE OF STE2 ACCUMULATION OR CONSTANT

The pathway will continuously be activated, and GFP accumulates without constraints, so the system is unstable.

Sensitivity analysis

Values can be tuned: initial concentration of Ste2, Vht2hd, VGa2ht, VGFPmRNAsyn