Evolutionary constraint and disease associations of post-translational modification sites in human genomes

Post-translational modification sites are present in hydrophilic cavities of alpha-synuclein, tau, FUS, and TDP-43 fibrils: A molecular dynamics study

Hydration plays a crucial role in the refolding of intrinsically disordered proteins into amyloid fibrils; however, the specific interactions between water and protein that may contribute to this process are still unknown. In our previous studies of alpha-synuclein (aSyn), we have shown that waters confined in fibril cavities are stabilizing features of this pathological fold; and that amino acids that hydrogen bond with these confined waters modulate primary and seeded aggregation. Here, we extend our aSyn molecular dynamics (MD) simulations with three new polymorphs and correlate MD trajectory information with known post-translational modifications (PTMs) and experimental data. We show that cavity residues are more evolutionarily conserved than non-cavity residues and are enriched with PTM sites. As expected, the confinement within hydrophilic cavities results in more stably hydrated amino acids. Interestingly, cavity PTM sites display the longest protein-water hydrogen bond lifetimes, three-fold greater than non-PTM cavity sites. Utilizing the deep mutational screen dataset by Newberry et al. and the Thioflavin T aggregation review by Pancoe et al. parsed using a fibril cavity/non-cavity definition, we show that hydrophobic changes to amino acids in cavities have a larger effect on fitness and aggregation rate than residues outside cavities, supporting our hypothesis that these sites are involved in the inhibition of aSyn toxic fibrillization. Finally, we expand our study to include analysis of fibril structures of tau, FUS, TDP-43, prion, and hnRNPA1; all of which contained hydrated cavities, with tau, FUS, and TDP-43 recapitulating our PTM results in aSyn fibril cavities.

Protein diversification through post-translational modifications, alternative splicing, and gene duplication

Proteins provide the basis for cellular function. Having multiple versions of the same protein within a single organism provides a way of regulating its activity or developing novel functions. Post-translational modifications of proteins, by means of adding/removing chemical groups to amino acids, allow for a well-regulated and controlled way of generating functionally distinct protein species. Alternative splicing is another method with which organisms possibly generate new isoforms. Additionally, gene duplication events throughout evolution generate multiple paralogs of the same genes, resulting in multiple versions of the same protein within an organism. In this review, we discuss recent advancements in the study of these three methods of protein diversification and provide illustrative examples of how they affect protein structure and function.

Predicting Deleterious Non-Synonymous Single Nucleotide Polymorphisms (nsSNPs) of HRAS Gene and In Silico Evaluation of Their Structural and Functional Consequences towards Diagnosis and Prognosis of Cancer

The Harvey rat sarcoma (HRAS) proto-oncogene belongs to the RAS family and is one of the pathogenic genes that cause cancer. Deleterious nsSNPs might have adverse consequences at the protein level. This study aimed to investigate deleterious nsSNPs in the HRAS gene in predicting structural alterations associated with mutants that disrupt normal protein-protein interactions. Functional and structural analysis was employed in analyzing the HRAS nsSNPs. Putative post-translational modification sites and the changes in protein-protein interactions, which included a variety of signal cascades, were also investigated. Five different bioinformatics tools predicted 33 nsSNPs as "pathogenic" or "harmful". Stability analysis predicted rs1554885139, rs770492627, rs1589792804, rs730880460, rs104894227, rs104894227, and rs121917759 as unstable. Protein-protein interaction analysis revealed that HRAS has a hub connecting three clusters consisting of 11 proteins, and changes in HRAS might cause signal cascades to dissociate. Furthermore, Kaplan-Meier bioinformatics analyses indicated that the HRAS gene deregulation affected the overall survival rate of patients with breast cancer, leading to prognostic significance. Thus, based on these analyses, our study suggests that the reported nsSNPs of HRAS may serve as potential targets for different proteomic studies, diagnoses, and therapeutic interventions focusing on cancer.

Use of viral motif mimicry improves the proteome-wide discovery of human linear motifs

Linear motifs have an integral role in dynamic cell functions, including cell signaling. However, due to their small size, low complexity, and frequent mutations, identifying novel functional motifs poses a challenge. Viruses rely extensively on the molecular mimicry of cellular linear motifs. In this study, we apply systematic motif prediction combined with functional filters to identify human linear motifs convergently evolved also in viral proteins. We observe an increase in the sensitivity of motif prediction and improved enrichment in known instances. We identify >7,300 non-redundant motif instances at various confidence levels, 99 of which are supported by all functional and structural filters. Overall, we provide a pipeline to improve the identification of functional linear motifs from interactomics datasets and a comprehensive catalog of putative human motifs that can contribute to our understanding of the human domain-linear motif code and the associated mechanisms of viral interference.

Human phospho-signaling networks of SARS-CoV-2 infection are rewired by population genetic variants

SARS-CoV-2 infection hijacks signaling pathways and induces protein-protein interactions between human and viral proteins. Human genetic variation may impact SARS-CoV-2 infection and COVID-19 pathology; however, the genetic variation in these signaling networks remains uncharacterized. Here, we studied human missense single nucleotide variants (SNVs) altering phosphorylation sites modulated by SARS-CoV-2 infection, using machine learning to identify amino acid substitutions altering kinase-bound sequence motifs. We found 2,033 infrequent phosphorylation-associated SNVs (pSNVs) that are enriched in sequence motif alterations, potentially reflecting the evolution of signaling networks regulating host defenses. Proteins with pSNVs are involved in viral life cycle and host responses, including RNA splicing, interferon response (TRIM28), and glucose homeostasis (TBC1D4) with potential associations with COVID-19 comorbidities. pSNVs disrupt CDK and MAPK substrate motifs and replace these with motifs of Tank Binding Kinase 1 (TBK1) involved in innate immune responses, indicating consistent rewiring of signaling networks. Several pSNVs associate with severe COVID-19 and hospitalization (STARD13, ARFGEF2). Our analysis highlights potential genetic factors contributing to inter-individual variation of SARS-CoV-2 infection and COVID-19 and suggests leads for mechanistic and translational studies.

Structural analysis of TrkA mutations in patients with congenital insensitivity to pain reveals PLCγ as an analgesic drug target

