A scalable solver for a stochastic, hybrid cellular automaton model of personalized breast cancer therapy
Open Access
- 30 October 2021
- journal article
- research article
- Published by Wiley in International Journal for Numerical Methods in Biomedical Engineering
- Vol. 38 (1) , e3542
Abstract
Mathematical modeling and simulation is a promising approach to personalized cancer medicine. Yet, the complexity, heterogeneity and multi‐scale nature of cancer pose significant computational challenges. Coupling discrete cell‐based models with continuous models using hybrid cellular automata (CA) is a powerful approach for mimicking biological complexity and describing the dynamical exchange of information across different scales. However, when clinically relevant cancer portions are taken into account, such models become computationally very expensive. While efficient parallelization techniques for continuous models exist, their coupling with discrete models, particularly CA, necessitates more elaborate solutions. Building upon FEniCS, a popular and powerful scientific computing platform for solving partial differential equations, we developed parallel algorithms to link stochastic CA with differential equations (https://bitbucket.org/HTasken/cansim). The algorithms minimize the communication between processes that share CA neighborhood values while also allowing for reproducibility during stochastic updates. We demonstrated the potential of our solution on a complex hybrid cellular automaton model of breast cancer treated with combination chemotherapy. On a single‐core processor, we obtained nearly linear scaling with an increasing problem size, whereas weak parallel scaling showed moderate growth in solving time relative to increase in problem size. Finally, we applied the algorithm to a problem that is 500 times larger than previous work, allowing us to run personalized therapy simulations based on heterogeneous cell density and tumor perfusion conditions estimated from magnetic resonance imaging data on an unprecedented scale.Keywords
All Related Versions
Funding Information
- H2020 European Research Council
- Horizon 2020 Framework Programme
- Norges Forskningsråd
- Universitetet i Oslo
This publication has 18 references indexed in Scilit:
- Biocellion : accelerating computer simulation of multicellular biological system modelsBioinformatics, 2014
- Morpheus: a user-friendly modeling environment for multiscale and multicellular systems biologyBioinformatics, 2014
- Chaste: An Open Source C++ Library for Computational Physiology and BiologyPLOS Computational Biology, 2013
- Modelling the effects of cell-cycle heterogeneity on the response of a solid tumour to chemotherapy: Biological insights from a hybrid multiscale cellular automaton modelJournal of Theoretical Biology, 2012
- Multiscale Modelling of Vascular Tumour Growth in 3D: The Roles of Domain Size and Boundary ConditionsPLOS ONE, 2011
- Mathematical Modeling Predicts Synergistic Antitumor Effects of Combining a Macrophage-Based, Hypoxia-Targeted Gene Therapy with ChemotherapyCancer Research, 2011
- A cell-based simulation software for multi-cellular systemsBioinformatics, 2010
- An evolutionary hybrid cellular automaton model of solid tumour growthJournal of Theoretical Biology, 2007
- A Multiple Scale Model for Tumor GrowthMultiscale Modeling & Simulation, 2005
- hypre: A Library of High Performance Preconditioners2002