Rule–based regulatory and metabolic model for Quorum sensing in P. aeruginosa
https://doi.org/10.1186/1752-0509-7-81
© Schaadt et al.; licensee BioMed Central Ltd. 2013
Received: 28 January 2013
Accepted: 9 August 2013
Published: 21 August 2013
Abstract
Background
In the pathogen P. aeruginosa, the formation of virulence factors is regulated via Quorum sensing signaling pathways. Due to the increasing number of strains that are resistant to antibiotics, there is a high interest to develop novel antiinfectives. In the combat of resistant bacteria, selective blockade of the bacterial cell–to–cell communication (Quorum sensing) has gained special interest as anti–virulence strategy. Here, we modeled the las, rhl, and pqs Quorum sensing systems by a multi–level logical approach to analyze how enzyme inhibitors and receptor antagonists effect the formation of autoinducers and virulence factors.
Results
Our rule–based simulations fulfill the behavior expected from literature considering the external level of autoinducers. In the presence of PqsBCD inhibitors, the external HHQ and PQS levels are indeed clearly reduced. The magnitude of this effect strongly depends on the inhibition level. However, it seems that the pyocyanin pathway is incomplete.
Conclusions
To match experimental observations we suggest a modified network topology in which PqsE and PqsR acts as receptors and an autoinducer as ligand that up–regulate pyocyanin in a concerted manner. While the PQS biosynthesis is more appropriate as target to inhibit the HHQ and PQS formation, blocking the receptor PqsR that regulates the biosynthesis reduces the pyocyanin level stronger.
Keywords
Background
Quorum sensing (QS) describes how the communication between bacteria is established. Thus, the regulation of genes is adapted to cell population density through the activity of a combined regulatory and metabolic network. In P. aeruginosa, QS is involved in the formation of biofilms and in the production of a large number of virulence factors, such as elastase, rhamnolipids, and pyocyanin. These virulence factors are responsible e.g. for tissue damage and disruption of tight junctions. P. aeruginosa usually infects patients with immune system deficiencies. Since an increasing number of infecting strains are resistant to most current antibiotics, there is a large interest in developing novel antibacterial strategies. It has been suggested that selectively targeting the QS machinery by signaling molecule inhibitors may be advantageous over antibiotics that target central metabolism or DNA replication with respect to the development of resistance mutations because the former strategies have no impact on bacterial viability delay [1, 2].
QS network of P.aeruginosa . Regulation of the virulence factors elastase (LasB), rhamnolipids (Rhm2), and pyocyanin by the hierarchically organized systems las (blue), rhl (green), and pqs (red). Colored balls represent signaling molecules, squares denote enzymes, and colored rectangles are symbols for receptors or other proteins.
The pqs system (in Figure 1 colored in red) uses the Pseudomonas quinolone signal (PQS) that is synthesized from HHQ by the enzyme PqsH. Both HHQ and PQS are able to form complexes with the receptor PqsR (in the following denoted as C5 and C3) that regulate many genes, such as the biosynthesis operon pqsA–E[17, 18]. As first step, the enzyme PqsA produces anthraniloyl–CoA using anthranilic acid. Then, PqsD forms HHQ from ACoA and β–ketodecanoic acid or its bioactivated thioesters [19]. The enzymes PqsB and PqsC are also required to synthesize HHQ whereby their rule is not completely understood [20, 21]. Finally, PqsE acts as negative regulator of the pqsA–E operon [20]. In this study, we do not include further regulators related to the QS machinery. For example, it was shown that QscR represses the transcription of lasI[22] and VqsR negatively regulates the expression of qscR[23].
The dynamic behavior of cellular regulatory networks can be analyzed with a wide range of theoretical approaches: Boolean networks, Kauffman networks, probabilistic Boolean networks, dynamic Bayesian networks, Petri Nets, as well as differential equation models [24–26]. Differential equation models have the advantage of enabling a quantitative analysis. However, for many biochemical reactions the actual rate constants are unknown. Therefore, a Boolean model was used, e.g. to understand the FA/BRCA pathway [27]. A generalized logical formalism with multi–levels was applied to the Gap–gene System in Drosophila[28] and to analyze the cell cycle of budding yeast [29].
The dynamic behavior of the QS systems was modeled so far either as simple feedback loop or considering the Pseudomonas las and rhl systems using ordinary as well as partial differential equations [30, 31] or regarding the lux system of Vibrio fischeri applying so–called P systems [32]. Anguige et al. included a LasR degradation drug in their differential equation approach of the las system [33]. Furthermore, the development of biofilms was studied using the las system [34] or a 3D growth model of a self–producing signaling molecule including inhibition [35].
In this work, we implemented a multi–level logical approach and compared the influence of enzyme inhibitors and that of receptor antagonists on the formation of autoinducers and virulence factors. Here, different levels of inhibition were considered. Additionally, we analyzed the topology of the network. For this purpose, we modeled the QS in P. aeruginosa comprising the las, rhl, and pqs systems as well as the virulence factors elastase, rhamnolipids, and pyocyanin [36, 37].
Methods
QS network as a Boolean topology. Black edge = threshold is 1; blue edge = state of underlined node must be at least 2; orange edge = state of underlined node must be at least 3; green and thin edge = state of underlined node must be at least 4; numbers denote possible states for a node; dotted arrows are reactions involved in a transport process; red and thick edge = happens after a certain number of time steps (degradation). Gray and dashed arrows denote reactions that occurs by chance with a certain probability. Nodes named C represent a complex between autoinducer and receptor, C:G is the complex bound to an operon.
Computational analysis
Flow chart. Iteration over updating of states, transport, and at certain time points growth (represented as generation of new cells). The updating of states is performed only for cells with a life time longer than a delay time of ten. Here, t denotes the time step and d the number of cell divisions.
whereby W=(w i j ) is a weighting matrix and contains all relations between the nodes. Its entry w i j is when node j is one of l nodes that are together required to activate node i. If node j deactivates node i then w i j is −k where k is the number of reactions that activate node i. Otherwise, w i j is zero. In contrast to some other Boolean models, the second case implies that one inhibiting reaction overrules the activating reactions of node i. Also, a node that is not explicitly activated at a certain time step is inactive. As our network only included deactivating nodes in cases where there is a good evidence that they inhibit certain activations, the weighting scheme makes sure that deactivations are stronger than activations. As the concentration of the autoinducers must reach a certain threshold ε to start the positive feedback loop, some reactions in our model are only activated if a critical threshold of certain metabolites is exceeded. Function K controls that the state of every node j that is related to node i has reached the required threshold at time t. Function H rounds the weighted sum to an integer valued output. For multi–level nodes with more than two possible states, i.e. autoinducers, RsaL, and virulence factors the level at the current time step x i (t) is added via function G. Function M ensures that no node adopts a higher value than its maximum possible value. The explicit forms of the functions M, H, K, and G are given in Additional file 4.
After updating all nodes, the states of multi–level nodes are decreased by the number of reactions in which they were used. Let us assume, that for example, HHQ forms a complex with PqsR and HHQ is additionally required to yield PQS. Then, in this step the level of HHQ is decreased by two. Those multi–level nodes are subject to a further reaction that decreases the state value after every 20 time steps to account for degradation (see Section Choice of parameters). This degradation is realized equally for all multi–level nodes independent of their maximum possible states. A time step represents either the whole formation of a protein including transcription and translation or a faster enzymatic reaction. Hence, a degradation happens infrequently. However, a node has a relatively small number of possible states, such that a single degradation process consumes a comparably large part of its concentration. For example, in the case of RsaL a degradation decreases the level by two–thirds. Transport reactions (dotted lines in Figure 2) are neither considered in the logical updating scheme nor in the decreasing step afterwards. Additionally, we included some random reactions for processes that take place either infrequently or very slowly in comparison to the other processes, for example, the transformation of HHQ to PQS [18]. For this, we used a random reaction that happens about every fourth step. For simplicity, the enzymes PqsB, PqsC, and PqsD are considered as one node.
The communication between different cells is based on the exchange of the signaling molecules. For this, we modeled a simple transport process that occurs after updating all states inside the cells. When the internal concentration of an autoinducer is higher than the transport threshold, it is transported out of the cell except for the case where the external concentration relative to the number of cells is higher than the internal concentration. Otherwise, for external levels larger than zero, an autoinducer is transported into the cell. In the case that there is no autoinducer outside, there is also no transport. We do not distinguish between diffusion and efflux pumps or active transports and apply usually a transport threshold of three, i.e. if the autoinducer concentration is high.
Knock–out mutants, enzyme inhibitors and receptor antagonists
Realization of inhibitors. The inhibitor that is initialized as being active blocks its target with a certain probability. The time point when the dashed and gray colored reaction happens is realized randomly.
Result and discussion
We constructed a combined regulatory and metabolic network for the three QS systems in P. aeruginosa. The realization as a dynamic rule–based computational model also enabled us to study the effects of gene knock–outs and those of enzyme inhibitors as well as receptor antagonists.
Choice of parameters
For the sake of simplicity and interpretability, the model was designed as a minimal system. A synchronous, parameter–free Boolean logical analysis was applied where nodes can be either “on” or “off” except for some nodes and edges that are discussed below in detail. In a Boolean network with a synchronous updating scheme, there is no time scale so that all reactions take place at the same time with the same rate. This is a valid approach unless this alters the system behavior in a crucial waycompared to either considering different reaction time scales, or if certain reactions depend on special conditions. We will show in the following that the first assumption (same time scale) is not crucial for the system properties studied here. The second point, however, required a generalization of the Boolean model by introducing multi-level logical functions at various points of our model.
We will now discuss the elements of the systems where more than two states were found necessary. (i) The regulator RsaL blocks the transcription of the autoinducer synthase LasI. Hence, three different RsaL levels were used, namely “not available” (0), “available at low concentrations” (1), and “sufficient to block LasI” (2) are used. (ii) Besides PQS, HHQ is also able to bind to PqsR, and to activate PqsR, but PQS has a 100–fold higher affinity than HHQ [18]. Therefore, PQS was set to inactivate the complex C5 of HHQ and PqsR with a threshold of two representing a large concentration of PQS. (In Figure 2, all reactions with a threshold of two are colored in blue). (iii) Since the activation of gene expression (activation of Ci:Gi by Ci) is concentration–dependent, the complexation between autoinducer and receptor occurs only when the autoinducer concentration has reached the lower threshold of two. (iv) Moreover, the autoinducer of the las system AI–1 is able to form a weak complex with the second receptor RhlR as long as the autoinducer of the rhl system AI–2 has not reached a very high concentration. Thus, formation of the complex C4 between AI–1 and RhlR required at least four AI–1 and its inactivation at least four AI–2 (edges colored in green in Figure 2). (v) Due to the accumulation of autoinducers and virulence factors, six different levels were considered for the internal signaling molecules AI–1, AI–2, HHQ, as well as PQS and the virulence factors. The external autoinducer levels were allowed to adopt arbitrarily large numbers.
Conversion frequency
Effect of the PQS production rate. Effect on the average levels of the autoinducers HHQ and PQS, as well as that of pyocyanin. The reported values are mean levels in the time interval 100 to 600 averaged over ten runs obtained with different random numbers and their standard deviations. The conversion frequency represents the probability that PQS is produced by PqsH. The degradation frequency was set to once per 20 time steps and the transport threshold to three.
Degradation
If no degradation steps were included, all metabolites would be present for the whole time after their production. Hence, the state levels of the virulence factors would be at maximal levels since they have no other out–going edges. This is not realistic because metabolites are generally degraded after a certain time in a cell or may be consumed by reactions not considered in our network. Therefore, we have also considered degradation steps in our model that occur after the states have been updated according to the producing reactions. After a degradation step, no node reaches its maximum level anymore. We found that if degradation takes place after every time step, the three QS systems become inactive except for the proteins Vfr and LasR due to the self–loop of Vfr and the relation between them and the applied initialization conditions.
Additional file 5: Figure S1 shows how different degradation frequencies between 1 and influence the levels of virulence factors. If degradation occurs after five or more time steps the degradation frequency generally has only a very small effect. We compared the virulence factor level with the theoretical possible maximum values considering the numbers of degradations in the simulation time considered. The smaller the frequency the larger is the effect of the autoinducer degradation on the virulence factor level. The pyocyanin level comes close to the maximal level for frequencies of or less and reaches this value at a degradation frequency of or less. Therefore, we chose a degradation frequency of for all following computations.
Transport threshold
In our model, transport is subject to concentration gradients between individual cells and the exterior. For transport out of a cell to occur, the autoinducer level must also overcome a certain threshold. Additional file 6: Figure S2 shows the influence of this threshold on the level of virulence factors. In general, the effect of different threshold values between 1 and 5 is negligible, in particular when considering the time interval from 100 to 600. As expected, the external autoinducer levels decrease with an increasing transport threshold, whereas the internal levels increase. For thresholds larger or equal to four, the levels of external HHQ and PQS are zero. If the transport threshold is set to one, the internal HHQ and PQS levels are so small that the pyocyanin level is slightly decreased. In the following, we used a transport threshold of three, which corresponds to a relatively high autoinducer concentration.
Initialization and starting behavior
In this study, we are primarily interested in identifying attractors in which the QS systems are active. Obviously, when all initial states of the network are set to zero (inactive or not available), the whole system will remain inactive. The las system becomes active when (1) one or both of the global activators GacA and Vfr are initially set to one and when (2) the complex C1 between AI-1 and LasR binds to the las operon (C1:G1 is set to one). The rhl system is activated with a certain time delay. This delay can be avoided when GacA, C1:G1, and additionally C2:G2 (complex C2 between AI–2 and RhlR binds to the rhl operon) are active in the beginning. If in the beginning C5:G3 and C3:G3 are activated as well, the pqs system becomes active after a short delay, too. Also, switching on PqsR avoids the delay. In general, the initialization of Cx:Gy (complex Cx bound to operon Gy) to one is necessary since basal transcription rates are not included in the logical approach used here. In our system, the basal transcriptions could be included either with the same rate as the activated gene expression or in a stochastic way as often as degradation. The second case does not result in an initiation of the system. Alternatively, it is possible to start with elevated autoinducer levels or to begin with externally available autoinducers that are able to diffuse or can be transported into the cells.
When more nodes are initially active, the initialization phase of the system takes a different route because it is no longer required to start up the whole system. However, at later times, the behavior and the levels are the same as using minimal initial conditions. Additional file 7: Figure S3 shows that, after about 100 time steps, the metabolite levels reached in a computation with all nodes initially activated are very similar to those reached with minimal initial conditions. For all other further considerations, we used minimal initial conditions in which only Vfr, C1:G1, C3:G3, and C5:G3 are initially set to one.
Additional file 8: Table S4 shows an example trajectory where the system has reached a stable regime after time 100. Therefore, we can consider all time steps after 100 for computing equilibrated average values. For statistical reasons, we carried out all simulations until a time of 600. The differences of the levels using different growth rates were found to be very small. We compared growth processes using different division frequencies at time steps 60, 120, and 240 and modified the simulation time such that in each case about 200 time steps in stationary phase are possible. For most internal nodes, the levels were the same. Due to the stochastically modeled formation of PQS, the levels differ from each other in the same range as between different runs that used the same growth parameters but different random numbers for the pqs system. The external levels also display only negligible variations due to degradations. For example, the external level of autoinducer AI–1 divided by the total number of cells is 3.1 for a growth rate of 60, 3.3 for a growth rate of 120, and 2.9 for a growth rate of 240 averaged over the last 100 time steps. When comparing the behavior of a single cell over 600 time steps to that of 64 cells the internal and external levels are slightly changed. However, the qualitative behavior remains unchanged as well. For a single cell, the external level of AI–1 averaged over the last 100 time steps is 2.6. Hence, the behavior of the model is practically independent of the growth processes. In the following, we used a growth rate of 60 time steps and 6 division processes.
Start behavior. Based on the complexes using minimal initial conditions. The first complex C1 between AI–1 and LasR is shown in red and the second complex C2 between AI–2 and RhlR in blue. The levels averaged over ten different runs for complexes of the pqs system are shown in orange (C3 between PQS and PqsR) and in green (C5 between HHQ and PqsR). Dotted lines denote the first peak of the same colored complex.
Behavior of wild type and knock–out mutants
Behavior of wild type. Lines denote the maximum value for bars with the same color. Shown are average levels in the time interval 100 to 600 averaged over ten runs obtained with different random numbers and their standard deviations.
Detailed view of the pqs system. Black edge = threshold is 1; blue edge = state of underlined node must be at least 2; orange edge = state of underlined node must be at least 3; numbers denote the number of possible states for a node; dotted arrows are reactions involved in a transport process; red and thick edge = happens after 20 time steps (degradation). gray and dashed edge denotes a reaction that occurs by chance with a certain probability.
Calculated PQS and pyocyanin levels for wild type and knock–out mutants. Displayed are average values in the time interval 100 to 600 averaged over ten runs obtained with different random numbers and their standard deviations.
Behavior of enzyme inhibitors and receptor antagonists
External HHQ and PQS levels of wild type and inhibitors with varying inhibition levels. Shown are average values in the time interval from 100 to 600 averaged over ten runs obtained with different random numbers and their standard deviations.
Enzyme inhibition. Inhibition of PqsBCD with varying inhibition levels. Shown are average values in the time interval from 100 to 600 averaged over ten runs obtained with different random numbers with and standard deviations. The dotted line represents the theoretical maximum possible values.
Receptor antagonists. Inhibition of PqsR with varying inhibition levels. Shown are average values in the time interval from 100 to 600 averaged over ten runs obtained with different random numbers and their standard deviations. The dotted line represents the theoretical maximum possible values.
Updated network
While the behavior of the HHQ and PQS levels mostly fulfills experimental observations, the discrepancy for pyocyanin between predicted and observed behavior suggests that pyocyanin production is co–regulated by further functionally unknown proteins that may or may not be connected to the pqs system. Here, we only considered the connected case. Thus, we proceeded by studying networks with modified connectivities. For simplicity, we did not consider negative regulations or varying conditions such as reaction rates or modifications in the number of states for certain nodes. Instead of the feed–forward link pqs–member → X → pyocyanin we used a direct link from pqs–member → pyocyanin such that we do not have to include new hypothetical proteins. We further assumed that relations are only possible between a protein acting as receptor and a corresponding signaling molecule.
Possible reactions to form pyocyanin
Notation | Reaction |
---|---|
R1 | PqsE → pyocyanin |
R2 | C3:G3 → pyocyanin |
R3 | C5:G3 → pyocyanin |
R4 | PqsR → pyocyanin |
R5 | PqsH → pyocyanin |
R6 | PqsA → pyocyanin |
R7 | PqsBCD → pyocyanin |
R8 | PqsR + DHQ → pyocyanin |
R9 | PqsH + HHQ → pyocyanin |
R10 | PqsH + PQS → pyocyanin |
R11 | PqsH + DHQ → pyocyanin |
R12 | PqsA + HHQ → pyocyanin |
R13 | PqsA + PQS → pyocyanin |
R14 | PqsA + DHQ → pyocyanin |
R15 | PqsBCD + HHQ → pyocyanin |
R16 | PqsBCD + PQS → pyocyanin |
R17 | PqsBCD + DHQ → pyocyanin |
R18 | PqsE + HHQ → pyocyanin |
R19 | PqsE + PQS → pyocyanin |
R20 | PqsE + DHQ → pyocyanin |
R21 | PqsR + PqsE + HHQ → pyocyanin |
R22 | PqsR + PqsE + PQS → pyocyanin |
R23 | PqsR + PqsE + DHQ → pyocyanin |
R24 | PqsR + PqsE → pyocyanin |
Behavior of updated networks
Network | Used reactions | Results | |||
---|---|---|---|---|---|
PqsA−– PqsBCD− | PqsE− | PqsR antagonists | PqsBCD inhibitors | ||
N1 | R1, R2 | - - - - - | Pyocyanin | Pyocyanin | - - - - - |
N2 | R18 | - - - - - | - - - - - | Pyocyanin | - - - - - |
N3 | R19 | - - - - - | - - - - - | Pyocyanin | PQS |
N4 | R20, R2 | - - - - - | Pyocyanin | Pyocyanin | - - - - - |
N5 | R1, R2, R4 | Pyocyanin | Pyocyanin | Pyocyanin | - - - - - |
N6 | R18, R4 | Pyocyanin | Pyocyanin | Pyocyanin | - - - - - |
N7 | R19, R4 | Pyocyanin | Pyocyanin | Pyocyanin | PQS |
N8 | R20, R2, R4 | Pyocyanin | Pyocyanin | Pyocyanin | - - - - - |
N9 | R21 | - - - - - | - - - - - | - - - - - | - - - - - |
N10 | R22 | - - - - - | - - - - - | - - - - - | PQS |
N11 | R23, R2 | - - - - - | Pyocyanin | Pyocyanin | - - - - - |
N12 | R24, R2 | - - - - - | Pyocyanin | Pyocyanin | - - - - - |
Modified network termed N3. Numbers denote possible states for a node; red colored elements denote the number of new nodes or edges, whereas blue dotted lines are edges of the original network removed in this network. Dashed lines represent reactions given in Figure 8.
In the networks N5, N6, N7, and N8, PqsA −–PqsBCD − and PqsE − mutants produce pyocyanin at maximum level due to the postulated reaction R4, see Table 2. According to this or–linkage, PqsR is essential for the biosynthesis of pyocyanin, so that its level is not decreased in PqsA −–PqsBCD − and PqsE − mutants. The levels of all other metabolites are the same as for the corresponding network without reaction R4, e.g. N5 behaves as N1. Besides N1 PqsE − mutants also produce pyocyanin at maximum level in networks N4, N11, and N12. In the case of N2 and N3, where pyocyanin is activated by a complex between PqsE and HHQ or PQS, no pyocyanin is formed by knock–out mutants. However, weak PqsR antagonists again do not decrease the pyocyanin level.
We consider the networks N9 and N10 as being the closest to literature. In our model, PqsBCD inhibitors are more suitable to reduce the HHQ and PQS level and PqsR antagonists reduce the pyocyanin level at higher rate. In general, the levels strongly decrease with increasing inhibition level and the level of HHQ is clearly stronger decreased than the PQS and pyocyanin levels for the same inhibition level. In the case of network N10, a PqsBCD inhibitor strongly reduces the external HHQ level, whereas the effect on the external PQS level is clearly not as emphasized as for comparable inhibitors of the other networks. In reality, reactions R21 and R22 may be represented in two ways. (1) Two complexes up–regulate the transcription of two different enzymes. The first enzyme catalyzes the formation of a metabolite that is used by the second enzyme to form pyocyanin. (2) A first complex up–regulates the transcription of an enzyme that catalyzes the formation of a certain metabolite and a second complex up–regulates a receptor protein that forms a third complex with this metabolite. Then, the third complex up–regulates the transcription of pyocyanin. In both ways, either the first or the second complex contains PqsR and the other PqsE and one of those receptors uses an autoinducer as substrate. The unknown proteins in our model may be directly or indirectly related to the proteins MvaT and MvaU that are required to form pyocyanin [45]. It was reported that the biosynthesis of pyocyanin is a complex pathway including two operons that share high sequence homology. These are responsible for forming a precursor as well as the enzymes PhzM and PhzS which are required to convert this precursor to pyocyanin [46]. Therefore, we suggest that PqsE might be a regulatory protein and that pyocyanin biosynthesis might be regulated by the two regulators PqsR and PqsE as well as an autoinducer as ligand.
Conclusions
We have introduced a novel rule–based logical model for QS in P. aeruginosa. The model incorporates the key steps of bacterial cells leading to the synthesis of autoinducer molecules and virulence factors. Despite its simplicity, the model enables almost quantitative predictions about the effect of inhibitors of PQS biosynthesis and antagonists of their corresponding receptor PqsR. Also, one can easily investigate the effect of reviewed network topologies. We were able to identify PqsR and PqsE as a key region of the network where the predicted effects of pyocyanin production deviate from experiment. We suggest additional regulatory interactions that can be validated experimentally. In the future, this approach can serve as basis for further integrating the effect of random mutations in various parts of the network.
Declarations
Acknowledgements
We thank Dr. Tihamér Geyer and Dr. Michael Hutter for discussions and careful reading of the manuscript.
Authors’ Affiliations
References
- Rasmussen T, Givskov M, et al.: Quorum-sensing inhibitors as anti-pathogenic drugs. IJMM. 2006, 296: 149-PubMedGoogle Scholar
- Hentzer M, Wu H, Andersen J, Riedel K, Rasmussen T, Bagge N, Kumar N, Schembri M, Song Z, Kristoffersen P, et al.: Attenuation of Pseudomonas aeruginosa virulence by quorum sensing inhibitors. EMBO. 2003, 22: 3803-3815. 10.1093/emboj/cdg366.View ArticleGoogle Scholar
- De Kievit T, Iglewski B: Bacterial quorum sensing in pathogenic relationships. Infect Immun. 2000, 68: 4839-4849. 10.1128/IAI.68.9.4839-4849.2000.PubMedPubMed CentralView ArticleGoogle Scholar
- Antunes L, Ferreira R, Buckner M, Finlay B: Quorum sensing in bacterial virulence. Microbiol. 2010, 156: 2271-2282. 10.1099/mic.0.038794-0.View ArticleGoogle Scholar
- Pearson J, Passador L, Iglewski B, Greenberg E: A second N-acylhomoserine lactone signal produced by Pseudomonas aeruginosa. Proc Natl Acad Sci. 1995, 92: 1490-10.1073/pnas.92.5.1490.PubMedPubMed CentralView ArticleGoogle Scholar
- Pesci E, Milbank J, Pearson J, McKnight S, Kende A, Greenberg E, Iglewski B: Quinolone signaling in the cell-to-cell communication system of Pseudomonas aeruginosa. Proc Natl Acad Sci. 1999, 96: 11229-10.1073/pnas.96.20.11229.PubMedPubMed CentralView ArticleGoogle Scholar
- Seed P, Passador L, Iglewski B: Activation of the Pseudomonas aeruginosa lasI gene by LasR and the Pseudomonas autoinducer PAI: an autoinduction regulatory hierarchy. J Bacteriol. 1995, 177: 654-659.PubMedPubMed CentralGoogle Scholar
- De Kievit T, Seed P, Nezezon J, Passador L, Iglewski B: RsaL, a novel repressor of virulence gene expression in Pseudomonas aeruginosa. J Bacteriol. 1999, 181: 2175-2184.PubMedPubMed CentralGoogle Scholar
- Passador L, Cook J, Gambello M, Rust L, Iglewski B: Expression of Pseudomonas aeruginosa virulence genes requires cell-to-cell communication. Science. 1993, 260: 1127-1130. 10.1126/science.8493556.PubMedView ArticleGoogle Scholar
- Ochsner U, Reiser J: Autoinducer-mediated regulation of rhamnolipid biosurfactant synthesis in Pseudomonas aeruginosa. Proc Natl Acad Sci. 1995, 92: 6424-10.1073/pnas.92.14.6424.PubMedPubMed CentralView ArticleGoogle Scholar
- Pesci E, Pearson J, Seed P, Iglewski B: Regulation of las and rhl quorum sensing in Pseudomonas aeruginosa. J Bacteriol. 1997, 179: 3127-3132.PubMedPubMed CentralGoogle Scholar
- Reimmann C, Beyeler M, Latifi A, Winteler H, Foglino M, Lazdunski A, Haas D: The global activator GacA of Pseudomonas aeruginosa PAO positively controls the production of the autoinducer N-butyryl-homoserine lactone and the formation of the virulence factors pyocyanin, cyanide, and lipase. Mol Microbiol. 1997, 24: 309-319. 10.1046/j.1365-2958.1997.3291701.x.PubMedView ArticleGoogle Scholar
- Albus A, Pesci E, Runyen-Janecky L, West S, Iglewski B: Vfr controls quorum sensing in Pseudomonas aeruginosa. J Bacteriol. 1997, 179: 3928-3935.PubMedPubMed CentralGoogle Scholar
- Ochsner U, Fiechter A, Reiser J: Isolation, characterization, and expression in Escherichia coli of the Pseudomonas aeruginosa rhlAB genes encoding a rhamnosyltransferase involved in rhamnolipid biosurfactant synthesis. J Biological Chem. 1994, 269: 19787-Google Scholar
- Medina G, Juárez K, Soberón-Chávez G: The Pseudomonas aeruginosa rhlAB operon is not expressed during the logarithmic phase of growth even in the presence of its activator RhlR and the autoinducer N-butyryl-homoserine lactone. J Bacteriol. 2003, 185: 377-380. 10.1128/JB.185.1.377-380.2003.PubMedPubMed CentralView ArticleGoogle Scholar
- Maier R, Soberon-Chavez G: Pseudomonas aeruginosa rhamnolipids: biosynthesis and potential applications. Appl Microbiol Biotechnol. 2000, 54: 625-633. 10.1007/s002530000443.PubMedView ArticleGoogle Scholar
- Déziel E, Lépine F, Milot S, He J, Mindrinos M, Tompkins R, Rahme L: Analysis of Pseudomonas aeruginosa 4-hydroxy-2-alkylquinolines (HAQs) reveals a role for 4-hydroxy-2-heptylquinoline in cell-to-cell communication. Proc Natl Acad Sci. 2004, 101: 1339-1344. 10.1073/pnas.0307694100.PubMedPubMed CentralView ArticleGoogle Scholar
- Xiao G, He J, Rahme L: Mutation analysis of the Pseudomonas aeruginosa mvfR and pqsABCDE gene promoters demonstrates complex quorum-sensing circuitry. Microbiol. 2006, 152: 1679-1686. 10.1099/mic.0.28605-0.View ArticleGoogle Scholar
- Pistorius D, Ullrich A, Lucas S, Hartmann R, Kazmaier U, Müller R: Biosynthesis of 2-Alkyl-4 (1H)-Quinolones in Pseudomonas aeruginosa: Potential for therapeutic interference with pathogenicity. ChemBioChem. 2011, 12: 850-853. 10.1002/cbic.201100014.PubMedView ArticleGoogle Scholar
- Rampioni G, Pustelny C, Fletcher M, Wright V, Bruce M, Rumbaugh K, Heeb S, Cámara M, Williams P: Transcriptomic analysis reveals a global alkyl-quinolone-independent regulatory role for PqsE in facilitating the environmental adaptation of Pseudomonas aeruginosa to plant and animal hosts. Environ Microbiol. 2010, 12: 1659-1673.PubMedPubMed CentralGoogle Scholar
- Heeb S, Fletcher M, Chhabra S, Diggle S, Williams P, Cámara M: Quinolones: from antibiotics to autoinducers. FEMS Microbiol Rev. 2011, 35: 247-274. 10.1111/j.1574-6976.2010.00247.x.PubMedPubMed CentralView ArticleGoogle Scholar
- Chugani SA, Whiteley M, Lee KM, D’Argenio D, Manoil C, Greenberg EP: QscR, a modulator of quorum-sensing signal synthesis and virulence in Pseudomonas aeruginosa. PNAS. 2001, 98: 2752-2757. 10.1073/pnas.051624298.PubMedPubMed CentralView ArticleGoogle Scholar
- Liang H, Deng X, Ji Q, Sun F, Shen T, He C: The Pseudomonas aeruginosa global regulator VqsR directly inhibits QscR to control quorum-sensing and virulence gene expression. J Bacteriol. 2012, 194: 3098-3108. 10.1128/JB.06679-11.PubMedPubMed CentralView ArticleGoogle Scholar
- Atlan H, Fogelman-Soulie F, Salomon J, Weisbuch G: Random boolean networks. Cybern Syst. 1981, 12: 103-121. 10.1080/01969728108927667.View ArticleGoogle Scholar
- Lähdesmäki H, Hautaniemi S, Shmulevich I, Yli-Harja O: Relationships between probabilistic Boolean networks and dynamic Bayesian networks as models of gene regulatory networks. Signal Process. 2006, 86: 814-834. 10.1016/j.sigpro.2005.06.008.View ArticleGoogle Scholar
- Steggles L, Banks R, Wipat A: Modelling and analysing genetic networks: From Boolean networks to Petri nets. Comp Meth Systems Biol. 2006, Springer, 127-141.4210 of Lecture Notes in Computer Science,View ArticleGoogle Scholar
- Rodríguez A, Sosa D, Torres L, Molina B, Frías S, Mendoza L: A Boolean network model of the FA/BRCA pathway. Bioinformatics. 2012, 28: 858-866. 10.1093/bioinformatics/bts036.PubMedView ArticleGoogle Scholar
- Sanchez L, Thieffry D: A logical analysis of the drosophila gap-gene system. J Theo Biol. 2001, 211: 115-141. 10.1006/jtbi.2001.2335.View ArticleGoogle Scholar
- Fauré A, Naldi A, Lopez F, Chaouiya C, Ciliberto A, Thieffry D: Modular logical modelling of the budding yeast cell cycle. Mol BioSyst. 2009, 5: 1787-1796. 10.1039/b910101m.PubMedView ArticleGoogle Scholar
- Dockery J, Keener J: A mathematical model for quorum sensing in Pseudomonas aeruginosa. Bull Math Biol. 2001, 63: 95-116. 10.1006/bulm.2000.0205.PubMedView ArticleGoogle Scholar
- Melke P, Sahlin P, Levchenko A, Jönsson H: A cell-based model for quorum sensing in heterogeneous bacterial colonies. PLoS Comp Biol. 2010, 6: e1000819-10.1371/journal.pcbi.1000819.View ArticleGoogle Scholar
- Bernardini F, Gheorghe M, Krasnogor N: Quorum sensing P systems. Theor Comput Sci. 2007, 371: 20-33. 10.1016/j.tcs.2006.10.012.View ArticleGoogle Scholar
- Anguige K, King J, Ward J, Williams P: Mathematical modelling of therapies targeted at bacterial quorum sensing. Math Biosci. 2004, 192: 39-83. 10.1016/j.mbs.2004.06.008.PubMedView ArticleGoogle Scholar
- Chopp D, Kirisits M, Moran B, Parsek M: A mathematical model of quorum sensing in a growing bacterial biofilm. J Ind Microbiol Biotech. 2002, 29: 339-346. 10.1038/sj.jim.7000316.View ArticleGoogle Scholar
- Fozard J, Lees M, King J, Logan B: Inhibition of quorum sensing in a computational biofilm simulation. Biosystems. 2012, 109: 105-114. 10.1016/j.biosystems.2012.02.002.PubMedView ArticleGoogle Scholar
- Gambello M, Iglewski B: Cloning and characterization of the Pseudomonas aeruginosa lasR gene, a transcriptional activator of elastase expression. J Bacteriol. 1991, 173: 3000-3009.PubMedPubMed CentralGoogle Scholar
- Wilder C, Diggle S, Schuster M: Cooperation and cheating in Pseudomonas aeruginosa: the roles of the las, rhl and pqs quorum-sensing systems. ISME J. 2011, 5: 1332-1343. 10.1038/ismej.2011.13.PubMedPubMed CentralView ArticleGoogle Scholar
- Kesarwani M, Hazan R, He J, Que Y, Apidianakis Y, Lesic B, Xiao G, Dekimpe V, Milot S, Deziel E, et al.: A Quorum Sensing Regulated Small Volatile Molecule Reduces Acute Virulence and Promotes Chronic Infection Phenotypes. PLoS Pathog. 2011, 7: e1002192-10.1371/journal.ppat.1002192.PubMedPubMed CentralView ArticleGoogle Scholar
- Gallagher L, McKnight S, Kuznetsova M, Pesci E, Manoil C: Functions required for extracellular quinolone signaling by Pseudomonas aeruginosa. J Bacteriol. 2002, 184: 6472-6480. 10.1128/JB.184.23.6472-6480.2002.PubMedPubMed CentralView ArticleGoogle Scholar
- Farrow III J, Sund Z, Ellison M, Wade D, Coleman J, Pesci E: PqsE functions independently of PqsR-Pseudomonas quinolone signal and enhances the rhl quorum-sensing system. J Bacteriol. 2008, 190: 7043-7051. 10.1128/JB.00753-08.View ArticleGoogle Scholar
- Storz M, Maurer C, Zimmer C, Wagner N, Brengel C, de Jong J, Lucas S, Müsken M, Häussler S, Steinbach A, Hartmann R: Validation of PqsD as anti-biofilm target in Pseudomonas aeruginosa by development of small molecule inhibitors. J Am Chem Soc. 2012, 134: 16143-16146. 10.1021/ja3072397.PubMedView ArticleGoogle Scholar
- Klein T, Henn C, de Jong, Zimmer C, Kirsch B, Maurer C, Pistorius D, Müller R, Steinbach A, Hartmann R: Identification of small-molecule antagonists of the Pseudomonas aeruginosa transcriptional regulator PqsR: Biophysically guided hit discovery and optimization. ACS Chem Biol. 2012, 7: 1496-1501. 10.1021/cb300208g.PubMedView ArticleGoogle Scholar
- Lu C, Kirsch B, Zimmer C, de Jong J, Henn C, Maurer C, Müsken M, Steinbach A, Hartmann R: Discovery of antagonists of PqsR, a key player in 2-Alkyl-4-quinolone-dependent quorum sensing in Pseudomonas aeruginosa. Chem Biol. 2012, 19: 381-390. 10.1016/j.chembiol.2012.01.015.PubMedView ArticleGoogle Scholar
- Yu S, Jensen V, Seeliger J, Feldmann I, Weber S, Schleicher E, Häussler S, Blankenfeldt W: Structure elucidation and preliminary assessment of hydrolase activity of PqsE, the Pseudomonas quinolone signal (PQS) response protein. Biochem. 2009, 48: 10298-10307. 10.1021/bi900123j.View ArticleGoogle Scholar
- Li C, Wally H, Miller S, Lu C: The multifaceted proteins MvaT and MvaU, members of the H-NS family, control arginine metabolism, pyocyanin synthesis, and prophage activation in Pseudomonas aeruginosa PAO1. J Bacteriol. 2009, 191: 6211-6218. 10.1128/JB.00888-09.PubMedPubMed CentralView ArticleGoogle Scholar
- Mavrodi D, Bonsall R, Delaney S, Soule M, Phillips G, Thomashow L: Functional analysis of genes for biosynthesis of pyocyanin and phenazine-1-carboxamide from Pseudomonas aeruginosa PAO1. J Bacteriol. 2001, 183: 6454-6465. 10.1128/JB.183.21.6454-6465.2001.PubMedPubMed CentralView ArticleGoogle Scholar
Copyright
This article is published under license to 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.