Home | About | Journals | Submit | Contact Us | Français |

**|**Cogn Neurodyn**|**v.3(3); 2009 September**|**PMC2727162

Formats

Article sections

- Abstract
- Introduction
- Model and analysis
- Stochastic stability condition of uncertain genetic networks with time-varying delays
- Illustrative examples
- Conclusions
- References

Authors

Related links

Cogn Neurodyn. 2009 September; 3(3): 271–280.

Published online 2009 July 30. doi: 10.1007/s11571-009-9077-1

PMCID: PMC2727162

Zhengxia Wang, Email: moc.361@gnawaixz.

Received 2008 June 14; Revised 2009 January 31; Accepted 2009 January 31.

Copyright © Springer Science+Business Media B.V. 2009

This article has been cited by other articles in PMC.

Gene regulation is an intrinsically noisy process, which is subject to intracellular and extracellular noise perturbations and environment fluctuations. In this paper, we consider the robust stability analysis problem of genetic regulatory networks with time-varying delays and stochastic perturbation. Different from other papers, the genetic regulate system considers not only stochastic perturbation but also parameter disturbances, it is in close proximity to the real gene regulation process than determinate model. Based on the Lyapunov functional theory, sufficient conditions are given to ensure the stability of the genetic regulatory networks. All the stability conditions are given in terms of LMIs which are easy to be verified. Illustrative examples are presented to show the effectiveness of the obtained results.

Development of genome sequencing and gene recognition has been accumulating a wealth of experimental data. This creates a new challenge to biology of understanding how genes and proteins work collectively, i.e. the analysis of these data. This challenge leads to a significant increase of computer applications for modeling and data interpretation methods. The aim is to develop computer simulations that mimic biological phenomena, data or patterns, such as complex biochemical reactions and genetic networks in cellular media (Turner et al. 2004). Real genetic systems are composed of a large number of reactions and reacting species. There are too many items to include them all in models. It is rather difficult to find an effective method to construct a complete model and analyze such a complex model. So we only consider the concentrations of mRNAs and proteins. Recently, there have been many efforts for modeling genetic regulatory networks using different classes of mathematical models (Hasty et al. 2001; De Jong 2002). Basically, there are two types of genetic network models. i.e. the Boolean models and the differential equation model (Smolen et al. 2000; Kobayashi et al. 2002; Bolouri and Davidson 2002; Wang et al. 2004; Benuskova and Kasabov 2008). Boolean network model interprets gene interactions as connections between genes. The state of gene expression is simplified as being either completely ON or completely OFF. The Boolean expression state converges to a terminal state via a series of state transitions that is determined by the designed Boolean rules. In the differential equation models, the variables describe the concentrations of gene products, such as mRNAs and proteins, as continuous values of the gene regulation systems.

Recently, the genetic regulatory networks in the form of differential equations have been well studied. Becskei and Serrano (2000) designed and constructed simple gene circuits consisting of a regulator and transcriptional repressor modules in *Escherichia coli* and they showed the gain of stability produced by negative feedback. In Chen and Aihara (2002), a model for genetic regulatory networks with time delays was proposed and nonlinear properties of the model in terms of local stability and bifurcation was analyzed. In Li et al. (2006a), a nonlinear model for genetic regulatory networks with SUM regulatory functions was presented. Genetic networks with delays and stochastic perturbations were studied and sufficient conditions of stability were derived in terms of linear matrix inequalities (LMI). Li et al. (2006b) provided a theoretical method for analyzing the synchronization of coupled nonidentical genetic oscillators. Sufficient conditions for the synchronization as well as the estimation of the bound of the synchronization error were also obtained. Authors Ren and Cao (2008) studied the robust stability of the genetic regulatory networks with time-delays, and present some sufficient conditions using Lyapunov functional theory and LMI technique. In Cao and Ren (2008), discrete-time versions of the continuous-time genetic regulatory networks with SUM regulatory functions are formulated and studied, and obtained sufficient conditions for exponential stability of the discrete-time genetic regulatory networks with delays.

