Towards personalized computer simulation of breast cancer treatment: a multi-scale pharmacokinetic and pharmacodynamic model informed by multi-type patient data

Mathematical modeling and simulation have emerged as a potentially powerful, time and cost effective approach to personalized cancer treatment. The usefulness of mechanistic models to disentangle complex multi-scale cancer processes such as treatment response has been widely acknowledged. However, a major barrier for multi-scale models to predict the outcomes of therapeutic regimens in a particular patient lies in their initialization and parameterization which need to reflect individual cancer characteristics accurately. In this study we use multi-type routinely acquired measurements on a single breast tumor, including histopathology, magnetic resonance imaging, and molecular profiling to personalize parts of a complex multi-scale model of breast cancer treated with chemotherapeutic and anti-angiogenic agents. We model the dynamics of drugs in tissue (pharmacokinetics) and the corresponding effects on their targets (pharmacodynamics). We developed a open-source computer program that simulates cross-sections of tumors under 12-week therapy regimes and use it to individually reproduce and elucidate treatment outcomes of four patients. For two of the tumors that did not respond to therapy, we used model simulations to suggest alternative regimes, depending on their individual characteristics, with improved outcomes. We found that more frequent doses of chemothereapy reduce tumor burden in a low proliferative tumor while lower doses of anti-angiogenic agents improve drug penetration in a poorly perfused tumor. In addition to bridge multi-type clinical data to shed light on individual treatment outcomes, our approach identified a few tumor-related aspects that need to be clinically portraited better to allow for future model-driven personalized cancer therapy.