Email updates

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

Open Access Research article

Noise regulation by quorum sensing in low mRNA copy number systems

Marc Weber and Javier Buceta*

Author Affiliations

Computer Simulation and Modelling (Co.S.Mo.) Lab, Parc Científic de Barcelona, C/Baldiri Reixac 10-12, Barcelona 08028, Spain

For all author emails, please log on.

BMC Systems Biology 2011, 5:11  doi:10.1186/1752-0509-5-11

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


Received:28 July 2010
Accepted:20 January 2011
Published:20 January 2011

© 2011 Weber and Buceta; 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

Cells must face the ubiquitous presence of noise at the level of signaling molecules. The latter constitutes a major challenge for the regulation of cellular functions including communication processes. In the context of prokaryotic communication, the so-called quorum sensing (QS) mechanism relies on small diffusive molecules that are produced and detected by cells. This poses the intriguing question of how bacteria cope with the fluctuations for setting up a reliable information exchange.

Results

We present a stochastic model of gene expression that accounts for the main biochemical processes that describe the QS mechanism close to its activation threshold. Within that framework we study, both numerically and analytically, the role that diffusion plays in the regulation of the dynamics and the fluctuations of signaling molecules. In addition, we unveil the contribution of different sources of noise, intrinsic and transcriptional, in the QS mechanism.

Conclusions

The interplay between noisy sources and the communication process produces a repertoire of dynamics that depends on the diffusion rate. Importantly, the total noise shows a non-monotonic behavior as a function of the diffusion rate. QS systems seems to avoid values of the diffusion that maximize the total noise. These results point towards the direction that bacteria have adapted their communication mechanisms in order to improve the signal-to-noise ratio.

Background

Gene regulation at the transcriptional level is one of the corner stones of molecular and cellular biology [1]. Recent studies in prokaryotes have revealed the existence of antisense and alternative transcripts and multiple regulators per gene that imply a highly dynamic transcriptome more similar to that of eukaryotes than first thought [2]. Still, prokaryotic gene regulation mainly relies on the binding of regulatory proteins that attach to DNA for either stimulating or repressing transcription. These binding/unbinding events are intrinsically probabilistic because of the significance of thermal fluctuations at that scale and the low number of molecules involved in the process. In this regard, over the past years a growing number of experiments have indeed characterized not only the levels of randomness in cellular biochemical processes but also their functionality [3-8].

Technical advances such as the use of fluorescent tags in single-cell experiments have allowed for quantitative measurements of the noise in protein concentration and have shed light on the mechanisms of gene expression that lead to cell-to-cell variability [9-11]. Moreover, the advent of experimental approaches that permit to count individual mRNA and protein molecules in single cells has further evidenced the role played by fluctuations and their characteristics [12,13]. Thus, in E. coli, the direct measurement of integer-valued numbers of mRNA as a function of time has revealed transcriptional bursts with Poissonian statistics [14]. The latter is in agreement with the two-state gene expression model where switching between the active and inactive transcriptional regimes occurs with constant probability [15]. It is worth noticing that these noisy sources, far for being a nuisance, have been recognized to play a constructive role in many gene regulatory processes. Examples in this direction include the efficiency of the phage lambda switch [16] or the differentiation into the competence state in bacteria [17]. All in all, it is now accepted that stochastic and non-linear approaches are required for understanding the randomness in the dynamics of biochemical reactions and the effects of fluctuations in gene regulatory networks [7,8].

While a lot of modeling studies have focused on the single cell level [18,19], few of them have addressed the role played by noise at the colony level [20-23]. Our recent contributions within this topic, that illustrate the constructive role of stochasticity, include the noise-induced coherence resonance phenomenon in multicellular circadian clocks [24] and the interplay between the stochasticity of the cell cycle duration and the protein expression noise in bacterial colonies [25]. One relevant question within this context is how cellular populations deal with stochasticity in communication processes. In particular, we focus on the simplest cell-to-cell communication mechanism in prokaryotes: the so-called quorum sensing (QS). The term QS generically refers to the mechanism that allow bacteria to count their number, i.e. the colony size, by producing, exporting/importing into/from the environment, and detecting a diffusive signaling molecule, namely, the autoinducer [26,27]. As the population of bacteria grows, the autoinducer accumulates both in the extracellular medium and inside the cells. When the autoinducer concentration surpasses a threshold the expression of QS-controlled genes starts. This mechanism ultimately results in a response of the colony in a cell density dependent manner. Importantly, QS has opened the door to the design of gene circuits using synthetic biology approaches that control cell populations at the collective level [28-31].

Recent studies have shown that diffusion in QS reduces the noise at the level of the autoinducer [23]. However, to the best of our knowledge, the role played by different sources of stochasticity and their contribution to the dynamics of the signaling molecule has not been characterized yet. Moreover, while in eukaryotes the diffusion seems to contribute for enhancing the precision of regulatory processes, similar effects have not been reported in the context of QS [32]. We point out that a deep understanding of these issues is key in order to design robust synthetic circuits based on such bacterial communication system. Herein, we address these problems by studying the interplay between the QS communication and the transcriptional noise in bacterial populations. First, we aim at understanding how that interaction determines the dynamics of the autoinducer. Second, we aim at shedding light on the mechanisms that confer robustness to noise in QS communication. Transcriptional noise is expected to be particularly relevant when transcription events are short and rare. Under these conditions two main sources of stochasticity naturally arise: the dichotomous fluctuating dynamics of mRNA and the intrinsic noise due to low copy number of species. Thus, we restrict ourselves to the study of the aforementioned problems near the QS activation threshold where we can assume that the transcription events produce basal constitutive levels of mRNA as low as one molecule per cell at a time. We note that such mRNA production processes have been experimentally validated in prokaryotes revealing that the statistics of proteins bursts originates from the translation of a single mRNA molecule [33].

Our main findings are threefold. First, we show how the diffusion process leads to a repertoire of dynamics in regards of the signaling molecule. Second, we demonstrate that, for a large range of diffusion rate values, the main contribution to the total noise of the autoinducer concentration is the mRNA fluctuations. Finally, we show that the total noise exhibits a non-monotomic behavior as a function of the diffusion rate in contrast to previous results [23].

The paper is organized as follows. In the Methods section we introduce our modeling approach, the analytical calculations, and the parameter values used in our in silico experiments. In the next section, Results, we present our findings and compare the results from stochastic simulations with those from analytical calculations. Finally, we discuss further implications of our results in the Discussion, and summarize our findings in the Conclusion.

Methods

Modeling Approach

