Institute of Fundamental Technological Research

Cells respond dynamically to pulsatile cytokine stimulation. Here we report that single, or well-spaced pulses of TNFα (>100 min apart) give a high probability of NF-κB activation. However, fewer cells respond to shorter pulse intervals (<100 min) suggesting a heterogeneous refractory state. This refractory state is established in the signal transduction network downstream of TNFR and upstream of IKK, and depends on the level of the NF-κB system negative feedback protein A20. If a second pulse within the refractory phase is IL-1β instead of TNFα, all of the cells respond. This suggests a mechanism by which two cytokines can synergistically activate an inflammatory response. Gene expression analyses show strong correlation between the cellular dynamic response and NF-κB-dependent target gene activation. These data suggest that refractory states in the NF-κB system constitute an inherent design motif of the inflammatory response and we suggest that this may avoid harmful homogenous cellular activation.

Dynamic cellular systems reprogram gene expression to ensure appropriate cellular fate responses to specific extracellular cues. Here we demonstrate that the dynamics of Nuclear Factor kappa B (NF-κB) signalling and the cell cycle are prioritised differently depending on the timing of an inflammatory signal. Using iterative experimental and computational analyses, we show physical and functional interactions between NF-κB and the E2 Factor 1 (E2F-1) and E2 Factor 4 (E2F-4) cell cycle regulators. These interactions modulate the NF-κB response. In S-phase, the NF-κB response was delayed or repressed, while cell cycle progression was unimpeded. By contrast, activation of NF-κB at the G1/S boundary resulted in a longer cell cycle and more synchronous initial NF-κB responses between cells. These data identify new mechanisms by which the cellular response to stress is differentially controlled at different stages of the cell cycle.

A central challenge in computational modeling of dynamic biological systems is parameter inference from experimental time course measurements. However, one would not only like to infer kinetic parameters but also study their variability from cell to cell. Here we focus on the case where single-cell fluorescent protein imaging time series data are available for a population of cells. Based on van Kampen’s linear noise approximation, we derive a dynamic state space model for molecular populations which is then extended to a hierarchical model. This model has potential to address the sources of variability relevant to single-cell data, namely, intrinsic noise due to the stochastic nature of the birth and death processes involved in reactions and extrinsic noise arising from the cell-to-cell variation of kinetic parameters. In order to infer such a model from experimental data, one must also quantify the measurement process where one has to allow for nonmeasurable molecular species as well as measurement noise of unknown level and variance. The availability of multiple single-cell time series data here provides a unique testbed to fit such a model and quantify these different sources of variation from experimental data.

Linear noise approximation, kinetic parameter estimation, intrinsic and extrinsic noise, state space model and Kalman filter, Bayesian hierarchical modeling

Heterogeneity between individual cells is a common feature of dynamic cellular processes, including signaling, transcription, and cell fate; yet the overall tissue level physiological phenotype needs to be carefully controlled to avoid fluctuations. Here we show that in the NF-κB signaling system, the precise timing of a dual-delayed negative feedback motif [involving stochastic transcription of inhibitor κB (IκB)-α and -ε] is optimized to induce heterogeneous timing of NF-κB oscillations between individual cells. We suggest that this dual-delayed negative feedback motif enables NF-κB signaling to generate robust single cell oscillations by reducing sensitivity to key parameter perturbations. Simultaneously, enhanced cell heterogeneity may represent a mechanism that controls the overall coordination and stability of cell population responses by decreasing temporal fluctuations of paracrine signaling. It has often been thought that dynamic biological systems may have evolved to maximize robustness through cell-to-cell coordination and homogeneity. Our analyses suggest in contrast, that this cellular variation might be advantageous and subject to evolutionary selection. Alternative types of therapy could perhaps be designed to modulate this cellular heterogeneity.

Gene expression in living cells is highly dynamic, but temporal patterns of gene expression in intact tissues are largely unknown. The mammalian pituitary gland comprises several intermingled cell types, organised as interdigitated networks that interact functionally to generate co-ordinated hormone secretion. Live-cell imaging was used to quantify patterns of reporter gene expression in dispersed lactotrophic cells or intact pituitary tissue from bacterial artificial chromosome (BAC) transgenic rats in which a large prolactin genomic fragment directed expression of luciferase or destabilised enhanced green fluorescent protein (d2EGFP). Prolactin promoter activity in transgenic pituitaries varied with time across different regions of the gland. Although amplitude of transcriptional responses differed, all regions of the gland displayed similar overall patterns of reporter gene expression over a 50-hour period, implying overall co-ordination of cellular behaviour. By contrast, enzymatically dispersed pituitary cell cultures showed unsynchronised fluctuations of promoter activity amongst different cells, suggesting that transcriptional patterns were constrained by tissue architecture. Short-term, high resolution, single cell analyses in prolactin-d2EGFP transgenic pituitary slice preparations showed varying transcriptional patterns with little correlation between adjacent cells. Together, these data suggest that pituitary tissue comprises a series of cell ensembles, which individually display a variety of patterns of short-term stochastic behaviour, but together yield long-range and long-term coordinated behaviour.

Live-cell, Microscopy, Pituitary, Prolactin, Transcription

Nuclear factor kappa B (NF-kappaB) signalling is activated by cellular stress and inflammation and regulates cytokine expression. We applied single-cell imaging to investigate dynamic responses to different doses of tumour necrosis factor alpha (TNFalpha). Lower doses activated fewer cells and those responding showed an increasingly variable delay in the initial NF-kappaB nuclear translocation and associated IkappaBalpha degradation. Robust 100 minute nuclear:cytoplasmic NF-kappaB oscillations were observed over a wide range of TNFalpha concentrations. The result is supported by computational analyses, which identified a limit cycle in the system with a stable 100 minute period over a range of stimuli, and indicated no co-operativity in the pathway activation. These results suggest that a stochastic threshold controls functional all-or-nothing responses in individual cells. Deterministic and stochastic models simulated the experimentally observed activation threshold and gave rise to new predictions about the structure of the system and open the way for better mechanistic understanding of physiological TNFalpha activation of inflammatory responses in cells and tissues.

NF-

The nuclear factor κB (NF-κB) transcription factor regulates cellular stress responses and the immune response to infection. NF-κB activation results in oscillations in nuclear NF-κB abundance. To define the function of these oscillations, we treated cells with repeated short pulses of tumor necrosis factor–α at various intervals to mimic pulsatile inflammatory signals. At all pulse intervals that were analyzed, we observed synchronous cycles of NF-κB nuclear translocation. Lower frequency stimulations gave repeated full-amplitude translocations, whereas higher frequency pulses gave reduced translocation, indicating a failure to reset. Deterministic and stochastic mathematical models predicted how negative feedback loops regulate both the resetting of the system and cellular heterogeneity. Altering the stimulation intervals gave different patterns of NF-κB–dependent gene expression, which supports the idea that oscillation frequency has a functional role.