Allosteric regulation of E2:E3 interactions promote a processive ubiquitination machine

RING finger proteins constitute the large majority of ubiquitin ligases (E3s) and function by interacting with ubiquitin‐conjugating enzymes (E2s) charged with ubiquitin. How low‐affinity RING–E2 interactions result in highly processive substrate ubiquitination is largely unknown. The RING E3, gp78, represents an excellent model to study this process. gp78 includes a high‐affinity secondary binding region for its cognate E2, Ube2g2, the G2BR. The G2BR allosterically enhances RING:Ube2g2 binding and ubiquitination. Structural analysis of the RING:Ube2g2:G2BR complex reveals that a G2BR‐induced conformational effect at the RING:Ube2g2 interface is necessary for enhanced binding of RING to Ube2g2 or Ube2g2 conjugated to Ub. This conformational effect and a key ternary interaction with conjugated ubiquitin are required for ubiquitin transfer. Moreover, RING:Ube2g2 binding induces a second allosteric effect, disrupting Ube2g2:G2BR contacts, decreasing affinity and facilitating E2 exchange. Thus, gp78 is a ubiquitination machine where multiple E2‐binding sites coordinately facilitate processive ubiquitination.     

Differential Ubiquitin Binding by the Acidic Loops of Ube2g1 and Ube2r1 Enzymes Distinguishes Their Lys-48-ubiquitylation Activities

The ubiquitin E2 enzymes, Ube2g1 and Ube2r1, are able to synthesize Lys-48-linked polyubiquitins without an E3 ligase but how that is accomplished has been unclear. Although both E2s contain essential acidic loops, only Ube2r1 requires an additional C-terminal extension (184–196) for efficient Lys-48-ubiquitylation activity. The presence of Tyr-102 and Tyr-104 in the Ube2g1 acidic loop enhanced both ubiquitin binding and Lys-48-ubiquitylation and distinguished Ube2g1 from the otherwise similar truncated Ube2r1^1–183 (Ube2r1C). Replacement of Gln-105–Ser-106–Gly-107 in the acidic loop of Ube2r1C (Ube2r1C^YGY) by the corresponding residues from Ube2g1 (Tyr-102–Gly-103–Tyr-104) increased Lys-48-ubiquitylation activity and ubiquitin binding. Two E2∼UB thioester mimics (oxyester and disulfide) were prepared to characterize the ubiquitin binding activity of the acidic loop. The oxyester but not the disulfide derivative was found to be a functional equivalent of the E2∼UB thioester. The ubiquitin moiety of the Ube2r1C^C93S-[¹⁵N]UB^K48R oxyester displayed two-state conformational exchange, whereas the Ube2r1C^C93S/YGY-[¹⁵N]UB^K48R oxyester showed predominantly one state. Together with NMR studies that compared UB^K48R oxyesters of the wild-type and the acidic loop mutant (Y102G/Y104G) forms of Ube2g1, in vitro ubiquitylation assays with various mutation forms of the E2s revealed how the intramolecular interaction between the acidic loop and the attached donor ubiquitin regulates Lys-48-ubiquitylation activity.     

Solution structure and dynamics of human ubiquitin conjugating enzyme Ube2g2

Ube2g2 is an E2 enzyme which functions as part of the endoplasmic reticulum‐associated degradation (ERAD) pathway responsible for identification and degradation of misfolded proteins in the endoplasmic reticulum. In tandem with a cognate E3 ligase, Ube2g2 assembles K48‐linked polyubiquitin chains and then transfers them to substrate, leading ultimately to proteasomal degradation of the polyubiquitin‐tagged substrate. We report here the solution structure and backbone dynamics of Ube2g2 solved by nuclear magnetic resonance spectroscopy. Although the solution structure agrees well with crystallographic structures for the E2 core, catalytically important loops (encompassing residues 95–107 and 130–135) flanking the active site cysteine are poorly defined. ¹⁵N spin relaxation and residual dipolar coupling analysis directly demonstrates that these two loops are highly dynamic in solution. These results suggest that Ube2g2 requires one or more of its protein partners, such as cognate E3, acceptor ubiquitin substrate or thiolester‐linked donor ubiquitin, to assume its catalytically relevant conformation. Within the NMR structural ensemble, interactions were observed between His94 and the highly mobile loop residues Asp98 and Asp99, supporting a possible role for His94 as a general base activated by the carboxylate side‐chains of Asp98 or Asp99. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Deducing hydration sites of a protein from molecular dynamics simulations

