Stochastic focusing. The first panel from the top compares the time evolution of product molecules P obtained with different methods: the black line represents the solution of ODE model (64) for P; the blue solid line and blue dashed lines are the solutions of the moment-based model (65) for μP and μP ± σP, respectively. The second panel indicates that the stochastic simulations without propensity adjustment (blue error bar) diverge from the prediction by the ODE model (64) (black line). The third panel shows the propensity adjustment function paf4(t) = exp(-fI fS cov [log I(t), log S(t)]) for the reaction . The bottom panel demonstrates that the propensity adjustment function paf4 achieves convergence between the stochastic simulation and the ODE model (64) (black line): the blue error bars were computed from 100 independent stochastic simulations with propensity adjustment paf4. The simulation parameters are (i(0), p(0), s(0), k1, k2, k3, k4, k5, k6, fI, fS) = (0, 10, 100, 104, 103, 1, 9.9 × 103, 104, 103, 1.1, 0.9); at t = 0.1, the value of k5 changes from 104 to 0.5 ×104.
Wu et al. BMC Systems Biology 2011 5:187 doi:10.1186/1752-0509-5-187