pSpatiocyte: A Parallel Stochastic Method for Particle Reaction-Diffusion Systems
Authors:
Atsushi Miyauchi,
Kazunari Iwamoto,
Satya Nanda Vel Arjunan,
Koichi Takahashi
Abstract:
Computational systems biology has provided plenty of insights into cell biology. Early on, the focus was on reaction networks between molecular species. Spatial distribution only began to be considered mostly within the last decade. However, calculations were restricted to small systems because of tremendously high computational workloads. To date, application to the cell of typical size with mole…
▽ More
Computational systems biology has provided plenty of insights into cell biology. Early on, the focus was on reaction networks between molecular species. Spatial distribution only began to be considered mostly within the last decade. However, calculations were restricted to small systems because of tremendously high computational workloads. To date, application to the cell of typical size with molecular resolution is still far from realization. In this article, we present a new parallel stochastic method for particle reaction-diffusion systems. The program called pSpatiocyte was created bearing in mind reaction networks in biological cells operating in crowded intracellular environments as the primary simulation target. pSpatiocyte employs unique discretization and parallelization algorithms based on a hexagonal close-packed lattice for efficient execution particularly on large distributed memory parallel computers. For two-level parallelization, we introduced isolated subdomain and tri-stage lockstep communication for process-level, and voxel-locking techniques for thread-level. We performed a series of parallel runs on RIKEN's K computer. For a fine lattice that had relatively low occupancy, pSpatiocyte achieved 7686 times speedup with 663552 cores relative to 64 cores from the viewpoint of strong scaling and exhibited 74\% parallel efficiency. As for weak scaling, efficiencies at least 60% were observed up to 663552 cores. In addition to computational performance, diffusion and reaction rates were validated by theory and another well-validated program and had good agreement. Lastly, as a preliminary example of real-world applications, we present a calculation of the MAPK model, a typical reaction network motif in cell signaling pathways.
△ Less
Submitted 12 May, 2016;
originally announced May 2016.
A computational framework for bioimaging simulation
Authors:
Masaki Watabe,
Satya N. V. Arjunan,
Seiya Fukushima,
Kazunari Iwamoto,
Jun Kozuka,
Satomi Matsuoka,
Yuki Shindo,
Masahiro Ueda,
Koichi Takahashi
Abstract:
Using bioimaging technology, biologists have attempted to identify and document analytical interpretations that underlie biological phenomena in biological cells. Theoretical biology aims at distilling those interpretations into knowledge in the mathematical form of biochemical reaction networks and understanding how higher level functions emerge from the combined action of biomolecules. However,…
▽ More
Using bioimaging technology, biologists have attempted to identify and document analytical interpretations that underlie biological phenomena in biological cells. Theoretical biology aims at distilling those interpretations into knowledge in the mathematical form of biochemical reaction networks and understanding how higher level functions emerge from the combined action of biomolecules. However, there still remain formidable challenges in bridging the gap between bioimaging and mathematical modeling. Generally, measurements using fluorescence microscopy systems are influenced by systematic effects that arise from stochastic nature of biological cells, the imaging apparatus, and optical physics. Such systematic effects are always present in all bioimaging systems and hinder quantitative comparison between the cell model and bioimages. Computational tools for such a comparison are still unavailable. Thus, in this work, we present a computational framework for handling the parameters of the cell models and the optical physics governing bioimaging systems. Simulation using this framework can generate digital images of cell simulation results after accounting for the systematic effects. We then demonstrate that such a framework enables comparison at the level of photon-counting units.
△ Less
Submitted 7 July, 2015; v1 submitted 5 November, 2014;
originally announced November 2014.