|
Advertisement | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
J. Biol. Chem., Vol. 279, Issue 12, 11926-11936, March 19, 2004
Specificity of the Interaction between Ubiquitin-associated Domains and Ubiquitin*![]() ![]() From the Department of Chemistry and Biochemistry, University of California, Los Angeles, California 90095-1569
Received for publication, November 25, 2003 , and in revised form, December 31, 2003.
Ubiquitin-associated (UBA) domains are found in a large number of proteins with diverse functions involved in ubiquitination, DNA repair, and signaling pathways. Recent studies have shown that several UBA domain proteins interact with ubiquitin (Ub), specifically p62, the phosphotyrosine-independent ligand of the SH2 domain of p56lck; HHR23A, a human nucleotide excision repair protein; and DDI1, another damage-inducible protein. NMR chemical shift mapping reveals that Ub binds specifically but weakly to a conserved hydrophobic epitope on HHR23A UBA(1) and UBA(2) and that the UBA domains bind on the hydrophobic patch on the surface of the five-stranded -sheet of Ub. Models of the UBA(1)-Ub and UBA(2)-Ub complexes obtained from de novo docking reveal different orientations of the UBA domains on the Ub surface compared with those obtained by homology modeling with the related CUE domains, which also bind Ub. Our results suggest that UBA domains may interact with Ub as well as other proteins in more than one way while utilizing the same binding surface.
The lifespan of proteins inside and outside a cell is tightly regulated by the ubiquitin-proteasome system, and numerous studies show that protein degradation is tightly interlocked with cell cycle progression and is therefore an integral part of transduction pathways and other cellular processes (1-4). For protein degradation by the Ub/proteasome1 system, the target proteins need to be tagged with a poly-Ub chain. These covalent complexes are then recognized and degraded by the 26 S proteasome (1, 2). The principle mechanism of this covalent modification has been identified: an enzyme cascade known as E1-E2-E3 is responsible for activation and transfer of Ub onto the target protein in a linkage-specific manner (1, 5).
The 26 S proteasome is formed by a 20 S cylindrical proteolytically active subunit and two 19 S regulatory subunits (1, 2, 6, 7). The 19 S particles represent the lid of the proteasome and regulate the access to the proteolysis (8). Although the polyubiquitinated substrate seems to be recognized by the S5a subunit in the 19 S particle (9-11), additional contacts between poly-Ub chains and parts of the 19 S regulatory subunit have been identified (12). Deletion studies indicate that other polyubiquitin-binding sites must exist (13). Monoubiquitination is not sufficient for targeting proteins to the proteasome, however; assembly of a poly-Ub chain of at least four Ub moieties is required to create a degradation signal (11, 14, 15). Although Ub contains seven lysine residues, they are not used with the same frequency in poly-Ub chain assembly. The predominant linkages observed are Lys48-Gly76 (16), Lys29-Gly76 (17), and Lys63-Gly76 (18), of which Lys48-linked chains appear to be the most frequent degradation signal. Poly-Ub assembly via Lys29 and Lys63 is less common, and chain formation via Lys63 seems to be involved in nondegradation signal events, e.g. DNA repair (18). Although key steps of Ub activation and transfer to a substrate as well as the structure of the 20 S subunit of the proteasome are known, the question of how proteins are targeted to the proteasome remains unanswered. It is not known whether there is an additional mechanism to regulate the time point of degradation. One possibility is that monoubiquitination leads to a "point of no return," which proceeds to substrate destruction in a defined time span. Recently, several groups reported that proteins containing a UBA motif can bind directly to Ub and/or poly-Ub, leading to an inhibition of the degradation of target substrates through the proteasome (19-21). UBA domains are a common motif in a variety of protein families involved in protein degradation, cell cycle control, or DNA repair. In vitro and in vivo assays have revealed that UBA domains of the DNA damage-inducible proteins, RAD23 (as well as the fission yeast homolog Rhp23) and DDI1, as well as proteins with no function in DNA repair, p62 and Mud1p, interact specifically with Ub. Furthermore, poly-Ub chain formation is inhibited by RAD23 in vitro in a concentration-dependent manner (20, 21). In addition, it has been shown that poly-Ub chain extension stops at a length of three Ub moieties and that the inhibition of chain extension of RAD23 is specific for Lys48-linked chains (22). As a consequence of the UBA-Ub interaction, Clarke et al. (23) concluded that RAD23 and DDI1 are involved in the checkpoint control of the cell cycle. They proposed that the UBA domains of RAD23 and DDI1 could bind to the nascent poly-Ub chain of the Pds1 substrate, inhibiting chain extension and thereby increasing the lifetime of Pds1, which would otherwise be rapidly degraded. In this report we provide a structural basis for the interaction of UBA domains with monomeric Ub based on an NMR chemical shift mapping study as well as Ub mutagenesis. Ub binds specifically to both UBA(1) and UBA(2) of HHR23A. The binding interface for both UBA domains is almost identical despite the low overall sequence similarity. Both UBA domains bind to the same region of monomeric Ub, which is also involved in the binding to the proteasome subunit S5a. Models for the UBA(1)-Ub and UBA(2)-Ub complexes were generated from the chemical shift mapping data by de novo docking as well as by homology modeling with the solution structure of the closely related CUE domain in complex with Ub (24). These models revealed very different orientations of the UBA domains on the surface of Ub. Our results suggest that UBA domains may interact with Ub as well as other proteins, e.g. the HHR23A-binding proteins HIV-1 Vpr (25), methyladenine DNA glycosylase (26), p300/cyclic AMP-responsive element-binding protein (27), and peptide:N-glycanase (Png1) (28), in more than one way while utilizing the same binding surface.
Preparation of ProteinsUBA(1) and UBA(2) of HHR23A were prepared as previously described (29, 30). 15N- and 15N-,13C-labeled proteins were prepared by growing cells on M9 minimal medium using 15NH4Cl and 13C6-glucose as the sole nitrogen and carbon sources. For the titration of the UBA domain proteins, bovine Ub (amino acid sequence is identical to human Ub) was purchased from Sigma and purified by gel filtration. Purity was subsequently checked by SDS-PAGE and analytical reversed phase HPLC. For the NMR chemical shift mapping of Ub, yeast Ub was prepared from the expression plasmid pET11c-Ub (gift from A. Varshavsky). The protein was purified by reversed phase HPLC. Alanine mutants of yeast Ub were generated by site-directed mutagenesis using QuikChange (Stratagene). The cDNA sequence was verified by sequencing using the DyeDeoxy-Terminator method (PerkinElmer). Ub mutant proteins were prepared similar to wild type Ub. Protein concentrations were determined using an extinction coefficient at A280 of 1280 M-1 for Ub and UBA(2) and 5120 M-1 cm-1 for UBA(1), based on their respective amino acid compositions.
Chemical Shift Mapping ExperimentsProteins for titration experiments were dialyzed against identical buffer (50 mM sodium phosphate, pH 6.5, 100 mM sodium chloride, 2 mM deuterated dithiothreitol (Cambridge Isotope Laboratories CIL) to avoid chemical shift changes because of differences in buffer conditions. NMR samples for titration studies contained 0.25-0.5 mM 15N-labeled protein. Unlabeled protein was added stepwise up to a final ratio of 1:10. Unlabeled Ub was concentrated as far as possible to minimize volume changes throughout the titration. The "reverse" titration of Ub with UBA proteins was performed similar to the experiment described above. In addition, the UBA domains were tested for homo- and heterodimerization with NMR mapping experiments. 15N-Labeled UBA(1) (0.7 mM) was mixed with unlabeled UBA(2). All of the measurements were performed at 27 °C using a Bruker DRX-500 or DRX-600. Two-dimensional 1H-15N HSQC experiments with Watergate water suppression (60) and water flip-back pulses were used for monitoring the chemical shift changes. All of the titration data sets were processed identical using the software XWINNMR (Bruker). 1H dimensions were referenced using external 3-trimethyl-2,2,3,3-tetradeuteropropionian. 15N nitrogen and 13C carbon frequencies were referenced using the gyromagnetic ratios 15N/1H = 0.101329118 and 13C/1H = 0.251449530. The changes of proton and nitrogen chemical shifts were averaged based on their gyromagnetic ratios using the following equation.
The chemical shifts of yeast Ub were reassigned (chemical shifts published for human Ub differ from yeast Ub because of three amino acid changes) using 15N,13C-labeled yeast Ub. A set of triple-resonance experiments (31) (CBCA(CO)NH, CBCANH, HBHA(CO)NH, H(C)(CO)NH-total correlation spectroscopy, CC(CO)NH-total correlation spectroscopy, and HCCH-correlation spectroscopy) was acquired to assign all backbone and side chain chemical shifts for yeast Ub under the conditions used for the NMR chemical shift perturbation studies.
Homology Modeling of the UBA-Ubiquitin ComplexesModel complexes of the UBA(1)/UBA(2) with Ub were built using the structure of the CUE domain of Cue2 and Ub as a template (24). Initial structures of the complexes were obtained by fitting the C De Novo Docking for the UBA-Ubiquitin InteractionThe program HADDOCK (32) was applied to dock Ub and the UBA domains of HHR23A. In this approach, the data obtained from chemical shift mapping are transformed into a set of ambiguous distance restraints that are used together with geometrical and electrostatic complementarity to dock the two molecules. Chemical shift mapping of UBA(1), UBA(2), and Ub and surface accessibility data calculated with the program NACCESS (33) were used to define ambiguous distance restraints as described. The target distance of these restraints was set to 2.0 Å; force constants for empirical and experimental restraints were used as suggested by the default settings. Initially 750 structures for a Ub-UBA complex were generated by docking Ub and UBA(1)/UBA(2) as rigid bodies using only ambiguous distance restraints, van der Waals' energy, and electrostatic energy terms using the program suite ARIA1.2 and CNS (34, 35). Of those, 200 structures with the lowest overall energy were subsequently refined, allowing the side chain conformations of residues within the binding interface to be flexible. Finally, 100 refined structures with lowest overall energy were then refined using a short molecular dynamics simulation in explicit solvent to allow for a correct implementation of the electrostatic contribution. The method was applied to two closely related complexes for which structures were available, the UIM-2 S5a-HHR23A Ubl and the Cue2 CUE domain-Ub complex, to test for the reliability of that procedure. For the simulation of the Cue2 CUE domain-Ub interaction, the distance restraint set was defined based on the residues buried in the complex structure; for the simulation of the HHR23A Ubl-UIM-2 S5a interaction the same criteria for chemical shift changes and residue accessibility as used for the Ub-UBA domain modeling were applied for generating the distance restraint set. The test calculations yielded model structures that were within 1.5 Å rmsd of the experimentally determined structures. Different cut-off values for the surface accessibility in the definition of the distance restraint set used for the computational docking did not lead to altered complex architectures, confirming that the computational results are stable. The influence of the template structures of the individual complex components was determined by using the Ub structure of the Cue2 CUE-Ub complex for docking to the UBA domains as well as the structure of free Ub (Protein Data Bank entry code 1UBI [PDB] ). Again no change in complex architecture could be observed.
Ubiquitin Binds Specifically to HHR23A UBA(1) and UBA(2) DomainsHigh resolution structures of both UBA domains of HHR23A have been determined showing that overall the three-helical bundles are largely identical, despite a low ( 20%) sequence identity (29, 30, 36). The main structural differences are slightly different packing of side chains in the hydrophobic core and the conformation of the N and C termini. Both domains have an unusually large hydrophobic surface patch that was predicted to be a protein-protein interface (29, 30). Based on the structures and sequence analysis of UBA domains, we predicted that a conserved portion of this hydrophobic patch would be the binding region for Ub (29). To map the binding site of Ub on UBA domains, NMR chemical shift perturbation experiments were performed on both UBA(1) and UBA(2) from HHR23A. Samples of 15N-labeled UBA(1) and UBA(2) were mixed with bovine Ub, and two-dimensional HSQC spectra were acquired to monitor the changes in the chemical shifts of the backbone amides induced by the binding to Ub (Fig. 1, a and b). Changes in chemical shifts are observed for several residues in both UBA(1) and UBA(2) starting at about 0.6 molar equivalents, indicating the formation of specific complexes in fast exchange on both the 500 and 600 MHz NMR time scale. Assuming a binary interaction between the UBA domains and Ub, a KD in the range of 500-600 µM can be determined from nonlinear fitting of the titrations with either UBA(1) or UBA(2) (see also Fig. 3, a and b). However, a binding constant of about 10 µM was reported for full-length Rad23 (19), which may indicate cooperative binding of both UBA domains (37). Plots of the chemical shift changes versus residue (Fig. 1, c and d) show that the absolute magnitudes of the chemical shift changes are similar for UBA(1) and UBA(2).
A recent chemical shift mapping study of the closely related HHR23B reported that the Ubl and UBA domains interact with each other with a reported KD of 2 mM, which is 10-fold higher than their calculated KD for the HHR23B-UBA interactions of 300 µM (38). The Ubl domain is structurally very similar to Ub (39, 40) and exhibits a high sequence identity. Weak binding between the HHR23A Ubl and UBA domains was also detected in the context of the full-length protein as well as isolated domains (41) under slightly different buffer conditions and at higher field strength, with a binding affinity at least 10-fold lower compared with the already weak interaction between the UBA domains and Ub. Ubiquitin Binds to a Conserved Surface Epitope of UBA DomainsThe results of the chemical shift perturbation study were mapped onto the structures of UBA(1) and UBA(2) (Fig. 1). Although UBA(1) and UBA(2) have a relatively low sequence identity, the binding interfaces revealed by the chemical shift mapping are remarkably similar. Residues exhibiting the largest chemical shift changes upon binding to Ub, Gly174 and His192 for UBA(1) and Leu330, Gly331, and Glu348 for UBA(2), respectively, are in similar locations. A large cluster of residues in UBA(1) involved in binding is located at the C-terminal end of the first helix, Ile170, Met171, Ser172, and the first three residues in the short and highly conserved loop 1, Met173, Gly174, and Tyr175 (Fig. 1, a, c, and e). Several positions on the third helix of UBA(1) also show large changes upon binding, specifically His192 and Arg193 at the N terminus of helix 3 and Tyr197 at the second turn of helix 3. Together, these residues form a consecutive patch of about 520 Å on the surface of UBA(1) (Fig. 1e). The residues with the largest chemical shift changes, His192, Met173, and Gly174, are in the center of the epitope. Residue Tyr197 also exhibits a significant change in chemical shift upon the addition of Ub but is located on the "back face" of the binding interface. The change for the amide proton and nitrogen frequencies of Tyr197 might be attributed to small structural changes in the hydrophobic core because of interactions of the side chain of Tyr197 with Arg193, which is part of the binding interface. Residues in the C terminus of UBA(1) do not exhibit chemical shift changes when Ub is added (Fig. 1c) and do not seem to be part of the binding interface. This was surprising, because the C terminus is close to the hydrophobic patch involved in binding and is of relatively rigid nature (29).
Analysis of the binding of UBA(2) to Ub reveals that the general location of the epitope remains the same compared with UBA(1) (Fig. 1). Identical positions in the C-terminal end of helix 1 (Arg326, Leu327, and Ala329) as well as the first three residues of the hydrophobic loop 1 (Leu330, Gly331, and Phe332) form a large part of the binding interface to Ub. In addition, residues at similar positions at the N terminus of helix 3 (Glu348 and Asn349) are also among the residues with the highest chemical shift changes upon binding. The structurally equivalent residue of Glu348 in UBA(1) is Pro191; therefore no information about whether or not Pro191 is involved in the binding of Ub can be deduced from this analysis. A larger number of residues at the C-terminal end of helix 3 of UBA(2) are involved in binding than observed for UBA(1). On the last turn of helix 3, residue Leu356 is affected by the binding, whereas residues of UBA(1), Leu199 and Thr200, do not exhibit significant changes. Overall, with a binding area of
UBA Domains Bind on the Five-stranded
Chemical shift mapping of the amide resonances of Ub as a function of added UBA(1) or UBA(2) up to a ratio of 1:10 Ub:UBA again revealed complexes in fast exchange on the NMR timescale. Upon the addition of UBA(1) to Ub, significant changes in chemical shift are observed for 23 residues (Fig. 2, a, c, and e). Most of these residues are located on the
These results suggest that both UBA domains are recognized by and bind to Ub using a similar binding epitope. The binding epitopes of UBA(1) and UBA(2) on Ub overlap perfectly, and the differences in several of the amino acids, e.g. His192 versus Glu348 and Arg193 versus Asn349, might account for the chemical shift differences observed for the titration of Ub with either UBA(1) or UBA(2). However, another possibility is that the two UBA domains are oriented differently on Ub, as discussed below. Differential Chemical Shift Mapping of UBA DomainsTo further address the question of how the UBA domains of HHR23A bind to Ub, especially whether the binding mechanism differs between the two UBA domains, we tried differential chemical shift mapping (42). Five Ub mutants (L8A, R42A, K48A, H68A, and R72A) with single amino acid substitutions were chosen for study. All of the mutated residues have side chains oriented toward the binding interface and show significant chemical shift changes upon binding to the UBA domains. Differences in the chemical shift changes between a Ub mutant and wild type Ub should be observed if the residue pair is in close proximity in the UBA-Ub complex. The binding affinities were determined by nonlinear fitting of the chemical shift changes of three residues located within the binding epitope (Ile170, Gly174, and His192 for UBA(1); Gly331, Glu348, and Leu356 for UBA(2)). Surprisingly, the mutant Ub L8A exhibits much smaller absolute chemical shift changes upon binding to UBA(1) and UBA(2), suggesting a lower binding affinity. However, quantitative analysis of the binding curves (Fig. 3, a and b) shows that the affinity of UBA(1) and Ub L8A is not changed significantly (490 and 570 µM, respectively). Similar small changes in KD (less than 2-fold) were observed for the other Ub mutants (K48A, H68A, and R72A). Interestingly, the binding affinities of the Ub mutant R42A for UBA(1) and UBA(2) seem to be increased 3.5- and 2-fold, respectively. Comparison of the UBA-Ub titration studies with those of the Ub mutants showed no significant differences in the course of the chemical shift changes for UBA(2). In contrast, for UBA(1), the chemical shift pattern of Ub mutants L8A and R42A is different compared with wild type Ub. For L8A significant changes are observed for residues Glu169 in helix 1 as well as Arg193, Ala194, and Glu196 in helix 3 (Fig. 3c), indicating that these residues are in close proximity to each other in the complex. For Ub mutant R42A, the largest differences are observed for residues Glu169 (helix 1) and Ala194, Glu196, Leu199 and Thr200 located in helix 3 (Fig. 3d). Because the side chains of Leu8 and Arg42 are separated by only 6 Å, similar residues can be affected by both Ub mutants. Model building of a UBA(1)-Ub complex based on just the two identified residue pairs is, however, not possible because two interaction points do not define the orientation of two rigid bodies unambiguously. These mutagenesis results are consistent with an analysis of the interaction of the CUE domain of Cue2 protein using two Ub mutants K48A and H68A (24). Although both mutations are located in the center of the interacting patch, they did not influence the interaction significantly. Apparently, the binding specificity between Ub and its interaction partners is achieved by the overall surface topology, which explains why a single point mutation is not able to destabilize the binding substantially.
Homology Modeling of the UBA-Ub InteractionRecently, the structures of a UBA homolog, the CUE domain in complex with Ub, was determined by NMR and x-ray crystallography (24, 43). Based on the structure of the CUE domain of Cue2 bound to Ub (24), we generated a model for the interaction of UBA(1) and UBA(2) with Ub (Fig. 4). The solution structure of the Cue2-Ub complex instead of the crystal structure of the CUE domain of Vps9 bound to Ub (43) was used for model building because of the higher similarity between the CUE domain structure of Cue2 and the UBA domains of HHR23A. The crystal structure of the CUE domain of Vps9p forms a domain-swapped intertwined dimer, and binding to Ub results in a large conformational change. Furthermore, the binding affinity of the Cue2 CUE domain to Ub (KD =
In the UBA-Ub models the axis of helix 1 is oriented at an angle of 45° to the -strand 5 of Ub, helix 3 is running at an angle of 50° across the 5, and the helical axis of helix 2 is running almost in parallel with 5 (Fig. 4a). Hydrophobic residues in loop 1 (Gly174-Tyr175 UBA(1) and Gly331-Phe332 UBA(2)) of the UBA domain interact with the 3- 4 loop residues of Ub (Fig. 4a). The N terminus of helix 3 of the UBA domains (Pro191-His192 UBA(1) and Glu348-Asn349 UBA(2)) is close to the C terminus of -strand 5 and the 1- 2 loop of Ub. The important "hydrophobic triad" of Ub, residues Leu8, Ile44, and Val70, is buried in the interface. In the Ub-UBA(1) complex, Leu8 of Ub is surrounded by UBA(1) residues Glu169, Ile170 (helix 1), Pro191, and His192 (loop 2) (Fig. 4a). In the complex of Ub-UBA(2), residue Leu8 of Ub has van der Waals' contacts to residues of UBA(2) occupying similar positions in the three-helical bundle (Ala323, Leu327 (helix 1), and Glu348 (loop 2)) (Fig. 4b). Ile44 of Ub is packed against Met173 (helix 1), Tyr175 (loop 1), and Val195 (helix 3) of UBA(1); in Ub-UBA(2) residues Leu330 (helix 1), Ala352, and Leu356 (helix 3) of UBA(2) are close to Ile44 of Ub. Val70 of Ub contacts either His192 (loop 2) of UBA(1) or residues Ala323, Leu327 (helix 1), and Glu348 (loop 2) of UBA(2).
The differences in the environment of Ub Val70 bound to UBA(1) versus UBA(2) are due to slightly different interhelical angles of the three-helical bundle, which lead to helix 1 of UBA(1) pointing further away from the binding interface. Although different residues of both UBA domains interact with Leu8, Ile44, and Val70 of Ub, the hydrophobic nature in the contact area is preserved. Only a very small number of possible intermolecular hydrogen bonds can be identified in the binding interfaces of the model complexes: five for Ub-UBA(1) (Ub-UBA: Lys6-Glu169, Gly47-Tyr175, Gly47-Met173, His68-Glu169, and Leu71-His192) and two for Ub-UBA(2) (Ub-UBA: Thr7-Arg326 and His68-Ala329). Similarly, in the complex of Cue2 CUE domain and Ub, only two hydrogen bonds are observed. The dearth of polar interactions in addition to the relatively small interface area ( De Novo Docking of the Ubiquitin-UBA Domain ComplexesIn addition to the homology modeling, we performed a de novo docking using the program HADDOCK and employing data from our chemical shift mapping (32). The chemical shift changes of UBA(1), UBA(2), and Ub were used together with surface accessibility data to define ambiguous distance restraints between the two molecules. The method relies on geometrical and electrostatic complementarity of the binding epitopes of the interacting molecules, and it has been applied successfully to other complexes (32). Much to our surprise, the results of the de novo docking of the Ub-UBA domain complexes revealed completely different complex structures compared with those obtained by homology modeling (Figs. 4 and 5). Additionally, the architecture of the complexes resulting from de novo docking for the Ub-UBA(1) and Ub-UBA(2) interaction (Fig. 5) is also different, indicating that the UBA domains of HHR23A might interact with Ub differently.
A structural alignment of the Ub molecules of the two models of UBA(2)-Ub (homology model versus de novo docking) shows the large difference. In the complex of UBA(2)-Ub obtained by de novo docking, helices 1 and 3 of UBA(2) run almost parallel to -strand 5 of Ub, and UBA(2) helix 2 and Ub -strand 5 are oriented at an angle of about 20° (Fig. 5b) relative to each other. In the homology model (and hence for the template complex Cue2 CUE domain-Ub), helices 1 and 3 of the UBA domain and -strand 5 of Ub share an angle of almost 45°, and the helical axis of helix 2 of UBA(2) runs parallel to -strand 5 of Ub (Fig. 4b). Consequently, the de novo docking model can be transformed into the homology model by a rotation of about 45° counterclockwise. A large positional movement is observed for the residues in the Gly-Phe-Pro loop of UBA(2), with the C positions of Leu330, Gly331, and Phe332 differing by 7-8 Å between the two modeling approaches. The positional change of the residues in loop 2 (Glu348 and Asn349) and helix 3 (Leu356) with respect to the Ub binding interface is smaller (distances for C positions: Glu348, 3 Å; Asn349, 4 Å; and Leu356, 4 Å). Consequently a different residue pairing is observed in the interface of the de novo docking model. Leu8 of Ub is surrounded by Arg326, Leu327, Leu330, and Glu348 of UBA(2), Ile44 of Ub is in van der Waals' contact with residues Ala352, Asn353, and Leu356 of UBA(2), and Val70 of Ub is in close proximity to Asn349 and Ala352 of UBA(2). In the complex of the Cue2-CUE domain and Ub (and thus in the homology model of UBA(2)Ub), the residues located at the N terminus of the helix 1 are part of the binding interface, which is not the case for the de novo docking. In the homology model, Ala323 and Arg326 both contact residues Leu8 and Thr9 of Ub, whereas for the de novo docking model, Ala323 shares no contacts with any residues of Ub.
De novo docking of UBA(1) and Ub resulted in a model complex that differs from the homology model to an even greater extent. The helical axes of helices 1, 2, and 3 of UBA(1) in the two models differ by about 55°, 70°, and 100°, respectively. To transform both models into each other, a rotation of almost 90° is required. In the UBA(2)-Ub models the largest positional discrepancies between the two different models are observed for the Gly-Phe-Pro loop, whereas for the UBA(1)-Ub models the major positional changes are for helix 3. Residue His192 is located on top of the
Because the results of the de novo docking for UBA(1)-Ub and UBA(2)-Ub deviate from the models obtained by homology modeling, we tested whether the differences might be an artifact of the docking procedure. Because all 100 refined structures of both complexes cluster into one single structure family exhibiting an rmsd within the cluster of less than 0.8 Å, the observed complex architecture obtained by the docking simulation seems to be very stable. Different sets of distance restraints for the docking did not lead to altered complex architectures. We also tested the docking procedure on two other complexes, the Cue2 CUE domain-Ub interaction (24) and the complex of HHR23A Ubl domain and UIM-2 of S5a (39). For both simulations, structure data, and for the latter, chemical shift mapping data, were available. The docking of the Cue2 CUE domain and Ub reproduced the experimental structure of Kang et al. (24) with an rmsd of less than 1 Å for the C
Is the Hydrophobic Knob a General Ubiquitin-binding Site for UBA Domains?The chemical shift mapping study presented here reveals a structural basis for the recognition and binding of UBA domains of HHR23A to Ub. The hydrophobic surface patches that were predicted to be the site of protein-protein interactions for the UBA domains (29) comprise a large portion of the binding epitope for Ub. In addition, some charged and polar residues appear to be important based on the chemical shift mapping. These results agree very well with a mapping study of HHR23B UBA domains published very recently by Choi and co-workers (38), except for the equivalent residue to HHR23A UBA(1) Tyr197. Despite the relatively low binding affinity of 500 µM, the interaction of the UBA domains with Ub is specific. Furthermore, binding between HHR23A/B UBA and Ubl domains is an order of magnitude weaker (38, 41).
The amino acid sequence identity of UBA domains is about 25%, which is too low to propose a common biological function or even a common binding mechanism. However, if the analysis is limited just to the binding interface determined in this study, the degree of similarity of the binding epitope is as high as for the hydrophobic core of the UBA domains (Fig. 6). The binding epitope of Ub on UBA(1) and UBA(2) is comprised of 9 and 12 residues, respectively. The high degree of conservation for the residues located on the surface (similarity for binding region >50%, similarity within hydrophobic core
The Five-stranded -Sheet of Ubiquitin Is a Universal Binding Site for Ubiquitin-interacting ProteinsDespite the differences in the binding epitopes determined for the UBA domains of HHR23A, the position of the interface on Ub is practically identical for both UBA domains. This region has been identified to interact with several other protein domains, indicating that the hydrophobic surface patch on the five-stranded -sheet of Ub is probably a general protein-protein interface that can facilitate various interactions. At least six ubiquitin-interacting domains have been described so far (44). The structures of five domains have been determined, and the site of their interaction with Ub has been mapped. The three-dimensional structures of these domains, UEV (ubiquitin E2 enzyme variant) (45), NZF (novel zinc finger domain) (46), UIM (ubiquitin-interacting motif) (39, 47, 48), UBA (29, 30, 36), and CUE (24, 43), vary greatly in architecture and size, being either a single helix (UIM), a pure -strand structure (NFZ), three helical bundles (CUE and UBA), or a mixed - structure (UEV). Despite this large variety, all domains interact with Ub via the same hydrophobic patch on the five-stranded -sheet of Ub. Comparing the binding to Ub for all domains reveals only very minor differences in the location of the binding sites. The center of the binding site, the hydrophobic patch around residue Ile44 of Ub, seems to be identical for all ubiquitin-interacting domains so far, although residues close to that patch have been mapped to different biological functions (49). Very little structural data for the interaction of monomeric Ub or ubiquitin-like domains in complex with an ubiquitin-interacting domain are available so far (24, 39, 43), probably because of the low to moderate binding affinities for such interactions (KD = 10-500 µM) (44). A comparison of the binding mechanism of a single -helix with an Ub-like domain (complex of HHR23A Ubl-UIM-2 of S5a) with that of a three-helical bundle bound to Ub (complexes of Cue2 CUE domain and Vps9p CUE with Ub) shows that complexes of Ub and Ub-interacting domains can adopt different architectures. The single -helix of the UIM motif of S5a binds on top of -strand 5 of HHR23A Ubl, with the axes of the -helix of the UIM motif and of -strand 5 of Ub running anti-parallel (39). In contrast, the three-helical bundle of the CUE domain of Cue2 binds via helices 1 and 3, which run across -strands 1, 3, and 5 at an angle of about 30° (24). In the case of the CUE domain of Vps9p, the interaction is even more complex. In the x-ray crystal structure of the complex of the CUE domain of Vps9p with Ub, the CUE domain forms a domain-swapped dimer in which helices 1 and 3 interact with Ub in a similar manner to that observed for the Cue2 CUE domain-Ub complex, but here helix 2 has additional contacts with residues of Ub (43). Despite the differences in the architecture of the complexes, the interacting amino acids are conserved, with only hydrophobic amino acids taking part in the interaction. The absence of hydrogen bonds in the center of the interface probably explains the limited specificity of Ub, because there is no requirement to maintain the geometry of hydrogen bonding acceptors or donors, and therefore specificity for a binding partner is only generated by geometrical restrictions for the interacting hydrophobic side chains. It is also interesting to note that the binding affinity between Ub and the Ub-interacting domains of known structures correlates with the size of the interface. The Cue2 CUE domain Ub interface measures 450 Å2 and has a binding affinity of about 150 µM (24), whereas the S5a UIM-2 HHR23A Ubl complex has a buried surface area of roughly 600 Å2 (39) and a KD of 10 µM. The interface between the CUE domain of Vps9p and Ub measures 520 Å2 (because of additional interacting residues in the second helix) resulting in an affinity of 20 µM (43).
A Model for the Interaction of UBA Domains and UbiquitinOne very surprising result of the modeling of the UBA-Ub interaction presented in this study is the large differences between the models obtained by homology modeling and de novo docking using the program HADDOCK. The homology models were built using the NMR structure of the Cue2 CUE-ubiquitin complex, with the CUE domain being replaced by the UBA domains of HHR23A. Residues of the Cue2 CUE domain that interact with residues of the five-stranded We note, in addition, that the differential chemical shift mapping employing Ub mutants showed differences for UBA(1) and UBA(2). These might indicate that both UBA domains bind to Ub by a different binding mechanism, supporting the results obtained by the de novo docking procedure. Analysis of the distribution of charged residues surrounding the hydrophobic patch, which is in the center of the binding interface for CUE and UBA domains, clearly shows that the electrostatic potentials are distinct for these domains. Biological data suggest differences in functions of the CUE and UBA domain (50-52), and sequence comparison clearly distinguishes between the CUE and UBA family (53). All of the CUE domains identified so far bind to monomeric Ub with low to moderate affinity (20-160 µM), and some but not all CUE domains seem to bind to poly-Ub in addition (52). It was not reported whether their binding affinity for poly-Ub is higher than for monomeric Ub; however, data from Shih et al. (52) suggest that at least the CUE domain of Vps9 has no preference for long poly-Ub chains over short poly-Ub chains. This binding preference is probably required for the maintenance of monoubiquitination (52), which is in turn a signal for trafficking and receptor endocytosis. However, for a more quantitative analysis the binding constants of several CUE domains for mono- and poly-Ub have to be determined. The biological function of the UBA domain is, on the other hand, still in debate (20, 22, 23, 37, 54). Several groups have reported that UBA domains bind to monomeric Ub as well as to poly-Ub, although the data are contradictory in some cases (19, 54). Binding to monomeric Ub was associated with inhibition of further extension of the nascent Ub chain, which results in an inhibition of the degradation of a substrate (20-23). The binding to poly-Ub was explained with a possible shuttle function for the transport of a substrate to the proteasome (37, 54-57). Although the physiological "Ub target" of the UBA domains of RAD23 and other proteins has not been determined, all in vitro binding experiments have yielded an at least 1000-fold greater affinity of the UBA domains for tetra-Ub compared with mono-Ub (22, 54). Therefore it is very likely that in the presence of polyubiquitinated substrates, the main binding partner of UBA domains might be these poly-Ub chains rather than monoubiquitinated substrates or free Ub, unless the concentration of polyubiquitinated substrates is very low. The molecular nature of the tighter binding of UBA to poly-Ub is not yet clear; however, structure analysis (58) as well as mutagenesis data on poly-Ub (59) suggest that they do not form plain linear poly-protein chains like pearls on a string but rather adopt globular structures. Hence additional contacts between the linked Ub moieties and possibly other epitope(s) of the UBA domain might lead to the increase in affinity in comparison with monomeric Ub. Because UBA domains are often associated together with Ub-like domains in modular proteins, the UBA domains could act as poly-Ub "receptors," whereas the Ubl domain might interact directly with the proteasome. However, such a shuttle mechanism has not yet been confirmed in vivo. The differences for UBA and CUE domains in their biological function as well as in the probable binding target, monomeric Ub versus poly-Ub, make it plausible that the binding mechanisms of CUE and UBA domains do not need to be identical.
Interestingly, our de novo model was recently at least partially confirmed by results by Walters et al. (41). Based on chemical shift differences between the isolated Ubl and UBA domains and the domains in the context of the full-length HHR23A protein, it was concluded that in full-length HHR23A the Ubl domain interacts in a dynamic fashion with the individual UBA domains, and this interaction has a 1:1 stoichiometry, i.e. Ubl is exchanging between one or the other UBA domains. No interdomain NOEs were observed, consistent with very weak binding. Using residual dipolar couplings measurements, Walters et al. (41) defined the relative orientation of the domains in the modular protein. Although the coordinates for their models of the interaction between Ubl and the UBA domains are not available, using the figures of their publication we find a striking similarity between their models for Ubl-UBA interaction and our Ub-UBA de novo models. Similar to the interaction with Ubl, an identical surface epitope of Ub is involved in the binding to both UBA(1) and UBA(2) domains. Helix 1 of either UBA(1) (residues 191-199) or UBA(2) (residues 348-356) contacts the five-stranded Finally, we note that RAD23 has been described as a binding partner through its UBA(2) domain for various proteins involved in DNA repair and cell cycle control, e.g. HIV-1 Vpr (25), methyladenine DNA glycosylase (26), p300/cyclic AMP-responsive element-binding protein (27), and Png1 (28). Because the UBA-binding epitopes of these proteins exhibit very different structures, the UBA domain of RAD23 must be able to bind to various structural architectures. Therefore the binding mechanism of UBA domains is also likely to vary for the diverse interacting proteins.
* This work was supported by National Institutes of Health Grants AI43190 (to I. S. Y. C. and J. F.). The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
1 The abbreviations used are: Ub, ubiquitin; UBA, ubiquitin-associated; E1, ubiquitin-activating enzyme; E2, ubiquitin carrier protein; E3, ubiquitin-protein isopeptide ligase; HPLC, high pressure liquid chromatography; rmsd, root mean square deviation; HSQC, heteronuclear single quantum correlation.
We thank Dr. Alexander Varshavsky for the gift of yeast ubiquitin expression plasmid pET11c-Ub; Dr. Carina Johansson, Darian Cash, and Nathan Cho for performing some of the NMR experiments; and Evan Feinstein for manuscript and figure preparation.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||