research papers
Structural insights into the histidine-containing phosphotransfer protein and receiver domain of sensor histidine kinase suggest a complex model in the two-component regulatory system in Pseudomonas aeruginosa
aDepartment of Biotechnology and Bioindustry Sciences, National Cheng Kung University, Tainan City 701, Taiwan, bLife Science Group, Scientific Research Division, National Synchrotron Radiation Research Center, Hsinchu 30076, Taiwan, cInstitute of Molecular Medicine, National Tsing Hua University, Hsinchu 30013, Taiwan, dInstitute of Bioinformatics and Structural Biology, National Tsing Hua University, Hsinchu 30013, Taiwan, eDepartment of Physics, National Tsing Hua University, Hsinchu 30013, Taiwan, and fDepartment of Biological Science and Technology, National Chiao Tung University, Hsinchu 30010, Taiwan
*Correspondence e-mail: cjchen@nsrrc.org.tw
In Pseudomonas aeruginosa, an important opportunistic pathogen that causes numerous acute and chronic infections, the hybrid two-component system (TCS) regulates the swarming ability and biofilm formation with a multistep phosphorelay, and consists of hybrid-sensor histidine kinase (HK), histidine-containing phosphotransfer protein (Hpt) and response regulator (RR). In this work, two crystal structures of HptB and the receiver domain of HK PA1611 (PA1611REC) of P. aeruginosa have been determined in order to elucidate their interactions for the transfer of the phosphoryl group. The structure of HptB folds into an elongated four-helix bundle – helices α2, α3, α4 and α5, covered by the short N-terminal helix α1. The imidazole side chain of the conserved active-site histidine residue His57, located near the middle of helix α3, protrudes from the bundle and is exposed to solvent. The structure of PA1611REC possesses a conventional (β/α)5 topology with five-stranded parallel β-sheets folded in the central region, surrounded by five α-helices. The divalent Mg2+ ion is located in the negatively charged active-site cleft and interacts with Asp522, Asp565 and Arg567. The HptB–PA1611REC complex is further modeled to analyze the binding surface and interactions between the two proteins. The model shows a shape complementarity between the convex surface of PA1611REC and the kidney-shaped HptB with fewer residues and a different network involved in interactions compared with other TCS complexes, such as SLN1-R1/YPD1 from Saccharomyces cerevisiae and AHK5RD/AHP1 from Arabidopsis thaliana. These structural results provide a better understanding of the TCS in P. aeruginosa and could potentially lead to the discovery of a new treatment for infection.
Keywords: two-component regulatory systems; sensor histidine kinase; histidine-containing phosphotransfer proteins; Pseudomonas aeruginosa.
PDB references: PA1611REC-Mg2+, 7c1j; PA1611REC-Ca2+, 7cfw; WT-HptB, 7c1i
1. Introduction
In aiding the adaptation of cells to environmental changes, two-component regulatory systems (TCS) are widely distributed in prokaryotes, whereas few are identified in lower eukaryotic organisms and plants. The basic mechanism of the TCS signaling transduction is the transfer of a phosphoryl group, which serves as a signaling molecule. The system is composed of a sensor histidine kinase (HK), which is autophosphorylated at a conserved histidine residue of its cytoplasmic transmitter domain after sensing outside signals with the extracellular input domain; and a response regulator (RR), which accepts the phosphoryl group at the conserved aspartate residue of the receiver domain (REC) from HK or a histidine-containing phosphotransfer protein (Hpt) (Stock et al., 1989; 2000). The RR protein initiates the response by binding to a specific target protein, such as a transcriptional factor or DNA element (Bell et al., 2010; Bauer et al., 2013). The chemistry of the phosphoramide and acylphosphate linkage is involved in this mechanism (Thomason & Kay, 2000). Biofilm formation and the implementation of drug resistance are typically regulated by TCS, which makes TCS an attractive target for research and drug intervention (Anjali et al., 2019).
The sensor–regulator protein pairs in a TCS can be classified into three major types based on the organization of functional domains of HK: classic, unorthodox and hybrid systems (Rodrigue et al., 2000; Lin et al., 2006). In the classic system, a phosphate is transferred directly from phosphorylated HK to RR (Skerker et al., 2008; Capra et al., 2010). In the unorthodox system, an extra receiver domain (REC) and a histidine-containing phosphotransfer (Hpt) domain are sequentially connected to the C-terminus of the HK transmitter domain. After the extracellular input domain of HK senses signals from the environment, the transmitter domain autophosphorylates; the phosphate is then relayed to the receiver domain and the Hpt domain of HK. Subsequently, the Hpt domain transphosphorylates the conserved aspartate residue at the receiver domain of RR (Posas et al., 1996; Thomason & Kay, 2000; Zhao et al., 2008; Kaserer et al., 2010). The hybrid system is similar to the unorthodox system, but, instead of an Hpt domain attached to HK, an independent Hpt protein transfers the phosphoryl group to the downstream RR (Rodrigue et al., 2000).
Pseudomonas aeruginosa (P. aeruginosa) is a versatile Gram-negative bacterium grown in soil, coastal marine habitats, and on plant and animal tissues (Stover et al., 2000). It is an elongated rod-shaped bacterium with only one-way flagellum mobility. Previous studies have shown that the metabolic versatility, intrinsic and acquired antibiotic resistance, biofilm formation and production of multiple virulence factors make P. aeruginosa a formidable pathogen causing numerous acute and chronic infections (Balasubramanian et al., 2013). It accounts for 10–20% of all hospital-acquired infections and is the leading cause of chronic pulmonary infections and mortality in cystic fibrosis (CF) patients and burn victims (Stover et al., 2000; National Nosocomial Infections Surveillance System, 2004; Lyczak et al., 2002). In P. aeruginosa, the Hpt-mediated hybrid TCS serves as a basic stimulus–response coupling mechanism to decrease the environmental harm and to ensure survival by activating downstream responses, including antibiotic susceptibility, swarming activity and biofilm formation (Stock et al., 2000; Wolanin et al., 2002). Many TCS components have been described for their key roles during the infection processes of P. aeruginosa (Rodrigue et al., 2000).
The three orphan Hpt proteins in P. aeruginosa are HptA, HptB and HptC. HptB, the target protein in this work, is closely related to biofilm formation, which could increase drug resistance and protect the cell from pinocytosis by the host. According to previous studies, after the sensor proteins PA1611, PA1976 and PA2824 are autophosphorylated, the phosphoryl group is transferred to HptB and subsequently to the downstream RR protein PA3346 as illustrated in Fig. 1. The results showed that the HptB-mediated multistep phosphorelay plays a central role in the infection process of P. aeruginosa (Lin et al., 2006; Hsu et al., 2008). Moreover, the previous structures of Hpt proteins/domains in complex with receiver domains of HK or RR, such as SLN1-R1/YPD1 of yeast (PDB entry 2r25) and AHK5RD/AHP1 of Arabidopsis thaliana (PDB entry 4euk; Bauer et al., 2013), have shown that the Hpt protein/domain would interact with the receiver domain. To elucidate the HptB-mediated TCS pathway of P. aeruginosa in depth, we report two crystal structures of HptB and the C-terminal receiver domain of HK PA1611 (PA1611REC). The structural details allow us to model the complex of HptB and PA1611REC, which reveals the possible binding conformation and the residues involved in the interaction between the two proteins. The result would further clarify the mechanism in the HptB-mediated signal transduction pathway of P. aeruginosa and potentially lead to the discovery of a new treatment for P. aeruginosa infection, such as novel small molecules which could potentially interfere with the interaction between HptB and the corresponding receiver domains.
2. Materials and methods
2.1. Cloning and expression of target proteins
Native HptB was cloned into pET23a (Novagen) using NdeI and XhoI restriction sites. For the mutant HptB, the only cysteine (C75) was mutated to alanine by . E. coli strain DH5α was used for the plasmid preparation and recombinant DNA manipulation.
(PCR) overlap extension. The C-terminal fragment (a.a. 507–651) including the receiver domain of HK PA1611 (PA1611REC) was amplified by PCR and cloned into pET28a (Novagen) using NheI and EcoRI restriction sites with T4 DNA ligase (Invitrogen). The primers used are listed in Table 1
|
After DNA sequencing confirmed the correctness of the insert, the expression constructs of both native HptB and PA1611REC were transformed into E. coli strain BL21(DE3) for the expression of recombinant proteins. The bacteria were cultured with shaking in Luria–Bertani media supplemented with specific antibiotics (100 µg ml−1) until the at 600 nm (OD600) reached 0.6. Overexpression of the target protein was induced on adding isopropyl-β-D-thiogalactoside (IPTG) to the final concentrations of 0.5 mM and 0.1 mM for native HptB and PA1611REC, respectively. The culture was grown overnight at 20°C with shaking and harvested using centrifugation at 8000g for 25 min at 4°C.
For the production of selenomethionine-labeled (SeMet) mutant HptB, E. coli strain BL21-Gold (DE3) harboring the expression vector served as the host cell. The culture was grown in M9 medium (6 mg l−1 Na2HPO4, 3 mg l−1 KH2PO4, 0.5 mg l−1 NaCl, 1 mg l−1 NH4Cl, 1 mM MgSO4, 0.1 mM CaCl2, 1 mg l−1 thiamine-HCl, 0.2% glucose) with L-lysine, L-threonine, L-phenylalanine (120 mg l−1), L-leucine, L-isoleucine, L-valine, L-selenomethionine (60 mg l−1) and ampicillin (100 µg ml−1) until the at 600 nm (OD600) reached 0.6. The following steps for the expression and purification of the SeMet-labeled protein were the same as the native HptB protein.
2.2. Purification of HptB protein
The cell pellet from a 1 l culture was suspended in a lysis buffer (35 ml) consisting of Tris–HCl (50 mM, pH 8.0) and subjected to cell disruption by ultrasonication using a pulsation cycle 2 s on and 3 s off with a total duration of 20 min at 40% energy on ice. The soluble protein extract was collected by centrifugation at 12 000g for 30 min at 4°C and passed through an open column packed with resin (Q Sepharose Fast Flow, 60 ml, GE Healthcare Bio-Sciences AB, Uppsala, Sweden), which was equilibrated with binding buffer (90 ml) consisting of Tris–HCl (50 mM, pH 8.0). The column was washed with binding buffer (90 ml) containing (NH4)2SO4 (40 mM). The target protein was eluted with portions (90 ml) of binding buffer containing (NH4)2SO4 at stepwise-increasing concentrations of 60 and 100 mM. The eluted fractions were concentrated by the centricon of a 3 kDa filter unit (Amicon Ultra, Merck KGaA, Darmstadt, Germany) and then loaded onto a column (HiLoad 16/60 Superdex 75 PG, GE Healthcare Bio-Sciences AB) connected to a (FPLC) system (ÄKTA, GE Healthcare) and equilibrated with a buffer composed of NaCl (300 mM) and Tris–HCl (20 mM, pH 8.0). The step yielded homogenous protein fractions.
2.3. Purification of PA1611REC protein
The procedures of cell disruption and collection of the soluble protein extract were the same as those for the purification of HptB except for the lysis buffer, which was composed of NaCl (500 mM) and Tris–HCl (50 mM, pH 8.0). The extract was then passed through an open column packed with resin (10 ml Ni-Sepharose 6 Fast Flow, GE Healthcare Bio-Sciences AB, Uppsala, Sweden), which was equilibrated with a binding buffer (50 ml) consisting of NaCl (500 mM) and Tris–HCl (50 mM, pH 8.0). The column was washed with a binding buffer (50 ml) containing imidazole (20 mM). The target protein was eluted with portions (50 ml) of binding buffer containing imidazole at stepwise-increasing concentrations of 80, 100, 300 and 500 mM.
2.4. Crystallization
After the eluted fractions were collected and verified with SDS–PAGE (15%) followed by staining (Coomassie Brilliant Blue R-250), the purified HptB protein was dialyzed against Tris–HCl (50 mM, pH 7.5) and concentrated to ∼8 mg ml−1 with a Centricon 3 kDa filter unit (Amicon Ultra, Merck KGaA, Darmstadt, Germany), whereas the purified PA1611REC protein was dialyzed against distilled water and concentrated to ∼10 mg ml−1 by the same method. Initial screening of crystallization conditions was performed by the hanging-drop vapor-diffusion method in 96-well plates with a liquid-handling robot (Mosquito, TTP Labtech). The hanging drops formed on mixing protein solution (0.1 µl) with reservoir solution (0.1 µl) and equilibrated against the reservoir solution (100 µl). To optimize the protein crystals, the hanging drops were made on mixing a protein solution (1 µl) with a reservoir solution (1 µl) and equilibrated against a reservoir solution (120 µl) in VDX48 plates (Hampton Research). All plates were placed in an incubator at 18°C.
Crystals formed under several conditions. The crystal of SeMet-labeled mutant HptB, grown in a reservoir solution containing PEG20000 [12%(w/v)] and MES monohydrate (0.1 M, pH 6.5), was picked for data collection. The structure of native HptB determined was from the crystals grown in a reservoir solution containing PEG6000 [10%(w/v)] and bicine/sodium hydroxide (100 mM, pH 9.0). For PA1611REC, the structure determined was from crystals grown in a reservoir solution containing PEG4000 [15%(w/v)], MgCl2 (0.2 M) and Tris–HCl (100 mM, pH 8.5).
2.5. Data collection and processing
The crystals were transferred from the hanging drops into a cryoprotectant solution (1 µl) made from adding glycerol [20%(v/v)] into the reservoir solution for a few seconds. The crystals were then picked with synthetic nylon loops (Hampton Research, Aliso Viejo, CA, USA) and flash-cooled in liquid nitrogen. The X-ray diffraction datasets were collected on beamlines TPS 05A, BL13B1 and BL15A1 at the National Synchrotron Radiation Research Center (NSRRC), Taiwan, and BL44XU and BL12B2 at SPring-8, Japan. For SeMet-labeled HptB, a multi-wavelength anomalous-dispersion (Se-MAD) experiment was performed. All datasets were indexed, integrated and scaled using the HKL2000 program suite (Otwinowski & Minor, 1997). The data collection statistics are listed in Table 2.
‡Rsym = , where Ii is the ith measurement and 〈I(h)〉 is the of all measurements of I(h). §Rwork = , where Fo and Fc are the observed and calculated amplitudes of reflection h, respectively. ¶Rfree is the same as Rwork, but calculated with 5% of randomly chosen reflections omitted from |
2.6. and refinement
The positions of nine selenium atoms of three SeMet-labeled mutant HptB–C75A molecules in the I4122 were automatically determined and used to calculate the multi-wavelength anomalous-diffraction phases with SOLVE (Terwilliger & Berendzen, 1999). The phases of wild-type HptB in P43212 and PA1611REC in P3121 were determined with the method using Phaser MR in the CCP4 program suite (Winn et al., 2011), with the structures of ab inito determined mutant HptB–C75A and the receiver domain from yeast histidine kinase SLN1 (PDB entry 2r25; B-chain) as the search models, respectively. Structural was performed using alternating rounds of model building with the program Coot (Emsley et al., 2010) and with Refmac5 in the CCP4 program suite (Murshudov et al., 2011) to improve the value of the R factor and Rfree. The are listed in Table 2. The secondary structure was assigned with the DSSP web server (Kabsch & Sander, 1983; Touw et al., 2015), and the structure figures were prepared using PyMOL (https://www.pymol.org/).
of2.7. Modeling of the HptB–PA1611REC complex
To establish the complex model of HptB and PA1611REC, docking was performed using the crystal structures of two individual proteins with the HADDOCK web server (version 2.2; van Zundert et al., 2016) in the easy interface mode. His57 of HptB and Asp565 of PA1611REC were considered the active residues because these two residues directly participate in the transfer of the phosphoryl group according to previous studies, which confirmed this transfer between PA1611 and HptB (Lin et al., 2006; Hsu et al., 2008). The passive residues, which potentially make contact in the complex, were automatically defined in the easy interface mode that included all residues within a radius of 6.5 Å from the active residue and on the surface with relative surface accessibility of either a main chain or a side chain above 15% as determined with NACCESS (Hubbard & Thornton, 1993). The divalent metal ions (Mg2+) in the active-site cleft of PA1611REC were also applied to the HADDOCK web server. The complex models were clustered according to the pairwise backbone root-mean-square deviation (RMSD) at the interface; the clusters were sorted according to a HADDOCK score calculated from energies of van der Waals, electrostatic interaction, restraint violation and buried surface area. The HADDOCK web server clustered 115 structures into eight clusters, which represents 57.5% of the water-refined models generated by HADDOCK for the HptB–PA1611REC complex with Mg2+. The cluster with the lowest HADDOCK score is the most likely solution; the first model in this cluster was chosen for the analysis of the interface with the PDBePISA web server (Krissinel & Henrick, 2007). The score of cluster 1 is −63.7 ± 2.6, whereas other clusters range from −48 to −7.6, indicating that the highest score is a unique solution.
3. Results
3.1. Overall structure of HptB
HptB consists of 116 amino acids with a total molecular mass of ∼13.2 kDa and shares a low-sequence identity with other Hpt proteins or domains [Fig. 2(a)]. Initially, we found crystallizing native HptB to be difficult. We suspected that Cys75 would introduce an intermolecular disulfide bond and interrupt the regular orientation of the protein in solution. We thus prepared SeMet-labeled mutant HptB–C75A for crystallization and phase determination. The crystals of SeMet-labeled mutant HptB–C75A belong to the tetragonal I4122 with three protein molecules in one Then, we successfully crystallized the native HptB in another tetragonal P43212, with six protein molecules in one The two pairs of trimeric proteins related to the non-crystallographic twofold symmetry exhibit similar structures [Fig. 3(a)], but the protein is monomeric in solution, as confirmed by The initial phases and structure of SeMet-labeled mutant HptB–C75A were first determined with the Se-MAD method at ∼2 Å. The of native HptB was subsequently determined with the mutant HptB structure as an initial model and refined to high resolution (1.58 Å). The structural differences between mutant and native HptB are small, with an RMSD of 0.15 Å. The structures of six native HptBs in one are essentially the same, with a small RMSD of 0.16 Å overall for the Cα atoms [Fig. 3(b)]. The final atomic model of native HptB comprises 115 residues, except for the first methionine. The entire main chain and most side chains are well defined with clear electron density.
The structure of HptB comprises five α-helices and folds into an elongated α-helical bundle with an up-and-down topology. The helices α2 (Tyr22–Gln45), α3 (Ala49–Asn65), α4 (Val69–Arg83) and α5 (Ala90–Arg112) are twisted around the central axis and form an antiparallel four-helix bundle. Residues Ala46–Asp48, Met66–Ala68 and Arg84–Arg89 form short turns between helices α2–α3, α3–α4 and α4–α5, respectively. The four-helix bundle is covered with the short N-terminal helix α1 (Asp8–Met18) which protrudes into the solvent environment and connects to helix α2 with a 310 helix (Glu19–Glu21) [Fig. 3(c)]. All these inter-helical regions are well ordered with some residues contributing to the stabilization of the helix-bundle through hydrogen bonds (Fig. 4).
The four-helix bundle is stabilized predominately by inter-helix hydrophobic contacts contributed by the side chains of residues Ala3, Pro4 and Leu6 at the N-terminal loop; Leu11, Leu14, Val17 and Met18 at helix α1; Tyr22, Pro23, Val24, Leu25, Leu26, Phe29, Val30, Leu37, Leu40, Ala43 and Leu44 at helix α2; Ala49, Ala51, Leu52, Ala56 and Phe59 at helix α3; Met66 and Ala68 at the turn between helices α3 and α4; Leu70, Leu71, Tyr74, Leu78 and Ala82 at helix α4; Leu87 at the turn between helices α4 and α5; Ala90, Pro91, Ile94, Met97, Phe101, Ile103, Val104, Ile106, Leu107 and Phe108 at helix α5; and Tyr115 at the C-terminal loop. Helix α1 sits on top of the bundle to form a cap and provides additional stabilization [Fig. 3(c)]. Aside from directly contributing the hydrophobic side chains to the inter-helix hydrophobic contact, helix α1 helps to stabilize the bundle with the side chains of non-hydrophobic residues, which are highly accessible on the surface. The non-hydrophobic residues on helix α1 shield the hydrophobic regions of helices α2 and α5, which are not covered by helices α3 and α4, and the hydrophobic turn connecting helices α3 and α4 from exposure to solvent.
3.2. The conserved Hpt fold and structural features in the active site
The protein family of Hpt proteins and domains could be further divided into enzymatically active and inactive orthologs. The active orthologs have a conserved histidine residue that serves as the phosphorelay active site (Ruszkowski et al., 2013). In HptB, this conserved histidine residue, His57, is located near the middle of helix α3. The imidazole side chain of His57 protrudes from the bundle and is exposed to solvent [Figs. 3(c) and 3(d)]. The surrounding residues are featured so as to maximize the accessibility of His57. The hydrogen-bonding linkages between Lys60 and Glu79 prohibit their side chains from shielding the imidazole ring of the active-site histidine. High solvent accessibility of the imidazole ring is achieved through the nearby small side-chain volume, which results from the highly conserved glycine residue (Gly61 in HptB) at position +4 from the active-site histidine (Fig. 5).
Some conserved residues are related to the stabilization of the structural integrity. Aside from several residues contributed to inter-helix hydrophobic contacts, the residues at positions −5 (Leu52 in HptB), −9 (Asp48), +6 (Ser63), +9 (Met66), +10 (Gly67) and +11 (Ala68) from the active-site histidine are conserved among Hpt proteins/domains with the structures determined listed in Fig. 2(a). The residues at positions −5 and −9 are involved in the hydrogen-bonding linkage between the first helix (α2 in HptB) and the second helix (α3), whereas the residues at positions +9, +10 and +11 are involved in the hydrogen-bonding linkage between the second helix (α3) and third helix (α4) of the bundle in the inter-helical region [Figs. 4(b) and 4(c)]. The hydrogen-bonding linkages include the hydrogen bonds between (i) the main-chain oxygen of Asp48 at position −9 and the main-chain nitrogen of Leu52 at position −5 [Fig. 4(b)]; (ii) the side-chain oxygen of Asp48 at position −9 and the main-chain nitrogen of Ala51 and Gln50 at positions −6 and −7, respectively [Fig. 4(b)]; (iii) the main-chain oxygen of Ser63 at position +6 and the main-chain nitrogen of Ala68 at position +11 [Fig. 4(c)]; (iv) the main-chain oxygen of Met66 at position +9 and the main-chain nitrogen of Ala68 at position +11 [Fig. 4(c)]; (v) the main-chain oxygen of Gly67 at position +10 and the main-chain nitrogen of Val69 at position +12 [Fig. 4(c)]; (vi) the main-chain oxygen of Ala68 at position +11 and the main-chain nitrogen of Ala72 at position +15 [Fig. 4(c)]. The N and C termini of the second helix (α3), at which the active-site histidine is located, are additionally stabilized due to these hydrogen bonds contributed by the conserved residues.
3.3. Overall structure of PA1611REC
The recombinant PA1611REC protein comprises amino acids 507–651 of hybrid sensor histidine kinase PA1611 and contains the C-terminal receiver domain. The molecular mass of recombinant PA1611REC with a 6×His fusion tag is ∼18 kDa. The crystals of PA1611REC with Mg2+ belong to trigonal P3121 with one protein molecule in the There are 122 amino-acid residues (a.a. 513–634) defined in the final atomic model; the electron-density map is clear enough to trace most parts of the main chains and side chains, except residues 507–512 and 635–651 were not seen, likely because of flexibility. The divalent magnesium ion (Mg2+) from the crystallization solution was found to be coordinated in the active-site cleft.
The overall structure of PA1611REC possesses a conventional (β/α)5 topology as observed in other receiver domains of response regulators despite the fact they share low-sequence identities [Fig. 2(b)]. In this (β/α)5 motif, the five-stranded parallel β-sheets, which are β1 (Glu516–Val520), β2 (Arg540–Val544), β3 (Ala561–Asp565), β4 (Ile591–Thr595) and β5 (Asp613–Ala616), are folded in the central core. The core comprises β1, β3 and β4 which contain conserved hydrophobic residues. The surrounding five α-helices are α1 (Pro524–Ser536), α2 (Gly547–Arg556), α3 (Gly573–Arg581), α4 (Gly601–Gln608) and α5 (Arg621–Trp631), among which two (helices α1 and α5) are located on one side, whereas the others (helices α2, α3 and α4) are on the opposite side. The β-strands and α-helices are connected with loops. Based on the assignment of the secondary structure from the DSSP web server, helix α4 is composed of two 310 helices (Gly601–Glu604 and Cys606–Gln608) and separated at Asn605 (Fig. 6). An inspection of the secondary structural geometry shows that the φ and ψ angles of Asn605 are −150.3 and −100.1°, respectively, which deviate from the normal range for a helical structure. Moreover, the adjacent residue Cys606 is positioned structurally closely to Cys566 (on the loop β3→α3) to form a disulfide bridge. Interestingly, both Cys606 and Cys566 are unique in P1611REC. Two 310 helices separated at Asn605 in helix α4 might be presumably affected by the disulfide formation of Cys606–566. This distinctive structural feature in PA1611REC, which does not exist in other receiver domains of response regulators, might contribute to the interaction with HptB through the residue on loop β4→α4 discussed later.
3.4. Metal binding in the active-site cleft of PA1611REC
In PA1611REC, the conserved aspartate residue, Asp565, responsible for phosphoryl group binding is located in the active-site cleft formed by loops β1→α1 (Glu521 to Asn523), β2→α2 (Ala545 to Asp546) and β3→α3 (Cys566 to Asp572) [Fig. 7(a)]. The hydrogen bond between the carbonyl oxygen of Cys566 and the amide nitrogen of Gly573 on helix α3 contributes to the stabilization of the active-site cleft. Previous studies have shown that the binding of divalent metal ions in the active-site cleft is essential to activate the receiver domain (Stock et al., 1993). The surface electrostatic potential around the active-site cleft of PA1611REC is highly negatively charged and suitable to accommodate metal ions [Fig. 7(a)]. The divalent magnesium ion (Mg2+) is hexa-coordinated, of which three positions coordinate with three water molecules, and three with the main-chain oxygen atom of Arg567 and the side-chain oxygen atoms of Asp522 and Asp565. The three water molecules coordinated to the magnesium ion form further hydrogen bonds with the main-chain oxygen atom of Arg567 and the side-chain oxygen atoms of Asp522 and Asp565. Moreover, one of these water molecules also interacts with Glu521 via hydrogen bonding [Fig. 7(b)].
3.5. The docking model of the HptB–PA1611REC complex
The key point in phosphorelay is the recognition between Hpt proteins/domains and receiver domains of HK or RR. The complex structures of SLN1-R1/YPD1 (PDB entry 2r25; Zhao et al., 2008) and AHK5RD/AHP1 (PDB entry 4euk) have highlighted that there is a general shape complementarity between the kidney-shaped Hpt protein/domain and the convex surface of the receiver domain (Zhao et al., 2008; Bauer et al., 2013). The transfer of the phosphoryl group between PA1611 and HptB was confirmed with an in vitro phosphorelay assay using [γ-32P]ATP in our previous work (Lin et al., 2006; Hsu et al., 2008); however, our attempts to obtain the of HptB in complex with the receiver domain of PA1611 were unsuccessful. We thus modeled the HptB–PA1611REC complex with the HADDOCK web server to elucidate the interactions between HptB and the cognate receiver domain. Previous studies have shown that the divalent metal ion is involved in the phosphotransfer reactions of the TCS transduction pathway; the relevant ion in vivo is presumably Mg2+ (Stock et al., 2000). Divalent Mg2+ in the structure of PA1611REC was hence taken into consideration in the docking process. Moreover, as mentioned earlier, the structures of six HptBs in one are essentially the same; thus, we selected one representative protomer for docking with PA1611REC.
According to the docking model, HptB and PA1611 REC would form a complex in a 1:1 ratio. Most of the interaction surface is composed of the second helix (helix α3) of the four-helix bundle of HptB, the N-terminus of helix α1 of PA1611REC, and the loops (β1→α1, β3→α3 and β4→α4) around the active-site cleft of PA1611REC [Fig. 8(a)]. According to an analysis calculated with the PDBePISA web server, there are 14 hydrogen bonds and salt bridges at the interface, in which 15 residues in HptB and 20 residues in PA1611REC are involved (Fig. 9, Table 3). Among the above 20 residues, the active-site residue Asp565 of PA1611REC also participates in the hydrogen-bonding network between HptB and PA1611REC. Although the active-site residue His57 of HptB is in the interface of the complex model, it is, however, not involved in the hydrogen bonds or salt bridges between the two proteins [Figs. 8(a) and 9]. There are few structural differences between the proteins of the apo form and the HptB–PA1611REC complex model, with an RMSD of 0.44 Å for HptB and 0.43 Å for PA1611REC. Furthermore, an inspection of of the residues of HptB involved in the interface of the complex model revealed that their side chains are positioned in consistent and similar orientations among the six HptB protomers in one suggesting that the docking model serves the representative complex.
|
4. Discussion
4.1. Comparison of HptB with other Hpt proteins/domains
HptB exhibits a low level of overall sequence similarity with other Hpt proteins/domains [Fig. 2(a)], but they all share a common structural feature as an elongated four-helix bundle, which is kidney-shaped with up-and-down topology and stabilized by inter-helix hydrophobic contacts and hydrogen-bond linkages in the inter-helical regions. The conserved residues are located mainly at the helices forming the four-helix bundle and surrounding the active-site histidine (His57 in HptB), indicating that they contribute to the stability of the bundle and maximize the accessibility of the active-site histidine [Fig. 2(a)]. The active-site histidine is located at a similar position, which is near the middle of the second helix (α3 in HptB) of the helix bundle. The major difference is located at the N-terminal helix (α1 in HptB) and C-terminus beyond the four-helix bundle [Fig. 10(a)]. The differences in the segments outside the four-helix bundle might influence the specific interactions between separate Hpt proteins/domains in each species and the cognate receiver domains.
A search for structural homologs with the DALI server (Holm, 2019) showed that HptB bears structural similarity not only to other Hpt proteins/domains but also to proteins without functional similarities, such as flagellin (PDB entry 5maw; Altegoer, et al., 2018) with a Z-score of 9.7, mannose-6-phosphate receptor binding protein 1 (PDB entry 1szi; Hickenbottom et al., 2004) with a Z-score of 8.7 and cytochrome c′ (PDB entry 5b3i; Fujii et al., 2017) with a Z-score of 8.5. The four-helix bundle motif is evidently suitable for diverse biological functions (Xu & West, 1999).
The orientation of helix α1 in HptB is similar to that of helix α2 in AHP1 of Arabidopsis. Helix α2 of AHP1 contributes additional hydrogen bonds and hydrophobic interactions to the interface with the receiver domain of AHK5 (AHK5RD). Compared with the AHK5RD–AHP1 complex, the CheA3HP1 of R. sphaeroides, which lacks extra helices outside the conserved Hpt bundle, has shown a smaller contact area and affinity with its partner, CheY6. The difference in affinity indicates that the extra helix at the N-terminus of the four-helix bundle might play an important role in the interaction with the cognate receiver domains (Bauer et al., 2013). Our docking model showed, however, that helix α1 of HptB is not involved in the interaction with the receiver domain of PA1611. The major role of this helix might therefore be to stabilize the four-helix bundle of HptB or to interact with other subjects, such as another domain of PA1611 or other proteins.
4.2. Structural features of HptB facilitates its interaction with receiver domains
The hydrogen-bond linkages around the side chains of the residues near the active-site histidine have been suggested to facilitate its accessibility to the phosphoryl group and to increase the nucleophilicity of the nitrogen atom at the imidazole ring to improve the efficiency of phosphotransfer (Xu & West, 1999). Moreover, the glycine residue (Gly61 in HptB) at position +4 from the active-site histidine (His57 in HptB) is located on its adjacent ridge in the second helix of the bundle (α3 in HptB) and seems to carve a space for the histidine to expose its imidazole side chain to solvent. It has been suggested that, if an amino acid with a larger side chain replaces glycine, the activity of the active-site histidine to accept and to transfer a phosphoryl group would be disturbed (Kato et al., 1997). In YPD1 and the Hpt domain of ArcB, the phosphorylation efficiencies were decreased when this glycine was replaced by Gln or Asp (Sugawara et al., 2005). Aside from the active-site histidine and the glycine residue at position +4, the lysine residue at position +3 (Lys60 in HptB) is conserved throughout all active Hpt orthologs. The conservation and position imply that this lysine residue is involved in the phosphorelay process. The phosphorylation efficiency of a mutant in which the +3 lysine was substituted with alanine was decreased in YPD1 (Janiak-Spens & West, 2000). In HptB, Lys60 contributes to the formation of the positively charged region around the active-site histidine, which may help to neutralize the negatively charged phosphoryl group as well as form hydrogen bonds to the oxygen atoms of phosphorylated histidine, providing an additional stabilizing force.
Although the sequence identity is low among Hpt proteins/domains, all Hpt proteins/domains share a common structural feature and active site, with some conserved residues that could be important for the structure and function described above. The key structural feature for the interactions between Hpt proteins/domains and receiver domains has been proposed to be the general shape complementarity between them, observed in several TCSs (Kato et al., 1997; Zhao et al., 2008; Bauer et al., 2013). In contrast, non-conserved residues might be important for the specific recognition of molecular partners (Xu & West, 1999; Xu et al., 2009). Although some Hpt proteins/domains could interact with non-cognate partners in vitro and in vivo, the rates of phosphotransfer in these cases might be substantially less than with cognate partners (Rogov et al., 2004; Xu et al., 2009). Because a different arrangement of amino acid residues would alter the affinity and specificity of the interaction between each Hpt and receiver domain, one reasonably suspects that the presence of Hpt proteins/domains could help to diminish or prevent the cross-interaction between distinct pathways in a TCS (Sugawara et al., 2005; Xu et al., 2009).
Aside from the general shape matching, other structural features might be involved in the interaction between Hpt proteins/domains and receiver domains. For example, some Hpt proteins/domains have a hydrophobic cavity in the proximity of the active-site histidine. The hydrophobic interaction has been proposed to contribute to the formation of the Hpt–REC complex; mutation of the hydrophobic residues in that cavity would decrease the phosphotransfer activity (Kato et al., 1999a; Rogov et al., 2004). The active-site histidine (His57) of HptB is surrounded largely by positively charged residues. We hence suggest that the electrostatic interactions and hydrogen bonding also help the binding of HptB and its partner in vivo. Based on an analysis of the complex model with the PDBePISA web server, the residues involved in the interface between HptB and PA1611REC contribute to the hydrogen bonds and salt bridges as well as hydrophobic and van der Waals interactions [Figs. 8(a) and 9, Table 3].
4.3. Comparison of PA1611REC with other receiver domains
The receiver domains typically exist in response regulators and also in some sensor histidine kinases in the hybrid TCS. Although the sequence similarity is not high, the receiver domain of PA1611 shares structural and functional homology with most bacterial receiver domains [Figs. 2(b) and 10(b)] (Lin et al., 2006; Hsu et al., 2008). The structure alignment of the receiver domains listed in Fig. 2(b) shows that the secondary structures are organized in a similar pattern with a typical (β/α)5 topology [Fig. 10(b)]. The central β strands and the regions around the active-site cleft are highly conserved [Fig. 2(b)]. In contrast, the region containing helix α4 and the loop β4→α4 displays a large variability, which causes an overall RMSD above 0.6 Å between PA1611REC and other receiver domains listed in Fig. 2(b). The RMSD values between PA1611REC-Mg2+ and RR468, VcCheY4, CtrAREC, RR02rec, PhoB REC and BfmRN are 0.92, 0.95, 0.67, 0.79, 0.78 and 0.98 Å, respectively.
The conserved residues of the receiver domains are located mainly on β-strands and α-helices. Aside from the aspartate residue that is directly phosphorylated (Asp565 in PA1611REC), highly conserved residues are involved in the process of phosphorelay. The aspartate residue at the end of β1 (Asp522 in PA1611REC) participates in the binding of the divalent metal ion. The structure of the CheY receiver domain with beryllofluoride (phosphoryl group analog) in the active-site cleft (PDB entry 1fqw; Lee et al., 2001) shows that the ligand is located close to the metal ion and is stabilized by Mg2+, conserved Thr87 on β4 with a hydrogen bond, and conserved Lys109 on β5 by salt-bridge interactions (Lee et al., 2001). We hence suggest that Thr595 (Thr87 in CheY) and Lys617 (Lys109 in CheY) in PA1611REC would also interact with the phosphoryl group and stabilize the conformation of the activated receiver domain during phosphorylation. The alanine residue (Ala596 in PA1611REC) following Thr595 was thought to improve access of the conserved Thr/Ser to the phosphorylation site (Bourret, 2010).
Although the α→β loops opposite the active-site cleft demonstrate greater variability in sequence and length, the β3→α3 loop near the active-site cleft is highly conserved. This loop has been suggested to form a universal recognition element across the receiver domain superfamily (Usher et al., 1998). Moreover, an inspection of temperature B-factors of PA1611REC reveals that the β3→α3 loop and the other two loops (β1→α1, β2→α2) surrounding the active-site cleft are relatively rigid, compared with the loops away from the active site (see Fig. S1 of the supporting information). In our complex model, this β3→α3 loop participates in the interaction with HptB. We hence suggest that the β3→α3 loop in PA1611REC would be involved in the binding with HptB in vivo.
4.4. Structural alterations in receiver domains
Conformational changes due to phosphorylation have been assumed to be common in proteins related to signal transduction. For example, in E. coli CheY, a series of structural changes occur near the active-site cleft during phosphorylation (Cho et al., 2000; Lee et al., 2001). The structures of Thermotoga maritima CheY (TMY) indicated that Ser82 and Phe101 might be sensitive to the phosphorylation status of the active-site aspartate and alter the conformation simultaneously. In the unphosphorylated state, the side chain of Phe101 is exposed to solvent; the side-chain hydroxyl of Ser82 points toward Phe101. After TMY is phosphorylated, the side-chain hydroxyl of Ser82 is oriented towards the active site; the aromatic ring of Phe101 is buried (Usher et al., 1998). Similar conformational changes have been reported for the residue substitutions in other receiver domains (Cho et al., 2000; Lee et al., 2001; Zhao et al., 2008). In our determined structure of PA1611REC, there is no phosphate group in the active-site cleft. The side chain of Tyr614 (Phe101 in TMY) is exposed to solvent; the side chain of Thr595 (Ser82 in TMY) points towards it. The results suggest that our PA1611REC is in an inactivated form; there might be conformational changes that occur in the side chain of Thr595 and Tyr614 after phosphorylation, similar to the residue substitutions in E. coli CheY and TMY.
Helix α4 and its adjacent regions display a great variability in several structures of RR receiver domains. In previous studies, these regions were shown to undergo conformational changes in the activated receiver domain (Cho et al., 2000; Lee et al., 2001; Zhao et al., 2008). Aside from helix α4, the N-terminus of helix α1 has been reported to be partially unfolded when CheY interacts with cognate Hpt proteins/domains (Kato et al., 1999a). The conformational alterations of activated REC are varied. For example, in the YPD1/SLN1-R1 complex with Mg2+ and BeF3−, conformational changes of a significant number were observed in comparison with the apo complex. In contrast, there are very few conformational changes in the Spo0B/Spo0F-Mg2+-BeF3− complex compared with the apo complex (Varughese et al., 2006; Zhao et al., 2008). Although the loops and some residues may undergo conformational changes during phosphorylation, there is no dramatic structural rearrangement in the main secondary-structure architecture of receiver domains (Lee et al., 2001). For example, the difference is small between the structures of BeF3−BfmRN and the apo form BfmRN, with an RMSD of 0.315 across 227 Cα atoms (Draughn et al., 2018). The structural differences between PA1611REC in the apo form and the complex model are also small, with an RMSD of 0.43 Å.
4.5. Interactions among the receiver domain, metal ion and Hpt
Previous studies have shown that a divalent metal cation helps to stabilize the phosphorylated active-site aspartate residue of the receiver domain. In the study of EL_LovR, a protein with only one receiver domain in Erythrobacter litoralis, its lack of divalent metal ion was suggested to create electrostatic repulsion among the negatively charged residues around the active-site cleft which destabilize the protein (Ocasio et al., 2015). The manners for metal binding are similar in most receiver domains, in which three coordination positions occur with the atoms of amino-acid residues and the others with solvent molecules. The structure of the CheY receiver domain with the phosphoryl group analog, beryllofluoride, in the active-site cleft (PDB entry1fqw; Lee et al., 2001) shows that beryllofluoride is located close to the metal ion; Oδ on the side chain of the conserved active-site aspartate (Asp57 in CheY) can interact with Mg2+ and the beryllium atom.
The studies of CheY have shown that Mn2+ has the same geometry as Mg2+ in the active-site and can also coordinate to Asp13, Asp57, Asn59 and BeF3− (Lee et al., 2001). The structure of PA1611REC with Ca2+ in the active-site cleft was also determined from a crystal grown in a crystallization solution containing Ca2+ (Fig. S2). The overall structures of PA1611REC-Mg2+ and PA1611REC-Ca2+ are similar with an RMSD of 0.23 Å. These results coincide with suggestions from previous studies that the divalent metal ion would not induce a conformational change of the REC active-site cleft, although it is required for the phosphorylation process (Lee et al., 2001).
The overall arrangement between Hpt proteins/domains and REC is similar in complex structures of SLN1-R1/YPD1, AHK5RD/AHP1 and our docking model of HptB–PA1611REC. However, compared with SLN1-R1/YPD1 and AHK5RD/AHP1 complexes, the slight deflection of the relative positions of HptB and PA1611REC results in fewer residues involved in the interface of our complex model according to an analysis with the PDBePISA web server (Fig. 8). The relative positions of the active-site residues in HptB (His57) and PA1611REC (Asp565) are also slightly different from those in SLN1-R1/YPD1 and AHK5RD/AHP1 complexes because of the minor deflection. We hence suggest the interactions between the two proteins in the HptB–PA1611REC complex model are weaker than those in SLN1-R1/YPD1 and AHK5RD/AHP1 complexes. Perhaps for this reason we could not directly observe the interaction between HptB and PA1611REC in the experiment using and hence could not successfully obtain suitable crystals of the HptB–PA1611REC complex. Because the transfer of the phosphoryl group between the cytosolic part of PA1611 (a.a. 202–652) and HptB has been confirmed with [γ-32P]ATP in our previous work (Lin et al., 2006; Hsu et al., 2008), we surmise that another part in addition to REC of PA1611 or some other molecule is necessary to enhance the interaction with HptB.
According to previous studies and an analysis with the PDBePISA web server, the active-site histidine residues of Hpt-related proteins in SLN1-R1/YPD1, AHK5RD/AHP1 and the HptB–PA1611REC model are involved in the interface of the complex but do not participate in the hydrogen-bond network and salt bridges of the binding surface (Figs. 8 and 9; Table 3). The conserved lysine residue at position +3 and glycine residue at position +4 from the active-site histidine also contribute to the interaction. These two residues in YPD1 and AHP1 are involved in the hydrogen-bond network or salt bridges of the interface, whereas the two residues (Lys60 and Gly61) in HptB are not involved in the hydrogen bonds and salt bridges between the two proteins in the complex model. Previous studies of the SLN1-R1/YPD1 complex described how the side chain of Gln86 in YPD1, which corresponds to Glu79 in HptB, forms a hydrogen bond with Gln1146 of SLN1-R1 (Xu et al., 2003). An analysis of the interface in the HptB–PA1611REC complex model with the PDBePISA web server showed that Glu79 of HptB is also involved in the hydrogen-bond networks between the two proteins [Figs. 8(a) and 9; Table 3]. Glu79 of HptB might contribute to the interaction with the cognate receiver domain in vivo.
The aspartate residue (Asp565) responsible for phosphoryl group binding in PA1611REC forms a hydrogen bond with His54 of HptB according to an analysis with the PDBePISA web server. In contrast, the conserved aspartate residue in AHK5RD is not involved in the interface of the complex. Previous studies of the SLN1-R1/YPD1 complex described that this aspartate residue in SLN1-R1 does not participate in the interface of the complex, although our analysis with the PDBePISA web server shows that it is involved in the interaction between SLN1-R1 and YPD1. The N-terminus of helix α1 and the loop β1→α1, which is one of the loops around the active-site cleft, are highly conserved in structures and sequences among SLN1-R1, AHK5RD and PA1611REC. These two segments contribute largely to the interface of the complex. Loop β3→α3 is structurally conserved and contributes largely to the interaction between HptB and PA1611REC, whereas in loop β3→α3 of SLN1-R1 and AHK5RD, only one non-conserved residue is involved in the binding surface of the complex. Loops β4→α4 of SLN1-R1, AHK5RD and PA1611REC also majorly contribute to the interface of the three complexes, but these loops display a low level of sequence similarity. According to these analyses, we suggest that the loops around the active-site cleft of REC are not only involved in the interaction with Hpt proteins/domains but also might influence the specificity of recognition.
The properties and distributions of the amino-acid residues involved in the interactions vary across different complexes. The HptB–PA1611REC complex model, SLN1-R1/YPD1 and AHK5RD/AHP1 are all well conserved in the structures of their components and the overall organization of the complexes. However, the amino-acid residues contributing to hydrogen bonds, salt bridges and non-polar contacts are distributed differently across the interfaces (Fig. 8) (Zhao et al., 2008; Bauer et al., 2013).
Supporting information
PDB references: PA1611REC-Mg2+, 7c1j; PA1611REC-Ca2+, 7cfw; WT-HptB, 7c1i
Supporting figures. DOI: https://doi.org/10.1107/S2052252520009665/lz5039sup1.pdf
Acknowledgements
We are indebted to the staff at beamlines TLS 13B1, 15A1 and TPS 05A at the National Synchrotron Radiation Research Center (NSRRC) in Taiwan; Masato Yoshimura and Eiki Yamashita at the BL44XU beamline and Taiwan beamline BL12B2 at SPring-8 in Japan for technical assistance.
Funding information
This work was supported in part by the Ministry of Science and Technology (MOST) (grant Nos. 105–2311-B-213–001-MY3, 107–2923-B-213–001-MY3 and 108–2311-B-213–001-MY3 awarded to CJC).
References
Altegoer, F., Mukherjee, S., Steinchen, W., Bedrunka, P., Linne, U., Kearns, D. B. & Bange, K. (2018). Sci. Rep. 8, 11552. CrossRef PubMed Google Scholar
Anjali, Y. B., Yanqi, L., Rakesh, P., Katherine, Y., Xiaoxue, M., Ayush, K. & Kangmin, D. (2019). Int. J. Mol. 20, 1781–1810. Google Scholar
Balasubramanian, D., Schneper, L., Kumari, H. & Mathee, K. (2013). Nucleic Acids Res. 41, 1–20. CrossRef CAS PubMed Google Scholar
Bauer, J., Reiss, K., Veerabagu, M., Heunemann, M., Harter, K. & Stehle, T. (2013). Mol. Plant. 6, 959–970. CrossRef CAS PubMed Google Scholar
Bell, C. H., Porter, S. L., Strawson, A., Stuart, D. I. & Armitage, J. P. (2010). PLoS Biol. 8, e1000306. CrossRef PubMed Google Scholar
Bent, C. J., Isaacs, N. W., Mitchell, T. J. & Riboldi-Tunnicliffe, A. (2004). J. Bacteriol. 186, 2872–2879. Web of Science CrossRef PubMed CAS Google Scholar
Biswas, M., Dey, S., Khamrui, S., Sen, U. & Dasgupta, J. (2013). PLoS One, 8, e73923–e73923. CrossRef CAS PubMed Google Scholar
Bourret, R. B. (2010). Curr. Opin. Microbiol. 13, 142–149. Web of Science CrossRef CAS PubMed Google Scholar
Capra, E. J., Perchuk, B. S., Lubin, E. A., Ashenberg, O., Skerker, J. M. & Laub, M. T. (2010). PLoS Genet. 6, e1001220. Web of Science CrossRef PubMed Google Scholar
Casino, P., Rubio, V. & Marina, A. (2009). Cell, 139, 325–336. Web of Science CrossRef PubMed CAS Google Scholar
Cho, H. S., Lee, S. Y., Yan, D., Pan, X., Parkinson, J. S., Kustu, S., Wemmer, D. E. & Pelton, J. G. (2000). J. Mol. Biol. 297, 543–551. Web of Science CrossRef PubMed CAS Google Scholar
Draughn, G. L., Milton, M. E., Feldmann, E. A., Bobay, B. G., Roth, B. M., Olson, A. L., Thompson, R. J., Actis, L. A., Davies, C. & Cavanagh, J. (2018). J. Mol. Biol. 430, 806–821. CrossRef CAS PubMed Google Scholar
Emsley, P., Lohkamp, B., Scott, W. G. & Cowtan, K. (2010). Acta Cryst. D66, 486–501. Web of Science CrossRef CAS IUCr Journals Google Scholar
Fujii, S., Oki, H., Kawahara, K., Yamane, D., Yamanaka, M., Maruno, T., Kobayashi, Y., Masanari, M., Wakai, S., Nishihara, H., Ohkubo, T. & Sambongi, Y. (2017). Protein Sci. 26, 737–748. CrossRef CAS PubMed Google Scholar
Hickenbottom, S. J., Kimmel, A. R., Londos, C. & Hurley, J. H. (2004). Structure, 12, 1199–1207. CrossRef PubMed CAS Google Scholar
Holm, L. (2019). Bioinformatics, 35, 5326–5327. Web of Science CrossRef CAS PubMed Google Scholar
Hsu, J. L., Chen, H. C., Peng, H. L. & Chang, H. Y. (2008). J. Biol. Chem. 283, 9933–9944. CrossRef PubMed CAS Google Scholar
Hubbard, S. J. & Thornton, J. M. (1993). NACCESS. Department of Biochemistry and Molecular Biology, University College, London, UK. Google Scholar
Janiak-Spens, F. & West, A. H. (2000). Mol. Microbiol. 37, 136–144. PubMed CAS Google Scholar
Kabsch, W. & Sander, C. (1983). Biopolymers, 22, 2577–2637. CrossRef CAS PubMed Web of Science Google Scholar
Kaserer, A. O., Andi, B., Cook, P. F. & West, A. H. (2010). Methods Enzymol. 471, 59–75. CrossRef CAS PubMed Google Scholar
Kato, M., Mizuno, T., Shimizu, T. & Hakoshima, T. (1997). Cell, 88, 717–723. CrossRef CAS PubMed Web of Science Google Scholar
Kato, M., Shimizu, T., Mizuno, T. & Hakoshima, T. (1999a). Acta Cryst. D55, 1257–1263. CrossRef CAS IUCr Journals Google Scholar
Kato, M., Shimizu, T., Mizuno, T. & Hakoshima, T. (1999b). Acta Cryst. D55, 1842–1849. CrossRef CAS IUCr Journals Google Scholar
Krissinel, E. & Henrick, K. (2007). J. Mol. Biol. 372, 774–797. Web of Science CrossRef PubMed CAS Google Scholar
Lee, S. Y., Cho, H. S., Pelton, J. G., Yan, D., Berry, E. A. & Wemmer, D. E. (2001). J. Biol. Chem. 276, 16425–16431. CrossRef PubMed CAS Google Scholar
Lin, C. T., Huang, Y. J., Chu, P. H., Hsu, J. L., Huang, C. H. & Peng, H. L. (2006). Res. Microbiol. 157, 169–175. CrossRef PubMed CAS Google Scholar
Lyczak, J. B., Cannon, C. L. & Pier, G. B. (2002). Clin. Microbiol. Rev. 15, 194–222. Web of Science CrossRef PubMed CAS Google Scholar
Murshudov, G. N., Skubák, P., Lebedev, A. A., Pannu, N. S., Steiner, R. A., Nicholls, R. A., Winn, M. D., Long, F. & Vagin, A. A. (2011). Acta Cryst. D67, 355–367. Web of Science CrossRef CAS IUCr Journals Google Scholar
National Nosocomial Infections Surveillance System (2004). Am. J. Infect. Control, 32, 470–485. CrossRef PubMed Google Scholar
Ocasio, V. J., Corrêa, F. & Gardner, K. H. (2015). Biochemistry, 54, 1353–1363. CrossRef CAS PubMed Google Scholar
Otwinowski, Z. & Minor, W. (1997). Methods Enzymol. 276, 307–326. CrossRef CAS PubMed Web of Science Google Scholar
Posas, F., Wurgler-Murphy, S. M., Maeda, T., Witten, E. A., Thai, T. C. & Saito, H. (1996). Cell, 86, 865–875. CrossRef CAS PubMed Web of Science Google Scholar
Rodrigue, A., Quentin, Y., Lazdunski, A., Méjean, V. & Foglino, M. (2000). Trends Microbiol. 8, 498–504. Web of Science CrossRef PubMed CAS Google Scholar
Rogov, V. V., Bernhard, F., Löhr, F. & Dötsch, V. (2004). J. Mol. Biol. 343, 1035–1048. CrossRef PubMed CAS Google Scholar
Ruszkowski, M., Brzezinski, K., Jedrzejczak, R., Dauter, M., Dauter, Z., Sikorski, M. & Jaskolski, M. (2013). FEBS J. 280, 3709–3720. CrossRef CAS PubMed Google Scholar
Skerker, J. M., Perchuk, B. S., Siryaporn, A., Lubin, E. A., Ashenberg, O., Goulian, M. & Laub, M. T. (2008). Cell, 133, 1043–1054. Web of Science CrossRef PubMed CAS Google Scholar
Sola, M., Gomis-Ruth, F. X., Serrano, L., Gonzalez, A. & Coll, M. (1999). J. Mol. Biol. 285, 675–687. Web of Science CrossRef CAS PubMed Google Scholar
Stock, A. M., Martinez-Hackert, E., Rasmussen, B. F., West, A. H., Stock, J. B., Ringe, D. & Petsko, G. A. (1993). Biochemistry, 32, 13375–13380. CrossRef CAS PubMed Web of Science Google Scholar
Stock, A. M., Mottonen, J. M., Stock, J. B. & Schutt, C. E. (1989). Nature, 337, 745–749. CrossRef CAS PubMed Google Scholar
Stock, A. M., Robinson, V. L. & Goudreau, P. N. (2000). Annu. Rev. Biochem. 69, 183–215. Web of Science CrossRef PubMed CAS Google Scholar
Stover, C. K., Pham, X. Q., Erwin, A. L., Mizoguchi, S. D., Warrener, P., Hickey, M. J., Brinkman, F. S., Hufnagle, W. O., Kowalik, D. J., Lagrou, M., Garber, R. L., Goltry, L., Tolentino, E., Westbrock-Wadman, S., Yuan, Y., Brody, L. L., Coulter, S. N., Folger, K. R., Kas, A., Larbig, K., Lim, R., Smith, K., Spencer, D., Wong, G. K., Wu, Z., Paulsen, I. T., Reizer, J., Saier, M. H., Hancock, R. E., Lory, S. & Olson, M. V. (2000). Nature, 406, 959–964. CrossRef PubMed CAS Google Scholar
Sugawara, H., Kawano, Y., Hatakeyama, T., Yamaya, T., Kamiya, N. & Sakakibara, H. (2005). Protein Sci. 14, 202–208. Web of Science CrossRef PubMed CAS Google Scholar
Terwilliger, T. C. & Berendzen, J. (1999). Acta Cryst. D55, 849–861. Web of Science CrossRef CAS IUCr Journals Google Scholar
Thomason, P. & Kay, R. (2000). J. Cell Sci. 113, 3141–3150. PubMed CAS Google Scholar
Touw, W. G., Baakman, C., Black, J., te Beek, T. A., Krieger, E., Joosten, R. P. & Vriend, G. (2015). Nucleic Acids Res. 43, D364–D368. Web of Science CrossRef CAS PubMed Google Scholar
Usher, K. C., De La Cruz, A. F., Dahlquist, F. W., James Remington, S., Swanson, R. V. & Simon, M. I. (1998). Protein Sci. 7, 403–412. CrossRef CAS PubMed Google Scholar
Varughese, K. I., Tsigelny, I. & Zhao, H. (2006). J. Bacteriol. 188, 4970–4977. Web of Science CrossRef PubMed CAS Google Scholar
Willett, J. W., Herrou, J., Briegel, A., Rotskoff, G. & Crosson, S. (2015). Proc. Natl Acad. Sci. USA, 112, E3709–E3718. CrossRef CAS PubMed Google Scholar
Winn, M. D., Ballard, C. C., Cowtan, K. D., Dodson, E. J., Emsley, P., Evans, P. R., Keegan, R. M., Krissinel, E. B., Leslie, A. G. W., McCoy, A., McNicholas, S. J., Murshudov, G. N., Pannu, N. S., Potterton, E. A., Powell, H. R., Read, R. J., Vagin, A. & Wilson, K. S. (2011). Acta Cryst. D67, 235–242. Web of Science CrossRef CAS IUCr Journals Google Scholar
Wolanin, P. M., Thomason, P. A. & Stock, J. B. (2002). Genome Biol. 3, reviews301. CrossRef Google Scholar
Xu, Q., Carlton, D., Miller, M. D., Elsliger, M. A., Krishna, S. S., Abdubek, P., Astakhova, T., Burra, P., Chiu, H. J., Clayton, T., Deller, M. C., Duan, L., Elias, Y., Feuerhelm, J., Grant, J. C., Grzechnik, A., Grzechnik, S. K., Han, G. W., Jaroszewski, L., Jin, K. K., Klock, H. E., Knuth, M. W., Kozbial, P., Kumar, A., Marciano, D., McMullan, D., Morse, A. T., Nigoghossian, E., Okach, L., Oommachen, S., Paulsen, J., Reyes, R., Rife, C. L., Sefcovic, N., Trame, C., Trout, C. V., van den Bedem, H., Weekes, D., Hodgson, K. O., Wooley, J., Deacon, A. M., Godzik, A., Lesley, S. A. & Wilson, I. A. (2009). J. Mol. Biol. 390, 686–698. CrossRef PubMed CAS Google Scholar
Xu, Q., Porter, S. W. & West, A. H. (2003). Structure, 11, 1569–1581. Web of Science CrossRef PubMed CAS Google Scholar
Xu, Q. & West, A. H. (1999). J. Mol. Biol. 292, 1039–1050. Web of Science CrossRef PubMed CAS Google Scholar
Zhao, X., Copeland, D. M., Soares, A. S. & West, A. H. (2008). J. Mol. Biol. 375, 1141–1151. CrossRef PubMed CAS Google Scholar
Zundert, G. C. P. van, Rodrigues, J. P. G. L. M., Trellet, M., Schmitz, C., Kastritis, P. L., Karaca, E., Melquiond, A. S. J., van Dijk, M., de Vries, S. J. & Bonvin, A. M. J. J. (2016). J. Mol. Biol. 428, 720–725. Web of Science PubMed Google Scholar
This is an open-access article distributed under the terms of the Creative Commons Attribution (CC-BY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.