The IκB-NF-κB Signaling Module: Temporal Control and Selective Gene Activation

See allHide authors and affiliations

Science  08 Nov 2002:
Vol. 298, Issue 5596, pp. 1241-1245
DOI: 10.1126/science.1071914

This article has a correction. Please see:


Nuclear localization of the transcriptional activator NF-κB (nuclear factor κB) is controlled in mammalian cells by three isoforms of NF-κB inhibitor protein: IκBα, -β, and -ɛ. Based on simplifying reductions of the IκB–NF-κB signaling module in knockout cell lines, we present a computational model that describes the temporal control of NF-κB activation by the coordinated degradation and synthesis of IκB proteins. The model demonstrates that IκBα is responsible for strong negative feedback that allows for a fast turn-off of the NF-κB response, whereas IκBβ and -ɛ function to reduce the system's oscillatory potential and stabilize NF-κB responses during longer stimulations. Bimodal signal-processing characteristics with respect to stimulus duration are revealed by the model and are shown to generate specificity in gene expression.

The transcription factor NF-κB regulates numerous genes that play important roles in inter- and intracellular signaling, cellular stress responses, cell growth, survival, and apoptosis (1–3). As such, the specificity and temporal control of gene expression are of crucial physiological interest. Furthermore, the realization of the potential of NF-κB as a drug target for chronic inflammatory diseases or within chemotherapy regimens (4, 5) is dependent on understanding the specificity mechanisms that govern NF-κB–responsive gene expression.

Five related mammalian gene products participate in NF-κB functions (RelA/p65, cRel, RelB, p50, p52), but the predominant species in many cell types is a p65:p50 heterodimer. Its activity is largely controlled by three IκB isoforms (IκBα, -β, and -ɛ) that bind to NF-κB, preventing its association with DNA and causing its localization to the cytoplasm. Signals from various stimuli are transduced to the IκB kinase (IKK) complex, which phosphorylates each IκB isoform, leading to its ubiquitination and proteolysis (6). IκB degradation allows NF-κB to translocate to the nucleus and bind DNA (Fig. 1A). The specific role of each IκB protein in regulating NF-κB is not understood. Mice with an engineered deletion of oneIκB gene show notable molecular compensation by the remaining IκB family members (7,3) and mild phenotypes (8, 9). TheIκBα–/– mouse, however, is perinatal lethal with multi-organ inflammation presumably caused by up-regulation of many NF-κB–responsive genes (10,11). This phenotype is largely rescued by placing the IκBβ coding region under transcriptional control of IκBα (12). IκBα synthesis is controlled by a highly NF-κB–responsive promoter generating autoregulation of NF-κB signaling (13).

Figure 1

Negative feedback and the IκB–NF-κB signaling module. (A) The IκB–NF-κB signaling module. NF-κB is held inactive in the cytoplasm by three IκB isoforms. Cell stimulation activates the IKK complex, leading to phosphorylation and degradation of IκB proteins. Free NF-κB translocates to the nucleus, activating genes, including IκBα. IκBβ and -ɛ are synthesized at a steady rate, allowing for complex temporal control of NF-κB activation involving a negative feedback. (B) A two-component system (x and y) with a negative feedback exhibits dynamic behavior that depends on the relative efficiency of the feedback regulation (β and γ) regulating oscillation persistence versus self-regulation (α and δ) causing oscillation damping. This relationship can be described mathematically as dx/dt = S – αx – βy and dy/dt = γx – δy, where S represents the stimulus. The output, y, ranges from persistent oscillations (green line, high feedback efficiency and no damping, α = δ = 0), to damped oscillations (red line, intermediate feedback efficiency and intermediate damping), to gradual rising to a plateau level (blue line, low feedback efficiency and high damping). In simulations corresponding to the red line, α and δ are 30% of the respective values used for the generation of the blue line. (C) EMSA for NF-κBn in TNFα-stimulated human T cells, human monocytes, and mouse fibroblasts. Nuclear extracts were prepared at the indicated times after the beginning of persistent stimulation with TNFα (10 ng/ml). Equal amounts of nuclear protein were reacted with a radioactively labeled double-stranded oligonucleotide containing a consensus κB site sequence (14). Arrows indicate specific nuclear NF-κB binding activity; asterisks indicate nonspecific DNA binding complexes.

The interactions of IKK, IκB isoforms, and NF-κB can be thought of as a negative feedback–containing signal-transduction module (Fig. 1A) that receives signals from pathways emanating from cell-surface receptors (input) and transmits signals to nuclear promoter-bound protein complexes regulating gene expression (output). In a minimal system with negative feedback (Fig. 1B), cross-regulation between the two signaling components determines the strength of the negative feedback (constants β and γ) and thus the propensity for oscillations, whereas self-regulation (constants α and δ) determines the degree of damping. Depending on the values of these constants, the response to persistent stimulation may thus vary from continual oscillations, to damped oscillations, to a virtually monotonic rise to a plateau level (Fig. 1B). We examined the dynamics of IκB–NF-κB signaling experimentally by measuring nuclear NF-κB (NF-κBn) activity with the electrophoretic mobility shift assay (EMSA) (14). Two phases of NF-κB activation were revealed in response to the stimulation of tumor necrosis factor–α (TNF-α) in various human and mouse cell lines (Fig. 1C). Although the degree and precise timing of postinduction attenuation are variable (15, 16), the output resembles strongly damped oscillations but is distinct from that predicted for a minimal two-component system (Fig. 1B). We conclude that the coordinated degradation, synthesis, and localization of all three IκB isoforms is required to generate the characteristic NF-κB activation profile.

To address their differential functions, we constructed a computational model based on ordinary differential equations. The model behavior depends on the specific values of various control parameters, including those describing (i) the synthesis of each IκB isoform (transcription, mRNA stability, and translation), (ii) the stability of free and NF-κB–bound IκB proteins, (iii) the formation of binary and tertiary IKK–IκB–NF-κB complexes, (iv) the enzymatic rate constants of IKK-containing complexes, and (v) the transport rates affecting localization of each of the components (IκBα, -β, and -ɛ; NF-κB; and derived complexes). Many of the 30 independent model parameters have been previously determined biochemically, and the values of others may be constrained by published data (14). We used reverse genetics to create three IκB–NF-κB signaling modules of reduced complexity and to provide further constraints for effective model parameter fitting.

We engineered mice deleted for IκBβ and -ɛ with the use of standard homologous recombination technology of embryonic stem cells (14). These mice, as well as existing IκBα gene–deleted mice (10), were inter-crossed appropriately to yield embryonic fibroblasts in which the nuclear localization of NF-κB was controlled by a single IκB isoform (14). TNFα stimulation of fibroblasts that contained only the IκBα isoform resulted in a highly oscillatory NF-κB response, with four equally spaced peaks over the course of the 6-hour experiment (Fig. 2A, top). In contrast, in cells harboring only IκBβ (Fig. 2A, center) or -ɛ (Fig. 2A, bottom), NF-κBn increased monotonically to a plateau at 1 hour with no notable subsequent repression.

Figure 2

A computational model based on genetically reduced systems. (A) Analysis of NF-κBn by EMSAs of nuclear extracts prepared at indicated times after stimulation with TNF-α (10 ng/ml) of fibroblasts of the indicated genotype. Arrows indicate specific nuclear NF-κB binding activity; asterisks indicate nonspecific DNA binding complexes. (B) The NF-κB–specific mobility shift in cells of the indicated genotype was quantitated by phosphoimager and normalized and graphed against a linear time scale. (C) Computational modeling of each genetically simplified signaling module results in characteristic kinetics of the NF-κBn response. Model-fitting allows previously undetermined biochemical parameters to be estimated. (D) Models of the simplified signaling modules are combined, with increasing IκBβ and -ɛ transcription rates, while keeping the IκBα transcription rate constant. Model behaviors are shown that result as the constitutive mRNA synthesis parameters for IκBβ and IκBɛ are increased fivefold (top to middle) and then sevenfold (middle to bottom). The bottom panel represents the NF-κBn output predicted by a model with mRNA synthesis parameters identical to those employed in the single IκB isoform models shown inFig. 2C. (E) Biochemical analysis of NF-κB and IκB isoforms in wild-type fibroblasts. NF-κBn (top) assayed by EMSA at the indicated times after persistent stimulation with TNF-α. The specific NF-κB–specific mobility shift was quantitated by phosphoimager and normalized and graphed at the indicated nonlinear time scale. Western blots of corresponding cytoplasmic fractions are probed with anti-bodies specific to IκBα and -β (bottom) and IκBɛ (above). (F) Verifica- tion of the computational model for wild-type cells. IκBα and -β mRNA synthesis parameters were determined by qualitative model fitting to yield the graphed outputs in response to persistent stimulation of NF-κBn (top) and total cellular concentrations of IκBα, -β, and -ɛ (lower panels as indicated).

