Detailed proteome analysis of growing cells of the planctomycete Rhodopirellula baltica SH1T

Rhodopirellula baltica SH1(T), which was isolated from the water column of the Kieler Bight, a bay in the southwestern Baltic Sea, is a marine aerobic, heterotrophic representative of the ubiquitous bacterial phylum Planctomycetes. We analyzed the R. baltica proteome by applying different preanalytical protein as well as peptide separation techniques (1-D and 2-DE, HPLC separation) prior to MS. That way, we could identify a total of 1115 nonredundant proteins from the intracellular proteome and from different cell wall protein fractions. With the contribution of 709 novel proteins resulting from this study, the current comprehensive R. baltica proteomic dataset consists of 1267 unique proteins (accounting for 17.3% of the total putative protein-coding ORFs), including 261 proteins with a predicted signal peptide. The identified proteins were functionally categorized using Clusters of Orthologous Groups (COGs), and their potential cellular locations were predicted by bioinformatic tools. A unique protein family that contains several YTV domains and is rich in cysteine and proline was found to be a component of the R. baltica proteinaceous cell wall. Based on this comprehensive proteome analysis a global schema of the major metabolic pathways of growing R. baltica cells was deduced.     

Hypoxia Is a Dominant Remodeler of the Effector T Cell Surface Proteome Relative to Activation and Regulatory T Cell Suppression

Immunosuppressive factors in the tumor microenvironment (TME) impair T cell function and limit the antitumor immune response. T cell surface receptors and surface proteins that influence interactions and function in the TME are proven targets for cancer immunotherapy. However, how the entire surface proteome remodels in primary human T cells in response to specific suppressive factors in the TME remains to be broadly and systematically characterized. Here, using a reductionist cell culture approach with primary human T cells and stable isotopic labeling with amino acids in cell culture–based quantitative cell surface capture glycoproteomics, we examined how two immunosuppressive TME factors, regulatory T cells (Tregs) and hypoxia, globally affect the activated CD8⁺ surface proteome (surfaceome). Surprisingly, coculturing primary CD8⁺ T cells with Tregs only modestly affected the CD8⁺ surfaceome but did partially reverse activation-induced surfaceomic changes. In contrast, hypoxia drastically altered the CD8⁺ surfaceome in a manner consistent with both metabolic reprogramming and induction of an immunosuppressed state. The CD4⁺ T cell surfaceome similarly responded to hypoxia, revealing a common hypoxia-induced surface receptor program. Our surfaceomics findings suggest that hypoxic environments create a challenge for T cell activation. These studies provide global insight into how Tregs and hypoxia remodel the T cell surfaceome and we believe represent a valuable resource to inform future therapeutic efforts to enhance T cell function.     

Hypoxia Is a Dominant Remodeler of the Effector T Cell Surface Proteome Relative to Activation and Regulatory T Cell Suppression

Immunosuppressive factors in the tumor microenvironment (TME) impair T cell function and limit the antitumor immune response. T cell surface receptors and surface proteins that influence interactions and function in the TME are proven targets for cancer immunotherapy. However, how the entire surface proteome remodels in primary human T cells in response to specific suppressive factors in the TME remains to be broadly and systematically characterized. Here, using a reductionist cell culture approach with primary human T cells and stable isotopic labeling with amino acids in cell culture–based quantitative cell surface capture glycoproteomics, we examined how two immunosuppressive TME factors, regulatory T cells (Tregs) and hypoxia, globally affect the activated CD8⁺ surface proteome (surfaceome). Surprisingly, coculturing primary CD8⁺ T cells with Tregs only modestly affected the CD8⁺ surfaceome but did partially reverse activation-induced surfaceomic changes. In contrast, hypoxia drastically altered the CD8⁺ surfaceome in a manner consistent with both metabolic reprogramming and induction of an immunosuppressed state. The CD4⁺ T cell surfaceome similarly responded to hypoxia, revealing a common hypoxia-induced surface receptor program. Our surfaceomics findings suggest that hypoxic environments create a challenge for T cell activation. These studies provide global insight into how Tregs and hypoxia remodel the T cell surfaceome and we believe represent a valuable resource to inform future therapeutic efforts to enhance T cell function.     

Towards the proteome of the marine bacterium Rhodopirellula baltica: mapping the soluble proteins

