Catalytic proficiency of ubiquitin conjugation enzymes: balancing pK(a) suppression, entropy, and electrostatics

The Lysine Deprotonation Mechanism in a Ubiquitin Conjugating Enzyme

Ubiquitination is a biochemical reaction in which a small protein, ubiquitin (Ub), is covalently linked to a lysine on a target protein. This type of post-translational modification can signal for protein degradation, DNA repair, or inflammation response. Ubiquitination is catalyzed by three families of enzymes: ubiquitin activating enzymes (E1), ubiquitin conjugating enzymes (E2), and ubiquitin ligases (E3). In this study, we focus on the chemical mechanism used by the E2 enzyme, Ubc13, which forms polyubiquitin chains by linking a substrate Ub to Lys63 on a target ubiquitin (Ub*). Initially, Ubc13 is covalently linked to the substrate Ub. Next, Lys63 in the Ub* is deprotonated, becomes an active nucleophile, and attacks the thioester bond in the Ubc13∼Ub conjugate. The deprotonation mechanism is not well understood. There are two, conserved nearby residues that may act as conjugate bases (Asp119 on Ubc13 and Glu64 on Ub*.) It is also hypothesized that the active site environment suppresses the lysine's pKa, favoring deprotonated lysine. We test these hypotheses by simulating both WT and mutant Ubc13 with constant pH molecular dynamics (CpHMD), which allows titratable residues to change their protonation states. In our simulations, we have five titratable residues, including Lys63, and we use these simulations to monitor the protonation states and to generate titration curves of lysine 63. We found that the pKa of Lys63 is highly dependent on its distance from the active site. Also, mutating Asp119 or Glu64 to Ala has little effect on the lysine pKa, indicating that neither residue acts as a generalized base. Finally, we note that mutating a structural residue (Asn79 to Ala) increases the lysine pKa, suggesting that alterations to the active site hydrogen bonding network can affect nucleophile activation.

Structural inscrutabilities of Histone (H2BK123) monoubiquitination: A systematic review

Histone H2B monoubiquitination in budding yeast is a highly conserved post-translational modification. It is involved in normal functions of the cells like DNA Repair, RNA Pol II activation, trans-histone H3K and H79K methylation, meiosis, vesicle budding, etc. Deregulation of H2BK123ub can lead to the activation of proto-oncogenes and is also linked to neurodegenerative and heart diseases. Recent discoveries have enhanced the mechanistic underpinnings of H2BK123ub. For the first time, the Rad6's acidic tail has been implicated in histone recognition and interaction with Bre1's RBD domain. The non-canonical backside of Rad6 showed inhibition in polyubiquitination activity. Bre1 domains RBD and RING play a role in site-specific ubiquitination. The role of single Alaline residue in Rad6 activity. Understanding the mechanism of ubiquitination before moving to therapeutic applications is important. Current advancements in this field indicate the creation of novel therapeutic approaches and a foundation for further study.

A ubiquitin-specific, proximity-based labeling approach for the identification of ubiquitin ligase substrates

Over 600 E3 ligases in humans execute ubiquitination of specific target proteins in a spatiotemporal manner to elicit desired signaling effects. Here, we developed a ubiquitin-specific proximity-based labeling method to selectively biotinylate substrates of a given ubiquitin ligase. By fusing the biotin ligase BirA and an Avi-tag variant to the candidate E3 ligase and ubiquitin, respectively, we were able to specifically enrich bona fide substrates of a ligase using a one-step streptavidin pulldown under denaturing conditions. We applied our method, which we named Ub-POD, to the really interesting new gene (RING) E3 ligase RAD18 and identified proliferating cell nuclear antigen and several other critical players in the DNA damage repair pathway. Furthermore, we successfully applied Ub-POD to the RING ubiquitin ligase tumor necrosis factor receptor-associated factor 6 and a U-box-type E3 ubiquitin ligase carboxyl terminus of Hsc70-interacting protein. We anticipate that our method could be widely adapted to all classes of ubiquitin ligases to identify substrates.

Near atmospheric carbon dioxide activates plant ubiquitin cross-linking

