Next Article in Journal
Study on the Discharge Process and Mechanism of Anti-Corrosion Pill Particles in the Oil and Gas Field Wellbore Casing Annulus Based on the Discrete Element Method
Next Article in Special Issue
Observer-Based Predefined-Time Attitude Control for Spacecraft Subject to Loss of Actuator Effectiveness
Previous Article in Journal
Fault Diagnosis of Rotating Equipment Bearing Based on EEMD and Improved Sparse Representation Algorithm
Previous Article in Special Issue
Fast Finite-Time Stability and Its Application in Adaptive Control of High-Order Stochastic Nonlinear Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Approximate Method of System Entropy in Discrete-Time Nonlinear Biological Networks

1
College of Mathematics and Systems, Shandong University of Science and Technology, Qingdao 266590, China
2
College of Electrical Engineering and Automation, Shandong University of Science and Technology, Qingdao 266590, China
3
College of Electronic and Information Engineering, Shandong University of Science and Technology, Qingdao 266590, China
4
School of Engineering, QuFu Normal University, Rizhao 276800, China
*
Author to whom correspondence should be addressed.
Processes 2022, 10(9), 1736; https://doi.org/10.3390/pr10091736
Submission received: 29 June 2022 / Revised: 17 August 2022 / Accepted: 23 August 2022 / Published: 1 September 2022
(This article belongs to the Special Issue Advances in Nonlinear and Stochastic System Control)

Abstract

:
This study discusses the calculation of entropy of discrete-time stochastic biological systems. First, measurement methods of the system entropy of discrete-time linear stochastic networks are introduced. The system entropy is found to be characterized by system matrices of the discrete-time biological systems. Secondly, the system entropy of nonlinear discrete-time stochastic biological systems is discussed and is calculated based on a global linearization method. The approximation of the values of system entropy of nonlinear stochastic systems needs to solve an optimization problem that is constrained by a kind of linear matrix inequality (LMI). Finally, a practical biochemical system is provided to verify the effectiveness of the proposed calculation method.

1. Introduction

In general, entropy is considered as a measure of the randomness or disorder of a physical or biological system under intrinsic random fluctuations and environmental disturbances [1,2,3,4,5]. According to the second law of thermodynamics, entropy is used to describe the dispersion of energy in a thermally isolated system, in which energy has a natural tendency to spontaneously change toward states with higher entropy [6,7,8]. In order to maintain life, biological systems need to exchange material and energy with their environment in a continual process, so they are open systems [9].
In this exchange process, the entropy of biological systems can maintain a dynamic balance [10]. Under such a background, the calculation of biological system entropy is particularly important. In the past few decades, the system entropy in biological networks has been extensively studied [11,12,13]. The discrete-time model plays an important role in numerical calculation, stochastic simulation and numerical analysis [14,15,16,17,18]. For continuous-time stochastic systems, the authors of [16] calculated the system entropy of biological systems from their system matrices by the global linearization technique. In this way, the measurement of the system entropy of nonlinear biological networks could be transformed to solve an optimization problem constrained by a set of LMIs. In this paper, we follow the line of [16] and extend the LMI method to the calculation of the discrete-time system entropy of stochastic biological networks. With the aid of the Matlab software package, we solve the corresponding LMI-constrained optimization problem to measure the discrete-time system entropy of a nonlinear biological network.
This paper is organized as follows. In Section 2, we discuss how to calculate the entropy of the discrete-time linear network. Section 3 gives the system entropy measurement of the discrete-time nonlinear random biological network. Section 4 presents how to calculate the system entropy in discrete-time nonlinear stochastic biological networks, which is approximated by the global linearization method. Lastly, in Section 5, an example is given to illustrate the measurement procedure and to validate the feasibility of the proposed system entropy measurement method.

2. System Entropy in Discrete-Time Linear Biological Networks

