Eukaryotic Proteomes Research Articles

The Functional Comparison of Eukaryotic Proteomes: Implications for Choosing an Appropriate Model Organism to Probe Human Biology.

Phenotypic differences between species are, in significant part, determined by their proteomic diversity. The link between proteomic and phenotypic diversity can be best understood in the context of the various pathways and biological processes in which proteins participate. While the conservation pattern for individual proteins across species is expected to follow the phylogenetic relationships among the species, the diversity patterns of individual pathways may not: certain pathways may be much more conserved among distantly related species than two closely related species, owing to the ecological histories of the species. Thus, a pathway-centric analysis of proteome conservation and diversity has important implications for the appropriate choice of a model organism when investigating specific aspects of human biology. Exploiting the complete genome sequences and protein-coding gene annotations, here we perform a comprehensive gene-set-centric analysis of proteomic diversity between humans and 54 eukaryotic organisms, resulting in a catalog of organisms that are most similar to humans in terms of specific pathways, processes, expression patterns, and diseases. We corroborate our findings using species-specific mass spectrometry data.Our analysis provides a general framework to identify conserved and unique pathways in a group of organisms and a resource to prioritize appropriate model systems to study a specific biological system in a reference organism such as humans.

Elaboration of the Homer1 recognition landscape reveals incomplete divergence of paralogous EVH1 domains.

Short sequences that mediate interactions with modular binding domains are ubiquitous throughout eukaryotic proteomes. Networks of short linear motifs (SLiMs) and their corresponding binding domains orchestrate many cellular processes, and the low mutational barrier to evolving novel interactions provides a way for biological systems to rapidly sample selectable phenotypes. Mapping SLiM binding specificity and the rules that govern SLiM evolution is fundamental to uncovering the pathways regulated by these networks and developing the tools to manipulate them. We used high-throughput screening of the human proteome to identify sequences that bind to the Enabled/VASP homology 1 (EVH1) domain of the postsynaptic density scaffolding protein Homer1. This expanded our understanding of the determinants of Homer EVH1 binding preferences and defined a new motif that can facilitate the discovery of additional Homer-mediated interactions. Interestingly, the Homer1 EVH1 domain preferentially binds to sequences containing an N-terminally overlapping motif that is bound by the paralogous family of Ena/VASP actin polymerases, and many of these sequences can bind to EVH1 domains from both protein families. We provide evidence from orthologous EVH1 domains in pre-metazoan organisms that the overlap in human Ena/VASP and Homer binding preferences corresponds to an incomplete divergence from a common Ena/VASP ancestor. Given this overlap in binding profiles, promiscuous sequences that can be recognized by both families either achieve specificity through extrinsic regulatory strategies or may provide functional benefits via multi-specificity. This may explain why these paralogs incompletely diverged despite the accessibility of further diverged isoforms.

Evolutionary analyses of intrinsically disordered regions reveal widespread signals of conservation.

Intrinsically disordered regions (IDRs) are segments of proteins without stable three-dimensional structures. As this flexibility allows them to interact with diverse binding partners, IDRs play key roles in cell signaling and gene expression. Despite the prevalence and importance of IDRs in eukaryotic proteomes and various biological processes, associating them with specific molecular functions remains a significant challenge due to their high rates of sequence evolution. However, by comparing the observed values of various IDR-associated properties against those generated under a simulated model of evolution, a recent study found most IDRs across the entire yeast proteome contain conserved features. Furthermore, it showed clusters of IDRs with common "evolutionary signatures," i.e. patterns of conserved features, were associated with specific biological functions. To determine if similar patterns of conservation are found in the IDRs of other systems, in this work we applied a series of phylogenetic models to over 7,500 orthologous IDRs identified in the Drosophila genome to dissect the forces driving their evolution. By comparing models of constrained and unconstrained continuous trait evolution using the Brownian motion and Ornstein-Uhlenbeck models, respectively, we identified signals of widespread constraint, indicating conservation of distributed features is mechanism of IDR evolution common to multiple biological systems. In contrast to the previous study in yeast, however, we observed limited evidence of IDR clusters with specific biological functions, which suggests a more complex relationship between evolutionary constraints and function in the IDRs of multicellular organisms.

RMATS-turbo: an efficient and flexible computational tool for alternative splicing analysis of large-scale RNA-seq data.

