Cells operate in ever changing 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 types.1 Cell-to-cell communication is primarily mediated by signaling molecules that form spatiotemporal concentration gradients, requiring cells to respond to a wide range of signal intensities. Here we use high-throughput microfluidic cell culture, quantitative gene expression analysis and mathematical modeling to investigate how single mammalian cells respond to different concentrations of the signaling molecule TNF-α, and relay information to the gene expression programs via the transcription 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 studies, 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 analog parameters to modulate the outcome; these parameters include NF-κB peak intensity, response time and number of oscillations. We developed a mathematical model that reproduces both the digital and analog dynamics as well as the most gene expression profiles at all measured conditions, constituting a broadly applicable model for TNF-α induced NF-κB signaling in various types of cells. These results highlight the value of high-throughput quantitative measurements at the single-cell level in understanding how biological systems operate.

No abstract

No abstract

The stochastic dynamics of T cell receptor (TCR) signaling are studied using a mathematical model intended to capture kinetic proofreading (sensitivity to ligand-receptor binding kinetics) and negative and positive feedback regulation mediated respectively by the phosphatase SHP1 and the MAP kinase ERK. The model incorporates protein-protein interactions involved in initiating TCR-mediated cellular responses and reproduces several experimental observations about the behavior of TCR signaling, including robust responses to as few as a handful of ligands (agonist peptide-MHC complexes on an antigen-presenting cell), distinct responses to ligands that bind TCR with different lifetimes, and antagonism. Analysis of the model indicates that TCR signaling dynamics are marked by significant stochastic fluctuations and bistability, which is caused by the competition between the positive and negative feedbacks. Stochastic fluctuations are such that single-cell trajectories differ qualitatively from the trajectory predicted in the deterministic approximation of the dynamics. Because of bistability, the average of single-cell trajectories differs markedly from the deterministic trajectory. Bistability combined with stochastic fluctuations allows for switch-like responses to signals, which may aid T cells in making committed cell-fate decisions.

No abstract

Living cells may be considered noisy or stochastic biochemical reactors. In eukaryotic cells, in which the number of protein or mRNA molecules is relatively large, the stochastic effects originate primarily in regulation of gene activity. Transcriptional activity of a gene can be initiated by transactivator molecules binding to the specific regulatory site(s) in the target gene. The stochasticity of activator binding and dissociation is amplified by transcription and translation, since target gene activation results in a burst of mRNAs molecules, and each copy of mRNA then serves as a template for numerous protein molecules. In this article, we reformulate our model of the NF-kappaB regulatory module to analyze a single cell regulation. Ordinary differential equations, used for description of fast reaction channels of processes involving a large number of molecules, are combined with a stochastic switch to account for the activity of the genes involved. The stochasticity in gene transcription causes simulated cells to exhibit large variability. Moreover, none of them behaves like an average cell. Although the average mRNA and protein levels remain constant before tumor necrosis factor (TNF) stimulation, and stabilize after a prolonged TNF stimulation, in any single cell these levels oscillate stochastically in the absence of TNF and keep oscillating under the prolonged TNF stimulation. However, in a short period of approximately 90 min, most cells are synchronized by the TNF signal, and exhibit similar kinetics. We hypothesize that this synchronization is crucial for proper activation of early genes controlling inflammation. Our theoretical predictions of single cell kinetics are supported by recent experimental studies of oscillations in NF-kappaB signaling made on single cells.

BackgroundThe NF-κB regulatory network controls innate immune response by transducing variety of pathogen-derived and cytokine stimuli into well defined single-cell gene regulatory events.ResultsWe analyze the network by means of the model combining a deterministic description for molecular species with large cellular concentrations with two classes of stochastic switches: cell-surface receptor activation by TNFα ligand, and IκBα and A20 genes activation by NF-κB molecules. Both stochastic switches are associated with amplification pathways capable of translating single molecular events into tens of thousands of synthesized or degraded proteins. Here, we show that at a low TNFα dose only a fraction of cells are activated, but in these activated cells the amplification mechanisms assure that the amplitude of NF-κB nuclear translocation remains above a threshold. Similarly, the lower nuclear NF-κB concentration only reduces the probability of gene activation, but does not reduce gene expression of those responding.ConclusionThese two effects provide a particular stochastic robustness in cell regulation, allowing cells to respond differently to the same stimuli, but causing their individual responses to be unequivocal. Both effects are likely to be crucial in the early immune response: Diversity in cell responses causes that the tissue defense is harder to overcome by relatively simple programs coded in viruses and other pathogens. The more focused single-cell responses help cells to choose their individual fates such as apoptosis or proliferation. The model supports the hypothesis that binding of single TNFα ligands is sufficient to induce massive NF-κB translocation and activation of NF-κB dependent genes.

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.DOI: http://dx.doi.org/10.7554/eLife.08931.001