## Dynamics of Delay Differential Equations with Its Applications 2014

View this Special IssueResearch Article | Open Access

# Quasisynchronization in Quorum Sensing Systems with Parameter Mismatches

**Academic Editor:**Chuangxia Huang

#### Abstract

The paper investigates quasisynchronization in a communication system, which consists of cells communicating through quorum sensing. With the help of Lyapunov function method and Lur’e system approach, some sufficient conditions for quasisynchronization are presented, and a bound on the synchronization errors is derived. The obtained theoretical results show that the synchronization quality is influenced by two parameters detrimentally: the error bound depends almost linearly on the mismatches between cells and depends sensitively on the diffusion rates of the signals inward the cell membrane. Numerical experiments are carried out to verify the theoretical results.

#### 1. Introduction

In the past decades, increasing interest has been shown to study the dynamics of coupled oscillator networks, which describes many complex systems in the field of nature and science. Due to the couplings among these oscillators, different types of synchronization can be realized in such systems.

In particular, a coupled oscillator network is called to be quasisynchronized, or weak synchronized, if the synchronization errors will be in some neighborhood of zero but will not tend to zero eventually. In other words, quasisynchronization means that the dynamical trajectories of each cell are similar but different from each other [1, 2]. By Lyapunov function method and a differential inequality method, two coupled identical oscillators with parameter mismatches were quasisynchronized via periodically intermittent control in [3]. Similar results were also obtained in a discontinuous master-response system with parameter mismatches [4]. It has also been shown that two coupled delayed oscillators with parameter mismatches can be lag quasisynchronized, and the error level is estimated by applying a generalized Halanary inequality and matrix measure [5].

To the best of our knowledge, most of researches on synchronization focused on two oscillators coupled directly [6–8]. And there are very few researches focused on the system composed of several oscillators coupled indirectly. However, many biological systems coupled indirectly can also exhibit synchronization such as a global cellular response. For an example of such mechanisms achieving synchronization, the unicellular bacteria are highly coupled through chemical signaling molecules. This process, termed quorum sensing [9–11], allows bacterial populations to exchange intercellular signals with their neighbor cells to coordinate gene expression and integrates cells to realize synchronization effectively. That is to say, each bacterium is connected with all the other bacteria through a mean field and coordinated precisely. As the result, all the bacteria form a “microsociety,” behaving synchronously and exhibiting various collective dynamics. In the past decade, synchronization induced by the intercell signaling mechanism has been widely investigated. For instance, synchronization induced by quorum sensing has been studied in networks composed by genetic relaxation oscillators [12], limit-cycle oscillators [13], and synthetic gene oscillators [14]. Another research has shown that a noisy community of such genetic oscillators can self-synchronize in a robust way and lead to a substantially global rhythmicity [15].

Due to the biological diversity, there are usually some parameter mismatches between the coupled oscillators in biological systems. Therefore, complete synchronization is hard to be achieved. Instead, quasisynchronization, which implies a state of synchronization with an error level, is more common in the biological systems. Motivated by the complexity and similarity of biological organisms and the potential applications, the paper studies quasisynchronization in quorum sensing systems with parameter mismatches. Through a new method different from many previous researches [3–5], the bound on the synchronization errors is estimated with help of Lur’e system, linear matrix inequalities, and Lyapunov function [16–19]. Both the theoretical results and the numerical simulations indicate that the synchronization errors stay in a neighborhood of zero, increase roughly linearly with the mismatches between individual cells, and depend sensitively on the diffusion rate of the signals inward the cell membrane.

The rest of the paper is organized as follows. In Section 2, we give some sufficient conditions for quasisynchronization in quorum sensing systems with parameter mismatches. The bound of the synchronization errors is also estimated. In Section 3, numerical examples are carried out to verify the theoretical results.

#### 2. Main Results

##### 2.1. Quorum Sensing Systems with Mismatches

Quorum sensing is a cell concentration dependent phenomenon in bacteria and fungi, which is mediated by small, diffusible signaling molecules that accumulate in the extracellular environment [9]. In such a multicell system, the individual oscillator in each cell is a network with three genes, , , and , the products of which inhibit the transcription of each other in a cyclic way. The gene expresses protein , which inhibits transcription of the gene . The product of inhibits transcription of the gene , the protein product of which in turn inhibits expression of , completing the cycle. These bacteria exhibit cell-to-cell communication through a mechanism that makes use of two proteins, the first one of which (LuxI) synthesizes a small molecule known as an autoinducer (AI), which can diffuse freely through the cell membrane. The principle of the phenomenon is that when a single bacterium releases autoinducers (AIs) into the environment, their concentration is too low to be detected. However, when sufficient bacteria are present, AI concentrations reach a threshold level that allows the bacteria to sense a critical cell mass and to activate target genes [9]. When a second protein (LuxR) binds to this molecule, the resulting complex activates transcription of various genes, including some coding for light-producing enzymes. The scheme of the network is shown in Figure 1. For further details, one is referred to previous articles [9–11].

