University of Cambridge > Talks.cam > Isaac Newton Institute Seminar Series > Closure Scheme for Chemical Master Equations - Is the Gibbs entropy maximum for stochastic reaction networks at steady state?

Closure Scheme for Chemical Master Equations - Is the Gibbs entropy maximum for stochastic reaction networks at steady state?

Add to your list(s) Download to your calendar using vCal

If you have a question about this talk, please contact INI IT.

SDBW03 - Advances in numerical and analytic approaches for the study of non-spatial stochastic dynamical systems in molecular biology

Stochasticity is a defining feature of biochemical reaction networks, with molecular fluctuations influencing cell physiology. In principle, master probability equations completely govern the dynamic and steady state behavior of stochastic reaction networks. In practice, a solution had been elusive for decades, when there are second or higher order reactions. A large community of scientists has then reverted to merely sampling the probability distribution of biological networks with stochastic simulation algorithms. Consequently, master equations, for all their promise, have not inspired biological discovery.

We recently presented a closure scheme that solves chemical master equations of nonlinear reaction networks [1]. The zero-information closure (ZI-closure) scheme is founded on the observation that although higher order probability moments are not numerically negligible, they contain little information to reconstruct the master probability [2]. Higher order moments are then related to lower order ones by maximizing the entropy of the network. Using several examples, we show that moment-closure techniques may afford the quick and accurate calculation of steady-state distributions of arbitrary reaction networks.

With the ZI-closure scheme, the stability of the systems around steady states can be quantitatively assessed computing eigenvalues of the moment Jacobian [3]. This is analogous to Lyapunov’s stability analysis of deterministic dynamics and it paves the way for a stability theory and the design of controllers of stochastic reacting systems [4, 5].

In this seminar, we will present the ZI-closure scheme, the calculation of steady state probability distributions, and discuss the stability of stochastic systems.

1. Smadbeck P, Kaznessis YN. A closure scheme for chemical master equations. Proc Natl Acad Sci U S A . 2013 Aug 27;110(35):14261-5.

2. Smadbeck P, Kaznessis YN. Efficient moment matrix generation for arbitrary chemical networks, Chem Eng Sci, 20

This talk is part of the Isaac Newton Institute Seminar Series series.

Tell a friend about this talk:

This talk is included in these lists:

Note that ex-directory lists are not shown.

 

© 2006-2024 Talks.cam, University of Cambridge. Contact Us | Help and Documentation | Privacy and Publicity