Skip to Main content Skip to Navigation
Journal articles

Push-forward method for piecewise deterministic biochemical simulations

Abstract : A biochemical network can be simulated by a set of ordinary differential equations (ODE) under well stirred reactor conditions, for large numbers of molecules, and frequent reactions. This is no longer a robust representation when some molecular species are in small numbers and reactions changing them are infrequent. In this case, discrete stochastic events trigger changes of the smooth deterministic dynamics of the biochemical network. Piecewise-deterministic Markov processes (PDMP) are well adapted for describing such situations. Although PDMP models are now well established in biology, these models remain computationally challenging. Previously we have introduced the push-forward method to compute how the probability measure is spread by the deterministic ODE flow of PDMPs, through the use of analytic expressions of the corresponding semigroup. In this paper we provide a more general simulation algorithm that works also for non-integrable systems. The method can be used for biochemical simulations with applications in fundamental biology, biotechnology and biocomputing.This work is an extended version of the work presented at the conference CMSB2019.
Document type :
Journal articles
Complete list of metadata
Contributor : Arnaud Debussche Connect in order to contact the contributor
Submitted on : Monday, February 15, 2021 - 2:47:56 PM
Last modification on : Friday, August 5, 2022 - 10:56:57 AM
Long-term archiving on: : Sunday, May 16, 2021 - 7:23:30 PM


Files produced by the author(s)



Guilherme C. P. Innocentini, Arran Hodgkinson, Fernando Antoneli, Arnaud Debussche, Ovidiu Radulescu. Push-forward method for piecewise deterministic biochemical simulations. Theoretical Computer Science, Elsevier, 2021, 893, pp.17-40. ⟨10.1016/j.tcs.2021.05.025⟩. ⟨hal-03137189⟩



Record views


Files downloads