Publications

Journal Articles


Mathematical Modeling of Tumor-Immune Interactions: Methods, Applications, and Future Perspectives

Published in CSIAM Transactions on Life Sciences, 2025

Mathematical oncology is a rapidly evolving interdisciplinary field that uses mathematical models to enhance our understanding of cancer dynamics, including tumor growth, metastasis, and treatment response. Tumor-immune interactions play a crucial role in cancer biology, influencing tumor progression and the effectiveness of immunotherapy and targeted treatments. However, studying tumor dynamics in isolation often fails to capture the complex interplay between cancer cells and the immune system, which is critical to disease progression and therapeutic efficacy. Mathematical models that incorporate tumor-immune interactions offer valuable insights into these processes, providing a framework for analyzing immune escape, treatment response, and resistance mechanisms. In this review, we provide an overview of mathematical models that describe tumor-immune dynamics, highlighting their applications in understanding tumor growth, evaluating treatment strategies, and predicting immune responses. We also discuss the strengths and limitations of current modeling approaches and propose future directions for the development of more comprehensive and predictive models of tumor-immune interactions. We aim to offer a comprehensive guide to the state of mathematical modeling in tumor immunology, emphasizing its potential to inform clinical decision-making and improve cancer therapies.

Recommended citation: C Li, and J Lei. Mathematical Modeling of Tumor-Immune Interactions: Methods, Applications, and Future Perspectives. CSIAM Transactions on Life Sciences. 2025. doi:10.1038/s41540-025-00513-1
Download Paper

Quantitative cancer-immunity cycle modeling for predicting disease progression in advanced metastatic colorectal cancer

Published in npj Systems Biology and Applications, 2025

Patients with advanced metastatic colorectal cancer (mCRC) typically exhibit significant interindividual differences in treatment responses and face poor survival outcomes. To systematically analyze the heterogeneous tumor progression and recurrence observed in advanced mCRC patients, we developed a quantitative cancer-immunity cycle (QCIC) model. The QCIC model employs differential equations to capture the biological mechanisms underlying the cancer-immunity cycle and predicts tumor evolution dynamics under various treatment strategies through stochastic computational methods. We introduce the treatment response index (TRI) to quantify disease progression in virtual clinical trials and the death probability function (DPF) to estimate overall survival. Additionally, we investigate the impact of predictive biomarkers on survival prognosis in advanced mCRC patients, identifying tumor-infiltrating CD8+ cytotoxic T lymphocytes (CTLs) as key predictors of disease progression and the tumor-infiltrating CD4+ Th1/Treg ratio as a significant determinant of survival outcomes. This study presents an approach that bridges the gap between diverse clinical data sources and the generation of virtual patient cohorts, providing valuable insights into interindividual treatment variability and survival forecasting in mCRC patients.

Recommended citation: Li C, Wei Y, Lei J. Quantitative cancer-immunity cycle modeling for predicting disease progression in advanced metastatic colorectal cancer. NPJ Syst Biol Appl. 2025. doi:10.1038/s41540-025-00513-1
Download Paper

Mathematical Modeling of Tumor Immune Interactions: The Role of Anti-FGFR and Anti-PD-1 in the Combination Therapy

Published in Bulletin of Mathematical Biology, 2024

Bladder cancer poses a significant global health burden with high incidence and recurrence rates. This study addresses the therapeutic challenges in advanced bladder cancer, focusing on the competitive mechanisms of ligand or drug binding to receptors. We developed a refined mathematical model that integrates the dynamics of tumor cells and immune responses, particularly targeting fibroblast growth factor receptor 3 (FGFR3) and immune checkpoint inhibitors (ICIs). This study contributes to understanding combination therapies by elucidating the competitive binding dynamics and quantifying the synergistic effects. The findings highlight the importance of personalized immunotherapeutic strategies, considering factors such as drug dosage, dosing schedules, and patient-specific parameters. Our model further reveals that ligand-independent activated-state receptors are the most essential drivers of tumor proliferation. Moreover, we found that PD-L1 expression rate was more important than PD-1 in driving the dynamic evolution of tumor and immune cells. The proposed mathematical model provides a comprehensive framework for unraveling the complexities of combination therapies in advanced bladder cancer. As research progresses, this multidisciplinary approach contributes valuable insights toward optimizing therapeutic strategies and advancing cancer treatment paradigms.

Recommended citation: Li, C., Ren, Z., Yang, G. et al. Mathematical Modeling of Tumor Immune Interactions: The Role of Anti-FGFR and Anti-PD-1 in the Combination Therapy. Bull Math Biol. 2024. doi.org/10.1007/s11538-024-01329-6.
Download Paper