# Coiled-Coil Antagonism Regulates Activity of Venus Flytrap-Domain-Containing Sensor Kinases of the BvgS Family

Elodie Lesne, Elian Dupré, Marc Lensink, Camille Locht, Rudy Antoine,
Françoise Jacob-Dubuisson

► To cite this version:

Elodie Lesne, Elian Dupré, Marc Lensink, Camille Locht, Rudy Antoine, et al. Coiled-Coil Antagonism Regulates Activity of Venus Flytrap-Domain-Containing Sensor Kinases of the BvgS Family. mBio, 2018, 9 (1), pp.e02052-17. 10.1128/mBio.02052-17. hal-02121569

HAL Id: hal-02121569

https://hal.science/hal-02121569

Submitted on 6 May 2019

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Coiled-Coil Antagonism Regulates Activity of Venus Flytrap-Domain-Containing Sensor Kinases of the BvgS Family

Elodie Lesne,ª,ᵇ,ᶜ,ᵈ,ᵉ Elian Dupré,ª,ᵇ,ᶜ,>d,>e* Marc F. Lensink,f Camille Locht,ª,ᵇ,ᶜ,>d,>e Rudy Antoine,ª,ᵇ,ᶜ,>d,>e

fFrançoise Jacob-Dubuissonª,ᵇ,ᶜ,>d,>e

**ABSTRACT** *Bordetella pertussis* controls the expression of its virulence regulon through the two-component system BvgAS. BvgS is a prototype for a family of multidomain sensor kinases. In BvgS, helical linkers connect periplasmic Venus flytrap (VFT) perception domains to a cytoplasmic Per-Ant-Sim (PAS) domain and the PAS domain to the dimerization/histidine phosphotransfer (DHp) domain of the kinase. The two linkers can adopt coiled-coil structures but cannot do so simultaneously. The first linker forms a coiled coil in the kinase mode and the second in the phosphatase mode, with the other linker in both cases showing an increase in dynamic behavior. The intervening PAS domain changes its quaternary structure between the two modes. In BvgS homologues without a PAS domain, a helical "X" linker directly connects the VFT and DHp domains. Here, we used BvgS as a platform to characterize regulation in members of the PAS-less subfamily. BvgS chimeras of homologues with natural X linkers displayed various regulation phenotypes. We identified two distinct coiled-coil registers in the N- and C-terminal portions of the X linkers. Stable coil formation in the C-terminal moiety determines the phosphatase mode, similarly to BvgS; in contrast, coil formation in the N-terminal moiety along the other register leads to the kinase mode. Thus, antagonism between two registers in the VFT-DHp linker forms the basis for activity regulation in the absence of the PAS domain. The N and C moieties of the X linker play roles similar to those played by the two independent linkers in sensor kinases with a PAS domain, providing a unified mechanism of regulation for the entire family.

**IMPORTANCE** The whooping cough agent *Bordetella pertussis* uses the BvgAS sensory transduction two-component system to regulate production of its virulence factors. BvgS serves as a model for a large family of multidomain bacterial sensor kinases. *B. pertussis* is virulent when BvgS functions as a kinase and avirulent when it switches to phosphatase activity in response to modulating signals. Understanding the molecular regulation of those proteins might lead to new antibacterial strategies. Here, we show that the linker regions between the perception and the enzymatic domains shift between distinct states of conformation in an alternating manner in response to signals and that their antagonistic changes control sensor kinase activity. These linker regions and mechanistic principles appear to be conserved among BvgS homologues, irrespective of the presence or absence of an intervening domain between the perception and the enzymatic domains. This work has thus uncovered general molecular mechanisms that regulate activity of sensor kinases in the BvgS family.

Received 1 December 2017
Accepted 24 January 2018
Published 27 February 2018

**Citation** Lesne E, Dupré E, Lensink MF, Locht C, Antoine R, Jacob-Dubuisson F. 2018. Coiled-coil antagonism regulates activity of Venus flytrap-domain-containing sensor kinases of the BvgS family. *mbio* 9:e02052-17. https://doi.org/10.1128/mbio.02052-17.

**Editor** Scott J. Hultgren, Washington University School of Medicine

**Copyright** © 2018 Lesne et al. This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

Address correspondence to Rudy Antoine, rudy.antoine@pasteur-lille.fr, or Françoise Jacob-Dubuisson, francoise.jacob@ibl.cnrs.fr.

* Present address: Elian Dupré, University of Lille, CNRS, UMR 8576, UGSF, Villeneuve d'Ascq, France.

<footnote>ªUniversity of Lille, Lille, France</footnote>

<footnote>ᵇCNRS UMR 8204, Lille, France</footnote>

<footnote>ᶜInserm U1019, Lille, France</footnote>

<footnote>ᵈCHU Lille, Lille, France</footnote>

<footnote>ᵉCentre d'Infection & d'Immunité de Lille, Institut Pasteur de Lille, Lille, France</footnote>

<footnote>fUniversity of Lille, CNRS, UMR 8576, UGSF—Unité de Glycobiologie Structurale & Fonctionnelle, Villeneuve d'Ascq, France</footnote>
**KEYWORDS** *Bordetella pertussis*, BvgS family, coiled coil, sensor kinases, two-component regulatory systems

Bordetella pertussis is the agent of whooping cough, an acute respiratory disease. Colonization of the human respiratory tract requires the production of virulence factors whose expression is coordinately regulated by a sensory transduction two-component system (TCS) called BvgAS. TCSs are widespread in bacteria and involved in regulating various cellular processes, including metabolism, mobility, bacterial development, and virulence. Perception of a specific signal by the sensor kinase leads to autophosphorylation of the kinase domain and transfer of the phosphoryl group to the response regulator (1-3).

BvgA is a classical response regulator and in its phosphorylated form activates the transcription of the virulence genes, thus setting the bacterium in the virulent, Bvg-positive (Bvg⁺) phase. Under standard culture conditions, i.e., at 37°C and without ligand, BvgS is in a kinase mode. In the laboratory, BvgS shifts to a phosphatase mode upon perception of modulators such as MgSO₄ and nicotinate or derivatives or under conditions of nutrient restriction or low temperature. This sets the bacteria in the avirulent, Bvg⁻ phase (4-7). An intermediate virulent phase, called Bvg¹ and characterized by the expression of only a subset of virulence factors, is encountered at intermediate concentrations of modulators (8).

BvgS is a dimeric sensor kinase that serves as a model for a large family of sensor kinases harboring periplasmic Venus flytrap (VFT) domains for signal perception (9, 10). Each monomer of BvgS is composed of two periplasmic VFT domains (Pfam SBP_bac_3) in tandem; a transmembrane (TM) segment; a cytoplasmic Per-Ant-Sim (PAS) domain; and a histidine kinase (HK), comprising a dimerization/histidine phosphotransfer (DHp) domain and a catalytic ATP-binding (CA) domain (10). They are followed by receiver and histidine phosphotransfer (HPt) domains, which form a phosphorelay (11).

Two linkers, called linker 1 and linker 2, connect the VFT domains to the PAS domain and the PAS domain to the DHp domain, respectively (Fig. 1A). Found in many sensory transduction proteins, PAS domains are involved in signal perception and/or in regulatory functions (12–14). We have thus far not found evidence that the PAS domain of BvgS perceives specific chemical or physical signals (15), suggesting that it participates

<figure><img src="image_2.png" /><figcaption>**FIG 1** Domain architecture of sensor kinases in the BvgS family. (Left panel) Schematic representation of the BvgS dimer comprising VFT, PAS, and histidine kinase (HK) domains, connected to one another with α-helical linkers 1 and 2. The receiver and HPt domains are omitted for clarity. (Right panel) Approximately 30% of BvgS homologues are devoid of the PAS domain and present an α-helical X linker between the VFT and HK domains. Linkers 1 and 2 are rather highly conserved in size (see Fig. 1 in reference 16 and Fig. 3 in reference 19), while the length of linker X is variable (see Fig. 3 in reference 19). Their sequence conservation is shown in Fig. 9.</figcaption></figure>
only in signal transduction. The PAS domains of BvgS appear to modify their dimeric interface in response to negative modulators. This change in quaternary structure amplifies small changes of conformation and dynamics of upstream linker 1 (16).

