Determinants of functionality in the ubiquitin conjugating enzyme family

A transposon insertion in the promoter of OsUBC12 enhances cold tolerance during japonica rice germination

Low-temperature germination (LTG) is an important agronomic trait for rice (Oryza sativa). Japonica rice generally has greater capacity for germination at low temperatures than the indica subpopulation. However, the genetic basis and molecular mechanisms underlying this complex trait are poorly understood. Here, we report that OsUBC12, encoding an E2 ubiquitin-conjugating enzyme, increases low-temperature germinability in japonica, owing to a transposon insertion in its promoter enhancing its expression. Natural variation analysis reveals that transposon insertion in the OsUBC12 promoter mainly occurs in the japonica lineage. The variation detected in eight representative two-line male sterile lines suggests the existence of this allele introgression by indica-japonica hybridization breeding, and varieties carrying the japonica OsUBC12 locus (transposon insertion) have higher low-temperature germinability than varieties without the locus. Further molecular analysis shows that OsUBC12 negatively regulate ABA signaling. OsUBC12-regulated seed germination and ABA signaling mainly depend on a conserved active site required for ubiquitin-conjugating enzyme activity. Furthermore, OsUBC12 directly associates with rice SUCROSE NON-FERMENTING 1-RELATED PROTEIN KINASE 1.1 (OsSnRK1.1), promoting its degradation. OsSnRK1.1 inhibits LTG by enhancing ABA signaling and acts downstream of OsUBC12. These findings shed light on the underlying mechanisms of UBC12 regulating LTG and provide genetic reference points for improving LTG in indica rice.

Characterization of the chimeric protein cUBC1 engineered by substituting the linker of E2-25K into UBC1 enzyme of Saccharomyces cerevisiae

Ubiquitination is an important posttranslational modification of proteins in eukaryotic cells, wherein ubiquitin molecules are conjugated to target proteins. Ubiquitination is catalyzed by the cascade of ubiquitin activating enzyme (E1), ubiquitin conjugating enzyme (E2), and ubiquitin ligase (E3). The number of E2s encoded in eukaryotes partly explains their contribution to the inherent specificity of the ubiquitin system. The ubiquitin conjugating enzyme UBC1 of Saccharomyces cerevisiae participates the degradation of short-lived and abnormal proteins. UBC1 consists of two well-defined domains separated by a long flexible linker. E2-25K, the human homolog of UBC1 is crucial to neurons and its failure leads to neurodegenerative disorders. The linker of UBC1 is of 22 amino acids, while that of E2-25K has 6 amino acids. To understand the importance of the linker, the chimeric protein, cUBC1 was constructed by substituting the linker of E2-25K in UBC1. cUBC1 shows minor changes in its secondary structure. cUBC1 expression in ubc1 deletion mutants showed no effect over growth, thermotolerance and resistance to antibiotic stress. However, survival under heat stress was enhanced with cUBC1. Western blot analysis of the enzymatic activity showed cUBC1 performed equally well as UBC1. Hence, cUBC1 demonstrates that the shorter linker increased the stability of UBC1.

Functional conservation and divergence of the helix-turn-helix motif of E2 ubiquitin-conjugating enzymes

Polyubiquitination by E2 and E3 enzymes is crucial to cell cycle control, epigenetic regulation, and development. The hallmark of the E2 family is the ubiquitin (Ub)-conjugating (UBC) domain that forms a dynamic thioester conjugate with ubiquitin (E2~Ub). Numerous studies have focused on E2 surfaces, such as the N-terminal and crossover helices, that directly interact with an E3 or the conjugated ubiquitin to stabilize the active, "closed" state of the E2~Ub. However, it remains unclear how other E2 surfaces regulate ubiquitin transfer. Here, we demonstrate the helix-turn-helix (HTH) motif of the UBC tunes the intrinsic polyubiquitination activity through distinct functions in different E2s. Interestingly, the E2^HTH motif is repurposed in UBE2S and UBE2R2 to interact with the conjugated or acceptor ubiquitin, respectively, modulating ubiquitin transfer. Furthermore, we propose that Anaphase-Promoting Complex/Cyclosome binding to the UBE2S^HTH reduces the conformational space of the flexible E2~Ub, demonstrating an atypical E3-dependent activation mechanism. Altogether, we postulate the E2^HTH motif evolved to provide new functionalities that can be harnessed by E3s and permits additional regulation to facilitate specific E2-E3-mediated polyubiquitination.

UBE2W down-regulation promotes cell apoptosis and correlates with hypospermatogenesis

Ubiquitin conjugating enzyme (E2) is crucial for mediating N-terminal ubiquitination. Recent study reports that UBE2W is involved in male infertility. However, the correlation between UBE2W expression and hypospermatogenesis is unclear. The present study is to explore the biological role of UBE2W and its association with hypospermatogenesis. Results showed that the sexpression levels of UBE2W in mouse testes were gradually elevated from 2 to 10 weeks, while were significantly deceased in the testes with hypospermatogenesis. When UBE2W expression was successfully down-regulated in spermatogenic cells, the rate of apoptosis was significantly increased and the P53/Bcl-2/caspase 6/caspase 9 signal pathways were activated. Thus, these data indicate that UBE2W down-regulation promotes cell apoptosis and correlates with hypospermatogenesis, which may be helpful for the diagnosis of male infertility.

Evolution and expression analysis of the sorghum ubiquitin-conjugating enzyme family

Ubiquitin-conjugating enzymes (UBCs), which catalyse the transfer of ubiquitin to substrate or E3 ligases, are key enzymes in ubiquitination modifications of target proteins. Current knowledge regarding the sorghum (Sorghum bicolor (L.) Moench) ubiquitin-conjugating enzyme (SbUBC) family remains very limited. We identified 53 UBC-encoding genes in the sorghum genome and divided these into 18 groups according to their phylogenetic relationship with Arabidopsis thaliana (L.) Heynh., which was further supported by conserved motif and gene structure analyses. Different expression levels under a variety of abiotic stresses suggested that these might participate in distinct signalling pathways and that they underwent functional divergence during evolution. Furthermore, several SbUBC genes responded to single treatments, and individual SbUBC genes responded to multiple treatments, suggesting that sorghum UBCs may mediate crosstalk among different signalling pathways. Overall, the results provide valuable information for better understanding the classification and putative functions of sorghum UBC-encoding genes.

Functional reconstruction of a eukaryotic-like E1/E2/(RING) E3 ubiquitylation cascade from an uncultured archaeon

