Research article

Cyclic response of a reinforced concrete frame: Comparison of experimental results with different hysteretic models

  • An accurate hysteresis model is fundamental to well capture the non-linearity phenomena occurring in structural and non-structural elements in building structures, that are usually made of reinforced concrete or steel materials. In this sense, this paper aims to numerically estimate through simplified non-linear analyses, the cyclic response of a reinforced concrete frame using different hysteretic models present in the literature. A commercial Finite Element Method package is used to carry out most of the simulations using polygonal hysteretic models and a fiber model, and additionally, a MATLAB script is developed to use a smooth hysteresis model. The experimental data is based on the experiments carried out in the Laboratório Nacional de Engenharia Civil, Portugal. The numerical outcomes are further compared with the experimental result to evaluate the accuracy of the simplified analysis based on the lumped plasticity or plastic hinge method for the reinforced concrete bare frame. Results show that the tetralinear Takeda's model fits closely the experimental hysteresis loops. The fiber model can well capture the hysteresis behavior, though it requires knowledge and expertise on parameter calibration. Sivaselvan and Reinhorn's smooth hysteresis model was able to satisfactorily reproduce the actual non-linear cyclic behavior of the RC frame structure in a global way.

    Citation: Pedro Folhento, Manuel Braz-César, Rui Barros. Cyclic response of a reinforced concrete frame: Comparison of experimental results with different hysteretic models[J]. AIMS Materials Science, 2021, 8(6): 917-931. doi: 10.3934/matersci.2021056

    Related Papers:

    [1] LanJiang Luo, Haihong Liu, Fang Yan . Dynamic behavior of P53-Mdm2-Wip1 gene regulatory network under the influence of time delay and noise. Mathematical Biosciences and Engineering, 2023, 20(2): 2321-2347. doi: 10.3934/mbe.2023109
    [2] Changyong Dai, Haihong Liu, Fang Yan . The role of time delays in P53 gene regulatory network stimulated by growth factor. Mathematical Biosciences and Engineering, 2020, 17(4): 3794-3835. doi: 10.3934/mbe.2020213
    [3] Tingzhe Sun, Dan Mu . Multi-scale modeling identifies the role of p53-Gys2 negative feedback loop in cellular homeostasis. Mathematical Biosciences and Engineering, 2020, 17(4): 3260-3273. doi: 10.3934/mbe.2020186
    [4] Moxuan Zhang, Quan Zhang, Jilin Bai, Zhiming Zhao, Jian Zhang . Transcriptome analysis revealed CENPF associated with glioma prognosis. Mathematical Biosciences and Engineering, 2021, 18(3): 2077-2096. doi: 10.3934/mbe.2021107
    [5] Menghan Chen, Haihong Liu, Ruiqi Wang . Dynamical behaviors of quorum sensing network mediated by combinatorial perturbation. Mathematical Biosciences and Engineering, 2022, 19(5): 4812-4840. doi: 10.3934/mbe.2022225
    [6] Rujing Zhao, Xiulan Lai . Evolutionary analysis of replicator dynamics about anti-cancer combination therapy. Mathematical Biosciences and Engineering, 2023, 20(1): 656-682. doi: 10.3934/mbe.2023030
    [7] Ze-Yu Zhang, Zhu-Jun Mao, Ye-ping Ruan, Xin Zhang . Computational identification of Shenshao Ningxin Yin as an effective treatment for novel coronavirus infection (COVID-19) with myocarditis. Mathematical Biosciences and Engineering, 2022, 19(6): 5772-5792. doi: 10.3934/mbe.2022270
    [8] Shunyi Li . Hopf bifurcation, stability switches and chaos in a prey-predator system with three stage structure and two time delays. Mathematical Biosciences and Engineering, 2019, 16(6): 6934-6961. doi: 10.3934/mbe.2019348
    [9] Xiao-Wei Ye, Xiao-Peng Zhang, Feng Liu . CSB modulates the competition between HIF-1 and p53 upon hypoxia. Mathematical Biosciences and Engineering, 2019, 16(5): 5274-5262. doi: 10.3934/mbe.2019262
    [10] Honghua Bin, Daifeng Duan, Junjie Wei . Bifurcation analysis of a reaction-diffusion-advection predator-prey system with delay. Mathematical Biosciences and Engineering, 2023, 20(7): 12194-12210. doi: 10.3934/mbe.2023543
  • An accurate hysteresis model is fundamental to well capture the non-linearity phenomena occurring in structural and non-structural elements in building structures, that are usually made of reinforced concrete or steel materials. In this sense, this paper aims to numerically estimate through simplified non-linear analyses, the cyclic response of a reinforced concrete frame using different hysteretic models present in the literature. A commercial Finite Element Method package is used to carry out most of the simulations using polygonal hysteretic models and a fiber model, and additionally, a MATLAB script is developed to use a smooth hysteresis model. The experimental data is based on the experiments carried out in the Laboratório Nacional de Engenharia Civil, Portugal. The numerical outcomes are further compared with the experimental result to evaluate the accuracy of the simplified analysis based on the lumped plasticity or plastic hinge method for the reinforced concrete bare frame. Results show that the tetralinear Takeda's model fits closely the experimental hysteresis loops. The fiber model can well capture the hysteresis behavior, though it requires knowledge and expertise on parameter calibration. Sivaselvan and Reinhorn's smooth hysteresis model was able to satisfactorily reproduce the actual non-linear cyclic behavior of the RC frame structure in a global way.



    p53 protein is one of the most important tumor inhibitors, whose mutation is contained in more than 50 percent of human cancer cells[1,2]. It has been proved that the wild-type p53 is able to block the activation of oncogenes in vitro and inhibit the development of tumor cells[3]. Current evidence has suggested that the p53 signaling pathway is one of the promising and effective strategies to treat cancer[4]. In response to stress stimulation or DNA damage, the key functions of p53 protein are the promotion of DNA repair, cell cycle arrest, apoptosis to prevent the damage inherited to daughter cells [1,5]. In addition to activating the transcription of promoters containing homologous binding sites, wild-type p53 can also inhibit the activities of various promoters, and whose expression is positively correlated with the genes of cell proliferation or increased malignant tumor [6]. There are three mechanisms for Mdm2 to inhibit p53. First, activate Mdm2 and trans p53, inhibit the transcriptional activity of p53. Second, the combination of Mdm2 and p53 contributes to the degradation of the proteasome. Third, Mdm2 acts as a p53 ubiquitin ligase (E3), thus it down-regulate the level of p53 protein [2,7,8,9]. The reason that we are concerned about Mdm2 is (a) Mdm2 deletion mutants are mortal in early mouse development, and (b) Mdm2 directly regulates the activity and stability of tumor suppressor p53 [2]. The p53 oscillation induced by undamped IR at the single-cell level consists of three subsystems: DNA damage repair module, ataxia-telangiectasia mutation (ATM) switch, and p53-Mdm2 oscillator [10]. Disrupting the interaction between p53 and Mdm2 is considered as a new strategy for the treatment of cancer without p53 mutation [11]. Low expression often occurs in tumors that carry wild-type p53 [12]. This phenomenon is caused by the over-expression of Mdm2, a major competitor of p53 [12]. For example, it has been observed in breast cancer, brain cancer, and lung cancer that when the p14 gene suppresses Mdm2 by isolating Mdm2 in the nucleus, low expression of p53 may occur in tumors that carry wild-type p53 [12]. Binding with viral proteins in infected cells also leads to low expression of p53 [12]. Tp53 mutates in half of all tumors, which is the most common mutation gene found in human tumors, and this also emphasizes the important anti-tumor effect of p53 [10]. Sometimes p53 is referred to as the gene "Guardian" [5,13].

    Nutlins are potent and selective small-molecule antagonists of Mdm2 that in the p53-binding pocket and activate the p53 pathway in cells with wild-type p53, thereby inhibiting the growth of cells and cause cell apoptosis or cell cycle arrest (including aging and quiescence)[14,15]. It has been proved to be effective against a variety of wild-type p53 tumor cells, including neuroblastoma, osteosarcoma, retinocytoma, and leukemia[15]. Although Nutlin is more effective than RiRNA, this binding affinity cannot change the autoubiquitination effect of Mdm2[4]. Nutlin has been taken as an activator of p53 in preclinical studies currently[11]. Some studies have shown that Nutlin can be activated independently of the p53 phosphorylation[14]. The structure of Nutlin is similar to SN15 peptide that is bound to Mdm2, which is a competitive inhibitor of p53-Mdm2 interaction[16]. Nutlin activates p53 by releasing p53 from negative control mediated by Mdm2, thus making up for the lack of the upstream signal pathway of p53[8]. Nutlin is a p53 activator that can induce G1-S and G2-M detection points to protect cell proliferation, but it cannot protect the cells with mutant p53 from the toxicity of steviol cells[17]. However, stable tetraploid can be isolated to clone the cells treated by Nutlin, and clone diploid with higher drug resistance to apoptotic induced by ionizing radiation and cisplatin[18]. Due to understanding the new uses of existing drugs can provide valuable information for the molecular mechanism and also ensured safety, drug relocation has become a new strategy in drug discovery[16]. Nutlin-3a can show the ability to activate p53 in cell culture at a certain concentration (about 510μM) and inhibit tumor growth when oral administration, which confirms the interaction between Nutlin and p53- Mdm2[19]. In conclusion, Nutlin has a significant anti-tumor effect regardless of the state of the p53 gene, which shows that it brings hope for the treatment of human tumors[15].

    Recently, researchers believe that the dynamic system with time delay has more complex behaviors. For example, the time delay can lead to the loss of stability, various oscillations, and periodic solutions induced by Hopf bifurcation[20]. Scientists have developed a variety of mathematical models, including continuous-time differential equations, discrete-time differential equations, time delays differential equations, and stochastic models[5]. It is worth noting that Lëv Bär-Or et al. established an ordinary differential equations (ODEs) model, in which the negative feedback loop between p53 and Mdm2 is taken into account[5]. If the damage can be repaired, the DNA damage will be gradually repaired with the oscillation of p53. It needs to be aware that the form of oscillation is different between single-cell and multi-cell. Concretely, the oscillation in single-cell is basically limit cycle oscillation, while that in multi-cell is damped oscillation[21]. A four-dimensional model[4] is proposed by Häseeb et al, which is mainly integrated by Hünziker's model[2] and Püszynski's model[12]. However, in fact, both protein synthesis and drug response are time-consuming processes that are neglected in the above model. Here, we focus on the collaborative effect of time delays, drug doses, and important biochemical reaction rates on the p53 drug model.

    To sum up, the innovation of this paper is mainly reflected in four aspects. First of all, the effect of time delays on the integrated model with Nutlin PBK and p53 pathway is considered. Secondly, the collaborative effects of drug dose of Nutlin and time delay on the dynamic behavior of p53 are studied. Thirdly, the stability of positive equilibrium and the Hopf bifurcation is studied theoretically and numerically by choosing time delays as the bifurcation parameters. Fourthly, according to the results of numerical simulation, the direction of Hopf bifurcation and the stability of the period are determined. Through numerical simulation, it is found that the time delay of protein synthesis can determine the period and amplitude of the oscillation. The dose of the drug and the time delay of Nutlin reaction also affect the p53 signal pathway.

    The integrated module of Nutlin PBK dynamics and p53 pathway dynamics is shown in Figure 1. The time delay of gene expression is marked in the diagram, with the red arrow representing degradation, the blue arrow representing generation and the orange arrow representing decomposition. Among them, Mdm2 is the main negative regulator of p53, which can inhibit the activity of p53 that can in turn induce the transcription and translation of Mdm2. In this way a negative feedback loop is formed, which can restrict p53 in an inactive state in the normal cells. Thus blocking-up the negative feedback loop between p53 and Mdm2 is a focus of treat disease. Upon the damage caused by extracellular stress, the activation of p53 plays a vital role by promoting cell cycle arrest or apoptosis of the damaged cell. Nutlin is a specific small molecule antagonist of Mdm2, which can bind to the binding site of Mdm2 and p53, so prevent the binding of Mdm2 and p53. It should be noted that the synthesis of Mdm2 protein is a complex and time-consuming process, which requires a certain amount of time to complete gene expression. Thus, there is always an inevitable time delay between the initiation of transcription and the emergence of fully functional Mdm2 proteins, represented by τ2(transcriptional delay) and τ3(translation delay), respectively. In addition, a certain amount of time is also need for the module to response to the drug Nutlin, which is denoted as τ1.

    Figure 1.  An integrated network diagram of Nutlin PBK and p53 pathway. The red arrows represent degradation, the blue arrows represent formation, and the orange arrows represent dissociation.

    According to the interactions among the components, Häseeb et al propsed the following differential equations (2.1) to design the drug dosage in order to revive p53 activity[4].

    {˙p(t)=σp(t)αkfp(t)m(t)+kbc(t)+c(t)γ,˙mm(t)=ktp2(t)mm(t)β,˙m(t)=ktlmm(t)kfp(t)m(t)+kbc(t)+c(t)δm(t)γka3nm(t),˙c(t)=kfp(t)m(t)kbc(t)c(t)δc(t)γ, (2.1)

    where p(t), mm(t), m(t) and c(t) are denote the concentration of nuclear-p53, Mdm2 protein, Mdm2 mRNA and p53-Mdm2 complex, respectively.

    In this part, based on the model (2.1), the response time delay of Nutlin drug and the protein synthesis time delay of mdm2 are included in our new model. The new model (2.2) can be shown as follows:

    {˙p(t)=σp(t)αkfp(t)m(t)+kbc(t)+c(t)γ,˙mm(t)=ktp2(tτ)mm(t)β,˙m(t)=ktlmm(t)kfp(t)m(t)+kbc(t)+c(t)δm(t)γka3nm(tτ1),˙c(t)=kfp(t)m(t)kbc(t)c(t)δc(t)γ, (2.2)

    The meaning and value of the parameters in our new model (2.2) are shown in Table 1.

    Table 1.  List of the meaning and value of parameters in the system (2.2).
    ParameterDescriptionvalueReference
    σProduction rate of p531000nM h1 [2]
    αMdm2 independent deactivation/degradation of p530.1 h1 [2]
    δMdm2 dependent deactivation/degradation of p5311 h1 [2]
    ktTranscription of Mdm20.03 nM1h1 [2]
    ktlTranslation of Mdm21.4 h1 [2]
    βDegradation rate of Mdm2 mRNA0.6 h1 [2]
    γMdm2 degradation/deactivation0.2 h1 [2]
    kbDissociation of Mdm2-p537200 h1 [2]
    kfDissociation constant of Mdm2-p535000 nM1h1 [2]
    ka3Nutlin-Mdm2 association rate0.00006 sec1 [12]
    τTranslational and transported time delays1025minestimate
    τ1The time delays of the reaction of adding Nutlin1030minestimate
    nThe actual amount of Nutlin present in the cell5 μM[4]

     | Show Table
    DownLoad: CSV

    The ubiquitous delay in protein synthesis and drug response has important consequences on dynamics of gene expression and consequently on downstream biological function[37]. In particular, such delay usually can result in oscillatory behavior and other more complex dynamics[29,30]. However, so far the study of the impact of such delay on the integrated module of Nutlin PBK and p53 pathway has not been reported both experimentally and theoretically. Therefore, we focus on the effect of varying protein synthesis delay and the response delay of drug Nutlin on the stability and amplitude and period of the p53 oscillations. Biologically speaking, the p53 protein is one of the most important tumor inhibitors, with mutations present in more than 50 percent of human cancer cells[1,2]. The oscillation of p53 can promote the repair of damage[31]. If the damage cannot be completely repaired, p53 can also promote programmed cell death[32]. In this way, cells can avoid inaccurate genetic information inherited to the next generation of cells[33,34,35]. In general, there is an average transcriptional delay of 10–20 minutes between the action of a transcription factor on a gene promoter and the presence of the corresponding mature mRNA in the cytoplasm[36]. Similarly, the synthesis of a typical protein from mRNA requires a translation delay of about 1–3 minutes[36]. Therefore, the total Mdm2 protein synthesis time delay τ is estimated about 11–23 minutes, which can be observed in a real experiment. In our manuscript, we set it among 10-25 minutes. In addition, the critical value of Nutlin response time delay is calculated theoretically to be about 24 minutes[19]. In our paper, we estimate the Nutlin response time delay τ1 is about 10–30 minutes that is within the biologically permissible range [14,16,17,18,19].

    p53, as a transcription factor, regulates the expression of the Mdm2 gene, which is inevitably associated with the protein synthesis time delay. Due to the fact that the process of transcription, translation and translocation is slow, time delay inevitably appears in the process of gene regulation. The appearance of time delays often causes great changes in the dynamic properties of the network system. On account of the limitation of signal transmission speed, the effect of time delays on the system is an important factor that has to be considered. In consequence, it is of great significance to analyze the stability of gene regulatory networks with time delays[22,23]. In this part, based on the results of numerical simulation, we will study how the time delays affect the dynamic behavior of the system (2.2). Special attention will be focused on the stability of positive equilibrium and the existence of Hopf bifurcation by analyzing the distribution of the corresponding characteristic equation root. In addition, the properties of Hopf bifurcation are discussed by giving the judgment formulas of three indexes based on the standard form theory and the central manifold theorem for differential equations with time delays developed by Hässard et al. The detailed theoretical calculation can be seen in the Supplementary.

    The numerical simulation will be carried out by using Mathematica 12.0 to verify the theoretical results. The parameters of system (2.2) are based on the previous experimental data[2,12] and estimated on the biochemical constraints, and they are shown in Table 1.

    {˙p(t)=10000.1p(t)5000p(t)m(t)+7200c(t)+0.2c(t),˙mm(t)=0.03p2(tτ)0.6mm(t),˙m(t)=1.4mm(t)5000p(t)m(t)+7200c(t)+11c(t)0.2m(t)0.2165m(tτ1),˙c(t)=5000p(t)m(t)7200c(t)11c(t)0.2c(t). (3.1)

    According to the parameters in Table 1, it is easy to calculate that the system (2.2) has a positive equilibrium E.

    E=(19.537,19.085,6.698,90.731). (3.2)

    To explore the effect of time delays on the dynamic behavior of network integrated by Nutlin PBK and p53 pathways, numerical simulation are performed under different cases of time delays via Mathematica 12.0, as shown in Figures 2 and 3. The two sets of figures show that both the protein synthesis time delay τ and the Nutlin drug response time delay τ1 can lead to oscillation of the Nutlin PBK and p53 pathway integration network. When there are no time delays or the time delays are too little, the system will trend to a stable state. Once the time delays exceed their thresholds, the system will lose stability and begin oscillation. Moreover, as the time delays continue to increase, the oscillation will be more obvious, which also can be observed from Figures 4 and 6. In addition, the amplitude and period increase with the increase of time delay, as shown in Figures 5 and 7.

    Figure 2.  The effect of time delay τ on dynamic behavior of network integrated by Nutlin PBK and p53 pathways at τ1=0h. The blue, red, yellow and green lines respectively represent the concentrations of p53, Mdm2, p53-Mdm2 complex, Mdm2 mRNA. (a) The time histories of the system (2.2) when τ=0h. (b) The time histories of the system (2.2) when τ=0.14h<τ0=0.214h. (c) The time histories of the system (2.2) when τ=0.5h>τ0=0.214h. (d) The time histories of the system (2.2) when τ=0.8h>τ0=0.214h.
    Figure 3.  The effect of time delay τ1 on dynamic behavior of network integrated by Nutlin PBK and p53 pathways at τ=0h. The blue, red, yellow and green lines respectively represent the concentrations of p53, Mdm2, p53-Mdm2 complex, Mdm2 mRNA. (a) The time histories of the system (2.2) when τ1=0h. (b) The time histories of the system (2.2) when τ1=0.3h<τ10=0.402h. (c) The time histories of the system (2.2) when τ1=0.5h>τ10=0.402h. (d) The time histories of the system (2.2) when τ1=0.8h>τ10=0.402h.
    Figure 4.  The effects of different time delays τ on the oscillation amplitude and period of the Nutlin PBK and p53 pathways at τ1=0h. The yellow, purple, blue, red and green lines indicate that the time delay is taken as τ=0,0.15,0.2,0.6,0.8h, respectively. (a) The time histories of the concentration levels of p53. (b) The time histories of the concentration levels of Mdm2. (c) The time histories of the concentration levels of p53-Mdm2 complex.(d) The time histories of the concentration levels of Mdm2 mRNA. The results show that as the time delay τ increases, the amplitude and period of the oscillation also increases.
    Figure 5.  When τ1=0h, the changing trend of the properties of the oscillation as the time delay τ increases. (a) The changing trend of the amplitude as the time delay τ increases. (b) The changing trend of the period as the time delay τ increases.
    Figure 6.  The effects of different time delays τ1 on the oscillation amplitude and period of the Nutlin PBK and 53 pathways at τ=0h. The yellow, purple, blue, red and green lines indicate that the time delay is taken as τ1=0,0.3,0.4,0.6,0.8h. (a) The time histories of the concentration levels of p53. (b) The time histories of the concentration levels of Mdm2. (c) The time histories of the concentration levels of p53-Mdm2 complex. (d) The time histories of the concentration levels of Mdm2 mRNA. The results show that as the time delay τ1 increases, the amplitude and period of the oscillation also increases.
    Figure 7.  When τ=0h, the changing trend of the properties of the oscillation as the time delay τ1 increases. (a) The changing trend of the amplitude as the time delay τ1 increases. (b) The changing trend of the period as the time delay τ1 increases.

    The above simulation results also can be verified. To do this, the first step is to consider the dynamic behavior of the system (2.2) without time delays. Based on the Routh-Hurwitz criterion, it is verified that the positive equilibrium point of the system (2.2) is asymptotically stable when τ=0h and τ1=0h.

    In the second step, there are a pair of pure imaginary roots λ=±ω0i as the characteristic roots of the system (2.2), where ω0=1.94754. Moreover, a set of critical time delays are obtained as follows.

    τj=0.214+0.3.266j,j=0,1,2,. (3.3)

    Define τ0=hour(τj)=0.214. According to the lemma in Supplementary (S1.2), the transversal condition is satisfied. Therefore, the positive balance E is gradually stable at τ<τ0=0.214h, as shown in Figure 2(a), (b). With the increase of time delay, when the time delay τ>τ0=0.214h, the positive equilibrium E would be unstable and Hopf bifurcation would occur, as shown in Figure 2(c), (d).

    The third step is to investigate the effect of time delays on the properties of Hopf bifurcation and the bifurcated oscillations. Combined with the above discussion and the calculation process in the Supplementary S2, we get three indexes to characteristic the properties of the Hopf bifurcation and the bifurcated oscillations.

    μ2=331451.6547490089>0,T2=2.132144488497094×106>0,β2=715228.8602363452<0. (3.4)

    Similarly, for the Nutlin drug response time delay τ1, there are a pair of pure imaginary roots λ=±ω10i as the characteristic roots of the system (2.2), where ω10=2.30206. Moreover, a set of critical values of time delay τ1 are found as follows.

    τj1=0.402+2.729j,j=0,1,2,. (3.5)

    Define τ10=hour(τj1)=0.402 and λ(τ10)=1.071+0.100i. The three indexes are also obtained.

    μ2=0.001168196819386769>0,T2=0.0055866866227185819>0,β2=0.0021632348152077793<0. (3.6)

    According to the theorem in Supplementary (S2.1), μ2 means that the system (2.2) has experienced a supercritical Hopf bifurcation at τ=τ0 or τ1=τ10. β2<0 means that the bifurcated-period solution from the positive equilibrium E=(19.537,19.085,6.698,90.731) is stable in the central prevalence. In addition, T2>0 means that the period of the periodic solution increases as τ or τ1 increases.

    In order to further study the effect of time delays in system (2.2), the dynamic behaviors under two cases including fixed biochemical rates with different time delays and fixed time delays with different biochemical rates are discussed. Here, we mainly focus on the effects of three parameters on p53 oscillations, namely Mdm2 degradation rate γ, Mdm2 dependent p53 degradation rate δ, and the actual dose of Nutlin in the cell n.

    As an effective inhibitor of p53, Mdm2 binds to the transcriptional activation domain of p53, and it also promotes the rapid degradation of p53 in other stable conditions[25]. To investigate the role of Mdm2-dependent degradation of p53 on the integrated network of Nutlin PBK and p53 pathways, the evolution processes of the system (2.2) are simulated under different Mdm2 degradation rates γ with time delay τ=0.8h(τ1=0.8h). As shown in Figure 8 (a)(d) (Figure 10 (a)(d)). When γ=0.01, the system in Figure 8 (a) maintains at a stable state, while the system in Figure 9 (a) is in a progressively stable state. With the increase of γ, oscillation occurs, and the amplitude of oscillation increases with the increase of γ. Due to the different meanings of the two time delays, the amplitude of oscillation is slightly different even under the same time delays and degradation rate. The results showed that Mdm2 degradation rate γ could cause Hopf bifurcation, and only a certain value of Mdm2 degradation rate could drive p53 oscillation. The changing trends of amplitude and period of the bifurcated oscillation are shown in Figures 9 and 11.

    Figure 8.  The evolution of the Nutlin PBK and p53 integrated network with time delays τ=0.8h and τ1=0h. The blue, red, yellow, green lines indicate the concentration levels of p53, Mdm2, p53-Mdm2, and Mdm2 mRNA, respectively. (a)–(d) The evolution processes of system (2.2) when the Mdm2 degradation rates γ=0.01,0.07,0.12,0.2, respectively.
    Figure 9.  When τ1=0h, the changing trend of the oscillation properties as the Mdm2 degradation rates γ increases. (a) The changing trend of the amplitude as the Mdm2 degradation rates γ increases. (b) The changing trend of the period as the Mdm2 degradation rates γ increases.
    Figure 10.  The evolution of the Nutlin PBK and p53 integrated network with time delays τ1=0.8h and τ=0h. The blue, red, yellow, green lines indicate the concentration levels of p53, Mdm2, p53-Mdm2, and Mdm2 mRNA, respectively. (a)–(d) The evolution processes of system (2.2) when the Mdm2 degradation rates γ=0.01,0.07,0.12,0.2, respectively.
    Figure 11.  When τ=0h, the changing trend of the oscillation properties as the Mdm2 degradation rates γ increases. (a) The changing trend of the amplitude as the Mdm2 degradation rates γ increases. (b) The changing trend of the period as the Mdm2 degradation rates γ increases.

    In the absence of time delays, as the degradation rate δ of p53 dependent on Mdm2 decreases, the system is always in an oscillating state[4]. In order to understand the collaborative effect of time delays and the degradation rate of p53 depending on Mdm2, the dynamic behaviors of the system (2.2) under different values of δ are studied. We select three groups of parameter levels, i.e. δ=11, δ=30, δ=60, then carry out the numerical simulation under three situations including without time delay and time delay is before and after the critical value, as shown in Figure 12. From (a)–(c), we can see that the system is always aesthetically stable regardless of δ=11, 30, and 60 when τ=0h. However, in the middle channel (d)–(f) of Figure 12, the time delay τ=0.2h that is less than the critical value. When δ=11, the system (2.2) is asymptotically stable, while δ=30 or δ=60, the system (2.2) is stable. In Figure 12 (g)(i), time delay is set as τ=0.7h that is bigger than the critical value. When δ=11, the system (2.2) begins to oscillate significantly and the amplitude of oscillation increases significantly. When δ=30 and δ=60, the amplitude and period of oscillation of system (2.2) decrease with the increase of degradation rate δ. Similarly, when τ1=0h, the system is always asymptotically stable regardless of δ=11, 30, 60 as shown in Figure 14 (a)(i). In Figure 14 (d)(f), time delay is set as τ1=0.4h that is less than the critical value. When δ=11, the system (2.2) is asymptotically stable. When δ=30 and δ=60, the system (2.2) is stable. However, when the time delay is larger than the critical value, for example τ1=0.7h as shown in Figure 14 (g)(i). When δ=11, the system (2.2) begins to oscillate significantly and the amplitude of oscillation increases significantly. When δ=30, the oscillation of system (2.2) still exists, but along with it, the amplitude is apparently fallen. When δ=60, the system (2.2) is in a state of damped oscillation. In other words, when the time delays are less than the critical value, whether the degradation rate is low, medium or high, the system (2.2) is asymptotically stable. When they are larger than their critical values, the amplitude and period of the oscillation are decreased obviously as the degradation rate increases. The changing trends of amplitude and period of the bifurcated oscillation are shown in Figures 13 and 15.

    Figure 12.  When τ1=0h, the influence of parameter δ and time delay τ on the dynamic behavior of the system (2.2). Blue, red, yellow, green lines indicate the concentration levels of p53, Mdm2, p53-Mdm2, and Mdm2 mRNA, respectively. (a)–(c) The evolution process of system (2.2) when τ=0h and δ takes 11, 30, 60 in turn. (d)–(f) The evolution process of system (2.2) when τ=0.2h and δ takes 11, 30, 60 in turn. (g)–(i) The evolution process of system (2.2) when τ=0.7h and δ takes 11, 30, 60 in turn.
    Figure 13.  When τ1=0h, the changing trend of the oscillation properties as the degradation rate δ of p53 dependent on Mdm2 increases. (a) The changing trend of the amplitude as δ increases. (b) The changing trend of the period as δ increases.
    Figure 14.  When τ=0h, the influence of parameter δ and time delay τ1 on the dynamic behavior of the system (2.2). Blue, red, yellow, green lines indicate the concentration levels of p53, Mdm2, p53-Mdm2, and Mdm2 mRNA, respectively. (a)–(c) The evolution process of system (2.2) when τ1=0h and δ takes 11, 30, 60 in turn. (d)–(f) The evolution process of system (2.2) when τ1=0.2h and δ takes 11, 30, 60 in turn. (g)–(i) The evolution process of system (2.2) when τ1=0.7h and δ takes 11, 30, 60 in turn.
    Figure 15.  When τ=0h, the changing trend of the oscillation properties as the degradation rate δ of p53 dependent on Mdm2 increases. (a) The changing trend of the amplitude as δ increases. (b) The changing trend of the period as δ increases.

    Blocking HDM2 by Nutlin may also restore the level of p53[26]. In order to further study how Nutlin recovers the level of p53 in the presence of time delays, the time evolution processes are displayed by changing the actual amount n of Nutlin present in the cell under two cases of time delay τ=0.2h and 0.5. As shown in Figure 16, (a)(c), when n takes 5, 8, 10 in turn, the system (2.2) is asymptotically stable when τ=0.2 regardless of the values of n. Figure 16 (d)(f) shows the evolution processes of the system (2.2) when τ=0.5h. When n=5, the system (2.2) oscillates. When n=8 and n=10, the amplitude and period of the oscillaiton is decreasing with the increase of n. In other words, when the time delay is less than the critical value, the system is asymptotically stable in spite of the value of n. When the time delay is larger than the critical value, the amplitude and period of the system decrease with the increase of n. As shown in Figure 18 (a)(c), we can see that when τ1=0.38h, the oscillation mode are gradually distinct when n increases from 5 and 8 to 10. Figure 18 (d)–(f) depicted that oscillation always exists when τ1=0.5h, n=5, n=8, n=10, and the amplitude of oscillation increases with the increase of n. The changing trends of amplitude and period of the bifurcated oscillation are shown in Figures 17 and 19.

    Figure 16.  When τ1=0h, the influence of parameter n and time delay τ on the dynamic behavior of the system (2.2). Blue, red, yellow, green lines indicate the concentration levels of p53, Mdm2, p53-Mdm2, and Mdm2 mRNA, respectively. (a)–(c) The evolution process of system (2.2) when τ=0.2h and n takes 5, 8, 10 in turn. (d)–(f) The evolution process of system (2.2) when τ=0.5h and n takes 5, 8, 10 in turn.
    Figure 17.  When τ1=0h, the changing trend of the oscillation properties as the actual amount n of Nutlin present in the cell increases. (a) The changing trend of the amplitude as n increases. (b) The changing trend of the period as n increases.
    Figure 18.  When τ=0h, the influence of parameter n and time delay τ1 on the dynamic behavior of the system (2.2). Blue, red, yellow, green lines indicate the concentration levels of p53, Mdm2, p53-Mdm2, and Mdm2 mRNA, respectively. (a)–(c) The evolution process of system (2.2) when τ1=0.38h and n takes 5, 8, 10 in turn. (d)–(f) The evolution process of system (2.2) when τ1=0.5h and n takes 5, 8, 10 in turn.
    Figure 19.  When τ=0h, the changing trend of the oscillation properties as the actual amount n of Nutlin present in the cell increases. (a) The changing trend of the amplitude as n increases. (b) The changing trend of the period as n increases.

    To sum up, when the time delay τ is less than the critical value, the system is asymptotically stable. When it is larger than the critical value, the oscillation occurs, and the amplitude and period decrease with the increase of n as shown in Figure 16. In contrast, when the time delay τ1 is less than its critical value, oscillation occurs with the increase of n. When it is larger than the critical value, the amplitude of system (2.2) increases with the increase of n. This also indicates that the time delay in protein synthesis, the time delay in Nutlin response, and the actual dose of Nutlin in the cell have effects on the dynamic behavior of the system, so it is reasonable and valuable to introduce the delay into the model.

    In this paper, based on the model proposed by Häseeb et al., a model with the time delay of protein synthesis and the time delay of Nutlin drug response is developed. The stability and oscillation of the system (2.2) with and without time delays are analyzed. The dynamic behavior of the Nutlin PBK and p53 integrated network is studied by Hopf bifurcation theory and numerical simulation. By choosing time delays as bifurcation parameter, we find that they play a key role in inducing Hopf bifurcation and determine the amplitude and period of oscillation. Both the Mdm2 protein synthesis time delay τ and drug Nutlin response time delay τ1 can control the appearance and disappearance of p53 oscillation. Moreover, the mode of the p53 oscillation, including the period and amplitude, is tightly relying on the Mdm2 protein synthesis time delay τ and drug Nutlin response time delay τ. We can infer spontaneously that the p53 oscillation is regulated in a cooperative mode of the two time delays. In fact, changing the number of introns in a gene can bring the changes of protein synthesis time delay[27].

    Biologically speaking, the p53 protein is one of the most important tumor inhibitors, with mutations present in more than 50 percent of human cancer cells [1,2]. The oscillation of p53 can promote the repair of damage [31]. If the damage cannot be completely repaired, p53 can also promote programmed cell death [32]. In this way, cells can avoid inaccurate genetic information inherited to the next generation of cells [33,34,35]. Therefore, cell fate decision can be regulated by the Mdm2 protein synthesis time delay τ and drug Nutlin response time delay τ1 through controlling the p53 oscillation. In addition to taking the time delays as the bifurcation parameter, we also consider the collaborative effect of several important parameters and time delays on the dynamic behavior of the Nutlin PBK and p53 pathways in the integrated network. p53 is a tumor suppressor, and Nutlin is a highly effective selective small molecule antagonist of Mdm2 found in clinical trials in recent years. Nutlin is a p53 activator[28]. The results in this paper can not only help us to understand the dynamic behavior of the Nutlin PBK and p53 integrated networks, but also can provide theoretical evidence to select proper drug Nutlin and set the expected time delays to for the treatment of human tumors.

    The authors express gratitude to the anonymous referee for his/her helpful suggestions and the partial supports of the National Natural Science Foundations of China (12062027/11762022) and Cultivating Plan Program for the Leader in Science and Technology of Yunnan Province (2019HB015).

    The authors declare that they have no conflict of interest.



    [1] Belleri A, Brunesi E, Nascimbene R, et al. (2015) Seismic performance of precast industrial facilities following major earthquakes in the Italian territory. J Perform Constru Fac 29: 04014135. doi: 10.1061/(ASCE)CF.1943-5509.0000617
    [2] Perrone D, Calvi P, Nascimbene R, et al. (2019) Seismic performance of non-structural elements during the 2016 central Italy earthquake. B Earthq Eng 17: 5655–5677. doi: 10.1007/s10518-018-0361-5
    [3] Brunesi E, Peloso S, Pinho R, et al. (2018) Cyclic testing of a full-scale two-storey reinforced precast concrete wall-slab-wall structure. B Earthq Eng 16: 5309–5339. doi: 10.1007/s10518-018-0359-z
    [4] Bianchi F, Nascimbene R, Pavese A (2017) Experimental vs. numerical aimulations: Seismic response of a half scale three-storey infilled RC building strengthened using FRP retrofit. TOCEJ 11: 1158–1169.
    [5] Nascimbene R (2015) Numerical model of a reinforced concrete building: Earthquake analysis and experimental validation. Period Polytech-Civ 59: 521–530. doi: 10.3311/PPci.8247
    [6] Rahnavard R, Rebelo C, Hélder C, et al. (2020) Numerical investigation of the cyclic performance of reinforced concrete frames equipped with a combination of a rubber core and a U-shaped metallic damper. Eng Struct 225: 111307. doi: 10.1016/j.engstruct.2020.111307
    [7] Visintin A (1994) Differential Models of Hysteresis, Germany: Springer-Verlag.
    [8] Sivaselvan M, Reinhorn A (2000) Hysteretic models for deteriorating inelastic structures. J Eng Mech 126: 633–640. doi: 10.1061/(ASCE)0733-9399(2000)126:6(633)
    [9] Pires F (1990) Influência das paredes de alvenaria no comportamento de estruturas reticuladas de betão armado sujeitas a acções horizontais [PhD Thesis]. LNEC, Lisboa.
    [10] Braz-César M, Oliveira D, Barros R (2008) Comparison of cyclic response of reinforced concrete infilled frames with experimental results. The 14th World Conference on Earthquake Engineering, Beijing, China.
    [11] Fardis M, Panagiotakos T (1997) Seismic design and response of bare and masonry-infilled reinforced concrete buildings. Part Ⅰ: Bare structures. J Earthq Eng 1: 219–256.
    [12] Paulay T, Priestley MJN (1992) Seismic Design of Reinforced Concrete and Masonry Buildings, New York: John Wiley & Sons. doi: 10.1002/9780470172841
    [13] Clough RW (1966) Effects of stiffness degradation on earthquake ductility requirement. UCB/SESM 1966/16, University of California, Berkeley, USA.
    [14] Takeda T, Sozen MA, Nielsen NN (1971) Reinforced concrete response to simulated earthquakes. OHBAYASHI-GUMI Technical Research Report 5, Tokyo, Japan.
    [15] Midas Inc. (2004) Analysis manual: Inelastic time history analysis, Korea.
    [16] Natick, Massachusetts, MATLAB R2019a 9.6.0.1072779, USA: The MathWorks, Inc.
    [17] Deng H, Chang Y, Lau D, et al. (2003) A simplified approach for nonlinear response analysis of composite structural members. International Workshop on Steel and Concrete Composite Constructions NCREE, Taiwan, 207–216.
    [18] Kent D, Park T (1971) Flexural members with confined concrete. J Struct Div 97: 1969–1990. doi: 10.1061/JSDEAG.0002957
    [19] Bouc R (1968) Forced vibration of mechanical systems with hysteresis, In: Rupakhety R, Olafsson S, Bessaon B, Proceedings of the Fourth Conference on Non-linear Oscillation, Prague: Academia.
    [20] Wen Y (1976) Method for random vibration of hysteretic systems. J Eng Mech-ASCE 102: 249–263.
    [21] Wen Y (1980) Equivalent linearization for hysteretic system under random excitation. J Appl Mech 47: 150–154. doi: 10.1115/1.3153594
    [22] D'Errico J (2021) fminsearchbnd, fminsearchcon. Available from: https://www.mathworks.com/matlabcentral/fileexchange/8277-fminsearchbnd-fminsearchcon.
  • This article has been cited by:

    1. Amina Yasin, Abdul Rehman Yasin, Muhammad Rizwan Azam, Ali Raza, Robina Nazir, Ridab Adlan Elamin Abdalla, Computational-Model-Based Biopharmaceutics for p53 Pathway Using Modern Control Techniques for Cancer Treatment, 2022, 12, 2076-3417, 5748, 10.3390/app12115748
  • Reader Comments
  • © 2021 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(2922) PDF downloads(227) Cited by(3)

Figures and Tables

Figures(12)  /  Tables(1)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog