Predator-Prey in Tumor-Immune Interactions: A Wrong Model or Just an Incomplete One?

The Hallmarks of Cancer as Eco-Evolutionary Processes

SIGNIFICANCE

Viewing the hallmarks as a sequence of adaptations captures the "why" behind the "how" of the molecular changes driving cancer. This eco-evolutionary view distils the complexity of cancer progression into logical steps, providing a framework for understanding all existing and emerging hallmarks of cancer and developing therapeutic interventions.

Sustainability of a Three-Species Predator-Prey Model in Tumor-Immune Dynamics with Periodic Treatment

Using a tumor-immune growth model, we investigate how immunotherapy affects its dynamical characteristics. Specifically, we extend the prey-predator model of tumor cells and immune cells by including periodic immunotherapy, the nonlinear damping of cancer cells, and the dynamics of a healthy cell population, and investigate the effects of the model parameters. The ideal value of immunotherapy, which promotes the growth of immune (and healthy) cells while contributing to the elimination or control of the cancer cells, is determined by using Fisher information as a measure of variability throughout our study.

Integration of Immune Cell-Target Cell Conjugate Dynamics Changes the Time Scale of Immune Control of Cancer

The human immune system can recognize, attack, and eliminate cancer cells, but cancers can escape this immune surveillance. Variants of ecological predator-prey models can capture the dynamics of such cancer control mechanisms by adaptive immune system cells. These dynamical systems describe, e.g., tumor cell-effector T cell conjugation, immune cell activation, cancer cell killing, and T cell exhaustion. Target (tumor) cell-T cell conjugation is integral to the adaptive immune system's cancer control and immunotherapy. However, whether conjugate dynamics should be explicitly included in mathematical models of cancer-immune interactions is incompletely understood. Here, we analyze the dynamics of a cancer-effector T cell system and focus on the impact of explicitly modeling the conjugate compartment to investigate the role of cellular conjugate dynamics. We formulate a deterministic modeling framework to compare possible equilibria and their stability, such as tumor extinction, tumor-immune coexistence (tumor control), or tumor escape. We also formulate the stochastic analog of this system to analyze the impact of demographic fluctuations that arise when cell populations are small. We find that explicit consideration of a conjugate compartment can (i) change long-term steady-state, (ii) critically change the time to reach an equilibrium, (iii) alter the probability of tumor escape, and (iv) lead to very different extinction time distributions. Thus, we demonstrate the importance of the conjugate compartment in defining tumor-effector T cell interactions. Accounting for transitionary compartments of cellular interactions may better capture the dynamics of tumor control and progression.

Spatial interactions modulate tumor growth and immune infiltration

Direct observation of tumor-immune interactions is unlikely in tumors with currently available technology, but computational simulations based on clinical data can provide insight to test hypotheses. It is hypothesized that patterns of collagen evolve as a mechanism of immune escape, but the exact nature of immune-collagen interactions is poorly understood. Spatial data quantifying collagen fiber alignment in squamous cell carcinomas indicates that late-stage disease is associated with highly aligned fibers. Our computational modeling framework discriminates between two hypotheses: immune cell migration that moves (1) parallel or (2) perpendicular to collagen fiber orientation. The modeling recapitulates immune-extracellular matrix interactions where collagen patterns provide immune protection, leading to an emergent inverse relationship between disease stage and immune coverage. Here, computational modeling provides important mechanistic insights by defining a kernel cell-cell interaction function that considers a spectrum of local (cell-scale) to global (tumor-scale) spatial interactions. Short-range interaction kernels provide a mechanism for tumor cell survival under conditions with strong Allee effects, while asymmetric tumor-immune interaction kernels lead to poor immune response. Thus, the length scale of tumor-immune interaction kernels drives tumor growth and infiltration.

Modeling tumors as complex ecosystems