A large class of gram-negative bacteria use acyl homoserine lactones (AHLs) as signaling molecules [27,34]. These autoinducer molecules are typically synthesized by enzymes of the LuxI family and can freely diffuse across the cell membrane, i.e. by means of passive diffusion. When the concentration of autoinducer surpasses a critical threshold, it binds to its receptor (a cytoplasmic protein of the LuxR family) which then activates the expression of target genes, e.g. in Vibrio fisheri species the luxICDABE operon, that is responsible for the expression of luxI and luciferase. In contrast, when the autoinducer concentration is below the activation threshold, the transcription of the luxI gene occurs at a low basal rate, thus producing low levels of the enzyme. In this regime, the feedback regulation of the luxI gene leading to autoinduction can be disregarded. As a matter of fact, a number of QS systems lack autoinduction [35]. Consequently, to neglect the feedback in those cases is a valid approximation even above the activation threshold. Herein, for the sake of simplicity, we focus on this situation and describe the dynamics of the autoinducer near, but below, the activation threshold when feedback loops can be neglected and the downstream QS genes are not activated. We stress that, generally speaking, our results are not applicable above the activation threshold since the dynamics of the autoinducer is obviously conditioned by feedback effects.

Following previous approaches we consider a two-stage model for gene expression/regulation [15,36-38]. Thus, we assume that during the transcription events a single mRNA molecule, e.g. a luxI transcript, is produced and its dynamics can be then described by means of a Markovian dichotomous process [39],

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

(1)

where <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M2">View MathML</a> stands for the number of mRNA molecules at cell i and α and β for the transition rates between these states; i.e. α and β account for the probabilities per unit of time of mRNA degradation and transcription frequency respectively. Notice that the stochastic alternation of the mRNA between the values 0 (no mRNA) and 1 (a mRNA molecule) is not memoryless, i.e. white. Once a mRNA molecule is produced, and until it becomes degraded, the cell keeps noticing its presence and keep producing the autoinducer. That is, the transcriptional noise is a colored noise, and its autocorrelation decays exponentially with a characteristic time scale τc = (α + β)-1 [39].

Once a mRNA molecule is produced the translational, and post-translational processes (if any), leads to the appearance of functional LuxI synthetases. Yet, our interest here focuses on the dynamics of the signaling molecule. It has been shown that the amount of the synthetase substrate is not a limiting factor for the production of the autoinducer [40,41]. As a consequence, the levels of the signaling molecule depends directly on the expression levels of the synthetase. Ignoring intermediate biochemical steps in the autoinducer synthesis reduces the number of noise sources and may even change, under some circumstances, the observed dynamics [42]. Still, it is a valid approximation in many situations, and here we assume that the translation of the synthetase and the subsequent synthesis of the autoinducer, A, can be effectively described by a single chemical step with rate k+. In addition, we consider that the autoinducer becomes degraded at a rate k-, that is,

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

(2)

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

(3)

Passive diffusion of the autoinducer can be implemented by considering a new species, Aext, that accounts for the number of signaling molecules in the extracellular medium such that,

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

(4)

where D stands for the diffusion rate and r = V/Vext represents the ratio of the volume of a cell to the total extracellular volume. We consider all cells to have the same value of the diffusion rate. In addition, we assume a well-stirred system where spatial effects can be neglected.

As the bacterial population grows the autoinducer accumulates in the media. In experiments, in order to keep the concentration of the autoinducer below the activation threshold, such growth is compensated by means of a dilution protocol. As detailed below (see Parameters) the latter constitutes the main source of effective degradation of the signaling molecule. Thus, hereinafter we assume the degradation rate of the signaling molecule to be the same inside and outside the cell,

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

(5)

Figure 1 schematically represents the biochemical processes considered in our approach. The set of reactions (1)-(5) characterizes the stochastic dynamics of the autoinducer and that of the mRNA. Their probabilistic description is given by the corresponding master equation that is exactly sampled by means of the Gillespie algorithm in a N-cells system [43].

thumbnailFigure 1. Scheme of a simplified biochemical network of QS systems near the activation threshold. Schematic representation of the biochemical processes considered in our approach for describing the dynamics of the signaling molecule, A, in cell i. The mRNA dynamics satisfies a dichotomous process characterized by the states M0,1 corresponding to zero and one molecules respectively. Once the autoinducer has been produced, it can diffuse in and out the cell leading to cell communication (see text).

Analytical Calculations: Null Intrinsic Noise Approximation

Further insight into the dynamics of the signaling molecule can be obtained by analytical means as follows. Two stochastic contributions drive the dynamics of A: the mRNA fluctuations due to the random switching (mRNA present or not) and the molecular, i.e. intrinsic, noise due to low copy number of the autoinducer. As for the latter, it can be neglected if over the course of time Ai/(Ai + 1) ≃ 1 (large number of autoinducer molecules). While in our system this approximation is not totally justified (see parameters values below), it is useful to implement it in order to discriminate between the effects caused by different stochastic contributions and to obtain analytical expressions. In this case, it is straightforward to demonstrate that the dynamics of the autoinducer, Eqs. (1)-(5), can be described by the following coupled stochastic equations,

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

(6)

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

(7)

where <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M9">View MathML</a>, <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M10">View MathML</a> and <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M11">View MathML</a> stand for the concentration of species A and <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M12">View MathML</a> at cell i and for species Aext at the extracellular medium respectively, N is the colony size (number of cells), and 〈·〉 represents the population average. In Eq. (6) the term <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M13">View MathML</a> accounts for a dichotomous stochastic process characterized by the rates and states (α, β) and (0, V −1) respectively, and describes the fluctuating mRNA dynamics. We point out that in case that D = 0, Eq. (6) has been proposed to study graded and binary responses in stochastic gene expression. Interestingly, it has been shown that despite its simplicity it can actually reproduce some gene expression phenomena [37,44].

We can further proceed with the analytical calculations by implementing, as in previous studies e.g. [45], a quasisteady approximation for the dynamics of the external autoinducer, i.e. <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M14">View MathML</a> so that,

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

(8)

By substituting Eq. (8) into Eq. (6) we obtain an equation for the concentration of the signaling molecule inside a given cell that depends on the average 〈cA〉(the index i has been dropped),

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

(9)

In the absence of diffusion, Eq. (9) reveals that the concentration of the signaling molecule reaches a maximum value of <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M17">View MathML</a> when <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M18">View MathML</a>. In terms of <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M19">View MathML</a> and the time scale tc = 1/k-, the typical lifetime of a signaling molecule, the dimensionless version of Eq. (9) reads

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

(10)

where <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M21">View MathML</a> and <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M22">View MathML</a> being a Markovian dichotomous noise characterized by the states <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M23','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M23">View MathML</a> and the rates <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M24','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M24">View MathML</a> and <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M25','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M25">View MathML</a> Equation (10) can be formally closed by invoking the following self-consistency condition:

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

(11)