The marine bacterium Rhodopirellula baltica, a member of the phylum Planctomycetes, has distinct morphological properties and contributes to remineralization of biomass in the natural environment. On the basis of its recently determined complete genome we investigated its proteome by 2-DE and established a reference 2-DE gel for the soluble protein fraction. Approximately 1000 protein spots were excised from a colloidal Coomassie-stained gel (pH 4-7), analyzed by MALDI-MS and identified by PMF. The non-redundant data set contained 626 distinct protein spots, corresponding to 558 different genes. The identified proteins were classified into role categories according to their predicted functions. The experimentally determined and the theoretically predicted proteomes were compared. Proteins, which were most abundant in 2-DE gels and the coding genes of which were also predicted to be highly expressed, could be linked mainly to housekeeping functions in glycolysis, tricarboxic acid cycle, amino acid biosynthesis, protein quality control and translation. Absence of predictable signal peptides indicated a localization of these proteins in the intracellular compartment, the pirellulosome. Among the identified proteins, 146 contained a predicted signal peptide suggesting their translocation. Some proteins were detected in more than one spot on the gel, indicating post-translational modification. In addition to identifying proteins present in the published sequence database for R. baltica, an alternative approach was used, in which the mass spectrometric data was searched against a maximal ORF set, allowing the identification of four previously unpredicted ORFs. The 2-DE reference map presented here will serve as framework for further experiments to study differential gene expression of R. baltica in response to external stimuli or cellular development and compartmentalization.     

Proteomic analysis of carbohydrate catabolism and regulation in the marine bacterium Rhodopirellula baltica

The marine bacterium Rhodopirellula baltica is a model organism for aerobic carbohydrate degradation in marine systems, where polysaccharides represent the dominant components of biomass. On the basis of the genome sequence and a 2-D map of soluble proteins, the central catabolic routes of R. baltica were reconstructed. Almost all enzymes of glycolysis and TCA cycle were identified. In addition, almost all enzymes of the oxidative branch of the pentose phosphate cycle were detected. This proteomic reconstruction was corroborated by determination of selected enzymatic activities. To study substrate-dependent regulation in R. baltica, cells were adapted to growth with eight different carbohydrates and profiled with 2-DE for changes in protein patterns. Relative abundances of regulated proteins were determined using the 2-D DIGE technology and protein identification was achieved by PMF. Most of the up-regulated proteins were either dehydrogenases/oxidoreductases or proteins of unknown function which are unique for R. baltica. For only some of the regulated proteins, the coding genes are located in a physiologically meaningful genomic context. e.g., a ribose-induced alcohol dehydrogenase is encoded within an operon-like structure together with genes coding for a ribose-specific ABC-transporter. However, most of the regulated genes are randomly distributed across the genome.     

Novel protein domains and motifs in the marine planctomycete Rhodopirellula baltica

The planctomycetes are a phylum of bacteria that have a unique cell compartmentalisation and yeast-like budding cell division and peptidoglycan-less proteinaceous cell walls. We wished to further our understanding of these unique organisms at the molecular level by searching for conserved amino acid sequence motifs and domains in the proteins encoded by Rhodopirellula baltica. Using BLAST and single-linkage clustering, we have discovered several new protein domains and sequence motifs in this planctomycete. R. baltica has multiple members of the newly discovered GEFGR protein family and the ASPIC C-terminal domain family, whilst most other organisms for which whole genome sequence is available have no more than one. Many of the domains and motifs appear to be restricted to the planctomycetes. It is possible that these protein domains and motifs may have been lost or replaced in other phyla, or they may have undergone multiple duplication events in the planctomycete lineage. One of the novel motifs probably represents a novel N-terminal export signal peptide. With their unique cell biology, it may be that the planctomycete cell compartmentalisation plan in particular needs special membrane transport mechanisms. The discovery of these new domains and motifs, many of which are associated with secretion and cell-surface functions, will help to stimulate experimental work and thus enhance further understanding of this fascinating group of organisms.     

Growth phase dependent regulation of protein composition in Rhodopirellula baltica

Growth phase dependent changes of protein composition in the marine bacterium Rhodopirellula baltica were quantitatively monitored by applying the two-dimensional difference gel electrophoresis (2D DIGE) technology. The number of regulated proteins (fold changes in protein abundance > absolute value(2)) increased from early (10) to late stationary growth phase (179), with fold changes reaching maximal values of 40. About 110 of these regulated protein spots were analysed by MALDI-TOF-MS and identified by mapping of peptide masses. Results indicate an opposing regulation of tricarboxylic acid cycle and oxidative pentose phosphate cycle, a downregulation of several enzymes involved in amino acid biosynthesis and an upregulation of the alternative sigma factor sigmaH in stationary phase. Interestingly, 26 proteins of unknown function were up- or downregulated in the stationary phase. Several proteins were specifically regulated during growth on solid surface (agar plates). These proteins could possibly be involved in the development of the different R. baltica morphotypes, i.e. motile swarmer cells and sessile cell aggregates (so-called rosettes).     

