Adaptive stochastic-deterministic chemical kinetic simulations.
BACKGROUND
Biochemical signaling pathways and genetic circuits often involve very small numbers of key signaling molecules. Computationally expensive stochastic methods are necessary to simulate such chemical situations. Single-molecule chemical events often co-exist with much larger numbers of signaling molecules where mass-action kinetics is a reasonable approximation. Here, we describe an adaptive stochastic method that dynamically chooses between deterministic and stochastic calculations depending on molecular count and propensity of forward reactions. The method is fixed timestep and has first order accuracy. We compare the efficiency of this method with exact stochastic methods.


RESULTS
We have implemented an adaptive stochastic-deterministic approximate simulation method for chemical kinetics. With an error margin of 5%, the method solves typical biologically constrained reaction schemes more rapidly than exact stochastic methods for reaction volumes >1-10 micro m(3). We have developed a test suite of reaction cases to test the accuracy of mixed simulation methods.


BACKGROUND
Simulation software used in the paper is freely available from http://www.ncbs.res.in/kinetikit/download.html
Vasudeva
Karan
K
National Centre for Biological Sciences, TIFR, GKVK Campus, Bangalore 560065, India.
Bhalla
Upinder S
US
D000465:Algorithms; D020650:Combinatorial Chemistry Techniques; D005246:Feedback; D007700:Kinetics; D008660:Metabolism / Q000502:physiology*; D008954:Models, Biological; D008956:Models, Chemical; D015233:Models, Statistical; D009097:Multienzyme Complexes / Q000737:chemistry / Q000378:metabolism; D015203:Reproducibility of Results; D012680:Sensitivity and Specificity; D015398:Signal Transduction / Q000502:physiology*; D013269:Stochastic Processes
D009097:Multienzyme Complexes