In this section, we consider a discrete-time linear network, which is described as follows:
x t + 1 = A x t + B v t , y t = C x t , t = 0 , 1 , 2 , T f
where x 0 = 0 , x t R n , v t R m , y t R l denotes the biological network’s state vector, random input and output, respectively. T f is the finite terminal time. A, B and C are matrices with proper dimensions with the following formats:
A = a 11 a 1 n a n 1 a n n , B = b 11 b 1 m b n 1 b n m , C = c 11 c 1 n c l 1 c l n .
The randomness of this system’s output y t can be measured by
r o = E [ 1 T f t = 0 T f y t T y t ]
while the randomness of input signals v t is denoted as
r i = E [ 1 T f t = 0 T f v t T v t ]
where E [ · ] denotes the expectation. Similar to the definitions of entropy in [16], the entropy of the input signal or output signal for a discrete-time system is also defined by
s i = log E [ 1 T f t = 0 T f v t T v t ]
and the entropy of y t is defined by
s o = log E [ 1 T f t = 0 T f y t T y t ] .
Thus, it is natural to define the net signal entropy of a biological system as the discrete-time system entropy, i.e.,
s = s i s o = log E [ 1 T f t = 0 T f y t T y t ] log E [ 1 T f t = 0 T f v t T v t ]
s = log E [ t = 0 T f y t T y t ] E [ t = 0 T f v t T v t ]
Thus, if the system randomness r of System (1) is defined as the following
r = E [ t = 0 T f y t T y t ] E [ t = 0 T f v t T v t ]
the system entropy is represented as
s = log r .
In order to calculate the system entropy s in (1), we have to calculate or approximate the system randomness first. Of course, it is not easy to approximate such randomness directly. Therefore, we need to estimate the system randomness indirectly as follows:
r = E [ t = 0 T f y t T y t ] E [ t = 0 T f v t T v t ] r ¯
which is equivalent to
E [ t = 0 T f y t T y t ] r ¯ E [ t = 0 T f v t T v t ]
Here, r ¯ denotes the upper bound of r.
We will decrease the upper bound r ¯ to be as small as possible, to approach the randomness r of the discrete-time biological network (1), which is suggested in [16].
Proposition 1.
Suppose that a positive definite matrix p > 0 and a positive real number r ¯ > 0 satisfy the following inequality:
A T P A + C T C P + A T P B ( r ¯ I B T P B ) 1 B T P A < 0 r ¯ I B T P B > 0
Then, r ¯ is an upper bound of the system randomness of network (1).
Proof. 
Choose the Lyapunov function V ( x ) = x T P x , then
V ( x t + 1 ) V ( x t ) = ( x t T A T + v t T B T ) P ( A x t + B v t ) = x t T ( A T P A P ) x t + 2 x t T A T P B v t + v t T B T P B v t
Taking summation, and then taking expectation on both sides, we have
E [ V ( x T f + 1 ) ] E [ V ( x 0 ) ] = E { t = 0 T f [ x t T ( A T P A P ) x t + 2 x t T A T P B v t + v t T B T P B v t ] } .
Recalling that V ( x ) 0 , x 0 = 0 and V ( 0 ) = 0 , we have
0 E t = 0 T f [ x t T ( A T P A P ) x t + 2 x t T A T P B v t + v t T B T P B v t ] .
Thus,
E [ t = 0 T f y t T y t ] r ¯ E [ t = 0 T f v t T v t ] E { t = 0 T f [ x t T ( A T P A + C T C P ) x t + 2 x t T A T P B v t v t T ( r ¯ I B T P B ) v t ] } .
Completing the square on the right side, we obtain
E [ t = 0 T f y t T y t ] r ¯ E [ t = 0 T f v t T v t ] E { t = 0 T f [ x t T ( A T P A + C T C P + A T P B × ( r ¯ I B T P B ) 1 B T P A ) x t v t ( r ¯ I B T P B ) 1 B T P A x t r ¯ I B T P B 2 ] } .
where the notation · A 2 denotes y A 2 = y T A y with A > 0 . By inequality (4), there exists
E [ t = 0 T f y t T y t ] r ¯ E [ t = 0 T f v t T v t ] 0 .
Thus, we have
E [ t = 0 T f y t T y t ] r ¯ E [ t = 0 T f v t T v t ]
This shows that the randomness r of System (1) has an upper bound r ¯ . □
By the Schur lemma, we know that the matrix inequality (4) equals the following LMI:
C T C + A T P A P A T P B B T P A B T P B r ¯ I < 0 .
Thus, we can obtain the following corollary described by LMI (5).
Corollary 1.
Suppose that there exists a positive definite matrix p > 0 and a positive real number r ¯ > 0 that satisfy the inequality (5); then, r ¯ is an upper bound of the system randomness of network (1).
Remark 1.
Compared with the results of Reference [16], the matrix inequality (4) that is obtained by the completing square method is different to the results of Proposition 2 of [16]. Moreover, the structure of LMI (5) for a discrete-time network is different to that of (12) for the continuous-time system discussed in [16].
From Equation (3), we can see that the upper bound of the discrete-time system randomness r is r ¯ , i.e., r r ¯ . Thus, the calculation of the system randomness of network (1) can be estimated by the following problem:
r = min r ¯
with the constraint of LMI (5), where p > 0.
Based on (6), with the aid of Matlab’s LMI toolbox, we can decrease r ¯ by the constraint of LMI in (5) until there is no positive matrix P appearing; then, the above LMI-constrained optimization problem can be solved, and the system randomness r can be estimated. Thus, the discrete-time system entropy of the linear stochastic System (1) can be obtained as s = log r . Moreover, the measurement of the system randomness r or system entropy s is dependent on matrices A, B and C in System (1) to some extent.

