ECOD domain classification of 48 whole proteomes from AlphaFold Structure Database using DPAM2.

Protein structure prediction has now been deployed widely across several different large protein sets. Large-scale domain annotation of these predictions can aid in the development of biological insights. Using our Evolutionary Classification of Protein Domains (ECOD) from experimental structures as a basis for classification, we describe the detection and cataloging of domains from 48 whole proteomes deposited in the AlphaFold Database. On average, we can provide positive classification (either of domains or other identifiable non-domain regions) for 90% of residues in all proteomes. We classified 746,349 domains from 536,808 proteins comprised of over 226,424,000 amino acid residues. We examine the varying populations of homologous groups in both eukaryotes and bacteria. In addition to containing a higher fraction of disordered regions and unassigned domains, eukaryotes show a higher proportion of repeated proteins, both globular and small repeats. We enumerate those highly populated domains that are shared in both eukaryotes and bacteria, such as the Rossmann domains, TIM barrels, and P-loop domains. Additionally, we compare the sampling of homologous groups from this whole proteome set against our stable ECOD reference and discuss groups that have been enriched by structure predictions. Finally, we discuss the implication of these results for protein target selection for future classification strategies for very large protein sets.

DisEnrich: database of enriched regions in human dark proteome.

MOTIVATION: Intrinsically disordered proteins (IDPs) are involved in numerous processes crucial for living organisms. Bias in amino acid composition of these proteins determines their unique biophysical and functional features. Distinct intrinsically disordered regions (IDRs) with compositional bias play different important roles in various biological processes. IDRs enriched in particular amino acids in human proteome have not been described consistently. RESULTS: We developed DisEnrich-the database of human proteome IDRs that are significantly enriched in particular amino acids. Each human protein is described using Gene Ontology (GO) function terms, disorder prediction for the full-length sequence using three methods, enriched IDR composition and ranks of human proteins with similar enriched IDRs. Distribution analysis of enriched IDRs among broad functional categories revealed significant overrepresentation of R- and Y-enriched IDRs in metabolic and enzymatic activities and F-enriched IDRs in transport. About 75% of functional categories contain IDPs with IDRs significantly enriched in hydrophobic residues that are important for protein-protein interactions. AVAILABILITY AND IMPLEMENTATION: The database is available at http://prodata.swmed.edu/DisEnrichDB/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics Advances online.

Functional analysis of Rossmann-like domains reveals convergent evolution of topology and reaction pathways.

Rossmann folds are ancient, frequently diverged domains found in many biological reaction pathways where they have adapted for different functions. Consequently, discernment and classification of their homologous relations and function can be complicated. We define a minimal Rossmann-like structure motif (RLM) that corresponds for the common core of known Rossmann domains and use this motif to identify all RLM domains in the Protein Data Bank (PDB), thus finding they constitute about 20% of all known 3D structures. The Evolutionary Classification of protein structure Domains (ECOD) classifies RLM domains in a number of groups that lack evidence for homology (X-groups), which suggests that they could have evolved independently multiple times. Closely related, homologous RLM enzyme families can diverge to bind different ligands using similar binding sites and to catalyze different reactions. Conversely, non-homologous RLM domains can converge to catalyze the same reactions or to bind the same ligand with alternate binding modes. We discuss a special case of such convergent evolution that is relevant to the polypharmacology paradigm, wherein the same drug (methotrexate) binds to multiple non-homologous RLM drug targets with different topologies. Finally, assigning proteins with RLM domain to the Enzyme Commission classification suggest that RLM enzymes function mainly in metabolism (and comprise 38% of reference metabolic pathways) and are overrepresented in extant pathways that represent ancient biosynthetic routes such as nucleotide metabolism, energy metabolism, and metabolism of amino acids. In fact, RLM enzymes take part in five out of eight enzymatic reactions of the Wood-Ljungdahl metabolic pathway thought to be used by the last universal common ancestor (LUCA). The prevalence of RLM domains in this ancient metabolism might explain their wide distribution among enzymes.

Molecular dynamics simulations of the Nip7 proteins from the marine deep- and shallow-water Pyrococcus species.

BACKGROUND: The identification of the mechanisms of adaptation of protein structures to extreme environmental conditions is a challenging task of structural biology. We performed molecular dynamics (MD) simulations of the Nip7 protein involved in RNA processing from the shallow-water (P. furiosus) and the deep-water (P. abyssi) marine hyperthermophylic archaea at different temperatures (300 and 373 K) and pressures (0.1, 50 and 100 MPa). The aim was to disclose similarities and differences between the deep- and shallow-sea protein models at different temperatures and pressures. RESULTS: The current results demonstrate that the 3D models of the two proteins at all the examined values of pressures and temperatures are compact, stable and similar to the known crystal structure of the P. abyssi Nip7. The structural deviations and fluctuations in the polypeptide chain during the MD simulations were the most pronounced in the loop regions, their magnitude being larger for the C-terminal domain in both proteins. A number of highly mobile segments the protein globule presumably involved in protein-protein interactions were identified. Regions of the polypeptide chain with significant difference in conformational dynamics between the deep- and shallow-water proteins were identified. CONCLUSIONS: The results of our analysis demonstrated that in the examined ranges of temperatures and pressures, increase in temperature has a stronger effect on change in the dynamic properties of the protein globule than the increase in pressure. The conformational changes of both the deep- and shallow-sea protein models under increasing temperature and pressure are non-uniform. Our current results indicate that amino acid substitutions between shallow- and deep-water proteins only slightly affect overall stability of two proteins. Rather, they may affect the interactions of the Nip7 protein with its protein or RNA partners.

DrugDomain: The evolutionary context of drugs and small molecules bound to domains.

Interactions between proteins and small organic compounds play a crucial role in regulating protein functions. These interactions can modulate various aspects of protein behavior, including enzymatic activity, signaling cascades, and structural stability. By binding to specific sites on proteins, small organic compounds can induce conformational changes, alter protein-protein interactions, or directly affect catalytic activity. Therefore, many drugs available on the market today are small molecules (72% of all approved drugs in the last 5 years). Proteins are composed of one or more domains: evolutionary units that convey function or fitness either singly or in concert with others. Understanding which domain(s) of the target protein binds to a drug can lead to additional opportunities for discovering novel targets. The evolutionary classification of protein domains (ECOD) classifies domains into an evolutionary hierarchy that focuses on distant homology. Previously, no structure-based protein domain classification existed that included information about both the interaction between small molecules or drugs and the structural domains of a target protein. This data is especially important for multidomain proteins and large complexes. Here, we present the DrugDomain database that reports the interaction between ECOD of human target proteins and DrugBank molecules and drugs. The pilot version of DrugDomain describes the interaction of 5160 DrugBank molecules associated with 2573 human proteins. It describes domains for all experimentally determined structures of these proteins and incorporates AlphaFold models when such structures are unavailable. The DrugDomain database is available online: http://prodata.swmed.edu/DrugDomain/.

Deep Learning for Subtypes Identification of Pure Seminoma of the Testis.

The most critical step in the clinical diagnosis workflow is the pathological evaluation of each tumor sample. Deep learning is a powerful approach that is widely used to enhance diagnostic accuracy and streamline the diagnosis process. In our previous study using omics data, we identified 2 distinct subtypes of pure seminoma. Seminoma is the most common histological type of testicular germ cell tumors (TGCTs). Here we developed a deep learning decision making tool for the identification of seminoma subtypes using histopathological slides. We used all available slides for pure seminoma samples from The Cancer Genome Atlas (TCGA). The developed model showed an area under the ROC curve of 0.896. Our model not only confirms the presence of 2 distinct subtypes within pure seminoma but also unveils the presence of morphological differences between them that are imperceptible to the human eye.

Structure classification of the proteins from Salmonella enterica pangenome revealed novel potential pathogenicity islands.

