# Regulation of signaling directionality revealed by 3D snapshots of a kinase:regulator complex in action

Felipe Trajtenberg, Juan A Imelio, Matías R Machado, Nicole Larrieux, Marcelo A Marti, Gonzalo Obal, Ariel E Mechaly, Alejandro Buschiazzo

► To cite this version:

Felipe Trajtenberg, Juan A Imelio, Matías R Machado, Nicole Larrieux, Marcelo A Marti, et al. Regulation of signaling directionality revealed by 3D snapshots of a kinase:regulator complex in action. eLife, 2016, 5, pp.e21422. 10.7554/eLife.21422. pasteur-02554309

HAL Id: pasteur-02554309

https://pasteur.hal.science/pasteur-02554309

Submitted on 25 Apr 2020

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
<figure><img src="image_2.png" /></figure>

# Regulation of signaling directionality revealed by 3D snapshots of a kinase: regulator complex in action

Felipe Trajtenberg¹, Juan A Imelio¹, Matías R Machado², Nicole Larrieux¹, Marcelo A Marti³, Gonzalo Obal⁴†, Ariel E Mechaly¹, Alejandro Buschiazzo¹,⁵*

¹Laboratory of Molecular and Structural Microbiology, Institut Pasteur de Montevideo, Montevideo, Uruguay; ²Biomolecular Simulations, Institut Pasteur de Montevideo, Montevideo, Uruguay; ³Departamento de Química Biológica e IQUIBICEN-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Buenos Aires, Argentina; ⁴Protein Biophysics Unit, Institut Pasteur de Montevideo, Montevideo, Uruguay; ⁵Département de Microbiologie, Institut Pasteur, Paris, France

**Abstract** Two-component systems (TCS) are protein machineries that enable cells to respond to input signals. Histidine kinases (HK) are the sensory component, transferring information toward downstream response regulators (RR). HKs transfer phosphoryl groups to their specific RRs, but also dephosphorylate them, overall ensuring proper signaling. The mechanisms by which HKs discriminate between such disparate directions, are yet unknown. We now disclose crystal structures of the HK:RR complex DesK:DesR from *Bacillus subtilis*, comprising snapshots of the phosphotransfer and the dephosphorylation reactions. The HK dictates the reactional outcome through conformational rearrangements that include the reactive histidine. The phosphotransfer center is asymmetric, poised for dissociative nucleophilic substitution. The structural bases of HK phosphatase/phosphotransferase control are uncovered, and the unexpected discovery of a dissociative reactional center, sheds light on the evolution of TCS phosphotransfer reversibility. Our findings should be applicable to a broad range of signaling systems and instrumental in synthetic TCS rewiring.

DOI: 10.7554/eLife.21422.001

## Introduction

Perception of environmental and intracellular cues is an essential feature of life. Signaling pathways enable cells to regulate genetic and biochemical programs for adaptation and survival. Two-component systems (TCSs) play a particularly important role among protein machineries that cells have evolved to carry out signaling, widely distributed in bacteria, fungi and plants. The simplest TCSs comprise a sensory histidine kinase (HK) and a response regulator (RR) component (*Gao and Stock, 2009*). The signal is transmitted downstream from HK to RR, through an orderly sequence of conformational rearrangements coupled to phosphoryl-transfer reactions. The HK is turned on or off after signal-dependent allosteric rearrangements, which control autophosphorylation on a conserved histidine residue. The phosphoryl group is then transferred from the phosphorylated kinase to a conserved aspartate on the RR receiver domain, activating it to effect a specific output response. When auto-kinase activity is turned off, HKs often act as a phosphatase of their cognate phosphorylated RR (P~RR) (*Goulian, 2010*), contributing to shutting down the pathway. Despite being paradoxical for a protein kinase, the phosphatase activity of HKs is physiologically relevant ensuring robust homeostatic responses (*Schramke et al., 2016*).

<footnote>*For correspondence: alebus@pasteur.edu.uy</footnote>

<footnote>**Present address:** †Molecular Mechanisms of Membrane Transport, Institut Pasteur, Paris, France</footnote>

<footnote>**Competing interests:** The authors declare that no competing interests exist.</footnote>

<footnote>**Funding:** See page 23</footnote>

<footnote>**Received:** 11 September 2016</footnote>

<footnote>**Accepted:** 09 December 2016</footnote>

<footnote>**Published:** 12 December 2016</footnote>

<footnote>**Reviewing editor:** Michael Laub, Massachusetts Institute of Technology, United States</footnote>

<footnote>© Copyright Trajtenberg et al. This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.</footnote>
HKs are homodimeric proteins, comprising an N-terminal sensor domain and a more conserved catalytic core. Trans-membrane HKs display the catalytic region in the cytoplasm and the sensory domain within the lipid bilayer or toward the extracellular/periplasmic space. The catalytic core typically includes a helical Dimerization and Histidine phosphotransfer (DHp) domain, followed by a Catalytic and ATP-binding (CA) domain. HKs are classified according to DHp sequence signatures in HisKA, HisKA_2, HisKA_3, HWE_HK and His_kinase families. CheA-like HKs instead harbor the reactive His within a Histidine Phosphotransfer (HPt) domain that is not involved in dimerization. Overall sequence conservation among HKs can be very low, but the entire catalytic region is remarkably conserved at the structural level (Zschiedrich et al., 2016), suggesting that functional mechanisms are shared. RRs comprise a receiver domain (REC), which may be adjacent to additional effector domains (DNA binding, enzymatic, etc.). REC phosphorylation stabilizes structural rearrangements associated with effector domain activation (Gao and Stock, 2010).

Phosphorelays are examples of more complex TCS pathways, involving additional intermediate phosphotransfer proteins, such as extra REC, HPt and/or modified DHp domains. Phosphorelays allow for more complex signaling circuits and sharper regulation (Goulian, 2010). HKs have evolved to catalyze irreversible phosphotransfer reactions (Porter et al., 2008; Potter et al., 2002) ensuring a unidirectional input→output information flow (Bourret and Stock, 2002). However, evolutionary pressure has also resulted in signaling pathways where reversible steps are simultaneously required, mostly in phosphorelays (Burbulys et al., 1991; Janiak-Spens et al., 2005), where dedicated HPt modules act as phospho-donors and -acceptors. The molecular bases underlying such differential reversibility patterns remain unknown, raising the question of how signaling pathway directionality is ensured, from input signal to adaptive response.

To address this question we have used the TCS DesK-DesR from B. subtilis as a model (de Mendoza, 2014). We had previously shown that the HK DesK, a member of the HisKA_3 family, undergoes important conformational changes to switch between phosphatase- and phosphotransfer-competent states (Albanesi et al., 2009). The cytoplasmic portion of DesK (DesKC) including its entire catalytic region, displays a symmetric and rigid structure in the phosphatase state. DesK activation implies its rearrangement into an auto-kinase competent form with substantially higher flexibility. Upon auto-phosphorylation DesKC adopts a strongly asymmetric conformation (Albanesi et al., 2009) able to transfer the phosphoryl group to DesR, its downstream RR partner. The DHp α-helices include a conserved membrane-proximal coiled-coil motif, upstream of the phosphorylatable His. Such DHp segment is critical for signal transmission, folding into a coiled-coil in the phosphatase state, whereas in the phosphotransferase form this coiled-coil breaks apart (Saita et al., 2015). We now report the crystal structures of the DesKC:DesR complex, trapped in the phosphatase and the phosphotransferase functional states. Extensive conformational rearrangements reveal how the two reaction centers are remodeled. Our data indicate that the relative orientation and distance of the reactive histidine with respect to the receiver aspartate is a molecular determinant controlling signal directionality.

## Results

### Two crystal structures of the DesKC:DesR complex represent snapshots of the dephosphorylation and the phosphotransfer reactions

To grasp the molecular determinants of unidirectional TCS signaling, the crystal structures of DesKC in complex with DesR were determined in two functional states (Figure 1, Table 1). One using the phosphatase-constitutive mutant DesKC<sub>STAB</sub>, and the other with DesKC<sub>H188E</sub> carrying a phosphomimetic substitution (Albanesi et al., 2009), both in complex with the REC domain of DesR. DesKC<sub>STAB</sub> includes amino acid replacements Ser150Ile, Ser153Leu and Arg157Ile in the DHp domain (Figure 1—figure supplement 1A), designed to stabilize the coiled-coil region (Saita et al., 2015).

### The DesKC<sub>STAB</sub>:DesR-REC complex

The asymmetric unit (ASU) of this crystal structure shows one full complex and half of a second one, which is completed through the crystallographic two-fold symmetry operator (Figure 1—figure supplement 1B). Several features indicate that this complex represents a snapshot of the phosphatase reaction, hereafter denominated the 'phosphatase complex'. Symmetric organizations are a hallmark
of HKs in the 'auto-kinase off / phosphatase on' states (Albanesi et al., 2009; Casino et al., 2009; Yamada et al., 2009). The phosphatase complex we have now crystallized indeed displays high symmetry, both in the way the REC domains associate to the HK dimer, as well as between HK monomers. The two independently refined phosphatase complexes show one dimer of DesKC<sub>STAB</sub> bound to two molecules of DesR-REC (Figure 1A and Figure 1—figure supplement 1C), with the latter occupying equivalent positions on either side of the HK core. The 2:2 HK:RR stoichiometry was confirmed in solution by isothermal titration calorimetry (ITC) (Figure 1B, Table 2), revealing an entropy-driven, endothermic association reaction. Size exclusion chromatography-coupled small-angle X ray scattering (SEC-SAXS) further supported the 2:2 stoichiometry (Figure 1C; Tables 3 and 4). There is also a high internal symmetry between both HK monomers, with CA domains rigidly fixed onto the central DHp and a resulting butterfly-like shape of the whole molecule. A second feature consistent with a phosphatase-competent configuration, is that one AMP-PCP (ATP analogue) moiety is bound to each CA domain, far (>27 Å) from the His<sub>188(HK)</sub> phosphorylation sites on the central DHp, precluding auto-kinase activity (subscripts HK and RR highlight the protein to which the indicated residues belong). Finally, the structural rearrangements triggered at high temperatures (Saita et al., 2015), driving DesK to its phosphatase state, have been linked to the stabilization of an N-terminal coiled-coil (Albanesi et al., 2009). The three point-mutations engineered within DesK helix α1, are indeed observed to stabilize a coiled-coil structure toward the N-terminus in the phosphatase complex (Figure 1—figure supplement 1D). Coiled-coil stabilization was devised according to previous structures of free DesKC<sub>H188V</sub> (Albanesi et al., 2009), a mutant version that maintains normal P~DesR-specific phosphatase activities in vitro and in vivo (Albanesi et al., 2004). Free DesKC<sub>H188V</sub> structures superimpose extremely well onto the HK partner of the phosphatase complex (Figure 1—figure supplement 1E). SEC-SAXS data further suggest that the coiled-coil conformation in the phosphatase complex was the most abundant species in solution (Figure 1C and Figure 1—figure supplement 1B; Table 4), involving a large α1:α1 dimerization interface of ~2100 Å² (Figure 1A and Figure 1—figure supplement 1D and E). Regarding the RR component within the phosphatase complex, the three DesR-REC molecules in the ASU display identical 'active-like' structures (Trajtenberg et al., 2014). BeF<sub>3</sub><sup>-</sup> moieties mimicking phosphoryl groups are observed bonded to each phosphorylatable Asp<sub>54(RR)</sub>, and coordinated to Mg<sup>2+</sup> cations, consistent with this crystal form being a snapshot of the phosphatase reaction at its pre-dephosphorylation step.

## The DesKC<sub>H188E</sub>:DesR-REC complex

Wild-type phosphorylated DesKC (wt P~DesKC) in complex with DesR variants resisted crystallization. The crystal structure was eventually solved by using the phosphomimetic DesKC<sub>H188E</sub> mutant in complex with DesR-REC. Three X ray diffraction datasets were obtained from DesKC<sub>H188E</sub>:DesR-REC crystals soaked with different concentrations of MgCl<sub>2</sub> and BeF<sub>3</sub><sup>-</sup> (Table 1). The three structures are similar, showing two DesKC<sub>H188E</sub>:DesR-REC complexes in the ASU (Figure 1D). A number of observations readily distinguish this complex from the phosphatase one, supporting it has captured a ground state of the P~HK→RR phosphoryl-transfer reaction, hereafter referred to as the 'phosphotransferase complex'. In the first place, the phosphotransferase complex is highly asymmetric, a typical feature of phosphorylated and kinase-active forms of HKs (Albanesi et al., 2009; Mechaly et al., 2014; Casino et al., 2014). One DesKC<sub>H188E</sub> dimer was observed bound to one DesR-REC monomer in the crystal. This 2:1 stoichiometry was further supported by ITC (Figure 1E and Table 2) and SEC-SAXS (Figure 1F; Tables 3 and 4) in solution. The symmetry is also broken within the HK dimer, with one of the two CAs interacting with the DHp helix α1 of the other protomer (Figure 1D). The second CA domain is free, visible in the structure due to fortuitous crystal packing contacts. Such arrangement leaves only one RR-binding site available on the HK in the phosphotransferase state, explaining the asymmetric stoichiometry. Secondly, structural superposition of the HK partner extracted from the phosphotransferase complex, with the available structure of free phosphorylated DesK (wt P~DesKC [PDB 3GIG]), reveals closely similar conformations (Figure 1—figure supplement 1F), confirming that the structure of DesKC<sub>H188E</sub> in complex is a reliable model of the wt phosphorylated species. Finally, BeF<sub>3</sub><sup>-</sup> was not observed bonded to Asp<sub>54(RR)</sub> in any of the DesKC<sub>H188E</sub>:DesR-REC complex structures, despite testing elevated soaking concentrations. Yet, the phosphorylation sites of the RR partners are very similar compared to the phosphatase complex, allowing for the Mg<sup>2+</sup> cation to bind alike (the cation was observed bound only in phosphotransferase complex structures
<figure><img src="image_3.png" /><figcaption>**Figure 1.** The phosphatase and the phosphotransferase complexes: crystal structures and stoichiometries in solution. (A) Cartoon representation of the DesKC<sub>STAB</sub>:DesR-REC (phosphatase) complex, along two orthogonal views (left and right subpanels). The two chains within the DesKC<sub>STAB</sub> dimer are depicted in green and yellow, and the two DesR-REC molecules, in magenta and orange. Solvent exposed surfaces are shown in transparent gray. (B) Isothermal titration calorimetry (ITC) of the phosphatase DesKC<sub>H188V</sub>:DesR-REC binding reaction, top panel showing the raw heat flow data; in the bottom, integrated heat exchange as a function of RR(monomer):HK(dimer) molar ratio. ITC was performed with the mutant DesKC<sub>H188V</sub>, equivalently trapped in the phosphatase state (Albanesi et al., 2009), because DesKC<sub>STAB</sub> was insoluble when not co-expressed with DesR-REC. (C) X ray scattering curve for the DesKC<sub>STAB</sub>:DesR-REC complex. Experimental data are plotted as gray dots, with theoretical curves overlaid as colored lines, revealing best fitting for the AB+ABCD mixture (see Table 4). Letters distinguish protomer chains: DesKC<sub>STAB</sub> dimers AB and EE' (primed labels distinguish crystallographically related partners) and DesR-REC monomers C, D, F and F'. (D) Cartoon representation of the DesKC<sub>H188E</sub>:DesR-REC (phosphotransferase) complex, along two orthogonal views (left and right subpanels). Coloring scheme and solvent exposed surface displayed as in (A). (E) ITC of the DesKC<sub>H188E</sub>:DesR-REC binding reaction, details as in (B). (F) SAXS curve for the DesKC<sub>H188E</sub>:DesR-REC complex, details as in (C). Note best fitting with a one dimer DesKC<sub>H188E</sub> (chains AB) to one monomer DesR-REC (chain C) model (see Table 4).</figcaption></figure>

DOI: 10.7554/eLife.21422.002

The following figure supplement is available for figure 1:

**Figure supplement 1.** Phosphatase and phosphotransferase complexes: structural details.

DOI: 10.7554/eLife.21422.003

soaked with high Mg²⁺ concentrations). The absence of BeF₃⁻ in the phosphotransferase state, even though Mg²⁺ was eventually bound, is consistent with phosphorylated HKs not interacting with phosphorylated RR species. The molecular basis for a decreased BeF₃⁻ affinity in the phosphotransferase became clear later by comparing both states' reaction centers (see below). Taken together, the evidence indicates the phosphotransferase complex represents a snapshot of the transfer reaction, prior to the P~His→Asp phosphoryl group migration.

