Hybrid modelling of stochastic chemical kinetics
53 mins 33 secs,
244.79 MB,
WebM
640x360,
29.97 fps,
44100 Hz,
624.13 kbits/sec
Share this media item:
Embed this media item:
Embed this media item:
About this item
| Description: |
Zygalakis, K (University of Edinburgh)
Wednesday 3rd February 2016 - 15:00 to 16:00 |
|---|
| Created: | 2016-02-18 15:03 |
|---|---|
| Collection: | Stochastic Dynamical Systems in Biology: Numerical Methods and Applications |
| Publisher: | Isaac Newton Institute |
| Copyright: | Zygalakis, K |
| Language: | eng (English) |
| Distribution: |
World
(downloadable)
|
| Explicit content: | No |
| Aspect Ratio: | 16:9 |
| Screencast: | No |
| Bumper: | UCS Default |
| Trailer: | UCS Default |
| Abstract: | 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. |
|---|---|
Available Formats
| Format | Quality | Bitrate | Size | |||
|---|---|---|---|---|---|---|
| MPEG-4 Video | 640x360 | 1.93 Mbits/sec | 778.82 MB | View | Download | |
| WebM * | 640x360 | 624.13 kbits/sec | 244.79 MB | View | Download | |
| iPod Video | 480x270 | 521.96 kbits/sec | 204.66 MB | View | Download | |
| MP3 | 44100 Hz | 249.78 kbits/sec | 98.03 MB | Listen | Download | |
| Auto | (Allows browser to choose a format it supports) | |||||