Many cancers resist therapeutic intervention. This is fundamentally related to intratumor heterogeneity: multiple cell populations, each with different phenotypic signatures, coexist within a tumor and its metastases. Like species in an ecosystem, cancer populations are intertwined in a complex network of ecological interactions. Most mathematical models of tumor ecology, however, cannot account for such phenotypic diversity or predict its consequences. Here, we propose that the generalized Lotka-Volterra model (GLV), a standard tool to describe species-rich ecological communities, provides a suitable framework to model the ecology of heterogeneous tumors. We develop a GLV model of tumor growth and discuss how its emerging properties provide a new understanding of the disease. We discuss potential extensions of the model and their application to phenotypic plasticity, cancer-immune interactions, and metastatic growth. Our work outlines a set of questions and a road map for further research in cancer ecology.

Understanding the role of B cells in CAR T-cell therapy in leukemia through a mathematical model

Chimeric antigen receptor T (CAR T) cell therapy has been proven to be successful against a variety of leukemias and lymphomas. This paper undertakes an analytical and numerical study of a mathematical model describing the competition of CAR T, leukemia, tumor, and B cells. Considering its significance in sustaining anti-CD19 CAR T-cell stimulation, a B-cell source term is integrated into the model. Through stability and bifurcation analyses, the potential for tumor eradication, contingent on the continuous influx of B cells, has been revealed, showing a transcritical bifurcation at a critical B-cell input. Additionally, an almost heteroclinic cycle between equilibrium points is identified, providing a theoretical basis for understanding disease relapse. Analyzing the oscillatory behavior of the system, the time-dependent dynamics of CAR T cells and leukemic cells can be approximated, shedding light on the impact of initial tumor burden on therapeutic outcomes. In conclusion, the study provides insights into CAR T-cell therapy dynamics for acute lymphoblastic leukemias, offering a theoretical foundation for clinical observations and suggesting avenues for future immunotherapy modeling research.

Spatial interactions modulate tumor growth and immune infiltration

Lenia, a cellular automata framework used in artificial life, provides a natural setting to implement mathematical models of cancer incorporating features such as morphogenesis, homeostasis, motility, reproduction, growth, stimuli response, evolvability, and adaptation. Historically, agent-based models of cancer progression have been constructed with rules that govern birth, death and migration, with attempts to map local rules to emergent global growth dynamics. In contrast, Lenia provides a flexible framework for considering a spectrum of local (cell-scale) to global (tumor-scale) dynamics by defining an interaction kernel governing density-dependent growth dynamics. Lenia can recapitulate a range of cancer model classifications including local or global, deterministic or stochastic, non-spatial or spatial, single or multi-population, and off or on-lattice. Lenia is subsequently used to develop data-informed models of 1) single-population growth dynamics, 2) multi-population cell-cell competition models, and 3) cell migration or chemotaxis. Mathematical modeling provides important mechanistic insights. First, short-range interaction kernels provide a mechanism for tumor cell survival under conditions with strong Allee effects. Next, we find that asymmetric interaction tumor-immune kernels lead to poor immune response. Finally, modeling recapitulates immune-ECM interactions where patterns of collagen formation provide immune protection, indicated by an emergent inverse relationship between disease stage and immune coverage.

The role of immune cells in resistance to oncolytic viral therapy

Resistance to treatment poses a major challenge for cancer therapy, and oncoviral treatment encounters the issue of viral resistance as well. In this investigation, we introduce deterministic differential equation models to explore the effect of resistance on oncolytic viral therapy. Specifically, we classify tumor cells into resistant, sensitive, or infected with respect to oncolytic viruses for our analysis. Immune cells can eliminate both tumor cells and viruses. Our research shows that the introduction of immune cells into the tumor-virus interaction prevents all tumor cells from becoming resistant in the absence of conversion from resistance to sensitivity, given that the proliferation rate of immune cells exceeds their death rate. The inclusion of immune cells leads to an additional virus-free equilibrium when the immune cell recruitment rate is sufficiently high. The total tumor burden at this virus-free equilibrium is smaller than that at the virus-free and immune-free equilibrium. Therefore, immune cells are capable of reducing the tumor load under the condition of sufficient immune strength. Numerical investigations reveal that the virus transmission rate and parameters related to the immune response significantly impact treatment outcomes. However, monotherapy alone is insufficient for eradicating tumor cells, necessitating the implementation of additional therapies. Further numerical simulation shows that combination therapy with chimeric antigen receptor (CAR T-cell) therapy can enhance the success of treatment.

