Single-cell NF-kappaB dynamics reveal digital activation and analogue information processing

Quantitative characterization and analysis of the dynamic NF-κB response in microglia

Background

Activation of the NF-κB transcription factor and its associated gene expression in microglia is a key component in the response to brain injury. Its activation is dynamic and is part of a network of biochemical species with multiple feedback regulatory mechanisms. Mathematical modeling, which has been instrumental for understanding the NF-κB response in other cell types, offers a valuable tool to investigate the regulation of NF-κB activation in microglia at a systems level.

Results

We quantify the dynamic response of NF-κB activation and activation of the upstream kinase IKK using ELISA measurements of a microglial cell line following treatment with the pro-inflammatory cytokine TNFα. A new mathematical model is developed based on these data sets using a modular procedure that exploits the feedback structure of the network. We show that the new model requires previously unmodeled dynamics involved in the stimulus-induced degradation of the inhibitor IκBα in order to properly describe microglial NF-κB activation in a statistically consistent manner. This suggests a more prominent role for the ubiquitin-proteasome system in regulating the activation of NF-κB to inflammatory stimuli. We also find that the introduction of nonlinearities in the kinetics of IKK activation and inactivation is essential for proper characterization of transient IKK activity and corresponds to known biological mechanisms. Numerical analyses of the model highlight key regulators of the microglial NF-κB response, as well as those governing IKK activation. Results illustrate the dynamic regulatory mechanisms and the robust yet fragile nature of the negative feedback regulated network.

Conclusions

We have developed a new mathematical model that incorporates previously unmodeled dynamics to characterize the dynamic response of the NF-κB signaling network in microglia. This model is the first of its kind for microglia and provides a tool for the quantitative, systems level study the dynamic cellular response to inflammatory stimuli.

Dynamic chromatin localization of Sirt6 shapes stress- and aging-related transcriptional networks

The sirtuin Sirt6 is a NAD-dependent histone deacetylase that is implicated in gene regulation and lifespan control. Sirt6 can interact with the stress-responsive transcription factor NF-κB and regulate some NF-κB target genes, but the full scope of Sirt6 target genes as well as dynamics of Sirt6 occupancy on chromatin are not known. Here we map Sirt6 occupancy on mouse promoters genome-wide and show that Sirt6 occupancy is highly dynamic in response to TNF-α. More than half of Sirt6 target genes are only revealed upon stress-signaling. The majority of genes bound by NF-κB subunit RelA recruit Sirt6, and dynamic Sirt6 relocalization is largely driven in a RelA-dependent manner. Integrative analysis with global gene expression patterns in wild-type, Sirt6−/−, and double Sirt6−/− RelA−/− cells reveals the epistatic relationships between Sirt6 and RelA in shaping diverse temporal patterns of gene expression. Genes under the direct joint control of Sirt6 and RelA include several with prominent roles in cell senescence and organismal aging. These data suggest dynamic chromatin relocalization of Sirt6 as a key output of NF-κB signaling in stress response and aging.

Sirtuins (Sirt) are a family of enzymes that modify chromatin and other proteins to affect gene activity. Loss of Sirt6 leads to a progeria-like phenotype in mice, but the target of SIRT6 action has been elusive. Here we show that Sirt6 binds to thousands of gene promoters in a stress-inducible fashion, guided by the stress-responsive transcription factor NF-κB. Both the departure and arrival of Sirt6 alter gene expression, shaping the temporal dynamics of NF-κB transcriptional response and directly controlling the expression of other key regulators of aging. These findings provide the first view of how an oscillatory transcription factor can drive a progression of chromatin changes over time.

Systems biology in immunology: a computational modeling perspective

Systems biology is an emerging discipline that combines high-content, multiplexed measurements with informatic and computational modeling methods to better understand biological function at various scales. Here we present a detailed review of the methods used to create computational models and to conduct simulations of immune function. We provide descriptions of the key data-gathering techniques employed to generate the quantitative and qualitative data required for such modeling and simulation and summarize the progress to date in applying these tools and techniques to questions of immunological interest, including infectious disease. We include comments on what insights modeling can provide that complement information obtained from the more familiar experimental discovery methods used by most investigators and the reasons why quantitative methods are needed to eventually produce a better understanding of immune system operation in health and disease.

Follicular helper CD4 T cells (TFH)

