Dynamic partitioning for hybrid simulation of the bistable HIV-1 transactivation network

Mark Griffith, Tod Courtney, Jean Peccoud, William H. Sanders

Research output: Contribution to journalArticlepeer-review

Abstract

Motivation: The stochastic kinetics of a well-mixed chemical system, governed by the chemical Master equation, can be simulated using the exact methods of Gillespie. However, these methods do not scale well as systems become more complex and larger models are built to include reactions with widely varying rates, since the computational burden of simulation increases with the number of reaction events. Continuous models may provide an approximate solution and are computationally less costly, but they fail to capture the stochastic behavior of small populations of macromolecules. Results: In this article we present a hybrid simulation algorithm that dynamically partitions the system into subsets of continuous and discrete reactions, approximates the continuous reactions deterministically as a system of ordinary differential equations (ODE) and uses a Monte Carlo method for generating discrete reaction events according to a time-dependent propensity. Our approach to partitioning is improved such that we dynamically partition the system of reactions, based on a threshold relative to the distribution of propensities in the discrete subset. We have implemented the hybrid algorithm in an extensible framework, utilizing two rigorous ODE solvers to approximate the continuous reactions, and use an example model to illustrate the accuracy and potential speedup of the algorithm when compared with exact stochastic simulation.

Original languageEnglish (US)
Pages (from-to)2782-2789
Number of pages8
JournalBioinformatics
Volume22
Issue number22
DOIs
StatePublished - Nov 15 2006

ASJC Scopus subject areas

  • Statistics and Probability
  • Biochemistry
  • Molecular Biology
  • Computer Science Applications
  • Computational Theory and Mathematics
  • Computational Mathematics

Fingerprint

Dive into the research topics of 'Dynamic partitioning for hybrid simulation of the bistable HIV-1 transactivation network'. Together they form a unique fingerprint.

Cite this