<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M27','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M27">View MathML</a> being the probability density solving Eq. (10) and <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M28','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M28">View MathML</a> its support [39]:

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

(12)

with

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

(13)

being the normalization constant. The condition (11) can be exactly solved and leads to the following value for the average concentration:

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

(14)

where <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M32','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M32">View MathML</a> is the average concentration of signaling molecules in the absence of diffusion. Thus, as expected, <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M33','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M33">View MathML</a>. For the sake of concision, on what follows we drop the argument term <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M34','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M34">View MathML</a> from the notation of <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M27','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M27">View MathML</a>. Note that <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M35','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M35">View MathML</a> has two states (barriers) that define its support. That is, the minimum and maximum values that the concentration of the autoinducer can reach as a function of the diffusion are:

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

(15)

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

(16)

It is easy to prove that the probability density <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M35','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M35">View MathML</a> shows a single extremum if,

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

(17)

Where the extremum is a maximum if <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M39','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M39">View MathML</a> and a minimum if <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M40','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M40">View MathML</a>. In the other cases the probability density does not display any extrema. Therefore, as a function of <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M41','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M41">View MathML</a> and <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M42">View MathML</a>, the probability density <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M35','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M35">View MathML</a> may show four different behaviours depending on the value of the diffusion coefficient as schematically represented in Figure 2A. However, a constraint in our modeling restricts the regions, i.e. behaviors, accessible to the autoinducer dynamics. We have assumed a low constitutive expression such that only a single mRNA molecule can be transcribed at a time. The latter implies that <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M43','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M43">View MathML</a>

thumbnailFigure 2. Probability densities of the signaling molecule and parameter space. Panel A: Schematic representation of the different probability densities of the autoinducer concentration depending on the value of <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M41','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M41">View MathML</a> and <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M42">View MathML</a> with respect to that of <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M44">View MathML</a>. Given a set of values (<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M41','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M41">View MathML</a>, <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M42">View MathML</a>) the dynamics of the autoinducer shows different behaviors depending on the value of the diffusion parameter since the transitions lines are located at <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M41','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M41">View MathML</a>, <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M42">View MathML</a> = 1 + <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M44">View MathML</a>. The constraints of our modelling in terms of the parameters values make the region on the top-left corner non-accessible (see text). Panel B: Parameter space diagram (<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M41','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M41">View MathML</a>, <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M42">View MathML</a>) indicating the sets of parameters used in the simulations (solid squares): γ1 = (8, 2), γ2 = (15, 5), γ3 = (8, 0.5), γ4 = (15, 0.5). The experimental values reported for the degradation rate of the mRNA leads to a biological meaningful range for <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M41','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M41">View MathML</a> (blue region). The low constitutive expression assumption is prescribed by the constraint <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M45','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M45">View MathML</a> (red colored region).

(the degradation rate of the mRNA is larger than the transcription rate) in order to assure that a maximum of one mRNA molecule is present in a cell at a given time. As a consequence, and independently of the diffusion value, the dynamics leading to the probability density shown at the top-left region of Figure 2A cannot be considered as physical in the context of our modeling approach. We stress that this constraint is not a fundamental ingredient for obtaining our results (see Discussion).

Finally, the noise of the autoinducer concentration reads,

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

(18)

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

Nonetheless we have disregarded the intrinsic noise in the analytical calculations, Eq. (18) will allow us to elucidate the contributions of different sources of noise as follows. By means of the numerical simulations (Gillespie) of the set of reactions (1)-(5) in a N-cell system, we can evaluate the total (intrinsic+transcriptional) noise. Hence, by substracting from that quantity the contribution of the transcriptional noise, i.e. Eq. (18), we obtain the levels of intrinsic noise (see Results).

Parameters

We are interested in the role played by the fluctuations of the signaling molecule, A, when its concentration is close to the activation of the QS switch. Therefore, we fix the mean concentration of the autoinducer and modulate the rest of the parameters in order to keep constant this value. Pai and You [46] have recently studied the core architecture of the QS mechanism for a comprehensive set of systems. These authors have estimated that the critical concentration of autoinducer needed for the activation of the QS genes ranges from 10 to 50 nM for most of the bacterial species. In our model, we set the average concentration of A to a typical value of <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M48','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M48">View MathML</a> = 25 nM. Yet, our results do not depend on the particular value we choose within that range. As shown below, see Results, this value fixes the level of intrinsic noise of the system. However, the interplay between diffusion and transcriptional noise does not depend on that. Moreover, by defining the so-called sensing potential, ν = (rN)-1 Pai and You estimated the range of critical cell densities for the QS activation. They concluded that its characteristic value is ν ~ 103 - 104. In our simulations we set this parameter to ν - 103. In the experimental setups the cells are typically present in a volume of a few milliliters and the total number of cells is of the order of 108 - 109. Therefore, the concentration of autoinducer in the medium is determined by the exchange of signaling molecules coming from many cells. In contrast, the behavior of QS systems with a very low number of cells can be significantly different, as shown by microfluidic confinement of cells in picoliter droplets [47-49]. In our study, in order to discard small system size effects, we choose a sufficiently large number of cells in the numerical simulations, N = 102. Since the typical volume of an E. coli cell is V = 1.5 μm3 then Vext = 105V (i.e. r = 10−5). We point out that keeping ν to a constant value necessarily requires an external dilution protocol for maintaining constant the cell density. In experiments, the control of the dilution rate is usually achieved by the use of chemostats or microfluidic devices [50]. The rate of dilution should compensate for the cell growth, ~ 2 · 10−2 min−1 (i.e. cell cycle duration ~50 min). In our modeling, by keeping constant the number of cells and the average concentration of the autoinducer, we tacitly assume a dilution protocol too. Importantly, the dilution rate effectively modifies the degradation rate of the signaling molecule. In this regard, while some bacteria species have hydrolytic enzymes that degrade AHLs, generally speaking, bacteria that synthetize AHLs do not degrade them enzymatically. In fact, AHLs are chemically stable species in aqueous solutions [51]. The degradation rate of the homoserine lactone 3-Oxo-C6-AHL has been measured in vitro revealing that this autoinducer is rather stable: ~ 3 · 10−4 min−1 [51]. Measurements of the degradation rate of other AHL autoinducers show similar results. Based on experimental data and mathematical modeling, the degradation rate of the signaling molecule in vivo has been also estimated [46]. Depending on the pH of the medium, the latter ranges from ~ 5 · 10−3 min−1 to ~ 2 · 10−2 min−1. Consequently, the dilution process constitutes the main source of effective degradation of A, both inside and outside the cell, and here we set k= 2 · 10−2 min−1.

