An Algorithm for the Stochastic Simulation of Gene Expression and Heterogeneous Population Dynamics

We present an algorithm for the stochastic simulation of gene expression and heterogeneous population dynamics. The algorithm combines an exact method to simulate molecular-level fluctuations in single cells and a constant-number Monte Carlo method to simulate time-dependent statistical characteristics of growing cell populations. To benchmark performance, we compare simulation results with steady-state and time-dependent analytical solutions for several scenarios, including steady-state and time-dependent gene expression, and the effects on population heterogeneity of cell growth, division, and DNA replication. This comparison demonstrates that the algorithm provides an efficient and accurate approach to simulate how complex biological features influence gene expression. We also use the algorithm to model gene expression dynamics within “bet-hedging” cell populations during their adaption to environmental stress. These simulations indicate that the algorithm provides a framework suitable for simulating and analyzing realistic models of heterogeneous population dynamics combining molecular-level stochastic reaction kinetics, relevant physiological details and phenotypic variability.

[1]  Daniel T Gillespie,et al.  Stochastic simulation of chemical kinetics. , 2007, Annual review of physical chemistry.

[2]  Themis Matsoukas,et al.  Constant-number Monte Carlo simulation of population balances , 1998 .

[3]  S. Stearns.  Life-History Tactics: A Review of the Ideas , 1976, The Quarterly Review of Biology.

[4]  W. Schaffer.  Optimal Reproductive Effort in Fluctuating Environments , 1974, The American Naturalist.

[5]  Vahid Shahrezaei,et al.  Analytical distributions for stochastic gene expression , 2008, Proceedings of the National Academy of Sciences.

[6]  D. Cohen.  Optimizing reproduction in a randomly varying environment. , 1966, Journal of theoretical biology.

[7]  D. Volfson,et al.  Origins of extrinsic variability in eukaryotic gene expression , 2006, Nature.

[8]  Uri Alon,et al.  A fluctuation method to quantify in vivo fluorescence data. , 2006, Biophysical journal.

[9]  D. Gillespie.  Exact Stochastic Simulation of Coupled Chemical Reactions , 1977 .

[10]  Nikos V Mantzaris,et al.  From single-cell genetic architecture to cell population dynamics: quantitatively decomposing the effects of different population heterogeneity sources for a genetic network with positive feedback architecture. , 2007, Biophysical journal.

[11]  D. Gillespie.  Approximate accelerated stochastic simulation of chemically reacting systems , 2001 .

[12]  Kangtaek Lee,et al.  Simultaneous coagulation and break-up using constant-N Monte Carlo , 2000 .

[13]  P. Swain,et al.  Intrinsic and extrinsic contributions to stochasticity in gene expression , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[14]  David McMillen,et al.  Biochemical Network Stochastic Simulator (BioNetS): software for stochastic modeling of biochemical networks , 2004, BMC Bioinformatics.

[15]  J Hasty,et al.  Cellular growth and division in the Gillespie algorithm. , 2004, Systems biology.

[16]  A. Arkin,et al.  From Fluctuations to Phenotypes: The Physiology of Noise , 2006, Science's STKE.

[17]  Kangtaek Lee,et al.  Solution of the population balance equation using constant-number Monte Carlo , 2002 .

[18]  Andre S. Ribeiro,et al.  CellLine, a stochastic cell lineage simulator , 2007, Bioinform..

[19]  Alexander van Oudenaarden,et al.  Stochastic Gene Expression: from Single Molecules to the Proteome This Review Comes from a Themed Issue on Chromosomes and Expression Mechanisms Edited Measuring Noise Mrna Fluctuations , 2022 .

[20]  D. Gillespie.  A General Method for Numerically Simulating the Stochastic Time Evolution of Coupled Chemical Reactions , 1976 .

[21]  Jerome T. Mettetal,et al.  Stochastic switching as a survival strategy in fluctuating environments , 2008, Nature Genetics.

[22]  Edward D. Lazowska,et al.  Speedup Versus Efficiency in Parallel Systems , 1989, IEEE Trans. Computers.

[23]  J. Tyson,et al.  Sloppy size control of the cell division cycle. , 1986, Journal of theoretical biology.

[24]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[25]  Mads Kærn,et al.  A chance at survival: gene expression noise and phenotypic diversification strategies , 2009, Molecular microbiology.

[26]  Andrzej M. Kierzek,et al.  STOCKS: STOChastic Kinetic Simulations of biochemical systems with Gillespie algorithm , 2002, Bioinform..

[27]  Margaritis Kostoglou,et al.  Evaluation of Zero Order Methods for Simulating Particle Coagulation , 1994 .

[28]  T. Elston,et al.  Stochasticity in gene expression: from theories to phenotypes , 2005, Nature Reviews Genetics.

[29]  W. L. Fangman,et al.  Cell cycle phases in the unequal mother/daughter cell cycles of Saccharomyces cerevisiae , 1984, Molecular and cellular biology.

[30]  T. Kepler,et al.  Stochasticity in transcriptional regulation: origins, consequences, and mathematical representations. , 2001, Biophysical journal.

[31]  Rui Zhu,et al.  Validation of an algorithm for delay stochastic simulation of transcription and translation in prokaryotic gene expression , 2006, Physical biology.

[32]  J. Paulsson.  Summing up the noise in gene networks , 2004, Nature.

[33]  E. O’Shea,et al.  Living with noisy genes: how cells function reliably with inherent variability in gene expression. , 2007, Annual review of biophysics and biomolecular structure.

[34]  D. Ramkrishna.  The Status of Population Balances , 1985 .

[35]  R. Levins.  Evolution in Changing Environments: Some Theoretical Explorations. (MPB-2) , 1968 .

[36]  Nikos V Mantzaris,et al.  Stochastic and deterministic simulations of heterogeneous cell population dynamics. , 2006, Journal of theoretical biology.