Email updates

Keep up to date with the latest news and content from BMC Systems Biology and BioMed Central.

Open Access Highly Accessed Research article

Robust synchronization control scheme of a population of nonlinear stochastic synthetic genetic oscillators under intrinsic and extrinsic molecular noise via quorum sensing

Bor-Sen Chen* and Chih-Yuan Hsu

Author Affiliations

Lab of Control and Systems Biology, Department of Electrical Engineering, National Tsing Hua University, Hsinchu, 30013, Taiwan

For all author emails, please log on.

BMC Systems Biology 2012, 6:136  doi:10.1186/1752-0509-6-136


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1752-0509/6/136


Received:26 July 2012
Accepted:10 October 2012
Published:26 October 2012

© 2012 Chen and Hsu; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

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.

Results

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.

Conclusion

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.

Background

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 [22]. Recently, several mechanisms of intercell coupling of synthetic genetic oscillators have been discussed, to enhance the oscillating response of the synthetic biological system [1]. 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 [1]. 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 [41], 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 [42]. 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 [43], and a wide variety of regular and irregular rhythms induced by periodic stimulation of squid giant axons [21]. 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 [9]. 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 [41].

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.

Methods

Stochastic models of nonlinear stochastic synthetic genetic oscillators under intrinsic kinetic parameter fluctuations and external molecular noise

Model description

Before discussion of synchronization of more general synthetic genetic oscillators, we here provide a design example of coupled repressilators to illustrate the interesting phenomenon of synchronization in coupled dynamic cells. Then model description, definition, and theoretical results of robust synchronization of more general coupled synthetic genetic networks will be introduced for further study in the sequel. The repressilator is a network of three genes, the products of which inhibit the transcription of each other in a cyclical manner [1]. The gene lacI (from E. coli.) codes for the protein LacI, which inhibits the transcription of the gene tetR. The product of the latter, TetR, inhibits the transcription of gene cI (from λ phage); the protein product CI in turn inhibits the expression of lacI, thus completing the cycle. Garcia-Ojalvo et al. proposed a modular addition to the repressilator, with the aim of coupling a population of cells containing this network [1]. The basic mechanism of communication among the cells is based on quorum sensing, which was first discovered in the bioluminescent bacteria, Vibrio fisceri. These bacteria exhibit collective behaviors, mainly using two proteins (see Figure 1). The first protein, LuxI, synthesizes a small molecule known as an auto-inducer (AI), which diffuses freely through the cell membrane. The second protein, LuxR, binds to the AI molecule to form a complex, which subsequently activates transcription of various genes, as shown in Figure 1 (a detailed description of this molecular pathway is provided by Garcia-Ojalvo et al. [1]). Recently, several studies have indicated that quorum sensing can be engineered and used as an intercellular signaling system in E. coli[1,7,9,15-17].

thumbnailFigure 1. Synchronization scheme of N coupled genetic oscillators distributed in different host cells by the quorum sensing mechanism.

Based on the synchronization scheme of coupled repressilators via quorum sensing mechanism (Figure 1), the mRNA dynamics in the cell i = 1, 2, …, N are governed by repressible transcription of all three genes of the repressilator plus transcription activation of the additional copy of the lacI gene, and by mRNA degradation [7].

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M1">View MathML</a>

(1)

where xai , xbi and xci are the concentration of mRNA transcribed from tetR, cI and lacI in cell i, respectively; concentrations of the corresponding proteins are represented by xAi , xBi and xCi , respectively. The concentration of AI inside each cell is denoted by xSi. α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 [7].

The dynamics of proteins TetR, CI and LacI are given respectively as [7]

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M2">View MathML</a>

(2)

where parameters βA, βB and βC are the translation rates of the proteins from their mRNA, and γP represents the dimensionless degradation rate of proteins TetR, CI, and LacI in the cell. The intercellular protein AI in cell i is synthesized by the protein LuxI and diffuses through the cell wall where it undergoes degradation, leading to the following equation [7].

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M3">View MathML</a>

(3)

where ηs measures the diffusion rate of AI across the cell membrane, βs is the synthesis rate of AI, and γs gives the rate of decay of AI. Consequently, the whole coupled synthetic system is expressed by (1)-(3), which can be represented by the following more generalized nonlinear dynamic equation

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M4">View MathML</a>

