Personalized Immunotherapy via Multiscale Tumor-Immune Modeling and Optimal Control

Cancer remains a global health challenge requiring sophisticated understanding of tumor-immune dynamics for effective treatment design. Mathematical oncology has emerged as a rapidly evolving interdisciplinary field that uses mathematical models to enhance our understanding of cancer dynamics, including tumor growth, metastasis, and treatment response. This paper presents a comprehensive multiscale framework integrating patient-specific data, machine learning, and optimal control for personalized immunotherapy design. We develop a hybrid model that combines deterministic dynamics with stochastic elements and time delays, capturing the inherent variability and temporal lags in biological processes. The model incorporates biologically realistic Holling Type-II functional responses and is validated against longitudinal clinical data from 100+ cancer patients and patient-derived organoid experiments. Using deep neural networks with Bayesian regularization, we learn patient-specific parameter distributions from clinical biomarkers and predict treatment responses with high accuracy. Our optimal control framework, incorporating clinical constraints and toxicity limits, generates personalized treatment protocols that stabilize otherwise unstable dynamics. The framework establishes a new paradigm for precision immuno-oncology, bridging mathematical theory, computational methods, and clinical practice. Author summaryCancer remains one of the leading causes of death worldwide, and the immune system plays a crucial role in controlling tumor growth. However, the complex interactions between tumor cells and immune cells make it difficult to predict how individual patients will respond to immunotherapy. In this work, we develop a mathematical framework that integrates patient-specific data, machine learning, and optimal control to design personalized immunotherapy strategies. Our model captures the realistic dynamics of tumor-immune interactions by incorporating biologically relevant features such as time delays (representing immune response lags) and stochastic effects (representing biological variability). Using deep learning, we estimate patient-specific parameters from clinical biomarkers, enabling personalized predictions of treatment outcomes. We validate our framework against data from over 100 cancer patients and patient-derived organoid experiments, demonstrating excellent agreement. Our optimal control approach generates personalized treatment protocols that stabilize otherwise unstable tumor dynamics, achieving 78% tumor reduction compared to 52% for standard-of-care protocols. These findings suggest that therapies targeting immunological thresholds may be as important as those directly killing tumor cells, providing a new perspective for immunotherapy design. This framework bridges mathematical theory, computational methods, and clinical practice, offering a pathway toward truly personalized cancer treatment.