Oscillatory dynamics of p38 activity with transcriptional and translational time delays

Computational analysis of synergism in small networks with different logic

Cell fate decision processes are regulated by networks which contain different molecules and interactions. Different network topologies may exhibit synergistic or antagonistic effects on cellular functions. Here, we analyze six most common small networks with regulatory logic AND or OR, trying to clarify the relationship between network topologies and synergism (or antagonism) related to cell fate decisions. We systematically examine the contribution of both network topologies and regulatory logic to the cell fate synergism by bifurcation and combinatorial perturbation analysis. Initially, under a single set of parameters, the synergism of three types of networks with AND and OR logic is compared. Furthermore, to consider whether these results depend on the choices of parameter values, statistics on the synergism of five hundred parameter sets is performed. It is shown that the results are not sensitive to parameter variations, indicating that the synergy or antagonism mainly depends on the network topologies rather than the choices of parameter values. The results indicate that the topology with "Dual Inhibition" shows good synergism, while the topology with "Dual Promotion" or "Hybrid" shows antagonism. The results presented here may help us to design synergistic networks based on network structure and regulation combinations, which has promising implications for cell fate decisions and drug combinations.

Dynamic behavior of P53-Mdm2-Wip1 gene regulatory network under the influence of time delay and noise

The tumor suppressor protein P53 can regulate the cell cycle, thereby preventing cell abnormalities. In this paper, we study the dynamic characteristics of the P53 network under the influence of time delay and noise, including stability and bifurcation. In order to study the influence of several factors on the concentration of P53, bifurcation analysis on several important parameters is conducted; the results show that the important parameters could induce P53 oscillations within an appropriate range. Then we study the stability of the system and the existing conditions of Hopf bifurcation by using Hopf bifurcation theory with time delays as the bifurcation parameter. It is found that time delay plays a key role in inducing Hopf bifurcation and regulating the period and amplitude of system oscillation. Meanwhile, the combination of time delays can not only promote the oscillation of the system but it also provides good robustness. Changing the parameter values appropriately can change the bifurcation critical point and even the stable state of the system. In addition, due to the low copy number of the molecules and the environmental fluctuations, the influence of noise on the system is also considered. Through numerical simulation, it is found that noise not only promotes system oscillation but it also induces system state switching. The above results may help us to further understand the regulation mechanism of the P53-Mdm2-Wip1 network in the cell cycle.

N-Acetylcysteine or Sodium Selenite Prevent the p38-Mediated Production of Proinflammatory Cytokines by Microglia during Exposure to Mercury (II)

Mercury (Hg) is known for its neurotoxicity and is reported to activate microglia cells at low exposure levels. Since mercury decreases the activity of the glutathione and thioredoxin systems, we hypothesize that Hg would, in turn, disrupt microglia homeostasis by interfering with redox regulation of signaling pathways. Thus, in this work, we analyzed the effect of exposure to Hg²⁺ on nuclear translocation and activation of NF-kB (p50) and p38 and pro-inflammatory gene transcription (IL-1ß; iNOS, TNF-alpha) considering the interaction of Hg with the glutathione system and thioredoxin systems in microglial cells. N9 (mouse) microglia cells were exposed to different concentrations of Hg²⁺ and the 24 h EC₅₀ for a reduction in viability was 42.1 ± 3.7 μM. Subsequent experiments showed that at sub-cytotoxic levels of Hg²⁺, there was a general increase in ROS (≈40%) accompanied by a significant depletion (60-90%) of glutathione (GSH) and thioredoxin reductase (TrxR) activity. Upon 6 h of exposure to Hg²⁺, p38 (but not p50) accumulated in the nucleus (50% higher than in control), which was accompanied by an increase in its phosphorylation. Transcript levels of both IL1-ß and iNOS were increased over two-fold relative to the control. Furthermore, pre-exposure of cells to the p38 inhibitor SB 239063 hindered the activation of cytokine transcription by Hg²⁺. These results show that disruption of redox systems by Hg²⁺ prompts the activation of p38 leading to transcription of pro-inflammatory genes in microglia cells. Treatment of N9 cells with NAC or sodium selenite-which caused an increase in basal GSH and TrxR levels, respectively, prevented the activation of p38 and the transcription of pro-inflammatory cytokines. This result demonstrates the importance of an adequate nutritional status to minimize the toxicity resulting from Hg exposure in human populations at risk.

