Characterizing and quantifying the effects of breast cancer therapy using mathematical modeling

We designed a mathematical model to describe and quantify the mechanisms and dynamics of tumor growth, cell-kill and resistance as they affect durations of benefit after cancer treatment. Our aim was to explore how treatment efficacy may be related to primary tumor characteristics, with the potential to guide future trial design and appropriate selection of therapy. Assuming a log-normal distribution of both resistant disease and tumor doubling times generates disease-free survival (DFS) or invasive DFS curves with specific shapes. Using a multivariate mathematical model, both treatment and tumor characteristics are related to quantified resistant disease and tumor regrowth rates by allowing different mean values for the influence of different treatments or clinical subtypes on these two log-normal distributions. Application of the model to the CALGB 9741 adjuvant breast cancer trial showed that dose-dense therapy was estimated to achieve an extra 3/4 log of cell-kill compared to standard therapy, but only in patients with more rapidly growing ER-negative tumors. Application of the model to the AZURE trial of adjuvant bisphosphonate treatment suggested that the 5-year duration of zoledronic acid was adequate for ER-negative tumors, but may not be so for ER-positive cases, with increased recurrences after ceasing the intervention. Mathematical models can identify different effects of treatment by subgroup and may aid in treatment design, trial analysis, and appropriate selection of therapy. They may provide a more appropriate and insightful tool than the conventional Cox model for the statistical analysis of response durations. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s10549-016-3684-4) contains supplementary material, which is available to authorized users.