De novo protein identification in mammalian sperm using in situ cryoelectron tomography and AlphaFold2 docking.

To understand the molecular mechanisms of cellular pathways, contemporary workflows typically require multiple techniques to identify proteins, track their localization, and determine their structures in vitro. Here, we combined cellular cryoelectron tomography (cryo-ET) and AlphaFold2 modeling to address these questions and understand how mammalian sperm are built in situ. Our cellular cryo-ET and subtomogram averaging provided 6.0-Å reconstructions of axonemal microtubule structures. The well-resolved tertiary structures allowed us to unbiasedly match sperm-specific densities with 21,615 AlphaFold2-predicted protein models of the mouse proteome. We identified Tektin 5, CCDC105, and SPACA9 as novel microtubule-associated proteins. These proteins form an extensive interaction network crosslinking the lumen of axonemal doublet microtubules, suggesting their roles in modulating the mechanical properties of the filaments. Indeed, Tekt5 -/- sperm possess more deformed flagella with 180° bends. Together, our studies presented a cellular visual proteomics workflow and shed light on the in vivo functions of Tektin 5.

SARS-CoV-2 variants evolve convergent strategies to remodel the host response.

SARS-CoV-2 variants of concern (VOCs) emerged during the COVID-19 pandemic. Here, we used unbiased systems approaches to study the host-selective forces driving VOC evolution. We discovered that VOCs evolved convergent strategies to remodel the host by modulating viral RNA and protein levels, altering viral and host protein phosphorylation, and rewiring virus-host protein-protein interactions. Integrative computational analyses revealed that although Alpha, Beta, Gamma, and Delta ultimately converged to suppress interferon-stimulated genes (ISGs), Omicron BA.1 did not. ISG suppression correlated with the expression of viral innate immune antagonist proteins, including Orf6, N, and Orf9b, which we mapped to specific mutations. Later Omicron subvariants BA.4 and BA.5 more potently suppressed innate immunity than early subvariant BA.1, which correlated with Orf6 levels, although muted in BA.4 by a mutation that disrupts the Orf6-nuclear pore interaction. Our findings suggest that SARS-CoV-2 convergent evolution overcame human adaptive and innate immune barriers, laying the groundwork to tackle future pandemics.

The Global Phosphorylation Landscape of SARS-CoV-2 Infection.

The causative agent of the coronavirus disease 2019 (COVID-19) pandemic, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has infected millions and killed hundreds of thousands of people worldwide, highlighting an urgent need to develop antiviral therapies. Here we present a quantitative mass spectrometry-based phosphoproteomics survey of SARS-CoV-2 infection in Vero E6 cells, revealing dramatic rewiring of phosphorylation on host and viral proteins. SARS-CoV-2 infection promoted casein kinase II (CK2) and p38 MAPK activation, production of diverse cytokines, and shutdown of mitotic kinases, resulting in cell cycle arrest. Infection also stimulated a marked induction of CK2-containing filopodial protrusions possessing budding viral particles. Eighty-seven drugs and compounds were identified by mapping global phosphorylation profiles to dysregulated kinases and pathways. We found pharmacologic inhibition of the p38, CK2, CDK, AXL, and PIKFYVE kinases to possess antiviral efficacy, representing potential COVID-19 therapies.

Brahma safeguards canalization of cardiac mesoderm differentiation.

