A2 : Due by 11:30 AM, Oct. 8 via Catalyst
- Consider the reaction network 2X + Y → X + 2Y. Assume the rate constant for the reaction is 1.
- Find the stoichiometric matrix.
- Assume the system is mass action, and find the kinetics vector.
- Write out the ODE model.
- Simulate the ODEs from a couple of initial conditions and show (nicely formatted) plots of the resulting trajectories for each species.
- Determine a mass vector or prove that the system is not conservative.
- Find the equilibria of the system. Which of the equilibra are stable. Which are not stable? Why?
- Extra credit: Linearize the ODE model around the equilibra and find the eigenvalues to determine stability. Note that if the system is conservative, you will need to reduce the dimension of the system by using the conservation law before you analyze stability.
- Repeat (1) with the model of transcription and translation from the lecture notes (Lec. 4).
- Repeat (1) with the following network
- X → 2X
- X+Y → 2Y
- Y → ∅
- Work through the gro tutorial number 1 (available on the main page of the gro website). Then work out a new example in gro that models both RNA and GFP with the following reactions:
- ∅ → RNA
- RNA → RNA + GFP
- RNA → ∅
- Derive the ODEs that model the rate of change of the concentration of RNA and the concentration of GFP, simulate them, and overlay them with the data obtained from gro, as was done with the simpler system in the tutorial.
- Extra credit: Save the trajectories of all the cells, and make a plot like that in the gro documentation under Tutorial->Molecules or in Fig. 3 of the gro paper. This plot should give you a better idea of the variance in the behaviors you can expect to see.
Turn in A2 here: https://catalyst.uw.edu/collectit/assignment/koishi/23996/95975.