Background Identifying CO₂-binding proteins is vital for our knowledge of CO₂-regulated molecular processes. The carbamate post-translational modification is a reversible CO₂-mediated adduct that can form on neutral N-terminal α-amino or lysine ε-amino groups. Methods We have developed triethyloxonium ion (TEO) as a chemical proteomics tool to trap the carbamate post-translational modification on protein covalently. We use ¹³C-NMR and TEO and identify ubiquitin as a plant CO₂-binding protein. Results We observe the carbamate post-translational modification on the Arabidopsis thaliana ubiquitin ε-amino groups of lysines 6, 33, and 48. We show that biologically relevant near atmospheric PCO₂ levels increase ubiquitin conjugation dependent on lysine 6. We further demonstrate that CO₂ increases the ubiquitin E2 ligase (AtUBC5) charging step via the transthioesterification reaction in which Ub is transferred from the E1 ligase active site to the E2 active site. Conclusions and general significance Therefore, plant ubiquitin is a CO₂-binding protein, and the carbamate post-translational modification represents a potential mechanism through which plant cells can respond to fluctuating PCO₂.

On the Possibility That Bond Strain Is the Mechanism of RING E3 Activation in the E2-Catalyzed Ubiquitination Reaction

Ubiquitination is a type of post-translational modification wherein the small protein ubiquitin (Ub) is covalently bound to a lysine on a target protein. Ubiquitination can signal for several regulatory pathways including protein degradation. Ubiquitination occurs by a series of reactions catalyzed by three types of enzymes: ubiquitin activating enzymes, E1; ubiquitin conjugating enzymes, E2; and ubiquitin ligases, E3. E2 enzymes directly catalyze the transfer of Ub to the target protein─the RING E3 improves the efficiency. Prior to its transfer, Ub is covalently linked to the E2 via a thioester bond and the Ub∼E2 conjugate forms a quaternary complex with the RING E3. It is hypothesized that the RING E3 improves the catalytic efficiency of ubiquitination by placing the E2∼Ub conjugate in a "closed" position, which tensions and weakens the thioester bond. We interrogate this hypothesis by analyzing the strain on the thioester during molecular dynamics simulations of both open and closed E2∼Ub/E3 complexes. Our data indicate that the thioester is strained when the E2∼Ub conjugate is in the closed position. We also show that the amount of strain is consistent with the experimental rate enhancement caused by the RING E3. Finally, our simulations show that the closed configuration increases the populations of key hydrogen bonds in the E2∼Ub active site. This is consistent with another hypothesis stating that the RING E3 enhances reaction rates by preorganizing the substrates.

Mutations of Rad6 E2 ubiquitin-conjugating enzymes at alanine-126 affect ubiquitination activity and decrease enzyme stability

Rad6, an E2 ubiquitin-conjugating enzyme conserved from yeast to humans, functions in transcription, genome maintenance, and proteostasis. The contributions of many conserved secondary structures of Rad6 and its human homologs UBE2A and UBE2B to their biological functions are not understood. A mutant RAD6 allele with a missense substitution at alanine-126 (A126) of helix-3 that causes defects in telomeric gene silencing, DNA repair, and protein degradation was reported over 2 decades ago. Here, using a combination of genetics, biochemical, biophysical, and computational approaches, we discovered that helix-3 A126 mutations compromise the ability of Rad6 to ubiquitinate target proteins without disrupting interactions with partner E3 ubiquitin-ligases that are required for their various biological functions in vivo. Explaining the defective in vitro or in vivo ubiquitination activities, molecular dynamics simulations and NMR showed that helix-3 A126 mutations cause local disorder of the catalytic pocket of Rad6 in addition to disorganizing the global structure of the protein to decrease its stability in vivo. We also show that helix-3 A126 mutations deform the structures of UBE2A and UBE2B, the human Rad6 homologs, and compromise the in vitro ubiquitination activity and folding of UBE2B. Providing insights into their ubiquitination defects, we determined helix-3 A126 mutations impair the initial ubiquitin charging and the final discharging steps during substrate ubiquitination by Rad6. In summary, our studies reveal that the conserved helix-3 is a crucial structural constituent that controls the organization of catalytic pockets, enzymatic activities, and biological functions of the Rad6-family E2 ubiquitin-conjugating enzymes.