The simplified computational models of β–/– ɛ–/–, α–/–β–/– , and α–/– ɛ–/– cells each contain 17 partially overlapping parameters that control the output, 9 of which have been biochemically determined previously. In addition, five transport parameters are constrained by ratios derived from steady state nuclear and cytoplasmic localization (17,18). Their absolute values, as well as IκB synthesis parameters, were determined by semiquantitative fitting of the model outputs to the experimental data.

Combining the three models of genetically reduced signaling modules results in a computational model of the IκB–NF-κB signaling module in wild-type cells. Varying their relative contributions revealed discrete functional roles for the mammalian IκB proteins in NF-κB regulation (Fig. 2D). IκBα mediates rapid NF-κB activation and strong negative feedback regulation, resulting in an oscillatory NF-κB activation profile. IκBβ and -ɛrespond more slowly to IKK activation and act to dampen the long-term oscillations of the NF-κB response. The interplay between these isoforms can result in remarkably rapid responses to the onset or cessation of stimulation and can allow a relatively stable NF-κB response during long-term stimulation (Fig. 2D, middle). In the absence of strong damping mechanisms, high negative feedback efficiency can lead to long-term oscillations in the output (Fig. 2D, top), as observed, for example, in p53-mdm2 cross-regulation (19). In wild-type NF-κB signaling, pronounced oscillations are absent and may be detrimental, asIκBβ–/– IκBɛ–/–females have a dramatically shortened fertility span (3). Our results also suggest that varying relative synthesis levels of IκBα, -β, and -ɛ may constitute a mechanism for altering the responsiveness of the NF-κB signal transduction pathway, a mechanism that cells or cell lineages may use in response to environmental or developmental cues.

To determine the individual contributions of the single-IκB isoform models in wild-type fibroblasts, we analyzed nuclear and cytoplasmic extracts during TNF-α stimulation by EMSA for NF-κBn (Fig. 2E, top) and by Western blot for IκB proteins (Fig. 2E, lower panels). NF-κBn first appeared within 5 minutes, as IκBα levels rapidly decreased, and was at its maximum level at 30 min, when cytoplasmic IκBβ and -ɛ disappeared. After 60 min, κB-binding activity was reduced concomitantly with a transient increase in IκBα protein levels. After 2 hours, NF-κBn returned to about maximum levels and then decreased to about half of maximum levels as IκBα levels stabilized and IκBβ and -ɛproteins reappeared. Model-fitting allowed the determination of IκBβ and -ɛ mRNA synthesis parameters that resulted in an optimal fit of the reconstituted model (Fig. 2F) with experimental wild-type responses. These parameters were about sevenfold lower than those determined for the respective single-IκB models, suggesting previously unrecognized cross-regulation in the expression of IκB genes. The outputs of the resulting “wild-type” model describe NF-κBn (Fig. 2F, top) and cellular IκB isoform levels (Fig. 2F, lower panels) in good qualitative and quantitative agreement with experimental data, thus justifying its use as a predictive tool in experimentation.

We then investigated how transient stimuli control NF-κB activation. Simulated short pulses led to transient NF-κB activation responses, whose duration was insensitive to the duration of the stimulus within the first hour (Fig. 3A). Measuring NF-κBn in transiently stimulated cells confirmed this prediction: TNF-α pulses of 5, 15, 30, and 60 min led to DNA binding activity profiles of similar duration, equivalent to the first peak of persistently stimulated NF-κB activity (Fig. 3B). To explore the signal-processing characteristics of the IκB–NF-κB signaling module, we considered the availability of NF-κBn for binding to an arbitrary κB-responsive promoter. We plotted the duration of NF-κBn availability above an arbitrarily chosen threshold concentration presumed to allow for efficient binding of the κB element of some promoters (Fig. 3C) as a function of the duration of TNF-α stimulation. For long stimulations, NF-κBn lasts as long as the stimulus. For stimulations of less than 1 hour, the duration of the response is largely invariant. Hence, the IκB–NF-κB signaling module has bimodal signal processing characteristics: One mode of signal processing ensures that even short stimulations result in substantial NF-κB responses; the other mode, operative at stimulations longer than 1 hour, generates responses proportional in duration to the stimulus.

Figure 3

