Abstract

This study discusses the relationship between the entropy and the dissipativity of stochastic systems under the background of biological systems. First, measurement methods of the system entropy and energy dissipativity of linear stochastic biological systems are introduced. We found that the system entropy is negatively proportional to the energy dissipativity in logarithmic scale. Some opposite effects between system entropy and energy dissipativity are also discussed and compared based on their measured values to get insight into the understanding of the system mechanisms and the system characteristics. We found that the intrinsic random fluctuation and the enhancement of the system robust stability both can increase the system entropy but decrease the system dissipativity. The system entropy and the energy dissipativity of nonlinear stochastic biological systems are also discussed and compared based on a global linearization method. Computation methods are also provided. Finally, two numerical examples are demonstrated to verify theoretical prediction.

1. Introduction

In general, entropy is considered as a measurement of randomness of thermodynamic systems [14]. According to the second law of thermodynamics, entropy is used to describe the dispersion of energy and the natural tendency of spontaneous change toward states with higher entropy in a thermally isolated system [58]. Since biological systems can exchange material and energy with their environment to maintain life, they are open systems, i.e., their entropy can be maintained through exchanging material and energy with their environment [915]. Recently, system entropy of biological system has been introduced as the net signal entropy between input signal entropy and output signal entropy [9, 10]. Further, biological systems need energy consumption to maintain their mechanisms [1622]. The power dissipation denotes the power exhaustion of biological system [2326]. Nonnegative dissipation corresponds to the passivity of biological system, and negative dissipation corresponds to the activity of biological system [23, 26]. In general, a biological system is passive (nonnegative dissipation), if the biological system has more energy dissipation (heat output) due to passive components than energy generation due to active components (i.e., net energy dissipation [26]). A biological system is with negative dissipation (i.e., active) if the biological system has more energy generation than energy dissipation (i.e., net energy generation) [13, 14, 26]. The largest dissipation of the biological system is called dissipativity [1, 26]. Obviously, there is some relationship between system entropy and dissipativity of biology systems, which is valuable to investigate from the systematic perspective.

In this study, we first measure the system entropy and energy dissipativity of a more general linear biological system. We could obtain system entropy by solving a linear matrix inequality (LMI), constrained minimization problem and energy dissipativity by solving a LMI, constrained maximization problem. We found that both system entropy and energy dissipativity are both system characteristics dependent on system parameters of biological system. By comparing two LMI constraints for system entropy and energy dissipativity and their corresponding Riccati–like inequalities, we found that system randomness and system energy dissipativity are approximately related by for some constant and system entropy and system energy dissipativity are related by as shown in Figure 1. We found that a more stable (with eigenvalues in the far left s–domain) biological system is with smaller system randomness and entropy, and with a nonnegative dissipativity. However, when biological system is less stable, the system randomness and entropy become larger and the energy dissipativity becomes negative dissipation (energy generation). When the coupling matrices between the biological system and external environment become weaker, the system randomness, system entropy, and system energy dissipation become smaller and vice versa. In nonlinear biological system, we need to solve two Hamilton Jacobi inequalities (HJI), constrained optimization problems for calculating system entropy and energy dissipativity, respectively [27]. In general it is very difficult to solve these two HJI–constrained optimization problems for system entropy and energy dissipativity of nonlinear biological systems except very simple systems [2834]. Moreover, such difficulties still exist in nonlinear system’s stability theories [3537] and other topics [3841]. In this study, the global linearization method is introduced to interpolate several linearized local biological systems to approximate the nonlinear biological system [23, 42]. In this way, the nonlinear HJI could be interpolated by a set of LMIs. In this situation, the HJI-constrained optimization problems for solving the system entropy and energy dissipativity of nonlinear biological system could be replaced by the corresponding LMI-constrained optimization problems, which could be efficiently solved with help of LMI toolbox in Matlab [23, 4345]. Based on those techniques, the relationship between system entropy and energy dissipativity can be easily investigated and discussed.

