Abstract
Gene regulatory networks are highly complex dynamical systems comprising biomolecular components which interact with each other and through those interactions determine gene expression levels, that is, determine the rate of gene transcription. In this paper, a particle filter with Markov Chain Monte Carlo move step is employed for the estimation of reaction rate constants in gene regulatory networks modeled by chemical Langevin equations. Simulation studies demonstrate that the proposed technique outperforms previously considered methods while being computationally more efficient. Dynamic behavior of gene regulatory networks averaged over a large number of cells can be modeled by ordinary differential equations. For this scenario, we compute an approximation to the Cramer-Rao lower bound on the mean-square error of estimating reaction rates and demonstrate that, when the number of unknown parameters is small, the proposed particle filter can be nearly optimal.
Highlights
Gene regulatory networks (GRN) are systems comprising biomolecular components that interact with each other and through those interactions determine gene expression levels, that is, determine the rate of gene transcription to mRNA [1,2,3]
We studied the problem of estimating reaction rates in a gene regulatory network modeled by a chemical Langevin equation, that is, a high-dimensional stochastic differential equation
We proposed a solution which employs a particle filtering algorithm with Markov Chain Monte Carlo move step
Summary
Gene regulatory networks (GRN) are systems comprising biomolecular components (genes, mRNA, proteins) that interact with each other and through those interactions determine gene expression levels, that is, determine the rate of gene transcription to mRNA [1,2,3]. Dynamic behavior of gene regulatory networks averaged over a large number of cells can be modeled by ordinary differential equations For this scenario, we compute an approximation to the CramerRao lower bound on the mean-square error of estimating reaction rates and demonstrate that, when the number of unknown parameters is small, the proposed particle filter can be nearly optimal. The function hm(X(t)) counts all possible combinations of individual molecules that may lead to a reaction in the mth channel.) The chemical master equation is often used to simulate the Markov process X(t) and enable computational studies of GRN To this end, one may employ various stochastic simulation algorithms, originally proposed by Gillespie [4].
Sign-in/Register to view highlights & summary Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callDisclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