Invariant water molecules that are of structural or functional importance to proteins are detected from their presence in the same location in different crystal structures of the same protein or closely related proteins. In this study we have investigated the location of invariant water molecules from MD simulations of ribonuclease A, HIV1-protease and Hen egg white lysozyme. Snapshots of MD trajectories represent the structure of a dynamic protein molecule in a solvated environment as opposed to the static picture provided by crystallography. The MD results are compared to an analysis on crystal structures. A good correlation is observed between the two methods with more than half the hydration sites identified as invariant from crystal structures featuring as invariant in the MD simulations which include most of the functionally or structurally important residues. It is also seen that the propensities of occupying the various hydration sites on a protein for structures obtained from MD and crystallographic studies are different. In general MD simulations can be used to predict invariant hydration sites when there is a paucity of crystallographic data or to complement crystallographic results.     

Comparison of computational model and X-ray crystal structure of human serotonin transporter: potential application for the pharmacology of human monoamine transporters

Abstract

The human serotonin transporter (hSERT) played a significant role in neurological process whose structural basis had been analysed for many years. Recently, the first homology model was constructed for hSERT based on the crystal structure of drosophila melanogaster dopamine transporter was published, and the inhibitory mechanism underlying the binding mode between hSERT and approved antidepressants was substantially investigated by molecular dynamics (MD) simulation. Right after this publication, the X-ray crystallographic structures of hSERT were reported, which provided a good opportunity to reassess the performance of previous simulation. In this study, the analyses of side-chain contact map, stereochemical quality and ligand-binding pocket were firstly conducted, which revealed that the constructed homology model of hSERT could successfully reproduce the reported crystal structure. Secondly, the approved antidepressant escitalopram was docked into the X-ray structure, and its binding pose was consistent with the reported docking pose in the homology model. Finally, MD simulation were performed based on the crystal structure of hSERT, and structural features revealed as critical for escitalopram-hSERT interaction by previous simulation were successfully recaptured. Thus, the newly reported X-ray crystal structure of hSERT was precisely predicted by computational model, which demonstrated its reliability in understanding the pharmacology of other human monoamine transporters whose 3-D structure remained unknown.     

A Correspondence Between Solution-State Dynamics of an Individual Protein and the Sequence and Conformational Diversity of its Family

Conformational ensembles are increasingly recognized as a useful representation to describe fundamental relationships between protein structure, dynamics and function. Here we present an ensemble of ubiquitin in solution that is created by sampling conformational space without experimental information using “Backrub” motions inspired by alternative conformations observed in sub-Angstrom resolution crystal structures. Backrub-generated structures are then selected to produce an ensemble that optimizes agreement with nuclear magnetic resonance (NMR) Residual Dipolar Couplings (RDCs). Using this ensemble, we probe two proposed relationships between properties of protein ensembles: (i) a link between native-state dynamics and the conformational heterogeneity observed in crystal structures, and (ii) a relation between dynamics of an individual protein and the conformational variability explored by its natural family. We show that the Backrub motional mechanism can simultaneously explore protein native-state dynamics measured by RDCs, encompass the conformational variability present in ubiquitin complex structures and facilitate sampling of conformational and sequence variability matching those occurring in the ubiquitin protein family. Our results thus support an overall relation between protein dynamics and conformational changes enabling sequence changes in evolution. More practically, the presented method can be applied to improve protein design predictions by accounting for intrinsic native-state dynamics.     

Computer modeling and molecular dynamics simulations of ligand bound complexes of bovine angiogenin: dinucleotide topology at the active site of RNase a family proteins

We have undertaken the modeling of substrate-bound structures of angiogenin. In our recent study, we modeled the dinucleotide ligand binding to human angiogenin. In the present study, the substrates CpG, UpG, and CpA were docked onto bovine angiogenin. This was achieved by overcoming the problem of an obstruction to the B1 site by the C-terminus and identifying residues that bind to the second base. The modeled complexes retain biochemically important interactions. The docked models were subjected to 1 ns of molecular dynamics, and structures from the simulation were refined by using simulated annealing. Our models explained the enzyme's specificity for both B1 and B2 bases as observed experimentally. The nature of binding of the dinucleotide substrate was compared with that of the mononucleotide product. The models of these complexes were also compared with those obtained earlier with human angiogenin. On the basis of the simulations and annealed structures, we came up with a consensus topology of dinucleotide ligands that binds to human and bovine angiogenins. This dinucleotide conformation can serve as a starting model for ligand-bound complex structures for RNase A family of proteins. We demonstrated this capability by generating the complex structure of CpA bound to eosinophil-derived neurotoxin (EDN) by fitting the consensus topology of CpA to the crystal structure of native EDN.     

