The Scientific World Journal

Volume 2014 (2014), Article ID 476912, 8 pages

http://dx.doi.org/10.1155/2014/476912

Research Article

Observer Design for a Core Circadian Rhythm Network

Yuhuan Zhang

School of Mathematics and Information Sciences, Henan University, Kaifeng 475004, China

Received 8 June 2014; Accepted 15 June 2014; Published 8 July 2014

Academic Editor: Guanghui Wen

Copyright © 2014 Yuhuan Zhang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The paper investigates the observer design for a core circadian rhythm network in Drosophila and Neurospora. Based on the constructed highly nonlinear differential equation model and the recently proposed graphical approach, we design a rather simple observer for the circadian rhythm oscillator, which can well track the state of the original system for various input signals. Numerical simulations show the effectiveness of the designed observer. Potential applications of the related investigations include the real-world control and experimental design of the related biological networks.

1. Introduction

Controllability, observability, and stability are typical problems of dynamical systems [1, 2]. Suppose we have a dynamical system with the following form:

where is the state vector of the system and is the input vector. System (1) is said to be observable, if we can find a function

from which one can determine all the state variables of system (1). Here, depends on , a set of the system’s state and the external input [1–6].

For with polynomial or rational expressions, existing results have reported that system (1) is observable if the Jacobian matrix has full rank [3, 4], where is the Lie derivative of the output function .

Recently, Liu et al. [4] proposed a graphical approach to reduce the observability problem to a property of the inference diagram of a system. The inference diagram is based on the dynamical equation of the system. If appears in ’s differential equation, then there is a link from to in the inference diagram. The inference diagram can be decomposed into some strongly connected component (SCC). Those SCCs without incoming edges are called root SCCs. Liu et al. [4] reported that the number of the root SCCs provides the lower bound of the monitored state variables in . Furthermore, for many nonlinear systems, they declared that the number of the root SCCs provides not only necessary but also sufficient numbers of state variables to realize observability. That is, if the inference diagram of a system has root SCCs, then one may only need to select state variables from different root SCCs, and the system will be observable through monitoring these state variables.

Biological systems are typical nonlinear systems [7–23]. Observer design for biological systems has important real-world implications. For example, through monitoring a few state variables of a complex biological system, if we can infer the state of the whole system, then lots of resources can be saved [4]. Circadian rhythms are typical biological phenomenon, which display endogenous, entrainable oscillations with a period that lasts approximately 24 hours. Circadian rhythms are widely in existence in various plants and animals [12], which are controlled by biomolecular networks. Circadian rhythms have been extensively investigated during the last decades [18–23]. For example, in 1995, Goldbeter established a mathematical model for the circadian rhythms in the Drosophila [11]. In 1999, Leloup et al. established a model for the circadian rhythms in the Drosophila and Neurospora [13]. In 2002, Gonze and coauthors [7] investigated the deterministic and stochastic dynamics in a core circadian rhythm network. They found that the core network can display roughly the same circadian oscillations under both deterministic and stochastic descriptions.

In this paper, based on the graphical approach introduced by Liu et al. [4], for the core circadian rhythm network [7] in Drosophila and Neurospora, we design some simple observers for the network. Based on the Lyapunov stability theory, we theoretically verify the correctness of the designed observers. Finally, numerical simulations show the effectiveness of the designed observers for various input signals. The rest of the paper is organized as follows. In Section 2, we briefly introduce the mathematical model for the core circadian rhythm network. Observers will be designed in Section 3. We perform numerical simulations in Section 4. Discussions and some concluding remarks will be in Section 5.

2. The Core Circadian Rhythm Model

The core circadian rhythm network is reported by Gonze et al. [7]; the detailed biochemical processes are shown in Figure 1(a). Figure 1(a) represents a prototype for the molecular mechanism of circadian oscillations based on negative autoregulation of gene expression. Real-world circuits corresponding to Figure 1(a) include the per mRNA and PER protein in Drosophila [11, 12] and frq mRNA and FRQ protein in Neurospora [13].

(a)

(b)

Figure 1: (a) A core model for circadian rhythms [7]. Parameters marked on the directed edges denote reaction rates. denotes the concentration of the mRNA; are the concentrations of four forms of the protein expressed by the gene; are cytosolic proteins; is nuclear protein. The circadian gene transcripts mRNA and then translates protein . Proteins can be reversibly phosphorylated successively. can be transported into the nucleus; the nucleus protein can negatively regulate the expression of the circadian gene. The mRNA and protein suffer from degradations. (b) The inference diagram [4] for the circadian rhythm model (3). If appears in ’s differential equation, then there will be a directed edge from to .

