Mathematical Modeling of Therapies for MCF7 Breast Cancer Cells

Estrogen receptor (ER)-positive breast cancer is responsive to a number of targeted therapies used clinically. Unfortunately, the continuous application of any targeted therapy often results in resistance to the therapy. Our ultimate goal is to use mathematical modelling to optimize alternating therapies that not only decrease proliferation but also stave off resistance. Toward this end, we measured levels of key proteins and proliferation over a 7-day time course in ER-positive MCF7 breast cancer cells. Treatments included endocrine therapy, either estrogen deprivation, which mimics the effects of an aromatase inhibitor, or fulvestrant, an ER degrader. These data were used to calibrate a mathematical model based on key interactions between ER signaling and the cell cycle. We show that the calibrated model is capable of predicting the combination treatment of fulvestrant and estrogen deprivation. Further, we show that we can add a new drug, palbociclib, to the model by measuring only two key proteins, c-Myc and hyperphosphorylated RB1, and adjusting only parameters associated with the drug. The model is then able to predict the combination treatment of estrogen deprivation and palbociclib. Then we added the dynamics of estrogen concentration in the medium into the model and extended the short-term model to a long-term model. The long-term model can simulate various mono- or combination treatments at different doses over 28 days. In addition to palbociclib, we add another Cdk4/6 inhibitor to the model, abemaciclib, which can induce apoptosis at high concentrations. Then the model can match the effects of abemaciclib treatment at two different doses and also capture the apoptosis effects induced by abemaciclib. After calibrating the model to these different treatment conditions, we used the model to explore the synergism among these different treatments. The mathematical model predicts a significant synergism between palbociclib or abemaciclib in combination with fulv