As for the value of the diffusion rate, the coefficient of passive diffusion has been estimated for the 3-Oxo-C6-AHL autoinducer based on the measure of the diffusion of glucose and lactose through the outer membrane of E. coli [23]. For a typical cell volume of 1.5 μm3 the estimated coefficient of diffusion is ~ 103 min−1. Under these conditions the typical value for the dimensionless parameter <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M44">View MathML</a> is of the order of 104. Yet, active transport mechanisms for the autoinducer leads to much smaller effective diffusion values (see Discussion) and we explore the role played by this parameter.

In regards to the mRNA dynamics, α, the degradation rate, depends on the cell degradative machinery. To this respect, the half-lives of all mRNAs of Staphylococcus aureus have been recently measured during the mid-exponential phase. Most of the transcripts, 90%, have half-lifes shorter than 5 minutes [52,53].

According to these studies we restrict the mRNA degradation rate to the range ln(2)/5 min−1 < α < ln(2)/2 min−1. Consequently, <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M41','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M41">View MathML</a> > 1. As for the frequency of the transcription events, β is determined by particular characteristics of the gene regulatory process under consideration, e.g. the affinity of the regulatory proteins to the operator site and the initiation rate of transcription. Due to the assumption of a low constitutive expression, we choose values of parameter β satisfying the relation α > β. In particular, we implement the more restrictive condition α > 2β. Figure 2B recapitulate these constraints and show the different sets of <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M41','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M41">View MathML</a> and <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M42">View MathML</a> values that we have used in our simulations and analytical calculations.

Summarizing, N, r, and kare kept fixed in our simulations and analytical calculations and we explore the parameter space α, β, and D within the ranges, and satisfying the constraints, mentioned above. In every particular situation, once a set of those parameters is prescribed, we set the value of k+ by using Eq. (14) in order to keep the average value of cA near its critical concentration value (25 nM).

Results

Comprehensive Study of the Autoinducer Dynamics as a Function of the Diffusion Rate

The distribution of cA at the steady-state is computed for the different parameter sets according to the ranges and constraints described above (section Methods). In order to explore the role of the diffusion in the dynamics of the signaling molecule we first study the case <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M49','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M49">View MathML</a>. According to the analytical calculations, see Eq. (17), in this case two possible distributions for the concentration of cA can be observed depending on the value of <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M42">View MathML</a>. Since <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M41','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M41">View MathML</a> > 1 we can expect a maximum only if <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M42">View MathML</a> > 1 (note that <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M50','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M50">View MathML</a> if <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M49','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M49">View MathML</a>), otherwise extrema are not expected. The results of the numerical simulations (Gillespie), Figure 3A, reveal that scenario. Note that in all cases the histogram obtained from the simulations fits fairly well to the expression (12) except for deviations due to the intrinsic noise that are not taken into account by the analytical approach. The differences among dynamics are evidenced by the trajectories, Figure 3B. Thus, for <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M42">View MathML</a> < 1 the dynamics of the autoinducer shows a burst-like behavior. If <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M42">View MathML</a> > 1 the frequency of bursts is high enough to maintain the concentration of signaling molecules near the average and a single-peak distribution develops.

thumbnailFigure 3. Distributions and dynamics of the signaling molecule in a diffusionless system. Panel A: Distributions of cA at steady-state for different sets of parameters (<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M41','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M41">View MathML</a>, <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M42">View MathML</a>) as indicated in Figure 2B. In all cases <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M44">View MathML</a> = 0. The histogram obtained in simulations (blue bars) compares well with the distribution from the analytical calculations (blue line). Yet, deviations are observed due to intrinsic noise (see text). Panel B: The dynamics of the autoinducer show different behaviors depending on the region of the parameters phase space (see Figure 2). Two typical trajectories are shown with a grey-shaded background indicating the presence of a mRNA molecule in the cell.

If <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M44">View MathML</a> > 0 we expect a more fruitful phenomenology since the transition lines between behaviors in the parameter space (<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M41','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M41">View MathML</a>, <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M42">View MathML</a>) shift as a function of the diffusion (see Figure 2A). According to the analytical calculations we can anticipate that, for a given parameter set and as <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M44">View MathML</a> increases, the system explores different dynamical regimes. By taking as a reference the case γ2, that is (<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M41','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M41">View MathML</a>, <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M42">View MathML</a>) = (15, 5), Figure 4 shows the effect of the diffusion on the distribution (left column) and dynamics (center column) of cA in a given cell. The system initially displays a single-peak distribution for <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M44">View MathML</a> = 1. By increasing the diffusion coefficient we observe transitions to the other phases (monotonically decreasing and double-peak distributions). The corresponding dynamics of cA (right panels) show how the diffusion, acting as an additional effective degradation on A, first increases the sharpness of the bursts of production. For <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M44">View MathML</a> = 10, the diffusion is large enough to remove signaling molecules between consecutive burst events, thus leading to a monotonically decreasing distribution. Increasing the diffusion rate to <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M44">View MathML</a> = 100 leads to the situation where both <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M41','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M41">View MathML</a> and <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M42">View MathML</a> becomes smaller than 1 + <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M44">View MathML</a> and a bistable dynamics develops. Under these circumstances the concentration of autoinducer alternates between two states that correspond to a low concentration, when there is no mRNA production, and a high concentration, following the mRNA synthesis. As the diffusion further increases, e.g. <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M44">View MathML</a> = 2 · 103, the autoinducer molecules diffusing from the external medium into the cell set a constitutive level of this species. The latter explains the presence of A molecules in the cell even if no mRNA is produced. Finally, at very large values of <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M44">View MathML</a>, e.g. <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M44">View MathML</a> = 5 · 104, the low constitutive concentration of the autoinducer increases due to the influx of molecules when no mRNA is present whereas the concentration of A that is internally produced decreases due to the efflux of molecules. In this case, the whole N-cells system can be considered as a single volume with no diffusive barriers between cells. Thus, the burst events average out and, as a consequence, a single effective peak appears.

