Understanding therapeutic tolerance through a mathematical model of drug-induced resistance

Jana L. Gevertz,James M. Greene,Samantha Prosperi,Natacha Comandante-Lou,Eduardo D. Sontag
DOI: https://doi.org/10.1101/2024.09.04.611211
2024-09-07
Abstract:There is growing recognition that phenotypic plasticity enables cancer cells to adapt to various environmental conditions. An example of this adaptability is the persistence of an initially sensitive population of cancer cells in the presence of therapeutic agents. Understanding the implications of this druginduced resistance is essential for predicting transient and long-term tumor tumor dynamics subject to treatment. This paper introduces a mathematical model of this phenomenon of drug-induced resistance which provides excellent fits to time-resolved experimental data. From observational data of total numbers of cells, the model unravels the relative proportions of sensitive and resistance subpopulations, and quantifies their dynamics as a function of drug dose. The predictions are then validated using data on drug doses which were not used when fitting parameters. The model is then used, in conjunction with optimal control techniques, in order to discover dosing strategies that might lead to better outcomes as quantified by lower total cell volume.
Systems Biology
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to understand the mechanism of drug - induced resistance in cancer treatment. Specifically, the author explains why cancer cells that were originally sensitive to drugs gradually develop resistance during chemotherapy or targeted therapy by establishing a mathematical model. This problem is of great significance for predicting the dynamic changes of tumors during treatment and designing more effective treatment regimens. ### Background and Problem Description of the Paper 1. **Drug - Induced Resistance** - Drug - induced resistance refers to the fact that during treatment, cancer cells that were originally sensitive to drugs gradually develop resistance to the drugs. This phenomenon weakens the treatment effect and even leads to treatment failure. - Understanding this phenomenon is crucial for predicting the dynamic changes of tumors at different doses and designing more effective treatment strategies. 2. **Research Motivation** - Traditional treatment methods are often based on fixed - dose administration regimens, but these regimens may not be able to effectively deal with drug - induced resistance. - By establishing a mathematical model, the mechanism of drug - induced resistance can be better understood and better dosing strategies can be explored. ### Construction of the Mathematical Model 1. **Model Overview** - The author constructed a mathematical model that includes sensitive cells (\( S \)), irreversibly resistant cells (\( R_n \)) and reversibly resistant cells (\( R_r \)). - The model describes the kinetic changes of these cell populations through a set of ordinary differential equations (ODEs). 2. **Model Equations** - The kinetic equation of sensitive cells \( S \): \[ \dot{S} = r_S \left(1 - \frac{V}{K}\right) S - (\epsilon_n + \epsilon_r) S - (\alpha_n + \alpha_r) u(t) S - d_s u(t) S + \eta R_r \] - The kinetic equation of irreversibly resistant cells \( R_n \): \[ \dot{R}_n = r_n \left(1 - \frac{V}{K}\right) R_n + \epsilon_n S + \alpha_n u(t) S - d_n u(t) R_n \] - The kinetic equation of reversibly resistant cells \( R_r \): \[ \dot{R}_r = r_r \left(1 - \frac{V}{K}\right) R_r + \epsilon_r S + \alpha_r u(t) S - d_r u(t) R_r - \eta R_r \] 3. **Parameter Explanation** - \( r_S, r_n, r_r \): The growth rates of sensitive cells, irreversibly resistant cells and reversibly resistant cells respectively. - \( \epsilon_n, \epsilon_r \): The conversion rates of spontaneously generated irreversibly resistant cells and reversibly resistant cells respectively. - \( \alpha_n, \alpha_r \): The conversion rates of drug - induced irreversibly resistant cells and reversibly resistant cells respectively. - \( d_s, d_n, d_r \): The killing rates of drugs on sensitive cells, irreversibly resistant cells and reversibly resistant cells respectively. - \( \eta \): The conversion rate at which reversibly resistant cells become sensitive cells again. - \( u(t) \): Represents the drug dose applied at time \( t \). - \( V \): The total number of tumor cells, that is, \( V = S + R_n + R_r \). - \( K \): The maximum carrying capacity of the environment. ### Experimental Data and Model Validation 1. **Experimental Data** - The author used in vitro experimental data, which showed that under different doses of vemurafenib treatment, COLO85