Mathematical modeling of intratumoral immunotherapy yields strategies to improve the treatment outcomes.

Intratumoral injection of immunotherapy aims to maximize its activity within the tumor. However, cytokines are cleared via tumor vessels and escape from the tumor periphery into the host-tissue, reducing efficacy and causing toxicity. Thus, understanding the determinants of the tumor and immune response to intratumoral immunotherapy should lead to better treatment outcomes. In this study, we developed a mechanistic mathematical model to determine the efficacy of intratumorally-injected conjugated-cytokines, accounting for properties of the tumor microenvironment and the conjugated-cytokines. The model explicitly incorporates i) the tumor vascular density and permeability and the tumor hydraulic conductivity, ii) conjugated-cytokines size and binding affinity as well as their clearance via the blood vessels and the surrounding tissue, and iii) immune cells-cancer cells interactions. Model simulations show how the properties of the tumor and of the conjugated-cytokines determine treatment outcomes and how selection of proper parameters can optimize therapy. A high tumor tissue hydraulic permeability allows for the uniform distribution of the cytokines into the tumor, whereas uniform tumor perfusion is required for sufficient access and activation of immune cells. The permeability of the tumor vessels affects the blood clearance of the cytokines and optimal values depend on the size of the conjugates. A size >5 nm in radius was found to be optimal, whereas the binding of conjugates should be high enough to prevent clearance from the tumor into the surrounding tissue. In conclusion, development of strategies to improve vessel perfusion and tissue hydraulic conductivity by reprogramming the microenvironment along with optimal design of conjugated-cytokines can enhance intratumoral immunotherapy.