In fact, for most genetic regulatory system, there are two types of reactions (De Jong 2002): fast reaction and slow reaction. Fast reaction, such as dimerization, binding reactions and other medical modification reaction, we can assume this reaction is immediately and time delay is reduced to zero. While transcription and translation involve a number of multi-stage reactions, there is a time lag in the peaks between mRNA molecules and proteins of gene. On the other side, mRNA and proteins may be synthesized at different locations (i.e. nucleus and cytoplasm, respectively), thus transportation or diffusion of mRNA and proteins between these two locations results in sizeable delays. That is, time delays exist in genetic regulatory networks, and possible effects of time delays have attracted some attentions (Chen and Aihara 2002; Li et al. 2006a; Ren and Cao 2008; He and Cao 2008).

Stochasticity is ubiquitous in biology. Noise in the form of random fluctuations arises in genetic regulatory network in one of two ways. Intrinsic noise is inherent in the biochemical reactions. Its magnitude is proportional to the inverse of the system size, and its origin is often thermal. On the other hand, external noise originates environment fluctuation (Hasty et al. 2000).

In the applications and designs of genetic networks, there are often some unavoidable uncertainties such as model errors, external perturbations, and parameter fluctuations, which can cause the networks to be unstable (Ren and Cao 2008). There are some papers have studied stability of neural networks with stochastic perturbations or parameter uncertainties (Huang and Feng 2007; Liao et al. 2001; Wang et al. 2006, 2007; Zhang et al. 2007), that give us some suggests for studying genetic regulatory networks. In this paper, we aim to analyze the stability of genetic networks in the forms of differential equations. We consider the delayed genetic regulatory networks not only with stochastic perturbations but also with parameter uncertainties. To our best knowledge, there are few paper to investigate it. By using Lyapunov functional theory and LMI technique, Novel criteria are derived to guarantee the asymptotic and robust stability of such genetic networks.

The rest of this paper is organized as follows. In section “Model and analysis”, problem formulation and preliminaries are given. In section “Stochastic stability condition of uncertain genetic networks with time-varying delays”, several sufficient criteria are derived for checking globally robust stability of the genetic regulatory networks with stochastic perturbations and time-varying delays. In section “Illustrative examples”, two examples are given to show the effectiveness of the proposed results. Finally, conclusions are given in section “Conclusions”.

**Notation** For convenience, some notations are introduced. For a real square matrix *X*, the notation *X* > 0(*X* < 0) means that *X* is symmetric and positive definite (negative definite). *I* is the identity matrix with appropriate dimension. The superscript ‘‘*T*’’ represents the transpose. For *τ* > 0, ([−*τ*, 0]; *R*^{n}) denotes the family of continuous functions from [−*τ*, 0] to *R*^{n} with the norm |||| = sup_{–τ ≤ ≤ 0}|()|. Let (Ω, *F*, {*F*_{t}}_{t 0}, *P*) be a complete probability space with a filtration {*F*_{t}}_{t ≥ 0} satisfying the usual conditions (i.e. it is right continuous and *F*_{0} contains all *P*-pull sets); the family of all *F*_{0}-measurable ([−*τ*, 0]; *R*^{n})-valued random variables *ξ* = {*ξ*(*θ*): −*τ* ≤ *θ* ≤ 0} such that sup_{–τ ≤ θ ≤ 0} *E*|*ξ*(*θ*)|^{p} < ∞ where *E* stands for the mathematical expectation operator with respect to the given probability measure *P*; ^{2,1}(*R*^{n} × *R*^{+}; *R*^{+}) the family of all nonnegative functions *V*(*x*, *t*) on *R*^{n} × *R*^{+} which are continuously twice differentiable in *x* and differentiable in *t*.

Authors considered a genetic regulatory network model (Li et al. 2006a):

1

where are the concentrations of mRNA and protein of the *i*th node. *a*_{i} and *c*_{i} are the rates of degradation of mRNA and protein, respectively; *d*_{i} is the translation rate, and *b*_{i} is the regulatory function of the *i*th gene, which is generally a nonlinear function of the variables (*p*_{1}(*t*), *p*_{2}(*t*), …, *p*_{n}(*t*)), and nonlinear function is monotonic with each variable.

In this paper, we considered the genetic regulatory networks with time delay described by the following differential equation:

2

