Structural analyses unravel the molecular mechanism of cyclic di-GMP regulation of bacterial chemotaxis via a PilZ adaptor protein

The bacterial second messenger cyclic di-GMP (c-di-GMP) has emerged as a prominent mediator of bacterial physiology, motility, and pathogenicity. c-di-GMP often regulates the function of its protein targets through a unique mechanism that involves a discrete PilZ adaptor protein. However, the molecular mechanism for PilZ protein–mediated protein regulation is unclear. Here, we present the structure of the PilZ adaptor protein MapZ cocrystallized in complex with c-di-GMP and its protein target CheR1, a chemotaxis-regulating methyltransferase in Pseudomonas aeruginosa. This cocrystal structure, together with the structure of free CheR1, revealed that the binding of c-di-GMP induces dramatic structural changes in MapZ that are crucial for CheR1 binding. Importantly, we found that restructuring and repositioning of two C-terminal helices enable MapZ to disrupt the CheR1 active site by dislodging a structural domain. The crystallographic observations are reinforced by protein–protein binding and single cell–based flagellar motor switching analyses. Our studies further suggest that the regulation of chemotaxis by c-di-GMP through MapZ orthologs/homologs is widespread in proteobacteria and that the use of allosterically regulated C-terminal motifs could be a common mechanism for PilZ adaptor proteins. Together, the findings provide detailed structural insights into how c-di-GMP controls the activity of an enzyme target indirectly through a PilZ adaptor protein.

