Constructing stochastic models from deterministic process equations by propensity adjustment

Abstract

BackgroundGillespie's stochastic simulation algorithm (SSA) for chemical reactions admits three kinds of elementary processes, namely, mass action reactions of 0th, 1st or 2nd order. All other types of reaction processes, for instance those containing non-integer kinetic orders or following other types of kinetic laws, are assumed to be convertible to one of the three elementary kinds, so that SSA can validly be applied. However, the conversion to elementary reactions is often difficult, if not impossible. Within deterministic contexts, a strategy of model reduction is often used. Such a reduction simplifies the actual system of reactions by merging or approximating intermediate steps and omitting reactants such as transient complexes. It would be valuable to adopt a similar reduction strategy to stochastic modelling. Indeed, efforts have been devoted to manipulating the chemical master equation (CME) in order to achieve a proper propensity function for a reduced stochastic system. However, manipulations of CME are almost always complicated, and successes have been limited to relative simple cases.ResultsWe propose a rather general strategy for converting a deterministic process model into a corresponding stochastic model and characterize the mathematical connections between the two. The deterministic framework is assumed to be a generalized mass action system and the stochastic analogue is in the format of the chemical master equation. The analysis identifies situations: where a direct conversion is valid; where internal noise affecting the system needs to be taken into account; and where the propensity function must be mathematically adjusted. The conversion from deterministic to stochastic models is illustrated with several representative examples, including reversible reactions with feedback controls, Michaelis-Menten enzyme kinetics, a genetic regulatory motif, and stochastic focusing.ConclusionsThe construction of a stochastic model for a biochemical network requires the utilization of information associated with an equation-based model. The conversion strategy proposed here guides a model design process that ensures a valid transition between deterministic and stochastic models.

Biochemical Systems Theory: A Review

“…Bajzer's [674][675][676] group tested fractal power-law representations with carefully designed experiments and found them to be valid and more appropriate than models with time-dependent reaction rates. Wu et al analyzed the same data with stochastic equivalents of GMA systems [176].…”

“…Wu and Voit embedded BST into hybrid functional Petri nets, which expanded the range of phenomena that could be modeled [172][173][174] (see also [175]). ey also considered stochastic analogues to GMA systems [176]; see also [177].…”

Global sensitivity analysis in stochastic simulators of uncertain reaction networks

“…Examples where these features would be of special importance include the analysis of system robustness, model reduction, and parameter calibration. 4,[34][35][36][37][38] The application of the proposed method to complex stochastic systems remains challenging because of the sampling effort. Even though the method is trivial to implement in parallel, future works should focus on the improvement of the methods to lower its computational complexity.…”

Optic tectum morphogenesis: A step‐by‐step model based on the temporal–spatial organization of the cell proliferation. Significance of deterministic and stochastic components subsumed in the spatial organization

“…This study also presents a biomathematical approach that could provide a conceptual framework to analyze a recurrent and controversial point in developmental biology: the problem about ''determinism'' and ''stochasticity'' and their possible role in developing organisms (Bronk et al, 1968;Saunders et al, 2002;Meng et al, 2004;Song et al, 2006;Raj et al, 2010;Zernicka-Goetz and Huang, 2010;Wu et al, 2011).…”