The contributions of this paper are described as follows: we investigate the relationship between two important systematic characters, i.e., system entropy and dissipativity of stochastic system. We found that the system entropy is negatively proportional to the energy dissipativity in logarithmic scale. In order to overcome the difficulty of the measurement of the system entropy and energy dissipativity of nonlinear biological systems given by the HJI-constrained optimization method, the global linearization technique is also employed to interpolate a nonlinear system by a set of local linearized systems such that system entropy and energy dissipativity could be measured by the LMIs-constrained optimization method. Then the relationship between system entropy and energy dissipativity could be also investigated by solving the corresponding LMIs. Finally, two in silico examples of a phosphorelay biological system and a predator-prey ecological system are given to illustrate how to measure the system entropy and energy dissipativity and to confirm their relationship investigated by the proposed method.

2. Systems Entropy and Energy Dissipativity of Linear Biological Systems

In this section, we consider a more general linear biological system with the following stoichiometric equation [9, 26]:where denotes the state vector of biological system with species, the random signal denotes the system input which can be drawn from any distribution with zero moment and finite variance for every , and denotes the system output. , , , and , respectively, denote the system interaction matrix, input state coupling matrix, state/output coupling matrix, and input/output coupling matrix of the linear biological system as follows:

2.1. System Entropy of Linear Biological System

The randomness of output and input signals can be considered as the dispersion of random signals; i.e., the randomness of input signal is represented by [9, 10]when the randomness of output signal is represented by [9, 10]where denotes the expectation of and denotes the period of the corresponding signal for trajectories with and . The randomness in (3) and (4) is an estimator of the variance of random signal.

The entropy of input signal is denoted as [911]and the entropy of output signal is denoted as [911]Then the system entropy of linear biological network in (1) is defined as the difference between input signal entropy and output signal entropy with zero initial state condition , i.e., the net signal entropy [9, 10]

The system entropy denotes the entropy increment between input signals and output signals. In this process, by the second law of thermodynamics, the dynamical system (1) should provide/absorb enough energy to achieve such transformation [46]. Similarly, when system entropy decreases/increases, the systems should absorb/release energy from/to the environment to maintain the system entropy at a steady state. For the biological systems, these processes are achieved by exchanging energy with external environment.

Let us denote the system randomness as [9, 10]Then the system entropy of linear biological system in (1) is the logarithm of the system randomness as [9, 10]

2.2. System Dissipativity of Linear Biological System

Suppose the dimensions of and are equivalent. The linear stochastic biological system in (1) is said to have energy dissipation if [23, 26]holds for every stochastic input and the corresponding process with and all .

The energy dissipation denotes the energy exhaustion of biological system. The largest dissipation of the biological system, i.e., the largest number such that (10) holds, will be called its dissipativity .

Remark 1. The biological system in (1) is said to be passive if energy solution with satisfies [23, 26] i.e., passivity corresponds to nonnegative dissipation (i.e., in (10)). The nonnegative dissipation () corresponds to the passivity of the biological system, and negative dissipation () corresponds to the activity of the biological system. In general, a biological system is passive () if the biological system has more energy dissipation due to passive components than energy generation due to active components (net energy dissipation). On the other hand, biological system is of negative dissipation () if the biological system has energy generation than energy dissipation (net energy generation). Obviously, there is some relationship between system entropy and dissipativity needful to be investigated for biological systems.

2.3. Measurement of System Entropy and Energy Dissipativity of Linear Biological Systems

From fractional functions in (7)-(9), it is difficult to calculate the system entropy for the linear biological system in (1) directly. An indirect method is proposed to measure system entropy in (9) through system matrices , , , and of the biological system. Let us denote the upper bound of system randomness as follows:or equivalentlyWe will first estimate the upper bound of system randomness in (12) and then decrease as small as possible to approach the system randomness of biological system.

Proposition 2. If the linear matrix inequality (LMI) has a positive definite solution then the system randomness of linear biological system has an upper bound .

Proof. See Appendix A.

Based on Proposition 2, the system randomness of biological system (1) could be estimated by solving the following LMI-constrained optimization problem:The above LMI-constrained optimization problem for system randomness could be easily solved by decreasing in (14) until there exists no positive definite symmetric , which is easily performed with the LMI toolbox available in Matlab software package. For further discussion, we prefer to reference [23]. After solving from (15), we could obtain system entropy from (9).

