Cooperative Control and Its Engineering Applications in Power SystemsView this Special Issue
Observer Design for a Core Circadian Rhythm Network
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.
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.  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.  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 . 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 , 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 . In 1999, Leloup et al. established a model for the circadian rhythms in the Drosophila and Neurospora . In 2002, Gonze and coauthors  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. , for the core circadian rhythm network  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. ; 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 .
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.  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 . 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.
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. , 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 , 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 . 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  proposed a general theorem for the observer design. However, due to the complexity of the biological model, the theorem obtained in  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.
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
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.
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 . 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 , 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.
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.
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.View at: Publisher Site | Google Scholar | MathSciNet
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.View at: Google Scholar
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.View at: Publisher Site | Google Scholar | MathSciNet
R. Rajamani, “Observers for Lipschitz nonlinear systems,” IEEE Transactions on Automatic Control, vol. 43, no. 3, pp. 397–401, 1998.View at: Publisher Site | Google Scholar | MathSciNet
J. Tsinias, “Observer design for nonlinear systems,” Systems and Control Letters, vol. 13, no. 2, pp. 135–142, 1989.View at: Publisher Site | Google Scholar | Zentralblatt MATH
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.View at: Publisher Site | Google Scholar
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.View at: Google Scholar
P. Wang and J. Lü, “Control of genetic regulatory networks: opportunities and challenges,” Acta Automatica Sinica, vol. 39, no. 12, pp. 1969–1979, 2013.View at: Google Scholar
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.View at: Publisher Site | Google Scholar
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.View at: Publisher Site | Google Scholar
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.View at: Publisher Site | Google Scholar
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.View at: Publisher Site | Google Scholar | MathSciNet
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.View at: Publisher Site | Google Scholar
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.View at: Google Scholar
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.View at: Publisher Site | Google Scholar
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.View at: Publisher Site | Google Scholar
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.View at: Publisher Site | Google Scholar
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.View at: Publisher Site | Google Scholar
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.View at: Publisher Site | Google Scholar | Zentralblatt MATH
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.View at: Publisher Site | Google Scholar
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.View at: Publisher Site | Google Scholar
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.View at: Publisher Site | Google Scholar | MathSciNet
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.View at: Publisher Site | Google Scholar | MathSciNet
J. Qiu and J. Cao, “Global synchronization of delay-coupled genetic oscillators,” Neurocomputing, vol. 72, no. 16–18, pp. 3845–3850, 2009.View at: Publisher Site | Google Scholar
T. Zhou, J. Zhang, Z. Yuan, and L. Chen, “Synchronization of genetic oscillators,” Chaos, vol. 18, no. 3, Article ID 037126, 2008.View at: Publisher Site | Google Scholar | MathSciNet
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.View at: Publisher Site | Google Scholar | MathSciNet