T cell help to B cells is a fundamental aspect of adaptive immunity and the generation of immunological memory. Follicular helper CD4 T (T(FH)) cells are the specialized providers of B cell help. T(FH) cells depend on expression of the master regulator transcription factor Bcl6. Distinguishing features of T(FH) cells are the expression of CXCR5, PD-1, SAP (SH2D1A), IL-21, and ICOS, among other molecules, and the absence of Blimp-1 (prdm1). T(FH) cells are important for the formation of germinal centers. Once germinal centers are formed, T(FH) cells are needed to maintain them and to regulate germinal center B cell differentiation into plasma cells and memory B cells. This review covers T(FH) differentiation, T(FH) functions, and human T(FH) cells, discussing recent progress and areas of uncertainty or disagreement in the literature, and it debates the developmental relationship between T(FH) cells and other CD4 T cell subsets (Th1, Th2, Th17, iTreg).

Principles and strategies for developing network models in cancer

The flood of genome-wide data generated by high-throughput technologies currently provides biologists with an unprecedented opportunity: to manipulate, query, and reconstruct functional molecular networks of cells. Here, we outline three underlying principles and six strategies to infer network models from genomic data. Then, using cancer as an example, we describe experimental and computational approaches to infer "differential" networks that can identify genes and processes driving disease phenotypes. In conclusion, we discuss how a network-level understanding of cancer can be used to predict drug response and guide therapeutics.

Impulse control: temporal dynamics in gene transcription

Regulatory circuits controlling gene expression constantly rewire to adapt to environmental stimuli, differentiation cues, and disease. We review our current understanding of the temporal dynamics of gene expression in eukaryotes and prokaryotes and the molecular mechanisms that shape them. We delineate several prototypical temporal patterns, including "impulse" (or single-pulse) patterns in response to transient environmental stimuli, sustained (or state-transitioning) patterns in response to developmental cues, and oscillating patterns. We focus on impulse responses and their higher-order temporal organization in regulons and cascades and describe how core protein circuits and cis-regulatory sequences in promoters integrate with chromatin architecture to generate these responses.

Inflammatory gene regulatory networks in amnion cells following cytokine stimulation: translational systems approach to modeling human parturition

A majority of the studies examining the molecular regulation of human labor have been conducted using single gene approaches. While the technology to produce multi-dimensional datasets is readily available, the means for facile analysis of such data are limited. The objective of this study was to develop a systems approach to infer regulatory mechanisms governing global gene expression in cytokine-challenged cells in vitro, and to apply these methods to predict gene regulatory networks (GRNs) in intrauterine tissues during term parturition. To this end, microarray analysis was applied to human amnion mesenchymal cells (AMCs) stimulated with interleukin-1β, and differentially expressed transcripts were subjected to hierarchical clustering, temporal expression profiling, and motif enrichment analysis, from which a GRN was constructed. These methods were then applied to fetal membrane specimens collected in the absence or presence of spontaneous term labor. Analysis of cytokine-responsive genes in AMCs revealed a sterile immune response signature, with promoters enriched in response elements for several inflammation-associated transcription factors. In comparison to the fetal membrane dataset, there were 34 genes commonly upregulated, many of which were part of an acute inflammation gene expression signature. Binding motifs for nuclear factor-κB were prominent in the gene interaction and regulatory networks for both datasets; however, we found little evidence to support the utilization of pathogen-associated molecular pattern (PAMP) signaling. The tissue specimens were also enriched for transcripts governed by hypoxia-inducible factor. The approach presented here provides an uncomplicated means to infer global relationships among gene clusters involved in cellular responses to labor-associated signals.

System Integration - A Major Step toward Lab on a Chip

Microfluidics holds great promise to revolutionize various areas of biological engineering, such as single cell analysis, environmental monitoring, regenerative medicine, and point-of-care diagnostics. Despite the fact that intensive efforts have been devoted into the field in the past decades, microfluidics has not yet been adopted widely. It is increasingly realized that an effective system integration strategy that is low cost and broadly applicable to various biological engineering situations is required to fully realize the potential of microfluidics. In this article, we review several promising system integration approaches for microfluidics and discuss their advantages, limitations, and applications. Future advancements of these microfluidic strategies will lead toward translational lab-on-a-chip systems for a wide spectrum of biological engineering applications.

Macroscopic law of conservation revealed in the population dynamics of Toll-like receptor signaling