Before we carry out the dynamics model for the -cell system described by differential equations, we make the following declaration throughout the paper. Let be the -dimensional Euclidean space, the set of all real matrices, the transpose of a square matrix , and the usual norm of a vector or the usual spectral norm of a square matrix. The notation (<0) is used to define a real symmetric positive (negative) definite matrix. If not explicitly stated, matrices are assumed to have compatible dimensions.

Then the dynamics model for the -cell system is built as follows:
where ; , , and are the concentrations of* mRNA* transcribed from genes , , and in cell , respectively; , , and are the concentrations of the corresponding proteins, respectively; and are concentrations of AI inside each cell and in the environment, respectively. The parameters and are the dimensionless degradation rates of the chemical molecules in cell ; are the dimensionless transcription rates in the absence of repressor; is the maximal contribution to the gene transcription in the presence of saturating amounts of AI; are the translation rates of the proteins from the* mRNAs*; is the synthesis rate of AI; is the Hill coefficient; and measure the diffusion rate of AI inward and outward the cell membrane.

We suppose that all the parameters mentioned above are time-varying in the vicinity of certain constants, which indicates that all the cells are similar but different from each other. For convenience, we decompose these parameters into two parts: a constant part that determines the values of the parameters and a time-varying part representing the parameter mismatches. For example, let the parameter , , which implies that the parameter is time-varying around the constant with the parameter mismatches , .

Denoting , , and for convenience, then muiticell system (1) can be rewritten as follows: where ; , with elements equal to zero except , , , and ; with elements equal to zero except , , and ; with elements equal to zero except ; ; ; and The matrices , , and describe the mismatches of the parameters , , and . Note the dimension of individual cells for multicell system (1). In fact, system (2) can describe the general model of quorum sensing mechanism, where the vector components of vector functions and are increasing functions and describes which components are coupled with the environment. Thus, the results obtained in the paper are also valid for any general multicell system based on quorum sensing.

##### 2.2. Sufficient Conditions for Quasisynchronization

Since the cells in realistic organisms are similar but different from each other, one can suppose that the mismatch matrices , , and are bounded as follows: Noticing that the concentrations of chemical molecules are bounded and and are monotonic functions of , we can conclude that

Since there exist parameter mismatches between different cells, multicell system (1) cannot be completely synchronized. Instead, we present another type of synchronization, which is defined as follows.

*Definition 1. *The multicell system (1) is said to reach quasisynchronization or weak synchronization, if, for any initial condition , there exists a positive constant such that

The regulatory functions and are both monotonic increasing functions; there exist two diagonal matrices and such that the following sector conditions are satisfied:
where , , and . Notice that the Lur’e system consists of a linear system feedback interconnected with a static nonlinearity that satisfies a sector condition [20]; multicell system (2) can be regarded as a Lur’e system. Consequently, with help of Lur’e system method in control theory and Lyapunov direct method, we obtain the following sufficient conditions for quasisynchronization of system (2).

Theorem 2. *If there exist symmetric matrices and , , and a constant such that the symmetric matrix
**
where is the unit matrix, multicell system (1) is quasisynchronized. The bound on the synchronization errors can be estimated by , where
*

*Proof. *Define a Lyapunov function with respect to multicell system (2) of the following form:
where is a positive definite matrix. According to [21], the Lyapunov function is equivalent to the following form:
where denotes the synchronization errors.

Based on Lyapunov direct method, if the time derivative of along the trajectories of (2) is negative outside of a neighborhood of the origin , multicell system (2) will achieve quasisynchronization with the errors staying in the neighborhood. Calculating the time derivative of along (2),
We have
Noticing the sector conditions (7) of and , we have
where . Then, for any and and any constant , there holds
where is defined in (8),

If matrix inequalities (8) and hold, one obtains that
According to Lyapunov direct method, all the trajectories of cells in system (2) will go closer to each other when stay outside of the neighborhood of the origin . Therefore, multicell system (2) realizes quasisynchronization with synchronization errors staying in the neighborhood of the origin .

Theorem 2 implies that the estimation of the bound on synchronization error is influenced by two factors, the mismatches between cells and the parameter . On one hand, quasisynchronization can be realized if the mismatches are small and the synchronization errors oscillate in a certain neighborhood of the origin. On the other hand, the parameter is determined by the dynamics of the individual cells and reflects the ability of synchronization of the inherent dynamics. Furthermore, the larger the parameter is, the smaller the bound on the synchronization errors is.

