A reaction-diffusion delay model for spatial tumor virotherapy with optimal control.

Abstract

Oncolytic virotherapy uses viruses that selectively infect and lyse tumor cells while promoting antitumor immunity. Because viral spread, tumor growth, and immune-cell migration are spatially heterogeneous, purely temporal models may miss important treatment dynamics. We develop a reaction-diffusion delay model for spatial tumor virotherapy with a virus-induced cytotoxic T-lymphocyte (CTL) response, thereby incorporating tumor growth, infection, lysis, immune-mediated killing, diffusion, and delayed immune activation. We prove positivity, local well-posedness, and global boundedness under explicit sufficient conditions, and analyze spatially homogeneous equilibria and modal stability. Then, we extend the model with a numerical optimal-control formulation for viral administration and immune stimulation, thus minimizing tumor burden, treatment cost, and excessive CTL proliferation. Simulations indicate that tumor-virus-immune dynamics can remain spatially heterogeneous and that the computed optimal control schedules improve tumor suppression for the chosen parameter set.