Next Article in Journal
Construction of a Family of Maximally Entangled Bases in d d′
Previous Article in Journal
A TOPSIS-Inspired Ranking Method Using Constrained Crowd Opinions for Urban Planning
Previous Article in Special Issue
A Fractional-Order Sinusoidal Discrete Map
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stability of Gene Regulatory Networks Modeled by Generalized Proportional Caputo Fractional Differential Equations

1
Center for Research and Development in Mathematics and Applications, Department of Mathematics, University of Aveiro, 3810-193 Aveiro, Portugal
2
Department of Mathematics, Texas A&M University-Kingsville, Kingsville, TX 78363, USA
3
Faculty of Mathematics and Informatics, Plovdiv University “P. Hilendarski”, 4000 Plovdiv, Bulgaria
4
School of Mathematical and Statistical Sciences, National University of Ireland, H91 TK33 Galway, Ireland
*
Author to whom correspondence should be addressed.
Entropy 2022, 24(3), 372; https://doi.org/10.3390/e24030372
Submission received: 7 February 2022 / Revised: 28 February 2022 / Accepted: 4 March 2022 / Published: 5 March 2022
(This article belongs to the Special Issue Fractional Order Systems: From Local Behavior to Network Dynamics)

Abstract

:
A model of gene regulatory networks with generalized proportional Caputo fractional derivatives is set up, and stability properties are studied. Initially, some properties of absolute value Lyapunov functions and quadratic Lyapunov functions are discussed, and also, their application to fractional order systems and the advantage of quadratic functions are pointed out. The equilibrium of the generalized proportional Caputo fractional model and its generalized exponential stability are defined, and sufficient conditions for the generalized exponential stability and asymptotic stability of the equilibrium are obtained. As a special case, the stability of the equilibrium of the Caputo fractional model is discussed. Several examples are provided to illustrate our theoretical results and the influence of the type of fractional derivative on the stability behavior of the equilibrium.

1. Introduction

Gene expression is the process where the hereditary code of a gene is used for synthesizing proteins and producing the structures of the cell. Genes that code for amino acid sequences are named ‘structural genes’. Gene expression processes include two main stages known as ‘Transcription and translations’. Transcription is the creating of messenger RNA (mRNA) by the enzyme RNA polymerase and the processing of the resulting mRNA molecule. A gene regulatory network consists of a number of genes interacting by proteins. Mathematical models of gene regulatory networks are described and studied in several papers (see, for example, [1,2], for fractional order [3,4,5,6], and with delays [7,8]).
Recently, fractional calculus, fractional derivatives, and fractional integrals of various types have been extensively studied and applied in mathematical modeling. The memory property of fractional derivatives makes them well suited in modeling and describing the complex nature of real-world problems, in comparison to local derivatives (see, for example [9,10,11]).
In this paper, a gene regulated model with the generalized proportional Caputo fractional derivative is set up, and the equilibrium is defined. The generalized exponential stability is introduced and studied via the application of Lyapunov functions and their generalized Caputo proportional fractional derivatives. Generalized proportional Caputo fractional derivatives were recently introduced (see [12,13]); this type of derivative is a generalization of the Caputo fractional derivative, and their application provides wider possibilities for modeling adequately the complexity of real-world problems. The stability of fractional order systems with a proportional Caputo fractional derivatives is quite recent (see, for example, [14,15]). In this paper, some properties of absolute values of Lyapunov functions and their fractional derivatives are discussed, and several examples are provided to illustrate the properties. The advantages of the application of the quadratic Lyapunov functions are considered, and sufficient conditions for generalized exponential stability and asymptotic stability are obtained. Several examples are provided to illustrate the theoretical results and the dependence of the fractional derivative on the behavior of the solutions.

2. Notes on Fractional Calculus

We recall the definitions needed in this paper, namely fractional integrals and derivatives (cf. [13]):
The generalized proportional fractional integral of a function u : [ a , ) R is defined by (as long as all integrals are well defined)
( a I q , ρ u ) ( t ) = 1 ρ α Γ ( α ) a t e ρ 1 ρ t s t s q 1 u ( s ) d s , t ( a , b ] , q 0 , ρ ( 0 , 1 ] .
The generalized Caputo proportional fractional derivative of a function u : [ a , ) R is defined by (as long as all integrals are well defined)
( a C D q , ρ u ) ( t ) = ( a I 1 q , ρ ( D 1 , ρ u ) ) ( t ) = 1 ρ 1 q Γ ( 1 q ) a t e ρ 1 ρ t s t s q ( D ρ u ) ( s ) d s , t ( a , b ] , q ( 0 , 1 ) , ρ ( 0 , 1 ] ,
where ( D ρ u ) ( t ) = ( 1 ρ ) u ( t ) + ρ u ( t ) .
Remark 1.
Note that the generalized proportional Caputo fractional derivative is defined for u C ( [ a , b ] , R n ) via component-wise.
Remark 2.
If ρ = 1 , then the generalized Caputo proportional fractional derivative reduces to the classical Caputo fractional derivative of order q ( 0 , 1 ) :
( a C D q , ρ u ) ( t ) = a C D q u ( t ) = 1 Γ ( 1 q ) a t t s q u ( s ) d s .
Definition 1.
We say u C q , ρ ( [ t 0 , T ] , R n ) if u ( · ) is differentiable and the generalized proportional Caputo fractional derivative ( a C D q , ρ u ) ( t ) exists for all t ( t 0 , T ] .
Lemma 1.
Let q , ρ ( 0 , 1 ) . Then, the generalized proportional fractional derivative of a constant K R is
( a C D q , ρ K ) ( t ) = ( 1 ρ ) K ρ 1 q Γ ( 1 q ) a t e ρ 1 ρ t s t s q d s = K ρ 1 q Γ ( 1 q ) 1 ( ρ 1 ) q Γ 1 + q , ρ 1 ρ ( t a ) Γ ( 1 + q ) , f o r t > a .
Proposition 1.
([13], Remark 3.2) The relation
( a C D q , ρ e ρ 1 ρ ( . ) ) ( t ) = 0 , f o r t > a
holds.
We will use the following property of the Mittag–Leffler function with one parameter, defined by E q ( z ) = k = 0 z k Γ ( 1 + k q ) with Γ ( ξ ) the gamma function.
Proposition 2.
([16], Theorem 1.2) For every q ( 0 , 1 ) , the function e t / q E q ( t q ) is completely monotonic.
Corollary 1.
[16] If q ( 0 , 1 ) , then E q ( t q ) < e t / q , t 0 .