Differentiation proceeds along a continuum of increasingly fate-restricted intermediates, referred to as canalization1,2. Canalization is essential for stabilizing cell fate, but the mechanisms that underlie robust canalization are unclear. Here we show that the BRG1/BRM-associated factor (BAF) chromatin-remodelling complex ATPase gene Brm safeguards cell identity during directed cardiogenesis of mouse embryonic stem cells. Despite the establishment of a well-differentiated precardiac mesoderm, Brm-/- cells predominantly became neural precursors, violating germ layer assignment. Trajectory inference showed a sudden acquisition of a non-mesodermal identity in Brm-/- cells. Mechanistically, the loss of Brm prevented de novo accessibility of primed cardiac enhancers while increasing the expression of neurogenic factor POU3F1, preventing the binding of the neural suppressor REST and shifting the composition of BRG1 complexes. The identity switch caused by the Brm mutation was overcome by increasing BMP4 levels during mesoderm induction. Mathematical modelling supports these observations and demonstrates that Brm deletion affects cell fate trajectory by modifying saddle-node bifurcations2. In the mouse embryo, Brm deletion exacerbated mesoderm-deleted Brg1-mutant phenotypes, severely compromising cardiogenesis, and reveals an in vivo role for Brm. Our results show that Brm is a compensable safeguard of the fidelity of mesoderm chromatin states, and support a model in which developmental canalization is not a rigid irreversible path, but a highly plastic trajectory.

Profiling the proximal proteome of the activated µ-opioid receptor.

The µ-opioid receptor (µOR) represents an important target of therapeutic and abused drugs. So far, most understanding of µOR activity has focused on a subset of known signal transducers and regulatory molecules. Yet µOR signaling is coordinated by additional proteins in the interaction network of the activated receptor, which have largely remained invisible given the lack of technologies to interrogate these networks systematically. Here we describe a proteomics and computational approach to map the proximal proteome of the activated µOR and to extract subcellular location, trafficking and functional partners of G-protein-coupled receptor (GPCR) activity. We demonstrate that distinct opioid agonists exert differences in the µOR proximal proteome mediated by endocytosis and endosomal sorting. Moreover, we identify two new µOR network components, EYA4 and KCTD12, which are recruited on the basis of receptor-triggered G-protein activation and might form a previously unrecognized buffering system for G-protein activity broadly modulating cellular GPCR signaling.

A SARS-CoV-2 protein interaction map reveals targets for drug repurposing.

A newly described coronavirus named severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which is the causative agent of coronavirus disease 2019 (COVID-19), has infected over 2.3 million people, led to the death of more than 160,000 individuals and caused worldwide social and economic disruption1,2. There are no antiviral drugs with proven clinical efficacy for the treatment of COVID-19, nor are there any vaccines that prevent infection with SARS-CoV-2, and efforts to develop drugs and vaccines are hampered by the limited knowledge of the molecular details of how SARS-CoV-2 infects cells. Here we cloned, tagged and expressed 26 of the 29 SARS-CoV-2 proteins in human cells and identified the human proteins that physically associated with each of the SARS-CoV-2 proteins using affinity-purification mass spectrometry, identifying 332 high-confidence protein-protein interactions between SARS-CoV-2 and human proteins. Among these, we identify 66 druggable human proteins or host factors targeted by 69 compounds (of which, 29 drugs are approved by the US Food and Drug Administration, 12 are in clinical trials and 28 are preclinical compounds). We screened a subset of these in multiple viral assays and found two sets of pharmacological agents that displayed antiviral activity: inhibitors of mRNA translation and predicted regulators of the sigma-1 and sigma-2 receptors. Further studies of these host-factor-targeting agents, including their combination with drugs that directly target viral enzymes, could lead to a therapeutic regimen to treat COVID-19.

Protein Interaction Map of APOBEC3 Enzyme Family Reveals Deamination-Independent Role in Cellular Function.