## DesK and DesR interact through a 'slippery' interface

Both phosphatase and phosphotransferase complexes show that the REC domain of DesR interacts with DesK through its α1(RR) helix and the β5α5(RR) loop, with a few additional contacts from the N-terminal portion of helix α5(RR) and the beginning of loop β4α4(RR). On the HK side, the main
<table><caption>**Table 1.** X ray diffraction data collection and refinement statistics.</caption><thead><tr><th></th><th>Phosphatase complex</th><th>Phosphotransferase complex I (low [Mg<sup>2+</sup>])</th><th>Phosphotransferase complex II (high [Mg<sup>2+</sup>])</th><th>Phosphotransferase complex III (high [Mg<sup>2+</sup>] + BeF<sub>3</sub><sup>-</sup>)</th><th>wt P-DesKC</th></tr></thead><tbody><tr><td><strong>Data collection</strong></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Space group</td><td>P3<sub>1</sub>2<sub>1</sub></td><td>P2<sub>1</sub></td><td>P2<sub>1</sub></td><td>P2<sub>1</sub></td><td>P3<sub>1</sub>2<sub>1</sub></td></tr><tr><td>Cell dimensions</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>a, b, c (Å)</td><td>94.3, 94.3, 239.9</td><td>87.8, 114.6, 91.6</td><td>87.9, 115.6, 91.6</td><td>88.2, 116.7, 91.9</td><td>94.4, 94.4, 161.8</td></tr><tr><td>α, β, γ (°)</td><td>90, 90, 120</td><td>90, 116.4, 90</td><td>90, 116.7, 90</td><td>90, 117.1, 90</td><td>90, 90, 120</td></tr><tr><td>Resolution (Å)</td><td>48.35–2.79 (2.94–2.79)*</td><td>66.70–3.2 (3.37–3.2)</td><td>66.78–2.9 (3.06–2.9)</td><td>38.59–3.15 (3.32–3.15)</td><td>36.5–3.16 (3.33–3.16)</td></tr><tr><td>Unique reflections</td><td>31510 (4385)</td><td>26369 (3906)</td><td>36132 (5279)</td><td>28033 (3944)</td><td>14829 (2109)</td></tr><tr><td>R<sub>meas</sub></td><td>0.1 (2.07)</td><td>0.14 (0.98)</td><td>0.13 (1.61)</td><td>0.1 (0.57)</td><td>0.07 (1.18)</td></tr><tr><td>R<sub>pim</sub></td><td>0.03 (0.63)</td><td>0.09 (0.62)</td><td>0.05 (0.65)</td><td>0.05 (0.29)</td><td>0.03 (0.43)</td></tr><tr><td>CC<sub>1/2</sub></td><td>1.00 (0.56)</td><td>0.99 (0.55)</td><td>0.99 (0.67)</td><td>1.00 (0.83)</td><td>1.00 (0.77)</td></tr><tr><td>I / σI</td><td>15.6 (1.5)</td><td>7.7 (1.7)</td><td>12.5 (2.4)</td><td>11.8 (2.7)</td><td>23.8 (1.9)</td></tr><tr><td>Completeness (%)</td><td>99.5 (96.7)</td><td>97.8 (98.8)</td><td>99.1 (99.2)</td><td>98 (95.6)</td><td>99.9 (99.6)</td></tr><tr><td>Redundancy</td><td>10.8 (10.3)</td><td>2.4 (2.4)</td><td>5.8 (6.0)</td><td>3.8 (3.8)</td><td>7.2 (7.2)</td></tr><tr><td><strong>Refinement</strong></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Resolution (Å)</td><td>48.35–2.79</td><td>66.70–3.2</td><td>39.27–2.9</td><td>38.59–3.15</td><td>36.5–3.16</td></tr><tr><td>Number of refls used (N in the free set)</td><td>31454 (1568)</td><td>26231 (1320)</td><td>35374 (1720)</td><td>28015 (1361)</td><td>14787 (791)</td></tr><tr><td>R<sub>work</sub> / R<sub>free</sub></td><td>0.214/0.249</td><td>0.187/0.24</td><td>0.195/0.239</td><td>0.19/0.231</td><td>0.257/0.295</td></tr><tr><td>Number of atoms</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Protein</td><td>7972</td><td>8869</td><td>8702</td><td>8709</td><td>3151</td></tr><tr><td>Ligands + ions</td><td>93 (AMP-PCP)/6 (Mg<sup>2+</sup>)/12 (MES)/12 (BeF<sub>3</sub>)/5 (SO<sub>4</sub><sup>2-</sup>)</td><td>124 (AMP-PCP)/4 (Mg<sup>2+</sup>)/2 (K<sup>+</sup>)</td><td>124 (AMP-PCP)/5 (Mg<sup>2+</sup>)/2 (K<sup>+</sup>)</td><td>124 (AMP-PCP)/6 (Mg<sup>2+</sup>)/2 (K<sup>+</sup>)</td><td>62 (AMP-PCP)/2 (Mg<sup>2+</sup>)/6 (glycerol)</td></tr><tr><td>Water</td><td>12</td><td>15</td><td>9</td><td>18</td><td>4</td></tr><tr><td>B-factors (Å<sup>2</sup>)</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Wilson plot</td><td>98.3</td><td>91.4</td><td>95.9</td><td>84.1</td><td>128.9</td></tr><tr><td>Mean (overall)</td><td>118.7 <sup>‡</sup></td><td>93.1 <sup>‡</sup></td><td>102.5 <sup>‡</sup></td><td>89.6 <sup>‡</sup></td><td>123.8 <sup>‡</sup></td></tr><tr><td>R.m.s. deviations</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Bond lengths (Å)</td><td>0.01</td><td>0.01</td><td>0.01</td><td>0.01</td><td>0.01</td></tr><tr><td>Bond angles (°)</td><td>1.2</td><td>1.3</td><td>1.25</td><td>1.25</td><td>1.3</td></tr><tr><td>Number of residues in Ramachandran plot <sup>§</sup> (favored / outliers)</td><td>1025/1</td><td>1099/3</td><td>1074/2</td><td>1073/5</td><td>382/4</td></tr><tr><td>PDB ID</td><td>5IUN</td><td>5IUJ</td><td>5IUK</td><td>5IUL</td><td>5IUM</td></tr></tbody></table>

element engaged in the interface is helix α₁(HK) on the DHp domain from one DesK protomer, and minor contacts with helix α₂(HK) from the second protomer. This interface broadly resembles the one found in other TCS complexes (Casino et al., 2009; Willett et al., 2015; Yamada et al., 2009; Zapf et al., 2000) (Figure 2A).

The DHp:REC interface engages van der Waals contacts almost exclusively, including several hydrophobic residues (Figure 2B), which clustered among the few HK:RR covariant residue pairs (Figure 2—figure supplement 1A) found to be conserved in other HK and RR families

<footnote>*Values in parentheses correspond to the highest-resolution shell.</footnote>

<footnote>‡ Including TLS contribution.</footnote>

<footnote>§ Calculated with Molprobity (Chen et al., 2010)</footnote>

<footnote>DOI: 10.7554/eLife.21422.004</footnote>
<table><caption>**Table 2.** Isothermal titration calorimetry parameters. For details on the two different titration procedures, see Materials and methods.</caption><thead><tr><th>HK:RR complexes</th><th>K<sub>a</sub> (M<sup>-1</sup>)</th><th>ΔG (x10<sup>3</sup> kcal.mol<sup>-1</sup>)</th><th>ΔH (x10<sup>4</sup> kcal.mol<sup>-1</sup>)</th><th>TΔS (x10<sup>4</sup> kcal.mol<sup>-1</sup>)</th></tr></thead><tbody><tr><td rowspan="2">DesKC<sub>H188V</sub>:DesR-REC</td><td>7.7×10<sup>5</sup></td><td>-7.8</td><td>1.5</td><td>2.3</td></tr><tr><td>9.5×10<sup>4</sup></td><td>-6.5</td><td>1.4</td><td>2.1</td></tr><tr><td rowspan="2">DesKC<sub>H188E</sub>:DesR-REC</td><td>1.7×10<sup>6</sup></td><td>-8.2</td><td>2.2</td><td>3.0</td></tr><tr><td>3.0×10<sup>4</sup></td><td>-5.9</td><td>2.1</td><td>2.7</td></tr></tbody></table>

DOI: 10.7554/eLife.21422.005

(Skerker et al., 2008). Leu<sub>200(HK)</sub> plays a central role by inserting its side chain into a hydrophobic pocket within DesR, delimited by the side chains of Leu<sub>13(RR)</sub>, Ala<sub>16(RR)</sub>, Leu<sub>17(RR)</sub> and Leu<sub>20(RR)</sub> on one side, and the main chain atoms of the loop spanning residues 102–106<sub>(RR)</sub>, on the other (Figure 2—figure supplement 1B). Ser<sub>196(HK)</sub> and Gln<sub>193(HK)</sub>, located on the same face of α1<sub>(HK)</sub> as Leu<sub>200(HK)</sub>, establish additional van der Waals contacts within a groove on the RR’s surface, extending the Leu<sub>200(HK)</sub>-lodging pocket all the way to the reaction center, which is rather open (Figure 2—figure supplement 1B). Very few polar contacts are observed surrounding the hydrophobic core, on both ends of the interface patch: Asp<sub>203(HK)</sub>/Arg<sub>206(HK)</sub> with Ser<sub>106(RR)</sub>/Glu<sub>107(RR)</sub>, and a salt bridge between Arg<sub>235(HK)</sub> and Asp<sub>103(RR)</sub> on the other end. The buried surface is ~1400 Å<sup>2</sup> in the phosphatase state, the DHp:REC interface being the major contributor (~1000 Å<sup>2</sup>) (Figure 2C and Figure 2—figure supplement 1C), although the CA domain adds a buried area of >400 Å<sup>2</sup>. On the other hand, the phosphotransferase complex shows a slightly smaller interface (~1000 Å<sup>2</sup>), dominated by the DHp: REC interaction (~900 Å<sup>2</sup>) (Figure 2D and Figure 2—figure supplement 1C).

Superposing the Cα atoms of the DHp region involved in DesR-binding (DesK residues 190–235) among the nine independently refined DesKC:DesR interfaces, reveals large variability in the relative positioning of the REC domains, with rmsd as high as 2.4 Å for the RECs’ 129 Cα atoms (Figure 2E, Table 5). These significant movements engage corresponding shifts in the positions of interacting residues, albeit preserving the contacts, reminiscent of a ‘slippery’ interface. The slippery nature of the interface is consistent with its planarity, a low number of mainly hydrophobic contacts, leading to low shape complementarity and unfilled cavities within the RR partner (Figure 2F).

## A single response regulator conformation is selected by the two functional states of the kinase

Comparing the phosphatase vs phosphotransferase complexes, the HK partner reveals substantial rearrangements. Concerning the DHp domain, a large rotational shift is observed on the membrane-proximal side of the phosphorylation site (Figure 3—figure supplement 1A), in contrast to the distal RR-binding site, which remains essentially invariant. The relative orientation of the CA domains is also dramatically changed (Figure 3—figure supplement 1B). On the other hand, the RR REC domain adopts a unique conformation in both complexes, with an average rmsd of 0.3 Å aligning all Cα atoms, and no significant contrast clustering phosphatase vs phosphotransferase structures (Table 5). Interestingly, this conformation adopted by the REC domain in complex with the kinase, combines structural features of both active and inactive states as seen in the free forms of RRs (Gao and Stock, 2009; Trajtenberg et al., 2014) (Figure 3). Namely, the β4α4 loop and β5 strand

<table><caption>**Table 3.** Small angle X ray scattering and derived molecular size parameters. Data derived from experiments performed by SEC-coupled SAXS. Figures for the phosphatase-trapped species DesKC<sub>H188V</sub>, are to be compared with those corresponding to HK:RR species. Full details in Materials and methods.</caption><thead><tr><th>Protein species</th><th>I<sub>(0)</sub></th><th>I<sub>(0)</sub> real</th><th>Rg (Guinier)</th><th>Rg (real)</th><th>Vc</th><th>MM<sub>Vc</sub> (Da)</th><th>MM<sub>seq</sub> (Da)</th><th>Dmax (Å)</th></tr></thead><tbody><tr><td>DesKC<sub>H188V</sub></td><td>4.47×10<sup>-2</sup></td><td>4.26×10<sup>-2</sup></td><td>31.25</td><td>31.21</td><td>442.9</td><td>46400</td><td>44780</td><td>100</td></tr><tr><td>DesKC<sub>STAB</sub>:DesR-REC</td><td>1.16×10<sup>-1</sup></td><td>1.11×10<sup>-1</sup></td><td>31.48</td><td>31.42</td><td>494.8</td><td>58400</td><td>77110</td><td>103</td></tr><tr><td>DesKC<sub>H188E</sub>:DesR-REC</td><td>1.05×10<sup>-1</sup></td><td>9.76×10<sup>-2</sup></td><td>32.84</td><td>32.18</td><td>486.2</td><td>51800</td><td>61790</td><td>103</td></tr></tbody></table>

DOI: 10.7554/eLife.21422.006
<table><caption>**Table 4.** Fitting figures (χ2) comparing alternative theoretical SAXS curves to the experimentally collected scattering data.</caption><thead><tr><th>Experimental Data</th><th>Structural models</th><th>χ<sup>2</sup></th></tr></thead><tbody><tr><td></td><td>AB+ABCD mixture</td><td>6.1</td></tr><tr><td></td><td>ABCD (folded coiled-coil)</td><td>8.9</td></tr><tr><td></td><td>EE'FF' (disrupted coiled-coil)</td><td>11.6</td></tr><tr><td>DesKC<sub>STAB</sub>:DesR-REC</td><td>ABC</td><td>8</td></tr><tr><td>DesKC<sub>H188E</sub>:DesR-REC</td><td>AB</td><td>24.9</td></tr></tbody></table>

DOI: 10.7554/eLife.21422.007

emulate the inactive RR, while the rest of the protein resembles the free phosphorylated species. The β4α4 loop is closed in phosphorylated RRs (e.g. PDB 4LE0), coupled to a fully wound first turn of helix α4 (Pro<sub>85(RR)</sub> displays a dihedral ψ angle near -40°, allowing for α-helicity). The REC molecules in both phosphatase and phosphotransferase complexes, show instead this β4α4 loop open (with Pro<sub>85(RR)</sub> ψ angles of respectively 126° and 130°, and the first turn of α4 unwound, as observed in the inactive state of the Mg<sup>2+</sup>-free protein [PDB 4LE1]). Thus acting as a 'phosphate lid', the β4α4 loop opens and closes to control phosphoryl-transfer chemistry on Asp<sub>54(RR)</sub> (Figure 3).

Regarding the positions of the phosphorylatable Asp, as well those of the β1α1 loop and contiguous helix α1, HK-binding triggers RR 'activation', as we had previously predicted based on functional readouts and in silico modeling (Trajtenberg et al., 2014). The shift of loop β1α1 locates key acidic residues in place to coordinate the essential Mg<sup>2+</sup> cation within the RR's reaction center.

