Agent based modelling (ABM) offers a powerful abstraction for scientific study in a broad range of domains. The use of agent based simulators encourages good software engineering design such as separation of concerns, that is, the uncoupling of the model description from its implementation detail. A major limitation in current approaches to ABM simulation is that of the trade off between simulator flexibility and performance. It is common that highly optimised simulations, such as those which target graphics processing units (GPU) hardware, are implemented as standalone software. This work presents a software framework (FLAME GPU 2) which balances flexibility with performance for general purpose ABM. Methods for ensuring high computational efficacy are demonstrated by, minimising data movement, and ensuring high device utilisation by exploiting opportunities for concurrent code execution within a model and through the use of ensembles of simulations. A novel hierarchical sub‐modelling approach is also presented which can be used to model certain types of recursive behaviours. This feature is shown to be essential in providing a mechanism to resolve competition for resources between agents within a parallel environment which would otherwise introduce race conditions. To understand the performance characteristics of the software, a benchmark model with millions of agents is used to explore the use of simulation ensembles and to parametrically investigate concurrent code execution within a model. Performance speedups are demonstrated of 3.5 and 10 respectively over a baseline GPU implementation. Our hierarchical sub‐modelling approach is used to demonstrate the implementation of a recursive algorithm to resolve competition of agent movement which occurs as a result of agent desire to simultaneously occupy discrete areas high in a ‘resource’. The algorithm is used to implement a classical socio‐economics model, Sugarscape, with populations of up to 16M agents.

This work presents an advanced agent-based model developed within the FLAMEGPU2 framework, aimed at simulating the intricate dynamics of cell microenvironments. Our primary objective is to showcase FLAMEGPU2's potential in modelling critical features such as cell-cell and cell-ECM interactions, species diffusion, vascularisation, cell migration, and/or cell cycling. By doing so, we provide a versatile template that serves as a foundational platform for researchers to model specific biological mechanisms or processes. We highlight the utility of our approach as a microscale component within multiscale frameworks. Through four example applications, we demonstrate the model's versatility in capturing phenomena such as strain-stiffening behaviour of hydrogels, cell migration patterns within hydrogels, spheroid formation and fibre reorientation, and the simulation of diffusion processes within a vascularised and deformable domain. This work aims to bridge the gap between computational efficiency and biological fidelity, offering a scalable and flexible platform to advance our understanding of tissue biology and engineering.

Neuroblastoma is a complex and aggressive type of cancer that affects children. Current treatments involve a combination of surgery, chemotherapy, radiotherapy, and stem cell transplantation. However, treatment outcomes vary due to the heterogeneous nature of the disease. Computational models have been used to analyse data, simulate biological processes, and predict disease progression and treatment outcomes. While continuum cancer models capture the overall behaviour of tumours, and agent-based models represent the complex behaviour of individual cells, multiscale models represent interactions at different organisational levels, providing a more comprehensive understanding of the system. In 2018, the PRIMAGE consortium was formed to build a cloud-based decision support system for neuroblastoma, including a multi-scale model for patient-specific simulations of disease progression. In this work we have developed this multi-scale model that includes data such as patient's tumour geometry, cellularity, vascularization, genetics and type of chemotherapy treatment, and integrated it into an online platform that runs the simulations on a high-performance computation cluster using Onedata and Kubernetes technologies. This infrastructure will allow clinicians to optimise treatment regimens and reduce the number of costly and time-consuming clinical trials. This manuscript outlines the challenging framework's model architecture, data workflow, hypothesis, and resources employed in its development.

Better software drives better research, a fundamental principle in the research software engineering community. Similarly, better architecture underpins better software, a core belief in the software engineering research community. Therefore, we advocate and emphasize the importance of designing robust architectures for research software to elevate the quality of research outcomes, and illustrate this with two case studies.

Neuroblastoma is the most common extra-cranial solid tumour in children. Over half of all high-risk cases are expected to succumb to the disease even after chemotherapy, surgery, and immunotherapy. Although the importance of MYCN amplification in this disease is indisputable, the mechanistic details remain enigmatic. Here, we present a multicellular model of neuroblastoma comprising a continuous automaton, discrete cell agents, and a centre-based mechanical model, as well as the simulation results we obtained with it. The continuous automaton represents the tumour microenvironment as a grid-like structure, where each voxel is associated with continuous variables such as the oxygen level therein. Each discrete cell agent is defined by several attributes, including its cell cycle position, mutations, gene expression pattern, and more with behaviours such as cell cycling and cell death being stochastically dependent on these attributes. The centre-based mechanical model represents the properties of these agents as physical objects, describing how they repel each other as soft spheres. By implementing a stochastic simulation algorithm on modern GPUs, we simulated the dynamics of over one million neuroblastoma cells over a period of months. Specifically, we set up 1200 heterogeneous tumours and tracked the MYCN-amplified clone’s dynamics in each, revealed the conditions that favour its growth, and tested its responses to 5000 drug combinations. Our results are in agreement with those reported in the literature and add new insights into how the MYCN-amplified clone’s reproductive advantage in a tumour, its gene expression profile, the tumour’s other clones (with different mutations), and the tumour’s microenvironment are inter-related. Based on the results, we formulated a hypothesis, which argues that there are two distinct populations of neuroblastoma cells in the tumour; the p53 protein is pro-survival in one and pro-apoptosis in the other. It follows that alternating between inhibiting MDM2 to restore p53 activity and inhibiting ARF to attenuate p53 activity is a promising, if unorthodox, therapeutic strategy. The multicellular model has the advantages of modularity, high resolution, and scalability, making it a potential foundation for creating digital twins of neuroblastoma patients.