Malfunctioning of adipocytes in obesity is linked to quantitative surfaceome changes

Increased triglyceride accumulation in adipocytes caused by a misbalance between energy intake and energy consumption, results in increased adipocyte size, excess adipose tissue, increased body weight and ultimately, obesity. It is well established that enlarged adipocytes exhibit malfunctions that contribute to whole body insulin resistance, a key factor for the development of type 2 diabetes. However, the underlying molecular cause for dysfunctional adipocyte behavior and signaling is poorly understood. Since the adipocyte cell surface proteome, or surfaceome, represents the cellular signaling gateway to the microenvironment, we studied the contribution of this subproteome to adipocyte malfunctions in obesity. By using the chemoproteomic Cell Surface Capture (CSC) technology, we established surfaceome maps of primary adipocytes derived from different mouse models for metabolic disorders. Relative quantitative comparison between these surfaceome maps revealed a set of cell surface glycoproteins with modulated location-specific abundance levels. RNAi mediated targeting of a subset of the detected obesity modulated cell surface glycoproteins in an in vitro model system provided functional evidence for their role in adiponectin secretion and the lipolytic activity of adipocytes. Thus, we conclude that the identified cell surface glycoproteins which exhibit obesity induced abundance changes and impact adipocyte function at the same time contribute to adipocyte malfunction in obesity. The regulation of their concerted activities could improve adipocyte function in obesity.     

Protein cell surface display in Gram-positive bacteria: from single protein to macromolecular protein structure

In the course of evolution, Gram-positive bacteria, defined here as prokaryotes from the domain Bacteria with a cell envelope composed of one biological membrane (monodermita) and a cell wall composed at least of peptidoglycan and covalently linked teichoic acids, have developed several mechanisms permitting to a cytoplasmic synthesized protein to be present on the bacterial cell surface. Four major types of cell surface displayed proteins are currently recognized: (i) transmembrane proteins, (ii) lipoproteins, (iii) LPXTG-like proteins and (iv) cell wall binding proteins. The subset of proteins exposed on the bacterial cell surface, and thus interacting with extracellular milieu, constitutes the surfaceome. Here, we review exhaustively the current molecular mechanisms involved in protein attachment within the cell envelope of Gram-positive bacteria, from single protein to macromolecular protein structure.     

A Mass Spectrometric-Derived Cell Surface Protein Atlas