Taken together, the RR adopts a single intermediate pre-activated conformation, instrumental for both dephosphorylation and phosphotransfer reactions, due to increased exposure of the phosphorylation site and stabilized configuration of the pocket poised for reaction.

## The phosphatase to phosphotransferase transition: reconfiguring the reaction center for catalysis

The crystal structure of wt P-DesKC was determined (Table 1), improving previous resolution limits (Albanesi et al., 2009). The reactive His<sub>188(HK)</sub> was observed to be phosphorylated on its Nε2 atom (Figure 4—figure supplement 1A), and superposition of this model onto the HK partner of the phosphotransferase complex (PDB 5IUK) revealed very similar structures, with 1.2 Å rmsd aligning all 1994 atoms that comprise the DHp domain and the fixed CA (Figure 4—figure supplement 1B). Substituting the DesKC<sub>H188E</sub> component in the crystal complex by the superposed wt P-DesKC, revealed surprisingly good geometry between P-DesKC and DesR-REC (Figure 4—figure supplement 1C). Energy minimization easily fixed minor interatomic bumping, rendering a model with P-His<sub>188(HK)</sub> well oriented toward Asp<sub>54(RR)</sub>, poised for phosphotransfer (Figure 4A). The P-His<sub>188(HK)</sub>Nε2, phosphorus and Asp<sub>54(RR)</sub>Oδ1 atoms describe a~180° angle, with the phosphoryl group at H-bonding distance to the side chain oxygens of Thr<sub>80(RR)</sub> and Thr<sub>81(RR)</sub>. The model proved to be stable in molecular dynamics simulations (Figure 4B), overall making chemical and biological sense. We shall thus use the minimized wt P-DesKC:DesR-REC model to perform further structural analyses of the phosphotransferase complex reaction center.

The phosphatase→phosphotransferase switch engages major rearrangements. An ~80° rotation of helix α1 is coupled to a cogwheel-like shift of helix α2, at the level of the phosphorylation site within the DHp domain (Albanesi et al., 2009). This rotation includes the phosphorylationable His<sub>188(HK)</sub> (Figure 4C), which had not been previously observed. Due to this shift, the side chain of His<sub>188(HK)</sub> ends up buried in the phosphatase state, establishing a H-bond with the same residue of the other protomer (Figure 4C and Figure 4—figure supplement 1D), displaced away from the reaction center as compared to the phosphotransferase (Figure 4D). Asp<sub>189(HK)</sub> right next to the phosphorylation site, follows a similarly dramatic rearrangement, well positioned only in the phosphatase complex to
<figure><img src="image_4.png" /><figcaption>**Figure 2.** The HK:RR interface and structural variability. (A) Structural superposition of the DesK:DesR phosphatase complex (this study, PDB 5IUN) with HK853:RR468 from *Thermotoga maritima* (PDB 3DGE). Only selected structural elements from both partners are shown as cartoons for clarity. The DesK:DesR phosphotransferase complex is overall similar (not shown here). DesK and DesR are illustrated with the same color code as in **Figure 1**. Superposed *T. maritima* RR468 helix α1 is shown in gray, and HK853 protomers in cyan and blue. Primed labels distinguish DHp helices from each HK protomer. (B) Close-up of the phosphatase complex DesK:DesR interface showing only a few selected interactions for clarity (see text for detailed description). Coloring scheme as in (A). Note Leu200(HK) inserted into a hydrophobic pocket of DesR, and several other hydrophobic residues completing the interface (not labeled, shown as sticks). Among the polar contacts surrounding the hydrophobic core, Arg235 (HK):Asp103(RR) and Gln193(HK):BeF3⁻-modified Asp54(RR), are highlighted. (C) Solvent accessible surface of DesKCSTAB indicating the interaction footprint (in magenta) with DesR in the phosphatase complex. DesK domains are highlighted in green (DHp) and blue (CAs); these domains participate differently in DesR interaction, depicted in light magenta (DHp) and dark magenta (CA). Note the ATP analogue (AMP-PCP, in sticks colored by atom) bound to the CA domains, visible in this view on the rightmost one. (D) Same as (C), but for the phosphotransferase complex, illustrating the DesR-binding surface of DesKC<sub>H188E</sub>. (E) The variable position of the REC<sub>(RR)</sub> domain with respect to the DHp<sub>(HK)</sub> is illustrated, superimposing the structurally invariant region of the DHp (including the nine Figure 2 continued on next page</figcaption></figure>
Figure 2 continued

independently refined DesK:DesR complexes). The coloring scheme is the same as in (A), with light and dark colors distinguishing phosphotransferase and phosphatase complexes, respectively. (F) The solid volume corresponds to DesR-REC, which is shown sliced to highlight its outline, revealing unfilled cavities at the protein: protein interface. The relative position of interfacing amino acids is roughly indicated with residue labels. DesK DHp helices are shown in cartoon representation, with its superposed molecular surface on top. Two key residues on the HK partner are shown as sticks.

DOI: 10.7554/eLife.21422.008

The following figure supplements are available for figure 2:

**Figure supplement 1.** Structural details of the DesK:DesR interface.

DOI: 10.7554/eLife.21422.009

**Figure supplement 2.** Molecular dynamics (MD) simulation of the phosphotransferase complex.

DOI: 10.7554/eLife.21422.010

make a salt bridge with Arg84(RR). Neighboring Phe82(RR) thus inserts its aromatic side chain into a hydrophobic pocket formed between DHp α1 of one HK protomer and α2' of the other (Figure 4D). That the reactive histidine moves away from the catalytic site in the phosphatase state is consistent with its marginal role in P~RR~ dephosphorylation catalysis in HisKA_3 HKs. Indeed, the phosphorylable His was proved nonessential for NarX-mediated P~NarL~ dephosphorylation (Huynh et al., 2010). His has been implicated in EnvZ-mediated phosphatase catalysis (Zhu et al., 2000) in which case its imidazole would act as a base assisting the attack of water on the P~Naspartyl~ group (Bhate et al., 2015). But contradictory results cast doubts on such a direct role (Hsing and Silhavy, 1997; Skarphol et al., 1997). Instead, a Gln in HisKA_3 (Huynh et al., 2010), corresponding to Gln193(HK) in DesK, or a Thr/Asn in HisKA HKs (Willett and Kirby, 2012), equivalently positioned one helical turn C-terminal to the phosphorylable His, play a key role in the phosphatase reaction, positioning a catalytic water molecule. Such water or hydroxide molecule is well located to perform the nucleophilic attack on the phosphoryl group, as observed in the P~CheY3~ phosphatase CheX (Pazy et al., 2010), a shared geometry among phosphatases and HKs as we now observe in the DesKC:DesR-REC phosphatase complex (Figure 4—figure supplement 2).

<table><caption>**Table 5.** Square matrix of all possible pair-wise structural superpositions. Superpositions were calculated using the nine independently refined DesKC:DesR-REC pairs, taken from the 1.5 complexes in the phosphatase asymmetric unit (named STAB1, STAB2 and STAB3) and the 2 complexes present in the ASU of each of the three different phosphotransferase structures (named as E188[1-6]). Below the main diagonal the structurally invariant region of DesK DHp domains (residues 190-234(HK)) were superposed, whereas above the diagonal, the DesR REC domains (residues 1-129(RR)) were superposed instead. Resulting root mean squared deviations (in Å) were calculated, for both matrix halves, between all Cα atoms of the REC domains (residues 1-129(RR)) after superposition, as indicated in each individual matrix cell. Font colors highlight the two different alignment procedures performed. The significant differences between corresponding blue- and red-colored values, indicate that the REC domains do not change within, but rather 'slip' with respect to the kinase DHp domain.</caption><thead><tr><th></th><th>STAB1</th><th>STAB2</th><th>STAB3</th><th>E<sub>188</sub>1</th><th>E<sub>188</sub>2</th><th>E<sub>188</sub>3</th><th>E<sub>188</sub>4</th><th>E<sub>188</sub>5</th><th>E<sub>188</sub>6</th></tr></thead><tbody><tr><th>STAB1</th><td></td><td>0.229</td><td>0.158</td><td>0.414</td><td>0.349</td><td>0.442</td><td>0.361</td><td>0.344</td><td>0.343</td></tr><tr><th>STAB2</th><td>0.826</td><td></td><td>0.188</td><td>0.496</td><td>0.459</td><td>0.512</td><td>0.334</td><td>0.343</td><td>0.333</td></tr><tr><th>STAB3</th><td>1.362</td><td>1.917</td><td></td><td>0.393</td><td>0.343</td><td>0.414</td><td>0.292</td><td>0.287</td><td>0.277</td></tr><tr><th>E<sub>188</sub>1</th><td>1.168</td><td>1.133</td><td>1.682</td><td></td><td>0.275</td><td>0.146</td><td>0.305</td><td>0.301</td><td>0.301</td></tr><tr><th>E<sub>188</sub>2</th><td>0.909</td><td>1.140</td><td>1.224</td><td>0.584</td><td></td><td>0.301</td><td>0.351</td><td>0.309</td><td>0.341</td></tr><tr><th>E<sub>188</sub>3</th><td>1.085</td><td>1.110</td><td>1.530</td><td>0.233</td><td>0.459</td><td></td><td>0.346</td><td>0.315</td><td>0.313</td></tr><tr><th>E<sub>188</sub>4</th><td>1.548</td><td>1.140</td><td>2.453</td><td>0.893</td><td>1.316</td><td>1.040</td><td></td><td>0.175</td><td>0.117</td></tr><tr><th>E<sub>188</sub>5</th><td>1.295</td><td>0.978</td><td>2.096</td><td>0.593</td><td>0.969</td><td>0.716</td><td>0.449</td><td></td><td>0.140</td></tr><tr><th>E<sub>188</sub>6</th><td>1.424</td><td>1.058</td><td>2.300</td><td>0.770</td><td>1.182</td><td>0.905</td><td>0.193</td><td>0.306</td><td></td></tr></tbody></table>

DOI: 10.7554/eLife.21422.011
<figure><img src="image_5.png" /></figure>