The core model involves gene transcription and transport of mRNA into the cytosol where it is translated into protein and degraded. Protein can be reversibly phosphorylated from the form into the forms and , successively. The phosphorylated protein is degraded or transported into the nucleus, and the nucleus protein can negatively regulate the expression of its gene. Based on the work from Gonze et al. [7] in 2002, the modified mathematical model for the core circadian model can be established as follows. Consider

where are state variables, which correspond to species concentrations of the mRNA, the four forms of proteins. represents the external input, which can be seen as the effect of the environment on the system. are Michaelis constants, first-order reaction rate constants, and maximum rates of protein degradation, transcription, and phosphorylation. It is noted that if we set , then system (3) degenerates into the model investigated in [7]. A typical set of parameter values for system (3) are shown in Table 1. Under the parameter values as shown in Table 1 and for , dynamical system (3) can display circadian rhythms with a period close to 24 hours.

Table 1: Parameter values.

Parameter
Value

0.5 nMh⁻¹

4

0.2 nM

6.0 nMh⁻¹

3.0 nMh⁻¹

6.0 nMh⁻¹

3.0 nMh⁻¹

1.5 nMh⁻¹

2.0 h⁻¹

2 nM

0.3 nMh⁻¹

2.0 h⁻¹

1.5 nM

2.0 nM

1.5 nM

2.0 nM

0.1 nM

1.0 h⁻¹

3. Observer Design for the Circadian Rhythm Model

For simplicity, in the following, we rewrite system (3) as the following form:

where denotes the nonlinear term. denotes the monitored output. is a constant matrix. denotes the linear term, with

For system (4), similar to the works in [4–6], our objective is to design the following observer, which can track the states of system (4):

where is a gain matrix, which is to be determined. The estimation error dynamics are then given by

where .

From Liu et al. [4], the observability of a dynamical system can be revealed by its inference diagram. The inference diagram of system (3) is shown in Figure 1(b), where the five nodes are strongly connected and consist of the single root SCC. From the conclusion in [4], system (3) is observable, and one should only monitor any one of the five nodes in Figure 1(b). In the following, for simplicity, we assume that only is monitored. The designed observer of system (3) is described as

where are the nonzero elements in matrix and of (6). For simplicity, we denote

Then,

The corresponding error dynamics are described as

Before we analyze the stability of system (11), we note that since represent species concentrations, they must be nonnegative and bounded [24]. Furthermore, since is continuous and differentiable, by the mean value theorem, there must exist between and , satisfying

For system (11), we construct the following Lyapunov function:

Based on (12), the derivative of (13) along system (11) is

Here,

which is a symmetrical matrix. are values between and . For appropriate gain , if , . System (11) will be globally asymptotically stable. In other words, system (3) can be observed through the observer (8).

Remark 1. If any one of the other variables is used to track the state of the original system (3), one should only slightly revise the observer (8). If more than one variable is measured to track the original system, the observer can be similarly designed. For example, if in (4) is , the control gain matrix is chosen as

Then, the observer is designed as

For appropriate parameters , one can easily prove that the original system is also observable from the observer (17). Obviously, the corresponding observer (17) is more complex than the observer (11).

Remark 2. There are many methods to prove the stability of a dynamical system. One can easily prove that the nonlinear terms on the right-hand side of system (11) are Lipschitz. For Lipschitz nonlinear systems, Rajamani [5] proposed a general theorem for the observer design. However, due to the complexity of the biological model, the theorem obtained in [5] fails to work for the error system (11). Therefore, we have used the mean value theorem to simplify the error system (11) and obtained a sufficient condition for the observer design. The matrix relies on the bound of system (3).

4. Numerical Simulations

Hereinafter, we numerically verify the effectiveness of the designed observers. Firstly, we assume ; the output function . For , all the state variables in the original system (3) can be tracked by the observer (8). Under randomly initial values, Figure 2 shows the state trajectories of systems (3) and (8), as well as the error dynamics of system (11). From Figure 2, one can see that the state variables of (3) oscillate with a period close to 24 hours. The observer system (8) can track the states of the original system. The error between the observer system and the original system quickly approximates to zero.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 2: (a)–(e) State trajectories of systems (3) and (8). (f) The error dynamics of system (11). Here, . Initial values are randomly chosen.

