Background: Ability to adapt to temperature changes trough the Heat Shock Response (HSR) pathways is one of the most fundamental and clinically relevant cellular response systems. Heat Shock (HS) affects the signalling and gene expression responses of the Nuclear Factor κB (NF-κB) transcription factor, a critical regulator of proliferation and inflammation, however, our quantitative understanding of how cells sense and adapt to temperature changes is limited.Methods: We used live-cell time-lapse microscopy and mathematical modelling to understand the signalling of the NF-κB system in the human MCF7 breast adenocarcinoma cells in response to pro-inflammatory Interleukin 1β (IL1β) and Tumour Necrosis Factor α (TNFα) cytokines, following exposure to a 37–43 °C range of physiological and clinical temperatures.Results: We show that exposure to 43 °C 1 h HS inhibits the immediate NF-κB signalling response to TNFα and IL1β stimulation although uptake of cytokines is not impaired. Within 4 h after HS treatment IL1β-induced NF-κB responses return to normal levels, but the recovery of the TNFα-induced responses is still affected. Using siRNA knock-down of Heat Shock Factor 1 (HSF1) we show that this stimulus-specificity is conferred via the Inhibitory κB kinase (IKK) signalosome where HSF1-dependent feedback regulates TNFα, but not IL1β-mediated IKK recovery post HS. Furthermore, we demonstrate that through the temperature-dependent denaturation and recovery of IKK, TNFα and IL1β-mediated signalling exhibit different temperature sensitivity and adaptation to repeated HS when exposed to a 37–43 °C temperature range. Specifically, IL1β-mediated NF-κB responses are more robust to temperature changes in comparison to those induced by TNFα treatment.Conclusions: We demonstrate that the kinetics of the NF-κB system following temperature stress is cytokine specific and exhibit differential adaptation to temperature changes. We propose that this differential temperature sensitivity is mediated via the IKK signalosome, which acts as a bona fide temperature sensor trough the HSR cross-talk. This novel quantitative understanding of NF-κB and HSR interactions is fundamentally important for the potential optimization of therapeutic hyperthermia protocols.

Elevated temperature induces the heat shock (HS) response, which modulates cell proliferation, apoptosis, the immune and inflammatory responses. However, specific mechanisms linking the HS response pathways to major cellular signaling systems are not fully understood. Here we used integrated computational and experimental approaches to quantitatively analyze the crosstalk mechanisms between the HS-response and a master regulator of inflammation, cell proliferation, and apoptosis the Nuclear Factor κB (NF-κB) system. We found that populations of human osteosarcoma cells, exposed to a clinically relevant 43°C HS had an attenuated NF-κB p65 response to Tumor Necrosis Factor α (TNFα) treatment. The degree of inhibition of the NF-κB response depended on the HS exposure time. Mathematical modeling of single cells indicated that individual crosstalk mechanisms differentially encode HS-mediated NF-κB responses while being consistent with the observed population-level responses. In particular “all-or-nothing” encoding mechanisms were involved in the HS-dependent regulation of the IKK activity and IκBα phosphorylation, while others involving transport were “analogue”. In order to discriminate between these mechanisms, we used live-cell imaging of nuclear translocations of the NF-κB p65 subunit. The single cell responses exhibited “all-or-nothing” encoding. While most cells did not respond to TNFα stimulation after a 60 min HS, 27% showed responses similar to those not receiving HS. We further demonstrated experimentally and theoretically that the predicted inhibition of IKK activity was consistent with the observed HS-dependent depletion of the IKKα and IKKβ subunits in whole cell lysates. However, a combination of “all-or-nothing” crosstalk mechanisms was required to completely recapitulate the single cell data. We postulate therefore that the heterogeneity of the single cell responses might be explained by the cell-intrinsic variability of HS-modulated IKK signaling. In summary, we show that high temperature modulates NF-κB responses in single cells in a complex and unintuitive manner, which needs to be considered in hyperthermia-based treatment strategies. Hyperthermia has been considered a promising strategy to sensitize cancer cells to treatment. As such, it might potentially increase treatment efficacy while reducing negative side effects. So far, this potential has not been fully realized. One of the major obstacles is a lack of quantitative understanding of crosstalk mechanisms involved in the heat shock and other cellular stress-factor systems. In this work, we used experimental population-level and live-cell imaging approaches to characterize the NF-κB signaling system response to physiologically high temperature. We then iteratively applied mathematical modeling of single cells to infer the underlying regulatory mechanisms. We showed experimentally that elevated temperature attenuated cytokine-induced NF-κB system responses, resulting in complex and unintuitive single cell behavior. Our analyses suggested that these responses are consistent with stochastic cell-intrinsic heat-shock-mediated regulation of the NF-κB system via the IKK signaling module. Considering the oncogenic associations of the NF-κB pathway, the crosstalk mechanisms described in this work are important for the understanding of chemo- and radio-resistance to treatment. However, more but similar, quantitative studies are required to ultimately change currently used therapy protocols.

