Hybrid modelling of stochastic chemical kinetics

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

• Andrew Duncan (University of Sussex; The Alan Turing Institute)
• Thursday 07 April 2016, 09:45-10:30
• Seminar Room 1, Newton Institute.

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

Co-authors: Radek Erban (University of Oxford), Kostantinos Zygalakis (University of Edinburgh)

It is well known that stochasticity can play a fundamental role in various biochemical processes, such as cell regulatory networks and enzyme cascades. Isothermal, well-mixed systems can be adequately modelled by Markov processes and, for such systems, methods such as Gillespie's algorithm are typically employed. While such schemes are easy to implement and are exact, the computational cost of simulating such systems can become prohibitive as the frequency of the reaction events increases. This has motivated numerous coarse grained schemes, where the “fast” reactions are approximated either using Langevin dynamics or deterministically. While such approaches provide a good approximation for systems where all reactants are present in large concentrations, the approximation breaks down when the fast chemical species exist in small concentrations, giving rise to significant errors in the simulation. This is particularly problematic when using such methods to compute statistics of extinction times for chemical species, as well as computing observables of cell cycle models. In this talk, we present a hybrid scheme for simulating well-mixed stochastic kinetics, using Gillepsie-type dynamics to simulate the network in regions of low reactant concentration, and chemical Langevin dynamics when the concentrations of all species is large. These two regimes are coupled via an intermediate region in which a “blended”' jump-diffusion model is introduced. Examples of gene regulatory networks involving reactions occurring at multiple scales, as well as a cell-cycle model are simulated, using the exact and hybrid scheme, and compared, both in terms weak error, as well as computational cost.

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

This talk is included in these lists:

• All CMS events
• Featured lists
• INI info aggregator
• Isaac Newton Institute Seminar Series
• School of Physical Sciences
• Seminar Room 1, Newton Institute
• bld31

Note that ex-directory lists are not shown.