Institute of Fundamental Technological Research

White's lab established that strong, continuous stimulation with tumour necrosis factor- () can induce sustained oscillations in the subcellular localisation of the transcription factor nuclear factor (NF-). But the intensity of the signal varies substantially, from picomolar in the blood plasma of healthy organisms to nanomolar in diseased states. We report on a systematic survey using computational bifurcation theory to explore the relationship between the intensity of stimulation and the existence of sustained NF- oscillations. Using a deterministic model developed by Ashall et al. in 2009, we find that the system's responses to are characterised by a supercritical Hopf bifurcation point: above a critical intensity of the system exhibits sustained oscillations in NF-kB localisation. For below this critical value, damped oscillations are observed. This picture depends, however, on the values of the model's other parameters. When the values of certain reaction rates are altered the response of the signalling pathway to stimulation changes: in addition to the sustained oscillations induced by high-dose stimulation, a second oscillatory regime appears at much lower doses. Finally, we define scores to quantify the sensitivity of the dynamics of the system to variation in its parameters and use these scores to establish that the qualitative dynamics are most sensitive to the details of NF- mediated gene transcription.

NF-kB signalling pathway, Parameter sensitivity, Bifurcation analysis, Oscillations

Background

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.
Results

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.
Conclusions

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.