Chronic pain is a major health issue, and the search for new analgesics has become increasingly important because of the addictive properties and unwanted side effects of opioids. To explore potentially new drug targets, we investigated mutations in the NTRK1 gene found in individuals with congenital insensitivity to pain with anhidrosis (CIPA). NTRK1 encodes tropomyosin receptor kinase A (TrkA), the receptor for nerve growth factor (NGF) and that contributes to nociception. Molecular modeling and biochemical analysis identified mutations that decreased the interaction between TrkA and one of its substrates and signaling effectors, phospholipase Cγ (PLCγ). We developed a cell-permeable phosphopeptide derived from TrkA (TAT-pQYP) that bound the Src homology domain 2 (SH2) of PLCγ. In HEK-293T cells, TAT-pQYP inhibited the binding of heterologously expressed TrkA to PLCγ and decreased NGF-induced, TrkA-mediated PLCγ activation and signaling. In mice, intraplantar administration of TAT-pQYP decreased mechanical sensitivity in an inflammatory pain model, suggesting that targeting this interaction may be analgesic. The findings demonstrate a strategy to identify new targets for pain relief by analyzing the signaling pathways that are perturbed in CIPA.

Computational and experimental methods for classifying variants of unknown clinical significance

The increase in sequencing capacity, reduction in costs, and national and international coordinated efforts have led to the widespread introduction of next-generation sequencing (NGS) technologies in patient care. More generally, human genetics and genomic medicine are gaining importance for more and more patients. Some communities are already discussing the prospect of sequencing each individual's genome at time of birth. Together with digital health records, this shall enable individualized treatments and preventive measures, so-called precision medicine. A central step in this process is the identification of disease causal mutations or variant combinations that make us more susceptible for diseases. Although various technological advances have improved the identification of genetic alterations, the interpretation and ranking of the identified variants remains a major challenge. Based on our knowledge of molecular processes or previously identified disease variants, we can identify potentially functional genetic variants and, using different lines of evidence, we are sometimes able to demonstrate their pathogenicity directly. However, the vast majority of variants are classified as variants of uncertain clinical significance (VUSs) with not enough experimental evidence to determine their pathogenicity. In these cases, computational methods may be used to improve the prioritization and an increasing toolbox of experimental methods is emerging that can be used to assay the molecular effects of VUSs. Here, we discuss how computational and experimental methods can be used to create catalogs of variant effects for a variety of molecular and cellular phenotypes. We discuss the prospects of integrating large-scale functional data with machine learning and clinical knowledge for the development of accurate pathogenicity predictions for clinical applications.

Glycomic and Glycoproteomic Techniques in Neurodegenerative Disorders and Neurotrauma: Towards Personalized Markers

The proteome represents all the proteins expressed by a genome, a cell, a tissue, or an organism at any given time under defined physiological or pathological circumstances. Proteomic analysis has provided unparalleled opportunities for the discovery of expression patterns of proteins in a biological system, yielding precise and inclusive data about the system. Advances in the proteomics field opened the door to wider knowledge of the mechanisms underlying various post-translational modifications (PTMs) of proteins, including glycosylation. As of yet, the role of most of these PTMs remains unidentified. In this state-of-the-art review, we present a synopsis of glycosylation processes and the pathophysiological conditions that might ensue secondary to glycosylation shortcomings. The dynamics of protein glycosylation, a crucial mechanism that allows gene and pathway regulation, is described. We also explain how-at a biomolecular level-mutations in glycosylation-related genes may lead to neuropsychiatric manifestations and neurodegenerative disorders. We then analyze the shortcomings of glycoproteomic studies, putting into perspective their downfalls and the different advanced enrichment techniques that emanated to overcome some of these challenges. Furthermore, we summarize studies tackling the association between glycosylation and neuropsychiatric disorders and explore glycoproteomic changes in neurodegenerative diseases, including Alzheimer's disease, Parkinson's disease, Huntington disease, multiple sclerosis, and amyotrophic lateral sclerosis. We finally conclude with the role of glycomics in the area of traumatic brain injury (TBI) and provide perspectives on the clinical application of glycoproteomics as potential diagnostic tools and their application in personalized medicine.

Enhancing the Discovery of Functional Post-Translational Modification Sites with Machine Learning Models - Development, Validation, and Interpretation

Protein posttranslational modifications (PTMs) are a rapidly expanding feature class of significant importance in cell biology. Due to a high burden of experimental proof, the number of functionals PTMs in the eukaryotic proteome is currently underestimated. Furthermore, not all PTMs are functionally equivalent. Computational approaches that can confidently recommend PTMs of probable function can improve the heuristics of PTM investigation and alleviate these problems. To address this need, we developed SAPH-ire: a multifeature heuristic neural network model that takes community wisdom into account by recommending experimental PTMs similar to those which have previously been established as having regulatory impact. Here, we describe the principle behind the SAPH-ire model, how it is developed, how we evaluate its performance, and important caveats to consider when building and interpreting such models. Finally, we discus current limitations of functional PTM prediction models and highlight potential mechanisms for their improvement.

Evolutionary and functional lessons from human-specific amino acid substitution matrices

Human genetic variation in coding regions is fundamental to the study of protein structure and function. Most methods for interpreting missense variants consider substitution measures derived from homologous proteins across different species. In this study, we introduce human-specific amino acid (AA) substitution matrices that are based on genetic variations in the modern human population. We analyzed the frequencies of >4.8M single nucleotide variants (SNVs) at codon and AA resolution and compiled human-centric substitution matrices that are fundamentally different from classic cross-species matrices (e.g. BLOSUM, PAM). Our matrices are asymmetric, with some AA replacements showing significant directional preference. Moreover, these AA matrices are only partly predicted by nucleotide substitution rates. We further test the utility of our matrices in exposing functional signals of experimentally-validated protein annotations. A significant reduction in AA transition frequencies was observed across nine post-translational modification (PTM) types and four ion-binding sites. Our results propose a purifying selection signal in the human proteome across a diverse set of functional protein annotations and provide an empirical baseline for interpreting human genetic variation in coding regions.

A global map of associations between types of protein posttranslational modifications and human genetic diseases