For different input signals and under appropriate gain , system (8) can always track the states of system (3). For example, when , Figure 3 shows the state trajectories of systems (3) and (8) as well as the error dynamics of system (11). From Figure 3, we can see that the states of system (3) can be observed by the observer (8). The error system (11) converges to zero quickly. When is a step signal, the observer (11) can also well monitor the states of the original (3). Figure 4 shows the case for , and the step input signal

(a)

(b)

Figure 3: (a) State trajectories of systems (3) and (8) with . (b) The error dynamics of system (11). Here, .

(a)

(b)

Figure 4: (a) State trajectories of systems (3) and (8) with step input signal . (b) The error dynamics of system (11). Here, for and for .

For the cases discussed in Remark 1, when we choose , and the following input signal,

Figure 5 shows the numerical simulation results for such case. From Figure 5, we see that the observer (17) can also well track the states of system (3). Additionally, combined with the simulation results as shown in Figures 2–5, we can conclude that the designed observers have good performance under various kinds of inputs, and the input signal can affect the period of the circadian oscillator.

(a)

(b)

Figure 5: (a) State trajectories of systems (3) and (17) with step input signal in (19). (b) The error dynamics between systems (3) and (17). Here .

5. Discussions and Conclusions

Biological systems are typical complex dynamical systems. To efficiently infer the state of a biological system, it is necessary to develop some simple observers via monitoring a few system variables. Based on the recently proposed graphical approach, we have designed some rather simple observers for a core circadian rhythm network. For various input signals and under appropriate control gains, the designed observer can well infer the states of the original system. The investigations in this paper further support the conclusions in [4]. Real-world applications of the related investigations on biological networks include the experimental design and control of the related biological systems.

We have considered three types of inputs, and it is intriguing to investigate the observer design problems for stochastic systems, since biological systems are inherent stochastic and perturbed by environment [14–16]. Another question that deserves to be further investigated is to develop some general theorems to guarantee the observability of the biological systems [25–28]. Finally, it is also intriguing to investigate the observability of large-scale biological networks [17], such as the yeast cell cycle network with boolean dynamical model or differential equation models [18–23]. These topics will be discussed in our future works.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China under Grant no. 61304151 and the Science Foundation of Henan University under Grants nos. 2012YBZR007 and 2013YBRW005.

References

B. C. Moore, “Principal component analysis in linear systems: controllability, observability, and model reduction,” IEEE Transactions on Automatic Control, vol. 26, no. 1, pp. 17–32, 1981.

J. P. Gauthier and I. Kupka, Deterministic Observation Theory and Applications, Cambridge University Press, 2001.

S. Diop and M. Fliess, “Nonlinear observability, identifiability, and persistent trajectories,” in Proceedings of the 30th IEEE Conference on Decision and Control, vol. 1, pp. 714–719, December 1991.

Y.-Y. Liu, J.-J. Slotine, and A.-L. Barabási, “Observability of complex systems,” Proceedings of the National Academy of Sciences of the United States of America, vol. 110, no. 7, pp. 2460–2465, 2013.

R. Rajamani, “Observers for Lipschitz nonlinear systems,” IEEE Transactions on Automatic Control, vol. 43, no. 3, pp. 397–401, 1998.

J. Tsinias, “Observer design for nonlinear systems,” Systems and Control Letters, vol. 13, no. 2, pp. 135–142, 1989.

D. Gonze, J. Halloy, and A. Goldbeter, “Deterministic versus stochastic models for circadian rhythms,” Journal of Biological Physics, vol. 28, no. 4, pp. 637–653, 2002.

P. Wang, X. Yu, and J. Lü, “Identification and evolution of structurally dominant nodes in protein protein interaction networks,” IEEE Transactions on Biomedical Circuits and Systems, vol. 8, no. 1, pp. 87–97, 2014.

P. Wang and J. Lü, “Control of genetic regulatory networks: opportunities and challenges,” Acta Automatica Sinica, vol. 39, no. 12, pp. 1969–1979, 2013.

P. Wang, J. Lü, and M. J. Ogorzalek, “Global relative parameter sensitivities of the feed-forward loops in genetic networks,” Neurocomputing, vol. 78, no. 1, pp. 155–165, 2012.