The paper is focused on application of sensitivity methods to analysis of signaling pathway models. Two basic methods are compared: local, based on standard sensitivity functions, and global, based on Sobol indices. Firstly, a general outline of modeling of signaling pathways by means of ordinary differential equations is briefly described. Afterwards, the methods of sensitivity analysis, known from literature, are introduced and illustrated with a simple example of a dynamical system of the second order. Subsequently, the analysis of the p53/Mdm2 regulatory module, which is a key element of any pathway involving p53 protein, is presented. The results of this analysis suggest that no single method should be chosen for investigation of any signaling pathway model but several of them should be applied to answer important questions about sources of heterogeneity in cells behavior, robustness of signaling pathways and possible molecular drug targets.

Mathematical modeling of signaling pathways and regulatory networks has been supporting experimental research for some time now. Sensitivity analysis, aimed at finding model parameters whose changes yield significantly altered cellular responses, is an important part of modeling work. However, sensitivity methods are often directly transplanted from analysis of technical systems, and thus, they may not serve the purposes of analysis of biological systems. This paper presents a novel sensitivity analysis method that is particularly suited to the task of searching for potential molecular drug targets in signaling pathways. Using two sample models of pathways, p53/Mdm2 regulatory module and IFN-$\beta$-induced JAK/STAT signaling pathway, we show that the method leads to biologically relevant conclusions, identifying processes suitable for targeted pharmacological inhibition, represented by the reduction of kinetic parameter values. That, in turn, facilitates subsequent search for active drug components.

The paper is focused on sensitivity analysis of large-scale models of biological systems that describe dynamics of the so called signaling pathways. These systems are continuous in time but their models are based on discrete-time measurements. Therefore, if sensitivity analysis is used as a tool supporting model development and evaluation of its quality, it should take this fact into account. Such models are usually very complex and include many parameters difficult to estimate in an experimental way. Changes of many of those parameters have little effect on model dynamics, and therefore they are called sloppy. In contrast, other parameters, when changed, lead to substantial changes in model responses and these are called stiff parameters. While this is a well-known fact, and there are methods to discern sloppy parameters from the stiff ones, they have not been utilized, so far, to create parameter rankings and quantify the influence of single parameter changes on system time responses. These single parameter changes are particularly important in analysis of signalling pathways, because they may pinpoint parameters, associated with the processes to be targeted at the molecular level in laboratory experiments. In the paper we present a new, original method of creating parameter rankings, based on an Hessian of a cost function which describes the fit of the model to a discrete experimental data. Its application is explained with simple dynamical systems, representing two typical dynamics exhibited by the signaling pathways.

This paper is concerned with analysis of two anticancer therapy models focused on sensitivity of therapy outcome with respect to model structure and parameters. Realistic periodic therapies are considered, combining cytotoxic and antiangiogenic agents, defined on a fixed time horizon. Tumor size at the end of therapy and average tumor size calculated over therapy horizon are chosen to represent therapy outcome. Sensitivity analysis has been performed numerically, concentrating on model parameters and structure at one hand, and on treatment protocol parameters, on the other. The results show that sensitivity of the therapy outcome highly depends on the model structure, helping to discern a good model. Moreover, it is possible to use this analysis to find a good protocol in case of heterogeneous tumors.