Abstract

In this paper, we investigate a class of fractional-order time-varying delays gene regulatory networks with structured uncertainties and controllers (DFGRNs). Our contributions lie in three aspects: first, a necessary and sufficient condition on the existence of the solution for the DFGRNs is given by using the properties of the Riemann–Liouville fractional derivative and Caputo’s fractional derivative; second, the unique solution of the DFGRNs is proved under given initial function and certain condition; third, some novel sufficient conditions on finite-time stability of the DFGRNs are established by using a generalized Gronwall inequality and norm technique, and some conclusions on the finite-time stability of the DFGRNs with memory state-feedback controllers are reached, and those conditions and conclusions depend on the fractional order of the DFGRNs. One of the most interesting findings is that the “estimated time” of the finite-time stability is indeed related to the structured uncertainties, state-feedback controllers, time delays, and the fractional order of the system.

1. Introduction

Genetic regulatory networks (GRNs), which describe the interaction functions in gene expressions between DNAs, RNAs, proteins, and small molecules in an organism, are fundamental and important biological networks. The analysis and control of GRNs involve two aspects: first, understanding the widespread phenomena in living organisms and providing potential routes to prolong life span, cure cancer and diabetes, and so on; second, potential application of GRNs in the development of related disciplines, such as synthetic biology, network medicine, and personalized medicine [15].

With the development of sequencing technology, more and more genes and their regulatory sites are discovered. The structures and functions of a large number of genes are confirmed through experimental techniques, and even the regulatory mechanisms of the gene expressions controlled by some proteins are also identified [2]. Some qualitative models were proposed for distinguishing these regulatory mechanisms, such as directed graphs, Boolean networks, generalized logical networks, and rule-based formalisms [2]. However, it is difficult to know the biochemical reaction mechanisms’ underlying regulatory interactions when qualitatively handling great quantities of experimental data. Therefore, studying GRNs needs more accurately quantitative model. Based on experimental data, some simple genetic networks have been constructed, for example, genetic repressilator network [5], negative feedback GRNs [6], and genetic switch network [7]. The results in these experiments show that quantitative mathematical modeling approaches on dynamical systems have been great tools in providing insights on the mechanisms underlying the structure and the behaviors of GRNs [812].

Using integer-order differential equation to model GRNs is a classical method. However, the fractional-order differential equations are more suitable for modeling the gene regulatory mechanism. Ji et al. [13] applied the particle swarm optimization technique in modeling the fractional-order GRNs with eight real target genes. The experimental results confirmed that the fractional-order model has achieved much lower fitting error on test data than integer-order model. Other studies also revealed that the fractional-order systems have excellent performance in describing the memory and hereditary properties of various processes in GRNs, which could be far better than the integer-order ones [811].

Due to slow biochemical processes such as gene transcription, translation, and transportation (the synthesis of mRNAs and proteins at nucleus and cytoplasm, respectively, in eukaryotic cells), time delays are omnipresent in GRNs [14]. Many nonlinear differential equations with time delays have been proposed to model general GRNs, and the important role of time delays in dynamics of GRNs is now widely accepted [1518]. Actually, time delays often degrade the system performance or destabilize the system [1, 19, 20]; even GRN models without time delay may generate wrong predictions [21]. As time delays often change with time and their precise measurement is difficult in real GRNs, the dynamics of fractional-order linear and nonlinear systems with time-varying delays has attracted increasing interest, and the results show that it is naturally of better practical significance than those with constant delays [2125].

In addition, in order to avoid undesirable states associated with disease, the control of GRNs is often regarded as developing therapeutic intervention strategies for some diseases [26, 27]. And many literatures focus on the research of control in the dynamic system [2833]. In [32], the authors obtained some stabilization results for neural networks with leakage delay by designing state-feedback controller. Ebihara et al. [33] discovered that exact robust control is indeed attained for discrete-time linear systems by designing periodically time-varying memory state-feedback controller. Therefore, it is necessary to consider the controller for the DFGRNs.

Since the modeling of GRNs is underlined with the real-world gene expression time-series data, some limitations of the current experimental techniques in GRNs make the modeling errors and parameter fluctuations unavoidable. Moreover, some point out that the system parameters identified with the experimental data may construct an unknown but bounded time-varying function, and this unknown nature is referred to as the structural uncertainty or the parametric uncertainty, also known as variation or fluctuation [34]. As is known, the structural uncertainties in GRNs may lead to the poor performance or even instability in real genetic networks [28, 3437]. In [28], the authors studied the robust stabilization and state-feedback controller design for a class of integer-order GRNs with time-varying delays (DGRNs) and structured uncertainties and established some delay-dependent stability results by using some matrix techniques. Therefore, taking into account the structural uncertainties while investigating the dynamical behaviors of DFGRNs is essential.

