A mathematical model for the dynamics of SARS-CoV-2 virus using the Caputo-Fabrizio operator

Tahir Khan; Roman Ullah; Gul Zaman; Jehad Alzabut; Tahir Khan; Roman Ullah; Gul Zaman; Jehad Alzabut

doi:10.3934/mbe.2021305

Mathematical Biosciences and Engineering

2021, Volume 18, Issue 5: 6095-6116. doi: 10.3934/mbe.2021305

Previous Article Next Article

Research article Special Issues

A mathematical model for the dynamics of SARS-CoV-2 virus using the Caputo-Fabrizio operator

1.
Department of Mathematics, University of Malakand Chakdara, Dir (L), Pakhtunkhwa, Pakistan
2.
Department of Computing, Muscat College, Muscat Oman
3.
Department of Mathematics and General Sciences, Prince Sultan, University, Riyadh, Saudi Arabia
4.
Department of Industrial Engineering, OSTIM Technical University, Ankara 06374, Turkey

Academic Editor:Higinio Ramos

Received: 12 May 2021 Accepted: 27 June 2021 Published: 12 July 2021

The pandemic of SARS-CoV-2 virus remains a pressing issue with unpredictable characteristics which spread worldwide through human interactions. The current study is focusing on the investigation and analysis of a fractional-order epidemic model that discusses the temporal dynamics of the SARS-CoV-2 virus in a community. It is well known that symptomatic and asymptomatic individuals have a major effect on the dynamics of the SARS-CoV-2 virus therefore, we divide the total population into susceptible, asymptomatic, symptomatic, and recovered groups of the population. Further, we assume that the vaccine confers permanent immunity because multiple vaccinations have commenced across the globe. The new fractional-order model for the transmission dynamics of SARS-CoV-2 virus is formulated via the Caputo-Fabrizio fractional-order approach with the maintenance of dimension during the process of fractionalization. The theory of fixed point will be used to show that the proposed model possesses a unique solution whereas the well-posedness (bounded-ness and positivity) of the fractional-order model solutions are discussed. The steady states of the model are analyzed and the sensitivity analysis of the basic reproductive number is explored. Moreover to parameterize the model a real data of SARS-CoV-2 virus reported in the Sultanate of Oman from January 1st, 2021 to May 23rd, 2021 are used. We then perform the large scale numerical findings to show the validity of the analytical work.

Keywords:

fractional-order SARS-CoV-2 epidemiological model,
existence and uniqueness,
bounded-ness and positivity,
steady state analysis,
sensitivity analysis,
numerical simulation

Citation: Tahir Khan, Roman Ullah, Gul Zaman, Jehad Alzabut. A mathematical model for the dynamics of SARS-CoV-2 virus using the Caputo-Fabrizio operator[J]. Mathematical Biosciences and Engineering, 2021, 18(5): 6095-6116. doi: 10.3934/mbe.2021305

Related Papers:

[1]	Kexin Li, Hu Chen, Shusen Xie . Error estimate of L1-ADI scheme for two-dimensional multi-term time fractional diffusion equation. Networks and Heterogeneous Media, 2023, 18(4): 1454-1470. doi: 10.3934/nhm.2023064
[2]	Xiongfa Mai, Ciwen Zhu, Libin Liu . An adaptive grid method for a singularly perturbed convection-diffusion equation with a discontinuous convection coefficient. Networks and Heterogeneous Media, 2023, 18(4): 1528-1538. doi: 10.3934/nhm.2023067
[3]	Huanhuan Li, Meiling Ding, Xianbing Luo, Shuwen Xiang . Convergence analysis of finite element approximations for a nonlinear second order hyperbolic optimal control problems. Networks and Heterogeneous Media, 2024, 19(2): 842-866. doi: 10.3934/nhm.2024038
[4]	Jinhu Zhao . Natural convection flow and heat transfer of generalized Maxwell fluid with distributed order time fractional derivatives embedded in the porous medium. Networks and Heterogeneous Media, 2024, 19(2): 753-770. doi: 10.3934/nhm.2024034
[5]	Panpan Xu, Yongbin Ge, Lin Zhang . High-order finite difference approximation of the Keller-Segel model with additional self- and cross-diffusion terms and a logistic source. Networks and Heterogeneous Media, 2023, 18(4): 1471-1492. doi: 10.3934/nhm.2023065
[6]	Zhe Pu, Maohua Ran, Hong Luo . Effective difference methods for solving the variable coefficient fourth-order fractional sub-diffusion equations. Networks and Heterogeneous Media, 2023, 18(1): 291-309. doi: 10.3934/nhm.2023011
[7]	Patrick Henning, Mario Ohlberger . The heterogeneous multiscale finite element method for advection-diffusion problems with rapidly oscillating coefficients and large expected drift. Networks and Heterogeneous Media, 2010, 5(4): 711-744. doi: 10.3934/nhm.2010.5.711
[8]	Yinlin Ye, Hongtao Fan, Yajing Li, Ao Huang, Weiheng He . An artificial neural network approach for a class of time-fractional diffusion and diffusion-wave equations. Networks and Heterogeneous Media, 2023, 18(3): 1083-1104. doi: 10.3934/nhm.2023047
[9]	Alexander Mielke, Sina Reichelt, Marita Thomas . Two-scale homogenization of nonlinear reaction-diffusion systems with slow diffusion. Networks and Heterogeneous Media, 2014, 9(2): 353-382. doi: 10.3934/nhm.2014.9.353
[10]	Yimamu Maimaiti, Zunyou Lv, Ahmadjan Muhammadhaji, Wang Zhang . Analyzing vegetation pattern formation through a time-ordered fractional vegetation-sand model: A spatiotemporal dynamic approach. Networks and Heterogeneous Media, 2024, 19(3): 1286-1308. doi: 10.3934/nhm.2024055