Single-cell-resolved measurement of enzyme activity at the tissue level using drop-on-demand microkits

Drop-on-demand microkits with a diameter of ∼20 μm are used to measure the activity of acetylcholinesterase (AChE) in a brain slice with single-cell resolution. The relative standard deviation from 25 cellular regions reached 73.3% exhibiting the difference of enzyme activity in the brain slice. Therefore, this approach utilizing the well-established kits provides an alternative single-cell-resolved strategy for the elucidation of enzymatic heterogeneity at the tissue level.

The Role of Conformational Dynamics in the Recognition and Regulation of Ubiquitination

The ubiquitination pathway is central to many cell signaling and regulatory events. One of the intriguing aspects of the pathway is the combinatorial sophistication of substrate recognition and ubiquitin chain building determinations. The abundant structural and biological data portray several characteristic protein folds among E2 and E3 proteins, and the understanding of the combinatorial complexity that enables interaction with much of the human proteome is a major goal to developing targeted and selective manipulation of the pathway. With the commonality of some folds, there are likely other aspects that can provide differentiation and recognition. These aspects involve allosteric effects and conformational dynamics that can direct recognition and chain building processes. In this review, we will describe the current state of the knowledge for conformational dynamics across a wide timescale, address the limitations of present approaches, and illustrate the potential to make new advances in connecting dynamics with ubiquitination regulation.

RYBP Is a K63-Ubiquitin-Chain-Binding Protein that Inhibits Homologous Recombination Repair

Ring1-YY1-binding protein (RYBP) is a member of the non-canonical polycomb repressive complex 1 (PRC1), and like other PRC1 members, it is best described as a transcriptional regulator. However, several PRC1 members were recently shown to function in DNA repair. Here, we report that RYBP preferentially binds K63-ubiquitin chains via its Npl4 zinc finger (NZF) domain. Since K63-linked ubiquitin chains are assembled at DNA double-strand breaks (DSBs), we examined the contribution of RYBP to DSB repair. Surprisingly, we find that RYBP is K48 polyubiquitylated by RNF8 and rapidly removed from chromatin upon DNA damage by the VCP/p97 segregase. High expression of RYBP competitively inhibits recruitment of BRCA1 repair complex to DSBs, reducing DNA end resection and homologous recombination (HR) repair. Moreover, breast cancer cell lines expressing high endogenous RYBP levels show increased sensitivity to DNA-damaging agents and poly ADP-ribose polymerase (PARP) inhibition. These data suggest that RYBP negatively regulates HR repair by competing for K63-ubiquitin chain binding.

Role of Polyubiquitin Chain Flexibility in the Catalytic Mechanism of Cullin-RING Ubiquitin Ligases

Cullin-RING ubiquitin ligases are a diverse family of ubiquitin ligases that catalyze the synthesis of K48-linked polyubiquitin (polyUb) chains on a variety of substrates, ultimately leading to their degradation by the proteasome. The cullin-RING enzyme scaffold processively attaches a Ub molecule to the distal end of a growing chain up to lengths of eight Ub monomers. However, the molecular mechanism governing how chains of increasing size are built using a scaffold of largely fixed dimensions is not clear. We developed coarse-grained molecular dynamics simulations to describe the dependence of k_cat for cullin-RING ligases on the length and flexibility of the K48-linked polyUb chain attached to the substrate protein, key factors that determine the rate of subsequent Ub attachment to the chain, and therefore, the ensuing biological outcomes of ubiquitination. The results suggest that a number of regulatory mechanisms may lead to variations in the rate of chain elongation for different cullin-RING ligases. Specifically, modulation of the distance between the target lysine and the phosphodegron sequence of the substrate, the distance between the substrate lysine and the active site cysteine of the Ub conjugation enzyme (E2) bound to the cullin-RING scaffold, and flexibility of the bound E2 can lead to significant differences in the processing of K48-linked chains on substrates, potentially leading to differences in biological outcomes.