For the energy dissipation in (10), we get the following result.

Proposition 3. If the LMI has a positive definite solution then the biological system in (1) has energy dissipation as shown in (10).

Proof. See Appendix B.

Since the energy dissipativity is the largest number of energy dissipation [23, 26], the energy dissipativity of linear biological system in (1) is measured by solving the following LMI-constrained optimization problem:In the constrained maximization of , to obtain energy dissipativity in (17), we could increase until there exists no in the LMI of (16), and we could obtain . The LMI-constrained maximization in (17) could be also solved with the help of LMI toolbox in Matlab software package.

Remark 4. (i) In the case of strict linear biological system (i.e., without direct coupling matrix ) [9] the LMI of system randomness of biological system in (18) is reduced from (14) to the following LMI [9]:In this situation, the system randomness of biological system in (18) is solved as follows [9]:Similarly, the LMI of energy dissipation of biological system in (18) is reduced from (16) to the following LMI:In this situation, the system energy dissipativity of biological system in (18) is solved as(ii) After solving from (15), by the Schur complement, the LMI of system randomness in (14) of biological system in (1) is equivalent toSimilarly, after solving from (17), the LMI of system dissipativity in (16) of biological system (1) is equivalent to In the case of strict biological system in (18), the inequality of system randomness in (24) is reduced to [9]and the inequality of energy dissipativity in (24) is reduced to

From the analysis in Remark 4, it is seen that the relationship between the system randomness and the energy dissipativity has the following proportion approximately:orfor some constant .

According to (9), the relationship between system entropy and the energy dissipativity is approximately given byfor some constant .

Remark 5. (i) In Proposition 2, in order to guarantee the negative definite of LMI in (14) all the diagonal components should be negative definite. In this situation, the lower bound of system randomness of the biological system in (1) is given byTherefore the constant has the following lower bound(ii) In Proposition 3, in order to guarantee the negative definite of LMI in (16), the upper bound of system dissipativity of the biological system in (1) is given by(iii) Since there are some extra-terms in the left hand side of (26) which do not appear in (25), the relation in (27) or (28) holds only if the equalities of (30) and (31) are satisfied for some positively definite matrix .

3. System Entropy and Energy Dissipativity of Nonlinear Biological Systems

Biological systems are always nonlinear stochastic systems. Therefore, the biological dynamical system in (1) should be modified as the following nonlinear stochastic system where denotes the nonlinear interaction vector among species within the biological system, denotes the nonlinear coupling vector between the biological dynamic system and environmental disturbance , denotes nonlinear outputs, and denotes the direct coupling vector from to .

3.1. System Entropy of Nonlinear Biological Systems

We assume the equilibrium point of interest (i.e., the phenotype to be discussed) of nonlinear biological system in (33) is at , i.e. . If the equilibrium point is of interest (see Figure 2), the origin of nonlinear biological system must be shifted to for the convenience of analysis. Based on the above assumptions, we get the following result of system randomness in (12) for nonlinear biological system in (33).

Proposition 6. The system randomness of nonlinear biological system in (33) has an upper bound if the following Hamilton Jacobi Inequality (HJI) has a positive Lyapunov solution ()

Proof. See Appendix C.

Based on the result of Proposition 6, the system randomness of nonlinear biological system in (33) could be estimated by solving the following HJI-constrained optimization problemAfter measuring system randomness from solving the above HJI-constrained optimization problem, we could estimate system entropy for the nonlinear biological system in (33).

3.2. Energy Dissipativity of Nonlinear Biological Systems

For the energy dissipation in (10) of nonlinear biological system in (33), we get the following result.

Proposition 7. If the HJI has a positive solution then the nonlinear biological system in (33) has energy dissipation as shown in (10).

Proof. See Appendix D.

Since the energy dissipativity is the largest number of energy dissipation , the energy dissipativity of nonlinear biological system in (33) could be measured by solving the following HJI-constrained optimization problem:

