Extracellular growth factors signal to transcription factors via a limited number of cytoplasmic kinase cascades. It remains unclear how such cascades encode ligand identities and concentrations. In this paper, we use live-cell imaging and statistical modeling to study FOXO3, a transcription factor regulating diverse aspects of cellular physiology that is under combinatorial control. We show that FOXO3 nuclear-to-cytosolic translocation has two temporally distinct phases varying in magnitude with growth factor identity and cell type. These phases comprise synchronous translocation soon after ligand addition followed by an extended back-and-forth shuttling; this shuttling is pulsatile and does not have a characteristic frequency, unlike a simple oscillator. Early and late dynamics are differentially regulated by Akt and ERK and have low mutual information, potentially allowing the two phases to encode different information. In cancer cells in which ERK and Akt are dysregulated by oncogenic mutation, the diversity of states is lower.

In most solid cancers, cells harboring oncogenic mutations represent only a sub-fraction of the entire population. Within this sub-fraction the expression level of mutated proteins can vary significantly due to cellular variability limiting the efficiency of targeted therapy. To address the causes of the heterogeneity, we performed a systematic analysis of one of the most frequently mutated pathways in cancer cells, the phosphatidylinositol 3 kinase (PI3K) signaling pathway. Among others PI3K signaling is activated by the hepatocyte growth factor (HGF) that regulates proliferation of hepatocytes during liver regeneration but also fosters tumor cell proliferation. HGF-mediated responses of PI3K signaling were monitored both at the single cell and cell population level in primary mouse hepatocytes and in the hepatoma cell line Hepa1_6. Interestingly, we observed that the HGF-mediated AKT responses at the level of individual cells is rather heterogeneous. However, the overall average behavior of the single cells strongly resembled the dynamics of AKT activation determined at the cell population level. To gain insights into the molecular cause for the observed heterogeneous behavior of individual cells, we employed dynamic mathematical modeling in a stochastic framework. Our analysis demonstrated that intrinsic noise was not sufficient to explain the observed kinetic behavior, but rather the importance of extrinsic noise has to be considered. Thus, distinct from gene expression in the examined signaling pathway fluctuations of the reaction rates has only a minor impact whereas variability in the concentration of the various signaling components even in a clonal cell population is a key determinant for the kinetic behavior.