where *τ*_{1}(*t*), *τ*_{2}(*t*) are transcriptional delays, for any single gene in the network, there are one output *p*_{i}(*t *− *τ*_{1}(*t*)) to other genes and multiple inputs *p*_{j}(*t* − *τ*_{1}(*t*)) (*j* = 1, 2, …, *n*) from other genes. The structure and regulation mechanism of the genetic work can be seen in Chen and Aihara (2002). As a monotonic increasing or decreasing regulatory function, *B*_{i} is usually the Michaelis–Menten or Hill form. In this paper we take *B*_{i}(*p*_{1}(*t*), *p*_{2}(*t*), …, *p*_{n}(*t*)) = ∑_{j}*B*_{ij}(*p*_{j}(*t*)), which is called SUM logic (Li et al. 2006a). That is each transcription factor acts additively to regulate the *i*th gene. The function *B*_{ij}(*p*_{j}(*t*)) is a monotonic function of the Hill form.

3

where *H* is the Hill coefficient, *β* is a positive constant, and *α*_{ij} is the dimensionless transcriptional rate of transcriptional factor *j* to gene *i*, which is a bounded constant. Note that

Hence, we can rewrite (2) as:

4

where *g*(*x*) = (*x*/*β*)^{H}/[1 + (*x*/*β*)^{H}] is a monotonically increasing function, *I*_{i} is defined as basal rate and *L*_{i} is the set of repressors of gene *i*. is defined as:

5

In compact matrix form, (4) can be rewritten as

6

where *m*(*t*) = [*m*_{1}(*t*), *m*_{2}(*t*), …, *m*_{n}(*t*)]^{T}, *p*(*t*) = [*p*_{1}(*t*), *p*_{2}(*t*), …, *p*_{n}(*t*)]^{T}, *I* = [*I*_{1}, *I*_{2}, …, *I*_{n}]^{T}, *A* = diag{*a*_{1}, *a*_{2, }…, *a*_{n}}, *C* = diag{*c*_{1}, *c*_{2, }…, *c*_{n}}, *D* = diag{*d*_{1}, *d*_{2}, …, *d*_{n}}. Let (*m**, *p**)^{T} be an equilibrium point of (6), which satisfies the following relationship:

7

For convenience, we will shift an intended equilibrium point (*m**, *p**) of the system (6) to the origin. Using *x*(*t*) = *m*(*t*) − *m**, *y*(*t*) = *p*(*t*) − *p**, we have

8

where *x*(*t*) = [*x*_{1}(*t*), *x*_{2}(*t*), …, *x*_{n}(*t*)]^{T}, *y*(*t*) = [*y*_{1}(*t*), *y*_{2}(*t*), …, *y*_{n}(*t*)]^{T}, *f*(*y*(*t*)) = [*f*_{1}(*y*_{1}(*t*)), *f*_{2}(*y*_{2}(*t*)), …, *f*_{n}(*y*_{n}(*t*))]^{T}, *f*_{i}(*y*_{i}(*t*)) = *g*_{i}(*y*_{i}(*t*) + *p*_{i}*) − *g*_{i}(*p*_{i}*), since *g*_{i}(·) is a monotonically increasing function with saturation, and if there exists matrix *G* = diag{*g*_{1}, *g*_{2}, …}, then it satisfies that

for all with *x* ≠ *y*. From the relationship of *f*(·) and *g*(·), *f*(·) satisfied that

9

which is equivalent to the following one:

10

In the following, we consider stability of delayed genetic regulatory networks with parameter uncertainties and stochastic perturbations:

11

where Δ*A*, Δ*B*, Δ*C*, Δ*D*, Δ*H*_{0} and Δ*H*_{1} are parameter uncertainties, *τ*_{1}(*t*), *τ*_{2}(*t*) are unknown time-varying delays satisfying 0 ≤ *τ*_{1}(*t*) ≤ *h*_{1}, 0 ≤ *τ*_{2}(*t*) ≤ *h*_{2} and where *h*_{1}, *h*_{2}, *d*_{1} and *d*_{2} are known constants. is a *m*-dimensional Brownian motion defined on a complete probability space (Ω, *F*, {*F*_{t}}_{t ≥ 0}, *P*), and Δ*A*, Δ*B*, Δ*C*, Δ*D*, Δ*H*_{0} and Δ*H*_{1} are defined as:

12

where *M*_{1}, *M*_{2}, *M*_{3}, *M*_{4}, *M*_{5}, *M*_{6}, *E*_{1}, *E*_{2}, *E*_{3}, *E*_{4}, *E*_{5} and *E*_{6} are known constant real matrices with appropriate dimensions. *F*(*t*) is unknown time-varying matrix satisfying