Salmonella enterica is a pathogenic bacterium known for causing severe typhoid fever in humans, making it important to study due to its potential health risks and significant impact on public health. This study provides evolutionary classification of proteins from Salmonella enterica pangenome. We classified 17,238 domains from 13,147 proteins from 79,758 Salmonella enterica strains and studied in detail domains of 272 proteins from 14 characterized Salmonella pathogenicity islands (SPIs). Among SPIs-related proteins, 90 proteins function in the secretion machinery. 41% domains of SPI proteins have no previous sequence annotation. By comparing clinical and environmental isolates, we identified 3682 proteins that are overrepresented in clinical group that we consider as potentially pathogenic. Among domains of potentially pathogenic proteins only 50% domains were annotated by sequence methods previously. Moreover, 36% (1330 out of 3682) of potentially pathogenic proteins cannot be classified into Evolutionary Classification of Protein Domains database (ECOD). Among classified domains of potentially pathogenic proteins the most populated homology groups include helix-turn-helix (HTH), Immunoglobulin-related, and P-loop domains-related. Functional analysis revealed overrepresentation of these protein in biological processes related to viral entry into host cell, antibiotic biosynthesis, DNA metabolism and conformation change, and underrepresentation in translational processes. Analysis of the potentially pathogenic proteins indicates that they form 119 clusters or novel potential pathogenicity islands (NPPIs) within the Salmonella genome, suggesting their potential contribution to the bacterium's virulence. One of the NPPIs revealed significant overrepresentation of potentially pathogenic proteins. Overall, our analysis revealed that identified potentially pathogenic proteins are poorly studied.

Seminoma subtypes differ in the organization and functional state of the immune microenvironment.

Seminoma is the most common type of testicular germ cell tumors (TGCTs) among 15-44 years old men. Seminoma treatments include orchiectomy, platinum-based chemotherapy and radiotherapy. These radical treatment methods cause up to 40 severe adverse long-term side effects including secondary cancers. Immunotherapy based on immune checkpoint inhibitors, which showed its efficiency for many types of cancer, can be important alternative to the platinum-based therapy for seminoma patients. However, five independent clinical trials evaluating the efficiency of immune checkpoint inhibitors for TGCTs treatment were shut down at the phase II due to lacking clinical efficacy and detailed mechanisms of this phenomena are yet to be discovered. Recently we identified two distinct seminoma subtypes based on transcriptomic data and here we focused on the analysis of seminoma microenvironment and its subtype-specific characteristics. Our analysis revealed that less differentiated subtype 1 of seminoma has immune microenvironment with significantly lower immune score and larger fraction of neutrophils. Both are features of the immune microenvironment at an early developmental stage. On the contrary, subtype 2 seminoma is characterized by the higher immune score and overexpression of 21 genes related to senescence-associated secretory phenotype. Seminoma single cell transcriptomic data showed that 9 out of 21 genes are predominantly expressed in immune cells. Therefore, we hypothesized that senescence of immune microenvironment can be one of the reasons for seminoma immunotherapy failure. Supplementary Information: The online version contains supplementary material available at 10.1007/s13205-023-03530-1.

Pan-cancer structurome reveals overrepresentation of beta sandwiches and underrepresentation of alpha helical domains.

The recent progress in the prediction of protein structures marked a historical milestone. AlphaFold predicted 200 million protein models with an accuracy comparable to experimental methods. Protein structures are widely used to understand evolution and to identify potential drug targets for the treatment of various diseases, including cancer. Thus, these recently predicted structures might convey previously unavailable information about cancer biology. Evolutionary classification of protein domains is challenging and different approaches exist. Recently our team presented a classification of domains from human protein models released by AlphaFold. Here we evaluated the pan-cancer structurome, domains from over and under expressed proteins in 21 cancer types, using the broadest levels of the ECOD classification: the architecture (A-groups) and possible homology (X-groups) levels. Our analysis reveals that AlphaFold has greatly increased the three-dimensional structural landscape for proteins that are differentially expressed in these 21 cancer types. We show that beta sandwich domains are significantly overrepresented and alpha helical domains are significantly underrepresented in the majority of cancer types. Our data suggest that the prevalence of the beta sandwiches is due to the high levels of immunoglobulins and immunoglobulin-like domains that arise during tumor development-related inflammation. On the other hand, proteins with exclusively alpha domains are important elements of homeostasis, apoptosis and transmembrane transport. Therefore cancer cells tend to reduce representation of these proteins to promote successful oncogeneses.

Pathogenic mutation hotspots in protein kinase domain structure.

