Most models of cancer cell population expansion assume exponential growth kinetics at low cell densities, with deviations to account for observed slowing of growth rate only at higher densities due to limited resources such as space and nutrients. However, recent preclinical and clinical observations of tumor initiation or recurrence indicate the presence of tumor growth kinetics in which growth rates scale positively with cell numbers. These observations are analogous to the cooperative behavior of species in an ecosystem described by the ecological principle of the Allee effect. In preclinical and clinical models, however, tumor growth data are limited by the lower limit of detection (i.e., a measurable lesion) and confounding variables, such as tumor microenvironment, and immune responses may cause and mask deviations from exponential growth models. In this work, we present alternative growth models to investigate the presence of an Allee effect in cancer cells seeded at low cell densities in a controlled in vitro setting. We propose a stochastic modeling framework to disentangle expected deviations due to small population size stochastic effects from cooperative growth and use the moment approach for stochastic parameter estimation to calibrate the observed growth trajectories. We validate the framework on simulated data and apply this approach to longitudinal cell proliferation data of BT-474 luminal B breast cancer cells. We find that cell population growth kinetics are best described by a model structure that considers the Allee effect, in that the birth rate of tumor cells increases with cell number in the regime of small population size. This indicates a potentially critical role of cooperative behavior among tumor cells at low cell densities with relevance to early stage growth patterns of emerging and relapsed tumors.

While acquired chemoresistance is recognized as a key challenge to treating many types of cancer, the dynamics with which drug sensitivity changes after exposure are poorly characterized. Most chemotherapeutic regimens call for repeated dosing at regular intervals, and if drug sensitivity changes on a similar time scale then the treatment interval could be optimized to improve treatment performance. Theoretical work suggests that such optimal schedules exist, but experimental confirmation has been obstructed by the difficulty of deconvolving the simultaneous processes of death, adaptation, and regrowth taking place in cancer cell populations. Here we present a method of optimizing drug schedules in vitro through iterative application of experimentally calibrated models, and demonstrate its ability to characterize dynamic changes in sensitivity to the chemotherapeutic doxorubicin in three breast cancer cell lines subjected to treatment schedules varying in concentration, interval between pulse treatments, and number of sequential pulse treatments. Cell populations are monitored longitudinally through automated imaging for 600–800 hours, and this data is used to calibrate a family of cancer growth models, each consisting of a system of ordinary differential equations, derived from the bi-exponential model which characterizes resistant and sensitive subpopulations. We identify a model incorporating both a period of growth arrest in surviving cells and a delay in the death of chemosensitive cells which outperforms the original bi-exponential growth model in Akaike Information Criterion based model selection, and use the calibrated model to quantify the performance of each drug schedule. We find that the inter-treatment interval is a key variable in determining the performance of sequential dosing schedules and identify an optimal retreatment time for each cell line which extends regrowth time by 40%-239%, demonstrating that the time scale of changes in chemosensitivity following doxorubicin exposure allows optimization of drug scheduling by varying this inter-treatment interval.