Since the expression of gene and mRNA-translated protein is accomplished in a much relatively short period, in recent decades, some scholars have paid more attention to the finite-time stability of GRNs [4, 38]. For example, Wu et al. [4] investigated the finite-time stability associated with a class of integer-order GRNs by designing adaptive controllers. Wang et al. [38] established some new sufficient conditions of the finite-time stability for a class of integer-order uncertain GRNs with time-varying delays. Lazarević [22] investigated the finite-time stability for fractional-order nonlinear differential equation with time-varying delays by using generalized Gronwall inequality and the classical Bellman–Gronwall inequality, respectively. Phat and Thanh [23] established some new sufficient conditions of robust finite-time stability for a class of nonlinear fractional-order differential systems with time-varying delays. Wang et al. [39] considered a class of nonlinear fractional-order systems with constant delays and studied the existence and uniqueness of the solution for this kind of systems by using relevant properties of the fractional derivative.

However, the discussions on the existence and uniqueness of the solutions and the finite-time stability results for the fractional-order uncertain GRNs with time-varying delays and controllers seem rare.

From above discussions, we focus on the existence and uniqueness of the solution and the finite-time stability for a class of DFGRNs with structured uncertainties and controllers. The remainder of this paper is organized as follows. In Section 2, we give the model description, some definitions, and related properties on fractional calculus. In Section 3, we discuss the existence and uniqueness of the solution and give some sufficient criteria on the finite-time stability for the DFGRNs. In Section 4, we perform some numerical simulations, which support our findings. In Section 5, we briefly review and summarize the main results.

2. Problem Description and Preliminaries

For any vector and matrix , denote

Let be the largest singular value of matrix,where are positive constants that satisfy the later assumptions (I) and (II), respectively.

We will focus on a class of DFGRNs with structured uncertainties and controllers, which is established as follows:wherein which represents Caputo’s fractional derivative and are the concentrations of mRNA and protein of the th node, respectively. The parameters and are the decay rates of mRNA and protein, respectively; are the translation rates; are the translation rates. Both and represent the feedback regulation of the protein on the transcription. Generally, each one of the two functions is a nonlinear function but has a form of monotonicity with its variable. As a monotonic increasing or decreasing regulatory function, and are usually of the Michaelis–Menten or Hill forms [21]. , where and are bounded constants which are, respectively, the dimensionless transcriptional rates of transcription factor to at time and , and , respectively, are the set of all the where the transcription factor is a repressor of gene at time and . are the coupling matrices of the gene network, which are defined as follows:

The transcriptional delay and translational delay are bounded continuous functions on with ; here is a positive constant. , are controller vectors, are dimensions coefficient matrices. , , , , , are norm-bounded unknown matrices with time-varying structured uncertainties.

The initial conditions for DFGRN (3) are as follows:where is the given initial function with and .(i)Assumption (I): the norm-bounded unknown matrices satisfy the following inequalities:where are positive constants.(ii)Assumption (II): the functions satisfy the following inequalities:where are positive constants.

Next, we give some definitions and lemmas.

Definition 1 (see [40]). The fractional integral of order for a function is defined aswhere . The gamma function is defined by the integral .

Definition 2 (see [40]). Caputo’s fractional derivative of order for a function is defined bywhere and is a positive integer such that .

Definition 3 (see [40]). The Riemann–Liouville fractional derivative of order for a function is defined aswhere and is a positive integer such that .
For convenience, we choose the notation , , .

Definition 4. A mild solution of DFGRN (3) with initial condition (6) is a vector composed of continuous functionssatisfying

Definition 5 (see [22]). The system given by (3) (when ) satisfying the initial condition (6) is finite-time stable with respect to , if and only if imply .

Definition 6 (see [22]). The system given by (3) satisfying the initial condition (6) is finite-time stable with respect to , if and only if and imply , where is a positive constant.

Lemma 1 (see [40]). If and , then(i).(ii).(iii).

Lemma 2 (see [41]). Suppose ; if (some ), is a locally integrable nonnegative function, is a nonnegative and nondecreasing continuous function, (constant), and is a nonnegative and locally integrable function with , then

In addition, if is a nondecreasing function, then , where is the Mittag-Leffler function defined by .

3. Main Results

3.1. The Existence and Uniqueness of the Mild Solution of DFGRNs

Theorem 1. Continuously differentiable functions form a mild solution to DFGRN (3) with initial condition (6) if and only if

Proof. We firstly give the sufficient condition of the existence of the mild solution to DFGRN (3).
When , is obvious. For , according to (15), applying and property (ii) of Lemma 1, we obtainAccording to the property (iii) of Lemma 1 and , we getFrom (16) and (17), we haveWe secondly give the necessary condition of the existence of the mild solution to DFGRN (3).
When , the solution of DFGRN (3) isIf , from DFGRN (3), we haveIn the case of , from property (i) of Lemma 1, we can obtainThe proof is completed.