A. Goldbeter, “A model for circadian oscillations in the Drosophila period protein (PER),” Proceedings of the Royal Society B: Biological Sciences, vol. 261, no. 1362, pp. 319–324, 1995.

A. Goldbeter, Biochemical Oscillations and Cellular Rhythms. The Molecular Bases of Periodic and Chaotic Behaviour, Cambridge University Press, Cambridge, UK, 1996.

J. Leloup, D. Gonze, and A. Goldbeter, “Limit cycle models for circadian rhythms based on transcriptional regulation in Drosophila and Neurospora,” Journal of Biological Rhythms, vol. 14, no. 6, pp. 433–448, 1999.

L. Wan, Q. Zhou, Z. Zhou, and P. Wang, “Dynamical behaviors of the stochastic Hopfield neural networks with mixed time delays,” Abstract and Applied Analysis, vol. 2013, Article ID 384981, 8 pages, 2013.

S. Wang, Y. Shen, C. Shi, T. Wang, Z. Wei, and H. Li, “Defining biological networks for noise buffering and signaling sensitivity using approximate Bayesian computation,” The Scientific World Journal, vol. 2014, Article ID 625754, 12 pages, 2014.

P. Wang, J. Lü, L. Wan, and Y. Chen, “A stochastic simulation algorithm for biochemical reactions with delays,” in Proceedings of the IEEE International Conference on Systems Biology, pp. 109–114, Huangshan, China, August 2013.

P. Wang, R. Lu, Y. Chen, and X. Wu, “Hybrid modelling of the general middle-sized genetic regulatory networks,” in Proceedings of the IEEE International Symposium on Circuits and Systems (ISCAS '13), pp. 2103–2106, Beijing, China, May 2013.

B. Novak and J. J. Tyson, “Modeling the control of DNA replication in fission yeast,” Proceedings of the National Academy of Sciences of the United States of America, vol. 94, no. 17, pp. 9147–9152, 1997.

A. Sveiczer, J. J. Tyson, and B. Novak, “Modelling the fission yeast cell cycle,” Briefings in Functional Genomics and Proteomics, vol. 2, no. 4, pp. 298–307, 2004.

F. Li, T. Long, Y. Lu, Q. Ouyang, and C. Tang, “The yeast cell-cycle network is robustly designed,” Proceedings of the National Academy of Sciences of the United States of America, vol. 101, no. 14, pp. 4781–4786, 2004.

B. Novak, Z. Pataki, A. Ciliberto, and J. J. Tyson, “Mathematical model of the cell division cycle of fission yeast,” Chaos, vol. 11, no. 1, pp. 277–286, 2001.

B. Novak and J. J. Tyson, “Quantitative analysis of a molecular model of mitotic control in fission yeast,” Journal of Theoretical Biology, vol. 173, no. 3, pp. 283–305, 1995.

B. Novak, K. C. Chen, and J. J. Tyson, “Systems biology of the yeast cell cycle engine,” in Systems Biology: Definitions and Perspectives, L. Alberghina and H. V. Westerhoff, Eds., vol. 13 of Topics in Current Genetics, pp. 305–324, Springer, Berlin, Germany, 2005.

P. Wang, D. Li, X. Wu, and J. Lü, “Ultimate bound estimation of a class of high dimensional quadratic autonomous dynamical systems,” International Journal of Bifurcation and Chaos, vol. 21, no. 9, pp. 2679–2694, 2011.

G. Wen, Z. Duan, G. Chen, and W. Yu, “Consensus tracking of multi-agent systems with Lipschitz-type node dynamics and switching topologies,” IEEE Transactions on Circuits and Systems I: Regular Papers, vol. 61, no. 2, pp. 499–511, 2014.

J. Qiu and J. Cao, “Global synchronization of delay-coupled genetic oscillators,” Neurocomputing, vol. 72, no. 16–18, pp. 3845–3850, 2009.

T. Zhou, J. Zhang, Z. Yuan, and L. Chen, “Synchronization of genetic oscillators,” Chaos, vol. 18, no. 3, Article ID 037126, 2008.

G. Wen, W. Yu, Y. Zhao, and J. Cao, “Pinning synchronisation in fixed and switching directed networks of Lorenz-type nodes,” IET Control Theory & Applications, vol. 7, no. 10, pp. 1387–1397, 2013.