Research article Special Issues

Global stability for a class of functional differential equations with distributed delay and non-monotone bistable nonlinearity

  • The present work is devoted to the global stability analysis for a class of functional differential equations with distributed delay and non-monotone bistable nonlinearity. First, we characterize some subsets of attraction basins of equilibria. Next, by Lyapunov functional and fluctuation method, we obtain a series of criteria for the global stability of equilibria. Finally, we illustrate our results by applying them to a problem with Allee effect.

    Citation: Toshikazu Kuniya, Tarik Mohammed Touaoula. Global stability for a class of functional differential equations with distributed delay and non-monotone bistable nonlinearity[J]. Mathematical Biosciences and Engineering, 2020, 17(6): 7332-7352. doi: 10.3934/mbe.2020375

    Related Papers:

    [1] Gheorghe Craciun, Baltazar Aguda, Avner Friedman . Mathematical Analysis Of A Modular Network Coordinating The Cell Cycle And Apoptosis. Mathematical Biosciences and Engineering, 2005, 2(3): 473-485. doi: 10.3934/mbe.2005.2.473
    [2] Floriane Lignet, Vincent Calvez, Emmanuel Grenier, Benjamin Ribba . A structural model of the VEGF signalling pathway: Emergence of robustness and redundancy properties. Mathematical Biosciences and Engineering, 2013, 10(1): 167-184. doi: 10.3934/mbe.2013.10.167
    [3] Cicely K. Macnamara, Mark A. J. Chaplain . Spatio-temporal models of synthetic genetic oscillators. Mathematical Biosciences and Engineering, 2017, 14(1): 249-262. doi: 10.3934/mbe.2017016
    [4] Alexis B. Cook, Daniel R. Ziazadeh, Jianfeng Lu, Trachette L. Jackson . An integrated cellular and sub-cellular model of cancer chemotherapy and therapies that target cell survival. Mathematical Biosciences and Engineering, 2015, 12(6): 1219-1235. doi: 10.3934/mbe.2015.12.1219
    [5] Bojie Yang, Zhuoqin Yang . Deterministic and stochastic approaches to a minimal model for the transition from autophagy to apoptosis. Mathematical Biosciences and Engineering, 2024, 21(2): 3207-3228. doi: 10.3934/mbe.2024142
    [6] Kewei Ni, Gaozhong Sun . The identification of key biomarkers in patients with lung adenocarcinoma based on bioinformatics. Mathematical Biosciences and Engineering, 2019, 16(6): 7671-7687. doi: 10.3934/mbe.2019384
    [7] Li-Bing Wang, Liang Feng, Jing He, Bo Liu, Jian-Guang Sun . MiR-125a-5p inhibits the proliferation and invasion of breast cancer cells and induces apoptosis by targeting GAB2. Mathematical Biosciences and Engineering, 2019, 16(6): 6923-6933. doi: 10.3934/mbe.2019347
    [8] Juenu Yang, Fang Yan, Haihong Liu . Dynamic behavior of the p53-Mdm2 core module under the action of drug Nutlin and dual delays. Mathematical Biosciences and Engineering, 2021, 18(4): 3448-3468. doi: 10.3934/mbe.2021173
    [9] Olga Vasilyeva, Tamer Oraby, Frithjof Lutscher . Aggregation and environmental transmission in chronic wasting disease. Mathematical Biosciences and Engineering, 2015, 12(1): 209-231. doi: 10.3934/mbe.2015.12.209
    [10] Tao Chen, Zhiming Li, Ge Zhang . Analysis of a COVID-19 model with media coverage and limited resources. Mathematical Biosciences and Engineering, 2024, 21(4): 5283-5307. doi: 10.3934/mbe.2024233
  • The present work is devoted to the global stability analysis for a class of functional differential equations with distributed delay and non-monotone bistable nonlinearity. First, we characterize some subsets of attraction basins of equilibria. Next, by Lyapunov functional and fluctuation method, we obtain a series of criteria for the global stability of equilibria. Finally, we illustrate our results by applying them to a problem with Allee effect.


    Apoptosis is a natural physiological process of programmed cell death. It helps an organism to eliminate cells which are no longer needed or irreparably damaged. Apoptosis plays a critical role in embryonic development and adult tissue homeostasis. Deregulation of apoptosis is an important aspect of cancer pathogenesis [1], and has been widely recognized as a hallmark of most types of cancer [2]. In general, apoptosis is strictly regulated by the body to maintain the stability of tissues and organs throughout the life cycle. However, when the regulation of apoptosis is out of balance, excessive proliferation or excessive apoptosis may occur, thus leading to the occurrence of related diseases such as ovarian cancer and endometriosis [3,4]. Therefore, deeper understanding the mechanisms of apoptosis contributes to the study of the pathogenesis of related diseases and plays an important role in the prevention and treatment of them.

    There are three major signal transduction pathways in cell apoptosis: the death receptor, the mitochondrial, and the endoplasmic reticulum pathways, of which the former belongs to the outer pathway and the latter two belong to the intrinsic pathways. Here we only study the mechanism of apoptosis related to the intrinsic pathway. Caspases are a class of cysteine-aspartyl proteases that are synthesized as inactive precursor enzymes or proenzymes. These proteases normally lie dormant in healthy cells. In response to diverse stimuli they become activated, by either proteolytic cleavage or by recruitment into large complexes. Once activated, caspases cleave their substrates typically after conserved aspartate residues and are responsible for most of biochemical and morphological features of apoptosis. The p53 tumor suppressor gene, which acts as a key conduit in cells between DNA damage detection and apoptosis initiation, is frequently inactivated by mutations in cancer. Inactivation of p53 enables cancer cells not only to bypass intrinsic apoptotic responses to their underlying genomic aberrations, but also to block apoptosis initiation in response to DNA damage they incur upon exposure to conventional therapies [5]. Therefore, agents that can restore apoptosis in cancer cells hold promise for therapy and have been the focus of many preclinical drug discovery studies. Developmental cues or severe cell stress, such as DNA or cytoskeletal damage, stimulate the intrinsic pathway by inducing transcription or post-translational activation of BH3-only proapoptotic BCL-2 family proteins, which, in turn, drive activation of other family members such as BAX and BAK. In healthy cells, antiapoptotic family members such as BCL-2, BCL-XL, or myeloid cell leukemia sequence 1 (MCL-1) bind to proapoptotic members such as BAX and BAK to prevent their activation. Upon stimulation, BH3-only proteins sequester the antiapoptotic proteins away, allowing BAX and BAK to act. BAX and BAK engage the mitochondria to release factors that promote caspase activation in the cytoplasm. One such factor is cytochrome c, which cooperates with another protein called Apaf-1 to activate caspase 9. This caspase, in turn, activates the effector caspases 3, 6, and 7, which carry out apoptosis [6].

    A mathematical model of the intrinsic pathway that addresses specifically the dynamical basis of the three physiological features of apoptosis was presented [7]. And they subdivided the intrinsic pathway into three modules (initiator, amplifier, and executioner modules). Based on prior work, we identify important regulatory processes in apoptosis through single-parameter sensitivity analysis and bifurcation theory. Parameter sensitivity analysis has been used by many researchers [8] and has been supported by some experiments. The work of [9] is encouraging in this regard. The authors studied the mitochondrial outer membrane permeabilization network that controls apoptosis, and found that the oncogenic mutation of each protein functional domain is strongly correlated with the parameter sensitivity of the bifurcation point. In [10], the authors determined the location of the bifurcation point rather than the concentration of some specified proteins by identifying parameters that have a large influence on the bifurcation points. The bifurcation-based sensitivity is a better measure for determining the function of the bio-network and plays an important role in the field of genetic mutation [8,9,10]. Fortunately, we found that modules or feedback loops but not specific parameters determine the network functions although similar bifurcation-based sensitivity analysis is performed.

    In previous studies, there were many models of apoptotic pathways. In this paper, we use the model described in Figures 1 and 2 to analyze. For more detailed information of the model, please refer to literature [7]. The model of the intrinsic pathway describes the dynamics of the three characteristic physiological features of apoptosis. In Figure 1, the intrinsic path is subdivided into three modules (initiators, amplifiers, and actuator modules). It also shows how the positive feedback in the initiator module is responsible for signal thresholds and time delays, and how the actuator module ensures commitment to cell death. The initiator module describes BAX activation by BH3 proteins. The amplifier module describes how BAXm causes the release of CytoC and SMAC. The executioner module describes how CytoC and SMAC cooperate to activate caspase. The schematic draft of the network is shown in Figure 2. Here, we mainly analyze three modules: initiator module, and executioner modules and the how p53 module governs cell fates after DNA damage.

    Figure 1.  Wiring diagrams for a three-module decomposition of the intrinsic apoptotic pathway. Initiator module in the green box, amplifier module in the red box, and executioner module in the blue box.
    Figure 2.  Wiring diagrams for processes regulating p53 and those regulated by p53.

    We perform bifurcation and sensitivity analysis to investigate dynamics of the system. Because the BAXm concentration, i.e., [BAXm], is an important measure of the initiator module (startup the executioner module), we use stress as the parameter and [BAXm] as the output, then we can get bifurcation diagram, as shown in Figure 3(a). As the stress increases from low to high, [BAXm] stays at low levels until the stress exceeds a threshold SN1, at which a saddle-node bifurcation occurs and the [BAXm] increases suddenly. Such a bifurcation corresponds to the conversion of the initiator module to the executioner module. Similarly, the other saddle-node bifurcation occurs at SN2 if the stress is decreased from high to low.

    Figure 3.  Bifurcation diagrams with stress as a control parameter. (a) Two solid lines indicate stable states and dashed line shows unstable states. SN1 is one transition point when the stress increases from low to high and SN2 is the other transition point as stress decreases instead. (b) The change of the locations of the bifurcation points SN1 with increasing the parameter kd3. Black, control; blue, 0.5-fold increase; green, 1-fold increase.

    At the high stress levels, the system starts producing lots of BH3 protein. At first, most [BH3] molecules bind with [BCL2], so they cannot activate [BAX]. Eventually, BH3 titrates out the BCL2 pool, and excess BH3 molecules are now free to convert [BAX] to [BAXm]. BAXm displaces additional BH3 from the BH3/BCL2 dimer pool, thereby speeding up its own production from inactive [BAX]. This positive feedback causes the [BAXm] level to rise quickly, which is sufficient to activate the amplifier module [7]. When the stress decreases from high to low, there is another saddle-node bifurcation point SN2. When the stress lies between SN1 and SN2, the module is bistable. The bistable states represent two different functional states. The higher level of [BAXm] stands for the initiator module and the lower level indicates the amplifier module.

    It is known that bifurcations exhibit qualitative changes in biological systems and are found in many biological functions [11,12]. In this system, the saddle-node bifurcation provides a mechanism of switching between modules, i.e., a switch mechanism that ultimately indicates life-to-death. Bifurcation point location has an important effect on apoptosis. However, some parameter changes may lead to large changes in positions of bifurcation points, resulting in abnormal system functions. Therefore, it is an important issue to determine which parameters have more significant influence on the positions of bifurcation points. In this paper, we determine the important parameters and therefore important regulatory processes in fate decisions through single parameter sensitivity analysis of the bifurcation point. For instance, the bifurcation point SN1 is clearly shifted to the right by increasing the parameter kd3, i.e., the BH3 degradation rate, by 0.5 fold, as shown in Figure 3(b). Increasing kd3 to 1 fold of the basal value, SN1 continues to move to the right, and high [BAXm] can be achieved at higher stress, which means that if BH3 accelerates self-degeneration due to mutations in certain genes, even if the stress is very high, apoptosis will be difficult to carry out. Sufficient stress has no way to induce apoptosis, providing the possibility of tumorigenesis. Similarly, we performed bifurcation analysis of the executor and the p53 modules. This conclusion is similar to the initiator module, as shown in Figures 4 and 5.

    Figure 4.  Bifurcation diagrams with [CytoC] as a control parameter. (a) Two solid lines indicate stable states and dashed line shows unstable states. SN1 is one transition point when [CytoC] increases from low to high and SN2 is the other transition point as [CytoC] decreases instead. (b) The change of the locations of the bifurcation points SN1 with increasing the parameter kdac9. Black, control; blue, 0.5-fold increase; green, 1-fold increase.
    Figure 5.  Bifurcation diagrams with [p53killer] as a control parameter. (a) Two solid lines indicate stable states and dashed line shows unstable states. SN1 is one transition point when [p53killer] increases from low to high and SN2 is the other transition point as [p53killer] decreases instead. (b) The change of the locations of the bifurcation points SN1 with increasing the parameter k2s3dr. Black, control; blue, 0.5-fold increase; green, 1-fold increase.

    As far as we know, the influence of parameter variations on dynamics is uneven in many biological systems [13,14]. And the effects of parameter variations on the positions of the bifurcation points have not been well investigated. There are many parameters in the model that can affect the positions of bifurcation points. In order to determine which parameters have more significant influences on the locations of bifurcation points in the model, we perform bifurcation analysis and record the percentage change in bifurcation points by increasing or decreasing each parameter by 10% off its base value. Thus, single parameter sensitivity of bifurcation points are derived and further the sensitivity spectrum of all parameters can be achieved.

    BCL2 is an anti-apoptotic gene that prevents apoptosis caused by various stimulating factors and promotes cell survival. When a stress signal is applied to the initiator module, the BH3 level increases. At first, these BH3 molecules are taken out of service by binding to the inhibitory protein, BCL2. If the stress signal is sustained, enough BH3 will accumulate to saturate the pool of anti-apoptotic [BCL2] molecules. Free and active BH3 induces conformational changes in BAX which facilitates its insertion into the mitochondrial outer membrane. The membrane-localized form of BAX, i.e., BAXm, binds strongly to BCL2, releasing additional free BH3 to promote the conversion of BAX to BAXm [7,15,16,17].

    The initiator module can be described by four ordinary differential equations (ODEs), once we assume constant total concentrations of BAX and BCL2 proteins. All numerical simulations and bifurcation analysis are performed based on these equations [7]. The parameter sensitivity analysis is shown in Figure 6. All the saddle-node bifurcation points are distributed under the standard parameter set and the perturbed parameter set with each parameter increased or decreased by 10% off its base value, as shown in Figure 6(a). Notably, one sees that the all SN2 are always positive, which means the switch is robustly reversible under the given perturbations.

    Figure 6.  Parameter sensitivity analysis of the bifurcation points. (a) The level of [BAXm] at the two saddle-node bifurcation points SN1-2 when each of 10 parameters in the feedback loops is increased or decreased by 10% with respect to the standard parameter set. (b) Percentage of changes in the levels of [BAXm] at SN1-2 compared to that under the standard parameter set.

    We aim to identify specific proteins through sensitive analysis, and then find out important feedback loops or modules in apoptosis. Due to very few parameters in the initiator module, we analyze the role of individual proteins. We then further calculate the percentage of changes in the two thresholds, as shown in Figure 6(b). Since the relative influence of the parameters on the bifurcation is basically unchanged, we distinguish the strong and weak influences according to the overall change of the parameters. Here we choose 10% as the standard. Most of the parameters in initiator module (approximately 70%) are insensitive, i.e., perturbations to these parameters result in small changes in the bifurcation. But some (approximately 30%) are sensitive and perturbation to them has a strong influence on the bifurcation, e.g., kd3 (decomposition rate of BH3), ks32 (stress stimulates the speed of BH3), and kb (the conversion rate from BAX to BAXm). Increasing parameter kd3 causes the SN1 point to move to the right, which requires a higher stress level to achieve a high level of BAXm than the basal rate. In other words, more stress is needed to start the amplifier module.

    As the free BAX increases, the oligomers form and open channels in the outer membrane of the granules. CytoC and SMAC leave the mitochondria and enter the cytoplasm through these channels, where they cooperate to activate the execution daughter caspase 3 and then eventually degenerate. Given the parameter settings in the model, the amplifier module only transfers the signal from the mitochondrial BAX to the caspase in the cytoplasm [7,18,19]. Therefore, we do not conduct in-depth research on the amplifier module.

    The executioner module describes caspase activation by [CytoC] and [SMAC]. In order to further study the executive module, we use parameter sensitivity analysis again. The executioner module uses caspase 3 activity as an indicator of cell death [7]. Therefore we use CytoC as the input variable and caspase 3 activity as the output. The results of parameter sensitivity are shown in Figure 7. To map cancer-associated mutations to network dynamics through protein-protein interaction dynamics, we have narrowed our research to the parameters of protein interactions. We separated the parameters into the following three clusters according to their impact on the bifurcation point locations: most sensitive (red and blue), sensitive (green), and insensitive (black). In other words, the thresholds for the switch are most sensitive to the parameters in both the aC9 autocrine expression and C9/XIAP feedback loop, whereas the thresholds for the switch are sensitive to the parameters in the C3/XIAP feedback loop. Caspase 9 plays a more important role in the intrinsic pathway of apoptosis. A large number of experimental articles have also shown that apoptosis can be promoted by inducing caspase 9. For example, BMP4 inhibits caspase 9 via the PI3K/AKT signaling pathway, thereby inhibiting apoptosis of GCs [20]. Caspase 9 and AIF play important roles in the elongation-induced apoptosis of myoblasts, and AIF participates in this process in a manner independent of caspase-9 [21].

    Figure 7.  Parameter sensitivity analysis. 10% decrease or increase of each parameter in the executioner module causes a percentage change in the positions of the bifurcation points SN1 and SN2. Most sensitive (red and blue): Red, aC9 autocrine expression, blue, C9/XIAP feedback loop; sensitive (green), C3/XIAP feedback loop; and insensitive (black), the rest proteins.

    After DNA damage, the p53 module dominates the cell fates. We want to explore which proteins play a decisive role in the process. From the study of the initiator module, we know that the level of [BAXm] is an important measure of the initiator module related to the executioner module, we use DNA damage as the input variable and [BAXm] as the output, then we can get the bifurcations diagram. Similar to previous analysis, we perform the bifurcation analysis and record the percentage change in bifurcation locations by increasing each parameter in the model by 10%. The sensitivity spectrum of all parameters is achieved in Figure 8. Similar results can be obtained for SN2 and the case of decreasing each parameter 10% off its base value.

    Figure 8.  Parameter sensitivity analysis. 10% increase of each parameter in the p53 module causes a percentage change in the position of the bifurcation points SN1. Sensitive (red, yellow and green): Red, BH3AC/BH3DR feedback loop, yellow, MDM2/ARF/ASPP feedback loop, green, p53helper/p53killer feedback loop; and insensitive (blue), the rest proteins.

    In this module we choose 4% as the standard. Of course, other criteria can be selected, but the relative influence of the parameters on the bifurcation point is substantially unchanged. From Figure 8, most of the parameters (approximately 82%) are insensitive: perturbations in these parameters result in small changes in the bifurcation. But some (approximately 18%) are sensitive: perturbing these parameters individually has a strong influence on the bifurcation. For sensitive parameters, we observed that the proteins represented by the sensitive parameters mainly focus on three feedback loops: Red, BH3AC/BH3DR feedback loop; yellow, MDM2/ARF/ASPP feedback loop; green, p53helper/p53killer feedback loop. In other words, these three feedback loops play a key role in the cell fate of p53-regulated after DNA damage. It has been experimentally proved that selection of a suitable combination of BH3 mimetic compounds can induce apoptosis of basophils and mast cells in mice [22]. UCHL1 activates the MDM2-ARF-p53 signaling pathway, which has two opposite enzymatic activities in the ubiquitin pathway, further leading to tumor cell apoptosis [23]. The results of the above experiments also indicate that BH3 plays an important role in inducing apoptosis, consistent with the conclusion of BH3AC/BH3DR feedback loop in this paper.

    To investigate the connection between core functions of biological networks and oncogenesis, we conducted sensitivity analysis of all the parameters in the apoptosis model. We found the sensitive parameters that have greater influence on the bifurcation locations. The positions of bifurcation point represent threshold of the initiator module triggering the amplifier module. Only when the stress signal exceeds a specific value, the system will induce the transition from the initiator module to the amplifier module. If BH3 accelerates self-degeneration due to mutations in certain genes, even if the stress is very high, apoptosis will be difficult to carry out. Sufficient stress has no way to induce apoptosis, providing the possibility of tumorigenesis.

    In this paper, we determine the important regulatory processes through single parameter sensitivity analysis of the bifurcation points. In executioner module, we found that aC9 autocrine expression and C9/XIAP feedback loop is more sensitive than the C3/XIAP feedback loop. In addition, we found that BH3AC/BH3DR feedback loop, MDM2/ARF/ASPP feedback loop, and p53helper/p53killer feedback play more important roles in cell fates after DNA damage. These results mean that modules or feedback loops but not detailed parameter values determine the network functions. Such results are partially supported by experimental findings. Our research provides a rational basis for the design of clinical trials of anticancer drugs.

    The authors thank the anonymous reviewer for his valuable comments and suggestions which are helpful to improve this paper. This work was supported by the National Science Foundation of Shanghai (Grant No. 17ZR1410800).

    The authors declare there is no conflict of interest.

    Initiator module

    We use x1, x2, x3, and x4 to denote [BAXm], [BAXm/BCL2], [BH3], and [BH3/BCL2], so we can obtain the dynamical equations as follows

    dx1dt=(kf1+kf2×bh3)×xfkb×x1,dx2dt=kam2×mf×bcl2kdm2×x2kb×x2,dx3dt=ks31+ks32×stresskd3×x3,dx4dt=ka32×bh3×bcl2kd32×x4kd3×x4,

    where xf=xtx1, mf=x1x2, bh3=x3x4, and bcl2=rtx2x4.

    The parameter values are chosen as

    stress=0.1, rt=80, xt=100, kf1=1

    kf2=3, kb=2, kam2=90, kdm2=0.05

    ka32=10, kd32=0.01, ks31=0.1, ks32=0.6, and kd3=0.01.

    The initial values are

    x1=33.4, x2=33.4, x3=16, and x4=16.

    Executioner module

    We use x5, x6, x7, x8, x9, and x10 to denote [C9], [aC9], [C3], [XIAP/C9], [XIAP/C3], and [SMAC/XIAP], respectively. Then we get the dynamics as follows

    dx5dt=ksc9×(c9tx5x8x6)×(cytoc)nkdc9×x5(k1ac9+k2ac9×(x7)n)×x5kas9x×x5×pf+kds9x×x8,dx6dt=(k1ac9+k2ac9×(x7)n)×x5kdac9×x6,dx7dt=(k1ac3+k2ac3×(x5)n+k3ac3×(x6)n)×(c3tx7x9)kdc3×x7kas3x×x7×pf+kds3x×x9,dx8dt=kas9x×x5×pfkds9x×x8kd9x×x8,dx9dt=kas3x×x7×pfkds3x×x9kd3x×x9,dx10dt=kassx×cf×pfkdssx×x10kdsx×x10,

    where cf=smaccytotx10, pf=ptx10x9x8.

    Parameter values are

    cytoc=0.1, smaccytot=2, k1ac9=0.001, k2ac9=0.5

    ksc9=0.001, c9t=1, kdc9=0.002, kdac9=0.003

    k1ac3=0.001, k2ac3=0.02, k3ac3=0.5, c3t=1

    kdc3=0.002, kas3x=0.2, kds3x=0.5, kd3x=0.1

    kas9x=0.1, kds9x=0.6, kd9x=0.2, kassx=2

    kdssx=0.01, kdsx=0.007, pt=6, and n=2.

    The initial values are

    x5=0, x6=0, x7=0, x8=0

    x9=0, and x10=0.

    P53 module

    We use x11, x12, x13, x14, x15, x16, x17, x18, x19, x20, x21, x22, x23, x24, x25, x26, x27, x28, x29, x30, and x31 to denote [DYRK2], [p53helper], [p53killer], [ASPP], [MDM2], [ARF], [ARF/MDM2], [CycE], [E2F1], [RBp], [p21], [p21/CycE], [BAX], [BAXm], [BCL2], [BH3AC], [BH3DR], [BAXm/BCL2], [BH3AC/BCL2], [BH3DR/BCL2], and [p73]. The dynamics of the P53 model can be expressed by the following equations

    dx11dt=kin×dna×(dyrtx11)kout×x11,dx12dt=ksp53(k1dp53+x151+0.1×dna)×x12kpp53×x11×x12jpp53+x12+kdpp53×pp×x13jdpp53+x13,dx13dt=kpp53×x11×x12jpp53+x12kdpp53×pp×x13jdpp53+x13(k1dp53+x151+0.1×dna)×x13,dx14dt=k1saspp+k2saspp×x19kdaspp×x14,dx15dt=k1smdm2+k2smdm2×p53t(k1dmdm2+k2dmdm2×dna)×x15,kasam×x16×x15+(kdsam+kdarf)×x17,dx16dt=k1sarf+k2sarf×x19kdarf×x16kasam×x16×x15,+kdsam×x17+(k1dmdm2+k2dmdm2×dna)×x17,dx17dt=kasam×x16×x15kdsam×x17(k1dmdm2+k2dmdm2×dna)×x17kdarf×x17,dx18dt=k1se+k2se×x19kde×x18kasp21e×x21×x18+(kdsp21e+kdp21)×x22,dx19dt=kasre×rb×x19+kdsre×rbe2f,dx20dt=kprb×x18×rbjprb+rbkdprb×pp×x20jdprb+x20,dx21dt=k1sp21+k2sp21×x12+k3sp21×x13kdp21×x21kasp21e×x21×x18+(kdsp21e+kde)×x22,dx22dt=kasp21e×x21×x18kdsp21e×x22kde×x22kdp21×x22,
    dx23dt=ksbaxkdbax×x23(kf1+kf2×x26)×x23+kb1×x24+kb2×x28,dx24dt=(kf1+kf2×x26)×x23kb1×x24kdbaxm×x24kasml×x24×x25+kdsml×x28,dx25dt=ksbcl2kdbcl2×x25kasml×x24×x25+kdsml×x28kas3acl×x26×x25,+kds3acl×x29kas3drl×x27×x25+kds3drl×x30+kb2×x28,dx26dt=ks3ackd3ac×x26kas3acl×x26×x25+kds3acl×x29,dx27dt=ks3drkd3dr×x27kas3drl×x27×x25+kds3drl×x30,dx28dt=kasml×x24×x25kdsml×x28kb2×x28kdml×x28,dx29dt=kas3acl×x26×x25kds3acl×x29kd3acl×x29,dx30dt=kas3drl×x27×x25kds3drl×x30kd3drl×x30,dx31dt=k1sp73+k2sp73×x19kdp73×x31,

    where

    p53t=x12+x13, rbe2f1=e2f1tx19, rb=rbtx20rbe2f1

    ks3ac=k1s3ac+k2s3ac×(1+x14)×x13+k3s3ac×x21+k4s3ac×(1+x14)×x12

    ks3dr=k1s3dr+k2s3dr×(1+x14)×x13+k3s3dr×x21+k4s3dr×(1+x14)×x12.

    Parameter values are

    dna=0.1, kin=0.00015, kout=0.0003, dyrt=2

    ksp53=0.5, k1dp53=0.1, kpp53=1, jpp53=0.1

    kdpp53=0.5, jdpp53=0.1, pp=1, k1saspp=0.01

    k2saspp=0.05, kdaspp=0.1, k1smdm2=0.02,

    k2smdm2=0.3, k1dmdm2=0.1, k2dmdm2=0.1, k1sarf=0.01

    k2sarf=0.3, kdarf=0.1, kasam=10, kdsam=2

    k1se=0.01, k2se=0.5, kde=0.12

    kasre=5, kdsre=1, kprb=1, jprb=0.1

    kdprb=0.5, jdprb=0.1, rbt=2, e2f1t=1

    k1sp21=0.03, k2sp21=0.3, k3sp21=0.01, kdp21=0.2

    kasp21e=10, kdsp21e=1, ksbax=0.97, kdbax=0.01

    kf1=0.1, kf2=3.4, kb1=4.2, kb2=4.2

    kdbaxm=0.01, kasml=87, kdsml=0.05, ksbcl2=0.57,

    kdbcl2=0.01, kas3acl=11.14, kds3acl=0.01, kas3drl=11.14

    kds3drl=0.01, kd3ac=0.01, kd3dr=0.01, kdml=0.01

    kd3acl=0.01, kd3drl=0.01, k1sp73=0.01, k2sp73=0.01

    kdp73=0.5, k1s3ac=0.01, k2s3ac=0.05, k3s3ac=0.1

    k4s3ac=0.02, k1s3dr=0.01, k2s3dr=0.1, k3s3dr=0.1, and k4s3dr=0.02.

    The initial values are

    x11=0, x12=0, x13=0.78, x14=0.55

    x15=2.53, x16=2.79, x17=1.99, x18=3.83

    x19=0.9, x20=1.88, x21=1.32, x22=1.25

    x23=95, x24=0, x25=47.49, x26=0

    x27=0, x28=2.26, x29=3.78, x30=3.78, and x31=0.04.



    [1] E. Braverman, D. Kinzebulatov, Nicholson's blowflies equation with distributed delay, Can. Appl. Math. Quart., 14 (2006), 107-128.
    [2] C. Foley, M. C. Mackey, Dynamics hematological disease: a review, J. Math. Biol., 58 (2009), 285-322. doi: 10.1007/s00285-008-0165-3
    [3] K. Gopalsamy, Stability and Oscillation in Delay Differential Equations of Population Dynamics, Springer, Dordrecht, 1992.
    [4] W. S. C. Gurney, S. P. Blythe, R. M. Nisbet, Nicholson's blowflies revisited, Nature, 287 (1980), 17-21. doi: 10.1038/287017a0
    [5] C. Huang, Z. Yang, T. Yi, X. Zou, On the basins of attraction for a class of delay differential equations with non-monotone bistable nonlinearities, J. Diff. Equat., 256 (2014), 2101-2114. doi: 10.1016/j.jde.2013.12.015
    [6] Y. Kuang, Delay Differential Equations: With Application in Population Dynamics, Academic Press, INC, 1993.
    [7] M. C. Mackey, Unified hypothesis for the origin of aplastic anemia and periodic hematopoiesis, Blood, 51 (1978), 941-956. doi: 10.1182/blood.V51.5.941.941
    [8] E. Liz, A. Ruis-Herrera, Delayed population models with Allee effects and exploitation, Math. Biosci. Eng., 12 (2015), 83-97. doi: 10.3934/mbe.2015.12.83
    [9] E. Liz, M. Pinto, V. Tkachenko, S. Tromichuk, A global stability criterion for a family of delayed population models, Quart. Appl. Math., 63 (2005), 56-70. doi: 10.1090/S0033-569X-05-00951-3
    [10] E. Liz, G. Röst, On the global attractor of delay differential equations with unimodal feedback, Disc. Cont. Dyn. Syst., 24 (2009), 1215-1224. doi: 10.3934/dcds.2009.24.1215
    [11] E. Liz, V. Tkachenko, S. Trofimchuk, Mackey-Glass type delay differential equations near the boundary of absolute stability, J. Math. Anal. Appl., 275 (2002), 747-760. doi: 10.1016/S0022-247X(02)00416-X
    [12] E. Liz, V. Tkachenko, S. Tromichuk, A global stability criterion for scalar functional differential equations, SIAM J. Math. Anal., 35 (2003), 596-622. doi: 10.1137/S0036141001399222
    [13] G. Röst, J. Wu, Domain-decomposition method for the global dynamics of delay differential equations with unimodal feedback, Proc. R. Soc. A, 463 (2007), 2655-2669. doi: 10.1098/rspa.2007.1890
    [14] Y. Yuan, J. Belair, Stability and Hopf bifurcation analysis for functional differential equation with distributed delay, SIAM J. Appl. Dyn. Syst., 10 (2011), 551-581. doi: 10.1137/100794493
    [15] Y. Yuan, X. Q. Zhao, Global stability for non monotone delay equations (with application to a model of blood cell production), J. Diff. Equat., 252 (2012), 2189-2209. doi: 10.1016/j.jde.2011.08.026
    [16] N. Bessonov, G. Bocharov, T. M. Touaoula, S. Trofimchuk, V. Volpert, Delay reaction-diffusion equation for infection dynamics, Disc. Cont. Dyn. Syst. Series B, 24 (2019), 2073-2091.
    [17] K. Deng, Y. Wu, On the diffusive Nicholso's blowflies equation with distributed delay, Appl. Math. Lett., 50 (2015) 126-132. doi: 10.1016/j.aml.2015.06.013
    [18] W. Hu, Y. Zhou, Stability of a nonlocal delayed reaction-diffusion equation with a non-monotone bistable nonlinearity, Appl. Math. Lett., 79 (2018), 20-26. doi: 10.1016/j.aml.2017.11.018
    [19] R. H. Martin, H. L. Smith, Abstract functional differential equations and reaction-diffusion systems, Trans. Amer. Math. Soc., 321 (1990), 1-44.
    [20] R. H. Martin, H. L. Smith, Reaction-diffusion systems with time delay: monotonicity, invariance, comparison and convergence, J. Reine Angew. Math., 413 (1991), 1-35.
    [21] T. M. Touaoula, Global dynamics for a class of reaction-diffusion equations with distributed delay and Neumann condition, Commun. Pur. Appl. Anal., 19 (2020), 2473-2490. doi: 10.3934/cpaa.2020108
    [22] T. M. Touaoula, M. N. Frioui, N. Bessonov, V. Volpert, Dynamics of solutions of a reactiondiffusion equation with delayed inhibition, Disc. Cont. Dyn. Syst. Series S, 13 (2020), 2425-2442.
    [23] E. Braverman, S. Zhukovskiy, Absolute and delay-dependent stability of equations with a distributed delay, Disc. Cont. Dyn. Syst., 32 (2012), 2041-2061. doi: 10.3934/dcds.2012.32.2041
    [24] L. Berezansky, E. Braverman, L. Idels, Nicholson's blowflies differential equations revisited: main results and open problems, Appl. Math. Model., 34 (2010), 1405-1417. doi: 10.1016/j.apm.2009.08.027
    [25] L. Berezansky, E. Braverman, L. Idels, Mackey-Glass model of hematopoiesis with non-monotone feedback: stability, oscillation and control, Appl. Math. Compt., 219 (2013), 6268-6283. doi: 10.1016/j.amc.2012.12.043
    [26] H. A. El-Morshedy, Global attractivity in a population model with nonlinear death rate and distributed delays, J. Math. Anal. Appl., 410 (2014), 642-658. doi: 10.1016/j.jmaa.2013.08.060
    [27] I. Gyori, S. Trofimchuk, Global attractivity in x'(t) = -δx(t) + pf (x(t - h)), Dynam. Syst. Appl., 8 (1999), 197-210.
    [28] H. L. Smith, H. R. Thieme, Dynamical Systems and Population Persistence, AMS, Providence, 2011.
    [29] H. L. Smith, Monotone Dynamical Systems: An Introduction to the Theory of Competitive and Cooperative Systems, AMS, Providence, 1995.
    [30] J. Hale, S. M. Verduyn Lunel, Introduction to Functional Differential Equations, Springer-Verlag, New York, 1993.
    [31] T. M. Touaoula, Global stability for a class of functional differential equations (Application to Nicholson's blowflies and Mackey-Glass models), Disc. Cont. Dyn. Syst., 38 (2018), 4391-4419. doi: 10.3934/dcds.2018191
    [32] H. R. Thieme, Mathematics in Population Biology, Princeton University Press, Princeton, 2003.
    [33] T. Yi, Q. Li, Unstable sets, heteroclinic orbits and generic quasi-convergence for essentially strongly order-preserving semiflows, Proc. Edinb. Math. Soc., 52 (2009), 797-807. doi: 10.1017/S0013091506000320
  • This article has been cited by:

    1. Michael J. Lippincott, Jenna Tomkinson, Dave Bunten, Milad Mohammadi, Johanna Kastl, Johannes Knop, Ralf Schwandner, Jiamin Huang, Grant Ongo, Nathaniel Robichaud, Milad Dagher, Andrés Mansilla-Soto, Cynthia Saravia-Estrada, Masafumi Tsuboi, Carla Basualto-Alarcón, Gregory P. Way, Assaf Zaritsky, A morphology and secretome map of pyroptosis, 2025, 36, 1059-1524, 10.1091/mbc.E25-03-0119
  • Reader Comments
  • © 2020 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(4234) PDF downloads(97) Cited by(8)

Figures and Tables

Figures(9)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog