Stochastic and deterministic simulations of heterogeneous cell population dynamics

A Monte Carlo algorithm, which can accurately simulate the dynamics of entire heterogeneous cell populations, was developed. The algorithm takes into account the random nature of cell division as well as unequal partitioning of cellular material at cell division. Moreover, it is general in the sense that it can accommodate a variety of single-cell, deterministic reaction kinetics as well as various stochastic division and partitioning mechanisms. The validity of the algorithm was assessed through comparison of its results with those of the corresponding deterministic cell population balance model in cases where stochastic behavior is expected to be quantitatively negligible. Both algorithms were applied to study: (a) linear intracellular kinetics and (b) the expression dynamics of a genetic network with positive feedback architecture, such as the lac operon. The effects of stochastic division as well as those of different division and partitioning mechanisms were assessed in these systems, while the comparison of the stochastic model with a continuum model elucidated the significance of cell population heterogeneity even in cases where only the prediction of average properties is of primary interest.