Cell surface proteins are major targets of biomedical research due to their utility as cellular markers and their extracellular accessibility for pharmacological intervention. However, information about the cell surface protein repertoire (the surfaceome) of individual cells is only sparsely available. Here, we applied the Cell Surface Capture (CSC) technology to 41 human and 31 mouse cell types to generate a mass-spectrometry derived Cell Surface Protein Atlas (CSPA) providing cellular surfaceome snapshots at high resolution. The CSPA is presented in form of an easy-to-navigate interactive database, a downloadable data matrix and with tools for targeted surfaceome rediscovery (http://wlab.ethz.ch/cspa). The cellular surfaceome snapshots of different cell types, including cancer cells, resulted in a combined dataset of 1492 human and 1296 mouse cell surface glycoproteins, providing experimental evidence for their cell surface expression on different cell types, including 136 G-protein coupled receptors and 75 membrane receptor tyrosine-protein kinases. Integrated analysis of the CSPA reveals that the concerted biological function of individual cell types is mainly guided by quantitative rather than qualitative surfaceome differences. The CSPA will be useful for the evaluation of drug targets, for the improved classification of cell types and for a better understanding of the surfaceome and its concerted biological functions in complex signaling microenvironments.     

Surfaceome of Exosomes Secreted from the Colorectal Cancer Cell Line SW480: Peripheral and Integral Membrane Proteins Analyzed by Proteolysis and TX114

Exosomes are important bidirectional cell–cell communicators in normal and pathological physiology. Although exosomal surface membrane proteins (surfaceome) enable target cell recognition and are an attractive source of disease marker, they are poorly understood. Here, a comprehensive surfaceome analysis of exosomes secreted by the colorectal cancer cell line SW480 is described. Sodium carbonate extraction/Triton X‐114 phase separation and mild proteolysis (proteinase K, PK) of intact exosomes is used in combination with label‐free quantitative mass spectrometry to identify 1025 exosomal proteins of which 208 are predicted to be integral membrane proteins (IMPs) according to TOPCONS and GRAVY scores. Interrogation of UniProt database‐annotated proteins reveals 124 predicted peripherally‐associated membrane proteins (PMPs). Surprisingly, 108 RNA‐binding proteins (RBPs)/RNA nucleoproteins (RNPs) are found in the carbonate/Triton X‐114 insoluble fraction. Mild PK treatment of SW480‐GFP labeled exosomes reveal 58 proteolytically cleaved IMPs and 14 exoplasmic PMPs (e.g., CLU/GANAB/LGALS3BP). Interestingly, 18 RBPs/RNPs (e.g., EIF3L/RPL6) appear bound to the outer exosome surface since they are sensitive to PK proteolysis. The finding that outer surface‐localized miRNA Let‐7a‐5p is RNase A–resistant, but degraded by a combination of RNase A/PK treatment suggests exosomal miRNA species also reside on the outer surface of exosomes bound to RBPs/RNPs.     

Functional analysis of the RING-type ubiquitin ligase family of Arabidopsis

Approximately 5% of the Arabidopsis (Arabidopsis thaliana) proteome is predicted to be involved in the ubiquitination/26S proteasome pathway. The majority of these predicted proteins have identity to conserved domains found in E3 ligases, of which there are multiple types. The RING-type E3 is characterized by the presence of a cysteine-rich domain that coordinates two zinc atoms. Database searches followed by extensive manual curation identified 469 predicted Arabidopsis RING domain-containing proteins. In addition to the two canonical RING types (C3H2C3 or C3HC4), additional types of modified RING domains, named RING-v, RING-D, RING-S/T, RING-G, and RING-C2, were identified. The modified RINGs differ in either the spacing between metal ligands or have substitutions at one or more of the metal ligand positions. The majority of the canonical and modified RING domain-containing proteins analyzed were active in in vitro ubiquitination assays, catalyzing polyubiquitination with the E2 AtUBC8. To help identity regions of the proteins that may interact with substrates, domain analyses of the amino acids outside the RING domain classified RING proteins into 30 different groups. Several characterized protein-protein interaction domains were identified, as well as additional conserved domains not described previously. The two largest classes of RING proteins contain either no identifiable domain or a transmembrane domain. The presence of such a large and diverse number of RING domain-containing proteins that function as ubiquitin E3 ligases suggests that target-specific proteolysis by these E3 ligases is a complex and important part of cellular regulation in Arabidopsis.     

Characterization of the mouse brain proteome using global proteomic analysis complemented with cysteinyl-peptide enrichment

We report a global proteomic approach for analyzing brain tissue and for the first time a comprehensive characterization of the whole mouse brain proteome. Preparation of the whole brain sample incorporated a highly efficient cysteinyl-peptide enrichment (CPE) technique to complement a global enzymatic digestion method. Both the global and the cysteinyl-enriched peptide samples were analyzed by SCX fractionation coupled with reversed phase LC-MS/MS analysis. A total of 48,328 different peptides were confidently identified (>98% confidence level), covering 7792 nonredundant proteins ( approximately 34% of the predicted mouse proteome). A total of 1564 and 1859 proteins were identified exclusively from the cysteinyl-peptide and the global peptide samples, respectively, corresponding to 25% and 31% improvements in proteome coverage compared to analysis of only the global peptide or cysteinyl-peptide samples. The identified proteins provide a broad representation of the mouse proteome with little bias evident due to protein pI, molecular weight, and/or cellular localization. Approximately 26% of the identified proteins with gene ontology (GO) annotations were membrane proteins, with 1447 proteins predicted to have transmembrane domains, and many of the membrane proteins were found to be involved in transport and cell signaling. The MS/MS spectrum count information for the identified proteins was used to provide a measure of relative protein abundances. The mouse brain peptide/protein database generated from this study represents the most comprehensive proteome coverage for the mammalian brain to date, and the basis for future quantitative brain proteomic studies using mouse models. The proteomic approach presented here may have broad applications for rapid proteomic analyses of various mouse models of human brain diseases.     

Analysis of Surface-Exposed Outer Membrane Proteins in Helicobacter pylori

More than 50 Helicobacter pylori genes are predicted to encode outer membrane proteins (OMPs), but there has been relatively little experimental investigation of the H. pylori cell surface proteome. In this study, we used selective biotinylation to label proteins localized to the surface of H. pylori, along with differential detergent extraction procedures to isolate proteins localized to the outer membrane. Proteins that met multiple criteria for surface-exposed outer membrane localization included known adhesins, as well as Cag proteins required for activity of the cag type IV secretion system, putative lipoproteins, and other proteins not previously recognized as cell surface components. We identified sites of nontryptic cleavage consistent with signal sequence cleavage, as well as C-terminal motifs that may be important for protein localization. A subset of surface-exposed proteins were highly susceptible to proteolysis when intact bacteria were treated with proteinase K. Most Hop and Hom OMPs were susceptible to proteolysis, whereas Hor and Hof proteins were relatively resistant. Most of the protease-susceptible OMPs contain a large protease-susceptible extracellular domain exported beyond the outer membrane and a protease-resistant domain at the C terminus with a predicted β-barrel structure. These features suggest that, similar to the secretion of the VacA passenger domain, the N-terminal domains of protease-susceptible OMPs are exported through an autotransporter pathway. Collectively, these results provide new insights into the repertoire of surface-exposed H. pylori proteins that may mediate bacterium-host interactions, as well as the cell surface topology of these proteins.     

Mapping and comprehensive analysis of the extracellular and cell surface proteome of the human pathogen Corynebacterium diphtheriae

Secreted proteins of the human pathogen Corynebacterium diphtheriae might be involved in important pathogen-host cell interactions. Here, we present the first systematic reference map of the extracellular and cell surface proteome fractions of the type strain C. diphtheriae C7s(-)tox-. The analysis window of 2-DE covered the pI range from 3 to 10 along with a MW range from 8 to 150 kDa. Computational analysis of the 2-D gels detected almost 150 protein spots in the extracellular proteome fraction and about 80 protein spots of the cell surface proteome. MALDI-TOF-MS and PMF with trypsin unambiguously identified 107 extracellular protein spots and 53 protein spots of the cell surface, representing in total 85 different proteins of C. diphtheriae C7s(-)tox-. Several of the identified proteins are encoded by pathogenicity islands and might represent virulence factors of C. diphtheriae. Additionally, four solute-binding proteins (HmuT, Irp6A, CiuA, and FrgD) of different iron ABC transporters were identified, with the hitherto uncharacterized FrgD protein being the most abundant one of the cell surface proteome of C. diphtheriae C7s(-)tox-.     

Expansion of divergent SEA domains in cell surface proteins and nucleoporin 54

SEA (s ea urchin sperm protein, e nterokinase, a grin) domains, many of which possess autoproteolysis activity, have been found in a number of cell surface and secreted proteins. Despite high sequence divergence, SEA domains were also proposed to be present in dystroglycan based on a conserved autoproteolysis motif and receptor‐type protein phosphatase IA‐2 based on structural similarity. The presence of a SEA domain adjacent to the transmembrane segment appears to be a recurring theme in quite a number of type I transmembrane proteins on the cell surface, such as MUC1, dystroglycan, IA‐2, and Notch receptors. By comparative sequence and structural analyses, we identified dystroglycan‐like proteins with SEA domains in Capsaspora owczarzaki of the Filasterea group, one of the closest single‐cell relatives of metazoans. We also detected novel and divergent SEA domains in a variety of cell surface proteins such as EpCAM, α/ε‐sarcoglycan, PTPRR, collectrin/Tmem27, amnionless, CD34, KIAA0319, fibrocystin‐like protein, and a number of cadherins. While these proteins are mostly from metazoans or their single cell relatives such as choanoflagellates and Filasterea, fibrocystin‐like proteins with SEA domains were found in several other eukaryotic lineages including green algae, Alveolata, Euglenozoa, and Haptophyta, suggesting an ancient evolutionary origin. In addition, the intracellular protein Nucleoporin 54 (Nup54) acquired a divergent SEA domain in choanoflagellates and metazoans.     

Cross-species analyses identify the BNIP-2 and Cdc42GAP homology (BCH) domain as a distinct functional subclass of the CRAL_TRIO/Sec14 superfamily

The CRAL_TRIO protein domain, which is unique to the Sec14 protein superfamily, binds to a diverse set of small lipophilic ligands. Similar domains are found in a range of different proteins including neurofibromatosis type-1, a Ras GTPase-activating Protein (RasGAP) and Rho guanine nucleotide exchange factors (RhoGEFs). Proteins containing this structural protein domain exhibit a low sequence similarity and ligand specificity while maintaining an overall characteristic three-dimensional structure. We have previously demonstrated that the BNIP-2 and Cdc42GAP Homology (BCH) protein domain, which shares a low sequence homology with the CRAL_TRIO domain, can serve as a regulatory scaffold that binds to Rho, RhoGEFs and RhoGAPs to control various cell signalling processes. In this work, we investigate 175 BCH domain-containing proteins from a wide range of different organisms. A phylogenetic analysis with ~100 CRAL_TRIO and similar domains from eight representative species indicates a clear distinction of BCH-containing proteins as a novel subclass within the CRAL_TRIO/Sec14 superfamily. BCH-containing proteins contain a hallmark sequence motif R(R/K)h(R/K)(R/K)NL(R/K)xhhhhHPs ('h' is large and hydrophobic residue and 's' is small and weekly polar residue) and can be further subdivided into three unique subtypes associated with BNIP-2-N, macro- and RhoGAP-type protein domains. A previously unknown group of genes encoding 'BCH-only' domains is also identified in plants and arthropod species. Based on an analysis of their gene-structure and their protein domain context we hypothesize that BCH domain-containing genes evolved through gene duplication, intron insertions and domain swapping events. Furthermore, we explore the point of divergence between BCH and CRAL-TRIO proteins in relation to their ability to bind small GTPases, GAPs and GEFs and lipid ligands. Our study suggests a need for a more extensive analysis of previously uncharacterized BCH, 'BCH-like' and CRAL_TRIO-containing proteins and their significance in regulating signaling events involving small GTPases.     

Proteome analysis on an early transformed human bronchial epithelial cell line,BEP2D,after alpha-particle irradiation

To probe the mechanism of carcinogenesis of lung cancer at the molecular level and to find potential protein markers involved in the early phase of tumorgenesis, differential proteome analysis on primary passage cell line R15H, and early transformed cell line R15H20 derived from (238)Pu alpha-particle irradiation of human papillomavirus (HPV) 18-immortalized human bronchial epithelial cell line (BEP2D), was carried out using two-dimensional electrophoresis (2-DE) and peptide mass fingerprinting (PMF) with matrix-assisted laser desorption/ionisation-time of flight mass spectrometry. Image analysis and Student's t-test (p < 0.05) showed that three protein spots were only expressed in R15H, intensities of 43 protein spots on the gels were altered between R15H and R15H20. Two of the three spots that were only expressed in R15H were identified as high mobility group protein 1. Two proteins decreased in abundance in R15H20 were identified as maspin precursor, a tumor suppressor and aminoacylase-1. Ornithine aminotransferase and peptidyl-prolyl cis-trans isomerase A that were increased in R15H20, were also identified. Relationships between these differentially expressed proteins and the carcinogenesis mechanism of lung cancer are discussed. The protein expression profile of the R15H cell line was also constructed during the study as a reference map for further comparative proteome analysis of the irradiation induced BEP2D cell line. Of the 90 spots analyzed with PMF in the 2-DE gel of R15H cell line, 50 proteins were identified by searching the nonredundant protein database SWISS-PROT/TrEMBL.     

Use of a Probabilistic Motif Search to Identify Histidine Phosphotransfer Domain-Containing Proteins

The wealth of newly obtained proteomic information affords researchers the possibility of searching for proteins of a given structure or function. Here we describe a general method for the detection of a protein domain of interest in any species for which a complete proteome exists. In particular, we apply this approach to identify histidine phosphotransfer (HPt) domain-containing proteins across a range of eukaryotic species. From the sequences of known HPt domains, we created an amino acid occurrence matrix which we then used to define a conserved, probabilistic motif. Examination of various organisms either known to contain (plant and fungal species) or believed to lack (mammals) HPt domains established criteria by which new HPt candidates were identified and ranked. Search results using a probabilistic motif matrix compare favorably with data to be found in several commonly used protein structure/function databases: our method identified all known HPt proteins in the Arabidopsis thaliana proteome, confirmed the absence of such motifs in mice and humans, and suggests new candidate HPts in several organisms. Moreover, probabilistic motif searching can be applied more generally, in a manner both readily customized and computationally compact, to other protein domains; this utility is demonstrated by our identification of histones in a range of eukaryotic organisms.