Perspective: Stochastic algorithms for chemical kinetics

Abstract

We outline our perspective on stochastic chemical kinetics, paying particular attention to numerical simulation algorithms. We first focus on dilute, well-mixed systems, whose description using ordinary differential equations has served as the basis for traditional chemical kinetics for the past 150 years. For such systems, we review the physical and mathematical rationale for a discretestochastic approach, and for the approximations that need to be made in order to regain the traditional continuous-deterministic description. We next take note of some of the more promising strategies for dealing stochastically with stiff systems, rare events, and sensitivity analysis. Finally, we review some recent efforts to adapt and extend the discrete-stochastic approach to systems that are not wellmixed. In that currently developing area, we focus mainly on the strategy of subdividing the system into well-mixed subvolumes, and then simulating diffusional transfers of reactant molecules between adjacent subvolumes together with chemical reactions inside the subvolumes. © 2013 AIP Publishing LLC.[http://dx

Relative frequencies of constrained events in stochastic processes: An analytical approach

“…The assumption is that the events belonging to each subset (j 1 , j 2 , .., j N ) are independent and identically distributed (i.i.d.). Such an assumption does not lead to loss of generality because the constraints used in building different subsets can contain the information about the dependencies between events [13]. Then if n(j 1 , j 2 , .., j N ) → +∞, the probability P i (j 1 , j 2 , .., j N ) for the i-th outcome to occur in the set (j 1 , j 2 , .., j N ) corresponds to the limit of frequency of the i-th outcome, i.e.…”

Lazy Updating of hubs can enable more realistic models by speeding up stochastic simulations

“…For example, reaction volume affects many propensities. 3,16,21 If the volume is affected in small steps by many actions in the system, then its impact on propensities is similar to that of hubs: frequent small changes in volume trigger frequent small updates to many propensities. Biological examples for such endogenously driven changes include active cellular nutrient import, 44 osmosis, 44 and export of molecules.…”

“…Continuous-time Markov chains [10][11][12][13][14][15] (CTMCs) have a track record of excellence in this area and provide a rigorous basis for transforming certain types of well-mixed stochastic models with infinitely large part numbers into equivalent ODE models. 3,14,15 CTMCs have no "memory," as they assume that the future state of a model only depends on its current state and not on its history. CTMCs also assume that the waiting times between actions (which lead to state transitions) are exponentially distributed, unless external noise affects the probabilities that determine how often these actions occur.…”

“…However, many processes in biology are discrete and involve very few key parts, such as switching on single-copy genes or single viral particles that cause dramatic growth but only if they infect a cell. Stochastic effects are important in these systems [3][4][5] and simulation models should respect the discreteness of molecules, individuals and unique events. Continuous-time Markov chains [10][11][12][13][14][15] (CTMCs) have a track record of excellence in this area and provide a rigorous basis for transforming certain types of well-mixed stochastic models with infinitely large part numbers into equivalent ODE models.…”

“…Time series can help by showing how the amounts of important entities in a network change over time. They are used in many areas of biology, including biochemistry, [1][2][3][4][5] ecology, evolution, 6,7 and systems biology. 1,2 Since simulation techniques for predicting time series are of general interest, we will discuss models in terms of "parts" and "actions" instead of "molecules" and "reactions."…”