AbstractAdaptation and survival of cancer cells to various stress and growth factor conditions is crucial for successful metastasis. A double-negative feedback loop between two serine/threonine kinases AMPK and Akt can regulate the adaptation of breast cancer cells to matrix-deprivation stress. This feedback loop can generate majorly two phenotypes or cell states: matrix detachment-triggered pAMPKhigh/ pAktlow state, and matrix (re)attachment-triggered pAkthigh/ pAMPKlow state. However, whether these two cell states can exhibit phenotypic plasticity and heterogeneity in a given cell population, i.e., whether they can co-exist and undergo spontaneous switching to generate the other subpopulation, remains unclear. Here, we develop a mechanism-based mathematical model that captures the set of experimentally reported interactions among AMPK and Akt. Our simulations suggest that the AMPK-Akt feedback loop can give rise to two co-existing phenotypes (pAkthigh/ pAMPKlow and pAMPKhigh/pAktlow) in specific parameter regimes. Next, to test the model predictions, we segregated these two subpopulations in MDA-MB-231 cells and observed that each of them was capable of switching to another in adherent conditions. Finally, the predicted trends are supported by clinical data analysis of TCGA breast cancer and pan-cancer cohorts that revealed negatively correlated pAMPK and pAkt protein levels. Overall, our integrated computational-experimental approach unravels that AMPK-Akt feedback loop can generate multistability and drive phenotypic switching and heterogeneity in a cancer cell population.
A double-negative feedback loop between EpCAM and ERK contributes to the regulation of epithelial–mesenchymal transition in cancer