Control of eukaryotic cellular function is heavily reliant on the phosphorylation of proteins at specific amino acid residues, such as serine, threonine, tyrosine, and histidine. Protein kinases that are responsible for this process comprise one of the largest families of evolutionarily related proteins. Dysregulation of protein kinase signaling pathways is a frequent cause of a large variety of human diseases including cancer, autoimmune, neurodegenerative, and cardiovascular disorders. In this study, we mapped all pathogenic mutations in 497 human protein kinase domains from the ClinVar database to the reference structure of Aurora kinase A (AURKA) and grouped them by the relevance to the disease type. Our study revealed that the majority of mutation hotspots associated with cancer are situated within the catalytic and activation loops of the kinase domain, whereas non-cancer-related hotspots tend to be located outside of these regions. Additionally, we identified a hotspot at residue R371 of the AURKA structure that has the highest number of exclusively non-cancer-related pathogenic mutations (21) and has not been previously discussed.

Integrated Molecular Analysis Reveals 2 Distinct Subtypes of Pure Seminoma of the Testis.

Objective: Testicular germ cell tumors (TGCT) are the most common solid malignancy in adolescent and young men, with a rising incidence over the past 20 years. Overall, TGCTs are second in terms of the average life years lost per person dying of cancer, and clinical therapeutics without adverse long-term side effects are lacking. Platinum-based regimens for TGCTs have heterogeneous outcomes even within the same histotype that frequently leads to under- and over-treatment. Understanding of molecular differences that lead to diverse outcomes of TGCT patients may improve current treatment approaches. Seminoma is the most common subtype of TGCTs, which can either be pure or present in combination with other histotypes. Methods: Here we conducted a computational study of 64 pure seminoma samples from The Cancer Genome Atlas, applied consensus clustering approach to their transcriptomic data and revealed 2 clinically relevant seminoma subtypes: seminoma subtype 1 and 2. Results: Our analysis identified significant differences in pluripotency stage, activity of double stranded DNA breaks repair mechanisms, rates of loss of heterozygosity, and expression of lncRNA responsible for cisplatin resistance between the subtypes. Seminoma subtype 1 is characterized by higher pluripotency state, while subtype 2 showed attributes of reprograming into non-seminomatous TGCT. The seminoma subtypes we identified may provide a molecular underpinning for variable responses to chemotherapy and radiation. Conclusion: Translating our findings into clinical care may help improve risk stratification of seminoma, decrease overtreatment rates, and increase long-term quality of life for TGCT survivors.

A Fifth of the Protein World: Rossmann-like Proteins as an Evolutionarily Successful Structural unit.

The Rossmann-like fold is the most prevalent and diversified doubly-wound superfold of ancient evolutionary origin. Rossmann-like domains are present in a variety of metabolic enzymes and are capable of binding diverse ligands. Discerning evolutionary relationships among these domains is challenging because of their diverse functions and ancient origin. We defined a minimal Rossmann-like structural motif (RLM), identified RLM-containing domains among known 3D structures (20%) and classified them according to their homologous relationships. New classifications were incorporated into our Evolutionary Classification of protein Domains (ECOD) database. We defined 156 homology groups (H-groups), which were further clustered into 123 possible homology groups (X-groups). Our analysis revealed that RLM-containing proteins constitute approximately 15% of the human proteome. We found that disease-causing mutations are more frequent within RLM domains than within non-RLM domains of these proteins, highlighting the importance of RLM-containing proteins for human health.

Completeness and Consistency in Structural Domain Classifications.

Domain classifications are a useful resource for computational analysis of the protein structure, but elements of their composition are often opaque to potential users. We perform a comparative analysis of our classification ECOD against the SCOPe, SCOP2, and CATH domain classifications with respect to their constituent domain boundaries and hierarchal organization. The coverage of these domain classifications with respect to ECOD and to the PDB was assessed by structure and by sequence. We also conducted domain pair analysis to determine broad differences in hierarchy between domains shared by ECOD and other classifications. Finally, we present domains from the major facilitator superfamily (MFS) of transporter proteins and provide evidence that supports their split into domains and for multiple conformations within these families. We find that the ECOD and CATH provide the most extensive structural coverage of the PDB. ECOD and SCOPe have the most consistent domain boundary conditions, whereas CATH and SCOP2 both differ significantly.

Recent advances suggest increased influence of selective pressure in allostery.

