Abstract

This paper is concerned with positive stability analysis and bio-circuits design for nonlinear biochemical networks. A fuzzy interpolation approach is employed to approximate nonlinear biochemical networks. Based on the Lyapunov stability theory, sufficient conditions are developed to guarantee the equilibrium points of nonlinear biochemical networks to be positive and asymptotically stable. In addition, a constrained bio-circuits design with positive control input is also considered. It is shown that the conditions can be formulated as a solution to a convex optimization problem, which can be easily facilitated by using the Matlab LMI control toolbox. Finally, a real biochemical network model is provided to illustrate the effectiveness and validity of the obtained results.

1. Introduction

In the past decades, biochemical networks, such as metabolic networks [1] and genetic networks [2], have received considerable attention and become a hot research topic [35]. A great number of results have been obtained, such as gene expression data modeling [69] and dynamic analysis of biochemical networks [10, 11].

It is not surprising that dynamical system theory plays a central role in understanding biological and physiological processes [4, 12] since it provides a powerful tool to quantitatively analyze these biochemical networks from a systematic viewpoint. In addition, dynamical system theory is important and useful for the development of synthetic and systems biology, which has a great potential in gene therapy and drug design [13].

It is worth noting that dynamical models of many biological and physiological processes, such as metabolic systems and endocrine systems [1], biochemical reactions [12], are derived from mass and energy balance considerations that take nonnegative chemical concentrations as dynamic states. Hence, state trajectories of such biochemical systems remain in the nonnegative orthant of the state space for arbitrary nonnegative initial conditions. Such systems are commonly referred to as nonnegative systems or positive systems in the literature [1417]. In this paper, we call them positive systems for convenience. A subclass of positive dynamical systems are compartmental systems [18], which involve dynamical models that are characterized by conservation laws (e.g., mass and energy) capturing exchange of materials between coupled macroscopic subsystems known as compartments. Each compartment is assumed to be kinetically homogeneous. That is, any material entering the compartment is instantaneously mixed with materials of the compartment. There have been some studies on this kind of systems with application to biochemical networks [4, 19].

Recently, there have been some results on stability analysis of gene networks with some special regulation functions [2025]. However, due to nonlinearity and complexity of biochemical networks, there does not exist a systematic approach to stability analysis of such nonlinear biochemical networks. It is noted that the fuzzy interpolation approach proposed in [26] can be seen as a promising way of dealing with nonlinear complex systems. In [27], the authors firstly used the fuzzy approximation method to investigate the robust stability of stochastic biochemical regulatory networks, where, however, the positive constraint of network states was ignored.

On the other hand, there have been some other results on bio-circuits design for biochemical networks. In [28], an external optimal control input was applied to stabilize a gene regulatory system. In [29], a simple robust circuit has been designed for the S-system model without considering stochastic noises. In [30], a robust engineering principle was proposed for stochastic biochemical regulatory networks with parameter uncertainties and disturbances. However, most of the aforementioned results paid little attention to the positive constraint of the states of biochemical networks, not to mention constrained bio-circuits with positive control input. In many cases, the states of those studied systems would become negative in simulation, which could be seen as the main drawbacks of these results.

In this paper, the T-S fuzzy system is employed to approximate nonlinear biochemical networks by interpolating several local linear systems. A positive stability condition and a bio-circuits design procedure will be developed for the nonlinear biochemical networks. In addition, a constrained bio-circuits design with positive control input will also be considered. Finally, a real biochemical network is given to illustrate the effectiveness of the obtained results.

The rest of the paper is organized as follows. In Section 2, some useful definitions and lemmas for positive systems are introduced. In Section 3, by using the fuzzy approximation approach, a sufficient condition for positive stability will be derived for nonlinear biochemical networks. In Section 4, bio-circuits design will be developed from a systematic point of view. In addition, a constrained bio-circuits design with positive control input will also be considered. In Section 5, a real biochemical network is given to illustrate the effectiveness of the obtained results. Finally, the paper will be closed with a conclusion.

2. Notation and Preliminaries

Notation. denotes the set of real numbers, stands for the vector space of all -tuples of real numbers, and is the space of matrices with real entries. For in , denotes the th component of . denotes the sets of all real matrices with nonnegative entries and . For a real matrix , () means that all its entries are nonnegative (positive). is the transpose of , and is the inverse of . The notation means that is symmetric and positive definite. The following notations of matrices are used throughout this paper: , .

Consider a general nonlinear system where is the state vector and is the control input. The nonlinear function satisfies and ; that is, has the second-order continuous derivative with respect to and .

Firstly, some definitions and useful lemmas for positive nonlinear systems are given as follows.

Definition 1. Given any positive initial condition , the unforced nonlinear system (1) is said to be positive if the corresponding trajectory for all .