Human APOBEC3 enzymes are a family of single-stranded (ss)DNA and RNA cytidine deaminases that act as part of the intrinsic immunity against viruses and retroelements. These enzymes deaminate cytosine to form uracil which can functionally inactivate or cause degradation of viral or retroelement genomes. In addition, APOBEC3s have deamination-independent antiviral activity through protein and nucleic acid interactions. If expression levels are misregulated, some APOBEC3 enzymes can access the human genome leading to deamination and mutagenesis, contributing to cancer initiation and evolution. While APOBEC3 enzymes are known to interact with large ribonucleoprotein complexes, the function and RNA dependence are not entirely understood. To further understand their cellular roles, we determined by affinity purification mass spectrometry (AP-MS) the protein interaction network for the human APOBEC3 enzymes and mapped a diverse set of protein-protein and protein-RNA mediated interactions. Our analysis identified novel RNA-mediated interactions between APOBEC3C, APOBEC3H Haplotype I and II, and APOBEC3G with spliceosome proteins, and APOBEC3G and APOBEC3H Haplotype I with proteins involved in tRNA methylation and ncRNA export from the nucleus. In addition, we identified RNA-independent protein-protein interactions with APOBEC3B, APOBEC3D, and APOBEC3F and the prefoldin family of protein-folding chaperones. Interaction between prefoldin 5 (PFD5) and APOBEC3B disrupted the ability of PFD5 to induce degradation of the oncogene cMyc, implicating the APOBEC3B protein interaction network in cancer. Altogether, the results uncover novel functions and interactions of the APOBEC3 family and suggest they may have fundamental roles in cellular RNA biology, their protein-protein interactions are not redundant, and there are protein-protein interactions with tumor suppressors, suggesting a role in cancer biology. Data are available via ProteomeXchange with the identifier PXD044275.

A proximity proteomics pipeline with improved reproducibility and throughput.

Proximity labeling (PL) via biotinylation coupled with mass spectrometry (MS) captures spatial proteomes in cells. Large-scale processing requires a workflow minimizing hands-on time and enhancing quantitative reproducibility. We introduced a scalable PL pipeline integrating automated enrichment of biotinylated proteins in a 96-well plate format. Combining this with optimized quantitative MS based on data-independent acquisition (DIA), we increased sample throughput and improved protein identification and quantification reproducibility. We applied this pipeline to delineate subcellular proteomes across various compartments. Using the 5HT2A serotonin receptor as a model, we studied temporal changes of proximal interaction networks induced by receptor activation. In addition, we modified the pipeline for reduced sample input to accommodate CRISPR-based gene knockout, assessing dynamics of the 5HT2A network in response to perturbation of selected interactors. This PL approach is universally applicable to PL proteomics using biotinylation-based PL enzymes, enhancing throughput and reproducibility of standard protocols.

Quantitative Proteomic Analysis Reveals apoE4-Dependent Phosphorylation of the Actin-Regulating Protein VASP.

Apolipoprotein (apo) E4 is the major genetic risk factor for Alzheimer's disease. While neurons generally produce a minority of the apoE in the central nervous system, neuronal expression of apoE increases dramatically in response to stress and is sufficient to drive pathology. Currently, the molecular mechanisms of how apoE4 expression may regulate pathology are not fully understood. Here, we expand upon our previous studies measuring the impact of apoE4 on protein abundance to include the analysis of protein phosphorylation and ubiquitylation signaling in isogenic Neuro-2a cells expressing apoE3 or apoE4. ApoE4 expression resulted in a dramatic increase in vasodilator-stimulated phosphoprotein (VASP) S235 phosphorylation in a protein kinase A (PKA)-dependent manner. This phosphorylation disrupted VASP interactions with numerous actin cytoskeletal and microtubular proteins. Reduction of VASP S235 phosphorylation via PKA inhibition resulted in a significant increase in filopodia formation and neurite outgrowth in apoE4-expressing cells, exceeding levels observed in apoE3-expressing cells. Our results highlight the pronounced and diverse impact of apoE4 on multiple modes of protein regulation and identify protein targets to restore apoE4-related cytoskeletal defects.

Data-Independent Acquisition Protease-Multiplexing Enables Increased Proteome Sequence Coverage Across Multiple Fragmentation Modes.