(4)

where xi(t) = (xai(t), xbi(t), xci(t), xAi(t), xBi(t), xCi(t), xSi(t))T ∈ Rm is the state vector of the ith synthetic oscillator; f(·) : Rm → Rm is a smooth nonlinear function that characterizes the behavior of the synthetic oscillator; g(·) : Rm → Rm is a smooth nonlinear inner-coupling function; and C = (cij)N × N ∈ RN × N is the coupling configuration matrix, where cij>0 means that the ith synthetic oscillator is coupled with the jth synthetic oscillator directly, otherwise cij=0. Assume that the diagonal elements of C satisfy <a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M5">View MathML</a>.

Synthetic genetic oscillators under intrinsic parameter fluctuations

In general, the synthetic genetic oscillators suffer kinetic parameter fluctuations from transcription control, alternative splicing, translation, genetic mutation and diffusion to biological modification of transcription factors [17,19,31,34,44]. In this situation, dynamic equations of coupled synthetic genetic oscillators in (1)-(3) are modified as

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M6">View MathML</a>

(5)

where Δαj, j ∈ {a, b, c, S}, Δβj, j ∈ {A, B, C, s}, Δγj, j ∈ {m, p, s} and Δηs denote the amplitudes of the kinetic parameter fluctuations and ni(t) is random white noise with zero-mean and unit variance, i.e., Δαj, Δβj, Δγj and Δηs denote the deterministic part of the stochastic kinetic parameter fluctuations Δαjni(t), Δβjni(t), Δγjni(t) and Δηsni(t), respectively, and ni(t) absorbs the stochastic property of random intrinsic kinetic parameter fluctuations. The covariances of the stochastic parameter fluctuations Δαjni(t), Δβjni(t), Δγjni(t) and Δηsni(t) are given as cov(ξni(t), ξni(t)) = ξ2Δt,τ, with ξ ∈ {Δαj, Δβj, Δγj, Δηs}, respectively, where Δt,τ denotes the delta function, that is, Δt,τ = 1 if t = τ and Δt,τ = 0 if t ≠ τ, i.e., Δαj, Δβj, Δγj and Δηs denote the corresponding standard deviations of the stochastic parameter fluctuations Δαjni(t), Δβjni(t), Δγjni(t) and Δηsni(t), respectively. Thus, the nonlinear stochastic coupled synthetic genetic oscillators under intrinsic kinetic parameter fluctuations in the host cell i can be represented by

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M7">View MathML</a>

(6)

where <a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M8">View MathML</a> denotes the intrinsic kinetic parameter fluctuations of synthetic genetic oscillator in the host cell i. For the convenience of analysis and control design of the synchronization in synthetic genetic oscillators inserted into different host cells, the nonlinear stochastic equation in (6) can be represented by the following Ito’s stochastic differential equation

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M9">View MathML</a>

(7)

where wi(t) is a standard Wiener process or Brownian motion with dwi(t) = ni(t)dt to represent the random kinetic fluctuations of synthetic gene circuits.

Synthetic genetic oscillators under external environmental molecular noise

In general, a synthetic genetic oscillator in vivo also suffers from extrinsic environmental molecular noise, such as transmitted noise from upstream and global noise affecting all cells. Therefore, the nonlinear stochastic equation of coupled oscillators in (7) should be modified to mimic realistic dynamic behavior, as follows:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M10">View MathML</a>

(8)

where the signal vector vi(t) denotes extrinsic environmental molecular noise from the environment and hi 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

The nonlinear stochastic coupled synthetic genetic oscillators in (8) are said to reach synchronization asymptotically if x1(t) = x2(t) = ⋯ = xN(t) = s(t) in probability as t → tf or

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M11">View MathML</a>

(9)

in which s(t) ∈ Rm is the synchronization solution satisfying

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M12">View MathML</a>

(10)

Let us denote the synchronization error signal for synthetic genetic oscillators as

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M13">View MathML</a>

(11)

According to (8), the synchronization error dynamics for cell i are then described by

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M14">View MathML</a>

(12)

for i = 1, 2, …, N which can be augmented as

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M15">View MathML</a>

(13)

where

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M16">View MathML</a>

