A solver for the stochastic master equation applied to gene regulatory networks

Markus Hegland*, Conrad Burden, Lucia Santoso, Shev MacNamara, Hilary Booth

*Corresponding author for this work

    Research output: Contribution to journalArticlepeer-review

    84 Citations (Scopus)

    Abstract

    An important driver of gene regulatory networks is noise arising from the stochastic nature of interactions of genes, their products and regulators. Thus, such systems are stochastic and can be modelled by the chemical master equations. A major challenge is the curse of dimensionality which occurs when one attempts to integrate these equations. While stochastic simulation techniques effectively address the curse, many repeated simulations are required to provide precise information about stationary points, bifurcation phenomena and other properties of the stochastic processes. An alternative way to address the curse of dimensionality is provided by sparse grid approximations. The sparse grid methodology is applied and the application demonstrated to work efficiently for up to 10 proteins. As sparse grid methods have been developed for the approximation of smooth functions, a variant for infinite sequences had to be developed together with a multiresolution analysis similar to Haar wavelets. Error bounds are provided which confirm the effectiveness of sparse grid approximations for smooth high-dimensional probability distributions.

    Original languageEnglish
    Pages (from-to)708-724
    Number of pages17
    JournalJournal of Computational and Applied Mathematics
    Volume205
    Issue number2
    DOIs
    Publication statusPublished - 15 Aug 2007

    Fingerprint

    Dive into the research topics of 'A solver for the stochastic master equation applied to gene regulatory networks'. Together they form a unique fingerprint.

    Cite this