A hybrid discrete-continuum modelling approach to explore the impact of T-cell infiltration on anti-tumour immune response

Luis Almeida,Chloe Audebert,Emma Leschiera,Tommaso Lorenzi

DOI: https://doi.org/10.48550/arXiv.2207.09110

2022-10-07

Abstract:We present a spatial hybrid discrete-continuum modelling framework for the interaction dynamics between tumour cells and cytotoxic T cells, which play a pivotal role in the immune response against tumours. In this framework, tumour cells and T cells are modelled as individual agents while chemokines that drive the chemotactic movement of T cells towards the tumour are modelled as a continuum. We formally derive the continuum counterpart of this model, which is given by a coupled system that comprises an integro-differential equation for the density of tumour cells, a partial differential equation for the density of T cells, and a partial differential equation for the concentration of chemokines. We report on computational results of the hybrid model and show that there is an excellent quantitative agreement between them and numerical solutions of the corresponding continuum model. These results shed light on the mechanisms that underlie the emergence of different levels of infiltration of T cells into the tumour and elucidate how T-cell infiltration shapes anti-tumour immune response. Moreover, to present a proof of concept for the idea that, exploiting the computational efficiency of the continuum model, extensive numerical simulations could be carried out, we investigate the impact of T-cell infiltration on the response of tumour cells to different types of anti-cancer immunotherapy.

Analysis of PDEs

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is to explore the impact of T - cell infiltration on the tumor immune response by constructing a hybrid discrete - continuous modeling framework. Specifically, the paper focuses on:

1. **Understanding the mechanism of T - cell infiltration**: By establishing a model that includes the interaction kinetics between tumor cells and cytotoxic T - cells, the paper explores how different levels of T - cell infiltration are formed. This helps to understand how T - cells move in the tumor microenvironment and affect tumor development.

2. **Quantifying the impact of T - cell infiltration**: Through numerical simulation, the paper quantitatively analyzes the specific impact of T - cell infiltration on the tumor immune response. These results help to reveal how T - cell infiltration shapes the anti - tumor immune response.

3. **Evaluating the effectiveness of immunotherapy**: The paper also explores the impact of T - cell infiltration on the effectiveness of different types of anti - cancer immunotherapies. By simulating tumor cell dynamics under different parameter settings, the paper aims to identify possible targets to improve the efficacy of anti - cancer treatment.

### Model overview

The model proposed in the paper is a hybrid discrete - continuous model, in which:

- **Tumor cells and T - cells** are modeled as individual agents, whose behaviors follow a series of rules, such as proliferation, death, and migration.
- **Chemokines** (such as CXCL9/10) are modeled as continuous variables, which describe the concentration of chemokines secreted by tumor cells. These chemokines drive T - cells to migrate towards the tumor area.

### Mathematical model

The core equations of the model are as follows:

1. **Evolution of tumor cell density**:
\[
\frac{\partial n}{\partial t}=\alpha_n n-\mu_n \rho_n(t) n-\zeta_n K(x, t) n
\]
where:
- \(\alpha_n\) is the proliferation rate of tumor cells.
- \(\mu_n \rho_n(t)\) is the mortality rate due to internal competition within the tumor.
- \(\zeta_n K(x, t)\) is the mortality rate due to the cytotoxic effect of T - cells.

2. **Evolution of T - cell density**:
\[
\frac{\partial c}{\partial t}-\nabla\cdot\left[\beta_c \psi(w)\nabla c-\gamma_c \psi(w) c\nabla\phi-\beta_c c\psi'(w)\nabla w\right]=-\mu_c \rho_c(t) c+\alpha_c r(x, t)
\]
where:
- \(\beta_c\) is the diffusion coefficient of T - cells.
- \(\gamma_c\) is the chemotactic sensitivity of T - cells.
- \(\psi(w)\) is the volume - filling effect function, which describes the reduction in cell migration ability at high cell densities.
- \(\mu_c \rho_c(t)\) is the mortality rate due to the homeostatic regulation mechanism.
- \(\alpha_c r(x, t)\) is the inflow rate of T - cells entering from blood vessels.

3. **Evolution of chemokine concentration**:
\[
\frac{\partial\phi}{\partial t}-\beta_\phi\Delta\phi=\alpha_\phi n-\kappa_\phi\phi
\]
where:
- \(\beta_\phi\) is the diffusion coefficient of chemokines.
- \(\alpha_\phi\) is the rate at which tumor cells secrete chemokines.
- \(\kappa_\phi\) is the natural decay rate of chemokines.

### Biological significance

Through these models, the paper provides the following biological insights:

- **Mechanism of T - cell infiltration**: The model reveals how T - cells migrate to the tumor area guided by chemokines and play a role in the tumor microenvironment.
- **Regulation of the immune response**