Robust synchronization control scheme of a population of nonlinear stochastic synthetic genetic oscillators under intrinsic and extrinsic molecular noise via quorum sensing
© Chen and Hsu; licensee BioMed Central Ltd. 2012
Received: 26 July 2012
Accepted: 10 October 2012
Published: 26 October 2012
Collective rhythms of gene regulatory networks have been a subject of considerable interest for biologists and theoreticians, in particular the synchronization of dynamic cells mediated by intercellular communication. Synchronization of a population of synthetic genetic oscillators is an important design in practical applications, because such a population distributed over different host cells needs to exploit molecular phenomena simultaneously in order to emerge a biological phenomenon. However, this synchronization may be corrupted by intrinsic kinetic parameter fluctuations and extrinsic environmental molecular noise. Therefore, robust synchronization is an important design topic in nonlinear stochastic coupled synthetic genetic oscillators with intrinsic kinetic parameter fluctuations and extrinsic molecular noise.
Initially, the condition for robust synchronization of synthetic genetic oscillators was derived based on Hamilton Jacobi inequality (HJI). We found that if the synchronization robustness can confer enough intrinsic robustness to tolerate intrinsic parameter fluctuation and extrinsic robustness to filter the environmental noise, then robust synchronization of coupled synthetic genetic oscillators is guaranteed. If the synchronization robustness of a population of nonlinear stochastic coupled synthetic genetic oscillators distributed over different host cells could not be maintained, then robust synchronization could be enhanced by external control input through quorum sensing molecules. In order to simplify the analysis and design of robust synchronization of nonlinear stochastic synthetic genetic oscillators, the fuzzy interpolation method was employed to interpolate several local linear stochastic coupled systems to approximate the nonlinear stochastic coupled system so that the HJI-based synchronization design problem could be replaced by a simple linear matrix inequality (LMI)-based design problem, which could be solved with the help of LMI toolbox in MATLAB easily.
If the synchronization robustness criterion, i.e. the synchronization robustness ≥ intrinsic robustness + extrinsic robustness, then the stochastic coupled synthetic oscillators can be robustly synchronized in spite of intrinsic parameter fluctuation and extrinsic noise. If the synchronization robustness criterion is violated, external control scheme by adding inducer can be designed to improve synchronization robustness of coupled synthetic genetic oscillators. The investigated robust synchronization criteria and proposed external control method are useful for a population of coupled synthetic networks with emergent synchronization behavior, especially for multi-cellular, engineered networks.
Many biochemically dynamical systems are controlled by intrinsic rhythms generated by specialized cellular clocks within the organism itself. These rhythm generators are composed of thousands of clock cells that are intrinsically diverse, but nevertheless manage to function in a coherent oscillatory state [1–6]. However, the synchronization mechanisms by which this collective behavior arises remains to be understood, even if individual clock cells are known to operate through biochemical networks comprising of multiple regulatory feedback loops [7–11]. The complexity of these cellular systems has hindered a complete understanding of natural genetic oscillators and their synchronization [12–16].
Recently designed synthetic genetic oscillators can offer an alternative approach, and provide a relatively well-controlled test bed in which the function and behavior of natural genetic oscillators can be isolated and characterized in detail [1, 3, 7, 8, 17]. As an example, a synthetic biological oscillator, termed the “repressilator,” was developed in Escherichia coli from a network of three transcriptional repressors that inhibit one another in a cyclic way [18–21]. Spontaneous oscillations were observed in individual cells within a growing culture, although substantial variability and noise were present among the different cells. Another synthetic oscillator was designed and built that exhibited damped oscillatory responses to perturbations in culture . Recently, several mechanisms of intercell coupling of synthetic genetic oscillators have been discussed, to enhance the oscillating response of the synthetic biological system . In general, coupling among oscillators is not sufficient to achieve synchronization, and many ensembles of coupled oscillators exhibit phase dispersion rather than a synchronized state, because the oscillators may actively resist oscillation or because the coupling is too small. Therefore, the synchronization of a population of nonlinear stochastic coupled oscillators must be analyzed carefully.
In previous studies [1, 12–16, 23–26], the collective behavior of synthetic genetic oscillators was discussed on the basis of cell-to-cell communication through quorum sensing, which may lead to synchronization in an ensemble of identical genetic oscillators. In general, intercellular communication is accomplished by transmitting individual cell reactions via intercellular signals to neighboring cells, which can generate a global cellular synchronization at the level of molecules, tissues, organs, or the body . The ability to communicate among cells is an absolute requisite to ensure appropriate and robust synchronization at all levels in organisms living in uncertain environments. Synchronization of coupled networks has been investigated intensively in past decades because of its biological implications and potential applications [3, 4, 20]. The global synchronization mechanism of oscillators via direct coupling has been derived based on the Lyapunov method and linearization technique [27–29]. Different synchronization mechanisms in a population of nonlinear stochastic genetic oscillators with noise and impulse control inputs have also been previously discussed [7–9].
Generally, biological systems or organisms are subject to time-varying uncertainties, assumed to be in the form of both internal noise resulting from the birth and death of biochemical molecules, and external noise resulting from environmental perturbations [15, 19, 30–34]. It has been shown that environmental molecular noise plays an important role in the stochastic behavioral phenomena of biological systems at various levels. In particular, gene regulation is an inherently noisy process [35–37], from transcriptional control, alternative splicing, translation, and diffusion to biological modification of transcription factors. Such stochastic noise cannot only significantly affect the dynamics of biological systems, but also may be exploited by living organisms to actively facilitate certain cellular functions, such as cellular communication and synchronization. In this study, the time-varying parameter fluctuations in gene regulation processes are modeled as stochastic intrinsic kinetic noise, whereas environmental molecular noise of the host cells is modeled as an extrinsic disturbance of synthetic genetic oscillators. A population of synthetic oscillators coupled by quorum sensing is modeled by a set of coupled nonlinear stochastic equations with intrinsic and extrinsic noise. Our main purpose is to discuss the robust synchronization mechanism of a population of nonlinear stochastic coupled synthetic genetic oscillators under intrinsic kinetic parameter fluctuations and extrinsic environmental molecular noise on the host cells.
Based on nonlinear stochastic equations of coupled synthetic genetic oscillators distributed in different host cells, the robust synchronization mechanism is discussed from the H ∞ noise-filtering perspective. The robust ability to tolerate stochastic kinetic parameter fluctuations and the filtering ability to attenuate extrinsic environmental molecular noise to maintain synchronization of nonlinear stochastic coupled synthetic genetic oscillations was measured from the nonlinear stochastic system theory point of view. In the case where robust synchronization cannot be maintained or is corrupted by intrinsic kinetic parameter fluctuations and extrinsic environmental molecular noise, some robust synchronization design methods are discussed to enhance synchronization. Both the physical insight into the robust synchronization mechanisms and the designs to improve these mechanisms require solving nonlinear HJIs, which cannot be easily achieved by any analytical or numerical method at present. In order to simplify this analysis and design, the Takagi-Sugeno (T-S) fuzzy model [38–40] is employed to interpolate several local linear stochastic genetic oscillators to approximate the nonlinear stochastic genetic oscillators. Consequently, the analysis and design of robust synchronization can be achieved by solving of a set of LMIs , via the help of the LMI toolbox in MATLAB.
If the robust synchronization of nonlinear stochastic coupled synthetic genetic oscillators could not be achieved spontaneously, an external control input was developed to synchronize the coupled synthetic genetic oscillators. External stimulation inputs are known to play an important role in the synchronization of biological rhythms. For instance, many organisms display a circadian rhythm of 24-hours periodicity entrained to the light–dark cycle . Other examples include physiological rhythms stimulated by regular or periodic inputs occurring in the context of medical devices, synchronization of electronic genetic networks by an external forcing, i.e., external voltage , and a wide variety of regular and irregular rhythms induced by periodic stimulation of squid giant axons . In general, physiological oscillations can be synchronized by appropriate external or internal stimuli. Recently, Wang et al. constructed an impulse control system to model the process of periodical injection of coupling substances with constant or random impulse control into a common extracellular medium and studied its effect on the dynamics of collective rhythms . Synchronization with time delays via direct coupling and linearization methods has also been discussed [29, 43]. In our study, a control input design method was developed to guarantee the robust synchronization of synthetic genetic oscillators under intrinsic kinetic parameter fluctuations and extrinsic environmental molecular noise. In order to avoid solving the complicated HJI, the fuzzy interpolation method was also employed to simplify the control design procedure by only solving a set of simple LMIs .
The contributions of this paper are fourfold: (1) A nonlinear stochastic system is introduced to model a population of nonlinear stochastic coupled synthetic genetic oscillators under random intrinsic kinetic parameter fluctuations and extrinsic molecule noise in vivo; (2) The conditions of robust synchronization are developed from the nonlinear stochastic dynamical system point of view, i.e. if the synchronization robustness ≥ intrinsic robustness + extrinsic robustness, then the intrinsic parameter fluctuation can be tolerated and the extrinsic noise can be buffered so that the robust synchronization of coupled oscillators can be guaranteed. Therefore, we could obtain better insight into synchronization mechanisms of coupled synthetic molecular systems distributed in host cells, and to provide further systematic analysis and control design to improve the synchronization robustness. Further, if the robust synchronization conditions cannot be guaranteed, some synchronization control schemes are also developed to improve the synchronization robustness of coupled synthetic gene networks; (3) The fuzzy interpolation method is introduced to simplify the analysis and design procedure of robust synchronization of coupled nonlinear stochastic synthetic gene networks; and (4) An external inducer input control design method is also developed to guarantee robust synchronization of the nonlinear stochastic synthetic genetic oscillators, when spontaneous synchronization cannot be achieved under intrinsic kinetic parameter fluctuations and extrinsic environmental molecular noise. Finally, a design example is provided in silico to illustrate the design procedure and to confirm the performance of the proposed design methods.
Stochastic models of nonlinear stochastic synthetic genetic oscillators under intrinsic kinetic parameter fluctuations and external molecular noise
where x ai , x bi and x ci are the concentration of mRNA transcribed from tetR, cI and lacI in cell i, respectively; concentrations of the corresponding proteins are represented by x Ai , x Bi and x Ci , respectively. The concentration of AI inside each cell is denoted by x Si . α a , α b , and α c are the dimensionless transcription rate in the absence of the repressor, and μ is the repression coefficient. α S is the maximal contribution to lacI transcription in the presence of saturating amounts of AI, and μ S is the activation coefficient. γ m is the respective dimensionless degradation rate of mRNA for tetR, cI and lacI in the cell, and n is the Hill coefficient .
where x i (t) = (x ai (t), x bi (t), x ci (t), x Ai (t), x Bi (t), x Ci (t), x Si (t)) T ∈ R m is the state vector of the i th synthetic oscillator; f(·) : R m → R m is a smooth nonlinear function that characterizes the behavior of the synthetic oscillator; g(·) : R m → R m is a smooth nonlinear inner-coupling function; and C = (c ij )N × N ∈ RN × N is the coupling configuration matrix, where c ij >0 means that the i th synthetic oscillator is coupled with the j th synthetic oscillator directly, otherwise c ij =0. Assume that the diagonal elements of C satisfy .
Synthetic genetic oscillators under intrinsic parameter fluctuations
where w i (t) is a standard Wiener process or Brownian motion with dw i (t) = n i (t)dt to represent the random kinetic fluctuations of synthetic gene circuits.
Synthetic genetic oscillators under external environmental molecular noise
where the signal vector v i (t) denotes extrinsic environmental molecular noise from the environment and h i denotes the noise-coupling matrix in cell i.
Robust synchronization control design for nonlinear synthetic genetic oscillators under intrinsic kinetic parameter fluctuations and external environmental molecular noise
for all possible environmental noises v(t), where ρ2 denotes the upper bound of the effect of v(t) on synchronization error from the mean energy point of view, i.e., the noise filtering level ρ in (14) denotes the upper bound of the noise-filtering ability ρ0 of coupled synthetic genetic oscillators, and R is a symmetric weighting matrix to be specified by designer.
The inequality in (14) means that the effect of extrinsic environmental molecular noise on the synchronization error is less than ρ from the mean energy point of view, i.e., the value of the noise-filtering ability ρ 0 is the lower bound of ρ. Because the statistics of extrinsic molecular noise may be unavailable or uncertain, it is very difficult to obtain the noise filtering ability ρ 0 for all possible extrinsic noises v(t) directly and only the upper bound ρ of the noise-filtering ability ρ 0 can be given in (14) at first. Then, we will decrease the upper bound ρ to as small a value as possible to approach its lower bound for the noise-filtering ability ρ 0, i.e., to get ρ 0 by minimizing ρ (or ρ 0 = min ρ) indirectly. If the noise filtering ability ρ 0 is small, it means that the environmental molecular noise has less influence on synchronization and vice versa. It will be further discussed in the sequel.
Based on the synchronization error dynamic equation in (13) and the H∞ noise filtering performance in (14) or (15), we obtain the following robust synchronization result for nonlinear stochastic coupled synthetic genetic oscillators under intrinsic kinetic parameter fluctuations and extrinsic environmental molecular noise.
then the stochastic intrinsic noise can be tolerated (i.e. the synchronization of coupled synthetic genetic networks cannot destroyed by intrinsic noise) and the influence of extrinsic environmental molecular noise v(t) on the synchronization of the nonlinear stochastic coupled synthetic oscillation systems in (8) is less than or equal to a prescribed filtering level ρ, i.e., the inequality in (14) holds.
Proof: see Additional file 1: Appendix A
subject to HJI in (16) with V(e) > 0
i.e., the noise-filtering ability ρ0 on v(t) for the synchronized synthetic oscillators can be evaluated by solving the constrained optimization in (17). The noise-filtering ability ρ0 of the synchronized synthetic oscillators in (17) can be obtained by decreasing ρ until no positive solution V(e) > 0 exists for HJI in (16) again.
subject to HIJ in (16) with V(e) > 0
Synchronization robustness: The ability of coupled synthetic genetic oscillators to resist both intrinsic noise and extrinsic noise so that the synchronization can be maintained.
Intrinsic robustness: The ability of coupled synthetic genetic oscillator to tolerate intrinsic parameter fluctuation to maintain synchronization.
Extrinsic robustness: The filtering ability to attenuate the effect of environmental noise on the synchronization of coupled synthetic genetic network.
The first term on the left hand side of (19) indicates the intrinsic robustness to tolerate the intrinsic parameter fluctuation in (13) because this term is induced by intrinsic noise (or random parameter fluctuation), the second and third term on the left hand side are due to the noise filtering in (14) and indicate the extrinsic robustness to filter the extrinsic noise with the noise filtering ability ρ0, and the term on the right hand side of (19) indicates the synchronization robustness of the coupled synthetic gene networks. The biological meaning of synchronization robustness criterion in (19) is that if the synchronization robustness can confer both the intrinsic robustness to tolerate intrinsic parameter fluctuation and extrinsic robustness to filter the environmental noise, then the coupled synthetic networks will synchronize with a noise filtering ability ρ0. If the synchronization robustness criterion in (19) is violated, then the synchronization of coupled synthetic gene networks may not be achieved due to the intrinsic parameter fluctuation and extrinsic noise.
In general, it is still very difficult to solve the second-order HJI in (16) with V(e) > 0 and V(0) = 0 to guarantee robust synchronization of nonlinear stochastic coupled synthetic genetic oscillators with a prescribed attenuation level ρ under intrinsic kinetic parameters fluctuations and extrinsic environmental molecular noise or to solve the constrained minimization in (18) for robust synchronization design to achieve the optimal molecular noise filtering of the synchronized coupled synthetic oscillators. Recently, the fuzzy dynamic model has been widely used to interpolate several local dynamic models to efficiently approximate a nonlinear dynamic system [32, 38, 39]. Hence, in this situation, we employ the T-S fuzzy model to interpolate several linear synthetic stochastic oscillators at different local operation points to efficiently and globally approximate the error dynamic in (13), so that the analysis and design procedure for robust synchronization of nonlinear stochastic coupled synthetic genetic oscillators can be simplified.
Robust synchronization design of synthetic genetic oscillators via t-s fuzzy methodology
In this study, the T-S fuzzy method is employed to simplify the analysis and design procedure for robust synchronization of nonlinear stochastic coupled synthetic oscillators under intrinsic kinetic parameter fluctuations and extrinsic environmental molecular noise. The T-S fuzzy model for the synchronization error dynamics is described by fuzzy if-then rules. The k th rule of the fuzzy model for the synchronization error dynamics for cell i in (12) is proposed in the following form [30, 32, 38, 39]:
where μ k (z) = diag(μk,1(z1), …, μk,N(z N )) and z = z1, …, z N T .
In , Takagi and Sugeno have proposed the systematic method to build T-S fuzzy model for nonlinear function approximation by the system identification tool, i.e. the local system matrix A k , B k , A Wk , and B Wk in (21) or (22) can be identified by least square estimation method. On the other hand, many studies have proved that the T-S fuzzy model can approximate a continuous function with any degree of accuracy. Actually, there is still some fuzzy approximation error in (22). In the robust synchronization control design, for simplicity, the fuzzy approximation error can be merged into the external noise, which could be efficiently attenuated by the proposed H∞ robust synchronization control design in the sequel.
After investigating the approximation of nonlinear stochastic coupled synthetic oscillators by the fuzzy interpolation method, in order to avoid solving the nonlinear constrained optimization problem in (18) for the robust synchronization design problem of coupled synthetic oscillators under intrinsic kinetic parameter fluctuation and extrinsic environmental molecular noise, the measurement procedure for the noise-filtering ability of synchronized synthetic genetic oscillators could also be simplified by the fuzzy approximation method. Then, we get the following result.
for k = 1, 2, …, L, then the noise-filtering level ρ in (14) holds, or the intrinsic kinetic parameter fluctuations are tolerated by the synchronized synthetic genetic oscillators, and the influence of extrinsic environmental molecular noise v(t) on the synchronized synthetic oscillation systems in (13) is less than or equal to a prescribed filtering level ρ.
Proof: See Additional file 1: Appendix B
subject to P>0 and LMIs in (23)
If the prescribed noise-filtering level ρ is prescribed by a biological engineer, a robust synchronization design would involve specifying the design parameters α a , α b , α c , β A , β B , β C , γ m, and γ p of the synthetic gene oscillators in A k , so that the LMIs in (23) have a positive solution P>0 in Proposition 2. If we want to achieve optimal filtering of extrinsic noise for the synchronized synthetic oscillators, some design parameters need to be specified for coupled synthetic oscillators to achieve the constrained optimization in (24).
In this study, the fuzzy approximation method in (21) or (22) is only employed to simplify the analysis and design procedure via solving P>0 for LMIs in (23) instead of solving V(e) > 0 for HJI in (16) directly. Further, based on the fuzzy interpolation of local linear systems, i.e., replacing F(x, s), G(x, s), F W (x, s) and G W (x, s) by the fuzzy approximations in (22), in Proposition 2, V(e) = e T Pe is employed to solve the HJI (16) in Proposition 1. The HJI in Proposition 1 is replaced with a set of LMIs in Proposition 2 and we only need to solve P>0 for LMIs to guarantee the coupled synthetic genetic oscillators have a noise filtering level ρ.
In general, the constrained optimization problems in (24) are called eigenvalue problem , which can be efficiently solved by the MATLAB LMI toolbox.
In addition to the robust oscillation synchronization, the proposed method can be applied to robust synchronization design of coupled synthetic gene networks with any kind of dynamic behavior.
- (v)In the fuzzy approximation case, the synchronization robustness criterion in (19) is equivalent to the following(25)
for k = 1, 2, …, L which is equivalent to (23) with ρ being replace by ρ0. The biological meaning of synchronization robustness criterion in (25) is that if the local synchronization robustness of local coupled synthetic genetic oscillators can confer local intrinsic robustness to tolerate local intrinsic parameter fluctuation and local extrinsic robustness to filter external noise, then the coupled synthetic genetic oscillators can be synchronized with a noise filtering ability ρ0. If the synchronization robustness criterion in (25) is violated, then the synchronization of coupled synthetic genetic oscillators may not be achieved due to intrinsic parameter fluctuation and extrinsic noise. In general, if the design parameters of coupled synthetic genetic oscillators are specified so that the eigenvalues of local coupled system matrix I N ⊗ A Wk + C ⊗ B Wk are far in the left hand side of complex s-domain (i.e. with more negative real part), then the coupled synthetic genetic networks are more easy to synchronize in spite of intrinsic parameter fluctuation and extrinsic noise.
Robust synchronization of synthetic genetic oscillators by external control input
where u e = Q represents an extracellular control input, which can be implemented via the injection of inducer AI.
where C e (Q) = (c eij (Q))N × N ∈ RN × N is the coupling configuration matrix, in which c eii (Q) = − η s (1 − N− 1)(Q e + Q) if i=j, otherwise c eij (Q) = η s N− 1(Q e + Q). Then, we can also obtain the robust synchronization control design of coupled oscillation systems under intrinsic kinetic parameter fluctuations and extrinsic environmental molecular noises as follows.
for a prescribed filtering level ρ, then the stochastic intrinsic kinetic noise can be robustly tolerated, and the influence of extrinsic environmental molecular noise v(t) on the synchronization of the nonlinear stochastic coupled synthetic oscillation systems in (27) is less than or equal to ρ, i.e., the inequality in (14) or (15) holds.
Proof: similar to the proof of proposition 1
subject to V(e) > 0 and HJI in (28)
Applying the fuzzy approximation method, the external signal control design can be obtained as described in the following corollary, for robust filtering of synchronized oscillation systems with intrinsic kinetic parameter fluctuations and extrinsic environmental molecular noise.
for k = 1, 2, ⋯, L, then intrinsic parametric noise can be tolerated and the effect of extrinsic molecular noise v(t) on the synchronization of nonlinear stochastic coupled oscillation systems is less than or equal to a prescribed filtering level ρ.
Proof: similar to the proof of proposition 2
The physical meaning of Corollary 4 is that if we can select a control parameter Q, such that the LMIs in (32) have a positive definite solution P > 0, then the robust synchronization with a prescribed noise filtering level ρ on extrinsic environmental molecular noise is guaranteed for the nonlinear stochastic coupled synthetic genetic oscillators. If we specify control parameter Q so that the eigenvalues of local system matrix of I N ⊗ A k + C e (Q) ⊗ B k of coupled synthetic gene oscillators have more negative real part (i.e. in far left hand complex s-domain), the coupled synthetic gene oscillators are with more robust synchronization to tolerate more intrinsic parameter fluctuations and to filter more extrinsic noise.
subject to V(e) > 0 and LMIs in (32)
The physical meaning of the constrained optimization in (33) is that if we can select a control parameter Q through the inducer concentration control method to solve the constrained optimization problem, we can achieve both robust synchronization against intrinsic kinetic parameter fluctuations and optimal filtering against external environmental molecular noise on the synchronization by using the external control signal in the coupled synthetic oscillation systems.
Consider a synthetic genetic network of N coupled oscillators with intrinsic kinetic parameter fluctuations and extrinsic environmental molecular noise.
Given the prescribed disturbance attenuation level ρ.
Represent the nonlinear stochastic synchronization error dynamic by the T-S fuzzy synchronization error dynamic model, using the interpolation of several local linear stochastic systems.
Specify Q to solve LMI in (32) with the help of LMI toolbox in MATLAB so that N coupled synthetic genetic oscillators can be synchronized with a prescribed noise filtering level ρ.
Design examples in silico for robust synchronization design in the genetic oscillation systems
In this section, we provide a simulated example to illustrate the design procedure of robust synchronization of the nonlinear stochastic coupled synthetic oscillation systems and to confirm the performance of the robust synchronization of proposed method against intrinsic kinetic parameter fluctuations and extrinsic environmental molecular noise.
The cell is the functional unit of all living things, either unicellular or multicellular. A cell can sense many different signals from the internal or external context and can respond to the constantly changing environment via appropriate cellular processes. Also, cells can interact with each other via cell-to-cell communication and achieve specific physiological functions essential for life in a cooperative manner. However, many fundamental questions remain regarding how cellular phenomena arise from the interactions between genes and proteins, what features make the cell operate reliably in diverse conditions, and how the cell is responsible for these operations. To gain insight into these questions, one can construct the underlying mechanisms that constitute the web of interactions. This idea is useful to separate a complicated network into many simpler ones, which can work independently but also cooperate with each other. It may not only enhance our understanding of collective behavior particularly via synchronization but may also establish a foundation to design robust implementation of coupled synthetic gene networks [47, 48].
In this paper, we consider a nonlinear stochastic coupled network with two or more coupled synthetic oscillators. By transforming nonlinear stochastic coupled network dynamics into synchronization error dynamics, we can use Lyapunov’s direct method to infer a sufficient condition required for robustness of the nonlinear synchronized network. Assuming that each synthetic oscillator suffers from intrinsic kinetic parameter fluctuations and extrinsic molecular noise, robust synchronization performance is defined as the effect of extrinsic molecular noise upon the synchronization error. Based on this definition, robust synchronization performance of a nonlinear coupled network can be calculated by solving an associated HJI-constrained optimization problem. We also show that nonlinear coupled networks with robust synchronization performance are also synchronizable. Based on this synchronization performance, we propose a procedure for designing or compensating a coupled network with two or more coupled synthetic oscillator through a given connected topology toward a desired robust synchronization performance. Using the proposed method, the coupled synthetic oscillators cannot only tolerate kinetic parameter variations but also attenuate the extrinsic molecular noise below a desired level to achieve a desired robust synchronization. However, the HJI-constrained optimization problem is difficult to solve directly by any analytical or numerical method because of the complexity of nonlinear synchronization error dynamics. Hence, we employ a T–S fuzzy model to solve the HJI easily and indirectly. The T–S fuzzy model has been widely applied to approximate nonlinear systems by interpolating several local linearized systems. Here, we use the T–S fuzzy model to approximate the nonlinear stochastic synchronization error dynamics. By using the T–S fuzzy model and choosing the appropriate Lyapunov function, the HJI-constrained optimization for calculating the robust synchronization performance of a nonlinear coupled network is reduced to an equivalent LMI-constrained optimization problem, which can be solved efficiently by MATLAB’s LMI toolbox.
Recently a simple synthetic device was engineered in a cell, and several cells were then combined, so that their connections allowed the construction of a more complex synthetic gene circuit, i.e., so-called multi-cellular engineered networks. This approach not only uses cellular consortia as an efficient way of engineering complex gene circuits, but also demonstrates the great potential for reutilization of small parts of the gene circuit. In such situations, our proposed evaluation framework may offer a possible guideline for the design or compensation of such coupled networks with a given connected topology toward a desired collective behavior.
In this study, several robust synchronization criteria and designs are proposed for a population of synthetic genetic oscillators in order to exploit an emergent synchronization phenomenon by quorum sensing molecules under intrinsic kinetic parameter fluctuations and extrinsic molecular noise. When the synchronization of nonlinear stochastic coupled synthetic genetic oscillators cannot be maintained, a robust H∞ control scheme is developed to enhance synchronization by adding external control to increase the cell-to-cell communication through quorum sensing. This study enhances our understanding in this area in the following ways: (a) nonlinear stochastic systems are employed to model the coupled synthetic genetic oscillators with intrinsic kinetic parameter fluctuations, extrinsic molecular noises on the host cells and quorum sensing molecules; (b) two robust synchronization criteria (19) and (25) of coupled synthetic oscillators are developed from the nonlinear stochastic filtering point of view, so that we can gain more insight into the robust synchronization mechanism from a systematic perspective. If these robust synchronization criterion cannot be guaranteed, robust synchronization control schemes via selecting adequate kinetic parameters and inducer concentration are also developed to improve the synchronization robustness of coupled synthetic genetic oscillators; and (c) the fuzzy approximation is employed to approximate the nonlinear stochastic synchronization error model by interpolating several linear stochastic systems, so that the powerful LMI toolbox in MATLAB can be used to significantly simplify the system analysis and design procedure for robust synchronization of coupled synthetic gene networks, which is very important for the emergent phenomenon of synthetic gene networks through quorum sensing at the molecular level. In addition, the proposed robust synchronization design and control scheme for nonlinear stochastic coupled genetic oscillators can be easily applied to the robust synchronization problem of other coupled genetic networks as a cellular consortium.
The work was supported by the National Science Council of Taiwan under grant NSC 100-2745-E-007-001-ASP and NSC 101-2745-E-007-001-ASP.
- Garcia-Ojalvo J, Elowitz MB, Strogatz SH: Modeling a synthetic multicellular clock: repressilators coupled by quorum sensing. Proc Natl Acad Sci USA. 2004, 101: 10955-10960. 10.1073/pnas.0307095101.View ArticleGoogle Scholar
- Gardner TS, Cantor CR, Collins JJ: Construction of a genetic toggle switch in E. coli. Nature. 2000, 403: 339-342. 10.1038/35002131.View ArticleGoogle Scholar
- Glass L: Synchronization and rhythmic processes in physiology. Nature. 2001, 410: 277-284. 10.1038/35065745.View ArticleGoogle Scholar
- Gonze D, Bernard S, Waltermann C, Kramer A, Herzel H: Spontaneous synchronization of coupled circadian oscillators. Biophys J. 2005, 89: 120-129. 10.1529/biophysj.104.058388.View ArticleGoogle Scholar
- Wolf J, Heinrich R: Dynamics of two-component biochemical systems in interacting cells; Synchronization and desynchronization of oscillations and multiple steady states. Biosystems. 1997, 43: 1-24. 10.1016/S0303-2647(97)01688-2.View ArticleGoogle Scholar
- Wolf J, Passarge J, Somsen OJG, Snoep JL, Heinrich R, Wasterhoff HV: Transduction of intracellular and intercellular dynamics in yeast glycolytic oscillations. Biophys J. 2000, 78: 1145-1153. 10.1016/S0006-3495(00)76672-0.View ArticleGoogle Scholar
- Wang RQ, Chen LN: Synchronizing genetic oscillators by signaling molecules. J Biol Rhythm. 2005, 20: 257-269. 10.1177/0748730405275653.View ArticleGoogle Scholar
- Wang RQ, Jing ZJ, Chen LN: Modelling periodic oscillation in gene regulatory networks by cyclic feedback systems. Bull Math Biol. 2005, 67: 339-367. 10.1016/j.bulm.2004.07.005.View ArticleGoogle Scholar
- Wang RQ, Chen LN, Aihara K: Synchronizing a multicellular system by external input: an artificial control strategy. Bioinformatics. 2006, 22: 1775-1781. 10.1093/bioinformatics/btl182.View ArticleGoogle Scholar
- Chen LN, Wang RS, Zhang XS: Biomolecular Networks: Methods and Application in Systems Biology. 2009, John Wiley&Sons, New YorkView ArticleGoogle Scholar
- Chen LN, Wang RS, Li C, Aihara K: Modeling Biomolecular Network in Cells: Structures and Dynamics. 2010, Springer, LondonView ArticleGoogle Scholar
- Csermely P: Creative elements: network-based predictions of active centres in proteins and social networks. Trends Biochem Sci. 2008, 33: 569-576. 10.1016/j.tibs.2008.09.006.View ArticleGoogle Scholar
- Csermely P: Weak links: The universal key to the ability of networks and complex systems. 2009, Springer, HeidelbergView ArticleGoogle Scholar
- Kovacs IA, Palotai R, Szalay MS, Csermely P: Community landscapes: an integrative approach to determine overlapping network module hierarchy, indentify key nodes and predict network dynamics. PLoS One. 2010, 5: e12528-10.1371/journal.pone.0012528.View ArticleGoogle Scholar
- Zhou TS, Chen LN, Aihara K: Molecular communication through stochastic synchronization induced by extracellular fluctuations. Phys Rev Lett. 2005, 95: 178103-View ArticleGoogle Scholar
- Zhou TS, Zhang JJ, Yuan ZJ, Chen LN: Synchronization of genetic oscillators. Chaos. 2008, 18: 037126-10.1063/1.2978183.View ArticleGoogle Scholar
- Chen LN, Wang RQ, Zhou TS, Aihara K: Noise-induced cooperative behavior in a multicell system. Bioinformatics. 2005, 21: 2722-2729. 10.1093/bioinformatics/bti392.View ArticleGoogle Scholar
- Elowitz MB, Leibler S: A synthetic oscillatory network of transcriptional regulators. Nature. 2000, 403: 335-338. 10.1038/35002125.View ArticleGoogle Scholar
- Elowitz MB, Levine AJ, Siggia ED, Swain PS: Stochastic gene expression in a single cell. Science. 2002, 297: 1183-1186. 10.1126/science.1070919.View ArticleGoogle Scholar
- Hasty J, McMillen D, Collins JJ: Engineered gene circuits. Nature. 2002, 420: 224-230. 10.1038/nature01257.View ArticleGoogle Scholar
- Kaplan DT, Clay JR, Manning T, Glass L, Guevara MR, Shrier A: Subthreshold dynamics in periodically stimulated squid giant axons. Phys Rev Lett. 1996, 76: 4074-4077. 10.1103/PhysRevLett.76.4074.View ArticleGoogle Scholar
- Atkinson MR, Savageau MA, Myers JT, Ninfa AJ: Development of genetic circuitry exhibiting toggle switch or oscillatory behavior in Escherichia coli. Cell. 2003, 113: 597-607. 10.1016/S0092-8674(03)00346-5.View ArticleGoogle Scholar
- You LC, Cox RS, Weiss R, Arnold FH: Programmed population control by cell-cell communication and regulated killing. Nature. 2004, 428: 868-871. 10.1038/nature02491.View ArticleGoogle Scholar
- Pai A, Tanouchi Y, Collins CH, You L: Engineering multicellular systems by cell–cell communication. Curr Opin Biotech. 2009, 20: 461-470. 10.1016/j.copbio.2009.08.006.View ArticleGoogle Scholar
- Gonze D: Coupling oscillations and switches in genetic networks. Biosystems. 2010, 99: 60-69. 10.1016/j.biosystems.2009.08.009.View ArticleGoogle Scholar
- Zhang J, Yuan Z, Zhou T: Synchronization and clustering of synthetic genetic networks: a role for cis-regulatory modules. Phys Rev E. 2009, 79: 041903-View ArticleGoogle Scholar
- Lu JH, Chen GR: A time-varying complex dynamical network model and its controlled synchronization criteria. IEEE Trans Auto Control. 2005, 50: 841-846.View ArticleGoogle Scholar
- Lu JH, Yu XH, Chen GR, Cheng DZ: Characterizing the synchronizability of small-world dynamical networks. IEEE Trans Circ Syst I. 2004, 51: 787-796. 10.1109/TCSI.2004.823672.View ArticleGoogle Scholar
- Lu WL, Chen TP: Synchronization of coupled connected neural networks with delays. IEEE Trans Circ Syst I. 2004, 51: 2491-2503. 10.1109/TCSI.2004.838308.View ArticleGoogle Scholar
- Chen BS, Zhang WH: Stochastic H-2/H∞ control with state-dependent noise. IEEE Trans Auto Control. 2004, 49: 45-57. 10.1109/TAC.2003.821400.View ArticleGoogle Scholar
- Chen BS, Wang YC: On the attenuation and amplification of molecular noise in genetic regulatory networks. BMC Bioinforma. 2006, 7: 52-10.1186/1471-2105-7-52.View ArticleGoogle Scholar
- Chen BS, Tseng CS, Uang HJ: Robustness design of nonlinear dynamic systems via fuzzy linear control. IEEE Trans Fuzzy Syst. 1999, 7: 571-585. 10.1109/91.797980.View ArticleGoogle Scholar
- Teramae J, Tanaka D: Robustness of the noise-induced phase synchronization in a general class of limit cycle oscillators. Phys Rev Lett. 2004, 93: 204103-View ArticleGoogle Scholar
- Thattai M, van Oudenaarden A: Intrinsic noise in gene regulatory networks. Proc Natl Acad Sci USA. 2001, 98: 8614-8619. 10.1073/pnas.151588598.View ArticleGoogle Scholar
- Blake WJ, Kærn M, Cantor CR, Collins JJ: Noise in enkaryotic gene expression. Nature. 2003, 422: 633-637. 10.1038/nature01546.View ArticleGoogle Scholar
- McAdams HH, Arkin A: It's a noisy business! Genetic regulation at the nanomolar scale. Trends Genet. 1999, 15: 65-69. 10.1016/S0168-9525(98)01659-X.View ArticleGoogle Scholar
- Li C, Chen L, Aihara K: Synchronization of coupled nonidentical genetic oscillators. Phys Biol. 2006, 3: 37-44. 10.1088/1478-3975/3/1/004.View ArticleGoogle Scholar
- Passino KM, Yurkovich S: Fuzzy control. 1998, Addison-Wesley, Menlo ParkGoogle Scholar
- Takagi T, Sugeno M: Fuzzy identification of systems and its applications to modeling and control. IEEE Trans Syst Man Cyber. 1985, 15: 116-132.View ArticleGoogle Scholar
- Tseng CS, Chen BS: H∞ decentralized fuzzy model reference tracking control design for nonlinear interconnected systems. IEEE Trans Fuzzy Syst. 2001, 9: 795-809. 10.1109/91.971729.View ArticleGoogle Scholar
- Boyd S, El Ghaoui L, Feron E, Balakrishnan V: Linear Matrix Inequalities in System and Control Theory. 1994, Society for Industrial & Applied Mathematics, PhiladelphiaView ArticleGoogle Scholar
- Schultz TF, Kay SA: Circadian clocks in daily and seasonal control of development. Science. 2003, 301: 326-328. 10.1126/science.1085935.View ArticleGoogle Scholar
- Wagemakers A, Buldu JM, Garcia-Ojalvo J, Sanjuan MAF: Synchronization of electronic genetic networks. Chaos. 2006, 16: 013127-10.1063/1.2173048.View ArticleGoogle Scholar
- Chen BS, Wu WS: Robust filtering circuit design for stochastic gene networks under intrinsic and extrinsic molecular noises. Math Biosciences. 2008, 211: 342-355. 10.1016/j.mbs.2007.11.002.View ArticleGoogle Scholar
- Zhang W, Chen BS, Tseng CS: Robust H∞ filtering for nonlinear stochastic systems. IEEE Trans Signal Proc. 2005, 53: 589-598.View ArticleGoogle Scholar
- Tseng CS, Chen BS, Uang HJ: Fuzzy tracking control design for nonlinear dynamic systems via T-S fuzzy model. IEEE Trans Fuzzy Syst. 2001, 9: 381-392. 10.1109/91.928735.View ArticleGoogle Scholar
- Endy D: Foundations for engineering biology. Nature. 2005, 438: 449-453. 10.1038/nature04342.View ArticleGoogle Scholar
- Andrianantoandro E, Basu S, Karig D, Weiss R: Synthetic biological: new engineering rules for an emerging discipline. Mol Syst Biol. 2006, 2: 2006.0028-View ArticleGoogle Scholar