Both linker 1 and linker 2 can form two-helix coiled coils. In canonical coiled coils, the two parallel α helices coil around one another in a left-handed manner. They are characterized by heptads of amino acid residues that form two helical turns, denoted by "abcdefg" letters (17). At the central positions of the coiled coil that form the interhelix interface, the "a" and "d" residues are generally hydrophobic and nonaromatic. The surface-exposed residues are preferentially polar and favor helix formation or mediate inter- and intrahelix interactions. Thus, both hydrophobic effects and ionic interactions stabilize coiled coils (17, 18).

Signal transduction in BvgS proceeds by changes of the conformation and the dynamics of the various segments of the protein. The kinase state is characterized by great dynamics of the membrane-distal VFT1 domains, linker 1 adopting a coiled-coil conformation, the PAS domains forming a dimer, and linker 2 being dynamic (16, 19). This dynamic behavior is compatible with kinase activity, which was shown in model sensor kinases to be linked to dynamic, asymmetrical conformations of the DHp and CA domains in the dimer (20–22). Binding of nicotinate to the VFT2 domains causes BvgS to shift to the phosphatase mode (7, 23). The periplasmic portion of BvgS then becomes less mobile and more compact, which induces a small symmetrical piston movement of the TM helices toward the periplasm (7, 16). This leads to the disruption of the linker 1 coiled coil and the splaying out of the last part of the linker 1 helices, which loosens the PAS domain interface. These changes enable linker 2 to adopt a stable coiled-coil conformation, setting the enzymatic moiety in the phosphatase mode (16, 19).

BvgS serves as a prototype for a family of more than 11,000 VFT-containing sensor kinases, some of which are found in major bacterial pathogens (9, 10). A sizeable proportion (~30%) of BvgS homologues are devoid of a PAS domain (19). In these proteins, α helices predicted to form coiled coils directly link the VFT and DHp domains (Fig. 1). In this work, we investigated the regulation of sensor kinases devoid of a PAS domain. We found that their activity is regulated by an antagonism between coiled coils at the N and C termini of the interdomain linker, very much like linkers 1 and 2 of sensor kinases with a PAS domain.

## RESULTS

**BvgS<sub>Δ65</sub>: an eccentric regulation.** Among the sensor kinases homologous to BvgS and comprising two VFT domains, approximately 30% of the predicted proteins have no PAS domain (19). These proteins harbor a single α-helical linker (here called linker X) directly connecting the TM helices to the DHp domain (Fig. 1). We previously generated four chimeric BvgS variants, BvgS<sub>ΔR1</sub>, BvgS<sub>ΔR2</sub>, BvgS<sub>ΔR3</sub>, and BvgS<sub>ΔR4</sub>, in which linker 1, the PAS domain, and linker 2 of BvgS were replaced with the X linkers of four PAS-less BvgS homologues whose lengths differed by multiples of 7 residues (19) (Fig. 2A). The four chimeras displayed kinase activity at the basal state, and the BvgS<sub>ΔR3</sub> and BvgS<sub>ΔR4</sub> chimeras shifted to the phosphatase mode in response to chemical modulation (19). Thus, as BvgS-like regulation can be obtained with linker X sequences, this system may serve as a platform to decipher the rules that govern regulation in the entire family.

Here, we constructed a new chimera called BvgS<sub>Δ65</sub>, whose natural X linker from the *Pseudomonas* sp. strain M1 BvgS homologue is 8 residues longer than that of BvgS<sub>ΔR2</sub> (Fig. 2A). We characterized the activity of BvgS<sub>Δ65</sub> using two reporter systems, the *ptx-lacZ* and *fhaB-lacZ* transcriptional fusions (24) positively regulated by BvgAS and downregulated by chemical modulation. Both promoters are activated at high concentrations of phosphorylated BvgA, but only the *fhaB* promoter can be transactivated at low concentrations that correspond to Bvg<sup>I</sup> conditions (25). At the basal state (i.e., under standard growth conditions), BvgS<sub>Δ65</sub> exhibited hardly any activity using the *ptx-lacZ* reporter and high β-galactoside (β-Gal) activity levels with *fhaB-lacZ*, indicating that the recombinant strain is in an intermediate, Bvg<sup>I</sup> phase (Fig. 2B and C). Addition
<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /><figcaption>**FIG 2** Chimeric BvgS variants with natural X linkers of homologs. (A) Sequences of the variants in which the region between the TM segment and the DHp domain of BvgS was replaced with natural linkers of PAS-less BvgS homologues. The species from which the sequences originated are *Pseudomonas brassicacearum* (accession no. WP_039012684; Pb), *Pseudomonas fragi* (accession no. WP_016782378; Pf), *Pseudomonas* sp. strain M1 (accession no. WP_009621264; Psp), *Pseudomonas putida* (GenBank accession no. AAN69015; Pp), and *Enterobacter cloacae* (accession no. WP_050860054; Ec). The variants obtained were BvgS<sub>ΔR1</sub>, BvgS<sub>ΔR2</sub>, BvgS<sub>Δ65</sub>, BvgS<sub>ΔR3</sub>, and BvgS<sub>ΔR4</sub>, respectively. (B and C) The *ptx-lacZ* (B) and *fhaB-lacZ* (C) reporters were used to determine the basal activities and the responses to 2 mM chloronicotinate (CN 2) or 50 mM MgSO<sub>4</sub> (Mg 50) for the BvgS<sub>Δ65</sub> variant compared to those of wt BvgS (WT). AU, arbitrary units. (D) The *ptx-lacZ* reporter was used to determine the effect of VFT1 domain closing on activity (BvgS<sub>Δ65</sub> VFT1–SS variant). Where indicated, the reducing agent Tris(2-carboxyethyl)phosphine hydrochloride (TCEP) was added to the cultures at 10 mM to reduce the S-S bond in VFT1. The means and standard errors of the means are given.</figcaption></figure>

of chloronicotinate, a potent modulator, to the bacterial cultures caused a shift to β-Gal activity levels similar to those of wild-type (wt) BvgS under standard conditions (i.e., in the kinase mode) for the *ptx-lacZ* reporter (Fig. 2B). No change of *fhaB-lacZ* activity was detected under modulating conditions, most likely because the maximum level of *fhaB* transcription was achieved at the basal state (Fig. 2C). Thus, BvgS<sub>Δ65</sub> regulates *ptx* expression in a manner opposite that seen with wt BvgS, as it shifts to a high-kinase mode of activity upon modulation. However, it is in a low-kinase (rather than phosphatase) mode under basal conditions.

Addition of MgSO₄, another chemical modulator (4, 26), caused the same increase of *ptx-lacZ* activity as chloronicotinate (Fig. 2B). We also tested the closed membrane-distal VFT1 domain of BvgS, as VFT1 may represent the natural ligand-binding sites *in vivo* (10). Disulfide (S-S) bond-mediated closing of VFT1 to mimic a liganded form shifts wt BvgS to phosphatase, similarly to chemical modulation (10). Introduction of the same S-S bond in BvgS<sub>Δ65</sub> restored kinase activity at the basal state with the *ptx-lacZ* reporter, and reduction of the S-S bond (i.e., allowing reopening of VFT1) by the use of Tris(2-carboxyethyl)phosphine hydrochloride (TCEP) markedly reduced the kinase activity (Fig. 2D). Thus, regulation of BvgS<sub>Δ65</sub> is globally inverted relative to that of BvgS. Taken together, our data suggest that the extra residue in linker X of BvgS<sub>Δ65</sub> relative to the heptad periodicity of the BvgS<sub>ΔR1</sub> to BvgS<sub>ΔR4</sub> chimeras displaces the conformational equilibrium away from the kinase mode under basal conditions. Chemical modulation or VFT1 closing shifts BvgS<sub>Δ65</sub> back to the kinase mode.