Theorem 2. If assumptions (I) and (II) hold, then DFGRN (3) with initial condition (6) has a unique mild solution.

Proof. Let and be any two different solutions to DFGRN (3) with initial condition (6); denote , , . According to Theorem 1, we know that both and satisfy condition (15).
If , thenFrom (22), by using the norm and assumptions (I) and (II), we can obtainFirst, when , . So, for .
Second, when , from (23), we haveIn accordance with (24), we can obtainFrom Lemma 2, we can getThus, . That is to say, . So, for .
Third, when , according to (23), we haveLet . Then, we obtainFrom Lemma 2, we can getThen, we have . So, for .
Summarizing the above three cases, we can obtain that , for . Due to the arbitrary nature of the solution and of DFGRN (3) and in accordance with Definition 4, we can conclude that DFGRN (3) has a unique mild solution. The proof is completed.

3.2. Finite-Time Stability of DFGRNs with Structured Uncertainties

Theorem 3. If assumptions (I) and (II) and hold, then the uncertain DFGRNs with controllers given by (3) with initial condition (6) are finite-time stable with respect to , where .

Proof. According to Theorem 1 and Theorem 2, we can know that DFGRN (3) has a mild solution and the solution satisfies the following integral equation:Using the norm , we can obtain the solution estimate of system (30):By applying norm to DFGRN (3) and combining assumptions (I) and (II), we can getLet . According to (31), (3), and (32), if , we haveLetThen, we know that is a nonnegative and nondecreasing function. By using Lemma 2 (the generalized Gronwall inequality), we haveIf , we haveBecause and , thenHence,The proof is completed.

Remark 1. If we adopt in DFGRN (3), we can obtain the following conclusion.
The uncertain DFGRN (3) satisfying the initial condition (6) is finite-time stable with respect to , if assumptions (I) and (II) hold and the following condition is satisfied:where .

Remark 2. In the proof of Theorem 3, if we use the “classical” Bellman–Gronwall inequality instead of the generalized Gronwall inequality, we can get the following result.
The uncertain DFGRN with controllers given by (3) satisfying the initial condition (6) is finite-time stable with respect to , if assumptions (I) and (II) hold and the following condition is satisfied:

Remark 3. If we take in system (3), the above results turn into the following conclusion.
The uncertain DFGRN (3) satisfying the initial condition (6) is finite-time stable with respect to . If assumptions (I) and (II) hold, the following condition is satisfied:

3.3. Finite-Time Stability of DFGRNs with Memory State-Feedback Controllers

We consider the following memory state-feedback controllers on DFGRN (3):where are the gain matrices of . Then, DFGRN (3) can be changed into

Theorem 4. If assumptions (I) and (II) andhold, then the uncertain DFGRN (3) with memory state-feedback controllers given by (43) satisfying the initial condition (6) is finite-time stable with respect to , .

Proof. Similar to Theorem 1 and Theorem 2, it is easy to prove that DFGRN (43) has a mild solution satisfying the following integral equation:Using the norm , we haveFrom (43), by using assumptions (I) and (II), we haveLet . From (46) and (47), we obtainHence,LetThen, we know is a nonnegative and nondecreasing function. From Lemma 2, we haveIf , we obtainFrom the condition of and , we can getTherefore,The proof is completed.

Remark 4. Similar to Remark 2, we can get the following result.
The uncertain DFGRN (3) with memory state-feedback controller given by (43) satisfying the initial condition (6) is finite-time stable with respect to , if assumptions (I) and (II) hold and the following condition is satisfied:

Remark 5. We can obtain the same conclusion as Theorem 3 and Theorem 4 if the inequalities in assumption (II) are

Remark 6. All the results in Remarks 14 are still new.

4. Numerical Examples

In this section, some numerical examples are given to illustrate the effectiveness of above theoretical results. In the following examples, the functions and are taken as the Hill form. And in the Adams–Bashforth–Moulton predictor-corrector scheme [42], the step length is .

Example 1. Consider the following DFGRNs of three mRNA and protein nodes with structured uncertainties and memory state-feedback controllers:LetThe memory state-feedback controllers are defined as follows:whereLet , , , . According to the notations in Section 2, we obtain , , . When , simple computation reveals thatFrom Theorem 4, system (57) is finite-time stable with respect to . Denote as the “estimated time” of finite-time stability. The transient states of the variable and of DFGRN (57) with and are shown in Figures 1(a) and 1(b), respectively.

Example 2. Consider the following DFGRNs of three mRNA and protein nodes with structured uncertainties and without controller:Using the same parameters in Example 1, we similarly get . When , we haveFrom Remark 1, system (62) is finite-time stable with respect to ; then, the “estimated time” of finite-time stability . The transient states of the variables and of DFGRN (62) with and are shown in Figures 2(a) and 2(b), respectively.
In Example 2, when , the case of infinite time, DFGRN (62) with structured uncertainties is unstable. The numerical simulations of the variables and of DFGRN (62) with and are shown in Figures 3(a) and 3(b), respectively.