Circadian Rhythm Regulated by Tumor Suppressor p53 and Time Delay in Unstressed Cells

Circadian function and p53 network are interconnected on the molecular level, but the dynamics induced by the interaction between the circadian factor Per2 and the tumor suppressor p53 remains poorly understood. Here, we constructed an integrative model composed of a circadian clock module and a p53-Mdm2 feedback module to study the dynamics of p53-Per2 network in unstressed cells. As expected, the model can accurately predict the circadian rhythm, which is consistent with diverse experimental observations. In addition, using a combination of theoretical analysis and numerical simulation, the results demonstrated that p53 expression enhances the phase advance of circadian rhythm and reduces the robustness of circadian rhythm. Furthermore, the time delay required for the transcription and translation of Per2 protein induces oscillations by undergoing a supercritical Hopf bifurcation, and improves the robustness of circadian rhythm. In summary, this work shows that the p53-Per2 interaction and the time delay are two essential factors for circadian functions.

Dynamical behaviors of quorum sensing network mediated by combinatorial perturbation

The dynamical behaviors of the quorum sensing (QS) system are closely related to the release drugs and control the PH value in microorganisms and plants. However, the effect of the main molecules AiiA, LuxI, H$ _2 $O$ _2 $, and time delayed individual and combinatorial perturbation on the QS system dynamics and the above-mentioned biological phenomena is still unclear, which are seen as a key consideration in our paper. This paper formulates a QS computational model by incorporating these several substances. First, for the protein production time delay, a critical value is given by Hopf bifurcation theory. It is found that a larger time delay can lead to a larger amplitude and a longer period. This indicates that the length of time for protein synthesis has a regulatory effect on the release of drugs from the bacterial population. Second, hen the concentrations of AiiA, LuxI, and H$ _2 $O$ _2 $ is modulated individually, the QS system undergoes periodic oscillation and bistable state. Meanwhile, oscillatory and bistable regions can be significantly affected by simultaneously perturbing any two parameters related to AiiA, LuxI, and H$ _2 $O$ _2 $. This means that the individual or simultaneous changes of the three intrinsic molecular concentrations can effectively control the drugs release and the PH value in microorganisms and plants. Finally, the sensitivity relationship between the critical value of the delay and AiiA, LuxI, H$ _2 $O$ _2 $ parameters is analyzed.

Dynamic behavior of the p53-Mdm2 core module under the action of drug Nutlin and dual delays

Nutlin is a family of p53-targeting drugs. It is able to bind to Mdm2, thereby accelerate the accumulation of p53 that is a prominent tumor suppressor. An integrated module of the Nutlin PBK and p53 pathway is composed of positive feedback mediated by Mdm2 mRNA as well as the drug Nutlin and negative feedback mediated by Mdm2 protein. The main research content of our paper is how the time delay of protein synthesis, response time delay of Nutlin drug, the degradation rate of Mdm2, the degradation rate of p53 depended on Mdm2 and the actual dose of Nutlin in the cell affect the oscillatory behavior caused by Hopf bifurcation in the integrated network system of Nutlin PBK and p53 pathways. The stability of the unique positive equilibrium point and the existence of Hopf bifurcation are studied by taking the time delays as the bifurcation parameters and applying bifurcation theory. Based on the normal form theory and central manifold theorem, explicit criteria to determine the Hopf bifurcation direction and stability of the bifurcated periodic solution are established. It is found that the time delays and key parameters in the integrated network system of Nutlin PBK and p53 pathways play an important role in the amplitude and period of p53 oscillation according to the results from the numerical simulation and theoretical calculation. These results may provide us with a better understanding of the biological functions of the p53 pathway and some clues for cancer treatment.

The common misuse of noise decomposition as applied to genetic systems