There are >200 types of protein posttranslational modifications (PTMs) described in eukaryotes, each with unique proteome coverage and functions. We hypothesized that some genetic diseases may be caused by the removal of a specific type of PTMs by genomic variants and the consequent deregulation of particular functions. We collected >320,000 human PTMs representing 59 types and crossed them with >4M nonsynonymous DNA variants annotated with predicted pathogenicity and disease associations. We report >1.74M PTM-variant co-occurrences that an enrichment analysis distributed into 215 pairwise associations between 18 PTM types and 148 genetic diseases. Of them, 42% were not previously described. Removal of lysine acetylation exerts the most pronounced effect, and less studied PTM types such as S-glutathionylation or S-nitrosylation show relevance. Using pathogenicity predictions, we identified PTM sites that may produce particular diseases if prevented. Our results provide evidence of a substantial impact of PTM-specific removal on the pathogenesis of genetic diseases and phenotypes.

•

There is an enrichment of disease-associated nsSNVs preventing certain types of PTMs

•

We report 215 pairwise associations between 18 PTM types and 148 genetic diseases

•

The removal of lysine acetylation exerts the most pronounced effect

•

We report a set of PTM sites that may produce particular diseases if prevented

ActiveDriverDB: Interpreting Genetic Variation in Human and Cancer Genomes Using Post-translational Modification Sites and Signaling Networks (2021 Update)

Deciphering the functional impact of genetic variation is required to understand phenotypic diversity and the molecular mechanisms of inherited disease and cancer. While millions of genetic variants are now mapped in genome sequencing projects, distinguishing functional variants remains a major challenge. Protein-coding variation can be interpreted using post-translational modification (PTM) sites that are core components of cellular signaling networks controlling molecular processes and pathways. ActiveDriverDB is an interactive proteo-genomics database that uses more than 260,000 experimentally detected PTM sites to predict the functional impact of genetic variation in disease, cancer and the human population. Using machine learning tools, we prioritize proteins and pathways with enriched PTM-specific amino acid substitutions that potentially rewire signaling networks via induced or disrupted short linear motifs of kinase binding. We then map these effects to site-specific protein interaction networks and drug targets. In the 2021 update, we increased the PTM datasets by nearly 50%, included glycosylation, sumoylation and succinylation as new types of PTMs, and updated the workflows to interpret inherited disease mutations. We added a recent phosphoproteomics dataset reflecting the cellular response to SARS-CoV-2 to predict the impact of human genetic variation on COVID-19 infection and disease course. Overall, we estimate that 16-21% of known amino acid substitutions affect PTM sites among pathogenic disease mutations, somatic mutations in cancer genomes and germline variants in the human population. These data underline the potential of interpreting genetic variation through the lens of PTMs and signaling networks. The open-source database is freely available at www.ActiveDriverDB.org.

The Burden of Post-Translational Modification (PTM)-Disrupting Mutations in the Tumor Matrisome

BACKGROUND

To evaluate the occurrence of mutations affecting post-translational modification (PTM) sites in matrisome genes across different tumor types, in light of their genomic and functional contexts and in comparison with the rest of the genome.

METHODS

This study spans 9075 tumor samples and 32 tumor types from The Cancer Genome Atlas (TCGA) Pan-Cancer cohort and identifies 151,088 non-silent mutations in the coding regions of the matrisome, of which 1811 affecting known sites of hydroxylation, phosphorylation, N- and O-glycosylation, acetylation, ubiquitylation, sumoylation and methylation PTM.

RESULTS

PTM-disruptive mutations (PTMmut) in the matrisome are less frequent than in the rest of the genome, seem independent of cell-of-origin patterns but show dependence on the nature of the matrisome protein affected and the background PTM types it generally harbors. Also, matrisome PTMmut are often found among structural and functional protein regions and in proteins involved in homo- and heterotypic interactions, suggesting potential disruption of matrisome functions.

CONCLUSIONS

Though quantitatively minoritarian in the spectrum of matrisome mutations, PTMmut show distinctive features and damaging potential which might concur to deregulated structural, functional, and signaling networks in the tumor microenvironment.

Non-Histone Protein Methylation: Biological Significance and Bioengineering Potential

Protein methylation is a key post-translational modification whose effects on gene expression have been intensively studied over the last two decades. Recently, renewed interest in non-histone protein methylation has gained momentum for its role in regulating important cellular processes and the activity of many proteins, including transcription factors, enzymes, and structural complexes. The extensive and dynamic role that protein methylation plays within the cell also highlights its potential for bioengineering applications. Indeed, while synthetic histone protein methylation has been extensively used to engineer gene expression, engineering of non-histone protein methylation has not been fully explored yet. Here, we report the latest findings, highlighting how non-histone protein methylation is fundamental for certain cellular functions and is implicated in disease, and review recent efforts in the engineering of protein methylation.

Comprehensive characterization of amino acid positions in protein structures reveals molecular effect of missense variants

Interpretation of the colossal number of genetic variants identified from sequencing applications is one of the major bottlenecks in clinical genetics, with the inference of the effect of amino acid-substituting missense variations on protein structure and function being especially challenging. Here we characterize the three-dimensional (3D) amino acid positions affected in pathogenic and population variants from 1,330 disease-associated genes using over 14,000 experimentally solved human protein structures. By measuring the statistical burden of variations (i.e., point mutations) from all genes on 40 3D protein features, accounting for the structural, chemical, and functional context of the variations' positions, we identify features that are generally associated with pathogenic and population missense variants. We then perform the same amino acid-level analysis individually for 24 protein functional classes, which reveals unique characteristics of the positions of the altered amino acids: We observe up to 46% divergence of the class-specific features from the general characteristics obtained by the analysis on all genes, which is consistent with the structural diversity of essential regions across different protein classes. We demonstrate that the function-specific 3D features of the variants match the readouts of mutagenesis experiments for BRCA1 and PTEN, and positively correlate with an independent set of clinically interpreted pathogenic and benign missense variants. Finally, we make our results available through a web server to foster accessibility and downstream research. Our findings represent a crucial step toward translational genetics, from highlighting the impact of mutations on protein structure to rationalizing the variants' pathogenicity in terms of the perturbed molecular mechanisms.

Post-translational Modification Crosstalk and Hotspots in Sirtuin Interactors Implicated in Cardiovascular Diseases