13

**Remark 1** It should be noted that, if let *H*_{0} = *H*_{1} = 0, the model (11) becomes the same form as Ren and Cao (2008). If we let Δ*A* = Δ*B* = Δ*C* = Δ*D* = 0, the model (11) is equivalent to the one investigated in Li et al. (2006a).

**Remark 2** The parameter uncertainty structure described in (12) and (13) has been widely exploited in robust control and robust filtering of uncertain systems. Many practical systems possess parameter uncertainties which can be either exactly modeled or over bounded by (13). The stochastic term, [(*H*_{0} + Δ*H*_{0})*y*(*t*) + (*H*_{1} + Δ*H*_{1})*y*(*t* − *τ*_{1}(*t*)]d*ω*(*t*), can be viewed as stochastic perturbations related to the gene current states and delayed states.

**Remark 3** Time delay is an important factor in considering the stability of real genetic system. It is known that the time delay is different at different reaction stages and at different time due to complexity and sensitiveness of biochemical reaction (De Jong 2002), so constant delays cannot reflect the real regulate process, accordingly, we use time-varying delays in this paper.

Next, we give the definition of global robust asymptotic stability for the uncertain stochastic genetic regulatory networks:

**Definitions 1** For the genetic regulatory network (11) and the trivial solution (equilibrium point) is robustly, globally, asymptotically stable in the mean square if, for all admissible uncertainties satisfying (12), the following holds:

14

We recall the following useful lemmas:

**Lemma 1** (Boyd et al. 1994). The following LMI

where *Q*(*x*) = *Q*^{T}(*x*), *R*(*x*) = *R*^{T}(*x*), and *S*(*x*) depend affinely on *x*, is equivalent to

**Lemma 2** (Gu 2000). For any constant matrix , *M* = *M*^{T} > 0, scalar *γ* > 0, vector function *ω*:[0, *γ*] → *R*^{n} such that the integrations are well defined, the following inequality holds:

**Lemma 3** (Liao et al. 2002) Given any real matrices Σ_{1}, Σ_{2}, Σ_{3} of appropriate dimensions and a scalar > 0 such that Then, the following inequality holds:

We first consider stochastic genetic regulatory networks without parameter uncertainties:

15

We have the following result:

**Theorem 1** System (15) is globally asymptotically stable in the mean square, if there exist symmetric positive definite matrices *P*_{1}, *P*_{2}, *Q*_{1}, *Q*_{2}, *S*_{1}, *S*_{2} and positive scalar *λ*, such that the following LMIs holds:

16

where

**Proof** Using schur complement, Ω_{2} < 0 implies that

Constructing a positive definite Lyapunov–Krasovskii functional as follows:

By *Itô’s* differential formula (Oksendal 2003), the stochastic derivative of *V*(*x*(*t*), *y*(*t*), *t*) along the trajectory of system (15)

From Lemma 3, the following inequalities hold:

17

18

From (9), we obtain the follow inequality easily,

19

Noticing that, for a scalar *λ* > 0, we have

20

Substituting (17), (18) and (20) into d*V*(*x*(*t*), *y*(*t*), *t*) and by use of Lemma 2, we have

21

where

Since the expectation of is equal to zero, when Ω_{1} < 0 and Ω_{2} < 0. Taking the mathematical expectation of both sides of (21), we have

22

where *E* is the mathematical expectation operator. It indicates that the genetic regulatory network (15) is asymptotically stable in mean square. This completes the proof.

**Remark 4** We introduce a new Lyapunov–Krasovskii functional. The Lyapunov–krasovskii functional not only dealt with the time-varying delays but also consider the upper bounds of time delays. That is, both time-varying delay and the upper bound of time delays have been brought into the final robust stability condition, and the effect factor of time delays is considered in genetic regulatory networks, the LMIs condition can reflect more characters of genetic regulatory networks.

**Remark 5** In spite of considerable variations and random perturbations of biochemical parameters, the genetic regulatory networks can keep homeostasis in metabolism and developmental programs of living cells (Yuh et al. 1998). i.e. despite the stochastic function of regulatory within cells, most cellular events are ordered and precisely regulated. That is, genetic networks system can reach stability by auto regulation. In fact, genetic regulatory network is an essentially continuous and complicated dynamical system. So the research of stability of genetic regulatory networks is necessary. The LMIs conditions provide one sufficient criteria of estimating the stability of genetic regulatory networks.

**Theorem 2** The genetic regulatory network (11) is asymptotically robustly stable if there exist positive matrices *P*_{i}, *Q*_{i}, *S*_{i} (*i* = 1, 2) and scalars *λ* > 0 and *l*_{i} > 0 (*i* = 1, 2, …, 6) such that the following LMIs hold:

23

where

and * denotes the symmetric term in a symmetric matrix.

**Proof** Taking the same Lyapunov–Krasovskii functional as that in the proof of Theorem 2, and replacing *A*, *B*, *C*, *D*, *H*_{0} and *H*_{1} by *A* + *M*_{1}*F*(*t*)*E*_{1}, *B* + *M*_{2}*F*(*t*)*E*_{2}, *C* + *M*_{3}*F*(*t*)*E*_{3}, *D* + *M*_{4}*F*(*t*)*E*_{4}, *H*_{0} + *M*_{5}*F*(*t*)*E*_{5}, *H*_{1} + *M*_{6}*F*(*t*)*E*_{6}, respectively.

Since we have

26

For all *x*(*t*), *y*(t) except for *x*(*t*) = *y*(*t*) = 0. Therefore, the genetic regulatory network (11) is asymptotically stable in mean square. This completes the proof.

**Remark 6** Noise always exist and need processing in process of data measurement. The gene expression data obtained from experiments is especially noisy. Robustness is a accuracy measurement of extracting the weight matrix of genetic regulatory network with the existence of noise

**Remark 7** In Li et al. (2006a), the authors studied stochastic stability of genetic networks with time-varying delays, however, the parameter uncertainty was not considered in the models, and The stochastic disturbance term in this paper we give is the form of Therefore, our results and those established in [11] are complementary each other.

**Remark 8** Gene regulation process is an inherent noisy process, the determinate model cannot commendably reflect the real gene regulation process, the stochastic model considering noise disturbances and parameter uncertainties in this paper is in close proximity to the real gene regulation process than determinate model. Ren and Cao studied the uncertain genetic networks with time-varying delays, and several LMI-based conditions were proposed to guarantee the stability of the equilibrium point of genetic networks. However, the stochastic term was not taken into account in the models. Therefore, our developed results in this paper are more general than those reported in Ren and Cao (2008).

**Remark 9** A cycle of most eukaryotes is composed of four stages (Chen and Wang 2004): G1(Gap) phase in which size of the cell is increased by constantly producing RNA and synthesizing protein. S phase in which the cell continuous to produce new proteins and grows in size. And M(mitosis) phase in which chromosomes segregate and cell division takes place. While regulation of gene expression at different stages of protein synthesis (De Jong 2002). Transcription and translation are important process, from the description of introduction, we known that there exist time-varying delays in process of transcription and translation, that is, the time-varying delays can also affect the cell cycle, which need much more further research.

In this section, we present two examples to show the effectiveness and correctness of our theoretical results.

** Example 1** In Elowitz and Leibler (2000), the dynamics of repressilator has been theoretically predicted and experimentally investigated in

where *i* = lacl, tetR, cl; *j* = cl, lacl, tetR, *m*_{i} and *p*_{i} are the concentrations of the three mRNA and repressor-proteins, and *β* > 0 denotes the ratio of the protein decay rate to mRNA decay rate.

We consider a three-node genetic regulatory networks (15) with stochastic disturbance and transcriptional time delay, where

and *f*(*x*) = *x*^{2}/(1 + *x*^{2}), i.e. the Hill coefficient is 2. It is easy to know that the maximal value of the derivative of *f*(*x*) is less than *g*_{i} = 0.65, Assume time delay *τ*_{1}(*t*) = 1 + 0.5 sin (*t*), *τ*_{1}(*t*) = 0.5 + 0.5 cos (*t*) and *τ*_{1} = 1.25, *τ*_{2} = 0.75, *d*_{1} = *d*_{2} = 0.25. According to Theorem 1, if the LMI (16) hold, then the genetic network is globally asymptotic stable in mean square. The simulation of trajectories and phase graph of *m*(*t*) and *p*(*t*) are show in Fig. 1. By using the MATLAB LMI Toolbox (Gahinet et al. 1995), we solve the LMI (16) and obtain

λ = 63.9422.

Therefore, the genetic regulatory network (15) with stochastic disturbance and time-varying delays is asymptotically stable.

** Example 2** Now, let us consider a five-node delayed stochastic genetic regulatory network (11) with norm bounded uncertainties. The network coefficients are given as follows:

The simulation of trajectories and phase graph of *m*(*t*) and *p*(*t*) with are show in Fig. 2. Again, by solving the LMI (23) for the same parameters, we have

*P*_{1}= diag(100.0412, 81.0543, 59.9681, 77.5957, 81.0543),*P*_{2}= diag(75.4597, 99.5427, 75.0213, 85.6906, 99.5427),*Q*_{1}= diag(294.2091, 291.9178, 294.2575, 293.3227, 291.9178),*Q*_{2}= diag(208.6455, 208.6287, 208.8624, 208.6488, 208.6287),*S*_{1}= diag(120.3644, 120.3778, 120.3997, 120.3833, 120.3778),*S*_{2}= diag(111.6799, 110.7737, 111.6990, 111.3293, 110.7737),- λ = 157.0317,
*l*_{1}= 150.5257,*l*_{2}= 150.6112,*l*_{3}= 150.4322, *l*_{4}= 150.5630,*l*_{5}= 151.2316,*l*_{6}= 151.2004

which indicates from Theorem 2 that the delayed uncertain stochastic genetic regulatory network (11) is robustly, globally, asymptotically stable in the mean square.

In this paper, we have dealt with the problem of global asymptotic and robust stability analysis for stochastic delayed genetic regulatory networks, which involve stochastic perturbations, parameter uncertainties and time-varying delays. There are exist time-varying delays in process of transcription and translation of gene expression, while process of transcription and translation occupy a majority of cell cycle. So the time-varying delays can also affect the cell cycle. Effect on gene regulate process of cell cycle may result in other dynamical behaviors, such as switch, oscillation and bifurcation, which need much more further research.

The authors would like to thank the reviewers and the editor for their valuable suggestions and comments which have led to a much improved paper. This work is partially supported by Natural Science Foundation Project of Chongqing CSTC (Grant No. 2007BB2395, 2008BB2366). Teaching and Research Program of Chongqing Education Committee (KJ090803). Fund for Young on Campus (No. 2005(35)).

Becskei A, Serrano L (2000) Engineering stability in gene networks by autoregulation. Nature 405:590–593. doi:10.1038/35014651 [PubMed]

Benuskova L, Kasabov N (2008) Modeling brain dynamics using computational neurogenetic approach. Cogn Neurodyn 2:319–334. doi:10.1007/s11571-008-9061-1 [PMC free article] [PubMed]

Bolouri H, Davidson EH (2002) Modelling transcriptional regulatory networks. BioEssay 24:1118–1129. doi:10.1002/bies.10189 [PubMed]

Boyd S, El Ghaoui L, Feron E, Balakrishnan V (1994) Linear matrix inequalities in system and control theory. SIAM, Philadelphia

Cao JD, Ren FL (2008) Exponential stability of discrete-time genetic regulatory networks with delays. IEEE Trans Neural Netw 19(3):520–523. doi:10.1109/TNN.2007.911748 [PubMed]

Chen LN, Aihara K (2002) Stability of genetic regulatory networks with time delay. IEEE Trans Circuits Syst I Fundam Theory Appl 49(5):602–608. doi:10.1109/TCSI.2002.1001949

Chen LN, Wang RQ et al (2004) Dynamics of gene regulatory networks with cell division cycle. Phys Rev E Stat Nonlin Soft Matter Phys 70:011909. doi:10.1103/PhysRevE.70.011909 [PubMed]

De Jong H (2002) Modeling and simulation of genetic regulatory systems: a literature review. J Comput Biol 9:67–103. doi:10.1089/10665270252833208 [PubMed]

Elowitz MB, Leibler S (2000) A synthetic oscillatory network of transcriptional regulators. Nature 403:335–338. doi:10.1038/35002125 [PubMed]

Gahinet P, Nemirovski A, Laub A, Chialali M (1995) LMI control toolbox user’s guide. The Mathworks, Natick- Gu K (2000) An integral inequality in the stability problem of time-delay systems. In: Proceedings of 39th IEEE conference on decision and control. Sydney, Australia, pp 2805–2810

Hasty J, Pradines J, Dolnik M, Collins J (2000) Noised-based switches and amplifiers for gene expressions. Proc Natl Acad Sci USA 97(5):2075–2080. doi:10.1073/pnas.040411297 [PubMed]

Hasty J, McMillen D, Isaacs F, Collins JJ (2001) Computational studies of gene regulatory networks: in numero molecular biology. Nat Rev Genet 2:268–279. doi:10.1038/35066056 [PubMed]

He WL, Cao JD (2008) Robust stability of genetic regulatory networks with distributed delay. Cogn Neurodyn 2:355–361. doi:10.1007/s11571-008-9062-0 [PMC free article] [PubMed]

Huang H, Feng G (2007) Delay-dependent stability for uncertain stochastic neural networks with time-varying delay. Physica A 381:93–103. doi:10.1016/j.physa.2007.04.020

Kobayashi T, Chen L, Aihara K (2002) Modeling genetic switches with positive feedback loops. J Theor Biol 221:379–399. doi:10.1006/jtbi.2003.3190 [PubMed]

Li CG, Chen LN, Aihara K (2006a) Stability of genetic networks with SUM regulatory logic: Lur’e system and LMI approach. IEEE Trans Circuits Syst I Regul Pap 53(11):2451–2458. doi:10.1109/TCSI.2006.883882

Li CG, Chen LN, Aihara K (2006b) Synchronization of coupled nonidentical genetic oscillators. Phys Biol 3:37–44. doi:10.1088/1478-3975/3/1/004 [PubMed]

Liao XF, Wong K-W, Wu Z, Chen G (2001) Novel robust stability criteria for interval-delayed Hopfield neural networks. IEEE Trans Circuits Syst I 48(11):1355–1359. doi:10.1109/81.964428

Liao XF, Chen GR, Sanche ZEN (2002) LMI-based approach for asymptotical stability analysis of delayed neural network. IEEE Tans Circuits Syst I 49(7):1033–1039. doi:10.1109/TCSI.2002.800842

Oksendal B (2003) Stochastic differential equations, 4th edn. Springer-Verlag, New York

Ren FL, Cao JD (2008) Asymptotic and robust stability of genetic regulatory networks with time-varying delays. Neurocomputing 71:834–842. doi:10.1016/j.neucom.2007.03.011

Smolen P, Baxter DA, Byrne JH (2000) Mathematical modeling of gene networks. Neuron 26:567–580. doi:10.1016/S0896-6273(00)81194-0 [PubMed]

Turner TE, Schnell S, Burrage K (2004) Stochastic approaches for modelling in vivo reactions. Comput Biol Chem 28:165–178. doi:10.1016/j.compbiolchem.2004.05.001 [PubMed]

Wang R, Zhou T, Jing Z, Chen L (2004) Modelling periodic oscillation of biological systems with multiple time scale networks. Syst Biol 1:71–84. doi:10.1049/sb:20045007 [PubMed]

Wang ZD, Shu H, Fang J, Liu X (2006) Robust stability for stochastic Hopfield neural networks with time delays. Nonlinear Anal Real World Appl 7(5):1119–1128. doi:10.1016/j.nonrwa.2005.10.004

Wang ZD, Lauria S, Fang J, Liu X (2007) Exponential stability of uncertain stochastic neural networks with mixed time-delays. Chaos, Solitons Fractals 32:62–72. doi:10.1016/j.chaos.2005.10.061

Yuh CH, Bolouri H, Davidson EH (1998) Genomic cis-regulatory logic: experimental and computational analysis of a sea urchin gene. Science 279:1896–1902. doi:10.1126/science.279.5358.1896 [PubMed]

Zhang JH, Shi P, Qiu JQ (2007) Novel roubust stability criteria for uncertain stochastic Hopfield neural networks with time-varying delays. Nonlinear Anal Real World Appl 8:1349–1357. doi:10.1016/j.nonrwa.2006.06.010

Articles from Cognitive Neurodynamics are provided here courtesy of **Springer Science+Business Media B.V.**

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |