Modeling tumor growth: a simple individual-based model and its analysis

Yuri Kozitsky,Krzysztof Pilorz

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

2020-03-20

Abstract:Initiation and development of a malignant tumor is a complex phenomenon that has critical stages determining its long time behavior. This phenomenon is mathematically described by means of various models: from simple heuristic models to those employing stochastic processes. In this chapter, we discuss some aspects of such modeling by analyzing a simple individual-based model, in which tumor cells are presented as point particles drifting in $\mathbf{R}_{+}:=[0,+\infty)$ towards the origin with unit speed. At the origin, each of them splits into two new particles that instantly appear in $\mathbf{R}_{+}$ at random positions. During their drift the particles are subject to a random death before splitting. In this model, trait $x\in \mathbf{R}_{+}$ of a given cell corresponds to time to its division and the death is caused by therapeutic factors. On its base we demonstrate how to derive a condition -- involving the therapy related death rate and cell cycle distribution parameters -- under which the tumor size remains bounded in time, which practically means combating the disease.

Populations and Evolution,Dynamical Systems

What problem does this paper attempt to address?

This paper aims to explore the mathematical modeling methods of tumor growth by analyzing a simple individual - based model. Specifically, in this model, tumor cells are regarded as point particles drifting towards the origin at a unit speed on \( \mathbb{R}^+ := [0, +\infty) \). When each particle reaches the origin, it will split into two new particles, and these new particles will then appear at random positions on \( \mathbb{R}^+ \). During the drifting process, particles may die randomly due to treatment factors and will not split.

### Core problems of the model
1. **Describe the dynamic behavior of tumor cells**:
- The life cycle of tumor cells (the time from birth to division) is random.
- Cells may die due to treatment factors before division.
- Daughter cells after division continue to drift and divide in the same way.

2. **Determine the behavior of tumor size changing over time**:
- By analyzing the model, derive a condition that involves treatment - related mortality rates and cell - cycle distribution parameters, so that the tumor size remains bounded in the long term.
- This condition actually means controlling or eliminating the disease.

### Key points of mathematical modeling
- **Individual - based model**: The life cycle and mortality rate of each cell are explicitly considered.
- **Stochastic process**: Cell division and death are random events, following specific probability distributions.
- **Treatment influence**: The cell mortality rate caused by treatment is a key parameter, which determines whether the tumor can be effectively controlled.

### Formulas and mathematical tools
- **Branching parameter**:
\[
a(\tau) = 2e^{-m\tau}
\]
where \( \tau \) is the life cycle of the cell, and \( m \) is the mortality rate caused by treatment.

- **Average branching parameter**:
\[
\langle a \rangle = 2 \int_0^{+\infty} e^{-m\tau} g(\tau) d\tau = 2 \hat{g}(m)
\]
where \( g(\tau) \) is the distribution density function of the cell life cycle, and \( \hat{g}(m) \) is its Laplace transform.

- **Extinction condition**:
\[
\hat{g}(m) < \frac{1}{2}
\]

- **Key parameter**:
\[
m^* = \frac{1}{\theta} \left( 2^{1/k} - 1 \right)
\]
where \( k \) and \( \theta \) are parameters of the cell life cycle distribution.

### Conclusion
This paper explores how to describe and control tumor growth by mathematical methods through establishing and analyzing a simple individual - based model. The key to the model is to determine the cell mortality rate caused by treatment, making it sufficient to keep the tumor size bounded in the long term, thereby achieving disease control. This research provides a new perspective for understanding the complex mechanisms of tumor growth and theoretical support for developing more effective treatment strategies.