Allosteric regulation of protein functions is ubiquitous in organismal biology, but the principles governing its evolution are not well understood. Here we discuss recent studies supporting the large-scale existence of latent allostery in ancestor proteins of superfamilies. As suggested, the evolution of allostery could be driven by the need for specificity in paralogs of slow evolving protein complexes with conserved active sites. The same slow evolution is displayed by purifying selection exhibited in allosteric proteins with somatic mutations involved in cancer, where disease-associated mutations are enriched in both orthosteric and allosteric sites. Consequently, disease-associated variants can be used to identify druggable allosteric sites that are specific for paralogs in protein superfamilies with otherwise similar functions.

ECOD: identification of distant homology among multidomain and transmembrane domain proteins.

The manual classification of protein domains is approaching its 20th anniversary. ECOD is our mixed manual-automatic domain classification. Over time, the types of proteins which require manual curation has changed. Depositions with complex multidomain and multichain arrangements are commonplace. Transmembrane domains are regularly classified. Repeatedly, domains which are initially believed to be novel are found to have homologous links to existing classified domains. Here we present a brief summary of recent manual curation efforts in ECOD generally combined with specific case studies of transmembrane and multidomain proteins wherein manual curation was useful for discovering new homologous relationships. We present a new taxonomy for the classification of ABC transporter transmembrane domains. We examine alternate topologies of the leucine-specific (LS) domain of Leucine tRNA-synthetase. Finally, we elaborate on a distant homologous links between two helical dimerization domains.

High temperature and pressure influence the interdomain orientation of Nip7 proteins from P. abyssi and P. furiosus: MD simulations.

Interactions between protein domains and their position and movement relative to each other are essential for the stability and normal functioning of a protein molecule. Features of the movement of domains may define the mechanism of enzymatic reactions. Therefore, the description of this motion is an important task in the analysis of the structures and functions of multidomain proteins. In the current work, we investigated the influence of pressure and temperature on changes in the movement of the two domains of the protein Nip7, expressed by deep-water (Pyrococcus abyssi) and shallow-water (Pyrococcus furiosus) archaea. The results of the present study show that the interdomain interfaces of the Nip7 proteins of P. abyssi and P. furiosus are formed by stable hydrophobic interactions. It was shown that high pressure and high temperature significantly changed the orientation of domains in Nip7 proteins which perhaps was connected with functional features of these domains. It was found that increasing the pressure significantly changed the angle of rotation of these domains, to a greater extent in the shallow-water protein, while an increase in temperature slightly reduced the angle of rotation of these domains. Moreover, the results suggest that the type of motion of the domains under study is similar to shear motion.

Functional and evolutionary analysis of viral proteins containing a Rossmann-like fold.

Viruses are the most abundant life form and infect practically all organisms. Consequently, these obligate parasites are a major cause of human suffering and economic loss. Rossmann-like fold is the most populated fold among α/ß-folds in the Protein Data Bank and proteins containing Rossmann-like fold constitute 22% of all known proteins 3D structures. Thus, analysis of viral proteins containing Rossmann-like domains could provide an understanding of viral biology and evolution as well as could propose possible targets for antiviral therapy. We provide functional and evolutionary analysis of viral proteins containing a Rossmann-like fold found in the evolutionary classification of protein domains (ECOD) database developed in our lab. We identified 81 protein families of bacterial, archeal, and eukaryotic viruses in light of their evolution-based ECOD classification and Pfam taxonomy. We defined their functional significance using enzymatic EC number assignments as well as domain-level family annotations.

Identification of residues of the archaeal RNA-binding Nip7 proteins specific to environmental conditions.

The understanding of biological and molecular mechanisms providing survival of cells under extreme temperatures and pressures will help to answer fundamental questions related to the origin of life and to design of biotechnologically important enzymes with new properties. Here, we analyze amino acid sequences of the Nip7 proteins from 35 archaeal species to identify positions containing mutations specific to the hydrostatic pressure and temperature of organism's habitat. The number of such positions related to pressure change is much lower than related to temperature change. The results suggest that adaptation to temperature changes of the Nip7 protein cause more pronounced modifications in sequence and structure, than to the pressure changes. Structural analysis of residues at these positions demonstrated their involvement in salt-bridge formation, which may reflect the importance of protein structure stabilization by salt-bridges at extreme environmental conditions.