The use of multiple proteases has been shown to increase protein sequence coverage in proteomics experiments; however, due to the additional analysis time required, it has not been widely adopted in routine data-dependent acquisition (DDA) proteomic workflows. Alternatively, data-independent acquisition (DIA) has the potential to analyze multiplexed samples from different protease digests, but has been primarily optimized for fragmenting tryptic peptides. Here we evaluate a DIA multiplexing approach that combines three proteolytic digests (Trypsin, AspN, and GluC) into a single sample. We first optimize data acquisition conditions for each protease individually with both the canonical DIA fragmentation mode (beam type CID), as well as resonance excitation CID, to determine optimal consensus conditions across proteases. Next, we demonstrate that application of these conditions to a protease-multiplexed sample of human peptides results in similar protein identifications and quantitative performance as compared to trypsin alone, but enables up to a 63% increase in peptide detections, and a 45% increase in nonredundant amino acid detections. Nontryptic peptides enabled noncanonical protein isoform determination and resulted in 100% sequence coverage for numerous proteins, suggesting the utility of this approach in applications where sequence coverage is critical, such as protein isoform analysis.

A global view of structure-function relationships in the tautomerase superfamily.

The tautomerase superfamily (TSF) consists of more than 11,000 nonredundant sequences present throughout the biosphere. Characterized members have attracted much attention because of the unusual and key catalytic role of an N-terminal proline. These few characterized members catalyze a diverse range of chemical reactions, but the full scale of their chemical capabilities and biological functions remains unknown. To gain new insight into TSF structure-function relationships, we performed a global analysis of similarities across the entire superfamily and computed a sequence similarity network to guide classification into distinct subgroups. Our results indicate that TSF members are found in all domains of life, with most being present in bacteria. The eukaryotic members of the cis-3-chloroacrylic acid dehalogenase subgroup are limited to fungal species, whereas the macrophage migration inhibitory factor subgroup has wide eukaryotic representation (including mammals). Unexpectedly, we found that 346 TSF sequences lack Pro-1, of which 85% are present in the malonate semialdehyde decarboxylase subgroup. The computed network also enabled the identification of similarity paths, namely sequences that link functionally diverse subgroups and exhibit transitional structural features that may help explain reaction divergence. A structure-guided comparison of these linker proteins identified conserved transitions between them, and kinetic analysis paralleled these observations. Phylogenetic reconstruction of the linker set was consistent with these findings. Our results also suggest that contemporary TSF members may have evolved from a short 4-oxalocrotonate tautomerase-like ancestor followed by gene duplication and fusion. Our new linker-guided strategy can be used to enrich the discovery of sequence/structure/function transitions in other enzyme superfamilies.

Organic and inorganic mercurials have distinct effects on cellular thiols, metal homeostasis, and Fe-binding proteins in Escherichia coli.

The protean chemical properties of the toxic metal mercury (Hg) have made it attractive in diverse applications since antiquity. However, growing public concern has led to an international agreement to decrease its impact on health and the environment. During a recent proteomics study of acute Hg exposure in E. coli, we also examined the effects of inorganic and organic Hg compounds on thiol and metal homeostases. On brief exposure, lower concentrations of divalent inorganic mercury Hg(II) blocked bulk cellular thiols and protein-associated thiols more completely than higher concentrations of monovalent organomercurials, phenylmercuric acetate (PMA) and merthiolate (MT). Cells bound Hg(II) and PMA in excess of their available thiol ligands; X-ray absorption spectroscopy indicated nitrogens as likely additional ligands. The mercurials released protein-bound iron (Fe) more effectively than common organic oxidants and all disturbed the Na(+)/K(+) electrolyte balance, but none provoked efflux of six essential transition metals including Fe. PMA and MT made stable cysteine monothiol adducts in many Fe-binding proteins, but stable Hg(II) adducts were only seen in CysXxx(n)Cys peptides. We conclude that on acute exposure: (a) the distinct effects of mercurials on thiol and Fe homeostases reflected their different uptake and valences; (b) their similar effects on essential metal and electrolyte homeostases reflected the energy dependence of these processes; and (c) peptide phenylmercury-adducts were more stable or detectable in mass spectrometry than Hg(II)-adducts. These first in vivo observations in a well-defined model organism reveal differences upon acute exposure to inorganic and organic mercurials that may underlie their distinct toxicology.