Abstract

1. Introduction

Biological systems are permeated with deterministic laws and randomness. Noise are inevitable at every level of biology, from the most essential molecular, sub-cellular processes to the kinetics of tissues, organs, organisms and populations ^[1]. According to different origination of noise, noise can be divided into two terms: intrinsic and extrinsic noise ^[2,3,4,5]. Intrinsic noise originates from the stochasticity of biochemical events, while extrinsic noise is caused by status differences in individual cells, e.g., ribosomes and metabolites, the concentration of RNA polymerases or local environmental conditions ^[6]. Some experimental and theoretical studies have shown that noise may lead to genetic diseases, gene damage, etc. ^[7,8]. However, over the past few decades, it has been revealed that noise may also play a positive role in various dynamic systems and dynamic processes ^[9,10], for example, stochastic transitions ^[11] and switch processes ^[3,12].

Experimental studies have exposed that fluctuations occur in various stages of gene expression, including transcription, degradation, translation, binding, and so on ^[6,13,14] and noise affects gene expression ^[15]. Smolen et al. ^[16] proposed a simple model of gene transcription system dynamics, which is a bistable model with a positive feedback loop. They pointed out that the simple model manifested multiple stable state and brief perturbations could switch the model between those states. Based on this model, Liu et al. showed that the switch process could be successfully induced by the fluctuation of the rate of degradation and the synthesis of transcription factors ^[17]. The state transition of a rapid response to environmental changes is a common and important phenomenon in biological systems. Noise can also be used to regulate the dynamics of bistable genetic regulatory systems ^[18]. Wang et al. ^[19] further considered the effect of a time delay on the switch process. The related studies have focused on white Gaussian noise and colored noise ^{[3,17,19,20,21,22,23,24]}. However, the range of the Gaussian noise is unbounded and large values can be obtained, but the nature of real physical quantities is always bounded. The use of Gaussian noise may lead to unreasonable results in some cases ^[25,26]. Therefore, researchers have explored the modelling of bounded noise in dynamic systems ^[27]. Among them, sine-Wiener noise is the most widely used form of bounded noise ^[28] that is generated by a sinusoidal function with a constant amplitude, constant frequency and random change of phase with the Wiener process. Sinusoidal functions are bounded, resulting in their noise values not exceeding a fixed amplitude. In recent years, the influence of cross-correlated sine-Wiener (CCSW) noises on nonlinear dynamic systems has attracted wide attention ^[29,30]. Furthermore, our previous studies ^{[31,32,33,34]} revealed temporal coherent resonance phenomenon induced by CCSW noises in regular and small-world neuronal networks and indicated the occurrence of CCSW noise-induced coherent resonance. The increasing degree of correlation of CCSW noises may enhance or impair temporal regularity, which is dependent on the cross-correlation time ^[34]. CCSW noises in the gene transcriptional regulatory system has not been reported to date. Therefore, this study focuses on the influence of CCSW noises on the switch process in the gene transcriptional regulatory system.

2. Model and algorithm

Smolen et al. ^[16] have developed a simple kinetic model to examine the complex dynamic activity of genetic regulatory systems. The dynamics of the system are governed by a single ordinary differential equation to determine the concentration of the transcription factor activator (TF-A):

$\frac{{dx}}{{dt}} = \frac{{{k_f}{x^2}}}{{{x^2} + {K_d}}} - {k_d}x + {R_{bas}},$

(1)

where x(t) denotes the concentration of TF-A. k_d and R_bas represent the decomposition rate and basic synthesis rate of TF-A, respectively. K_d describes the dissociation concentration of the TF-A dimer from the responsive elements, TF-Res. Under the following parameters:

${\left[ { - {{\left( {\frac{{{k_f} + {R_{bas}}}}{{3{k_d}}}} \right)}^3} + \frac{{{K_d}\left( {{k_f} + {R_{bas}}} \right)}}{{6{k_d}}} - \frac{{{K_d}{R_{bas}}}}{{2{k_d}}}} \right]^2} + {\left[ {\frac{{{K_d}}}{3} - {{\left( {\frac{{{k_f} + {R_{bas}}}}{{3{k_d}}}} \right)}^2}} \right]^3} \lt 0,$

(2)

the potential

${U_0}{\rm{(}}x{\rm{)}} = {k_f}\sqrt {{K_d}} {\rm{arctan}}\frac{x}{{\sqrt {{K_d}} }} + \frac{{{k_d}}}{2}{x^2} - {\rm{(}}{R_{bas}}{\rm{ + }}{k_f}{\rm{)}}x.$

(3)

Figure 1(a) shows that when K_d = 10, k_d = 1, R_bas = 0.4, and k_f is in the interval [5.45, 6.68], each value of k_f corresponds to three TF-A, representing two steady states and one unstable state. The corresponding bistable potential U₀ is shown in Figure 1b. When k_f = 6, the stable steady states are x_- ≈ 0.62685 and x₊ ≈ 4.28343 and the unstable steady state is x_u ≈ 1.48971, as shown in Figure 1. To simulate the stochastic effects of the biochemical reaction rates R_bas and k_d, it is assumed that the stochasticity is added to the reaction rates as k_d→k_d + η₁(t) and R_bas→R_bas +η₂(t). Here, multiplicative and additive cross-correlated sine-Wiener noises η₁(t) and η₂(t) are introduced into the original deterministic kinetics. Thus,

$\frac{{dx}}{{dt}} = \frac{{{k_f}{x^2}}}{{{x^2} + {K_d}}} - \left[ {{k_d} + {\eta _1}\left( t \right)} \right]x + {R_{bas}} + {\eta _2}\left( t \right).$

