Distinct Protein Classes in Human Red Cell Proteome Revealed by Similarity of Phylogenetic Profiles

The minimal set of proteins necessary to maintain a vertebrate cell forms an interesting core of cellular machinery. The known proteome of human red blood cell consists of about 1400 proteins. We treated this protein complement of one of the simplest human cells as a model and asked the questions on its function and origins. The proteome was mapped onto phylogenetic profiles, i.e. vectors of species possessing homologues of human proteins. A novel clustering approach was devised, utilising similarity in the phylogenetic spread of homologues as distance measure. The clustering based on phylogenetic profiles yielded several distinct protein classes differing in phylogenetic taxonomic spread, presumed evolutionary history and functional properties. Notably, small clusters of proteins common to vertebrates or Metazoa and other multicellular eukaryotes involve biological functions specific to multicellular organisms, such as apoptosis or cell-cell signaling, respectively. Also, a eukaryote-specific cluster is identified, featuring GTP-ase signalling and ubiquitination. Another cluster, made up of proteins found in most organisms, including bacteria and archaea, involves basic molecular functions such as oxidation-reduction and glycolysis. Approximately one third of erythrocyte proteins do not fall in any of the clusters, reflecting the complexity of protein evolution in comparison to our simple model. Basically, the clustering obtained divides the proteome into old and new parts, the former originating from bacterial ancestors, the latter from inventions within multicellular eukaryotes. Thus, the model human cell proteome appears to be made up of protein sets distinct in their history and biological roles. The current work shows that phylogenetic profiles concept allows protein clustering in a way relevant both to biological function and evolutionary history.     

Comparing function and structure between entire proteomes

More than 30 organisms have been sequenced entirely. Here, we applied a variety of simple bioinformatics tools to analyze 29 proteomes for representatives from all three kingdoms: eukaryotes, prokaryotes, and archaebacteria. We confirmed that eukaryotes have relatively more long proteins than prokaryotes and archaes, and that the overall amino acid composition is similar among the three. We predicted that approximately 15%-30% of all proteins contained transmembrane helices. We could not find a correlation between the content of membrane proteins and the complexity of the organism. In particular, we did not find significantly higher percentages of helical membrane proteins in eukaryotes than in prokaryotes or archae. However, we found more proteins with seven transmembrane helices in eukaryotes and more with six and 12 transmembrane helices in prokaryotes. We found twice as many coiled-coil proteins in eukaryotes (10%) as in prokaryotes and archaes (4%-5%), and we predicted approximately 15%-25% of all proteins to be secreted by most eukaryotes and prokaryotes. Every tenth protein had no known homolog in current databases, and 30%-40% of the proteins fell into structural families with >100 members. A classification by cellular function verified that eukaryotes have a higher proportion of proteins for communication with the environment. Finally, we found at least one homolog of experimentally known structure for approximately 20%-45% of all proteins; the regions with structural homology covered 20%-30% of all residues. These numbers may or may not suggest that there are 1200-2600 folds in the universe of protein structures. All predictions are available at http://cubic.bioc.columbia.edu/genomes.     

Single eubacterial origin of eukaryotic sulfide:quinone oxidoreductase,a mitochondrial enzyme conserved from the early evolution of eukaryotes during anoxic and sulfidic times

Mitochondria occur as aerobic, facultatively anaerobic, and, in the case of hydrogenosomes, strictly anaerobic forms. This physiological diversity of mitochondrial oxygen requirement is paralleled by that of free-living alpha-proteobacteria, the group of eubacteria from which mitochondria arose, many of which are facultative anaerobes. Although ATP synthesis in mitochondria usually involves the oxidation of reduced carbon compounds, many alpha-proteobacteria and some mitochondria are known to use sulfide (H2S) as an electron donor for the respiratory chain and its associated ATP synthesis. In many eubacteria, the oxidation of sulfide involves the enzyme sulfide:quinone oxidoreductase (SQR). Nuclear-encoded homologs of SQR are found in several eukaryotic genomes. Here we show that eukaryotic SQR genes characterized to date can be traced to a single acquisition from a eubacterial donor in the common ancestor of animals and fungi. Yet, SQR is not a well-conserved protein, and our analyses suggest that the SQR gene has furthermore undergone some lateral transfer among prokaryotes during evolution, leaving the precise eubacterial lineage from which eukaryotes obtained their SQR difficult to discern with phylogenetic methods. Newer geochemical data and microfossil evidence indicate that major phases of early eukaryotic diversification occurred during a period of the Earth's history from 1 to 2 billion years before present in which the subsurface ocean waters contained almost no oxygen but contained high concentrations of sulfide, suggesting that the ability to deal with sulfide was essential for prokaryotes and eukaryotes during that time. Notwithstanding poor resolution in deep SQR phylogeny and lack of a specifically alpha-protebacterial branch for the eukaryotic enzyme on the basis of current lineage sampling, a single eubacterial origin of eukaryotic SQR and the evident need of ancient eukaryotes to deal with sulfide, a process today germane to mitochondrial quinone reduction, are compatible with the view that eukaryotic SQR was an acquisition from the mitochondrial endosymbiont.     

