A Mathematical Model of Breast Tumor Progression Based on Immune Infiltration

Modeling cancer progression: an integrated workflow extending data-driven kinetic models to bio-mechanical PDE models

Computational modeling of cancer can help unveil dynamics and interactions that are hard to replicate experimentally. Thanks to the advancement in cancer databases and data analysis technologies, these models have become more robust than ever. There are many mathematical models which investigate cancer through different approaches, from sub-cellular to tissue scale, and from treatment to diagnostic points of view. In this study, we lay out a step-by-step methodology for a data-driven mechanistic model of the tumor microenvironment. We discuss data acquisition strategies, data preparation, parameter estimation, and sensitivity analysis techniques. Furthermore, we propose a possible approach to extend mechanistic ordinary differential equation models to PDE models coupled with mechanical growth. The workflow discussed in this article can help understand the complex temporal and spatial interactions between cells and cytokines in the tumor microenvironment and their effect on tumor growth.

A new treatment for breast cancer using a combination of two drugs: AZD9496 and palbociclib

In this study, we examined a mathematical model of breast cancer (BC) treatment that combines an oral oestrogen receptor inhibitor, AZD9496 with Palbociclib, a selective inhibitor of cyclin- dependent kinases CDK4 and CDK6. Treatment is described by analytical functions that enable us to control the dosage and time interval of the treatment, thus personalising the treatment for each patient. Initially, we investigated the effect of each treatment separately, and finally, we investigated the combination of both treatments. By applying numerical simulations, we confirmed that the combination of AZD9496 with palbociclib was the optimal treatment for BC. The dosage of AZD9496 increased and decreased throughout the treatment period, while the intervals were constant between treatments. Palbociclib changed almost cyclically, whereas the time intervals remained constant. To investigate the mathematical model, we applied the singularly perturbed homotopy analysis method, which is a numerical algorithm. The significant advantage of this method is that the mathematical model does not have to contain a small parameter (as is standard in perturbation theory). However, it is possible to artificially introduce a small parameter into the system of equations, making it possible to study the model using asymptotic methods.

Global stability and parameter analysis reinforce therapeutic targets of PD-L1-PD-1 and MDSCs for glioblastoma

Glioblastoma (GBM) is an aggressive primary brain cancer that currently has minimally effective treatments. Like other cancers, immunosuppression by the PD-L1-PD-1 immune checkpoint complex is a prominent axis by which glioma cells evade the immune system. Myeloid-derived suppressor cells (MDSCs), which are recruited to the glioma microenviroment, also contribute to the immunosuppressed GBM microenvironment by suppressing T cell functions. In this paper, we propose a GBM-specific tumor-immune ordinary differential equations model of glioma cells, T cells, and MDSCs to provide theoretical insights into the interactions between these cells. Equilibrium and stability analysis indicates that there are unique tumorous and tumor-free equilibria which are locally stable under certain conditions. Further, the tumor-free equilibrium is globally stable when T cell activation and the tumor kill rate by T cells overcome tumor growth, T cell inhibition by PD-L1-PD-1 and MDSCs, and the T cell death rate. Bifurcation analysis suggests that a treatment plan that includes surgical resection and therapeutics targeting immune suppression caused by the PD-L1-PD1 complex and MDSCs results in the system tending to the tumor-free equilibrium. Using a set of preclinical experimental data, we implement the approximate Bayesian computation (ABC) rejection method to construct probability density distributions that estimate model parameters. These distributions inform an appropriate search curve for global sensitivity analysis using the extended fourier amplitude sensitivity test. Sensitivity results combined with the ABC method suggest that parameter interaction is occurring between the drivers of tumor burden, which are the tumor growth rate and carrying capacity as well as the tumor kill rate by T cells, and the two modeled forms of immunosuppression, PD-L1-PD-1 immune checkpoint and MDSC suppression of T cells. Thus, treatment with an immune checkpoint inhibitor in combination with a therapeutic targeting the inhibitory mechanisms of MDSCs should be explored.

Mathematical and Machine Learning Models of Renal Cell Carcinoma: A Review

This review explores the multifaceted landscape of renal cell carcinoma (RCC) by delving into both mechanistic and machine learning models. While machine learning models leverage patients' gene expression and clinical data through a variety of techniques to predict patients' outcomes, mechanistic models focus on investigating cells' and molecules' interactions within RCC tumors. These interactions are notably centered around immune cells, cytokines, tumor cells, and the development of lung metastases. The insights gained from both machine learning and mechanistic models encompass critical aspects such as signature gene identification, sensitive interactions in the tumors' microenvironments, metastasis development in other organs, and the assessment of survival probabilities. By reviewing the models of RCC, this study aims to shed light on opportunities for the integration of machine learning and mechanistic modeling approaches for treatment optimization and the identification of specific targets, all of which are essential for enhancing patient outcomes.

Global stability and parameter analysis reinforce therapeutic targets of PD-L1-PD-1 and MDSCs for glioblastoma