(4)

Figure 1. (a) Bifurcation plot for the steady state of TF-A on the control parameter of the transcription rate k_f. (b) The bistable potential of Eq (3). k_f = 6, K_d = 10, k_d = 1, and R_bas = 0.4.

DownLoad: Full-Size Img PowerPoint

The explicit definitions of CCSW noises are represented as follows:

${\eta _1}\left( t \right) = A\sin \left( {\sqrt {\frac{2}{{{\tau _1}}}} {\omega _1}\left( t \right)}\right) ,$

(5)

${\eta _2}\left( t \right) = B\sin \left( {\sqrt {\frac{2}{{{\tau _2}}}} {\omega _2}\left( t \right)} \right).$

(6)

The noise strength and self-correlation times of η₁(t) and η₂(t) are indicated as A, B, τ₁ and τ₂. ω₁ and ω₂ are two cross-correlated standard Wiener processes, i.e., $\left\langle {{\omega _1}\left(t \right){\omega _2}\left({t'} \right)} \right\rangle = \left\langle {{\omega _1}\left({t'} \right){\omega _2}\left(t \right)} \right\rangle = \lambda \cdot \min \left({t, t'} \right)$ . $\left\langle { \cdot \cdot \cdot } \right\rangle$ denotes an ensemble average and $\min \left({t, t'} \right)$ denotes taking the smaller value between t and $t'$ . As previously mentioned, the statistical properties of ${\eta _1}\left(t \right)$ and ${\eta _2}\left(t \right)$ are given by:

$\left\langle {{\eta _1}\left( t \right)} \right\rangle = \left\langle {{\eta _2}\left( t \right)} \right\rangle = 0,$

(7)

