Institute of Fundamental Technological Research

Digital signaling enhances robustness of cellular decisions in noisy environments, but it is unclear how digital systems transmit temporal information about a stimulus. To understand how temporal input information is encoded and decoded by the NF-κB system, we studied transcription factor dynamics and gene regulation under dose- and duration-modulated inflammatory inputs. Mathematical modeling predicted and microfluidic single-cell experiments confirmed that integral of the stimulus (or area, concentration × duration) controls the fraction of cells that activate NF-κB in the population. However, stimulus temporal profile determined NF-κB dynamics, cell-to-cell variability, and gene expression phenotype. A sustained, weak stimulation lead to heterogeneous activation and delayed timing that is transmitted to gene expression. In contrast, a transient, strong stimulus with the same area caused rapid and uniform dynamics. These results show that digital NF-κB signaling enables multidimensional control of cellular phenotype via input profile, allowing parallel and independent control of single-cell activation probability and population heterogeneity.

NF-κB is a key transcription factor that regulates innate immune response. Its activity is tightly controlled by numerous feedback loops, including two negative loops mediated by NF-κB inducible inhibitors, IκBα and A20, which assure oscillatory responses, and by positive feedback loops arising due to the paracrine and autocrine regulation via TNFα, IL-1 and other cytokines. We study the NF-κB system of interlinked negative and positive feedback loops, combining bifurcation analysis of the deterministic approximation with stochastic numerical modeling. Positive feedback assures the existence of limit cycle oscillations in unstimulated wild-type cells and introduces bistability in A20-deficient cells. We demonstrated that cells of significant autocrine potential, i.e., cells characterized by high secretion of TNFα and its receptor TNFR1, may exhibit sustained cytoplasmic–nuclear NF-κB oscillations which start spontaneously due to stochastic fluctuations. In A20-deficient cells even a small TNFα expression rate qualitatively influences system kinetics, leading to long-lasting NF-κB activation in response to a short-pulsed TNFα stimulation. As a consequence, cells with impaired A20 expression or increased TNFα secretion rate are expected to have elevated NF-κB activity even in the absence of stimulation. This may lead to chronic inflammation and promote cancer due to the persistent activation of antiapoptotic genes induced by NF-κB. There is growing evidence that A20 mutations correlate with several types of lymphomas and elevated TNFα secretion is characteristic of many cancers. Interestingly, A20 loss or dysfunction also leaves the organism vulnerable to septic shock and massive apoptosis triggered by the uncontrolled TNFα secretion, which at high levels overcomes the antiapoptotic action of NF-κB. It is thus tempting to speculate that some cancers of deregulated NF-κB signaling may be prone to the pathogen-induced apoptosis.

Cells operate in dynamic environments using extraordinary communication capabilities that emerge from the interactions of genetic circuitry. The mammalian immune response is a striking example of the coordination of different cell types1. Cell-to-cell communication is primarily mediated by signalling molecules that form spatiotemporal concentration gradients, requiring cells to respond to a wide range of signal intensities2. Here we use high-throughput microfluidic cell culture3 and fluorescence microscopy, quantitative gene expression analysis and mathematical modelling to investigate how single mammalian cells respond to different concentrations of the signalling molecule tumour-necrosis factor (TNF)-α, and relay information to the gene expression programs by means of the transcription factor nuclear factor (NF)-κB. We measured NF-κB activity in thousands of live cells under TNF-α doses covering four orders of magnitude. We find, in contrast to population-level studies with bulk assays2, that the activation is heterogeneous and is a digital process at the single-cell level with fewer cells responding at lower doses. Cells also encode a subtle set of analogue parameters to modulate the outcome; these parameters include NF-κB peak intensity, response time and number of oscillations. We developed a stochastic mathematical model that reproduces both the digital and analogue dynamics as well as most gene expression profiles at all measured conditions, constituting a broadly applicable model for TNF-α-induced NF-κB signalling in various types of cells. These results highlight the value of high-throughput quantitative measurements with single-cell resolution in understanding how biological systems operate.

Cell biology, Biophysics, Immunology, Genetics, Genomics