The noise-decomposition technique is applied in several fields, including genetic systems, optical images, recording, and navigation. In genetic systems, noise decomposition is usually achieved by using two reporters [Elowitz M.B., Levine A.J., Siggia E.D., Swain P·S., 2002. Stochastic gene expression in a single cell. Science 297, 1183-6.]. A reporter is a protein with fluorescence, an RNA hybridized with a fluorescent probe, or any other detectable intracellular component. If a reporter is constructed in addition to the original reporter, the system's stochasticity may change. Such phenomena became severe for genes in plasmids with a high copy number. By SSA (stochastic simulation algorithm), we observed an approximately 50% increment in the coefficient of variation while introducing additional reporters. Besides, if two reporters respond to the upstream element at a different time, the trunk noise (or extrinsic noise) cannot be accurately determined. This is because the "calculative trunk noise" changes along with the delay, though the real trunk noise does not. For RNA reporters, a 5-min transcriptional delay caused a calculative trunk noise that was 90% less than the real trunk noise. Fortunately, this problem is negligible when the degradation rate constant is low, and it is usually true in the case of the protein reporters. One can check the lifespan of the reporter before applying the noise-decomposition technique.

Combinatorial dynamics of protein synthesis time delay and negative feedback loop in NF-κB signalling pathway

The transcription factor NF-κB links immune response and inflammatory reaction and its different oscillation patterns determine different cell fates. In this study, a mathematical model with IκBα protein synthesis time delay is developed based on the experimental evidences. The results show that time delay has the ability to drive oscillation of NF-κB via Hopf bifurcation. Meanwhile, the amplitude and period are sensitive to the time delay. Moreover, the time delay threshold is a function of four parameters characterising the negative feedback loop. Likewise, the parameters also have effects on the amplitude and period of NF-κB oscillation induced by time delay. Therefore, the oscillation patterns of NF-κB are collaborative results of time delay coupled with the negative feedback loop. These results not only enhance the understanding of NF-κB biological oscillation but also provide clues for the development of anti-inflammatory or anti-cancer drugs.

Bifurcation and oscillatory dynamics of delayed CDK1-APC feedback loop

Extensive experimental evidence has been demonstrated that the dynamics of CDK1-APC feedback loop play crucial roles in regulating cell cycle processes, but the dynamical mechanisms underlying the regulation of this loop are still not completely understood. Here, the authors systematically investigated the stability and bifurcation criteria for a delayed CDK1-APC feedback loop. They showed that the maximum reaction rate of CDK1 inactivation by APC can drive sustained oscillations of CDK1 activity ([inline-formula removed]) and APC activity ([inline-formula removed]), and the amplitude of these oscillations is increasing with the increase of the reaction rate over a wide range; a certain range of the self-activation rate for CDK1 is also significant for generating these oscillations, for too high or too low rates the oscillations cannot be generated. Moreover, they derived the sufficient conditions to determine the stability and Hopf bifurcations, and found that the sum of time delays required for activating CDK1 and APC can induce [inline-formula removed] and [inline-formula removed] to be oscillatory, even when the [inline-formula removed] and [inline-formula removed] settle in a definite stable steady state. Furthermore, they presented an explicit algorithm for the properties of periodic oscillations. Finally, numerical simulations have been presented to justify the validity of theoretical analysis.

Oscillatory Dynamics of p53-Mdm2 Circuit in Response to DNA Damage Caused by Ionizing Radiation

Although the dynamical behavior of the p53-Mdm2 loop has been extensively studied, the understanding of the mechanism underlying the regulation of this pathway still remains limited. Herein, we developed an integrated model with five basic components and three ubiquitous time delays for the p53-Mdm2 interaction in response to DNA damage following ionizing radiation (IR). We showed that a sufficient amount of activated ATM level can initiate the p53 oscillations with nearly the same amplitude over a wide range of the ATM level; a proper range of p53 level is also required for generating the oscillations, for too high or too low levels it would fail to generate the oscillations; and increased Mdm2 level leads to decreased amplitude of the p53 oscillation and reduced expression of the p53 activity. Moreover, we found that the negative feedback loop formed between p53 and nuclear Mdm2 plays a dominant role in determining the p53 dynamics, whereas when interaction strength of the negative feedback loop becomes weaker, the positive feedback loop formed between p53 and cytoplasmatic Mdm2 can induce different types of dynamics. Furthermore, we demonstrated that the total time delay required for protein production and nuclear translocation of Mdm2 can induce p53 oscillations even when the p53 level is at a certain stable high steady state or at a certain stable low steady state. In addition, the two important features of the oscillatory dynamics-amplitude and period-can be controlled by such time delay. These results are in agreement with multiple experimental observations and may enrich our understanding of the dynamics of the p53 network.

Modelling and analysing biological oscillations in quorum sensing networks

Recent experiments have shown that the biological oscillation of quorum sensing (QS) system play a vital role not only in the process of bacterial synthesis but also in the treatment of cancer by releasing drugs. As known, these five substances TetR, CI, LacI, AiiA and AI are the core components of the QS system. However, the effects of AiiA and protein synthesis time delay on QS system are often ignored in the theoretical model, which is taken as a priority in the proposed research. Therefore, the authors developed a new mathematical model to explore the effects of AiiA and time delay on the dynamical behaviour of QS system theoretically and numerically. The results show that time delay can induce oscillation of QS system. Concretely, there exists a time delay threshold [inline-formula removed]. When time delay is less than [inline-formula removed], the system is stable. With the increasing of time delay and once it passes [inline-formula removed], oscillation behaviour occurs. Moreover, the length of time delay determines the amplitude and period of the QS oscillation. In addition, the value of [inline-formula removed] is sensitive to AiiA. These results may enhance the understanding of QS oscillations and provide new insights for bacterial release drugs to treat cancer.

Oscillatory behavior of p53-Mdm2 system driven by transcriptional and translational time delays

Biological experiments clarify that p53-Mdm2 module is the core of tumor network and p53 oscillation plays an important role in determining the tumor cell fate. In this paper, we investigate the effect of time delay on the oscillatory behavior induced by Hopf bifurcation in p53-Mdm2 system. First, the stability of the unique positive equilibrium point and the existence of Hopf bifurcation are investigated by using the time delay as the bifurcation parameter and by applying the bifurcation theory. Second, the explicit criteria determining the direction of Hopf bifurcation and the stability of bifurcating periodic solutions are developed based on the normal form theory and the center manifold theorem. In addition, the combination of numerical simulation results and theoretical calculation results indicates that time delays in p53-Mdm2 system are critical for p53 oscillations. The results may help us to better understand the biological functions of p53 pathway and provide clues for treatment of cancer.

The role of time delays in P53 gene regulatory network stimulated by growth factor

In this paper, a delayed mathematical model for the P53-Mdm2 network is developed. The P53-Mdm2 network we study is triggered by growth factor instead of DNA damage and the amount of DNA damage is regarded as zero. We study the influences of time delays, growth factor and other important chemical reaction rates on the dynamic behaviors in the system. It is shown that the time delay is a critical factor and its length determines the period, amplitude and stability of the P53 oscillation. Furthermore, as for some important chemical reaction rates, we also obtain some interesting results through numerical simulation. Especially, S (growth factor), k₃ (rate constant for Mdm2_p dephosphorylation), k₁₀ (basal expression of PTEN) and k₁₄ (Rate constant for PTEN-induced Akt dephosphorylation) could undermine the dynamic behavior of the system in different degree. These findings are expected to understand the mechanisms of action of several carcinogenic and tumor suppressor factors in humans under normal conditions.

Stability and Oscillation Analysis of a Gene Regulatory Network With Multiple Time Delays and Diffusion Rate

In genetic regulatory networks (GRNs), the diffusion rate of mRNA and protein play a key role in regulatory mechanisms of gene expression, especially in translation and transcription. However, the influence of diffusion rate on oscillatory gene expression is not well understood. In this paper, by considering the diffusion rate of mRNA and protein, a novel GRN is proposed. Then, two basic problems of such network, i.e. stability and oscillation, are solved in detail. Moreover, the properties of oscillation are also investigated. it is found that the total biochemistry reaction time can affect the stability of the positive equilibrium and give rise to the oscillation. The diffusion rate of mRNA and proteins have a major impact on the oscillation properties. Finally, two examples not only verify the theoretical results, but also show that a slight diffusion rate increasing may lead to huge change in oscillatory gene expressions.

