Optimal control of combination immunotherapy in a glioblastoma-immune dynamics model

Abstract The immune checkpoint inhibitor anti-PD-1, commonly used in cancer immunotherapy, has not been successful as a monotherapy for the highly aggressive brain cancer glioblastoma. However, when used in conjunction with a CC-chemokine receptor-2 (CCR2) antagonist, anti-PD-1 has shown efficacy in preclinical studies. In this paper, we aim to optimize treatment regimens for this combination immunotherapy using optimal control theory. We extend a treatment-free glioblastoma-immune dynamics ODE model to include interventions with anti-PD-1 and the CCR2 antagonist. An optimized regimen increases the survival of an average mouse from 32 days post-tumor implantation without treatment to 111 days with treatment. We scale this approach to a virtual murine cohort to evaluate mortality and quality of life concerns during treatment, and predict survival, tumor recurrence, or death after treatment. A parameter identifiability analysis identifies five parameters suitable for personalizing treatment within the virtual cohort. Sampling from these five practically identifiable parameters for the virtual murine cohort reveals that personalized, optimized regimens enhance survival: 84% of the virtual mice survive to day 100, compared to 60% survival in a previously studied experimental regimen. Subjects with high tumor growth rates and low T cell kill rates are identified as more likely to die during and after treatment due to their compromised immune systems and more aggressive tumors. Notably, the MDSC death rate emerges as a long-term predictor of either disease-free survival or death. Highlights A mathematical model of glioma-immune dynamics integrates combination immunotherapy. An optimized regimen extends survival in an average virtual mouse by 79 days. Quality of life and survival outcomes were evaluated for a virtual murine cohort. A high death rate of myeloid-derived suppressor cells predicts long-term survival.