The evolution of the genetic code took place in an anaerobic environment

We have compared orthologous proteins from an aerobic organism, Cytophaga hutchinsonii, and from an obligate anaerobe, Bacteroides thetaiotaomicron. This comparison allows us to define the oxyphobic ranks of amino acids, i.e. a scale of the relative sensitivity to oxygen of the amino acid residues. The oxyphobic index (OI), which can be simply obtained from the amino acids' oxyphobic ranks, can be associated to any protein and therefore to the genetic code, if the number of synonymous codons attributed to the amino acids in the code is assumed to be the frequency with which the amino acids appeared in ancestral proteins. Sampling of the OI variable from the proteins of obligate anaerobes and aerobes has established that the OI value of the genetic code is not significantly different from the mean OI value of anaerobe proteins, while it is different from that of aerobe proteins. This observation would seem to suggest that the terminal phases of the evolution of genetic code organization took place in an anaerobic environment. This result is discussed in the framework of hypotheses suggested to explain the timing of the evolutionary appearance of the aerobic metabolism.     

Variation among species in proteomic sulphur content is related to environmental conditions

The elemental composition of proteins influences the quantities of different elements required by organisms. Here, we considered variation in the sulphur content of whole proteomes among 19 Archaea, 122 Eubacteria and 10 eukaryotes whose genomes have been fully sequenced. We found that different species vary greatly in the sulphur content of their proteins, and that average sulphur content of proteomes and genome base composition are related. Forces contributing to variation in proteomic sulphur content appear to operate quite uniformly across the proteins of different species. In particular, the sulphur content of orthologous proteins was frequently correlated with mean proteomic sulphur contents. Among prokaryotes, proteomic sulphur content tended to be greater in anaerobes, relative to non-anaerobes. Thermophiles tended to have lower proteomic sulphur content than non-thermophiles, consistent with the thermolability of cysteine and methionine residues. This work suggests that persistent environmental growth conditions can influence the evolution of elemental composition of whole proteomes in a manner that may have important implications for the amount of sulphur used by living organisms to build proteins. It extends previous studies that demonstrated links between transient changes in environmental conditions and the elemental composition of subsets of proteins expressed under these conditions.     

The mechanism of superoxide scavenging by Archaeoglobus fulgidus neelaredoxin

Neelaredoxin is a mononuclear iron protein widespread among prokaryotic anaerobes and facultative aerobes, including human pathogens. It has superoxide scavenging activity, but the exact mechanism by which this process occurs has been controversial. In this report, we present the study of the reaction of superoxide with the reduced form of neelaredoxin from the hyperthermophilic archaeon Archaeoglobus fulgidus by pulse radiolysis. This protein reduces superoxide very efficiently (k = 1.5 x 10(9) m(-1)s(-1)), and the dismutation activity is rate-limited, in steady-state conditions, by the much slower superoxide oxidation step. These data show unambiguously that the superfamily of neelaredoxin-like proteins (including desulfoferrodoxin) presents a novel type of reactivity toward superoxide, a result of particular relevance for the understanding of both oxygen stress response mechanisms and, in particular, how pathogens may respond to the oxidative burst produced by the defense cells in eukaryotes. The actual in vivo functioning of these enzymes will depend strongly on the cell redox status. Further insight on the catalytic mechanism was obtained by the detection of a transient intermediate ferric species upon oxidation of neelaredoxin by superoxide, detectable by visible spectroscopy with an absorption maximum at 610 nm, blue-shifted approximately 50 nm from the absorption of the resting ferric state. The role of the iron sixth ligand, glutamate-12, in the reactivity of neelaredoxin toward superoxide was assessed by studying two site-directed mutants: E12Q and E12V.     

Orderly order in protein intrinsic disorder distribution: disorder in 3500 proteomes from viruses and the three domains of life

Intrinsically disordered proteins and intrinsically disordered protein regions are highly abundant in nature. However, the quantitative and qualitative measures of protein intrinsic disorder in species with known genomes are still not available. Furthermore, although the correlation between high fraction of disordered residues and advanced species has been reported, the details of this correlation and the connection between the disorder content and proteome complexity have not been reported as of yet. To fill this gap, we analysed entire proteomes of 3484 species from three domains of life (archaea, bacteria and eukaryotes) and from viruses. Our analysis revealed that the evolution process is characterized by distinctive patterns of changes in the protein intrinsic disorder content. We are showing here that viruses are characterized by the widest spread of the proteome disorder content (the percentage of disordered residues ranges from 7.3% in human coronavirus NL63 to 77.3% in Avian carcinoma virus). For several organisms, a clear correlation is seen between their disorder contents and habitats. In multicellular eukaryotes, there is a weak correlation between the complexity of an organism (evaluated as a number of different cell types) and its overall disorder content. For both the prokaryotes and eukaryotes, the disorder content is generally independent of the proteome size. However, disorder shows a sharp increase associated with the transition from prokaryotic to eukaryotic cells. This suggests that the increased disorder content in eukaryotic proteomes might be used by nature to deal with the increased cell complexity due to the appearance of the various cellular compartments.     

Acidophilic character of yeast PID261/BUD32, a putative ancestor of eukaryotic protein kinases

Yeast piD261/Bud32 and its homologues are present in eukaryotes and in archaea but not in bacteria and are believed to make up a primordial branch of the eukaryotic protein kinase superfamily. Here, we show that, at variance with the majority of Ser/Thr protein kinases which recognize phosphoacceptor sites specified by basic and/or proline residues, piD261 phosphorylates in vitro a number of acidic proteins and peptides, and it recognizes seryl residues specified by carboxylic side chains. These data suggest that recognition of acidic sites might have been a primordial trait of protein kinases, which was modified during evolution to cope with the increasing complexity of protein phosphorylation in eukaryotes.     

Direct Detection of Alternative Open Reading Frames Translation Products in Human Significantly Expands the Proteome

A fully mature mRNA is usually associated to a reference open reading frame encoding a single protein. Yet, mature mRNAs contain unconventional alternative open reading frames (AltORFs) located in untranslated regions (UTRs) or overlapping the reference ORFs (RefORFs) in non-canonical +2 and +3 reading frames. Although recent ribosome profiling and footprinting approaches have suggested the significant use of unconventional translation initiation sites in mammals, direct evidence of large-scale alternative protein expression at the proteome level is still lacking. To determine the contribution of alternative proteins to the human proteome, we generated a database of predicted human AltORFs revealing a new proteome mainly composed of small proteins with a median length of 57 amino acids, compared to 344 amino acids for the reference proteome. We experimentally detected a total of 1,259 alternative proteins by mass spectrometry analyses of human cell lines, tissues and fluids. In plasma and serum, alternative proteins represent up to 55% of the proteome and may be a potential unsuspected new source for biomarkers. We observed constitutive co-expression of RefORFs and AltORFs from endogenous genes and from transfected cDNAs, including tumor suppressor p53, and provide evidence that out-of-frame clones representing AltORFs are mistakenly rejected as false positive in cDNAs screening assays. Functional importance of alternative proteins is strongly supported by significant evolutionary conservation in vertebrates, invertebrates, and yeast. Our results imply that coding of multiple proteins in a single gene by the use of AltORFs may be a common feature in eukaryotes, and confirm that translation of unconventional ORFs generates an as yet unexplored proteome.     

Energy metabolism among eukaryotic anaerobes in light of Proterozoic ocean chemistry

Recent years have witnessed major upheavals in views about early eukaryotic evolution. One very significant finding was that mitochondria, including hydrogenosomes and the newly discovered mitosomes, are just as ubiquitous and defining among eukaryotes as the nucleus itself. A second important advance concerns the readjustment, still in progress, about phylogenetic relationships among eukaryotic groups and the roughly six new eukaryotic supergroups that are currently at the focus of much attention. From the standpoint of energy metabolism (the biochemical means through which eukaryotes gain their ATP, thereby enabling any and all evolution of other traits), understanding of mitochondria among eukaryotic anaerobes has improved. The mainstream formulations of endosymbiotic theory did not predict the ubiquity of mitochondria among anaerobic eukaryotes, while an alternative hypothesis that specifically addressed the evolutionary origin of energy metabolism among eukaryotic anaerobes did. Those developments in biology have been paralleled by a similar upheaval in the Earth sciences regarding views about the prevalence of oxygen in the oceans during the Proterozoic (the time from ca 2.5 to 0.6 Ga ago). The new model of Proterozoic ocean chemistry indicates that the oceans were anoxic and sulphidic during most of the Proterozoic. Its proponents suggest the underlying geochemical mechanism to entail the weathering of continental sulphides by atmospheric oxygen to sulphate, which was carried into the oceans as sulphate, fueling marine sulphate reducers (anaerobic, hydrogen sulphide-producing prokaryotes) on a global scale. Taken together, these two mutually compatible developments in biology and geology underscore the evolutionary significance of oxygen-independent ATP-generating pathways in mitochondria, including those of various metazoan groups, as a watermark of the environments within which eukaryotes arose and diversified into their major lineages.     