and

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M17">View MathML</a>

Suppose the influence of extrinsic environmental molecular noise on the synchronization error can be bounded by the following noise-filtering level ρ of coupled synthetic oscillators [30,31,44]

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M18">View MathML</a>

(14)

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.

Remark 1

(i) 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.

(ii) If the extrinsic environmental molecular noise v is deterministic, then the expectation on v(t) in (14) should be neglected. If the initial condition e(0) is considered, then the noise-filtering level in (14) should be modified as follows [15,16]:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M19">View MathML</a>

(15)

For some Lyapunov function V(e(0)), i.e., the energy due to the initial condition e(0) should be considered in the influence of noise on synchronization [30,45].

Results

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.

Proposition 1

If there exists a positive function V(e) > 0 with V(0) > 0 solving the following HJI,

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M20">View MathML</a>

(16)

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

Since ρ denotes an upper bound of the effect of v(t) on synchronization, the real effect can be obtained by minimizing ρ to as small a value as possible. Therefore, the noise-filtering ability of synchronized oscillators on v(t) can be obtained by solving the following constrained optimization:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M21">View MathML</a>

(17)

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.

If the noise-filtering ability ρ0 cannot satisfy the designer’s specification, in order to enhance the noise filtering of extrinsic noise, we need to specify the design parameters of nonlinear stochastic coupled synthetic genetic oscillators, for example, the kinetic parameters αa, αb, αc, βA, βB, βC, γm, and γp in (5) to solve the constrained optimization in (17) to enhance the noise-filtering ability, i.e.,

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M22">View MathML</a>

(18)

subject to HIJ in (16) with V(e) > 0

Before the discussion on the synchronization robustness criterion of coupled synthetic genetic oscillators, some definitions on synchronization robustness, intrinsic robustness and extrinsic robustness are given as follows:

(1) Synchronization robustness: The ability of coupled synthetic genetic oscillators to resist both intrinsic noise and extrinsic noise so that the synchronization can be maintained.

(2) Intrinsic robustness: The ability of coupled synthetic genetic oscillator to tolerate intrinsic parameter fluctuation to maintain synchronization.

(3) Extrinsic robustness: The filtering ability to attenuate the effect of environmental noise on the synchronization of coupled synthetic genetic network.

Remark 2

Substituting the noise-filtering ability ρ0 of (17) into (16) in Proposition 1, we get the following equivalent synchronization robustness criterion

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M23','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M23">View MathML</a>

(19)

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 kth 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]:

Rulek: If z1,i(t) is Fk1 and z2,i(t) is Fk2 … and zg,i(t) is Fkg, then

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M24','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M24">View MathML</a>

(20)

for k = 1, 2, ⋯, L, where zg,i is the element of premise variables of the ith coupled oscillation system, i.e., zi = z1,i, …, zg,iT; Fkg is the fuzzy set; Ak, Bk, AWk, and BWk are the fuzzy system matrices; L is the number of if-then rules; and g is the number of premise variables. The physical meaning of fuzzy rule k is that if the premise variables z1,i(t), z2,i(t), ⋯, zg,i(t) are with the fuzzy sets Fk1, Fk2, ⋯, Fkg, then the synchronization error dynamics in (12) can be represented by interpolating the linearized synchronization error dynamics in (20) via the fuzzy basis. The fuzzy synchronization error dynamics in (20) is referred as follows [32,40,46]:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M25','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M25">View MathML</a>

(21)

where <a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M26','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M26">View MathML</a>, Fkj(zj,i) is the grade of membership of zj,i(t) in Fkj or the possibility function of zj,i(t) in Fkj, and μk(zi) is called fuzzy basis function for k = 1, 2, …, L. The denominator or <a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M27','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M27">View MathML</a> in the above fuzzy basis function is only for normalization, so that the total sum of fuzzy basis is <a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M28','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M28">View MathML</a>. The physical meaning of (21) is that the fuzzy stochastic system interpolates L local linear stochastic systems through nonlinear basis μk(zi) to approximate the nonlinear stochastic system in (13). In this situation, the nonlinear stochastic coupled oscillation systems in (13) can be represented by the fuzzy interpolation system as follows:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M29','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M29">View MathML</a>

(22)

