Introduction: A defining hallmark of cancer is aberrant cell proliferation. Efforts to understand the generative properties of cancer cells span all biological scales: from genetic deviations and alterations of metabolic pathways, to physical stresses due to overcrowding, as well as the effects of therapeutics and the immune system. While these factors have long been studied in the laboratory, mathematical and computational techniques are being increasingly applied to help understand and forecast tumor growth and treatment response. Advantages of mathematical modeling of proliferation include the ability to simulate and predict the spatiotemporal development of tumors across multiple experimental scales. Central to proliferation modeling is the incorporation of available biological data and validation with experimental data. Areas covered: We present an overview of past and current mathematical strategies directed at understanding tumor cell proliferation. We identify areas for mathematical development as motivated by available experimental and clinical evidence, with a particular emphasis on emerging, non-invasive imaging technologies. Expert commentary: The data required to legitimize mathematical models are often difficult or (currently) impossible to obtain. We suggest areas for further investigation to establish mathematical models that more effectively utilize available data to make informed predictions on tumor cell proliferation.

Optimal control theory is branch of mathematics that aims to optimize a solution to a dynamical system. While the concept of using optimal control theory to improve treatment regimens in oncology is not novel, many of the early applications of this mathematical technique were not designed to work with routinely available data or produce results that can eventually be translated to the clinical setting. The purpose of this review is to discuss clinically relevant considerations for formulating and solving optimal control problems for treating cancer patients. Our review focuses on two of the most widely used cancer treatments, radiation therapy and systemic therapy, as they naturally lend themselves to optimal control theory as a means to personalize therapeutic plans in a rigorous fashion. To provide context for optimal control theory to address either of these two modalities, we first discuss the major limitations and difficulties oncologists face when considering alternate regimens for their patients. We then provide a brief introduction to optimal control theory before formulating the optimal control problem in the context of radiation and systemic therapy. We also summarize examples from the literature that illustrate these concepts. Finally, we present both challenges and opportunities for dramatically improving patient outcomes via the integration of clinically relevant, patient-specific, mathematical models and optimal control theory.

Triple-negative breast cancer (TNBC) is persistently refractory to therapy, and methods to improve targeting and evaluation of responses to therapy in this disease are needed. Here, we integrate quantitative MRI data with biologically based mathematical modeling to accurately predict the response of TNBC to neoadjuvant systemic therapy (NAST) on an individual basis. Specifically, 56 patients with TNBC enrolled in the ARTEMIS trial (NCT02276443) underwent standard-of-care doxorubicin/cyclophosphamide (A/C) and then paclitaxel for NAST, where dynamic contrast-enhanced MRI and diffusion-weighted MRI were acquired before treatment and after two and four cycles of A/C. A biologically based model was established to characterize tumor cell movement, proliferation, and treatment-induced cell death. Two evaluation frameworks were investigated using: (i) images acquired before and after two cycles of A/C for calibration and predicting tumor status after A/C, and (ii) images acquired before, after two cycles, and after four cycles of A/C for calibration and predicting response following NAST. For Framework 1, the concordance correlation coefficients between the predicted and measured patient-specific, post-A/C changes in tumor cellularity and volume were 0.95 and 0.94, respectively. For Framework 2, the biologically based model achieved an area under the receiver operator characteristic curve of 0.89 (sensitivity/specificity = 0.72/0.95) for differentiating pathological complete response (pCR) from non-pCR, which is statistically superior (P < 0.05) to the value of 0.78 (sensitivity/specificity = 0.72/0.79) achieved by tumor volume measured after four cycles of A/C. Overall, this model successfully captured patient-specific, spatiotemporal dynamics of TNBC response to NAST, providing highly accurate predictions of NAST response. Significance: Integrating MRI data with biologically based mathematical modeling successfully predicts breast cancer response to chemotherapy, suggesting digital twins could facilitate a paradigm shift from simply assessing response to predicting and optimizing therapeutic efficacy.