Remark 8. (i) After solving the system randomness of nonlinear biological system (33) from (35), then the HJI of system randomness in (34) should be modified as follows:Similarly, after solving the energy dissipativity of nonlinear biological system (33) from (37), the HJI of energy dissipativity in (36) should be modified as follows:(ii) By comparing two HJIs of system randomness and energy dissipativity in (38) and (39), we found that the relationships among system randomness, system entropy, and energy dissipativity in (27), (28), and (29) are still approximately true near the equilibrium point .
(iii) If the nonlinear biological system is more stable at the equilibrium in Figure 2 (i.e., the basin of the equilibrium point is more deep and sharp), in this situation, the term in (38) is more negative and could be smaller to maintain in the negative inequality of (38), i.e., a more stable nonlinear biological system is with a smaller system randomness. Further large couplings , , and in (38) will lead to large system randomness and vice versa.
(iv) From the inequality of energy dissipativity in (39), if nonlinear biological system is more stable at (i.e., the basin of in Figure 2 is more deep and sharp and is more negative at ), the energy dissipativity is more nonnegative. When the biological system is less stable at , the energy dissipativity become less nonnegative or positive (i.e., active with energy generation).

3.3. Measurement of System Entropy and Energy Dissipativity of Nonlinear Biological Systems

In general, it is very difficult to solve the HJI-constrained optimization problems in (35) and (37) for system randomness and energy dissipativity of nonlinear biological system in (33), respectively. In this study, the global linearization technique will be employed to interpolate several local biological systems to approximate a nonlinear biological system so that the HJI in (34) or (36) could be interpolated by a set of LMIs to simplify the solution of HJI-constrained optimization problem in (30) or (37).

If the global linearized systems of nonlinear biological system in (33) are bounded by a polytope consisting of the vertices [23, 42]then the state trajectory of nonlinear biological system in (33) can be represented by the convex combination of state trajectories of the following local linearized biological systems derived from the vertices of the polytope in (40): According to the global linearization theory in [23, 42], every trajectory of the nonlinear biological system in (33) can be interpolated by the following local linearized biological systems: where the interpolation functionswith interpolation nodes , satisfy and .

Remark 9. The linear biological system in (1) could be considered as a special linearization case of nonlinear biological system (33) at the origin . Therefore the trajectories of the nonlinear biological system in (33) can be represented by the global linearization biological system in (42), i.e., the convex combination of local linearized biological systems at the vertices of polytope in (40).
After the nonlinear biological system in (33) is represented by the global linearization system in (42), we get following results for nonlinear biological system in (33) or (42).

Proposition 10. If the linear matrix inequalities (LMIs) have a positive definite solution then the system randomness of nonlinear biological system in (42) or (33) has an upper bound .

Proof. See Appendix E.

Therefore, according to Proposition 10, the system randomness of nonlinear biological system in (42) or (33) could be calculated by solving the following LMIs-constrained optimization problem:

Similarly, a systematic characteristic of energy dissipativity in (10) of nonlinear biological system in (33) is derived as follows.

Proposition 11. If the LMIs have a positive definite solution then the nonlinear biological system in (33) or (42) has energy dissipation .

Proof. See Appendix F.

Therefore, the energy dissipativity of nonlinear biological system in (33) or (42) could be calculated by solving the following LMIs-constrained optimization method:After solving the LMIs-constrained optimization problems in (45) and (47) for system randomness and energy dissipativity for nonlinear biological system in (33) or (42), LMIs in (44) for system randomness and LMIs in (46) for energy dissipativity are given in the following, respectively: which are equivalent to the following Riccati-like inequalities, respectively,which could be considered as local linearizations of the HJIs in (34) and (36), respectively; i.e., the interpolation of local Riccati-like inequalities in (50) and (51) could approach HJIs in (34) and (36), respectively.

Remark 12. (i) Comparing (48) and (50) with (14) and (24), it is seen that the system randomness of nonlinear biological system in (33) or (42) should satisfy LMIs (48) or Riccati-like inequalities (50) of local linearized biological systems in (40) or (42). Similarly, comparing (49) and (51) with (24) and (24), the energy dissipativity of nonlinear biological system in (33) or (42) should satisfy LMIs in (49) and Riccati-like inequalities in (51) of local linear biological systems in (40) or (42).
(ii) Comparing (48) and (50) with (49) and (51), it is seen that the relationships between system random, system entropy, and energy dissipativity in (27), (28), and (29) are also true near the equilibrium point .