Stimulating the receptors of a single cell generates stochastic intracellular signaling. The fluctuating response has been attributed to the low abundance of signaling molecules and the spatio-temporal effects of diffusion and crowding. At population level, however, cells are able to execute well-defined deterministic biological processes such as growth, division, differentiation and immune response. These data reflect biology as a system possessing microscopic and macroscopic dynamics. This commentary discusses the average population response of the Toll-like receptor (TLR) 3 and 4 signaling. Without requiring detailed experimental data, linear response equations together with the fundamental law of information conservation have been used to decipher novel network features such as unknown intermediates, processes and cross-talk mechanisms. For single cell response, however, such simplicity seems far from reality. Thus, as observed in any other complex systems, biology can be considered to possess order and disorder, inheriting a mixture of predictable population level and unpredictable single cell outcomes.

Rapid prototyping of arrayed microfluidic systems in polystyrene for cell-based assays

Microfluidic cell-based systems have enabled the study of cellular phenomena with improved spatiotemporal control of the microenvironment and at increased throughput. While poly(dimethylsiloxane) (PDMS) has emerged as the most popular material in microfluidics research, it has specific limitations that prevent microfluidic platforms from achieving their full potential. We present here a complete process, ranging from mold design to embossing and bonding, that describes the fabrication of polystyrene (PS) microfluidic devices with similar cost and time expenditures as PDMS-based devices. Emphasis was placed on creating methods that can compete with PDMS fabrication methods in terms of robustness, complexity, and time requirements. To achieve this goal, several improvements were made to remove critical bottlenecks in existing PS embossing methods. First, traditional lithographic techniques were adapted to fabricate bulk epoxy molds capable of resisting high temperatures and pressures. Second, a method was developed to emboss through-holes in a PS layer, enabling creation of large arrays of independent microfluidic systems on a single device without need to manually create access ports. Third, thermal bonding of PS layers was optimized in order to achieve quality bonding over large arrays of microsystems. The choice of materials and methods was validated for biological function in two different cell-based applications to demonstrate the versatility of our streamlined fabrication process.

Multiscale computational modeling reveals a critical role for TNF-α receptor 1 dynamics in tuberculosis granuloma formation

Multiple immune factors control host responses to Mycobacterium tuberculosis infection, including the formation of granulomas, which are aggregates of immune cells whose function may reflect success or failure of the host to contain infection. One such factor is TNF-α. TNF-α has been experimentally characterized to have the following activities in M. tuberculosis infection: macrophage activation, apoptosis, and chemokine and cytokine production. Availability of TNF-α within a granuloma has been proposed to play a critical role in immunity to M. tuberculosis. However, in vivo measurement of a TNF-α concentration gradient and activities within a granuloma are not experimentally feasible. Further, processes that control TNF-α concentration and activities in a granuloma remain unknown. We developed a multiscale computational model that includes molecular, cellular, and tissue scale events that occur during granuloma formation and maintenance in lung. We use our model to identify processes that regulate TNF-α concentration and cellular behaviors and thus influence the outcome of infection within a granuloma. Our model predicts that TNF-αR1 internalization kinetics play a critical role in infection control within a granuloma, controlling whether there is clearance of bacteria, excessive inflammation, containment of bacteria within a stable granuloma, or uncontrolled growth of bacteria. Our results suggest that there is an interplay between TNF-α and bacterial levels in a granuloma that is controlled by the combined effects of both molecular and cellular scale processes. Finally, our model elucidates processes involved in immunity to M. tuberculosis that may be new targets for therapy.

Illuminating intracellular signaling and molecules for single cell analysis

Fluorescent and bioluminescent proteins are now widely used for detection of small molecules and various intracellular events ranging from protein conformational change to cell death in living cells. To analyze the dynamics of molecular processes in real time at the level of single cells, engineered protein-based probes with higher sensitivity and selectivity are required. The probes can be entirely genetically encoded and can comprise fusions of different proteins or domains. This review specifically examines basic concepts of designing genetically encoded fluorescent and bioluminescent probes developed in the past decade, highlighting some potential applications for basic research and for drug discovery.

Origins of regulated cell-to-cell variability

Single-cell measurements and lineage-tracing experiments are revealing that phenotypic cell-to-cell variability is often the result of deterministic processes, despite the existence of intrinsic noise in molecular networks. In most cases, this determinism represents largely uncharacterized molecular regulatory mechanisms, which places the study of cell-to-cell variability in the realm of molecular cell biology. Further research in the field will be important to advance quantitative cell biology because it will provide new insights into the mechanisms by which cells coordinate their intracellular activities in the spatiotemporal context of the multicellular environment.

Light control of plasma membrane recruitment using the Phy-PIF system