**Two coiled-coil registers in linker X of BvgS<sub>Δ65</sub>.** Predictions of coiled coils in linker X of BvgS<sub>Δ65</sub> indicated two distinct registers (Fig. 3A). One of them (in gray letters below the sequence in the figure) is best predicted immediately before the DHp domain of the kinase moiety. The “a” and “d” interfacial positions of that register harbor the Leu and Ala residues of the coiled-coil motif “LxxxLxxAxxxAxxA” (here denoted “LLAAA”) found in the linkers 2 of BvgS and its homologs (19). Increased coiled-coil stability resulting from changing the LLAAA motif to LLLLA in linker 2 of BvgS or in the
<figure><img src="image_5.png" /><figcaption>**FIG 3** Organization of the X linker of BvgS<sub>Δ65</sub>. (A) Amino acid sequences of the X linkers of BvgS<sub>Δ65</sub> and the 5-Ala variant. The "a" and "d" positions of the two coiled-coil registers (P and K) are indicated below the sequence. The central "LLAAA" motif is shown in bold letters. (B) Representation of the X linker coiled coils of BvgS<sub>Δ65</sub> along registers P and K using helical wheel diagrams prepared with DrawCoil 1.0 (http://www.grigoryanlab.org/drawcoil/). Hydrophobic residues are in black, polar residues in yellow, negatively charged residues in red, and positively charged residues in blue. Interfacial electrostatic bridges and core-to-interface electrostatic bridges are indicated with dashed lines, blue for favorable ionic interactions and red for unfavorable interactions. (C) The *ptx-lacZ* reporter system was used to determine the activities of the BvgS<sub>Δ65</sub> variants. The means and standard errors of the means are given.</figcaption></figure>

last portion of linker X of the PAS-less BvgS<sub>ΔR1</sub> variant leads to the phosphatase mode (19). Therefore, the coiled-coil register corresponding to the LLAAA motif is denoted here "register P," for the register inducing the phosphatase mode.

In BvgS<sub>Δ65</sub>, another coiled-coil register (in blue letters in Fig. 3A) is predicted in the first part of linker X, but it appears to be less favorable in the last portion. We called it "register K," for the register inducing the kinase mode. Given that the kinase mode for wt BvgS is characterized by linker 2 being dynamic (19), we hypothesized that for BvgS<sub>Δ65</sub>, the modulated state is characterized by coiled-coil formation along register K in the first portion of linker X while keeping the linker region immediately preceding the DHp domain dynamic.

In BvgS<sub>Δ65</sub>, repulsive interhelical interactions are predicted in the first portion of linker X due to the density of positively charged residues (Fig. 3B). Replacing five consecutive positively charged residues (Arg<sup>569</sup> to Arg<sup>573</sup>) with Ala generated BvgS<sub>Δ65-5A</sub>, which is locked in the kinase mode but with a lower maximum activity than that of modulated BvgS<sub>Δ65</sub> (Fig. 3A and C). Among those five residues, the central Lys<sup>571</sup> residue is in the "d" position of register K, which should affect coiled-coil stability in that register (17, 27). Its replacement yielded a variant locked in a mode of high kinase activity (Fig. 3C). Thus, amino acid residues at key coiled-coil positions determine the balance between the kinase and phosphatase modes and their interconversion.

To validate the coiled-coil predictions for linker X of BvgS<sub>Δ65</sub>, cysteine cross-linking-mediated ("Cys-scanning") analyses were performed after the positively charged region. Residues 574 to 595 were individually replaced by Cys, and an oxidative treatment was performed to favor *in vivo* formation of intradimer S-S bonds (Fig. 4A; see also Fig. S1 in the supplemental material). High proportions of S-S bond-mediated cross-links were observed between residues predicted to adopt interfacial "a" and "d" positions of register P, especially in the portion comprising residues 585 to 595. Coiled-coil formation in register P immediately before the DHp domain is consistent with the close-to-
<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /><figcaption>**FIG 4** Topology and dynamics of the X linker of BvgS<sub>Δ65</sub>. (A) Cys-scanning analyses of BvgS<sub>Δ65</sub> were performed under basal conditions (0; black curve) or after the addition of 5 mM chloronicotinate (CN 5; red curve). The proportions of dimers from two independent experiments are graphed, and the means and standard errors of the means are indicated. The “a” and “d” positions of the two coiled-coil registers (Reg P and Reg K) are indicated below amino acid positions. (B and C) The activities of the BvgS variants were determined using the *ptx-lacZ* (B) or *fhaB-lacZ* (C) reporters. The means and standard errors of the means are indicated.</figcaption></figure>

phosphatase (low-kinase) mode of BvgS<sub>Δ65</sub> under basal conditions. In the region encompassing residues 579 to 584, S-S bond patterns indicated a coiled-coil irregularity that might be the site of interconversion between the two registers. Surprisingly, the addition of a modulator did not modify the proportions of interhelical cross-links
formed at most positions. This result might be related to the fact that the protein is not
fully in the phosphatase mode under basal conditions and, therefore, that the cross-
links represent the superimposition of two conformations, rather than a completely
displaced equilibrium. In addition, many Cys substitutions affected BvgSΔ65 activity,
which might also explain the limited effects of modulation on S-S-mediated cross-
linking (Fig. 4B and C). While the regulation trends were at many positions similar to
that seen with BvgSΔ65, i.e., the kinase activity increased upon modulation, the levels of
activity were lower. In addition, several Cys variants appeared to be locked in the kinase
mode (positions 578, 581, 582, 585, and 589) or in the phosphatase mode (positions 593
and 594). We checked the production and stability of the latter two in B. pertussis
membrane extracts to dismiss the possibility that the proteins were proteolyzed *in vivo*.
Both were properly produced (Fig. S2A), suggesting that charged residues at positions
593 and 594 are important for BvgSΔ65 function. Of note, charged residues were also
found at the corresponding positions in variants BvgSΔR1 to BvgSΔR4 (Fig. 2A), as is often
the case in the members of that family (see below). The observation that replacements
of interfacial and noninterfacial residues affect regulation indicates that coiled-coil
stability in each portion of linker X, which is dictated by the nature and the arrange-
ment of both interfacial and surface residues, is critical for function.

**Effect of linker length on regulation.** To restore the heptad periodicity to linker X of BvgSΔ65 as found in chimeras BvgSΔR1 to BvgSΔR4, we deleted one residue immediately after the TM segment, yielding BvgSΔ64. We similarly constructed size variants BvgSΔ62, BvgSΔ63, and BvgSΔ66 and tested the effects of those modifications on activity (Fig. 5). BvgSΔ64 exhibited regulation similar to that seen with wt BvgS; i.e., it was in the kinase and phosphatase modes of activity under basal conditions and in response to modulation, respectively (Fig. 5B). BvgSΔ63 was locked in the kinase mode, while both BvgSΔ62 and BvgSΔ66 were locked in the phosphatase mode, despite normal BvgS production levels (Fig. 5B and C; see also Fig. S2B). This indicates that the length and the organization of the helices are critical for function of the chimeras. In the context of BvgS, the only functional lengths are those found in BvgSΔ64 and BvgSΔ65, while the other chimeras are locked in a single state.

**Effect of linker composition on regulation.** We constructed new chimeras har-
boring natural linkers of the same lengths as that of BvgSΔ64 or BvgSΔ65 (Fig. 5A). The
BvgSΔ64 *Ec* variant was locked in the kinase mode (Fig. 5B). It notably features a
degenerate IQSAA interfacial motif in register P, unlikely to favor coiled-coil formation
before the DHp domain. The BvgSΔ65 *Se* variant exhibited low kinase activity at the basal
state, which increased upon modulation, although only with the *fhaB-lacZ* reporter
(Fig. 5B and C). It is thus regulated, in an inverted manner, between a phosphatase
mode under basal conditions and an intermediate mode of activity under modulating
conditions. Another chimera with the same linker size, BvgSΔ65 *Sy*, was locked in an
intermediate mode of activity (Fig. 5B and C). The X linkers of BvgSΔ65 *Sy* and BvgSΔ65 *Se*
harbored interfacial motif RLAAA in register P before the DHp domain (Fig. 5A).

For BvgS<sub>Δ65</sub> Se and BvgS<sub>Δ65</sub> Sy, we deleted one residue to generate variants with linker lengths similar to that seen in BvgS<sub>Δ64</sub>. This deletion locked BvgS<sub>Δ64</sub> Se in kinase mode, but BvgS<sub>Δ64</sub> Sy remained in the same intermediate state as BvgS<sub>Δ65</sub> Sy (Fig. 5B and C). Thus, not all linker X sequences from homologues function properly when inserted between the TM segment and the DHp domain of BvgS. A notable difference between BvgS<sub>Δ65</sub> and the BvgS<sub>Δ65</sub> Se and BvgS<sub>Δ65</sub> Sy variants is that the α helices of the former might be stabilized by a number of potential intrahelix ionic bonds between i and i+3 residues or between i and i+4 residues at positions “b,” “c,” and “f” of the coil (28, 29) (Fig. S3). In contrast, the X linkers of the other two chimeras have fewer predicted intrahelix interactions.

**Search for functional variants with longer linkers.** We constructed chimeras with X linkers longer than that of BvgSΔ65 and screened for functional sequences. In addition, we modified their length to obtain the same heptad periodicities as in BvgSΔ65 and BvgSΔ64 (Fig. 6A).
A

<figure><img src="image_9.png" /></figure>

B

<figure><img src="image_10.png" /></figure>

C

<figure><img src="image_11.png" /><figcaption>**FIG 5** Effect of linker length and composition on activity. (A) Sequences of the X linkers of the size variants of BvgS<sub>Δ65</sub> (BvgS<sub>Δ62r</sub>, BvgS<sub>Δ63</sub>, BvgS<sub>Δ64</sub>, and BvgS<sub>Δ66</sub>) and of other variants are indicated. The species from which the sequences originated were *Enterobacter cancerogenus* (accession no. WP_058610751; Ec), *Selenomonas* sp. strain ND2010 (accession no. WP_081823318; Se), and *Synergistes* sp. strain 3_1_syn1 (accession no. WP_083829553; Sy). The "a" and "d" positions of the two coiled-coil registers (Reg P and Reg K) are indicated below the alignment. (B and C) The activities of the BvgS variants were determined using the *ptx-lacZ* (B) and *fhaB-lacZ* (C) reporters. The means and standard errors of the means are indicated.</figcaption></figure>

BvgS<sub>Δ70A</sub>, BvgS<sub>Δ71A</sub>, and BvgS<sub>Δ72A</sub> were all locked in the phosphatase mode, although they were all produced at levels similar to the BvgS level (Fig. 6B and C; see also Fig. S2C). BvgS<sub>Δ70B</sub>, although produced at normal levels, was also locked in the phosphatase mode, but the longer BvgS<sub>Δ71B</sub> and BvgS<sub>Δ72B</sub> variants were locked in the kinase mode (Fig. 6B and C; see also Fig. S2C). Thus, the linker length affects the phenotype, but modulation cannot regulate the activity of any of the three chimeras.

In contrast, BvgS<sub>Δ72C</sub>, with an additional full heptad relative to BvgS<sub>Δ65</sub>, showed moderate kinase activity under basal conditions that increased after addition of the modulator. The shorter counterparts, BvgS<sub>Δ70C</sub> and BvgS<sub>Δ71C</sub>, displayed kinase-locked phenotypes. The partially inverted regulation of BvgS<sub>Δ72C</sub> is somewhat similar to that of BvgS<sub>Δ65</sub>. Specific features of this sequence include the presence of several predicted
A

<figure><img src="image_12.png" /></figure>

B

<figure><img src="image_13.png" /></figure>

C

<figure><img src="image_14.png" /><figcaption>**FIG 6** Activity of chimeras with longer linker. (A) Sequences of the X linkers of chimeric variants. The species from which the sequences originated were *Treponema primitia* (accession no. WP_015709136; Δ70A), *Citrobacter* sp. strain MGH104 (accession no. WP_048232642; Δ70B), and *Butytrivibrio* sp. strain VCB2001 (accession no. WP_081669071; Δ70C). The “a” and “d” positions of the two coiled-coil registers (Reg P and Reg K) are indicated below the alignment. (B and C) The activities of the BvgS<sub>chim</sub> variants were determined using the *ptx-lacZ* (B) and *fhaB-lacZ* (C) reporters. The means and standard errors of the means are indicated.</figcaption></figure>

interhelical ionic bonds, in particular, in the first portion of the linker in register K and
in an interfacial motif, KLAAA, in register P (Fig. 6A).

**Molecular dynamics simulations.** We generated *in silico* models of the linkers described above along the two registers and performed molecular dynamics (MD) simulations for each. The simulations were analyzed to determine the features of the linkers that are functional or partially functional in the context of our BvgS platform, i.e., that cause shifts of activity upon modulation. For the BvgS<sub>Δ65</sub>, BvgS<sub>Δ65 Se'</sub> and BvgS<sub>Δ70C</sub> linker sequences, the levels of helical stability notably differed between the two sets of simulations (Fig. 7). Thus, the α helices along register K remained mostly stable over time, while the linkers tended to form P<sub>1</sub> helices in register P, mostly at their N-terminal side. For the phenotypically locked BvgS<sub>65-5A'</sub>, BvgS<sub>Δ70A'</sub>, and BvgS<sub>Δ70B</sub> variants, in contrast, the MD simulation profiles were almost identical to one another along both registers. MD simulations of the linker of the BvgS<sub>Δ65 Sy</sub> variant, locked in the intermediate phase, showed strong disruption of helicity over the entire length of the coiled-coil structure in register P and a somewhat less severe disruption in register K. From these data, it appears that the stability of the α helices that form the X linkers contributes to function. In addition, differential levels of α-helical stability between registers K and P enable the interconversion between two states. In contrast, locked phenotypes appear to correspond to low overall helical stability or to similar helical stabilities in the two registers.

**Surface and core residues determine regulation.** On the basis of all the information gathered on those systems, we selected a chimera to use for sequence modifications aimed at altering its regulation properties in a controlled manner. BvgS<sub>Δ72C</sub> was chosen because its length matches that of BvgS<sub>Δ65</sub> but with an additional full heptad. In addition, it appears to be weakly upmodulated by chloronicotinate (Fig. 6B). We thus tested our hypotheses of regulation by modifying its composition in a rational manner.
<figure><img src="image_15.png" /><figcaption>**FIG 7** Helix length as calculated from the MD simulations. Averages and standard deviations were calculated for both helices of the coiled coil and for all extracted time frames. Values are shown for each sequence in its assigned registers, namely, registers K (gray) and P (brown). The square data points corresponding to the value of the averages are sized proportionally to the size of the primary (i.e., largest) structural cluster, the largest value being 96.9% for BvgS<sub>Δ70C</sub> in register K and the smallest value being 24.7% for BvgS<sub>Δ65S</sub> in register P. The representative (i.e., "center") structures of these clusters are shown above the graph, colored red and green for registers K and P, respectively.</figcaption></figure>

To favor the phosphatase mode, we modified the KLAAA motif by consecutively replacing the Lys and first Ala residues by Leu residues in "d" position along register P to increase coiled-coil stability before the DHp domain (Fig. 8A). BvgS<sub>Δ72C LLAAA</sub> was in the phosphatase mode at the basal state (Fig. 8B and C), and addition of the modulator increased its kinase activity with the *fhaB-lacZ* reporter system (Fig. 8C). Further strengthening the coiled coil by generating the LLLAA motif yielded a bona fide "inverted-regulation" variant, BvgS<sub>Δ72C LLLAA</sub>, whose phenotype paralleled that of BvgS<sub>Δ65</sub> (Fig. 8B and C). Its linker length corresponded to that of the latter plus one full heptad.

<figure><img src="image_16.png" /></figure>

<figure><img src="image_17.png" /><figcaption>**FIG 8** Modifications of the effects of coiled-coil registers P and K on regulation. (A) Sequences of BvgS<sub>Δ72C</sub> variants where Leu substitutions were progressively introduced to increase coiled-coil stability along register P and of BvgS<sub>Δ72C LLLAA</sub> variants in which specific amino acids in the "e" and "g" positions or the "b," "c," and "f" positions were replaced (denoted in green) to disrupt interhelix interactions or to favor intrahelix interactions. The "a" and "d" positions of the two coiled-coil registers (Reg P and Reg K) are indicated below the alignment. (B and C) The activities of the BvgS<sub>Δ72C</sub> variants were determined using the *ptx-lacZ* (B) and *fhaB-lacZ* (C) reporters. The means and standard errors of the means are indicated.</figcaption></figure>
Thus, one can restore regulation by strengthening register P. We tried to obtain the same effect by altering the other register, starting from $BvgS_{Δ72C LLAAA}$. In particular, we modified the inter- or intrahelix interactions in the first segment of its linker X by introducing mutations that exclusively affect residues at noninterfacial positions b, c, e, f, and g of register K. The combined substitutions $K_{571}E$ and $E_{578}R$ in $BvgS_{Δ72C LLAAA}$, designed to weaken interhelix interactions in register K (Fig. 8A; see also Fig. S4), yielded a variant, $BvgS_{Δ72C LLAAA Inter}$, locked in the phosphatase mode (Fig. 8C). It was produced at normal levels (Fig. S2D). Thus, as we predicted, decreasing coiled-coil stability in register K by altering noninterfacial positions favors register P, with the LLAAA motif in the latter register taking precedence.

The substitutions $S_{575}E$, $A_{579}R$, $K_{582}E$, $I_{589}E$, and $I_{593}R$ were combined to strengthen the α helices by increasing the number of potential intrahelix ionic interactions (Fig. 8A; see also Fig. S4). The resulting $BvgS_{Δ72C LLAAA Intra}$ variant had low kinase activity at the basal state with the *ptx-lacZ* reporter that markedly increased upon modulation (Fig. 8B). The *fhaB-lacZ* reporter system indicated that this variant was in the BvgI phase at the basal state (Fig. 8C). By increasing helix stability, one could thus obtain the same regulation phenotype as would be obtained by strengthening the coiled coil in register P. This demonstrates that surface-exposed, noninterfacial residues of linker X also contribute to function by determining the intrinsic stability of the helices.

**Conserved elements involved in sensor kinase regulation.** Last, we analyzed a large set of linker sequences encompassing the variety found in the BvgS family to determine whether our findings could be generalized. More than 8,000 sequences of two-component sensor kinases with two VFT domains were collected, including 2,479 proteins without a PAS domain and 5,009 containing a single PAS domain. For the PAS-less proteins, 740 unique sequences of linker X were retrieved and clustered at 80% sequence identity. The first 30 clusters, which are the most highly populated and comprise 393 sequences, were further analyzed. Clusters were grouped depending on heptad composition, yielding four "superclusters" with 214, 53, 39, and 31 sequences, respectively. The sizes of the X linkers range from two to six heptads in this set. Occurrences of longer linkers (up to 13 heptads) in the initial set of sequences were very rare. Linker X of $BvgS_{Δ65}$ belonged to the largest supercluster.

Sequence alignments showed a high level of conservation in the two heptads that precede the DHp domain (Fig. 9). The LLAAA motif was present in all superclusters, with the exception that the first Leu was not found in the linker sequences harboring only three or four heptads. In addition, a few charged or polar residues were relatively well conserved at non-"a" and non-"d" positions in the last heptads. The conserved features are most likely important for the control of the enzymatic moiety that follows immediately.

In the first heptads of the X linkers, a high density of charged residues was generally conserved, as were a few residues at both interfacial and noninterfacial positions of register K. In particular, a positively charged residue ($K^{571}$ in $BvgS_{Δ65}$ and $R^{570}$ in wt BvgS) shown to be critical for regulation was very often found at the "d" position in the second heptad of the X linkers composed of at least 4 heptads. Overall, conserved features of linker X in the family support our model that its composition generates the coiled-coil antagonism that regulates the output activity of the system.

As for the PAS domain-containing sensor kinases, most harbored four heptads in linker 1 and three in linker 2 (Fig. 9). The sequences of the linkers 2 aligned quite well with those of the last three heptads of the X linkers, except for very short linker sequences (19). Similarly, the sequences of the first three heptads of the linkers 1 aligned relatively well with those of the long X linkers. In contrast, the fourth heptads of the linkers 1 had conserved sequences with no corresponding heptad in the X linkers. This last heptad is most likely necessary to control the PAS domain.

## DISCUSSION

BvgS serves as the prototype for a family of TCS sensor kinases with VFT domains. In BvgS, a cytoplasmic PAS domain flanked by two helical linkers connects the VFT and DHp domains. In a sizeable fraction of BvgS homologues, however, the PAS domain is
<figure><img src="image_18.png" /><figcaption>**FIG 9** Heptad composition of the linkers in the BvgS family. The sequence alignments of the region between the predicted TM segment and the DHp domain of sensor kinases of the BvgS family are displayed. Each line represents the consensus sequence for one cluster, with the heptads boxed in colors. The first four blocks of sequences (4 superclusters) represent 43% (322 of 740 sequences) of the X linkers of PAS-less proteins. The last block (11 clusters) represents 66.5% (510 of 767 sequences) of linkers 1 and 2 of sensor kinases with one PAS domain. The two registers (registers K and P) are shown in two lines between the two groups of proteins. A large gap, which corresponds to the PAS domain in the second group, was inserted between the first three and the last three heptads of the PAS-less proteins. Note that the fourth heptad (boxed in black) of the proteins with a PAS domain has no corresponding heptad in the PAS-less proteins. BvgS belongs to the sequence cluster indicated by an asterisk.</figcaption></figure>

replaced with a single helical linker between the periplasmic and cytoplasmic parts. In previous work, we identified rules for BvgS regulation (16, 19). Briefly, under basal conditions, linker 1 forms a coiled coil whereas linker 2 is dynamic and BvgS is in the kinase mode. Upon modulation, disruption of the linker 1 coiled coil affects the PAS domain interface, which enables linker 2 to adopt a stable coiled-coil conformation leading to the phosphatase mode (16, 19). In this work, we extended our analyses to the subset of homologues devoid of the PAS domain. We found that linker X between the TM and DHp domains harbors two antagonistic coiled-coil registers. Regulation of activity appears to involve coiled-coil interconversion between these two registers. Therefore, function is governed by similar mechanical principles in the entire family, with specific sequence adaptations in each type of sensor kinase.

In linker X, coiled-coil register P, with the characteristic LLAAA motif also found in linker 2 of PAS-containing family members, has been extensively predicted to be present just before the DHp domain. Coiled-coil formation along that register leads to phosphatase activity, for both PAS-containing and PAS-less sensor kinases (19; this work). In the first moiety of linker X, a moderately stable coiled coil is predicted along register K, with a suboptimal interfacial motif that is compensated for by favorable
inter- and intrahelical ionic bonds. Removing interhelical ionic interactions along register K in that region favors the phosphatase mode. Conversely, enhancing coiled-coil formation along register K by strengthening the interface favors the kinase mode, as shown by the kinase-locked phenotypes of the Lys<sup>571</sup>-to-Ala variant in linker X of BvgS<sub>Δ65</sub> (this work) or the Arg<sup>570</sup>-to-Ala variant in linker 1 of BvgS (19). A charged residue at this position is highly conserved in the family.

We propose that the antagonism between the two coils in the PAS-less BvgS homologues forms the basis for regulation: coiled-coil formation along register K prevents coiled-coil formation along register P and vice versa, and signal perception triggers interconversion. The dynamics of the last portion of linker X is conducive to kinase activity, which has been proposed to require dynamic asymmetry of the DHp helices in model sensor kinases, while the phosphatase state corresponds to a more rigid state (21, 22). By affecting the TM helices, modulation displaces the equilibrium in one direction or the other in a manner dependent on the initial state of the system, which is itself determined by the length of linker X. Such a mechanism of interconversion between two states of coiled coils triggered by external signals has been reported for Nek2, a eukaryotic kinase involved in cell cycle regulation (30). Other bacterial His kinases appear to function in a similar manner (31). Successive domains are connected by two-helix coiled-coil linkers, and a register mismatch between the coiled coil of the linker and that of the preceding or following domain allows conversion between states. Similarly, in the Tar and Tsr chemoreceptors, a "control cable" modulates the intensity of the register mismatch between the TM2 helix and the following HAMP domain (32).

The BvgS<sub>Δ65</sub> chimera displays inverted regulation; i.e., modulators cause a strong increase of kinase activity. In wt BvgS, modulator binding to the periplasmic domain causes a small piston motion that pulls one residue of linker 1 into the membrane, thus disrupting the first coiled coil (16). In BvgS<sub>Δ65</sub>, modulation is expected to cause a similar motion that shortens the cytoplasmic length of linker X by one residue, in this case leading to the kinase mode. Another way to bring about the kinase mode is to artificially remove one residue from linker X of BvgS<sub>Δ65</sub>, such as in BvgS<sub>Δ64</sub>. Not coincidentally, the length of BvgS<sub>Δ64</sub> linker X differs by full heptad counts from those of the chimeras of BvgS<sub>ΔR1</sub> to BvgS<sub>ΔR4</sub>, which are also in kinase mode under basal conditions (19).

The length of linker X is thus a critical determinant of function. Thus, one-residue additions or deletions relative to the functional lengths affect coiled-coil packing, with functional consequences. In addition, the ability to undergo interconversion is determined by the composition of the linker X helices, at both interfacial and noninterfacial coiled-coil positions. Local helical propensities and charges have recently been shown to modulate coil stability in a model system (33). The α helices of BvgS<sub>Δ65</sub> linker X harbor a number of potential intrahelical ionic bonds that contribute to their stability. Similarly, we were able to recover regulatory activity of BvgS<sub>Δ72C</sub> LLAAA by increasing helix stability. Taking the data together, responsiveness to modulation thus requires marginally stable coils in the two registers and sufficiently pronounced helical stability to avoid the status of being trapped in a single state. According to this model, phosphatase-locked phenotypes might be explained by a loss of helicity, causing the register P coil in the C-terminal portion to prevail. Intermediate (i.e., low-kinase) locked phenotypes might be due to the accommodation of the two coiled-coil registers by loose α helices (34), resulting in suboptimal conformations.

The X linkers that caused a loss of function in BvgS are most likely functional in their natural sensor kinases. For those sensor kinases, signal perception might trigger conformational changes of the VFT and TM domains that are larger than those seen in BvgS, thus allowing register interconversion and activity regulation. It is also possible that the composition of the cytoplasmic membrane in their host organisms differs from that of *B. pertussis*. Differences in membrane thickness or in lipid composition are likely to affect the length of the TM segment, its motions, and its interactions with the polar phospholipid heads, with a possible impact on regulation.

The major finding of this work is that the N-terminal and C-terminal moieties of linker X functionally correspond to linker 1 and linker 2 of sensor kinases with PAS
domains, respectively. The first heptads of the X linkers, similarly to the results seen with linker 1, directly follow the TM segment; they "process" the input signal, which consists of changes of the conformation and dynamics of the VFT domains and the ensuing motions of the TM helices. The C-terminal heptads of the X linkers, as seen with linker 2, are involved in the output of the system, i.e., in the control of enzymatic activity. In sensor kinases without a PAS domain, a balanced antagonism between the upstream and downstream registers of linker X is necessary for interconversion between the two states. Sensor kinases with a PAS domain regulate enzymatic activity with similar antagonism between the coiled coils formed by linkers 1 and 2 (16, 19). The sequences of the fourth heptad of linker 1 of the PAS-containing sensor kinases are relatively well conserved in that subfamily, but there is no corresponding heptad in the PAS-less proteins. It is likely that this heptad directly preceding the PAS domain serves to control the quaternary structure of the latter in response to signals. The PAS domain thus arbitrates between the two linkers, serving as a signal amplifier. Therefore, globally similar principles of mechanical signal transduction are most likely valid for the entire family, with specific adaptations for the presence or the absence of an intervening PAS domain.

## MATERIALS AND METHODS

**Strains, plasmids, and culture conditions.** *B. pertussis* was grown on Bordet-Gengou agar plates for 2 days at 37°C and then cultured in modified stainer Scholte (SS) liquid medium at 37°C under conditions of agitation. To construct BvgS chimeras devoid of the PAS domain, synthetic gene portions (Genecust, Luxembourg) coding for the junction between the TM segment and the DHp domain of BvgS homologues were introduced by BglII-XbaI cassette exchange in pUC19mpla (10) to replace the wt BvgS fragment. Then, the EcoRI-HindIII fragments of the resulting plasmids were transferred into pBBR1-MCS4, yielding pBBRmpla variants as previously described (19). To introduce the two Cys codons used to generate an S-S bond between the VFT1 lobes, the Spel-BamHI cassette of the pUC19mpla derivative coding for BvgSΔΔ65 was replaced by that extracted from pUC19mosaic containing the Glu¹¹³Cys and Asn¹¹⁷Cys mutations (10). Variants with short deletions or insertions were constructed by site-directed mutagenesis on a pUC19-derived plasmid containing the relevant DNA region, followed by cassette exchange in pUC19mpla and then in pBBRmpla. The recombinant *B. pertussis* strains were obtained by introducing the pBBRmpla variants by conjugation into BPSM<sub>newΔAS</sub> carrying chromosomal *ptx-lacZ* or *fhaB-lacZ* transcriptional fusions (10). The *bvgSΔ65* variants used for Cys-scanning analyses in *Escherichia coli* were constructed by mutagenesis on a pUC19-derived plasmid carrying the sequence coding for the transmembrane segment of BvgS, the linker of *Pseudomonas* sp. strain M1 (GenBank accession no. WP_009621264), and the first part of the BvgS kinase. Cassette exchanges were performed in pPORVPH (19). To check the effect of Cys substitutions on BvgS activity, the pBBRmpla<sub>Δ65</sub> Cys variants were prepared by cassette exchange as previously described (19).

**Cys-scanning analyses.** The disulfide (S-S)-mediated cysteine cross-linking-mediated (Cys-scanning) analyses were performed in *E. coli* UT5600 carrying pPOR<sub>Δ65</sub> variants that code for truncated BvgSΔ65 Cys variants (called BvgS<sub>Δ65</sub><sup>f</sup>), encompassing the VFT, TM, linker X, and HK domains as previously described (19). In this case, however, the variants harbored their native Cys<sup>BB1</sup>, since this residue was not involved in S-S bond formation. The Cys-scanning analyses were performed at least two times for each position, and the results were reproducible. Full-length BvgSΔ65 Cys variants (called BvgS<sub>Δ65</sub><sup>f</sup>) were used to determine the effect of Cys substitution on BvgS activity.

**Molecular dynamics simulations.** Models of the coiled coils were prepared in registers K and P using CCBUILDER (35). The models were processed for molecular dynamics (MD) simulations using the Gromacs suite of programs, version 5.0.4, and the Gromos 43a2 force field (36). C-terminal and N-terminal helix ends were neutralized by capping them with an acetyl group and NH₂, respectively. The systems were solvated with simple point charge (SPC) water (37) in a rhombic dodecahedron box with a minimum of 1 nm between any atom of the peptides and the box edge and were neutralized by adding the appropriate amount of Na⁺ or Cl⁻ counterions. The systems were then relaxed by performing 2,000 steps of steepest descent energy minimization followed by a 100-ps MD simulation with positional restraints on the heavy atoms, using an integration time step of 2 fs and a neighbor list update performed every 5 steps. For the MD simulations, a temperature bath of 300 K with a coupling constant of 0.1 ps was employed (38), using independent coupling for the coiled coil and the solvent. Pressure was maintained at 1 bar using isotropic pressure coupling with a coupling constant of 1 ps. Bond lengths were constrained to their equilibrium values with LinCS (39). van der Waals interactions were cut off at a distance of 1 nm, and electrostatic interactions were calculated with the particle mesh Ewald method (40). Equations of motion for the water molecules were solved with the SETTLE algorithm (41). The production runs, consisting of 3 independent runs of 500 ns each, used the same settings but employed a time step of 4 fs, enabled by the use of virtual interaction sites (42). Data extraction was done from the concatenated production runs at an interval of 100 ps. Helix length was calculated from the values corresponding to the instantaneous helical rise per residue and averaged over all time frames. Clustering of structures was performed on the basis of the root mean square deviation (RMSD) matrix, using an algorithm described in reference 43 and an RMSD cutoff of 0.3 nm. The center structure of the most highly populated cluster was used for the figure.
<figure><img src="image_19.png" /></figure>

**Protein sequence analyses.** hmmsearch output files consisting of all predicted proteins identified from the NCBI nonredundant database (release of September 2017) containing both PF00497 (Bacterial extracellular solute-binding protein, family 3) and PF00512 (His kinase A domain) domains were parsed, considering only hits with an E value of <0.001 for both domains. Using a customized Python script, we selected the sequences containing two VFT domains, which resulted in selection of 8,218 TCS sequences. They were divided in two groups, those without PAS domain (2,479 sequences) and those with a single PAS domain (5,009 sequences), using modified versions of the script mentioned above. TM domain predictions were done for all those sequences, and a few sequences were discarded because of the absence of a clearly predicted TM segment.

Sequences devoid of a PAS domain were then reduced to a fragment starting from the end of the predicted TM segment and extending to the end of the HisKA motif, yielding 740 unique sequences of linker X due to high redundancy of this part of the sequences. Sequences containing a single PAS domain were similarly reduced to the end of the HisKA motif, yielding 767 unique sequences, after discarding atypical sequences in which linker 1 was longer than 34 residues. From this set, the linker 1 sequences were extracted for the alignment from the end of the predicted TM segment to the last residue before the PNP motif that initiates the PAS domain. Similarly, the linker 2 sequences included the residues found after the DIT (Asp-Ile-Thr) motif that terminates the PAS domain and before the beginning of the DHp domain.

The two sets of sequences were clustered using CD-HIT (http://weizhongli-lab.org/cdhit s_uite), with an 80% cutoff for the linker X sequences and a 50% cutoff for the linker 1 sequences. The 30 most highly populated clusters of linker X sequences and the 15 most highly populated clusters of linker 1 sequences were further analyzed. Sequences were divided into heptads along register P based on a prior coiled-coil prediction using PCOILS (https://toolkit.tuebingen.mpg.de). Heptad-based manual alignment of different clusters was carried out using Jalview (44), generating four superclusters for the linker X sequences and one for the linker 1 and linker 2 sequences.

**Other methods.** β-Galactosidase assays were performed as described previously (10) with three different clones at different times, and the means and standard errors of the means were determined. The chimeras that gave intermediate-phase phenotypes, i.e., high reporter activity with the *fhaB-lacZ* fusion but no activity with the *ptx-lacZ* fusion, were described as having low kinase activity, while the chimeras yielding no reporter activity with either fusion were described as being in the phosphatase state.

For the detection of inactive BvgS variants in *B. pertussis*, the bacteria were lysed, and the membrane proteins were harvested by ultracentrifugation before immunodetection of BvgS was performed as previously described (19). For the Cys variants, to avoid S-S bond formation during sample handling, 10 mM N-ethyl maleimide was added to the resuspended pellet before lysis. Electrophoresis and immunoblotting were performed as previously described (19).

## SUPPLEMENTAL MATERIAL

Supplemental material for this article may be found at https://doi.org/10.1128/mBio .02052-17.

FIG S1, TIF file, 0.8 MB.

FIG S2, TIF file, 0.8 MB.

FIG S3, TIF file, 1.5 MB.

FIG S4, TIF file, 2.3 MB.

## ACKNOWLEDGMENTS

We thank Eva-Maria Krammer for her initial modeling work.

E.L. received a doctoral fellowship from the Region Nord-Pas-de-Calais and INSERM. This work was supported by the Agence Nationale de la Recherche (grant number ANR-13-BSV8-0002-01 to F.J.-D.). The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

## REFERENCES

1. Gao R, Stock AM. 2009. Biological insights from structures of two-component proteins. Annu Rev Microbiol 63:133-154. https://doi.org/10.1146/annurev.micro.091208.073214.

2. Casino P, Rubio V, Marina A. 2010. The mechanism of signal transduction by two-component systems. Curr Opin Struct Biol 20:763-771. https:// doi.org/10.1016/j.sbi.2010.09.010.

3. Zschiedrich CP, Keidel V, Szurmant H. 2016. Molecular mechanisms of two-component signal transduction. J Mol Biol 428:3752-3775. https:// doi.org/10.1016/j.jmb.2016.08.003.

4. Lacey BW. 1960. Antigenic modulation of *Bordetella pertussis*. J Hyg 8:57-93. https://doi.org/10.1017/S0022172400038134.

5. Melton AR, Weiss AA. 1993. Characterization of environmental regulators of *Bordetella pertussis*. Infect Immun 61:807-815.

6. Nakamura MM, Liew SY, Cummings CA, Brinig MM, Dieterich C, Relman DA.

2006. Growth phase- and nutrient limitation-associated transcript abundance regulation in Bordetella pertussis. Infect Immun 74:5537-5548. https://doi.org/10.1128/IAI.00781-06.

7. Dupré E, Lesne E, Guérin J, Lensink MF, Verger A, de Ruyck J, Brysbaert G, Vezin H, Locht C, Antoine R, Jacob-Dubuisson F. 2015. Signal transduction by BvgS sensor-kinase: binding of modulator nicotinate affects conformation and dynamics of entire periplasmic moiety. J Biol Chem 290:23307-23319. https://doi.org/10.1074/jbc.M115.655720.

8. Deora R, Bootsma HJ, Miller JF, Cotter PA. 2001. Diversity in the Bordetella virulence regulon: transcriptional control of a Bvg-intermediate phase gene. Mol Microbiol 40:669-683. https://doi.org/10.1046/j.1365-2958.2001.02415.x.

9. Jacob-Dubuisson F, Wintjens R, Herrou J, Dupré E, Antoine R. 2012. BvgS of pathogenic Bordetellae: a paradigm for sensor kinase with Venus
Flytrap perception domains, p 57–83. In *Gross R, Beier D (ed), Two-component system in bacteria*. Caister Academic Press, Norfolk, United Kingdom.

10. Dupré E, Herrou J, Lensink MF, Wintjens R, Vagin A, Lebedev A, Crosson S, Villeret V, Locht C, Antoine R, Jacob-Dubuisson F. 2015. Virulence regulation with Venus flytrap domains: structure and function of the periplasmic moiety of the sensor-kinase BvgS. PLoS Pathog 11:e1004700. https://doi.org/10.1371/journal.ppat.1004700.

11. Hoch JA. 2000. Two-component and phosphorelay signal transduction. Curr Opin Microbiol 3:165–170. https://doi.org/10.1016/S1369-5274(00)00070-9.

12. Taylor BL, Zhulin IB. 1999. PAS domains: internal sensors of oxygen, redox potential, and light. Microbiol Mol Biol Rev 63:479–506.

13. Möglich A, Ayers RA, Moffat K. 2009. Structure and signaling mechanism of Per-ARNT-Sim domains. Structure 17:1282–1294. https://doi.org/10.1016/j.str.2009.08.011.

14. Henry JT, Crosson S. 2011. Ligand-binding PAS domains in a genomic, cellular, and structural context. Annu Rev Microbiol 65:261–286. https://doi.org/10.1146/annurev-micro-121809-151631.

15. Dupré E, Wohlkonig A, Herrou J, Locht C, Jacob-Dubuisson F, Antoine R. 2013. Characterization of the PAS domain in the sensor-kinase BvgS: mechanical role in signal transmission. BMC Microbiol 13:172. https://doi.org/10.1186/1471-2180-13-172.

16. Lesne E, Dupré E, Locht C, Antoine R, Jacob-Dubuisson F. 2017. Conformational changes of an interdomain linker mediate mechanical signal transmission in sensor kinase BvgS. J Bacteriol 199:e00114-17. https://doi.org/10.1128/JB.00114-17.

17. Lupas AN, Gruber M. 2005. The structure of alpha-helical coiled coils. Adv Protein Chem 70:37–78. https://doi.org/10.1016/S0065-3233(05)70003-6.

18. Mason JM, Arndt KM. 2004. Coiled coil domains: stability, specificity, and biological implications. Chembiochem 5:170–176. https://doi.org/10.1002/cbic.200300781.

19. Lesne E, Krammer EM, Dupre E, Locht C, Lensink MF, Antoine R, Jacob-Dubuisson F. 2016. Balance between coiled-coil stability and dynamics regulates activity of BvgS sensor kinase in Bordetella. mBio 7:e02089. https://doi.org/10.1128/mBio.02089-15.

20. Ferris HU, Coles M, Lupas AN, Hartmann MD. 2014. Crystallographic snapshot of the Escherichia coli EnvZ histidine kinase in an active conformation. J Struct Biol 186:376–379. https://doi.org/10.1016/j.jsb.2014.03.014.

21. Mechaly AE, Sassoon N, Betton JM, Alzari PM. 2014. Segmental helical motions and dynamical asymmetry modulate histidine kinase autophosphorylation. PLoS Biol 12:e1001776. https://doi.org/10.1371/journal.pbio.1001776.

22. Mechaly AE, Soto Diaz S, Sassoon N, Buschiazzo A, Betton JM, Alzari PM. 2017. Structural coupling between autokinase and phosphotransferase reactions in a bacterial histidine kinase. Structure 25:939–944.e3. https://doi.org/10.1016/j.str.2017.04.011.

23. Herrou J, Bompard C, Wintjens R, Dupré E, Willery E, Villeret V, Locht C, Antoine R, Jacob-Dubuisson F. 2010. Periplasmic domain of the sensor-kinase BvgS reveals a new paradigm for the Venus flytrap mechanism. Proc Natl Acad Sci U S A 107:17351–17355. https://doi.org/10.1073/pnas.1006267107.

24. Antoine R, Alonso S, Raze D, Coutte L, Lesjean S, Willery E, Locht C, Jacob-Dubuisson F. 2000. New virulence-activated and virulence-repressed genes identified by systematic gene inactivation and generation of transcriptional fusions in *Bordetella pertussis*. J Bacteriol 182: 5902–5905. https://doi.org/10.1128/JB.182.20.5902-5905.2000.

25. Jones AM, Boucher PE, Williams CL, Stibitz S, Cotter PA. 2005. Role of BvgA phosphorylation and DNA binding affinity in control of Bvg-mediated phenotypic phase transition in Bordetella pertussis. Mol Microbiol 58:700–713. https://doi.org/10.1111/j.1365-2958.2005.04875.x.

26. Boulanger A, Chen Q, Hinton DM, Stibitz S. 2013. In vivo phosphorylation dynamics of the Bordetella pertussis virulence-controlling response reg-

ulator BvgA. Mol Microbiol 88:156–172. https://doi.org/10.1111/mmi.12177.

27. Cohen C, Parry DAD. 1986. Alpha-helical coiled coil-a widespread motif in proteins. Trends Biochem Sci 11:245–248. https://doi.org/10.1016/0968-0004(86)90186-6.

28. Burkhard P, Ivaninskii S, Lustig A. 2002. Improving coiled-coil stability by optimizing ionic interactions. J Mol Biol 318:901–910. https://doi.org/10.1016/S0022-2836(02)00114-6.

29. Meier M, Burkhard P. 2006. Statistical analysis of intrahelical ionic interactions in alpha-helices and coiled coils. J Struct Biol 155:116–129. https://doi.org/10.1016/j.jsb.2006.02.019.

30. Croasdale R, Ivins FJ, Muskett F, Daviter T, Scott DJ, Hardy T, Smerdon SJ, Fry AM, Pfuhl M. 2011. An undecided coiled coil: the leucine zipper of Nek2 kinase exhibits atypical conformational exchange dynamics. J Biol Chem 286:27537–27547. https://doi.org/10.1074/jbc.M110.196972.

31. Stewart V, Chen LL. 2010. The S helix mediates signal transmission as a HAMP domain coiled-coil extension in the NarX nitrate sensor from *Escherichia coli* K-12. J Bacteriol 192:734–745. https://doi.org/10.1128/JB.00172-09.

32. Ames P, Hunter S, Parkinson JS. 2016. Evidence for a helix-clutch mechanism of transmembrane signaling in a bacterial chemoreceptor. J Mol Biol 428:3776–3788. https://doi.org/10.1016/j.jmb.2016.03.017.

33. Drobnak I, Gradišar H, Ljubetić A, Merljak E, Jerala R. 2017. Modulation of coiled-coil dimer stability through surface residues while preserving pairing specificity. J Am Chem Soc 139:8229–8236. https://doi.org/10.1021/jacs.7b01690.

34. Schmidt NW, Grigoryan G, DeGrado WF. 2017. The accommodation index measures the perturbation associated with insertions and deletions in coiled-coils. Application to understand signaling in histidine kinases. Protein Sci 26:414–435. https://doi.org/10.1002/pro.3095.

35. Wood CW, Bruning M, Ibarra AA, Bartlett GJ, Thomson AR, Sessions RB, Brady RL, Woolfson DN. 2014. CCBuilder: an interactive web-based tool for building, designing and assessing coiled-coil protein assemblies. Bioinformatics 30:3029–3035. https://doi.org/10.1093/bioinformatics/btu502.

36. Van Der Spoel D, Lindahl E, Hess B, Groenhof G, Mark AE, Berendsen HJ. 2005. GROMACS: fast, flexible, and free. J Comput Chem 26:1701–1718. https://doi.org/10.1002/jcc.20291.

37. Berendsen HJC, Postma JPM, Van Gunsteren WF, Hermans J. 1991. Interaction models for water in relation to protein hydration, p 331–342. *In* Pullman B (ed), Intermolecular forces. Reidel, Dordrecht, The Netherlands.

38. Bussi G, Donadio D, Parrinello M. 2007. Canonical sampling through velocity rescaling. J Chim Phys 126:014101. https://doi.org/10.1063/1.2408420.

39. Hess B. 2008. P-LINCS: a parallel linear constraint solver for molecular simulation. J Chem Theory Comput 4:116–122. https://doi.org/10.1021/ct700200b.

40. Essmann U, Perera L, Berkowitz ML, Darden T, Lee H, Pedersen LG. 1995. A smooth particle mesh Ewald method. J Chem Phys 103:8577–8593. https://doi.org/10.1063/1.470117.

41. Miyamoto S, Kollman PA. 1992. Settle: an analytical version of the SHAKE and RATTLE algorithm for rigid water models. J Comput Chem 13: 952–962. https://doi.org/10.1002/jcc.540130805.

42. Feenstra KA, Hess B, Berendsen HJC. 1999. Improving efficiency of large time scale molecular dynamics simulations of hydrogen-rich systems. J Comput Chem 20:786–798. https://doi.org/10.1002/(SICI)1096-987X (199906)20:8<786::AID-JCC5>3.0.CO;2-B.

43. Daura X, Gademann K, Jaun B, Seebach D, van Gunsteren WF, Mark AE. 1999. Peptide folding: when simulation meets experiment. Angew Chem Int Ed 38:236–240. https://doi.org/10.1002/(SICI)1521-3773(19990115)38:1/2<236::AID-ANIE236>3.0.CO;2-M.

44. Waterhouse AM, Procter JB, Martin DM, Clamp M, Barton GJ (ed). 2009. Jalview Version 2—a multiple sequence alignment editor and analysis workbench. Bioinformatics 25:1189–1191. https://doi.org/10.1093/bioinformatics/btp033.