skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Coupled Feedback Loops Involving PAGE4, EMT and Notch Signaling Can Give Rise to Non-Genetic Heterogeneity in Prostate Cancer Cells
Non-genetic heterogeneity is emerging as a crucial factor underlying therapy resistance in multiple cancers. However, the design principles of regulatory networks underlying non-genetic heterogeneity in cancer remain poorly understood. Here, we investigate the coupled dynamics of feedback loops involving (a) oscillations in androgen receptor (AR) signaling mediated through an intrinsically disordered protein PAGE4, (b) multistability in epithelial–mesenchymal transition (EMT), and (c) Notch–Delta–Jagged signaling mediated cell-cell communication, each of which can generate non-genetic heterogeneity through multistability and/or oscillations. Our results show how different coupling strengths between AR and EMT signaling can lead to monostability, bistability, or oscillations in the levels of AR, as well as propagation of oscillations to EMT dynamics. These results reveal the emergent dynamics of coupled oscillatory and multi-stable systems and unravel mechanisms by which non-genetic heterogeneity in AR levels can be generated, which can act as a barrier to most existing therapies for prostate cancer patients.  more » « less
Award ID(s):
2019745
PAR ID:
10233838
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Entropy
Volume:
23
Issue:
3
ISSN:
1099-4300
Page Range / eLocation ID:
288
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    The epithelial-mesenchymal transition (EMT) and the corresponding reverse process, mesenchymalepithelial transition (MET), are dynamic and reversible cellular programs orchestrated by many changes at both biochemical and morphological levels. A recent surge in identifying the molecular mechanisms underlying EMT/MET has led to the development of various mathematical models that have contributed to our improved understanding of dynamics at single-cell and population levels: (a) multi-stability—how many phenotypes can cells attain during an EMT/MET?, (b) reversibility/irreversibility—what time and/or concentration of an EMT inducer marks the “tipping point” when cells induced to undergo EMT cannot revert?, (c) symmetry in EMT/MET—do cells take the same path when reverting as theytook during the induction of EMT?, and (d) non-cell autonomous mechanisms—how does a cell undergoing EMT alter the tendency of its neighbors to undergo EMT? These dynamical traits may facilitate a heterogenous response within a cell population undergoing EMT/MET. Here, we present a few examples of designing different mathematical models that can contribute to decoding EMT/MET dynamics. Key words Mathematical modeling, Epithelial-mesenchymal plasticity, Nongenetic heterogeneity,Multi-stability, Epithelial-mesenchymal heterogeneity 
    more » « less
  2. The transcriptional coactivator YAP1 (yes-associated protein 1) is a critical nuclear effector of the Hippo pathway. The serine/threonine protein kinases STK3/4 and LATS1/2, core components of the Hippo pathway, phosphorylate and inhibit YAP1 nuclear localization. Previously, we reported that the interaction of nuclear YAP1 with androgen receptor (AR) might play a critical role in prostate cancer progression and therapeutic relapse (Kuser-Abali et al., Nat. Commun. 2015). Here, we investigated the regulation of YAP1 by androgens in isogenic, androgen-responsive LNCaP and androgen non-responsive C4-2 prostate cancer cell models. We demonstrated that androgen suppressed the inhibitory phospho-Ser127 site on YAP1 in LNCaP cells, but the effects of androgen on phospho-Ser127 was modest in C4-2 cells. In agreement with this observation, androgen increased the presence of nuclear YAP1 in LNCaP cells, whereas regardless of androgen exposure the YAP1 protein was primarily expressed in C4-2 cell nuclei. We also demonstrated that androgen exposure suppressed the levels of phospho-Ser127 induced by okadaic acid, which is a potent inhibitor of the Ser/Thr phosphatases PP1 and PP2A. Moreover, the pharmacological inhibition of androgen receptor (AR) signaling by enzalutamide reversed the inhibitory effects of androgen on phospho-Ser127, which coincided with the inhibition of YAP1 nuclear localization. Similarly, the genetic inhibition of AR signaling by small interfering RNA (siRNA) reduced phospho-Ser127 levels. Additionally, the silencing of the STK3/4 and LATS1/2 signaling by siRNA resulted in increases in YAP1 protein levels. Furthermore, our analysis of the TCGA (The Cancer Genome Atlas) prostate adenocarcinoma data set indicates that the levels of YAP1 and AR mRNA expression were positively correlated in prostate cancer clinical samples. These observations suggest that AR signaling promotes YAP1 nuclear localization by suppressing phospho-Ser127, possibly through the protein phosphatases PP1 and PP2A, and supporting a new mechanism of YAP1 regulation and YAP1-mediated cancer cell growth and survival. 
    more » « less
  3. The transitions between epithelial (E) and mesenchymal (M) cell phenotypes are essential in many biological processes like tissue development and cancer metastasis. Previous studies, both modeling and experimental, suggested that in addition to E and M states, the network responsible for these phenotypes exhibits intermediate phenotypes between E and M states. The number and importance of such states is subject to intense discussion in the epithelial-mesenchymal transition (EMT) community. Previous modeling efforts used traditional bifurcation analysis to explore the number of the steady states that correspond to E, M and intermediate states by varying one or two parameters at a time. Since the system has dozens of parameters that are largely unknown, it remains a challenging problem to fully describe the potential set of states and their relationship across all parameters. We use the computational tool DSGRN (Dynamic Signatures Generated by Regulatory Networks) to explore the intermediate states of an EMT model network by computing summaries of the dynamics across all of parameter space. We find that the only attractors in the system are equilibria, that E and M states dominate across parameter space, but that bistability and multistability are common. Even at extreme levels of some of the known inducers of the transition, there is a certain proportion of the parameter space at which an E or an M state co-exists with other stable steady states. Our results suggest that the multistability is broadly present in the EMT network across parameters and thus response of cells to signals may strongly depend on the particular cell line and genetic background. 
    more » « less
  4. null (Ed.)
    Epithelial-to-mesenchymal transition (EMT) plays an important role in many biological processes during development and cancer. The advent of single-cell transcriptome sequencing techniques allows the dissection of dynamical details underlying EMT with unprecedented resolution. Despite several single-cell data analysis on EMT, how cell communicates and regulates dynamics along the EMT trajectory remains elusive. Using single-cell transcriptomic datasets, here we infer the cell–cell communications and the multilayer gene–gene regulation networks to analyze and visualize the complex cellular crosstalk and the underlying gene regulatory dynamics along EMT. Combining with trajectory analysis, our approach reveals the existence of multiple intermediate cell states (ICSs) with hybrid epithelial and mesenchymal features. Analyses on the time-series datasets from cancer cell lines with different inducing factors show that the induced EMTs are context-specific: the EMT induced by transforming growth factor B1 (TGFB1) is synchronous, whereas the EMTs induced by epidermal growth factor and tumor necrosis factor are asynchronous, and the responses of TGF-β pathway in terms of gene expression regulations are heterogeneous under different treatments or among various cell states. Meanwhile, network topology analysis suggests that the ICSs during EMT serve as the signaling in cellular communication under different conditions. Interestingly, our analysis of a mouse skin squamous cell carcinoma dataset also suggests regardless of the significant discrepancy in concrete genes between in vitro and in vivo EMT systems, the ICSs play dominant role in the TGF-β signaling crosstalk. Overall, our approach reveals the multiscale mechanisms coupling cell–cell communications and gene–gene regulations responsible for complex cell-state transitions. 
    more » « less
  5. null (Ed.)
    The epithelial-mesenchymal transition (EMT) plays a critical role in cancer progression, being responsible in many cases for the onset of the metastatic cascade and being integral in the ability of cells to resist drug treatment. Most studies of EMT focus on its induction via chemical signals such as TGF-β or Notch ligands, but it has become increasingly clear that biomechanical features of the microenvironment such as extracellular matrix (ECM) stiffness can be equally important. Here, we introduce a coupled feedback loop connecting stiffness to the EMT transcription factor ZEB1, which acts via increasing the secretion of LOXL2 that leads to increased cross-linking of collagen fibers in the ECM. This increased cross-linking can effectively increase ECM stiffness and increase ZEB1 levels, thus setting a positive feedback loop between ZEB1 and ECM stiffness. To investigate the impact of this non-cell-autonomous effect, we introduce a computational approach capable of connecting LOXL2 concentration to increased stiffness and thereby to higher ZEB1 levels. Our results indicate that this positive feedback loop, once activated, can effectively lock the cells in a mesenchymal state. The spatial-temporal heterogeneity of the LOXL2 concentration and thus the mechanical stiffness also has direct implications for migrating cells that attempt to escape the primary tumor. 
    more » « less