Simulation results of the community effect model. (A, B) Numerical simulations of the deterministic rate equations Eqs.7-14. (A) is the plot for [Ap] and (B) for [Bp
]. Simulation results for different community size are shown. With 100 cells, very little gene expression occurs at steady state (not shown) as 100 cells are close to n
≈ 97. (C) Average number of Ap at steady state (10000 min) as a function of community size (solid lines). Dotted curves indicate number of Ap at steady state ([Ap]*) obtained by Eqs.35 in additional file 1. Plots are shown for ε = 2.31 × 10-6 and 5.78 × 10-7. (D, E) Time series of [Ap] and [Bp
] as overlays of 100 stochastic simulation results (temperature map) for community size of 300 cells and ε = 5.78 × 10 -7. Solid lines show a typical simulation result. (F) Probability distributions of [Ap] at steady state (t = 10000 min in stochastic simulations) for community size n = 140, 200, and 500 cells. ε = 5.78 × 10 -7. All simulations in this figure are with one gene copy each for gene A and gene B.