In fact, the proof of the theorem shows that the mismatches that go against synchronization could be compensated by the linear function . Both the error dynamics caused by the mismatches and could be compensated by the linear function , which is determined by the diffusion rate of the signals inward the cell membrane . Consequently, the parameter measures the quasisynchronization ability of the inherent dynamics of the cells. Therefore, if the two factors mentioned above can be controlled, the synchronization errors can be controlled.

As a special case, if there are no parameter mismatches between different cells, which implies that all the cells are identical, then it is easy to conclude that multicell system (2) realizes complete synchronization. In such a case, multicell systems (1) and (2) could be rewritten as where the parameters , , , , and are defined in multicell system (2). Then we obtain the following corollary.

Corollary 3. *If there exist symmetric matrices and , , and a constant such that the symmetric matrix
**
where is the unit matrix, multicell system (18) realizes complete synchronization.*

#### 3. Numerical Simulations

In order to demonstrate the effectiveness of our theoretical analysis, we give numerical examples based on multicell system (1) consisting of 6 cells. Set the parameters as follows: Suppose that the initial mismatches , , and are taken randomly in the open intervals , , and , respectively. The evolutions of the synchronization errors are shown to converge to a neighborhood of the origin in Figure 2. As can be seen, the smaller the bound on mismatches between the oscillators is, the smaller the synchronization errors are. And complete synchronization can be realized if .

**(a)**

**(b)**

In order to verify the relationship between the synchronization errors and the mismatches between cells, the figure of the synchronization errors transition for increasing the mismatches is plotted in Figure 3. It is obvious that quasisynchronization is realized and the synchronization errors increase with the mismatches roughly linearly.

Numerical simulations are also carried out to verify the relationship between the synchronization errors and the diffusion rate of the signals inward the cell membrane in the region . Figure 4 shows that the error bounds depend sensitively on the diffusion rates , which play the role of the coupling strength. The increase of the coupling strength makes for the decrease of the synchronization errors, but if the parameter is too large, the dynamics of oscillation of multicell systems (1) breaks. Therefore, all the parameters should be taken appropriately to ensure the dynamics of oscillation.

From Figures 3 and 4, one can see that the error concentrations of gene products are much smaller than those of genes; specifically, the concentrations of biosignals AI remain quasisynchronized even when the mismatches are big. During the time course of achieving synchronization, the synchronous states of AI decrease the mismatches between cells through quorum sensing. And the strong synchronizability of AI makes it act as biosignals to synchronize other chemical molecules. On the other hand, the figures also imply that the error concentrations of genes decrease greatly when the mismatches are decreased. The genes’ strong sensitivity to regulations (the synchronous states of the concentrations of gene products, especially AI) is the reason why it can be led to synchronization by AI. The two aspects make up of the mechanisms of collective behavior caused by quorum sensing during transcription, translation, translocation, and signal-transduction.

#### 4. Conclusions

Many previous researches studied the collective behavior of biological systems by using complete synchronization. However, in many real biological systems, individual organisms are similar but different from each other, and their behaviors are not completely identical either. Therefore, it is meaningful to carry out researches on quasisynchronization instead of complete synchronization. Our results on quasisynchronization in multicell systems coupled by quorum sensing indicate that mismatches between cells can lead to quasisynchronization, and the synchronization errors depend heavily on the parameter mismatches. Theoretical analysis shows that the synchronization errors will decrease if the coupling strength increases. All the results agree well with numerical simulations and biological phenomena in practice.

#### Conflict of Interests

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

#### Acknowledgments

This research was supported by Zhejiang Provincial NSF of China (nos. LQ12A01003 and LQ12A01002), NNSF of China (nos. 61203155, and 11001069), Guangxi Key Laboratory of Trusted Software (no. kx201417), and the Science Foundation of Guangxi Province (no. 2013GXNSFAA019006).

#### References