Sirtuins are protein deacetylases that play a protective role in cardiovascular diseases (CVDs), as well as many other diseases. Absence of sirtuins can lead to hyperacetylation of both nuclear and mitochondrial proteins leading to metabolic dysregulation. The protein post-translational modifications (PTMs) are known to crosstalk among each other to bring about complex phenotypic outcomes. Various PTM types such as acetylation, ubiquitination, and phosphorylation, and so on, drive transcriptional regulation and metabolism, but such crosstalks are poorly understood. We integrated protein–protein interactions (PPI) and PTMs from several databases to integrate information on 1,251 sirtuin-interacting proteins, of which 544 are associated with cardiac diseases. Based on the ∼100,000 PTM sites obtained for sirtuin interactors, we observed that the frequency of PTM sites (83 per protein), as well as PTM types (five per protein), is higher than the global average for human proteome. We found that ∼60–70% PTM sites fall into ordered regions. Approximately 83% of the sirtuin interactors contained at least one competitive crosstalk (in situ) site, with half of the sites occurring in CVD-associated proteins. A large proportion of identified crosstalk sites were observed for acetylation and ubiquitination competition. We identified 614 proteins containing PTM hotspots (≥5 PTM sites) and 133 proteins containing crosstalk hotspots (≥3 crosstalk sites). We observed that a large proportion of disease-associated sequence variants were found in PTM motifs of CVD proteins. We identified seven proteins (TP53, LMNA, MAPT, ATP2A2, NCL, APEX1, and HIST1H3A) containing disease-associated variants in PTM and crosstalk hotspots. This is the first comprehensive bioinformatics analysis on sirtuin interactors with respect to PTMs and their crosstalks. This study forms a platform for generating interesting hypotheses that can be tested for a deeper mechanistic understanding gained or derived from big-data analytics.

Parallel selection on ecologically relevant gene functions in the transcriptomes of highly diversifying salmonids

Background

Salmonid fishes are characterised by a very high level of variation in trophic, ecological, physiological, and life history adaptations. Some salmonid taxa show exceptional potential for fast, within-lake diversification into morphologically and ecologically distinct variants, often in parallel; these are the lake-resident charr and whitefish (several species in the genera Salvelinus and Coregonus). To identify selection on genes and gene categories associated with such predictable diversifications, we analysed 2702 orthogroups (4.82 Mbp total; average 4.77 genes/orthogroup; average 1783 bp/orthogroup). We did so in two charr and two whitefish species and compared to five other salmonid lineages, which do not evolve in such ecologically predictable ways, and one non-salmonid outgroup.

Results

All selection analyses are based on Coregonus and Salvelinus compared to non-diversifying taxa. We found more orthogroups were affected by relaxed selection than intensified selection. Of those, 122 were under significant relaxed selection, with trends of an overrepresentation of serine family amino acid metabolism and transcriptional regulation, and significant enrichment of behaviour-associated gene functions. Seventy-eight orthogroups were under significant intensified selection and were enriched for signalling process and transcriptional regulation gene ontology terms and actin filament and lipid metabolism gene sets. Ninety-two orthogroups were under diversifying/positive selection. These were enriched for signal transduction, transmembrane transport, and pyruvate metabolism gene ontology terms and often contained genes involved in transcriptional regulation and development. Several orthogroups showed signs of multiple types of selection. For example, orthogroups under relaxed and diversifying selection contained genes such as ap1m2, involved in immunity and development, and slc6a8, playing an important role in muscle and brain creatine uptake. Orthogroups under intensified and diversifying selection were also found, such as genes syn3, with a role in neural processes, and ctsk, involved in bone remodelling.

Conclusions

Our approach pinpointed relevant genomic targets by distinguishing among different kinds of selection. We found that relaxed, intensified, and diversifying selection affect orthogroups and gene functions of ecological relevance in salmonids. Because they were found consistently and robustly across charr and whitefish and not other salmonid lineages, we propose these genes have a potential role in the replicated ecological diversifications.

Precision Profiling of the Cardiovascular Post-Translationally Modified Proteome: Where There Is a Will, There Is a Way

There is an exponential increase in biological complexity as initial gene transcripts are spliced, translated into amino acid sequence, and post-translationally modified. Each protein can exist as multiple chemical or sequence-specific proteoforms, and each has the potential to be a critical mediator of a physiological or pathophysiological signaling cascade. Here, we provide an overview of how different proteoforms come about in biological systems and how they are most commonly measured using mass spectrometry-based proteomics and bioinformatics. Our goal is to present this information at a level accessible to every scientist interested in mass spectrometry and its application to proteome profiling. We will specifically discuss recent data linking various protein post-translational modifications to cardiovascular disease and conclude with a discussion for enablement and democratization of proteomics across the cardiovascular and scientific community. The aim is to inform and inspire the readership to explore a larger breadth of proteoform, particularity post-translational modifications, related to their particular areas of expertise in cardiovascular physiology.

ActiveDriverDB: human disease mutations and genome variation in post-translational modification sites of proteins

Interpretation of genetic variation is needed for deciphering genotype-phenotype associations, mechanisms of inherited disease, and cancer driver mutations. Millions of single nucleotide variants (SNVs) in human genomes are known and thousands are associated with disease. An estimated 21% of disease-associated amino acid substitutions corresponding to missense SNVs are located in protein sites of post-translational modifications (PTMs), chemical modifications of amino acids that extend protein function. ActiveDriverDB is a comprehensive human proteo-genomics database that annotates disease mutations and population variants through the lens of PTMs. We integrated >385,000 published PTM sites with ∼3.6 million substitutions from The Cancer Genome Atlas (TCGA), the ClinVar database of disease genes, and human genome sequencing projects. The database includes site-specific interaction networks of proteins, upstream enzymes such as kinases, and drugs targeting these enzymes. We also predicted network-rewiring impact of mutations by analyzing gains and losses of kinase-bound sequence motifs. ActiveDriverDB provides detailed visualization, filtering, browsing and searching options for studying PTM-associated mutations. Users can upload mutation datasets interactively and use our application programming interface in pipelines. Integrative analysis of mutations and PTMs may help decipher molecular mechanisms of phenotypes and disease, as exemplified by case studies of TP53, BRCA2 and VHL. The open-source database is available at https://www.ActiveDriverDB.org.

Minimotif Miner 4: a million peptide minimotifs and counting