Protein interaction map of APOBEC3 enzyme family reveals deamination-independent role in cellular function.

Human APOBEC3 enzymes are a family of single-stranded (ss)DNA and RNA cytidine deaminases that act as part of the intrinsic immunity against viruses and retroelements. These enzymes deaminate cytosine to form uracil which can functionally inactivate or cause degradation of viral or retroelement genomes. In addition, APOBEC3s have deamination independent antiviral activity through protein and nucleic acid interactions. If expression levels are misregulated, some APOBEC3 enzymes can access the human genome leading to deamination and mutagenesis, contributing to cancer initiation and evolution. While APOBEC3 enzymes are known to interact with large ribonucleoprotein complexes, the function and RNA dependence is not entirely understood. To further understand their cellular roles, we determined by affinity purification mass spectrometry (AP-MS) the protein interaction network for the human APOBEC3 enzymes and map a diverse set of protein-protein and protein-RNA mediated interactions. Our analysis identified novel RNA-mediated interactions between APOBEC3C, APOBEC3H Haplotype I and II, and APOBEC3G with spliceosome proteins, and APOBEC3G and APOBEC3H Haplotype I with proteins involved in tRNA methylation and ncRNA export from the nucleus. In addition, we identified RNA-independent protein-protein interactions with APOBEC3B, APOBEC3D, and APOBEC3F and the prefoldin family of protein folding chaperones. Interaction between prefoldin 5 (PFD5) and APOBEC3B disrupted the ability of PFD5 to induce degradation of the oncogene cMyc, implicating the APOBEC3B protein interaction network in cancer. Altogether, the results uncover novel functions and interactions of the APOBEC3 family and suggest they may have fundamental roles in cellular RNA biology, their protein-protein interactions are not redundant, and there are protein-protein interactions with tumor suppressors, suggesting a role in cancer biology.

A foundational atlas of autism protein interactions reveals molecular convergence.

Translating high-confidence (hc) autism spectrum disorder (ASD) genes into viable treatment targets remains elusive. We constructed a foundational protein-protein interaction (PPI) network in HEK293T cells involving 100 hcASD risk genes, revealing over 1,800 PPIs (87% novel). Interactors, expressed in the human brain and enriched for ASD but not schizophrenia genetic risk, converged on protein complexes involved in neurogenesis, tubulin biology, transcriptional regulation, and chromatin modification. A PPI map of 54 patient-derived missense variants identified differential physical interactions, and we leveraged AlphaFold-Multimer predictions to prioritize direct PPIs and specific variants for interrogation in Xenopus tropicalis and human forebrain organoids. A mutation in the transcription factor FOXP1 led to reconfiguration of DNA binding sites and altered development of deep cortical layer neurons in forebrain organoids. This work offers new insights into molecular mechanisms underlying ASD and describes a powerful platform to develop and test therapeutic strategies for many genetically-defined conditions.

Discovering mercury protein modifications in whole proteomes using natural isotope distributions observed in liquid chromatography-tandem mass spectrometry.

