DSpace Repository

Moment Propagation Methods for Stochastic Simulation of Complex Biochemical Systems

Show simple item record

dc.contributor.advisor Petar M. Djuric. en_US
dc.contributor.author Mane, Vibha en_US
dc.contributor.other Department of Electrical Engineering en_US
dc.date.accessioned 2012-05-17T12:21:39Z
dc.date.accessioned 2015-04-24T14:48:24Z
dc.date.available 2012-05-17T12:21:39Z
dc.date.available 2015-04-24T14:48:24Z
dc.date.issued 2011-05-01 en_US
dc.identifier Mane_grad.sunysb_0771E_10426.pdf en_US
dc.identifier.uri http://hdl.handle.net/1951/56069 en_US
dc.identifier.uri http://hdl.handle.net/11401/71654 en_US
dc.description.abstract We 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. en_US
dc.description.sponsorship This work is sponsored by the Stony Brook University Graduate School in compliance with the requirements for completion of degree. en_US
dc.format Monograph en_US
dc.format.medium Electronic Resource en_US
dc.language.iso en_US en_US
dc.publisher The Graduate School, Stony Brook University: Stony Brook, NY. en_US
dc.subject.lcsh Electrical Engineering -- Bioinformatics en_US
dc.subject.other biochemical processes, moment closure, moment propagation methods, stochastic simulation en_US
dc.title Moment Propagation Methods for Stochastic Simulation of Complex Biochemical Systems en_US
dc.type Dissertation en_US
dc.mimetype Application/PDF en_US
dc.contributor.committeemember Murali Subbarao en_US
dc.contributor.committeemember Monica F. Bugallo en_US
dc.contributor.committeemember Sanford Simon. en_US

Files in this item

This item appears in the following Collection(s)

Show simple item record

Search DSpace

Advanced Search


My Account