Minimotif Miner (MnM) is a database and web system for analyzing short functional peptide motifs, termed minimotifs. We present an update to MnM growing the database from ∼300 000 to >1 000 000 minimotif consensus sequences and instances. This growth comes largely from updating data from existing databases and annotation of articles with high-throughput approaches analyzing different types of post-translational modifications. Another update is mapping human proteins and their minimotifs to know human variants from the dbSNP, build 150. Now MnM 4 can be used to generate mechanistic hypotheses about how human genetic variation affect minimotifs and outcomes. One example of the utility of the combined minimotif/SNP tool identifies a loss of function missense SNP in a ubiquitylation minimotif encoded in the excision repair cross-complementing 2 (ERCC2) nucleotide excision repair gene. This SNP reaches genome wide significance for many types of cancer and the variant identified with MnM 4 reveals a more detailed mechanistic hypothesis concerning the role of ERCC2 in cancer. Other updates to the web system include a new architecture with migration of the web system and database to Docker containers for better performance and management. Weblinks:minimotifminer.org and mnm.engr.uconn.edu

Homogeneous peptide-break assay for luminescent detection of enzymatic protein post-translational modification activity utilizing charged peptides

We have developed a rapid and sensitive universal peptide-based time-resolved luminescence assay for detection of enzymatic post-translational modifications (PTMs). PTMs play essential roles in intracellular signaling and cell regulation, thus providing functional protein diversity in cell. Due this, impaired PTM patterns have been linked to multiple disease states. Clear link between PTMs and pathological conditions have also driven assay development further, but still today most of the methodologies are based on single-specificity or group-specific PTM-recognition. We have previously introduced leuzine-zipper based peptide-break technology as a viable option for universal PTM detection. Here, we introduce peptide-break technology utilizing single-label homogeneous quenching resonance energy transfer (QRET) and charge-based peptide-peptide interaction. We demonstrate the functionality of the new assay concept in phosphorylation, deacetylation, and citrullination. In a comparable study between previously introduced leucine-zipper and the novel charge-based approach, we found equal PTM detection performance and sensitivity, but the peptide design for new targets is simplified with the charged peptides. The new concept allows the use of short <20 amino acid peptides without limitations rising from the leucine-zipper coiled-coil structure. Introduced methodology enables wash-free PTM detection in a 384-well plate format, using low nanomolar enzyme concentrations. Potentially, the peptide-break technique using charged peptides may be applicable for natural peptide sequences directly obtained from the target protein.

Improved measures for evolutionary conservation that exploit taxonomy distances

Selective pressures on protein-coding regions that provide fitness advantages can lead to the regions' fixation and conservation in genome duplications and speciation events. Consequently, conservation analyses relying on sequence similarities are exploited by a myriad of applications across all biosciences to identify functionally important protein regions. While very potent, existing conservation measures based on multiple sequence alignments are so pervasive that improvements to solutions of many problems have become incremental. We introduce a new framework for evolutionary conservation with measures that exploit taxonomy distances across species. Results show that our taxonomy-based framework comfortably outperforms existing conservation measures in identifying deleterious variants observed in the human population, including variants located in non-abundant sequence domains such as intrinsically disordered regions. The predictive power of our approach emphasizes that the phenotypic effects of sequence variants can be taxonomy-level specific and thus, conservation needs to be interpreted accordingly.

Information on protein sequence variability and conservation can be leveraged to identify functionally important regions. Here, the authors develop new conservation measures that exploit taxonomy distances and LIST, a tool for predicting deleteriousness of human variants.

Analyzing Biomolecular Ensembles

Several techniques are available to generate conformational ensembles of proteins and other biomolecules either experimentally or computationally. These methods produce a large amount of data that need to be analyzed to identify structure-dynamics-function relationship. In this chapter, we will cover different tools to unveil the information hidden in conformational ensemble data and to guide toward the rationalization of the data. We included routinely used approaches such as dimensionality reduction, as well as new methods inspired by high-order statistics and graph theory.

A resource of variant effect predictions of single nucleotide variants in model organisms

The effect of single nucleotide variants (SNVs) in coding and noncoding regions is of great interest in genetics. Although many computational methods aim to elucidate the effects of SNVs on cellular mechanisms, it is not straightforward to comprehensively cover different molecular effects. To address this, we compiled and benchmarked sequence and structure-based variant effect predictors and we computed the impact of nearly all possible amino acid and nucleotide variants in the reference genomes of Homo sapiens, Saccharomyces cerevisiae and Escherichia coli Studied mechanisms include protein stability, interaction interfaces, post-translational modifications and transcription factor binding sites. We apply this resource to the study of natural and disease coding variants. We also show how variant effects can be aggregated to generate protein complex burden scores that uncover protein complex to phenotype associations based on a set of newly generated growth profiles of 93 sequenced S. cerevisiae strains in 43 conditions. This resource is available through mutfunc (www.mutfunc.com), a tool by which users can query precomputed predictions by providing amino acid or nucleotide-level variants.

Maximizing the Utility of Cancer Transcriptomic Data

Transcriptomic profiling has been applied to large numbers of cancer samples, by large-scale consortia, including The Cancer Genome Atlas, International Cancer Genome Consortium, and Cancer Cell Line Encyclopedia. Advances in mining cancer transcriptomic data enable us to understand the endless complexity of the cancer transcriptome and thereby to discover new biomarkers and therapeutic targets. In this paper, we review computational resources for deep mining of transcriptomic data to identify, quantify, and determine the functional effects and clinical utility of transcriptomic events, including noncoding RNAs, post-transcriptional regulation, exogenous RNAs, and transcribed genetic variants. These approaches can be applied to other complex diseases, thereby greatly leveraging the impact of this work.

A Methodological Assessment and Characterization of Genetically-Driven Variation in Three Human Phosphoproteomes

