BioDynaMo: a modular platform for high-performance agent-based simulation

Computer simulation is an indispensable tool for studying complex biological systems. In particular, agent-based modeling is an attractive method to describe biophysical dynamics. However, two barriers limit faster progress. First, simulators do not always take full advantage of parallel and heterogeneous hardware. Second, many agent-based simulators are written with a specific research problem in mind and lack a flexible software design. Consequently, researchers have to spend an unnecessarily long time implementing their simulation and have to compromise either on model resolution or system size. We present a novel simulation platform called BioDynaMo that alleviates both of these problems researchers face in computer simulation of complex biological systems. BioDynaMo features a general-purpose and high-performance simulation engine. The engine simulates cellular elements, their interactions within a 3D physical environment, and their cell-internal genetic dynamics. Cell-internal dynamics can be described in C++ code or using system biology markup language (SBML). We demonstrate BioDynaMo’s wide range of application with three example use cases: soma clustering, neural development, and tumor spheroid growth. We validate our results with experimental data, and evaluate the performance of the simulation engine. We compare BioDynaMo’s performance with a state-of-the-art baseline, and analyze its scalability. We observe a speedup of 20–124× over the state-of-the-art baseline using one CPU core and a parallel speedup between 67× and 76× using 72 physical CPU cores with hyperthreading enabled. Combining these two results, we conclude that, on our test system, BioDynaMo is at least three orders of magnitude faster than the state-of-the-art serial baseline. These improvements make it feasible to simulate neural development with 1.24 billion agents on a single server with 1TB memory, and 12 million agents on a laptop with 16GB memory. BioDynaMo is an open-source project under the Apache 2.0 license and is available at www.biodynamo.org. Author summary Computer simulations of biological systems are crucial to gain insights into complex processes of living organisms. However, the development of increasingly large and complex simulations is a difficult task, partly because a strong background in biology as well as computer science is required. In this paper, we introduce BioDynaMo, an agent-based simulation platform with which life scientists can create simulations that are three orders of magnitude faster than the state-of-the-art baseline. By taking advantage of the latest developments in computing hardware, we build a platform that is highly optimized. This enables the simulation of 1.24 billion agents on a single server and 12 million agents on a laptop. BioDynaMo places a lot of focus on hiding computational complexity and providing an easy-to-use interface, such that the life scientist can concentrate on biological aspects, rather than computational. BioDynaMo helps scientists to translate an idea quickly into a simulation by providing common building blocks, and a modular and extensible software design. We analyze the performance of the platform and demonstrate the capabilities with three example use cases: soma clustering, neural development, and tumor spheroid growth. The results support the view that BioDynaMo will help open up new research opportunities for large-scale biological simulations.

[1]  Daniela M. Romano,et al.  High performance cellular level agent-based simulation with FLAME for the GPU , 2010, Briefings Bioinform..

[2]  Elisabete C Costa,et al.  3D tumor spheroids as in vitro models to mimic in vivo human solid tumors resistance to therapeutic drugs , 2018, Biotechnology and bioengineering.

[3]  Randy Heiland,et al.  PhysiCell: An open source physics-based cell simulator for 3-D multicellular systems , 2017, bioRxiv.

[4]  Andrew Phillips,et al.  Computational modeling of synthetic microbial biofilms. , 2012, ACS synthetic biology.

[5]  Philipp M. Altrock,et al.  The mathematics of cancer: integrating quantitative models , 2015, Nature Reviews Cancer.

[6]  Steven F. Railsback,et al.  Agent-Based and Individual-Based Modeling: A Practical Introduction , 2011 .

[7]  J. Jefferys,et al.  Specific cytoarchitectureal changes in hippocampal subareas in daDREAM mice , 2016, Molecular Brain.

[8]  Cristian Picioreanu,et al.  iDynoMiCS: next-generation individual-based modelling of biofilms. , 2011, Environmental microbiology.

[9]  Somchai Pattana,et al.  Division d'un milieu cellulaire sous contraintes mécaniques : utilisation de la mécanique des matériaux granulaires , 2006 .

[10]  J. Gurdon,et al.  Morphogen gradient interpretation , 2001, Nature.

[11]  C. C. Law,et al.  ParaView: An End-User Tool for Large-Data Visualization , 2005, The Visualization Handbook.

[12]  Greg Miller,et al.  A Scientist's Nightmare: Software Problem Leads to Five Retractions , 2006, Science.

[13]  O. Sten-Knudsen,et al.  Passive Transport Processes , 1978 .

[14]  Thomas E. Yankeelov,et al.  Multi-scale Modeling in Clinical Oncology: Opportunities and Barriers to Success , 2016, Annals of Biomedical Engineering.

[15]  Alexander G. Fletcher,et al.  Chaste: An Open Source C++ Library for Computational Physiology and Biology , 2013, PLoS Comput. Biol..

[16]  Frederico A. C. Azevedo,et al.  Equal numbers of neuronal and nonneuronal cells make the human brain an isometrically scaled‐up primate brain , 2009, The Journal of comparative neurology.

[17]  Kent L. Beck,et al.  Extreme programming explained - embrace change , 1990 .

[18]  Maciej Cytowski,et al.  Large-Scale Parallel Simulations of 3D Cell Colony Dynamics , 2014, Computing in Science & Engineering.

[19]  Erik De Schutter,et al.  Context-aware modeling of neuronal morphologies , 2014, Front. Neuroanat..

[20]  Hiroaki Kitano,et al.  The systems biology markup language (SBML): a medium for representation and exchange of biochemical network models , 2003, Bioinform..

[21]  Fons Rademakers,et al.  ROOT — An object oriented data analysis framework , 1997 .

[22]  L. Levin,et al.  Biodiversity on the Rocks: Macrofauna Inhabiting Authigenic Carbonate at Costa Rica Methane Seeps , 2015, PloS one.

[23]  Tianlin Liu,et al.  Generation of Multicellular Tumor Spheroids with Microwell-Based Agarose Scaffolds for Drug Testing , 2015, PloS one.

[24]  Marco Manca,et al.  An in silico hybrid continuum-/agent-based procedure to modelling cancer development: interrogating the interplay amongst glioma invasion, vascularity and necrosis. , 2020, Methods.

[25]  Jacky L. Snoep,et al.  BioModels Database: a free, centralized database of curated, published, quantitative kinetic models of biochemical and cellular systems , 2005, Nucleic Acids Res..

[26]  Zhiwei Ji,et al.  Mathematical and Computational Modeling in Complex Biological Systems , 2017, BioMed research international.

[27]  K. Basler,et al.  Dpp controls growth and patterning in Drosophila wing precursors through distinct modes of action , 2017, eLife.

[28]  Vittorio Cristini,et al.  Patient-calibrated agent-based modelling of ductal carcinoma in situ (DCIS): from microscopic measurements to macroscopic predictions of clinical progression. , 2012, Journal of theoretical biology.

[29]  Randal A. Koene,et al.  NETMORPH: A Framework for the Stochastic Generation of Large Scale Neuronal Networks With Realistic Neuron Morphologies , 2009, Neuroinformatics.

[30]  Nicholas S. Flann,et al.  Biocellion: accelerating computer simulation of multicellular biological system models , 2014, Bioinform..

[31]  I W Taylor,et al.  Cell proliferation kinetics of MCF-7 human mammary carcinoma cells in culture and effects of tamoxifen on exponentially growing and plateau-phase cells. , 1983, Cancer research.

[32]  Mario di Bernardo,et al.  BSim 2.0: An Advanced Agent-Based Cell Simulator. , 2017, ACS synthetic biology.

[33]  G. Amdhal,et al.  Validity of the single processor approach to achieving large scale computing capabilities , 1967, AFIPS '67 (Spring).

[34]  K. Perez Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment , 2014 .

[35]  G.E. Moore,et al.  Cramming More Components Onto Integrated Circuits , 1998, Proceedings of the IEEE.

[36]  Michael J. North,et al.  AgentCell: a digital single-cell assay for bacterial chemotaxis , 2005, Bioinform..

[37]  Rodney J. Douglas,et al.  Frontiers in Computational Neuroscience , 2022 .

[38]  R.H. Dennard,et al.  Design Of Ion-implanted MOSFET's with Very Small Physical Dimensions , 1974, Proceedings of the IEEE.