Spatial interactions modulate tumor growth and immune infiltration

Direct observation of immune cell trafficking patterns and tumor-immune interactions is unlikely in human tumors with currently available technology, but computational simulations based on clinical data can provide insight to test hypotheses. It is hypothesized that patterns of collagen formation evolve as a mechanism of immune escape, but the exact nature of the interaction between immune cells and collagen is poorly understood. Spatial data quantifying the degree of collagen fiber alignment in squamous cell carcinomas indicates that late stage disease is associated with highly aligned fibers. Here, we introduce a computational modeling framework (called Lenia) to discriminate between two hypotheses: immune cell migration that moves 1) parallel or 2) perpendicular to collagen fiber orientation. The modeling recapitulates immune-ECM interactions where collagen patterns provide immune protection, leading to an emergent inverse relationship between disease stage and immune coverage. We also illustrate the capabilities of Lenia to model the evolution of tumor progression and immune predation. Lenia provides a flexible framework for considering a spectrum of local (cell-scale) to global (tumor-scale) dynamics by defining a kernel cell-cell interaction function that governs tumor growth dynamics under immune predation with immune cell migration. Mathematical modeling provides important mechanistic insights into cell interactions. Short-range interaction kernels provide a mechanism for tumor cell survival under conditions with strong Allee effects, while asymmetric tumor-immune interaction kernels lead to poor immune response. Thus, the length scale of tumor-immune interactions drives tumor growth and infiltration.

Onco-Breastomics: An Eco-Evo-Devo Holistic Approach

Known as a diverse collection of neoplastic diseases, breast cancer (BC) can be hyperbolically characterized as a dynamic pseudo-organ, a living organism able to build a complex, open, hierarchically organized, self-sustainable, and self-renewable tumor system, a population, a species, a local community, a biocenosis, or an evolving dynamical ecosystem (i.e., immune or metabolic ecosystem) that emphasizes both developmental continuity and spatio-temporal change. Moreover, a cancer cell community, also known as an oncobiota, has been described as non-sexually reproducing species, as well as a migratory or invasive species that expresses intelligent behavior, or an endangered or parasite species that fights to survive, to optimize its features inside the host's ecosystem, or that is able to exploit or to disrupt its host circadian cycle for improving the own proliferation and spreading. BC tumorigenesis has also been compared with the early embryo and placenta development that may suggest new strategies for research and therapy. Furthermore, BC has also been characterized as an environmental disease or as an ecological disorder. Many mechanisms of cancer progression have been explained by principles of ecology, developmental biology, and evolutionary paradigms. Many authors have discussed ecological, developmental, and evolutionary strategies for more successful anti-cancer therapies, or for understanding the ecological, developmental, and evolutionary bases of BC exploitable vulnerabilities. Herein, we used the integrated framework of three well known ecological theories: the Bronfenbrenner's theory of human development, the Vannote's River Continuum Concept (RCC), and the Ecological Evolutionary Developmental Biology (Eco-Evo-Devo) theory, to explain and understand several eco-evo-devo-based principles that govern BC progression. Multi-omics fields, taken together as onco-breastomics, offer better opportunities to integrate, analyze, and interpret large amounts of complex heterogeneous data, such as various and big-omics data obtained by multiple investigative modalities, for understanding the eco-evo-devo-based principles that drive BC progression and treatment. These integrative eco-evo-devo theories can help clinicians better diagnose and treat BC, for example, by using non-invasive biomarkers in liquid-biopsies that have emerged from integrated omics-based data that accurately reflect the biomolecular landscape of the primary tumor in order to avoid mutilating preventive surgery, like bilateral mastectomy. From the perspective of preventive, personalized, and participatory medicine, these hypotheses may help patients to think about this disease as a process governed by natural rules, to understand the possible causes of the disease, and to gain control on their own health.

Reuniting philosophy and science to advance cancer research