Phosphorylation of proteins on serine, threonine, and tyrosine residues is a ubiquitous post-translational modification that plays a key part of essentially every cell signaling process. It is reasonable to assume that inter-individual variation in protein phosphorylation may underlie phenotypic differences, as has been observed for practically any other molecular regulatory phenotype. However, we do not know much about the extent of inter-individual variation in phosphorylation because it is quite challenging to perform a quantitative high throughput study to assess inter-individual variation in any post-translational modification. To test our ability to address this challenge with SILAC-based mass spectrometry, we quantified phosphorylation levels for three genotyped human cell lines within a nested experimental framework, and found that genetic background is the primary determinant of phosphoproteome variation. We uncovered multiple functional, biophysical, and genetic associations with germline driven phosphopeptide variation. Variants affecting protein levels or structure were among these associations, with the latter presenting, on average, a stronger effect. Interestingly, we found evidence that is consistent with a phosphopeptide variability buffering effect endowed from properties enriched within longer proteins. Because the small sample size in this 'pilot' study may limit the applicability of our genetic observations, we also undertook a thorough technical assessment of our experimental workflow to aid further efforts. Taken together, these results provide the foundation for future work to characterize inter-individual variation in post-translational modification levels and reveal novel insights into the nature of inter-individual variation in phosphorylation.

Minimotifs dysfunction is pervasive in neurodegenerative disorders

Minimotifs are modular contiguous peptide sequences in proteins that are important for posttranslational modifications, binding to other molecules, and trafficking to specific subcellular compartments. Some molecular functions of proteins in cellular pathways can be predicted from minimotif consensus sequences identified through experimentation. While a role for minimotifs in regulating signal transduction and gene regulation during disease pathogenesis (such as infectious diseases and cancer) is established, the therapeutic use of minimotif mimetic drugs is limited. In this review, we discuss a general theme identifying a pervasive role of minimotifs in the pathomechanism of neurodegenerative diseases. Beyond their longstanding history in the genetics of familial neurodegeneration, minimotifs are also major players in neurotoxic protein aggregation, aberrant protein trafficking, and epigenetic regulation. Generalizing the importance of minimotifs in neurodegenerative diseases offers a new perspective for the future study of neurodegenerative mechanisms and the investigation of new therapeutics.

Proteomic analysis and prediction of amino acid variations that influence protein posttranslational modifications

Accumulative studies have indicated that amino acid variations through changing the type of residues of the target sites or key flanking residues could directly or indirectly influence protein posttranslational modifications (PTMs) and bring about a detrimental effect on protein function. Computational mutation analysis can greatly narrow down the efforts on experimental work. To increase the utilization of current computational resources, we first provide an overview of computational prediction of amino acid variations that influence protein PTMs and their functional analysis. We also discuss the challenges that are faced while developing novel in silico approaches in the future. The development of better methods for mutation analysis-related protein PTMs will help to facilitate the development of personalized precision medicine.

Disruption of protein function by pathogenic mutations: common and uncommon mechanisms 1

Mutations in protein-coding regions underlie almost all Mendelian disorders, drive tumorigenesis, and contribute to susceptibility to common diseases. Despite the great diversity of diseases that are caused by coding mutations, the cellular processes that affect, and are affected by, pathogenic variants at the molecular level are fundamentally conserved. Experimental and computational approaches have revealed that a substantial fraction of disease mutations are not simple loss-of-function alleles. Rather, these pathogenic variants disrupt protein function in more subtle ways by tuning protein folding pathways, altering subcellular trafficking, interrupting signaling cascades, and rewiring highly connected interaction networks. Focusing mainly on Mendelian disorders, this review discusses the common mechanisms by which deleterious mutations disrupt protein function and how these disruptions can be exploited in the development of novel therapies.

Pathogenic Germline Variants in 10,389 Adult Cancers

We conducted the largest investigation of predisposition variants in cancer to date, discovering 853 pathogenic or likely pathogenic variants in 8% of 10,389 cases from 33 cancer types. Twenty-one genes showed single or cross-cancer associations, including novel associations of SDHA in melanoma and PALB2 in stomach adenocarcinoma. The 659 predisposition variants and 18 additional large deletions in tumor suppressors, including ATM, BRCA1, and NF1, showed low gene expression and frequent (43%) loss of heterozygosity or biallelic two-hit events. We also discovered 33 such variants in oncogenes, including missenses in MET, RET, and PTPN11 associated with high gene expression. We nominated 47 additional predisposition variants from prioritized VUSs supported by multiple evidences involving case-control frequency, loss of heterozygosity, expression effect, and co-localization with mutations and modified residues. Our integrative approach links rare predisposition variants to functional consequences, informing future guidelines of variant classification and germline genetic testing in cancer.

Prediction and interpretation of deleterious coding variants in terms of protein structural stability

The classification of human genetic variants into deleterious and neutral is a challenging issue, whose complexity is rooted in the large variety of biophysical mechanisms that can be responsible for disease conditions. For non-synonymous mutations in structured proteins, one of these is the protein stability change, which can lead to loss of protein structure or function. We developed a stability-driven knowledge-based classifier that uses protein structure, artificial neural networks and solvent accessibility-dependent combinations of statistical potentials to predict whether destabilizing or stabilizing mutations are disease-causing. Our predictor yields a balanced accuracy of 71% in cross validation. As expected, it has a very high positive predictive value of 89%: it predicts with high accuracy the subset of mutations that are deleterious because of stability issues, but is by construction unable of classifying variants that are deleterious for other reasons. Its combination with an evolutionary-based predictor increases the balanced accuracy up to 75%, and allowed predicting more than 1/4 of the variants with 95% positive predictive value. Our method, called SNPMuSiC, can be used with both experimental and modeled structures and compares favorably with other prediction tools on several independent test sets. It constitutes a step towards interpreting variant effects at the molecular scale. SNPMuSiC is freely available at https://soft.dezyme.com/.

Comprehensive Analysis of Constraint on the Spatial Distribution of Missense Variants in Human Protein Structures

The spatial distribution of genetic variation within proteins is shaped by evolutionary constraint and provides insight into the functional importance of protein regions and the potential pathogenicity of protein alterations. Here, we comprehensively evaluate the 3D spatial patterns of human germline and somatic variation in 6,604 experimentally derived protein structures and 33,144 computationally derived homology models covering 77% of all human proteins. Using a systematic approach, we quantify differences in the spatial distributions of neutral germline variants, disease-causing germline variants, and recurrent somatic variants. Neutral missense variants exhibit a general trend toward spatial dispersion, which is driven by constraint on core residues. In contrast, germline disease-causing variants are generally clustered in protein structures and form clusters more frequently than recurrent somatic variants identified from tumor sequencing. In total, we identify 215 proteins with significant spatial constraints on the distribution of disease-causing missense variants in experimentally derived protein structures, only 65 (30%) of which have been previously reported. This analysis identifies many clusters not detectable from sequence information alone; only 12% of proteins with significant clustering in 3D were identified from similar analyses of linear protein sequence. Furthermore, spatial analyses of mutations in homology-based structural models are highly correlated with those from experimentally derived structures, supporting the use of computationally derived models. Our approach highlights significant differences in the spatial constraints on different classes of mutations in protein structure and identifies regions of potential function within individual proteins.