thumbnailFigure 4. Distributions and dynamics of the signaling molecule in a system with diffusion. Distributions (left column) and dynamics (center column) of cA at steady-state for different values of <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M44">View MathML</a>. The right most column stands for a density plot of the distribution of <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M51','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M51">View MathML</a> as a function of <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M52','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M52">View MathML</a> for discerning a putative increase in the molecular noise (see text). In all cases the parameters set (<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M41','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M41">View MathML</a>, <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M42">View MathML</a>) is γ2 (see Figure 2B). The production rate <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M53','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M53">View MathML</a> is modulated as a function of (<a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M41','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M41">View MathML</a>, <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M42">View MathML</a>, <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M44">View MathML</a>) in order to maintain constant the average 〈cA〉 = 25 nM. The histograms obtained in the stochastic simulations (blue bars, left column) are in qualitative agreement with the probability densities from the analytical calculations (blue line, left column). When increasing the diffusion coefficient, the system explores different dynamics as revealed by the trajectories shown in the center column. The grey-shaded background shown in the trajectories of cA indicates the presence of a mRNA molecule in the cell. The density plots (right column) reveals that the diffusion does not contribute to an increase of the intrinsic noise since the spreading of the distributions in a direction perpendicular to the diagonal does not grow.

Diffusion and Intrinsic Noise

Figure 4 shows that the theoretical distribution captures the essential features of the dynamics obtained in numerical simulations (Gillespie). The noticeable deviations are due to the intrinsic noise that are not considered in the theoretical analysis. Notice that as the diffusion increases those deviations seem to be larger. We stress that in our simulations we keep constant the average concentration of the autoinducer by modulating the effective production rate (see Eq. (14)). That is, as <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M44">View MathML</a> becomes larger, we increase the production rate k+ so that the average number of autoinducer molecules per cell remains the same. Consequently, the deviations between simulations and the theoretical analysis cannot be ascribed to a putative intracellular decrease of the number of A molecules (i.e. to an increase of intrinsic noise). Moreover, the deviations cannot be attributed either to a failure of the quasisteady approximation introduced in Eq. (7) because the larger the diffusion, the more accurate that approximation is. As indicated by equation (18), the transcriptional noise behaves as <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M54','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M54">View MathML</a>. Thus, for large values of the diffusion rate, the transcriptional noise level decreases. Therefore, we must conclude that the deviations between the theoretical and the numerical approaches become accentuated as the diffusion increases because there is a drop of the fluctuations related to the mRNA dynamics. This also indicates that for large enough diffusion rate, the intrinsic noise constitutes the main source of stochasticity.

In order to ensure that the intrinsic fluctuations are not actually increasing due to diffusion we perform the following in silico experiment. We consider a modification of our system such that a single mRNA molecule transcript leads to two autoinducer molecules that are considered to be distinguishable. The latter can be experimentally achieved by placing two consecutive copies of the encoding sequence of the autoinducer synthetase labeled with different fluorescent tags in the operon. Thus, we double the set of equations (2)-(5) in order to account for <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M55','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M55">View MathML</a> and <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M56','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M56">View MathML</a> molecules synthesis at cell i due to a single mRNA transcript <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M57','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M57">View MathML</a>. Following [9], by plotting the distribution of <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M51','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M51">View MathML</a> as a function of <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M52','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M52">View MathML</a> a putative increase of the intrinsic fluctuations can be discerned by these means. Right column of Figure 4 displays the results in this regard. The width of the distribution in a direction perpendicular to the diagonal is a measure of the intrinsic fluctuations (see [9] for details). As shown, as the diffusion increases there is no amplification of this quantity.

Diffusion and Total Noise

It is interesting to place the previous result in the context of the total noise of the autoinducer concentration. Figure 5 reveals that <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M58','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M58">View MathML</a> shows a non-monotonic behavior. As a function of <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M44">View MathML</a> the total noise first increases, reaches a maximum, and then decreases as the diffusion becomes larger. While Figure 5 represents data for the γ2 parameter set, this behavior applies for all the values of <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M41','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M41">View MathML</a> and <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M44">View MathML</a> explored in our simulations (data not shown). We point out that if <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M58','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M58">View MathML</a> > 1 then the dispersion is larger than the mean. Under this circumstance the fluctuations can lead to catastrophic events. E.g., fluctuations are able to remove all the signaling molecules within the cell. Note that the analytical calculations, that just account for the transcriptional noise, are in agreement with the numerical simulations, that account for both the transcriptional and the intrinsic noise, for a large range of <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M44">View MathML</a> values. This indicates that the main contribution to the total fluctuations within a large range of diffusion values is the transcriptional noise. Yet, as mentioned above, the latter diminishes as the diffusion increases while the intrinsic fluctuations remain constant. Consequently, the contribution of the intrinsic noise must become more relevant than the mRNA stochasticity beyond some value of <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M44">View MathML</a>.

thumbnailFigure 5. Noise of the signaling molecule as a function of the diffusion coecient. Noise <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M58','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M58">View MathML</a> as a function of diffusion coefficient <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M44">View MathML</a> for the set of parameters γ2 (see Figure 2B): stochastic simulations (circles) and analytical expression. Eq. (18), (solid line). By using the decomposition <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M59','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M59">View MathML</a> the differences between the computational and the theoretical distributions quantifies the amount of intrinsic noise (squares). As evidenced by the linear regression (blue short-dashed line) the later remains constant and is the main contribution to the total noise only for large diffusion values, <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M44">View MathML</a> > 104 (see text).

We address this point quantitatively by calculating the relative importance of these noisy sources. To this end, we make use of the decomposition <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M59','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M59">View MathML</a> where <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M60','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M60">View MathML</a> and <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M61','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M61">View MathML</a>, stand respectively for the intrinsic and the transcriptional contributions to the total noise [54]. Thus, by subtracting the analytical expression given by Eq. (18) to the total noise obtained in the numerical simulations, we are able to compute the intrinsic noise as a function of the diffusion. By performing a linear regression of the points that corresponds to the intrinsic noise we obtain that the slope of the curve is indeed zero in practical terms, ~ 2 · 10−7. Therefore, in agreement with the qualitative results obtained in Figure 4 (right column), the intrinsic noise remains constant as the diffusion increases, <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M60','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M60">View MathML</a> = 0.054 ± 0.003, and is the main stochastic component if <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M44">View MathML</a> > 104.

Discussion

As a matter of discussion, our modeling consider passive diffusion as the mechanism for the transport of the signaling molecule. This is indeed true in many QS systems. However, in other cases the autoinducer is actively transported in and out of the cell. For example, in the bacterial species Pseudomonas aeruginosa, C4-HSL can freely diffuse but C12-HSL, a larger signaling molecule, is subjected to active influx and efflux at rates of ~ 10−2 min−1 and ~ 10−1 min−1 respectively [55]. Other example corresponds to the AI-2 signaling molecule. The latter is present in many Gram-positive and Gram-negative species and it is believed to allow for interspecies communication [56]. In E. coli and Salmonella enterica extracellular AI-2 accumulates during the exponential phase, but then decreases drastically upon entry into the stationary phase. This reduction is due to the import and processing of AI-2 by the Lsr transporter [56,57]. Moreover, excretion from the cell of this autoinducer also appears to be an active process involving the putative transport protein YdgG (or alternatively named TqsA) [58]. In the case of E. coli these rates have been estimated by computational and experimental means: Dout ≃ 10-1 min−1 and Din ≃ 10-3 - 10-2 min−1 [59]. All in all, the diffusion rates when driven by active processes are four orders of magnitude smaller than the diffusion rates of small molecules through the membrane. Hence, the diffusion rates in QS systems can be categorized into two main, well separated, classes: small diffusion rates due to active transport, and large diffusion rates due to passive mechanisms.

In principle our model does not account for active diffusion processes, but transport driven by concentration differences. Still, it can be demonstrated that if the following condition holds, r Din/Dout, then our simple model fairly describes the active diffusion with D = Dout. However, that regime is not accessible in our simulations near the critical sensing potential since that would imply a non-physical situation (N < 1). If in any case we consider the rates of active transport of these QS systems, they would fit in our model with a normalized effective diffusion coefficient in the range <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M44">View MathML</a> ∈ [10-1, 10]. Note that, either driven by passive diffusion or by active transport, the region that maximizes the total noise, <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M58','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M58">View MathML</a>, is not accessible: <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M44">View MathML</a> ~ 5 · 101 − 102. Notice also that this range of diffusion values corresponds to the dynamics that produce separated peaks in the distribution of the autoinducer (see Figure 4). While our modeling is certainly very simple and the derived consequences should be carefully taken, this observation suggest that bacteria have developed mechanisms for coping with the noise and keep their functional QS regime away from the region where <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M58','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M58">View MathML</a> > 1. In this regard, let us point out that in every form of information exchange the precision is key. If the precision of the information is fuzzy then the related biological function lacks robustness. Thus, our results point towards the direction that bacteria have adapted their communication mechanisms in order to improve the signal to noise ratio.

Diffusion can lessen the effects of fluctuations, both in eukaryotes [32] and prokaryotes [23,60]. This behavior is certainly obtained in our model: if the diffusion rate is larger than a given amount, then the noise decreases as the diffusion increases. However, to the best of our knowledge, the reverse effect, i.e. the noise increases as the diffusion increases, had not been reported. These opposed behaviors are responsible of the non-monotonic comportment of the noise described herein. One may wonder the reason underlying this phenomenology. By manipulating Eq. (18), it is easy to demonstrate that a) the slope of the noise at <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M44">View MathML</a> = 0 is always positive and b) in the white noise limit, τc → 0, the slope becomes null and the transcriptional noise decreases monotonously as the diffusion increases. Consequently, the observed behavior is due to the colored character of the transcriptional noise, i.e. due to a competition of temporal scales: those of the diffusion and the noise correlation time. This result opens the possibility of finding a similar phenomenology in other systems subjected to relevant levels of transcriptional noise and diffusive signals. In addition, it points up the relevance of non-memoryless noisy sources in biological systems [11]. Moreover, it shows that our results do not depend on details as the precise number of transcripts (one in our case). As long as there are distinct transcriptional phases inside cells with characteristic time scales, that is, presence versus absence of transcripts, the observed phenomenology is, qualitatively, the same. Yet, those phases are prone to appear when the autoinducer is produced at constitutive levels and the number of transcripts is small. Herein we restrict ourselves to the case of one transcript simply because we are able to obtain analytical results in that situation.