The ability to control the activity of intracellular signaling processes in live cells would be an extraordinarily powerful tool. Ideally, such an intracellular input would be (i) genetically encoded, (ii) able to be turned on and off in defined temporal or spatial patterns, (iii) fast to switch between on and off states, and (iv) orthogonal to other cellular processes. The light-gated interaction between fragments of two plant proteins--termed Phy and PIF--satisfies each of these constraints. In this system, Phy can be switched between two conformations using red and infrared light, while PIF only binds one of these states. This chapter describes known constraints for designing genetic constructs using Phy and PIF and provides protocols for expressing these constructs in mammalian cells, purifying the small molecule chromophore required for the system's light responsivity, and measuring light-gated binding by microscopy.

Reporting from the field: genetically encoded fluorescent reporters uncover signaling dynamics in living biological systems

Real-time visualization of a wide range of biochemical processes in living systems is being made possible through the development and application of genetically encoded fluorescent reporters. These versatile biosensors have proven themselves tailor-made to the study of signal transduction, and in this review, we discuss some of the unique insights that they continue to provide regarding the spatial organization and dynamic regulation of intracellular signaling networks. In addition, we explore the more recent push to expand the scope of biological phenomena that can be monitored using these reporters, while also considering the potential to integrate this highly adaptable technology with a number of emerging techniques that may significantly broaden our view of how networks of biochemical processes shape larger biological phenomena.

FGF-1 induces ATP release from spinal astrocytes in culture and opens pannexin and connexin hemichannels

Spinal astrocytes are coupled by connexin (Cx) gap junctions and express pannexin 1 (Px1) and purinergic receptors. Fibroblast growth factor 1 (FGF-1), which is released in spinal cord injury, activated spinal astrocytes in culture, induced secretion of ATP, and permeabilized them to relatively large fluorescent tracers [ethidium (Etd) and lucifer yellow (LY)] through "hemichannels" (HCs). HCs can be formed by connexins or pannexins; they can open to extracellular space or can form gap junction (GJ) channels, one HC from each cell. (Pannexins may not form gap junctions in mammalian tissues, but they do in invertebrates). HC types were differentiated pharmacologically and by Px1 knockdown with siRNA and by use of astrocytes from Cx43 knockout mice. Permeabilization was reduced by apyrase (APY), an ATPase, and by P2X(7) receptor antagonists, implicating secretion of ATP and autocrine and/or paracrine action. Increased permeability of cells exposed to FGF-1 or ATP for 2 h was mediated largely by Px1 HCs activated by P2X(7) receptors. After a 7-h treatment, the permeability was mediated by both Cx43 and Px1 HCs. FGF-1 also caused reduction in gap junctional communication. Botulinum neurotoxin A, a blocker of vesicular release, reduced permeabilization when given 30 min before FGF-1 application, but not when given 1 h after FGF-1. We infer that ATP is initially released from vesicles and then it mediates continued release by action on P2X(7) receptors and opening of HCs. These changes in HCs and gap junction channels may promote inflammation and deprive neurons of astrocyte-mediated protection in spinal cord trauma and neurodegenerative disease.

Theoretical modeling techniques and their impact on tumor immunology

Currently, cancer is one of the leading causes of death in industrial nations. While conventional cancer treatment usually results in the patient suffering from severe side effects, immunotherapy is a promising alternative. Nevertheless, some questions remain unanswered with regard to using immunotherapy to treat cancer hindering it from being widely established. To help rectify this deficit in knowledge, experimental data, accumulated from a huge number of different studies, can be integrated into theoretical models of the tumor-immune system interaction. Many complex mechanisms in immunology and oncology cannot be measured in experiments, but can be analyzed by mathematical simulations. Using theoretical modeling techniques, general principles of tumor-immune system interactions can be explored and clinical treatment schedules optimized to lower both tumor burden and side effects. In this paper, we aim to explain the main mathematical and computational modeling techniques used in tumor immunology to experimental researchers and clinicians. In addition, we review relevant published work and provide an overview of its impact to the field.

Becoming self-aware: the thymic education of regulatory T cells

The generation of Foxp3(+) regulatory T (Treg) cells in the thymus is essential for immune homeostasis. In the past several years, substantial progress has been made in understanding the mechanisms by which a minor portion of developing thymocytes are selected to become Treg cells. Although previously controversial, recent data support the importance of TCR specificity as a primary determinant for selecting self-reactive thymocytes to become Treg cells in a multi-step process involving cytokines, co-stimulatory molecules, and a variety of antigen-presenting cells. Importantly, the antigenic niche for Treg cell development appears to be typically quite small, implying the recognition of tissue-specific, rather than ubiquitous, self-antigens. Finally, it appears that an NF-κB transcription factor, c-Rel, may be the link between TCR recognition and the induction of Foxp3 expression, which is required for the function and stability of the natural Treg cell population.