Pre-mRNA alternative splicing is a prevalent mechanism for diversifying eukaryotic transcriptomes and proteomes. Regulated alternative splicing plays a role in many biological processes, and dysregulated alternative splicing is a feature of many human diseases. Short-read RNA sequencing (RNA-seq) is now the standard approach for transcriptome-wide analysis of alternative splicing. Since 2011, our laboratory has developed and maintained Replicate Multivariate Analysis of Transcript Splicing (rMATS), a computational tool for discovering and quantifying alternative splicing events from RNA-seq data. Here we provide a protocol for the contemporary version of rMATS, rMATS-turbo, a fast and scalable re-implementation that maintains the statistical framework and user interface of the original rMATS software, while incorporating a revamped computational workflow with a substantial improvement in speed and data storage efficiency. The rMATS-turbo software scales up to massive RNA-seq datasets with tens of thousands of samples. To illustrate the utility of rMATS-turbo, we describe two representative application scenarios. First, we describe a broadly applicable two-group comparison to identify differential alternative splicing events between two sample groups, including both annotated and novel alternative splicing events. Second, we describe a quantitative analysis of alternative splicing in a large-scale RNA-seq dataset (~1,000 samples), including the discovery of alternative splicing events associated with distinct cell states. We detail the workflow and features of rMATS-turbo that enable efficient parallel processing and analysis of large-scale RNA-seq datasets on a compute cluster. We anticipate that this protocol will help the broad user base of rMATS-turbo make the best use of this software for studying alternative splicing in diverse biological systems.

An atlas of protein homo-oligomerization across domains of life

Protein structures are essential to understanding cellular processes in molecular detail. While advances in artificial intelligence revealed the tertiary structure of proteins at scale, their quaternary structure remains mostly unknown. We devise a scalable strategy based on AlphaFold2 to predict homo-oligomeric assemblies across four proteomes spanning the tree of life. Our results suggest that approximately 45% of an archaeal proteome and a bacterial proteome and 20% of two eukaryotic proteomes form homomers. Our predictions accurately capture protein homo-oligomerization, recapitulate megadalton complexes, and unveil hundreds of homo-oligomer types, including three confirmed experimentally by structure determination. Integrating these datasets with omics information suggests that a majority of known protein complexes are symmetric. Finally, these datasets provide a structural context for interpreting disease mutations and reveal coiled-coil regions as major enablers of quaternary structure evolution in human. Our strategy is applicable to any organism and provides a comprehensive view of homo-oligomerization in proteomes.

Tracking the effects of PLGA-based nanoparticles on protein expression in living cells through quantitative proteomics.

Mass spectrometry (MS)-based proteomics can identify and quantify the differential abundance of expressed proteins in parallel, and bottom-up proteomic approaches are even approaching comprehensive coverage of the complex eukaryotic proteome. Protein-nanoparticle (NP) interactions have been extensively studied owing to their importance in biological applications and nanotoxicology. However, the proteome-level effects of NPs on cells have received little attention, although changes in protein abundance can reflect the direct effects of nanocarriers on protein expression. Herein, we investigated the effect of PLGA-based NPs on protein expression in HepG2 cells using a label-free quantitative proteomics approach with data independent acquisition (DIA). The percentage of two-fold change in the protein expression of cells treated with PLGA-based NPs was less than 10.15% during a 6 hour observation period. Among the changed proteins, we found that dynamic proteins involved in cell division, localization, and transport are more likely to be more susceptible to PLGA-based NPs.

Mechanistic Insights into Alternative Gene Splicing in Oxidative Stress and Tissue Injury.

Significance: Oxidative stress (OS) and inflammation are inducers of tissue injury. Alternative splicing (AS) is an essential regulatory step for diversifying the eukaryotic proteome. Human diseases link AS to OS; however, the underlying mechanisms must be better understood. Recent Advances: Genome‑wide profiling studies identify new differentially expressed genes induced by OS-dependent ischemia/reperfusion injury. Overexpression of RNA-binding protein RBFOX1 protects against inflammation. Hypoxia-inducible factor-1α directs polypyrimidine tract binding protein 1 to regulate mouse carcinoembryonic antigen-related cell adhesion molecule 1 (Ceacam1) AS under OS conditions. Heterogeneous nuclear ribonucleoprotein L variant 1 contains an RGG/RG motif that coordinates with transcription factors to influence human CEACAM1 AS. Hypoxia intervention involving short interfering RNAs directed to long-noncoding RNA 260 polarizes M2 macrophages toward an anti-inflammatory phenotype and alleviates OS by inhibiting IL-28RA gene AS. Critical Issues: Protective mechanisms that eliminate reactive oxygen species (ROS) are important for resolving imbalances that lead to chronic inflammation. Defects in AS can cause ROS generation, cell death regulation, and the activation of innate and adaptive immune factors. We propose that AS pathways link redox regulation to the activation or suppression of the inflammatory response during cellular stress. Future Directions: Emergent studies using molecule-mediated RNA splicing are being conducted to exploit the immunogenicity of AS protein products. Deciphering the mechanisms that connect misspliced OS and pathologies should remain a priority. Controlled release of RNA directly into cells with clinical applications is needed as the demand for innovative nucleic acid delivery systems continues to be demonstrated.