- R. Femat and J. Alvarez-Ram{\'{\i}}rez, “Synchronization of a class of strictly different chaotic oscillators,”
*Physics Letters. A*, vol. 236, no. 4, pp. 307–313, 1997. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - J. A. K. Suykens, P. F. Curran, J. Vandewalle, and L. O. Chua, “Robust nonlinear ${H}_{\infty}$ synchronization of chaotic Lur’e systems,”
*IEEE Transactions on Circuits and Systems I. Fundamental Theory and Applications*, vol. 44, no. 10, pp. 891–904, 1997. View at: Publisher Site | Google Scholar | MathSciNet - W. Zhang, J. Huang, and P. Wei, “Weak synchronization of chaotic neural networks with parameter mismatch via periodically intermittent control,”
*Applied Mathematical Modelling*, vol. 35, no. 2, pp. 612–620, 2011. View at: Publisher Site | Google Scholar | MathSciNet - X. Liu, T. Chen, J. Cao, and W. Lu, “Dissipativity and quasi-synchronization for neural networks with discontinuous activations and parameter mismatches,”
*Neural Networks*, vol. 24, no. 10, pp. 1013–1021, 2011. View at: Publisher Site | Google Scholar | Zentralblatt MATH - W. He, F. Qian, Q. Han, and J. Cao, “Lag quasi-synchronization of coupled delayed systems with parameter mismatch,”
*IEEE Transactions on Circuits and Systems I: Regular Papers*, vol. 58, no. 6, pp. 1345–1357, 2011. View at: Publisher Site | Google Scholar | MathSciNet - J. Zhang, Z. Ma, and G. Zhang, “Cluster synchronization induced by one-node clusters in networks with asymmetric negative couplings,”
*Chaos*, vol. 23, no. 4, Article ID 043128, 2013. View at: Google Scholar - X. Yang, C. Huang, and J. Cao, “An LMI approach for exponential synchronization of switched stochastic competitive neural networks with mixed delays,”
*Neural Computing and Applications*, vol. 21, no. 8, pp. 2033–2047, 2012. View at: Publisher Site | Google Scholar - J. Zhang, Z. Ma, and J. Cao, “Full synchronization studied by a set of partitions connected together,”
*Abstract and Applied Analysis*, vol. 2013, Article ID 352826, 8 pages, 2013. View at: Publisher Site | Google Scholar | MathSciNet - M. B. Elowitz and S. Leibier, “A synthetic oscillatory network of transcriptional regulators,”
*Nature*, vol. 403, no. 6767, pp. 335–338, 2000. View at: Publisher Site | Google Scholar - R. Wang and L. Chen, “Synchronizing genetic oscillators by signaling molecules,”
*Journal of Biological Rhythms*, vol. 20, no. 3, pp. 257–269, 2005. View at: Publisher Site | Google Scholar - J. Zamora-Munt, C. Masoller, J. Garcia-Ojalvo, and R. Roy, “Crowd synchrony and quorum sensing in delay-coupled lasers,”
*Physical Review Letters*, vol. 105, no. 26, Article ID 264101, 2010. View at: Publisher Site | Google Scholar - A. F. Taylor, M. R. Tinsley, F. Wang, Z. Huang, and K. Showalter, “Dynamical quorum sensing and synchronization in large populations of chemical oscillators,”
*Science*, vol. 323, no. 5914, pp. 614–617, 2009. View at: Publisher Site | Google Scholar - D. J. Schwab, A. Baetica, and P. Mehta, “Dynamical quorum-sensing in oscillators coupled through an external medium,”
*Physica D: Nonlinear Phenomena*, vol. 241, no. 21, pp. 1782–1788, 2012. View at: Publisher Site | Google Scholar - D. McMillen, N. Kopell, J. Hasty, and J. J. Collins, “Synchronizing genetic relaxation oscillators by intercell signaling,”
*Proceedings of the National Academy of Sciences of the United States of America*, vol. 99, no. 2, pp. 679–684, 2002. View at: Publisher Site | Google Scholar - J. Garcia-Ojalvo, M. B. Elowitz, and S. H. Strogatz, “Modeling a synthetic multicellular clock: repressilators coupled by quorum sensing,”
*Proceedings of the National Academy of Sciences of the United States of America*, vol. 101, no. 30, pp. 10955–10960, 2004. View at: Publisher Site | Google Scholar | MathSciNet - S. P. Boyd, L. E. I. Ghaoui, E. Feron, and V. BalaKrishnan, “Linear matrix inequalities in system and control theory,”
*SIAM Studies in Applied Mathematics*, vol. 15, pp. 119–126, 1994. View at: Google Scholar - C. Huang, H. Kuang, X. Chen et al., “An LMI approach for dynamics of switched cellular neural networks with mixed delays,”
*Abstract and Applied Analysis*, vol. 2013, Article ID 870486, 8 pages, 2013. View at: Publisher Site | Google Scholar | MathSciNet - C. Li, L. Chen, and K. Aihara, “Synchronization of coupled nonidentical genetic oscillators,”
*Physical Biology*, vol. 3, no. 1, pp. 37–44, 2006. View at: Publisher Site | Google Scholar - C. Huang and J. Cao, “Convergence dynamics of stochastic Cohen-Grossberg neural networks with unbounded distributed delays,”
*IEEE Transactions on Neural Networks*, vol. 22, no. 4, pp. 561–572, 2011. View at: Publisher Site | Google Scholar - M. Vidyasagar,
*Nonlinear Systems Analysis*, Prentice-Hall, Englewood Cliffs, NJ, USA, 2nd edition, 1993. - C. W. Wu,
*Synchronization in Coupled Chaotic Circuits and Systems*, World Scientific, Singapore, 2002. View at: Publisher Site | MathSciNet

#### Copyright

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