- Niels Bohr Institute, University of Copenhagen, Copenhagen, Denmark
The transcription factor NF-κB plays a vital role in the control of the immune system, and following stimulation with TNF-α its nuclear concentration shows oscillatory behaviour. How environmental factors, in particular temperature, can control the oscillations and thereby affect gene stimulation is still remains to be resolved question. In this work, we reveal that the period of the oscillations decreases with increasing temperature. We investigate this using a mathematical model, and by applying results from statistical physics, we introduce temperature dependency to all rates, resulting in a remarkable correspondence between model and experiments. Our model predicts how temperature affects downstream protein production and find a crossover, where high affinity genes upregulates at high temperatures. Finally, we show how or that oscillatory temperatures can entrain NF-κB oscillations and lead to chaotic dynamics presenting a simple path to chaotic conditions in cellular biology.
1 Introduction
The fine-tuned regulation of protein production is fundamental to all living organisms. This production is complex and includes a number of components, but central to the stimulation of genes is the concentration of transcription factors inside the nucleus. In the past two decades, it has been revealed that the nuclear concentration of a number of central transcription factors, can be highly dynamic and it is expected that this dynamics might be an important element in the complex gene regulation of cells. One example of such a transcription factor, is the p53 tumour suppressor that has a period of ≈5 h (Lahav et al., 2004; Chen et al., 2016; Heltberg et al., 2019a) and another, which is the scope of the present work, is NF-κB, which oscillates with a period ≈1.5 h and has been shown to control the production of a number of proteins related to the immune response (Gonze et al., 2002; Hoffmann et al., 2002; Nelson et al., 2004; Krishna et al., 2006; Levine et al., 2013; Kellogg and Tay, 2015). It is at present debated what the functional role (if any) of these oscillations is, but it seems certain that the downstream genes are affected by this dynamics (Hoffmann et al., 2002; Nelson et al., 2004; Mengel et al., 2010; Tay et al., 2010; Heltberg et al., 2021).
The NF-κB signaling pathway is one of the most essential signaling pathways in eukaryotic cells. Among other functions, it has a role in cancer, inflammation, ageing, and in the immune defence, and moreover, it is serving as a general stress response (Wang and Tarbell, 1995; Mercurio and Manning, 1999; Hoffmann et al., 2002; Nelson et al., 2004; Courtois and Gilmore, 2006; O’Dea and Hoffmann, 2009; Sugioka et al., 2010; Kellogg and Tay, 2015; Liu et al., 2015; Mitchell et al., 2016; Heltberg et al., 2021). Since the regulation of temperature is also predicted to be a fundamental part of especially the immune response, we were interested in studying the interplay between oscillatory NF-κB and variations in temperature. A heat shock protein-dependent mechanism has been proposed as a link between the NF-κB signaling pathway and heat (Liu et al., 2015), and recently a model introducing a delay of the A20 signaling protein has been suggested (Harper et al., 2018), but we wanted to test this further, and establish the results for different levels of TNFα induction. Furthermore, little is known about the mechanisms of how the affected NF-κB oscillations could affect downstream protein production which we also set out to investigate in this study.
The theory of how reaction rates were affected by temperature, was pioneered by the work of the polish physicist Marian Smoluchowski (Smoluchowski, 1916), who calculated the reaction rate between two spherical particles diffusing in a potential. Even though this theory is more than 100 years old, and is a part of every physics curriculum, it has rarely been applied to systems that show dynamical behaviour such as limit cycles.
In this paper, in parallel to the experimental investigation, we apply the theory of temperature dependent reactions, to predict how the oscillatory dynamics of the transcription factor NF-κB is affected. In the experiments, we use both single and double additions of the ligand TNF-α as well as applying this through a custom made flow chamber thereby achieving a constant concentration in TNF-α, as opposed to a declining concentration as is the case if only added once [as done, e.g., in Ref. (Harper et al., 2018)] The latter condition enables us to detect clear oscillations in the nuclear concentration of TNF-α, and we succeed in varying the temperature from 32°C to 41.5°C in the flow chamber, where oscillations are still maintained. Through the analysis of these data we find a clear dependency of the period on the external temperature, so increasing the temperature leads to faster oscillations. By using a well-tested model of the dynamics of NF-κB, we find a striking correspondence between the model and the experimental findings and our model additionally predicts that lowering the temperature can lower the threshold level of the Hopf Bifurcation which defines the onset of oscillations in NF-κB. Furthermore we use the model to predict how the change in temperature might lead to a significantly different downstream production level. Finally, we simulate how an oscillatory temperature can affect the NF-κB dynamics and we find, that this can lead to entrainment even for small temperature oscillations. Interestingly, for amplitudes of
2 Results
2.1 Experimental Set-Up and Quantification of Nuclear NF-κB Concentration Dynamics at Different Temperatures
We investigate how temperature affects oscillations of NF-κB in mouse embryonic fibroblast (MEF) cells, through three experiments:
• TNF-α was added once to the MEF cells,
• TNF-α was added twice with 40 min in between
• MEF cells were exposed to a constant concentration of TNF-α in a flow chamber.
Each type of experiment was conducted at 32.0°C, 34.5°C, 37.0°C, and 39.5°C while a temperature of 41.5°C was additinally used in the flow experiments. The sample was imaged in a fluorescence time-lapse microscope, where DsRed-labeled p65 would give information about the relative nuclear to cytoplasmic concentration of NF-κB in the MEF cells.
In the single addition experiments, TNF-α was added to the MEF cells 2 min prior to initiation of the experiments and the concentration was increased from 0 ng TNF-α/ml to 10 ng TNF-α/ml. At t = 0 min the fluorescence time-lapse microscope would start imaging at 10–12 different locations every 10 min. Experiments were terminated after 20–48 h. In the double addition experiments, TNF-α was added to the MEF cells increasing the concentration from 0 ng TNF-α/ml to 10 ng TNF-α/ml 2 min prior to initiation of the experiments, and then after ≈38 min the TNF-α concentration was increased to 18 ng TNF-α/ml. These experiments were also terminated after 20–48 h. The reason we systematically tested different doses in the single- and double addition experiments, was to map out the dependency (if any) on TNF-α addition, and by testing it though different protocols, also through the constant flow experiments, is to explore the robustness and origin of the oscillations. In the flow experiments, the TNF-α concentration converged towards a value of 1 ng/ml. These experiments were terminated between 40 h and 7 days after initialization (See Methods for details).
The flow set-up is presented schematically in Figure 1A and a more detailed illustration is shown in Figure 5 in the methods section. For these experiments, all tubing in the system would be filled up with TNF-α-free medium and the syringe connected to the inlet would thereafter inject medium with TNF-α, meaning that the TNF-α level in the flow chamber would increase. It is known from earlier experiments that the degradation of TNF-α is slow, but it is not clear whetherthe declining concentration seen both in the double, but in particular in the single-addition experiment, would affect the results. By comparing to the flow experiments, where the concentration of TNF-α was surely constant, we could control any effect of gradually decreasing concentrations of TNF-α.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/910738/fcell-10-910738-HTML-r1/image_m/fcell-10-910738-g001.jpg)
FIGURE 1. Experiments on NF-κB oscillations conducted on mouse embryonic fibroblast cells. (A) Schematic drawing of the flow experiment with a flow chamber containing NF-κB oscillating cells. Through a syringe, a medium containing TNF-α is injected into the flow chamber, which initiates NF-κB oscillations. (B) Representative example of experiment where NF-κB is primarily located in the cytoplasms. (C) Same as B, but with NF-κB is primarily located in the nuclei. (D) The top row images show a time-lapse of a single fibroblast cell with p65 fluorescently labeled DsRed. At temperature 32.0°C the cell is exposed to 10 ng TNF-α/ml at t = −2 min and a new concentration increment, to a total of 18 ng TNF-α/ml, at t = 38 min. Each image corresponds to the time points in the plot below. The graph below, shows IN/IC ratio of the first 200 min of the data series. Data is extracted in 10 min intervals, with image examples every 30 or 40 min. (E) Same as D, but for temperature at 39.5°C. (F) Power spectrum of the IN/IC ratio of the data presented in (D). (G) Power spectrum of the IN/IC ratio of the data presented in (E). (H) Ratio (IN/IC) versus time, measured after adding 10 ng TNF-α/ml to μ-wells at t = −2 min. (I) Ratio IN/IC versus time where 10 ng TNF-α/ml was added at t = −2 min and the concentration was again increased at t = 38 min to a total of 18 ng TNF-α/ml. (J) Ratio IN/IC for the flow experiments where the TNF-α concentration is converging towards 1 ng/ml during the entire range plotted.
In Figure 1B an example of a fluorescent image of an MEF cell at 37°C is presented at t = 0 min, where TNF-α was added to a concentration of 10 ng/ml at t = −2 min. The emitted light is from the DsRed-p65 complex and hence the image taken at t = 0 is a signature of the NF-κB accumulation in the cytoplasm in all the cells. Additionally, we also see from Figure 1B that all nuclei are depleted from NF-κB, thus showing near complete translocation of NF-κB from the nucleus to the cytoplasm. This is also shown in Figure 1C, but at time t = 30 min. Here the cytoplasm has only a fraction of the NF-κB compared to Figure 1B, and instead the NF-κB has translocated into the nucleus.
To analyse the dynamics of the nuclear NF-κB concentration, two different methods were used. One by doing statistics on a large number of cells (see Period-extraction method in Methods) and one for the visualization of the oscillations (see Trace-vizualisation method). With these two methods combined, we made sure that the oscillations qualitatively existed, by applying the Trace-visualization, and extracted the correct frequency of the oscillations by using the Period-extraction method. Based on the images, we were able to quantify the relative nuclear concentration of NF-κB and plot this as a time series, where we at this temperature (32°C), observe approximately two periods in this time interval (Figure 1D, below). Here, the image series shows how the cytoplasm is bright and the nucleus is dark at t = 0 meaning that NF-κB is accumulated in the cytoplasm (Figure 1D, above). By inspection for a higher temperature (39.5°C), we find approximately three periods, which can be seen both in the images and in the plot of IN/IC ratio (Figure 1E).
We quantified the frequency of these periods, by calculating the power spectrum (see Methods), and by applying this to the time series above, we could extract frequency of 8.26 ⋅ 10−3min−1 for the time series at 32°C, which corresponds to a period of 121 min (Figure 1F). Similarly, we could do the same for the time series at 39.5°C, finding a frequency of 13.2 ⋅ 10–3 min−1 corresponding to a period of 74.1 min (Figure 1G). We further visualized the data (see Methods) of a single addition, a double addition and a flow experiment at 37°C, by showing the nuclear to cytoplasm intensity ratio, IN/IC, plotted versus time (Figures 1H–J). We note that when the TNF-α concentration is increased only once (Figure 1H), a transient peak in the NF-κB concentration appears, followed by periodic oscillations with lower amplitude (Figure 1H). However, when TNF-α is added at t = −2 min and again at t = 38 min, the transient peak is followed by another peak that has higher concentration than the average amplitudes, which again is followed by oscillations with lower, and slightly decaying amplitudes (Figure 1I). Based on this, we note that when the TNF-α concentration is abruptly increased, the following NF-κB peak has an larger amplitude, which is the case for both addition of one and two subsequent doses of TNF-α. These initial peaks after TNF-α addition is expected (Kellogg and Tay, 2015; Zambrano et al., 2016), but here we reveal that two successive additions will create two initial peaks with higher amplitude.
Our goal is to identify the periods of the oscillations in order to measure the difference between the three types of experiments and to finally identify how the temperature affected the period of NF-κB oscillations.
We first compare the oscillations at 37.0°C, by calculating the power spectrum for each of the three types of experiments (see Methods). Here we find the periods to be.
• 105.3 min ± 21.3 min (for the single addition with n37single = 12)
• 114.7 ± 29.7 min (for the double addition with n37double = 10)
• 100.8 ± 10.9 min (for the flow experiments with n37flow = 10)
From this, it is concluded that there is no significant difference across the three types of experiments at this physiologically relevant temperature, with Student’s t-test resulting in p-values
Next, we applied this method to all experiments at different temperatures. Here we measured the dynamics for temperatures down to 32°C for all three experiments. The maximum temperatures where NF-κB oscillations were possible to extract, were 39.5°C for single- and double addition experiments and 41.5°C for flow experiments, indicating that the flow setup might be more robust in order to measure the dynamics during external stresses. By analysis of these time series, we find that the NF-κB oscillation period decreases as a function of increasing temperature (Figure 2A). This relation was found in all three experimental conditions, emphasizing the generality if this result.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/910738/fcell-10-910738-HTML-r1/image_m/fcell-10-910738-g002.jpg)
FIGURE 2. Experimental observations (A) and model data (B–F) of the temperature dependency in NF-κB frequency. (A) Oscillation period of NF-κB in fibroblast cells vs. temperature for the three different types of experiments: Single addition experiment, where 10 ng TNF-α/ml was added at t = −2 min, double addition experiment, where 10 ng TNF-α/ml was added at t = −2 min and followed by another increase of 8 ng TNF-α/ml at t = 38 min, and lastly, flow experiments, where cells were exposed to flow in a flow chamber where the TNF-α concentration was converging towards 1 ng/ml. All types of experiments were performed at 32.0°C, 34.5°C,37.0°C, and 39.5°C, and the flow experiments were in addition performed at 41.5°C (a dot corresponds to a measurement for one cell). (B) Schematic figure, showing the components in the NF-κB network. (C) Time traces at three different temperatures. (D) Phase space for NF-κB, IRNA and IkB at three different temperatures. (E) Power spectrum of the oscillations. The rightt most peak corresponds to the observed frequency. (F) NF-κB oscillation period as a function of the external temperature shown for three different levels of TNF. Points in black corresponds to the experimental findings and with bars. (G) Absolute value of the derivative hog NF-κB traces, calculated as a mean over the time series, as a function of the applied TNF level. Note that these curves rise at the onset of oscillations and thereby locates the point of the Hopf Bifurcation.
The flow experiments did result in a more smooth variation than in the experiments with single and double addition of TNF-alpha. When grouping the assays as a weighted average the periods are 160.4 ± 32.4 min and 106.8 ± 22.0 min for the 32.0 and 37.0°C, respectively. In spite of the similar results across the different types of experiments, the flow experiments resulted in more stable results, and so, for the calculation of the change of period per change of temperature, the flow experiments were used.
From these experiments, the period is constant throughout different TNF-α concentrations and throughout different methods TNF-α exposure, however, varying the temperature +4.5°C or −5.0°C compared to body temperature significantly changes the period, and is found that Δperiod/Δtemperature = −11.9 ± 2.8 min/°C.
2.2 Inclusion of Temperature Dependency in Mathematical Model Predicts Changes in Periods and Hopf Bifurcation
To gain insight into the mechanism behind the temperature dependency of the oscillations, we tested whether these experimental findings could be explained by the classical theory of temperature dependency of reaction rates. To describe the dynamics of the NF-κB concentration, we used a mathematical model, that has been simplified in order to reduce the number of parameters, and that has previously been applied to model the dynamics of NF-κB (Jensen and Krishna, 2012; Heltberg et al., 2016; Heltberg et al., 2019b; Heltberg and Jensen, 2019). In this model, we consider the NF-κB inside the nucleus (Nn), acting as a transcription factor for many proteins, including Iκ-B. The equations take the following form:
Here, Nn is the nuclear NF-κB concentration, Im is the IkB mRNA level, and I is the concentration of cytoplasmic I-κB protein. Iκ-B kinase (IKK) is activated by the external signal TNF and causes eventual targeted degradation of Iκ-B when it is bound to NF-κB. IKK has the forms: neutral (n), active (a) and in-active (i). The network is schematized in Figures 2B, and a more detailed description can be found in the methos section, where all the parameters are listed as well (note that in the model section we use the notation TNF instead of TNF-α for simplicity).
We now include the temperature dependency into the simulations. In the model, we have 9 rates, and these were all made temperature dependent (See Table 1 in Methods). We assumed that the fastest ones: α, kNin and kt followed a Smoluchowski dependency (i.e., they are diffusion limited, k+ ∝ D(T)), whereas the others followed an Arrhenius dependency (i.e., they are reaction limited
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/910738/fcell-10-910738-HTML-r1/image_m/fcell-10-910738-t001.jpg)
TABLE 1. Default values of parameters in the model. Here we assume that the slowest rates are governed by the Arrhenius equation, whereas the fastest rates will follow the temperature dependency of the Smoluchowski rate (Diffusion limited). The temperature dependencies and their activation energies are listed when applicable.
With this set-up we were ready to simulate the dynamics of NF-κB at different temperatures. By varying the temperature we found that the NF-κB oscillations were highly affected by the change in the temperature level and that low temperatures led to large periods and amplitudes (Figure 2C). It should be noted that it is possible to change the period of oscillations, without affecting the amplitudes by simply rescaling the time dependent parameters. However, since the temperature dependency acts both on the diffusion limited and the reaction limited parameters, we do not obtain this scaling and as a result the amplitudes are affected by this temperature variation. We visualized this further in the three-dimensional phase space spanned by the NF-κB, IRNA and IkB. Here we note that the change in temperature affects the entire size of the limit cycle and for high temperatures the limit cycle shrinks, leading to faster oscillations (Figure 2D). To quantify these oscillations, we calculated the power spectrum, by applying the FFT algorithm, finding the leading frequency of the time series. We note that multiples of this oscillation appear as well, but taking the maximal value of the power spectrum we find the correct frequency of the oscillations (Figure 2E). We now used this algorithm to calculate the period as a function of the applied temperature. Here we found that the curve decreases monotonically for increasing temperature, which is what we would expect by observing the time series, and by comparing these results to the experimental observations we find a striking compliance (Figure 2F). We tested this for different values of TNF baseline levels, and here we found similar patterns indicating that this result is quite robust and not sensitive to our initial choice of external TNF level (Figure 2F). We also conclude that this scaling with temperature is not dependent on the followed protocol, which solidifies the robustness of this study and can be interpreted as TNF-α is extremely slowly degraded in these environments. Since the temperature affected the period of the oscillations, we were also interested if this meant that changes in the temperature could also induce oscillations and thereby affect the point of the onset of oscillations (i.e. the Hopf bifurcation). Here we found to our surprise that decreasing the temperature, would lead to oscillations for smaller values of the external TNF levels (Figure 2G). This means that it would be possible to use regulations of temperature to induce or stop the oscillations in living organisms.
2.3 Control of Downstream Protein Production by Changing the Temperature
As we have established how the dynamics of NF-κB could vary with different temperature levels, we wanted to investigate how this could affect the downstream production of proteins stimulated by NF-κB. Here we used a previously suggested model (Heltberg et al., 2016), where all genes are divided into groups based on their affinity and cooperativity from stimulation with NF-κB. We assume that NF-κB can bind to an enhancer or operator region, and can form complexes to bind the RNA polymerase, with different affinity, depending on the gene (schematically shown in Figure 3A). We describe the transcription and translation of each gene, labelled i = 1, 2, 3, … , using the differential equations [see (Heltberg et al., 2019b)]:
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/910738/fcell-10-910738-HTML-r1/image_m/fcell-10-910738-g003.jpg)
FIGURE 3. Temperature dependency of the protein production in downstream genes from NF − κB oscillations. (A) Schematic figure showing the stimulation of a high affinity gene (left) and a low affinity gene (right). (B) Protein production of Protein 1 (from a high affinity gene) at three different temperatures. (C) Same as (B) but for Protein 2 (from a low affinity gene). (D) Relative steady state production of Protein 1 and Protein 2 respectively as a function of temperature.
Here, the mi represent the mRNA level transcribed from gene i, and Pi represents the concentration of proteins produced from the corresponding mRNA. The first term in the equation for the mRNA is known as a Hill function; the canonical way to describe the protein production for genes governed by transcription factors where each gene has a specific Hill coefficient and effective affinity (Werner et al., 2007; Kaplan et al., 2008; Maienschein-Cline et al., 2010; Mengel et al., 2010; Sneppen et al., 2010).
The effective affinity Ki is a parameter that combines the strength of binding of the transcription factor to the operator/enhancer region, the strength of binding of RNA polymerase to the promoter and transcription factor, as well as the effect of DNA looping that may be needed to bring the enhancer/operator close to the promoter region. Operationally, Ki sets the concentration of NF-κB that results in 50% of maximal gene expression enhancement.
With this set-up, we simulated the model of NF-κB with varying the temperature as shown in Figure 2, and measuring the produced proteins. For simplicity we will only consider two proteins, Protein 1 and Protein 2, being stimulated from a High and Low affinity gene respectively. From the model we see that as we increase the temperature, the steady state level of Protein 1 is enhanced (Figure 3B). However as we assessed the steady state level of Protein 2, we realised that this was significantly reduced and thus these types of proteins would be up-regulated if the cell could lower the temperature (Figure 3C). The reason for this change in the protein level, is predominantly due to the variation in the amplitudes of NF-κB. The low affinity genes (here Protein 2) will rarely be expressed, unless very high levels of transcription factor available, which happens during the large transient pulses. However, due to the spiky oscillations in this regime, NF-κB remains very low for long periods of time which will cause a decrease in the high affinity genes that are otherwise always expressed. In this way the tuning of oscillations might alter the overall protein production. With this information we simulated the steady state protein level of both proteins as a function of temperature, and here we found a very interesting crossover effect (Figure 3D) which indicates that proteins from High affinity can be up-regulated as one increase the temperature whereas proteins from Low affinity genes are monotonically decreased for increasing temperature. This result highlights the fascinating prospect, that one can use temperature as a regulator for the downstream production of proteins.
2.4 Synchronization and Chaotic Dynamics Emerges From Temperature Oscillations
At this point, we had established that temperature variations affect the oscillatory properties of NF-κB. Therefore, we hypothesized that a periodically varying temperature could lead to highly complex dynamics in the nuclear NF-κB. Mathematically, the introduction of periodically varying temperature variations turns the system into two coupled oscillators (Jensen et al., 1983; Jensen et al., 1984; Heltberg et al., 2016; Heltberg et al., 2021), where NF-κB is an internal oscillator, stimulated by an external temperature oscillator. In the experiments, we found that cells would survive and remain oscillatory for temperature variations of ± 5°C, and in order to remain close to the experimental observations, we allow the amplitude of the temperature oscillations to be maximally 5°C.
To a start, we oscillated the temperature with an amplitude of 1°C and by varying the frequency we observed that different entrainment modes emerged. First we found the 1/1 coupling, which means that one full period of the temperature corresponds to one period of the NF-κB system, where the phases are locked, Figure 4A [see e.g., (Heltberg et al., 2016; Heltberg et al., 2021)]. By varying the period, it became clear that this temperature dependency could lead to entrainment for different rational numbers, for instance with a 5/3 coupling (Figure 4B) and strong 2/1 coupling (Figure 4C). A 5/3 coupling corresponding to a synchronized situation where the NF-κB signal performs 3 cycles while the temperature performs 5 cycles (see Figure 4B). Similarly, the 2/1 coupling corresponds to a state with one NF-κB cycle for each two temperature cycles (Figure 4C). Based on these first observations, we varied the frequency and measured the rotation number of the system. This rotation number we define as the (externally fixed) frequency of the temperature divided by the (measured) NF-κB frequency. By doing this we were able to extract the resulting plot known as a Devil’s staircase (Figures 4D,E) for two different values of the external amplitude. Here we find entrainment plateaus (horizontal regions), were the NF-κB frequency is entirely determined by the temperature frequency. This means that inside these regions, one can completely control of behaviour of the NF-κB oscillations and either speed up or slow down the dynamics. We note by comparing Figures 4D,E that the dominating plateaus grow in range, as we increase the amplitude from 0.5°C to 1°C. We therefore wanted to measure the width of these entrainment regions especially for the dominating ones. First we tested how this was affected by the external level of TNF, and we found that even though small variations occurred, these entrainment regions were quite stable and robust to changes in the TNF level (Figure 4F). This is a promising observation, since it allows future experiments in vivo and in vitro to focus on the changes in temperature, without worrying about the small differences in the levels of TNF. Next we increased the temperature amplitudes, and here we found that all the dominating regions were growing, whereas the smaller rational regions (such as 5/3) loose their stability as they are being “squeezed” out by the dominating ones (Figure 4G).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/910738/fcell-10-910738-HTML-r1/image_m/fcell-10-910738-g004.jpg)
FIGURE 4. Temperature oscillations lead to entrainment and chaotic dynamics for NF-κB. (A) Oscillations in NF-κB (blue, left axis) coupled to temperature (red, right axis). Temperature amplitude = 1°C and frequency (ftmp) = 0.01[ min−1]. (B) Same as A, but with temperature frequency = 0.018[ min−1] resulting in a 5/3 coupling. (C) Same as A, but with temperature frequency = 0.02[ min−1] resulting in a 2/1 coupling. (D) Rotation number measured as ftmp/fNF−κB as a function of the oscillation frequency for the temperature. Amplitude of temperature = 0.5°C. (E) Same as D, but with amplitude of temperature = 1°C. (F) Width of the dominating steps in the figures D + E as a function of the external TNF level. Amplitude of the temperature oscillations = 1°C. (G) Width of the dominating steps in the figures D + E as a function of the amplitude of temperature oscillations. TNF level = 0.2. (H) Dynamics of NF-κB for different initial conditions (indicated by colours) as a function of time (above) and in the phase space spanned as IRNAκB and IκB (I) Same as H, but for two initial conditions separated by 10–4%. (J) Mean distance of the trajectories, initially separated by 10–4% as a function of the amplitude of the temperature oscillations. Curves shown for three different frequencies of temperature oscillations.
It has previously been observed that for large amplitudes of oscillatory TNF, the dynamics of NF-κB, could show “modehopping” which corresponds to transitions between two stable limit cycles (Heltberg et al., 2016; Heltberg and Jensen, 2019). We were interested if this was also the case, if we applied temperature oscillations, and we realised that already for amplitudes of 2°C, multistability occured, and we were even able to find three stable limit cycles by varying the initial conditions of the simulation (Figure 4H). We note that for some of these limit cycles, the dynamics could entrain in different phases even though they were part of the same attractor in the phase space (see IC2+IC4 in the time series above in Figure 4H). This observation also indicated that the applied amplitude was above the “critical value”, and due to a theorem of Poincare, above this line various complex dynamics could emerge and in particular chaotic dynamics (Jensen et al., 1984). Therefore, we investigated this by increasing the external amplitude, and found that chaotic dynamics could emerge for amplitudes around 2.5°C. We visualized the chaotic dynamics, by simulating the system with identical parameter values, and initial conditions only separated by 10–4%. By studying the time series, it is clear that even though they show completely similar dynamics for a long time, then the time series evolve completely differently (Figure 4I above). We also visualized this in the three-dimensional phase space, and could see that the trajectories moved on a strange attractor (Figure 4I below). Finally we were interested in studying for which values of the temperature oscillations, one expect to see the chaotic transition. To quantify this, we measured the average distance between two trajectories only separated by 10–4% in the initial conditions. Using this measure we observed that for relatively large frequencies of the temperature oscillations (periods of 30–50 min), we found a chaotic transition for
3 Discussion
Temperature is known to vary in living organisms as a response to external stresses, but a mechanistic understanding of how this can affect the protein production and control transcription factor dynamics has been scarcely investigated. In this work, we have shown how the transcription factor NF-κB can change it is properties of oscillations by varying the temperature in a realistic range where the cells can still survive and maintain oscillations. We use theories of statistical physics and in particular the Smoluchowski equation to predict how the individual rates are affected by a changing temperature, and show that these are sufficient to describe the results of the experimental data. We further use this calibrated model, to predict how different families of downstream genes are affected by a change in temperature and we point out a crossover effect with strong regulatory implications. Finally, we show that temperature can be used as a strong source to induce more complex dynamics to the NF-κB and we suggest that this can be used as a key measure to investigate the role of dynamics in transcription factors on the downstream gene production.
The role of temperature dependency on the NF-κB oscillations has previously been studied (Harper et al., 2018), where an increasing frequency for increasing temperature was also observed, even though this effect was only investigated for a smaller temperature range than investigated here. One of the findings in (Harper et al., 2018), was an A20 dependent mechanism that was a key element in the perturbation of NF-κB frequency. Here, we apply a model where A20 is assumed to be constant, which is originally based on the findings of Ashall et al. (2009), in which they argue that there is a range of constitutive A20 expression values that can functionally replace A20 negative feedback. Therefore, the mechanism of the temperature dependency may not be linked directly to the dynamics of A20, however, two of the most impactful parameters (ki and kp, see Table 1) both affect the part of the network that A20 also acts on. Therefore, our results are generally in agreeing with those of Harper et al., however, they suggest that on a fundamental level, it is the effective turnover of IKK that mediates the temperature dependency in the oscillations of NF-κB.
Our findings reveal a fundamental temperature dependency in the oscillations of NF-κB and it is therefore intriguing to speculate what the functional role of this might be. We hypothesize that the alterations in the oscillations might affect and stimulate different groups of genes. However, the most studied oscillator in biological systems, the circadian clock, is known to be robust to changes in temperature (Leloup and Goldbeter, 1997; Hong et al., 2007), and therefore it could be argued that it would be beneficial to biological oscillators to adopt this robustness. We believe that a key difference between these biological oscillators is that the NF-κB oscillations should be regarded as a response and therefore it would be advantageous to be able to tune this response to different stimuli and the temperature in particular. On the other hand, the circadian oscillator is basically a very advanced clockwork, where robustness is the key to maintain time keeping. Finally, in the tissues where the NF-κB oscillator is important, the temperature is well-regulated and thereby changes in temperature would be a parameter to trigger a response where the NF-κB oscillator would be part of it.
The finding that oscillations emerge by lowering the temperature while keeping TNF at a constant low level, combined with the fact that the low affinity genes are in general enhanced in the low-temperature regime, suggests a potential functional role of the temperature stimulation on NF-κB. Since NF-κB stimulates several hundreds of genes downstream, it is potentially important to stimulate a subset of genes, while lowering other families of genes. The mechanism suggested here can create this type of cross-over effect by simply changing the temperature for the cells. In this way an entire cascade can potentially be initiated by lowering the temperature, which will lead to an up-regulation of specific groups of genes that are otherwise not expressed.
Even though chaotic dynamics has been studied theoretically and experimentally for more than 50 years, it has still not been shown to exist and play a role in cellular dynamics. Complex phenomena such as synchronization has been shown to exist (Danino et al., 2010; Kellogg and Tay, 2015; Heltberg et al., 2016), and even modehopping that reveals the existence of multistable cycles, when the amplitude of the external oscillator has been sufficiently increased. Theoretically this should also guide a way to induce chaotic dynamics (Jensen et al., 1983; Jensen et al., 1984; Heltberg and Jensen, 2019), however often cells have trouble surviving the transiently high concentrations of TNF-α, and therefore it has been difficult to predict how chaotic dynamics might be investigated for such systems. Our results predict that temperature oscillations is a simple and effective way to induce highly complex dynamics when the temperature is oscillated externally. Since our experimental results reveal that cells can survive and maintain oscillations under temperature variations of
Temperature presents a fundamental, physical property with the potential to control and regulate the dynamical properties of protein concentration in cells. It is our hope that this work will inspire theoretical and experimental explorations these prospects of transcription factors, as the system is affected by either constant or dynamically varying temperatures in living cells.
4 Methods
4.1 Cell Culturing
All cells used in this article were stably transduced 3T3 mouse fibroblast cells obtained from Tay et al. (2010). Using lentivirus, DsRed was introduced into the genome of the cells, resulting in the expression of the fusion protein DsRed-p65. This allowed for tracking of p65 when shuttling between cytoplasm and nucleus of the cells. These cells also had the nuclear marker H2B-GFP which was not utilized in the presented data of this article.
Cells were incubated in 5.0% CO2 at 37.0°C in a HERA CELL VIOS 160i CO2 incubator. For culturing Gibco™ DMEM culture medium with high glucose, L-Glutamine, phenol red, and no Sodium Pyruvate and no HEPES were used. While performing experiments this medium was substituted with Gibco™ DMEM culture medium with high glucose, HEPES, L-Glutamine, no phenol red, and no Sodium Pyruvate. To both media +10% Fetal bovine serum (FBS) and +1% Penicillin-Streptomycin (PS) was added.
4.2 Single Well Experiments
48 h prior to the experiments fibroblast cells were seeded on collagen-coated, γ-radiated, 35 mm, No. 1.5 glass bottom microwell dishes (Mat Tek, 2021). The microwells were brought from the 37.0°C, 5% CO2, high humidity incubator into the incubation chamber of the microscope with an atmosphere of 5% CO2, high humidity, and a specific target temperature. Approximately 15 min after, when the temperature had stabilized, the 10 ng TNF-α/ml was added to a total of 40 ng TNF-α. One to two minutes thereafter the fluorescence microscope time-lapse was initiated. This initiation was defined as time = 0 min in the experiments. Timelapse images were obtained every 10 min where 10–12 positions were captured at each time point throughout the entire timelapse series. In the experiments where TNF-α was added twice, the first addition was once again 10 ng TNF-α/ml to a total of 40 ng TNF-α, and the second addition took place one to two minutes before the fourth set of timelapse images at t = ≈ 38 min, and consist of a total increase of 8 ng TNF-α/ml. At the second addition 50 ng TNF-α was added and the total TNF-α added to these experiments were 90 ng TNF-α.
4.3 Flow Experiments
When performing flow experiments illustrated in Figure 5, a major task is to prevent gas formation in the perfusion chamber where cells grow. The media utilized in the experiments have to be stored at room temperature but the temperature of the perfusion chamber is approximately 10–20°C higher than room temperature, and since Gas is extruded from a liquid when it is heated up, gas formation in the perfusion chamber is an inherent problem to flow chamber (MBR) flow experiments. If the perfusion chamber is filled up with gas, even temporarily, in the order of minutes, cells will be injured. They will not be able to fully recover and likely cells will undergo necrosis. If cells survive, however, our pilot experiments showed that nuclear translocation of NF-κB shuttling was prevented. Following, the extensive measures required to prevent gas formation are described.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/910738/fcell-10-910738-HTML-r1/image_m/fcell-10-910738-g005.jpg)
FIGURE 5. Schematic drawing of the flow system with numbers added, indicating how the medium flows through the system, chronologically. The solid line connecting the devices is a gas permeable silicone tubing. 1) Pump injecting medium into the system at a rate controlled by a computer via Labview software. 2) Medium flows into the microscope enclosure. Medium changes from room temperature to target temperature −2.0°C and changes from the starting height to the same level as all devices inside the microscope enclosure (ΔHeight ∼ −40 cm). 3) 1.0 m of extra tubing lets the medium reach target temperature −2.0°C before entering the bubble trap. 4) Two bubble traps capturing excess gas in the tubing. 5) Medium enters incubation chamber and the atmosphere changes from room atmospheric CO2 and room atmospheric humidity to 5% CO2 and high humidity while also the temperature increases with 2.0°C to the target temperature. 6) Extra tubing allows the medium to adjust to CO2 and temperature. 7) Perfusion chamber where the fibroblast cells grow and the microscopy data is collected. 8) Medium exits incubation chamber and CO2 and humidity changes to atmospheric levels. Temperature changes to target temperature −2.0°C. 9) Medium exits microscope enclosure and temperature changes from target temperature −2.0°C to room temperature. Height of medium is increased (ΔHeight ∼ +30 cm). 10) Medium flows out of the flow system and into a glass beaker.
Prior to the initialization of the experiments, the medium was degassed for 90 min. Subsequently, the medium was brought to a flow hood and filtered using a 0.2 μm filter to remove biological contaminants from the degassing process. The medium was separated into two portions. One where TNF-α was added and one that was left for filling up the flow system before starting the experiment. From degassing, the medium was cooled down, so before injecting it into the flow system5, consisting of tubing, bubble traps, and perfusion chamber, the medium was heated up to room temperature in a sealed falcon tube to prevent new gas from being dissolved in the medium while the medium was heating up. The entire flow system was then filled up with medium in steps and assembled inside the flow hood. In each bubble trap, 3 ml of medium was injected, leaving 2 ml in each bubble trap for air. The entire flow system was then brought to the microscope set-up and inserted into the microscope enclosure. The chamber and 1 m of tubing were inserted into the incubation chamber, an inlet was connected to the controllable pump outside the enclosure and an outlet was inserted into a glass beaker.
The inlet and outlet of the system were elevated compared to the other devices in the system. This was to create increased pressure in the perfusion chamber and tubing situated in the incubation chamber in order to minimize gas formation. In between the injection and the perfusion chamber, the temperature increases in two steps. The first increment is close incubation chamber, where the temperature is reached, which serves two purposes: 1) It stabilizes the temperature surrounding the incubation chamber, which will minimize temperature fluctuations. 2) When gas forms in the medium due to increased temperature, the temperature has already increased significantly before entering the bubble trap and most gas will be trapped herein. In between the microscope enclosure and the bubble trap, extra tubing is inserted to make sure the temperature of the medium is increased and to give the medium time to extrude gas before entering the bubble trap. After the bubble trap, the medium flows into the incubation chamber, where the medium and perfusion chamber is exposed to the target temperature, 5% CO2, and high humidity.
4.4 Microscopy
A Nikon inverted fluorescence microscope was utilized throughout these experiments (Nikon, 2021). A Lumencor Sola Light Engine (Lumencor, 2021a) solid-state illumination provides fluorescent light combined with excitation- and emission filters in the microscope at 575 and 641 nm, respectively. Bright-field images are provided by the build-in lamp in the microscope. Images are captured with an Andor Neo sCMOS model DC-152Q-COO-FI camera (Lumencor, 2021b).
4.5 Temperature Control
Temperature control inside the microscope enclosure was provided by Oko lab incubation unit model H201-T-0016 (Oko Cage Incubator, 2021). The temperature of the incubation chamber, provided by the Warner Duel Automatic Temperature Controller TC-344B (Warner, 2021), was 2.0°C higher than the temperature in the microscope enclosure. Each of its two controllers had two heaters, a control thermometer, and a monitor thermometer. The two heaters connected to one of the controllers were mounted to the bottom stage of the outsides of the incubation chamber and two heaters connected to the other controller were mounted to the inside of the incubation chamber. The heaters outside the incubation chamber were constantly exposed to 6 V from the controller, while the heaters inside were continuously controlled by a loop, meaning the temperature inside the incubation chamber could be set to a constant value. The control thermometer of the inside heaters was taped on top of the flow slide in the case of flow experiments and taped onto the bottom of the inside of the incubation chamber in the case of micro dish experiments. During flow experiments, the temperature of the cells inside the perfusion chamber was assumed to be equal to the temperature on top of the flow slide, where the control thermometer was attached. In micro dish calibration experiment it was shown that the temperature inside a micro dish containing 4 ml of water was 1.0°C lower than a thermometer attached to the bottom of the inside of the incubation chamber, so in these experiments, the temperature of the Warner temperature controller was set to 1.0°C higher than the target temperature.
4.6 Data Analysis
A semi-automatic Matlab program was built for the analysis, where nuclear light intensity (IN) and cytoplasm intensity (IC) were measured and an oscillation period was extracted for a single cell at the time, by taking the ratio of IN/IC. Prior to being analysed in the Matlab program, cells with clear oscillation patterns were chosen from the visual image analysis in Fiji. In the Matlab program, the average light intensity of a circular area for each of the nucleus and cytoplasm was selected for each time frame, corresponding to a change in time of 10 min. Depending on the size of the cell, the measured area would be between 3.1 and 8.7 μm2, corresponding to a circle with a diameter between 2 μ m and 3.3 μm, and corresponding to between 27 pixels and 79 pixels for each measured nuclear or cytoplasm intensity.
This approach was used throughout all data analysis, but it was found to be an advantage to both have a Trace-vizualisation method and the Period-extraction method for the total data analysis. The Trace-vizualisation method would visualize the oscillations to give an understanding of the shape of the oscillations as well as give an understanding of how well the system behaved. This method was very time-consuming, so another method, namely the and Period-extraction method was used to extract statistics on oscillations of a high number of cells in a more time-efficient manner.
4.7 Trace-Vizualisation Method
With the analysis in Fiji described above, a single area from the cytoplasm and a single area from the nucleus were extracted. Each cell was analysed in a time span of between 90 and 770 min, equivalent to be between 2 and 11 oscillations.
The Matlab program would then do a power spectral analysis of IN/IC—mean(IN/IC) and extract the frequency value of the highest peak, see Figures 1D–G). The inverse of this frequency value would then correspond to the most significant period in the data series. However, in some cases, a peak that did not correspond to the cell’s oscillations was the highest peak. In these cases, the high peak near the visually observed oscillation was chosen as the frequency for extracting a period.
4.8 Period-Extraction Method
In this analysis, four areas of the same size as described above were chosen for both the cytoplasm and for the nucleus. This made the total areas for the calculations of the IN/IC ratio larger resulting in local intensity variations being minimized. The IN/IC ratio was then smoothed with the Matlab Smooth Function, with a span = 4. From experimental errors, the IN/IC data had intensity variations larger than the typical oscillation period. These variations were filtered out by taking the Matlab Smooth Function with span = 13. The Matlab Smooth function with span = 13 was then subtracted to the Matlab Smooth Function with span = 4, which resulted in a data series well suited for qualitative understanding of the oscillations.
4.9 Temperature Dependency
With the mathematical input based on Smoluchowski equations (see Supplementary Materials) we define how we expect the rates should change according to a change in temperature. As standard we set
Now differentiating with respect to temperature gives:
Assuming linearity around the value of T, we can rewrite the differential equation into a difference equation:
The constant A, is known as the exponential prefactor in the Arrhenius equation, and this is a free parameter that we do not have knowledge about. In the equations, we obtain the scaling factor
Next consider the case of diffusion limited reactions. This is the other limit in the equations above, where we also assume there is no potential. Here we obtain the maximal rate as:
If we apply the Einstein-Stokes relation we find the temperature dependency in D0
This means that the rate would simply follow the temperature according to:
4.10 The NF-κB System
In this model, we consider the NF-κB inside the nucleus (Nn), acting as a transcription factor for many proteins, including I-κB. The equations are repeated here with the addition of an oscillating TNF value:
Here, Nn is the nuclear NF-κB concentration, Im is the IkB mRNA level, and I is the concentration of cytoplasmic I-κB protein.
All the parameters used in the NF-κB model are found in the table below.
We outline briefly the biological correspondences of the different terms in the model:
• In the equation for
• The equation for IRNA describes the NF-κB activated transcription of IκB mRNA and the spontaneous degradation of the mRNA with a half-life of ln(2)/γm.
• The first term in the equation for IκB models translation of IκB mRNA into IκB protein in the cytoplasm, and the second term models the TNF-triggered degradation of IκB in the cytoplasm when it is bound to NF-κB.
• The triggering stimulus TNF, acts by changing the level of active IκB kinase, [IKKa], which phosphorylates IκB, resulting eventually in its degradation. This degradation rate is set by the parameter α in the model. It is thus only this protein complex with IKK that can phosphorylate the NF-κB - I-κB complex and make NF-κB active again.
This model assumes that there is a constant amount of IKK (IKKtot), which can be in three states: active (IKKa), inactive(IKKi) and neutral (IKKtot − IKKa − IKKi). TNF increases the rate at which neutral IKK is made active, and decreases the rate at which inactive IKK is made neutral.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://github.com/Mathiasheltberg/TemperaturePrpject. Note that the original videos are too large to be placed in this folder and are available upon request to heltberg@nbi.ku.dk.
Author Contributions
MH performed the mathematical analysis and performed the simulations. MvB performed the experiments and the data analysis. All authors wrote and reviewed the paper.
Funding
We acknowledge support from the Independent Research Fund Denmark (grant number 9040-00116B), the Novo Nordisk Foundation (grant number NNF20OC0064978) and the Danish National Research Foundation through StemPhys Center of Excellence (grant number DNRF116). We acknowledge financial support from the Novo Nordisk Foundation (NNF) (grant agreement number NNF20OC0061176). MH acknowledges MH acknowledges the Carlsberg Foundation grant CF20-0621 and the the Lundbeck Foundation grant R347-2020-2250. MHJ acknowledges support from the Independent Re- search Fund Denmark grant number 9040-00116B.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We are grateful to Markus Covert for providing the cells, to Sandeep Krishna for knowledge about the TNF model and to Axel Boisen for assistance with constructing the custom-made flow chamber.
References
Amon, A., and Lefranc, M. (2004). Topological Signature of Deterministic Chaos in Short Nonstationary Signals from an Optical Parametric Oscillator. Phys. Rev. Lett. 92 (9), 094101. doi:10.1103/PhysRevLett.92.094101
Ashall, L., Horton, C. A., Nelson, D. E., Paszek, P., Harper, C. V., Sillitoe, K., et al. (2009). Pulsatile Stimulation Determines Timing and Specificity of NF-κb-dependent Transcription. Science 324 (5924), 242–246. doi:10.1126/science.1164860
Chen, S.-h., Forrester, W., and Lahav, G. (2016). Schedule-dependent Interaction between Anticancer Treatments. Science 351 (6278), 1204–1208. doi:10.1126/science.aac5610
Courtois, G., and Gilmore, T. D. (2006). Mutations in the NF-Κb Signaling Pathway: Implications for Human Disease. Oncogene 25 (51), 6831–6843. doi:10.1038/sj.onc.1209939
Danino, T., Mondragón-Palomino, O., Tsimring, L., and Hasty, J. (2010). A Synchronized Quorum of Genetic Clocks. Nature 463 (7279), 326–330. doi:10.1038/nature08753
Gonze, D., Halloy, J., and Goldbeter, A. (2002). Robustness of Circadian Rhythms with Respect to Molecular Noise. Proc. Natl. Acad. Sci. U.S.A. 99, 673–678. doi:10.1073/pnas.022628299
Harper, C. V., Woodcock, D. J., Lam, C., Garcia-Albornoz, M., Adamson, A., Ashall, L., et al. (2018). Temperature Regulates NF-Κb Dynamics and Function through Timing of A20 Transcription. Proc. Natl. Acad. Sci. U. S. A. 115 (22), E5243–E5249. doi:10.1073/pnas.1803609115
Heltberg, M., Kellogg, R. A., Krishna, S., Tay, S., and Jensen, M. H. (2016). Noise Induces Hopping between NF-Κb Entrainment Modes. Cell. Syst. 3 (6), 532–539. doi:10.1016/j.cels.2016.11.014
Heltberg, M. L., Krishna, S., and Jensen, M. H. (2019). On Chaotic Dynamics in Transcription Factors and the Associated Effects in Differential Gene Regulation. Nat. Commun. 10 (1), 71–10. doi:10.1038/s41467-018-07932-1
Heltberg, M. L., Chen, S.-h., Jiménez, A., Jambhekar, A., Jensen, M. H., and Lahav, G. (2019). Inferring Leading Interactions in the p53/Mdm2/Mdmx Circuit through Live-Cell Imaging and Modeling. Cell. Syst. 9 (6), 548–558. doi:10.1016/j.cels.2019.10.010
Heltberg, M. L., and Jensen, M. H. (2019). Locked Body Clocks. Nat. Phys. 15 (10), 989–990. doi:10.1038/s41567-019-0617-2
Heltberg, M. L., Krishna, S., Kadanoff, L. P., and Jensen, M. H. (2021). A Tale of Two Rhythms: Locked Clocks and Chaos in Biology. Cell. Syst. 12 (4), 291–303. doi:10.1016/j.cels.2021.03.003
Hoffmann, A., Levchenko, A., Scott, M. L., and Baltimore, D. (2002). The IκB-NF-Κb Signaling Module: Temporal Control and Selective Gene Activation. Science 298, 1241–1245. doi:10.1126/science.1071914
Hong, C. I., Conrad, E. D., and Tyson, J. J. (2007). A Proposal for Robust Temperature Compensation of Circadian Rhythms. Proc. Natl. Acad. Sci. U.S.A. 104 (4), 1195–1200. doi:10.1073/pnas.0601378104
Jensen, M. H., Bak, P., and Bohr, T. (1983). Complete Devil's Staircase, Fractal Dimension, and Universality of Mode- Locking Structure in the Circle Map. Phys. Rev. Lett. 50, 1637–1639. doi:10.1103/physrevlett.50.1637
Jensen, M. H., Bak, P., and Bohr, T. (1984). Transition to Chaos by Interaction of Resonances in Dissipative Systems. I. Circle Maps. Phys. Rev. A 30 (4), 1960–1969. doi:10.1103/physreva.30.1960
Jensen, M. H., and Krishna, S. (2012). Inducing Phase-Locking and Chaos in Cellular Oscillators by Modulating the Driving Stimuli. FEBS Lett. 586, 1664–1668. doi:10.1016/j.febslet.2012.04.044
Kaplan, S., Bren, A., Zaslaver, A., Dekel, E., and Alon, U. (2008). Diverse Two-Dimensional Input Functions Control Bacterial Sugar Genes. Mol. Cell. 29 (6), 786–792. doi:10.1016/j.molcel.2008.01.021
Kellogg, R. A., and Tay, S. (2015). Noise Facilitates Transcriptional Control under Dynamic Inputs. Cell. 160, 381–392. doi:10.1016/j.cell.2015.01.013
Krishna, S., Jensen, M. H., and Sneppen, K. (2006). Minimal Model of Spiky Oscillations in NF-Κb Signaling. Proc. Natl. Acad. Sci. U.S.A. 103, 10840–10845. doi:10.1073/pnas.0604085103
Lahav, G., Rosenfeld, N., Sigal, A., Geva-Zatorsky, N., Levine, A. J., Elowitz, M. B., et al. (2004). Dynamics of the P53-Mdm2 Feedback Loop in Individual Cells. Nat. Genet. 36 (2), 147–150. doi:10.1038/ng1293
Leloup, J.-C., and Goldbeter, A. (1997). Temperature Compensation of Circadian Rhythms: Control of the Period in a Model for Circadian Oscillations of the PER Protein in Drosophila. Chronobiology Int. 14 (5), 511–520. doi:10.3109/07420529709001472
Levine, J. H., Lin, Y., and Elowitz, M. B. (2013). Functional Roles of Pulsing in Genetic Circuits. Science 342, 1193–1200. doi:10.1126/science.1239999
Liu, Y., Zhou, G., Wang, Z., Guo, X., Xu, Q., Huang, Q., et al. (2015). NF-κB Signaling Is Essential for Resistance to Heat Stress-Induced Early Stage Apoptosis in Human Umbilical Vein Endothelial Cells. Sci. Rep. 5 (1), 13547. doi:10.1038/srep13547
Lumencor (2021b). Andor Neo sCMOS. Available at: https://andor.oxinst.com/products/scmos-camera-series/neo-5-5-scmos (Accessed: 01 14, 2021).
Lumencor (2021a). Lumencor Sola Light Engine. Available at: https://lumencor.com/products/sola-light-engine/(Accessed: 01 14, 2021).
Maienschein-Cline, M., Warmflash, A., and Dinner, A. R. (2010). Defining Cooperativity in Gene Regulation Locally through Intrinsic Noise. IET Syst. Biol. 4, 379–392. doi:10.1049/iet-syb.2009.0070
Mat Tek (2021). MatTek Microwell Dishes. Available at: https://www.mattek.com/products/glass-bottom-dishes/(Accessed 01 13, 2021).
Mengel, B., Hunziker, A., Pedersen, L., Trusina, A., Jensen, M. H., and Krishna, S. (2010). Modeling Oscillatory Control in NF-Κb, P53 and Wnt Signaling. Curr. Opin. Genet. Dev. 20, 656–664. doi:10.1016/j.gde.2010.08.008
Mercurio, F., and Manning, A. M. (1999). NF-κB as a Primary Regulator of the Stress Response. Oncogene 18 (45), 6163–6171. doi:10.1038/sj.onc.1203174
Mitchell, S., Vargas, J., and Hoffmann, A. (2016). Signaling via the NFκB System. WIREs Mech. Dis. 8 (3), 227–241. doi:10.1002/wsbm.1331
Nelson, D. E., Ihekwaba, A. E. C., Elliott, M., Johnson, J. R., Gibney, C. A., Foreman, B. E., et al. (2004). Oscillations in NF-Κb Signaling Control the Dynamics of Gene Expression. Science 306, 704–708. doi:10.1126/science.1099962
Nikon (2021). Nikon Microscopes. Available at: https://www.microscope.healthcare.nikon.com/(Accessed 01 14, 2021).
O’Dea, E., and Hoffmann, A. (2009). NF-κB Signaling. Wiley Interdiscip. Rev. Syst. Biol. Med. 1 (1), 107–115.
Oko Cage Incubator (2021). Oko Cage Incubator. Available at: http://www.oko-lab.com/live-cell-imaging/cage-incubator (Accessed 01 14, 2021).
Smoluchowski, M. V. (1916). Über Brownsche Molekularbewegung unter Einwirkung äußerer Kräfte und deren Zusammenhang mit der verallgemeinerten Diffusionsgleichung. Ann. Phys. 353 (24), 1103–1112. doi:10.1002/andp.19163532408
Sneppen, K., Krishna, S., and Semsey, S. (2010). Simplified Models of Biological Networks. Annu. Rev. Biophys. 39, 43–59. doi:10.1146/annurev.biophys.093008.131241
Sugioka, K., Nakagawa, K., Murata, R., Ochiai, N., Sasho, T., Arai, M., et al. (2010). Radial Shock Waves Effectively Introduced Nf-Kappa B Decoy into Rat Achilles Tendon Cells In Vitro. J. Orthop. Res. 28 (8), 1078–1083. doi:10.1002/jor.21081
Tay, S., Hughey, J. J., Lee, T. K., Lipniacki, T., Quake, S. R., and Covert, M. W. (2010). Single-cell NF-Κb Dynamics Reveal Digital Activation and Analogue Information Processing. Nature 466 (7303), 267–271. doi:10.1038/nature09145
Wang, D. M., and Tarbell, J. M. (1995). Modeling Interstitial Flow in an Artery Wall Allows Estimation of Wall Shear Stress on Smooth Muscle Cells. J. Biomech. Eng. 117 (3), 358–363. doi:10.1115/1.2794192
Warner (2021). Warner Temperature Controller. Available at: https://www.warneronline.com/dual-channel-temperature-controller-tc-344c (Accessed 01 14, 2021).
Werner, M., Zhu, L., and Aurell, E. (2007). Cooperative Action in Eukaryotic Gene Regulation: Physical Properties of a Viral Example. Phys. Rev. E Stat. Nonlin Soft Matter Phys. 76, 061909. doi:10.1103/PhysRevE.76.061909
Keywords: dynamical systems, oscillations, gene regulation, controllability, chaos, transcription factors, NF-kB
Citation: Heltberg M, von Borries M, Bendix PM, Oddershede LB and Jensen MH (2022) Temperature Controls Onset and Period of NF-κB Oscillations and can Lead to Chaotic Dynamics. Front. Cell Dev. Biol. 10:910738. doi: 10.3389/fcell.2022.910738
Received: 04 April 2022; Accepted: 18 May 2022;
Published: 20 June 2022.
Edited by:
Fernando Antunes, University of Lisbon, PortugalReviewed by:
Didier Gonze, Université libre de Bruxelles, BelgiumMarc Lefranc, Université de Lille, France
Copyright © 2022 Heltberg, von Borries, Bendix, Oddershede and Jensen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Lene B. Oddershede, oddershede@nbi.ku.dk; Mogens H. Jensen, mhjensen@nbi.dk
†These authors share first authorship