Interactions among oscillatory pathways in NF-kappa B signaling
© Wang et al; licensee BioMed Central Ltd. 2011
Received: 26 April 2010
Accepted: 3 February 2011
Published: 3 February 2011
Sustained stimulation with tumour necrosis factor alpha (TNF-alpha) induces substantial oscillations—observed at both the single cell and population levels—in the nuclear factor kappa B (NF-kappa B) system. Although the mechanism has not yet been elucidated fully, a core system has been identified consisting of a negative feedback loop involving NF-kappa B (RelA:p50 hetero-dimer) and its inhibitor I-kappa B-alpha. Many authors have suggested that this core oscillator should couple to other oscillatory pathways.
First we analyse single-cell data from experiments in which the NF-kappa B system is forced by short trains of strong pulses of TNF-alpha. Power spectra of the ratio of nuclear-to-cytoplasmic concentration of NF-kappa B suggest that the cells' responses are entrained by the pulsing frequency. Using a recent model of the NF-kappa B system due to Caroline Horton, we carried out extensive numerical simulations to analyze the response frequencies induced by trains of pulses of TNF-alpha stimulation having a wide range of frequencies and amplitudes. These studies suggest that for sufficiently weak stimulation, various nonlinear resonances should be observable. To explore further the possibility of probing alternative feedback mechanisms, we also coupled the model to sinusoidal signals with a wide range of strengths and frequencies. Our results show that, at least in simulation, frequencies other than those of the forcing and the main NF-kappa B oscillator can be excited via sub- and superharmonic resonance, producing quasiperiodic and even chaotic dynamics.
Our numerical results suggest that the entrainment phenomena observed in pulse-stimulated experiments is a consequence of the high intensity of the stimulation. Computational studies based on current models suggest that resonant interactions between periodic pulsatile forcing and the system's natural frequencies may become evident for sufficiently weak stimulation. Further simulations suggest that the nonlinearities of the NF-kappa B feedback oscillator mean that even sinusoidally modulated forcing can induce a rich variety of nonlinear interactions.
Nuclear factor kappa B (NF-κ B) transcription factors are critical to the control of response to cellular stress and are also involved in the regulation of cell-cycle/growth, survival, apoptosis, inflammation and immunity [1–5]. They are dimeric molecules, composed of either homo-or hetero-dimers, with the most common form being the RelA:p50 hetero-dimer.
In resting cells, NF-κ B dimers are sequestered in the cytoplasm by members of a family of molecules, the inhibitors of κ B or Iκ Bs. When the cells are stimulated with tumour necrosis factor alpha (TNFα) certain kinases, the IKKs, are activated and phosphorylate Iκ Bs, marking them for degradation. The degraded Iκ Bs then release NF-κ B dimers which are free to translocate into the nucleus, where they bind to specific sequences in the promoter or enhancer regions of target genes, including those for Iκ Bα and the zinc finger protein A20 [6, 7]. Newly synthesized Iκ Bα migrates to the nucleus, binds to NF-κ B dimers and removes them from the nucleus, while A20 protein stays in cytoplasm and represses the activity of TNFα receptors . Hence the NF-κ B system includes at least two negative feedback loops, one involving cytoplasmic sequestration mediated by Iκ Bα and another involving A20.
Coupled negative feedback loops are known to have the potential to support oscillations [9, 10]. Biological evidence suggesting NF-κ B oscillations was reported by Hoffmann et al. , who did population-level studies of Iκ Bα-/-embryonic fibroblasts using electro-mobility shift assays, and also observed by Nelson et al.  in the cellular concentration and nuclear: cytoplasmic localisation of Iκ Bα and RelA-fluorescent fusion proteins in single SK-N-AS cells. Single cell time-lapse imaging data showed persistent cycling of NF-κ B localisation between the cytoplasm and nucleus of SK-N-AS cells in response to continual TNFα stimulation [12, 13].
The expression of genes regulated by NF-κ B is tightly coordinated with the activities of many other signalling and transcription-factor pathways [14–18] including the p53 signalling pathway. Though the canonical NF-κ B signalling pathway has been studied extensively, the existence and mechanisms of the interactions between the NF-κ B pathway and other signalling pathways are unclear.
Ashall et al.  reported that different patterns of pulsatile stimulation of the NF-κ B system lead to different patterns of NF-κ B dependent gene expression, supporting the view that the frequency of the oscillations may have a functional role. It is therefore important to understand how the system's response frequencies are influenced by interactions with other oscillatory pathways. If both the NF-κ B network and its couplings to other other oscillatory pathways were purely linear, then it would be straightforward to use the machinery of transfer functions to characterise their interactions. In particular, one would expect the power spectrum of a periodically stimulated system to have its power concentrated at the forcing signal's frequency and its harmonics. But the NF-κ B network is a highly nonlinear system of coupled chemical reactions and hence, as we will demonstrate in modelling studies below, its response to periodic stimulation can include complex interactions between the intrinsic negative feedback oscillator and the stimulus. In what follows we work with a deterministic model first described in  and examine the power spectral densities of the time courses of NF-κ B localisation when the model is subjected to two types of periodic signals: trains of rectangular pulses and sinusoidal signals. Both sorts of stimulus are intended as a proxies for the influence of other oscillatory pathways on the core NF-κ B feedback loop and we find that such interactions can produce a rich variety of nonlinear dynamical behaviour. Our results on sinusoidal forcing are in broad agreement with those of Fonslet et al. , who applied a related family of sinusoidal stimulation protocols to a reduced model of the core NF-κ B oscillator originally developed in .
Results and Discussion
Experimental phenomena and time series analysis
The bottom row of Figure 1 shows the corresponding power spectral densities and makes the point that the intrinsic oscillation—the one induced by constant stimulation—does not appear to have been excited by these pulsed stimulation protocols: in all cases the peaks in the power spectral density appear at multiples of the stimulus frequency. In mathematical terms, the cells respond as though they have been perturbed away from a stable resting state and do not offer evidence of any more complex internal dynamics. Nonetheless, these data provided crucial input to the development of the models in . In particular, note that the amplitude of the response to the second and third pulses of stimulation varies as a function of the inter-pulse period, being somewhat reduced in panels (a) and (b), but not in (c): this constrains the rate at which the model should relax toward equilibrium.
One possibility is that the pulsing signal is too strong and suppresses the response at the natural frequency. In the following section we report in numerical experiments, based on the deterministic model in , through which we investigate this possibility by simulating the system's response to a periodic train of square pulses with strength varying from 0 to 10 ng/ml and inter-pulse intervals ranging from 20 to 295 minutes.
Model Introduction and Bifurcation Analysis
Here we have suppressed the time dependence of the concentrations of the various forms of IKK and of the A20 protein. The k* are fixed parameters and T R is a parameter indicating the presence or absence TNFα stimulation: T R = 1 when the system is being stimulated with 10 ng/ml TNFα and T R = 0 otherwise. This is a simplification: TNFα does not act directly on IKK, but rather produces its effect through a cascade of chemical reactions that begins when TNFα binds its receptor at the cell surface.
Following Ashall et al., we treat only the final stages of this chain. Thus when we model strong stimulation by setting T R = 1, we abstract away a great deal of transduction machinery. In the studies that follow we will want to explore the consequences of weaker stimulation (TNFα concentrations of 10 ng/ml are orders of magnitude higher than physiological levels: Matalka et al.  report measurements from various tissues in healthy mice and found concentrations in the range 1-5 pg/ml while Prabha et al.  studied the plasma of healthy human subjects and found TNFα concentrations on the order of 100 pg/ml.) and it is natural to generalize the role of T R , allowing it to vary across the interval 0 ≤ T R ≤ 1. In light of the modelling assumptions discussed above, one should not imagine that T R depends linearly on the concentration of TNFα, but only that T R increases monotonically with dose.
All the numerical experiments reported here begin from the same initial condition, whose preparation is described in the Methods section. We simulated the consequences of two TNFα stimulation protocols: pulse-like stimulation similar to that used in Nelson's experiments and sinusoidally modulated stimulation.
From the experiments and simulations in , we know that stimulation with a finite train of three strong pulses tends to entrain the cell's response. In this section, we first show that when subjected to a long periodic train of strong pulses, the model's response frequencies are also entrained, but that when the model is driven by sufficiently weak periodic pulse trains, various resonance phenomena connected to the underlying natural oscillations become observable.
When T R ≥ 0.2 the power spectral density of the response shows power only at integer multiples of the pulsing frequency, indicating that the response is fully entrained by the forcing. Given that the pulse-strength used in the experiments corresponds to T R = 1, these modelling results are in qualitative agreement with the experimental data in Figure 1: the stimulation was so strong that we should have expected it to have entrained the responses completely.
Sinusoidally modulated stimulation
The case with η = 1 is the straightforward substitution of pulse trains with sinusoidal waves having peak-to-trough amplitude 2ε. Our motivations for this formulation are twofold: firstly, provided η < 1 we have T R (t) > 0 at all times and so expect the system more readily to exhibit oscillations similar to those induced by constant stimulation. Secondly, periodic forcing of the form (2) converts our model into a sinusoidally forced nonlinear oscillator, a class of systems that has been studied very extensively (see, for example, Pikovsky et al. , Wiggins  or Nayfeh and Mook ).
where p and q are integers.
This idea allows us to give precise definitions for the terms sub-and superharmonic resonance used above: the former corresponds to resonances where p = 1 and q > 1 in (4), while the latter corresponds to q = 1 and p > 1. Finally, when both η and the difference (ν -ν0) are large, the system's response can become chaotic, so that the features in its power spectral density bear no simple relationship to the frequencies ν and ν0.
which gives rise to the network of lines with slope of ±1. Finally, Figure 7 also exhibits sub-harmonic resonance: vertical strips for which ν ≈ qν0 (with q a whole number) show power concentrated at the forcing frequency, but also at frequencies f = ν/p, where p is a whole number. This is perhaps clearest in the strip ν ≈ 2ν0, where the strongest spectral feature lies along the line f = ν/2 ≈ ν0.
where p and q are whole numbers. Generally speaking these periodic responses—which are a nonlinear generalization of the familiar phenomenon of resonance in linear systems—are easiest to excite when the numbers p and q are small: the tongues in Figure 9 are labelled by ratios q:p where p and q are as in (5). Although Figure 9 illustrates this story for the specific family of modulations with ε = 0.5, the qualitative picture is essentially the same for all values in the range of 0.366 < ε < 0.5: in all cases there will be narrow tongues of parameter combinations (ν, η) for which the response is periodic with a period given by a resonance relation like (5). For modulations whose parameters lie outside the tongues the qualitative behaviour of the response will be more complex, having a power spectrum qualitatively similar to those in Figure 7, with features at harmonics of the modulation frequency, at harmonics of the natural frequency ν0 (which varies with ε) and at near-resonant combinations (3) of the two.
When TR < 0.366, the unforced system is a damped oscillator and so, when subjected to periodic forcing, may exhibit resonant phenomena. By using the same numerical simulation and power spectral density analyses as we did for the case ε = 0.5, we find that the responses of the pathway to forcing with with T R < 0.366 can be divided into two groups: those in which harmonic, sub-and superharmonic resonances can be observed (in the region 0 < ε ≤ 0.01), and those for which only harmonic and subharmonic resonance can be observed.
We have focussed here on periodic and quasiperiodic oscillatory interactions amenable to power-spectral analyses of the sort illustrated in Figures 4, 5, 6, 7, 8 and, for single-cell recordings, in the lower panels of Figure 1. But there is every reason to expect a much richer range of dynamical behaviour: Fonslet et al. , who applied forcing of the form (2) with η = 1 to a simplified NF-κ B model, saw evidence of a period-doubling cascade as well as chaos and strange attractors. Complete analysis of these more complex dynamical regimes requires tools beyond the power spectrum, and so we defer their exploration to a future paper.
Given that intrinsically nonlinear chemical kinetics underpin cell-signalling networks, one shouldn't expect these systems to be linear and any analysis of experimental time series, whether in the time or frequency domain, must take this intrinsic nonlinearity into account. Our modelling studies suggest that coupling even the simplest, sinusoidal signal into the NF-κ B network can give rise to a host nonlinear phenomena, including harmonic, subharmonic and superharmonic resonances as well as quasiperiodic and even chaotic behaviour.
The simulation studies reported here used an externally-imposed oscillatory forcing as a proxy for interactions between the core NF-κ B feedback loop and other oscillatory networks. Our results suggest that interactions between the NF-κ B oscillator and other oscillatory pathways can give rise to extremely rich temporal signalling programs and so, perhaps, to many distinct patterns of expression for target genes.
while all others start at zero. We then integrated the ODEs through 2000 minutes of simulated time and used the result as an initial condition for the various TNFα-stimulation studies.
The bifurcation diagram in Figure 3 and the plot of the Arnol'd tongues in Figure 9 were prepared using Bard Ermentrout's XPPAUT , a front end for the powerful numerical bifurcation package AUTO . All our other numerical work used MATLAB© to integrate the ODEs, compute the power spectra and plot the figures.
Appendix: circle maps
Here we discuss a standard mathematical tool, the circle map, used to study the response of a nonlinear oscillator subjected to periodic forcing. The idea is to consider sufficiently weak forcing that the response is close to that of the unforced system, and then define a phase angle ϕ that is close to the phase of the corresponding unforced oscillator.
where τ = 2π/ω is the period of the forcing, ε is a measure of its amplitude and G(ϕ) is a phase-dependent function that characterizes the response.
That is, rational rotation numbers correspond to periodic dynamics: the orbit of ϕ0 is called a (p, q) cycle. On the other hand, irrational rotation numbers correspond to quasiperiodic dynamics [27, 31, 32].
Where . This implies that forced oscillators whose corresponding circle maps have irrational rotation number never repeat periodically.
Results of Arnold  show that the qualitative behaviour of the circle map (6) is stable against arbitrary small perturbations if and only if the rotation number is rational. We can thus expect that if F has a rational rotation number , there exists a region of parameter values (α, ε) such that the all the forced systems whose parameters lie in this region share the same rational rotation number: that is, they all have the same sorts of periodic response. Such regions of parameter values are called Arnol'd tongues.
where T = 2π/ω is the period of the forcing. Although the function G' in (9) is not as simple as the the G (ϕ) in (6), the map F' still has (p, q) periodic responses and a corresponding system of Arnol'd tongues: Figure 9 shows examples.
which has some solutions ϕ ∈ [0, 2π] provided that |2πβ| ≤ ε or |β| ≤ ε/2π.
CAH and YW were supported by BBSRC grant BBD0088081. PP and MRHW were supported by MRC Grant G0500346. PP, MRM, DSB and MRHW were supported by BBSRC grants BBF005381 & BBF00561X1.
- Hayden M, Ghosh S: Signaling to NF- κ B. Genes & Development 2004, 18: 2195-2224. [http://dx.doi.org/10.1101/gad.1228704]View ArticleGoogle Scholar
- Bonizzi G, Karin M: The two NF- κ B activation pathways and their role in innate and adaptive immunity. Trends in Immunology 2004,25(6):280-288. 10.1016/j.it.2004.03.008View ArticlePubMedGoogle Scholar
- Pasparakis M, Luedde T, Schmidt-Supprian M: Dissection of the NF- κ B signalling cascade in transgenic and knockout mice. Cell Death & Differentiation 2006, 13: 861-872.View ArticleGoogle Scholar
- Pahl HL: Activators and target genes of Rel/NF- κ B transcription factors. Oncogene 1999,18(49):6853-6866. [http://www.nature.com/onc/journal/v18/n49/abs/1203239a.html] 10.1038/sj.onc.1203239View ArticlePubMedGoogle Scholar
- Ghosh S, May M, Kopp E: NF- κ B and rel proteins: Evolutionarily conserved mediators of immune responses. Annual Review of Immunology 1998, 16: 225-260. 10.1146/annurev.immunol.16.1.225View ArticlePubMedGoogle Scholar
- Scott ML, Fujita T, Liou HC, Nolan GP, Baltimore D: The p65 subunit of NF-kappa B regulates I kappa B by two distinct mechanisms. Genes & Development 1993,7(7a):1266-1276. [http://genesdev.cshlp.org/content/7/7a/1266.abstract]View ArticleGoogle Scholar
- Song HY, Rothe M, Goeddel DV: The tumor necrosis factor-inducible zinc finger protein A20 interacts with TRAF1/TRAF2 and inhibits NF-kappaB activation. Proceedings of the National Academy of Sciences of the United States of America 1996,93(13):6721-6725. 10.1073/pnas.93.13.6721PubMed CentralView ArticlePubMedGoogle Scholar
- Wertz IE, O'Rourke KM, Zhou H, Eby M, Aravind L, Seshagiri S, Wu P, Wiesmann C, Baker R, Boone DL, Ma A, Koonin EV, Dixit VM: De-ubiquitination and ubiquitin ligase domains of A20 downregulate NF- κ B signalling. Nature 2004, 430: 694-699. 10.1038/nature02794View ArticlePubMedGoogle Scholar
- Pigolotti S, Krishna S, Jensen MH: Oscillation patterns in negative feedback loops. Proceedings of the National Academy of Sciences 2007,104(16):6533-6537. 10.1073/pnas.0610759104View ArticleGoogle Scholar
- Thomas R: Logical analysis of systems comprising feedback loops. Journal of Theoretical Biology 1978,73(4):631-656. 10.1016/0022-5193(78)90127-3View ArticlePubMedGoogle Scholar
- Hoffmann A, Levchenko A, Scott M, Baltimore D: The I kappa B-NF-kappa B signaling module: Temporal control and selective gene activation. Science 2002,298(5596):1241-1245. 10.1126/science.1071914View ArticlePubMedGoogle Scholar
- Nelson DE, Ihekwaba AEC, Elliott M, Johnson JR, Gibney CA, Foreman BE, Nelson G, See V, Horton CA, Spiller DG, Edwards SW, McDowell HP, Unitt JF, Sullivan E, Grimley R, Benson N, Broomhead DS, Kell DB, White MRH: Oscillations in NF-kappa B signaling control the dynamics of gene expression. Science 2004,306(5696):704-708. 10.1126/science.1099962View ArticlePubMedGoogle Scholar
- Horton CA: Computational Modelling of Cell Signalling Pathway Dynamics. PhD thesis. School of Biological Sciences, University of Liverpool; 2006.Google Scholar
- Garside H, Stevens A, Farrow S, Normand C, Houle B, Berry A, Maschera B, Ray D: Glucocorticoid ligands specify different interactions with NF-kappa B by allosteric effects on the glucocorticoid receptor DNA binding domain. Journal of Biological Chemistry 2004,279(48):50050-50059. 10.1074/jbc.M407309200View ArticlePubMedGoogle Scholar
- Ikeda A, Sun X, Li Y, Zhang Y, Eckner R, Doi T, Takahashi T, Obata Y, Yoshioka K, Yamamoto K: p300/CBP-dependent and -independent transcriptional interference between NF-kappa B RelA and p53. Biochemical and Biophysical Research Communications 2000,272(2):375-379. 10.1006/bbrc.2000.2786View ArticlePubMedGoogle Scholar
- Perkins ND: Integrating cell-signalling pathways with NF- κk B and IKK function. Nature Reviews Molecular Cell Biology 2007, 8: 49-62. 10.1038/nrm2083View ArticlePubMedGoogle Scholar
- Salminen A, Ojala J, Huuskonen J, Kauppinen A, Suuronen T, Kaarniranta K: Interaction of aging-associated signaling cascades: Inhibition of NF-kappa B signaling by longevity factors FoxOs and SIRT1. Cellular And Molecular Life Sciences 2008,65(7-8):1049-1058. 10.1007/s00018-008-7461-3View ArticlePubMedGoogle Scholar
- Wadgaonkar R, Phelps K, Haque Z, Williams A, Silverman E, Collins T: CREB-binding protein is a nuclear integrator of nuclear factor-kappa B and p53 signaling. Journal of Biological Chemistry 1999,274(4):1879-1882. 10.1074/jbc.274.4.1879View ArticlePubMedGoogle Scholar
- Ashall L, Horton CA, Nelson DE, Paszek P, Harper CV, Sillitoe K, Ryan S, Spiller DG, Unitt JF, Broomhead DS, Kell DB, Rand DA, See V, White MRH: Pulsatile stimulation determines timing and specificity of NF-kappaB-dependent transcription. Science 2009,324(5924):242-246. 10.1126/science.1164860PubMed CentralView ArticlePubMedGoogle Scholar
- Fonslet J, Rud-Petersen K, Krishna S, Jensen MH: Pulses and Chaos: Dynamical Response in a Simple Genetic Oscillator. International Journal of Modern Physics B 2007,21(23-24):4083-4090. 10.1142/S0217979207045256View ArticleGoogle Scholar
- Krishna S, Jensen MH, Sneppen K: Minimal model of spiky oscillations in NF- κ B signaling. PNAS 2006,103(29):10840-10845. 10.1073/pnas.0604085103PubMed CentralView ArticlePubMedGoogle Scholar
- Nelson DE: Charaterisation and Manipulation of NF- κ B Signal Oscillations in Single Living Cells. PhD thesis. School of Biological Sciences, University of Liverpool; 2005.Google Scholar
- Matalka KZ, Tutunji MF, Abu-Baker M, Abu Baker Y: Measurement of protein cytokines in tissue extracts by enzyme-linked immunosorbent assays: Application to lipopolysaccharide-induced differential milieu of cytokines. Neuroendocrinology Letters 2005,26(3):231-236.PubMedGoogle Scholar
- Prabha C, Jalapathy KV, Matsa RP, Das SD: Role of TNF-alpha in host immune response in tuberculous pleuritis. Current Science 2003,85(5):639-642.Google Scholar
- Wang Y, Paszek P, Horton CA, Hong Y, White MRH, Kell DB, Muldoon MR, Broomhead DS: A systematic survey of the response of a model NF- κ B signalling pathway to TNF α stimulation. MIMS EPrint 2009.104 Manchester Institute for Mathematical Sciences, The University of Manchester; 2010. [http://eprints.ma.man.ac.uk/1404/]Google Scholar
- Wiggins S: Introduction to Applied Nonlinear Dynamical Systems and Chaos. New York; London: Springer Verlag; 1990.View ArticleGoogle Scholar
- Pikovsky A, Rosenblum M, Kurths J: Synchronization : A Universal Concept in Nonlinear Sciences. Cambridge, UK: Cambridge University Press; 2003.Google Scholar
- Nayfeh AH, Mook DT: Nonlinear Oscillators. New York: John Wiley & Sons; 1979.Google Scholar
- Ermentrout B: Simulating, Analyzing, and Animating Dynamical Systems: A Guide to XPPAUT for Researchers and Students. Philadelphia: SIAM; 2002.View ArticleGoogle Scholar
- Doedel EJ, Paffenroth RC, Champneys AR, Fairgrieve TF, Kuznetsov YA, Oldeman BE, Sandstede B, Wang XJ: AUTO2000 : Software for continuation and bifurcation problems in ordinary differential equations. In Tech rep. California Institute of Technology, Department of Applied and Computational Mathematics, Pasadena CA 91125; 2000.Google Scholar
- Arnold VI: Geometrical Methods in the Theory of Ordinary Differential Equations. 2nd edition. New York; London: Springer Verlag; 1988.View ArticleGoogle Scholar
- Arrowsmith DK, Place CM: An introduction to Dynamical Systems. Cambridge, UK: Cambridge University Press; 1990.Google Scholar
- Denjoy A: Sur les courbes definies par les equations differentielles a la surface du tore. J Math 1932,17(IV):333-375.Google Scholar