**Figure 3.** A shared RR REC domain conformation in the phosphatase and phosphotransferase complexes, has active- and inactive-like features. The phosphotransferase complex (this study, PDB 5IUK) is shown in cartoon representation, following the same coloring scheme as in **Figure 2**. REC domains from reported free DesR species are shown superposed, in gray the inactive Mg²⁺-free RR (PDB 4LE1) and the active species in cyan (PDB 4LE0). Note that the complexed form of the regulator (in magenta) is mainly similar to the active species, including in the β1α1 loop and α1 helix (that play key roles in the structuring of the RR's active site), even though the RR is not phosphorylated on the reactive Asp₅₄(RR). Dashed outlines highlight the regions that most resemble each functional state. Note that the inactive-like region essentially involves the end of β4 and the β4α4 loop, which closes in as a phosphate lid in P~DesR. Phe₈₂(RR) (shown in transparent sticks for clarity) interacts with the HK, playing a role in phosphate lid opening in both phosphatase and phosphotransferase complexes. Tyr₉₉(RR) on strand β5 (the strand is not visible below α4), typically forming a H-bond with the main chain O of Phe₈₂(RR) in activated RRs, is observed farther away in both DesK:DesR complexes, yet another feature of inactive states. Inset: a close-up of the N-terminal tip of helix α4 is shown. The dihedral ψ angle of Pro₈₅(RR) is linked to the active/inactive shift of the β4α4 loop.

DOI: 10.7554/eLife.21422.012

The following figure supplement is available for figure 3:

**Figure supplement 1.** Major structural differences of the HK partner, comparing the phosphatase and phosphotransferase complexes.

DOI: 10.7554/eLife.21422.013

The rotational phosphatase→phosphotransferase rearrangement disrupts the Asp₁₈₉(HK):Arg₈₄(RR) ionic bond, and coupled to the α₁(HK) and α₂'(HK) movements, the position of Phe₈₂(RR) shifts, propagating along the rest of the phosphate lid β4α4 loop (**Figure 4D**). This shift accounts for the capacity of Thr₈₀(RR) to be at H-bonding distance from the P~His₁₈₈(HK) phosphoryl moiety, accompanying its transfer, and maintaining the H-bond to the phosphate on P~Asp₅₄(RR) (**Figure 4D**). The side chain of Phe₈₂(RR) is also contacting P~His₁₈₈(HK), which could increase the pKa of its imidazole Nδ1 atom due to cation-π interactions (**Loewenthal et al., 1992**). Indeed, P~DesKC appeared to catalyze slower phosphoryl-transfer to DesR-REC₁₈₂A compared to wt DesR-REC (**Figure 4E**), consistent with a
poorer Nδ1 protonation, and hence a stronger P-N phosphoramidate bond (Attwood et al., 2007). Systematic mutagenesis at RR position 82 and functional assays (catalytic activities and HK:RR affinities) will be needed to reach a definitive conclusion.

Thus, the phosphatase and phosphotransferase catalytic activities depend upon modifications in the position of a subset of residues in both partners. Both complexes share a common RR intermediate conformation, and the reaction outcome is dictated by the HK switch of helices α1 and α2.

## Dissociative phosphotransfer is linked to unidirectional signaling

DesK:DesR-catalyzed phosphotransfer occurred essentially in the forward, P~DesK→DesR direction (Figure 5A), extending earlier suggestions (Albanesi et al., 2004) into quantitative figures, and indicating that back-transfer is minimal. How is immediate P~RR dephosphorylation prevented once the phosphate has been transferred? The side chain of Thr<sub>80(RR)</sub> at the end of β4 strand, appears to drive the whole phosphate lid to a closed configuration, escorting the phosphoryl group during P~His→Asp migration (Figures 3 and 4D). To test whether a closed phosphate lid indeed plays a role in minimizing P~RR dephosphorylation, as a means to ensuring phosphotransfer unidirectionality, Phe<sub>82(RR)</sub> was substituted by alanine. Phe<sub>82(RR)</sub> covers the P~Asp<sub>54(RR)</sub> in active free DesR, while exposing it to bulk solvent through phosphate lid-opening when bound to DesK (Figure 4—figure supplement 2). The half-life of P~DesR<sub>F82A</sub> was significantly reduced when compared to wt P~DesR (Figure 5B), confirming that the phosphate lid is a key structural element to avoid otherwise futile P~RR dephosphorylation after phosphotransfer.

Then how is P~RR→HK back-transfer inhibited? The phosphate lid could also play a role. Yet, comparing wt P~DesR with point-mutants P~DesR<sub>F82A</sub> or P~DesR<sub>RR84A</sub>, no change was detected in their marginal abilities to transfer the phosphoryl group back to DesK (data not shown). There is however a parameter that was found to discriminate irreversible phosphotransfer systems from others that reveal appreciable reversibility: the reactional His-Asp distance. The distance between P~His<sub>188(HK)</sub>Nε2 and Asp<sub>54(RR)</sub>Oδ1 in the phosphotransferase complex is ~7.6 Å (Figures 4A and 5C), strongly suggesting that the phosphotransfer reaction occurs through a loose (or predominantly dissociative) nucleophilic substitution mechanism. This is in contrast to what has been observed in reversible complexes (Figure 5C), which efficiently catalyze phosphotransfer also in the opposite, P~Asp→His direction. The latter display significantly shorter Nε2-to-Oδ1 distances, consistent with an associative nucleophilic mechanism (Table 6). Examples include several phosphorelay protein pairs (Bauer et al., 2013; Zapf et al., 2000; Zhao et al., 2008) as well as Escherichia coli CheY variants in complex with imidazole (Page et al., 2015), or yet two autophosphorylating HKs where the phosphoryl group is poised to migrate from ATP to the reactive histidine (Casino et al., 2014; Mechaly et al., 2014). In contrast, HK:RR TCSs with unidirectional P~His→Asp transfer, such as B. subtilis DesK:DesR (this study), Thermotoga maritima HK853:RR469 (Casino et al., 2009) and Rhodobacter sphaeroides CheA3:CheY6 (Bell et al., 2010), consistently show >6.6 Å HisNε2-AspOδ1 distances and correlated large reaction coordinate distances when these could be measured (Table 6). It will be interesting to test whether the phosphorelay pair CbpT:CtrA from Brucella abortus (Willett et al., 2015) displays reversible transfer, as anticipated on the basis of the reported 5.9 Å HisNε2-AspOδ1 distance.

In a loose reaction center like that of DesK:DesR, a largely dissociated meta-phosphate intermediate would be stabilized migrating toward the positive charges of the Mg<sup>2+</sup> cation and the conserved Lys<sub>102(RR)</sub> in the RR partner. The position of these positive charges, next to the RR phosphorylation site, introduces an intrinsic asymmetry within the reaction center, which must translate into unequal likelihood for phosphoryl-transfer directions. As the Mg<sup>2+</sup> position becomes more asymmetric, the P~Asp→His back-transfer would be less favorable, explaining unidirectionality. The ratio between the Mg<sup>2+</sup>-Asp<sub>(Oδ1)</sub> and the Mg<sup>2+</sup>-His<sub>(Nε2)</sub> distances, can be used to quantitate such asymmetry within reaction centers. This ratio was calculated for different TCS complexes (Table 6), indeed clustering those complexes that catalyze reversible P~Asp→His/P~His→Asp reactions, apart from the ones that carry out P~His→Asp unidirectional transfer (Figure 5D). Signaling pathways appear to have evolved associative or dissociative mechanisms, corresponding respectively to more reversible or irreversible reactions, such that they are fit to drive a defined directionality in the flow of information.
<figure><img src="image_6.png" /><figcaption>**Figure 4.** At the reaction center: the phosphatase to phosphotransferase transition. (A) Cartoon model of the phosphotransferase reaction center prior to DesR phosphorylation (see model construction details the Materials and methods section). DesK protomers (orange and blue) and DesR (cyan, transparent for clarity) are displayed with selected key residues in sticks colored by atom. Mg²⁺ (green sphere) is already in place, coordinated by the two shown Asp residues and water molecules (not included). Phosphoryl moiety interactions with DesR, and the reactive Asp-His distance are indicated. (B) Evolution of the atomic coordinates of the phosphotransferase complex along molecular dynamics (MD) calculations. Selected HK or RR domains were structurally aligned (marked with orange subscripts on each curve's label on top), to thereafter compare the evolving MD model with the initial experimental structure (calculating rmsds of chosen domains as marked in black fonts on each curve's label on top). Resulting rmsds for all Cα atoms of chosen domains are plotted (colored curves) as a function of time. Note that the time lapse is enough to detect large CA mobility (gray curve) whereas the DHp:REC complex remains attached and stable (pink curve). (C) Cartoon illustration of the HK DHp domains of the phosphatase complex (P), with its two HK protomers in green and yellow, superposed onto the phosphotransferase (PT) in orange and blue. Residue His₁₈₈(HK) (in sticks) reveals the rotational rearrangement between both states. (D) Similar phosphatase vs phosphotransferase illustration as in (C), along a different view. The RR partner is now shown (magenta for the phosphatase complex [P], cyan for the phosphotransferase [PT]). The Arg₈₄(RR):Asp₁₈₉(HK) salt bridge is disrupted in the PT complex due to the DHp rotational shift (red arrow). Note the shift of the RR β4α4 loop, including Phe₈₂(RR), propagating toward Thr₈₀(RR). The latter is positioned at H-bonding distance to the phosphoryl group either on the P-His or on the P-Asp residue (the BeF₃⁻ moiety in the P complex is transparent), with a shift of 1.5 Å of its side chain O atom. (E) Phosphotransfer kinetics comparing wt DesR-REC (top panel) with phosphate lid mutant DesR<sub>F82A</sub>-REC (bottom panel), revealed by Phos-tag SDS-PAGE.</figcaption></figure>

DOI: 10.7554/eLife.21422.014

The following figure supplements are available for figure 4:

**Figure supplement 1.** Electron density maps and structural similarity of wt P~DesKC compared to DesKC<sub>H188E</sub> in the phosphotransferase complex.

DOI: 10.7554/eLife.21422.015

**Figure supplement 2.** Structural details of the reaction centers.

DOI: 10.7554/eLife.21422.016

**Figure supplement 3.** Phosphate lid opening for RR dephosphorylation.

DOI: 10.7554/eLife.21422.017
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Dissociative phosphotransfer and asymmetric position of the divalent cation. (A) Phosphoryl-transfer reactions were analyzed in vitro using Phos-tag SDS-PAGE. Band intensities corresponding to unphosphorylated and phosphorylated species were quantified by densitometry. Forward- and back-transfer directions were compared, with bars representing the amount of phosphotransfer in each direction, plotted as the percentage of phosphorylated DesKC with respect to its total initial amount (100%). Forward phosphotransfer (black bar) was measured by incubating pre-phosphorylated P~DesKC and DesR-REC. Back-transfer (gray bar) was assayed by incubating pre-phosphorylated P~DesR and DesKC. The standard deviation for triplicate replicas are shown as error bars. (B) Stabilization of the phosphorylated DesR species by the phosphorylation lid. The intrinsic dephosphorylation velocities of P~DesR-REC and P~DesR<sub>F82A</sub>-REC were compared using Phos-tag SDS-PAGE. One of three similar replicas is shown. (C) Structural superposition of the DesK:DesR phosphotransferase complex (colored in green:magenta) onto the *B. subtilis* Spo0B:Spo0F complex (PDB 2FTK, depicted in cyan:white), aligning the DesR and Spo0F REC domains. Selected key residues are labeled and numbered according to DesR's sequence. Note the strictly conserved position of the divalent cation (green sphere). The relative shift in the position of the phosphoryatable histidine (His<sub>188</sub> in DesK corresponds to His<sub>30</sub> in Spo0B), is highlighted by depicting the His<sub>(Ne2)</sub>-Asp<sub>(O51)</sub> distances, not including the phosphoryl group for clarity. Note that besides the His shift, the overall DHP:REC binding interface is otherwise preserved. (D) Plot of His<sub>(Ne2)</sub>-Asp<sub>(O51)</sub> distance vs degree of asymmetry in the position of the Mg<sup>2+</sup> cation, comparing independently refined TCS complexes (color-labeled). The inset shows the two interatomic distances used to calculate the ratio plotted in the x-axis. Note the clustering in two groups, which discriminate complexes catalyzing reversible (yellow-to-brown range of points) vs irreversible (blue range of points) phosphotransfer reactions.</figcaption></figure>

DOI: 10.7554/eLife.21422.018

Discussion

Signaling requires specificity, to ensure that a given stimulus is linked to a defined adaptive response while minimizing cross-talk; efficiency, to avoid wasting cellular energy through futile cycles; and directionality, to guarantee the information is directed from the stimulus to the output response. Two component systems constitute a particularly interesting biological model to understand signaling, considering that so many histidine kinases catalyze both kinase and phosphatase reactions according to their functional status, implying exquisite regulation mechanisms at play.

By solving the crystal structures of a bona fide HK in complex with its cognate RR in two distinct functional states, such molecular mechanisms of TCS unidirectional signal transmission are being uncovered. These structures correspond to snapshots of the phosphotransferase reaction just prior of actual phosphoryl migration, and of the phosphatase reaction in a pre-dephosphorylation complex.
<table><caption>**Table 6.** Summary of reactive His(HK)-to-Asp(RR) distances in reported TCS protein complexes.</caption><thead><tr><th></th><th>Group (based on distance discriminants*)</th><th>Complex partners (N° complexes per ASU†)</th><th>Distance his(Nε2)-Asp(Oδ1) (Å)‡</th><th>Reaction coordinate distance (Å)§</th><th>Distance ratio [Mg<sup>2+</sup>-Asp(Oδ1)/Mg<sup>2+</sup>-His(Nε2)]</th><th>Pdb (ID)</th><th>Phosphotransfer reaction direction</th><th>Remarks</th></tr></thead><tbody><tr><td rowspan="4">Phosphorelays</td><td rowspan="4">I</td><td>Spo0F: Spo0B¶ (4)</td><td>5.7; 5.3; 5.1; 5.1</td><td>4.1; 3.8; 3.6; 3.5</td><td>0.82; 0.74; 0.65; 0.67</td><td>2FTK</td><td>P-Asp(Spo0F)→His(Spo0B)</td><td>Spo0B bears a DHp domain</td></tr><tr><td>SLN1: YPD1** (1)</td><td>4.9</td><td>3.2</td><td>0.67</td><td>2R25</td><td>P-Asp(SLN1)→His(YPD1)</td><td>YPD1 bears an HPt domain</td></tr><tr><td>AHK5: AHP1†† (1)</td><td>6.4</td><td>ND[8]</td><td>0.72</td><td>4EUK</td><td>P-Asp(AHK5-REC)→His(AHP1)</td><td>AHP1 bears an HPt domain</td></tr><tr><td>ChpT: CtrA§§ (2)</td><td>5.9¶¶</td><td>ND††</td><td>0.64***</td><td>4QPJ</td><td>P-His(CptA)→Asp(CtrA)</td><td>ChpT bears a DHp domain P-Asp(CtrA)→His(CptA) remains to be confirmed</td></tr><tr><td rowspan="2">HKs undergoing auto-phosphorylation</td><td rowspan="2">I</td><td>CpxA: ATP††† (1)</td><td>5.0†††</td><td>3.4</td><td>0.82</td><td>4BIW</td><td>ATP→His(CpxA)</td><td>AMP-PCP was used as non-hydrolizable ATP analogs (between phosphates β and γ there is a C atom substituting the O)</td></tr><tr><td>HK853/ EnvZ chimera: ATP<sup>14</sup> (1)</td><td>5.3†††</td><td>3.6</td><td>0.72</td><td>4KP4</td><td>ATP→His(HK853)</td><td></td></tr><tr><td>RR-mediated phosphorylation of imidazole</td><td>I</td><td>Imidazole: CheY¶¶¶ (2)</td><td>5.1; 5.1</td><td>3.3; 3.3</td><td>0.71; 0.73</td><td>5DKF</td><td>P-Asp(CheY)→imidazole</td><td>In vitro engineering of a phosphorelaying system from CheY to PhoR, using imidazole as a rudimentary HPt</td></tr><tr><td rowspan="3">HK:RR complexes</td><td rowspan="3">II</td><td>HK853: RR468**** (2)</td><td>7.6; 8.0††††</td><td>5.7; 6.1††††</td><td>0.51; 0.48††††</td><td>3DGE</td><td>P-His(HK853)→Asp(RR468)</td><td>HK:RR complex, snapshot of the phosphatase state according to the authors.</td></tr><tr><td>DesK: DesR†††† (2)</td><td>7.6§§§§</td><td>5.8§§§§</td><td>0.47§§§§</td><td>5IUK</td><td>P-His(DesK)→Asp(DesR)</td><td>HK:RR complex in the phosphotransferase state</td></tr><tr><td>CheA3: CheY6¶¶¶¶ (1)</td><td>7.3</td><td>ND††</td><td>0.52*****</td><td>3KYJ</td><td>P-His(CheA3)→Asp(CheY6)</td><td>3KYJ is not phosphorylated on either partner; 3KYJ with P-His on CheA3 displays an unproductive P-His rotamer, otherwise confirming 3KYJ's 3D organization</td></tr></tbody></table>

<footnote>* Distance between reactive His(Nε2)-Asp(Oδ1) less than or greater than 6.5Å, and distance ratio [Mg2+-Asp(Oδ1) / Mg2+-His(Nε2)] less than or greater than 0.6</footnote>

<footnote>† ASU = asymmetric unit</footnote>

<footnote>‡ Each distance corresponds to the one measured in each one of the independently refined complexes in the ASU.</footnote>

<footnote>§ Distances are reported in the same order as in the previous column, with correspondence among same individual complexes.</footnote>

<footnote>¶ Bacillus subtilis (Varughese et al., 2006)</footnote>

<footnote>** Saccharomyces cerevisiae (Zhao et al., 2008)</footnote>

<footnote>††Arabidopsis thaliana (Bauer et al., 2013)</footnote>

<footnote>‡‡ Not determined: no phosphoryl group or phosphoryl-mimetic present in the structure.</footnote>

<footnote>§§ Brucella abortus (Willett et al., 2015)</footnote>

<footnote>¶¶ Only one of the two complexes in the ASU (ChpTchainA:CtrAchainC, display the reactive His and Asp properly oriented poised for reacting. Only one distance is thus recorded.</footnote>

<footnote>*** No Mg2+ cation was actually bound on CtrA. A metal atom was modeled by superimposing the structure of RR468 with bound BeF3 (PDB 3GL9), one of the top ranking structures in multiple structural alignments with CtrA (DALI Z score 12.9, rmsd 0.8Å superimposing 98 αCs).</footnote>

<footnote>††† Escherichia coli (Mechaly et al., 2014)</footnote>

<footnote>‡‡‡ In the autophosphorylation complexes the distance is recorded between His(Nε2) and the position of the O between phosphates β and γ.</footnote>

<footnote>§§§ Thermotoga maritima (HK853) / E. coli (EnvZ) (Casino et al., 2014)</footnote>

<footnote>¶¶¶ E. coli (Page et al., 2015)</footnote>

<footnote>**** Thermotoga maritima (Casino et al., 2009)</footnote>
†††† These distances have been calculated by superimposing the phosphorylation-mimetic structure of RR468 alone with bound BeF3 (PDB 3GL9), onto 3DGE, in order to use better estimations of the receiver Asp position as well as of the phosphorus atom.

††††B. subtilis, this report.

§§§§ Glu188 present in the crystal structure was substituted by wt P-His following superposition of wt phosphorylated DesKC (PDB 5IUM) onto one of the phosphotransferase DesKCH188E:DesR-REC complexes in the ASU (chains A-B:E). Distances are reported after energy minimization. See Materials and methods.

¹¹¹¹Rhodobacter sphaeroides (Bell et al., 2010)

**** The position of Mg2+ cation was modeled using PDB 4TMY as template.

DOI: 10.7554/eLife.21422.019

## Specificity determinants within a loose and slippery interface

The small number of interactions and low surface complementarity that we observe in the several DesKC:DesR complexes (Figure 2F) are consistent with the observed 'slippery' nature of the protein:protein interface, displaying significant shifts in the relative positions of both partners. Further supporting that this ensemble of conformations is functionally relevant, molecular dynamics simulations of the phosphotransferase complex (Figure 2—figure supplement 2) recapitulated the gliding shifts of one partner with respect to the other, reaching the different arrangements observed in the crystal structures at different moments of the simulations (Figure 2—figure supplement 2B). Loose HK:RR association architectures can actually be observed in other TCS-related complexes (Zapf et al., 2000; Casino et al., 2009; Willett et al., 2015), suggesting it is a universal feature. This explains the considerable promiscuity reported in different TCSs in in vitro assays, if long enough incubation times are allowed for a given HK to phosphotransfer to different RRs, including surrogate partners (Skerker et al., 2005). HK:RR interactions thus ensure high specificity, but also avoiding exceedingly high stability, which would hinder proper functional modulation and, ultimately, signal transmission altogether. Evolutionarily selected to be loose, slippery interfaces are inherently tolerant for substantial sequence variation of protein:protein contacting residues. Such tolerance is consistent with the extensive degeneracy found in TCS pairs (Podgornaia and Laub, 2015), likely a general feature of HK:RR interactions. However, a permissive tendency for substitutions would appear to be contradictory with the high protein:protein specificity requirements of these systems. The DesK:DesR structures now provide with molecular details explaining why epistasis is found as a key property of HK:RR interfaces (Podgornaia and Laub, 2015). Epistatic amino acid changes are those that depend on the presence of other substitutions, leading to different effects in combination than individually. While interfacing residue replacements on any one of the TCS components can easily be accommodated due to the loose character of the protein:protein association, epistatic substitutions in the same component and/or in the other partner, are called to play essential stabilizing or destabilizing compensatory roles. Epistatic substitutions are thus selected, such that partner recognition and non-cognate discrimination are maintained, but without leading to overstable, non-functional complexes.

## A novel conformation of the response regulator when complexed with the histidine kinase

Our data indicate that the HK is largely unable to discriminate binding to RR or P-RR on the basis of structural determinants. The structures show the RR REC domain adopting identical conformations in both phosphatase and phosphotransferase complexes, presenting active- and inactive-like features in different regions of the protein, arranging the RR's active site for phosphoryl group migration (Figure 3). The two complexes also show equivalent HK:RR interfacing surfaces, with all major contacts preserved. Among the few differential contacts, they mainly implicate the HK's CA domain, but unlikely to grant discrimination ability given the highly similar conformation of the RR in those regions.

Calorimetry data are consistent with the crystallographic evidence, revealing only modest differences between DesR-REC:DesKC<sub>H188E</sub> and DesR-REC:DesKC<sub>H188E</sub> binding affinities (Figure 1 and Table 2), further supporting that both HK functional states interact with DesR-REC with similar strengths. Collecting calorimetric data by titrating DesKC with P-DesR-REC would be interesting to directly compare with the unphosphorylated form, but binding confounds with simultaneous P-RR dephosphorylation, and the addition of BeF<sub>3</sub><sup>-</sup>/Mg<sup>2+</sup> to mimic RR phosphorylation resulted in protein
aggregation. We did succeed in engineering a different phosphatase-trapped DesKC mutant, by insertion of additional heptad repeats in the coiled-coil segment (to be published elsewhere), indeed displaying equivalent DesR-REC-binding $K_D$ figures compared to DesKC<sub>H188E</sub>. HK:RR binding data in other TCSs are scarce, particularly exploring differential affinities according to the proteins' phosphorylation status. The few cases dealing with HKs that switch between phosphatase/phosphotransferase-competent states, do provide further evidence in support of non-discrimination. EnvZ binds OmpR and P-OmpR with equivalent affinities (Yoshida et al., 2002), as well as PhoQ binding to PhoP or P-PhoP (Castelli et al., 2003). Future studies shall uncover the relevance of additional mechanisms at the cell level that might favor particular HK:RR pairs to form, such as absolute and relative TCS protein concentrations, stability of non-phosphorylated HK:RR species in complex, subcellular spatial dynamics of each partner, and precise time-courses for HK and RR phosphorylated species to appear/disappear.

## A HK:RR concerted switch controls efficiency along the signaling pathway

Considering all the evidence, we can now propose a conceptual model explaining the molecular workings of the DesK/DesR signaling pathway (Figure 6). When the pathway is turned off, HK autophosphorylation is inhibited by: 1) hindering CA domains mobility, holding ATP substrates far from the His<sub>188(HK)</sub> phosphorylation site; and, 2) DHp helical rotation, making the His<sub>188(HK)</sub> side chains inaccessible (Figures 4D and 6A, Video 1). This 'open' autokinase-off configuration is linked to a folded N-terminal coiled-coil, in which DesK can bind P-DesR and trigger phosphatase catalysis. Linked to P-DesR phosphate lid opening, Asp<sub>189(HK)</sub> binds to Arg<sub>84(RR)</sub>, and Phe<sub>82(RR)</sub> becomes inserted into a DHp pocket (Figure 6A, Video 1). Phe<sub>82(RR)</sub> is thus removed from the closed position, which would otherwise clash with an incoming nucleophilic water molecule (Figure 4—figure supplements 2 and 3) that brings about dephosphorylation (Pazy et al., 2010). DesK Gln<sub>193(HK)</sub> is well located to position the catalytic water to perform dephosphorylation, consistent with the essential role of the homologous Gln residues in the phosphatase reactions mediated by the phosphatase CheZ (Zhao et al., 2002) or by the HK NarX (Huynh et al., 2010). The rotational movement of His<sub>188(HK)</sub>, placing it at >18 Å away from the reactive Asp<sub>54(RR)</sub>, makes phosphoryl back-transfer P-DesR→DesK impossible (Figure 4E). Although sequence divergence makes the DHp rotational shift unlikely in HisKA HKs (Marina et al., 2005), DHp kinking by means of a conserved proline (Casino et al., 2014; Mechaly et al., 2014) could regulate the separation of the phosphorylatable His from the reaction center, analogously to the HisKA_3 rotational switch. Additional mechanisms that stabilize an 'open' state of HisKA HKs via extensive CA-DHp interfaces, have been reported to favor the phosphatase reaction (Dubey et al., 2016) further supporting our hypotheses.