Conclusions

Herein we have explored the role played by cell-cell communication and transcriptional noise in QS systems near the activation threshold. Within this context, we have shown that the dynamics of the signaling molecule exhibits different behaviors depending on the diffusion coefficient. When increasing the rate of diffusion, the probability distribution of the autoinducer changes from single-peak distribution (sustained bursts dynamics), to monotonically decreasing distribution (bursts dynamics), to double-peak distribution (bistable dynamics), and finally to narrow single-peak distribution (diffusion-averaged dynamics).

In addition, we have shown that the mRNA dynamics plays a crucial role for regulating the total amount of molecular noise of the signaling molecules. Transcriptional noise is the main contribution to the total noise for a large range of diffusion values, <a onClick="popup('http://www.biomedcentral.com/1752-0509/5/11/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/5/11/mathml/M44">View MathML</a> < 104. Only for very large values of the diffusion the intrinsic noise is the major source of stochasticity. Due to a competition of temporal scales, the total noise shows a non-monotonic behavior as a function of the diffusion rate. For large values of the diffusion coefficient, the total noise decreases as the diffusion rate increases. In this regard, our results are to be compared to previously reported noise reduction mechanisms, as for example in the case of bistable genetic switches coupled by QS communication [60]. On the other hand, when the diffusion rate is small enough compared to the characteristic rate of transcription events, the total noise increases as the diffusion becomes larger. The values of the diffusion rates in QS systems fall into two distinctive categories: either large values corresponding to passive transport mechanism or small values when an active transport mechanism applies. Surprisingly, these two QS classes avoid diffusion rates that maximize the total noise. According to these results, we conjecture that bacteria have engineered the communication mechanism for reducing the signal-to-noise ratio and produce a more reliable information exchange.

Our final comment refers to the possibility of considering other sources of stochasticity. Cell-to-cell variability and extrinsic noise have been proved to act as an important contribution in many cell processes [16,25,61,62]. In the context of the problem studied herein, we can envision that variability, either at the level of the mRNA dynamics or at the level of the diffusion rate, can effectively lead to significant changes in the reported phenomenology. In addition, by considering additional steps in the synthesis of the autoinducer, the levels of intrinsic noise would increase. Whether or not this extra level of fluctuations, coupled with feedback regulation, may generate new effects in the framework of QS is not known. Work in those directions is in progress.

Authors' contributions

MW and JB contributed equally to construct the model, perform the analytical calculations and write the paper. Numerical simulations were carried out by MW. All authors read and approved the final manuscript.

Acknowledgements

We thank Oriol Canela Xandri and Pieter Rein ten Wolde for fruitful comments. Financial support was provided by the Spanish MICINN under grants FIS2009-11104 and BFU2010-21847-C02-01/BMC and by DURSI through project 2009-SGR/01055. M.W. acknowledges the support of the Spanish MICINN through a doctoral fellowship (FPU AP2008-03272) and the partial support of the Catalan government by means of a FI doctoral fellowship (2009FI_B 00433).

