Block-tridiagonal state-space realization of Chemical Master Equations: A tool to compute explicit solutions

JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS

卷:296

Article

Borri, A.¹  Carravetta, F.¹  Mavelli, G.¹  Palumbo, P.^1,2 

[1] Consiglio Nazl Ric IASI CNR, Ist Anal Sistemi & Informat A Ruberti, Rome, Italy

[2] SYSBIO Ctr Syst Biol, Milan, Italy

关键词: Chemical Master Equation;    Markov processes;    Systems biology;   

DOI  :  10.1016/j.cam.2015.10.008

来源: Elsevier

【 摘 要 】

Chemical Master Equations (CMEs) provide a comprehensive way to model the probabilistic behavior in biochemical networks. Despite their widespread diffusion in systems biology, the explicit computation of their solution is often avoided in favor of purely statistic Monte Carlo methods, due to the dramatically high dimension of the CME system. In this work, we investigate some structural properties of CMEs and their solutions, focusing on the efficient computation of the stationary distribution. We introduce a generalized notion of one-step process, which results in a sparse dynamic matrix describing the collection of the scalar CMEs, showing a recursive block-tridiagonal structure as well. Further properties are inferred by means of a graph-theoretical interpretation of the reaction network. We exploit this structure by proposing different methods, including a dedicated LU decomposition, to compute the explicit solution. Examples are included to illustrate the introduced concepts and to show the effectiveness of the proposed approach. (C) 2015 Elsevier B.V. All rights reserved.

【 授权许可】

Free