Dimeric Ube2g2 simultaneously engages donor and acceptor ubiquitins to form Lys48‐linked ubiquitin chains

Cellular adaptation to proteotoxic stress at the endoplasmic reticulum (ER) depends on Lys48‐linked polyubiquitination by ER‐associated ubiquitin ligases (E3s) and subsequent elimination of ubiquitinated retrotranslocation products by the proteasome. The ER‐associated E3 gp78 ubiquitinates misfolded proteins by transferring preformed Lys48‐linked ubiquitin chains from the cognate E2 Ube2g2 to substrates. Here we demonstrate that Ube2g2 synthesizes linkage specific ubiquitin chains by forming an unprecedented homodimer: The dimerization of Ube2g2, mediated primarily by electrostatic interactions between two Ube2g2s, is also facilitated by the charged ubiquitin molecules. Mutagenesis studies show that Ube2g2 dimerization is required for ER‐associated degradation (ERAD). In addition to E2 dimerization, we show that a highly conserved arginine residue in the donor Ube2g2 senses the presence of an aspartate in the acceptor ubiquitin to position only Lys48 of ubiquitin in proximity to the donor E2 active site. These results reveal an unanticipated mode of E2 self‐association that allows the E2 to effectively engage two ubiquitins to specifically synthesize Lys48‐linked ubiquitin chains.     

The mechanism of linkage-specific ubiquitin chain elongation by a single-subunit E2

Ubiquitin chains of different topologies trigger distinct functional consequences, including protein degradation and reorganization of complexes. The assembly of most ubiquitin chains is promoted by E2s, yet how these enzymes achieve linkage specificity is poorly understood. We have discovered that the K11-specific Ube2S orients the donor ubiquitin through an essential noncovalent interaction that occurs in addition to the thioester bond at the E2 active site. The E2-donor ubiquitin complex transiently recognizes the acceptor ubiquitin, primarily through electrostatic interactions. The recognition of the acceptor ubiquitin surface around Lys11, but not around other lysines, generates a catalytically competent active site, which is composed of residues of both Ube2S and ubiquitin. Our studies suggest that monomeric E2s promote linkage-specific ubiquitin chain formation through substrate-assisted catalysis.     

Ligand-based 3-D pharmacophore generation and molecular docking of mTOR kinase inhibitors

The 3-D structure of the human mTOR kinase domain was modeled based on the crystal structure of PI3Kγ using comparative modeling methods, and the ATP-binding site of mTOR was characterized. The mTOR kinase 3-D model structure is similar to the structure of the PI3Kγ kinase domain, and exhibits great similarity to PI3Kγ at the active site of the kinase. Pharmacophore generation, the docking of mTOR inhibitors, and molecular dynamics (MD) simulations of mTOR–inhibitor docked complexes were carried out in this work. The best pharmacophore model generated from 27 ATP-competitive mTOR inhibitors comprised two hydrogen-bond acceptors, one aromatic ring, and one hydrophobic feature. These 27 inhibitors were docked into the ATP-binding site comprising the DFG motif, and the interactions in each protein–inhibitor complex were characterized. Mapping the pharmacophore model onto the docked inhibitors explained the specificity of the features of the pharmacophore and how they were arranged inside the active site of mTOR kinase. MD studies revealed important structural features, such as the large hydrophobic pocket “HP” and hydrophilic pocket “A1,” and the solvent-exposed hydrophilic pocket “A2” at the active site of mTOR. Our results provide structural models of mTOR–inhibitor complexes and clues that should aid in the design of better mTOR kinase inhibitors.     

Mulan E3 ubiquitin ligase interacts with multiple E2 conjugating enzymes and participates in mitophagy by recruiting GABARAP

