Abstract
HER2-positive (HER2+) metastatic breast cancer (mBC) is highly aggressive and a major threat to human health. Despite the significant improvement in patients' prognosis given the drug development efforts during the past several decades, many clinical questions still remain to be addressed such as efficacy when combining different therapeutic modalities, best treatment sequences, interindividual variability as well as resistance and potential coping strategies. To better answer these questions, we developed a mechanistic quantitative systems pharmacology model of the pathophysiology of HER2+ mBC that was extensively calibrated and validated against multiscale data to quantitatively predict and characterize the signal transduction and preclinical tumor growth kinetics under different therapeutic interventions. Focusing on the second-line treatment for HER2+ mBC, e.g., antibody-drug conjugates (ADC), small molecule inhibitors/TKI and chemotherapy, the model accurately predicted the efficacy of various drug combinations and dosing regimens at the in vitro and in vivo levels. Sensitivity analyses and subsequent heterogeneous phenotype simulations revealed important insights into the design of new drug combinations to effectively overcome various resistance scenarios in HER2+ mBC treatments. In addition, the model predicted a better efficacy of the new TKI plus ADC combination which can potentially reduce drug dosage and toxicity, while it also shed light on the optimal treatment ordering of ADC versus TKI plus capecitabine regimens, and these findings were validated by new in vivo experiments. Our model is the first that mechanistically integrates multiple key drug modalities in HER2+ mBC research and it can serve as a high-throughput computational platform to guide future model-informed drug development and clinical translation.