Systems modelling of the EGFR-PYK2-c-Met interaction network predicts and prioritizes synergistic drug combinations for triple-negative breast cancer

Prediction of drug combinations that effectively target cancer cells is a critical challenge for cancer therapy, in particular for triple-negative breast cancer (TNBC), a highly aggressive breast cancer subtype with no effective targeted treatment. As signalling pathway networks critically control cancer cell behaviour, analysis of signalling network activity and crosstalk can help predict potent drug combinations and rational stratification of patients, thus bringing therapeutic and prognostic values. We have previously showed that the non-receptor tyrosine kinase PYK2 is a downstream effector of EGFR and c-Met and demonstrated their crosstalk signalling in basal-like TNBC. Here we applied a systems modelling approach and developed a mechanistic model of the integrated EGFR-PYK2-c-Met signalling network to identify and prioritize potent drug combinations for TNBC. Model predictions validated by experimental data revealed that among six potential combinations of drug pairs targeting the central nodes of the network, including EGFR, c-Met, PYK2 and STAT3, co-targeting of EGFR and PYK2 and to a lesser extent of EGFR and c-Met yielded strongest synergistic effect. Importantly, the synergy in co-targeting EGFR and PYK2 was linked to switch-like cell proliferation-associated responses. Moreover, simulations of patient-specific models using public gene expression data of TNBC patients led to predictive stratification of patients into subgroups displaying distinct susceptibility to specific drug combinations. These results suggest that mechanistic systems modelling is a powerful approach for the rational design, prediction and prioritization of potent combination therapies for individual patients, thus providing a concrete step towards personalized treatment for TNBC and other tumour types.Author summary: We applied a systems modelling approach combining mechanistic modelling and biological experimentation to identify effective drug combinations for triple-negative breast cancer (TNBC), an aggressive subtype of breast cancer with no approved targeted treatment. The model predicted and prioritized the synergistic combinations as confirmed by experimental data, demonstrating the power of this approach. Moreover, analysis of clinical data of TNBC patients and patient-specific modelling simulation enabled us to stratify the patients into subgroups with distinct susceptibility to specific drug combinations, and thus defined a subset of patient that could benefit from the combined treatments.