The identification of peptides that result from post-translational modifications is critical for understanding normal pathways of cellular regulation as well as identifying damage from, or exposures to xenobiotics, i.e. the exposome. However, because of their low abundance in proteomes, effective detection of modified peptides by mass spectrometry (MS) typically requires enrichment to eliminate false identifications. We present a new method for confidently identifying peptides with mercury (Hg)-containing adducts that is based on the influence of mercury's seven stable isotopes on peptide isotope distributions detected by high-resolution MS. Using a pure protein and E. coli cultures exposed to phenyl mercuric acetate, we show the pattern of peak heights in isotope distributions from primary MS single scans efficiently identified Hg adducts in data from chromatographic separation coupled with tandem mass spectrometry with sensitivity and specificity greater than 90%. Isotope distributions are independent of peptide identifications based on peptide fragmentation (e.g. by SEQUEST), so both methods can be combined to eliminate false positives. Summing peptide isotope distributions across multiple scans improved specificity to 99.4% and sensitivity above 95%, affording identification of an unexpected Hg modification. We also illustrate the theoretical applicability of the method for detection of several less common elements including the essential element, selenium, as selenocysteine in peptides.

An automated proximity proteomics pipeline for subcellular proteome and protein interaction mapping.

Proximity labeling (PL) coupled with mass spectrometry has emerged as a powerful technique to map proximal protein interactions in living cells. Large-scale sample processing for proximity proteomics necessitates a high-throughput workflow to reduce hands-on time and increase quantitative reproducibility. To address this issue, we developed a scalable and automated PL pipeline, including generation and characterization of monoclonal cell lines, automated enrichment of biotinylated proteins in a 96-well format, and optimization of the quantitative mass spectrometry (MS) acquisition method. Combined with data-independent acquisition (DIA) MS, our pipeline outperforms manual enrichment and data-dependent acquisition (DDA) MS regarding reproducibility of protein identification and quantification. We apply the pipeline to map subcellular proteomes for endosomes, late endosomes/lysosomes, the Golgi apparatus, and the plasma membrane. Moreover, using serotonin receptor (5HT2A) as a model, we investigated agonist-induced dynamics in protein-protein interactions. Importantly, the approach presented here is universally applicable for PL proteomics using all biotinylation-based PL enzymes, increasing both throughput and reproducibility of standard protocols.

Impact of SARS-CoV-2 ORF6 and its variant polymorphisms on host responses and viral pathogenesis.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) encodes several proteins that inhibit host interferon responses. Among these, ORF6 antagonizes interferon signaling by disrupting nucleocytoplasmic trafficking through interactions with the nuclear pore complex components Nup98-Rae1. However, the roles and contributions of ORF6 during physiological infection remain unexplored. We assessed the role of ORF6 during infection using recombinant viruses carrying a deletion or loss-of-function (LoF) mutation in ORF6. ORF6 plays key roles in interferon antagonism and viral pathogenesis by interfering with nuclear import and specifically the translocation of IRF and STAT transcription factors. Additionally, ORF6 inhibits cellular mRNA export, resulting in the remodeling of the host cell proteome, and regulates viral protein expression. Interestingly, the ORF6:D61L mutation that emerged in the Omicron BA.2 and BA.4 variants exhibits reduced interactions with Nup98-Rae1 and consequently impairs immune evasion. Our findings highlight the role of ORF6 in antagonizing innate immunity and emphasize the importance of studying the immune evasion strategies of SARS-CoV-2.

Proteomic and genetic analyses of influenza A viruses identify pan-viral host targets.

Influenza A Virus (IAV) is a recurring respiratory virus with limited availability of antiviral therapies. Understanding host proteins essential for IAV infection can identify targets for alternative host-directed therapies (HDTs). Using affinity purification-mass spectrometry and global phosphoproteomic and protein abundance analyses using three IAV strains (pH1N1, H3N2, H5N1) in three human cell types (A549, NHBE, THP-1), we map 332 IAV-human protein-protein interactions and identify 13 IAV-modulated kinases. Whole exome sequencing of patients who experienced severe influenza reveals several genes, including scaffold protein AHNAK, with predicted loss-of-function variants that are also identified in our proteomic analyses. Of our identified host factors, 54 significantly alter IAV infection upon siRNA knockdown, and two factors, AHNAK and coatomer subunit COPB1, are also essential for productive infection by SARS-CoV-2. Finally, 16 compounds targeting our identified host factors suppress IAV replication, with two targeting CDK2 and FLT3 showing pan-antiviral activity across influenza and coronavirus families. This study provides a comprehensive network model of IAV infection in human cells, identifying functional host targets for pan-viral HDT.