Glioblastoma (GBM) is an aggressive primary brain cancer that currently has minimally effective treatments. Like other cancers, immunosuppression by the PD-L1-PD-1 immune checkpoint complex is a prominent axis by which glioma cells evade the immune system. Myeloid-derived suppressor cells (MDSCs), which are recruited to the glioma microenviroment, also contribute to the immunosuppressed GBM microenvironment by suppressing T cell functions. In this paper, we propose a GBM-specific tumor-immune ordinary differential equations model of glioma cells, T cells, and MDSCs to provide theoretical insights into the interactions between these cells. Equilibrium and stability analysis indicates that there are unique tumorous and tumor-free equilibria which are locally stable under certain conditions. Further, the tumor-free equilibrium is globally stable when T cell activation and the tumor kill rate by T cells overcome tumor growth, T cell inhibition by PD-L1-PD-1 and MDSCs, and the T cell death rate. Bifurcation analysis suggests that a treatment plan that includes surgical resection and therapeutics targeting immune suppression caused by the PD-L1-PD1 complex and MDSCs results in the system tending to the tumor-free equilibrium. Using a set of preclinical experimental data, we implement the Approximate Bayesian Computation (ABC) rejection method to construct probability density distributions that estimate model parameters. These distributions inform an appropriate search curve for global sensitivity analysis using the extended Fourier Amplitude Sensitivity Test (eFAST). Sensitivity results combined with the ABC method suggest that parameter interaction is occurring between the drivers of tumor burden, which are the tumor growth rate and carrying capacity as well as the tumor kill rate by T cells, and the two modeled forms of immunosuppression, PD-L1-PD-1 immune checkpoint and MDSC suppression of T cells. Thus, treatment with an immune checkpoint inhibitor in combination with a therapeutic targeting the inhibitory mechanisms of MDSCs should be explored.

Investigating the spatial interaction of immune cells in colon cancer

The intricate network of interactions between cells and molecules in the tumor microenvironment creates a heterogeneous ecosystem. The proximity of the cells and molecules to their activators and inhibitors is essential in the progression of tumors. Here, we develop a system of partial differential equations coupled with linear elasticity to investigate the effects of spatial interactions on the tumor microenvironment. We observe interesting cell and cytokine distribution patterns, which are heavily affected by macrophages. We also see that cytotoxic T cells get recruited and suppressed at the site of macrophages. Moreover, we observe that anti-tumor macrophages reorganize the patterns in favor of a more spatially restricted cancer and necrotic core. Furthermore, the adjoint-based sensitivity analysis indicates that the most sensitive model's parameters are directly related to macrophages. The results emphasize the widely acknowledged effect of macrophages in controlling cancer cells population and spatially arranging cells in the tumor microenvironment.

Analytical Models of Intra- and Extratumoral Cell Interactions at Avascular Stage of Growth in the Presence of Targeted Chemotherapy

In this study, we propose a set of nonlinear differential equations to model the dynamic growth of avascular stage tumors, considering nutrient supply from underlying tissue, innate immune response, contact inhibition of cell migration, and interactions with a chemotherapeutic agent. The model has been validated against available experimental data from the literature for tumor growth. We assume that the size of the modeled tumor is already detectable, and it represents all clinically observed existent cell populations; initial conditions are selected accordingly. Numerical results indicate that the tumor size and regression significantly depend on the strength of the host immune system. The effect of chemotherapy is investigated, not only within the malignancy, but also in terms of the responding immune cells and healthy tissue in the vicinity of a tumor.

Patient-Specific Mathematical Model of the Clear Cell Renal Cell Carcinoma Microenvironment

The interactions between cells and molecules in the tumor microenvironment can give insight into the initiation and progression of tumors and their optimal treatment options. In this paper, we developed an ordinary differential equation (ODE) mathematical model of the interaction network of key players in the clear cell renal cell carcinoma (ccRCC) microenvironment. We then performed a global gradient-based sensitivity analysis to investigate the effects of the most sensitive parameters of the model on the number of cancer cells. The results indicate that parameters related to IL-6 have high a impact on cancer cell growth, such that decreasing the level of IL-6 can remarkably slow the tumor's growth.

Exploring approaches for predictive cancer patient digital twins: Opportunities for collaboration and innovation

We are rapidly approaching a future in which cancer patient digital twins will reach their potential to predict cancer prevention, diagnosis, and treatment in individual patients. This will be realized based on advances in high performance computing, computational modeling, and an expanding repertoire of observational data across multiple scales and modalities. In 2020, the US National Cancer Institute, and the US Department of Energy, through a trans-disciplinary research community at the intersection of advanced computing and cancer research, initiated team science collaborative projects to explore the development and implementation of predictive Cancer Patient Digital Twins. Several diverse pilot projects were launched to provide key insights into important features of this emerging landscape and to determine the requirements for the development and adoption of cancer patient digital twins. Projects included exploring approaches to using a large cohort of digital twins to perform deep phenotyping and plan treatments at the individual level, prototyping self-learning digital twin platforms, using adaptive digital twin approaches to monitor treatment response and resistance, developing methods to integrate and fuse data and observations across multiple scales, and personalizing treatment based on cancer type. Collectively these efforts have yielded increased insights into the opportunities and challenges facing cancer patient digital twin approaches and helped define a path forward. Given the rapidly growing interest in patient digital twins, this manuscript provides a valuable early progress report of several CPDT pilot projects commenced in common, their overall aims, early progress, lessons learned and future directions that will increasingly involve the broader research community.

3D brain slice classification and feature extraction using Deformable Hierarchical Heuristic Model

Brain tumors are the most frequently occurring and severe type of cancer, with a life expectancy of only a few months in most advanced stages. As a result, planning the best course of therapy is critical to improve a patient's ability to fight cancer and their quality of life. Various imaging modalities, such as computed tomography (CT), magnetic resonance imaging (MRI) and ultrasound imaging, are commonly employed to assess a brain tumor. This research proposes a novel technique for extracting and classifying tumor features in 3D brain slice images. After input images are processed for noise removal, resizing, and smoothening, features of brain tumor are extracted using Volume of Interest (VOI). The extracted features are then classified using the Deformable Hierarchical Heuristic Model-Deep Deconvolutional Residual Network (DHHM-DDRN) based on surfaces, curves, and geometric patterns. Experimental results show that proposed approach obtained an accuracy of 95%, DSC of 83%, precision of 80%, recall of 85%, and F1 score of 55% for classifying brain cancer features.

A Dynamically Consistent Nonstandard Difference Scheme for a Discrete-Time Immunogenic Tumors Model

This manuscript deals with the qualitative study of certain properties of an immunogenic tumors model. Mainly, we obtain a dynamically consistent discrete-time immunogenic tumors model using a nonstandard difference scheme. The existence of fixed points and their stability are discussed. It is shown that a continuous system experiences Hopf bifurcation at one and only one positive fixed point, whereas its discrete-time counterpart experiences Neimark–Sacker bifurcation at one and only one positive fixed point. It is shown that there is no chance of period-doubling bifurcation in our discrete-time system. Additionally, numerical simulations are carried out in support of our theoretical discussion.

Cancer systems epidemiology: Overcoming misconceptions and integrating systems approaches into cancer research

Patricia Mabry and coauthors discuss application of systems approaches in cancer research.

A PDE Model of Breast Tumor Progression in MMTV-PyMT Mice

The evolution of breast tumors greatly depends on the interaction network among different cell types, including immune cells and cancer cells in the tumor. This study takes advantage of newly collected rich spatio-temporal mouse data to develop a data-driven mathematical model of breast tumors that considers cells' location and key interactions in the tumor. The results show that cancer cells have a minor presence in the area with the most overall immune cells, and the number of activated immune cells in the tumor is depleted over time when there is no influx of immune cells. Interestingly, in the case of the influx of immune cells, the highest concentrations of both T cells and cancer cells are in the boundary of the tumor, as we use the Robin boundary condition to model the influx of immune cells. In other words, the influx of immune cells causes a dominant outward advection for cancer cells. We also investigate the effect of cells' diffusion and immune cells' influx rates in the dynamics of cells in the tumor micro-environment. Sensitivity analyses indicate that cancer cells and adipocytes' diffusion rates are the most sensitive parameters, followed by influx and diffusion rates of cytotoxic T cells, implying that targeting them is a possible treatment strategy for breast cancer.

Investigating key cell types and molecules dynamics in PyMT mice model of breast cancer through a mathematical model

The most common kind of cancer among women is breast cancer. Understanding the tumor microenvironment and the interactions between individual cells and cytokines assists us in arriving at more effective treatments. Here, we develop a data-driven mathematical model to investigate the dynamics of key cell types and cytokines involved in breast cancer development. We use time-course gene expression profiles of a mouse model to estimate the relative abundance of cells and cytokines. We then employ a least-squares optimization method to evaluate the model’s parameters based on the mice data. The resulting dynamics of the cells and cytokines obtained from the optimal set of parameters exhibit a decent agreement between the data and predictions. We perform a sensitivity analysis to identify the crucial parameters of the model and then perform a local bifurcation on them. The results reveal a strong connection between adipocytes, IL6, and the cancer population, suggesting them as potential targets for therapies.

One of the outstanding challenges of the mathematical modeling of cancer progression is the existence of many unknown parameters. In this work, we develop a data-driven mathematical model of breast cancer progression by deriving a system of ordinary differential equations for the interaction networks of key cell types and molecules in breast tumors. To overcome the limitations of unknown parameters, we utilize a time course data of a PyMT mice model of breast cancer and estimate parameters using an optimization method. Although the predicted dynamics of cancer and necrotic cells using the obtained values of parameters are in good agreement with the data, the predicted values for a few other variables do not match the data. This might indicate that there are some other key interactions that have not been modeled, and/or there is a noise in the data. The sensitivity and bifurcation analyses show that the most important parameters in controlling the cancer cells population are the proliferation and death rates of cancer cells and adipocytes. These results are in agreement with some biological and clinical studies of breast cancer, which have reported a link between adipocytes and breast cancer progression.