Abstract
Oncogenic FGFR4 signaling represents a potential therapeutic target in various cancer types, including triple-negative breast cancer and hepatocellular carcinoma. However, resistance to FGFR4 single-agent therapy remains a major challenge, emphasizing the need for effective combinatorial treatments. Our study sought to develop a comprehensive computational model of FGFR4 signaling and to provide network-level insights into resistance mechanisms driven by signaling dynamics. An integrated approach, combining computational network modeling with experimental validation, uncovered potent AKT reactivation following FGFR4 targeting in triple-negative breast cancer cells. Analyzing the effects of cotargeting specific network nodes by systematically simulating the model predicted synergy of cotargeting FGFR4 and AKT or specific ErbB kinases, which was subsequently confirmed through experimental validation; however, cotargeting FGFR4 and PI3K was not synergistic. Protein expression data from hundreds of cancer cell lines was incorporated to adapt the model to diverse cellular contexts. This revealed that although AKT rebound was common, it was not a general phenomenon. For example, ERK reactivation occurred in certain cell types, including an FGFR4-driven hepatocellular carcinoma cell line, in which there is a synergistic effect of cotargeting FGFR4 and MEK but not AKT. In summary, this study offers key insights into drug-induced network remodeling and the role of protein expression heterogeneity in targeted therapy responses. These findings underscore the utility of computational network modeling for designing cell type–selective combination therapies and enhancing precision cancer treatment.
Significance: Computational predictive modeling of signaling networks can decipher mechanisms of cancer cell resistance to targeted therapies and enable identification of more effective cancer type–specific combination treatment strategies.