Active Site Gate Dynamics Modulate the Catalytic Activity of the Ubiquitination Enzyme E2-25K

The ubiquitin proteasome system (UPS) signals for degradation of proteins through attachment of K48-linked polyubiquitin chains, or alterations in protein-protein recognition through attachment of K63-linked chains. Target proteins are ubiquitinated in three sequential chemical steps by a three-component enzyme system. Ubiquitination, or E2 enzymes, catalyze the central step by facilitating reaction of a target protein lysine with the C-terminus of Ub that is attached to the active site cysteine of the E2 through a thioester bond. E2 reactivity is modulated by dynamics of an active site gate, whose central residue packs against the active site cysteine in a closed conformation. Interestingly, for the E2 Ubc13, which specifically catalyzes K63-linked ubiquitination, the central gate residue adopts an open conformation. We set out to determine if active site gate dynamics play a role in catalysis for E2-25K, which adopts the canonical, closed gate conformation, and which selectively synthesizes K48-linked ubiquitin chains. Gate dynamics were characterized using mutagenesis of key residues, combined with enzyme kinetics measurements, and main chain NMR relaxation. The experimental data were interpreted with all atom MD simulations. The data indicate that active site gate opening and closing rates for E2-25K are precisely balanced.

Ubc13: the Lys63 ubiquitin chain building machine

Ubc13 is an ubiquitin E2 conjugating enzyme that participates with many different E3 ligases to form lysine 63-linked (Lys63) ubiquitin chains that are critical to signaling in inflammatory and DNA damage response pathways. Recent studies have suggested Ubc13 as a potential therapeutic target for intervention in various human diseases including several different cancers, alleviation of anti-cancer drug resistance, chronic inflammation, and viral infections. Understanding a potential therapeutic target from different angles is important to assess its usefulness and potential pitfalls. Here we present a global review of Ubc13 from its structure, function, and cellular activities, to its natural and chemical inhibition. The aim of this article is to review the literature that directly implicates Ubc13 in a biological function, and to integrate structural and mechanistic insights into the larger role of this critical E2 enzyme. We discuss observations of multiple Ubc13 structures that suggest a novel mechanism for activation of Ubc13 that involves conformational change of the active site loop.

Ubiquitin phosphorylated at Ser57 hyper-activates parkin

Malfunction of the ubiquitin (Ub) E3 ligase, parkin, leads to defects in mitophagy and protein quality control linked to Parkinson's disease. Parkin activity is stimulated by phosphorylation of Ub at Ser65 (pUb^S65). Since the upstream kinase is only known for Ser65 (PINK1), the biochemical function of other phosphorylation sites on Ub remain largely unknown. We used fluorescently labelled and site-specifically phosphorylated Ub substrates to quantitatively relate the position and stoichiometry of Ub phosphorylation to parkin activation. Fluorescence measurements show that pUb^S65-stimulated parkin is 5-fold more active than auto-inhibited and un-stimulated parkin, which catalyzes a basal level of auto-ubiquitination. We consistently observed a low but detectable level of parkin activity with pUb^S12. Strikingly, pUb^S57 hyper-activates parkin, and our data demonstrate that parkin is able to selectively synthesize poly-pUb^S57 chains, even when 90% of the Ub in the reaction is un-phosphorylated. We further found that parkin ubiquitinates its physiological substrate Miro-1 with chains solely composed of pUb^S65 and more efficiently with pUb^S57 chains. Parkin hyper-activation by pUb^S57 demonstrates the first PINK1-independent route to active parkin, revealing the roles of multiple ubiquitin phosphorylation sites in governing parkin stimulation and catalytic activity. This article is part of a Special Issue entitled "Biochemistry of Synthetic Biology - Recent Developments" Guest Editor: Dr. Ilka Heinemann and Dr. Patrick O'Donoghue.

Molecular Basis for K63-Linked Ubiquitination Processes in Double-Strand DNA Break Repair: A Focus on Kinetics and Dynamics