NF-κB responses to transient stimulation of varying duration. (A) Modeling the NF-κB response to stimuli of various durations. Temporal profile of NF-κBn predicted by the computational model of the pathway in the wild-type cells in response to stimulations of the same intensity but varying durations as indicated. (B) Experimental NF-κBn data for transient TNF stimulation regimes. Each panel shows the results from EMSAs with nuclear extracts after the onset of a transient stimulation with TNF-α for 5, 15, 30, and 60 min. (C) Graph of the duration of above-threshold (20 nM) NF-κBn as a function of the duration of the transient stimulus as predicted by the computational model.

One implication of the above-described temporal regulatory switch is that some NF-κB–responsive genes might be efficiently activated by a stimulation pulse as short as 15 min (Fig. 4A, blue stippled line). However, there may be other genes that require longer (>1 hour) exposure to NF-κB. We tested this hypothesis by examining the behavior of two NF-κB–regulated genes that display different transcriptional activation profiles (3). One gene encodes the chemokine IP-10, which displays detectable mRNA levels within 30 min of TNF-α stimulation (Fig. 4B) (14). A 15-min pulse of TNF-α was sufficient to activate IP-10 transcription, although, as might be expected, persistent stimulation yielded greater amounts of IP-10 mRNA at later time points. In contrast, the chemokine gene RANTES required at least 2 hours of stimulation for detectable expression and was not induced by a shorter, transient stimulation (Fig. 4B). The model predicts that in cells lacking IκBα, NF-κB will have a longer nuclear lifetime after a short TNF-α pulse (Fig. 4A, red lines) than in wild-type cells. Indeed, transient stimulation led to an extended peak of κB-binding activity, similar to responses resulting from longer stimulations (compare upper panels in Fig. 4, B and C). In turn, RANTES gene transcription was induced not only with persistent stimulation but also with transient stimulations as short as 15 min in IκBα–/– cells (Fig. 4C).

Figure 4

Temporal control of NF-κB has qualitative effects on gene regulation. (A) Modeling NF-κBn concentrations in response to transient (15 min, blue) and persistent (red) stimulation in wild-type (dashed lines) and IκBα–/– (solid lines) cells. (B andC) Transcriptional NF-κB responses to persistent or transient stimulation with TNF-α. EMSAs (top) monitor NF-κBn after the onset of persistent (left) or transient (right) stimulation of wild-type (B) or IκBα–/– (C) fibroblasts with TNF-α. Ribonuclease protection assays (bottom) monitor the transcript levels of chemokine genes RANTES and IP-10 as well as the housekeeping gene glyceraldehyde phosphate dehydrogenase (GAPDH) at indicated times (in min and hours) after the onset of persistent (left) or transient (right) stimulation of wild-type (B) or IκBα–/– (C) fibroblasts with TNF-α.

Thus, the bimodal temporal signal-processing characteristics of the IκB–NF-κB module result in not only quantitative but also qualitative regulation of gene expression defining two classes of genes: those that require persistent NF-κB activation and those that do not. The latter genes (e.g., IP-10) undergo a standard activation program irrespective of the precise duration of a transient TNF-α stimulus. IκB–NF-κB signal processing ensures that even very short stimulations can produce easily detectable transcriptional activation of such genes. Conversely, the former genes (e.g., RANTES) may or may not require expression of NF-κB–induced transcription factors with which NF-κB must synergize for gene induction. RANTES induction is protein synthesis–independent (3), and although the mechanism for the apparent delay in RANTES induction is unknown, it is noteworthy that histone H4 acetylation in macrophages was shown to be lipopolysaccharide-inducible on the RANTES promoter while it was constitutive on many other chemokine promoters (20). If H4 acetylation, a marker for chromatin accessibility, precedes NF-κB binding to the endogenous promoter, it may be induced by a NF-κB–independent pathway. Our results imply that specificity in gene expression can be achieved by using two signal transduction pathways in temporally distinct ways. Gene induction only occurs when the two pathways are temporally coordinated.

Our analysis has revealed the IκB–NF-κB signaling module as a biological system that regulates cellular behavior through the control of system dynamics. The generation of genetically reduced systems enabled the computational analysis of complex dynamic behavior. Exploration of the computational model, in turn, provided insights into the physiologically relevant differential functions of heretofore seemingly redundant system components. Their distinct but coordinated regulation in synthesis and degradation allows for a transcriptional response system with signal-processing characteristics that exhibit both rapid signal responsiveness and stable long-term responses.

Supporting Online Material

Materials and Methods

SOM Text

  • * These authors contributed equally to this work.

  • Present address: Biogen, Incorporated, Cambridge, MA 02142, USA.

  • To whom correspondence should be addressed. E-mail: baltimo{at}


View Abstract

Stay Connected to Science

Navigate This Article