Oscillatory control of signalling molecules

The emergence of biological function from the dynamic control of cellular signalling molecules is a fundamental process in biology. Key questions include: How do cells decipher noisy environmental cues, encode these signals to control fate decisions and propagate information through tissues? Recent advances in systems biology, and molecular and cellular biology, exemplified by analyses of signalling via the transcription factor Nuclear Factor kappaB (NF-κB), reveal a critical role of oscillatory control in the regulation of these biological functions. The emerging view is that the oscillatory dynamics of signalling molecules and the epigenetically regulated specificity for target genes contribute to robust regulation of biological function at different levels of cellular organisation through frequency-dependent information encoding.

Systems biology approaches to dissect mammalian innate immunity

Advances in experimental tools have allowed for the systematic identification of components and biological processes as well as quantification of their activities over time. Together with computational analysis, these measurement and perturbation technologies have given rise to the field of systems biology, which seeks to discover, analyze and model the interactions of physical components in a biological system. Although in its infancy, recent application of this approach has resulted in novel insights into the machinery that regulates and modifies innate immune cell functions. Here, we summarize contributions that have been made through the unbiased interrogation of the mammalian innate immune system, emphasizing the importance of integrating orthogonal datasets into models. To enable application of approaches more broadly, however, a concerted effort across the immunology community to develop reagent and tool platforms will be required.

Understanding the temporal codes of intra-cellular signals

The health of organisms and cells depends on appropriate responses to diverse internal and external cues, stimuli, or challenges, such as changes in hormone or cytokine levels, or exposure to a pathogen. Cellular responses must be tailored to the identity and intensity of the stimulus and therefore intra-cellular signals must carry information about both. However, signaling mediators often form intricate networks that react to multiple stimuli yet manage to produce stimulus-specific responses. The multi-functionality ('functional pleiotropism') of signaling nodes suggests that biological networks have evolved ways of passing physiologically relevant stimulus information through shared channels. Increasing evidence supports the notion that this is achieved in part through temporal regulation of signaling mediators' activities. The present challenge is to identify the features of temporal activity profile that represent information about a given stimulus and understand how cells read the temporal codes to control their responses.

Cutting edge: TCR ligation triggers digital activation of NF-kappaB

TCR-mediated activation of the transcription factor NF-κB is required for T cell proliferation, survival, and effector differentiation. Although this pathway is the subject of intense study, it is not known whether TCR signaling to NF-κB is digital (switch-like) or analog in nature. Through analysis of the phosphorylation and degradation of IκBα and the nuclear translocation and phosphorylation of the NF-κB subunit RelA, we show that TCR-directed NF-κB activation is digital. Furthermore, digitization occurs well upstream of the IκB kinase complex, as protein kinase C translocation to the immunologic synapse and activation-associated aggregation of Bcl10 and Malt1 also demonstrate both digital behavior and high correlation with RelA nuclear translocation. Thus, similar to the TCR-to-MAPK signaling cascade, analog Ag inputs are converted to digital activation outputs to NF-κB at an early step downstream of TCR ligation.

High-throughput, single-cell NF-κB dynamics

Single cells in a population often respond differently to perturbations in the environment. Live-cell microscopy has enabled scientists to observe these differences at the single-cell level. Some advantages of live-cell imaging over population-based methods include better time resolution, higher sensitivity, automation, and richer datasets. One specific area where live-cell microscopy has made a significant impact is the field of NF-κB signaling dynamics, and recent efforts have focused on making live-cell imaging of these dynamics more high-throughput. We highlight the major aspects of increasing throughput and describe a current system that can monitor, image and analyze the NF-κB activation of thousands of single cells in parallel.

NF-kappaB: much learned, much to learn

Although NF-kappaB (nuclear factor kappaB) was discovered less than 25 years ago, it is one of the most well-studied transcription factors. It plays critical roles in various immune and cellular processes, and its dysregulation is a major driver of cancer and inflammatory pathology. A new book published by Cold Spring Harbor Laboratory Press, NF-kappaB: A Network Hub Controlling Immunity, Inflammation, and Cancer, provides a series of excellent reviews on the architecture, regulation, and pathophysiological roles of NF-kappaB signaling by some of the leading investigators in the field.