c-di-GMP 4 is a widespread second messenger that controls a multitude of cellular functions and bacterial behaviors (1)(2)(3). In human pathogenic bacteria such as Pseudomonas aeruginosa, Salmonella enterica, and Vibrio cholerae, c-di-GMP plays a prominent role regulating the motile-sessile transition and formation of antibiotic-resistant biofilm communities. c-di-GMP exerts its control on cellular functions by binding to a variety of protein and riboswitch targets (4,5). Among the known c-di-GMP-binding proteins, discrete single-domain PilZ proteins are widespread and by far the most prevalent (6 -12). As the single-domain PilZ proteins only contain 80 -130 amino acids and do not exhibit enzymatic or other activities, their physiological roles often remain elusive. Emerging evidence suggests that the single-domain PilZ proteins function as adaptor proteins to enable c-di-GMP to control its protein targets indirectly (8,9,13,14). However, the detailed molecular mechanism of how c-di-GMP regulates the activity of its protein targets through the PilZ adaptor proteins remains unknown. How specificity in c-di-GMP signaling, which is crucial for preventing cross-talk among different pathways, is achieved through the PilZ adaptor proteins also remains a key question.
We recently discovered that c-di-GMP controls the enzymatic activity of the methyltransferase CheR1 in P. aeruginosa through the PilZ adaptor protein MapZ (14). High cellular c-di-GMP concentrations lead to the formation of the MapZ-c-di-GMP-CheR1 ternary complex and subsequently inhibit the methyltransferase activity of CheR1. Inhibition of CheR1 results in a decrease in the methylation level of methyl-accepting chemoreceptor proteins (MCPs), which ultimately leads to the suppression of flagellar motor switching and a decrease of swimming reversal frequency. Because cellular c-di-GMP concentration fluctuates with metabolic and extracellular stimuli, P. aeruginosa can fine-tune its chemotactic response and motility in different environments using the MapZ-dependent regulatory mechanism. Although it is known that the protein-protein interaction between CheR1 and MapZ is highly specific and dependent on c-di-GMP, the structural basis for the c-di-GMP-mediated specific protein-protein interaction is unknown. Here, we report the crystal structures of the free CheR1 and the MapZ-c-di-GMP-CheR1 ternary complex along with other experimental data to support the physiological relevance of the structural observations. The results reveal an . The authors declare that they have no conflicts of interest with the contents of this article. This article contains Figs. S1-S9, Tables S1 and S2, and Movies S1-S9. The atomic coordinates and structure factors (codes 5Y4R and 5Y4S) have been deposited in the Protein Data Bank (http://wwpdb.org/). 1 Both authors contributed equally to this work. 2 To whom correspondence may be addressed. E-mail: zxliang@ntu.edu.sg. 3 To whom correspondence may be addressed. E-mail: ygao@ntu.edu.sg. 4 The abbreviations used are: c-di-GMP, cyclic di-GMP; MCP, methyl-accepting chemoreceptor protein; aa, amino acids; Bs, Bacillus subtilis; AdoMet, S-adenosylmethionine; SAH, S-adenosyl-L-homocysteine; Bis-Tris, 2-[bis(2hydroxyethyl)amino]-2-(hydroxymethyl)propane-1,3-diol; RNAP-␤, RNA polymerase ␤-subunit; ITC, isothermal titration calorimetry.
cro ARTICLE exquisite allosteric mechanism that is potentially conserved in many PilZ adaptor proteins.

Overall structure of the MapZ-c-di-GMP-CheR1 ternary complex
We determined the crystal structure of CheR1 in complex with MapZ and c-di-GMP at 2.3-Å resolution ( Fig. 1A and Fig.  S1) with the crystallographic statistics listed in Table 1. The difference Fourier map (F o Ϫ F c ) clearly shows the ligand c-di-GMP (Fig. S2). The c-di-GMP-bound MapZ is composed of a ␤-barrel core consisting of eight antiparallel ␤-strands, an N-terminal loop, and two C-terminal ␣-helices as seen in other PilZ proteins (15)(16)(17)(18)(19) (Fig. 1B). Observed as an intercalated dimer in the electron density map, c-di-GMP is held in position by residues stemming from the N-terminal loop and strands ␤4, ␤5, ␤7, and ␤8. The methyltransferase CheR1 is composed of two domains (Fig. 1C). The N-terminal domain of CheR1 is a relatively small domain consisting of four perpendicularly packed ␣-helices, whereas the C-terminal domain contains a large ␣/␤ fold core and a small ␤-subdomain (Fig. 1C). The core of the C-terminal domain adopts a Rossmann fold consisting of a mixed seven-stranded ␤-sheet flanked by ␣-helices. The architecture of the C-terminal domain core is similar to those of other class I methyltransferases (20). In addition to the ␣/␤ fold core, the C-terminal domain also contains a ␤-subdomain consisting of two short helices and three antiparallel strands (Fig.  1C). The ␤-subdomain, which is inserted between ␣7 and ␣10, is unique to chemotaxis methyltransferases (21-23) (Fig. S3). The N-and C-terminal domains are connected by a flexible domain linker (aa 70 -80) with a cleft formed between the two domains. In the ternary complex structure, the ␤-barrel core of MapZ sits on the top of CheR1 with the two C-terminal helices ␣1 and ␣2 wedged into the central cleft of CheR1, resulting in extensive contacts between the two proteins that are likely to account for the formation of a stable MapZ-c-di-GMP-CheR1 ternary complex.

Specific interactions between MapZ and CheR1
In the ternary complex, MapZ makes direct contacts with CheR1 through its C-terminal motifs, the ␤6-␤7 loop, and a small segment of the N-terminal loop. The C-terminal motifs of MapZ encompass the two helices ␣1 (aa 92-106) and ␣2 (aa 109 -121) and the adjacent loops. Helix ␣1 seems to play a cen-

PilZ adaptor-mediated c-di-GMP signaling
tral role in mediating the MapZ-CheR1 interaction by situating in the central cleft of CheR1 with several highly specific polar interactions formed between ␣1 and CheR1. In particular, Arg 99 interacts with Asp 67 and Thr 70 of CheR1 through salt bridges and hydrogen bonds, and Arg 100 forms a salt bridge with Asp 144 (Fig. 2A). The side chains of Asp 93 , Ser 96 , His 97 , and Asn 105 and the main chain oxygen of Leu 104 also interact with various CheR1 residues through hydrogen bonds. The side chains of several residues from ␣1, including Leu 92 , Ile 95 , Leu 101 , and Leu 104 , as well as Arg 99 , Arg 100 , and Glu 103 interact with various CheR1 residues through hydrophobic forces. In addtition, ␣2 seems to interact with CheR1 mainly through hydrophobic interaction via residues Leu 112 , Val 116 , Ala 117 , Leu 119 , and Val 120 (Fig. 2B). The only noticeable direct polar interaction that involves ␣2 seems to be the hydrogen bond formed between Glu 109 (␣2) and Ser 146 (CheR1). Apart from the residues in ␣1 and ␣2, the main chain groups of Leu 92 and Ile 90 from the loop prior to ␣1 form three hydrogen bonds with Val 22 and Gly 20 of CheR1 (Fig. 2C), whereas His 123 and Asp 124 from the loop downstream of ␣2 interact with Thr 18 and Arg 62 of CheR1 through a hydrogen bond and salt bridge, respectively.
In addition to the residues from the C-terminal motifs, several residues from other non-C-terminal regions of MapZ are also involved in direct binding of CheR1. Asp 66 and Val 67 from the ␤6-␤7 loop interact with Arg 78 , Ser 250 , and Ile 221 of CheR1 through hydrogen bond and hydrophobic interactions, whereas Phe 16 and Asp 17 from the N-terminal loop interact with CheR1 residues through hydrophobic and hydrogen bond interactions, respectively (Fig. 2C).

Binding affinity measurements reinforce the key role of the MapZ C-terminal residues in binding CheR1
The extensive contacts between the two proteins are consistent with the high thermodynamic binding affinity (K d ϭ 28.6 Ϯ 17.7 nM) between the two proteins. c-di-GMP is essential for MapZ-CheR1 interaction because no binding could be detected in the absence of c-di-GMP. To evaluate the contribution of the aforementioned residues of ␣1 and ␣2 to CheR1 binding, we measured the binding affinity of several MapZ mutants for CheR1. The results (Fig. 3, A and B) show that Arg 99 is the most crucial residue as the substitution of Arg 99 by Ala reduced the binding affinity to below the detection limit, consistent with the structural observation that the residue Arg 99 engages in multiple interactions with CheR1 residues (Fig. 2A). The substitution of other ␣1 residues, His 97 , Arg 100 , and Leu 104 , by Ala reduced the K d by 7.1-, 6-, and 2.4-fold, respectively. Replacement of Asp 93 , which positions His 97 for interacting with CheR1 residues, caused a 10.4-fold decrease in K d . In comparison, removal of the entire ␣2 helix (MapZ(⌬␣2)) decreased the K d by 14.6-fold. These observations support the relevance of the crystal structure and reinforce the central role of ␣1 in mediating the MapZ-CheR1 interaction.

The C-terminal binding residues are important for the function of MapZ in suppressing flagellar motor switching
To corroborate the structural and binding data and further validate the importance of the two C-terminal helices for the in vivo function of MapZ, we investigated how mutations in MapZ impact flagellar motor switching in P. aeruginosa cells by cell tethering assay (14,24). We previously observed that the overexpression of MapZ in P. aeruginosa cells dramatically suppressed flagellar motor switching, which is similar to the effect caused by the deletion of cheR1 (14). Most of the MapZ-overexpressing cells cannot switch motor direction and exhibit unidirectional flagellar rotation (14). The suppression is attributed to the ability of MapZ to inhibit CheR1 and decrease MCP methylation at high c-di-GMP levels. Hence, it was expected that overexpression of defective MapZ mutants that lack key residues for c-di-GMP or CheR1 binding would exhibit reduced suppression on flagellar motor switching. Indeed, the cell tethering assays showed that overexpression of MapZ mutants that lack binding residues for CheR1 suppressed motor switching to a lesser extent than MapZ. MapZ(R99A) did not seem to suppress motor switching at all, whereas MapZ(D93A), MapZ(H97A), MapZ(R110A), MapZ(L104A), and MapZ(⌬␣2) exhibited less suppression than MapZ to various extents ( Fig.  4A and Movies S1-S8). Overexpression of MapZ(R13A), which is defective in c-di-GMP binding (14), also exhibited a much weaker suppression effect. As supported by Western blotting assays, the different suppression effects observed for the mutant proteins are not caused by a reduction in protein expression level ( Fig. 4B and Fig. S4). Taken together, the results provide strong support for the physiological relevance of the MapZ-c-di-GMP-CheR1 ternary structure and establish the pivotal role of the C-terminal motifs in mediating MapZ-CheR1 interaction and suppressing flagellar motor switching.

PilZ adaptor-mediated c-di-GMP signaling Structural rearrangement in MapZ induced by the binding of c-di-GMP and CheR1
The structures of MapZ with and without c-di-GMP bound were determined previously by NMR spectroscopy (25,26). It was noted that MapZ undergoes large structural rearrangements upon the binding of c-di-GMP. Comparing the ternary complex structure with the NMR structures (Protein Data Bank code 1YWU (26)) reveals how the structural rearrangement induced by c-di-GMP binding prepares MapZ for CheR1 binding. Upon the binding of c-di-GMP, the N-terminal loop becomes organized by wrapping itself around the intercalated c-di-GMP dimer (Fig. 5A). Without the c-di-GMP-induced reorganization, the N-terminal loop would hinder the binding of CheR1 sterically (Fig. 5A). In the C-terminal region, one of the most notable changes is that the binding of c-di-GMP results in the straightening of the initially bent ␣2 and concomitant dissociation of ␣2 from the ␤-barrel core (Fig. 5A). The restructuring and displacement of ␣2 is crucial for setting up MapZ for the interaction with CheR1.
Unexpectedly, the c-di-GMP binding mode in the ternary complex differs to some extent from that in the c-di-GMP-MapZ binary complex (Protein Data Bank code 2L74) (25). The N-terminal loop adopts a rather dissimilar conformation as characterized by the different positions and orientations of several loop residues (e.g. Glu 7 , Phe 11 , and His 12 ) in the ternary complex (Fig. 5B). Also, noticeably, the helix ␣1 from the C terminus is tilted toward the ␤-barrel core by ϳ2 Å, and ␣2 is pushed toward ␣1 by as much as 12 Å in the ternary complex. Although residues Arg 8 , Arg 9 , Arg 10 , Arg 13 , Asp 35 , and Trp 77 are seen to interact with c-di-GMP directly in both the binary and ternary complexes, three additional residues (Glu 7 , Phe 11 , and Arg 79 ) are directly involved in the binding of c-di-GMP only in the ternary complex (Fig. 5, C and D). It should be noted that although Glu 7 and Phe 11 are conserved in most single-domain PilZ proteins, their functions remain unknown. The above observations suggest that the two residues are conserved because they play a direct role in c-di-GMP binding. On the contrary, His 12 , which is not a conserved residue and seen to form a hydrogen bond with c-di-GMP in the binary complex, does not participate in c-di-GMP binding in the ternary complex. The discrepancies in c-di-GMP binding between the binary and ternary complexes highlight the malleability of the c-di-GMP-binding pocket and call for caution when interpreting the structures of PilZ proteins determined in the absence of their protein partners.

Binding of MapZ to CheR1 resulting in disruption of the S-adenosylmethionine (AdoMet)-binding pocket and inhibition of methyltransferase activity
We also determined the crystal structure of free CheR1 protein to determine how the binding of MapZ impacts the structure of CheR1. Crystals of free CheR1 belong to the tetragonal space group I4 1 with 10 copies of the molecule (chains A-J) in an asymmetric unit. Superposition of crystallographic indepen- PilZ adaptor-mediated c-di-GMP signaling dent chains B-J onto chain A gives an average root mean square deviation of 0.76 Å over main-chain atoms, indicating almost identical conformations and a limited influence of crystal pack-ing on the individual chain. The structure was refined to 3.4-Å resolution, and the simulated annealing 2F o Ϫ F c composite omit map clearly demonstrates that the two domains and the domain linker are nicely modeled (Fig. S5). Comparison of the structures of the free CheR1 and the ternary complex shows that although the C-terminal domains are almost superimposable the N-terminal domains adopt very different conformations in the two structures (Fig. 5E). The spatial positioning of the N-terminal domain relative to the C-terminal domain is also different from that of BsCheR (29) (Fig. S6) despite the fact that the structures of the N-terminal domains are highly similar (Fig. S7). The N-terminal domain and the domain linker of CheR1 have relatively higher B-factors (Fig. S8), suggesting that they are intrinsically flexible. The most important implication from the comparison of the structures of free CheR1 and the ternary complex is that the binding of MapZ induces a relocation of the N-terminal domain by a translocation of more than 10 Å and a rotation of ϳ105°as well as a noticeable conformational change in the domain linker ( Fig.  5E and see Movie S9).
Based on the studies on other chemotaxis methyltransferases (22,23), the C-terminal domain and the domain linker are likely involved in the binding of the cofactor AdoMet. In the free CheR1 structure, an AdoMet-binding pocket can be readily identified by overlaying the structure with that of the S-adenosyl-L-homocysteine (SAH)-bound BsCheR (Protein Data Bank code 5FTW; Fig. 5F). Two highly conserved catalytic residues (Asp 144 and Tyr 222 ) that are essential for AdoMet binding can also be identified. In the ternary complex, the AdoMet-binding pocket is completely occupied by MapZ upon dislodging of the N-terminal domain and conformational changes in the domain linker. The AdoMet-binding pocket is filled by residues from the ␣1 of MapZ with Asp 144 and Tyr 222 engaging in multiple interactions with residues from ␣1 ( Figs. 2A and 5F). These observations strongly suggest that the C-terminal motifs of MapZ, particularly ␣1, have evolved to interact specifically with the conserved residues in the CheR1 active site to inhibit its enzymatic activity. MapZ only interacts with CheR1 but not the other three chemotaxis methyltransferases of P. aeruginosa (27,28). Sequence alignment of the four methyltransferases suggests that this is likely because some of the MapZ-binding residues in CheR1 are absent in the other three CheR proteins (Fig. S9).

MapZ-mediated regulation of chemotaxis by c-di-GMP is conserved in many proteobacteria
Together with our previous study (14), the current study firmly establishes the regulation of chemotaxis by c-di-GMP through a MapZ-mediated mechanism in P. aeruginosa. We further determined whether this regulatory mechanism is conserved in the genus Pseudomonas and other non-Pseudomonas species by searching MapZ orthologs/homologs encoded by Figure 3. Assessment of the contribution of the residues from helices ␣1 and ␣2 to protein-protein interaction by binding affinity measurement. A, binding isotherm (top panels) and data-fitting curves (bottom panels) obtained from ITC assay. Solution containing CheR1 was titrated into solution containing MapZ (or a mutant) and saturating c-di-GMP (250 M). The experimental conditions are described under "Experimental procedures." The binding measurements were performed in two independent experiments, and the representative data set is presented here. DP, differential power. B, apparent dissociation constant (K d ) values obtained by fitting the isotherm data to a 1:1 binding model.

PilZ adaptor-mediated c-di-GMP signaling
bacterial genomes. We first identified and collected the sequences of a total of 11,781 single-domain PilZ proteins from the InterPro protein database with most sequences belonging to the IPR009875 family (29). From this collection, we identified 1,063 potential functional orthologs/homologs of MapZ according to the following two criteria. First, the proteins The two methyltransferase (BsCheR in gray) structures are aligned according to their C-terminal domains. SAH is shown as a stick with carbon, oxygen, nitrogen, and sulfur atoms colored yellow, red, blue, and brown, respectively. In the enlarged view of the cofactor-binding site, it can be seen that Arg 100 and Leu 104 from ␣1 of MapZ occupy the SAH-binding site to prevent the binding of the cofactor AdoMet. The highly conserved catalytic residue Asp 144 /Asp 130 (CheR1/BsCheR) that interacts with Arg 100 (MapZ) is also highlighted as a stick.
PilZ adaptor-mediated c-di-GMP signaling must share high sequence similarity (Ͼ60% similarity) with MapZ. Second, the proteins must contain two highly similar C-terminal helices with the CheR1-binding residues conserved. As shown by the sequence logo (Fig. 6A), the critical CheR1-interacting residues are highly conserved in the putative orthologs/homologs.
A taxonomic and phylogenetic relationship analysis for MapZ and the 1,063 proteins suggests that most MapZ orthologs/homologs come from Pseudomonas, Vibrio, and other proteobacteria, in particular ␥-proteobacteria (Fig. 6B). A small percentage of the orthologs/homologs were encoded by genes from unclassified bacterial species uncovered by metagenomics projects. PlzA, one of the PilZ proteins from the pathogenic bacterium V. cholerae (9) with unknown cellular function was also identified as an ortholog/homolog of MapZ. This hints that PlzA is likely to play a similar role in V. cholerae by binding to a chemotactic methyltransferase to regulate chemotaxis. Furthermore, for the MapZ homolog/ortholog-containing bacterial species whose genomes have been sequenced, we could identify genes encoding chemotaxis methyltransferases that are highly similar to CheR1. Altogether, the bioinformatics analysis suggests that the MapZ-mediated mechanism for fine-tuning chemotactic response and bacterial motility is not a unique mechanism restricted to P. aeruginosa but a widely distributed mechanism adapted by many proteobacteria.
In summary, the studies described here establish the structural basis for the regulation of the methyltransferase CheR1 by c-di-GMP through the PilZ protein MapZ. The data support a model whereby MapZ, at high cellular c-di-GMP concentrations, undergoes dramatic conformational changes and occupies the CheR1 active site to inhibit the enzyme activity (Fig. 7). Inhibition of CheR1 decreases the methylation level of MCPs and ultimately results in the suppression of flagellar motor switching. Further analysis of the sequences of about 12,000 single-domain PilZ proteins encoded by sequenced bacterial genomes revealed that most of them contain C-terminal motifs that vary in sequence length and secondary structure. Although the variation in sequence and structure may reflect the diverse protein partners with which the PilZ proteins interact, the mechanism of allosteric regulation of C-terminal motifs for controlling protein-protein interaction observed in MapZ could be a common theme for the PilZ adaptor proteins.

Materials
All bacterial strains, plasmids, and primers used in this study are listed in Tables S1 and S2.

Protein expression and purification
For the recombinant proteins used for crystallization, the CheR1-and MapZ-encoding genes PA3348 and PA4608 were amplified by polymerase chain reaction (PCR) using the genomic DNA of P. aeruginosa PAO1 as template. The PCR fragments were cloned into pACYCDuet-1 expression vector (Novagen), which provides a hexahistidine tag at the N terminus of recombinant CheR1 and MapZ, respectively. In brief, the expression plasmid carrying cheR1 or mapZ gene was transformed into Escherichia coli BL21/DE3 strain using a heat-shock method. Protein expression was induced by 0.2 mM isopropyl ␤-D-thiogalactoside at 16°C overnight. Cells were harvested by centrifugation, and the pellet was resuspended and homogenized five times using an EmulsiFlex-C3 (Avestin). The lysate was centrifuged at 40,000 rpm for 1 h to remove insoluble cell debris. The supernatant was applied to a nickelnitrilotriacetic acid affinity column, sequentially washed, and eluted by a gradient of increasing imidazole concentration in buffer. The pooled fractions were treated with tobacco etch virus protease and dialyzed against a buffer containing 20 mM Tris-HCl, pH 7.5, 100 mM NaCl. After removal of tobacco etch virus protease and uncleaved protein by the second affinity column, the protein samples were further purified by size-exclusion chromatography. The target protein was concentrated to ϳ15 mg/ml for subsequent crystallization and biochemical analysis. The mutant MapZ proteins were generated using a QuikChange kit (Stratagene) and prepared following the aforementioned procedure except the final buffer contained20 mM Tris-HCl, pH 8.0, 100 mM NaCl, 1 mM tris(2-carboxyethyl) phosphine hydrochloride (Sigma).

Crystallization and data collection
The initial crystal of the isolated CheR1 was achieved at 20°C by mixing 1.0 l of protein with 1.0 l of reservoir solution containing 1.5 M sodium phosphate monobasic/potassium phosphate dibasic, pH 5.0. For further optimization, the crystallization condition was supplemented with 1% (v/v) glycerol and 3% (v/v) dextran sulfate sodium salt (M r 5,000). The crystals were flash frozen (100 K) in the reservoir solution supplemented with 30% (v/v) glycerol. The native data set (total of 3,600 images with 0.2 oscillations) was collected at PXIII, Swiss Light Source.
For the MapZ-c-di-GMP-CheR1 complex, the molar ratio between the two proteins was determined to be 1:1 by analytical gel filtration. The MapZ-c-di-GMP-CheR1 complex was prepared by incubating 0.3 mM CheR1 with 0.3 mM MapZ and 1 mM c-di-GMP (Sigma-Aldrich). Crystals of this ternary complex were grown at 20°C by mixing 2.0 l of protein with 1.0 l of reservoir solution containing 0.1 mM Bis-Tris, pH 5.5, 19% PEG 3350, and 0.2 M lithium sulfate. The crystals were flash frozen (100 K) in the reservoir solution with a final concentration of 40% (v/v) PEG 3350. The data set (total of 2,000 images with 0.1 oscillations) for MapZ-c-di-GMP-CheR1 complex was collected at the Swiss Light Source.

Data processing, structure determination, and refinement
All data were processed using the XDS package (30). The structure of MapZ-c-di-GMP-CheR1 ternary complex was determined by molecular replacement with Phenix combined with Phaser (31) using BsCheR (Protein Data Bank code 5FTW) and MapZ (Protein Data Bank code 2L74) as search models. The structure of apoCheR1 was determined by molecular replacement (Phaser/Phenix) (31) using the CheR1 model from our ternary complex as a search model. The models were improved by manual model building with Coot (32) and refined using Phenix (31). Crystallographic statistics of the structures are listed in Table 1, and all structural figures were generated with PyMOL (Schrödinger, LLC). Figure 6. MapZ homologs are widespread in proteobacteria. A, sequence log for MapZ and its 1,063 functional homologs. The two essential c-di-GMP-binding motifs (RXXXR and DXSXXG) and the C-terminal ␣1 and ␣2 helices are highlighted by bars. The residues from ␣1 and ␣2 that are involved in direct binding of CheR1 are highlighted by red residue numbers. B, phylogenetic tree analysis of MapZ homologs. The clades for the two major Pseudomonadales and Vibrionales orders are colored. Most of the remaining homologs belong to other orders from the Proteobacteria phylum, but a small number of homologs belong to unclassified species whose sequences were obtained by metagenomics projects. A FASTA file that contains the 1,063 sequences can be found in the supporting information.

Preparation of MapZ overexpression mutant strains and Western blotting assay
The DNA fragments of MapZ mutants, including MapZ (D93A), MapZ(H97A), MapZ(R99A), MapZ(R100A), and MapZ(L104A) were synthesized by Integrated DNA Technologies, Inc. Synthesized DNA fragments were digested and ligated into the pUCP18 overexpression vector (Novagen) (33). The constructs were verified by sequencing and transformed into P. aeruginosa mPAO1 strain following the established procedure (34). Single colonies of overexpression strains were used to inoculate 10 ml of LB medium that was kept at 37°C with shaking overnight. One milliliter of culture was centrifuged at 11,325 ϫ g for 10 min at 4°C. Cells were resuspended with 100 l of 2ϫ loading dye and heated at 95°C for 10 min. Samples were analyzed by SDS-PAGE and Western blotting using THE TM HA tag antibody (HRP) (GenScript) and RNA polymerase ␤-subunit (RNAP-␤) (loading control) antibody (a kind gift from Dr. Wang).

Binding affinity measurement by isothermal titration calorimetry (ITC)
ITC measurements were performed at 25°C using a Micro-Cal ITC200 (MicroCal Inc.). Protein samples were dialyzed into a buffer containing 20 mM Tris-HCl, pH 8.0, 100 mM NaCl, 1 mM tris(2-carboxyethyl)phosphine hydrochloride. A sample syringe with stirring speed of 800 rpm was used to titrate the CheR1 (250 -320 M) into a cell containing 25-30 M MapZ or MapZ mutant in the presence of c-di-GMP (250 M), which was prepared enzymatically as described previously (35,36). The titration comprised 19 injections of 2 l each separated by 150-s equilibration time. The isothermal data were analyzed using the Origin 7.0 program, fitting to a single-site binding model.

Cell tethering assay for the analysis of the flagellar rotation
Single colonies of P. aeruginosa PAO1 MapZ and mutants overexpression strains were used to inoculate 10 ml of LB broth (BD Biosciences) that was kept at 37°C with shaking at 250 rpm overnight. Cultures were diluted to an A 600 of 0.2 using M9 medium and grown at 37°C with shaking at 250 rpm until the cells reached the late-exponential growth phase at A 600 of 0.8. Cell cultures were then diluted 10-fold before the microfluidic experiments. Cell tethering assays were performed by using a microfluidic stagnation flow device precoated with flagellar antibodies (37). Flagella were sheared off by passing the bacterial cells through a 34-gauge blunt-end needle four times. Cells were loaded into the cell chamber of the flow device, and nontethered cells were washed away using M9 medium. Cells were PilZ adaptor-mediated c-di-GMP signaling visualized using an inverted microscope (Nikon TE2000-U) under a 40ϫ objective. Videos of tethered bacteria were taken at 40 frames/s for 1-5 min using a complementary metal-oxide semiconductor camera (Thorlabs DCC1645). The rotation of the cells was monitored and recorded using a procedure reported previously (38). The time spent by each cell in clockwise, counterclockwise, or pause phases was tallied in 20-s intervals. Image processing of the tethered cells was carried out as described previously (39) with videos of the tethered cells converted to grayscale and the contrast of the cells adjusted (saturated pixels, 0.4; histogram stack) using NIH ImageJ. Images of the cells were binarized to isolate them from the background for every frame in the video. The center of mass coordinates of the rotating cell bodies of each cell were measured from the obtained binarized image stack and imported into our in-house analysis program BTAP (bacterial tethering analysis program) in MATLAB. Student's two-tailed t test was used to test the statistical significance of the data. The normality assumption was considered to be sound because, by the central limit theorem, sample means of moderately large samples are well-approximated by a normal distribution.

Genome mining and phylogenetic relationship analysis
Single-domain PilZ proteins were identified by mining PilZ protein sequences deposited in the InterPro database (http:// www.ebi.ac.uk/interpro/) 5 (29). A total of 11,781 sequences of single-domain PilZ proteins were then submitted to the EFI/EST sequence similarity network server (http://efi.igb. illinois.edu/efi-est/tutorial_startscreen.php) 5 (40) for analysis and generation of a sequence similarity network by using an alignment score of 16. The clusters formed in the network were analyzed and visualized using Cytoscape v3.5. MapZ was found to be clustered with a group of 1,023 proteins in the network. The sequences of the 1,023 proteins were curated manually to verify whether they contain the crucial C-terminal helices and key residues for CheR1 binding.
For the generation of the sequence logo, the 1,024 protein sequences contained in a FASTA file were first aligned using ClustalW. The FASTA file containing the aligned sequences is included in the supporting information and can be visualized using the free software UGENE and MEGA. ClustalW was also used to generate the phylogenetic tree file (.ph file) to allow visualization of the tree using the web-based tool interactive Tree Of Life (iTOL) (41).