Interleukin-6 is a pleiotropic cytokine with high clinical relevance and an important mediator of cellular communication, orchestrating both pro- and anti-inflammatory processes. Interleukin-6-induced signalling is initiated by binding of IL-6 to the IL-6 receptor α and subsequent binding to the signal transducing receptor subunit gp130. This active receptor complex initiates signalling through the Janus kinase/signal transducer and activator of transcription pathway. Of note, IL-6 receptor α exists in a soluble and a transmembrane form. Binding of IL-6 to membrane-bound IL-6 receptor α induces anti-inflammatory classic signalling, whereas binding of IL-6 to soluble IL-6 receptor α induces pro-inflammatory trans-signalling. Trans-signalling has been described to be markedly stronger than classic signalling. Understanding the molecular mechanisms that drive differences between trans- and classic signalling is important for the design of trans-signalling-specific therapies. These differences will be addressed here using a combination of dynamic mathematical modelling and molecular biology. We apply an iterative systems biology approach using set-based modelling and validation approaches combined with quantitative biochemical and cell biological analyses. The combination of experimental analyses and dynamic modelling allows to relate the observed differences between IL-6-induced trans- and classic signalling to cell-type specific differences in the expression and ratios of the individual subunits of the IL-6 receptor complex. Canonical intracellular Jak/STAT signalling is indifferent in IL-6-induced trans- and classic signalling. This study contributes to the understanding of molecular mechanisms of IL-6 signal transduction and underlines the power of combined dynamical modelling, model-based validation and biological experiments. The opposing pro- and anti-inflammatory responses initiated by IL-6 trans- and classic signalling depend solely on the expression ratios of the subunits of the entire receptor complex. By pointing out the importance of the receptor expression ratio for the strength of IL-6 signalling this study lays a foundation for future precision medicine approaches that aim to selectively block pro-inflammatory trans-signalling. Furthermore, the derived models can be used for future therapy design.

Cellular communication via intracellular signalling pathways is crucial. Expression and activation of signalling proteins is heterogenous between isogenic cells of the same cell-type. However, mechanisms evolved to enable sufficient communication and to ensure cellular functions. We use information theory to clarify mechanisms facilitating IL-6-induced JAK/STAT signalling despite cell-to-cell variability. We show that different mechanisms enabling robustness against variability complement each other. Early STAT3 activation is robust as long as cytokine concentrations are low. Robustness at high cytokine concentrations is ensured by high STAT3 expression or serine phosphorylation. Later the feedback-inhibitor SOCS3 increases robustness. Channel Capacity of JAK/STAT signalling is limited by cell-to-cell variability in STAT3 expression and is affected by the same mechanisms governing robustness. Increasing STAT3 amount increases Channel Capacity and robustness, whereas increasing STAT3 tyrosine phosphorylation reduces robustness but increases Channel Capacity. In summary, we elucidate mechanisms preventing dysregulated signalling by enabling reliable JAK/STAT signalling despite cell-to-cell heterogeneity.

Background: Cell-to-cell heterogeneity is an inherent feature of multicellular organisms and is central in all physiological and pathophysiological processes including cellular signal transduction. The cytokine IL-6 is an essential mediator of pro- and anti-inflammatory processes. Dysregulated IL-6-induced intracellular JAK/STAT signalling is associated with severe inflammatory and proliferative diseases. Under physiological conditions JAK/STAT signalling is rigorously controlled and timely orchestrated by regulatory mechanisms such as expression of the feedback-inhibitor SOCS3 and activation of the protein-tyrosine phosphatase SHP2 (PTPN11). Interestingly, the function of negative regulators seems not to be restricted to controlling the strength and timely orchestration of IL-6-induced STAT3 activation. Exemplarily, SOCS3 increases robustness of late IL-6-induced STAT3 activation against heterogenous STAT3 expression and reduces the amount of information transferred through JAK/STAT signalling. Methods: Here we use multiplexed single-cell analyses and information theoretic approaches to clarify whether also SHP2 contributes to robustness of STAT3 activation and whether SHP2 affects the amount of information transferred through IL-6-induced JAK/STAT signalling. Results: SHP2 increases robustness of both basal, cytokine-independent STAT3 activation and early IL-6-induced STAT3 activation against differential STAT3 expression. However, SHP2 does not affect robustness of late IL-6-induced STAT3 activation. In contrast to SOCS3, SHP2 increases the amount of information transferred through IL-6-induced JAK/STAT signalling, probably by reducing cytokine-independent STAT3 activation and thereby increasing sensitivity of the cells. These effects are independent of SHP2-dependent MAPK activation. Conclusion: In summary, the results of this study extend our knowledge of the functions of SHP2 in IL-6-induced JAK/STAT signalling. SHP2 is not only a repressor of basal and cytokine-induced STAT3 activity, but also ensures robustness and transmission of information. Plain English summary Cells within a multicellular organism communicate with each other to exchange information about the environment. Communication between cells is facilitated by soluble molecules that transmit information from one cell to the other. Cytokines such as interleukin-6 are important soluble mediators that are secreted when an organism is faced with infections or inflammation. Secreted cytokines bind to receptors within the membrane of their target cells. This binding induces activation of an intracellular cascade of reactions called signal transduction, which leads to cellular responses. An important example of intracellular signal transduction is JAK/STAT signalling. In healthy organisms signalling is controlled and timed by regulatory mechanisms, whose activation results in a controlled shutdown of signalling pathways. Interestingly, not all cells within an organism are identical. They differ in the amount of proteins involved in signal transduction, such as STAT3. These differences shape cellular communication and responses to intracellular signalling. Here, we show that an important negative regulatory protein called SHP2 (or PTPN11) is not only responsible for shutting down signalling, but also for steering signalling in heterogeneous cell populations. SHP2 increases robustness of STAT3 activation against variable STAT3 amounts in individual cells. Additionally, it increases the amount of information transferred through JAK/STAT signalling by increasing the dynamic range of pathway activation in heterogeneous cell populations. This is an amazing new function of negative regulatory proteins that contributes to communication in heterogeneous multicellular organisms in health and disease.

Background: Modeling of biological pathways is a key issue in systems biology. When constructing a model, it is tempting to incorporate all known interactions of pathway species, which results in models with a large number of unknown parameters. Fortunately, unknown parameters need not necessarily be measured directly, but some parameter values can be estimated indirectly by fitting the model to experimental data. However, parameter fitting, or, more precisely, maximum likelihood parameter estimation, only provides valid results, if the complexity of the model is in balance with the amount and quality of the experimental data. If this is the case the model is said to be identifiable for the given data. If a model turns out to be unidentifiable, two steps can be taken. Either additional experiments need to be conducted, or the model has to be simplified. Results: We propose a systematic procedure for model simplification, which consists of the following steps: estimate the parameters of the model, create an identifiability ranking for the estimated parameters, and simplify the model based on the identifiability analysis results. These steps need to be applied iteratively until the resulting model is identifiable, or equivalently, until parameter variances are small. We choose parameter variances as stopping criterion, since they are concise and easy to interpret. For both, the parameter estimation and the calculation of parameter variances, multi-start parameter estimations are run on a parallel cluster. In contrast to related work in systems biology, we do not suggest simplifying a model by fixing some of its parameters, but change the structure of the model. Conclusions: We apply the proposed approach to a model of early signaling events in the JAK-STAT pathway. The resulting model is not only identifiable with small parameter variances, but also shows the best trade-off between goodness of fit and model complexity.