On the design and stability of cancer adaptive therapy cycles: deterministic and stochastic models

Yuri G. Vilela,Artur C. Fassoni,Armando G. M. Neves

DOI: https://doi.org/10.1101/2024.09.10.612338

2024-09-14

Abstract:Adaptive therapy is a promising paradigm for treating cancers, that exploits competitive interactions between drug-sensitive and drug-resistant cells, thereby avoiding or delaying treatment failure due to evolution of drug resistance within the tumor. Previous studies have shown the mathematical possibility of building cyclic schemes of drug administration which restore tumor composition to its exact initial value in deterministic models. However, algorithms for cycle design, the conditions on which such algorithms are certain to work, as well as conditions for cycle stability remain elusive. Here, we state biologically motivated hypotheses that guarantee existence of such cycles in two deterministic classes of mathematical models already considered in the literature: Lotka-Volterra and adjusted replicator dynamics. We stress that not only existence of cyclic schemes, but also stability of such cycles is a relevant feature for applications in real clinical scenarios. We also analyze stochastic versions of the above deterministic models, a necessary step if we want to take into account that real tumors are composed by a finite population of cells subject to randomness, a relevant feature in the context of low tumor burden. We argue that the stability of the deterministic cycles is also relevant for the stochastic version of the models. In fact, Dua, Ma and Newton [Cancers (2021)] and Park and Newton [Phys. Rev. E (2023)] observed breakdown of deterministic cycles in a stochastic model (Moran process) for a tumor. Our findings indicate that the breakdown phenomenon is not due to stochasticity itself, but to the deterministic instability inherent in the cycles of the referenced papers. We then illustrate how stable deterministic cycles avoid for very large times the breakdown of cyclic treatments in stochastic tumor models.

Cancer Biology

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to design and stabilize treatment cycles in cancer adaptive therapy. Specifically, the author focuses on how to construct and ensure the existence and stability of these treatment cycles in deterministic mathematical models such as the Lotka - Volterra model and the adjusted replicator dynamics model. In addition, the author also explores the stochastic versions of these deterministic models to take into account the situation where the number of cells in actual tumors is limited and affected by random fluctuations.

### Main research questions:
1. **Design of periodic treatment regimens**:
- How to design periodic treatment regimens in deterministic mathematical models that can restore the tumor composition to its initial state?
- What are the existence conditions of these periodic treatment regimens?

2. **Stability of periodic treatment regimens**:
- How to ensure the stability of these periodic treatment regimens in practical applications?
- What is the importance of stability for clinical applications?

3. **Periodic stability in stochastic models**:
- How to evaluate the stability of periodic treatment regimens in the case of considering random fluctuations?
- Is the periodic collapse phenomenon in stochastic models related to the instability in deterministic models?

### Specific research content:
- **Deterministic models**:
- **Lotka - Volterra model**: Explore the existence and stability of periodic treatment regimens by analyzing the system of differential equations in different treatment stages.
- **Adjusted replicator dynamics model**: Study the design and stability of periodic treatment regimens through the adjusted replicator dynamics equation.

- **Stochastic models**:
- **Moran process**: Explore the stability of periodic treatment regimens in stochastic models and analyze the impact of random fluctuations on periodic treatment.

### Key conclusions:
- **Existence of cycles**: The author proposes biologically meaningful assumptions to ensure the existence of periodic treatment regimens in the Lotka - Volterra model and the adjusted replicator dynamics model.
- **Stability of cycles**: Not only is existence important, but the stability of cycles is also a key factor in practical applications. Stable periodic treatment regimens can maintain the control of tumor composition over a long period.
- **Stability in stochastic models**: The stability of periodic treatment regimens in stochastic models is closely related to their stability in deterministic models. If the cycle in the deterministic model is unstable, then the cycle collapse phenomenon will also occur in the stochastic model.

### Formula presentation:
- **Lotka - Volterra model**:
\[
\begin{cases}
x_1' = r_{1,i} x_1 \left(1 - \frac{x_1 + a_{12,i} x_2}{k_{1,i}}\right) \\
x_2' = r_{2,i} x_2 \left(1 - \frac{a_{21,i} x_1 + x_2}{k_{2,i}}\right)
\end{cases}
\]
where \( r_{j,i} \) and \( k_{j,i} \) are the Malthusian growth rate and carrying capacity of the \( j \)-th type of cell in the \( i \)-th treatment stage respectively, and \( a_{jl,i} \) is the aggressiveness of the \( l \)-th type of cell against the \( j \)-th type of cell.

- **Adjusted replicator dynamics model**:
\[
x_i' = x_i \left(f_i(x) - \Phi(x)\right) \quad \text{where} \quad f_i(x) = (Bx)_i \quad \text{and} \quad \Phi(x) = \sum_{i = 1}^3 x_i f_i(x)
\]
where \( B \) is the payoff matrix, representing the interactions between different types of cells.

Through these models and analyses, the author hopes to provide a theoretical basis and practical guidance for cancer adaptive therapy.