Figure 4.

Noise statistics for cooperative two-subunit enzyme network. Plots showing the Fano factor multiplied by the volume, ΩFFS, (a) and the coefficient of variation squared, <a onClick="popup('http://www.biomedcentral.com/1752-0509/6/39/mathml/M88','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/39/mathml/M88">View MathML</a>, (b) for the substrate fluctuations as a function of the non-dimensional fraction Θ, in steady-state conditions. The latter is a measure of enzyme saturation. The solid lines are the ssLNA predictions, the dashed lines are the hLNA predictions, the solid circles are obtained from stochastic simulations of the full network and the open circles are obtained from stochastic simulations of the coarse-grained network using the SSA with heuristic propensities. The color coding indicates different values of the bimolecular rate constant k1: 5×10−3(yellow), 5×10−5(purple), and 5×10−7(blue). The remaining parameters are given by <a onClick="popup('http://www.biomedcentral.com/1752-0509/6/39/mathml/M89','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/39/mathml/M89">View MathML</a>, k2=1000, k−1=k−2=100, k3=k4=1. Note that in (a) the black dashed line indicates the hLNA prediction for all three different values of k1, which are indistinguishable in this figure. The stochastic simulations were carried out for a volume Ω=100. In (c) sample paths of the SSA for the full network (gray), the slow scale Langevin equation (red) as given by equation (18) and the SSA with heuristic propensities (blue) are compared for Θ=0.5. The slow scale Langevin equation is numerically solved using the Euler-Maruyama method with timestep δt=0.1. Note that in all cases, the chosen parameters guarantee timescale separation (validity of the deterministic QSSA) (see Figure 5)

Thomas et al. BMC Systems Biology 2012 6:39   doi:10.1186/1752-0509-6-39
Download authors' original image