skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: An In Silico Glioblastoma Microenvironment Model Dissects the Immunological Mechanisms of Resistance to PD‐1 Checkpoint Blockade Immunotherapy
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
Award ID(s):
1701322
PAR ID:
10450762
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Small Methods
Volume:
5
Issue:
6
ISSN:
2366-9608
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Programmed cell death protein-1 (PD-1) checkpoint immunotherapy efficacy remains unpredictable in glioblastoma (GBM) patients due to the genetic heterogeneity and immunosuppressive tumor microenvironments. Here, we report a microfluidics-based, patient-specific ‘GBM-on-a-Chip’ microphysiological system to dissect the heterogeneity of immunosuppressive tumor microenvironments and optimize anti-PD-1 immunotherapy for different GBM subtypes. Our clinical and experimental analyses demonstrated that molecularly distinct GBM subtypes have distinct epigenetic and immune signatures that may lead to different immunosuppressive mechanisms. The real-time analysis in GBM-on-a-Chip showed that mesenchymal GBM niche attracted low number of allogeneic CD154+CD8+ T-cells but abundant CD163+ tumor-associated macrophages (TAMs), and expressed elevated PD-1/PD-L1 immune checkpoints and TGF-β1, IL-10, and CSF-1 cytokines compared to proneural GBM. To enhance PD-1 inhibitor nivolumab efficacy, we co-administered a CSF-1R inhibitor BLZ945 to ablate CD163+ M2-TAMs and strengthened CD154+CD8+ T-cell functionality and GBM apoptosis on-chip. Our ex vivo patient-specific GBM-on-a-Chip provides an avenue for a personalized screening of immunotherapies for GBM patients. 
    more » « less
  2. 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
  3. Immune checkpoint therapies such as PD-1 blockade have vastly improved the treatment of numerous cancers, including basal cell carcinoma (BCC). However, patients afflicted with pancreatic ductal carcinoma (PDAC), one of the deadliest malignancies, overwhelmingly exhibit negative responses to checkpoint therapy. We sought to combine data analysis and machine learning to differentiate the putative mechanisms of BCC and PDAC non-response. We discover that increased MHC-I expression in malignant cells and suppression of MHC and PD-1/PD-L expression in CD8 + T cells is associated with nonresponse to treatment. Furthermore, we leverage machine learning to predict response to PD-1 blockade on a cellular level. We confirm divergent resistance mechanisms between BCC, PDAC, and melanoma and highlight the potential for rapid and affordable testing of gene expression in BCC patients to accurately predict response to checkpoint therapies. Our findings present an optimistic outlook for the use of quantitative cross-cancer analyses in characterizing immune responses and predicting immunotherapy outcomes. 
    more » « less
  4. null (Ed.)
    The advent of immune checkpoint therapy for metastatic skin cancer has greatly improved patient survival. However, most skin cancer patients are refractory to checkpoint therapy, and furthermore, the intra-immune cell signaling driving response to checkpoint therapy remains uncharacterized. When comparing the immune transcriptome in the tumor microenvironment of melanoma and basal cell carcinoma (BCC), we found that the presence of memory B cells and macrophages negatively correlate in both cancers when stratifying patients by their response, with memory B cells more present in responders. Moreover, inhibitory immune signaling mostly decreases in melanoma responders and increases in BCC responders. We further explored the relationships between macrophages, B cells and response to checkpoint therapy by developing a stochastic differential equation model which qualitatively agrees with the data analysis. Our model predicts BCC to be more refractory to checkpoint therapy than melanoma and predicts the best qualitative ratio of memory B cells and macrophages for successful treatment. 
    more » « less
  5. We present a mechanistic mathematical model of immune checkpoint inhibitor therapy to address the oncological need for early, broadly applicable readouts (biomarkers) of patient response to immunotherapy. The model is built upon the complex biological and physical interactions between the immune system and cancer, and is informed using only standard-of-care CT. We have retrospectively applied the model to 245 patients from multiple clinical trials treated with anti–CTLA-4 or anti–PD-1/PD-L1 antibodies. We found that model parameters distinctly identified patients with common ( n = 18) and rare ( n = 10) malignancy types who benefited and did not benefit from these monotherapies with accuracy as high as 88% at first restaging (median 53 days). Further, the parameters successfully differentiated pseudo-progression from true progression, providing previously unidentified insights into the unique biophysical characteristics of pseudo-progression. Our mathematical model offers a clinically relevant tool for personalized oncology and for engineering immunotherapy regimens. 
    more » « less