Dataset Information

Stochastic reaction-diffusion kinetics in the microscopic limit.

ABSTRACT: Quantitative analysis of biochemical networks often requires consideration of both spatial and stochastic aspects of chemical processes. Despite significant progress in the field, it is still computationally prohibitive to simulate systems involving many reactants or complex geometries using a microscopic framework that includes the finest length and time scales of diffusion-limited molecular interactions. For this reason, spatially or temporally discretized simulations schemes are commonly used when modeling intracellular reaction networks. The challenge in defining such coarse-grained models is to calculate the correct probabilities of reaction given the microscopic parameters and the uncertainty in the molecular positions introduced by the spatial or temporal discretization. In this paper we have solved this problem for the spatially discretized Reaction-Diffusion Master Equation; this enables a seamless and physically consistent transition from the microscopic to the macroscopic frameworks of reaction-diffusion kinetics. We exemplify the use of the methods by showing that a phosphorylation-dephosphorylation motif, commonly observed in eukaryotic signaling pathways, is predicted to display fluctuations that depend on the geometry of the system.

SUBMITTER: Fange D

PROVIDER: S-EPMC2993376 | biostudies-literature | 2010 Nov

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Stochastic reaction-diffusion kinetics in the microscopic limit.

Fange David D Berg Otto G OG Sjöberg Paul P Elf Johan J

Proceedings of the National Academy of Sciences of the United States of America 20101101 46

Quantitative analysis of biochemical networks often requires consideration of both spatial and stochastic aspects of chemical processes. Despite significant progress in the field, it is still computationally prohibitive to simulate systems involving many reactants or complex geometries using a microscopic framework that includes the finest length and time scales of diffusion-limited molecular interactions. For this reason, spatially or temporally discretized simulations schemes are commonly used ...[more]

PMID: 21041672

Similar Datasets

Project description:BACKGROUND: Models of cellular molecular systems are built from components such as biochemical reactions (including interactions between ligands and membrane-bound proteins), conformational changes and active and passive transport. A discrete, stochastic description of the kinetics is often essential to capture the behavior of the system accurately. Where spatial effects play a prominent role the complex morphology of cells may have to be represented, along with aspects such as chemical localization and diffusion. This high level of detail makes efficiency a particularly important consideration for software that is designed to simulate such systems. RESULTS: We describe STEPS, a stochastic reaction-diffusion simulator developed with an emphasis on simulating biochemical signaling pathways accurately and efficiently. STEPS supports all the above-mentioned features, and well-validated support for SBML allows many existing biochemical models to be imported reliably. Complex boundaries can be represented accurately in externally generated 3D tetrahedral meshes imported by STEPS. The powerful Python interface facilitates model construction and simulation control. STEPS implements the composition and rejection method, a variation of the Gillespie SSA, supporting diffusion between tetrahedral elements within an efficient search and update engine. Additional support for well-mixed conditions and for deterministic model solution is implemented. Solver accuracy is confirmed with an original and extensive validation set consisting of isolated reaction, diffusion and reaction-diffusion systems. Accuracy imposes upper and lower limits on tetrahedron sizes, which are described in detail. By comparing to Smoldyn, we show how the voxel-based approach in STEPS is often faster than particle-based methods, with increasing advantage in larger systems, and by comparing to MesoRD we show the efficiency of the STEPS implementation. CONCLUSION: STEPS simulates models of cellular reaction-diffusion systems with complex boundaries with high accuracy and high performance in C/C++, controlled by a powerful and user-friendly Python interface. STEPS is free for use and is available at http://steps.sourceforge.net/

Dataset Information

Stochastic reaction-diffusion kinetics in the microscopic limit.

Publications

Stochastic reaction-diffusion kinetics in the microscopic limit.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets