Elucidating the crosstalk mechanism between IFN-gamma and IL-6 via mathematical modelling
© Qi et al; licensee BioMed Central Ltd. 2013
Received: 19 September 2012
Accepted: 3 February 2013
Published: 6 February 2013
Interferon-gamma (IFN-gamma) and interleukin-6 (IL-6) are multifunctional cytokines that regulate immune responses, cell proliferation, and tumour development and progression, which frequently have functionally opposing roles. The cellular responses to both cytokines are activated via the Janus kinase/signal transducer and activator of transcription (JAK/STAT) pathway. During the past 10 years, the crosstalk mechanism between the IFN-gamma and IL-6 pathways has been studied widely and several biological hypotheses have been proposed, but the kinetics and detailed crosstalk mechanism remain unclear.
Using established mathematical models and new experimental observations of the crosstalk between the IFN-gamma and IL-6 pathways, we constructed a new crosstalk model that considers three possible crosstalk levels: (1) the competition between STAT1 and STAT3 for common receptor docking sites; (2) the mutual negative regulation between SOCS1 and SOCS3; and (3) the negative regulatory effects of the formation of STAT1/3 heterodimers. A number of simulations were tested to explore the consequences of cross-regulation between the two pathways. The simulation results agreed well with the experimental data, thereby demonstrating the effectiveness and correctness of the model.
In this study, we developed a crosstalk model of the IFN-gamma and IL-6 pathways to theoretically investigate their cross-regulation mechanism. The simulation experiments showed the importance of the three crosstalk levels between the two pathways. In particular, the unbalanced competition between STAT1 and STAT3 for IFNR and gp130 led to preferential activation of IFN-gamma and IL-6, while at the same time the formation of STAT1/3 heterodimers enhanced preferential signal transduction by sequestering a fraction of the activated STATs. The model provided a good explanation of the experimental observations and provided insights that may inform further research to facilitate a better understanding of the cross-regulation mechanism between the two pathways.
The Janus kinase/signal transducer and activator of transcription (JAK/STAT) pathway can be activated by a wide range of cytokines. The binding of cytokines to their receptors drives receptor dimerization and phosphorylation, which leads to the recruitment, activation and dimerization of STATs [1-5]. STAT homodimers are translocated to the nucleus, where they regulate the expression of target genes [6, 7]. IFN-gamma was first identified in a mitogen-activated lymphocyte supernatant based on its distinctive antiviral activity. It is widely involved in protection against tumour development and cancer immunoediting . IFN-gamma triggers the prolonged activation of the transcription factor STAT1 via the IFN-gamma receptor (IFNR) and JAK, which induces target gene expression by binding to the gamma-activated sequences (GAS) in the promoters of these genes . Thus, the activation of the IFN-gamma/JAK/STAT1 pathway can prevent the expansion of many different normal and neoplastic cells . However, IL-6 also has important roles in triggering the acute phase response of the body to injury or inflammation . The receptor involved with the recognition of IL-6 can be subdivided into the non-signalling part, i.e., glycoprotein 80 (gp80), and the signalling part, i.e., glycoprotein 130 (gp130). IL-6 attaches to gp80 and drives IL6-gp80 complex binding to the gp130-JAK complex, which forms the IL6-gp80-gp130-JAK complex. gp130 associates with JAK and it is tyrosine-phosphorylated in response to IL-6 stimulation, which leads to the activation of the JAK/STAT3 pathway and MAPK (mitogen-activated protein kinase) cascades [2, 9]. IL-6 activates target genes involved with differentiation, survival, apoptosis and proliferation, and it plays important roles in pro- and anti-inflammatory functions, acute phase and immune responses in the organism, and tumour progression [2, 6]. Three types of negative regulators are involved in the regulation of the IFN-gamma and IL-6 pathways: the suppressor of cytokine signalling (SOCS), SH2 domain-containing tyrosine phosphates 2 (SHP-2) and various cytoplasmic and nuclear protein tyrosine phosphates (PPs) [10, 11]. SOCS1 and SOCS3 are induced by JAK/STAT pathway and they bind to the activated receptors of IFN-gamma and IL-6, respectively, which negatively regulate the signal transduction of IFN-gamma and IL-6 [12, 13]. SHP-2 acts as a phosphatase in activated receptor complexes of IFN-gamma and IL-6 and it negatively regulates the activation of STATs [14, 15]. PP1 dephosphorylates STAT*s (* represents species-activation states in this study) in the cytoplasm, while PP2 dephosphorylates STAT*s in the nucleus, which results in STATs being returned to the cytosol, thereby influence the activation of STATs [16, 17].
Previous studies have indicated that IFN-gamma and IL-6 have opposing roles in cell proliferation, apoptotic death and inflammation, which are closely related to the specific patterns and duration of STAT activation after their stimulation . IFN-gamma mainly phosphorylates STAT1, which has many pro-inflammatory effects [8, 18, 19]. By contrast, IL-6 is a potent activator of STAT3, which contributes to its anti-inflammatory functions [2, 20, 21]. However, the detailed molecular mechanism leading to the unbalanced activation of STATs after IFN-gamma and IL-6 stimulation remains unclear. Qing et al. suggested that tyrosine 419 in the IFN-gamma receptor subunit 1(IFNGR1) is required for the activation of both STAT1 and STAT3 . In response to IL-6, STAT3 binds to phosphorylated (p)YXXQ motifs (Y767,Y814,Y905 and Y915) of gp130, whereas STAT1 is recruited to a more restricted consensus sequence pYXPQ (Y905 and Y915) in gp130 . In addition, researchers have provided some interesting experimental results using STAT-deficient cells. Qing et al. showed that the activation of STAT3 in response to IFN-gamma was much stronger and more prolonged in STAT1-null cells than wild-type cells . Costa-Pereira et al. showed that IL-6 mediated an IFN-gamma-like response in mouse embryo fibroblasts lacking STAT3, including the prolonged activation of STAT1 and it promoted the induction of multiple IFN-gamma-inducible genes . However, Regis et al. reported that the activation of STAT1 in human neoplastic T lymphocytes after IFN-gamma stimulation was generally unaffected by STAT3 silencing . Ho et al. argued that STAT3 did not affect the activation of STAT1 . Bluyssen et al. reported that pre-treatment of EC with IFN-gamma significantly reduced the activation of STAT3 after induction by IL-6, but without affecting the total amounts of STAT3 . However, Kaur et al. reported that the activation of STAT1 by IFN-gamma was mainly unaffected after pre-treatment with IL-6 or other gp130-related cytokines in SH-SY5Y human neuroblastoma cells . These experimental results indicate that STAT1 and STAT3 may have common binding sites within IFN-gamma and IL-6 receptors, while the activation of STAT3 may depend on the concentration of STAT1 and vice versa. Moreover, the interactions between IFN-gamma and IL-6 signals are not symmetric.
Thyrell et al. also reported that IFN-alpha could influence the signal response of IL-6 in multiple myeloma, which resulted in a decrease in STAT3 homodimer DNA-binding activity and a shift from STAT3 homodimers to STAT1/3 heterodimers . Herrero et al. observed that pre-treatment with IFN-gamma could affect the signal response of IL-10 in macrophages, which caused the IL-10 mediated STAT activation pattern to switch from STAT3 homodimers to STAT1/3 heterodimers . Therefore, changes between STAT homodimers and STAT1/3 heterodimers may represent a biologically relevant approach to determining the crosstalk between IFN-gamma and IL-6 pathways. However, how the formation of STAT1/3 heterodimers regulates the interactions between IFN-gamma and IL-6 signals is not fully understood.
Systems biology modelling generally aims to find fairly plausible mechanistic models that include all of the relevant key processes in a biochemical system. It is considered to be a powerful analytical approach to understanding the essential mechanisms of the physiological functions of normal tissues and pathological progression during complex diseases . The first systems biology model of the IFN-gamma/JAK/STAT1 pathway was developed by Yamada et al. . They modelled activation of the JAK/STAT1 pathway in response to IFN-gamma and analysed the effects of the feedback loop regulated by SOCS1. Later, Zi et al. conducted a multi-parametric sensitivity analysis of the model produced by Yamada et al. and indicated that the concentrations of SOCS1, nuclear phosphatase PP2 and cytoplasmic STAT1, as well as some of the reaction steps that affected those concentrations, were the most sensitive to perturbation . The first model of the IL-6/JAK/STAT3 pathway was produced by Singh et al. . Recently, Moya et al. proposed an updated model of IL-6 and IL-10 signalling via JAK/STAT and ERK-C/EBPβ activation . The model was used to investigate dynamical features of the system such as the activity ratio of JAK/STAT and ERK-C/EBPβ with different stimulation levels of IL-6 and IL-10. The dynamic behaviours of some individual molecules, such as STATs and SOCSs, in the IFN-gamma and IL-6 pathways were investigated in previous studies [28, 30, 31], but signalling of the crosstalk during signal transduction by IFN-gamma and IL-6 has still not been modelled.
In this study, we developed a crosstalk model of the IFN-gamma and IL-6 pathways by combining previously established mathematical models and by comprehensively analyzing the interactions between the two pathways. The model considered three possible levels of crosstalk between the two pathways: (1) the competition between STAT1 and STAT3 for IFNR and gp130; (2) the mutual negative regulation between IFN-gamma and IL-6 via the regulators SOCS1 and SOCS3; and (3) the restrictive effects of the formation of STAT1/3 heterodimers on the activation of the transcription factors STAT1 and STAT3. We considered a number of protocols where cells were stimulated by IFN-gamma and/or IL-6. The simulation results showed that the model provided a good explanation of the experimental observations and it provided new insights that could inform further research to facilitate a better understanding of the cross-regulation between the IFN-gamma and IL-6 pathways.
Based on the model of the IFN-gamma/JAK/STAT1 pathway produced by Yamada et al.  and the model of the IL-6/10/JAK/STAT3 pathway produced by Moya et al. , we established a crosstalk model of the IFN-gamma and IL-6 pathways. A schematic diagram of the model is shown in Additional file 1: Figure S1. In this model, the components of the two previous mathematical models, their structures and most of the parameters were left unchanged. For simplicity, we specified that SHP-2 could repress the activated receptors of IFN-gamma and IL-6, while PP1 and PP2 (PPX and PPN in Yamada et al.’s work) could dephosphorylate STAT1 and STAT3 in the cytoplasm and the nucleus, respectively. We removed any reactions and components that were not connected with IFN-gamma and IL-6 signalling, such as IL-10. Sixteen new reactions were added based on the possible mechanisms of cross-regulation between IFN-gamma and IL-6.
The structure of the STAT1 and STAT3 proteins contains an oligomerization domain, a coiled-coil domain, a DNA-binding domain, a linker domain, an SH2 domain and a transactivation domain [32, 33]. The recruitment of STAT1 and STAT3 to the activated receptor complexes is known to be mediated by their SH2 domains and phosphorylation of the receptor tyrosine motifs is required [34-36]. After receptor binding, the STATs are phosphorylated on a single tyrosine residue (Tyr 701 in STAT1 and Tyr705 in STAT3) [37, 38]. Many experimental observations have shown that STAT1 and STAT3 may combine with the same docking sites in IFNR and gp130 [2, 19]. Thus, we hypothesized that STAT1 and STAT3 might compete for the same phosphorylated docking sites in IFNR and gp130, via their SH2 domains in our model. After STAT1 and STAT3 combine with the activated receptors complexes via IFN-gamma and/or IL-6, they are phosphorylated and disassociate from the receptors. Based on these considerations, the new biochemical reactions (N1) - (N3) and (N4) - (N6) (see “New biochemical reactions added to the crosstalk model” in Additional file 1) were added to our model to simulate the activation of STAT3 after IFN-gamma stimulation and the activation of STAT1 after IL-6 stimulation, respectively.
The kinetic parameters of these new reactions are important because they reflect the properties of the biological system. Wiederkehr-Adam et al. indicated that the SH2 domain of STAT1 had a much higher affinity for the phosphotyrosine 419 motif in IFNGR1 than that in STAT3 . After IL-6 stimulation, STAT3 binds to the flexible pYXXQ motifs in gp130, whereas STAT1 is recruited to the more restricted consensus sequence of pYXPQ in gp130 . Based on these observations, we hypothesized the unbalanced competitive binding of STAT1 and STAT3 with IFNR and gp130 after IFN-gamma and IL-6 stimulation, respectively. Additional file 1: Tables S1-S3 show that the main effector of IL-6 signalling, STAT3, had a higher affinity for gp130 than STAT1. Similarly, the main effector of IFN-gamma signalling, STAT1, had a higher affinity for IFNR than STAT3.
SHP-2 and SOCS combine to regulate signal transduction by IFN-gamma and IL-6 [1, 12, 40]. SOCS1 inhibits the JAK/STAT pathway by binding to the activation loop of JAK via its SH2 domain . SOCS3 can also bind to JAK . SOCS1 and SHP-2 combine with different sites in the receptor complexes of IFN-gamma . However, SOCS3 and SHP-2 may have similar binding specificities. Experiments have suggested that SOCS3 and SHP2 may compete for same site (tyrosine Y759) in gp130 after IL-6 stimulation [42, 43]. In our model, SOCS1 and SHP-2 were capable of binding to the receptor complex of IFN-gamma without mutual interference, whereas SOCS3 and SHP-2 could competitively bind to the receptor complex. Distinct genes belonging to the SOCS family are induced as immediate early genes (IEGs) downstream of different STATs and they can inhibit STAT activation in a classical negative feedback loop . It is generally recognized that SOCS1 has an important role in modulating IFN-gamma signalling , whereas SOCS3 mainly affects IL-6 signalling . It is also well known that the STAT1 and STAT3 homodimers are direct transcription factors of the JAK/STAT pathway, which play important roles in signal transduction during IFN and gp130 receptor signalling . However, Bluyssen et al. reported that SOCS3 could also be induced after IFN-gamma stimulation in EC and that it could inhibit signal transduction by IL-6 . Qing et al. observed that the STAT1 and STAT3 homodimers could both be induced after IFN-gamma stimulation in MEFs, which both bound to the same GAS element in the SOCS3 promoter . The sequence of this GAS element is conserved in mice, rats and humans . It was shown that STAT3 activation was much stronger and more prolonged in STAT1-null cells, and that SOCS3 was strongly induced in wild-type and STAT1-null cells, while the levels of SOCS3 mRNA were greatly increased in STAT3-null cells . Thus, it is speculated that STAT1 homodimers might also promote the transcription of SOCS3 in the same way as STAT3 homodimers. However, no experimental evidence indicates that STAT3 homodimers can combine with the promoter region of SOCS1. Thus, our model does not regard STAT3 homodimers as an efficient transcription factor for SOCS1. We added equation (N7) to our model to simulate the transcription of SOCS3 mRNA after its induction by STAT1 homodimers in the nucleus, which is represented as (STAT1N*)2.
Thyrell et al. reported that IFN-alpha could affect the signal response of IL-6 in multiple myeloma, which resulted in a decrease in STAT3 homodimer DNA-binding activity and a shift from STAT3 homodimers to STAT1/3 heterodimers . Herrero et al. showed that pre-treatment with IFN-gamma could affect the signal response of IL-10 in macrophages, causing the IL-10-mediated activation pattern to switch from STAT3 homodimers to STAT1/3 heterodimers . These experimental results showed that STAT1/3 heterodimers play important roles in the crosstalk between different cytokines. The activation of STATs after cytokine stimulation led to the formation of STAT homo- and heterodimers [2, 6].
Haan et al. reported that IL-6 stimulation of primary human macrophages led to a different distribution of STAT dimer species in the cytosol and nucleus. In particular, they showed that STAT1/3 heterodimers were present in the cytosol and nucleus . The size of STATs exceeds 90 kDa, which is far beyond the exclusion limit of the nuclear pores, so STATs need to be translocated actively into the nucleus [2, 6]. Tyrosine phosphorylation is not necessarily required for STAT nuclear translocation . Basic residues (Lys410 and Arg413) are known to contribute to the nuclear localization signals (NLSs) of dimeric STAT1 . PP1 dephosphorylates STAT*s in the cytoplasm, while PP2 dephosphorylates STAT*s in the nucleus, which leads to STATs being returned to the cytosol [16, 17, 51]. It was postulated that the nuclear export signals (NESs) in the DNA-binding domain of STAT1 is comprised of the residues 399-410 . The formation and dephosphorylation of STAT1 and STAT3 homodimers in the cytosol and nucleus were modelled by Yamada et al. and Moya et al., respectively. Their models also described the translocation of STAT1 and STAT3 homodimers from the cytosol to the nucleus, and the export of dephosphorylated STAT monomer from the nucleus to the cytoplasm [28, 31]. In our model, biochemical reactions (N8) - (N9) were added to simulate the formation of STAT1/3 heterodimers in the cytoplasm and nuclei. For simplicity, it was supposed that only STAT1* and STAT3* could form STAT1/3 heterodimers. Biochemical reaction (N10) was added to simulate the STAT1/3 heterodimer translocation process from the cytoplasm to the nucleus based on the translocation of STAT homodimers. It was also assumed that STAT1/3 heterodimers could be dephosphorylated by PP1 and PP2, which resulted in STAT export from the nucleus to the cytoplasm. Biochemical reactions (N11)-(N16) were added to simulate this process. The specific biological role of STAT1/3 heterodimers remains obscure , so we did not regard the STAT1/3 heterodimers as efficient transcription factors in our model.
To construct the crosstalk model, we merged the common components from previous models, such as SHP-2, PP1 and PP2. JAK1 and JAK2 are two species in the JAK family, which play important roles in the signal responses of IFN-gamma and IL-6 . In previous studies, JAK1 and JAK2 were treated as JAK for simplicity [28, 31]. JAK1 and JAK2 can combine with the receptors of IFN-gamma and IL-6, but the signal transduction activities of IFN-gamma and IL-6 may relate to specific types of JAKs. The internal membrane proximal regions of JAK1 and JAK2,which responds to the IFN-gamma signal, bind the IFN-gamma receptor subunits IFNGR1 and IFNGR2 . During IL-6 signalling, JAK1 and JAK2 are activated via the conserved membrane proximal-binding domain of the receptors, and JAK1 plays a major role in the signal response to IL-6 . Moreover, Guschin et al. suggested that, although JAK2 was activated, it could not mediate the efficient activation of STAT1 and STAT3 after IL-6 stimulation in the absence of JAK1 . Thus, JAK2 may share redundant functions with JAK1. In our model, we used two different JAK species for IFN-gamma and IL-6 receptors, respectively.
Our model contains two main components: IL-6 signalling via the JAK/STAT3 pathway and IFN-gamma signalling via the JAK/STAT1 pathway. There are multi-level interactions between the two pathways. The model contains 108 species, 192 kinetic parameters and 119 reactions, of which 103 reactions are based on previous models whereas 16 reactions were new. To make the experimental results easier to compare, if not specified otherwise, the concentrations of STAT1* and STAT3* were the sum of the concentrations of all species containing activated STAT1 and STAT3, respectively, including their monomers and dimers.
Responses of the crosstalk model to separate IFN-gamma and IL-6 stimulation
SHP-2 knockout simulations were also performed to characterize the effects of SHP-2. As shown in Figure 2A-B, knocking out SHP-2 enhanced the signal responses of IFN-gamma and IL-6, which agreed with previous experimental observations. You et al. showed that IFN-gamma could induce a higher signal response in SHP-2 null cells . Schapter et al. also reported that over-expression of an inactive SHP-2 mutant in HepG2 cells enhanced STAT activation after IL-6 stimulation . After IFN-gamma or IL-6 stimulation, however, the JAK/STAT pathway exhibited different features to those when knocking out SHP-2. Without SHP-2, IFN-gamma stimulation induced higher levels of STAT1* and STAT3* (about 150%) than that in normal conditions (Figure 2A). By contrast, IL-6 stimulation induced rapid increases in STAT1* and STAT3*, the maximum concentrations of which reached 900 nM and 500 nM, respectively, which was about three times higher than that in normal conditions. After IL-6 stimulation, we also observed that SOCS3 reached a peak value about 75 nM at 2 h, which inhibited signal transduction by IL-6 and quickly caused the concentrations of STAT1* and STAT3* to drop to normal levels after 3 h (Figure 2B). Knockout simulations were also performed for SOCS1 and SOCS3. As shown in Figure 2C-D, knocking out SOCS1 enhanced the activation of STAT1 after IFN-gamma stimulation, while knocking out SOCS3 enhanced the activation of STAT3 after IL-6 stimulation. Our simulation results agreed with previous experimental observations. Brysha et al. demonstrated that in vitro and in vivo hepatocytes lacking SOCS-1 exhibited prolonged activation of STAT1 after IFN-gamma stimulation, which correlated with the dramatically increased sensitivity to the toxic effects of IFN-gamma . Niwa et al. reported that inhibition of SOCS3 expression enhanced the activation of STAT3 and cell growth . After IFN-gamma or IL-6 stimulation, however, the JAK/STAT pathway exhibited different features when knocking out SOCS1 or SOCS3. Without SOCS1, IFN-gamma stimulation induced higher levels of STAT1* and STAT3* (about 250%) compared with those in normal conditions (Figure 2C). Without SOCS3, however, IL-6 stimulation induced increases in STAT1* and STAT3*, the maximum concentrations of which reached 700 nM and 300 nM, respectively, which were about double those in normal conditions. After IL-6 stimulation, we also observed that SHP-2 dropped to a low level of about 80 nM at 1 h, which attenuated signal transduction by IL-6 and caused the concentrations of STAT1* and STAT3* to fall slowly after 3 h (Figure 2D). Thus, SOCS1 and SOCS3 had different interaction patterns with SHP-2. SOCS1 and SHP-2 synergistically regulated signal transduction by IFN-gamma. Knocking out SOCS1 or SHP-2 enhanced the integral activation of STAT1 induced by IFN-gamma stimulation. By contrast, SOCS3 and SHP-2 regulated signal transduction by IL-6 in a more complementary manner. Knocking out SHP-2 alone enhanced the rapid response of the IL-6 signal, due to a compensatory increase in SOCS3. Knocking out SOCS3 also led to lower levels of SHP-2, which caused a slow decline in STAT1* and STAT3* 3 h after IL-6 stimulation. Simulations of the combined knockout of SHP-2 and SOCS were performed to characterize their joint effects on IFN-gamma and IL-6 stimulations. First, we stimulated the SHP-2 and SOCS1 combined knockout model with IFN-gamma (0.1 nM) for 12 h and found that STAT1* reached its maximum concentration (960 nM) within about 2 h while STAT3* reached its maximum concentration (900 nM) within about 3 h (Figure 2E). IFN-gamma stimulation induced a similar strength of STAT1* and STAT3* in SHP-2 and SOCS1 combined knockout conditions. We then stimulated the SHP-2 and SOCS3 combined knockout model using IL-6 (0.1 nM) for 12 h and found that STAT1* and STAT3* both quickly reached their maximal concentration of 950 nM and 980 nM, respectively, within about 1 h (Figure 2F). IL-6 stimulation also induced similar strengths of STAT1* and STAT3* in SHP-2 and SOCS3 combined knockout conditions. Thus, the combined knockout of SHP-2 and SOCSs abolished the preferential activation of IFN-gamma and IL-6. The unbalance competition between STAT1 and STAT3 was not related directly to SHP-2 and SOCSs, but SHP-2 and SOCSs combined with the activated receptor complexes and inhibited signal transduction via the JAK/STAT pathway [2, 6, 12]. Therefore, we deduced that SHP-2 and SOCSs could limit the concentration of active receptor complexes, which indirectly affected the preferential activation of IFN-gamma and IL-6. Therefore, we investigated the signal transduction profiles of the activated receptor complexes in response to IFN-gamma and IL-6 with and without knocking out SHP-2 and/or SOCSs. Without any knockout, (IFN-R-JAK)2* reached its maximum concentration (0.1 nM) in about 0.5 h after IFN-gamma, before decreasing rapidly. After knocking out SHP-2, the level of (IFN-R-JAK)2* increased by about 120% compared with that in normal conditions. Without SOCS1, (IFN-R-JAK)2* increased rapidly and reached a new steady state (0.35 nM) after 2 h, whereas the combined knockout of SHP-2 and SOCS1 caused the level of (IFN-R-JAK)2* to increase significantly, reaching 4.5 nM in 12 h, which was about forty time as high as that in normal conditions (Figure 2G). Without any knockout and with IL-6, (IL6-gp80-gp130-JAK)2* reached its maximum concentration (0.15 nM) within about 0.25 h, before decreasing rapidly. After knocking out SHP-2, (IL6-gp80-gp130-JAK)2* increased rapidly and it reached its maximum concentration of 1.4 nM in 0.25 h, which was about nine times that in normal conditions, although it quickly returned to a normal level after 0.5 h. With SOCS3 knock out, the (IL6-gp80-gp130-JAK)2* level increased and reached a new steady state (0.15 nM) after 1 h. With the combined knockout of SHP-2 and SOCS3, the levels of (IL6-gp80-gp130-JAK)2* increased significantly and reached a new steady state (5 nM) after 1 h, which was about 35 times that in normal conditions (Figure 2H). The simulation results demonstrated that with the SHP-2 and SOCSs combined knockout, the levels of (IFN-R-JAK)2* and (IL6-gp80-gp130-JAK)2* increased significantly after IFN-gamma and IL-6 stimulation. STAT1 and STAT3 competed for the same motifs in IFNR and gp130, but there was sufficient (IFN-R-JAK)2* and (IL6-gp80-gp130-JAK)2*, so the preferential activations of IFN-gamma and IL-6 were abolished. These simulated observations still await further experimental verification.
Responses of the crosstalk model after disrupting STAT1 and STAT3
After IFN-gamma and IL-6 stimulation, the recruitment of STAT1 and STAT3 to the activated receptor complexes directly affected their phosphorylation, which has important roles in signal transduction by IFN-gamma and IL-6 [34-36]. In our unbalanced competition model, STAT1 and STAT3 had different affinities for IFNR and gp130, so we deduced that disrupting STAT1 and STAT3 may have different effects on the recruitment of STAT1 and STAT3. Next, we investigated the effect of changing the initial concentration of STAT1 and STAT3 on the associations of STATs with activated receptor complexes in response to IFN-gamma and IL-6. Our simulation results demonstrated that changing the concentration of STAT3 had little effect on the formation of (IFN-R-JAK)2*-STAT1 (Figure 3E), whereas altering the STAT1 level significantly affected the formation of (IFN-R-JAK)2*-STAT3 (Figure 3F) after IFN-gamma stimulation. Our simulations also showed that the formation of (IL6-gp80-gp130-JAK)2*-STAT3 was almost independent of STAT1 disruption after IL-6 stimulation, although changing the STAT3 level significantly affected the formation of (IL6-gp80-gp130-JAK)2*-STAT1(data not shown). Therefore, the unbalanced competition between STAT1 and STAT3 for IFNR and gp130 was not only the pivotal mechanism for the preferential activation of IFN-gamma and IL-6, but it also determined the recruitment of STAT1 and STAT3 to the activated receptor complexes.
Responses of the crosstalk model to combined stimulation with IFN-gamma and IL-6
Signal transduction via the JAK/STAT pathway depended on the formation of STAT homodimers , which are regarded as the main transcription factors during IFN-gamma and IL-6 signalling [2, 8]. We further investigated whether combined stimulation with IFN-gamma and IL-6 could induce higher (STAT1N)*2 and (STAT3N)*2 than separate treatments. Figure 4A shows that 12 h after combined stimulation, the STAT homodimers were not induced at a higher level than the separate treatments. However, STAT1N*-STAT3N* reached their maximum concentration (200 nM) within about 1 h, which was about 3 time higher than the individual treatment. The formation of STAT1N*-STAT3N* greatly restricted the formation of STATs homodimers. After we abolished the formation of STAT1/3 heterodimers, the maximum concentrations of (STAT3N)*2 and (STAT1N)*2 increased to about 100 nM with combined stimulation (Figure 4E). Combined stimulation with IFN-gamma and IL-6 led to greater activation of both STAT1 and STAT3, but the formation of STAT1/3 heterodimers played an important role in preventing mutual strengths between IFN-gamma and IL-6 signalling.
Responses of the crosstalk model to successive IFN-gamma and IL-6 stimulation
We analyzed previous studies that focused on the interactions between IFN-gamma and IL-6 signalling and found that their interactions were asymmetric. Bluyssen et al. reported that pre-treatment of EC with IFN-gamma significantly decreased STAT3* induction by IL-6 without affecting the total amount of STAT3 . By contrast, Kaur et al. reported that STAT1 activation induced by IFN-gamma was mainly unchanged after pre-treatment IL-6 or other gp130-related cytokines in SH-SY5Y human neuroblastoma cells . We tried to provide a reasonable explanation for the asymmetric interactions between IFN-gamma and IL-6 using simulation experiments with our model.
The formation of STAT1/3 heterodimers also contributed to the asymmetric interactions between IFN-gamma and IL-6 signalling. As previously described, the formation of STAT1/3 heterodimers enhanced the preferential signal transduction of IFN-gamma and IL-6 by sequestering a fraction of STAT1* and STAT3*. After we abolished the formation of STAT1/3 heterodimers, the maximum concentration of SOCS1 induced by IL-6 increased to 4.2 nM within about 1.5 h, and IL-6 exhibited a greater capacity for inhibiting IFN-gamma signalling (Figure 5F). Abolishing the formation of STAT1/3 heterodimers also enhanced the inhibition from IFN-gamma to IL-6 (Figure 5G). In addition, the mechanism of (STAT1N*)2 inducing SOCS3 also played an important role in the asymmetric interactions. The concentration of (STAT3N*)2 induced by IFN-gamma stimulation was very low due to the sequestering effect of STAT1/3 heterodimers. Therefore, we deduced that SOCS3 induction by (STAT3N*)2 was not sufficient to achieve the inhibition from IFN-gamma to IL-6. Indeed, when we abolished the (STAT1N)*2 induction of SOCS3, the inhibition from IFN-gamma to IL-6 was clearly mitigated (Figure 5H). IFN-gamma only slightly reduced the activation of STAT3 induction by IL-6, which did not agree with previous experimental observations . Therefore, (STAT1N*)2 induction of SOCS3 was essential for the inhibition from IFN-gamma to IL-6 during signal transduction.
Argument for the non-competitive model
Ho et al. reported that STAT3 did not affect the activation of STAT1 during type I interferon responses . They argued that STAT3 did not suppress STAT1 activation via tyrosine phosporylation and they excluded inhibitory mechanisms such as the competition for docking sites and the inhibition of signalling events upstream of STAT1 tyrosine phosphorylation. Based on their argument, we proposed an alternative hypothesis for the combined pattern between STATs and the receptors for IFN-gamma and IL-6 in addition to our previous competition model, which we refer to as the non-competitive model. In the non-competitive model, STAT1 and STAT3 do not compete for the same docking sites of IFNR and gp130, so we considered that STAT1 and STAT3 could bind different phosphorylated docking sites in IFNR and gp130 via their SH2 domain. Based on these considerations, we added the new biochemical reactions (N17) - (N20) (see “non-competitive model” in Additional file 1) to simulate the mutually independent combination of STAT1 and STAT3 with the receptor complexes. STAT3 did not suppress the activation of STAT1 via tyrosine phosphorylation, so we considered that the combination of STAT3 with receptor complexes did not suppress the phosphorylation processes of STAT1 and vice versa. Based on these considerations, we added the new biochemical reactions (N21) - (N24) to simulate the mutually independent phosphorylation of STAT1 and STAT3 within receptor complexes. In addition, we estimated the parameters for reactions (N17)-(N24) and re-estimated the parameters for reactions (N1)-(N6) to fit the non-competitive model (Additional file 1: Table S5). Detailed descriptions of the non-competitive model written in COPSI (using the filename extension “.cps”) are available in Additional file 2.
Next, we executed the same series simulation experiments described in the previous sections using the non-competitive model. We found that the non-competitive model approximately simulated the crosstalk between IFN-gamma and IL-6 signalling in all conditions (data not shown), with the exception of the disrupted STAT1 and STAT3 conditions. Figure 6E shows that when we fixed the concentration of STAT1 at 1000 nM and changed the concentration of STAT3 from 0, 500, 1000, to 2000 nM, the states of STAT1* were obviously dependent on the level of STAT3 after IFN-gamma stimulation. In particular, STAT1 was more strongly phosphorylated, and for longer, with higher levels of STAT3, which were completely different from the experimental observations reported by Regis et al. . Similarly, when we fixed the concentration of STAT3 at 1000 nM and changed the concentration of STAT1 from 0, 500, 1000, to 2000 nM, the states of STAT3* were also obviously dependent on the level of STAT1 after IL-6 stimulation (Figure 6F). Moreover, STAT3 was also more strongly phosphorylated and for longer with higher levels of STAT1, which did not agree with the experimental observations of Dimberg et al. . Indeed, the activation of STATs was determined by the swing between phosphorylation and dephosphorylation processes . In the non-competitive model, STAT1 and STAT3 did not suppress the phosphorylation of each other, although they share the phosphatases PP1 and PP2, so the activation of STAT1 prevented STAT3 from being dephosphorylated and vice versa. The disruption of the dynamic balance between the phosphorylation and dephosphorylation of STAT1 and STAT3 meant the responses of the non-competitive model were hypersensitive to the initial concentrations of STAT1 and STAT3. Therefore, the non-competitive model may not reflect the crosstalk mechanisms between IFN-gamma and IL-6 signalling in a physiological environment. Thus, alternative hypotheses should be proposed based on new experimental observations. Here, we provide an integrated platform that facilitates the verification of possible hypotheses related to the crosstalk between IFN-gamma and IL-6 signalling in future work.
Sensitivity analysis of the competition model
In this study, we used computational simulation to demonstrate the effectiveness of the competition model. Next, we determined the dynamic characteristics of our model by sensitivity analysis. The calculations used in the sensitivity analysis are shown in the “Methods” section. The upstream key regulatory molecules of (STAT1N*)2 and (STAT3N*)2 comprise cytokines, receptors, JAK, STAT1C, STAT3C, PP1, PP2, and SHP2, and their initial concentrations were considered in the sensitivity analysis. First, we applied a sensitivity analysis to determine the relative sensitivity of these components using IFN-gamma stimulation as the input.
Next, we performed a sensitivity analysis using IL-6 stimulation as the input. The sensitivities of the initial concentrations and the critical kinetic parameters are shown in Figure 7C and Figure 7D. The results indicated that the reactions with PP2 negative regulation and JAK combination were highly sensitive to IL-6 stimulation. As described in previous sections, the concentration of STAT1 did not significantly affect the activation of STAT3 after IL-6 stimulation, which was also confirmed by the sensitivity analysis. Similarly, we found that the SOCS3 synthesis, degradation, and negative regulation processes were highly sensitive to IL-6 stimulation, which confirmed the important role of the negative feedback factor SOCS3 during IL-6signal transduction. In particular, we found that the phosphorylation of STAT3 in STAT1/3 heterdimers in the nucleus also had a high impact on the state of (STAT3N*)2. Overall, the sensitivity analysis determined the sensitive components and parameters during JAK/STAT signal transduction. These results provide valuable information that may inform further investigations of regulation and drug target identification in aberrant pathways. The detailed results of the sensitivity analysis of the competition model are shown in Additional file 1: Tables S6-S7.
Sensitivity analysis of the non-competitive model
The results were affected by the limitations of our model. It is important to consider these limitations because they form the basis for future improvements. First, we only considered three possible levels of crosstalk between IFNγ and IL-6 systems. Haan et al. showed that, after IL-6 stimulation, STAT1 phosphorylation was entirely dependent on JAK1 kinase activity whereas STAT3 activation was not . Qing et al. showed that, in response to IFN-gamma, SCR-family kinases had to activate STAT3 (but not STAT1) via tyrosine phosphorylation . Some mechanisms that contribute to the specific signal responses of IFN-gamma and IL-6 were not considered in this study. Second, our computational model was a simplification of a biological process and it might not reflect the true regulatory mechanism for the signal interactions. Third, the simulation results in this study were affected by assumptions made based on our experience, e.g., the STAT1/3 heterodimers could be dephosphorylated by PP1 and PP2 in the cytoplasm and nucleus, respectively. Fourth, our model had very rich dynamics parameters that were not fully explored in this study. For simplicity, protein turnover, receptor recycling, de novo synthesis and the degradation of transcription factors were not included in the current model. The ultimate aim of our research is to construct a universal model that accurately reflects the crosstalk between IFN-gamma and IL-6 signals. However, the signal responses to cytokines may be cell type-dependent. Caldenhoven et al. reported that the activation of STAT3 induced by IFN-gamma was lineage-specific in human neutrophils . Zhang et al. showed that IL-6 stimulation could induce STAT1 phosphorylation in a dose- and time-dependent manner in M1, R2 and U937 cells, although it had little effect on STAT1 phosphorylation in 7TD1 and TF1 cells . Bluyssen et al. argued that IL-6 did not activate STAT1 in EC . Furthermore, the concentrations of molecules within the JAK/STAT pathway, such as STATs, are cell type-dependent . The experimental data are limited and not systematic, so we had to create our model based on experimental observations of many different cell types, while we neglected certain contradictory experimental observations. Therefore, it is not possible to expect that the dynamics predicted by our model will apply universally to all types of cells. Thus, structure and parameters of the model may need some adjustment to reflect signal transduction by IFN-gamma and IL-6 in certain cell types. Finally, our crosstalk model was based on experimentally established interactions, but further experimental verification and improvements are required.
Our simulation results showed that STAT1/3 heterodimers have three important functions during signal transductions from IFN-gamma to IL-6. First, the formation of STAT1/3 heterodimers enhances the preferential signal transduction by IFN-gamma and IL-6 because it sequesters a fraction of STAT1* and STAT3*. Second, the formation of STAT1/3 heterodimers prevents mutual reinforcement between IFN-gamma and IL-6 signalling. Finally, the formation of STAT1/3 heterodimers limits the reciprocal inhibition of IFN-gamma and IL-6 signalling. In our simulations, therefore, the formation of STAT1/3 heterodimers dramatically affected the interaction between the IFN-gamma and IL-6 systems, which suggests that STAT1/3 heterodimers may be a potential target for rectifying abnormal signal transduction by IFN-gamma and IL-6. The functional interference of STAT3 homodimers using STAT3 transcription decoys or small molecules in structure-activity relationship (SAR) studies could successfully inhibit the growth of tumour cells [59, 60]. However, the therapeutic potential of altering the formation of STAT1/3 heterodimers has not been fully investigated. Thus, further research is still required.
This is the first effort to construct a mathematical model of the crosstalk between IFN and IL-6 signalling. Moreover, our simulation results and theoretical findings provide new insights into the dynamical integration of IFN and IL-6 signals. The lack of experimental data and our current superficial understanding of signal transduction mean there are still many defects in our crosstalk model. We will follow current studies related to this issue and improve the quality and veracity of model in our future research.
Based on previous models and new experimental observations, we developed the first crosstalk model of IFN-gamma and IL-6 signalling. This theoretical study successfully reproduced key experimental findings and reached some definitive conclusions. First, the unbalanced competition between STAT1 and STAT3 for IFNR and gp130 led to preferential activation of IFN-gamma and IL-6. At the same time, the formation of STAT1/3 heterodimers enhanced preferential signal transduction by sequestering a fraction of STAT1* and STAT3*. Moreover, SOCSs with SHP-2 limited the concentration of the activated receptor complexes of IFN-gamma and IL-6, which also contributed to the preferential activation of IFN-gamma and IL-6. Second, the unbalanced competition between STAT3 and STAT1 was the pivotal mechanism during the mutual switch between IFN-gamma and IL-6 signals after knocking out STAT1 or STAT3. Finally, the formation of STAT1/3 heterodimers prevented the mutual reinforcement of IFN-gamma and IL-6 signalling, and also limited the reciprocal inhibition between IFN-gamma and IL-6 signalling. Moreover, the process of STAT1 homodimer induction of SOCS3 also contributed to asymmetric interactions between IFN-gamma and IL-6 signalling.
where A C represents a molecule in the cytoplasm; A N is A C translocated from the cytoplasm to nuclei; B is the catalyst; and K f is the rate of the reaction.
New biochemical reactions added to the crosstalk model
where R represents the receptor; C represents a species within the cytoplasm; N represents a species within nuclei; 2 represents dimers; and * represents the activation state. mRNA-SOCS3N represents the SOCS3 mRNA in the nucleus; (STAT1N*)2 represents phosphorylated STAT1 homodimers in the nucleus; Vmax represents the maximum reaction velocity of the SOCS3-mRNA transcription; Km represents half of the substrate concentration when the enzymatic reaction reaches its maximum reaction velocity.
Using this reaction scenario, we established systematic models that were described using ordinary differential equations (ODEs). COPASI (http://www.copasi.org/) and CellDesigner (http://www.celldesigner.org/) were used to calculate the model states and solve the ODE models [61, 62]. Detailed descriptions of the model written in COPSI (using the filename extension “.cps”) are available in Additional file 2.
The sensitivity analysis determined the change in the output after a tiny perturbation in the initial state or the reaction parameters. Thus, a perturbation Δx in the parameter set induced an overall state change Δy and the sensitivity was defined as Δy/Δx, which was normalized with respect to the actual values of the phenomenon investigated and the parameter x/y. Experimental studies have shown that phosphorylated STAT dimers in the nucleus are critical transcription factors for target gene activation in the JAK/STAT pathway. Therefore, our sensitivity analysis considered (STAT1N*)2 and (STAT3N*)2 as the outputs of IFN-gamma and IL-6 signalling, respectively. In this study, Δy was defined as the change in the concentration of (STAT1N*)2 or (STAT3N*)2 after 1.5 h, which adequately represented the signal transduction strength of IFN-gamma or IL-6. The perturbation ranges of each parameter were set as 30% of the actual absolute value to avoid errors in the sensitivity estimates. COPASI and CellDesigner were used to calculate the sensitivity.
This work was supported by the National Natural Science Foundation of China (No. 61172183), the Natural Science Foundation of Jilin Province (No. 20101503), the National High Technology Research and Development Program of China (No.2012AA02A407), the Research Foundation of Jilin Provincial Science & Technology Development (No. 20110711), and the Fundamental Research Funds for the Central Universities (Grant No. 11QNJJ021).
- Schindler C, Levy DE, Decker T: JAK-STAT signaling: from interferons to cytokines. J Biol Chem 2007,282(28):20059-20063. 10.1074/jbc.R700016200View ArticlePubMedGoogle Scholar
- Heinrich PC, Behrmann I, Haan S, Hermanns HM, Muller-Newen G, Schaper F: Principles of interleukin (IL)-6-type cytokine signalling and its regulation. Biochem J 2003,374(Pt 1):1-20.PubMed CentralView ArticlePubMedGoogle Scholar
- Shuai K: Modulation of STAT signaling by STAT-interacting proteins. Oncogene 2000,19(21):2638-2644. 10.1038/sj.onc.1203522View ArticlePubMedGoogle Scholar
- Ikeda H, Old LJ, Schreiber RD: The roles of IFN gamma in protection against tumor development and cancer immunoediting. Cytokine Growth Factor Rev 2002,13(2):95-109. 10.1016/S1359-6101(01)00038-7View ArticlePubMedGoogle Scholar
- Radaeva S, Jaruga B, Kim WH, Heller T, Liang TJ, Gao B: Interferon-gamma inhibits interferon-alpha signalling in hepatic cells: evidence for the involvement of STAT1 induction and hyperexpression of STAT1 in chronic hepatitis C. Biochem J 2004,379(Pt 1):199-208.PubMed CentralView ArticlePubMedGoogle Scholar
- Regis G, Pensa S, Boselli D, Novelli F, Poli V: Ups and downs: the STAT1:STAT3 seesaw of Interferon and gp130 receptor signalling. Semin Cell Dev Biol 2008,19(4):351-359. 10.1016/j.semcdb.2008.06.004View ArticlePubMedGoogle Scholar
- Koster M, Hauser H: Dynamic redistribution of STAT1 protein in IFN signaling visualized by GFP fusion proteins. Eur J Biochem 1999,260(1):137-144. 10.1046/j.1432-1327.1999.00149.xView ArticlePubMedGoogle Scholar
- Schroder K, Hertzog PJ, Ravasi T, Hume DA: Interferon-gamma: an overview of signals, mechanisms and functions. J Leukoc Biol 2004,75(2):163-189.View ArticlePubMedGoogle Scholar
- Eulenfeld R, Dittrich A, Khouri C, Muller PJ, Mutze B, Wolf A, Schaper F: Interleukin-6 signalling: More than Jaks and STATs. Eur J Cell Biol 2011,91(6-7):486-495.View ArticlePubMedGoogle Scholar
- O’Shea JJ, Gadina M, Schreiber RD: Cytokine signaling in 2002: new surprises in the Jak/Stat pathway. Cell 2002,109(Suppl):S121-131.View ArticlePubMedGoogle Scholar
- Shuai K, Liu B: Regulation of JAK-STAT signalling in the immune system. Nat Rev Immunol 2003,3(11):900-911. 10.1038/nri1226View ArticlePubMedGoogle Scholar
- Cooney RN: Suppressors of cytokine signaling (SOCS): inhibitors of the JAK/STAT pathway. Shock 2002,17(2):83-90. 10.1097/00024382-200202000-00001View ArticlePubMedGoogle Scholar
- Brysha M, Zhang JG, Bertolino P, Corbin JE, Alexander WS, Nicola NA, Hilton DJ, Starr R: Suppressor of cytokine signaling-1 attenuates the duration of interferon gamma signal transduction in vitro and in vivo. J Biol Chem 2001,276(25):22086-22089. 10.1074/jbc.M102737200View ArticlePubMedGoogle Scholar
- You M, Yu DH, Feng GS: Shp-2 tyrosine phosphatase functions as a negative regulator of the interferon-stimulated Jak/STAT pathway. Mol Cell Biol 1999,19(3):2416-2424.PubMed CentralView ArticlePubMedGoogle Scholar
- Schaper F, Gendo C, Eck M, Schmitz J, Grimm C, Anhuf D, Kerr IM, Heinrich PC: Activation of the protein tyrosine phosphatase SHP2 via the interleukin-6 signal transducing receptor protein gp130 requires tyrosine kinase Jak1 and limits acute-phase protein expression. Biochem J 1998,335(Pt 3):557-565.PubMed CentralView ArticlePubMedGoogle Scholar
- ten Hoeve J, de Jesus Ibarra-Sanchez M, Fu Y, Zhu W, Tremblay M, David M, Shuai K: Identification of a nuclear Stat1 protein tyrosine phosphatase. Mol Cell Biol 2002,22(16):5662-5668. 10.1128/MCB.22.16.5662-5668.2002PubMed CentralView ArticlePubMedGoogle Scholar
- Levy DE, Darnell JE Jr: Stats: transcriptional control and biological impact. Nat Rev Mol Cell Biol 2002,3(9):651-662.View ArticlePubMedGoogle Scholar
- Bluyssen HA, Rastmanesh MM, Tilburgs C, Jie K, Wesseling S, Goumans MJ, Boer P, Joles JA, Braam B: IFN gamma-dependent SOCS3 expression inhibits IL-6-induced STAT3 phosphorylation and differentially affects IL-6 mediated transcriptional responses in endothelial cells. Am J Physiol Cell Physiol 2010,299(2):C354-362. 10.1152/ajpcell.00513.2009View ArticlePubMedGoogle Scholar
- Qing Y, Stark GR: Alternative activation of STAT1 and STAT3 in response to interferon-gamma. J Biol Chem 2004,279(40):41679-41685. 10.1074/jbc.M406413200View ArticlePubMedGoogle Scholar
- Costa-Pereira AP, Tininini S, Strobl B, Alonzi T, Schlaak JF, Is’harc H, Gesualdo I, Newman SJ, Kerr IM, Poli V: Mutational switch of an IL-6 response to an interferon-gamma-like response. Proc Natl Acad Sci USA 2002,99(12):8043-8047. 10.1073/pnas.122236099PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang JY, Li Y, Ren YF, Shen BF: Comparison of IL-6-induced STAT1 activation in different responsive cells. Xi Bao Yu Fen Zi Mian Yi Xue Za Zhi 2003,19(2):106-108.PubMedGoogle Scholar
- Regis G, Icardi L, Conti L, Chiarle R, Piva R, Giovarelli M, Poli V, Novelli F: IL-6, but not IFN-gamma, triggers apoptosis and inhibits in vivo growth of human malignant T cells on STAT3 silencing. Leukemia 2009,23(11):2102-2108. 10.1038/leu.2009.139View ArticlePubMedGoogle Scholar
- Ho HH, Ivashkiv LB: Role of STAT3 in type I interferon responses. Negative regulation of STAT1-dependent inflammatory gene activation. J Biol Chem 2006,281(20):14111-14118. 10.1074/jbc.M511797200View ArticlePubMedGoogle Scholar
- Kaur N, Kim IJ, Higgins D, Halvorsen SW: Induction of an interferon-gamma Stat3 response in nerve cells by pre-treatment with gp130 cytokines. J Neurochem 2003,87(2):437-447. 10.1046/j.1471-4159.2003.02012.xView ArticlePubMedGoogle Scholar
- Thyrell L, Arulampalam V, Hjortsberg L, Farnebo M, Grander D, Pokrovskaja Tamm K: Interferon alpha induces cell death through interference with interleukin 6 signaling and inhibition of STAT3 activity. Exp Cell Res 2007,313(19):4015-4024. 10.1016/j.yexcr.2007.08.007View ArticlePubMedGoogle Scholar
- Herrero C, Hu X, Li WP, Samuels S, Sharif MN, Kotenko S, Ivashkiv LB: Reprogramming of IL-10 activity and signaling by IFN-gamma. J Immunol 2003,171(10):5034-5041.View ArticlePubMedGoogle Scholar
- Chen L, Wu J: Systems biology for complex diseases. J Mol Cell Biol 2012,4(3):125-126. 10.1093/jmcb/mjs022View ArticlePubMedGoogle Scholar
- Yamada S, Shiono S, Joo A, Yoshimura A: Control mechanism of JAK/STAT signal transduction pathway. FEBS Lett 2003,534(1-3):190-196.View ArticlePubMedGoogle Scholar
- Zi Z, Cho KH, Sung MH, Xia X, Zheng J, Sun Z: In silico identification of the key components and steps in IFN-gamma induced JAK-STAT signaling pathway. FEBS Lett 2005,579(5):1101-1108. 10.1016/j.febslet.2005.01.009View ArticlePubMedGoogle Scholar
- Singh A, Jayaraman A, Hahn J: Modeling regulatory mechanisms in IL-6 signal transduction in hepatocytes. Biotechnol Bioeng 2006,95(5):850-862. 10.1002/bit.21026View ArticlePubMedGoogle Scholar
- Moya C, Huang Z, Cheng P, Jayaraman A, Hahn J: Investigation of IL-6 and IL-10 signalling via mathematical modelling. IET Syst Biol 2011,5(1):15. 10.1049/iet-syb.2009.0060View ArticlePubMedGoogle Scholar
- Chen X, Vinkemeier U, Zhao Y, Jeruzalmi D, Darnell JE Jr, Kuriyan J: Crystal structure of a tyrosine phosphorylated STAT-1 dimer bound to DNA. Cell 1998,93(5):827-839. 10.1016/S0092-8674(00)81443-9View ArticlePubMedGoogle Scholar
- Becker S, Groner B, Muller CW: Three-dimensional structure of the Stat3beta homodimer bound to DNA. Nature 1998,394(6689):145-151. 10.1038/28101View ArticlePubMedGoogle Scholar
- Greenlund AC, Farrar MA, Viviano BL, Schreiber RD: Ligand-induced IFN gamma receptor tyrosine phosphorylation couples the receptor to its signal transduction system (p91). EMBO J 1994,13(7):1591-1600.PubMed CentralPubMedGoogle Scholar
- Heim MH, Kerr IM, Stark GR, Darnell JE Jr: Contribution of STAT SH2 groups to specific interferon signaling by the Jak-STAT pathway. Science 1995,267(5202):1347-1349. 10.1126/science.7871432View ArticlePubMedGoogle Scholar
- Hemmann U, Gerhartz C, Heesel B, Sasse J, Kurapkat G, Grotzinger J, Wollmer A, Zhong Z, Darnell JE Jr, Graeve L: Differential activation of acute phase response factor/Stat3 and Stat1 via the cytoplasmic domain of the interleukin 6 signal transducer gp130. II. Src homology SH2 domains define the specificity of stat factor activation. J Biol Chem 1996,271(22):12999-13007. 10.1074/jbc.271.22.12999View ArticlePubMedGoogle Scholar
- Shuai K, Ziemiecki A, Wilks AF, Harpur AG, Sadowski HB, Gilman MZ, Darnell JE: Polypeptide signalling to the nucleus through tyrosine phosphorylation of Jak and Stat proteins. Nature 1993,366(6455):580-583. 10.1038/366580a0View ArticlePubMedGoogle Scholar
- Kaptein A, Paillard V, Saunders M: Dominant negative stat3 mutant inhibits interleukin-6-induced Jak-STAT signal transduction. J Biol Chem 1996,271(11):5961-5964. 10.1074/jbc.271.11.5961View ArticlePubMedGoogle Scholar
- Wiederkehr-Adam M, Ernst P, Muller K, Bieck E, Gombert FO, Ottl J, Graff P, Grossmuller F, Heim MH: Characterization of phosphopeptide motifs specific for the Src homology 2 domains of signal transducer and activator of transcription 1 (STAT1) and STAT3. J Biol Chem 2003,278(18):16117-16128. 10.1074/jbc.M300261200View ArticlePubMedGoogle Scholar
- Pfeifer AC, Timmer J, Klingmuller U: Systems biology of JAK/STAT signalling. Essays Biochem 2008, 45: 109-120. 10.1042/BSE0450109View ArticlePubMedGoogle Scholar
- Sasaki A, Yasukawa H, Suzuki A, Kamizono S, Syoda T, Kinjyo I, Sasaki M, Johnston JA, Yoshimura A: Cytokine-inducible SH2 protein-3 (CIS3/SOCS3) inhibits Janus tyrosine kinase by binding through the N-terminal kinase inhibitory region as well as SH2 domain. Genes Cells 1999,4(6):339-351. 10.1046/j.1365-2443.1999.00263.xView ArticlePubMedGoogle Scholar
- Lehmann U, Schmitz J, Weissenbach M, Sobota RM, Hortner M, Friederichs K, Behrmann I, Tsiaris W, Sasaki A, Schneider-Mergener J: SHP2 and SOCS3 contribute to Tyr-759-dependent attenuation of interleukin-6 signaling through gp130. J Biol Chem 2003,278(1):661-671.View ArticlePubMedGoogle Scholar
- De Souza D, Fabri LJ, Nash A, Hilton DJ, Nicola NA, Baca M: SH2 domains from suppressor of cytokine signaling-3 and protein tyrosine phosphatase SHP-2 have similar binding specificities. Biochemistry 2002,41(29):9229-9236. 10.1021/bi0259507View ArticlePubMedGoogle Scholar
- Chen XP, Losman JA, Rothman P: SOCS proteins, regulators of intracellular signaling. Immunity 2000,13(3):287-290. 10.1016/S1074-7613(00)00028-5View ArticlePubMedGoogle Scholar
- Alexander WS, Starr R, Fenner JE, Scott CL, Handman E, Sprigg NS, Corbin JE, Cornish AL, Darwiche R, Owczarek CM: SOCS1 is a critical inhibitor of interferon gamma signaling and prevents the potentially fatal neonatal actions of this cytokine. Cell 1999,98(5):597-608. 10.1016/S0092-8674(00)80047-1View ArticlePubMedGoogle Scholar
- Croker BA, Krebs DL, Zhang JG, Wormald S, Willson TA, Stanley EG, Robb L, Greenhalgh CJ, Forster I, Clausen BE: SOCS3 negatively regulates IL-6 signaling in vivo. Nat Immunol 2003,4(6):540-545. 10.1038/ni931View ArticlePubMedGoogle Scholar
- He B, You L, Uematsu K, Matsangou M, Xu Z, He M, McCormick F, Jablons DM: Cloning and characterization of a functional promoter of the human SOCS-3 gene. Biochem Biophys Res Commun 2003,301(2):386-391. 10.1016/S0006-291X(02)03071-1View ArticlePubMedGoogle Scholar
- Haan S, Keller JF, Behrmann I, Heinrich PC, Haan C: Multiple reasons for an inefficient STAT1 response upon IL-6-type cytokine stimulation. Cell Signal 2005,17(12):1542-1550. 10.1016/j.cellsig.2005.03.010View ArticlePubMedGoogle Scholar
- Milocco LH, Haslam JA, Rosen J, Seidel HM: Design of conditionally active STATs: insights into STAT activation and gene regulatory function. Mol Cell Biol 1999,19(4):2913-2920.PubMed CentralView ArticlePubMedGoogle Scholar
- Meyer T, Begitt A, Lodige I, van Rossum M, Vinkemeier U: Constitutive and IFN-gamma-induced nuclear import of STAT1 proceed through independent pathways. EMBO J 2002,21(3):344-354. 10.1093/emboj/21.3.344PubMed CentralView ArticlePubMedGoogle Scholar
- Haspel RL, Darnell JE Jr: A nuclear protein tyrosine phosphatase is required for the inactivation of Stat1. Proc Natl Acad Sci USA 1999,96(18):10188-10193. 10.1073/pnas.96.18.10188PubMed CentralView ArticlePubMedGoogle Scholar
- McBride KM, McDonald C, Reich NC: Nuclear export signal located within theDNA-binding domain of the STAT1transcription factor. EMBO J 2000,19(22):6196-6206. 10.1093/emboj/19.22.6196PubMed CentralView ArticlePubMedGoogle Scholar
- Pestka S, Kotenko SV, Muthukumaran G, Izotova LS, Cook JR, Garotta G: The interferon gamma (IFN-gamma) receptor: a paradigm for the multichain cytokine receptor. Cytokine Growth Factor Rev 1997,8(3):189-206. 10.1016/S1359-6101(97)00009-9View ArticlePubMedGoogle Scholar
- Stahl N, Farruggella TJ, Boulton TG, Zhong Z, Darnell JE Jr, Yancopoulos GD: Choice of STATs and other substrates specified by modular tyrosine-based motifs in cytokine receptors. Science 1995,267(5202):1349-1353. 10.1126/science.7871433View ArticlePubMedGoogle Scholar
- Guschin D, Rogers N, Briscoe J, Witthuhn B, Watling D, Horn F, Pellegrini S, Yasukawa K, Heinrich P, Stark GR: A major role for the protein tyrosine kinase JAK1 in the JAK/STAT signal transduction pathway in response to interleukin-6. EMBO J 1995,14(7):1421-1429.PubMed CentralPubMedGoogle Scholar
- Niwa Y, Kanda H, Shikauchi Y, Saiura A, Matsubara K, Kitagawa T, Yamamoto J, Kubo T, Yoshikawa H: Methylation silencing of SOCS-3 promotes cell growth and migration by enhancing JAK/STAT and FAK signalings in human hepatocellular carcinoma. Oncogene 2005,24(42):6406-6417.PubMedGoogle Scholar
- Dimberg LY, Dimberg A, Ivarsson K, Fryknas M, Rickardson L, Tobin G, Ekman S, Larsson R, Gullberg U, Nilsson K: Stat1 activation attenuates IL-6 induced Stat3 activity but does not alter apoptosis sensitivity in multiple myeloma. BMC Cancer 2012, 12: 318. 10.1186/1471-2407-12-318PubMed CentralView ArticlePubMedGoogle Scholar
- Caldenhoven E, Buitenhuis M, van Dijk TB, Raaijmakers JA, Lammers JW, Koenderman L, de Groot RP: Lineage-specific activation of STAT3 by interferon-gamma in human neutrophils. J Leukoc Biol 1999,65(3):391-396.PubMedGoogle Scholar
- Fletcher S, Page BD, Zhang X, Yue P, Li ZH, Sharmeen S, Singh J, Zhao W, Schimmer AD, Trudel S: Antagonism of the Stat3-Stat3 protein dimer with salicylic acid based small molecules. ChemMedChem 2011,6(8):1459-1470. 10.1002/cmdc.201100194PubMed CentralView ArticlePubMedGoogle Scholar
- Souissi I, Ladam P, Cognet JA, Le Coquil S, Varin-Blank N, Baran-Marszak F, Metelev V, Fagard R: A STAT3-inhibitory hairpin decoy oligodeoxynucleotide discriminates between STAT1 and STAT3 and induces death in a human colon carcinoma cell line. Mol Cancer 2012, 11: 12. 10.1186/1476-4598-11-12PubMed CentralView ArticlePubMedGoogle Scholar
- Hoops S, Sahle S, Gauges R, Lee C, Pahle J, Simus N, Singhal M, Xu L, Mendes P, Kummer U: COPASI-a COmplex PAthway SImulator. Bioinformatics 2006,22(24):3067-3074. 10.1093/bioinformatics/btl485View ArticlePubMedGoogle Scholar
- Funahashi A, Jouraku A, Matsuoka Y, Kitano H: Integration of cell designer and SABIO-RK. In Silico Biol 2007,7(2 Suppl):S81-90.PubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.