3. Some Comments on Properties of the Fractional Derivatives of Lyapunov Functions

One of the most applicable Lyapunov functions is the absolute values Lyapunov function. In connection with this, we will give and discuss some results about their fractional derivatives.
In [17], the following result is proved:
Proposition 3.
([17], Lemma 12) If u C 1 ( [ 0 , + ) , R ) is a continuously differentiable function, and the following relation holds almost everywhere:
t 0 C D q | u ( t ) | = ( s i g n u ( t ) ) t 0 C D q u ( t ) , 0 < q < 1 .
This result is applied by many authors to study the stability of various types of Caputo fractional differential equations and models. For example, this equality is applied in the proof of global Mittag–Leffler stability for fractional-order gene regulatory networks in [3], and to study global uniform asymptotical stability for fractional-order gene regulatory networks with delays in [4,18]. Unfortunately, this equality is not satisfied for all continuously differentiable functions, and we will demonstrate this with an example.
Example 1.
Let u ( t ) = 0.5 t , t 0 .
For any t ( 0 , 0.5 ) we have u ( t ) > 0 , u ( t ) = 1 , s i g n ( 0.5 t ) = 1 , | u ( t ) | = 0.5 t , | u ( t ) | = 1 , and
0 C D q | u ( t ) | = 1 Γ ( 1 0.5 ) 0 t t s 0.5 d s = 2 t 0.5 Γ ( 0.5 ) = 1 Γ ( 0.5 ) 0 t t s 0.5 ( 0.5 s ) d s = 0 C D q u ( t ) = ( s i g n u ( t ) ) 0 C D q u ( t ) .
Let t > 0.5 . Then u ( t ) < 0 , s i g n ( 0.5 t ) = 1 , u ( t ) = 1 , | u ( t ) | = 0.5 + t , | u ( t ) | = 1 , and
0 C D q | u ( t ) | = 1 Γ ( 0.5 ) 0 t t s 0.5 ( | u ( s ) | ) d s = 1 Γ ( 0.5 ) 0 0.5 t s 0.5 ( | u ( s ) | ) d s + 1 Γ ( 0.5 ) 0.5 t t s 0.5 ( | u ( s ) | ) d s = 1 Γ ( 0.5 ) 0 0.5 t s 0.5 d s + 1 Γ ( 0.5 ) 0.5 t t s 0.5 d s = 1 Γ ( 0.5 ) ( 4 ( t 0.5 ) 0.5 2 t 0.5 ) 2 t 0.5 Γ ( 0.5 ) = 1 Γ ( 0.5 ) 0 t t s 0.5 d s = 1 Γ ( 0.5 ) 0 t t s 0.5 ( 0.5 s ) d s = 1 Γ ( 0.5 ) 0 t t s 0.5 ( u ( s ) ) d s = ( s i g n u ( t ) ) 0 C D q u ( t ) .
Therefore, for the Caputo fractional derivative, the equality
0 C D q | u ( t ) | = ( s i g n u ( t ) ) 0 C D q u ( t )
is not true for all t [ t 0 , T ] and any function u C 1 ( [ 0 , ) , R ) .
Note that in the proof of ([5], Theorem 1), the inequality of the type u ( 0 ) u ( t ) 0 is applied to prove the equality 0 C D q | u ( t ) | = ( s i g n u ( t ) ) 0 C D q u ( t ) . Unfortunately, this inequality is not true for all functions. We will illustrate this with an example:
Example 2.
Consider the Caputo fractional differential equation
0 C D 0.5 u ( t ) = 2 u ( t ) f ( u ( t ) ) , t > 0
with u ( 0 ) = c , c R and f ( u ) = sin ( u ) + 1.5 , which is Lipschitz with a constant L = 1 and 0.5 f ( u ) = sin ( u ) 1.5 2.5 . This initial value problem is equivalent to the integral equation
u ( t ) = c + 1 Γ ( 0.5 ) f ( u ( s ) ) ( t s ) 0.5 d s .
Use
1 Γ ( 0.5 ) 0 t t s 0.5 d s = 2 t 0.5
and obtain the bounds for the solution:
c t 0.5 x ( t ) c 5 t 0.5 .
For c = 1 , we have u ( t ) < 0 for t 2 (see Figure 1), i.e., u ( 0 ) u ( t ) is not non-negative for all t 0 .
Note a similar situation occurs when the generalized proportional Caputo fractional derivative is applied. We will illustrate this with an example.
Example 3.
Let u ( t ) = 0.5 t , t 0 , q ( 0 , 1 ) , and ρ ( 0 , 1 ) .
Case 1.1. Let t ( 0 , 0.5 ) . Then, we get
( D ρ | u ( · ) | ) ( t ) = 0.5 1.5 ρ ( 1 ρ ) t = ( s i g n u ( t ) ) ( D ρ u ( · ) ) ( t ) = ( D ρ u ( · ) ) ( t )
and
( 0 C D q , ρ | u ( · ) | ) ( t ) = 1 ρ 1 q Γ ( 1 q ) 0 t e ρ 1 ρ t s t s q ( D ρ u ( · ) ) ( s ) d s = ( 0 C D q , ρ u ( · ) ) ( t ) = ( s i g n u ( t ) ) ( 0 C D q , ρ u ( · ) ) ( t ) , t ( 0 , 0.5 ) .
Case 1.2. Let t > 0.5 . Then, we get
( D ρ | u ( · ) | ) ( t ) = ( 1 ρ ) ( t 0.5 ) + ρ = ( D ρ u ( · ) ) ( t ) = ( s i g n u ( t ) ) ( D ρ u ( · ) ) ( t )
and
( 0 C D q , ρ | u ( · ) | ) ( t ) = 1 ρ 1 q Γ ( 1 q ) 0 t e ρ 1 ρ t s t s q ( s i g n u ( s ) ) ( D ρ u ( · ) ) ( s ) d s = 1 ρ 1 q Γ ( 1 q ) 0 0.5 e ρ 1 ρ t s t s q ( D ρ u ( · ) ) ( s ) d s 1 ρ 1 q Γ ( 1 q ) 0.5 t e ρ 1 ρ t s t s q ( D ρ u ( · ) ) ( s ) d s 1 ρ 1 q Γ ( 1 q ) 0 t e ρ 1 ρ t s t s q ( D ρ u ( · ) ) ( s ) d s = ( s i g n u ( t ) ) ( 0 C D q , ρ u ( · ) ) ( t ) .
Therefore, for the generalized proportional Caputo fractional derivative, the equality
0 C D q , ρ | u ( t ) | = ( s i g n u ( t ) ) 0 C D q , ρ u ( t )
is not true for all t [ t 0 , T ] and any function u C q , ρ ( [ t 0 , T ] , R ) .
We will now prove the correct result. To be general, we will consider the generalized proportional Caputo fractional derivative:
Lemma 2.
Let q ( 0 , 1 ) , ρ ( 0 , 1 ] , u C q , ρ ( [ t 0 , T ] , R ) , and suppose that the sign of u ( · ) is not changeable in [ t 0 , T ] . Then, for any t [ t 0 , T ] , the equality
( t 0 C D q , ρ | u ( · ) | ) ( t ) = ( s i g n u ( t ) ) ( t 0 C D q , ρ u ) ( t )
holds.
Proof. 
For any t [ t 0 , T ] we get
( D ρ | u ( · ) | ) ( t ) = ( 1 ρ ) ( s i g n u ( t ) ) u ( t ) + ρ ( s i g n u ( t ) ) u ( t ) = ( s i g n u ( t ) ) ( D ρ u ( · ) ) ( t )
and
( t 0 C D q , ρ | u ( · ) | ) ( t ) = 1 ρ 1 q Γ ( 1 q ) t 0 t e ρ 1 ρ t s t s q ( s i g n u ( s ) ) ( D ρ u ( · ) ) ( s ) d s = ( s i g n u ( t ) ) 1 ρ 1 q Γ ( 1 q ) t 0 t e ρ 1 ρ t s t s q ( D ρ u ( · ) ) ( s ) d s = ( s i g n u ( t ) ) ( 0 C D q , ρ u ( · ) ) ( t ) .
In the case of the Caputo fractional derivative, we obtain the following result:
Corollary 2.
Let u C q ( [ t 0 , T ] , R ) and suppose that the sign of u ( · ) is not changeable in [ t 0 , T ] . Then, for any t [ t 0 , T ] , the equality
t 0 C D q | u ( t ) | = ( s i g n u ( t ) ) ( t 0 C D q u ( t ) )
holds.
The proof of Corollary 2 follows from Lemma 2 with ρ = 1 .
Remark 3.
If the function u ( · ) changes its sign in the interval [ t 0 , T ] , then because of the memory property of the fractional derivatives (different to integer order derivatives), the equalities (7) and (8) are not true for all points t [ t 0 , T ] (see Example 1).
When Lyapunov functions are applied to differential equations, including the absolute values Lyapunov function, the type of derivatives of Lyapunov functions is very important and it depends on the type of derivatives in the differential equations.
Let u ( · ) be a solution of the scalar fractional differential equation 0 C D q u ( t ) = f ( t , u ( t ) ) . In the literature, several types of fractional derivatives of Lyapunov function V ( t , u ) are applied. Consider the following derivative, which is called by some authors the upper right-hand derivative in Caputo’s sense or Caputo-type fractional derivative:
C D q V ( u ( t ) ) = lim h 0 + sup V ( u ( t ) ) V ( u ( t ) h q f ( t , u ( t ) ) ) h q .
We will illustrate some of properties of the absolute values Lyapunov function and its derivative given by the above definition with an example.
Example 4.
Let V ( x ) = | x | . Then, its Caputo-type fractional derivative is
C D q | u ( t ) | = lim h 0 + sup | u ( t ) | | u ( t ) h q f ( t , u ( t ) ) | h q = lim ε 0 + sup | u ( t ) | | u ( t ) ε f ( t , u ( t ) ) | ε = ( s i g n ( u ( t ) ) f ( t , u ( t ) ) = ( s i g n ( u ( t ) ) 0 C D q u ( t ) ,
where ε = h q .
However, the derivative C D q V ( u ( t ) ) has no memory and C D q V ( u ( t ) ) 0 C D q V ( u ( t ) ) , so we could not conclude that C D q | u ( t ) | = ( s i g n ( u ( t ) ) C D q u ( t ) .
Remark 4.
From Example 1, Example 2, and Example 4, it could be seen that in the case of Caputo fractional differential equations, for the Lyapunov function
V ( t , x ) = k = 1 n | x i | , x = ( x 1 , x 2 , , x n ) ,
we have:
-
If its Caputo fractional derivative is applied, then the inequality
0 C D q V ( t , x ( t ) ) i = 1 n ( s i g n x i ( t ) ) 0 C D q x i ( t )
holds in the general case (see Example 1–Case 2). According to Lemma 2, the equality is true only in a particular case;
-
If its Caputo-type fractional derivative C D q | u ( t ) | is applied, then in the general case
C D q V ( u ( t ) ) i = 1 n ( s i g n x i ( t ) ) C D q x i ( t ) .
Before the application of the equality (8), one needs to prove that the solution has a constant sign (see Corollary 2).
Remark 5.
The situation mentioned in Remark 4 is true also for generalized proportional Caputo fractional differential equations and the absolute values Lyapunov function.
According to the above discussions, in this section, the Lyapunov function of the type
V ( t , x ) = k = 1 n x i 2 , x = ( x 1 , x 2 , , x n )
is appropriate to apply to fractional differential equations using the following result:
Lemma 3.
[14] Let q ( 0 , 1 ) , ρ ( 0 , 1 ] , and u C q , ρ ( [ t 0 , ) , R n ) . Then, for any t t 0 , the inequality
( t 0 C D q , ρ u T ( · ) u ( · ) ( t ) 2 x T ( t ) ( t 0 C D q , ρ u ) ( t )
holds.
Remark 6.
Note that in the special case ρ = 1 of Lemma 3, i.e., the application of Caputo fractional derivative, the result is proved in [19,20].
Remark 7.
Note that quadratic Lyapunov functions for Caputo fractional order time-delayed gene regulatory networks are applied in [7].

4. Statement of the Problem

In this paper, we will consider a class of fractional order gene regulatory networks modeled by a generalized proportional Caputo fractional derivative for 0 < q < 1 , ρ ( 0 , 1 ] :
( 0 C D q , ρ x i ) ( t ) = d i x i ( t ) + k = 1 N a i k h k ( y k ( t ) ) + I i , t > 0 , i = 1 , 2 , , N , ( 0 C D q , ρ y i ) ( t ) = b i y i ( t ) + c i x i ( t ) , t > 0 , i = 1 , 2 , , N , x i ( 0 ) = x i 0 , y i ( 0 ) = y i 0 , i = 1 , 2 , , N ,
where x j 0 , y j 0 R , x j ( t ) , y j ( t ) , j = 1 , 2 , , N , denote the concentrations of messenger ribonucleic acid (mRNA) and protein of the j-th node at time t, respectively, d j and b j are degradation velocities of mRNA and protein, respectively, c j > 0 is the translation rate, the functions h k C ( R , R ) , k = 1 , 2 , , N , represent the activator initiates of protein of mRNA, and the coupling matrix of the network A = ( a j k ) R N × N is described by
a j k = γ j k k is a repressor of gene j 0 k does not regulate gene j γ j k k is a initiator of gene j ,
and I j = k J a j k , where J is the set of all repressors of gene j.
Remark 8.
Commonly, the activator functions h k ( · ) , k = 1 , 2 , , N , are indicated in the Hill form h k ( s ) = s β k α k β k + s β k , s R , where β k are the Hill coefficients and α k 0 are constants.
Remark 9.
Note that the model (9) is studied in the case of the Caputo fractional derivative and the absolute value Lyapunov function is applied (see Remark 4).
Introduce the following assumptions:
(A1) The activator functions h k ( · ) , k = 1 , 2 , , N are increasing and there exist constants γ k > 0 such that for any u , v R with u v , the inequalities
0 h k ( u ) h k ( v ) u v γ k , k = 1 , 2 , , N
hold.
(A2) There exist positive constants μ k , μ N + k , k = 1 , 2 , , N such that the coefficients in (9) satisfy the inequalities
j = 1 N γ j | a k j | + μ N + k μ k c k < 2 d k , k = 1 , 2 , , N , 2 b k c k + γ k j = 1 N μ j μ N + k | a j k | < 2 b k , k = 1 , 2 , , N .
Remark 10.
From Assumption (A1), it follows that Lemma 2 is applicable to the solutions of (9) and equality (7) holds; i.e., the absolute value Lyapunov function is applicable to (9).
From Lemma 1, it follows that the generalized proportional Caputo fractional derivative of a nonzero constant is not zero, and applying Corollary 1, we introduce the following definition.
Definition 2.
The couple of functions
( x * ( t ) , y * ( t ) ) = C e ρ 1 ρ t , Q e ρ 1 ρ t ,
with C = ( C 1 , C 2 , , C N ) R N , Q = ( Q 1 , Q 2 , , Q N ) R N , C i , Q i = c o n s t , is called an equilibrium of (9) if
d j C j e ρ 1 ρ t = k = 1 N a j k h k Q k e ρ 1 ρ ( t ) + I j , t > t 0 , i = 1 , 2 , , N , b j Q j = c j C j , t > t 0 , i = 1 , 2 , , N .
Remark 11.
Note that in the case of Caputo fractional derivative ( ρ = 1 ), the defined equilibrium in Definition 2 coincides with the one known in the literature (see, for example, [3]).
Definition 3.
The equilibrium X * ( t ) = ( x * ( t ) , y * ( t ) ) of the model (9) is generalized exponentially stable if there exist constants M , λ > 0 such that
X ( t ) X * ( t ) M X 0 X * ( 0 ) e ρ 1 2 ρ t E q λ t ρ q , t 0 ,
where X 0 = ( x 0 , y 0 ) is the solution of (9) with initial values X ( t ) = ( x ( t ) , y ( t ) ) .
Use the transformations
u j ( t ) = x j ( t ) C j e ρ 1 ρ t , v j ( t ) = y j ( t ) Q j e ρ 1 ρ t .
Then, (9) can be written in the form
( 0 C D q , ρ u i ) ( t ) = d i u i ( t ) + k = 1 N a i k H k ( t , v k ( t ) ) , t > t 0 , i = 1 , 2 , , N , ( 0 C D q , ρ v j ) ( t ) = b j v j ( t ) + c j u j ( t ) , t > t 0 , j = 1 , 2 , , N , u j ( 0 ) = x j 0 C j , v j ( 0 ) = y j 0 Q j , j = 1 , 2 , , N ,
where
H j ( t , v k ) = h k v k + Q k e ρ 1 ρ t h k Q k e ρ 1 ρ t .
The system (11) has a zero equilibrium.
The goal of our paper is to study the exponential and asymptotic stability of the equilibrium of (9); equivalently, we also study the stability properties of the zero solution of the IVP for FrDE (11).
We will apply quadratic Lyapunov functions, and in connection with this, we will use the following result:
Lemma 4.
([14], Lemma 2) Let the function U ( t ) = ( u ( t ) , v ( t ) ) C q , ρ ( [ 0 , ) , R 2 N ) , with u ( · ) C q , ρ ( [ 0 , ) , R N ) , v ( · ) C q , ρ ( [ 0 , ) , R N ) be a solution of (11), and suppose that, for any t 0 , the inequality
0 C D q , ρ U ( t ) 2 0
holds. Then,
U ( t ) U 0 e ρ 1 2 ρ t , f o r a l l t 0 ,
where
U 0 = ( x 1 0 C 1 , x 2 0 C 2 , , x N 0 C N , y 1 0 Q 1 , , y N 0 Q N ) .
Lemma 5.
([14], Lemma 3) Let the function U ( t ) = ( u ( t ) , v ( t ) ) C q , ρ ( [ 0 , ) , R 2 N ) , with u ( · ) C q , ρ ( [ 0 , ) , R N ) , v ( · ) C q , ρ ( [ 0 , ) , R N ) , be a solution of (11) and suppose that, for any t 0 , the inequality
0 C D q , ρ U ( t ) 2 K U ( t ) 2
holds, where K > 0 is a constant. Then,
U ( t ) U 0 e ρ 1 2 ρ t E q K t ρ q , f o r a l l t 0 ,
where
U 0 = ( x 1 0 C 1 , x 2 0 C 2 , , x N 0 C N , y 1 0 Q 1 , , y N 0 Q N ) .
Theorem 1.
Let the assumptions (A1) and (A2) be satisfied, and assume that there exists an equilibrium X * ( t ) = ( x * ( t ) , y * ( t ) ) of the model (9). Then, the equilibrium of the model (9) is generalized exponentially stable.
Proof. 
The generalized exponentially stability of equilibrium of the model (9) is equivalent to the generalized exponential stability of the zero solution of (11).
Consider the Lyapunov function
V ( U ) = M U 2 = k = 1 N μ k u k 2 + k = 1 N μ N + k v k 2 ,
where
M = ( μ 1 , μ 2 , , μ N , μ N + 1 , , μ 2 N ) .
Let U ( · ) be a solution of (11). According to Lemma 3, we obtain
0 C D q , ρ V ( U ( t ) ) = k = 1 N μ k 0 C D q , ρ u k 2 ( t ) + k = 1 N μ N + k 0 C D q , ρ v k 2 ( t ) 2 k = 1 N u k ( t ) 0 C D q , ρ u k ( t ) + 2 k = 1 N v k ( t ) 0 C D q , ρ v k ( t ) = 2 k = 1 N μ k d k u k 2 ( t ) + 2 k = 1 N j = 1 N μ k a k j u k ( t ) H j ( t , v j ( t ) ) 2 k = 1 N μ N + k b k v k 2 ( t ) + 2 k = 1 N μ N + k c k v k ( t ) u k ( t ) .
From assumption (A2), it follows that
H k ( t , v k ( t ) ) v k ( t ) = h k v k ( t ) + Q k e ρ 1 ρ t h k Q k e ρ 1 ρ t v k ( t ) γ k , k = 1 , 2 , , N ,
and thus
0 C D q , ρ V ( U ( t ) ) 2 k = 1 N μ k d k u k 2 ( t ) + 2 k = 1 N j = 1 N μ k | a k j | | u k ( t ) | γ j | v j ( t ) | 2 k = 1 N μ N + k b k v k 2 ( t ) + 2 k = 1 N μ N + k c k | v k ( t ) | | u k ( t ) | 2 k = 1 N μ k d k u k 2 ( t ) + k = 1 N μ k j = 1 N | a k j | γ j u k 2 ( t ) + k = 1 N j = 1 N μ k | a k j | γ j v j 2 ( t ) 2 k = 1 N μ N + k b k v k 2 ( t ) + k = 1 N μ N + k c k u k 2 ( t ) + k = 1 N μ N + k c k v k 2 ( t ) k = 1 N 2 d k + j = 1 N γ j | a k j | + μ N + k μ k c k μ k u k 2 ( t ) k = 1 N 2 b k + c k + γ k j = 1 N μ j μ N + k | a j k | μ N + k v k 2 ( t ) λ V ( U ( t ) ) ,
where
λ = max k = 1 , 2 , , N 2 d k j = 1 N γ j | a k j | μ N + k μ k c k , 2 b k c k γ k j = 1 N μ j μ N + k | a j k | .
According to Lemma 5, the inequality
μ m i n U ( t ) = μ m i n k = 1 N u k 2 ( t ) + k = 1 N v k 2 ( t ) k = 1 N μ k u k 2 ( t ) + k = 1 N μ N + k u k 2 ( t ) = M U ( t ) M U 0 e ρ 1 2 ρ t E q λ t ρ q μ m a x U 0 e ρ 1 2 ρ t E q λ t ρ q
holds, where
μ m i n = min min k = 1 , 2 , , N μ k , min k = 1 , 2 , , N μ N + k
and
μ m a x = max max k = 1 , 2 , , N μ k , max k = 1 , 2 , , N μ N + k ,
or
U ( t ) μ m a x μ m i n U 0 e ρ 1 2 ρ t E q λ t ρ q .
Corollary 3.
Let the conditions of Theorem 1 be satisfied. Then, the equilibrium of the model (9) is asymptotically stable, i.e.,
lim t x ( t ) x * ( t ) = 0 a n d lim t y ( t ) y * ( t ) = 0 .

5. Applications

Application 1
We will consider the model of three repressor-protein concentrations, p i , and their corresponding mRNA concentrations, m i , i = 1 , 2 , 3 , which are defined and studied in [21] when the kinetics of the system are determined by ordinary differential equations. To have a more appropriate model, we will adopt this model and use generalized proportional Caputo fractional derivatives; i.e., we will consider the model
( 0 C D q , ρ m i ) ( t ) = m i + α ( 1 + p i n ) + α 0 ( 0 C D q , ρ p i ) ( t ) = β ( p i m i ) , t > 0 , i = 1 , 2 , 3 ,
where (see [21]):
-
The number of protein copies per cell produced from a given promoter type during continuous growth is α 0 in the presence of saturating amounts of repressor and α + α 0 in its absence;
-
β is the ratio of the protein decay rate to the mRNA decay rate;
-
n is a Hill coefficient.
System (19) is similar to (9) with d i = 1 , b i = c i = β , h k ( u ) = 1 1 + u n , a i i = α and a i k = 0 for k i , I i = α 0 .
Take n = 2 , α = 1.3 , and β = 1 . Thus, γ k = 0.649519 , m i = m 3 + i = 1 , and
2 b k + c k + γ k j = 1 N μ j μ N + k | a j k | = β + ( 0.649519 ) · ( 1.3 ) μ i μ 3 + i = 1 + 0.844375 = 0.155625 < 0
and
2 d k + j = 1 N γ j | a k j | + μ N + k μ k c k = 2 + ( 0.649519 ) · ( 1.3 ) + 1 = 0.155625 < 0 ,
i.e., assumptions (A1) and (A2) are satisfied. According to Theorem 1, if there exists an equilibrium E * = ( m 1 ˜ , m 2 ˜ , m 3 ˜ , p 1 ˜ , p 2 ˜ , p 3 ˜ ) ) of (19), then it is generalized exponential stable.
Case 1. Caputo fractional derivative, i.e., ρ = 1 . The equilibrium m i ˜ , p i ˜ , i = 1 , 2 , 3 is a solution of the system
m i ˜ = p i ˜ , m i ˜ = α ( 1 + p i ˜ n ) + α 0 .
The system (20) has a solution for every value of α and α 0 .
Consider a particular case of α = 1.3 and α 0 = 0 . Then, the equilibrium is m i ˜ = p i ˜ = 0.795876 , i = 1 , 2 , 3 .
The solutions m i ( · ) , i = 1 , 2 , 3 are given in Figure 2 (left), and the solutions p i ( · ) , i = 1 , 2 , 3 , are given in Figure 2 (right). It could be seen that all components of the solution approach the equilibrium 0.795876 .
Note that problem (19) is considered in [3] with α = 2.5 , α 0 = 0 . However, in this case, the equilibrium is m i ˜ = p i ˜ = 1.11475 , which does not correspond to the provided graphs.
Let α = 1.3 , α 0 = 2 . Then, the equilibrium is m i ˜ = p i ˜ = 2.21938 , i = 1 , 2 , 3 .
The solutions m i ( · ) , i = 1 , 2 , 3 are given in Figure 3 (left) and the solutions p i ( · ) , i = 1 , 2 , 3 , are given in Figure 3 (right). It could be seen that all components of the solution approach 2.21938 .
Case 2. Generalized proportional Caputo fractional derivative, i.e., ρ ( 0 , 1 ) .
Since 0 C D q , ρ 0 0 and 0 C D q , ρ e ρ 1 ρ t = 0 , t 0 , the equilibrium m i ˜ = p i ˜ = C j e ρ 1 ρ t , i = 1 , 2 , 3 , is a solution of the system
C j e ρ 1 ρ t = α 1 + C j 2 e 2 ρ 1 ρ t + α 0 , t > 0 .
Case 2.1. Let α 0 = 1.3 . Then, the system (21) has zero solution w.r.t. C j and the system (19) have a zero equilibrium. The solutions m i ( · ) , i = 1 , 2 , 3 are given in Figure 4 (left) and the solutions p i ( · ) , i = 1 , 2 , 3 are given in Figure 4 (right). It could be seen that all components approach the equilibrium 0.
Case 2.2. Let α 0 1.3 . Then, the system (21) has no solution w.r.t. C j , j = 1 , 2 , 3 , and the system (19) has no equilibrium, and we could not apply Theorem 1.
Application 2
Consider the model of three repressor-protein concentrations, p i , and their corresponding mRNA concentrations, m i , i = 1 , 2 , 3 , (19) with the activator functions h k ( s ) = s β k α k β k + s β k , s R ; i.e., consider
( 0 C D q , ρ m i ) ( t ) = m i + α p i n ( 1 + p i n ) + α 0 , ( 0 C D q , ρ p i ) ( t ) = β ( p i m i ) , i = 1 , 2 , 3 , t > 0 ,
with a 0 = 0 . The system (8) has a zero equilibrium. Take n = 2 , α = 1.3 , and β = 1 . Thus, γ k = 0.649519 , m k = m 3 + k = 1 , k = 1 , 2 , 3 , and assumptions (A1) and (A2) are satisfied. According to Theorem 1, the zero equilibrium is generalized exponential stable. The graphs of the solutions m i ( · ) and p i ( · ) , i = 1 , 2 , 3 , of system (22) are given in Figure 5 (left) and Figure 5 (right), respectively, with n = 2 , ρ = 0.8 , q = 0.3 , α = 1.3 , α 0 = 0 , β = 1 , with initial values m 1 ( 0 ) = 1 , m 2 ( 0 ) = 3 , m 3 ( 0 ) = 5 , p 1 ( 0 ) = 2 , p 2 ( 0 ) = 4 , p 3 ( 0 ) = 6 .
Application 3
Consider the general model describing the dynamics of the interacting defects in the genome and in the proteome with the generalized proportional fractional derivative:
( 0 C D q , ρ u ) ( t ) = p v ( t ) α u ( t ) + f ( t ) , ( 0 C D q , ρ v ) ( t ) = β G K u ( t ) δ v ( t ) + G g ( t ) , t > 0 ,
where β is the coupling rate constant characterizing the regulation of gene expression by the proteins, K is the average number of genes regulated by any single protein and represents a simple measure of the overall connectivity of the genetic network, c reflects the combined efficiency of proteolysis and heat shock response systems, mediating the degradation and refolding of misfolded proteins, respectively, whereas δ characterizes the DNA repair rate, the “force” terms, f ( · ) and g ( · ) characterize the proteome and genome damage rates, respectively, and G is the total genome size.
Let f ( t ) = A e ρ 1 ρ t and g ( t ) = B e ρ 1 ρ t . Then, the model (23) has equilibrium
( u ˜ , v ˜ ) = G B α + A β G K α δ β G K p e ρ 1 ρ t , A δ + G B p α δ β G K p e ρ 1 ρ t , t > 0 .
Model (23) is in the form of (9) with N = 2 , h 1 ( s ) = h 2 ( s ) = s , d 1 = α , d 2 = δ , b = c = 0 . Thus, γ 1 = γ 2 = 1 , μ 1 = μ 2 = 1 and
2 d k + j = 1 N γ j | a k j | + μ N + k μ k c k = 2 α + p < 0 , k = 1 ,
and
2 d k + j = 1 N γ j | a k j | + μ N + k μ k c k = 2 δ + β G K < 0 , k = 2 ,
i.e., assumptions (A1) and (A2) are satisfied if 0.5 p < α and 0.5 β G K < δ . In other words, the DNA repair rate δ and the expressome (proteome, metabolome) turnover rate, c, have to be large enough. In Figure 6, the graphs of the solution ( u , v ) are given with q = 0.3 , ρ = 0.5 , p = 1 , α = 0.9 , A = 1 , β = 2 , K = 1 , δ = 0.8 , G = 0.5 , B = 2 , and the initial values u ( 0 ) = 0.5 , v ( 0 ) = 1.1 . Then, the equilibrium is ( u ˜ , v ˜ ) = 1.9 0.32 e t , 1.8 0.32 e t ) = ( 5.9375 e t , 5.625 e t ) .
Note that the model (23) in the case of ordinary derivatives is studied in [2] with the more restrictive assumption K p β G < c δ .

6. Conclusions

A new gene-regulated model is set up. The dynamics is decribed by the generalized proportional Caputo fractional derivative. The equilibrium is defined in an appropriate way. In the general case, the classical definition of the equilibrium differs. The generalized exponential stability is introduced and studied via the application of Lyapunov functions and their generalized Caputo proportional fractional derivatives. In connection with the application of Lyapunov functions to fractional type models, some properties of absolute values Lyapunov functions and their fractional derivatives are discussed. Several examples are provided to illustrate the properties. The advantages of the application of the quadratic Lyapunov functions are considered, and sufficient conditions for generalized exponential stability are obtained. Some examples illustrate the theoretical results and the dependence of the fractional derivative on the behavior of the solutions.

Author Contributions

Conceptualization, R.A., R.P.A., S.H. and D.O.; methodology, R.A., R.P.A., S.H. and D.O.; validation, R.A., R.P.A., S.H. and D.O.; formal analysis, R.A., R.P.A., S.H. and D.O.; writing—original draft preparation, R.A., R.P.A., S.H. and D.O.; writing—review and editing, R.A., R.P.A., S.H. and D.O.; funding acquisition, S.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Bulgarian National Science Fund under Project KP-06-N32/7.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

S. Hristova is partially supported by the Bulgarian National Science Fund under Project KP-06-N32/7. R. Almeida is supported by Portuguese funds through the CIDMA—Center for Research and Development in Mathematics and Applications, and the Portuguese Foundation for Science and Technology (FCT-Fundação para a Ciência e a Tecnologia), within project UIDB/04106/2020.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jin, Y.; Lindsey, M. Stability analysis of genetic regulatory network with additive noises. BMC Genom. 2008, 9, S21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Kogan, V.; Molodtsov, I.; Menshikov, L.I.; Reis, R.J.S.; Fedichev, P. Stability analysis of a model gene network links aging, stress resistance, and negligible senescence. Sci. Rep. 2015, 5, 13589. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Ren, F.; Cao, F.; Cao, J. Mittag-Leffler stability and generalized Mittag-Leffler stability of fractional-order gene regulatory networks. Neurocomputing 2015, 160, 185–190. [Google Scholar] [CrossRef]
  4. Wu, Z.; Wang, Z.; Zhou, T. Global uniform asymptotical stability for fractional-order gene regulatory networks with time-varying delays and structured uncertainties. Adv. Differ. Equ. 2021, 2021, 93. [Google Scholar] [CrossRef]
  5. Yu, J.; Jiang, C.H.H. α-stability and α-synchronization for fractional-oredr neural networks. Neural Netw. 2012, 35, 82–87. [Google Scholar] [CrossRef]
  6. Ren, F.; Cao, J. Asymptotic and robust stability of genetic regulatory networks with time-varying delays. Neurocomputing 2008, 71, 834–842. [Google Scholar] [CrossRef]
  7. Anbalagan, P.; Hincal, E.; Ramachandran, R.; Baleanu, D.; Cao, J.; Niezabitowski, M. A Razumikhin approach to stability and synchronization criteria for fractional order time delayed gene regulatory networks. AIMS Math. 2021, 6, 4526–4555. [Google Scholar] [CrossRef]
  8. Arjunan, M.M.; Abdeljawad, T.; Anbalagan, P. Impulsive effects on fractional order time delayed gene regulatory networks: Asymptotic stability analysis. Chaos Sol. Fract. 2022, 154, 111634. [Google Scholar] [CrossRef]
  9. Phong, T.T.; Long, L.D. Well-posed results for nonlocal fractional parabolic equation involving Caputo-Fabrizio operator. J. Math. Comput. Sci. 2022, 26, 357–367. [Google Scholar] [CrossRef]
  10. Khan, A.; Trikha, P.; Khan, T. Dislocated dual HPS between integer and complex fractional order chaotic systems using tracking controllers with application. Int. J. Appl. Nonlinear Sci. 2021, 3, 89–118. [Google Scholar] [CrossRef]
  11. Gardner, T.; Cantor, C.; Collins, J. Construction of a genetic toggle switch in Escherichia coli. Nature 2000, 403, 339–342. [Google Scholar] [CrossRef] [PubMed]
  12. Jarad, F.; Abdeljawad, T. Generalized fractional derivatives and Laplace transform. Discret. Contin. Dyn. Syst. Ser. S 2020, 13, 709–722. [Google Scholar] [CrossRef] [Green Version]
  13. Jarad, F.; Abdeljawad, T.; Alzabut, J. Generalized fractional derivatives generated by a class of local proportional derivatives. Eur. Phys. J. Spec. Top. 2017, 226, 3457–3471. [Google Scholar] [CrossRef]
  14. Almeida, R.; Agarwal, R.P.; Hristova, S.; O’Regan, D. Quadratic Lyapunov functions for stability of generalized proportional fractional differential equations with applications to neural networks. Axioms 2021, 10, 322. [Google Scholar] [CrossRef]
  15. Agarwal, R.; O’Regan, D.; Hristova, S. Stability of generalized proportional Caputo fractional differential equations by Lyapunov functions. Fractal Fract. 2022, 6, 34. [Google Scholar] [CrossRef]
  16. Simon, T. Mittag-Leffler functions and complete monotonicity. Integral Transf. Spec. Funct. 2015, 26, 36–50. [Google Scholar] [CrossRef]
  17. Zhang, S.; Yu, Y.; Hu, W. Robust Stability Analysis of Fractional-Order Hopfield Neural Networks with Parameter Uncertainties. Math. Probl. Eng. 2014, 2014, 302702. [Google Scholar] [CrossRef]
  18. Wu, Z.; Wang, Z.; Zhou, T. Global stability analysis of fractional-order gene regulatory networks with time delay. Int. J. Biomath. 2019, 12, 1950067. [Google Scholar] [CrossRef]
  19. Aguila-Camacho, N.; Duarte-Mermoud, M.A.; Gallegos, J.A. Lyapunov functions for fractional order systems. Commun. Nonlinear Sci. Numer. Simul. 2014, 19, 2951–2957. [Google Scholar] [CrossRef]
  20. Duarte-Mermoud, M.A.; Aguila-Camacho, N.; Gallegos, J.A.; Castro-Linares, R. Using general quadratic Lyapunov functions to prove Lyapunov uniform stability for fractional order systems. Commun. Nonlinear Sci. Numer. Simul. 2015, 22, 650–659. [Google Scholar] [CrossRef]
  21. Elowitz, M.B.; Leibler, S. A synthetic oscillatory network of transcriptional regulators. Nature 2000, 403, 335–338. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Graph of the bounds of the solution.
Figure 1. Graph of the bounds of the solution.
Entropy 24 00372 g001
Figure 2. Convergence of the solutions of (19) with q = 0.3 , ρ = 1 to the equilibrium 0.795876 .
Figure 2. Convergence of the solutions of (19) with q = 0.3 , ρ = 1 to the equilibrium 0.795876 .
Entropy 24 00372 g002
Figure 3. Convergence of the solution of (19) with q = 0.3 , ρ = 1 to the equilibrium 2.21938 .
Figure 3. Convergence of the solution of (19) with q = 0.3 , ρ = 1 to the equilibrium 2.21938 .
Entropy 24 00372 g003
Figure 4. Solution of system (19) with q = 0.3 , ρ = 0.8 .
Figure 4. Solution of system (19) with q = 0.3 , ρ = 0.8 .
Entropy 24 00372 g004
Figure 5. Solution of system (22) with q = 0.3 , ρ = 0.8 .
Figure 5. Solution of system (22) with q = 0.3 , ρ = 0.8 .
Entropy 24 00372 g005
Figure 6. Convergence of the solution of (23) to the equilibrium ( 5.9375 e t , 5.625 e t ) .
Figure 6. Convergence of the solution of (23) to the equilibrium ( 5.9375 e t , 5.625 e t ) .
Entropy 24 00372 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Almeida, R.; Agarwal, R.P.; Hristova, S.; O’Regan, D. Stability of Gene Regulatory Networks Modeled by Generalized Proportional Caputo Fractional Differential Equations. Entropy 2022, 24, 372. https://doi.org/10.3390/e24030372

AMA Style

Almeida R, Agarwal RP, Hristova S, O’Regan D. Stability of Gene Regulatory Networks Modeled by Generalized Proportional Caputo Fractional Differential Equations. Entropy. 2022; 24(3):372. https://doi.org/10.3390/e24030372

Chicago/Turabian Style

Almeida, Ricardo, Ravi P. Agarwal, Snezhana Hristova, and Donal O’Regan. 2022. "Stability of Gene Regulatory Networks Modeled by Generalized Proportional Caputo Fractional Differential Equations" Entropy 24, no. 3: 372. https://doi.org/10.3390/e24030372

APA Style

Almeida, R., Agarwal, R. P., Hristova, S., & O’Regan, D. (2022). Stability of Gene Regulatory Networks Modeled by Generalized Proportional Caputo Fractional Differential Equations. Entropy, 24(3), 372. https://doi.org/10.3390/e24030372

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop