Show simple item record

dc.identifier.urihttp://hdl.handle.net/1951/56069
dc.identifier.urihttp://hdl.handle.net/11401/71654
dc.description.sponsorshipThis work is sponsored by the Stony Brook University Graduate School in compliance with the requirements for completion of degree.en_US
dc.formatMonograph
dc.format.mediumElectronic Resourceen_US
dc.language.isoen_US
dc.publisherThe Graduate School, Stony Brook University: Stony Brook, NY.
dc.typeDissertation
dcterms.abstractWe are interested in predicting the time dependent behavior of biochemical networks such as interaction between proteins. These networks are represented by a system of chemical reactions. In the forward problem, we have knowledge of all the reactions and the associated rate constants. We want to determine the joint probability density of the populations of all the molecular species at any time instant. The given biochemical system is a discrete state, continuous time Markov process, and the time evolution of its probability density function is described by the chemical master equation (CME). We want to obtain the solution of the master equation in complex biochemical systems with a large number of species and reactions. Analytical solutions of the chemical master equation for first and second order reactions have only been obtained for select cases with a few species and reactions. Another approach to solving the problem is to approximate the CME with the Fokker-Planck equations. But this would require solving partial differential equations with a large number of variables. The present state of the art approaches for stochastic simulation of such systems are based on Monte Carlo methods. One such popular method is the Stochastic Simulation Algorithm (SSA) derived by Gillespie in 1976. Several authors have developed accelerated versions of the SSA such as the Next Reaction Method and time leaping methods, in order to reduce the computation time of SSA. The Monte Carlo methods provide approximations of the complete distribution, but they require simulations of many realizations of the Markov process and many time steps. Hence their computation times are prohibitively long for very complex systems. Methods for modeling the biochemical networks based on moment propagation is a relatively unexplored area. We propose a new method for propagating the first two moments of the joint probability distribution of the number of molecules. In many systems, the distribution can be approximated as Gaussians and therefore computing the first two moments is sufficient. Simulation results show that our method yields accurate results for first order and second order reactions. Compared with the Monte Carlo methods, our method yields significant savings in computation time. Compared with other moment propagation methods, the recursive expressions in our method can be implemented by specifying rate constants and stoichiometries, without having to derive or solve any differential equations. Whereas other moment propagation methods with similar accuracy have been demonstrated for a few species, we demonstrate our method for complex biochemical systems with hundreds of species.
dcterms.available2012-05-17T12:21:39Z
dcterms.available2015-04-24T14:48:24Z
dcterms.contributorPetar M. Djuric.en_US
dcterms.contributorMurali Subbaraoen_US
dcterms.contributorMonica F. Bugalloen_US
dcterms.contributorSanford Simon.en_US
dcterms.creatorMane, Vibha
dcterms.dateAccepted2012-05-17T12:21:39Z
dcterms.dateAccepted2015-04-24T14:48:24Z
dcterms.dateSubmitted2012-05-17T12:21:39Z
dcterms.dateSubmitted2015-04-24T14:48:24Z
dcterms.descriptionDepartment of Electrical Engineeringen_US
dcterms.formatApplication/PDFen_US
dcterms.formatMonograph
dcterms.identifierMane_grad.sunysb_0771E_10426.pdfen_US
dcterms.identifierhttp://hdl.handle.net/1951/56069
dcterms.identifierhttp://hdl.handle.net/11401/71654
dcterms.issued2011-05-01
dcterms.languageen_US
dcterms.provenanceMade available in DSpace on 2012-05-17T12:21:39Z (GMT). No. of bitstreams: 1 Mane_grad.sunysb_0771E_10426.pdf: 1356206 bytes, checksum: 3c7744bd4389e1d09aeaff228c95dcfb (MD5) Previous issue date: 1en
dcterms.provenanceMade available in DSpace on 2015-04-24T14:48:24Z (GMT). No. of bitstreams: 3 Mane_grad.sunysb_0771E_10426.pdf.jpg: 1894 bytes, checksum: a6009c46e6ec8251b348085684cba80d (MD5) Mane_grad.sunysb_0771E_10426.pdf: 1356206 bytes, checksum: 3c7744bd4389e1d09aeaff228c95dcfb (MD5) Mane_grad.sunysb_0771E_10426.pdf.txt: 139336 bytes, checksum: ee570d9fbcaab96fa4fd3aa026493690 (MD5) Previous issue date: 1en
dcterms.publisherThe Graduate School, Stony Brook University: Stony Brook, NY.
dcterms.subjectElectrical Engineering -- Bioinformatics
dcterms.subjectbiochemical processes, moment closure, moment propagation methods, stochastic simulation
dcterms.titleMoment Propagation Methods for Stochastic Simulation of Complex Biochemical Systems
dcterms.typeDissertation


Files in this item

Thumbnail

This item appears in the following Collection(s)

Show simple item record