This DHp helical switch is thus coupled to signal sensing, such that sensory domain rearrangements lead to intra-cytoplasmic HK coiled-coil disruption. Possibly implying trans-membrane helix tilting, or other signal-triggered helical reconfigurations, a direct connection of sensory to intracytoplasmic α-helices can be posited. This switch can be readily understood as an order-disorder transition, which has been observed in HKs and other modular proteins, to work as an intraprotein signaling mechanism (Schultz and Natarajan, 2013). Helical rotational shifts have also been uncovered in structural elements within chemotaxis and TCSs, such as in HAMP domains, likely revealing common underlying transduction mechanisms eventually controlling their physiologic directionality (Airola et al., 2013).

HK activation leads to autophosphorylation and subsequent switching to its phosphotransferase state (Figure 6B). Comparing the conformation of the phosphorylatable histidine between HKs performing auto-phosphorylation (Casino et al., 2014; Mechaly et al., 2014), and the phosphotransferase complex (this study), the His residue shows ~180° flipping of its imidazole side chain (Figure 6B). During auto-phosphorylation, His<sub>Nδ1</sub> has been observed to H-bond with the ultra-conserved acidic His + 1 residue, which corresponds to Asp<sub>189(HK)</sub> in DesK. Such H-bond enhances His Nε2 nucleophilicity (Quezada et al., 2005), and indeed a carboxylate at His + 1 is critical for autophosphorylation in a number of HKs (Willett and Kirby, 2012). Once phosphorylated, the His disrupts its contact with this neighbor carboxylate group, stabilizing the P-HK species and leading to His χ-2 dihedral angle flipping (Figure 6C). Subsequent re-protonation of His<sub>188</sub>Nδ1, triggered by RR binding, is consistent with P-N bond weakening, as a first step towards phosphotransfer (Figure 4E). This
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Conceptual model of TCS unidirectional signaling. (A) P~DesR~ dephosphorylation catalysis is favored by the opening of the RR phosphorylation lid, via Phe<sub>82(RR)</sub> insertion into an HK hydrophobic pocket, and an ionic interaction between Asp<sub>189(HK)</sub> and Arg<sub>84(RR)</sub>. The Phe<sub>82(RR)</sub> movement allows for entry of a nucleophilic water (yellow sphere) positioned by Gln<sub>193(HK)</sub>. In the phosphatase state the phosphorylatable His<sub>188(HK)</sub> is buried within the DHp domain core, avoiding autophosphorylation and phosphoryl back-transfer reactions. (B) Signal-mediated pathway activation triggers the cogwheel rotation of the HK DHp α-helices, via coiled-coil disruption. The His<sub>188(HK)</sub> becomes exposed to the solvent, CA domains are released, and autophosphorylation thus enabled. Binding to unphosphorylated RR positions the P~His<sub>188(HK)</sub> into a suitable orientation for phosphoryl-transfer to occur through a loose (dissociative) nucleophilic substitution. Movement of Thr<sub>80(RR)</sub> is coupled with phosphoryl group migration, triggering RR's phosphate lid closure and P~RR release. (C) Reaction scheme summarizing the autophosphorylation and phosphotransfer reactions. Stabilization of a protonated Nδ1 tautomer of the reactive histidine, is required in autophosphorylation. In keeping with the imidazole aromaticity, the nucleophilicity of Nε2 is thus finely regulated, allowing for phospho-acceptor/donor roles of the His.</figcaption></figure>

DOI: 10.7554/eLife.21422.020

mechanism might be critical in dissociative P~HK→RR reactions, which lack the simultaneous nucleophilic attack that the receiver Asp<sub>(RR)</sub> carries out in associative ones.

Our model predicts an important role of the highly conserved Thr<sub>80(RR)</sub> in stabilizing an intermediary metaphosphate group, ideally positioned to accompany the migrating phosphoryl moiety. This threonine has been substituted in CheY (Ganguli et al., 1995; Appleby and Bourret, 1998) and DosR (Gautam et al., 2011), suggesting it is not directly engaged in catalysis, but does play important roles in phosphoryl-transfer in vitro and proper in vivo signaling. Systematic mutagenesis studies substituting this threonine in different TCS pairs are needed to be conclusive about its role in phosphotransfer. Once the phosphotransferase reaction is complete, immediate dephosphorylation of
<figure><img src="image_9.png" /></figure>

**Video 1.** Animation illustrating the phosphotransferase to phosphatase transition, within the model of TCS-mediated phosphorylation and dephosphorylation of the RR component. As an animated support to main Figure 6, the fixed ground states correspond to crystal structures, and the moving intermediates in between obtained by linearly interpolated morphing, as implemented in Pymol v 1.8.2

(Schrodinger, LLC, 2015). This animation does not include the kinase's auto-phosphorylation step, shown at the beginning of Figure 6B. Structures of phosphorylated ('active') and non-phosphorylated ('inactive') DesR, when shown separate from the HK DesK, correspond respectively to PDB structures 4LE0 and 4LE1 (Trajtenberg et al., 2014). On the other hand, the fixed positions of P~DesK in complex with DesR (phosphotransferase state), and of DesK in complex with P~DesR (phosphatase state), correspond to the structures reported in this work (PDBs 5IUK and 5IUN, respectively). The reactive His<sub>188(HK)</sub> on DesK is depicted as spheres, as well as the active Asp<sub>54(RR)</sub> and Phe<sub>82(RR)</sub> on DesR. The rest of the proteins are rendered as cartoons. Chains A and B of DesK are colored in green and yellow, respectively. DesR is blue when not phosphorylated, and magenta once activated by phosphorylation. The P~His was modeled into DesK as explained in the text and Materials and methods. P~DesR corresponds to the phosphomimetic BeF<sub>3</sub><sup>-</sup> bonded to Asp<sub>54(RR)</sub>.

DOI: 10.7554/eLife.21422.021

the P~Asp<sub>54(RR)</sub>, as well as P~DesR→DesK back-transfer, are minimized. The β4α4 loop in DesR plays a key role in these controls, remodeling its configuration in concert with DesK-association. This loop, acting as a phosphate lid, is open in the phosphotransferase complex, enabling the transfer reaction (Video 1). Yet, the phosphate lid is expected to close in, following the phosphoryl migration during transfer, due to the strong phosphoryl-Thr<sub>80(RR)</sub> interaction. Phosphotransfer thus results in a fully activated RR configuration, with a closed phosphate lid that ensures the release of a stable P~DesR species (Figure 6B and Video 1), as observed in active free RR structures (Trajtenberg et al., 2014).

The β4α4 loop amino acid sequence is not well conserved among different RRs and RR families (Page et al., 2016), despite its direct implication in forming the HK:RR interface. Conspicuous examples of variable residues can be mapped to Thr<sub>81(RR)</sub>, well positioned to interact with the phosphoryl group of P~DesK, or yet Phe<sub>82(RR)</sub> and Arg<sub>84(RR)</sub>, key players in the open/closure switch of DesR's phosphate lid. Some of these positions had already been pinpointed as highly variable, yet functionally relevant (Thomas et al., 2013; Immormino et al., 2016; Page et al., 2016), recognizing positions T + 1 and T + 2 (starting from T, the conserved position Thr<sub>80(RR)</sub>) among others, as important modulators of RRs' catalytic activities. Particular sequence signatures of these variable motifs, have been found to correlate with RR families, in turn connected to the distinct physiological constraints of each family's signaling functions (Page et al., 2016). Our data now illustrate an additional selective pressure dimension to such variable positions, linked to the evolution of the HK:RR interface. Beyond sequence variability, the RR β4α4 loop is bound to play functional roles, as it is always located at the HK:RR interface in available structures. Muta-

genesis studies with functional readouts including HK-mediated activities are needed, together with a larger number of 3D structures of HK:RR complexes from different TCS families, in order to delimit universal TCS mechanisms and system-dependent variations.

## Associative vs dissociative phosphotransfer: evolution of signaling directionality

Regulating the direction of information flow is central to signaling performance. The energies of P-N and P-O bonds engaged in phosphoryl-transfer reactions among histidine and aspartate, appear to be equivalent, explaining a few known reversible systems (Burbulys et al., 1991). However, and except for such systems where P~Asp→His phosphorylation has been selected (such as in phosphoryl- elays), phosphoryl back-transfer (P~RR→HK) is typically minimized. Alternative sources of reactional asymmetry, other than bonding energetics, must be at play. Taking into account that phosphotransfer is a nucleophilic substitution at phosphorus, the reaction coordinate distance (i.e. between phosphorus and the entering atom, before the transfer reaction has begun) strongly suggests that P~DesK→DesR phosphotransfer occurs through a largely dissociative, or 'loose' mechanism
(Lassila et al., 2011): a reaction coordinate distance of ≥4.9 Å provides room for a fully dissociated metaphosphate intermediate (Mildvan, 1997). This distance is directly correlated with the one separating side chain N and O atoms, respectively on the reactive His and Asp residues. Considering van der Waals radii, full dissociation of the phosphate intermediate is predicted for N-O distances of ≥6.6 Å. This threshold distance classifying phosphotransfer reactions into loose vs tight nucleophilic substitutions, nicely discriminates signaling complexes involved in, respectively, phosphorylations from the histidine imidazole group, or towards it (Figure 5D, Table 6). The His-Asp distance thus correlates with associative vs dissociative phosphotransfer, but also with the symmetry of the Mg²⁺ cation with respect to the entering and leaving atoms (Figure 5C and D). Increased back-transfer has been observed when the cation is removed (Albanesi et al., 2004; Lukat et al., 1990; Shi et al., 1999), further supporting that Mg²⁺ asymmetry favors irreversible P~His→Asp transfer reactions. Furthermore, mutations that likely increase the distance between TCS partners, without impeding their association, have been reported to affect differentially the forward- and back-transfer directions in otherwise reversible phosphoryl-transfer steps (Janiak-Spens et al., 2005).

The DesK:DesR structures disclose the catalytic centers for TCS-mediated dephosphorylation and phosphotransfer reactions, built with contributing residues from both partners. Moreover, the HK serves as the driving component to modulate the construction of a given active site, selecting a conformation of the RR that is prone for transfer reactions on its catalytic aspartate. Except for the phosphate lid, the REC domain bound to the HK largely resembles the RR's active form, shifting the position of the β1α1 loop and helix α1, now able to recruit the Mg²⁺ cation essential for phosphoryl-transfer catalysis. The phosphatase/phosphotransferase transition, switches the distance and orientation of the HK histidine side chain with respect to the RR aspartate, enabling for the distinct reaction mechanisms to take place, and ultimately ensuring the unidirectional flow of information.

# Materials and methods

## Cloning and mutagenesis

