Steady-state simulation of metastable stochastic chemical systems

Onderzoeksoutput: ArticleAcademicpeer review

5 Citaten (Scopus)
65 Downloads (Pure)

Samenvatting

We address the problem of steady-state simulation for metastable continuous-time Markov chains with application to stochastic chemical kinetics. Such systems are characterized by the existence of two or more pseudo-equilibrium states and very slow convergence towards global equilibrium. Approximation of the stationary distribution of these systems by direct application of the Stochastic Simulation Algorithm (SSA) is known to be very inefficient. In this paper, we propose a new method for steady-state simulation of metastable Markov chains that is centered around the concept of stochastic complementation. The use of this mathematical device along with SSA results in an algorithm with much better convergence properties, that facilitates the analysis of rarely switching stochastic biochemical systems. The efficiency of our method is demonstrated by its application to two genetic toggle switch models.
Originele taal-2English
Artikelnummer184109
Aantal pagina's11
TijdschriftJournal of Chemical Physics
Volume138
Nummer van het tijdschrift18
DOI's
StatusPublished - 14-mei-2013
Extern gepubliceerdJa

Vingerafdruk

Duik in de onderzoeksthema's van 'Steady-state simulation of metastable stochastic chemical systems'. Samen vormen ze een unieke vingerafdruk.

Citeer dit