Cells are exposed to thousands of DNA damage events on a daily basis. This damage must be repaired to preserve genetic information and prevent development of disease. The most deleterious damage is a double-strand break (DSB), which is detected and repaired by mechanisms known as non-homologous end-joining (NHEJ) and homologous recombination (HR), which are components of the DNA damage response system. NHEJ is an error-prone first line of defense, whereas HR invokes error-free repair and is the focus of this review. The functions of the protein components of HR-driven DNA repair are regulated by the coordinated action of post-translational modifications including lysine acetylation, phosphorylation, ubiquitination, and SUMOylation. The latter two mechanisms are fundamental for recognition of DSBs and reorganizing chromatin to facilitate repair. We focus on the structures and molecular mechanisms for the protein components underlying synthesis, recognition, and cleavage of K63-linked ubiquitin chains, which are abundant at damage sites and obligatory for DSB repair. The forward flux of the K63-linked ubiquitination cascade is driven by the combined activity of E1 enzyme, the heterodimeric E2 Mms2-Ubc13, and its cognate E3 ligases RNF8 and RNF168, which is balanced through the binding and cleavage of chains by the deubiquitinase BRCC36, and the proteasome, and through the binding of chains by recognition modules on repair proteins such as RAP80. We highlight a number of aspects regarding our current understanding for the role of kinetics and dynamics in determining the function of the enzymes and chain recognition modules that drive K63 ubiquitination.

Protein dynamics and function from solution state NMR spectroscopy

It is well-established that dynamics are central to protein function; their importance is implicitly acknowledged in the principles of the Monod, Wyman and Changeux model of binding cooperativity, which was originally proposed in 1965. Nowadays the concept of protein dynamics is formulated in terms of the energy landscape theory, which can be used to understand protein folding and conformational changes in proteins. Because protein dynamics are so important, a key to understanding protein function at the molecular level is to design experiments that allow their quantitative analysis. Nuclear magnetic resonance (NMR) spectroscopy is uniquely suited for this purpose because major advances in theory, hardware, and experimental methods have made it possible to characterize protein dynamics at an unprecedented level of detail. Unique features of NMR include the ability to quantify dynamics (i) under equilibrium conditions without external perturbations, (ii) using many probes simultaneously, and (iii) over large time intervals. Here we review NMR techniques for quantifying protein dynamics on fast (ps-ns), slow (μs-ms), and very slow (s-min) time scales. These techniques are discussed with reference to some major discoveries in protein science that have been made possible by NMR spectroscopy.

Quantitative analysis of protein-ligand interactions by NMR

Protein-ligand interactions have been commonly studied through static structures of the protein-ligand complex. Recently, however, there has been increasing interest in investigating the dynamics of protein-ligand interactions both for fundamental understanding of the underlying mechanisms and for drug development. NMR is a versatile and powerful tool, especially because it provides site-specific quantitative information. NMR has widely been used to determine the dissociation constant (KD), in particular, for relatively weak interactions. The simplest NMR method is a chemical-shift titration experiment, in which the chemical-shift changes of a protein in response to ligand titration are measured. There are other quantitative NMR methods, but they mostly apply only to interactions in the fast-exchange regime. These methods derive the dissociation constant from population-averaged NMR quantities of the free and bound states of a protein or ligand. In contrast, the recent advent of new relaxation-based experiments, including R2 relaxation dispersion and ZZ-exchange, has enabled us to obtain kinetic information on protein-ligand interactions in the intermediate- and slow-exchange regimes. Based on R2 dispersion or ZZ-exchange, methods that can determine the association rate, kon, dissociation rate, koff, and KD have been developed. In these approaches, R2 dispersion or ZZ-exchange curves are measured for multiple samples with different protein and/or ligand concentration ratios, and the relaxation data are fitted to theoretical kinetic models. It is critical to choose an appropriate kinetic model, such as the two- or three-state exchange model, to derive the correct kinetic information. The R2 dispersion and ZZ-exchange methods are suitable for the analysis of protein-ligand interactions with a micromolar or sub-micromolar dissociation constant but not for very weak interactions, which are typical in very fast exchange. This contrasts with the NMR methods that are used to analyze population-averaged NMR quantities. Essentially, to apply NMR successfully, both the type of experiment and equation to fit the data must be carefully and specifically chosen for the protein-ligand interaction under analysis. In this review, we first explain the exchange regimes and kinetic models of protein-ligand interactions, and then describe the NMR methods that quantitatively analyze these specific interactions.