Redefining the phenotypic spectrum of de novo heterozygous CDK13 variants: Three patients without cardiac defects

Recently, 7 patients with de novo constitutional non-synonymous mutations in the CDK13 gene were ascertained through a trio exome analysis of a large cohort of 610 patients with congenital cardiac diseases. Despite another report describing 9 additional patients, the clinical spectrum of this condition has yet to be defined. Herein, we report 3 patients with heterozygous constitutional CDK13 mutations, who were ascertained through exome analysis of children with intellectual disability and minor anomalies, who lacked cardiac anomalies. Two patients had a c.2149G > A, p.Gly717Arg mutation, and one had a c.2525A > G, p. Asn842Ser mutation. A review of the previously described patients and those described herein has enabled the following points to be clarified. First, congenital heart diseases are not an essential feature (13/19). Second, nasal features may help syndromic recognition (14/16). Third, widely spaced and peg-shaped teeth may represent a previously unappreciated diagnostic clue for this newly identified syndrome. Here, we show that p.Gly717Arg represents a hotspot in addition to p.Asn842Ser. We suggest that this CDK13-related disorder may represent a clinically recognizable syndrome.

Posttranslational Modification as a Critical Determinant of Cytoplasmic Innate Immune Recognition

Cell surface innate immune receptors can directly detect a variety of extracellular pathogens to which cytoplasmic innate immune sensors are rarely exposed. Instead, within the cytoplasm, the environment is rife with cellular machinery and signaling pathways that are indirectly perturbed by pathogenic microbes to activate intracellular sensors, such as pyrin, NLRP1, NLRP3, or NLRC4. Therefore, subtle changes in key intracellular processes such as phosphorylation, ubiquitination, and other pathways leading to posttranslational protein modification are key determinants of innate immune recognition in the cytoplasm. This concept is critical to establish the “guard hypothesis” whereby otherwise homeostatic pathways that keep innate immune sensors at bay are released in response to alterations in their posttranslational modification status. Originally identified in plants, evidence that a similar guardlike mechanism exists in humans has recently been identified, whereby a mutation that prevents phosphorylation of the innate immune sensor pyrin triggers a dominantly inherited autoinflammatory disease. It is also noteworthy that even when a cytoplasmic innate immune sensor has a direct ligand, such as bacterial peptidoglycan (NOD1 or NOD2), RNA (RIG-I or MDA5), or DNA (cGAS or IFI16), it can still be influenced by posttranslational modification to dramatically alter its response. Therefore, due to their existence in the cytoplasmic milieu, posttranslational modification is a key determinant of intracellular innate immune receptor functionality.

Genome annotation for clinical genomic diagnostics: strengths and weaknesses

The Human Genome Project and advances in DNA sequencing technologies have revolutionized the identification of genetic disorders through the use of clinical exome sequencing. However, in a considerable number of patients, the genetic basis remains unclear. As clinicians begin to consider whole-genome sequencing, an understanding of the processes and tools involved and the factors to consider in the annotation of the structure and function of genomic elements that might influence variant identification is crucial. Here, we discuss and illustrate the strengths and weaknesses of approaches for the annotation and classification of important elements of protein-coding genes, other genomic elements such as pseudogenes and the non-coding genome, comparative-genomic approaches for inferring gene function, and new technologies for aiding genome annotation, as a practical guide for clinicians when considering pathogenic sequence variation. Complete and accurate annotation of structure and function of genome features has the potential to reduce both false-negative (from missing annotation) and false-positive (from incorrect annotation) errors in causal variant identification in exome and genome sequences. Re-analysis of unsolved cases will be necessary as newer technology improves genome annotation, potentially improving the rate of diagnosis.

Gene and Variant Annotation for Mendelian Disorders in the Era of Advanced Sequencing Technologies

Comprehensive annotations of genetic and noncoding regions and corresponding accurate variant classification for Mendelian diseases are the next big challenge in the new genomic era of personalized medicine. Progress in the development of faster and more accurate pipelines for genome annotation and variant classification will lead to the discovery of more novel disease associations and candidate therapeutic targets. This ultimately will facilitate better patient recruitment in clinical trials. In this review, we describe the trends in research at the intersection of basic and clinical genomics that aims to increase understanding of overall genomic complexity, complex inheritance patterns of disease, and patient-phenotype-specific genomic associations. We describe the emerging field of translational functional genomics, which integrates other functional "-omics" approaches that support next-generation sequencing genomic data in order to facilitate personalized diagnostics, disease management, biomarker discovery, and medicine. We also discuss the utility of this integrated approach for diagnostic clinics and medical databases and its role in the future of personalized medicine.

Racial Differences in Cancer Susceptibility and Survival: More Than the Color of the Skin?

Epidemiological studies point to race as a determining factor in cancer susceptibility. In US registries recording cancer incidence and survival by race (distinguishing 'black versus white'), individuals of African ancestry have a globally increased risk of malignancies compared with Caucasians and Asian Americans. Differences in socioeconomic status and health-care access play a key role. However, the lesser disease susceptibility of Hispanic populations with comparable lifestyles and socioeconomic status as African Americans (Hispanic paradox) points to the concomitant importance of genetic determinants. Here, we overview the molecular basis of racial disparity in cancer susceptibility ranging from genetic polymorphisms and cancer-driver gene mutations to obesity, chronic inflammation, and immune responses. We discuss implications for race-adapted cancer screening programs and clinical trials to reduce disparities in cancer burden.

Glycosylation and other PTMs alterations in neurodegenerative diseases: Current status and future role in neurotrauma