3. System Entropy in Discrete-Time Nonlinear Network

Nonlinear dynamic systems play an important role in biological networks, which causes the difficulty of estimating the entropy of such systems. Under this situation, the global linearized method is suggested [19,20,21,22,23,24,25,26,27], which is an interpolation method of local linearized systems of a nonlinear biological network. Suppose that the biological systems are described by the following discrete-time nonlinear stochastic biological network:
x t + 1 = f ( x t ) + g ( x t ) v t , y t = h ( x t ) , t = 0 , 1 , 2 , T f
where x 0 = 0 , f is a function, v t denotes m external input signal and g ( x t ) denotes m nonlinear couplings between the biological network and environment. h ( x t ) denotes l nonlinear outputs.
Based on the ideas of the discrete-time system randomness of (3), we obtain the following proposition.
Proposition 2.
Suppose that there exists a positive definite matrix p > 0 and a positive real number r ¯ > 0 that satisfy the following HJI:
f T ( x ) P f ( x ) x T P x + h T ( x ) h ( x ) + f T ( x ) P g ( x ) ( r ¯ I g T ( x ) P g ( x ) ) 1 g T ( x ) P f ( x ) 0 r ¯ I g T ( x ) P g ( x ) > 0 , x R n
Then, r ¯ is an upper bound of the system randomness of network (7).
Proof. 
Let V ( x ) = x T P x , then
V ( x t + 1 ) V ( x t ) = f T ( x t ) P f ( x t ) x t T P x t T + 2 f T ( x t ) P g ( x t ) v t + v t T g T ( x t ) P g ( x t ) v t .
Taking summation, and then taking expectation on both sides, we obtain
E [ V ( x T f + 1 ) ] E [ V ( x 0 ) ] = E { t = 0 T f [ f T ( x t ) P f ( x t ) x t T P x t T + 2 f T ( x t ) P g ( x t ) v t + v t T g T ( x t ) P g ( x t ) v t ] } .
Recalling that V ( x ) 0 , x 0 = 0 and V ( 0 ) = 0 , we have
0 E t = 0 T f f T ( x t ) P f ( x t ) x t T P x t T + 2 f T ( x t ) P g ( x t ) v t + v t T g T ( x t ) P g ( x t ) v t .
Thus,
E [ t = 0 T f y t T y t ] r ¯ E [ t = 0 T f v t T v t ] E { t = 0 T f [ f T ( x t ) P f ( x t ) x t T P x t T + h T ( x ) h ( x ) + 2 f T ( x t ) P g ( x t ) v t v t T ( r ¯ I g T ( x t ) P g ( x t ) ) v t ] } .
Completing the square on the right side, we obtain
E [ t = 0 T f y t T y t ] r ¯ E [ t = 0 T f v t T v t ] E { t = 0 T f [ f T ( x t ) P f ( x t ) x t T P x t T + h T ( x ) h ( x )                 + f T ( x ) P g ( x ) ( r ¯ I g T ( x ) P g ( x ) ) 1 g T ( x ) P f ( x ) v t ( r ¯ I g T ( x ) P g ( x ) ) 1 g T ( x ) P f ( x ) r ¯ I g T ( x ) P g ( x ) 2 ] } .
where the notation · A 2 denotes y A 2 = y T A y with A > 0 . By inequality (8), there exists
E [ t = 0 T f y t T y t ] r ¯ E [ t = 0 T f v t T v t ] 0 .
Thus, we have
E [ t = 0 T f y t T y t ] r ¯ E [ t = 0 T f v t T v t ]
This ends the proof. □
Remark 2.
Compared with the results of Proposition 4 in Reference [16], the HJI (8) in this paper does not depend on the input variables v ( t ) , but the HJIs in Proposition 4 of [16] include v ( t ) . Thus, the system randomness of network (7) can be obtained only by the coefficients f ( x ) , g ( x ) and h ( x ) , which is defined on the state space.
By Proposition 2, the system randomness r of network (7) can be approximated to solve the following optimization constrained by HJI:
r = min r ¯
It is constrained by HJI in (8), where p > 0. Denote
V = V ( x ) | V ( x ) = x T P x , x R n , P > 0 .
The system randomness r of network (7) could be estimated by solving the following optimization:
r = min V ( x ) V r ¯
with the constraint of HJI in (8).
Based on the above results, it is easy to obtain the system entropy as
s = log r .
We can obtain the system randomness r from the HJI-constrained optimization problem in (10) or (11). However, at present, there exists no efficient method to solve the HJI in (10) or (11) analytically or numerically. In this study, the global linearization method in [21,27] will be employed to interpolate several local linearized systems at the M matrices of the convex hull of the globalization systems to approach the nonlinear discrete-time biological systems in (7), to transform the difficult HJI-constrained optimization problem in (10) or (11) to an equivalent LMIs-constrained optimization problem for the calculation of system randomness in the following section.