Cancers rely on multiple, heterogeneous processes at different scales, pertaining to many biomedical fields. Therefore, understanding cancer is necessarily an interdisciplinary task that requires placing specialised experimental and clinical research into a broader conceptual, theoretical, and methodological framework. Without such a framework, oncology will collect piecemeal results, with scant dialogue between the different scientific communities studying cancer. We argue that one important way forward in service of a more successful dialogue is through greater integration of applied sciences (experimental and clinical) with conceptual and theoretical approaches, informed by philosophical methods. By way of illustration, we explore six central themes: (i) the role of mutations in cancer; (ii) the clonal evolution of cancer cells; (iii) the relationship between cancer and multicellularity; (iv) the tumour microenvironment; (v) the immune system; and (vi) stem cells. In each case, we examine open questions in the scientific literature through a philosophical methodology and show the benefit of such a synergy for the scientific and medical understanding of cancer.

Modeling the Role of Immune Cell Conversion in the Tumor-Immune Microenvironment

Tumors develop in a complex physical, biochemical, and cellular milieu, referred to as the tumor microenvironment. Of special interest is the set of immune cells that reciprocally interact with the tumor, the tumor-immune microenvironment (TIME). The diversity of cell types and cell-cell interactions in the TIME has led researchers to apply concepts from ecology to describe the dynamics. However, while tumor cells are known to induce immune cells to switch from anti-tumor to pro-tumor phenotypes, this type of ecological interaction has been largely overlooked. To address this gap in cancer modeling, we develop a minimal, ecological model of the TIME with immune cell conversion, to highlight this important interaction and explore its consequences. A key finding is that immune conversion increases the range of parameters supporting a co-existence phase in which the immune system and the tumor reach a stalemate. Our results suggest that further investigation of the consequences of immune cell conversion, using detailed, data-driven models, will be critical for greater understanding of TIME dynamics.

Data driven model discovery and interpretation for CAR T-cell killing using sparse identification and latent variables

In the development of cell-based cancer therapies, quantitative mathematical models of cellular interactions are instrumental in understanding treatment efficacy. Efforts to validate and interpret mathematical models of cancer cell growth and death hinge first on proposing a precise mathematical model, then analyzing experimental data in the context of the chosen model. In this work, we present the first application of the sparse identification of non-linear dynamics (SINDy) algorithm to a real biological system in order discover cell-cell interaction dynamics in in vitro experimental data, using chimeric antigen receptor (CAR) T-cells and patient-derived glioblastoma cells. By combining the techniques of latent variable analysis and SINDy, we infer key aspects of the interaction dynamics of CAR T-cell populations and cancer. Importantly, we show how the model terms can be interpreted biologically in relation to different CAR T-cell functional responses, single or double CAR T-cell-cancer cell binding models, and density-dependent growth dynamics in either of the CAR T-cell or cancer cell populations. We show how this data-driven model-discovery based approach provides unique insight into CAR T-cell dynamics when compared to an established model-first approach. These results demonstrate the potential for SINDy to improve the implementation and efficacy of CAR T-cell therapy in the clinic through an improved understanding of CAR T-cell dynamics.

Opportunities for improving brain cancer treatment outcomes through imaging-based mathematical modeling of the delivery of radiotherapy and immunotherapy

Immunotherapy has become a fourth pillar in the treatment of brain tumors and, when combined with radiation therapy, may improve patient outcomes and reduce the neurotoxicity. As with other combination therapies, the identification of a treatment schedule that maximizes the synergistic effect of radiation- and immune-therapy is a fundamental challenge. Mechanism-based mathematical modeling is one promising approach to systematically investigate therapeutic combinations to maximize positive outcomes within a rigorous framework. However, successful clinical translation of model-generated combinations of treatment requires patient-specific data to allow the models to be meaningfully initialized and parameterized. Quantitative imaging techniques have emerged as a promising source of high quality, spatially and temporally resolved data for the development and validation of mathematical models. In this review, we will present approaches to personalize mechanism-based modeling frameworks with patient data, and then discuss how these techniques could be leveraged to improve brain cancer outcomes through patient-specific modeling and optimization of treatment strategies.

Cellular kinetics: A clinical and computational review of CAR-T cell pharmacology

To the extent that pharmacokinetics influence the effectiveness of nonliving therapeutics, so too do cellular kinetics influence the efficacy of Chimeric Antigen Receptor (CAR) -T cell therapy. Like conventional therapeutics, CAR-T cell therapies undergo a distribution phase upon administration. Unlike other therapeutics, however, this distribution phase is followed by subsequent phases of expansion, contraction, and persistence. The magnitude and duration of these phases unequivocally influence clinical outcomes. Furthermore, the "pharmacodynamics" of CAR-T cells is truly dynamic, as cells can rapidly become exhausted and lose their therapeutic efficacy. Mathematical models are among the translational tools commonly applied to assess, characterize, and predict the complex cellular kinetics and dynamics of CAR-T cells. Here, we provide a focused review of the cellular kinetics of CAR-T cells, the mechanisms underpinning their complexity, and the mathematical modeling approaches used to interrogate them.

Tumour immunotherapy: lessons from predator-prey theory

With the burgeoning use of immune-based treatments for cancer, never has there been a greater need to understand the tumour microenvironment within which immune cells function and how it can be perturbed to inhibit tumour growth. Yet, current challenges in identifying optimal combinations of immunotherapies and engineering new cell-based therapies highlight the limitations of conventional paradigms for the study of the tumour microenvironment. Ecology has a rich history of studying predator-prey dynamics to discern factors that drive prey to extinction. Here, we describe the basic tenets of predator-prey theory as applied to 'predation' by immune cells and the 'extinction' of cancer cells. Our synthesis reveals fundamental mechanisms by which antitumour immunity might fail in sometimes counterintuitive ways and provides a fresh yet evidence-based framework to better understand and therapeutically target the immune-cancer interface.

Dose-dependent thresholds of dexamethasone destabilize CAR T-cell treatment efficacy

Chimeric antigen receptor (CAR) T-cell therapy is potentially an effective targeted immunotherapy for glioblastoma, yet there is presently little known about the efficacy of CAR T-cell treatment when combined with the widely used anti-inflammatory and immunosuppressant glucocorticoid, dexamethasone. Here we present a mathematical model-based analysis of three patient-derived glioblastoma cell lines treated in vitro with CAR T-cells and dexamethasone. Advanced in vitro experimental cell killing assay technologies allow for highly resolved temporal dynamics of tumor cells treated with CAR T-cells and dexamethasone, making this a valuable model system for studying the rich dynamics of nonlinear biological processes with translational applications. We model the system as a nonautonomous, two-species predator-prey interaction of tumor cells and CAR T-cells, with explicit time-dependence in the clearance rate of dexamethasone. Using time as a bifurcation parameter, we show that (1) dexamethasone destabilizes coexistence equilibria between CAR T-cells and tumor cells in a dose-dependent manner and (2) as dexamethasone is cleared from the system, a stable coexistence equilibrium returns in the form of a Hopf bifurcation. With the model fit to experimental data, we demonstrate that high concentrations of dexamethasone antagonizes CAR T-cell efficacy by exhausting, or reducing the activity of CAR T-cells, and by promoting tumor cell growth. Finally, we identify a critical threshold in the ratio of CAR T-cell death to CAR T-cell proliferation rates that predicts eventual treatment success or failure that may be used to guide the dose and timing of CAR T-cell therapy in the presence of dexamethasone in patients.

Bioengineering and gene-editing technologies have paved the way for advance immunotherapies that can target patient-specific tumor cells. One of these therapies, chimeric antigen receptor (CAR) T-cell therapy has recently shown promise in treating glioblastoma, an aggressive brain cancer often with poor patient prognosis. Dexamethasone is a commonly prescribed anti-inflammatory medication due to the health complications of tumor associated swelling in the brain. However, the immunosuppressant effects of dexamethasone on the immunotherapeutic CAR T-cells are not well understood. To address this issue, we use mathematical modeling to study in vitro dynamics of dexamethasone and CAR T-cells in three patient-derived glioblastoma cell lines. We find that in each cell line studied there is a threshold of tolerable dexamethasone concentration. Below this threshold, CAR T-cells are successful at eliminating the cancer cells, while above this threshold, dexamethasone critically inhibits CAR T-cell efficacy. Our modeling suggests that in the presence of high dexamethasone reduced CAR T-cell efficacy, or increased exhaustion, can occur and result in CAR T-cell treatment failure.