Mixed immunotherapy and chemotherapy of tumors: modeling, applications and biological interpretations

We develop and analyze a mathematical model, in the form of a system of ordinary differential equations (ODEs), governing cancer growth on a cell population level with combination immune, vaccine and chemotherapy treatments. We characterize the ODE system dynamics by locating equilibrium points, determining stability properties, performing a bifurcation analysis, and identifying basins of attraction. These system characteristics are useful not only to gain a broad understanding of the specific system dynamics, but also to help guide the development of combination therapies. Numerical simulations of mixed chemo-immuno and vaccine therapy using both mouse and human parameters are presented. We illustrate situations for which neither chemotherapy nor immunotherapy alone are sufficient to control tumor growth, but in combination the therapies are able to eliminate the entire tumor.