Systematic identification of conditionally folded intrinsically disordered regions by AlphaFold2.

The AlphaFold Protein Structure Database contains predicted structures for millions of proteins. For the majority of human proteins that contain intrinsically disordered regions (IDRs), which do not adopt a stable structure, it is generally assumed that these regions have low AlphaFold2 confidence scores that reflect low-confidence structural predictions. Here, we show that AlphaFold2 assigns confident structures to nearly 15% of human IDRs. By comparison to experimental NMR data for a subset of IDRs that are known to conditionally fold (i.e., upon binding or under other specific conditions), we find that AlphaFold2 often predicts the structure of the conditionally folded state. Based on databases of IDRs that are known to conditionally fold, we estimate that AlphaFold2 can identify conditionally folding IDRs at a precision as high as 88% at a 10% false positive rate, which is remarkable considering that conditionally folded IDR structures were minimally represented in its training data. We find that human disease mutations are nearly fivefold enriched in conditionally folded IDRs over IDRs in general and that up to 80% of IDRs in prokaryotes are predicted to conditionally fold, compared to less than 20% of eukaryotic IDRs. These results indicate that a large majority of IDRs in the proteomes of human and other eukaryotes function in the absence of conditional folding, but the regions that do acquire folds are more sensitive to mutations. We emphasize that the AlphaFold2 predictions do not reveal functionally relevant structural plasticity within IDRs and cannot offer realistic ensemble representations of conditionally folded IDRs.

Probing excited state 1Hα chemical shifts in intrinsically disordered proteins with a triple resonance-based CEST experiment: Application to a disorder-to-order switch.

Over 40% of eukaryotic proteomes and 15% of bacterial proteomes are predicted to be intrinsically disordered based on their amino acid sequence. Intrinsically disordered proteins (IDPs) exist as heterogeneous ensembles of interconverting conformations and pose a challenge to the structure-function paradigm by apparently functioning without possessing stable structural elements. IDPs play a prominent role in biological processes involving extensive intermolecular interaction networks and their inherently dynamic nature facilitates their promiscuous interaction with multiple structurally diverse partner molecules. NMR spectroscopy has made pivotal contributions to our understanding of IDPs because of its unique ability to characterize heterogeneity at atomic resolution. NMR methods such as Chemical Exchange Saturation Transfer (CEST) and relaxation dispersion have enabled the detection of 'invisible' excited states in biomolecules which are transiently and sparsely populated, yet central for function. Here, we develop a 1Hα CEST pulse sequence which overcomes the resonance overlap problem in the 1Hα-13Cα plane of IDPs by taking advantage of the superior resolution in the 1H-15N correlation spectrum. In this sequence, magnetization is transferred after 1H CEST using a triple resonance coherence transfer pathway from 1Hα (i) to 1HN(i+1) during which the 15N(t1) and 1HN(t2) are frequency labelled. This approach is integrated with spin state-selective CEST for eliminating spurious dips in CEST profiles resulting from dipolar cross-relaxation. We apply this sequence to determine the excited state 1Hα chemical shifts of the intrinsically disordered DNA binding domain (CytRN) of the bacterial cytidine repressor (CytR), which transiently acquires a functional globally folded conformation. The structure of the excited state, calculated using 1Hα chemical shifts in conjunction with other excited state NMR restraints, is a three-helix bundle incorporating a helix-turn-helix motif that is vital for binding DNA.

Fluorosulfate as a Latent Sulfate in Peptides and Proteins.

Sulfation widely exists in the eukaryotic proteome. However, understanding the biological functions of sulfation in peptides and proteins has been hampered by the lack of methods to control its spatial or temporal distribution in the proteome. Herein, we report that fluorosulfate can serve as a latent precursor of sulfate in peptides and proteins, which can be efficiently converted to sulfate by hydroxamic acid reagents under physiologically relevant conditions. Photocaging the hydroxamic acid reagents further allowed for the light-controlled activation of functional sulfopeptides. This work provides a valuable tool for probing the functional roles of sulfation in peptides and proteins.

Strategies for rapid production of crystallization quality coatomer WD40 domains

