Theoretical modeling of collaterally sensitive drug cycles: shaping heterogeneity to allow adaptive therapy

In previous work, we focused on the optimal therapeutic strategy with a pair of drugs which are collaterally sensitive to each other, that is, a situation in which evolution of resistance to one drug induces sensitivity to the other, and vice versa. Yoona (Bull Math Biol 8:1-34,Yoon et al. 2018) Here, we have extended this exploration to the optimal strategy with a collaterally sensitive drug sequence of an arbitrary length, N. To explore this, we have developed a dynamical model of sequential drug therapies with N drugs. In this model, tumor cells are classified as one of N subpopulations represented as {R-i vertical bar i = 1, 2, ... , N}. Each subpopulation, R-i, is resistant to 'Drug i' and each subpopulation, Ri-1 (or R-N, if i = 1), is sensitive to it, so that R-i increases under 'Drug i' as it is resistant to it, and after drug-switching, decreases under 'Drug i + 1' as it is sensitive to that drug(s). Similar to our previous work examining optimal therapy with two drugs, we found that there is an initial period of time in which the tumor is 'shaped' into a specific makeup of each subpopulation, at which time all the drugs are equally effective (R*). After this shaping period, all the drugs are quickly switched with duration relative to their efficacy in order to maintain each subpopulation, consistent with the ideas underlying adaptive therapy. West(Canver Res 80(7):578-589Gatenby et al. 2009) and Gatenby (Cancer Res 67(11):4894-4903West et al. 2020). Additionally, we have developed methodologies to administer the optimal regimen under clinical or experimental situations in which no drug parameters and limited information of trackable populations data (all the subpopulations or only total population) are known. The therapy simulation based on these methodologies showed consistency with the theoretical effect of optimal therapy .