Carbonic Anhydrase IX Suppression Shifts Partial Response to Checkpoint Inhibitors into Complete Tumor Eradication: Model-Based Investigation
Abstract
:1. Introduction
2. Results
2.1. Phase Portrait Analysis Reveals Two Possible Steady States
2.2. Synergistic Combination Therapy with Anti-CAIX and Immune Checkpoint Inhibitors
2.3. Model Fitting
3. Discussion
4. Materials and Methods
4.1. Differential Equation Model Formulation
- Tumor growth Tumor growth is logistic with the carrying capacity K. The model differentiates between cancer non-stem cells (denoted in the manuscript as CCs) and cancer stem cells ( denoted as CSCs). CSCs can only be killed by immune cells, whereas CCs experience apoptosis with the rate n [22]. CSCs divide asymmetrically with rate and symmetrically otherwise.
- Tumor–immune interactions T-cells’ infiltration is proportional to the tumor volume, and their number decreases exponentially due to cell death. They attack and kill cancer cells at a rate proportional to their density in the tumor, as proposed in [23]. Notably, this is a spin on the classical Kuznetsov-type interactions as presented in [24], where tumor cell killing is proportional to the product of the number of tumor cells and T-cells. We believe that our modification suits our needs better than the original interaction term. In particular, let us consider the scenario in which we compare two tumors consisting of the same number of CCs and T-cells. Let us assume further that one of the tumors also has a large population of CSCs, while the second has none. If we used the original Kuznetsov-type term, the CC decay due to interactions with T-cells would be the same in both tumors. On the other hand, in our model, the decay of CCs is smaller in the tumor with CSCs, since the T-cell density is smaller in this tumor. This seems more plausible, as the lymphocytes are then more likely to attack CSCs instead of just CCs. However, our model tacitly assumes that the tumor infiltration by lymphocytes is not over-saturated, i.e., the interactions between T-cells and cancer cells are not limited by a lack of cancer cells. In particular, T-cell decay due to interactions with cancer cells depends only on the number of T-cells. Finally, only cancer cells expressing MHC class I on their surface are recognized and attacked by T-cells. Moreover, the immune response is higher for tumors with a higher tumor mutational burden. Therefore, the rate of tumor cell killing is equal to , where a indicates the interaction rate between tumor and immune cells, m the fraction of cancer cells expressing MHC class I, and quantifies TMB, as proposed in [25].
- PD-1-PD-L1 pathway Tumor cell killing by immune cells is inhibited via the binding of PD-1 and PD-L1, which induces T-cell anergy. We assume that the fraction of PD-1-expressing cells is constant and equal to p. The expression of PD-L1, however, can be either constitutive or adaptive, i.e., induced by as a way of escaping the immune response [26]. We assume that the fraction of cancer cells with constitutive PD-L1 expression is constant and equal to . Adaptive PD-L1 expression is dynamic and bounded from above by the parameter .
- Substances in the TME IFN- is produced by active lymphocytes with rate r and decays naturally with rate . Protons are produced due to cancer cell metabolism with rate and due to CAIX expression with rate q. Outside of the tumor, we assume a physiological pH. The flux of protons into and out of the TME is proportional to the difference between the pH in the TME and the physiological pH. Immune cells that are exposed to acidosis die. The lower the pH, the greater the induced death rate.
- Anti-CAIX suppresses CAIX expression by the fraction .
- Anti-PD-1 suppresses PD-1 expression by the fraction .
- Anti-CTLA-4 is mainly responsible for reinvigorating early T-cell activation in the lymph nodes, which we include in our model by increasing lymphocyte influx by the rate .
4.2. Model Calibration
5. Conclusions
Author Contributions
Funding
Institutional Review Board Statement
Informed Consent Statement
Data Availability Statement
Conflicts of Interest
Abbreviations
CAIX | carbonic anhydrase IX |
PD-1 | programmed cell death protein 1 |
CTLA-4 | cytotoxic T-lymphocyte-associated protein 4 |
CC | cancer cells |
CSC | cancer stem cells |
TME | tumor microenvironment |
KO | knock-out |
NSCLC | non-small-cell lung cancer |
References
- Topalian, S.L.; Hodi, F.S.; Brahmer, J.R.; Gettinger, S.N.; Smith, D.C.; McDermott, D.F.; Powderly, J.D.; Carvajal, R.D.; Sosman, J.A.; Atkins, M.B.; et al. Safety, Activity, and Immune Correlates of Anti–PD-1 Antibody in Cancer. N. Engl. J. Med. 2012, 366, 2443–2454. [Google Scholar] [CrossRef] [Green Version]
- Seidel, J.A.; Otsuka, A.; Kabashima, K. Anti-PD-1 and Anti-CTLA-4 Therapies in Cancer: Mechanisms of Action, Efficacy, and Limitations. Front. Oncol. 2018, 8, 86. [Google Scholar] [CrossRef]
- Bagchi, S.; Yuan, R.; Engleman, E.G. Immune Checkpoint Inhibitors for the Treatment of Cancer: Clinical Impact and Mechanisms of Response and Resistance. Annu. Rev. Pathol. Mech. Dis. 2021, 16, 223–249. [Google Scholar] [CrossRef]
- Ciccolini, J.; Benzekry, S.; Barlesi, F. Deciphering the response and resistance to immune-checkpoint inhibitors in lung cancer with artificial intelligence-based analysis: When PIONeeR meets QUANTIC. Br. J. Cancer 2020, 123, 337–338. [Google Scholar] [CrossRef] [PubMed]
- Havel, J.J.; Chowell, D.; Chan, T.A. The evolving landscape of biomarkers for checkpoint inhibitor immunotherapy. Nat. Rev. Cancer 2019, 19, 133–150. [Google Scholar] [CrossRef] [PubMed]
- Brahmer, J.R. Harnessing the immune system for the treatment of non-small-cell lung cancer. J. Clin. Oncol. 2013, 31, 1021–1028. [Google Scholar] [CrossRef] [PubMed]
- Herzberg, B.; Campo, M.J.; Gainor, J.F. Immune Checkpoint Inhibitors in Non-Small Cell Lung Cancer. Oncologist 2017, 22, 81–88. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Corbet, C.; Feron, O. Tumour acidosis: From the passenger to the driver’s seat. Nat. Rev. Cancer 2017, 17, 577–593. [Google Scholar] [CrossRef] [PubMed]
- Lee, S.H.; McIntyre, D.; Honess, D.; Hulikova, A.; Pacheco-Torres, J.; Cerdán, S.; Swietach, P.; Harris, A.L.; Griffiths, J.R. Carbonic anhydrase IX is a pH-stat that sets an acidic tumour extracellular pH in vivo. Br. J. Cancer 2018, 119, 622–630. [Google Scholar] [CrossRef] [Green Version]
- Chiche, J.; Ilc, K.; Laferrière, J.; Trottier, E.; Dayan, F.; Mazure, N.M.; Brahimi-Horn, M.C.; Pouysségur, J. Hypoxia-Inducible Carbonic Anhydrase IX and XII Promote Tumor Cell Growth by Counteracting Acidosis through the Regulation of the Intracellular pH. Cancer Res. 2009, 69, 358–368. [Google Scholar] [CrossRef] [Green Version]
- Giatromanolaki, A.; Koukourakis, M.I.; Sivridis, E.; Pastorek, J.; Wykoff, C.C.; Gatter, K.C.; Harris, A.L. Expression of Hypoxia-inducible Carbonic Anhydrase-9 Relates to Angiogenic Pathways and Independently to Poor Outcome in Non-Small Cell Lung Cancer. Cancer Res. 2001, 61, 7992–7998. Available online: http://xxx.lanl.gov/abs/https://aacrjournals.org/cancerres/article-pdf/61/21/7992/2489028/ch2101007992.pdf (accessed on 22 April 2023).
- Neri, D.; Supuran, C.T. Interfering with pH regulation in tumours as a therapeutic strategy. Nat. Rev. Drug Discov. 2011, 10, 767–777. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- McDonald, P.C.; Chia, S.; Bedard, P.L.; Chu, Q.; Lyle, M.; Tang, L.; Singh, M.; Zhang, Z.; Supuran, C.T.; Renouf, D.J.; et al. A Phase 1 Study of SLC-0111, a Novel Inhibitor of Carbonic Anhydrase IX, in Patients with Advanced Solid Tumors. Am. J. Clin. Oncol. Cancer Clin. Trials 2020, 43, 484–490. [Google Scholar] [CrossRef] [PubMed]
- De Campos, N.S.P.; de Oliveira Beserra, A.; Pereira, P.H.B.; Chaves, A.S.; Fonseca, F.L.A.; da Silva Medina, T.; dos Santos, T.G.; Wang, Y.; Marasco, W.A.; Suarez, E.R. Immune Checkpoint Blockade via PD-L1 Potentiates More CD28-Based than 4-1BB-Based Anti-Carbonic Anhydrase IX Chimeric Antigen Receptor T Cells. Int. J. Mol. Sci. 2022, 23, 5448. [Google Scholar] [CrossRef]
- Suarez, E.R.; Chang, D.K.; Sun, J.; Sui, J.; Freeman, G.J.; Signoretti, S.; Zhu, Q.; Marasco, W.A.; Suarez, E.R.; Chang, D.K.; et al. Chimeric antigen receptor T cells secreting anti-PD-L1 antibodies more effectively regress renal cell carcinoma in a humanized mouse model. Oncotarget 2016, 7, 34341–34355. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Chafe, S.C.; McDonald, P.C.; Saberi, S.; Nemirovsky, O.; Venkateswaran, G.; Burugu, S.; Gao, D.; Delaidelli, A.; Kyle, A.H.; Baker, J.H.; et al. Targeting Hypoxia-Induced Carbonic Anhydrase IX Enhances Immune-Checkpoint Blockade Locally and Systemically. Cancer Immunol. Res. 2019, 7, 1064–1078. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Grajek, J.; Kather, J.N.; Poleszczuk, J. An in silico model to study the impact of carbonic anhydrase IX expression on tumour growth and anti-PD-1 therapy. J. R. Soc. Interface 2023, 20, 20220654. [Google Scholar] [CrossRef] [PubMed]
- Feng, L.; Dong, Z.; Tao, D.; Zhang, Y.; Liu, Z. The acidic tumor microenvironment: A target for smart cancer nano-theranostics. Natl. Sci. Rev. 2018, 5, 269–286. [Google Scholar] [CrossRef] [Green Version]
- Cortes-Dericks, L.; Galetta, D. Impact of Cancer Stem Cells and Cancer Stem Cell-Driven Drug Resiliency in Lung Tumor: Options in Sight. Cancers 2022, 14, 267. [Google Scholar] [CrossRef]
- Lock, F.E.; McDonald, P.C.; Lou, Y.; Serrano, I.; Chafe, S.C.; Ostlund, C.; Aparicio, S.; Winum, J.Y.; Supuran, C.T.; Dedhar, S. Targeting carbonic anhydrase IX depletes breast cancer stem cells within the hypoxic niche. Oncogene 2012, 32, 5210–5219. [Google Scholar] [CrossRef]
- Kravaris, C.; Hahn, J.; Chu, Y. Advances and selected recent developments in state and parameter estimation. Comput. Chem. Eng. 2013, 51, 111–123. [Google Scholar] [CrossRef]
- Poleszczuk, J.; Macklin, P.; Enderling, H. Agent-Based Modeling of Cancer Stem Cell Driven Solid Tumor Growth. Methods Mol. Biol. 2016, 1516, 335. [Google Scholar] [CrossRef] [PubMed]
- Poleszczuk, J.; Enderling, H. The Optimal Radiation Dose to Induce Robust Systemic Anti-Tumor Immunity. Int. J. Mol. Sci. 2018, 19, 3377. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Kuznetsov, V.A.; Makalkin, I.A.; Taylor, M.A.; Perelson, A.S. Nonlinear dynamics of immunogenic tumors: Parameter estimation and global bifurcation analysis. Bull. Math. Biol. 1994, 56, 295–321. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Valentinuzzi, D.; Simončič, U.; Uršič, K.; Vrankar, M.; Turk, M.; Jeraj, R. Predicting tumour response to anti-PD-1 immunotherapy with computational modelling. Phys. Med. Biol. 2019, 64, 025017. [Google Scholar] [CrossRef] [Green Version]
- Cottrell, T.R.; Taube, J.M. PD-L1 and Emerging Biomarkers in PD-1/PD-L1 Blockade Therapy. Cancer J. 2018, 24, 41. [Google Scholar] [CrossRef]
- Danciu, C.; Falamas, A.; Dehelean, C.; Soica, C.; Radeke, H.; Barbu-Tudoran, L.; Bojin, F.; Pînzaru, S.C.; Munteanu, M.F. A characterization of four B16 murine melanoma cell sublines molecular fingerprint and proliferation behavior. Cancer Cell Int. 2013, 13, 75. [Google Scholar] [CrossRef] [Green Version]
Parameter Name | Value | Lower Bound | Upper Bound |
---|---|---|---|
0 | |||
0 | |||
eta | 0 | ||
n | 0 | 10 | |
q | |||
0 | 1 | ||
0 | 20 | ||
0 | 1 | ||
initCSCrat | 0 | 1 | |
0 | 1 | ||
0 | 1 | ||
0 | 1 | ||
0 | 1 | ||
0 | 1 | ||
0 | 1 |
Par | Interpretation | Value | Unit | Source |
---|---|---|---|---|
maximal rate of tumor cell growth | day | [27] | ||
K | carrying capacity for tumor cells | 1200 | mm | permitted tumor volume limit [16] |
probability of asymmetric division | - | [17] | ||
a | interaction rate between tumor cells and TILs | - | day | free parameter |
volume of one tumor cell | [25] | |||
m | mean MHC class I expression | - | [25] | |
antigenicity strength (single nucleotide variations) | 908 | - | [25] | |
p | mean PD-1 expression by TILs | - | [25] | |
n | tumor cell apoptosis rate | - | day | free parameter |
b | infiltration rate of T-cells into TME | - | free parameter | |
d | apoptosis rate of T-cells | day | [25] | |
rate of T-cell death due to acidosis | - | day | free parameter | |
r | rate of production | [17] | ||
rate of decay | day | [17] | ||
rate of proton production due to tumor cell metabolism | assumption to yield realistic pH values | |||
q | rate of proton production due to CAIX expression | - | free parameter | |
v | rate of proton flux into and out of the TME | 5 | day | assumption to yield realistic pH values |
proton concentration at physiological pH | 3.98 × 1 × | |||
proton concentration equivalent to pH = 6.7 | 2 × 1 × | [17] | ||
constitutive PD-L1 expression | 0.1 | - | [17] | |
rate of adaptive PD-L1 expression | - | assumption, | ||
saturation constant | assumption | |||
effect of anti-PD-1 | - | - | free parameter | |
effect of anti-CTLA-4 | - | - | free parameter | |
effect of anti-CAIX | - | - | free parameter | |
CSCrat | ratio of CSC at inoculation | - | - | free parameter |
inocCells | ratio of inoculated cells that initiates the tumor | - | - | free parameter |
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Grajek, J.; Poleszczuk, J. Carbonic Anhydrase IX Suppression Shifts Partial Response to Checkpoint Inhibitors into Complete Tumor Eradication: Model-Based Investigation. Int. J. Mol. Sci. 2023, 24, 10068. https://doi.org/10.3390/ijms241210068
Grajek J, Poleszczuk J. Carbonic Anhydrase IX Suppression Shifts Partial Response to Checkpoint Inhibitors into Complete Tumor Eradication: Model-Based Investigation. International Journal of Molecular Sciences. 2023; 24(12):10068. https://doi.org/10.3390/ijms241210068
Chicago/Turabian StyleGrajek, Julia, and Jan Poleszczuk. 2023. "Carbonic Anhydrase IX Suppression Shifts Partial Response to Checkpoint Inhibitors into Complete Tumor Eradication: Model-Based Investigation" International Journal of Molecular Sciences 24, no. 12: 10068. https://doi.org/10.3390/ijms241210068
APA StyleGrajek, J., & Poleszczuk, J. (2023). Carbonic Anhydrase IX Suppression Shifts Partial Response to Checkpoint Inhibitors into Complete Tumor Eradication: Model-Based Investigation. International Journal of Molecular Sciences, 24(12), 10068. https://doi.org/10.3390/ijms241210068