Traumatic brain injuries (TBIs) present a chief public health threat affecting nations worldwide. As numbers of patients afflicted by TBI are expected to rise, the necessity to increase our understanding of the pathophysiological mechanism(s) as a result of TBI mounts. TBI is known to augment the risk of developing a number of neurodegenerative diseases (NDs) such as Alzheimer's disease (AD) and Parkinson's disease (PD). Hence, it is rational to assume that a common mechanistic ground links the pathophysiology of NDs to that of TBIs. Through this review, we aim to identify the protein-protein interactions, differential proteins expression, and PTMs, mainly glycosylation, that are involved in the pathogenesis of both ND and TBI. OVID and PubMed have been rigorously searched to identify studies that utilized advanced proteomic platforms (MS based) and systems biology tools to unfold the mechanism(s) behind ND in an attempt to unveil the mysterious biological processes that occur postinjury. Various PTMs have been found to be common between TBI and AD, whereas no similarities have been found between TBI and PD. Phosphorylated tau protein, glycosylated amyloid precursor protein, and many other modifications appear to be common in both TBI and AD. PTMs, differential protein profiles, and altered biological pathways appear to have critical roles in ND processes by interfering with their pathological condition in a manner similar to TBI. Advancement in glycoproteomic studies pertaining to ND and TBI is urgently needed in order to develop better diagnostic tools, therapies, and more favorable prognoses.

Frequent mutations in acetylation and ubiquitination sites suggest novel driver mechanisms of cancer

Background

Discovery of cancer drivers is a major goal of cancer research. Driver genes and pathways are often predicted using mutation frequency, assuming that statistically significant recurrence of specific somatic mutations across independent samples indicates their importance in cancer. However, many mutations, including known cancer drivers, are not observed at high frequency. Fortunately, abundant information is available about functional “active sites” in proteins that can be integrated with mutations to predict cancer driver genes, even based on low frequency mutations. Further, considering active site information predicts detailed biochemical mechanisms impacted by the mutations. Post-translational modifications (PTMs) are active sites that are regulatory switches in proteins and pathways. We analyzed acetylation and ubiquitination, two important PTM types often involved in chromatin organization and protein degradation, to find proteins that are significantly affected by tumor somatic mutations.

Methods

We performed computational analyses of acetylation and ubiquitination sites in a pan-cancer dataset of 3200 tumor samples from The Cancer Genome Atlas (TCGA). These analyses were targeted at different levels of biological organization including individual genes, pathway annotated gene sets, and protein-protein interaction networks.

Results

Acetylation and ubiquitination site mutations are enriched in cancer with significantly stronger evolutionary conservation and accumulation in protein domains. Gene-focused analysis with the ActiveDriver method reveals significant co-occurrences of acetylation and ubiquitination PTMs and mutation hotspots in known oncoproteins (TP53, AKT1, IDH1) and highlights candidate cancer driver genes with PTM-related mechanisms (e.g. several histone proteins and the splicing factor SF3B1). Pathway analysis shows that PTM mutations in acetylation and ubiquitination sites accumulate in cancer-related processes such as cell cycle, apoptosis, chromatin regulation, and metabolism. Integrated mutation analysis of clinical information and protein interaction networks suggests that many PTM-specific mutations associate with decreased patient survival.

Conclusions

Mutation analysis of acetylation and ubiquitination PTM sites reveals their importance in cancer. As PTM networks are increasingly mapped and related enzymes are often druggable, deeper investigation of specific associated mutations may lead to the discovery of treatment-relevant cellular mechanisms.

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-016-0311-2) contains supplementary material, which is available to authorized users.

g:Profiler-a web server for functional interpretation of gene lists (2016 update)

Functional enrichment analysis is a key step in interpreting gene lists discovered in diverse high-throughput experiments. g:Profiler studies flat and ranked gene lists and finds statistically significant Gene Ontology terms, pathways and other gene function related terms. Translation of hundreds of gene identifiers is another core feature of g:Profiler. Since its first publication in 2007, our web server has become a popular tool of choice among basic and translational researchers. Timeliness is a major advantage of g:Profiler as genome and pathway information is synchronized with the Ensembl database in quarterly updates. g:Profiler supports 213 species including mammals and other vertebrates, plants, insects and fungi. The 2016 update of g:Profiler introduces several novel features. We have added further functional datasets to interpret gene lists, including transcription factor binding site predictions, Mendelian disease annotations, information about protein expression and complexes and gene mappings of human genetic polymorphisms. Besides the interactive web interface, g:Profiler can be accessed in computational pipelines using our R package, Python interface and BioJS component. g:Profiler is freely available at http://biit.cs.ut.ee/gprofiler/.

Prioritizing functional phosphorylation sites based on multiple feature integration

Protein phosphorylation is an important type of post-translational modification that is involved in a variety of biological activities. Most phosphorylation events occur on serine, threonine and tyrosine residues in eukaryotes. In recent years, many phosphorylation sites have been identified as a result of advances in mass-spectrometric techniques. However, a large percentage of phosphorylation sites may be non-functional. Systematically prioritizing functional sites from a large number of phosphorylation sites will be increasingly important for the study of their biological roles. This study focused on exploring the intrinsic features of functional phosphorylation sites to predict whether a phosphosite is likely to be functional. We found significant differences in the distribution of evolutionary conservation, kinase association, disorder score, and secondary structure between known functional and background phosphorylation datasets. We built four different types of classifiers based on the most representative features and found that their performances were similar. We also prioritized 213,837 human phosphorylation sites from a variety of phosphorylation databases, which will be helpful for subsequent functional studies. All predicted results are available for query and download on our website (Predict Functional Phosphosites, PFP, http://pfp.biosino.org/).

g2pDB: A Database Mapping Protein Post-Translational Modifications to Genomic Coordinates

Large scale proteomics have made it possible to broadly screen samples for the presence of many types of post-translational modifications, such as phosphorylation, acetylation, and ubiquitination. This type of data has allowed the localization of these modifications to either a specific site on a proteolytically generated peptide or to within a small domain on the peptide. The resulting modification acceptor sites can then be mapped onto the appropriate protein sequences and the information archived. This paper describes the usage of a very large archive of experimental observations of human post-translational modifications to create a map of the most reproducible modification observations onto the complete set of human protein sequences. This set of modification acceptor sites was then directly translated into the genomic coordinates for the codons for the residues at those sites. We constructed the database g2pDB using this protein-to-codon site mapping information. The information in g2pDB has been made available through a RESTful-style API, allowing researchers to determine which specific protein modifications would be perturbed by a set of observed nucleotide variants determined by high throughput DNA or RNA sequencing.

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.