A seven-step guide to spatial, agent-based modelling of tumour evolution.

Spatial agent-based models are frequently used to investigate the evolution of solid tumours subject to localized cell-cell interactions and microenvironmental heterogeneity. As spatial genomic, transcriptomic and proteomic technologies gain traction, spatial computational models are predicted to become ever more necessary for making sense of complex clinical and experimental data sets, for predicting clinical outcomes, and for optimizing treatment strategies. Here we present a non-technical step by step guide to developing such a model from first principles. Stressing the importance of tailoring the model structure to that of the biological system, we describe methods of increasing complexity, from the basic Eden growth model up to off-lattice simulations with diffusible factors. We examine choices that unavoidably arise in model design, such as implementation, parameterization, visualization and reproducibility. Each topic is illustrated with examples drawn from recent research studies and state of the art modelling platforms. We emphasize the benefits of simpler models that aim to match the complexity of the phenomena of interest, rather than that of the entire biological system. Our guide is aimed at both aspiring modellers and other biologists and oncologists who wish to understand the assumptions and limitations of the models on which major cancer studies now so often depend.

Dissecting the Spatially Restricted Effects of Microenvironment-Mediated Resistance on Targeted Therapy Responses.

The response of tumors to anti-cancer therapies is defined not only by cell-intrinsic therapy sensitivities but also by local interactions with the tumor microenvironment. Fibroblasts that make tumor stroma have been shown to produce paracrine factors that can strongly reduce the sensitivity of tumor cells to many types of targeted therapies. Moreover, a high stroma/tumor ratio is generally associated with poor survival and reduced therapy responses. However, in contrast to advanced knowledge of the molecular mechanisms responsible for stroma-mediated resistance, its effect on the ability of tumors to escape therapeutic eradication remains poorly understood. To a large extent, this gap of knowledge reflects the challenge of accounting for the spatial aspects of microenvironmental resistance, especially over longer time frames. To address this problem, we integrated spatial inferences of proliferation-death dynamics from an experimental animal model of targeted therapy responses with spatial mathematical modeling. With this approach, we dissected the impact of tumor/stroma distribution, magnitude and distance of stromal effects. While all of the tested parameters affected the ability of tumor cells to resist elimination, spatial patterns of stroma distribution within tumor tissue had a particularly strong impact.

The bone ecosystem facilitates multiple myeloma relapse and the evolution of heterogeneous drug resistant disease.

Multiple myeloma (MM) is an osteolytic malignancy that is incurable due to the emergence of treatment resistant disease. Defining how, when and where myeloma cell intrinsic and extrinsic bone microenvironmental mechanisms cause relapse is challenging with current biological approaches. Here, we report a biology-driven spatiotemporal hybrid agent-based model of the MM-bone microenvironment. Results indicate MM intrinsic mechanisms drive the evolution of treatment resistant disease but that the protective effects of bone microenvironment mediated drug resistance (EMDR) significantly enhances the probability and heterogeneity of resistant clones arising under treatment. Further, the model predicts that targeting of EMDR deepens therapy response by eliminating sensitive clones proximal to stroma and bone, a finding supported by in vivo studies. Altogether, our model allows for the study of MM clonal evolution over time in the bone microenvironment and will be beneficial for optimizing treatment efficacy so as to significantly delay disease relapse.

Peristromal niches protect lung cancers from targeted therapies through a combined effect of multiple molecular mediators.

Targeted therapies directed against oncogenic signaling addictions, such as inhibitors of ALK in ALK+ NSCLC often induce strong and durable clinical responses. However, they are not curative in metastatic cancers, as some tumor cells persist through therapy, eventually developing resistance. Therapy sensitivity can reflect not only cell-intrinsic mechanisms but also inputs from stromal microenvironment. Yet, the contribution of tumor stroma to therapeutic responses in vivo remains poorly defined. To address this gap of knowledge, we assessed the contribution of stroma-mediated resistance to therapeutic responses to the frontline ALK inhibitor alectinib in xenograft models of ALK+ NSCLC. We found that stroma-proximal tumor cells are partially protected against cytostatic effects of alectinib. This effect is observed not only in remission, but also during relapse, indicating the strong contribution of stroma-mediated resistance to both persistence and resistance. This therapy-protective effect of the stromal niche reflects a combined action of multiple mechanisms, including growth factors and extracellular matrix components. Consequently, despite improving alectinib responses, suppression of any individual resistance mechanism was insufficient to fully overcome the protective effect of stroma. Focusing on shared collateral sensitivity of persisters offered a superior therapeutic benefit, especially when using an antibody-drug conjugate with bystander effect to limit therapeutic escape. These findings indicate that stroma-mediated resistance might be the major contributor to both residual and progressing disease and highlight the limitation of focusing on suppressing a single resistance mechanism at a time.