The vesicular secretion of soluble cargo proteins from the endoplasmic reticulum (ER) is accompanied by the export of ER-resident membrane proteins that are co-packaged in secretory vesicles. The cytosolic coatomer protein complex I (COPI) utilizes the N-terminal WD40 domains of α-COPI and β′-COPI subunits to bind these membrane protein “clients” for ER retrieval. These “αWD40” and “β'WD40” domains are structural homologs that demonstrate distinct selectivity for client proteins. However, elucidation of the atomic-level principles of coatomer-client interactions has been challenging due to the tendency of αWD40 domain to undergo aggregation during expression and purification. Here we describe a rapid recombinant production strategy from E. coli, which substantially enhances the quality of the purified αWD40 domain. The αWD40 purification and crystallization are completed within one day, which minimizes aggregation losses and yields a 1.9 Å resolution crystal structure. We demonstrate the versatility of this strategy by applying it to purify the β′WD40 domain, which yields crystal structures in the 1.2–1.3 Å resolution range. As an alternate recombinant production system, we develop a cost-effective strategy for αWD40 production in human Expi293 cells. Finally, we suggest a roadmap to simplify these protocols further, which is of significance for the production of WD40 mutants prone to rapid aggregation. The WD40 production strategies presented here are likely to have broad applications because the WD40 domain represents one of the largest families of biomolecular interaction modules in the eukaryotic proteome and is critical for trafficking of host as well as viral proteins such as the SARS-CoV-2 spike protein.

Origin and Early Diversification of the Papain Family of Cysteine Peptidases.

Peptidases of the papain family play a key role in protein degradation, regulated proteolysis, and the host-pathogen arms race. Although the papain family has been the subject of many studies, knowledge about its diversity, origin, and evolution in Eukaryota, Bacteria, and Archaea is limited; thus, we aimed to address these long-standing knowledge gaps. We traced the origin and expansion of the papain family with a phylogenomic analysis, using sequence data from numerous prokaryotic and eukaryotic proteomes, transcriptomes, and genomes. We identified the full complement of the papain family in all prokaryotic and eukaryotic lineages. Analysis of the papain family provided strong evidence for its early diversification in the ancestor of eukaryotes. We found that the papain family has undergone complex and dynamic evolution through numerous gene duplications, which produced eight eukaryotic ancestral paralogous C1A lineages during eukaryogenesis. Different evolutionary forces operated on C1A peptidases, including gene duplication, horizontal gene transfer, and gene loss. This study challenges the current understanding of the origin and evolution of the papain family and provides valuable insights into their early diversification. The findings of this comprehensive study provide guidelines for future structural and functional studies of the papain family.

Chemoselective Late-Stage Functionalization of Peptides via Photocatalytic C2-Alkylation of Tryptophan.

Across eukaryotic proteomes, tryptophan is the least abundant of the 20 canonical amino acids, which makes it an ideal chemical handle for the late-stage functionalization of peptide and protein scaffolds with minimal production of undesired isoforms. Herein, we report the photocatalytic C2-alkylation of tryptophan using bromodifluoroacetate/acetamide-derived radical precursors. This rapid visible-light-mediated reaction is additive-free, operationally simple, and tolerates diverse functionality. We demonstrate the late-stage modification of a variety of complex peptides, including examples of biological significance.

On the Unknown Proteins of Eukaryotic Proteomes.

To study unknown proteins on a large scale, a reference system has been set up for the three better studied eukaryotic kingdoms, built with 36 proteomes as taxonomically diverse as possible. Proteins from 362 other eukaryotic proteomes with no known homologue in this set were then analyzed, focusing noteworthy on singletons, that is, on such proteins with no known homologue in their own proteome. Consistently, for a given species, no more than 12% of the singletons thus found are known at the protein level, according to Uniprot. In addition, since they rely on the information found in the alignment of homologous sequences, predictions of AlphaFold2 for their tridimensional structure are poor. In the case of metazoan species, the number of singletons rarely exceeds 1000 for the species the closest to the reference system (divergence times below 75 Myr). Interestingly, in the cases of viridiplantae and fungi, larger amounts of singletons are found for such species, as if the timescale on which singletons are added to proteomes were different in metazoa and in other eukaryotic kingdoms. In order to confirm this phenomenon, further studies of proteomes closer to those of the reference system are, however, needed.

Alternative Splicing: From Abiotic Stress Tolerance to Evolutionary Genomics

The post-transcriptional regulation of gene expression, in particular alternative splicing (AS) events, substantially contributes to the complexity of eukaryotic transcriptomes and proteomes [...].

Dynamics-Based Regulatory Switches of Type II Antitoxins: Insights into New Antimicrobial Discovery.