Tumor-associated vasculature is responsible for the delivery of nutrients, removal of waste, and allowing growth beyond 2–3 mm³. Additionally, the vascular network, which is changing in both space and time, fundamentally influences tumor response to both systemic and radiation therapy. Thus, a robust understanding of vascular dynamics is necessary to accurately predict tumor growth, as well as establish optimal treatment protocols to achieve optimal tumor control. Such a goal requires the intimate integration of both theory and experiment. Quantitative and time-resolved imaging methods have emerged as technologies able to visualize and characterize tumor vascular properties before and during therapy at the tissue and cell scale. Parallel to, but separate from those developments, mathematical modeling techniques have been developed to enable in silico investigations into theoretical tumor and vascular dynamics. In particular, recent efforts have sought to integrate both theory and experiment to enable data-driven mathematical modeling. Such mathematical models are calibrated by data obtained from individual tumor-vascular systems to predict future vascular growth, delivery of systemic agents, and response to radiotherapy. In this review, we discuss experimental techniques for visualizing and quantifying vascular dynamics including magnetic resonance imaging, microfluidic devices, and confocal microscopy. We then focus on the integration of these experimental measures with biologically based mathematical models to generate testable predictions.

The ability to accurately predict response and then rigorously optimize a therapeutic regimen on a patient-specific basis, would transform oncology. Toward this end, we have developed an experimental-mathematical framework that integrates quantitative magnetic resonance imaging (MRI) data into a biophysical model to predict patient-specific treatment response of locally advanced breast cancer to neoadjuvant therapy. Diffusion-weighted and dynamic contrast-enhanced MRI data is collected prior to therapy, after 1 cycle of therapy, and at the completion of the first therapeutic regimen. The model is initialized and calibrated with the first 2 patient-specific MRI data sets to predict response at the third, which is then compared to patient outcomes (N = 18). The model's predictions for total cellularity, total volume, and the longest axis at the completion of the regimen are significant within expected measurement precision (P< 0.05) and strongly correlated with measured response (P < 0.01). Further, we use the model to investigate, in silico, a range of (practical) alternative treatment plans to achieve the greatest possible tumor control for each individual in a subgroup of patients (N = 13). The model identifies alternative dosing strategies predicted to achieve greater tumor control compared to the standard of care for 12 of 13 patients (P < 0.01). In summary, a predictive, mechanism-based mathematical model has demonstrated the ability to identify alternative treatment regimens that are forecasted to outperform the therapeutic regimens the patients clinically. This has important implications for clinical trial design with the opportunity to alter oncology care in the future.

Clinical methods for assessing tumor response to therapy are largely rudimentary, monitoring only temporal changes in tumor size. Our goal is to predict the response of breast tumors to therapy using a mathematical model that utilizes magnetic resonance imaging (MRI) data obtained non-invasively from individual patients. We extended a previously established, mechanically coupled, reaction-diffusion model for predicting tumor response initialized with patient-specific diffusion weighted MRI (DW-MRI) data by including the effects of chemotherapy drug delivery, which is estimated using dynamic contrast-enhanced (DCE-) MRI data. The extended, drug incorporated, model is initialized using patient-specific DW-MRI and DCE-MRI data. Data sets from five breast cancer patients were used—obtained before, after one cycle, and at mid-point of neoadjuvant chemotherapy. The DCE-MRI data was used to estimate spatiotemporal variations in tumor perfusion with the extended Kety–Tofts model. The physiological parameters derived from DCE-MRI were used to model changes in delivery of therapy drugs within the tumor for incorporation in the extended model. We simulated the original model and the extended model in both 2D and 3D and compare the results for this five-patient cohort. Preliminary results show reductions in the error of model predicted tumor cellularity and size compared to the experimentally-measured results for the third MRI scan when therapy was incorporated. Comparing the two models for agreement between the predicted total cellularity and the calculated total cellularity (from the DW-MRI data) reveals an increased concordance correlation coefficient from 0.81 to 0.98 for the 2D analysis and 0.85 to 0.99 for the 3D analysis (p < 0.01 for each) when the extended model was used in place of the original model. This study demonstrates the plausibility of using DCE-MRI data as a means to estimate drug delivery on a patient-specific basis in predictive models and represents a step toward the goal of achieving individualized prediction of tumor response to therapy.