Impact of SARS-CoV-2 ORF6 and its variant polymorphisms on host responses and viral pathogenesis.

We and others have previously shown that the SARS-CoV-2 accessory protein ORF6 is a powerful antagonist of the interferon (IFN) signaling pathway by directly interacting with Nup98-Rae1 at the nuclear pore complex (NPC) and disrupting bidirectional nucleo-cytoplasmic trafficking. In this study, we further assessed the role of ORF6 during infection using recombinant SARS-CoV-2 viruses carrying either a deletion or a well characterized M58R loss-of-function mutation in ORF6. We show that ORF6 plays a key role in the antagonism of IFN signaling and in viral pathogenesis by interfering with karyopherin(importin)-mediated nuclear import during SARS-CoV-2 infection both in vitro , and in the Syrian golden hamster model in vivo . In addition, we found that ORF6-Nup98 interaction also contributes to inhibition of cellular mRNA export during SARS-CoV-2 infection. As a result, ORF6 expression significantly remodels the host cell proteome upon infection. Importantly, we also unravel a previously unrecognized function of ORF6 in the modulation of viral protein expression, which is independent of its function at the nuclear pore. Lastly, we characterized the ORF6 D61L mutation that recently emerged in Omicron BA.2 and BA.4 and demonstrated that it is able to disrupt ORF6 protein functions at the NPC and to impair SARS-CoV-2 innate immune evasion strategies. Importantly, the now more abundant Omicron BA.5 lacks this loss-of-function polymorphism in ORF6. Altogether, our findings not only further highlight the key role of ORF6 in the antagonism of the antiviral innate immune response, but also emphasize the importance of studying the role of non-spike mutations to better understand the mechanisms governing differential pathogenicity and immune evasion strategies of SARS-CoV-2 and its evolving variants. ONE SENTENCE SUMMARY: SARS-CoV-2 ORF6 subverts bidirectional nucleo-cytoplasmic trafficking to inhibit host gene expression and contribute to viral pathogenesis.

A strategy for large-scale comparison of evolutionary- and reaction-based classifications of enzyme function.

Determining the molecular function of enzymes discovered by genome sequencing represents a primary foundation for understanding many aspects of biology. Historically, classification of enzyme reactions has used the enzyme nomenclature system developed to describe the overall reactions performed by biochemically characterized enzymes, irrespective of their associated sequences. In contrast, functional classification and assignment for the millions of protein sequences of unknown function now available is largely done in two computational steps, first by similarity-based assignment of newly obtained sequences to homologous groups, followed by transferring to them the known functions of similar biochemically characterized homologs. Due to the fundamental differences in their etiologies and practice, `how' these chemistry- and evolution-centric functional classification systems relate to each other has been difficult to explore on a large scale. To investigate this issue in a new way, we integrated two published ontologies that had previously described each of these classification systems independently. The resulting infrastructure was then used to compare the functional assignments obtained from each classification system for the well-studied and functionally diverse enolase superfamily. Mapping these function assignments to protein structure and reaction similarity networks shows a profound and complex disconnect between the homology- and chemistry-based classification systems. This conclusion mirrors previous observations suggesting that except for closely related sequences, facile annotation transfer from small numbers of characterized enzymes to the huge number uncharacterized homologs to which they are related is problematic. Our extension of these comparisons to large enzyme superfamilies in a computationally intelligent manner provides a foundation for new directions in protein function prediction for the huge proportion of sequences of unknown function represented in major databases. Interactive sequence, reaction, substrate and product similarity networks computed for this work for the enolase and two other superfamilies are freely available for download from the Structure Function Linkage Database Archive (http://sfld.rbvi.ucsf.edu).