Type II toxin-antitoxin (TA) modules are prevalent in prokaryotes and are involved in cell maintenance and survival under harsh environmental conditions, including nutrient deficiency, antibiotic treatment, and human immune responses. Typically, the type II TA system consists of two protein components: a toxin that inhibits an essential cellular process and an antitoxin that neutralizes its toxicity. Antitoxins of type II TA modules typically contain the structured DNA-binding domain responsible for TA transcription repression and an intrinsically disordered region (IDR) at the C-terminus that directly binds to and neutralizes the toxin. Recently accumulated data have suggested that the antitoxin's IDRs exhibit variable degrees of preexisting helical conformations that stabilize upon binding to the corresponding toxin or operator DNA and function as a central hub in regulatory protein interaction networks of the type II TA system. However, the biological and pathogenic functions of the antitoxin's IDRs have not been well discussed compared with those of IDRs from the eukaryotic proteome. Here, we focus on the current state of knowledge about the versatile roles of IDRs of type II antitoxins in TA regulation and provide insights into the discovery of new antibiotic candidates that induce toxin activation/reactivation and cell death by modulating the regulatory dynamics or allostery of the antitoxin.

Fuzzy Drug Targets: Disordered Proteins in the Drug-Discovery Realm.

Intrinsically disordered proteins (IDPs) and regions (IDRs) form a large part of the eukaryotic proteome. Contrary to the structure-function paradigm, the disordered proteins perform a myriad of functions in vivo. Consequently, they are involved in various disease pathways and are plausible drug targets. Unlike folded proteins, that have a defined structure and well carved out drug-binding pockets that can guide lead molecule selection, the disordered proteins require alternative drug-development methodologies that are based on an acceptable picture of their conformational ensemble. In this review, we discuss various experimental and computational techniques that contribute toward understanding IDP "structure" and describe representative pursuances toward IDP-targeting drug development. We also discuss ideas on developing rational drug design protocols targeting IDPs.

Recombinant Expression and Chemical Amidation of Isotopically Labeled Native Melittin.

Post-translational modifications are ubiquitous in the eukaryotic proteome. However, these modifications are rarely incorporated in NMR studies of eukaryotic proteins, which are typically produced through recombinant expression in E. coli. Melittin is the primary peptide in honey bee venom. Its native C-terminal amide significantly affects its equilibrium structure and dynamics in solution and is thus a prerequisite for studying its native structure and function. Here, we present a method for producing triply isotopically labeled (2H, 13C, and 15N) native melittin through recombinant expression followed by chemical amidation. We then show that structural models produced with AlphaFold-Multimer are in even better agreement with experimental residual dipolar couplings than the 2.0 Å resolution X-ray crystal structure for residues G3-K23.

Proteome-Wide Detection and Annotation of Receptor Tyrosine Kinases (RTKs): RTK-PRED and the TyReK Database.

Receptor tyrosine kinases (RTKs) form a highly important group of protein receptors of the eukaryotic cell membrane. They control many vital cellular functions and are involved in the regulation of complex signaling networks. Mutations in RTKs have been associated with different types of cancers and other diseases. Although they are very important for proper cell function, they have been experimentally studied in a limited range of eukaryotic species. Currently, there is no available database for RTKs providing information about their function, expression, and interactions. Therefore, the identification of RTKs in multiple organisms, the documentation of their characteristics, and the collection of related information would be very useful. In this paper, we present a novel RTK detection pipeline (RTK-PRED) and the Receptor Tyrosine Kinases Database (TyReK-DB). RTK-PRED combines profile HMMs with transmembrane topology prediction to identify and classify potential RTKs. Proteins of all eukaryotic reference proteomes of the UniProt database were used as input in RTK-PRED leading to a filtered dataset of 20,478 RTKs. Based on the information collected for these RTKs from multiple databases, the relational TyReK database was created.

Lectin Affinity Chromatography.

Glycosylation is a non-template-driven posttranslational modification during which linked-sugars and glycans are added to the nascent polypeptide. Over 70% of the eukaryotic proteome is thought to be glycosylated. It is now known that correct glycosylation is essential for the correct folding, solubility, stability, and immunogenicity of proteins. In this chapter, we describe the technique of lectin affinity chromatography (LAC), a procedure that has the ability to distinguish different glycans, which are attached to proteins or lipids, termed glycoproteins or glycolipids, respectively. This method utilizes different immobilized lectins that have affinity for specific sugar substrates, to separate a wide range of glycan-attached complexes (Ambrosi et al., Org Biomol Chem 3:1593-1608, 2005). To further enhance the specificity of LAC, a corresponding free sugar may be used to produce a specific elution. In general, the conditions under which lectin affinity chromatography operates are relatively mild resulting in good biological recoveries of the glycoproteins.