A PQL (protein quantity loci) analysis of mature pea seed proteins identifies loci determining seed protein composition

Legume seeds are a major source of dietary proteins for humans and animals. Deciphering the genetic control of their accumulation is thus of primary significance towards their improvement. At first, we analysed the genetic variability of the pea seed proteome of three genotypes over 3 years of cultivation. This revealed that seed protein composition variability was under predominant genetic control, with as much as 60% of the spots varying quantitatively among the three genotypes. Then, by combining proteomic and quantitative trait loci (QTL) mapping approaches, we uncovered the genetic architecture of seed proteome variability. Protein quantity loci (PQL) were searched for 525 spots detected on 2-D gels obtained for 157 recombinant inbred lines. Most protein quantity loci mapped in clusters, suggesting that the accumulation of the major storage protein families was under the control of a limited number of loci. While convicilin accumulation was mainly under the control of cis-regulatory regions, vicilins and legumins were controlled by both cis- and trans-regulatory regions. Some loci controlled both seed protein composition and protein content and a locus on LGIIa appears to be a major regulator of protein composition and of protein in vitro digestibility.     

Characterization of bud emergence 46 (BEM46) protein: Sequence,structural, phylogenetic and subcellular localization analyses

The bud emergence 46 (BEM46) protein from Neurospora crassa belongs to the α/β-hydrolase superfamily. Recently, we have reported that the BEM46 protein is localized in the perinuclear ER and also forms spots close by the plasma membrane. The protein appears to be required for cell type-specific polarity formation in N. crassa. Furthermore, initial studies suggested that the BEM46 amino acid sequence is conserved in eukaryotes and is considered to be one of the widespread conserved “known unknown” eukaryotic genes. This warrants for a comprehensive phylogenetic analysis of this superfamily to unravel origin and molecular evolution of these genes in different eukaryotes. Herein, we observe that all eukaryotes have at least a single copy of a bem46 ortholog. Upon scanning of these proteins in various genomes, we find that there are expansions leading into several paralogs in vertebrates. Usingcomparative genomic analyses, we identified insertion/deletions (indels) in the conserved domain of BEM46 protein, which allow to differentiate fungal classes such as ascomycetes from basidiomycetes. We also find that exonic indels are able to differentiate BEM46 homologs of different eukaryotic lineage. Furthermore, we unravel that BEM46 protein from N. crassa possess a novel endoplasmic-retention signal (PEKK) using GFP-fusion tagging experiments. We propose that three residues namely a serine 188S, a histidine 292H and an aspartic acid 262D are most critical residues, forming a catalytic triad in BEM46 protein from N. crassa. We carried out a comprehensive study on bem46 genes from a molecular evolution perspective with combination of functional analyses. The evolutionary history of BEM46 proteins is characterized by exonic indels in lineage specific manner.     

Families of membranous proteins can be characterized by the amino acid composition of their transmembrane domains

SUMMARY: In eukaryotes, membranous proteins account for 20-30% of the proteome. Most of these proteins contain one or more transmembrane (TM) domains. These are short segments that transverse the bilayer lipid membrane. Various properties of the TM domains, such as their number, their topology and their arrangement within the membrane, are closely related to the protein's cellular functions. The properties of the TM domains also determine the cellular targeting and localization of these proteins. It is not known, however, whether the information encoded by TM domains suffices for the purpose of classifying proteins into their functional families. This is the question we address here. We introduce an algorithm that creates a profile of each functional family of membranous proteins based only on the amino acid composition of their TM domains. This is complemented by a classifier program for each such family (to determine whether a given protein belongs to it or not). We find that in most instances TM domains contain enough information to allow an accurate discrimination of approximately 80% sensitivity and approximately 90% specificity among unrelated polytopic functional families with the same number of TM domains. SUPPLEMENTARY INFORMATION: Available at www.protonet.cs.huji.ac.il/TM/     

