INTRODUCTION
The TRAF family member-associated NF-kappa-B activator (TANK)-binding kinase 1 (TBK1, also referred to as NAK or T2K) was initially discovered by yeast 2-hybrid screening by using the TANK protein as a bait and was found to be a nuclear factor (NF)-κB-activating kinase (NAK) homologous to the IκB kinases IKKα and IKKβ, which are well known as the key kinases involved in the regulation of NF-κB activation (
Pomerantz and Baltimore, 1999). However, the functions of TBK1 and its close homolog IκB kinase epsilon (IKKε, also known as IKKi) were obscure until it was reported that the TBK1/IKKi complex plays an important role in the induction of type I interferon (IFN) expression (
Fitzgerald et al., 2003). There are at least 2 pathways that can activate the TBK1/IKKi complex. One is through the regulatory proteins toll-like receptor 3/4 (TLR3/TLR4), TIR domain-containing adapter molecule 1 (TICAM-1 or TRIF), and TNF receptor-associated factor 3 (TRAF3). The second is through retinoid-inducible gene 1 (RIG-1)-like receptors (RLRs), interferon β promoter stimulator protein 1 (IPS1; also known as MAVS, VISA, Cardif), mediator of IRF3 activation (MITA; also known as STING) and TRAF3 (
Baccala et al., 2007). The TLR3/4 and RLRs belong to the pattern recognition receptor (PRR) family, which recognizes pathogen-associated molecular patterns (PAMPs) and mediates innate immune responses against invading microorganisms (
Chen et al., 2009). Upon activation, TBK1/IKKi phosphorylates the IFN regulatory factors IRF3 and IRF7 (IRF3/7), leading to their translocation into the nucleus and activation of type I IFN expression (
Tamura et al., 2008). The IFNs and their regulated downstream gene products mediate the inhibition of viral replication, clearance of virus-infected cells, and enhancement of immune responses (
Zhong et al., 2008).
Because TBK1 is the key player in this cellular antiviral effect, some viruses have developed defensive mechanisms by antagonizing TBK1 function. For example, γ
134.5 protein of herpes simplex virus type 1 associates with TBK1 and inhibits TBK1-mediated phosphorylation of IRF3, thereby replicating more efficiently in host cells (
Verpooten et al., 2009). In addition to the induction of IFNs, TBK1 possesses various additional functions such as promoting antibacterial autophagy by phosphorylation of the autophagy receptor optineurin (
Wild et al., 2011), attenuating retroviral budding by phosphorylation of the endosomal sorting complex required for transport-I (ESCRT-I) subunit VPS37C (
Da et al., 2011), and involvement in malignant transformation (
Shen and Hahn, 2011). Thus, the structural determination of this molecule would greatly assist in understanding its functional mechanisms.
The IKKα/IKKβ and IKK-related kinases TBK1/IKKi are both members of the IKK kinase family. While IKKα and IKKβ share a 52% overall homology with one another, TBK1 and IKKi share a 61% overall homology with one another. However, the homology between classical IKK kinases and IKK-related kinases is limited, with 27% overall homology and 33% homology within the catalytic kinase domain (
Hiscott, 2007). Hence, substantial structural difference may exist between the two families. Recently, Xu et al. reported the crystal structure of
Xenopus laevis IKKβ (xIKKβ) without the C-terminal NEMO-binding domain (
Xu et al., 2011). The xIKKβ has a 72% sequence identity with human IKKβ. It has a trimodular architecture comprising a kinase domain, a ubiquitin-like domain (ULD), and an elongated α-helical scaffold/dimerization domain (SDD).
Bioinformatic analysis suggests that TBK1 folds into a 3-domain structure similar to that of IKKβ: an N-terminal Ser/Thr kinase domain, a middle ULD, and C-terminal elongated helices (Fig. 1A). The ULDs are reported to have a variety of functions, including sulfur transfer (
Hochstrasser, 2009) and signal transduction (
Welchman et al., 2005). Moreover, many proteins in the ubiquitination/proteasomal delivery system and from the deubiquitinase family possess one or more ULDs. One of the well-defined functions of ULDs is to mediate protein-protein interactions in signal transduction and proteasomal delivery (
Welchman et al., 2005). It has been reported that the double mutation L352A/I353A in the ULD prevents TBK1 from activating IFN-β or RANTES promoters (
Ikeda et al., 2007). In all these IKKs and IKK-related kinases, deletion of the ULD results in the loss of the kinase activity (
May et al., 2004;
Ikeda et al., 2007;
Xu et al., 2011), suggesting that ULD is absolutely required for the kinase domain to form an enzymatically active conformation. Hence, ULD is an indispensable part of TBK1; however, the underlying mechanism remains unknown.
In this study, we describe the crystal structure of the ULD of human TBK1. We compared this domain with some structurally similar domains, identified several specific and conserved residues by multiple sequence alignment, and analyzed the surface properties of this domain. Our structure provides detailed information about ULD-mediated potential intra- and inter-molecular interactions between TBK1 and other IKK-related kinases.
RESULTS
Overall structure
We overexpressed human TBK1 ULD covering amino acids 302–383 and finally obtained its crystals (Fig. 1) (see MATERIALS AND METHODS). The crystal adopts a P2(1) space group, with two molecules in an asymmetric unit; its structure was determined at 1.8 Å resolution with the single-wavelength anomalous diffraction (SAD) method by using a selenomethionyl derivative crystal. A summary of crystallographic data and refinement statistics data is shown in Table 1. The structure of ULD contains 5 β strands and 1 α helix. A Dali search (
Holm and Rosenström, 2010) from Protein Data Bank (PDB, http://www.pdb.org) identified hundreds of structures, predominantly the ULDs, similar to TBK1 ULD. The TBK1 ULD structure has a main chain root mean squared deviation of 2.3 Å with human ubiquitin (PDB: 1ubq) (Fig. 1B). Similar to ubiquitin, the ULD of human TBK1 adopts a β-grasp fold, with 5 β strands packed against an α helix (
Burroughs et al., 2007). The β1 (R308-S315) and β5 (I379-S383) strands are adjacent and parallel to one another, forming the inner 2 strands. The β2 (M319-H327) strand flanks β1, and the β3 (Q350-I353) strand flanks β5 on either side. Together with β4 (R357-L359), these structures form the outer strands. These 5 strands grasp a 12-residue α1 helix (A332-T343), which forms the conserved hydrophobic core of ULD. A long and primarily unstructured region (Y354-P378) connecting β3 and β5 covers the bottom part of this motif (Fig. 1B).
It has been suggested that TBK1 and IKKi can form both homodimers and heterodimers (
Chau et al., 2008). During protein purification of TBK1 ULD with gel filtration, we observed that the elution volume of TBK1 ULD was the same as that of ubiquitin, which exists as a monomer in solution (
Bolton et al., 2001) (Fig. 1C). Therefore, we propose that TBK1 ULD also exists as a monomer in solution. Although the ULDs have been reported to mediate protein-protein interactions in signal transduction, this observation indicates that the TBK1 ULD cannot form homodimers in solution. Therefore, we suggest that the formation of TBK1 homodimer is not mediated by the ULD alone.
The conserved Ile44 hydrophobic patch
A Dali search (
Holm and Rosenström, 2010) from PDB showed that the TBK1 ULD was similar to ULDs that broadly exist in different functional proteins. The ULDs can be classified into 2 categories: small ubiquitin-like modifiers and ULD containing multi-domain proteins. There are at least 10 types of ubiquitin-like modifiers such as Nedd8, SUMO, Fat10, and ISG15 (
Hochstrasser, 2009). The integral ULDs can be further divided into several subfamilies: UBX, a module possessed by some proteasomal shuttle factors; PB1 (Phox and Bem1), which heterodimerize by electrostatic force to form a functional signaling complex; Ras-binding domain, which binds and regulates Ras through its β1 and β2 strands; DCX of doublecortin protein, which binds microtubules; and classic ULD, such as hPLIC_N (human homologs of the yeast ubiquitin-like Dsk2 protein), Rad23N, and NIRF_N (amino-terminal ULD of NP95 and NIRF) (
Buchberger et al., 2001;
Lytle et al., 2004;
Grabbe and Dikic, 2009).
Of these subfamilies, we selected 5 representative structures that were structurally similar to TBK1 ULD and analyzed their sequence and structural characteristics. These structures included human ubiquitin (PDB 1ubq), Nedd8 (PDB 1r4m), the ULD of hHR23A (PDB 1qze), the ULD of NP95 (PDB 2faz), and the ULD of yeast Dsk2 (PDB 2bwf). Four of these form complex structures with their binding partners, which may provide further information about protein interaction sites.
Structure-based sequence alignment of these proteins reveals several conserved residues of TBK1 ULD, including Ile311, Val313, Leu316, Ile324, Ile326, Val339, Ile345, Gln350, Leu352, Ile353, Ile379, Val381, and Val382 (Fig. 2A). Alignment of TBK1 and IKKi in various species showed additional conserved residues, including Phe314, Phe335, and Pro378 (data not shown). Of these conserved residues, Ile311, Val313, Ile324, Ile326, Phe335, Val339, Ile345, Leu352, Pro378, Ile379, and Val381 show very low percentages of solvent-accessible surface areas, suggesting that these residues are primarily structural in nature and contribute to the hydrophobic core of this ubiquitin-like fold (Fig. 2B).
The surface residues Leu316, Ile353, and Val382 are distinct from the hydrophobic core residues and may be involved in protein-protein interactions between ULDs and their binding partners. Interestingly, Leu316, Ile353, and Val382 of TBK1 ULD correspond to the hydrophobic patch formed by Leu8, Ile44, and Val70 in ubiquitin that is involved in interaction with ubiquitin-binding partners (
Hurley et al., 2006) (Fig. 2A). In the case of ubiquitin, the majority of ubiquitin-binding domains—e.g. UIM, CUE, UBA, GAT, NZF, and UEV domains—bind to ubiquitin through the hydrophobic patch centered on Ile44, known as the Ile44 hydrophobic patch (
Hurley et al., 2006). Moreover, deubiquitinases (DUBs) interact with ubiquitin through the Ile44 patch (
Komander et al., 2009). Of the 5 identified homologous structures, hHR23A, yDsk2, and Nedd8 possess the same Ile44 hydrophobic surface patch as that present in ubiquitin and use this region to interact with their corresponding binding partners S5a-UIM, Ufd2, and UBA3, respectively (
Mueller and Feigon, 2003;
Walden et al., 2003;
Hanzelmann et al., 2010). Structural alignment shows that these 3 residues (Leu, Ile, and Val) have similar distributions on the surface of ULDs (Fig. 2C). Based on these results, we suggest that the Ile44 hydrophobic patch, including Leu316, Ile353, and Val382, on TBK1 ULD, should be considered potential protein-binding site.
The Ile44 hydrophobic patch of TBK1 ULD may mediate the ULD-SDD interaction
Recently, the crystal structure of
X. laevis IKKβ (xIKKβ) was reported (
Xu et al, 2011). The structure comprises 3 domains: a kinase domain, a ULD domain, and an SDD domain. The ULD maintains extensive interactions with SDD, specifically using the Ile44 hydrophobic patch to interact with the α helices of SDD (
Xu et al, 2011). We aligned the ULD sequences from human TBK1, IKKα, IKKβ, and IKKi with
X. laevis IKKβ and found that most residues involved in the ULD-SDD interactions of xIKKβ, especially the residues corresponding to the Ile44 hydrophobic patch, were also conserved in the protein homologs (Fig. 3A). In TBK1, residues F314 and F380 lie adjacent to Leu316, Ile353, and Val382. Together, they form a hydrophobic area comparable to that in xIKKβ ULD. It has been reported that the L352A/I353A double mutant of TBK1 ULD fails to interact with the full-length TBK1 and cannot activate the IFNβ promoter (
Ikeda et al., 2007). Considering the significant homology between IKK and IKK-related kinase, we speculate that residues Leu316, Ile353, and Val382 comprising the Ile44 hydrophobic patch of TBK1 ULD also play a role in the ULD-SDD interaction in the TKB1 molecule (Fig. 3B).
Surface charge comparison of ULD between IKK and IKK-related kinase
Because the ubiquitin-like domain is usually a protein-protein interaction motif (
May et al., 2004;
Ikeda et al., 2007;
Xu et al., 2011) and the IKK and IKK-related kinase possess different binding partners in the innate immune pathways, we compared their ULDs to examine the structural mechanism behind their specificities. Because human IKKi ULD shares a 50% sequence identity with TBK1, we performed homology modeling of IKKi ULD by using our structure as the template in the Swiss-Model web server (
Guex and Peitsch, 1997;
Schwede et al., 2003;
Arnold et al., 2006). We also modeled the ULD structures of human IKKβ and IKKα by using the xIKKβ ULD structure, because xIKKβ ULD has 57% and 32% sequence identity with the ULD of human IKKβ and IKKα, respectively. Subsequently, we analyzed and compared the different features of electrostatic surfaces between these ULDs (Fig. 4).
We observed a prominent positively charged area on one of the surfaces of TBK1/IKKi (Fig. 4I & J, highlighted in black circles), whereas the homologous region of IKKα/IKKβ was quite different (Fig. 4G & H). A conserved Arg (Arg308 in TBK1 and Arg308 in IKKi) from different species of TBK1/IKKi on the β1 strand contributes to this positively charged area. In contrast, IKKα/IKKβ shows prominent negative charges on the loop connecting the β2 strand and the α1 helix where several negatively charged Glu and Asp residues are situated (Asp 331 and Glu 332 in IKKα and Glu 329, Asp330, and Glu331 in IKKβ) (Fig. 4G & H, highlighted in black circles).
We also observed a negatively charged surface patch in all these ULDs when viewing the structure from the C-terminus of xIKKβ SDD (Fig. 4K–O, highlighted in black circles). This patch corresponds to the amino acid sequence surrounding β3 where Glu (Glu351 in TBK1) contributes to a portion of this negative charge. Because the amino acid sequence “350QELI353” in the β3 strand of TBK1 ULD is conserved in TBK1, IKKi, IKKβ, and IKKα (Fig. 3A), we suggest that this negative charge may have additional, unknown functions.
DISCUSSION
In this study, we present the crystal structure of the TBK1 ULD. Based on sequence and structural analysis, we showed that TBK1 ULD has a ubiquitin-like structure and an Ile44 hydrophobic patch, which is conserved among ULDs and IKK and IKK-related proteins. Analysis of the published complex structures in PDB showed that the exposed face of the β-sheets in the Ile44 patch mediates a greater number of interactions across the β-grasp fold (
Burroughs et al., 2007). The results of our sequence and structural alignment with xIKKβ suggest that this hydrophobic patch is involved in ULD-SDD interactions in TBK1 and other IKK and IKK-related proteins. Because many protein domains UIM, UBA, CUE, GAT (
Hurley et al., 2006) and NEMO (IKKγ) use their α helical regions to bind to ubiquitin (
Rahighi et al., 2009), it is not surprising that the ULDs tend to bind the α-helical SDD.
To date, research on the use of ULD as an interaction module is supported by the following findings: (1) IKKβ ULD binds weakly to its specific substrate IκB (
Xu et al., 2011), (2) ULD of human IKKβ alone can associate with Homer 3 coiled-coil/leucine zipper domain (
Yatherajam et al., 2010), and (3) TBK1 ULD interacts with its substrate (IRF3) through the IAD domain of IRF3 (
Ikeda et al., 2007). Hence, there must exist additional binding surfaces apart from the conserved Ile44 hydrophobic patch described above. We analyzed the electrostatic surfaces of the ULDs in IKK and IKK-related proteins by performing a structural alignment (Fig. 4) and observed several conserved and differential features that may be relevant to IKK specificities for their interaction partners. Moreover, after an extensive investigation of the known structures of ULD complexes with their binding partners in PDB, 3 distinct regions on ULD emerge as hotspots for protein interaction (Fig. 5) in addition to the Ile44 patch on the exposed face of β sheets. These 3 areas correspond to the 3 modes of interaction that are compatible with TBK1 ULD. In the first type of interaction, ULD uses its β1 strand to form a continuous intermolecular β sheet with its binding partner (Fig. 5A), represented by RBD-Ras and SUMO-SIM complexes (
Nassar et al., 1995;
Namanja et al., 2012). In the second type of interaction, the binding site is almost exactly opposite the Ile44 hydrophobic patch, comprising the α1 and β2-α1 loop (Fig. 5B). This is observed in the yHub1-HIND structure, in which a negatively charged Asp residue located in the β2-α1 loop forms a salt bridge with HIND (
Mishra et al., 2011). Interestingly, IKKα/IKKβ also display prominent negative charges in the same region, while TBK1/IKKi do not. In the third type of interaction, ubiquitin interacts with an A20 zinc-finger domain (A20 ZnF) through the region between β4 and β5 (
Lee et al., 2006) (Fig. 5C). This region is primarily an unstructured loop, with a length that is different between IKK and IKK-related family members. In the xIKKβ structure, several residues are missing from the electron-density map, suggesting the unordered nature of this region, which may be involved in interaction with its binding partners. The proposed binding sites are supported by the results from an online web server for the prediction of protein interaction sites PredUs and meta-PPISP (
Qin and Zhou, 2007;
Zhang et al., 2011). Taken together, we proposed 3 potential intermolecular modes of interaction for TBK1 ULD. Although this hypothesis requires further work, our structure and analysis may pave the way for further efforts to identify and characterize the ULD-partner interaction.
In conclusion, we provide a high-resolution crystal structure of human TBK1 ULD, a pivotal kinase in IFN production, and analyze its special features for its potential interaction. These results would be important for understanding the structural and functional differences of IKK and IKK-related kinase.
MATERIALS AND METHODS
Protein expression and purification
Recombinant human TBK1 ULD (amino acids 302–383) was expressed in E. coli as a GST fusion protein using a pGEX-6p-1 vector (GE Healthcare). After induction with 0.2 mmol/L IPTG at 16°C for 20 h, the cells were harvested and lysed in PBS buffer (pH 7.4) by sonication. The debris and the insoluble fraction were removed by centrifugation at 20,000 g for 30 min and the supernatant was loaded onto a glutathione affinity column (GE Healthcare). After extensive washing with PBS buffer, the GST fusion protein was cleaved in situ by PreScission Protease digestion in a buffer containing 50 mmol/L Tris-HCl, 150 mmol/L NaCl, 1 mol/L urea, and 5% glycerol (pH 8.0). The eluted protein was further purified by Q ion exchange chromatography followed by Superdex-75 gel filtration chromatography (GE Healthcare). Selenomethionine-derivatized protein was purified in a similar manner. The purified protein contains the additional vector-derived amino acid sequence GPLGS at the N terminus and LERPHRD at the C terminus.
The budding yeast ubiquitin was expressed in E. coli as a his-tagged protein using a pET-28a vector. Yeast ubiquitin was purified by nickel affinity column, and applied to Superdex-75 gel filtration chromatography.
Crystallization
Purified samples were concentrated to approximately 10 mg/mL in 50 mmol/L Tris-HCl, 150 mmol/L NaCl, 1 mol/L urea, 5% glycerol (pH 8.0). Crystals were grown by the hanging-drop vapor diffusion method at 16°C by mixing 1 μL protein with 1 μL reservoir solution. The best diffracting crystals were obtained in a reservoir solution containing 0.1 mol/L BIS-TRIS, 0.2 mol/L NaCl, and 21% PEG3350 (pH 5.3). For selenomethionine-derivatized protein, the best quality crystals were grown in a condition containing 0.1 mol/L BIS-TRIS, 0.2 mol/L lithium sulfate monohydrate, and 25% PEG3350 (pH 5.5).
Data collection and structure determination
Native and SeMet MAD data sets were collected on beamline BL5A (KEK-PF, Japan; λ = 0.978 Å) and BL17U (SSRF, China; λ = 0.979 Å), respectively. Data were processed using the HKL2000 program suite (
Minor, 1997). The structure of TBK1 ULD was solved by the SAD method. Heavy atom positions were located by SHELXD (
Sheldrick, 2008). An initial model was automatically built by Phenix Autosol (
Adams et al., 2010), manually modified with Coot (
Emsley and Cowtan, 2004), and refined with Phenix Refine. The final model contains no missing fragment and has an
Rwork of 17.0% and an
Rfree of 20.4%. Data scaling, refinement, and validation statistics are shown in Table 1. The atomic coordinates and diffraction data have been deposited in the Protein Data Bank (4EFO). Structural figures were prepared using PyMOL (http://www.pymol.org).
Higher Education Press and Springer-Verlag Berlin Heidelberg 2012