Note to users. If you're seeing this message, it means that your browser cannot find this page's style/presentation instructions -- or possibly that you are using a browser that does not support current Web standards. Find out more about why this message is appearing, and what you can do to make your experience of our site the best it can be.
Response to Comment on "Oscillations in NF-B Signaling Control the Dynamics of Gene Expression"
Single-cell oscillations of signaling proteins can only be detectedby time-lapse imaging of fluorescent proteins. We discussedthe functionality of our NF-B (RelA) and IB fusion constructsin Section A of the supplementary material in (1), and it isclearly important to control for and to minimize overexpression.We previously showed that the IB:RelA ratio affects the timingof translocation responses (2) and that sensitive parameterssuch as NF-Bdependent IB transcription significantlyaffect the system (1, 3). The ability to relate expression levelto behavior of the system [see figure 3 in (1)] allowed us toidentify when ectopic expression perturbs the system and providesinformation to improve computational models.
Using the Hoffmann model (4), we concluded by sensitivity analysis(3) that RelA overexpression would minimally perturb the system;therefore, we used RelA localization as our principal output.We estimated that average overexpression of RelA fusion proteinswas 3 to 5 times that of endogenous RelA levels in transfectedcells, with a distribution in the cell population. Using simulations,Barken et al. (5) suggest that RelA, IB, or RelA and IB expression(within the range suggested by our data) maintains oscillationsbut alters oscillation frequency and amplitude. They suggest[figure 1 in (5)], that simulated 4-fold overexpression of RelAdelays the second peak of nuclear localization by 3.5 hours(normal 2-hour peak delayed to 5.5 hours). This does not fitour experimental data, because in both cell lines used, thesecond peak of nuclear localization with RelA overexpressionalone was at around 3 hours [see figure 1 and supplemental materialin (1)]. Reanalysis of our data from SK-N-AS cells (Fig. 1)and HeLa and Swiss 3T3 cells (data not shown) demonstrates nocorrelation between RelA-DsRed expression level and successivepeak timing (also of amplitude, not shown) with cellular fluorescentlevels that vary up to 20 fold. This clear inconsistency betweentheir in silico results and our experimental data suggests thattheir computational model cannot faithfully reproduce all aspectsof the system. Altering the computational model (1, 4) by replacingthe second-order term for NF-Binduced IB synthesis witha linear expression [reaction 28 in table S1 in (1)] resultsin a new model that reproduces all of the fundamental characteristicsof our experimental data. This model shows much reduced sensitivityof period to RelA concentration (Fig. 2), which demonstratesthat continued refinement of the Hoffmann model (6, 7) may wellbring the simulations closer to observed biological phenomena.
Fig. 1. Correlation of RelA-DsRed expression level with peak timing in single SK-N-AS cells. SK-N-AS cells were transfected with a RelA-DsRed expression vector together with a control EGFP-N1 construct (Clontech). At 24 hours after transfection, cells were treated with 10 ng/ml TNF and imaged by time-lapse confocal microscopy, as described (1). Peak timings (expressed as the time difference between successive peaks) from 30 cells from three separate experiments were plotted against the RelA-DsRed average fluorescence intensity per pixel at the start of the experiment. The number of cells between graphs decreases primarily due to cells that divided during the course of the experiments. These data showed no statistically significant correlation (by linear regression) between initial RelA-DsRed expression level and peak timing. Analysis of absolute peak timing (rather than the difference between peaks) also showed no observable correlation.
[View Larger Version of this Image (18K GIF file)]
Fig. 2. Dependence of oscillation period on RelA concentration using a modified computational model. Simulated nuclear NF-B activity over time normalized to peak 1. NF-B concentrations were increased by factors of 1 (black), 1.5 (red), 2 (green), and 4 (blue) in the model. Simulations were performed using the computational model as in (1) with reduction of the term describing the induced synthesis of IB by NF-B to a linear form [reaction 28 in table S1 in (1)]. [This second-order term has been removed in other variations of the model (5, 6) and is unlikely to be relevant at higher RelA concentrations.] To produce oscillatory behavior, which was representative of the experimental data, appropriate fitting of a single parameter was performed, resulting in parameter tr2 set to 0.0582 min1. All other parameters were unchanged.
[View Larger Version of this Image (21K GIF file)]
Barken et al. (5) show immunocytochemistry (ICC) analysis ofserum-starved 3T3 fibroblasts after TNF stimulation [figure2 in (5)]. We have already described ICC analysis of localizationof endogenous RelA in HeLa cells after TNF stimulation [figureS3 in (1)]. These data showed that the nuclear (N) to cytoplasmic(C) ratio was entirely consistent with the oscillations thatwe saw using fluorescent protein imaging in the same cells.Barken et al. interpret their ICC data as suggesting that synchronousand highly damped oscillations occur in their fibroblast cellsand that the RelA N:C average ratios become less variable. Ourdata suggested that oscillations in the SK-N-AS cells were alsoheavily damped. Surprisingly, Barken et al. comment that "[i]fthe responses of individual cells are asynchronous and showundamped oscillations, one would expect higher variance at latetimes than at early times and lower averages due to prolongedphases between peaks." However, both our previous data and theirnew data clearly show damped oscillations. The anticipated changein variance is therefore irrelevant. It is also not appropriateto directly compare data from different cell lines, becausedamping may differ [e.g., stronger damping in HeLa comparedwith SK-N-AS cells, as shown in (1)].
Barken et al. show experimental data [figure 3 in (5)] fromgenetically engineered cells: IBß//and IB// knockout mouse embryonicfibroblasts expressing only IB or IBß, respectively.Cells expressing only IB showed synchronous oscillations (bynuclear electromobility shift assay), whereas cells in whichthis negative feedback loop was removed (expressing only IBß)showed a single extended phase of NF-BDNA binding. UsingRNAse protection analysis of five NF-Bregulated genes,they showed that IBß-expressing cells exhibit similarkinetics of transcription as IB-expressing cells over an 8-hourperiod. However, these cells lack two IB isoforms, which leadsto a significantly altered NF-B signaling network. The classicalNF-B heterodimer (p65/p50) is predominantly, although not exclusively,regulated by IB. Mice lacking IB show runting and skin defectsand die within 8 days after birth, which demonstrates that NF-Bsignaling is far from normal in these animals (8, 9). The IB-deficientmice display increased, but not constitutive, p65/p50 DNA bindingactivity, which suggests that other mechanisms regulate theactivity of this dimer. It is possible that compensatory changeshave occurred to allow the continued growth of these mutantcell lines, and we cannot exclude the possibility that the activitiesof some of their NF-Bregulating kinases/phosphatasesmay be altered. Redundancy in the NF-B pathway might allow othermembers of the family that are not regulated through IB to counteractthe mutation. Therefore, as with our experiments using leptomycinB to block nuclear export (thus blocking oscillations), theexperiments of Barken et al. are an experimental compromise.Other approaches will help to further elucidate the functionalrole of oscillations.
We do not claim that oscillations are essential for gene expressionat all NF-Bregulated promoters in response to all NF-Bactivatingstimuli. In the same way that calcium oscillations are not aprerequisite for functional calcium signaling, we propose thatoscillations are one feature of NF-B signaling that may be functionallyimportant. One of the major unanswered questions in the fieldis how an apparently identical signal can lead to very differenttranscriptional and cell-fate responses. Oscillations may leadto differential control at different promoters, and this mayvary depending on the rate of NF-B inactivation (e.g., by nucleardephosphorylation). A lowered rate of nuclear inactivation mightmean that oscillations could be less important. We find thatthe rate of dephosphorylation of RelA Ser536 is markedly differentbetween HeLa and SK-N-AS cells [supplementary material in (1)].The presumption that bulk translocation is absolutely linkedto transcription is disputed, and there are examples of normaltranslocation of NF-B without transcriptional up-regulation(10, 11). If the modification of NF-B is the key activationstep, then if NF-B is activated, but IB is not degraded, a lowlevel of endogenous shuttling into the nucleus without bulktranslocation could activate transcription. Alternatively, IBdegradation without NF-B activation could lead to movement ofnonfunctional (or differentially functional) NF-B into the nucleus(12).
As discussed by Lahav (13), the negative feedback intrinsicto this system seems wasteful compared with leaving the NF-Bin the nucleus for as long as the TNF is present. One possibleexplanation is that a message could reside in the oscillationfrequency or oscillations might allow robust independent samplingof the state of the transcription factor at several time points(13). We do consistently observe that the period is more robustthan the amplitude. The use of fluorescent-protein imaging hasbeen a very important tool for the analysis of the dynamicsof protein translocation, oscillations, and transcription incells (1416). Only such single-cell measurements canquantify the period, amplitude, and phase of oscillatory signals.When these are used, we find that oscillations in NF-B expressionare substantial and appear to exert a controlling influenceon the dynamics of gene expression.
D. E. Nelson C. A. Horton V. See J. R. Johnson G. Nelson D. G. Spiller
Centre for Cell Imaging School of Biological Sciences University of Liverpool Bioscience Research Building, Crown Street Liverpool, L69 7ZB, UK
D. B. Kell
Department of Chemistry University of Manchester Post Office Box 88 Sackville Street Manchester, M60 1QD, UK
M. R. H. White
Centre for Cell Imaging School of Biological Sciences University of Liverpool
To whom correspondence should be addressed. E-mail: mwhite{at}liv.ac.uk
Received for publication 28 December 2004. Accepted for publication 8 March 2005.
The editors suggest the following Related Resources on Science sites:
In Science Magazine
TECHNICAL COMMENTS
Derren Barken, Chiaochun Joanne Wang, Jeff Kearns, Raymond Cheong, Alexander Hoffmann, and Andre Levchenko (1 April 2005) Science308 (5718), 52a.
[DOI: 10.1126/science.1107904] |Full Text »|PDF »
REPORTS
D. E. Nelson, A. E. C. Ihekwaba, M. Elliott, J. R. Johnson, C. A. Gibney, B. E. Foreman, G. Nelson, V. See, C. A. Horton, D. G. Spiller, S. W. Edwards, H. P. McDowell, J. F. Unitt, E. Sullivan, R. Grimley, N. Benson, D. Broomhead, D. B. Kell, and M. R. H. White (22 October 2004) Science306 (5696), 704.
[DOI: 10.1126/science.1099962] |Abstract »|Full Text »|PDF »|Supporting Online Material »
THIS ARTICLE HAS BEEN CITED BY OTHER ARTICLES:
Pulsatile Stimulation Determines Timing and Specificity of NF-{kappa}B-Dependent Transcription.
L. Ashall, C. A. Horton, D. E. Nelson, P. Paszek, C. V. Harper, K. Sillitoe, S. Ryan, D. G. Spiller, J. F. Unitt, D. S. Broomhead, et al. (2009)
Science
324, 242-246
|Abstract »|Full Text »|PDF »
Dynamic Effect of Bortezomib on Nuclear Factor-{kappa}B Activity and Gene Expression in Tumor Cells.
M.-H. Sung, L. Bagain, Z. Chen, T. Karpova, X. Yang, C. Silvin, T. C. Voss, J. G. McNally, C. Van Waes, and G. L. Hager (2008)
Mol. Pharmacol.
74, 1215-1222
|Abstract »|Full Text »|PDF »