4. The Global Linearization Method to Estimate System Entropy for Nonlinear Networks

In this section, the global linearization technique is suggested to help in estimating the system entropy of nonlinear stochastic systems in (7). The main idea of this method is described as follows: we convert the nonlinear system into a set of interpolated locally linearized networks in which the linear system’s entropy is easy to calculate and approximate [16,21,27]. Therefore, the estimation of system randomness and entropy can be transformed to solve the HJI-constrained optimization problems (10) efficiently.
We prefer the detailed theory of the global linearization method to Reference [21]. According to this method, the global linearized systems are constructed by the convex hull of M vertices defined in Equation (12) as follows:
f ( x ) x g ( x ) x h ( x ) x C 0 A 1 B 1 C 1 A i B i C i A M B M C M , x t
Then, the state x t in the discrete-time nonlinear System (7) can be represented by those states of local linearized biological networks with (12) as follows:
x t + 1 = A i x t + B i v t , y t = C i x t , t = 0 , 1 , 2 , T
Thus, the combination of linearized systems in (13) can be represented as:
x t + 1 = i = 1 M α i ( x t ) ( A i x t + B i v t ) , y t = i = 1 M α i ( x ) C i x t , t = 0 , 1 , 2 , T
where the interpolation function (for x t x i )
α i ( x t ) = 1 x i x t 2 2 1 = 1 M 1 x i x t 2 2
and α i ( x t ) = 1 for some x t = x i satisfy 0 α i ( x ) 1 and i = 1 M α i ( x t ) = 1 , while x i denotes the ith local operation point with local linearization [27], i.e., the trajectory of the nonlinear System (7) can be represented by the trajectories of the interpolated biological network in (14).
If the nonlinear biological network in Equation (7) could be approximated by the global linearization system in Equation (14), then we obtain the following result.
Proposition 3.
Suppose that a positive definite matrix p > 0 and a real number r ¯ > 0 satisfy the following LIMs:
A i T P A i + C i T C i P A i T P B i B i T P A i B i T P B i r ¯ I 0 , f o r i = 1 , , M
Then, r ¯ is an upper bound of the system randomness of network (7).
Proof. 
Let Lyapunov function V ( x ) = x t T P x t , and the approximation of f ( x t ) and g ( x t ) are
f ( x t ) = i = 1 M α i ( x t ) A i x t
and
g ( x t ) = i = 1 M α i ( x t ) B i .
Then, we have
E [ V ( x 0 ) ] E [ V ( x t ) ] + E [ t = 0 T f ( y t T y t r ¯ v t T v t + Δ V ( x t ) ) ] 0
where
Δ V ( x t ) = i = 1 M α i ( x t ) [ x t T A i T + v t T B i T P ( A i x t + B i v t ) ] x t T P x t T
Thus, we obtain the following result:
E i = 1 M α i ( x t ) x t T v t T п A i T P B i B i T P A i Ξ x t v t 0
where п = A i T P A i + C i T C i P , Ξ = B i T P B i r ¯ I , i = 1 , , M .
This ends the proof. □
Based on this technique, the system randomness r of network (7) or (14) can be approximated by solving the following optimization problem:
r = min P > 0 r ¯
with the constraint of LMIs in (16) and p > 0. The system randomness r in (18) could be obtained by decreasing its upper bound r ¯ until there exists no p > 0.

