Quantitative analysis of time-resolved data in primary erythroid progenitor cells reveals that a dual negative transcriptional feedback mechanism underlies the ability of STAT5 to respond to the broad spectrum of physiologically relevant Epo concentrations.
A mathematical dual feedback model of the Epo-induced JAK2/STAT5 signaling pathway was calibrated with extensive time-resolved quantitative data sets from immunoblotting, mass spectrometry and qRT–PCR experiments in primary erythroid progenitor cells.We show that the amount of nuclear phosphorylated STAT5 integrated for 60 min post Epo stimulation directly correlates with the fraction of surviving cells 24 h later.CIS and SOCS3 were identified as the most relevant transcriptional feedback regulators of JAK2/STAT5 signaling in primary erythroid progenitor cells. Applying the model, we revealed that CIS-mediated inhibitory effects are most important at low ligand concentrations, whereas SOCS3 inhibition is more effective at high ligand doses.The distinct modes of inhibition of CIS and SOCS3 at various Epo concentrations provide a strategy for achieving control of JAK2/STAT5 signaling over the entire range of physiological Epo concentrations.
Cells interpret information encoded by extracellular stimuli through the activation of intracellular signaling networks and translate this information into cellular decisions. A prime example for a system that is exposed to extremely variable ligand concentrations is the erythroid lineage. The key regulator Erythropoietin (Epo) facilitates continuous renewal of erythrocytes at low basal levels but also secures compensation in case of, e.g., blood loss through an up to 1000-fold increase in hormone concentration. The Epo receptor (EpoR) is expressed on erythroid progenitor cells at the colony forming unit erythroid (CFU-E) stage. Stimulation of these cells with Epo leads to rapid but transient activation of receptor and JAK2 phosphorylation followed by phosphorylation of the latent transcription factor STAT5. Although STAT5 is known to be an essential regulator of survival and differentiation of erythroid progenitor cells, a quantitative link between the dynamic properties of STAT5 signaling and survival decisions remained unknown. STAT5-mediated responses in CFU-E cells are modulated by multiple attenuation mechanisms that operate on different time scales. Fast-acting mechanisms such as depletion of Epo by rapid receptor turnover and recruitment of the phosphatase SHP-1 control the initial signal amplitude at the receptor level. Transcriptional feedback regulators such as suppressor of cytokine signaling (SOCS) family members CIS and SOCS3 operate at a slower time scale. Despite the ample knowledge of the individual components involved, only little is known about the specific contributions of these regulators in controlling dynamic properties of STAT5 in response to a broad range of input signals. Therefore, dynamic pathway modeling is required to understand the complex regulatory network of feedback regulators.
To address these questions, we established a dual negative feedback model of JAK2/STAT5 signaling in primary erythroid progenitor cells isolated from mouse fetal livers. We provide a large data set of JAK2/STAT5 signaling dynamics employing quantitative immunoblotting, mass spectrometry and quantitative RT–PCR measured under different perturbation conditions to calibrate our model (Figure 3). The structure of our model was constructed to comprise the minimal number of parameters necessary to explain the data. Thereby, we aimed at a model with fully identifiable parameters that are essential to obtain high predictive power. Parameter identifiability was analyzed by the profile likelihood approach. Applying this method, we could establish a dual negative feedback model of JAK2-STAT5 signaling with structurally and in most cases practically identifiable parameters.
A major bottle-neck in combining signal transduction events with cellular phenotypes is the discrepancy in the time scale and stimuli concentrations that are applied in the different experiments. The sensitivity of biochemical assays to determine phosphorylation events within minutes or hours after stimulation is usually lower than the threshold of sensitivity in assays to determine the physiological response after one or more days. Facilitated by the model, we were able to compute the integrated response of JAK2/STAT5 signaling components for experimentally unaddressable Epo concentrations. Our results demonstrate that the integrated response of pSTAT5 in the nucleus accurately correlates with the experimentally determined survival of CFU-E cells. This provides a quantitative link of the dependency of primary CFU-E cells on pSTAT5 activation dynamics. By correlation analysis, we could identify the early signaling phase (⩽1 h) of STAT5 to be the most predictive for the fraction of surviving cells, which was determined ∼24 h later. Thus, we hypothesize that as a general principle in apoptotic decisions, ligand concentrations translated into kinetic-encoded information of early signaling events downstream of receptors can be predictive for survival decisions 24 h later.
After the first hour of stimulation, it is important to constrain signaling to a residual steady-state level. Constitutive phosphorylation of the JAK2/STAT5 pathway has a crucial role in the onset of polycythemia vera (PV), a disease associated with Epo-independent erythroid differentiation. The two identified transcriptional feedback proteins, CIS and SOCS3, are responsible for adjusting the phosphorylation level of STAT5 after 1 h of stimulation. Since the Epo input signal can vary over a broad range of ligand concentrations, we asked how CIS and SOCS3 can facilitate control of STAT5 long-term phosphorylation levels over the entire physiological relevant hormone concentrations. By using model simulations, we revealed that the two feedbacks are most effective at different Epo concentration ranges. Predicted by our mathematical model, the major role of CIS in modulating STAT5 phosphorylation levels is at low, basal Epo concentrations, whereas SOCS3 is essential to control the STAT5 phosphorylation levels at high Epo doses (Figure 6). As a potential molecular mechanism of this dose-dependent inhibitory effect, we could identify the quantity of pJAK2 relative to pEpoR that increases with higher Epo concentrations. Since SOCS3 can inhibit JAK2 directly via its KIR domain to attenuate downstream STAT5 activation, SOCS3 becomes more effective with the relative increase of JAK2 activation. Hence, CIS and SOCS3 act in a concerted manner to ensure tight regulation of STAT5 responses over the broad physiological range of Epo concentrations.
In summary, our mathematical approach provided new insights into the specific function of feedback regulation in STAT5-mediated life or death decisions of primary erythroid cells. We dissected the roles of the transcriptionally induced proteins CIS and SOCS3 that operate as dual feedback with divided function thereby facilitating the control of STAT5 activation levels over the entire range of physiological Epo concentrations. The detailed understanding of the molecular processes and control distribution of Epo-induced JAK/STAT signaling can be further applied to gain insights into alterations promoting malignant hematopoietic diseases.
Cellular signal transduction is governed by multiple feedback mechanisms to elicit robust cellular decisions. The specific contributions of individual feedback regulators, however, remain unclear. Based on extensive time-resolved data sets in primary erythroid progenitor cells, we established a dynamic pathway model to dissect the roles of the two transcriptional negative feedback regulators of the suppressor of cytokine signaling (SOCS) family, CIS and SOCS3, in JAK2/STAT5 signaling. Facilitated by the model, we calculated the STAT5 response for experimentally unobservable Epo concentrations and provide a quantitative link between cell survival and the integrated response of STAT5 in the nucleus. Model predictions show that the two feedbacks CIS and SOCS3 are most effective at different ligand concentration ranges due to their distinct inhibitory mechanisms. This divided function of dual feedback regulation enables control of STAT5 responses for Epo concentrations that can vary 1000-fold in vivo. Our modeling approach reveals dose-dependent feedback control as key property to regulate STAT5-mediated survival decisions over a broad range of ligand concentrations.