The covalent modification of protein substrates by ubiquitin regulates a diverse range of critical biological functions. Although it has been established that ubiquitin-like modifiers evolved from prokaryotic sulphur transfer proteins it is less clear how complex eukaryotic ubiquitylation system arose and diversified from these prokaryotic antecedents. The discovery of ubiquitin, E1-like, E2-like and small-RING finger (srfp) protein components in the Aigarchaeota and the Asgard archaea superphyla has provided a substantive step toward addressing this evolutionary question. Encoded in operons, these components are likely representative of the progenitor apparatus that founded the modern eukaryotic ubiquitin modification systems. Here we report that these proteins from the archaeon Candidatus ‘Caldiarchaeum subterraneum’ operate together as a bona fide ubiquitin modification system, mediating a sequential ubiquitylation cascade reminiscent of the eukaryotic process. Our observations support the hypothesis that complex eukaryotic ubiquitylation signalling pathways have developed from compact systems originally inherited from an archaeal ancestor.

In eukaryotic cells, the ubiquitylation system regulates several cellular processes central to protein homoeostasis. Here the authors demonstrate the existence of an eukaryotic-like ubiquitylation cascade requiring E1, E2 and E3-like enzymes in the archaeon C. subterraneum, shedding light on the evolution of the ubiquitin-proteasome system.

Characterization and expression analysis of genes encoding ubiquitin conjugating domain-containing enzymes in Carica papaya

Background

Ripening affects the quality and nutritional contents of fleshy fruits and is a crucial process of fruit development. Although several studies have suggested that ubiquitin-conjugating enzyme (E2s or UBC enzymes) are involved in the regulation of fruit ripening, little is known about the function of E2s in papaya (Carica papaya).

Methodology/Principal findings

In the present study, we searched the papaya genome and identified 34 putative UBC genes, which were clustered into 17 phylogenetic subgroups. We also analyzed the nucleotide sequences of the papaya UBC (CpUBC) genes and found that both exon-intron junctions and sequence motifs were highly conserved among the phylogenetic subgroups. Using real-time PCR analysis, we also found that all the CpUBC genes were expressed in roots, stems, leaves, male and female flowers, and mature fruit, although the expression of some of the genes was increased or decreased in one or several specific organs. We also found that the expression of 13 and two CpUBC genes were incresesd or decreased during one and two ripening stages, respectively. Expression analyses indicates possible E2s playing a more significant role in fruit ripening for further studies.

Conclusions

To the best of our knowledge, this is the first reported genome-wide analysis of the papaya UBC gene family, and the results will facilitate further investigation of the roles of UBC genes in fruit ripening and will aide in the functional validation of UBC genes in papaya.

The banana E2 gene family: Genomic identification, characterization, expression profiling analysis

The E2 is at the center of a cascade of Ub1 transfers, and it links activation of the Ub1 by E1 to its eventual E3-catalyzed attachment to substrate. Although the genome-wide analysis of this family has been performed in some species, little is known about analysis of E2 genes in banana. In this study, 74 E2 genes of banana were identified and phylogenetically clustered into thirteen subgroups. The predicted banana E2 genes were distributed across all 11 chromosomes at different densities. Additionally, the E2 domain, gene structure and motif compositions were analyzed. The expression of all of the banana E2 genes was analyzed in the root, stem, leaf, flower organs, five stages of fruit development and under abiotic stresses. All of the banana E2 genes, with the exception of few genes in each group, were expressed in at least one of the organs and fruit developments, which indicated that the E2 genes might involve in various aspects of the physiological and developmental processes of the banana. Quantitative RT-PCR (qRT-PCR) analysis identified that 45 E2s under drought and 33 E2s under salt were induced. To the best of our knowledge, this report describes the first genome-wide analysis of the banana E2 gene family, and the results should provide valuable information for understanding the classification, cloning and putative functions of this family.

Genome-Wide Identification, Phylogenetic and Expression Analyses of the Ubiquitin-Conjugating Enzyme Gene Family in Maize

Background

Ubiquitination is a post-translation modification where ubiquitin is attached to a substrate. Ubiquitin-conjugating enzymes (E2s) play a major role in the ubiquitin transfer pathway, as well as a variety of functions in plant biological processes. To date, no genome-wide characterization of this gene family has been conducted in maize (Zea mays).

Methodology/Principal Findings

In the present study, a total of 75 putative ZmUBC genes have been identified and located in the maize genome. Phylogenetic analysis revealed that ZmUBC proteins could be divided into 15 subfamilies, which include 13 ubiquitin-conjugating enzymes (ZmE2s) and two independent ubiquitin-conjugating enzyme variant (UEV) groups. The predicted ZmUBC genes were distributed across 10 chromosomes at different densities. In addition, analysis of exon-intron junctions and sequence motifs in each candidate gene has revealed high levels of conservation within and between phylogenetic groups. Tissue expression analysis indicated that most ZmUBC genes were expressed in at least one of the tissues, indicating that these are involved in various physiological and developmental processes in maize. Moreover, expression profile analyses of ZmUBC genes under different stress treatments (4°C, 20% PEG6000, and 200 mM NaCl) and various expression patterns indicated that these may play crucial roles in the response of plants to stress.

Conclusions

Genome-wide identification, chromosome organization, gene structure, evolutionary and expression analyses of ZmUBC genes have facilitated in the characterization of this gene family, as well as determined its potential involvement in growth, development, and stress responses. This study provides valuable information for better understanding the classification and putative functions of the UBC-encoding genes of maize.

Biochemical and structural characterization of a novel ubiquitin-conjugating enzyme E2 from Agrocybe aegeria reveals Ube2w family-specific properties

Ubiquitination is a post-translational modification that is involved in myriad cellar regulation and disease pathways. The ubiquitin-conjugating enzyme (E2) is an important player in the ubiquitin transfer pathway. Although many E2 structures are available, not all E2 families have known structures, and three-dimensional structures from fungal organisms other than yeast are lacking. We report here the crystal structure of UbcA1, which is a novel ubiquitin-conjugating enzyme identified from the edible and medicinal mushroom Agrocybe aegerita and displays potential antitumor properties. The protein belongs to the Ube2w family and shows similar biochemical characteristics to human Ube2w, including monomer-dimer equilibrium in solution, α-NH2 ubiquitin-transfer activity and a mechanism to recognize backbone atoms of intrinsically disordered N-termini in substrates. Its structure displays a unique C-terminal conformation with an orientation of helix α3 that is completely different from the reported E2 structures but similar to a recently reported NMR ensemble of Ube2w. A mutagenesis study on this novel enzyme revealed that an intact C-terminus is significant for protein dimerization and enzymatic activity. As the first crystallized full-length protein of this family, UbcA1 may supersede the truncated X-ray structure of Ube2w (PDB entry 2A7L) as the representative structure of the Ube2w family.

Selective multifaceted E3 ubiquitin ligases barricade extreme defense: Potential therapeutic targets for neurodegeneration and ageing