5. Example and Simulation

In this section, we consider a phosphorelay system in yeast, discussed in [16,23]; see Figure 1. This signal transduction pathway includes seven state variables: Sln1, Sln1H-P, Sln1D-P, Ypd1, Yod1-P, Ssk1, Ssk1-P. We prefer to References [16,23,25], for detailed information.
The state x t is denoted by the following vector:
x t = x 1 ( t ) x 2 ( t ) x 3 ( t ) x 4 ( t ) x 5 ( t ) x 6 ( t ) x 7 ( t ) = S l n 1 ( t ) S l n 1 H P ( t ) S l n 1 D P ( t ) Y p d 1 ( t ) Y p d 1 P ( t ) S s k 1 ( t ) S s k P ( t )
Suppose that the dynamic behavior of this system can be represented by discrete-time difference equations, which are seen as the discrete-time type of phosphorelay system discussed in [25]:
x 1 ( t + 1 ) = ( 1 k 1 ) x 1 ( t ) + k 3 x 3 ( t ) x 4 ( t ) + k 0 v t x 2 ( t + 1 ) = k 1 x 1 ( t ) + ( 1 k 2 ) x 2 ( t ) x 3 ( t + 1 ) = k 2 x 2 ( t ) + ( 1 k 3 x 4 ( t ) ) x 3 ( t ) x 4 ( t + 1 ) = k 4 x 5 ( t ) x 6 ( t ) + ( 1 k 3 x 3 ( t ) ) x 4 ( t ) x 5 ( t + 1 ) = ( 1 k 4 x 6 ( t ) ) x 5 ( t ) + k 3 x 3 ( t ) x 4 ( t ) x 6 ( t + 1 ) = k 5 x 7 ( t ) + ( 1 x 4 ( t ) x 5 ( t ) ) x 6 ( t ) x 7 ( t + 1 ) = ( 1 k 5 ) x 7 ( t ) + x 4 ( t ) x 5 ( t ) x 6 ( t ) y = x 6 ( t )
where k 0 , k 1 , k 2 , k 3 , k 4 and k 5 are the systematic characteristics, and v t is the random fluctuation. Figure 2 shows the trajectories of x with the systematic characteristics k 0 = 0.5 , k 1 = 0.4 , k 2 = 0.1 , k 3 = 50 , k 4 = 50 , k 5 = 0.5 and v t , which is the standard Gaussian white noise with zero mean.
In order to show the effects of systematic characteristics on system entropy, we take the three different systematic characteristics given in Table 1.
Due to the global linearization in (14), the approximation of discrete-time nonlinear network (20) can be presented by:
X t + 1 = t = 1 M α i ( x t ) ( A i X t + B i v t )
where α ( x t ) denotes the the interpolation functions. X t denotes the M states of local linearization systems [26].
By solving the LMI-constrained optimization problem in Equation (18) for the above three system characteristic cases in Table 1, the positive definite matrices P 1 , P 2 and P 3 are obtained as follows:
P 1 = 2.7010 0.4458 0.0177 0.1448 0.0034 0 0 0.4458 2.8184 0.0063 0.0032 0.0006 0 0 0.0177 0.0063 3.0983 0.0006 0.0249 0 0 0.1448 0.0032 0.0006 3.0790 0.0352 0 0 0.0034 0.0006 0.0249 0.0352 3.0920 0 0 0 0 0 0 0 3.9651 0.0001 0 0 0 0 0 0.0001 5.2752 ,
P 2 = 16.2540 0.1773 0.0102 0.0034 0.0000 0 0 0.1773 17.5868 0.0025 0.0001 0.0000 0 0 0.0102 0.0025 16.2441 0.0000 0.0001 0 0 0.0034 0.0001 0.0000 16.2445 0.0001 0 0 0.0000 0.0000 0.0001 0.0001 16.2445 0 0 0 0 0 0 0 17.2443 0.0006 0 0 0 0 0 0.0006 35.0822 ,
P 3 = 119.4032 0.0239 0.0001 0.0005 0.0000 0 0 0.0239 119.6182 0.0007 0.0000 0.0000 0 0 0.0001 0.0007 119.4032 0.0000 0.0240 0 0 0.0005 0.0000 0.0000 119.4032 0.00000 0 0 0.0000 0.0000 0.0000 0.0000 119.4032 0 0 0 0 0 0 0 120.40 0.00 0 0 0 0 0 0.0000 119.78 .
Corresponding system entropy is calculated, and the result is shown in following Table 2.
By the estimation results of Table 2, we see that the different systematic characteristics in (20) can affect the system entropy of the biological network.