4. Example of Calculating System Entropy and Energy Dissipativity of Biological Networks

In this section two examples of calculating the system entropy and the energy dissipativity of biological networks are given to illustrate the relationships between them.

4.1. Example 1: The System Entropy and Energy Dissipativity of a Phosphorelay System

A phosphorelay system from high osmolarity glycerol (HOG) signal transduction pathway in yeast is shown in Figure 3. This signal transduction pathway involves a transmembrane protein Sln1, which is a dimmer. Under normal conditions, the signal transaction pathway is active due to continuous autophosphorylation at a histine residue, sln1H-p, under the consumption of ATP. Then, this phosphate group is transferred to an aspartate group of Sln1, i.e., Sln1D-p, consequently to a histidine residue of Ypd1, and to an aspartate residue of Ssk1. Finally, the Ssk1-p is continuously dephosphorylated by phosphatase in the signal transduction pathway [20].

Let us denote the state vector of the signal transduction pathway as Then the dynamic system of the Sln1-phosphorelay signal pathway in yeast can be represented by the following dynamic equation: where is input signal, i.e., high osmolarity.

In order to investigate the relationship between system entropy and energy dissipativity of the phosphorelay system in (28) and (29), the following three systematic characteristics around the equilibrium point are considered to measure their system entropy and dissipativity:(a)Case 1: , , , , , and ;(b)Case 2: , , , , , and ;(c)Case 3: , , , , , and .

The temporal profiles of the phosphorelay system at the above three systematic characteristic cases are given in Figure 4.

Based on the proposed calculation procedure, we could measure the system entropy and energy dissipativity of the phosphorelay system. According to the global linearization in (42), the phosphorelay system in (53) could be represented by the following global linearization system:where the local linearized system matrices and for Case 1, Case 2, and Case 3 are given in Appendix G.

By solving the LMI-constrained optimization problem in (45) at the above three system characteristic cases with we could obtain the system randomness and the corresponding system entropy at three systematic characteristic cases as follows:

Similarly, by solving the LMIs-constrained optimization problem in (47) at three system characteristic cases with we could obtain the system energy dissipativity at three systematic characteristic cases as follows:

According to the above measurement results of system entropy and energy dissipativity at three systematic characteristic cases, with some interpolation technique, the relationship between of system entropy and system energy dissipativity of the phosphorelay system in Equation (53) is shown in Figure 5 which could confirm the relationship in Figure 1(b) between the system entropy and the energy dissipativity near the equilibrium state of nonlinear biological system. The simulations of the output for Case 1, Case 2, and Case 3 are shown in Figure 6. Obviously, the measurements of phosphorelay system can confirm the predicted relationship of system entropy and system energy dissipativity of biological system in (29) and in Figure 1.

4.2. Example 2: The System Entropy and Energy Dissipativity of Predator-Prey Ecological System

In the predator-prey ecological system [21, 4649] where and are external inputs, from (42), the nonlinear predator-prey system in (59) was approximated by the following global linearization system:where and , for three systematic characteristic cases of , , and given in Appendix H. The simulations of the predatory-prey ecological system in (59) with three systematic characteristic cases are shown in Figure 7, solving the LMIs-constrained optimization problems in (45) and in (47) with The system entropy and energy dissipativity of predatory-prey ecological system in (59) at the three systematic characteristic cases of , , and are obtained, respectively, The relationship between system entropy and dissipativity based on three systematic characteristic cases is shown in Figure 8 which could also confirm the relationship in Figure 1(b) between the system entropy and the energy dissipativity near the equilibrium state of biological system. The simulations of the first component of the output for three systematic characteristic cases , , and are shown in Figure 9. Obviously the measured results of predator-prey ecological system can also confirm the predicted relationship of system entropy and system energy dissipativity of biological system in (29) and Figure 1.

5. Conclusion

