An application of the Krylov-FSP-SSA method to parameter fitting with maximum likelihood

Physical Biology
Khanh N Dinh, Roger B Sidje

Abstract

Monte Carlo methods such as the stochastic simulation algorithm (SSA) have traditionally been employed in gene regulation problems. However, there has been increasing interest to directly obtain the probability distribution of the molecules involved by solving the chemical master equation (CME). This requires addressing the curse of dimensionality that is inherent in most gene regulation problems. The finite state projection (FSP) seeks to address the challenge and there have been variants that further reduce the size of the projection or that accelerate the resulting matrix exponential. The Krylov-FSP-SSA variant has proved numerically efficient by combining, on one hand, the SSA to adaptively drive the FSP, and on the other hand, adaptive Krylov techniques to evaluate the matrix exponential. Here we apply this Krylov-FSP-SSA to a mutual inhibitory gene network synthetically engineered in Saccharomyces cerevisiae, in which bimodality arises. We show numerically that the approach can efficiently approximate the transient probability distribution, and this has important implications for parameter fitting, where the CME has to be solved for many different parameter sets. The fitting scheme amounts to an optimization problem of fi...Continue Reading

References

Mar 31, 1999·Trends in Genetics : TIG·H H McAdams, Adam P Arkin
Feb 5, 2000·Nature·M B Elowitz, S Leibler
Feb 5, 2000·Nature·T S GardnerJ J Collins
Aug 31, 2004·The Journal of Chemical Physics·Yang CaoLinda Petzold
Jan 11, 2005·The Journal of Chemical Physics·Abhijit ChatterjeeMarkos A Katsoulakis
Jan 11, 2005·The Journal of Chemical Physics·Yang CaoLinda Petzold
Feb 11, 2005·Proceedings of the National Academy of Sciences of the United States of America·Michael SamoilovAdam P Arkin
Aug 20, 2005·The Journal of Chemical Physics·Yang CaoLinda Petzold
Sep 20, 2005·Nature·Alejandro Colman-LernerRoger Brent
Dec 3, 2005·Computational Biology and Chemistry·James M McCollumNagiza F Samatova
Dec 20, 2005·Cell·Ido GoldingEdward C Cox
Feb 8, 2006·The Journal of Chemical Physics·Brian Munsky, Mustafa Khammash
Feb 8, 2006·The Journal of Chemical Physics·Yang CaoLinda Petzold
May 19, 2006·Journal of Computational Biology : a Journal of Computational Molecular Cell Biology·Andrew Golightly, Darren J Wilkinson
Jul 19, 2006·PLoS Computational Biology·David Fange, Johan Elf
Sep 22, 2006·Systems Biology·S ReinkerJ Timmer
Dec 6, 2006·The Journal of Chemical Physics·Slaven PelesMustafa Khammash
May 13, 1983·Science·S KirkpatrickM P Vecchi
Jun 3, 2008·The Journal of Chemical Physics·Alexander SlepoySteven J Plimpton
Jan 14, 2009·Nature Reviews. Genetics·Darren J Wilkinson
Jul 2, 2009·The Journal of Chemical Physics·Rajesh RamaswamyIvo F Sbalzarini
Apr 10, 2010·BMC Systems Biology·Verena WolfThomas A Henzinger
Aug 10, 2010·BMC Bioinformatics·Suresh Kumar Poovathingal, Rudiyanto Gunawan
Nov 16, 2010·IET Systems Biology·B Munsky, Mustafa Khammash
May 3, 2012·BMC Bioinformatics·Bernie J DaigleJarad Niemi
Feb 2, 2013·Science·Gregor NeuertAlexander van Oudenaarden
Jun 12, 2013·Proceedings of the National Academy of Sciences of the United States of America·Min WuXiao Wang
Aug 22, 2016·The Journal of Chemical Physics·Zachary FoxBrian Munsky
Sep 27, 2016·Physical Review Letters·Heng XuIdo Golding
Aug 3, 2017·The Journal of Chemical Physics·Huy D Vo, Roger B Sidje

Related Concepts

Probability
Saccharomyces cerevisiae
Systems Biology
Gene Modules
Synthetic Biology
Extracellular Matrix
Fibrin split products
Chemicals
Adaptation
Size

Related Feeds

Biomaterial Engineering

Advances in biomaterial engineering have permitted the development of sophisticated drug-releasing materials with a biomimetic 3D support that allow a better control of the microenvironment of transplanted cells. Here is the latest research.