E2 superfamily of ubiquitin-conjugating enzymes: constitutively active or activated through phosphorylation in the catalytic cleft

Protein phosphorylation is a modification that offers a dynamic and reversible mechanism to regulate the majority of cellular processes. Numerous diseases are associated with aberrant regulation of phosphorylation-induced switches. Phosphorylation is emerging as a mechanism to modulate ubiquitination by regulating key enzymes in this pathway. The molecular mechanisms underpinning how phosphorylation regulates ubiquitinating enzymes, however, are elusive. Here, we show the high conservation of a functional site in E2 ubiquitin-conjugating enzymes. In catalytically active E2s, this site contains aspartate or a phosphorylatable serine and we refer to it as the conserved E2 serine/aspartate (CES/D) site. Molecular simulations of substrate-bound and -unbound forms of wild type, mutant and phosphorylated E2s, provide atomistic insight into the role of the CES/D residue for optimal E2 activity. Both the size and charge of the side group at the site play a central role in aligning the substrate lysine toward E2 catalytic cysteine to control ubiquitination efficiency. The CES/D site contributes to the fingerprint of the E2 superfamily. We propose that E2 enzymes can be divided into constitutively active or regulated families. E2s characterized by an aspartate at the CES/D site signify constitutively active E2s, whereas those containing a serine can be regulated by phosphorylation.

Covalent Inhibition of Ubc13 Affects Ubiquitin Signaling and Reveals Active Site Elements Important for Targeting

Ubc13 is an E2 ubiquitin conjugating enzyme that functions in nuclear DNA damage signaling and cytoplasmic NF-κB signaling. Here, we present the structures of complexes of Ubc13 with two inhibitors, NSC697923 and BAY 11-7082, which inhibit DNA damage and NF-κB signaling in human cells. NSC697923 and BAY 11-7082 both inhibit Ubc13 by covalent adduct formation through a Michael addition at the Ubc13 active site cysteine. The resulting adducts of both compounds exploit a binding groove unique to Ubc13. We developed a Ubc13 mutant which resists NSC697923 inhibition and, using this mutant, we show that the inhibition of cellular DNA damage and NF-κB signaling by NSC697923 is largely due to specific Ubc13 inhibition. We propose that unique structural features near the Ubc13 active site could provide a basis for the rational development and design of specific Ubc13 inhibitors.

The HIP2~ubiquitin conjugate forms a non-compact monomeric thioester during di-ubiquitin synthesis

Polyubiquitination is a post-translational event used to control the degradation of damaged or unwanted proteins by modifying the target protein with a chain of ubiquitin molecules. One potential mechanism for the assembly of polyubiquitin chains involves the dimerization of an E2 conjugating enzyme allowing conjugated ubiquitin molecules to be put into close proximity to assist reactivity. HIP2 (UBE2K) and Ubc1 (yeast homolog of UBE2K) are unique E2 conjugating enzymes that each contain a C-terminal UBA domain attached to their catalytic domains, and they have basal E3-independent polyubiquitination activity. Although the isolated enzymes are monomeric, polyubiquitin formation activity assays show that both can act as ubiquitin donors or ubiquitin acceptors when in the activated thioester conjugate suggesting dimerization of the E2-ubiquitin conjugates. Stable disulfide complexes, analytical ultracentrifugation and small angle x-ray scattering were used to show that the HIP2-Ub and Ubc1-Ub thioester complexes remain predominantly monomeric in solution. Models of the HIP2-Ub complex derived from SAXS data show the complex is not compact but instead forms an open or backbent conformation similar to UbcH5b~Ub or Ubc13~Ub where the UBA domain and covalently attached ubiquitin reside on opposite ends of the catalytic domain. Activity assays showed that full length HIP2 exhibited a five-fold increase in the formation rate of di-ubiquitin compared to a HIP2 lacking the UBA domain. This difference was not observed for Ubc1 and may be attributed to the closer proximity of the UBA domain in HIP2 to the catalytic core than for Ubc1.

Lineage-specific interface proteins match up the cell cycle and differentiation in embryo stem cells

The shortage of molecular information on cell cycle changes along embryonic stem cell (ESC) differentiation prompts an in silico approach, which may provide a novel way to identify candidate genes or mechanisms acting in coordinating the two programs. We analyzed germ layer specific gene expression changes during the cell cycle and ESC differentiation by combining four human cell cycle transcriptome profiles with thirteen in vitro human ESC differentiation studies. To detect cross-talk mechanisms we then integrated the transcriptome data that displayed differential regulation with protein interaction data. A new class of non-transcriptionally regulated genes was identified, encoding proteins which interact systematically with proteins corresponding to genes regulated during the cell cycle or cell differentiation, and which therefore can be seen as interface proteins coordinating the two programs. Functional analysis gathered insights in fate-specific candidates of interface functionalities. The non-transcriptionally regulated interface proteins were found to be highly regulated by post-translational ubiquitylation modification, which may synchronize the transition between cell proliferation and differentiation in ESCs.

Exploring the RING-catalyzed ubiquitin transfer mechanism by MD and QM/MM calculations

Ubiquitylation is a universal mechanism for controlling cellular functions. A large family of ubiquitin E3 ligases (E3) mediates Ubiquitin (Ub) modification. To facilitate Ub transfer, RING E3 ligases bind both the substrate and ubiquitin E2 conjugating enzyme (E2) linked to Ub via a thioester bond to form a catalytic complex. The mechanism of Ub transfer catalyzed by RING E3 remains elusive. By employing a combined computational approach including molecular modeling, molecular dynamics (MD) simulations, and quantum mechanics/molecular mechanics (QM/MM) calculations, we characterized this catalytic mechanism in detail. The three-dimensional model of dimeric RING E3 ligase RNF4 RING, E2 ligase UbcH5A, Ub and the substrate SUMO2 shows close contact between the substrate and Ub transfer catalytic center. Deprotonation of the substrate lysine by D117 on UbcH5A occurs with almost no energy barrier as calculated by MD and QM/MM calculations. Then, the side chain of the activated lysine gets close to the thioester bond via a conformation change. The Ub transfer pathway begins with a nucleophilic addition that forms an oxyanion intermediate of a 4.23 kcal/mol energy barrier followed by nucleophilic elimination, resulting in a Ub modified substrate by a 5.65 kcal/mol energy barrier. These results provide insight into the mechanism of RING-catalyzed Ub transfer guiding the discovery of Ub system inhibitors.

The ubiquitin-conjugating enzyme, UbcM2, is restricted to monoubiquitylation by a two-fold mechanism that involves backside residues of E2 and Lys48 of ubiquitin

Proteins can be modified on lysines (K) with a single ubiquitin (Ub) or with polymers of Ub (polyUb). These different configurations and their respective topologies are primary factors for determining whether substrates are targeted to the proteasome for degradation or directed to nonproteolytic outcomes. We report here on the intrinsic ubiquitylation properties of UbcM2 (UBE2E3/UbcH9), a conserved Ub-conjugating enzyme linked to cell proliferation, development, and the cellular antioxidant defense system. Using a fully recombinant ubiquitylation assay, we show that UbcM2 is severely limited in its ability to synthesize polyUb chains with wild-type Ub. Restriction to monoubiquitylation is governed by multiple residues on the backside of the enzyme, far removed from its active site, and by lysine 48 of Ub. UbcM2 with mutated backside residues can synthesize K63-linked polyUb chains and to a lesser extent K6- and K48-linked chains. Additionally, we identified a single residue on the backside of the enzyme that promotes monoubiquitylation. Together, these findings reveal that a combination of noncatalytic residues within the Ubc catalytic core domain of UbcM2 as well as a lysine(s) within Ub can relegate a Ub-conjugating enzyme to monoubiquitylate its cognate targets despite having the latent capacity to construct polyUb chains. The two-fold mechanism for restricting activity to monoubiquitylation provides added insurance that UbcM2 will not build polyUb chains on its substrates, even under conditions of high local Ub concentrations.