Efficient and regular performance of Ubiquitin Proteasome System and Autophagy continuously eliminate deleterious accumulation of nonnative protiens. In cellular quality control system, E3 ubiquitin ligases are significant employees for defense mechanism against abnormal toxic proteins. Few findings indicate that lack of functions of E3 ubiquitin ligases can be a causative factor of neurodevelopmental disorders, neurodegeneration, cancer and ageing. However, the detailed molecular pathomechanism implying E3 ubiquitin ligases in cellular functions in multifactorial disease conditions are not well understood. This article systematically represents the unique characteristics, molecular nature, and recent developments in the knowledge of neurobiological functions of few crucial E3 ubiquitin ligases. Here, we review recent literature on the roles of E6-AP, HRD1 and ITCH E3 ubiquitin ligases in the neuro-pathobiological mechanisms, with precise focus on the processes of neurodegeneration, and thereby propose new lines of potential targets for therapeutic interventions.

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 (GABAA 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.

Development of inhibitors in the ubiquitination cascade

The ubiquitin proteasome system (UPS) is essential in regulating myriad aspects of protein functions. It is therefore a fundamentally important regulatory mechanism that impacts most if not all aspects of cellular processes. Indeed, malfunction of UPS components is implicated in human diseases such as neurodegenerative and immunological disorders and many cancers. The success of proteasome inhibitors in cancer therapy suggests that modulating enzymes in the ubiquitination cascade would be clinically important for therapeutic benefits. In this review, we summarize advances in developing inhibitors of a variety of UPS components. In particular, we highlight recent work done on the protein engineering of ubiquitin as modulators of the UPS, a novel approach that may shed light on innovative drug discovery in the future.

Novel targets in drug design: enzymes in the protein ubiquitylation pathway

Protein ubiquitylation is a pathway by which many proteins are selectively degraded. Its role has been shown in processes such as cell division and differentiation, oncogenesis, apoptosis, DNA repair, membrane transport and the removal of abnormal proteins. The ubiquitylation pathway enzymes are an insufficiently researched area for drug development. A genetic method has been developed (supported by computational biology) to identify potentially useful small molecules that will have a positive impact on our battle against cancer and other diseases. In silico screening is used for initial selection of drug-like compounds. This method is based on docking three-dimensional chemical libraries onto the target enzyme's functional site for initial screens using a computational scheme, followed by genetic and in vivo methods for hit optimisation. Focus has been on using the ubiquitin conjugation pathway as target for therapeutic intervention against cancer and potent inhibitors of ubiquitylation subpathways have been obtained (including those that are vital for the survival of aggressive cancer cells/tumours). Leads from the development of in vitro inhibitors provided a direction for the development of in vivo inhibitors as investigational tools, and as promising therapeutic agents.

Tripartite motif ligases catalyze polyubiquitin chain formation through a cooperative allosteric mechanism

Ligation of polyubiquitin chains to proteins is a fundamental post-translational modification, often resulting in targeted degradation of conjugated proteins. Attachment of polyubiquitin chains requires the activities of an E1 activating enzyme, an E2 carrier protein, and an E3 ligase. The mechanism by which polyubiquitin chains are formed remains largely speculative, especially for RING-based ligases. The tripartite motif (TRIM) superfamily of ligases functions in many cellular processes including innate immunity, cellular localization, development and differentiation, signaling, and cancer progression. The present results show that TRIM ligases catalyze polyubiquitin chain formation in the absence of substrate, the rates of which can be used as a functional readout of enzyme function. Initial rate studies under biochemically defined conditions show that TRIM32 and TRIM25 are specific for the Ubc5 family of E2-conjugating proteins and, along with TRIM5α, exhibit cooperative kinetics with respect to Ubc5 concentration, with submicromolar [S]0.5 and Hill coefficients of 3-5, suggesting they possess multiple binding sites for their cognate E2-ubiquitin thioester. Mutation studies reveal a second, non-canonical binding site encompassing the C-terminal Ubc5α-helix. Polyubiquitin chain formation requires TRIM subunit oligomerization through the conserved coiled-coil domain, but can be partially replaced by fusing the catalytic domain to GST to promote dimerization. Other results suggest that TRIM32 assembles polyubiquitin chains as a Ubc5-linked thioester intermediate. These results represent the first detailed mechanistic study of TRIM ligase activity and provide a functional context for oligomerization observed in the superfamily.

Dynamic control of selectivity in the ubiquitination pathway revealed by an ASP to GLU substitution in an intra-molecular salt-bridge network

Ubiquitination relies on a subtle balance between selectivity and promiscuity achieved through specific interactions between ubiquitin-conjugating enzymes (E2s) and ubiquitin ligases (E3s). Here, we report how a single aspartic to glutamic acid substitution acts as a dynamic switch to tip the selectivity balance of human E2s for interaction toward E3 RING-finger domains. By combining molecular dynamic simulations, experimental yeast-two-hybrid screen of E2-E3 (RING) interactions and mutagenesis, we reveal how the dynamics of an internal salt-bridge network at the rim of the E2-E3 interaction surface controls the balance between an "open", binding competent, and a "closed", binding incompetent state. The molecular dynamic simulations shed light on the fine mechanism of this molecular switch and allowed us to identify its components, namely an aspartate/glutamate pair, a lysine acting as the central switch and a remote aspartate. Perturbations of single residues in this network, both inside and outside the interaction surface, are sufficient to switch the global E2 interaction selectivity as demonstrated experimentally. Taken together, our results indicate a new mechanism to control E2-E3 interaction selectivity at an atomic level, highlighting how minimal changes in amino acid side-chain affecting the dynamics of intramolecular salt-bridges can be crucial for protein-protein interactions. These findings indicate that the widely accepted sequence-structure-function paradigm should be extended to sequence-structure-dynamics-function relationship and open new possibilities for control and fine-tuning of protein interaction selectivity.

The interaction properties of the human Rab GTPase family--comparative analysis reveals determinants of molecular binding selectivity

Background

Rab GTPases constitute the largest subfamily of the Ras protein superfamily. Rab proteins regulate organelle biogenesis and transport, and display distinct binding preferences for effector and activator proteins, many of which have not been elucidated yet. The underlying molecular recognition motifs, binding partner preferences and selectivities are not well understood.

Methodology/Principal Findings

Comparative analysis of the amino acid sequences and the three-dimensional electrostatic and hydrophobic molecular interaction fields of 62 human Rab proteins revealed a wide range of binding properties with large differences between some Rab proteins. This analysis assists the functional annotation of Rab proteins 12, 14, 26, 37 and 41 and provided an explanation for the shared function of Rab3 and 27. Rab7a and 7b have very different electrostatic potentials, indicating that they may bind to different effector proteins and thus, exert different functions. The subfamily V Rab GTPases which are associated with endosome differ subtly in the interaction properties of their switch regions, and this may explain exchange factor specificity and exchange kinetics.

Conclusions/Significance

We have analysed conservation of sequence and of molecular interaction fields to cluster and annotate the human Rab proteins. The analysis of three dimensional molecular interaction fields provides detailed insight that is not available from a sequence-based approach alone. Based on our results, we predict novel functions for some Rab proteins and provide insights into their divergent functions and the determinants of their binding partner selectivity.

Human proteome-scale structural modeling of E2-E3 interactions exploiting interface motifs

Ubiquitination is crucial for many cellular processes such as protein degradation, DNA repair, transcription regulation, and cell signaling. Ubiquitin attachment takes place via a sequential enzymatic cascade involving ubiquitin activation (by E1 enzymes), ubiquitin conjugation (by E2 enzymes), and ubiquitin substrate tagging (by E3 enzymes). E3 ligases mediate ubiquitin transfer from E2s to substrates and as such confer substrate specificity. Although E3s can interact and function with numerous E2s, it is still unclear how they choose which E2 to use. Identifying all E2 partners of an E3 is essential for inferring the principles guiding E2 selection by an E3. Here we model the interactions of E3 and E2 proteins in a large, proteome-scale strategy based on interface structural motifs, which allows elucidation of (1) which E3s interact with which E2s in the human ubiquitination pathway and (2) how they interact with each other. Interface analysis of E2-E3 complexes reveals that loop L1 of E2s is critical for binding; the residue in the sixth position in loop L1 is widely utilized as an interface hot spot and appears indispensible for E2 interactions. Other loop L1 residues also confer specificity on the E2-E3 interactions: HECT E3s are in contact with the residue in the second position in loop L1 of E2s, but this is not the case for the RING finger type E3s. Our modeled E2-E3 complexes illuminate how slight sequence variations in E2 residues may contribute to specificity in E3 binding. These findings may be important for discovering drug candidates targeting E3s, which have been implicated in many diseases.

E1-E2 interactions in ubiquitin and Nedd8 ligation pathways

Initial rates of E1-catalyzed E2 transthiolation have been used as a reporter function to probe the mechanism of 125I-ubiquitin transfer between activation and ligation half-reactions of ubiquitin conjugation. A functional survey of 11 representative human E2 paralogs reveals similar Km for binding to human Uba1 ternary complex (Km(ave)=121±72 nm) and kcat for ubiquitin transfer (kcat(ave)=4.0±1.2 s(-1)), suggesting that they possess a conserved binding site and transition state geometry and that they compete for charging through differences in intracellular concentration. Sequence analysis and mutagenesis localize this binding motif to three basic residues within Helix 1 of the E2 core domain, confirmed by transthiolation kinetics. Partial conservation of the motif among E2 paralogs not recognized by Uba1 suggests that another factor(s) account for the absolute specificity of cognate E2 binding. Truncation of the Uba1 carboxyl-terminal β-grasp domain reduces cognate Ubc2b binding by 31-fold and kcat by 3.5×10(4)-fold, indicating contributions to E2 binding and transition state stabilization. Truncation of the paralogous domain from the Nedd8 activating enzyme has negligible effect on cognate Ubc12 transthiolation but abrogates E2 specificity toward non-cognate carrier proteins. Exchange of the β-grasp domains between ubiquitin and Nedd8 activating enzymes fails to reverse the effect of truncation. Thus, the conserved Helix 1 binding motif and the β-grasp domain direct general E2 binding, whereas the latter additionally serves as a specificity filter to exclude charging of non-cognate E2 paralogs in order to maintain the fidelity of downstream signaling.

Phosphorylation of threonine 190 is essential for nuclear localization and endocytosis of the FTS (Fused Toes Homolog) protein

Fused Toes Homolog (FTS) is a member of a group of proteins termed as E2 variants and this group of proteins lacks an active cysteine residue that is required for ubiquitin transfer. We have identified the expression of this protein in early neoplastic stages of cervical cancer and its translocation into nucleus from cytoplasm upon irradiation. Here we have reported that a threonine residue at position 190 is essential for its nucleocytoplasmic shuttling and function. Upon LMB treatment we found that FTS was located in the nucleus and it suggests that direct role of nuclear export signal (NES) is required for the binding to CRM1 and facilitates nuclear export. The threonine residue was phosphorylated and promoted the phosphorylation of EGFR, p38 and JNK facilitating vesicular trafficking of early to late endosomes. Mutational change of the threonine into alanine resulted in the cytoplasmic localization of FTS and failed to phosphorylate EGFR and its downstream effector proteins. In addition the mutation also reduced the number of early endosomes formed and also resulted in the clustering of late endosomes around the perinuclear region. These data suggest that threonine residue of FTS at position 190 is not only essential for its function but also for the formation, maturation and trafficking of early endosomes to late endosome/lysosome, as well as we speculate that FTS may function at a connection point in the vesicle tethering.

Atg3-mediated lipidation of Atg8 is involved in encystation of Acanthamoeba

Autophagy is a catabolic process involved in the degradation of a cell's own components for cell growth, development, homeostasis, and the recycling of cellular products. Autophagosome is an essential component in the protozoan parasite during differentiation and encystation. The present study identified and characterized autophagy-related protein (Atg) 3, a member of Atg8 conjugation system, in Acanthamoeba castellanii (AcAtg3). AcAtg3 encoding a 304 amino acid protein showed high similarity with the catalytic cysteine site of other E2 like enzymes of ubiquitin system. Predicted 3D structure of AcAtg3 revealed a hammer-like shape, which is the characteristic structure of E2-like enzymes. The expression level of AcAtg3 did not increase during encystation. However, the formation of mature cysts was significantly reduced in Atg3-siRNA transfected cells in which the production of Atg8-phosphatidylethanolamine conjugate was inhibited. Fluorescent microscopic analysis revealed that dispersed AcAtg3-EGFP fusion protein gathered around autophagosomal membranes during encystation. These results provide important information for understanding autophagic machinery through the lipidation reaction mediated by Atg3 in Acanthamoeba.

Evolutionary history of the poly(ADP-ribose) polymerase gene family in eukaryotes

BACKGROUND

The poly(ADP-ribose) polymerase (PARP) superfamily was originally identified as enzymes that catalyze the attachment of ADP-ribose subunits to target proteins using NAD+ as a substrate. The family is characterized by the catalytic site, termed the PARP signature. While these proteins can be found in a range of eukaryotes, they have been best studied in mammals. In these organisms, PARPs have key functions in DNA repair, genome integrity and epigenetic regulation. More recently it has been found that proteins within the PARP superfamily have altered catalytic sites, and have mono(ADP-ribose) transferase (mART) activity or are enzymatically inactive. These findings suggest that the PARP signature has a broader range of functions that initially predicted. In this study, we investigate the evolutionary history of PARP genes across the eukaryotes.

RESULTS

We identified in silico 236 PARP proteins from 77 species across five of the six eukaryotic supergroups. We performed extensive phylogenetic analyses of the identified PARPs. They are found in all eukaryotic supergroups for which sequence is available, but some individual lineages within supergroups have independently lost these genes. The PARP superfamily can be subdivided into six clades. Two of these clades were likely found in the last common eukaryotic ancestor. In addition, we have identified PARPs in organisms in which they have not previously been described.

CONCLUSIONS

Three main conclusions can be drawn from our study. First, the broad distribution and pattern of representation of PARP genes indicates that the ancestor of all extant eukaryotes encoded proteins of this type. Second, the ancestral PARP proteins had different functions and activities. One of these proteins was similar to human PARP1 and likely functioned in DNA damage response. The second of the ancestral PARPs had already evolved differences in its catalytic domain that suggest that these proteins may not have possessed poly(ADP-ribosyl)ation activity. Third, the diversity of the PARP superfamily is larger than previously documented, suggesting as more eukaryotic genomes become available, this gene family will grow in both number and type.

The family of ubiquitin-conjugating enzymes (E2s): deciding between life and death of proteins

The family of ubiquitin-conjugating (E2) enzymes is characterized by the presence of a highly conserved ubiquitin-conjugating (UBC) domain. These domains accommodate the ATP-activated ubiquitin (Ub) or ubiquitin-like (UBL) protein via a covalently linked thioester onto its active-site residue. E2 enzymes act via selective protein-protein interactions with the E1 and E3 enzymes and connect activation to covalent modification. By doing so, E2s differentiate effects on downstream substrates, either with a single Ub/UBL molecule or as a chain. While E3s are involved in substrate selection, E2s are the main determinants for selection of the lysine to construct ubiquitin chains, which thereby directly control the cellular fate of the substrate. In humans, 35 active E2 enzymes have been identified so far, while other eukaryotic genomes harbor 16 to 35 E2 family members. Some E2s possess N- and/or C-terminal extensions that mediate E2-specific processes. During the past two decades, strong support has led to the control of E2 enzymes in decisions concerning the life or death of a protein. Here, we summarize current knowledge and recent developments on E2 enzymes with respect to structural characteristics and functions. From this we propose a shell-like model to rationalize the selectivity of these key enzymes in directing Ub/UBL-conjugation pathways.-Van Wijk, S. J. L., Timmers, H. T. M. The family of ubiquitin-conjugating enzymes (E2s): deciding between life and death of proteins.

Cross-species analysis of the glycolytic pathway by comparison of molecular interaction fields

The electrostatic potential of an enzyme is a key determinant of its substrate interactions and catalytic turnover. Here we invoke comparative analysis of protein electrostatic potentials, along with sequence and structural analysis, to classify and characterize all the enzymes in an entire pathway across a set of different organisms. The electrostatic potentials of the enzymes from the glycolytic pathway of 11 eukaryotes were analyzed by qPIPSA (quantitative protein interaction property similarity analysis). The comparison allows the functional assignment of neuron-specific isoforms of triosephosphate isomerase from zebrafish, the identification of unusual protein surface interaction properties of the mosquito glucose-6-phosphate isomerase and the functional annotation of ATP-dependent phosphofructokinases and cofactor-dependent phosphoglycerate mutases from plants. We here show that plants possess two parallel pathways to convert glucose. One is similar to glycolysis in humans, the other is specialized to let plants adapt to their environmental conditions. We use differences in electrostatic potentials to estimate kinetic parameters for the triosephosphate isomerases from nine species for which published parameters are not available. Along the core glycolytic pathway, phosphoglycerate mutase displays the most conserved electrostatic potential. The largest cross-species variations are found for glucose-6-phosphate isomerase, enolase and fructose-1,6-bisphosphate aldolase. The extent of conservation of electrostatic potentials along the pathway is consistent with the absence of a single rate-limiting step in glycolysis.

Origin and evolution of ubiquitin-conjugating enzymes from Guillardia theta nucleomorph to hominoid

The origin of eukaryotic ubiquitin-conjugating enzymes (E2s) can be traced back to the Guillardia theta nucleomorph about 2500 million years ago (Mya). E2s are largely vertically inherited over eukaryotic evolution [Lespinet, O., Wolf, Y.I., Koonin, E.V., Aravind, L., 2002. The role of lineage-specific gene family expansion in the evolution of eukaryotes. Genome Res. 1048-1059], while mammal E2s experienced evolution of multigene families by gene duplications which have been accompanied by the increase in the species complexity. Because of alternatively splicing, primate-specific expansions of E2s happened once again at a transcriptional level. Both of them resulted in increasing genomic complexity and diversity of primate E2 proteomic function. The evolutionary processes of human E2 gene structure during expansions were accompanied by exon duplication and exonization of intronic sequences. Exonizations of Transposable Elements (TEs) in UBE2D3, UBE2L3 and UBE2V1 genes from primates indicate that exaptation of TEs also plays important roles in the structural innovation of primate-specific E2s and may create alternative splicing isoforms at a transcriptional level. Estimates for the ratio of dN/dS suggest that a strong purifying selection had acted upon protein-coding sequences of their orthologous UBE2D2, UBE2A, UBE2N, UBE2I and Rbx1 genes from animals, plants and fungi. The similar rates of synonymous substitutions are in accordance with the neutral mutation-random drift hypothesis of molecular evolution. Systematic detection of the origin and evolution of E2s, analyzing the evolution of E2 multigene families by gene duplications and the evolutionary processes of E2s during expansions, and testing its evolutionary force using E2s from distant phylogenetic lineages may advance our distinguishing of ancestral E2s from created E2s, and reveal previously unknown relationships between E2s and metazoan complexity. Analysis of these conserved proteins provides strong support for a close relationship between social amoeba and eukaryote, choanoflagellate and metazoan, and for the central roles of social amoeba and choanoflagellate in the origin and evolution of eukaryote and metazoan. Retracing the different stages of primate E2 exonization by monitoring genomic events over 63 Myr of primate evolution will advance our understanding of how TEs dynamically modified primate transcriptome and proteome in the past, and continue to do so.

What was the set of ubiquitin and ubiquitin-like conjugating enzymes in the eukaryote common ancestor?

Ubiquitin (Ub)-conjugating enzymes (E2) are key enzymes in ubiquitination or Ub-like modifications of proteins. We searched for all proteins belonging to the E2 enzyme super-family in seven species (Homo sapiens, Mus musculus, Drosophila melanogaster, Caenorhabditis elegans, Schizosaccharomyces pombe, Saccharomyces cerevisiae, and Arabidopsis thaliana) to identify families and to reconstruct each family’s phylogeny. Our phylogenetic analysis of 207 genes led us to define 17 E2 families, with 37 E2 genes, in the human genome. The subdivision of E2 into four classes did not correspond to the phylogenetic tree. The sequence signature HPN (histidine–proline–asparagine), followed by a tryptophan residue at 16 (up to 29) amino acids, was highly conserved. When present, the active cysteine was found 7 to 8 amino acids from the C-terminal end of HPN. The secondary structures were characterized by a canonical alpha/beta fold. Only family 10 deviated from the common organization because the proteins were devoid of enzymatic activity. Family 7 had an insertion between beta strands 1 and 2; families 3, 5 and 14 had an insertion between the active cysteine and the conserved tryptophan. The three-dimensional data of these proteins highlight a strong structural conservation of the core domain. Our analysis shows that the primitive eukaryote ancestor possessed a diversified set of E2 enzymes, thus emphasizing the importance of the Ub pathway. This comprehensive overview of E2 enzymes emphasizes the diversity and evolution of this superfamily and helps clarify the nomenclature and true orthologies. A better understanding of the functions of these enzymes is necessary to decipher several human diseases.

Rapid comparison of properties on protein surface

The mapping of physicochemical characteristics onto the surface of a protein provides crucial insights into its function and evolution. This information can be further used in the characterization and identification of similarities within protein surface regions. We propose a novel method which quantitatively compares global and local properties on the protein surface. We have tested the method on comparison of electrostatic potentials and hydrophobicity. The method is based on 3D Zernike descriptors, which provides a compact representation of a given property defined on a protein surface. Compactness and rotational invariance of this descriptor enable fast comparison suitable for database searches. The usefulness of this method is exemplified by studying several protein families including globins, thermophilic and mesophilic proteins, and active sites of TIM beta/alpha barrel proteins. In all the cases studied, the descriptor is able to cluster proteins into functionally relevant groups. The proposed approach can also be easily extended to other surface properties. This protein surface-based approach will add a new way of viewing and comparing proteins to conventional methods, which compare proteins in terms of their primary sequence or tertiary structure.

The ubiquitin system, disease, and drug discovery

The ubiquitin system of protein modification has emerged as a crucial mechanism involved in the regulation of a wide array of cellular processes. As our knowledge of the pathways in this system has grown, so have the ties between the protein ubiquitin and human disease. The power of the ubiquitin system for therapeutic benefit blossomed with the approval of the proteasome inhibitor Velcade in 2003 by the FDA. Current drug discovery activities in the ubiquitin system seek to (i) expand the development of new proteasome inhibitors with distinct mechanisms of action and improved bioavailability, and (ii) validate new targets. This review summarizes our current understanding of the role of the ubiquitin system in various human diseases ranging from cancer, viral infection and neurodegenerative disorders to muscle wasting, diabetes and inflammation. I provide an introduction to the ubiquitin system, highlight some emerging relationships between the ubiquitin system and disease, and discuss current and future efforts to harness aspects of this potentially powerful system for improving human health.

Republished from Current BioData's Targeted Proteins database (TPdb; ).

An FTS/Hook/p107(FHIP) complex interacts with and promotes endosomal clustering by the homotypic vacuolar protein sorting complex

Fused Toes (FTS) is a member of a small group of inactive variant E2 ubiquitin-conjugating enzyme domain-containing proteins of unknown function. Through proteomic analysis of FTS complexes purified from human embryonic kidney 293T cells, we identified a new multiprotein complex, the FHF complex, containing FTS, members of the microtubule-binding Hook family of coiled-coil proteins (Hook1, Hook2, and Hook3), and a previously uncharacterized 107-kDa protein, FTS and Hook Interacting Protein (FHIP). FTS associated with a conserved C-terminal motif in Hook proteins in the yeast two-hybrid system and in tissue culture cells, and Hook proteins were found to form homo- and heterodimers. The approximately 500-kDa FHF complex contained all three Hook proteins, and small interfering RNA depletion experiments suggest that Hook proteins can interact interchangeably within this complex. Hook proteins as well as FTS interact with members of both the class B and class C components of the homotypic vesicular protein sorting (HOPS) complex. Depletion of FTS by RNA interference affects both the trafficking of epidermal growth factor from early-to-late endosome/lysosomes and the efficiency by which overexpression of the HOPS component Vps18 promotes clustering of lysosomal-associated membrane protein 1-positive endosome/lysosomes. These data suggest that the FTS/Hook/FHIP complex functions to promote vesicle trafficking and/or fusion via the HOPS complex.

The basis for selective E1-E2 interactions in the ISG15 conjugation system

E1 and E2 enzymes coordinate the first steps in conjugation of ubiquitin (Ub) and ubiquitin-like proteins (Ubls). ISG15 is an interferon-alpha/beta-induced Ubl, and the E1 and E2 enzymes for ISG15 conjugation are Ube1L and UbcH8, respectively. UbcH7 is the most closely related E2 to UbcH8, yet it does not function in ISG15 conjugation in vivo, while both UbcH7 and UbcH8 have been reported to function in Ub conjugation. Kinetic analyses of wild-type and chimeric E2s were performed to determine the basis for preferential activation of UbcH8 by Ube1L and to determine whether UbcH8 is activated equally well by Ube1L and E1(Ub) (Ube1). K(m) determinations confirmed the strong preference of Ube1L for UbcH8 over UbcH7 (a 29-fold K(m) difference), similar to the preference of E1(Ub) for UbcH7 over UbcH8 (a 36-fold K(m) difference). Thioester assays of chimeric E2s identified two structural elements within residues 1-39 of UbcH8 that play a major role in defining Ube1L-UbcH8 specificity: the alpha1-helix and the beta1-beta2 region. The C-terminal ubiquitin fold domain (UFD) of Ube1L was required for transfer of ISG15 to UbcH8 and for binding of Ube1L to UbcH8. Replacement of the Ube1L UFD with that from E1(Ub) resulted in preferential transfer of ISG15 to UbcH7. Together, these results indicate that Ube1L discriminates between UbcH8 and closely related Ub E2s based on specific interactions between the Ube1L UFD and determinants within the N-terminal region of UbcH8.

Interactions between the quality control ubiquitin ligase CHIP and ubiquitin conjugating enzymes

Background

Ubiquitin (E3) ligases interact with specific ubiquitin conjugating (E2) enzymes to ubiquitinate particular substrate proteins. As the combination of E2 and E3 dictates the type and biological consequence of ubiquitination, it is important to understand the basis of specificity in E2:E3 interactions. The E3 ligase CHIP interacts with Hsp70 and Hsp90 and ubiquitinates client proteins that are chaperoned by these heat shock proteins. CHIP interacts with two types of E2 enzymes, UbcH5 and Ubc13-Uev1a. It is unclear, however, why CHIP binds these E2 enzymes rather than others, and whether CHIP interacts preferentially with UbcH5 or Ubc13-Uev1a, which form different types of polyubiquitin chains.

Results

The 2.9 Å crystal structure of the CHIP U-box domain complexed with UbcH5a shows that CHIP binds to UbcH5 and Ubc13 through similar specificity determinants, including a key S-P-A motif on the E2 enzymes. The determinants make different relative contributions to the overall interactions between CHIP and the two E2 enzymes. CHIP undergoes auto-ubiquitination by UbcH5 but not by Ubc13-Uev1a. Instead, CHIP drives the formation of unanchored polyubiquitin by Ubc13-Uev1a. CHIP also interacts productively with the class III E2 enzyme Ube2e2, in which the UbcH5- and Ubc13-binding specificity determinants are highly conserved.

Conclusion

The CHIP:UbcH5a structure emphasizes the importance of specificity determinants located on the long loops and central helix of the CHIP U-box, and on the N-terminal helix and loops L4 and L7 of its cognate E2 enzymes. The S-P-A motif and other specificity determinants define the set of cognate E2 enzymes for CHIP, which likely includes several Class III E2 enzymes. CHIP's interactions with UbcH5, Ube2e2 and Ubc13-Uev1a are consistent with the notion that Ubc13-Uev1a may work sequentially with other E2 enzymes to carry out K63-linked polyubiquitination of CHIP substrates.

webPIPSA: a web server for the comparison of protein interaction properties

Protein molecular interaction fields are key determinants of protein functionality. PIPSA (Protein Interaction Property Similarity Analysis) is a procedure to compare and analyze protein molecular interaction fields, such as the electrostatic potential. PIPSA may assist in protein functional assignment, classification of proteins, the comparison of binding properties and the estimation of enzyme kinetic parameters. webPIPSA is a web server that enables the use of PIPSA to compare and analyze protein electrostatic potentials. While PIPSA can be run with downloadable software (see http://projects.eml.org/mcm/software/pipsa), webPIPSA extends and simplifies a PIPSA run. This allows non-expert users to perform PIPSA for their protein datasets. With input protein coordinates, the superposition of protein structures, as well as the computation and analysis of electrostatic potentials, is automated. The results are provided as electrostatic similarity matrices from an all-pairwise comparison of the proteins which can be subjected to clustering and visualized as epograms (tree-like diagrams showing electrostatic potential differences) or heat maps. webPIPSA is freely available at: http://pipsa.eml.org.

qPIPSA: relating enzymatic kinetic parameters and interaction fields

Background

The simulation of metabolic networks in quantitative systems biology requires the assignment of enzymatic kinetic parameters. Experimentally determined values are often not available and therefore computational methods to estimate these parameters are needed. It is possible to use the three-dimensional structure of an enzyme to perform simulations of a reaction and derive kinetic parameters. However, this is computationally demanding and requires detailed knowledge of the enzyme mechanism. We have therefore sought to develop a general, simple and computationally efficient procedure to relate protein structural information to enzymatic kinetic parameters that allows consistency between the kinetic and structural information to be checked and estimation of kinetic constants for structurally and mechanistically similar enzymes.

Results

We describe qPIPSA: quantitative Protein Interaction Property Similarity Analysis. In this analysis, molecular interaction fields, for example, electrostatic potentials, are computed from the enzyme structures. Differences in molecular interaction fields between enzymes are then related to the ratios of their kinetic parameters. This procedure can be used to estimate unknown kinetic parameters when enzyme structural information is available and kinetic parameters have been measured for related enzymes or were obtained under different conditions. The detailed interaction of the enzyme with substrate or cofactors is not modeled and is assumed to be similar for all the proteins compared. The protein structure modeling protocol employed ensures that differences between models reflect genuine differences between the protein sequences, rather than random fluctuations in protein structure.

Conclusion

Provided that the experimental conditions and the protein structural models refer to the same protein state or conformation, correlations between interaction fields and kinetic parameters can be established for sets of related enzymes. Outliers may arise due to variation in the importance of different contributions to the kinetic parameters, such as protein stability and conformational changes. The qPIPSA approach can assist in the validation as well as estimation of kinetic parameters, and provide insights into enzyme mechanism.

Sequence determinants of E2-E6AP binding affinity and specificity

The conjugation of ubiquitin to substrates requires a series of enzymatic reactions consisting of an activating enzyme (E1), conjugating enzymes (E2) and ligases (E3). Tagging the appropriate substrate with ubiquitin is achieved by specific E2-E3 and E3-substrate interactions. E6AP, a member of the HECT family of E3s, has been previously shown to bind and function with the E2s UbcH7 and UbcH8. To decipher the sequence determinants of this specificity we have developed a quantitative E2-E3 binding assay based on fluorescence polarization and used this assay to measure the affinity of wild-type and mutant E2-E6AP interactions. Alanine scanning of the E6AP-UbcH7 binding interface identified four side-chains on UbcH7 and six side-chains on E6AP that contribute more than 1 kcal/mol to the binding free energy. Two of the hot spot residues from UbcH7 (K96 and K100) are conserved in UbcH8 but vary across other E2s. To determine if these are key specificity determining residues, we attempted to induce a tighter association between the E2 UbcH5b and E6AP by mutating the corresponding positions in UbcH5b to lysine residues. Surprisingly, the mutations had little effect, but rather a mutation at UbcH7 position 4, which is not at a hot spot on the UbcH7-E6AP interface, significantly strengthened UbcH5bs affinity for E6AP. This result indicates that E2-E3 binding specificities are a function of both favorable interactions that promote binding, and unfavorable interactions that prevent binding with unwanted partners.

Structure, interactions, and dynamics of the RING domain from human TRAF6

A key step in the signaling cascade responsible for activation of the transcription factor NF-kappaB involves Lys63-linked polyubiquitination of TRAF6. Covalent attachment of ubiquitin (Ub) to TRAF6, and subsequent poly(Ub) chain synthesis, is catalyzed by the hUev1a-hUbc13 heterodimer. hUbc13 is a catalytically competent E2 enzyme, and hUev1a is an E2-like protein that binds substrate Ub. The hUev1a-hUbc13 heterodimer is targeted to TRAF6 through interactions between hUbc13 and the N-terminal RING domain from TRAF6. Nuclear magnetic resonance (NMR) spectroscopy was used to determine the solution state structure of the RING domain from human TRAF6, and the interaction between hUbc13 and TRAF6 was characterized using NMR chemical shift mapping. The main-chain dynamics of the RING domain from TRAF6 were studied using (15)N NMR relaxation. Analysis of the main-chain dynamics data indicates that residues within the alpha-helix and beta-sheet of the RING domain are as rigid as regions of canonical secondary structure in larger proteins, consistent with the biological role of RING-domain E3 proteins, which requires that the E3 contain a recognition site for recruitment of E2 ubiquitin conjugation enzymes.

The crystal structure of Atg3, an autophagy-related ubiquitin carrier protein (E2) enzyme that mediates Atg8 lipidation

Atg3 is an E2-like enzyme that catalyzes the conjugation of Atg8 and phosphatidylethanolamine (PE). The Atg8-PE conjugate is essential for autophagy, which is the bulk degradation process of cytoplasmic components by the vacuolar/lysosomal system. We report here the crystal structure of Saccharomyces cerevisiae Atg3 at 2.5-A resolution. Atg3 has an alpha/beta-fold, and its core region is topologically similar to canonical E2 enzymes. Atg3 has two regions inserted in the core region, one of which consists of approximately 80 residues and has a random coil structure in solution and another with a long alpha-helical structure that protrudes from the core region as far as 30 A. In vivo and in vitro analyses suggested that the former region is responsible for binding Atg7, an E1-like enzyme, and that the latter is responsible for binding Atg8. A sulfate ion was bound near the catalytic cysteine of Atg3, suggesting a possible binding site for the phosphate moiety of PE. The structure of Atg3 provides a molecular basis for understanding the unique lipidation reaction that Atg3 carries out.

Issues in high-throughput comparative modelling: a case study using the ubiquitin E2 conjugating enzymes

Sequences of the ubiquitin-conjugating enzyme (UBC or E2) family were used as a test set to investigate issues associated with the high-throughput comparative modelling of protein structures. A semi-automatic method was initially developed with particular emphasis on producing models of a quality suitable for structural comparison. Structural and sequence features of the E2 family were used to improve the sequence alignment and the quality of the structural templates. Initially, failure to correct for subtle structural inconsistencies between templates lead to problems in the comparative analysis of the UBC electrostatic potentials. Modelling of known UBC structures using Modeller 4.0 showed that multiple templates produced, on average, no better models than the use of just one template, as judged by the root-mean-squared deviation between the comparative model and crystal structure backbones. Using four different quality-checking methods, for a given target sequence, it was not possible to distinguish the model most similar to the experimental structure. The UBC models were thus finally modelled using only the crystal structure template with the highest sequence identity to the target to be modelled, and producing only one model solution. Quality checking was used to reject models with obvious structural anomalies (e.g., bad side-chain packing). The resulting models have been used for a comparison of UBC structural features and of their electrostatic potentials. The work was extended through the development of a fully automated pipeline that identifies E2 sequences in the sequence databases, aligns and models them, and calculates the associated electrostatic potential.

Immunophysical properties and prediction of activities for vaccinia virus complement control protein and smallpox inhibitor of complement enzymes using molecular dynamics and electrostatics

We present immunophysical modeling for VCP, SPICE, and three mutants using MD simulations and Poisson-Boltzmann-type electrostatic calculations. VCP and SPICE are homologous viral proteins that control the complement system by imitating, structurally and functionally, natural regulators of complement activation. VCP and SPICE consist of four CCP modules connected with short flexible loops. MD simulations demonstrate that the rather complex modules of VCP/SPICE and their mutants exhibit a high degree of intermodular spatial mobility, which is affected by surface mutations. Electrostatic calculations using snapshots from the MD trajectories demonstrate variable spatial distribution of the electrostatic potentials, which suggests dynamic binding properties. We use covariance analysis to identify correlated modular oscillations. We also use electrostatic similarity indices to cluster proteins with common electrostatic properties. Our results are compared with experimental data to form correlations between the overall positive electrostatic potential of VCP/SPICE with binding and activity. We show how these correlations can be used to predict binding and activity properties. This work is expected to be useful for understanding the function of native CCP-containing regulators of complement activation and receptors and for the design of antiviral therapeutics and complement inhibitors.

Application of new multi-resolution methods for the comparison of biomolecular electrostatic properties in the absence of global structural similarity

In this paper we present a method for the multi-resolution comparison of biomolecular electrostatic potentials without the need for global structural alignment of the biomolecules. The underlying computational geometry algorithm uses multi-resolution attributed contour trees (MACTs) to compare the topological features of volumetric scalar fields. We apply the MACTs to compute electrostatic similarity metrics for a large set of protein chains with varying degrees of sequence, structure, and function similarity. For calibration, we also compute similarity metrics for these chains by a more traditional approach based upon 3D structural alignment and analysis of Carbo similarity indices. Moreover, because the MACT approach does not rely upon pairwise structural alignment, its accuracy and efficiency promises to perform well on future large-scale classification efforts across groups of structurally-diverse proteins. The MACT method discriminates between protein chains at a level comparable to the Carbo similarity index method; i.e., it is able to accurately cluster proteins into functionally-relevant groups which demonstrate strong dependence on ligand binding sites. The results of the analyses are available from the linked web databases http://ccvweb.cres.utexas.edu/MolSignature/ and http://agave.wustl.edu/similarity/. The MACT analysis tools are available as part of the public domain library of the Topological Analysis and Quantitative Tools (TAQT) from the Center of Computational Visualization, at the University of Texas at Austin (http://ccvweb.csres.utexas.edu/software). The Carbo software is available for download with the open-source APBS software package at http://apbs.sf.net/.

Structural basis for recruitment of Ubc12 by an E2 binding domain in NEDD8's E1

E2 conjugating enzymes play a central role in ubiquitin and ubiquitin-like protein (ublp) transfer cascades: the E2 accepts the ublp from the E1 enzyme and then the E2 often interacts with an E3 enzyme to promote ublp transfer to the target. We report here the crystal structure of a complex between the C-terminal domain from NEDD8's heterodimeric E1 (APPBP1-UBA3) and the catalytic core domain of NEDD8's E2 (Ubc12). The structure and associated mutational analyses reveal molecular details of Ubc12 recruitment by NEDD8's E1. Interestingly, the E1's Ubc12 binding domain resembles ubiquitin and recruits Ubc12 in a manner mimicking ubiquitin's interactions with ubiquitin binding domains. Structural comparison with E2-E3 complexes indicates that the E1 and E3 binding sites on Ubc12 may overlap and raises the possibility that crosstalk between E1 and E3 interacting with an E2 could influence the specificity and processivity of ublp transfer.