where μk(z) = diag(μk,1(z1), …, μk,N(zN)) and z = z1, …, zNT.

Remark 3

In [39], 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 Ak, Bk, AWk, and BWk 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.

Proposition 2

If there exists a positive definite symmetric matrix P>0 solving the following LMIs,

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M30','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M30">View MathML</a>

(23)

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

Therefore, the optimal noise-filtering design of synchronized oscillation systems obtained by solving the HJI-constrained optimization problem in (18) could be replaced by solving the following constrained optimizations, respectively

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M31','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M31">View MathML</a>

(24)

subject to P>0 and LMIs in (23)

Additional file 1 Supplementary Appendix. This file contains appendices A, B and C.

Format: PDF Size: 60KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Remark 4

(i)  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 Ak, 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).

(ii)  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), FW(x, s) and GW(x, s) by the fuzzy approximations in (22), in Proposition 2, V(e) = eTPe 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 ρ.

(iii)  In general, the constrained optimization problems in (24) are called eigenvalue problem [41], which can be efficiently solved by the MATLAB LMI toolbox.

(iv)  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

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M32','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M32">View MathML</a>

(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 IN ⊗ AWk + C ⊗ BWk 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

If robust synchronizations of coupled synthetic genetic oscillators cannot be achieved spontaneously via the parameter design in the above sections, then a control strategy is needed from external stimulation inputs to improve the robust synchronization of coupled synthetic genetic oscillators. External stimulation inputs are known to play an important role in the synchronization of biological rhythms. Recently, several methods of periodic stimulation for synchronization of nonlinear oscillators have been introduced [7-9,17]. However, even simple methods may show enormous complexity in the control scheme for synchronization of nonlinear stochastic coupled oscillators. In this study, based on nonlinear H stochastic control theory, an input control strategy is introduced to enhance the robust synchronization. If AI is injected into a common medium to increase the average concentration of AI protein in the exteracellular environment, which in turn increases the cellular communication of coupled oscillation systems, then the dynamics of the signaling molecule AI in the cellular environment, as shown in (3), should be modified as

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M33','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M33">View MathML</a>

(26)

where ue = Q represents an extracellular control input, which can be implemented via the injection of inducer AI.

For the simplicity of control design, suppose that the following control input ue = Q is employed to improve the robust synchronization of the nonlinear stochastic coupled synthetic oscillation systems. In this situation, the synchronization error dynamics in (13) should be modified as follows:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M34','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M34">View MathML</a>

(27)

where Ce(Q) = (ceij(Q))N × N ∈ RN × N is the coupling configuration matrix, in which ceii(Q) = − ηs(1 − N− 1)(Qe + Q) if i=j, otherwise ceij(Q) = ηsN− 1(Qe + 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.

Corollary 3

For the nonlinear stochastic coupled synthetic genetic oscillators with an extracellular control input ue=Q in the terms of Ce(Q) in (27), if there exists a positive solution V(e) > 0 with V(0) = 0 to the following HJI

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M35','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M35">View MathML</a>

(28)

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

The inequality (28) is equivalent to synchronization robustness criterion

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M36">View MathML</a>

(29)

The physical meaning of synchronization robustness criterion in (29) is that if we can specify control parameter Q to improve the synchronization robustness to provide more intrinsic robustness and more extrinsic robustness to tolerate more intrinsic parameter fluctuation and filter more extrinsic noise, then the robust synchronization of the nonlinear stochastic coupled synthetic genetic oscillators in (27) can be guaranteed. Similarly, the optimal noise-filtering design of synchronized oscillation systems by the extracellular control input in (27) can be achieved by solving the following constrained optimization problem:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M37','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M37">View MathML</a>

(30)

subject to V(e) > 0 and HJI in (28)

In general, it is still very difficult to specify the control parameter ue=Q to solve the HJI-constrained optimization in (30) for achieving the optimal noise filtering for synchronized synthetic genetic oscillators. Therefore, the fuzzy approximation method is again employed to simplify the control design procedure. Based on the fuzzy approximation method, the following fuzzy interpolation system is employed to approach the nonlinear stochastic coupled oscillation systems in (27):

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M38','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M38">View MathML</a>

(31)

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.

Corollary 4

For stochastic synchronized oscillation systems, if there exists a symmetric solution P > 0 to the following LMIs for a prescribed noise-filtering level ρ

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M39','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M39">View MathML</a>

(32)

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 IN ⊗ Ak + Ce(Q) ⊗ Bk 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.

Similarly, based on the fuzzy approximation method, an optimal noise-filtering design of synchronized oscillation systems by using the extracellular control input in (30) can be achieved by solving the following constrained optimization problem:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/6/136/mathml/M40','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/6/136/mathml/M40">View MathML</a>

(33)

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.

The design procedure of external inducer control for robust synchronization of the coupled network is summarized as follows.

(1) Consider a synthetic genetic network of N coupled oscillators with intrinsic kinetic parameter fluctuations and extrinsic environmental molecular noise.

(2) Given the prescribed disturbance attenuation level ρ.

(3) 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.

(4) 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 purpose of this example is to demonstrate the effectiveness of the theoretical synchronization result of synthetic gene oscillators in mimicking real biological oscillator systems. We consider a synthetic genetic network of N=10 coupled synthetic genetic oscillators with intrinsic kinetic parameter fluctuations and extrinsic environmental molecular noise in (5). The simulation results are shown in Figure 2. It can be seen that the parameter set of the coupled synthetic oscillator network, as listed below the figure, cannot make the whole network synchronize spontaneously. Suppose we want to specify a control parameter Q (which is proportional to the density of inducer AI) in (26) to compensate for the inefficiency of coupling between the synthetic genetic oscillators from the quasi-steady-state point of view. The design procedure first begins with representing the nonlinear stochastic synchronization error dynamic in (22) by the T-S fuzzy synchronization error dynamic model in (31), using the interpolation of several linear stochastic systems as presented in Additional file 1: Appendix C. According to the fuzzy approximation and Corollary 4, our control design problem is how to specify Q (i.e. the corresponding density of inducer AI), so that the ten coupled synthetic genetic oscillators have a positive solution P > 0 with a prescribed noise filtering level ρ = 0.56 to guarantee robust synchronization under intrinsic kinetic parameter fluctuations and extrinsic environmental molecular noise. The LMI toolbox in MATLAB can then be used to significantly simplify the system analysis and design procedure. With Q = 0.66 solved from LMIs in (32), the outputs of the coupled gene network of ten synthetic oscillators under intrinsic parametric fluctuations and extrinsic noise are shown in Figure 3. It can be seen that the coupled synthetic genetic oscillators have robust synchronizability to achieve the synchronous behavior despite the effect of uncertain initial state, intrinsic kinetic parameter fluctuations, and extrinsic environmental molecular noise on the host cell. According to a Monte Carlo simulation with 100 runs, the noise-filtering level of the coupled gene network is given by (E ∫ 0100eT(t)Re(t)dt)/(E ∫ 0100vT(t)v(t)dt) ≈ 0.192 < 0.562. It can be clearly seen that based on our proposed design method, the coupled gene network cannot only tolerate kinetic parameter variations but also attenuate the extrinsic molecular noise below a desired level to achieve a robust synchronization.

thumbnailFigure 2. Ten coupled genetic oscillators. The parameter values in (1), (2), and (3) are set as follows: αa = αb = αc = 216, αS = 20, μ = 1.2, μS = 1, n = 2, γS = 1, ηS = 2, βS = 0.1, βA = βB = βC = 1, γm = 6.9315, γp = 1.1552 and Qe = 0.09 [1]. Suppose the nonlinear stochastic coupled synthetic oscillators suffer from stochastic parameter fluctuations as shown in (8) with Δαa = Δαb = Δαc = 2.16, ΔαS = 0.2, ΔβA = ΔβB = ΔβC = 0.01, ΔβS = 0.001, ΔηS = 0.02, Δγm = 0.06, Δγp = 0.01, and ΔγS = 0.01. For the convenience of simulation, we assume that the extrinsic molecular noise v1~v10 is independent Gaussian white noise with a mean of zero and standard deviation of 0.02. It can be seen that coupled synthetic oscillators cannot achieve synchronization under these intrinsic kinetic parameter fluctuations and extrinsic molecular noise.

thumbnailFigure 3. The robust synchronization result of ten coupled synthetic oscillators in Figure2, by external control with Q=0.66. Based on a Monte Carlo simulation with 100 runs, the noise filtering level is given by (E ∫ 0100eT(t)Re(t)dt)/(E ∫ 0100vT(t)v(t)dt) ≈ 0.192 < 0.562.

Discussion

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.

Conclusions

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.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

BSC formulated the research topic and carried out some robust synchronization control designs. CYH participated in some control designs and performed simulation. Both authors read and approved the final manuscript.

Acknowledgements

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.

References

  1. 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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  2. Gardner TS, Cantor CR, Collins JJ: Construction of a genetic toggle switch in E. coli.

    Nature 2000, 403:339-342. PubMed Abstract | Publisher Full Text OpenURL

  3. Glass L: Synchronization and rhythmic processes in physiology.

    Nature 2001, 410:277-284. PubMed Abstract | Publisher Full Text OpenURL

  4. Gonze D, Bernard S, Waltermann C, Kramer A, Herzel H: Spontaneous synchronization of coupled circadian oscillators.

    Biophys J 2005, 89:120-129. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  5. 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. PubMed Abstract | Publisher Full Text OpenURL

  6. 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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  7. Wang RQ, Chen LN: Synchronizing genetic oscillators by signaling molecules.

    J Biol Rhythm 2005, 20:257-269. Publisher Full Text OpenURL

  8. Wang RQ, Jing ZJ, Chen LN: Modelling periodic oscillation in gene regulatory networks by cyclic feedback systems.

    Bull Math Biol 2005, 67:339-367. PubMed Abstract | Publisher Full Text OpenURL

  9. Wang RQ, Chen LN, Aihara K: Synchronizing a multicellular system by external input: an artificial control strategy.

    Bioinformatics 2006, 22:1775-1781. PubMed Abstract | Publisher Full Text OpenURL

  10. Chen LN, Wang RS, Zhang XS: Biomolecular Networks: Methods and Application in Systems Biology. John Wiley&Sons, New York; 2009. OpenURL

  11. Chen LN, Wang RS, Li C, Aihara K: Modeling Biomolecular Network in Cells: Structures and Dynamics. Springer, London; 2010. OpenURL

  12. Csermely P: Creative elements: network-based predictions of active centres in proteins and social networks.

    Trends Biochem Sci 2008, 33:569-576. PubMed Abstract | Publisher Full Text OpenURL

  13. Csermely P: Weak links: The universal key to the ability of networks and complex systems. Springer, Heidelberg; 2009. OpenURL

  14. 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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  15. Zhou TS, Chen LN, Aihara K: Molecular communication through stochastic synchronization induced by extracellular fluctuations.

    Phys Rev Lett 2005, 95:178103. PubMed Abstract | Publisher Full Text OpenURL

  16. Zhou TS, Zhang JJ, Yuan ZJ, Chen LN: Synchronization of genetic oscillators.

    Chaos 2008, 18:037126. PubMed Abstract | Publisher Full Text OpenURL

  17. Chen LN, Wang RQ, Zhou TS, Aihara K: Noise-induced cooperative behavior in a multicell system.

    Bioinformatics 2005, 21:2722-2729. PubMed Abstract | Publisher Full Text OpenURL

  18. Elowitz MB, Leibler S: A synthetic oscillatory network of transcriptional regulators.

    Nature 2000, 403:335-338. PubMed Abstract | Publisher Full Text OpenURL

  19. Elowitz MB, Levine AJ, Siggia ED, Swain PS: Stochastic gene expression in a single cell.

    Science 2002, 297:1183-1186. PubMed Abstract | Publisher Full Text OpenURL

  20. Hasty J, McMillen D, Collins JJ: Engineered gene circuits.

    Nature 2002, 420:224-230. PubMed Abstract | Publisher Full Text OpenURL

  21. 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. PubMed Abstract | Publisher Full Text OpenURL

  22. 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. PubMed Abstract | Publisher Full Text OpenURL

  23. You LC, Cox RS, Weiss R, Arnold FH: Programmed population control by cell-cell communication and regulated killing.

    Nature 2004, 428:868-871. PubMed Abstract | Publisher Full Text OpenURL

  24. Pai A, Tanouchi Y, Collins CH, You L: Engineering multicellular systems by cell–cell communication.

    Curr Opin Biotech 2009, 20:461-470. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. Gonze D: Coupling oscillations and switches in genetic networks.

    Biosystems 2010, 99:60-69. PubMed Abstract | Publisher Full Text OpenURL

  26. 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. OpenURL

  27. Lu JH, Chen GR: A time-varying complex dynamical network model and its controlled synchronization criteria.

    IEEE Trans Auto Control 2005, 50:841-846. OpenURL

  28. 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. Publisher Full Text OpenURL

  29. Lu WL, Chen TP: Synchronization of coupled connected neural networks with delays.

    IEEE Trans Circ Syst I 2004, 51:2491-2503. Publisher Full Text OpenURL

  30. Chen BS, Zhang WH: Stochastic H-2/H∞ control with state-dependent noise.

    IEEE Trans Auto Control 2004, 49:45-57. Publisher Full Text OpenURL

  31. Chen BS, Wang YC: On the attenuation and amplification of molecular noise in genetic regulatory networks.

    BMC Bioinforma 2006, 7:52. BioMed Central Full Text OpenURL

  32. Chen BS, Tseng CS, Uang HJ: Robustness design of nonlinear dynamic systems via fuzzy linear control.

    IEEE Trans Fuzzy Syst 1999, 7:571-585. Publisher Full Text OpenURL

  33. 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. PubMed Abstract | Publisher Full Text OpenURL

  34. Thattai M, van Oudenaarden A: Intrinsic noise in gene regulatory networks.

    Proc Natl Acad Sci USA 2001, 98:8614-8619. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  35. Blake WJ, Kærn M, Cantor CR, Collins JJ: Noise in enkaryotic gene expression.

    Nature 2003, 422:633-637. PubMed Abstract | Publisher Full Text OpenURL

  36. McAdams HH, Arkin A: It's a noisy business! Genetic regulation at the nanomolar scale.

    Trends Genet 1999, 15:65-69. PubMed Abstract | Publisher Full Text OpenURL

  37. Li C, Chen L, Aihara K: Synchronization of coupled nonidentical genetic oscillators.

    Phys Biol 2006, 3:37-44. PubMed Abstract | Publisher Full Text OpenURL

  38. Passino KM, Yurkovich S: Fuzzy control. Addison-Wesley, Menlo Park; 1998. OpenURL

  39. Takagi T, Sugeno M: Fuzzy identification of systems and its applications to modeling and control.

    IEEE Trans Syst Man Cyber 1985, 15:116-132. OpenURL

  40. Tseng CS, Chen BS: H∞ decentralized fuzzy model reference tracking control design for nonlinear interconnected systems.

    IEEE Trans Fuzzy Syst 2001, 9:795-809. Publisher Full Text OpenURL

  41. Boyd S, El Ghaoui L, Feron E, Balakrishnan V: Linear Matrix Inequalities in System and Control Theory. Society for Industrial & Applied Mathematics, Philadelphia; 1994. OpenURL

  42. Schultz TF, Kay SA: Circadian clocks in daily and seasonal control of development.

    Science 2003, 301:326-328. PubMed Abstract | Publisher Full Text OpenURL

  43. Wagemakers A, Buldu JM, Garcia-Ojalvo J, Sanjuan MAF: Synchronization of electronic genetic networks.

    Chaos 2006, 16:013127. PubMed Abstract | Publisher Full Text OpenURL

  44. Chen BS, Wu WS: Robust filtering circuit design for stochastic gene networks under intrinsic and extrinsic molecular noises.

    Math Biosciences 2008, 211:342-355. Publisher Full Text OpenURL

  45. Zhang W, Chen BS, Tseng CS: Robust H∞ filtering for nonlinear stochastic systems.

    IEEE Trans Signal Proc 2005, 53:589-598. OpenURL

  46. 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. Publisher Full Text OpenURL

  47. Endy D: Foundations for engineering biology.

    Nature 2005, 438:449-453. PubMed Abstract | Publisher Full Text OpenURL

  48. Andrianantoandro E, Basu S, Karig D, Weiss R: Synthetic biological: new engineering rules for an emerging discipline.

    Mol Syst Biol 2006, 2:2006.0028. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL