Skip to main content

Genetic analysis, structural modeling, and direct coupling analysis suggest a mechanism for phosphate signaling in Escherichia coli



Proper phosphate signaling is essential for robust growth of Escherichia coli and many other bacteria. The phosphate signal is mediated by a classic two component signal system composed of PhoR and PhoB. The PhoR histidine kinase is responsible for phosphorylating/dephosphorylating the response regulator, PhoB, which controls the expression of genes that aid growth in low phosphate conditions. The mechanism by which PhoR receives a signal of environmental phosphate levels has remained elusive. A transporter complex composed of the PstS, PstC, PstA, and PstB proteins as well as a negative regulator, PhoU, have been implicated in signaling environmental phosphate to PhoR.


This work confirms that PhoU and the PstSCAB complex are necessary for proper signaling of high environmental phosphate. Also, we identify residues important in PhoU/PhoR interaction with genetic analysis. Using protein modeling and docking methods, we show an interaction model that points to a potential mechanism for PhoU mediated signaling to PhoR to modify its activity. This model is tested with direct coupling analysis.


These bioinformatics tools, in combination with genetic and biochemical analysis, help to identify and test a model for phosphate signaling and may be applicable to several other systems.


Adapting to changes in the environment is one of the hallmarks of life. For all life, phosphate is an essential nutrient. Bacteria have several mechanisms to scavenge phosphate that are only expressed when the level of available environmental phosphate is limited: including a phosphate specific ABC transporter complex (PstSCAB) and a periplasmic phosphate scavenging enzyme (alkaline phosphatase (AP); the product of the phoA gene) [1]. Expression control of these genes is essential for optimal growth and has been implicated in the regulation of pathogenesis in several organisms [2, 3].

In Escherichia coli, a classic two-component signal transduction system, composed of the PhoR histidine kinase and the PhoB response regulator, is responsible for expression control of a group of genes called the Pho regulon. PhoR consists of an N-terminal membrane-spanning region, as well as cytoplasmic PAS, DHp and CA domains. The PAS domain was named for the Drosophila Per, Arnt, and Sim proteins, in which this domain was originally described and has been found in many signaling proteins. Many PAS domains bind cofactors such as heme [4]. The DHp domain is conserved in histidine kinases and functions in dimerization and contains the site of histidine phosphorylation. The CA domain is the catalytic, ATP-binding part of the protein. PhoB consists of an N-terminal phosphorylation domain that receives a phosphoryl group from PhoR and a C-terminal DNA binding domain. In low phosphate conditions, PhoR acts as a PhoB kinase. Upon phosphorylation, PhoB recruits RNA polymerase to promoters of the Pho regulon that contain a Pho box. In high phosphate conditions, PhoR acts as a phospho-PhoB phosphatase and removes the phosphate from PhoB to keep the expression level of Pho regulon genes very low. One unanswered question with this system is how PhoR perceives external phosphate concentrations. PhoR lacks a significant periplasmic domain that could detect phosphate abundance outside the cell. Past work has shown that the PstSCAB transporter and the PhoU protein play important roles in phosphate signaling to PhoR [1]. The mechanism of this signal has not been fully elucidated.

A deletion mutation of the phoU gene leads to poor growth and the frequent development of compensatory mutations in the other Pho regulon expression control genes, pstSCAB, phoR, and phoB [5]. The poor growth phenotype is likely due to overexpression and under-regulation of a functional PstSCAB transporter [6], which leads to phosphate poisoning when cells are grown in high phosphate environments [7]. Reference [7] proposed that PhoU modulates phosphate transport through the PstSCAB complex by inhibiting transport when internal phosphate levels are too high.

Recently, PhoU was shown to directly interact with PstB and PhoR [8]. This observation suggested a model that PhoU interacts with PstB to sense environmental phosphate levels and that it passes that signal along to PhoR to modulate its kinase/phosphatase activities (Figure 1). To further characterize these interactions, this work analyzes several phoU mutants for signaling activity, interactions with PhoR, and interactions with PstB. A scanning mutagenesis screen of the PAS domain of PhoR identified potential residues important for interaction with PhoU. We modeled potential docking interactions between PhoR and PhoU. One of the docking models showed a close proximity of identified residues in PhoU with a predicted interaction loop of the PAS domain of PhoR. To validate this model, we performed a Direct Coupling Analysis (DCA) between PhoR and PhoU using sequences from the gammaproteobacteria. The DCA results are consistent with the proposed docking model and may point to a mechanism of action for PhoU in controlling the opposing kinase and phosphatase activities of PhoR.

Figure 1
figure 1

Pho regulon expression control. Proper control of Pho regulon expression involves the histidine kinase, PhoR that has a Per-Arnt-Sim (PAS), dimerization/histidine phosphorylation (DHp), and catalytic (CA) domains; the response regulator PhoB (with a conserved aspartate (D) residue that is phosphorylated); the phosphate specific ABC transporter composed of PstS (Periplasmic phosphate binding protein), PstC and PstA (proteins that form the pore in the inner membrane), and PstB (the ATP binding portion of the transporter); and PhoU (a negative regulator of PhoR).


Strains, plasmids, and reconstructing of signaling system

Plasmids that were used include pKG116[9], p116U2[7] (a pKG116 with the phoU gene under salicylate expression control), p116A147E, p116U2 A147K, p116U2 R148A, p116U2 R148E, p116U2 D150A (all the mutant p116U2 plasmids were constructed with QuikChange site-directed mutagenesis kit from Agilent Technologies and verified by DNA sequence analysis as described in [8] (primers listed in Table S1 found in additional file 1), pRR48[10], and p48SCAB (a pRR48 with the pstSCAB genes inserted also as described in [8] (primers are listed in Table S1 found in additional file 1). For signaling analysis, combinations of p48SCAB and a pKG116 derived plasmid were introduced into the E. coli strain BW26337 (having a chromosomal deletion of pstSCABphoU constructed from the wildtype strain BW25113 [11]). The cultures were grown at 37°C with shaking overnight in morpholinepropanesulfonic acid (MOPS) defined medium with 0.06% glucose and 2.0 mM phosphate (a high phosphate minimal media) [8].

Growth of PhoU A147E mutant

STAC and STACΔphoU [8] were used with pKG116, p116U2, or p116U2 A147E (a p116U2 plasmid mutated with QuikChange site-directed mutagenesis kit from Agilent Technologies and verified by DNA sequence analysis (primers listed in Table S1 found in additional file 1)). Triplicate overnight cultures were grown in LB media with 100 µM IPTG and the OD600 of a 1:5 dilution of culture was measured. Values were adjusted for dilution and averages are reported with error bars representing the standard deviation.

Alkaline phosphatase assays

We based our Alkaline Phosphatase (AP) Assays on a 96 well plate β-Galactosidase Activity assay previously described [8]. Cultures were grown overnight at 37°C with shaking. The OD600 was read on a 1:4 dilution of the culture. 1 ml of culture was collected and resuspended in 900 µl of 1 M Tris-HCl pH 8.2. 1 drop of 0.1% sodium dodecyl sulfate and 2 drops of chloroform were added to each tube and tubes were vortexed vigorously for 15 sec. Tubes were then spun for 1 min at 16,000 × g and 200 µl of each sample was loaded into a well of a 96 well flat bottomed plate. 40 µl of 20 mM p-Nitrophenyl phosphate in 1 M Tris-HCl pH 8.2 was added. The OD420 values were read once a minute for 20 min at 37°C. Units of activity were calculated as (1000 × slope of a line fit to OD420 in mOD/min)/(4 × OD600 of 1:4 dilution of the overnight culture). Cultures were grown in triplicate and the average is reported with error bars representing the standard deviation.

BACTH and β-galactosidase assays for scanning mutagenesis

BACTH analysis and β-Galactosidase Assays were performed as described in [8]. Briefly, using the QuikChange Lightning Site-Directed Mutagenesis Kit from Agilent Technologies, the p18CRN-PAS plasmid was mutated to change every two amino acids of the PhoR PAS domain to code for alanine-cysteine (primers listed in Table S1 found in additional file 1). For example, p18CRN-PAS 109 had PhoR D109A and A110C changes. These residues were chosen because they are not predicted to cause major secondary structural changes. Each p18CRN-PAS plasmid of the library of alanine-cysteine mutants was transformed with pKT25phoU into the tester strain, BTH101 from EuroMedex. Cultures were grown in triplicate in LB and assayed as described in [8]. The percent of PhoU/PhoR Interaction was found by dividing the average activity of each sample with the activity of an unmutated p18CRN-PAS control and multiplied by 100.

Protein structure modeling and protein docking modeling

We used the ClusPro webserver [1215] to dock structures of a dimer of the cytoplasmic portion from E. coli PhoR (structure modeled from the structure for VicK from Streptococcus mutans [16]), and PhoU (modeled from Thermatoga maratima PhoU [17]) using Protein Homology/analogY Recognition Engine V 2.0 (Phyre2) [18].

Direct coupling analysis

We used PhoU and PhoR sequences from 196 species of bacteria from the gammaproteobacteria group. Sequences for proteins were collected using the Kegg webserver ( We used a list of PhoR and PhoB orthologs to identify species where both annotated PhoR PhoB encoding genes were found on the chromosome. The PhoR and PhoU sequences were collected manually using the predicated protein domains and the genomic context of the gene to select genes that were most likely correctly annotated. For example, for this analysis we were not interested in histidine kinases annotated as PhoR from species that did not have a phoU gene on the chromosome. Sequences were concatenated starting with PhoU, followed by PhoR from the same species. Twenty alanine residues were artificially placed between the two protein sequences to aid in differentiating the PhoU from the PhoR after the alignment. Using MAFFT version 7 (website is, we were able to create an alignment file to use in DCA analysis. A MATLAB implementation of direct coupling analysis reported in [19] distinguished differences between the direct information and mutual information between all protein residues. The command line argument used was matlab-nodisplay-nojvm-nosplash-r "dca $INPUT_ALIGNED_FASTA $OUT_FILE". The DCA implementation used the input aligned fasta file in its calculations, outputting a file with N(N-1)/2 (N = length of the sequences) rows and four columns: residue i (column 1), residue j (column 2), MI(i,j) (Mutual Information between i and j), and DI(i,j) (Direct Information between i and j). All inserts columns were removed from the alignment during the analysis. Concluding this analysis, a simple python script was used to remove all rows in which residues i and j were within five residue numbers of each other. This was done to eliminate false positive residue interactions caused by their proximity to other amino acids. Mutual information was later screened to signify relatedness between two residues.

Results and discussion

Mutations in any of the pstSCAB transporter genes or phoU lead to loss of phosphate signaling. To determine whether PhoU could act independently of the PstSCAB transporter, we expressed these two genetic entities from separate plasmids in E. coli strain BW26337, which contains a deletion of the pstSCABphoU operon and tested strains for control of the Pho regulon. We cloned the pstSCAB genes into the pRR48 plasmid (AmpR) [10] and the phoU gene into the compatible pKG116 plasmid (CamR) [9]. Control of the Pho regulon was analyzed by determining Alkaline Phosphatase (AP) expression. Figure 2 shows that AP expression was unregulated when cells were grown under high phosphate conditions in strains expressing neither protein (containing two empty vectors), or expressing each protein individually. However, when both proteins were expressed in the same cell, AP expression levels were significantly reduced. These results show that phoU expression is necessary but not sufficient for proper phosphate signaling and are consistent with a model that these proteins act together in signal transduction.

Figure 2
figure 2

Signaling necessity and sufficiency of PhoU and PstSCAB. Triplicate cultures of BW26337 (a ΔpstSCABphoU strain) cells with pKG116 and pRR48 (Empty vectors), pKG116 and pRR48SCAB (PstSCAB), p116U2 and pRR48 (PhoU), p116U2 and pRR48SCAB (PstSCAB PhoU) and BW25113 (a wild-type strain) were grown in LB and Bacterial Alkaline Phosphatase activity was assayed. Error bars represent the standard deviation. * = p-value of < 0.05 compared to the Empty vectors with a two-tailed T-test.

Early studies of phosphate homeostasis in E. coli isolated mutants that constitutively expressed AP [20]. One mutant, named C4 [21] was characterized [22] and was later named phoU35 because it was independent of the phosphate transport genes [23]. When the genes of the pstSCABphoU operon were sequenced and the phoU35 mutant was analyzed, they found this mutation encoded a change from alanine at position 147 to glutamic acid (A147E) [24]. Since a phoU deletion mutation results in loss of phosphate signaling and causes a severe growth phenotype, but the phoU35 allele was only reported to cause a loss of signaling without the accompanying poor growth, we hypothesized that the phoU35 mutation may disrupt PhoU's interaction with PhoR, preventing the signal for the switch to PhoR phosphatase activity, but may maintain the interaction with PstB, limiting excess transport of phosphate into the cell during phosphate replete conditions.

We wanted to confirm that strains expressing the phoU35 allele do not show a severe growth phenotype. We employed a phoU deletion strain in which the normal pstS promoter was replaced by the tac promoter to uncouple expression of this operon from mutations in any of its genes (STACΔphoU) [8]. By growing this strain in the absence of IPTG the accumulation of compensatory mutations was avoided. Figure 3 shows the results of characterizing the growth yield of various cultures grown in a high-phosphate medium with IPTG overnight with shaking. The STAC strain with a wildtype copy of phoU still in the chromosome grew normally. We observed a significant growth defect when the phoU knockout strain containing an empty vector, pKG116, was grown in identical conditions. When the p116U2 plasmid was introduced into this strain, expressing wild-type phoU, and grown under identical conditions, no growth defect was observed. Moreover, when the phoUA147E allele was introduced into the strain, the growth phenotype was similar to that of wild-type phoU. These results confirm the lack of growth defect in the phoU35 mutant and are consistent with the hypothesis that the A147 residue of PhoU is important for interactions with PhoR.

Figure 3
figure 3

Growth yield of overnight cultures. Triplicate cultures were grown overnight in LB and the absorbance of 600nm wavelength of light for each culture was recorded. Error bars represent standard deviations. * = p-value of < 0.01 compared to STAC with a two-tailed T-test.

To verify that the phoU35 allele caused a signaling defect when expressed from a multicopy plasmid, we tested AP expression in the strains constructed for the previous experiments. Figure 4 shows that high AP levels were observed in the STACΔphoU pKG116 strain and the STACΔphoU p116U2 A47E strain grown under high phosphate conditions, but that reduced AP levels were observed in the STACΔphoU p116U2 strain. This showed that PhoUA147E could not regulate PhoR and was consistent with the model presented. We wondered if other mutations that caused changes in the phoU protein in the vicinity of A147 would also display a signaling defect. We therefore introduced additional mutations by site-directed mutagenesis in the region near A147 (A147, R148, and D150) to create versions of PhoU with different charges and sizes of amino acid side chains in this region. The mutations that we created were A147K, R148E, R148A, and D150A. The rightmost four bars in Figure 4 show that the A147K, R148E and R148A mutants lost phosphate signaling activity as they expressed elevated AP levels, but that the D150A mutation retained signaling activity. These results suggest that both A147 and R148 are important for interactions with PhoR.

Figure 4
figure 4

Signaling activity of various PhoU mutants. Triplicate cultures were grown in LB and assayed for Bacterial Alkaline Phosphatase activity. BW25113 is a wildtype control. Other strains were STACΔphoU with p116U2 (wildtype phoU), pKG116 (empty parent plasmid), and various p116U2 plasmids with specific mutations in phoU (A147E, A147K, R148A, R148E, and D150A). Error bars represent standard deviations. Values significantly greater than p116U2 are labeled with * = p-value of < 0.001 compared to p116U2 with a two-tailed T-test.

The combination of no growth inhibition and disrupted signaling are the expected phenotypes of a phoU mutant that can interact with PstB to limit phosphate transport but can no longer interact with PhoR to signal the phosphate level.

The combination of no growth inhibition and disrupted signaling are the expected phenotypes of a phoU mutant that can interact with PstB to limit phosphate transport but can no longer interact with PhoR to signal the phosphate level. We further examined this hypothesis by using a bacterial adenylate cyclase two-hybrid (BACTH) system to test for interactions between mutant versions of PhoU and PstB and PhoR (Figure 5A and 5B).

Figure 5
figure 5

PhoU mutant interactions with PhoR and PstB. Values shown are the average of triplicate samples with the error bars representing standard deviation. A. BACTH analysis of PhoU/PhoR interaction. Various combinations of PhoU and mutants or an empty vector negative control were tested for interaction with PhoR. Activity of β-Galactosidase correlates with protein/protein interaction. * = p-value < 0.01 compared to PhoU/PhoR with a two-tailed T-test. B. This chart shows the interaction of PstB with various PhoU mutants and an empty vector negative control. * = p < 0.05 compared to Empty Vector with a two-tailed T-test.

The BACTH system employs the separable T18 and T25 domains of adenylate cyclase from Bordetella pertussis. When these two domains are in close proximity they create an active enzyme that produces cAMP. By creating gene fusions in which protein domains are connected to the T18 and T25 fragments, cAMP production is a measure of whether the fused proteins interact. Since cAMP binds to the CRP protein, cAMP can be measured indirectly by assaying β-galactosidase. We previously used this method to show interactions between the wild-type version of PhoU and PstB and various domains of PhoR [8].

Figure 5A shows the results of β-galactosidase activity assays of various BACTH strains. Each sample is BTH101 strain containing one plasmid that expresses the T18 domain fused to PhoR and another plasmid that expresses either the T25 domain alone (Empty vector), T25 fused to PhoU (PhoU/PhoR), or various PhoU mutant proteins (A147E, A147K, A148A, R148E, and D150A). The A147 and R148 mutants of phoU had a significantly weaker interaction with PhoR (as represented by low β-galactosidase activities) and the D150A mutant retained PhoR interaction. We also tested these plasmids expressing the T25-PhoU fusions for interaction with a T18-PstB fusion protein. Interestingly, the A147 and R148 mutant proteins maintained interactions with PstB (Figure 5B). This implies that the loss of interaction with PhoR is not due to protein instability or radical protein misfolding.

Previous work showed that the PhoR/PhoU interaction was dependent on the PAS domain of PhoR [8]. We used scanning mutagenesis of a plasmid expressing a T18-PhoR N-PAS (the portion of the PhoR protein from the C terminus through the PAS domain, but without the CA and DHp domains) domain fusion protein to identify the residues within this domain that are important for the interaction. Every two amino acids of the PhoR PAS domain were mutated to alanine and cysteine residues by sequentially replacing the six bases encoding adjacent codons with a SphI restriction site. The mutant versions were then introduced into the appropriate tester strain expressing wild-type phoU fused to the T25 fragment of adenylate cyclase and β-galactosidase assays were performed. Many mutants lost the ability to interact with PhoU as indicated by significantly reduced β-galactosidase levels (shown as blue bars in Figure 6B). Among the mutations there were several regions where neighbouring residues showed loss of interaction (for example, residues 141-146, 157-162 and 169-176). There were also several single mutations (residues 111, 115, 189) that reduced protein interactions.

Figure 6
figure 6

PAS domain scanning mutagenesis. A. Predicted secondary structure of the PhoR PAS domain with alpha helixes labeled with yellow arrows and β-sheets labeled with red cylinders. B. Scanning mutagenesis of the PhoR PAS domain used BACTH to identify regions of the protein that are essential for interaction with PhoU. Every two amino acids were changed to code for alanine and cysteine. Each construct was tested in triplicate for interaction with PhoU. Blue bars represent samples that had less than 40% of the activity of unmutated PhoR.

We modeled the cytoplasmic portion of the E. coli PhoR structure based upon the structure of the VicK protein from Streptococcus mutans using the Phyre2 webserver (Figure 7A). This modeled structure predicts that the 141-146 and the 157-162 regions of PhoR have residues that are exposed to the surface and face outwards, consistent with an interaction with another protein (Figure 7A, the residues colored red and purple respectively.) The 111, 115, 189, and the 165-176 mutations all map to regions of PhoR that are predicted to lie along a β-sheet that has been shown in other histidine kinases to be important for PAS domain dimerization contacts. For this reason, these mutations may not be involved in PhoU/PhoR interactions, but for maintaining a proper protein conformation.

Figure 7
figure 7

Structures of modeled PhoR and PhoU. A. Interaction sites identified from scanning mutagenesis that are found on the surface of a modeled PhoR structure, 141-146 are red and the 157-162 residues are purple. B. Shows the protein surface predictions of the PhoU/PhoR docking model. PhoR is colored yellow and PhoU is colored red. C. A cartoon of protein backbones with PhoU A147 and R148 shown as spheres in blue and PhoR K157-S162 shown as green spheres.

We docked a modeled E. coli PhoU protein onto the modeled PhoR dimer. The server returned ten potential PhoR/PhoU models with most of the models being unreasonable because they docked PhoU onto the part of PhoR that interacts with the cytoplasmic membrane. One of the reasonable models showed an interesting interaction between PhoR and PhoU (Figure 7B). In this model the PhoU residues that were identified through genetic analysis as important for interactions with PhoR were in close proximity to regions of the PhoR PAS domain that were identified by scanning mutagenesis. Figure 7C shows that some of the PhoR 157-162 residues (shown in green) appear to interact with the A147 and R148 residues of PhoU (shown in blue).

Support for this interaction model was obtained through a bioinformatic method, called direct coupling analysis (DCA) that identifies covariance between residues. DCA identifies which residues from a sequence tend to co-evolve with any other residue by measuring the predictive power of one residue on another. This can identify direct interactions between residues, such as residues that are directly involved in protein/protein interactions. Also, residues that are involved in indirect interaction (for example residues that play a role in proper structural conformation) and mechanistic interactions (as would be found in residues associated with an enzyme active site) are identified with DCA [25]. Mutual information (MI), which identifies direct and indirect coupling of residue pairs, has been used to identify both intradomain and interdomain interactions [2632]. Direct information (DI) is a value that attempts to eliminate the indirect coupling of MI due to neighbouring residues. Biologically, MI may identify residues that play a role in interaction through indirect means that may not be identified with DI. We used MI to sort our DCA results to identify the residues that play an important role in PhoU/PhoR interaction,

We analyzed sequences of PhoU and PhoR from 196 gammaproteobacteria. We concatenated the phoU and phoR sequences from each species with 20 ala residues in between the two genes to aid in alignment. We aligned the sequences and then compared all the residues with DCA.

Our previous work identified PhoU A147 and R148 residues as being involved in interaction with PhoR. To isolate residues from the PhoR PAS domain that show covariance with PhoU A147 and R148, we sorted the top co-varying residues of the PhoR PAS domain by MI for both A147 and R148 positions of PhoU. The top ten hits are shown in Table 1 and Table 2. Several of these co-varying residues fall in the same region of the PAS domain of PhoR identified by the BACTH analysis. When the top five residues for A147 and R148 are highlighted in the docking model, we see that many cluster near the predicted interaction surface. Given that some of the residues for A147 overlap with the R148 residues we colored residues that co-vary with A147 green, residues that co-vary with R148 cyan, and residues that co-vary with both purple (Figure 8A).

Table 1 DCA of PhoU A147 vs.PhoR PAS domain.
Table 2 DCA of PhoU R148 vs.PhoR PAS domain.
Figure 8
figure 8

Residues identified from DCA. A. PhoR in yellow and PhoU in red with A147 and R148 highlighted in blue and the top five DCA residues of PhoR PAS domain highlighted as spheres (residues that co-vary with A147 are green, R148 are cyan, and residues that co-vary with both are purple). B. PhoU E121 highlighted in blue and the top seven residues in the DCA analysis of the PhoR CA domain highlighted in green spheres. C. PhoU R148 potential interactions with the side chain of PhoR PAS R163, P164, and N166 residues (distances labeled are Å).

With the docking model it appears that the loops on the opposite side of PhoU from the PAS interacting surface are in close proximity to the CA domain of the opposite PhoR of the dimer. To observe any covariance between PhoU and PhoR at these surfaces, we compared the PhoU E121 residue to the CA domain with DCA. We see that the top PhoR CA domain residues (shown as green spheres) that co-vary with E121(shown as blue spheres) cluster at the surface of PhoR nearest to PhoU E121 (Table 3 Figure 8B).

Table 3 DCA of PhoU E121 vs.PhoR CA domain

For a more complete analysis, we also sorted the potential interacting pairs by DI and found that several of the same residues were found in the top pairs (Table 4 and Table 5). In the tables, residues that were also found in the top ten residues of MI sorted data are bolded. Several other residues sorted by MI show relatively large DI values (Table 1 and Table 2). Also, it is interesting that many of the interacting pairs sorted by DI fall into the loop regions identified by the PAS domain scanning mutagenesis (Figure 6).

Table 4 DCA of PhoU A147 sorted by Direct Information
Table 5 DCA of PhoU R148 sorted by Direct Information

We noticed that there are few residues that appear to interact based on the docking model and DCA that did not show loss of function in our scanning mutagenesis (~PhoR163-166). To better understand these results we looked closer at the residues in our docking model structures. When we highlight the side chain of the PhoU R148 residue and the side chains of residues in positions 163, 164, and 166 of the PhoR PAS domain, it appears that the R148 may interact with the amino acid backbone for these sites and not the side chain residues as they appear to point away from the PhoU structure (Figure 8C). Changing the side chains of these amino acids in our scanning mutagenesis of the PAS domain may not disrupt PhoU/PhoR interactions, which explains our scanning mutagenesis results.

Additionally, when distances between predicted interacting pairs is measured, several residues are found near R148 (within ~8Å) but few are found that near to A147. One explanation for this is that PhoU A147 itself does not directly interact with PhoR, but mutations of PhoU A147 may affect the ability of R148 to interact with PhoR. For example, the A147E mutation places a large acidic side chain right next to the R148 basic side chain and may disrupt proper PhoU/PhoR interaction and phosphate signaling. Using MI to sort our DCA results allows for identifying residues like this that are important for proper interaction but may not be directly involved. This phenomenon may also explain some of the PhoR PAS residues identified with DCA, like PhoR G118. The PhoR G118 residue is not found very near to PhoU R148 or A147 in our docking model (Table 1 and Table 2). However, G118 is on the surface of the PhoR PAS domain on a loop between α-helices that appears to form the top of the PhoU binding pocket. One would expect that mutations of a highly flexible glycine at this position could be associated with compensating mutations of PhoU in the R148 region.


Our results confirm that PhoU is necessary for proper signaling in high phosphate, but that the PstSCAB transporter is also necessary. We have identified the A148 and R149 residues of PhoU as being important for interaction with PhoR. Using scanning mutagenesis, we identified residues of the PhoR PAS domain essential for interaction with PhoU. A docking model was identified and tested with DCA. The docking model shows PhoU interacting with both the PAS domain and the CA domain of PhoR, pointing to a potential mechanism for PhoU mediated switching of PhoR's kinase to phosphatase activity.

The switch from kinase to phosphatase activity for histidine kinases is dependent upon the conformation between the CA and DHp domains [33, 34]. The CA domain must remain flexible for proper kinase function. Also, the DHp domain was previously shown biochemically to be the portion of PhoR that has phospho-PhoB phosphatase activity [35]. It is possible that when PhoU interacts with PhoR, that interaction may constrain the CA domain to inhibit kinase activity and expose the DHp domain to allow phosphatase activity of PhoR. These results are the first to point to a specific molecular mechanism for PhoU mediated modification of PhoR activity.

These studies have combined genetic, mutagenesis, computer modeling, and DCA analyses in studying the molecular interaction of PhoU and PhoR. It will be interesting to apply these techniques to identify residues involved in other signaling pathways and develop interaction models. For example, PhoU interaction with PstB appears to be weaker than PhoU interaction with PhoR based on our BACTH (Figure 6A and 6B), making it more difficult to directly identify residues essential for interaction. Using modeling and DCA, potential interacting residues may be identified. Efforts to isolate a complete signaling complex of PhoR, PhoU, and PstSCAB have been unsuccessful in the past. However, if the sites of interaction were identified, perhaps an entire signaling complex could be isolated and characterized using methods directed toward the specific interaction sites.


  1. Wanner BL: Phosphorus assimilation and control of the phosphate regulon. Escherichia coli and Salmonella: cellular and molecular biology. Edited by: Neidhardt FC, III RC, Ingraham JL, Lin ECC, Low KB, Magasanik B, Reznikoff WS, Riley M, Schaechter M, Umbarger HE. 1996, Washington D. C.: American Society for Microbiology, 1357-1381.

    Google Scholar 

  2. Chekabab SM, Harel J, Dozois CM: Interplay between genetic regulation of phosphate homeostasis and bacterial virulence. Virulence. 2014, 5 (5):

  3. Lamarche MG, Wanner BL, Crepin S, Harel J: The phosphate regulon and bacterial virulence: a regulatory network connecting phosphate homeostasis and pathogenesis. FEMS Microbiol Rev. 2008, 32 (3): 461-473. 10.1111/j.1574-6976.2008.00101.x.

    Article  CAS  PubMed  Google Scholar 

  4. Taylor BL, Zhulin IB: PAS domains: internal sensors of oxygen, redox potential, and light. Microbiology and molecular biology reviews : MMBR. 1999, 63 (2): 479-506.

    PubMed Central  CAS  PubMed  Google Scholar 

  5. Steed PM, Wanner BL: Use of the rep technique for allele replacement to construct mutants with deletions of the pstSCAB-phoU operon: evidence of a new role for the PhoU protein in the phosphate regulon. J Bacteriol. 1993, 175 (21): 6797-6809.

    PubMed Central  CAS  PubMed  Google Scholar 

  6. Haldimann A, Daniels LL, Wanner BL: Use of new methods for construction of tightly regulated arabinose and rhamnose promoter fusions in studies of the Escherichia coli phosphate regulon. J Bacteriol. 1998, 180 (5): 1277-1286.

    PubMed Central  CAS  PubMed  Google Scholar 

  7. Rice CD, Pollard JE, Lewis ZT, McCleary WR: Employment of a promoter-swapping technique shows that PhoU modulates the activity of the PstSCAB2 ABC transporter in Escherichia coli. Appl Environ Microbiol. 2009, 75 (3): 573-582. 10.1128/AEM.01046-08.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Gardner SG, Johns KD, Tanner R, McCleary WR: The PhoU protein from Escherichia coli interacts with PhoR, PstB, and metals to form a phosphate-signaling complex at the membrane. J Bacteriol. 2014, 196 (9): 1741-1752. 10.1128/JB.00029-14.

    Article  PubMed Central  PubMed  Google Scholar 

  9. Buron-Barral MC, Gosink KK, Parkinson JS: Loss- and gain-of-function mutations in the F1-HAMP region of the Escherichia coli aerotaxis transducer Aer. J Bacteriol. 2006, 188 (10): 3477-3486. 10.1128/JB.188.10.3477-3486.2006.

    Article  CAS  PubMed  Google Scholar 

  10. Studdert CA, Parkinson JS: Insights into the organization and dynamics of bacterial chemoreceptor clusters through in vivo crosslinking studies. Proceedings of the National Academy of Sciences of the United States of America. 2005, 102 (43): 15623-15628. 10.1073/pnas.0506040102.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  11. Datsenko KA, Wanner BL: One-step inactivation of chromosomal genes in Escherichia coli K-12 using PCR products. Proc Natl Acad Sci USA. 2000, 97: 6640-6645. 10.1073/pnas.120163297.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  12. Kozakov D, Brenke R, Comeau SR, Vajda S: PIPER: an FFT-based protein docking program with pairwise potentials. Proteins. 2006, 65 (2): 392-406. 10.1002/prot.21117.

    Article  CAS  PubMed  Google Scholar 

  13. Comeau SR, Gatchell DW, Vajda S, Camacho CJ: ClusPro: an automated docking and discrimination method for the prediction of protein complexes. Bioinformatics. 2004, 20 (1): 45-50. 10.1093/bioinformatics/btg371.

    Article  CAS  PubMed  Google Scholar 

  14. Comeau SR, Gatchell DW, Vajda S, Camacho CJ: ClusPro: a fully automated algorithm for protein-protein docking. Nucleic acids research. 2004, 32 (Web Server): W96-99. 10.1093/nar/gkh354.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Kozakov D, Beglov D, Bohnuud T, Mottarella S, Xia B, Hall DR, Vajda S: How good is automated protein docking?. Proteins: Structure, Function, and Bioinformatics. 2013, 81 (12): 2159-2166. 10.1002/prot.24403.

    Article  CAS  Google Scholar 

  16. Wang C, Sang J, Wang J, Su M, Downey JS, Wu Q, Wang S, Cai Y, Xu X, Wu J, et al: Mechanistic insights revealed by the crystal structure of a histidine kinase with signal transducer and sensor domains. PLoS biology. 2013, 11 (2): e1001493-10.1371/journal.pbio.1001493.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. Madej T, Addess KJ, Fong JH, Geer LY, Geer RC, Lanczycki CJ, Liu C, Lu S, Marchler-Bauer A, Panchenko AR, et al: MMDB: 3D structures and macromolecular interactions. Nucleic acids research. 2012, 40 (Database): D461-464.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  18. Kelley LA, Sternberg MJ: Protein structure prediction on the Web: a case study using the Phyre server. Nature protocols. 2009, 4 (3): 363-371. 10.1038/nprot.2009.2.

    Article  CAS  PubMed  Google Scholar 

  19. Morcos F, Pagnani A, Lunt B, Bertolino A, Marks DS, Sander C, Zecchina R, Onuchic JN, Hwa T, Weigt M: Direct-coupling analysis of residue coevolution captures native contacts across many protein families. Proceedings of the National Academy of Sciences of the United States of America. 2011, 108 (49): E1293-1301. 10.1073/pnas.1111471108.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  20. Torriani A, Rothman F: Mutants of Escherichia coli constitutive for alkaline phosphatase. J Bacteriol. 1961, 81: 835-836.

    PubMed Central  CAS  PubMed  Google Scholar 

  21. Garen A, Otsuji N: Isolation of a Protein Specified by a Regulator Gene. J Mol Biol. 1964, 8: 841-852. 10.1016/S0022-2836(64)80165-0.

    Article  CAS  PubMed  Google Scholar 

  22. Zuckier G, Torriani A: Genetic and physiological tests of three phosphate-specific transport mutants of Escherichia coli. J Bacteriol. 1981, 145 (3): 1249-1256.

    PubMed Central  CAS  PubMed  Google Scholar 

  23. Amemura M, Shinagawa H, Makino K, Otsuji N, Nakata A: Cloning of and complementation tests with alkaline phosphatase regulatory genes (phoS and phoT) of Escherichia coli. J Bacteriol. 1982, 152 (2): 692-701.

    PubMed Central  CAS  PubMed  Google Scholar 

  24. Muda M, Rao NN, Torriani A: Role of PhoU in phosphate transport and alkaline phosphatase regulation. J Bacteriol. 1992, 174 (24): 8057-8064.

    PubMed Central  CAS  PubMed  Google Scholar 

  25. Dwyer RS, Ricci DP, Colwell LJ, Silhavy TJ, Wingreen NS: Predicting functionally informative mutations in Escherichia coli BamA using evolutionary covariance analysis. Genetics. 2013, 195 (2): 443-455. 10.1534/genetics.113.155861.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  26. Liu Y, Eyal E, Bahar I: Analysis of correlated mutations in HIV-1 protease using spectral clustering. Bioinformatics. 2008, 24 (10): 1243-1250. 10.1093/bioinformatics/btn110.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  27. Liu Y, Bahar I: Sequence evolution correlates with structural dynamics. Mol Biol Evol. 2012, 29 (9): 2253-2263. 10.1093/molbev/mss097.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  28. Clark GW, Ackerman SH, Tillier ER, Gatti DL: Multidimensional mutual information methods for the analysis of covariation in multiple sequence alignments. BMC bioinformatics. 2014, 15: 157-10.1186/1471-2105-15-157.

    Article  PubMed Central  PubMed  Google Scholar 

  29. Simonetti FL, Teppa E, Chernomoretz A, Nielsen M, Marino Buslje C: MISTIC: Mutual information server to infer coevolution. Nucleic acids research. 2013, 41 (Web Server): W8-14.

    Article  PubMed Central  PubMed  Google Scholar 

  30. Mahony S, Auron PE, Benos PV: Inferring protein-DNA dependencies using motif alignments and mutual information. Bioinformatics. 2007, 23 (13): i297-304. 10.1093/bioinformatics/btm215.

    Article  CAS  PubMed  Google Scholar 

  31. Weil P, Hoffgaard F, Hamacher K: Estimating sufficient statistics in co-evolutionary analysis by mutual information. Computational biology and chemistry. 2009, 33 (6): 440-444. 10.1016/j.compbiolchem.2009.10.003.

    Article  CAS  PubMed  Google Scholar 

  32. Zhang S, McDonald PW, Thompson TA, Dennis AF, Akopov A, Kirkness EF, Patton JT, McDonald SM: Analysis of human rotaviruses from a single location over an 18-year time span suggests that protein coadaption influences gene constellations. Journal of virology. 2014, 88 (17): 9842-9863. 10.1128/JVI.01562-14.

    Article  PubMed Central  PubMed  Google Scholar 

  33. Dago AE, Schug A, Procaccini A, Hoch JA, Weigt M, Szurmant H: Structural basis of histidine kinase autophosphorylation deduced by integrating genomics, molecular dynamics, and mutagenesis. Proceedings of the National Academy of Sciences of the United States of America. 2012, 109 (26): E1733-1742. 10.1073/pnas.1201301109.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  34. Albanesi D, Martin M, Trajtenberg F, Mansilla MC, Haouz A, Alzari PM, de Mendoza D, Buschiazzo A: Structural plasticity and catalysis regulation of a thermosensor histidine kinase. Proceedings of the National Academy of Sciences of the United States of America. 2009, 106 (38): 16185-16190. 10.1073/pnas.0906699106.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. Carmany DO, McCleary WR: Genetic and biochemical studies of phosphate activity of PhoR. J Bacteriol. 2003, 185: 1112-1115. 10.1128/JB.185.3.1112-1115.2003.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references


We thank Kathryn Hanks, Alex Cummock, Evan Christensen, Bethany Evans, Gregory Bowden, and Michael Barrus for help with sequence collection and mutant characterization. This work was supported by Public Health Service grant R15GM96222 from the National Institute of General Medical Sciences.


The publication costs for this article were funded by the Department of Microbiology and Molecular Biology at Brigham Young University, the College of Life Sciences at Brigham Young University, and the Public Health Service grant R15GM96222 from the National Institute of General Medical Sciences.

This article has been published as part of BMC Genetics Volume 16 Supplement 2, 2015: Selected articles from The 11th Annual Biotechnology and Bioinformatics Symposium (BIOT-2014): Genetics. The full contents of the supplement are available online at

Author information

Authors and Affiliations


Corresponding author

Correspondence to William R McCleary.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

SGG planned, performed experiments, prepared, and edited the manuscript. JM performed DCA, helped with analysis of data, and aided in preparing the manuscript. TD, TR, and ME helped in the scanning mutagenesis. PR helped with bioinformatic experimental planning and analysis. WRM designed experiments, conducted experiments, aided in data analysis, and aided in preparing and editing the manuscript.

Electronic supplementary material

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit

The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Gardner, S.G., Miller, J.B., Dean, T. et al. Genetic analysis, structural modeling, and direct coupling analysis suggest a mechanism for phosphate signaling in Escherichia coli. BMC Genet 16 (Suppl 2), S2 (2015).

Download citation

  • Published:

  • DOI: