An official website of the United States government
Abstract 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.
more »
« less
Modeling and Analyzing Homogeneous Tumor Growth under Virotherapy
We present a mathematical model based on ordinary differential equations to investigate the spatially homogeneous state of tumor growth under virotherapy. The model emphasizes the interaction among the tumor cells, the oncolytic viruses, and the host immune system that generates both innate and adaptive immune responses. We conduct a rigorous equilibrium analysis and derive threshold conditions that determine the growth or decay of the tumor under various scenarios. Numerical simulation results verify our analytical predictions and provide additional insight into the tumor growth dynamics.
more »
« less
- Award ID(s):
- 1951345
- PAR ID:
- 10430158
- Date Published:
- Journal Name:
- Mathematics
- Volume:
- 11
- Issue:
- 2
- ISSN:
- 2227-7390
- Page Range / eLocation ID:
- 360
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Myeloid-derived suppressor cells (MDSCs) play a prominent role in the tumor microenvironment. A quantitative understanding of the tumor–MDSC interactions that influence disease progression is critical, and currently lacking. We developed a mathematical model of metastatic growth and progression in immune-rich tumor microenvironments. We modeled the tumor–immune dynamics with stochastic delay differential equations and studied the impact of delays in MDSC activation/recruitment on tumor growth outcomes. In the lung environment, when the circulating level of MDSCs was low, the MDSC delay had a pronounced impact on the probability of new metastatic establishment: blocking MDSC recruitment could reduce the probability of metastasis by as much as 50%. To predict patient-specific MDSC responses we fit to the model individual tumors treated with immune checkpoint inhibitors via Bayesian parameter inference. We reveal that control of the inhibition rate of natural killer (NK) cells by MDSCs had a larger influence on tumor outcomes than controlling the tumor growth rate directly. Posterior classification of tumor outcomes demonstrates that incorporating knowledge of the MDSC responses improved predictive accuracy from 63% to 82%. Investigation of the MDSC dynamics in an environment low in NK cells and abundant in cytotoxic T cells revealed, in contrast, that small MDSC delays no longer impacted metastatic growth dynamics. Our results illustrate the importance of MDSC dynamics in the tumor microenvironment overall and predict interventions promoting shifts toward less immune-suppressed states. We propose that there is a pressing need to consider MDSCs more often in analyses of tumor microenvironments.more » « less
-
Abstract Recent years have seen a tremendous growth of interest in understanding the role that the adaptive immune system could play in interdicting tumor progression. In this context, it has been shown that the density of adaptive immune cells inside a solid tumor serves as a favorable prognostic marker across different types of cancer. The exact mechanisms underlying the degree of immune cell infiltration is largely unknown. Here, we quantify the temporal dynamics of the density profile of activated immune cells around a solid tumor spheroid. We propose a computational model incorporating immune cells with active, persistent movement and a proliferation rate that depends on the presence of cancer cells, and show that the model able to reproduce semi-quantitatively the experimentally measured infiltration profile. Studying the density distribution of immune cells inside a solid tumor can help us better understand immune trafficking in the tumor micro-environment, hopefully leading towards novel immunotherapeutic strategies.more » « less
-
Abstract During progression from carcinoma in situ to an invasive tumor, the immune system is engaged in complex sets of interactions with various tumor cells. Tumor cell plasticity alters disease trajectories via epithelial-to-mesenchymal transition (EMT). Several of the same pathways that regulate EMT are involved in tumor-immune interactions, yet little is known about the mechanisms and consequences of crosstalk between these regulatory processes. Here we introduce a multiscale evolutionary model to describe tumor-immune-EMT interactions and their impact on epithelial cancer progression from in situ to invasive disease. Through simulation of patient cohorts in silico, the model predicts that a controllable region maximizes invasion-free survival. This controllable region depends on properties of the mesenchymal tumor cell phenotype: its growth rate and its immune-evasiveness. In light of the model predictions, we analyze EMT-inflammation-associated data from The Cancer Genome Atlas, and find that association with EMT worsens invasion-free survival probabilities. This result supports the predictions of the model, and leads to the identification of genes that influence outcomes in bladder and uterine cancer, including FGF pathway members. These results suggest new means to delay disease progression, and demonstrate the importance of studying cancer-immune interactions in light of EMT.more » « less
-
Abstract The PD‐1 immune checkpoint‐based therapy has emerged as a promising therapy strategy for treating the malignant brain tumor glioblastoma (GBM). However, patient response varies in clinical trials, mainly due to the tumor heterogeneity and immunological resistance in the tumor microenvironment. To further understand how mechanistically the niche interplay and competition drive anti‐PD‐1 resistance, an in silico model is established to quantitatively describe the biological rationale of critical GBM‐immune interactions, such as tumor growth and apoptosis, T cell activation and cytotoxicity, and tumor‐associated macrophage (TAM) mediated immunosuppression. Such an in silico experimentation and predictive model, based on the in vitro microfluidic chip‐measured end‐point data and patient‐specific immunological characteristics, allows for a comprehensive and dynamic analysis of multiple TAM‐associated immunosuppression mechanisms against the anti‐PD‐1 immunotherapy. The computational model demonstrates that the TAM‐associated immunosuppression varies in severity across different GBM subtypes, which results in distinct tumor responses. The prediction results indicate that a combination therapy by co‐targeting of PD‐1 checkpoint and TAM‐associated CSF‐1R signaling can enhance the immune responses of GBM patients, especially those patients with mesenchymal GBM who are irresponsive to the single anti‐PD‐1 therapy. The development of a patient‐specific in silico–in vitro GBM model will help navigate and personalize immunotherapies for GBM patients.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.3390/math11020360
