On minimising tumoural growth under treatment resistance

Abstract

Drug resistance is a major challenge for curative cancer treatment, representing the main reason of death in patients. Evolutionary biology suggests pauses between treatment rounds as a way to delay or even avoid resistance emergence. Indeed, this approach has already shown promising preclinical and early clinical results, and stimulated the development of mathematical models for finding optimal treatment protocols. Due to their complexity, however, these models do not lend themself to a rigorous mathematical analysis, hence so far clinical recommendations generally relied on numerical simulations and ad-hoc heuristics. Here, we derive two mathematical models describing tumour growth under genetic and epigenetic treatment resistance, respectively, which are simple enough for a complete analytical investigation. First, we find key differences in response to treatment protocols between the two modes of resistance. Second, we identify the optimal treatment protocol which leads to the largest possible tumour shrinkage rate. Third, we fit the “epigenetic model” to previously published xenograft experiment data, finding excellent agreement, underscoring the biological validity of our approach. Finally, we use the fitted model to calculate the optimal treatment protocol for this specific experiment, which we demonstrate to cause curative treatment, making it superior to previous approaches which generally aimed at stabilising tumour burden. Overall, our approach underscores the usefulness of simple mathematical models and their analytical examination, and we anticipate our findings to guide future preclinical and, ultimately, clinical research in optimising treatment regimes.