For protein crystallization purposes, recombinant protein expression plasmid pACYC-DesKC<sub>STAB</sub>: DesR-REC was generated by sub-cloning DesKC<sub>STAB</sub> from pHPKS/Pxyl-desKSTA (Saita et al., 2015) into pACYC-DesKC<sub>H188E</sub>:DesR-REC (Trajtenberg et al., 2014) using RF cloning (Unger et al., 2010) with primers STAB_F (5'-CCTGTATTTTCAGGGATCC GGTATTATAAAACTTCGCAAG-3') and STAB_R (5'-GTCAGACACTGTAATCACAACTTCCTTCCAG-3'). The two expression cassettes present in the pACYC-DesKC<sub>STAB</sub>:DesR-REC encode the two recombinant proteins each fused to a His-tag and a TEV protease cleavage site.

For phosphotransfer assays, two DesR point mutants Phe82Ala and Arg84Ala were designed by mutagenesis of pQE80-DesR-REC (Trajtenberg et al., 2014), using mutagenic primers DesR_Fw (5'-GATTAGTATTTATTGCAGAAGATCAGCAAATGC-3'), F82A_Rev (5'-CGGGTCTGGCGGCGGTGG TTAAG ATGATAATTTTG-3') and DesR_R84A_Rev (5'-GAAAGTAACCGGGTGCGGCGAAGGTGG TTAAG-3'), resulting in plasmids pQE80-DesR-REC<sub>F82A</sub> and pQE80-DesR-REC<sub>R84A</sub>. Cloning procedures, as well as protein expression, crystallization and structure solution methods, are described in more detail at Bio-protocol (Imelio et al., 2017).

## Protein expression, purification and crystallization

Recombinant proteins to be used in crystallization were co-expressed from pACYC-DesKC<sub>STAB</sub>: DesR-REC and pACYC-DesKC<sub>H188E</sub>:DesR-REC plasmids, transformed in *Escherichia coli* strain BL21 (DE3) (Novagen). For functional assays proteins were expressed from plasmids pQE32-DesKC (Albanesi et al., 2009), pQE80-DesR-REC (Trajtenberg et al., 2014), pQE80-DesR-REC<sub>F82A</sub> and pQE80-DesR-REC<sub>R84A</sub>, in *E. coli* TOP10F' (Invitrogen). Proteins were expressed and purified as previously described (Albanesi et al., 2009; Trajtenberg et al., 2014). For the DesK-DesR complexes the last size exclusion chromatography step (HiLoad 16/60 Superdex 75 preparation grade column; GE Healthcare) was used to select those fractions corresponding to the complex only.

Hanging-drop crystallizations were performed at 20°C in Linbro plates. The DesKC<sub>STAB</sub>:DesR-REC complex (10 mg/mL protein prepared in 5 mM AMP-PCP, 10 mM MgCl₂) crystallized in 1 mL mother liquor [30% (w/v) PEG 4000, 0.1 M Tris.HCl pH 8.5, 0.2 M Li₂SO₄]. Protein drops were set up by mixing 0.8 µL mother liquor + 2 µL protein solution + 1.2 µL additive solution [27% (v/v) PEG 600, 0.1 M
MES pH 6.5. 0.15 M MgSO₄, 5% (v/v) glycerol. Initial crystals were optimized by microseeding. Cryo-protection was achieved by slowly adding 4 µl of cryoprotection solution [32% (w/v) PEG 4000, 0.1 M Tris.HCl pH 8, 0.2 M Li₂SO₄, 20 mM MgCl₂, 18 mM BeF₃⁻, 5 mM AMP-PCP, 15% (v/v) glycerol], then soaked in 100% cryoprotection solution and frozen in liquid N₂. The DesKC<sub>H188E</sub>:DesR-REC complex (8.3 mg/mL protein prepared in 5 mM AMP-PCP, 20 mM MgCl₂) crystallized in mother liquor 18% (w/v) PEG 3350, 0.3 M tri-potassium citrate. Protein drops were set up by mixing 2 µL protein + 2 µL mother liquor. Cryo-protection was achieved by quick soaking in 20% (w/v) PEG 3350, 0.3 M tri-potassium citrate, 5 mM AMP-PCP, 25% (v/v) glycerol, and 20 to 150 mM MgCl₂ + 0 to 5 mM BeF₃⁻. wt P~DesKC was crystallized as described (Albanesi et al., 2009), except that 5 mM ATP was used throughout instead of AMP-PCP.

## X ray diffraction data processing, crystal structure solution, refinement and analysis

Single crystal X ray diffraction was performed with a copper rotating anode home-source (Protein Crystallography Facility, Institut Pasteur de Montevideo) or synchrotron radiation (Soleil, France). Data processing was performed with autoPROC (Vonrhein et al., 2011). Structures were solved by molecular replacement (McCoy et al., 2007), using an in silico-generated model of a DesKC:DesR complex (Trajtenberg et al., 2014) as search probe. The other domains were then manually located in the electron density maps, and complete models were rebuilt using Coot (Emsley et al., 2010) and refined with Buster (Bricogne et al., 2009). Validation was done throughout and towards the end of refinement, using MolProbity tools (Chen et al., 2010). DesKC:DesR-REC surface complementarity was calculated according to (Lawrence and Colman, 1993), resulting in figures of ~0.6. Visualization of protein models and structural analyses, figure rendering and morphing for animation were performed with Pymol (Schrodinger, LLC, 2015). Software for data processing, structure determination and analysis was provided by the SBGrid Consortium (Morin et al., 2013).

## Small angle X-ray scattering data acquisition and analysis

SEC-SAXS experiments were carried out at beamline SWING (Soleil synchrotron, France). Purified DesK-DesR complexes samples were brought to 10 mg/mL in 50 mM Tris.HCl pH 8.0, 300 mM NaCl, and injected into a Superdex 75 5/150 GL column (GE Healthcare Biosciences) equilibrated in the same buffer. Samples were eluted at constant flow (0.15 mL/min) and loaded into capillary cell for X-ray exposure. Recorded frames were processed with Foxtrot (David and Pérez, 2009) following standard procedures. Subsequent analyses of the scattering data were performed with tools from the ATSAS (Konarev et al., 2006). Forward scattering $I_{(0)}$ and the particle's radius of gyration ($R_g$) were estimated using the Guinier approximation. $I_{(0)}$ real was calculated as the extrapolated intensity at zero scattering angle. Pair-distance distribution functions were calculated from the scattering patterns with GNOM, which also provides maximum particle dimension ($D_{max}$) and $R_g$ (real) values. Molecular mass (MM) of particles was estimated from the volume-of-correlation ($V_c$) values (Rambo and Tainer, 2013). Theoretical scattering patterns from atomic models were calculated and fitted to experimental curves using CRY SOL. Mixture analysis was performed with OLIGOMER.

## Isothermal titration calorimetry

Isothermal titration calorimetry (ITC) assays were performed on a VP-ITC (MicroCal Inc., Northampton, MA). Titrations consisted of an initial injection (1 µL), followed by 25–30 injections of 10 µL of DesR-REC (ligand) on the cell containing the different DesKC variants (partner). Assays were carried out at 15°C in a buffer containing 20 mM Tris.HCl pH 8, 0.3 M NaCl, 10 mM MgCl₂ and 0.5 mM AMP-PNP. The heat released by the dilution of the ligand was determined injecting the ligand, on the cell containing the working solution without partner, using the same sequence of injections. The concentrations used for the experiments were between 25–30 µM for DesKC and 350–400 µM of DesR-REC and each titration was done in duplicate. The data were analyzed with MicroCal Origin version 7 software (MicroCal Software Inc.), after manual baseline correction, and subtraction of heat due to ligand dilution. Binding isotherms were fitted to a two independent sequential site model.
## Sequence-based direct coupling analysis of the HK:RR complex

Sequences of HK and RR pairs, with the same architecture as DesK and DesR (families HisKA_3 and NarL, respectively), belonging to the same operon, were selected from the Uniprot database using hmmsearch (*Eddy, 1998*). Hidden Markov model-profiles for each domain were obtained from Pfam (*Finn et al., 2016*). Redundancy was filtered using a cutoff of 90% (*Li and Godzik, 2006*) resulting in a total of 3318 sequences including HKs (DHp+CA) and RRs (REC+HTH). The concatenated HK:RR sequences were aligned with hmmalign (*Eddy, 1998*) and manually curated, removing sequences with large insertions. DI (direct information) of each pair of residues of the alignment was calculated using mfDCA matlab script (*Morcos et al., 2011*) and we only take in consideration coevolving pairs of residues farther than five positions apart in sequence. The same procedure was followed for the HisKA and PhoB families for comparison purposes. A similar trend in equivalent 3D positions for covariant residue pairs among different HK and RR families was thus confirmed (*Podgornaia et al., 2013; Skerker et al., 2008; Weigt et al., 2009*).

## In silico modeling and molecular dynamics simulations

The phosphotransferase complex was built by replacing the DesKC<sub>H188E</sub> dimer of 5IUK with wt P~DesKC (PDB 5IUM). This replacement resulted in model with barely no clashes, except for P~His<sub>188(HK)</sub> slightly bumping into Phe<sub>82(RR)</sub> and Arg<sub>84(RR)</sub>. Energy minimization was then performed on such constructed complex with all-atom constraints in Rosetta (*Das and Baker, 2008*), strongly restraining shifts to maintain the experimental coordinates. A final wt P~DesKC:DesR-REC model was thus obtained, with optimal stereochemical geometry and no clashes (0.09 rmsd between the energy minimized model and the pdb 5IUM experimental structure, aligning all 1930 DesK atoms). To further evaluate the minimized wt P~DesKC:DesR-REC model (see *Figure 4A*), molecular dynamics simulations were performed, observing that its conformation is stable throughout the trajectory (see *Figure 4B*), especially considering the positions of His<sub>188(HK)</sub>, Asp<sub>54(RR)</sub> and Thr<sub>80(RR)</sub> in contact with the phosphate.

For molecular dynamics (MD) calculations, the missing ATP lid loop within the CA domain in chain A (residues 328 to 335) were reconstructed by using a previously reported high resolution structure (PDB 3EHG) as template. The other CA domain (in chain B, residues 241–367) was deleted. The missing N-terminal portion of chain B helix α1, was reconstructed and extended up to residue Lys<sub>155(HK)</sub> using 5IUK as a template.

For computational efficiency only DesK residues Lys<sub>155(HK)</sub> to Asn<sub>368(HK)</sub> in chain A and Lys<sub>155(HK)</sub> to Ser<sub>239(HK)</sub> in chain B were considered. Only His<sub>188(HK)</sub> in chain B was set to be phosphorylated. His<sub>335(HK)</sub> was protonated at position Nδ1, while other histidines were protonated at Nε2, to preserve the interaction network. The ATP moiety and Mg<sup>2+</sup> cation present within the CA domain were both kept in the model. The structure of DesR-REC spanned residues Ser<sub>0(RR)</sub> to Leu<sub>131(RR)</sub>. The Mg<sup>2+</sup> cation in the active site of DesR-REC and the three crystallographic water molecules of the metal's coordination sphere, were included in the model.

Acetyl and N-methylamide capping groups were added to N- and C-terminal residues, which did not correspond to the real terminals in the full protein sequences. The protein complex was solvated within an octahedral box of 15 Å from the solute, and sodium counterions were added to preserve the electroneutrality of the system. The AMBER force field ff14SB (*Maier et al., 2015*) was used to represent the amino acids. The parameters for His<sub>188(HK)</sub> phosphorylated at NE2 with an unprotonated phosphate group were taken from (*Homeyer et al., 2006*). The ATP molecule was described with reported parameters (*Meagher et al., 2003*). The TIP3P model (*Jorgensen et al., 1983*) was used for water molecules. Monovalent ions were treated with previously described parameters (*Joung and Cheatham, 2008*), while Mg<sup>2+</sup> ions were modeled using the compromise set of parameters for the 12–6 non-bonded potential in TIP3P (*Li et al., 2013*).

All calculations were performed with the GPU version of AMBER14 (*Salomon-Ferrer et al., 2013*). Initially, the whole system was relaxed by energy minimization, then it was equilibrated for 0.2 ns in the NVT ensemble imposing harmonic positional restrains of 10 kcal mol<sup>-1</sup> Å<sup>-2</sup> on the protein atoms. A reference temperature of 300 K was set by coupling the system to the Langevin thermostat (*Pastor et al., 1988; Wu and Brooks, 2003*) with a friction constant of 50 ps<sup>-1</sup>, which approximates the physical collision frequency for liquid water (*Izaguirre et al., 2001*). A 10 Å cut-off was used for non-bonded interactions, while long-range electrostatics were evaluated using Particle
Mesh Ewald (PME) (Darden et al., 1993; Essmann et al., 1995). A time step of 2 fs was used and all bonds involving hydrogen atoms were restrained using the SHAKE algorithm (Miyamoto and Kollman, 1992; Ryckaert et al., 1977). Production simulations were performed in the NPT ensemble. The pressure was kept at 1 atm by means of the Berendsen barostat (Berendsen et al., 1984). Snapshots were recorded every 5 ps for analysis.

## Autodephosphorylation and phosphotransfer assays

To purify the phosphorylated species of DesR-REC and DesR-REC<sub>F82A</sub>, 600 µM of both recombinant proteins were auto-phosphorylated using 50 mM acetyl phosphate in a buffer containing 25 mM Tris.HCl pH 8, 300 mM NaCl and 30 mM MgCl<sub>2</sub>, at room temperature. Reactions were stopped by adding EDTA to a final concentration of 50 mM, and buffer exchanged by using a PD Minitrap G-25 (GE Healthcare) desalting column. Auto-dephosphorylation assays of P~DesR-REC variants were performed at 30 µM of protein concentration and incubated in the presence of 30 mM MgCl<sub>2</sub>. At different time points the reactions were stopped by adding SDS-PAGE sample buffer. For each time point DTT was added to a final concentration of 25 mM and incubated for 5 min. DTT in excess was blocked with 40 mM iodo-acetamide and loaded in a Phos-tag acrylamide SDS-PAGE, as described before (Trajtenberg et al., 2014), Coomassie blue-stained gels were scanned and quantification of all reactions was done by densitometry using ImageJ (Schneider et al., 2012).

To perform phosphotransfer assays (Figure 4E), DesKC was phosphorylated by incubation with 10 mM ATP and 20 mM MgCl<sub>2</sub> at 100 µM of protein, for 1 hr at room temperature. P~DesKC was then purified by size exclusion chromatography S75 10/300 (GE Healthcare), equilibrated in 20 mM Tris.HCl pH 8, 300 mM NaCl. Phosphotransfer reactions were performed in triplicate by combining 30 µM P~DesKC with equimolar concentrations of DesR-REC (both wild type and DesR-REC<sub>F82A</sub>), at 25°C in 20 mM Tris.HCl pH 8, 300 mM NaCl, 30 mM MgCl<sub>2</sub>. Reactions were stopped at different time points as described above. Samples were then separated by Phos-tag SDS-PAGE.

In order to assess unidirectionality (Figure 5A), the reverse and forward reactions were analyzed. The forward reaction was prepared as described above using P~DesKC as phosphodonor for 1 min. On the other hand, the reverse of the reaction was measured by incubating 30 µM of P~DesR-REC (wild type or point mutants DesR-REC<sub>F82A</sub> or DesR-REC<sub>R84A</sub>) with 30 µM DesKC, 50 mM acetyl phosphate and 30 mM MgCl<sub>2</sub>, for 1 min at room temperature. Reactions were stopped as described above and samples analyzed by Phos-tag SDS-PAGE and quantified by densitometry.

## Accession numbers

The X ray structures presented have been deposited in the wwPDB with accession codes 5IUJ (DesK-DesR complex in the phosphotransfer state with low Mg<sup>2+</sup> [20 mM]), 5IUK (DesK-DesR complex in the phosphotransfer state with high Mg<sup>2+</sup> [150 mM]), 5IUL (DesK-DesR complex in the phosphotransfer state with high Mg<sup>2+</sup> [150 mM] and BeF<sub>3</sub>'), 5IUM (phosphorylated wild type DesKC) and 5IUN (DesK-DesR complex in the phosphatase state).

Raw X ray diffraction data corresponding to each one of these structures are publicly available at SBGrid Data Bank (http://data.sbgrid.org) as dataset entries 399, 400, 401, 407 and 408.

## Acknowledgements

We want to thank the staffs at synchrotron beamlines Proxima 1 and 2, Soleil (France), especially William Shepard; and Daniela Albanesi for kindly providing plasmid pHPKS/Pxyl-desKSTA. The Institut Pasteur International Network is acknowledged for support through the IMiZA International Joint Unit. We acknowledge funding from Agencia Nacional de Investigación e Innovación (ANII), Uruguay (FCE2009_1_2679; FCE2007_219); Agence Nationale de la Recherche (ANR), France (PCV06_138918); Centro de Biología Estructural del Mercosur (www.cebem-lat.org) and Fondo para la Convergencia Estructural del MERCOSUR (COF 03/11).
# Additional information

## Funding

<table><thead><tr><td>Funder</td><td>Grant reference number</td><td>Author</td></tr></thead><tbody><tr><td>Agencia Nacional de Investigación e Innovación</td><td>FCE2009_1_2679</td><td>Felipe Trajtenberg<br>Alejandro Buschiazzo</td></tr><tr><td>Agencia Nationale de la Recherche</td><td>PCV06_138918</td><td>Alejandro Buschiazzo</td></tr><tr><td>FOCEM (MERCOSUR Structural Convergence Fund)</td><td>COF 03/11</td><td>Alejandro Buschiazzo</td></tr><tr><td>Centro de Biologia Estructural del Mercosur</td><td></td><td>Alejandro Buschiazzo</td></tr><tr><td>Agencia Nacional de Investigación e Innovación</td><td>FCE2007_219</td><td>Felipe Trajtenberg<br>Alejandro Buschiazzo</td></tr></tbody></table>

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

## Author contributions

FT, AB, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article; JAI, GO, Acquisition of data, Analysis and interpretation of data; MRM, Acquisition of data, Analysis and interpretation of data, Contributed unpublished essential data or reagents; NL, Acquisition of data, Contributed unpublished essential data or reagents; MAM, Analysis and interpretation of data, Drafting or revising the article, Contributed unpublished essential data or reagents; AEM, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article

## Author ORCIDs

Alejandro Buschiazzo, http://orcid.org/0000-0002-2509-6526

# Additional files

## Major datasets

The following datasets were generated:

<table><thead><tr><td>Author(s)</td><td>Year</td><td>Dataset title</td><td>Dataset URL</td><td>Database, license, and accessibility information</td></tr></thead><tbody><tr><td>Trajtenberg F, Buschiazzo A</td><td>2016</td><td>Crystal structure of phosphorylated DesKC</td><td>http://www.rcsb.org/pdb/search/structid-Search.do?structureId=5IUM</td><td>Publicly available at the RCSB Protein Data Bank (accession no: 5IUM)</td></tr><tr><td>Trajtenberg F, Imelio JA, Larrieux N, Buschiazzo A</td><td>2016</td><td>Crystal structure of the DesK-DesR complex in the phosphatase state</td><td>http://www.rcsb.org/pdb/search/structid-Search.do?structureId=5IUN</td><td>Publicly available at the RCSB Protein Data Bank (accession no: 5IUN)</td></tr><tr><td>Trajtenberg F, Imelio JA, Larrieux N, Buschiazzo A</td><td>2016</td><td>Crystal structure of the DesK-DesR complex in the phosphotransfer state with low Mg2+ (20 mM)</td><td>http://www.rcsb.org/pdb/search/structid-Search.do?structureId=5IUJ</td><td>Publicly available at the RCSB Protein Data Bank (accession no: 5IUJ)</td></tr><tr><td>Trajtenberg F, Imelio JA, Larrieux N, Buschiazzo A</td><td>2016</td><td>Crystal structure of the DesK-DesR complex in the phosphotransfer state with high Mg2+ (150 mM)</td><td>http://www.rcsb.org/pdb/search/structid-Search.do?structureId=5IUK</td><td>Publicly available at the RCSB Protein Data Bank (accession no: 5IUK)</td></tr><tr><td>Trajtenberg F, Imelio JA, Larrieux N, Buschiazzo A</td><td>2016</td><td>Crystal structure of the DesK-DesR complex in the phosphotransfer state with high Mg2+ (150 mM) and BeF3</td><td>http://www.rcsb.org/pdb/search/structid-Search.do?structureId=5IUL</td><td>Publicly available at the RCSB Protein Data Bank (accession no: 5IUL)</td></tr><tr><td>Trajtenberg F,</td><td>2016</td><td>Complex DesKC:DesR-REC (B.</td><td>http://dx.doi.org/10.</td><td>Publicly available at</td></tr></tbody></table>
<table><tr><td>Buschiazzo A</td><td>subtilis) - phosphotransferase state, low Mg2+</td><td>15785/SBGRID/399</td><td>the SBGrid database (dataset no. 399; X-ray diffraction raw data for PDB 51UJ)</td></tr><tr><td>Trajtenberg F, Buschiazzo A</td><td>2016 Complex DesKC:DesR-REC (B. subtilis) - phosphatase state</td><td>http://dx.doi.org/10.15785/SBGRID/400</td><td>Publicly available at the SBGrid database (dataset no. 400; X-ray diffraction raw data for PDB 51UN)</td></tr><tr><td>Trajtenberg F, Buschiazzo A</td><td>2016 Complex DesKC:DesR-REC (B. subtilis) - phosphotransferase state, high Mg2+</td><td>http://dx.doi.org/10.15785/SBGRID/401</td><td>Publicly available at the SBGrid database (dataset no. 401; X-ray diffraction raw data for PDB 51UK)</td></tr><tr><td>Buschiazzo A</td><td>2016 Phosphorylated DesKC (B. subtilis)</td><td>http://dx.doi.org/10.15785/SBGRID/407</td><td>Publicly available at the SBGrid database (dataset no. 407; X-ray diffraction raw data for PDB 51UM)</td></tr><tr><td>Trajtenberg F, Buschiazzo A</td><td>2016 Complex DesKC:DesR-REC (B. subtilis) - phosphotransferase state, high Mg2+ and BeF3-</td><td>http://dx.doi.org/10.15785/SBGRID/408</td><td>Publicly available at the SBGrid database (dataset no. 408; X-ray diffraction raw data for PDB 51UL)</td></tr></table>

## References

* Airola MV, Sukomon N, Samanta D, Borbat PP, Freed JH, Watts KJ, Crane BR. 2013. HAMP domain conformers that propagate opposite signals in bacterial chemoreceptors. *PLoS Biology* **11:e1001479. doi: 10.1371/journal.pbio.1001479, PMID: 23424282**

* Albanesi D, Mansilla MC, de Mendoza D. 2004. The membrane fluidity sensor DesK of *Bacillus subtilis* controls the signal decay of its cognate response regulator. *Journal of Bacteriology* **186:2655-2663. doi: 10.1128/JB.186.9.2655-2663.2004, PMID: 15090506**

* Albanesi D, Martín M, Trajtenberg F, Mansilla MC, Haouz A, Alzari PM, de Mendoza D, Buschiazzo A. 2009. Structural plasticity and catalysis regulation of a thermosensor histidine kinase. *PNAS* **106:16185-16190. doi: 10.1073/pnas.0906699106, PMID: 19805278**

* Appleby JL, Bourret RB. 1998. Proposed signal transduction role for conserved CheY residue Thr87, a member of the response regulator active-site quintet. *Journal of bacteriology* **180:3563-3569. PMID: 9657998**

* Attwood PV, Piggott MJ, Zu XL, Besant PG. 2007. Focus on phosphohistidine. *Amino Acids* **32:145-156. doi: 10.1007/s00726-006-0443-6, PMID: 17103118**

* Bauer J, Reiss K, Veerabagu M, Heunemann M, Harter K, Stehle T. 2013. Structure-function analysis of *Arabidopsis thaliana* histidine kinase AHK5 bound to its cognate phosphotransfer protein AHP1. *Molecular Plant* **6:959-970. doi: 10.1093/mp/sss126, PMID: 23132142**

* Bell CH, Porter SL, Strawson A, Stuart DI, Armitage JP. 2010. Using structural information to change the phosphotransfer specificity of a two-component chemotaxis signalling complex. *PLoS Biology* **8:e1000306. doi: 10.1371/journal.pbio.1000306, PMID: 20161720**

* Berendsen HJC, Postma JPM, van Gunsteren WF, DiNola A, Haak JR. 1984. Molecular dynamics with coupling to an external bath. *The Journal of Chemical Physics* **81:3684-3690. doi: 10.1063/1.448118**

* Bhate MP, Molnar KS, Goulian M, DeGrado WF. 2015. Signal transduction in histidine kinases: insights from new structures. *Structure* **23:981-994. doi: 10.1016/j.str.2015.04.002, PMID: 25982528**

* Bourret RB, Stock AM. 2002. Molecular information processing: lessons from bacterial chemotaxis. *Journal of Biological Chemistry* **277:9625-9628. doi: 10.1074/jbc.R100066200, PMID: 11779877**

* Bricogne G, Blanc E, Brandl M, Flensburg C, Keller P, Paciorek W, Roversi P, Smart OS, Vonrhein C, Womack TO. 2009. BUSTER. 2.8.0. Cambridge, United Kingdom: Global Phasing Ltd.

* Burbulys D, Trach KA, Hoch JA. 1991. Initiation of sporulation in *B. subtilis* is controlled by a multicomponent phosphorelay. *Cell* **64:545-552. doi: 10.1016/0092-8674(91)90238-T, PMID: 1846779**

* Casino P, Miguel-Romero L, Marina A. 2014. Visualizing autophosphorylation in histidine kinases. *Nature Communications* **5:3258. doi: 10.1038/ncomms4258, PMID: 24500224**

* Casino P, Rubio V, Marina A. 2009. Structural insight into partner specificity and phosphoryl transfer in two-component signal transduction. *Cell* **139:325-336. doi: 10.1016/j.cell.2009.08.032, PMID: 19800110**

* Castelli ME, Cauerhff A, Amongero M, Soncini FC, Vescovi EG. 2003. The H box-harboring domain is key to the function of the *Salmonella enterica* PhoQ Mg2+-sensor in the recognition of its partner PhoP. *Journal of Biological Chemistry* **278:23579-23585. doi: 10.1074/jbc.M303042200, PMID: 12702718**

* Chen VB, Arendall WB, Headd JJ, Keedy DA, Immormino RM, Kapral GJ, Murray LW, Richardson JS, Richardson DC. 2010. MolProbity: all-atom structure validation for macromolecular crystallography. *Acta Crystallographica Section D Biological Crystallography* **66:12-21. doi: 10.1107/S0907444909042073, PMID: 20057044**
Darden T, York D, Pedersen L. 1993. Particle mesh Ewald: An N-log(N) method for Ewald sums in large systems. *The Journal of Chemical Physics* **98**:10089-10092. doi: 10.1063/1.464397

Das R, Baker D. 2008. Macromolecular modeling with rosetta. *Annual Review of Biochemistry* **77**:363-382. doi: 10.1146/annurev.biochem.77.062906.171838, PMID: 18410248

David G, Pérez J. 2009. Combined sampler robot and high-performance liquid chromatography: a fully automated system for biological small-angle X-ray scattering experiments at the Synchrotron SOLEIL SWING beamline. *Journal of Applied Crystallography* **42**:892-900. doi: 10.1107/S0021889809029288

de Mendoza D. 2014. Temperature sensing by membranes. *Annual Review of Microbiology* **68**:101-116. doi: 10.1146/annurev-micro-091313-103612, PMID: 24819366

Dubey BN, Lori C, Ozaki S, Fucile G, Plaza-Menacho I, Jenal U, Schirmer T. 2016. Cyclic di-GMP mediates a histidine kinase/phosphatase switch by noncovalent domain cross-linking. *Science Advances* **2**:e1600823. doi: 10.1126/sciadv.1600823, PMID: 27652341

Eddy SR. 1998. Profile hidden Markov models. *Bioinformatics* **14**:755-763. doi: 10.1093/bioinformatics/14.9.755, PMID: 9918945

Emsley P, Lohkamp B, Scott WG, Cowtan K. 2010. Features and development of coot. *Acta Crystallographica Section D Biological Crystallography* **66**:486-501. doi: 10.1107/S0907444910007493, PMID: 20383002

Essmann U, Perera L, Berkowitz ML, Darden T, Lee H, Pedersen LG. 1995. A smooth particle mesh ewald method. *The Journal of Chemical Physics* **103**:8577-8593. doi: 10.1063/1.470117

Finn RD, Coggill P, Eberhardt RY, Eddy SR, Mistry J, Mitchell AL, Potter SC, Punta M, Qureshi M, Sangrador-Vegas A, Salazar GA, Tate J, Bateman A. 2016. The Pfam protein families database: towards a more sustainable future. *Nucleic Acids Research* **44**:D279-D285. doi: 10.1093/nar/gkv1344, PMID: 26673716

Ganguli S, Wang H, Matsumura P, Volz K. 1995. Uncoupled phosphorylation and activation in bacterial chemotaxis. the 2.1-A structure of a threonine to isoleucine mutant at position 87 of CheY. *Journal of Biological Chemistry* **270**:17386-17393. doi: 10.1074/jbc.270.29.17386, PMID: 7615544

Gao R, Stock AM. 2009. Biological insights from structures of two-component proteins. *Annual Review of Microbiology* **63**:133-154. doi: 10.1146/annurev.micro.091208.073214, PMID: 19575571

Gao R, Stock AM. 2010. Molecular strategies for phosphorylation-mediated regulation of response regulator activity. *Current Opinion in Microbiology* **13**:160-167. doi: 10.1016/j.mib.2009.12.009, PMID: 20080056

Gautam US, Sikri K, Tyagi JS. 2011. The residue threonine 82 of DevR (DosR) is essential for DevR activation and function in *Mycobacterium tuberculosis* despite its atypical location. *Journal of Bacteriology* **193**:4849-4858. doi: 10.1128/JB.05051-11, PMID: 21764934

Goulian M. 2010. Two-component signaling circuit structure and properties. *Current Opinion in Microbiology* **13**: 184-189. doi: 10.1016/j.mib.2010.01.009, PMID: 20149717

Homeyer N, Horn AH, Lanig H, Sticht H. 2006. AMBER force-field parameters for phosphorylated amino acids in different protonation states: phosphoserine, phosphothreonine, phosphotyrosine, and phosphohistidine. *Journal of Molecular Modeling* **12**:281-289. doi: 10.1007/s00894-005-0028-4, PMID: 16240095

Hsing W, Silhavy TJ. 1997. Function of conserved histidine-243 in phosphatase activity of EnvZ, the sensor for porin osmoregulation in *Escherichia coli*. *Journal of Bacteriology* **179**:3729-3735. doi: 10.1128/jb.179.11.3729-3735.1997, PMID: 9171423

Huynh TN, Noriega CE, Stewart V. 2010. Conserved mechanism for sensor phosphatase control of two-component signaling revealed in the nitrate sensor NarX. *PNAS* **107**:21140-21145. doi: 10.1073/pnas.1013081107, PMID: 21078995

Imelio J, Larrieux N, Mechaly A, Trajtenberg F, Buschiazzo A. 2017. Snapshots of the Signaling Complex DesK: DesR in Different Functional States Using Rational Mutagenesis and X-ray Crystallography. *BIO-PROTOCOL* **7**: e2510. doi: 10.21769/BioProtoc.2510

Immormino RM, Silversmith RE, Bourret RB. 2016. A variable active site residue influences the kinetics of response regulator phosphorylation and dephosphorylation. *Biochemistry* **55**:5595-5609. doi: 10.1021/acs.biochem.6b00645, PMID: 27589219

Izaguirre Jesús A., Catarello DP, Wozniak JM, Skeel RD. 2001. Langevin stabilization of molecular dynamics. *The Journal of Chemical Physics* **114**:2090-2098. doi: 10.1063/1.1332996

Janiak-Spens F, Cook PF, West AH. 2005. Kinetic analysis of YPD1-dependent phosphotransfer reactions in the yeast osmoregulatory phosphorelay system. *Biochemistry* **44**:377-386. doi: 10.1021/bi048433s, PMID: 1562 8880

Jorgensen WL, Chandrasekhar J, Madura JD, Impey RW, Klein ML. 1983. Comparison of simple potential functions for simulating liquid water. *The Journal of Chemical Physics* **79**:926-935. doi: 10.1063/1.445869

Joung IS, Cheatham TE. 2008. Determination of alkali and halide monovalent ion parameters for use in explicitly solvated biomolecular simulations. *The Journal of Physical Chemistry B* **112**:9020-9041. doi: 10.1021/jp8001614, PMID: 18593145

Konarev PV, Petoukhov MV, Volkov VV, Svergun DI. 2006. ATSAS 2.1, a program package for small-angle scattering data analysis. *Journal of Applied Crystallography* **39**:277-286. doi: 10.1107/S0021889806004699

Lassila JK, Zalatan JG, Herschlag D. 2011. Biological phosphoryl-transfer reactions: understanding mechanism and catalysis. *Annual Review of Biochemistry* **80**:669-702. doi: 10.1146/annurev-biochem-060409-092741, PMID: 21513457

Lawrence MC, Colman PM. 1993. Shape complementarity at protein/protein interfaces. *Journal of Molecular Biology* **234**:946-950. doi: 10.1007/jmbi.1993.1648, PMID: 8263940
Li P, Roberts BP, Chakravorty DK, Merz KM. 2013. Rational design of particle mesh Ewald compatible Lennard-Jones parameters for +2 metal cations in explicit solvent. *Journal of Chemical Theory and Computation* **9**: 2733-2748. doi: 10.1021/ct400146w, PMID: 23914143

Li W, Godzik A. 2006. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. *Bioinformatics* **22**:1658-1659. doi: 10.1093/bioinformatics/btl158, PMID: 16731699

Loewenthal R, Sancho J, Fersht AR. 1992. Histidine-aromatic interactions in barnase : elevation of histidine pKa and contribution to protein stability. *Journal of Molecular Biology* **224**:759-770. doi: 10.1016/0022-2836(92) 90560-7, PMID: 1569555

Lukat GS, Stock AM, Stock JB. 1990. Divalent metal ion binding to the CheY protein and its significance to phosphotransfer in bacterial chemotaxis. *Biochemistry* **29**:5436-5442. doi: 10.1021/bi00475a004, PMID: 2201404

Maier JA, Martinez C, Kasavajhala K, Wickstrom L, Hauser KE, Simmerling C. 2015. ff14SB: Improving the accuracy of protein side chain and backbone parameters from ff99SB. *Journal of Chemical Theory and Computation* **11**:3696-3713. doi: 10.1021/acs.jctc.5b00255, PMID: 26574453

Marina A, Waldburger CD, Hendrickson WA. 2005. Structure of the entire cytoplasmic portion of a sensor histidine-kinase protein. *The EMBO Journal* **24**:4247-4259. doi: 10.1038/sj.emboj.7600886, PMID: 16319927

McCoy AJ, Grosse-Kunstleve RW, Adams PD, Winn MD, Storoni LC, Read RJ. 2007. Phaser crystallographic software. *Journal of Applied Crystallography* **40**:658-674. doi: 10.1107/S0021889807021206, PMID: 19461840

Meagher KL, Redman LT, Carlson HA. 2003. Development of polyphosphate parameters for use with the AMBER force field. *Journal of Computational Chemistry* **24**:1016-1025. doi: 10.1002/jcc.10262, PMID: 1275 9902

Mechaly AE, Sassoon N, Betton JM, Alzari PM. 2014. Segmental helical motions and dynamical asymmetry modulate histidine kinase autophosphorylation. *PLoS Biology* **12**:e1001776. doi: 10.1371/journal.pbio.1001776, PMID: 24492262

Mildvan AS. 1997. Mechanisms of signaling and related enzymes. *Proteins: Structure, Function, and Genetics* **29**: 401-416. doi: 10.1002/(SICI)1097-0134(199712)29:4<401::AID-PROT1>3.0.CO;2-B, PMID: 9408938

Miyamoto S, Kollman PA. 1992. Settle: An analytical version of the SHAKE and RATTLE algorithm for rigid water models. *Journal of Computational Chemistry* **13**:952-962. doi: 10.1002/jcc.540130805

Morcos F, Pagnani A, Lunt B, Bertolino A, Marks DS, Sander C, Zecchina R, Onuchic JN, Hwa T, Weigt M. 2011. Direct-coupling analysis of residue coevolution captures native contacts across many protein families. *PNAS* **108**:E1293-1301. doi: 10.1073/pnas.1111471108, PMID: 22106262

Morin A, Eisenbraun B, Key J, Sanschagrin PC, Timony MA, Ottaviano M, Sliz P. 2013. Collaboration gets the most out of software. *eLife* **2**:e01456. doi: 10.7554/eLife.01456, PMID: 24040512

Page SC, Immormino RM, Miller TH, Bourret RB. 2016. Experimental analysis of functional variation within protein families: Receiver domain autodephosphorylation kinetics. *Journal of Bacteriology* **198**:2483-2493. doi: 10.1128/JB.00853-15, PMID: 27381915

Page SC, Silversmith RE, Collins EJ, Bourret RB. 2015. Imidazole as a small molecule analogue in Two-Component signal transduction. *Biochemistry* **54**:7248-7260. doi: 10.1021/acs.biochem.5b01082, PMID: 2656 9142

Pastor RW, Brooks BR, Szabo A. 1988. An analysis of the accuracy of Langevin and molecular dynamics algorithms. *Molecular Physics* **65**:1409-1419. doi: 10.1080/00268978800101881

Pazy Y, Motaleb MA, Guarnieri MT, Charon NW, Zhao R, Silversmith RE. 2010. Identical phosphatase mechanisms achieved through distinct modes of binding phosphoprotein substrate. *PNAS* **107**:1924-1929. doi: 10.1073/pnas.0911185107

Podgornaia AI, Casino P, Marina A, Laub MT. 2013. Structural basis of a rationally rewired protein-protein interface critical to bacterial signaling. *Structure* **21**:1636-1647. doi: 10.1016/j.str.2013.07.005, PMID: 23954504

Podgornaia AI, Laub MT. 2015. Protein evolution. Pervasive degeneracy and epistasis in a protein-protein interface. *Science* **347**:673-677. doi: 10.1126/science.1257360, PMID: 25657251

Porter SL, Roberts MA, Manning CS, Armitage JP. 2008. A bifunctional kinase-phosphatase in bacterial chemotaxis. *PNAS* **105**:18531-18536. doi: 10.1073/pnas.0808010105, PMID: 19020080

Potter CA, Ward A, Laguri C, Williamson MP, Henderson PJ, Phillips-Jones MK. 2002. Expression, purification and characterisation of full-length histidine protein kinase RegB from *Rhodobacter sphaeroides*. *Journal of Molecular Biology* **320**:201-213. doi: 10.1016/S0022-2836(02)00424-2, PMID: 12079379

Quezada CM, Hamel DJ, Gradinaru C, Bilwes AM, Dahlquist FW, Crane BR, Simon MI. 2005. Structural and chemical requirements for histidine phosphorylation by the chemotaxis kinase CheA. *Journal of Biological Chemistry* **280**:30581-30585. doi: 10.1074/jbc.M505316200, PMID: 15994328

Rambo RP, Tainer JA. 2013. Accurate assessment of mass, models and resolution by small-angle scattering. *Nature* **496**:477-481. doi: 10.1038/nature12070, PMID: 23619693

Ryckaert J-P, Ciccotti G, Berendsen HJC. 1977. Numerical integration of the cartesian equations of motion of a system with constraints: molecular dynamics of n-alkanes. *Journal of Computational Physics* **23**:327-341. doi: 10.1016/0021-9991(77)90098-5

Saita E, Abriata LA, Tsai YT, Trajtenberg F, Lemmin T, Buschiazzo A, Dal Peraro M, de Mendoza D, Albanesi D. 2015. A coiled coil switch mediates cold sensing by the thermosensory protein desk. *Molecular Microbiology* **98**:258-271. doi: 10.1111/mmi.13118, PMID: 26172072

Salomon-Ferrer R, Case DA, Walker RC. 2013. An overview of the amber biomolecular simulation package. *Wiley Interdisciplinary Reviews: Computational Molecular Science* **3**:198-210. doi: 10.1002/wcms.1121
Schneider CA, Rasband WS, Eliceiri KW. 2012. NIH image to ImageJ: 25 years of image analysis. *Nature Methods* **9**:671–675. doi: 10.1038/nmeth.2089, PMID: 22930834

Schramke H, Tostevin F, Heermann R, Gerland U, Jung K. 2016. A Dual-Sensing receptor confers robust cellular homeostasis. *Cell Reports* **16**:1–9. doi: 10.1016/j.celrep.2016.05.081, PMID: 27320909

Schrodinger, LLC. 2015. The PyMOL Molecular Graphics System. Version 1.8.

Schultz JE, Natarajan J. 2013. Regulated unfolding: a basic principle of intraprotein signaling in modular proteins. *Trends in Biochemical Sciences* **38**:538–545. doi: 10.1016/j.tibs.2013.08.005, PMID: 24051046

Shi L, Liu W, Hulett FM. 1999. Decay of activated *Bacillus subtilis* pho response regulator, PhoP approximately P, involves the PhoR approximately P intermediate. *Biochemistry* **38**:10119–10125. doi: 10.1021/bi990658t, PMID: 10433720

Skarphol K, Waukau J, Forst SA. 1997. Role of His243 in the phosphatase activity of EnvZ in *Escherichia coli*. *Journal of Bacteriology* **179**:1413–1416. doi: 10.1128/jb.179.4.1413-1416.1997, PMID: 9023231

Skerker JM, Perchuk BS, Siryaporn A, Lubin EA, Ashenberg O, Goulian M, Laub MT. 2008. Rewiring the specificity of two-component signal transduction systems. *Cell* **133**:1043–1054. doi: 10.1016/j.cell.2008.04.040, PMID: 18555780

Skerker JM, Prasol MS, Perchuk BS, Biondi EG, Laub MT. 2005. Two-component signal transduction pathways regulating growth and cell cycle progression in a bacterium: a system-level analysis. *PLoS Biology* **3**:e334. doi: 10.1371/journal.pbio.0030334, PMID: 16176121

Thomas SA, Immormino RM, Bourret RB, Silversmith RE. 2013. Nonconserved active site residues modulate CheY autophosphorylation kinetics and phosphodonor preference. *Biochemistry* **52**:2262–2273. doi: 10.1021/ bi301654m, PMID: 23458124

Trajtenberg F, Albanesi D, Ruétalo N, Botti H, Mechaly AE, Nieves M, Aguilar PS, Cybulski L, Larrieux N, de Mendoza D, Buschiazzo A. 2014. Allosteric activation of bacterial response regulators: the role of the cognate histidine kinase beyond phosphorylation. *mBio* **5**:e02105. doi: 10.1128/mBio.02105-14, PMID: 25406381

Unger T, Jacobovitch Y, Dantes A, Bernheim R, Peleg Y. 2010. Applications of the restriction free (RF) cloning procedure for molecular manipulations and protein expression. *Journal of Structural Biology* **172**:34–44. doi: 10.1016/j.jsb.2010.06.016, PMID: 20600952

Varughese KI, Tsigelny I, Zhao H. 2006. The crystal structure of beryllofluoride Spo0F in complex with the phosphotransferase Spo0B represents a phosphotransfer pretransition state. *Journal of Bacteriology* **188**:4970–4977. doi: 10.1128/JB.00160-06, PMID: 16788205

Vonrhein C, Flensburg C, Keller P, Sharff A, Smart O, Paciorek W, Womack T, Bricogne G. 2011. Data processing and analysis with the autoPROC toolbox. *Acta Crystallographica Section D Biological Crystallography* **67**:293–302. doi: 10.1107/S0907444911007773, PMID: 21460447

Weigt M, White RA, Szurmant H, Hoch JA, Hwa T. 2009. Identification of direct residue contacts in protein-protein interaction by message passing. *PNAS* **106**:67–72. doi: 10.1073/pnas.0805923106, PMID: 19116270

Willett JW, Herrou J, Briegel A, Rotskoff G, Crosson S. 2015. Structural asymmetry in a conserved signaling system that regulates division, replication, and virulence of an intracellular pathogen. *PNAS* **112**:E3709–E3718. doi: 10.1073/pnas.1503118112, PMID: 26124143

Willett JW, Kirby JR. 2012. Genetic and biochemical dissection of a HisKA domain identifies residues required exclusively for kinase and phosphatase activities. *PLoS Genetics* **8**:e1003084. doi: 10.1371/journal.pgen.1003084, PMID: 23226719

Wu X, Brooks BR. 2003. Self-guided Langevin dynamics simulation method. *Chemical Physics Letters* **381**:512–518. doi: 10.1016/j.cplett.2003.10.013

Yamada S, Sugimoto H, Kobayashi M, Ohno A, Nakamura H, Shiro Y. 2009. Structure of PAS-linked histidine kinase and the response regulator complex. *Structure* **17**:1333–1344. doi: 10.1016/j.str.2009.07.016, PMID: 19836334

Yoshida T, Cai S, Inouye M. 2002. Interaction of EnvZ, a sensory histidine kinase, with phosphorylated OmpR, the cognate response regulator. *Molecular Microbiology* **46**:1283–1294. doi: 10.1046/j.1365-2958.2002.03240.x, PMID: 12453215

Zapf J, Sen U, Madhusudan, Hoch JA, Varughese KI. 2000. A transient interaction between two phosphorelay proteins trapped in a crystal lattice reveals the mechanism of molecular recognition and phosphotransfer in signal transduction. *Structure* **8**:851–862. doi: 10.1016/S0969-2126(00)00174-X, PMID: 10997904

Zhao R, Collins EJ, Bourret RB, Silversmith RE. 2002. Structure and catalytic mechanism of the *E. coli* chemotaxis phosphatase CheZ. *Nature Structural Biology* **9**:570–575. doi: 10.1038/nsb816, PMID: 12080332

Zhao X, Copeland DM, Soares AS, West AH. 2008. Crystal structure of a complex between the phosphorelay protein YPD1 and the response regulator domain of SLN1 bound to a phosphoryl analog. *Journal of Molecular Biology* **375**:1141–1151. doi: 10.1016/j.jmb.2007.11.045, PMID: 18076904

Zhu Y, Qin L, Yoshida T, Inouye M. 2000. Phosphatase activity of histidine kinase EnvZ without kinase catalytic domain. *PNAS* **97**:7808–7813. doi: 10.1073/pnas.97.14.7808, PMID: 10884412

Zschiedrich CP, Keidel V, Szurmant H. 2016. Molecular mechanisms of two-component signal transduction. *Journal of Molecular Biology* **428**:3752–3775. doi: 10.1016/j.jmb.2016.08.003, PMID: 27519796