6. Conclusions

In this paper, the system entropy measurement of discrete-time nonlinear biological system is discussed. In order to overcome the nonlinear Hamilton-Jacobi inequality(HJI) in the measurement precedure, we extend the global linearization method in continuous-time system to the discrete-time system, so that the HJI-constrained optimization for the measurement of system entropy of discrete-time nonlinear biological system can be transformed to LMIs-constrained optimization problem to efficiently calculate the system entropy easily with the help of LMI Toobox in MATLAB. Moreover, the calculation methods of system entropy of more complex systems such as the nonlinear system with intrinsic randomness, stochastic systems driven by Markov processes are worth further study.

Author Contributions

Conceptualization, X.L.; methodology, X.L. and W.Z.; software, X.L., R.Z. and X.W.; validation, X.L. and R.Z.; formal analysis, X.L. and R.Z.; investigation, W.Z. and C.T.; writing—original draft preparation, X.L., X.W. and R.Z.; supervision, X.L., W.Z. and R.Z.; project administration, X.L. and R.Z.; funding acquisition, X.L., W.Z. and C.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant numbers 61973198, 62173206; the Research Fund for the Taishan Scholar Project of Shandong Province of China; the Natural Science Foundation of Shandong Province of China, grant numbers ZR2020MF062, ZR2019QF005, ZR2021ZD13; the China Postdoctoral Science Foundation, grant number 2021M691849; the SDUST Research Fund, grant number 2015TDJH105.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to thank the anonymous reviewers for their constructive suggestions to improve the quality of this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lebiedz, D. Entropy-related extremum principles for model reduction of dissipative dynamical systems. Entropy 2010, 12, 706–719. [Google Scholar] [CrossRef]
  2. Baierlein, R. Thermal Physics; Cambridge University Press: Cambridge, UK, 1999. [Google Scholar]
  3. Mettetal, J.T.; Oudenaarden, A. Microbiology Necessary noise. Science 2007, 317, 463–464. [Google Scholar] [CrossRef]
  4. Pedraza, J.M.; Oudenaarden, A. Noise propagation in gene networks. Science 2005, 307, 1965–1969. [Google Scholar] [CrossRef]
  5. Mettetal, J.T.; Muzzey, D.; Pedraza, J.M.; Ozbudak, E.M.; Oudenaarden, A. Predicting stochastic gene expression dynamics in single cells. Proc. Natl. Acad. Sci. USA 2006, 103, 7304–7309. [Google Scholar] [CrossRef]
  6. Meirovitch, H. Methods for calculating the absolute entropy and free energy of biological systems based on ideas from polymer physics. J. Mol. Recogn. JMR 2010, 23, 153–172. [Google Scholar] [CrossRef]
  7. Salamon, P.; Ghochani, M.; Nulton, J.; Terrence, G. Shape Entropy and the Time Scales for Thermodynamics in Biological Systems. Biophys. J. 2012, 102, 505. [Google Scholar] [CrossRef]
  8. Mitrokhin, Y. Two faces of entropy and information in biological systems. J. Theor. Biol. 2014, 359, 192–198. [Google Scholar] [CrossRef]
  9. Chen, B.S.; Lin, X.Y.; Zhang, W.H. On the System Entropy and Energy Dissipativity of Stochastic Systems and Their Application in Biological Systems. Complexity 2018, 2018, 1628472. [Google Scholar] [CrossRef]
  10. Baez, J.; Pollard, B. Relative Entropy in Biological Systems. Entropy 2016, 18, 46. [Google Scholar] [CrossRef]
  11. Cofre, R.; Herzog, R.; Corcoran, D.; Rosas, F.E. A Comparison of the Maximum Entropy Principle Across Biological Spatial Scales. Entropy 2019, 21, 1009. [Google Scholar] [CrossRef] [Green Version]
  12. Ana, J.L.; Rigoberto, P.S. Entropy Application for Forecasting. Entropy 2020, 22, 604. [Google Scholar]
  13. Lucia, U. Irreversible entropy variation and the problem of the trend to equilibrium. Physics 2007, 376, 289–292. [Google Scholar] [CrossRef]
  14. Lucia, U. Irreversibility, entropy and incomplete information. Physics 2009, 388, 4025–4033. [Google Scholar] [CrossRef]
  15. Lucia, U. Maximum entropy generation and kappa-exponential model. Physics 2010, 389, 4558–4563. [Google Scholar]
  16. Chen, B.S.; Wong, S.W.; Li, C.W. On the calculation of systrm entropy in nonlinear stochastic biological networks. Entropy 2015, 17, 6801–6833. [Google Scholar] [CrossRef]
  17. Zhang, T.; Deng, F.; Sun, Y.; Shi, P. Fault estimation and fault-tolerant control for linear discrete time-varying stochastic systems. Sci. China Inf. Sci. 2021, 64, 200201. [Google Scholar] [CrossRef]
  18. Jiang, X.; Zhao, D. Event-triggered fault detection for nonlinear discrete-time switched stochastic systems: A convex function method. Sci. China Inf. Sci. 2021, 64, 200204. [Google Scholar] [CrossRef]
  19. Johansson, R. System Modeling and Identification; Springer: London, UK, 1993. [Google Scholar]
  20. Chen, B.S.; Li, C.W. On the Interplay between Entropy and Robustness of Gene Regulatory Networks. Entropy 2010, 12, 1071–1101. [Google Scholar] [CrossRef]
  21. Boyd, S.P.; Ghaoui, L.E.; Feron, E.; Balakrishnan, V. Linear Matrix Inequalities in System and Control Theory; SIAM: Philadelphia, PA, USA, 1994. [Google Scholar]
  22. Chen, B.S.; Wang, Y.C. On the attenuation and amplification of molecular noise in genetic regulatory networks. BMC Bioinform. 2006, 7, 52. [Google Scholar] [CrossRef]
  23. Chen, B.S.; Wang, Y.C.; Wu, W.S.; Li, W.H. A new measure of the robustness of biochemical networks. Bioinformatics 2005, 21, 2698–2705. [Google Scholar] [CrossRef]
  24. Chen, B.S.; Wu, W.S. Robust filtering circuit design for stochastic gene networks under intrinsic and extrinsic molecular noises. Math. Biosci. 2008, 211, 342–355. [Google Scholar] [CrossRef] [PubMed]
  25. Klipp, E.; Herwig, R.; Kowald, A.; Wierling, C.; Lehrach, H. Systems Biology in Practice: Concepts, Implementation and Application; Wiley-VCH: Hoboken, NJ, USA, 2005; pp. 1–465. [Google Scholar]
  26. Chen, B.S.; Chen, P.W. On the estimation of robustness and filtering ability of dynamic biochemical networks under process delays, internal parametric perturbations and external disturbances. Math. Biosci. 2009, 222, 92–108. [Google Scholar] [CrossRef] [PubMed]
  27. Chen, B.S.; Chen, W.H.; Wu, H.L. Robust H2/H global linearization filter design for nonlinear stochastic systems. IEEE Trans. Circ. Syst. 2009, 56, 1441–1454. [Google Scholar]