Accuracy and precision of protein-ligand interaction kinetics determined from chemical shift titrations

NMR-monitored chemical shift titrations for the study of weak protein-ligand interactions represent a rich source of information regarding thermodynamic parameters such as dissociation constants (K ( D )) in the micro- to millimolar range, populations for the free and ligand-bound states, and the kinetics of interconversion between states, which are typically within the fast exchange regime on the NMR timescale. We recently developed two chemical shift titration methods wherein co-variation of the total protein and ligand concentrations gives increased precision for the K ( D ) value of a 1:1 protein-ligand interaction (Markin and Spyracopoulos in J Biomol NMR 53: 125-138, 2012). In this study, we demonstrate that classical line shape analysis applied to a single set of (1)H-(15)N 2D HSQC NMR spectra acquired using precise protein-ligand chemical shift titration methods we developed, produces accurate and precise kinetic parameters such as the off-rate (k ( off )). For experimentally determined kinetics in the fast exchange regime on the NMR timescale, k ( off ) ~ 3,000 s(-1) in this work, the accuracy of classical line shape analysis was determined to be better than 5 % by conducting quantum mechanical NMR simulations of the chemical shift titration methods with the magnetic resonance toolkit GAMMA. Using Monte Carlo simulations, the experimental precision for k ( off ) from line shape analysis of NMR spectra was determined to be 13 %, in agreement with the theoretical precision of 12 % from line shape analysis of the GAMMA simulations in the presence of noise and protein concentration errors. In addition, GAMMA simulations were employed to demonstrate that line shape analysis has the potential to provide reasonably accurate and precise k ( off ) values over a wide range, from 100 to 15,000 s(-1). The validity of line shape analysis for k ( off ) values approaching intermediate exchange (~100 s(-1)), may be facilitated by more accurate K ( D ) measurements from NMR-monitored chemical shift titrations, for which the dependence of K ( D ) on the chemical shift difference (Δω) between free and bound states is extrapolated to Δω = 0. The demonstrated accuracy and precision for k ( off ) will be valuable for the interpretation of biological kinetics in weakly interacting protein-protein networks, where a small change in the magnitude of the underlying kinetics of a given pathway may lead to large changes in the associated downstream signaling cascade.

Increased precision for analysis of protein-ligand dissociation constants determined from chemical shift titrations

NMR is ideally suited for the analysis of protein-protein and protein ligand interactions with dissociation constants ranging from ~2 μM to ~1 mM, and with kinetics in the fast exchange regime on the NMR timescale. For the determination of dissociation constants (K ( D )) of 1:1 protein-protein or protein-ligand interactions using NMR, the protein and ligand concentrations must necessarily be similar in magnitude to the K ( D ), and nonlinear least squares analysis of chemical shift changes as a function of ligand concentration is employed to determine estimates for the parameters K ( D ) and the maximum chemical shift change (Δδ(max)). During a typical NMR titration, the initial protein concentration, [P (0)], is held nearly constant. For this condition, to determine the most accurate parameters for K ( D ) and Δδ(max) from nonlinear least squares analyses requires initial protein concentrations that are ~0.5 × K ( D ), and a maximum concentration for the ligand, or titrant, of ~10 × [P (0)]. From a practical standpoint, these requirements are often difficult to achieve. Using Monte Carlo simulations, we demonstrate that co-variation of the ligand and protein concentrations during a titration leads to an increase in the precision of the fitted K ( D ) and Δδ(max) values when [P (0)] > K ( D ). Importantly, judicious choice of protein and ligand concentrations for a given NMR titration, combined with nonlinear least squares analyses using two independent variables (ligand and protein concentrations) and two parameters (K ( D ) and Δδ(max)) is a straightforward approach to increasing the accuracy of measured dissociation constants for 1:1 protein-ligand interactions.