Definition 2. Let the nonlinear function f is positive (or essentially nonnegative in [19]) if , for all and .

Lemma 3 (see [19]). Consider the unforced nonlinear system (1). If is positive and continuously differentiable in and , then is positive (essentially nonnegative).

Lemma 3 implies that if a nonlinear system is positive, then its linearization is also positive.

Theorem 4. For the unforced system (1), if the nonlinear function is positive, and if there exists a Lyapunov function and satisfying the following inequality: for all nonzero , then the equilibrium point of the nonlinear system (1) is asymptotically stable.

As a special case of the nonlinear system (1), the following linear system is considered where and . Some useful results from [14] are presented as follows.

Definition 5. A real matrix is called a Metzler matrix if its off-diagonal elements are nonnegative, that is

Lemma 6. The unforced linear system (3) is positive if and only if is a Metzler matrix.

Lemma 7. The unforced positive linear system (3) is asymptotically stable if and only if there exists a positive definite diagonal matrix such that

Remark 8. It follows from the physical consideration that all the states of biochemical networks should stay in the positive orthant. Thus, positive systems are suitable for quantitatively studying biochemical networks, such as gene microarray data modeling [68] and bio-circuits design for biochemical networks [3, 4].

3. Positive Stability Analysis of Nonlinear Biochemical Networks

As pointed out in Introduction, many applications in biochemical processes give rise to nonlinear dynamical systems, such as genetic networks, metabolic pathways and membrane transports, to cite just a few examples. Consider the nonlinear system (1) for representation of biochemical networks, which describes complex interactions between molecules. It is noted that the states of the system denote the concentrations of the molecules, the nonlinear function is the regulation function and the control input can be seen as the external sources, such as drugs, proteins, or other chemical complexes.

For a general nonlinear biochemical network (1), it is very difficult if not impossible to find a suitable Lyapunov function such that the condition (2) is satisfied, especially when the positive constraint of the states should be also maintained. However, the T-S fuzzy interpolation approach provides a way to approximate the nonlinear biochemical network and, thus, potentially provides a simplified method for positive stability analysis and bio-circuits design.

Consider a nonlinear biochemical network described by a T-S fuzzy system where and is the number of fuzzy rules; are the premise variables and are the fuzzy sets; and are the control input and output, respectively; , and are the known matrices of appropriate dimensions.

By using a center average defuzzifier, product inference, and a singleton fuzzifier, the global dynamics of the T-S fuzzy system (6) cab be described by where ’s are the so-called normalized activation functions in relation to the th submodel such that

Now, we are in the position to develop the global positive stability results for the unforced biochemical network (7); that is, .

Theorem 9. If there exists a diagonal matrix such that the following LMI conditions hold for , then the unforced biochemical network (7) is positive and asymptotically stable.

Proof. Construct the following quadratic Lyapunov function candidate for the unforced biochemical network (7) where the diagonal matrix is to be determined.
Taking the derivative along the trajectory of (7), one can readily get
Then, it follows immediately from condition (9) that for all nonzero , and if and only if . Hence, the unforced biochemical network (7) is asymptotically stable.
Furthermore, since is a positive definite diagonal matrix, that is, , it can be deduced from the LMI condition (10) that matrix is a Metzler matrix for every local linear model of (7). Moreover, it follows from condition (8) that the membership function satisfies which together with the fact that is Metzler can guarantee the unforced biochemical network (7) to be positive. Hence, together with condition (13), it can be concluded that the equilibrium point of the unforced biochemical network (7) is positive and asymptotically stable. The proof is thus completed.

Remark 10. It is noted that the obtained sufficient condition for positive stability is in the form of linear matrix inequalities (LMIs), which could be efficiently solved by using the Matlab LMI control toolbox. Compared with the results developed in [2025], the regulation function is no longer needed to satisfy a sector condition; thus, the regulation function considered here is more general.

4. Bio-Circuits Design for Nonlinear Biochemical Networks

If the equilibrium point of the unforced nonlinear biochemical network (7) is not stable, bio-circuits design would become necessary for these complex biological systems to work properly, which would be useful for drug design and gene therapy. In this paper, the following smooth controller is employed to stabilize the biochemical network (7): which can be rewritten as where are the feedback gains of the bio-circuits to be determined.

The controlled biochemical network (7) can be described as follows:

Then, the following bio-circuits design results would guarantee the positivity and asymptotic stability of the nonlinear biochemical network (17).

Theorem 11. There exists a smooth control scheme such that the biochemical network (17) is positive and asymptotically stable, if there exist a diagonal matrix and matrices such that the following LMI conditions hold for .

Moreover, the feedback gains can be computed as