Figure 1. Schematic representation of phosphorelay system.
Figure 1. Schematic representation of phosphorelay system.
Processes 10 01736 g001
Figure 2. Temporal profiles of the state variables of Sln-1-phosphorelay system in (20) at the systematic characteristic case: k 0 = 0.5 , k 1 = 0.4 , k 2 = 0.1 , k 3 = 50 , k 4 = 50 , k 5 = 0.5 with v t N ( 0 , 0 . 01 2 ) .
Figure 2. Temporal profiles of the state variables of Sln-1-phosphorelay system in (20) at the systematic characteristic case: k 0 = 0.5 , k 1 = 0.4 , k 2 = 0.1 , k 3 = 50 , k 4 = 50 , k 5 = 0.5 with v t N ( 0 , 0 . 01 2 ) .
Processes 10 01736 g002
Table 1. Systematic characteristics in three cases.
Table 1. Systematic characteristics in three cases.
Characteristics k 0 k 1 k 2 k 3 k 4 k 5
   Case 10.50.40.150500.5
   Case 20.30.10.2110.6
   Case 30.010.020.0310100.04
Table 2. System entropy of systems with three different systematic characteristics.
Table 2. System entropy of systems with three different systematic characteristics.
CharacteristicsPr s = log r
Case 1 P 1 3.64041.2921
Case 2 P 2 20.3077 3.0110
Case 3 P 3 149.25395.0056
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lin, X.; Wang, X.; Zhang, W.; Zhang, R.; Tan, C. An Approximate Method of System Entropy in Discrete-Time Nonlinear Biological Networks. Processes 2022, 10, 1736. https://doi.org/10.3390/pr10091736

AMA Style

Lin X, Wang X, Zhang W, Zhang R, Tan C. An Approximate Method of System Entropy in Discrete-Time Nonlinear Biological Networks. Processes. 2022; 10(9):1736. https://doi.org/10.3390/pr10091736

Chicago/Turabian Style

Lin, Xiangyun, Xinrui Wang, Weihai Zhang, Rui Zhang, and Cheng Tan. 2022. "An Approximate Method of System Entropy in Discrete-Time Nonlinear Biological Networks" Processes 10, no. 9: 1736. https://doi.org/10.3390/pr10091736

APA Style

Lin, X., Wang, X., Zhang, W., Zhang, R., & Tan, C. (2022). An Approximate Method of System Entropy in Discrete-Time Nonlinear Biological Networks. Processes, 10(9), 1736. https://doi.org/10.3390/pr10091736

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