title: Numerical solution of the chemical master equation uniqueness and stability of the stationary distribution for chemical networks, and mRNA bursting in a gene network with negative feedback regulation.
abstract: In this work, we introduce a couple of algorithms to compute the stationary probability distribution for the chemical master equation (CME) of arbitrary chemical networks. We further find the conditions guaranteeing the algorithms' convergence and the unity and stability of the stationary distribution. Next, we employ these algorithms to study the mRNA and protein probability distributions in a gene regulatory network subject to negative feedback regulation. In particular, we analyze the influence of the promoter activation/deactivation speed on the shape of such distributions. We find that a reduction of the promoter activation/deactivation speed modifies the shape of those distributions in a way consistent with the phenomenon known as mRNA (or transcription) bursting.
authors.affiliation: Centro de Investigación y de Estudios Avanzados del IPN, Departamento de Matemáticas, Av. Instituto Politécnico Nacional 2508, México DF, México.
authors.forename: E S
authors.initials: ES
authors.lastname: Zeron
authors.forename: M
authors.initials: M
authors.lastname: Santillán
mesh_terms: D000465:Algorithms; D053263:Gene Regulatory Networks; D008957:Models, Genetic; D011336:Probability; D012333:RNA, Messenger / Q000737:chemistry*
chemical_list: D012333:RNA, Messenger