Models of cancer cell population expansion assume exponential growth kinetics at low cell densities, with deviations from exponential growth only at higher densities due to limited resources such as space and nutrients. However, recent pre-clinical and clinical observations of tumor initiation or recurrence indicate the presence of tumor growth kinetics in which growth rates scale with cell numbers. These observations are analogous to the cooperative behavior of species in an ecosystem described by the ecological principle of the Allee effect. In preclinical and clinical models however, tumor growth data is limited by the lower limit of detection (i.e. a measurable lesion) and confounding variables, such as tumor microenvironment and immune responses may cause and mask deviations from exponential growth models. In this work, we present alternative growth models to investigate the presence of an Allee effect in cancer cells seeded at low cell densities in a controlled in vitro setting. We propose a stochastic modeling framework to consider the small number of cells in this low-density regime and use the moment approach for stochastic parameter estimation to calibrate the stochastic growth trajectories. We validate the framework on simulated data and apply this approach to longitudinal cell proliferation data of BT-474 luminal B breast cancer cells. We find that cell population growth kinetics are best described by a model structure that considers the Allee effect, in that the birth rate of tumor cells depends on cell number. This indicates a potentially critical role of cooperative behavior among tumor cells at low cell densities with relevance to early stage growth patterns of emerging tumors and relapse. Author SummaryThe growth kinetics of cancer cells at very low cell densities are of utmost clinical importance as the ability of a small number of newly transformed or surviving cells to grow exponentially and thus, to "take off" underlies tumor formation and relapse after treatment. Mathematical models of stochastic tumor cell growth typically assume a stochastic birth-death process of cells impacted by limited nutrients and space when cells reach high density, resulting in the widely accepted logistic growth model. Here we present an in-depth investigation of alternate growth models adopted from ecology to describe potential deviations from a simple cell autonomous birth-death model at low cell densities. We show that our stochastic modeling framework is robust and can be used to identify the underlying structure of stochastic growth trajectories from both simulated and experimental data taken from a controlled in vitro setting in which we can capture data from the relevant low cell density regime. This work suggests that the assumption of cell autonomous proliferation via a constant exponential growth rate at low cell densities may not be appropriate for all cancer cell growth dynamics. Consideration of cooperative behavior amongst tumor cells in this regime is critical for elucidating strategies for controll...

The development of resistance to chemotherapy is a major cause of treatment failure in breast cancer. While mathematical models describing the dynamics of resistant cancer cell subpopulations have been proposed, experimental validation has been difficult due to the complex nature of resistance that limits the ability of a single phenotypic marker to sufficiently identify the drug resistant subpopulations. We address this problem with a coupled experimental/modeling approach to reveal the composition of drug resistant subpopulations changing in time following drug exposure. We calibrate time-resolved drug sensitivity assays to three mathematical models to interrogate the models’ ability to capture drug response dynamics. The Akaike information criterion was employed to evaluate the three models, and it identified a multi-state model incorporating the role of population heterogeneity and cellular plasticity as the optimal model. To validate the model’s ability to identify subpopulation composition, we mixed different proportions of wild-type MCF-7 and MCF-7/ADR resistant cells and evaluated the corresponding model output. Our blinded two-state model was able to estimate the proportions of cell types with an R-squared value of 0.857. To the best of our knowledge, this is the first work to combine experimental time-resolved drug sensitivity data with a mathematical model of resistance development.

While acquired chemoresistance is recognized as a key challenge to treating many types of cancer, the dynamics with which drug sensitivity changes after exposure are poorly characterized.  Most chemotherapeutic regimens call for repeated dosing at regular intervals, and if drug sensitivity changes on a similar time scale then the treatment interval could be optimized to improve treatment performance.  Theoretical work suggests that such optimal schedules exist, but experimental confirmation has been obstructed by the difficulty of deconvolving the simultaneous processes of death, adaptation, and regrowth taking place in cancer cell populations.  Here we present work characterizing dynamic changes in sensitivity to the chemotherapeutic doxorubicin in three breast cancer cell lines subjected to treatment schedules varying in concentration, interval between pulse treatments, and number of sequential pulse treatments. Cell populations are monitored longitudinally through automated imaging for 600-800 hours, and this data is used to calibrate a family of cancer growth models derived from the bi-exponential model which characterizes resistant and sensitive subpopulations.  We identify a model incorporating both a period of growth arrest in surviving cells and a delay in the death of chemosensitive cells which outperforms the original bi-exponential growth model in Akaike Information Criterion based model selection, and use the calibrated model to quantify the performance of each drug schedule.  We find that the inter-treatment interval is a key variable in determining the performance of sequential dosing schedules and identify an optimal retreatment time for each cell line which extends regrowth time by 40%-106%, demonstrating that the time scale of changes in chemosensitivity following doxorubicin exposure allows optimization of drug scheduling by varying this inter-treatment interval.