Mulan is an E3 ubiquitin ligase embedded in the outer mitochondrial membrane (OMM) with its RING finger facing the cytoplasm and a large domain located in the intermembrane space (IMS). Mulan is known to have an important role in cell growth, cell death, and more recently in mitophagy. The mechanism of its function is poorly understood; but as an E3 ligase it is expected to interact with specific E2 ubiquitin conjugating enzymes and these complexes will bind and ubiquitinate specific substrates. The unique topology of Mulan can provide a direct link of communicating mitochondrial signals to the cytoplasm. Our studies identified four different E2 conjugating enzymes (Ube2E2, Ube2E3, Ube2G2 and Ube2L3) as specific interactors of Mulan. Each of these E2 conjugating enzymes was fused to the RING finger domain of Mulan and used in a modified yeast two-hybrid screen. Several unique interactors for each Mulan–E2 complex were isolated. One such specific interactor of Mulan–Ube2E3 was the GABARAP (GABA_A receptor-associated protein). GABARAP is a member of the Atg8 family of proteins that plays a major role in autophagy/mitophagy. The interaction of GABARAP with Mulan–Ube2E3 required an LC3-interacting region (LIR) located in the RING finger domain of Mulan as well as the presence of Ube2E3. The isolation of four different E2 conjugating enzymes, as specific partners of Mulan E3 ligase, suggests that Mulan is involved in multiple biological pathways. In addition, the interaction of GABARAP with Mulan–Ube2E3 supports the role of Mulan as an important regulator of mitophagy and provides a plausible mechanism for its function in this process.     

The UBA domain of conjugating enzyme Ubc1/Ube2K facilitates assembly of K48/K63‐branched ubiquitin chains

The assembly of a specific polymeric ubiquitin chain on a target protein is a key event in the regulation of numerous cellular processes. Yet, the mechanisms that govern the selective synthesis of particular polyubiquitin signals remain enigmatic. The homologous ubiquitin‐conjugating (E2) enzymes Ubc1 (budding yeast) and Ube2K (mammals) exclusively generate polyubiquitin linked through lysine 48 (K48). Uniquely among E2 enzymes, Ubc1 and Ube2K harbor a ubiquitin‐binding UBA domain with unknown function. We found that this UBA domain preferentially interacts with ubiquitin chains linked through lysine 63 (K63). Based on structural modeling, in vitro ubiquitination experiments, and NMR studies, we propose that the UBA domain aligns Ubc1 with K63‐linked polyubiquitin and facilitates the selective assembly of K48/K63‐branched ubiquitin conjugates. Genetic and proteomics experiments link the activity of the UBA domain, and hence the formation of this unusual ubiquitin chain topology, to the maintenance of cellular proteostasis.     

NMR Spectroscopic and Bioinformatic Analyses of the LTBP1 C-Terminus Reveal a Highly Dynamic Domain Organisation

Proteins from the LTBP/fibrillin family perform key structural and functional roles in connective tissues. LTBP1 forms the large latent complex with TGFβ and its propeptide LAP, and sequesters the latent growth factor to the extracellular matrix. Bioinformatics studies suggest the main structural features of the LTBP1 C-terminus are conserved through evolution. NMR studies were carried out on three overlapping C-terminal fragments of LTBP1, comprising four domains with characterised homologues, cbEGF14, TB3, EGF3 and cbEGF15, and three regions with no homology to known structures. The NMR data reveal that the four domains adopt canonical folds, but largely lack the interdomain interactions observed with homologous fibrillin domains; the exception is the EGF3-cbEGF15 domain pair which has a well-defined interdomain interface. ¹⁵N relaxation studies further demonstrate that the three interdomain regions act as flexible linkers, allowing a wide range of motion between the well-structured domains. This work is consistent with the LTBP1 C-terminus adopting a flexible “knotted rope” structure, which may facilitate cell matrix interactions, and the accessibility to proteases or other factors that could contribute to TGFβ activation.     

The ubiquitin E1 enzyme Ube1 mediates NEDD8 activation under diverse stress conditions