A Computational Algorithm for Functional Clustering of Proteome Dynamics During Development

Phenotypic traits, such as seed development, are a consequence of complex biochemical interactions among genes, proteins and metabolites, but the underlying mechanisms that operate in a coordinated and sequential manner remain elusive. Here, we address this issue by developing a computational algorithm to monitor proteome changes during the course of trait development. The algorithm is built within the mixture-model framework in which each mixture component is modeled by a specific group of proteins that display a similar temporal pattern of expression in trait development. A nonparametric approach based on Legendre orthogonal polynomials was used to fit dynamic changes of protein expression, increasing the power and flexibility of protein clustering. By analyzing a dataset of proteomic dynamics during early embryogenesis of the Chinese fir, the algorithm has successfully identified several distinct types of proteins that coordinate with each other to determine seed development in this forest tree commercially and environmentally important to China. The algorithm will find its immediate applications for the characterization of mechanistic underpinnings for any other biological processes in which protein abundance plays a key role.     

What is the total number of protein molecules per cell volume? A call to rethink some published values

Novel methods such as mass‐spectrometry enable a view of the proteomes of cells in unprecedented detail. Recently, these efforts have culminated in quantitative measurements of the number of copies per cell for most expressed proteins in organisms ranging from bacteria to mammalian cells. Here, we estimate the expected total number of proteins per unit of cell volume using known parameters related to the composition of cells such as the fraction of cell mass that is protein, and the average protein length. Using simple arguments, we estimate a range of 2–4 million proteins per cubic micron (i.e. 1 fL) in bacteria, yeast, and mammalian cells. Interestingly, we find that measured values that are reported for fission yeast and mammalian cells are often about 3–10 times lower. We discuss this apparent discrepancy and how to use the estimate as benchmark to recalibrate proteome‐wide quantitative censuses or to revisit assumptions about cell composition.     

Lysine biotinylation and methionine oxidation in the heat shock protein HSP60 synergize in the elimination of reactive oxygen species in human cell cultures

Previous studies suggest that the number of proteins containing covalently bound biotin is larger than previously thought. Here, we report the identity of some of these proteins. Using mass spectrometry, we discovered 108 novel biotinylation sites in the human embryonic kidney HEK293 cell proteome; members of the heat shock protein (HSP) superfamily were overrepresented among the novel biotinylated proteins. About half of the biotinylated proteins also displayed various degrees of methionine oxidation, which is known to play an important role in the defense against reactive oxygen species; for biotinylated HSPs, the percent of methionine sulfoxidation approached 100%. Protein structure analysis suggests that methionine sulfoxides localize in close physical proximity to the biotinylated lysines on the protein surface. Mass spectrometric analysis revealed that between one and five of the methionine residues in the C-terminal KEEKDPGMGAMGGMGGGMGGGMF motif are oxidized in HSP60. The likelihood of methionine sulfoxidation is higher if one of the adjacent lysine residues is biotinylated. Knockdown of HSP60 caused a 60% increase in the level of reactive oxygen species in fibroblasts cultured in biotin-sufficient medium. When HEK293 cells were transferred from biotin-sufficient medium to biotin-free medium, the level of reactive oxygen species increased by >9 times compared with baseline controls and a time-response relationship was evident. High levels of methionine sulfoxidation coincided with cell cycle arrest in the G0/G1 and S phases in biotin-depleted cells. We conclude that biotinylation of lysines synergizes with sulfoxidation of methionines in heat shock proteins such as HSP60 in the defense against reactive oxygen species.     

Predicting rates of interspecific interaction from phylogenetic trees

Integrating phylogenetic information can potentially improve our ability to explain species' traits, patterns of community assembly, the network structure of communities, and ecosystem function. In this study, we use mathematical models to explore the ecological and evolutionary factors that modulate the explanatory power of phylogenetic information for communities of species that interact within a single trophic level. We find that phylogenetic relationships among species can influence trait evolution and rates of interaction among species, but only under particular models of species interaction. For example, when interactions within communities are mediated by a mechanism of phenotype matching, phylogenetic trees make specific predictions about trait evolution and rates of interaction. In contrast, if interactions within a community depend on a mechanism of phenotype differences, phylogenetic information has little, if any, predictive power for trait evolution and interaction rate. Together, these results make clear and testable predictions for when and how evolutionary history is expected to influence contemporary rates of species interaction.