$\left\langle {{\eta _1}\left( t \right){\eta _1}\left( {t'} \right)} \right\rangle = \frac{{A_{}^2}}{2}\exp \left( { - \frac{{t - t'}}{{{\tau _1}}}} \right)\left[ {1 - \exp \left( { - \frac{{4t'}}{{{\tau _1}}}} \right)} \right],$

(8)

$\left\langle {{\eta _2}\left( t \right){\eta _2}\left( {t'} \right)} \right\rangle = \frac{{B_{}^2}}{2}\exp \left( { - \frac{{t - t'}}{{{\tau _2}}}} \right)\left[ {1 - \exp \left( { - \frac{{4t'}}{{{\tau _2}}}} \right)} \right],$

(9)

$\left\langle {{\eta _1}\left( t \right){\eta _2}\left( {t'} \right)} \right\rangle = \left\langle {{\eta _1}\left( {t'} \right){\eta _2}\left( t \right)} \right\rangle = \lambda \frac{{AB}}{2}\exp \left( { - \frac{{t - t'}}{\tau }} \right)\left[ {1 - \exp \left( { - \frac{{4t'}}{\tau }} \right)} \right], \;t \gt t', {\tau _1} = {\tau _2} = \tau.$

(10)

Here, $\lambda$ and τ stand for the cross-correlation intensity (0 ≤ $\lambda$ ≤ 1) and cross-correlation time between η₁(t) and η₂(t), respectively. The coupling between CCSW noises results in a difficulty of directly integrating the dynamic equation (4). However, CCSW noises can be decoupled and transformed as below ^[32,35,36]:

${\xi _1}\left( t \right) = A\sin \left( {\sqrt {\frac{2}{\tau }} {\omega _\alpha }\left( t \right)} \right),$

(11)

${\xi _2}\left( t \right) = B\lambda \sin \left( {\sqrt {\frac{2}{\tau }} {\omega _\alpha }\left( t \right)} \right) + B\sqrt {1 - {\lambda ^2}} \sin \left( {\sqrt {\frac{2}{\tau }} {\omega _\beta }\left( t \right)} \right).$

(12)

Here, ω_α and ω_β denote two independent standard Wiener processes. Note that the transformations do not change the statistical properties of Eqs 7–10. After the above-mentioned decoupling operations, the dynamic equation (4) is rewritten as follows:

$\frac{{dx}}{{dt}} = \frac{{{k_f}{x^2}}}{{{x^2} + {K_d}}} - \left[ {{k_d} + {\xi _1}\left( t \right)} \right]x + {R_{bas}} + {\xi _2}\left( t \right).$

(13)

According to the Euler forward procedure, Eq (13) is discretized in our numerical simulations as follows ^[37]:

$\begin{gathered} x\left( {t + \Delta t} \right) = x\left( t \right) + \left[ {\frac{{{k_f}{x^2}}}{{{x^2} + {K_d}}} - {k_d}x + {R_{bas}}} \right]\;\Delta t \\ - Ax\left( t \right){X_1}\left( t \right)\Delta t + B\left[ {\lambda {X_1}\left( t \right) + \sqrt {1 - {\lambda ^2}} {X_2}\left( t \right)} \right]\;\Delta t \\ - \frac{1}{2}A{X_1}\left( t \right)\;\left\{ {\; - Ax\left( t \right){X_1}\left( t \right) + B\left[ {\lambda {X_1}\left( t \right) + \sqrt {1 - {\lambda ^2}} {X_2}\left( t \right)} \right]\;} \right\}\;\Delta {t^2}, \\ \end{gathered}$

(14)

with

${X_1}\left( t \right) = \sin \left( {\sqrt {{2 \mathord{\left/ {\vphantom {2 \tau }} \right. } \tau }} {\omega _\alpha }\left( t \right)} \right),$

(15)

${X_2}\left( t \right) = \sin \left( {\sqrt {{2 \mathord{\left/ {\vphantom {2 \tau }} \right. } \tau }} {\omega _\beta }\left( t \right)} \right),$

(16)

${\omega _\alpha }\left( t \right) = {\omega _\alpha }\left( {t - \Delta t} \right) + \sqrt { - 2\Delta t\ln {\chi _1}} \cos \left( {2{\rm{ \mathsf{ π} }}{\chi _2}} \right),$

(17)

${\omega _\beta }\left( t \right) = {\omega _\beta }\left( {t - \Delta t} \right) + \sqrt { - 2\Delta t\ln {\chi _3}} \cos \left( {2{\rm{ \mathsf{ π} }}{\chi _4}} \right),$

(18)

Here, ${\chi _1}$ , ${\chi _2}$ , ${\chi _3}$ and ${\chi _4}$ are four independent random numbers uniformly distributed on the unit interval. Δt = 0.001 is the time step.

3. Results and discussion

3.1. Steady state characteristics

We first investigated the effect of the cross-correlation time τ on the steady-state probability distribution (SPD) and switching process. For a given time series of a concentration of TF-A, the total number of time points is N. The observed number n_i of x within the distance bin from x_i to x_i +^△x is counted for a long x (t) sequence that discards transient processing, 0 < x < 10, ^△x = 0.1. The probability n_i/N of each distance is calculated, and thus, the steady-state probability distribution (SPD) is obtained. As shown in Figure 2, the novel phenomenon is that the high concentration state of SPD undergoes a transition from unimodal to bimodal, and therefore, the overall SPD has a three-peak structure with a larger noise cross-correlation time τ. With the increase of the noise cross-correlation time τ, the peak of the probability distribution in the low steady state increases, which corresponds to a decrease of the TF-A monomer concentration and the occurrence of a switch process (i.e., "on"→"off"). That is, the cross-correlation time τ not only causes the on-off switching process but also leads to the formation of double peaks of a high steady state and a triple peaks structure in the overall SPD. Tristability itself may be valuable; for example, a genetic tristate can be found or used in various biochemical networks. This triple peaks structure has only ever been observed in an experiment in synthetic biology by Rui Ma et al. ^[38]. They noted that the triple steady state phenomenon resulted from the discrete and fluctuating nature of a small system. Stem cells and other kinds of multipotent cells might differentiate into two states, and then, one of them redifferentiates into three states according to Waddington's epigenetic landscape ^[39]. Our work also provides a possible mechanism for the formation of the triple peaks structure.

Figure 2. (a) Probability distribution of x(t) for different cross-correlation time τ with A = 0.1, B = 0.1 and λ = 0.5. (b). Sample paths of x(t) with τ = 0.5 and τ = 20. The other parameter values are the same as those in Figure 1.

DownLoad: Full-Size Img PowerPoint

Based on the different effects of the noise cross-correlation time τ on the system, the influence of the other parameters of noise is considered when different values of the cross-correlation time τ are selected. First, we consider the effect of multiplicative noise, which is shown in Figure 3. It can be seen that the increase of the multiplicative noise intensity A leads the system to be in a low steady state, corresponding to the decrease of the TF-A monomer concentration and the switch to an off position when τ is smaller and larger. However, multiplicative noise has a strong influence on the bimodal structure of a high steady state. With the increase of the multiplicative noise intensity A, the distance between the two peaks increases. Finally, the left peak of the two peaks weakens and disappears, and the whole peak structure returns to the double peak structure. The extreme value of the peak of the high steady state deviates from the steady-state point to the right.

Figure 3. Probability distribution and sample paths of x(t) for different multiplicative noise intensities A. The additive noise intensity B = 0.2, the cross-correlation intensity λ = 0.3, the cross-correlation intensity τ = 0.1 and τ = 20. The other parameter values are the same as those in Figure 1.

DownLoad: Full-Size Img PowerPoint

Then, we investigated the effect of additive noise, which is shown in Figure 4. When τ is small, the increase of the additive noise intensity B leads to a decrease of the low concentration and an increase of the high concentration (Figure 4a). Correspondingly, the time series diagrams of the TF-A monomer concentrations for additive noise intensities B = 0.1 and B = 0.6 are given in Figure 4b. It can be seen that the increase of the additive noise intensity B drives the system to the high steady state, which corresponds to an increase of the TF-A monomer concentration, and the system switches to the "on" position. However, when τ is large, the distance between the two peaks decreases, and they become a single peak with the increase of additive noise B, as shown in Figure 4c-d. At the same time, switching of the process from "on" to "off" is induced.

Figure 4. Probability distribution and sample paths of x(t) for different additive noise intensities B. (a-b) The multiplicative noise intensity A = 0.6, the cross-correlation intensity λ = 0.1, the cross-correlation intensity τ = 0.1. (c-d) the additive noise intensity A = 0.1, the cross-correlation intensity λ = 0.3, the cross-correlation intensity τ = 20. The other parameter values are the same as those in Figure 1.

DownLoad: Full-Size Img PowerPoint

Figure 5a shows the effect of the cross-correlation intensity on SPD when τ is small. It can be found that an increase of the noise cross-correlation intensity λ causes the SPD to shift from a high steady state to a low concentration state. However, it is interesting to note that the shift between the different concentration states does not mean the process switches from "on" → "off". As can be seen from Figure 5b, the mean first passage time (MFPT), which is the mean time spent in one state before switching to another, monotonously increases with the increase of the cross-correlation intensity λ. The mean first passage time reflects the switching time between on and off states. That is, the smaller MFPT is, the more frequent the jump between the two steady states. On the contrary, the longer MFPT is, the longer the system maintains the current state. The corresponding time series of the TF-A monomer concentration, as shown in Figure 5c, indicates that a large cross-correlation intensity λ can maintain the TF-A monomer in the current concentration state, that is, keep in the "on" or "off" state. When τ is large, the increase of the cross-correlation intensity λ makes the high-concentration binary states more distinct (Figure 5d-e) and also causes the process to switch from "off" to "on". Contrary to the case of a smaller τ, for a large cross-correlation intensity λ, the system cannot maintain the current state when τ is larger.

Figure 5. (a, d) Probability distribution, (b) MFPT and (c, e) sample paths of x(t) for different cross-correlation intensities λ. (a) A = 0.6, B = 0.2, τ = 0.1. (d) A = 0.1, B = 0.1 and τ = 10. The other parameter values are the same as those in Figure 1.

DownLoad: Full-Size Img PowerPoint

In addition, there is feedback delay in the real system, so equation (1) becomes:

$\frac{{dx\left( t \right)}}{{dt}} = \frac{{{k_f}x{{\left( {t - {t_d}} \right)}^2}}}{{x{{\left( {t - {t_d}} \right)}^2} + {K_d}}} - \left[ {{k_d} + {\eta _1}\left( t \right)} \right]x\left( t \right) + {R_{bas}} + {\eta _2}\left( t \right),$

(19)

${t_d}$ is the delay time, and the other parameters are the same as in Eq (1).

Figure 6 shows that the feedback time delay plays an important role in random systems. With the increase of the delay time, the peak of the probability distribution at a low steady state decreases, while those of the high steady state increase, corresponding to an increase of the TF-A monomer concentration, and the system switches from "off" to "on".

Figure 6. Probability distribution and sample paths of x(t) for different delay times t_d. A = 0.8, B = 0.1, λ = 0.5 and τ = 0.1. The other parameter values are the same as those in Figure 1.

DownLoad: Full-Size Img PowerPoint

3.2. Mean first passage time (MFPT)

When the system is stochastic bistable, the quantity of interest is the time from one state to another. This time is a random variable, often referred to as the first passage time. Here, we consider the mean first passage time (MFPT) ^[40]. We record the time of the process x(t) to reach the low concentration state x_- from the initial condition x (t = 0) = x₊ (the high concentration state) for the first time and repeat the process 1000 times to obtain the mean first pass time (MFPT). The MFPT can quantify the effect of noise on the transition between two steady states. The main results are shown in Figures 7-10.

Figure 7. The MFPT as a function of the cross-correlation τ (a) for different cross-correlation intensities, λ = 0.1, 0.3, 0.5, 0.7, 0.9 with A = 0.6, B = 0.2; (b) for different multiplicative noise intensities, A = 0.6, 0.7, 0.8, 0.9 with B = 0.2, λ = 0.5; (c) for different additive noise intensities, B = 0.1, 0.3, 0.5, 0.7, 0.9 with A = 0.6, λ = 0.5.

DownLoad: Full-Size Img PowerPoint

Figure 8. The MFPT as a function of the multiplicative noise intensity A. (a) for different cross-correlation intensities, λ = 0.1, 0.3, 0.5, 0.7, 0.9 with B = 0.2 and τ = 0.1; (b) for different additive noise intensities, B = 0.1, 0.3, 0.5, 0.7, 0.9 with λ = 0.5 and τ = 0.1.

DownLoad: Full-Size Img PowerPoint

Figure 9. The MFPT as a function of the additive noise intensity B; (a) for different cross-correlation intensities, λ = 0.1, 0.3, 0.5, 0.7, 0.9 with A = 0.6 and τ = 0.1; (b) for different multiplicative noise intensities, A = 0.6, 0.7, 0.8, 0.9 with λ = 0.5 and τ = 0.1.

DownLoad: Full-Size Img PowerPoint

Figure 10. The MFPT as a function of the cross-correlation intensity λ; (a) for different additive noise intensities, B = 0.1, 0.3, 0.5, 0.7, 0.9 with A = 0.6 and τ = 0.1; (b) for different multiplicative noise intensities, A = 0.6, 0.7, 0.8, 0.9 with λ = 0.5 and τ = 0.1.

DownLoad: Full-Size Img PowerPoint

The MFPT as a function of the cross-correlation time τ for different cross-correlation intensities λ, multiplicative noise intensities A and additive noise intensities B is shown in Figure 7. It can be seen that the MFPT always decreases monotonously with the increase of the cross-correlation time τ. These results indicate that the switch of the process from "on" to "off" becomes easier when the cross-correlation time τ is larger. In addition, the MFPT increases with the increase of the cross-correlation intensity λ, as illustrated in Figures 7a and 10, while the MFPT decreases with the increase of the multiplicative noise intensity A, as shown in Figures 7b and 8. Then, the MFPT varies nonmonotonically with the increase of the additive noise intensity B when λ = 0.5, as shown in Figures 7c, 9 and 10a.

Figure 8 shows the MFPT as a function of multiplicative noise A for different cross-correlation intensities λ and additive noise intensities B. The MFPT always decreases monotonously with the increase of the multiplicative noise intensity A, which is similar to the effect of the cross-correlation time τ. In addition, there are both monotonic and non-monotonic changes with the increase of λ or B when A is fixed in the MFPT as shown in Figure 8a, b. Further results are presented in Figures 9 and 10.

Then, the MFPT as a function of the additive noise intensity B for different cross-correlation intensities and multiplicative noise intensities is shown in Figure 9. The MFPT slightly decreases when the cross-correlation intensity λ is small (λ = 0.1). Moreover, the MFPT first increases and then decreases when λ reaches a moderate value (λ = 0.5). Furthermore, the MFPT increases monotonously and rapidly with a large λ (λ = 0.9). The MFPT showed different trends with the enhancement of B for different values of λ (monotonic decrease, non-monotonic, monotonic increase), which is also shown in Figure 10a. The larger λ is, the more TF-A is able to maintain its current state and the larger MFPT is (Figure 5b). Therefore, it can be inferred that the monotonous decrease of the black dotted line in Figure 9a is due to the role of B when λ is small. This result suggests that the increasing additive noise intensity B also makes the switching process easier. When the value of λ is large, the rapid increase of the MFPT is caused by the enhancement of the ability of λ to maintain the current state. Therefore, when the value of λ is set moderately, the enhancement of B reduces the MFPT, while a moderate value of λ allows the MFPT to be increased, which results in an increase and then a decrease of the MFPT. Moreover, Figure 9b shows that the above nonmonotonicity is destroyed with the increase of A when the value of λ is moderate (λ = 0.5), and the MFPT then undergoes a monotonous increase.

Figure 10 shows the MFPT as a function of the cross-correlation intensity λ for different additive noise intensities B and multiplicative noise intensities A. The MFPT always increases monotonously when an increase of the cross-correlation intensity λ occurs. The results imply that the switch of the process from "on" to "off" becomes difficult. The MFPT decreases monotonically with the increase of B in the region of 0 < λ < 0.1. The MFPT increases monotonically with the increase of B in the region of 0.6 < λ < 1 and presents non-monotonicity in the region of 0.1 < λ < 0.6 (Figure 10a). In Figure 10b, when A increases, the trend of the MFPT increases slowly with the increase of λ. That is, the role of A is opposite to that of λ.

4. Conclusion

In this study, we investigated switch processes in a transcriptional regulatory system driven by CCSW noises. It was found that the increase of the multiplicative noise intensity A and cross-correlation time τ could reduce the concentration of the TF-A monomer and switch to an "off" state. In addition, the effect of the cross-correlation time on SPD presents a novel phenomenon. When the cross-correlation time τ is small, the enhancement of the cross-correlation intensity λ presents a seemingly contradictory result, that is, an increase of the MFPT and an increase of the low steady state, which indicates that the role of the cross-correlation intensity λ is to maintain the current state. The combined effects of the additive noise intensity B and cross-correlation intensity λ show that when λ is small, the effect of maintaining the current state is weak. The enhancement of B leads the SPD to shift from a high steady state to a low concentration one, a switch to the "on" state, and a monotonous decrease of MFPT. Next, when λ is large, the ability to maintain the current state is enhanced and the MFPT increases rapidly. Third, when λ is moderate, the enhancement of B decreases the MFPT and the increase of cross-correlation intensity λ increases the MFPT, which results in the non-monotonicity of the MFPT (first increases and then decreases).

When the cross-correlation time τ is large, the single peak of the high steady state splits into two peaks and the SPD has a three-peak structure. The effect of the cross-correlation intensity λ is different from that of the small cross-correlation time τ. The enhancement of the cross-correlation intensity λ leads to a decrease of the low steady state and an increase of the high steady state, and the bimodal structure of the high steady state is sharper. The SPD is shifted from a high steady state to a low one when the additive noise intensity is decreased. However, the effect of multiplicative noise and additive noise is opposite to that of the double peaks in the high concentration state. The increase of the multiplicative noise intensity leads to an increase of the distance between two peaks in the high concentration state, while the increase of the intensity of additive noise presents an inverse process, that is, the distance between the two peaks decreases until they transition into a single peak. The above results show different characteristics from white Gaussian noise, which lays a foundation for further research on the influence of CCSW noises on nonlinear systems.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant No.11675060, 91730301, 11804111), the Huazhong Agricultural University Scientific and the Fundamental Research Funds for the Central University (Grant No. 26622018JC017).

Conflict of Interest

The authors declare that they have no competing interests.

References

[1]	S. G. Samko, Fractional Integrals and Derivatives, Theory and Applications, Minsk, Nauka I Tekhnika, 1987.
[2]	I. Podlubny, Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications, Elsevier, 1999.
[3]	D. Baleanu, Z. B. Güvenç, J. T. Machado, New Trends in Nanotechnology and Fractional Calculus Applications, Springer, New York, 2010.
[4]	D. Baleanu, J. A. T. Machado, A. C. Luo, Fractional Dynamics and Control, Springer Science and Business Media, 2011.
[5]	G. Zaman, Y. H. Kang, I. H. Jung, Stability analysis and optimal vaccination of an SIR epidemic model, BioSystems, 93 (2008), 240-249. doi: 10.1016/j.biosystems.2008.05.004
[6]	Y. Wang, J. Cao, Global dynamics of a network epidemic model for waterborne diseases spread, Appl. Math. Comput., 237 (2014), 474–488.
[7]	H. Abboubakar, J. C. Kamgang, D. Tieudjo, Backward bifurcation and control in transmission dynamics of arboviral diseases, Math. Biosci., 278 (2016), 100-129. doi: 10.1016/j.mbs.2016.06.002
[8]	T. Khan, G. Zaman, M. I. Chohan, The transmission dynamic and optimal control of acute and chronic hepatitis B, J. Biol. Dyn., 11 (2017), 172-189. doi: 10.1080/17513758.2016.1256441
[9]	Y. Tu, S. Gao, Y. Liu, D. Chen, Y. Xu, Transmission dynamics and optimal control of stage-structured HLB model, J. Biol. Dyn., 16 (2019), 5180.
[10]	A. Akgül, S. H. Khoshnaw, A. S. Abdalrahman, Mathematical modeling for enzyme inhibitors with slow and fast subsystems, Arab J. Basic Appl. Sci., 27 (2020), 442–449. doi: 10.1080/25765299.2020.1844369
[11]	A. Atangana, Modelling the spread of COVID-19 with new fractal-fractional operators: can the lockdown save mankind before vaccination, Chaos, Solitons Fractals, 136 (2020), 109860. doi: 10.1016/j.chaos.2020.109860
[12]	B. Ghanbari, A. Atangana, A new application of fractional Atangana–Baleanu derivatives: designing ABC-fractional masks in image processing, Physica A: Stat. Mech. Appl., 542 (2020), 123516. doi: 10.1016/j.physa.2019.123516
[13]	A. G. M. Selvam, J. Alzabut, D. A. Vianny, M. Jacintha, F. B. Yousef, Modeling and stability analysis of the spread of novel coronavirus disease COVID-19, Int. J. Biomath., 2021 (2021), 2150035.
[14]	F. Bozkurt, A. Yousef, D. Baleanu, J. Alzabut, A mathematical model of the evolution and spread of pathogenic coronaviruses from natural host to human host, Chaos, Solitons Fractals, 138 (2020), 109931. doi: 10.1016/j.chaos.2020.109931
[15]	K. M. Owolabi, A. Atangana, A. Akgul, Modelling and analysis of fractal-fractional partial differential equations: application to reaction-diffusion model, Alexandria Eng. J., 59 (2020), 2477–2490. doi: 10.1016/j.aej.2020.03.022
[16]	H. Mohammadi, S. Rezapour, A. Jajarmi, On the fractional SIRD mathematical model and control for the transmission of COVID-19: the first and the second waves of the disease in Iran and Japan, ISA Trans., 2021.
[17]	A. A. Kilbas, H. M. Srivastava, J. J. Trujillo, Theory and Applications of Fractional Differential Equations, Elsevier, 204 (2006).
[18]	H. M. Srivastava, K. M. Saad, J. F. Gómez-Aguilar, A. A. Almadiy, Some new mathematical models of the fractional-order system of human immune against IAV infection, Math. Biosci. Eng., 17 (2020), 4942–4969. doi: 10.3934/mbe.2020268
[19]	N. H. Tuan, H. Mohammadi, S. Rezapour, A mathematical model for COVID-19 transmission by using the Caputo fractional derivative, Chaos, Solitons Fractals, 140 (2020), 110107. doi: 10.1016/j.chaos.2020.110107
[20]	L. F. Ávalos-Ruiz, J. F. Gomez-Aguilar, A. Atangana, K. M. Owolabi, On the dynamics of fractional maps with power-law, exponential decay and Mittag–Leffler memory, Chaos, Solitons Fractals, 127 (2019), 364–388. doi: 10.1016/j.chaos.2019.07.010
[21]	K. M. Owolabi, A. Atangana, Numerical approximation of nonlinear fractional parabolic differential equations with Caputo–Fabrizio derivative in Riemann–Liouville sense, Chaos, Solitons Fractals, 99 (2017), 171–179. doi: 10.1016/j.chaos.2017.04.008
[22]	O. A. Arqub, Application of residual power series method for the solution of time-fractional Schrödinger equations in one-dimensional space, Fundam. Informaticae, 166 (2019), 87–110. doi: 10.3233/FI-2019-1795
[23]	O. A. Arqub, Numerical algorithm for the solutions of fractional order systems of Dirichlet function types with comparative analysis, Fundam. Informaticae, 166 (2019), 111–137. doi: 10.3233/FI-2019-1796
[24]	S. Rezapour, H. Mohammadi, M. E. Samei, SEIR epidemic model for COVID-19 transmission by Caputo derivative of fractional order. Advances in difference equations, Adv. Differ. Equations, 1 (2020), 1–9.
[25]	M. Caputo, Michele, M. Fabrizio, A new definition of fractional derivative without singular kernel, Progr. Fract. Differ. Appl., 1 (2015), 1–13.
[26]	H. A. A. El-Saka, The fractional-order SIS epidemic model with variable population size, J. Egypt. Math. Soc., 22 (2014), 50–54. doi: 10.1016/j.joems.2013.06.006
[27]	R. Toledo-Hernandez, V. Rico-Ramirez, G. A. Iglesias-Silva, U. M. Diwekar, A fractional calculus approach to the dynamic optimization of biological reactive systems. Part I: Fractional models for biological reactions, Chem. Eng. Sci., 117 (2014), 217–228. doi: 10.1016/j.ces.2014.06.034
[28]	K. M. Owolabi, Computational analysis of different Pseudoplatystoma species patterns the Caputo-Fabrizio derivative, Chaos, Solitons Fractals, 144 (2021), 110675. doi: 10.1016/j.chaos.2021.110675
[29]	D. Baleanu, A. Mousalou, S. Rezapour, A new method for investigating approximate solutions of some fractional integro-differential equations involving the Caputo-Fabrizio derivative, Adv. Differ. Equations, 2017 (2017), 1–12. doi: 10.1186/s13662-016-1057-2
[30]	M. S. Aydogan, D. Baleanu, A. Mousalou, S. Rezapour, On high order fractional integro-differential equations including the Caputo-Fabrizio derivative, Boundary Value Probl., 2018 (2018), 1–15. doi: 10.1186/s13661-017-0918-2
[31]	D. Baleanu, S. Rezapour, Z. Saberpour, On fractional integro-differential inclusions via the extended fractional Caputo-Fabrizio derivation, Boundary Value Probl., 2019 (2019), 1–17. doi: 10.1186/s13661-018-1115-7
[32]	D. Baleanu, A. Jajarmi, H. Mohammadi, S. Rezapour, A new study on the mathematical modelling of human liver with Caputo–Fabrizio fractional derivative, Chaos, Solitons Fractals, 134 (2020), 109705. doi: 10.1016/j.chaos.2020.109705
[33]	S. Rezapour, H. Mohammadi, A study on the AH1N1/09 influenza transmission model with the fractional Caputo-Fabrizio derivative, Adv. Differ. Equations, 2020 (2020), 1–15. doi: 10.1186/s13662-019-2438-0
[34]	S. Rezapour, H. Mohammadi, A. Jajarmi, A new mathematical model for Zika virus transmission, Adv. Differ. Equations, 2020 (2020), 1–15. doi: 10.1186/s13662-019-2438-0
[35]	D. Baleanu, H. Mohammadi, S. Rezapour, A mathematical theoretical study of a particular system of Caputo-Fabrizio fractional differential equations for the Rubella disease model, Adv. Differ. Equations, 2020 (2020), 1–19. doi: 10.1186/s13662-019-2438-0
[36]	N. M. Ferguson, D. Laydon, G. Nedjati-Gilani, N. Imai, K. Ainslie, M. Baguelin, et al., Impact of non-pharmaceutical interventions (NPIs) to reduce COVID-19 mortality and healthcare demand. Imperial College COVID-19 Response Team, Imp. Coll. COVID-19 Response Team, 2020.
[37]	A. J. Kucharski, T. W. Russell, C. Diamond, Y. Liu, J. Edmunds, S. Funk, et al., Early dynamics of transmission and control of COVID-19: a mathematical modelling study, Lancet Infect. Dis., 20 (2020), 553–558. doi: 10.1016/S1473-3099(20)30144-4
[38]	T. Kuniya, Prediction of the epidemic peak of coronavirus disease in Japan, 2020, J. Clin. Med., 9 (2020), 789. doi: 10.3390/jcm9030789
[39]	R. O. Stutt, R. Retkute, M. Bradley, C. A. Gilligan, J. Colvin, A modelling framework to assess the likely effectiveness of facemasks in combination with 'lock-down'in managing the COVID-19 pandemic, Proc. R. Soc. A, 476 (2020), 2238.
[40]	Z. Tang, X. Li, H. Li, Prediction of new coronavirus infection based on a modified SEIR model, medRxiv, 2020.
[41]	Z. Tang, X. Li, H. Li, New applications related to Covid-19, Results Phys., 20 (2021), 103663. doi: 10.1016/j.rinp.2020.103663
[42]	A. Atangana, S. İ. Araz, Nonlinear equations with global differential and integral operators: existence, uniqueness with application to epidemiology, Results Phys., 20 (2021), 103593. doi: 10.1016/j.rinp.2020.103593
[43]	M. Rahman, S. Ahmad, R. T. Matoog, N. A. Alshehri, T. Khan, Study on the mathematical modelling of COVID-19 with Caputo-Fabrizio operator, Chaos, Solitons Fractals, 2021 (2021), 111121.
[44]	H. M. Dobrovolny, Modeling the role of asymptomatics in infection spread with application to SARS-CoV-2, Plos one, 15 (2020), e0236976. doi: 10.1371/journal.pone.0236976
[45]	H. M. Dobrovolny, Quantifying the effect of remdesivir in rhesus macaques infected with SARS-CoV-2, Virology, 550 (2020), 61–69. doi: 10.1016/j.virol.2020.07.015
[46]	S. Mandal, T. Bhatnagar, N. Arinaminpathy, A. Agarwal, A. Chowdhury, M. Murhekar, et al., Prudent public health intervention strategies to control the coronavirus disease 2019 transmission in India: A mathematical model-based approach, Indian J. Med. Res., 151 (2020), 2–3.
[47]	R. F. Reis, B. de Melo Quintela, J. de Oliveira Campos, J. M. Gomes, B. M. Rocha, M. Lobosco, et al., Characterization of the COVID-19 pandemic and the impact of uncertainties, mitigation strategies, and underreporting of cases in South Korea, Italy, and Brazil, Chaos, Solitons Fractals, 136 (2020), 109888. doi: 10.1016/j.chaos.2020.109888
[48]	X. Yang, L. Chen, J. Chen, Permanence and positive periodic solution for the single-species nonautonomous delay diffusive models, Comput. Math. Appl., 32 (1996), 109–116.
[49]	H. Ramos, Z. Kalogiratou, T. Monovasilis, T. E. Simos, An optimized two-step hybrid block method for solving general second order initial-value problems, Numer. Algorithms, 72 (2016), 1089–1102. doi: 10.1007/s11075-015-0081-8
[50]	J. Vigo-Aguiar, H. Ramos, C. Clavero, A first approach in solving initial-value problems in ODEs by elliptic fitting methods, J. Comput. Appl. Math., 318 (2017), 599–603. doi: 10.1016/j.cam.2015.12.005
[51]	C. Li, F. Zeng, Numerical Methods for Fractional Calculus, Chapman and Hall/CRC, 2019.
[52]	A. Jajarmi, D. Baleanu, A new fractional analysis on the interaction of HIV with CD4+ T-cells, Chaos, Solitons Fractals, 113 (2018), 221–229. doi: 10.1016/j.chaos.2018.06.009
[53]	A. Jajarmi, D. Baleanu, On the nonlinear dynamical systems within the generalized fractional derivatives with Mittag-Leffler kernel, Nonlinear Dyn., 94 (2018), 397–414. doi: 10.1007/s11071-018-4367-y

This article has been cited by:

1.	Yubing Jiang, Hu Chen, Local convergence analysis of L1-ADI scheme for two-dimensional reaction-subdiffusion equation, 2024, 70, 1598-5865, 1953, 10.1007/s12190-024-02037-z
2.	Aniruddha Seal, Srinivasan Natesan, Suayip Toprakseven, A Dimensional-Splitting Weak Galerkin Finite Element Method for 2D Time-Fractional Diffusion Equation, 2024, 98, 0885-7474, 10.1007/s10915-023-02448-3

Reader Comments

Your name:*

Email:*
© 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)