Spatial cumulant models enable spatially informed treatment strategies and analysis of local interactions in cancer systems

Theoretical and applied cancer studies that use individual-based models (IBMs) have been limited by the lack of a mathematical formulation that enables rigorous analysis of these models. However, spatial cumulant models (SCMs), which have arisen from theoretical ecology, describe population dynamics generated by a specific family of IBMs, namely spatio-temporal point processes (STPPs). SCMs are spatially resolved population models formulated by a system of differential equations that approximate the dynamics of two STPP-generated summary statistics: first-order spatial cumulants (densities), and second-order spatial cumulants (spatial covariances). We exemplify how SCMs can be used in mathematical oncology by modelling theoretical cancer cell populations comprising interacting growth factor-producing and non-producing cells. To formulate model equations, we use computational tools that enable the generation of STPPs, SCMs and mean-field population models (MFPMs) from user-defined model descriptions (Cornell et al., 2019). To calculate and compare STPP, SCM and MFPM-generated summary statistics, we develop an application-agnostic computational pipeline. Our results demonstrate that SCMs can capture STPP-generated population density dynamics, even when MFPMs fail to do so. From both MFPM and SCM equations, we derive treatment-induced death rates required to achieve non-growing cell populations. When testing these treatment strategies in STPP-generated cell populations, our results demonstrate that SCM-informed strategies outperform MFPM-informed strategies in terms of inhibiting population growths. We thus demonstrate that SCMs provide a new framework in which to study cell-cell interactions, and can be used to describe and perturb STPP-generated cell population dynamics. We, therefore, argue that SCMs can be used to increase IBMs’ applicability in cancer research. Statements and Declarations The authors have no competing interests to declare that are relevant to the content of this article.

[1]  Jacob G. Scott,et al.  Measuring competitive exclusion in non–small cell lung cancer , 2020, bioRxiv.

[2]  Joshua M. Weiss,et al.  Spatially resolved transcriptomics reveals the architecture of the tumor-microenvironment interface , 2021, Nature Communications.

[3]  P. Altrock,et al.  Autocrine signaling can explain the emergence of Allee effects in cancer cell populations , 2021, bioRxiv.

[4]  Thomas O. McDonald,et al.  Identification of optimal dosing schedules of dacomitinib and osimertinib for a phase I/II trial in advanced EGFR-mutant non-small cell lung cancer , 2021, Nature Communications.

[5]  H. Mistry On the reporting and analysis of a cancer evolutionary adaptive dosing trial , 2021, Nature communications.

[6]  M. Chaplain,et al.  Targeting Cellular DNA Damage Responses in Cancer: An In Vitro-Calibrated Agent-Based Model Simulating Monolayer and Spheroid Treatment Responses to ATR-Inhibiting Drugs , 2019, Bulletin of Mathematical Biology.

[7]  Panu Somervuo,et al.  A general mathematical method for predicting spatio-temporal correlations emerging from agent-based models , 2020, Journal of the Royal Society Interface.

[8]  M. Oli,et al.  Eco‐oncology: Applying ecological principles to understand and manage cancer , 2020, Ecology and evolution.

[9]  G. Nowicka,et al.  Communication in the Cancer Microenvironment as a Target for Therapeutic Interventions , 2020, Cancers.

[10]  A. Boddy,et al.  Agent‐based modelling reveals strategies to reduce the fitness and metastatic potential of circulating tumour cell clusters , 2020, Evolutionary applications.

[11]  Joel s. Brown,et al.  Eradicating metastatic cancer and the eco-evolutionary dynamics of Anthropocene extinctions. , 2019, Cancer research.

[12]  Panu Somervuo,et al.  A unified framework for analysis of individual-based models in ecology and beyond , 2019, Nature Communications.

[13]  Katarzyna A Rejniak,et al.  Hybrid modeling frameworks of tumor development and treatment , 2019, Wiley interdisciplinary reviews. Systems biology and medicine.

[14]  Michael Hinczewski,et al.  The 2019 mathematical oncology roadmap , 2019, Physical biology.

[15]  Zhihui Wang,et al.  Dynamic Targeting in Cancer Treatment , 2019, Front. Physiol..

[16]  Sara Hamis,et al.  Blackboard to Bedside: A Mathematical Modeling Bottom-Up Approach Toward Personalized Cancer Treatments. , 2019, JCO clinical cancer informatics.

[17]  Randy Heiland,et al.  A Review of Cell-Based Computational Modeling in Cancer Biology , 2019, JCO clinical cancer informatics.

[18]  P. Newton,et al.  Cellular interactions constrain tumor growth , 2019, Proceedings of the National Academy of Sciences.

[19]  Yucheng Dong,et al.  A Unified Framework , 2018, Linguistic Decision Making.

[20]  Jan Poleszczuk,et al.  The Optimal Radiation Dose to Induce Robust Systemic Anti-Tumor Immunity , 2018, International journal of molecular sciences.

[21]  Anudeep Surendran,et al.  Spatial Moment Description of Birth–Death–Movement Processes Incorporating the Effects of Crowding and Obstacles , 2018, Bulletin of Mathematical Biology.

[22]  Adelle C F Coster,et al.  Mathematical Modelling of the Interaction Between Cancer Cells and an Oncolytic Virus: Insights into the Effects of Treatment Protocols , 2018, Bulletin of mathematical biology.

[23]  David Basanta,et al.  Fibroblasts and Alectinib switch the evolutionary games played by non-small cell lung cancer , 2017, Nature Ecology & Evolution.

[24]  Joel s. Brown,et al.  Integrating evolutionary dynamics into treatment of metastatic castrate-resistant prostate cancer , 2017, Nature Communications.

[25]  Matthew J Simpson,et al.  Quantifying rates of cell migration and cell proliferation in co-culture barrier assays reveals how skin and melanoma cells interact during melanoma spreading and invasion. , 2017, Journal of theoretical biology.

[26]  Ursula Klingmüller,et al.  An individual-based model for collective cancer cell migration explains speed dynamics and phenotype variability in response to growth factors , 2017, npj Systems Biology and Applications.

[27]  A. Schäffer,et al.  The evolution of tumour phylogenetics: principles and practice , 2017, Nature Reviews Genetics.

[28]  B. Perthame,et al.  On interfaces between cell populations with different mobilities , 2016 .

[29]  Alex James,et al.  Collective Cell Behaviour with Neighbour-Dependent Proliferation, Death and Directional Bias , 2016, Bulletin of mathematical biology.

[30]  Martin A. Nowak,et al.  A spatial model predicts that dispersal and cell turnover limit intratumour heterogeneity , 2015, Nature.

[31]  M. Archetti,et al.  Heterogeneity for IGF-II production maintained by public goods dynamics in neuroendocrine pancreatic cancer , 2015, Proceedings of the National Academy of Sciences.

[32]  Deborah C Markham,et al.  Choosing an Appropriate Modelling Framework for Analysing Multispecies Co-culture Cell Biology Experiments , 2014, bioRxiv.

[33]  I. Jamall,et al.  Cell-Cell Communication in the Tumor Microenvironment, Carcinogenesis, and Anticancer Treatment , 2014, Cellular Physiology and Biochemistry.

[34]  Otso Ovaskainen,et al.  A general mathematical framework for the analysis of spatiotemporal point processes , 2013, Theoretical Ecology.

[35]  Philip Hahnfeldt,et al.  The Tumor Growth Paradox and Immune System-Mediated Selection for Cancer Stem Cells , 2013, Bulletin of mathematical biology.

[36]  Dmitri Finkelshtein,et al.  Semigroup approach to birth-and-death stochastic dynamics in continuum , 2011, 1109.5094.

[37]  Matthew J Simpson,et al.  Correcting mean-field approximations for birth-death-movement processes. , 2010, Physical review. E, Statistical, nonlinear, and soft matter physics.

[38]  Jason M. Graham,et al.  A spatial model of tumor-host interaction: application of chemotherapy. , 2008, Mathematical biosciences and engineering : MBE.

[39]  Yuri G. Kondratiev,et al.  Markov evolutions and hierarchical equations in the continuum. I: one-component systems , 2007, 0707.0619.

[40]  Yuri Kondratiev,et al.  On non-equilibrium stochastic dynamics for interacting particle systems in continuum , 2008 .

[41]  Xiaodong Cai,et al.  Exact stochastic simulation of coupled chemical reactions with delays. , 2007, The Journal of chemical physics.

[42]  C. Maley,et al.  Cancer is a disease of clonal evolution within the body1–3. This has profound clinical implications for neoplastic progression, cancer prevention and cancer therapy. Although the idea of cancer as an evolutionary problem , 2006 .

[43]  Otso Ovaskainen,et al.  Space and stochasticity in population dynamics , 2006, Proceedings of the National Academy of Sciences.

[44]  Anatoli V. Skorokhod,et al.  ON CONTACT PROCESSES IN CONTINUUM , 2006 .

[45]  U. Dieckmann,et al.  POPULATION GROWTH IN SPACE AND TIME: SPATIAL LOGISTIC EQUATIONS , 2003 .

[46]  B. Bolker,et al.  Using Moment Equations to Understand Stochastically Driven Spatial Pattern Formation in Ecological Systems , 1997, Theoretical population biology.

[47]  R. Durrett,et al.  The Importance of Being Discrete (and Spatial) , 1994 .