{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,3,5]],"date-time":"2026-03-05T19:13:35Z","timestamp":1772738015174,"version":"3.50.1"},"update-to":[{"DOI":"10.1371\/journal.pcbi.1009214","type":"new_version","label":"New version","source":"publisher","updated":{"date-parts":[[2021,8,9]],"date-time":"2021-08-09T00:00:00Z","timestamp":1628467200000}}],"reference-count":49,"publisher":"Public Library of Science (PLoS)","issue":"7","license":[{"start":{"date-parts":[[2021,7,28]],"date-time":"2021-07-28T00:00:00Z","timestamp":1627430400000},"content-version":"vor","delay-in-days":0,"URL":"http:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"name":"institut national de recherche en informatique et en automatique","award":["Project Lab CoSy"],"award-info":[{"award-number":["Project Lab CoSy"]}]},{"name":"institut national de recherche en informatique et en automatique","award":["Project Lab CoSy"],"award-info":[{"award-number":["Project Lab CoSy"]}]},{"name":"institut national de recherche en informatique et en automatique","award":["Project Lab CoSy"],"award-info":[{"award-number":["Project Lab CoSy"]}]},{"name":"ANR","award":["ANR-18-CE91-0002"],"award-info":[{"award-number":["ANR-18-CE91-0002"]}]},{"name":"ANR","award":["ANR-16-CE33-0018"],"award-info":[{"award-number":["ANR-16-CE33-0018"]}]}],"content-domain":{"domain":["www.ploscompbiol.org"],"crossmark-restriction":false},"short-container-title":["PLoS Comput Biol"],"abstract":"<jats:p>The chemical master equation and its continuum approximations are indispensable tools in the modeling of chemical reaction networks. These are routinely used to capture complex nonlinear phenomena such as multimodality as well as transient events such as first-passage times, that accurately characterise a plethora of biological and chemical processes. However, some mechanisms, such as heterogeneous cellular growth or phenotypic selection at the population level, cannot be represented by the master equation and thus have been tackled separately. In this work, we propose a unifying framework that augments the chemical master equation to capture such auxiliary dynamics, and we develop and analyse a numerical solver that accurately simulates the system dynamics. We showcase these contributions by casting a diverse array of examples from the literature within this framework and applying the solver to both match and extend previous studies. Analytical calculations performed for each example validate our numerical results and benchmark the solver implementation.<\/jats:p>","DOI":"10.1371\/journal.pcbi.1009214","type":"journal-article","created":{"date-parts":[[2021,7,28]],"date-time":"2021-07-28T17:35:42Z","timestamp":1627493742000},"page":"e1009214","update-policy":"https:\/\/doi.org\/10.1371\/journal.pcbi.corrections_policy","source":"Crossref","is-referenced-by-count":13,"title":["Beyond the chemical master equation: Stochastic chemical kinetics coupled with auxiliary processes"],"prefix":"10.1371","volume":"17","author":[{"ORCID":"https:\/\/orcid.org\/0000-0002-1603-0840","authenticated-orcid":true,"given":"Davin","family":"Lunz","sequence":"first","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0003-1697-139X","authenticated-orcid":true,"given":"Gregory","family":"Batt","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0003-1615-3282","authenticated-orcid":true,"given":"Jakob","family":"Ruess","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"ORCID":"https:\/\/orcid.org\/0000-0003-0988-9865","authenticated-orcid":true,"given":"J. Fr\u00e9d\u00e9ric","family":"Bonnans","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"340","published-online":{"date-parts":[[2021,7,28]]},"reference":[{"issue":"1","key":"pcbi.1009214.ref001","doi-asserted-by":"crossref","first-page":"404","DOI":"10.1016\/0378-4371(92)90283-V","article-title":"A rigorous derivation of the chemical master equation","volume":"188","author":"DT Gillespie","year":"1992","journal-title":"Physica A"},{"issue":"1","key":"pcbi.1009214.ref002","doi-asserted-by":"crossref","first-page":"297","DOI":"10.1063\/1.481811","article-title":"The chemical Langevin equation","volume":"113","author":"DT Gillespie","year":"2000","journal-title":"J Chem Phys"},{"key":"pcbi.1009214.ref003","volume-title":"Concepts and models of a quantitative sociology: the dynamics of interacting populations","author":"W Weidlich","year":"2012"},{"issue":"6","key":"pcbi.1009214.ref004","doi-asserted-by":"crossref","first-page":"1333","DOI":"10.1111\/j.1365-2958.2009.06605.x","article-title":"A chance at survival: gene expression noise and phenotypic diversification strategies","volume":"71","author":"D Fraser","year":"2009","journal-title":"Mol Microbiol"},{"issue":"7312","key":"pcbi.1009214.ref005","doi-asserted-by":"crossref","first-page":"167","DOI":"10.1038\/nature09326","article-title":"Functional roles for noise in genetic circuits","volume":"467","author":"A Eldar","year":"2010","journal-title":"Nature"},{"issue":"10","key":"pcbi.1009214.ref006","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1371\/journal.pcbi.1008366","article-title":"A Bayesian non-parametric mixed-effects model of microbial growth curves","volume":"16","author":"PD Tonner","year":"2020","journal-title":"PLOS Comput Biol"},{"issue":"1","key":"pcbi.1009214.ref007","doi-asserted-by":"crossref","first-page":"1","DOI":"10.1038\/s41598-018-35927-x","article-title":"Intrinsic and extrinsic noise of gene expression in lineage trees","volume":"9","author":"P Thomas","year":"2019","journal-title":"Sci Rep"},{"issue":"8","key":"pcbi.1009214.ref008","first-page":"1","article-title":"Bistability versus Bimodal Distributions in Gene Regulatory Processes from Population Balance","volume":"7","author":"CC Shu","year":"2011","journal-title":"PLOS Comput Biol"},{"issue":"42","key":"pcbi.1009214.ref009","doi-asserted-by":"crossref","first-page":"13194","DOI":"10.1021\/jp403231f","article-title":"Effect of Phenotypic Selection on Stochastic Gene Expression","volume":"117","author":"T Mora","year":"2013","journal-title":"J Phys Chem B"},{"issue":"1","key":"pcbi.1009214.ref010","doi-asserted-by":"crossref","first-page":"74","DOI":"10.1088\/1478-3975\/3\/1\/008","article-title":"On the distribution of state values of reproducing cells","volume":"3","author":"K Sato","year":"2006","journal-title":"Phys Biol"},{"key":"pcbi.1009214.ref011","article-title":"Stochastic reaction networks in dynamic compartment populations","author":"L Duso","year":"2020","journal-title":"Proc Natl Acad Sci USA"},{"key":"pcbi.1009214.ref012","doi-asserted-by":"crossref","unstructured":"Perthame B. Transport equations in biology. Birkh\u00e4user Basel; 2006.","DOI":"10.1007\/978-3-7643-7842-4"},{"key":"pcbi.1009214.ref013","doi-asserted-by":"crossref","first-page":"168302","DOI":"10.1103\/PhysRevLett.97.168302","article-title":"Linking Stochastic Dynamics to Population Distribution: An Analytical Framework of Gene Expression","volume":"97","author":"N Friedman","year":"2006","journal-title":"Phys Rev Lett"},{"key":"pcbi.1009214.ref014","doi-asserted-by":"crossref","DOI":"10.1214\/EJP.v18-1720","article-title":"Total variation estimates for the TCP process","volume":"18","author":"JB Bardet","year":"2013","journal-title":"Electron J Probab"},{"issue":"4","key":"pcbi.1009214.ref015","doi-asserted-by":"crossref","first-page":"693","DOI":"10.1073\/pnas.1609012114","article-title":"First-passage time approach to controlling noise in the timing of intracellular events","volume":"114","author":"KR Ghusinga","year":"2017","journal-title":"Proc Natl Acad Sci USA"},{"key":"pcbi.1009214.ref016","unstructured":"Lunz D. Flips solver library. Version 0.0.2 [software]; 2021. [Available from]: https:\/\/gitlab.inria.fr\/dlunz\/flips."},{"key":"pcbi.1009214.ref017","volume-title":"Springer Series in Synergetics","author":"C Gardiner","year":"2009","edition":"3"},{"key":"pcbi.1009214.ref018","doi-asserted-by":"crossref","DOI":"10.1007\/978-3-642-61544-3_4","volume-title":"Fokker\u2013Planck Equation","author":"H Risken","year":"1996"},{"key":"pcbi.1009214.ref019","doi-asserted-by":"crossref","first-page":"3","DOI":"10.2307\/3318611","article-title":"Deterministic and Stochastic Models for Coalescence (Aggregation, Coagulation): a Review of the Mean-Field Theory for Probabilists","volume":"5","author":"DJ Aldous","year":"1997","journal-title":"Bernoulli"},{"issue":"12","key":"pcbi.1009214.ref020","doi-asserted-by":"crossref","first-page":"1059","DOI":"10.1038\/nchembio.2205","article-title":"A photoactivatable Cre\u2013loxP recombination system for optogenetic genome engineering","volume":"12","author":"F Kawano","year":"2016","journal-title":"Nat Chem Biol"},{"issue":"4","key":"pcbi.1009214.ref021","doi-asserted-by":"crossref","first-page":"403","DOI":"10.1016\/0021-9991(76)90041-3","article-title":"A general method for numerically simulating the stochastic time evolution of coupled chemical reactions","volume":"22","author":"DT Gillespie","year":"1976","journal-title":"J Comput Phys"},{"issue":"25","key":"pcbi.1009214.ref022","doi-asserted-by":"crossref","first-page":"2340","DOI":"10.1021\/j100540a008","article-title":"Exact stochastic simulation of coupled chemical reactions","volume":"81","author":"DT Gillespie","year":"1977","journal-title":"J Phys Chem"},{"issue":"9","key":"pcbi.1009214.ref023","doi-asserted-by":"crossref","first-page":"1876","DOI":"10.1021\/jp993732q","article-title":"Efficient Exact Stochastic Simulation of Chemical Systems with Many Species and Many Channels","volume":"104","author":"MA Gibson","year":"2000","journal-title":"J Phys Chem A"},{"issue":"4","key":"pcbi.1009214.ref024","doi-asserted-by":"crossref","first-page":"1716","DOI":"10.1063\/1.1378322","article-title":"Approximate accelerated stochastic simulation of chemically reacting systems","volume":"115","author":"DT Gillespie","year":"2001","journal-title":"J Chem Phys"},{"issue":"1","key":"pcbi.1009214.ref025","doi-asserted-by":"crossref","first-page":"100","DOI":"10.1016\/j.jcp.2007.07.020","article-title":"Hybrid method for the chemical master equation","volume":"227","author":"A Hellander","year":"2007","journal-title":"J Comput Phys"},{"issue":"1","key":"pcbi.1009214.ref026","doi-asserted-by":"crossref","first-page":"014116","DOI":"10.1063\/1.1824902","article-title":"The slow-scale stochastic simulation algorithm","volume":"122","author":"Y Cao","year":"2005","journal-title":"J Chem Phys"},{"issue":"4","key":"pcbi.1009214.ref027","doi-asserted-by":"crossref","first-page":"044104","DOI":"10.1063\/1.2145882","article-title":"The finite state projection algorithm for the solution of the chemical master equation","volume":"124","author":"B Munsky","year":"2006","journal-title":"J Chem Phys"},{"issue":"20","key":"pcbi.1009214.ref028","doi-asserted-by":"crossref","first-page":"204104","DOI":"10.1063\/1.2397685","article-title":"Reduction and solution of the chemical master equation using time scale separation and finite state projection","volume":"125","author":"S Pele\u0161","year":"2006","journal-title":"J Chem Phys"},{"issue":"13","key":"pcbi.1009214.ref029","doi-asserted-by":"crossref","first-page":"2018","DOI":"10.1021\/es970967b","article-title":"Condensing Complex Atmospheric Chemistry Mechanisms. 1. The Direct Constrained Approximate Lumping (DCAL) Method Applied to Alkane Photochemistry","volume":"32","author":"SW Wang","year":"1998","journal-title":"Environ Sci Technol"},{"issue":"1","key":"pcbi.1009214.ref030","doi-asserted-by":"crossref","first-page":"208","DOI":"10.1016\/j.cam.2008.10.029","article-title":"Spectral approximation of solutions to the chemical master equation","volume":"229","author":"S Engblom","year":"2009","journal-title":"J Comput Appl Math"},{"issue":"16","key":"pcbi.1009214.ref031","doi-asserted-by":"crossref","first-page":"5724","DOI":"10.1016\/j.jcp.2010.04.015","article-title":"Solving chemical master equations by adaptive wavelet compression","volume":"229","author":"T Jahnke","year":"2010","journal-title":"J Comput Phys"},{"issue":"3","key":"pcbi.1009214.ref032","doi-asserted-by":"crossref","first-page":"034118","DOI":"10.1063\/1.4905196","article-title":"Adaptive hybrid simulations for multiscale stochastic reaction networks","volume":"142","author":"B Hepp","year":"2015","journal-title":"J Chem Phys"},{"issue":"1","key":"pcbi.1009214.ref033","doi-asserted-by":"crossref","first-page":"37","DOI":"10.1007\/s00791-006-0045-6","article-title":"Fokker\u2013Planck approximation of the master equation in molecular biology","volume":"12","author":"P Sj\u00f6berg","year":"2009","journal-title":"Comput Vis Sci"},{"issue":"1","key":"pcbi.1009214.ref034","doi-asserted-by":"crossref","first-page":"61","DOI":"10.1007\/s10543-006-0082-z","article-title":"Conservative solution of the Fokker\u2013Planck equation for stochastic chemical reactions","volume":"46","author":"L Ferm","year":"2006","journal-title":"BIT"},{"issue":"15","key":"pcbi.1009214.ref035","doi-asserted-by":"crossref","first-page":"6959","DOI":"10.1063\/1.1505860","article-title":"Approximate simulation of coupled fast and slow reactions for stochastic chemical kinetics","volume":"117","author":"EL Haseltine","year":"2002","journal-title":"J Chem Phys"},{"issue":"5","key":"pcbi.1009214.ref036","doi-asserted-by":"crossref","first-page":"054103","DOI":"10.1063\/1.1835951","article-title":"Accurate hybrid stochastic simulation of a system of coupled chemical or biochemical reactions","volume":"122","author":"H Salis","year":"2005","journal-title":"J Chem Phys"},{"key":"pcbi.1009214.ref037","doi-asserted-by":"crossref","first-page":"177","DOI":"10.1016\/j.jcp.2014.10.026","article-title":"Hybrid discrete\/continuum algorithms for stochastic reaction networks","volume":"281","author":"C Safta","year":"2015","journal-title":"J Comput Phys"},{"key":"pcbi.1009214.ref038","first-page":"2007","volume-title":"PDE and Monte Carlo Approaches to Solving the Master Equation Applied to Gene Regulation","author":"P Sj\u00f6berg","year":"2007"},{"issue":"1","key":"pcbi.1009214.ref039","doi-asserted-by":"crossref","first-page":"241","DOI":"10.1006\/jcph.2000.6459","article-title":"New high-resolution central schemes for nonlinear conservation laws and convection\u2013diffusion equations","volume":"160","author":"A Kurganov","year":"2000","journal-title":"J Comput Phys"},{"key":"pcbi.1009214.ref040","doi-asserted-by":"crossref","first-page":"022409","DOI":"10.1103\/PhysRevE.93.022409","article-title":"Gene expression dynamics with stochastic bursts: Construction and exact results for a coarse-grained model","volume":"93","author":"YT Lin","year":"2016","journal-title":"Phys Rev E"},{"issue":"2","key":"pcbi.1009214.ref041","doi-asserted-by":"crossref","first-page":"367","DOI":"10.1073\/pnas.1808080116","article-title":"Predicting growth rate from gene expression","volume":"116","author":"TP Wytock","year":"2019","journal-title":"Proc Natl Acad Sci USA"},{"issue":"1","key":"pcbi.1009214.ref042","doi-asserted-by":"crossref","first-page":"294","DOI":"10.1137\/20M1332293","article-title":"On Continuum Approximations of Discrete-State Markov Processes of Large System Size","volume":"19","author":"D Lunz","year":"2021","journal-title":"Multiscale Modeling & Simulation"},{"issue":"6","key":"pcbi.1009214.ref043","doi-asserted-by":"crossref","first-page":"972","DOI":"10.1016\/j.cub.2018.02.006","article-title":"Size Laws and Division Ring Dynamics in Filamentous Escherichia coli cells","volume":"28","author":"M Wehrens","year":"2018","journal-title":"Current Biology"},{"issue":"1","key":"pcbi.1009214.ref044","doi-asserted-by":"crossref","first-page":"119","DOI":"10.1007\/s10665-018-9968-4","article-title":"Dynamics of particle chopping in blenders and food processors","volume":"112","author":"CM Rooney","year":"2018","journal-title":"J Eng Math"},{"issue":"4","key":"pcbi.1009214.ref045","doi-asserted-by":"crossref","first-page":"334","DOI":"10.1016\/j.matpur.2011.01.003","article-title":"Rate of convergence to an asymptotic profile for the self-similar fragmentation and growth-fragmentation equations","volume":"96","author":"MJ C\u00e1ceres","year":"2011","journal-title":"J Math Pure Appl"},{"issue":"5","key":"pcbi.1009214.ref046","doi-asserted-by":"crossref","first-page":"401","DOI":"10.1038\/sj.gt.3301102","article-title":"Cell cycle dependence of gene transfer by lipoplex, polyplex and recombinant adenovirus","volume":"7","author":"S Brunner","year":"2000","journal-title":"Gene ther"},{"issue":"1","key":"pcbi.1009214.ref047","doi-asserted-by":"crossref","first-page":"53","DOI":"10.1016\/S0167-4781(99)00039-1","article-title":"Mitosis enhances transgene expression of plasmid delivered by cationic liposomes","volume":"1445","author":"WC Tseng","year":"1999","journal-title":"BBA-Gene Struct Expr"},{"issue":"7","key":"pcbi.1009214.ref048","first-page":"1","article-title":"Cell-Cycle Dependence of Transcription Dominates Noise in Gene Expression","volume":"9","author":"CJ Zopf","year":"2013","journal-title":"PLOS Comput Biol"},{"issue":"34","key":"pcbi.1009214.ref049","doi-asserted-by":"crossref","first-page":"8505","DOI":"10.1073\/pnas.1718942115","article-title":"Solving high-dimensional partial differential equations using deep learning","volume":"115","author":"J Han","year":"2018","journal-title":"Proc Natl Acad Sci USA"}],"updated-by":[{"DOI":"10.1371\/journal.pcbi.1009214","type":"new_version","label":"New version","source":"publisher","updated":{"date-parts":[[2021,8,9]],"date-time":"2021-08-09T00:00:00Z","timestamp":1628467200000}}],"container-title":["PLOS Computational Biology"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/dx.plos.org\/10.1371\/journal.pcbi.1009214","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2021,8,9]],"date-time":"2021-08-09T19:08:02Z","timestamp":1628536082000},"score":1,"resource":{"primary":{"URL":"https:\/\/dx.plos.org\/10.1371\/journal.pcbi.1009214"}},"subtitle":[],"editor":[{"given":"Vassily","family":"Hatzimanikatis","sequence":"first","affiliation":[],"role":[{"role":"editor","vocabulary":"crossref"}]}],"short-title":[],"issued":{"date-parts":[[2021,7,28]]},"references-count":49,"journal-issue":{"issue":"7","published-online":{"date-parts":[[2021,7,28]]}},"URL":"https:\/\/doi.org\/10.1371\/journal.pcbi.1009214","relation":{"new_version":[{"id-type":"doi","id":"10.1371\/journal.pcbi.1009214","asserted-by":"object"}]},"ISSN":["1553-7358"],"issn-type":[{"value":"1553-7358","type":"electronic"}],"subject":[],"published":{"date-parts":[[2021,7,28]]}}}