In this study, the measurement methods of system entropy and energy dissipativity of more general biological system are introduced and the relationship between system entropy and energy dissipativity is also investigated by as shown in Figure 1(b). In the nonlinear biological systems, based on the global linearization technique, we could find that the same relationship between system entropy and energy dissipativity is also maintained at the phenotype near the equilibrium point . From the two simulation results, the relationship between the system entropy and energy dissipativity is confirmed by two nonlinear biological systems. Our theories verify that when biological systems exchange material and energy with their environment to maintain life, their entropy and energy dissipativity can be maintained through exchanging material and energy with their environment. Since the evolution of biological system is also related to the intrinsic genetic variations and environmental disturbances, our future work will be to investigate the relationship between the system entropy and evolution of biological systems. Furthermore, the experimental results in literature for the validation of the proposed theories are also important which also can be considered in the future studies.

Appendix

A. Proof of Proposition 2

Proof. From (13), we getor equivalentlywhere the Lyapunov function for some positive symmetric matrix . Since we assume and , the following inequality implies (A.1) and (A.2):By in (1), we getTherefore, if the LMI in (14) holds, then the inequality (A.4) holds and (A.1) also holds.

B. Proof of Proposition 3

Proof. From (10), we getor equivalentlyBy the fact and and , the following inequality implies the inequality in (B.2):By the fact in (1), we getBy the fact and , (B.4) is equivalent toTherefore, the LMI in (16) implies (B.5) and (B.1).

C. Proof of Proposition 6

Proof. From (12), we getor equivalently, for some Lyapunov energy function , when By the fact , and , the following inequality implies (C.2):Therefore, the HJI in (34) implies (C.3), and (C.3) implies (C.1); i.e., HJI in (34) implies that the system randomness of nonlinear biological system in (33) has an upper bound in (12).

D. Proof of Proposition 7

Proof. From (10), we getor equivalently, for some Lyapunov energy function , when By the fact and and , the following inequality implies the inequality in (D.2) or (D.1):The HJI in (36) could imply the inequality in (D.3) and then implies (D.2) and (D.1). Therefore the HJI in (36) implies the energy dissipation in (10).

E. Proof of Proposition 10

Proof. From (13), we getor equivalentlywhere the Lyapunov function for some positive symmetric matrix satisfying LMIs (46). Since and , the following inequality implies (E.1) and (E.2):Keeping (42) in mind, we get orTherefore, if the LMI in (44) holds, then the inequality (E.5) holds and (E.1) also holds.

F. Proof of Proposition 11

Proof. From (10), there existsor equivalentlyTaking and by the fact and and , the following inequality implies the inequality in (F.2):By (42), we getBy the fact and , (F.4) is equivalent toTherefore, the LMI in (46) implies (F.5) and (F.1).

G. System Matrices and of Local Linearized Systems in (40) Used in Example 1

In system (53), the corresponding nonlinear coefficients and are presented as In order to linearize system (53), let , and as the approximation variables. DenoteThen with and when takes the fixed approximating values.

The matrices and for Case 1, Case 2, and Case 3 in Example 1 are calculated based on (G.2) where takes values in , , and , respectively, .

H. System Matrices and of Local Linearized Systems in (40) Used in Example 2

In system (59), the corresponding nonlinear coefficients and are presented as In order to linearize system (59), let and as the approximation variables. DenoteThen with and when takes the fixed approximating values. The matrices and for , , and used in Example 2 are calculated based on (H.2), where takes values in , .

Denote , , and , and letwhere and takes values in with , takes values in with . Then the global linearization system (60) for the nonlinear system (59) is presented as

Figure 10 illustrates the relationships between and with different , ; i.e., can also be seen as the function of ( is omitted) and is also written as . It shows that, for each , the closer the distance between and the bigger the values of are, inversely, the further the distance between and the smaller the values of are.

Figure 11 illustrates the positional relationships among . It shows that for every given there exist several with different values corresponding with it and, by (H.3), they also satisfy .

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (no. 61573227), the Research Fund for the Taishan Scholar Project of Shandong Province of China, the Natural Science Foundation of Shandong Province of China (no. ZR2016FM48), SDUST Research Fund (no. 2015TDJH105), and the Ministry of Science and Technology of Taiwan (under contract no. MOST-106-2221-E-007-010-MY2).