Proof. Construct the following quadratic Lyapunov function for the biochemical network (17) where the diagonal matrix is to be determined.
Taking the derivative along the trajectory of system (17), one has It is noted that which is equivalent to Therefore, it follows immediately from conditions (18) and (20) that for all the nonzero , and if and only if . Hence, the controlled biochemical network (17) is asymptotically stable.
On the other hand, since is a positive definite diagonal matrix, that is, , it follows from conditions (19) and (20) that the off-diagonal elements of matrix are nonnegative; that is, This implies that for any , is a Metzler matrix for every linear local model. Similar to the proof of Theorem 9, from condition (14), it can be concluded that the controlled biochemical network (17) is positive. Therefore, together with condition (25), one can conclude that the controlled biochemical network (17) is positive and asymptotically stable. The proof is, thus, completed.

It is worth pointing out that in many practical applications, drugs or chemical complexes delivered to the human body are often taken as control inputs to biochemical systems, and in this case, control signals have to be nonnegative. Therefore, it is necessary to consider the positive constraints of control inputs when designing bio-circuits. We have the following results.

Theorem 12. There exists a positive smooth control scheme , that is, , such that the biochemical network (17) is positive and asymptotically stable, if there exist a diagonal matrix and matrices such that the following LMI conditions hold for .

Moreover, the feedback gains can be computed as

Remark 13. It is noted that the positive feedback control law is only available in some special cases. When is a Metzler matrix and , it is impossible to design a positive feedback bio-circuit.

Remark 14. In many practical applications, the concentrations of some materials should be rigorously kept below or above a certain level, otherwise they may have side effects on blood or other chemical complexes. In other words the states of biochemical networks should be subject to some kinds of constraints. How to deal with such scenarios is one of our future research interests.

5. Design Example In Silico for the Proposed Method

Consider a two-compartment model [31], which describes the kinetics of a drug in the human body. The drug is injected into the blood where it exchanges linearly with the tissues; the drug is irreversibly removed with a nonlinear saturative characteristic from the blood and with a linear one from the tissue. The model can be expressed by the following nonlinear differential equations: where are the drug masses in blood and tissues, respectively; is the drug input is the measured drug output in the blood; , and are the constant rate parameters; and are the Michaelis-Menten parameters; and are the input and output parameters, respectively.

Let the premise variable , then the membership functions can be chosen as By using and , the biochemical network (29) can be expressed by the following T-S fuzzy model where and

Let , and ; Figure 1 shows that the state trajectories of biochemical network (29) are divergent under with the initial condition . For convenience, we use the same initial conditions in the following simulations.

5.1. Bio-Circuits Design without Constraints

In this case study, we will consider the bio-circuits design for the biochemical network (31) without constraints. Solving the LMI conditions (18)-(19) leads to the feasible solutions as follows: Then the feedback gains can be calculated as , which could guarantee the controlled biochemical network (31) to be positive and asymptotically stable. The dynamic response of the controlled biochemical network can be seen in Figure 2, and the evolution of the control input can be seen in Figure 3. It can be observed that the control input is negative in some stage.

5.2. Bio-Circuits Design with Positive Control Constraints

In most situations, the control input is the drug or the mixture of some biochemical complexes, which should be constrained to be positive. From Theorem 12, the elements of nonsymmetric matrixes should satisfy conditions (27). By using the constrained LMI algorithm, the following feasible solutions are obtained: Then, the feedback gains can be calculated as . The state response of the biochemical network under this constrained control can be seen in Figure 4, and the evolution of the positive control input can be seen in Figure 5, where the control signal remains to be positive all the time.

Remark 15. By using the fuzzy interpolation approach, bio-circuits can be easily implemented for the nonlinear biochemical network (29), which could guarantee it to be positive and asymptotically stable. Although the stability conditions for this biochemical network can also be derived by the method proposed in [29], the positivity of the states cannot be guaranteed. This would greatly reduce the significance of those results in real application such as drug delivery.

6. Discussion

In this paper, the fuzzy interpolation approach has been employed to approximate nonlinear biochemical networks for positive stability analysis and bio-circuits design. A few results on sufficient conditions for positivity and asymptotic stability of the network have been obtained in terms of a number of linear matrix inequalities. In addition, the positive constraint on control input is also considered for bio-circuits design. Finally, a real biochemical network model was provided to illustrate the effectiveness and validity of the obtained results.

Due to the transcription, translation, diffusion, and translocation processes of genes, time delays are inevitable in describing the dynamics of biochemical networks [7]. In addition, biochemical networks are often subject to intrinsic and extrinsic perturbations such as gene expression noises, mutation and disturbance from uncertain environment, and the fractal and chaotic features of systems [32]. Therefore, how to design constrained robust bio-circuits for such biochemical networks will be an interesting and challenging task.

Acknowledgment

This work was supported in part by the National Natural Science Foundation of China under Grants 61104045, 51107032, and 51277052, and in part by the Fundamental Research Funds for the Central Universities of China under Grant 2012B03514.