Remark 7. It is worthy to note that in a special case of DFGRN (62) without structured uncertainties, it is proved that in the sense of infinite stability, (62) is globally asymptotically stable [16].

Example 3. Consider the following DFGRNs of three mRNA and protein nodes with memory state-feedback controllers and without structured uncertainties:Using the same parameters in Example 1, we similarly obtain . When , we can getFrom Theorem 4, system (64) is finite-time stable with respect to ; then, the “estimated time” of finite-time stability . The transient states of the variables and of DFGRN (64) with and are shown in Figures 4(a) and 4(b), respectively.

Example 4. Consider the following DFGRNs of three mRNA and protein nodes without structured uncertainties or controller:Using the same parameters in Example 1, we also obtain the “estimated time” of finite-time stability for system (66) as . The transient states of the variables and of DFGRN (66) with and are shown in Figures 5(a) and 5(b), respectively.
If we adopt constant time-delay and in DFGRN (66), then system (66) is finite-time stable, and the “estimated time” of finite-time stability is 0.0315. The transient states of the variables and of DFGRN (66) with are shown in Figure 6.

Remark 8. If and in DFGRN (66), then system (66) converts to system (4.1) in [16]. When , it is proved that system (4.1) is unstable in the sense of infinite-time stability [16], which means that the finite-time stability is different from the infinite-time stability of DFGRNs.
If we take and in DFGRNs (57), (62), (64), and (66) systems (57), (62), (64), and (66) convert to the corresponding fractional-order gene regulatory networks without time delays (FGRNs).
In order to investigate the effects of structured uncertainties, controllers and time delays on the stability of the DFGRNs, we calculate the “estimated time” of finite-time stability for above four examples and the corresponding FGRNs with different fractional-order ; the results are shown in Tables 1 and 2, respectively.
From Table 1 or Table 2, we have the following conclusions:(i)The effect of the controllers: comparing column 2 with 3 (or column 4 with 5), we can know that the controllers can shorten the “estimated time” of finite-time stability under the same conditions of fractional-order q and structured uncertainties.(ii)The effect of the structured uncertainties: comparing column 3 with 5, we can know that the structured uncertainties can shorten the “estimated time” of finite-time stability under the same fractional-order it q.(iii)The difference between the structured uncertainties and the controllers: comparing column 3 with 4, we can know that the size of “estimated time” of finite-time stability for DFGRN (62) with structured uncertainties is longer than DFGRN (64) with controllers under the same fractional-order q.(iv)The effect of the fractional-order : in the same column, we can know that decreasing the fractional-order will be useful to decrease the “estimated time” of finite-time stability for DFGRNs or FGRNs.(v)The effect of time delays: comparing Table 1 with Table 2, we can know that the “estimated time” of finite-time stability is reduced under the same fractional-order when considering time delays.

5. Concluding Remarks

This paper deals with the existence and uniqueness of the solution and the finite-time stability for a class of DFGRNs with structured uncertainties and controllers. In particular, we design the memory state-feedback controllers for DFGRNs with structured uncertainties and give the sufficient conditions for the system to achieve the finite-time stability.

It should be pointed out that the conditions of finite-time stability in the present paper are dependent on the fractional-order , which is more different from the previous stability results for the case of integer order, i.e., the finite-time stability is independent of the integer order.

In addition, from the numerical results, we find that all of the controllers, uncertain terms, fractional-order and time delays can affect the “estimated time” of finite-time stability. Particularly, (i) the size of “estimated time” of finite-time stability with controllers is shorter than the case without controller but only with structured uncertainties, which means that the controllers are more beneficial for controlling the “estimated time” than the structured uncertainties; (ii) the size of “estimated time” of finite-time stability with time delays is shorter than the case without time delays, which means that time delays degrade the GRN performance.

If we take and controllers terms , meanwhile, in the special case constant time delay, system (3) convert to (2.2) in [16], and we find that numerically: as , DFGRN (62) in this paper is unstable; however, DFGRN (4.1) in [16] is globally asymptotically stable, which means that the structured uncertainty can change the stability of DFGRNs. Furthermore, from Remark 8, we know that DFGRN (66) is finite-time stable, while the corresponding system (4.1) in [16] is infinite-time unstable, which means that an infinite-time unstable system can change to a finite-time stable one under extra conditions. The analytical study on above questions is desirable in the future.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This study was supported by the Hunan Provincial Natural Science Foundation (nos. 2019JJ50222 and 13JJ4065) and the Scientific Research Fund of Hunan Provincial Education Department (no. 19C0911).