Modification of proteins with ubiquitin and ubiquitin-like molecules is involved in the regulation of almost every biological process. Historically, each conjugation pathway has its unique set of E1, E2 and E3 enzymes that lead to activation and conjugation of their cognate molecules. Here, we present the unexpected finding that under stress conditions, the ubiquitin E1 enzyme Ube1 mediates conjugation of the ubiquitin-like molecule NEDD8. Inhibition of the 26S proteasome, heat shock and oxidative stress cause a global increase in NEDDylation. Surprisingly, this does not depend on the NEDD8 E1-activating enzyme, but rather on Ube1. A common event in the tested stress conditions is the depletion of “free” ubiquitin. A decrease in “free” ubiquitin levels in the absence of additional stress is sufficient to stimulate NEDDylation through Ube1. Further analysis on the NEDD8 proteome shows that the modified NEDDylated proteins are simultaneously ubiquitinated. Mass spectrometry on the complex proteome under stress reveals the existence of mixed chains between NEDD8 and ubiquitin. We further show that NEDDylation of the p53 tumor suppressor upon stress is mediated mainly through Ube1. Our studies reveal an unprecedented interplay between NEDD8 and ubiquitin pathways operating in diverse cellular stress conditions.     

How does a small molecule bind at a cryptic binding site?

Protein-protein interactions (PPIs) are ubiquitous biomolecular processes that are central to virtually all aspects of cellular function. Identifying small molecules that modulate specific disease-related PPIs is a strategy with enormous promise for drug discovery. The design of drugs to disrupt PPIs is challenging, however, because many potential drug-binding sites at PPI interfaces are “cryptic”: When unoccupied by a ligand, cryptic sites are often flat and featureless, and thus not readily recognizable in crystal structures, with the geometric and chemical characteristics of typical small-molecule binding sites only emerging upon ligand binding. The rational design of small molecules to inhibit specific PPIs would benefit from a better understanding of how such molecules bind at PPI interfaces. To this end, we have conducted unbiased, all-atom MD simulations of the binding of four small-molecule inhibitors (SP4206 and three SP4206 analogs) to interleukin 2 (IL2)—which performs its function by forming a PPI with its receptor—without incorporating any prior structural information about the ligands’ binding. In multiple binding events, a small molecule settled into a stable binding pose at the PPI interface of IL2, resulting in a protein–small-molecule binding site and pose virtually identical to that observed in an existing crystal structure of the IL2-SP4206 complex. Binding of the small molecule stabilized the IL2 binding groove, which when the small molecule was not bound emerged only transiently and incompletely. Moreover, free energy perturbation (FEP) calculations successfully distinguished between the native and non-native IL2–small-molecule binding poses found in the simulations, suggesting that binding simulations in combination with FEP may provide an effective tool for identifying cryptic binding sites and determining the binding poses of small molecules designed to disrupt PPI interfaces by binding to such sites.     

Ensemble MD simulations restrained via crystallographic data: Accurate structure leads to accurate dynamics

Currently, the best existing molecular dynamics (MD) force fields cannot accurately reproduce the global free‐energy minimum which realizes the experimental protein structure. As a result, long MD trajectories tend to drift away from the starting coordinates (e.g., crystallographic structures). To address this problem, we have devised a new simulation strategy aimed at protein crystals. An MD simulation of protein crystal is essentially an ensemble simulation involving multiple protein molecules in a crystal unit cell (or a block of unit cells). To ensure that average protein coordinates remain correct during the simulation, we introduced crystallography‐based restraints into the MD protocol. Because these restraints are aimed at the ensemble‐average structure, they have only minimal impact on conformational dynamics of the individual protein molecules. So long as the average structure remains reasonable, the proteins move in a native‐like fashion as dictated by the original force field. To validate this approach, we have used the data from solid‐state NMR spectroscopy, which is the orthogonal experimental technique uniquely sensitive to protein local dynamics. The new method has been tested on the well‐established model protein, ubiquitin. The ensemble‐restrained MD simulations produced lower crystallographic R factors than conventional simulations; they also led to more accurate predictions for crystallographic temperature factors, solid‐state chemical shifts, and backbone order parameters. The predictions for ¹⁵N relaxation rates are at least as accurate as those obtained from conventional simulations. Taken together, these results suggest that the presented trajectories may be among the most realistic protein MD simulations ever reported. In this context, the ensemble restraints based on high‐resolution crystallographic data can be viewed as protein‐specific empirical corrections to the standard force fields.     

SnugDock: Paratope Structural Optimization during Antibody-Antigen Docking Compensates for Errors in Antibody Homology Models

High resolution structures of antibody-antigen complexes are useful for analyzing the binding interface and to make rational choices for antibody engineering. When a crystallographic structure of a complex is unavailable, the structure must be predicted using computational tools. In this work, we illustrate a novel approach, named SnugDock, to predict high-resolution antibody-antigen complex structures by simultaneously structurally optimizing the antibody-antigen rigid-body positions, the relative orientation of the antibody light and heavy chains, and the conformations of the six complementarity determining region loops. This approach is especially useful when the crystal structure of the antibody is not available, requiring allowances for inaccuracies in an antibody homology model which would otherwise frustrate rigid-backbone docking predictions. Local docking using SnugDock with the lowest-energy RosettaAntibody homology model produced more accurate predictions than standard rigid-body docking. SnugDock can be combined with ensemble docking to mimic conformer selection and induced fit resulting in increased sampling of diverse antibody conformations. The combined algorithm produced four medium (Critical Assessment of PRediction of Interactions-CAPRI rating) and seven acceptable lowest-interface-energy predictions in a test set of fifteen complexes. Structural analysis shows that diverse paratope conformations are sampled, but docked paratope backbones are not necessarily closer to the crystal structure conformations than the starting homology models. The accuracy of SnugDock predictions suggests a new genre of general docking algorithms with flexible binding interfaces targeted towards making homology models useful for further high-resolution predictions.     

Scalable Nanohelices for Predictive Studies and Enhanced 3D Visualization

Spring-like materials are ubiquitous in nature and of interest in nanotechnology for energy harvesting, hydrogen storage, and biological sensing applications.  For predictive simulations, it has become increasingly important to be able to model the structure of nanohelices accurately.  To study the effect of local structure on the properties of these complex geometries one must develop realistic models.  To date, software packages are rather limited in creating atomistic helical models.  This work focuses on producing atomistic models of silica glass (SiO₂) nanoribbons and nanosprings for molecular dynamics (MD) simulations. Using an MD model of “bulk” silica glass, two computational procedures to precisely create the shape of nanoribbons and nanosprings are presented.  The first method employs the AWK programming language and open-source software to effectively carve various shapes of silica nanoribbons from the initial bulk model, using desired dimensions and parametric equations to define a helix.  With this method, accurate atomistic silica nanoribbons can be generated for a range of pitch values and dimensions.  The second method involves a more robust code which allows flexibility in modeling nanohelical structures.  This approach utilizes a C++ code particularly written to implement pre-screening methods as well as the mathematical equations for a helix, resulting in greater precision and efficiency when creating nanospring models.  Using these codes, well-defined and scalable nanoribbons and nanosprings suited for atomistic simulations can be effectively created.  An added value in both open-source codes is that they can be adapted to reproduce different helical structures, independent of material.  In addition, a MATLAB graphical user interface (GUI) is used to enhance learning through visualization and interaction for a general user with the atomistic helical structures.  One application of these methods is the recent study of nanohelices via MD simulations for mechanical energy harvesting purposes.     

Mechanism of polyubiquitin chain recognition by the human ubiquitin conjugating enzyme Ube2g2

Ube2g2 is a human ubiquitin conjugating (E2) enzyme involved in the endoplasmic reticulum-associated degradation pathway, which is responsible for the identification and degradation of unfolded and misfolded proteins in the endoplasmic reticulum compartment. The Ube2g2-specific role is the assembly of Lys-48-linked polyubiquitin chains, which constitutes a signal for proteasomal degradation when attached to a substrate protein. NMR chemical shift perturbation and paramagnetic relaxation enhancement approaches were employed to characterize the binding interaction between Ube2g2 and ubiquitin, Lys-48-linked diubiquitin, and Lys-63-linked diubiquitin. Results demonstrate that ubiquitin binds to Ube2g2 with an affinity of 90 μM in two different orientations that are rotated by 180° in models generated by the RosettaDock modeling suite. The binding of Ube2g2 to Lys-48- and Lys-63-linked diubiquitin is primarily driven by interactions with individual ubiquitin subunits, with a clear preference for the subunit containing the free Lys-48 or Lys-63 side chain (i.e. the distal subunit). This preference is particularly striking in the case of Lys-48-linked diubiquitin, which exhibits an ～3-fold difference in affinities between the two ubiquitin subunits. This difference can be attributed to the partial steric occlusion of the subunit whose Lys-48 side chain is involved in the isopeptide linkage. As such, these results suggest that Lys-48-linked polyubiquitin chains may be designed to bind certain proteins like Ube2g2 such that the terminal ubiquitin subunit carrying the reactive Lys-48 side chain can be positioned properly for chain elongation regardless of chain length.