models with androgen dynamics [15]. d) Allee effect: cell populations that grow more slowly
(per capita) at low populations [16] effectively introduce an element of cooperation, e) Phenotypic
plasticity: rapid changes in cell phenotypes can generate resistance far more quickly than mutation
or population dynamics [17].
A key result from our work on the basic model is a tradeoff curve between time for resistant cells
to emerge and the mean cancer burden [18]. This tradeoff holds for both adaptive (dose skipping)
and intermittent therapies, and is robust across all model extensions except for the Allee effect
and cell plasticity. With an Allee effect, results are quite different. Aggressive therapy can drive
cells below the threshold, and prevent both resistant cells and total cells from reaching their upper
thresholds. Adaptive therapy, by backing off early to avoid favoring resistance, can behave quite
poorly, leading to escape times nearly as short as those with no therapy and with a high total cell
burden.
With the exception of the success of high dose therapy with a strong Allee effect, no universal
therapy can achieve all three objectives of lowering average dose, delaying time to emergence of
resistant cells, and reducing total tumor burden. All strategies show a tradeoff between delaying
emergence of resistant cells and a high cancer burden. Choosing the appropriate treatment requires
assessing the individual patients and specific cancers, and include factors often not included in
models, such as therapy toxicity [19]. Phenotypic plasticity, where resistance is induced by therapy
rather than arising from mutations or pre-existing variants [17], makes resistance much more difficult
to suppress [20]. Reversible behaviors can create complex responses to therapeutic timing [21].
Effective adaptive therapies require fitting data on individual patients, and data may lack the
resolution to distinguish among alternative models. In the case of PSA in prostate cancer, a simple
model [21], a more complex model with basic androgen dynamics [22], and a detailed model
of androgen dynamics [15] all fit data on a set of patients reasonably well, although with some
exceptions [23]. If models can be fit to the dynamics, adaptive therapies may be more robust
to patient variability than prescribed timing of intermittent therapy. Although data may lack the
resolution to identify specific mechanisms of interaction, such as the strength of competition
between different cancer cell phenotypes, simple models may have the greatest potential to capture
dynamics and guide therapy. The ideal combination will be patient-specific models combined with
in vivo data, perhaps with mouse PDX models [24], or in vitro data on patient derived cells that can
reveal mechanisms of interaction in different treatment environments.
1.2 How competitive are treatment resistant phenotypes?
Rob Noble: Adaptive therapy aims to exploit competition between treatment-sensitive and resistant
cells. Key questions remain largely unanswered. First, what is the nature of this competition?
Mathematical modellers typically assume that the fitness of resistant cells is a simple function
of their relative abundance and/or the total tumor size (reviewed in [25]). But frequency- or
density-dependent mathematical functions only approximate average population dynamics. Actual
clonal growth rates depend on the spatial arrangement of cells, their interaction ranges, and local
levels of shared resources, all of which vary both within and between tumors [26,27,28]. For
example, if a tumor grows mainly at its boundary then spatial constraints alone could suffice to
contain rare resistant clones, but only if they are located away from the boundary [29,30]. A
corollary is that the effectiveness of adaptive therapy may vary between cancer types due to different
tumor architectures [26]. Although spatially-structured computational and experimental models
can account for some important factors – such as competition for space and oxygen – the ability
4