References

  1. Elf J, Li GW, Xie XS: Probing transcription factor dynamics at the single-molecule level in a living cell.

    Science 2007, 316(5828):1191-4. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  2. Güell M, van Noort V, Yus E, Chen WH, Leigh-Bell J, Michalodimitrakis K, Yamada T, Arumugam M, Doerks T, Kühner S, Rode M, Suyama M, Schmidt S, Gavin AC, Bork P, Serrano L: Transcriptome complexity in a genome-reduced bacterium.

    Science (New York, NY) 2009, 326(5957):1268-71. OpenURL

  3. Kæ rn M, Elston T, Blake W, Collins J: Stochasticity in gene expression: from theories to phenotypes. [http://www.nature.com/nrg/journal/v6/n6/abs/nrg1615.html] webcite

    Nature Reviews Genetics 2005, 6(6):451-464. PubMed Abstract | Publisher Full Text OpenURL

  4. Süel GM, Garcia-Ojalvo J, Liberman LM, Elowitz MB: An excitable gene regulatory circuit induces transient cellular differentiation.

    Nature 2006, 440(7083):545-50. PubMed Abstract | Publisher Full Text OpenURL

  5. Raj A, van Oudenaarden A: Nature, nurture, or chance: stochastic gene expression and its consequences.

    Cell 2008, 135(2):216-26. PubMed Abstract | Publisher Full Text OpenURL

  6. Brock A, Chang H, Huang S: Non-genetic heterogeneity - a mutation-independent driving force for the somatic evolution of tumours.

    Nature reviews. Genetics 2009, 10(5):336-42. PubMed Abstract | Publisher Full Text OpenURL

  7. Ao P: Global view of bionetwork dynamics: adaptive landscape.

    J Genet Genomics 2009, 36(2):63-73. PubMed Abstract | Publisher Full Text OpenURL

  8. Eldar A, Elowitz MB: Functional roles for noise in genetic circuits.

    Nature 2010, 467(7312):167-73. PubMed Abstract | Publisher Full Text OpenURL

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

    Science (New York, NY) 2002, 297(5584):1183-6. Publisher Full Text OpenURL

  10. Ozbudak EM, Thattai M, Kurtser I, Grossman AD, van Oudenaarden A: Regulation of noise in the expression of a single gene.

    Nature genetics 2002, 31:69-73. PubMed Abstract | Publisher Full Text OpenURL

  11. Rosenfeld N, Young JW, Alon U, Swain PS, Elowitz MB: Gene regulation at the single-cell level.

    Science (New York, NY) 2005, 307(5717):1962-5. Publisher Full Text OpenURL

  12. Raj A, van Oudenaarden A: Single-molecule approaches to stochastic gene expression.

    Annual review of biophysics 2009, 38:255-70. PubMed Abstract | Publisher Full Text OpenURL

  13. Larson DR, Singer RH, Zenklusen D: A single molecule view of gene expression.

    Trends in cell biology 2009, 19(11):630-7. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. Golding I, Paulsson J, Zawilski SM, Cox EC: Real-time kinetics of gene activity in individual bacteria.

    Cell 2005, 123(6):1025-36. PubMed Abstract | Publisher Full Text OpenURL

  15. Shahrezaei V, Swain PS: Analytical distributions for stochastic gene expression.

    Proceedings of the National Academy of Sciences of the United States of America 2008, 105(45):17256-61. PubMed Abstract | Publisher Full Text OpenURL

  16. Zhu XM, Yin L, Hood L, Ao P: Robustness, stability and efficiency of phage lambda genetic switch: dynamical structure analysis.

    J Bioinform Comput Biol 2004, 2(4):785-817. PubMed Abstract | Publisher Full Text OpenURL

  17. Cağatay T, Turcotte M, Elowitz MB, Garcia-Ojalvo J, Süel GM: Architecture-dependent noise discriminates functionally analogous differentiation circuits.

    Cell 2009, 139(3):512-22. PubMed Abstract | Publisher Full Text OpenURL

  18. Kepler T, Elston T: Stochasticity in transcriptional regulation: origins, consequences, and mathematical representations.

    Biophysical Journal 2001, 81(6):3116-3136. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  19. Wilkinson DJ: Stochastic modelling for quantitative description of heterogeneous biological systems.

    Nature reviews. Genetics 2009, 10(2):122-33. PubMed Abstract | Publisher Full Text OpenURL

  20. Garcia-Ojalvo J, Elowitz M, Strogatz S: Modeling a synthetic multicellular clock: Repressilators coupled by quorum sensing.

    Proceedings of the National Academy of Sciences of the United States of America 2004, 101(30):10955. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  21. Goryachev AB, Toh DJ, Wee KB, Zhang HB, Zhang LH, Lee T: Transition to quorum sensing in an Agrobacterium population: A stochastic model.

    PLoS computational biology 2005, 1(4):e37. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Zhou T, Chen L, Aihara K: Molecular Communication through Stochastic Synchronization Induced by Extracellular Fluctuations.

    Physical Review Letters 2005, 95(17):2-5. Publisher Full Text OpenURL

  23. Tanouchi Y, Tu D, Kim J, You L: Noise reduction by diffusional dissipation in a minimal quorum sensing motif.

    PLoS computational biology 2008, 4(8):e1000167. PubMed Abstract | PubMed Central Full Text OpenURL

  24. Ullner E, Buceta J, Díez-Noguera A, García-Ojalvo J: Noise-Induced Coherence in Multicellular Circadian Clocks.

    Biophysical Journal 2009, 96(May):3573-3581. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. Canela-Xandri O, Sagués F, Buceta J: Interplay between Intrinsic Noise and the Stochasticity of the Cell Cycle in Bacterial Colonies.

    Biophysical journal 2010., 98(11) PubMed Abstract | Publisher Full Text OpenURL

  26. Bassler BL, Losick R: Bacterially speaking.

    Cell 2006, 125(2):237-46. PubMed Abstract | Publisher Full Text OpenURL

  27. Ng WL, Bassler BL: Bacterial quorum-sensing network architectures.

    Annual review of genetics 2009, 43:197-222. PubMed Abstract | Publisher Full Text OpenURL

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

    Nature 2004, 428(6985):868-71. PubMed Abstract | Publisher Full Text OpenURL

  29. Basu S, Gerchman Y, Collins C, Arnold F, Weiss R: A synthetic multicellular system for programmed pattern formation.

    Nature 2005, 434(7037):1130-1134. PubMed Abstract | Publisher Full Text OpenURL

  30. Balagaddé FK, Song H, Ozaki J, Collins CH, Barnet M, Arnold FH, Quake SR, You L: A synthetic Escherichia coli predator-prey ecosystem.

    Molecular systems biology 2008, 4(187):187. PubMed Abstract | PubMed Central Full Text OpenURL

  31. Danino T, Mondragón-Palomino O, Tsimring L, Hasty J: A synchronized quorum of genetic clocks.

    Nature 2010, 463(7279):326-330. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  32. Erdmann T, Howard M, ten Wolde PR: Role of spatial averaging in the precision of gene expression patterns.

    Phys Rev Lett 2009, 103(25):258101. PubMed Abstract | Publisher Full Text OpenURL

  33. Yu J, Xiao J, Ren X, Lao K, Xie XS: Probing gene expression in live cells, one protein molecule at a time.

    Science (New York, NY) 2006, 311(5767):1600-3. Publisher Full Text OpenURL

  34. Boyer M, Wisniewski-Dyé F: Cell-cell signalling in bacteria: not simply a matter of quorum.

    FEMS microbiology ecology 2009, 70:1-19. PubMed Abstract | Publisher Full Text OpenURL

  35. Ravn L, Christensen A, Molin S, Givskov M, Gram L: Methods for detecting acylated homoserine lactones produced by Gram-negative bacteria and their application in studies of AHL-production kinetics.

    Journal of Microbiological Methods 2001, 44(3):239-251. PubMed Abstract | Publisher Full Text OpenURL

  36. Peccoud J, Ycart B: Markovian modeling of gene-product synthesis.

    Theoretical Population Biology 1995, 48(2):222-234. Publisher Full Text OpenURL

  37. Karmakar R, Bose I: Graded and binary responses in stochastic gene expression.

    Physical biology 2004, 1(3-4):197-204. PubMed Abstract | Publisher Full Text OpenURL

  38. Iyer-Biswas S, Hayot F, Jayaprakash C: Stochasticity of gene products from transcriptional pulsing.

    Physical Review E 2009, 79(3):1-9. Publisher Full Text OpenURL

  39. Horsthemke W, Lefever R: Noise-induced transitions: theory and applications in physics, chemistry and biology. Berlin Heidelberg: Springer-Verlag; 1984. OpenURL

  40. More MI, Finger LD, Stryker JL, Fuqua C, Eberhard A, Winans SC: Enzymatic Synthesis of a Quorum-Sensing Autoinducer Through Use of Defined Substrates.

    Science 1996, 272(5268):1655-1658. PubMed Abstract | Publisher Full Text OpenURL

  41. Parsek MR, Val DL, Hanzelka BL, Cronan JE, Greenberg EP: Acyl homoserine-lactone quorum-sensing signal generation.

    Proceedings of the National Academy of Sciences of the United States of America 1999, 96(8):4360-5. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  42. Shahrezaei V, Swain PS: Analytical distributions for stochastic gene expression.

    Proc Natl Acad Sci USA 2008, 105(45):17256-61. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  43. Gillespie D: Exact stochastic simulation of coupled chemical reactions.

    The journal of physical chemistry 1977, 81(25):2340-2361. Publisher Full Text OpenURL

  44. Tsimring LS, Volfson D, Hasty J: Stochastically driven genetic circuits.

    Chaos (Woodbury, NY) 2006, 16(2):026103. Publisher Full Text OpenURL

  45. McMillen D, Kopell N, Hasty J, Collins J: Synchronizing genetic relaxation oscillators by intercell signaling.

    Proceedings of the National Academy of Sciences of the United States of America 2002, 99(2):679. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  46. Pai A, You L: Optimal tuning of bacterial sensing potential.

    Molecular systems biology 2009, 5(286):286. PubMed Abstract | PubMed Central Full Text OpenURL

  47. Boedicker JQ, Vincent ME, Ismagilov RF: Microfluidic confinement of single cells of bacteria in small volumes initiates high-density behavior of quorum sensing and growth and reveals its variability.

    Angewandte Chemie (International ed in English) 2009, 48(32):5908-11. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  48. Carnes EC, Lopez DM, Donegan NP, Cheung A, Gresham H, Timmins GS, Brinker CJ: Confinement-induced quorum sensing of individual Staphylococcus aureus bacteria.

    Nature chemical biology 2010, 6:41-5. PubMed Abstract | Publisher Full Text OpenURL

  49. Hagen SJ, Son M, Weiss JT, Young JH: Bacterium in a box: sensing of quorum and environment by the LuxI/LuxR gene regulatory circuit.

    Journal of Biological Physics 2010, 36(3):317-327. Publisher Full Text OpenURL

  50. Bennett MR, Hasty J: Microfluidic devices for measuring gene network dynamics in single cells.

    Nature reviews. Genetics 2009, 10(9):628-38. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  51. Kaufmann G, Sartorio R, Lee S, Rogers C, Meijler M, Moss J, Clapham B, Brogan A, Dickerson T, Janda K: Revisiting quorum sensing: discovery of additional chemical and biological functions for 3-oxo-N-acylhomoserine lactones.

    Proceedings of the National Academy of Sciences 2005, 102(2):309. Publisher Full Text OpenURL

  52. Roberts C, Anderson KL, Murphy E, Projan SJ, Mounts W, Hurlburt B, Smeltzer M, Overbeek R, Disz T, Dunman PM: Characterizing the effect of the Staphylococcus aureus virulence factor regulator, SarA, on log-phase mRNA half-lives.

    Journal of bacteriology 2006, 188(7):2593-603. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  53. Anderson KL, Dunman PM: Messenger RNA Turnover Processes in Escherichia coli, Bacillus subtilis, and Emerging Studies in Staphylococcus aureus.

    International journal of microbiology 2009, 2009:525491. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  54. Swain P, Elowitz M, Siggia E: Intrinsic and extrinsic contributions to stochasticity in gene expression.

    Proceedings of the National Academy of Sciences 2002, 99(20):12795. Publisher Full Text OpenURL

  55. Pearson J, Van Delden C, Iglewski B: Active efflux and diffusion are involved in transport of Pseudomonas aeruginosa cell-to-cell signals.

    Journal of bacteriology 1999, 181(4):1203. PubMed Abstract | PubMed Central Full Text OpenURL

  56. Xavier K, Bassler B: Regulation of uptake and processing of the quorum-sensing autoinducer AI-2 in Escherichia coli.

    Journal of bacteriology 2005, 187:238. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  57. Wang L, Hashimoto Y, Tsao C, Valdes J, Bentley W: Cyclic AMP (cAMP) and cAMP receptor protein influence both synthesis and uptake of extracellular autoinducer 2 in Escherichia coli.

    Journal of bacteriology 2005, 187(6):2066. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  58. Herzberg M, Kaye I, Peti W, Wood T: YdgG (TqsA) controls biofilm formation in Escherichia coli K-12 through autoinducer 2 transport.

    Journal of bacteriology 2006, 188(2):587. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  59. Li J, Wang L, Hashimoto Y, Tsao CY, Wood TK, Valdes JJ, Zafiriou E, Bentley WE: A stochastic model of Escherichia coli AI-2 quorum signal circuit reveals alternative synthesis pathways.

    Molecular systems biology 2006, 2:67. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  60. Koseska A, Zaikin A, Kurths J, García-Ojalvo J: Timing cellular decision making under noise via cell-cell communication.

    PloS one 2009, 4(3):e4872. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  61. Lu T, Shen T, Bennett MR, Wolynes PG, Hasty J: Phenotypic variability of growing cellular populations.

    Proceedings of the National Academy of Sciences of the United States of America 2007, 104(48):18982-7. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  62. Shahrezaei V, Ollivier J, Swain P: Colored extrinsic fluctuations and stochastic gene expression.

    Molecular Systems Biology 2008, 4:196. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL