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

MOTIVATION: Agent-based modeling is an indispensable tool for studying complex biological systems. However, existing simulation platforms do not always take full advantage of modern hardware and often have a field-specific software design. RESULTS: We present a novel simulation platform called BioDynaMo that alleviates both of these problems. BioDynaMo features a modular and high-performance simulation engine. We demonstrate that BioDynaMo can be used to simulate use cases in: neuroscience, oncology and epidemiology. For each use case, we validate our findings with experimental data or an analytical solution. Our performance results show that BioDynaMo performs up to three orders of magnitude faster than the state-of-the-art baselines. This improvement makes it feasible to simulate each use case with one billion agents on a single server, showcasing the potential BioDynaMo has for computational biology research. AVAILABILITY AND IMPLEMENTATION: BioDynaMo is an open-source project under the Apache 2.0 license and is available at www.biodynamo.org. Instructions to reproduce the results are available in the supplementary information. SUPPLEMENTARY INFORMATION: Available at https://doi.org/10.5281/zenodo.5121618.

An in silico hybrid continuum-/agent-based procedure to modelling cancer development: Interrogating the interplay amongst glioma invasion, vascularity and necrosis.

This paper develops a three-dimensional in silico hybrid model of cancer, which describes the multi-variate phenotypic behaviour of tumour and host cells. The model encompasses the role of cell migration and adhesion, the influence of the extracellular matrix, the effects of oxygen and nutrient availability, and the signalling triggered by chemical cues and growth factors. The proposed in silico hybrid modelling framework combines successfully the advantages of continuum-based and discrete methods, namely the finite element and agent-based method respectively. The framework is thus used to realistically model cancer mechano-biology in a multiscale fashion while maintaining the resolution power of each method in a computationally cost-effective manner. The model is tailored to simulate glioma progression, and is subsequently used to interrogate the balance between the host cells and small sized gliomas, while the go-or-grow phenotype characteristic in glioblastomas is also investigated. Also, cell-cell and cell-matrix interactions are examined with respect to their effect in (macroscopic) tumour growth, brain tissue perfusion and tumour necrosis. Finally, we use the in silico framework to assess differences between low-grade and high-grade glioma growth, demonstrating significant differences in the distribution of cancer as well as host cells, in accordance with reported experimental findings.

Retinal self-organization: a model of retinal ganglion cells and starburst amacrine cells mosaic formation.

Individual retinal cell types exhibit semi-regular spatial patterns called retinal mosaics. Retinal ganglion cells (RGCs) and starburst amacrine cells (SACs) are known to exhibit such layouts. Mechanisms responsible for the formation of mosaics are not well understood but follow three main principles: (i) homotypic cells prevent nearby cells from adopting the same type, (ii) cell tangential migration and (iii) cell death. Alongside experiments in mouse, we use BioDynaMo, an agent-based simulation framework, to build a detailed and mechanistic model of mosaic formation. We investigate the implications of the three theories for RGC's mosaic formation. We report that the cell migration mechanism yields the most regular mosaics. In addition, we propose that low-density RGC type mosaics exhibit on average low regularities, and thus we question the relevance of regular spacing as a criterion for a group of RGCs to form a RGC type. We investigate SAC mosaics formation and interactions between the ganglion cell layer (GCL) and inner nuclear layer (INL) populations. We propose that homotypic interactions between the GCL and INL populations during mosaics creation are required to reproduce the observed SAC mosaics' characteristics. This suggests that the GCL and INL populations of SACs might not be independent during retinal development.