Oscillation induced by Hopf bifurcation in the p53-Mdm2 feedback module

This study develops an integrated model of the p53-Mdm2 interaction composed of five basic components and time delay in the DNA damage response based on the existing research work. Some critical factors, including time delay, system parameters, and their interactions in the p53-Mdm2 system are investigated to examine their effects on the oscillatory behaviour induced by Hopf bifurcation. It is shown that the positive feedback formed between p53 and the activity of Mdm2 in the cytoplasm can cause a slight decrease in the amplitude of the p53 oscillation. The length of the time delay plays an important role in determining the amplitude and period of the oscillation and can significantly extend the parameter range for the system to demonstrate oscillatory behaviour. The numerical simulation results are found to be in good agreement with the published experimental observation. It is expected that the results of this research would be helpful to better understand the biological functions of p53 pathway and provide some clues in the treatment of cancer.

Oscillatory dynamics mechanism induced by protein synthesis time delay in quorum-sensing system

Recent experimental evidence reports that the oscillatory behavior of quorum sensing plays an extremely important role in the process of bacterial synthesis and release drug to fight cancer. As we know, the six substances AiiA, LuxI, internal AHL, external AHL, AHL substrate, and H_{2}O_{2} are the core parts of the quorum-sensing system. Here, the effects of several important factors, including time delay, variable H_{2}O_{2}, AHL synthesis rate induced by LuxI, and AHL degradation rate induced by AiiA on the oscillatory behavior of the quorum-sensing system are studied theoretically based on a part of mathematical model describing the interaction of the above six substances proposed by Prindle et al. [Nature 508, 387 (2014)10.1038/nature13238]. The results show that the time delay is a prerequisite for inducing oscillation of the quorum-sensing system. Furthermore, the length of time delay can determine the amplitude and period of oscillation. As a further matter, the change of H_{2}O_{2} concentration can induce the oscillatory behavior of the quorum-sensing system. In addition, under the regulation of H_{2}O_{2}, the period robustness of the quorum-sensing system is increased. Similarly, the quorum-sensing system exhibits periodic oscillation when AHL synthesis rate induced by LuxI less than a certain critical value, unless it displays a steady state. Additionally, a too-high or too-low level of AHL degradation rate induced by AiiA will fail to generate oscillation of the quorum-sensing system, only the intermediate level will cause oscillation. Finally, the two and three parameter regions in which the quorum-sensing system exhibits oscillation behavior are generated.

Contribution of time delays to p53 oscillation in DNA damage response

Although the oscillatory dynamics of the p53 network have been extensively studied, the understanding of the mechanism of delay-induced oscillations is still limited. In this paper, a comprehensive mathematical model of p53 network is studied, which contains two delayed negative feedback loops. By studying the model with and without explicit delays, the results indicate that the time delay of Mdm2 protein synthesis can well control the pulse shape but cannot induce p53 oscillation alone, while the time delay required for Wip1 protein synthesis induces a Hopf bifurcation to drive p53 oscillation. In addition, the synergy of the two delays will cause the p53 network to oscillate in advance, indicating that p53 begins the repair process earlier in the damaged cell. Furthermore, the stability and bifurcation of the model are addressed, which may highlight the role of time delay in p53 oscillations.

Theoretical study on the oscillation mechanism of p53-Mdm2 network

In this paper, a delayed mathematical model was developed based on experimental data to understand how the time delays required for transcription and translation in Mdm2 gene expression affect the kinetic behavior of the p53-Mdm2 network. Taking the time delays as the main research parameters, the stability of the system at the positive equilibrium was studied by using the theoretical method of delay differential equation. We found that such delays can induce oscillations by undergoing a supercritical Hopf bifurcation. Then, we used the normal form theory and the center manifold reduction to study the direction and stability of the bifurcation in detail. Furthermore, we also studied the effects of the length of time delays and the model parameters by numerical simulations. We found that time delays in Mdm2 synthesis are required for p53 oscillations and the length of such delays can determine the amplitude and period of the oscillations. In addition, the model parameters can also change the stability of the system. These results illustrate that the repair process after DNA damage can be regulated by varying time delays and the model parameters.