Investigating the synergistic effects of immunotherapy and normalization treatment in modulating tumor microenvironment and enhancing treatment efficacy.

We developed a comprehensive mathematical model of cancer immunotherapy that takes into account: i) Immune checkpoint blockers (ICBs) and the interactions between cancer cells and the immune system, ii) characteristics of the tumor microenvironment, such as the tumor hydraulic conductivity, interstitial fluid pressure, and vascular permeability, iii) spatial and temporal variations of the modeled components within the tumor and the surrounding host tissue, iv) the transport of modeled components through the vasculature and between the tumor-host tissue with convection and diffusion, and v) modeling of the tumor draining lymph nodes were the antigen presentation and the development of cytotoxic immune cells take place. Our model successfully reproduced experimental data from various murine tumor types and predicted immune system profiling, which is challenging to achieve experimentally. It showed that combination of ICB therapy and normalization treatments, that aim to improve tumor perfusion, decreases interstitial fluid pressure and increases the concentration of both innate and adaptive immune cells at the tumor center rather than the periphery. Furthermore, using the model, we investigated the impact of modeled components on treatment outcomes. The analysis found that the number of functional vessels inside the tumor region and the ICB dose administered have the largest impact on treatment outcomes.

Dissecting the Impact of the Gut Microbiome on Cancer Immunotherapy.

The gut microbiome has emerged as a key regulator of response to cancer immunotherapy. However, there is a gap in our understanding of the underlying mechanisms by which the microbiome influences immunotherapy. To this end, we developed a mathematical model based on i) gut microbiome data derived from preclinical studies on melanomas after fecal microbiota transplant, ii) mechanistic modeling of antitumor immune response, and iii) robust association analysis of murine and human microbiome profiles with model-predicted immune profiles. Using our model, we could distill the complexity of these murine and human studies on microbiome modulation in terms of just two model parameters: the activation and killing rate constants of immune cells. We further investigated associations between specific bacterial taxonomies and antitumor immunity and immunotherapy efficacy. This model can guide the design of studies to refine and validate mechanistic links between the microbiome and immune system.

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.

Inducing Biomechanical Heterogeneity in Brain Tumor Modeling by MR Elastography: Effects on Tumor Growth, Vascular Density and Delivery of Therapeutics.

The purpose of this study is to develop a methodology that incorporates a more accurate assessment of tissue mechanical properties compared to current mathematical modeling by use of biomechanical data from magnetic resonance elastography. The elastography data were derived from five glioblastoma patients and a healthy subject and used in a model that simulates tumor growth, vascular changes due to mechanical stresses and delivery of therapeutic agents. The model investigates the effect of tumor-specific biomechanical properties on tumor anisotropic growth, vascular density heterogeneity and chemotherapy delivery. The results showed that including elastography data provides a more realistic distribution of the mechanical stresses in the tumor and induces anisotropic tumor growth. Solid stress distribution differs among patients, which, in turn, induces a distinct functional vascular density distribution-owing to the compression of tumor vessels-and intratumoral drug distribution for each patient. In conclusion, incorporating elastography data results in a more accurate calculation of intratumoral mechanical stresses and enables a better mathematical description of subsequent events, such as the heterogeneous development of the tumor vasculature and intrapatient variations in tumor perfusion and delivery of drugs.