Classification of 29 families of secondary transport proteins into a single structural class using hydropathy profile analysis

A classification scheme for membrane proteins is proposed that clusters families of proteins into structural classes based on hydropathy profile analysis. The averaged hydropathy profiles of protein families are taken as fingerprints of the 3D structure of the proteins and, therefore, are able to detect more distant evolutionary relationships than amino acid sequences. A procedure was developed in which hydropathy profile analysis is used initially as a filter in a BLAST search of the NCBI protein database. The strength of the procedure is demonstrated by the classification of 29 families of secondary transporters into a single structural class, termed ST[3]. An exhaustive search of the database revealed that the 29 families contain 568 unique sequences. The proteins are predominantly from prokaryotic origin and most of the characterized transporters in ST[3] transport organic and inorganic anions and a smaller number are Na(+)/H(+) antiporters. All modes of energy coupling (symport, antiport, uniport) are found in structural class ST[3]. The relevance of the classification for structure/function prediction of uncharacterised transporters in the class is discussed.     

Profiling human brain proteome by multi-dimensional separations coupled with MS

In our initial attempt to analyze the human brain proteome, we applied multi-dimensional protein separation and identification techniques using a combination of sample fractionation, 1-D SDS-PAGE, and MS analysis. The complexity of human brain proteome requires multiple fractionation strategies to extend the range and total number of proteins identified. According to the method of Klose (Methods Mol. Biol. 1999, 112, 67), proteins of the temporal lobe of human brain were fractionated into (i) cytoplasmic and nucleoplasmic, (ii) membrane and other structural, and (iii) DNA-binding proteins. Each fraction was then separated by SDS-PAGE, and the resulting gel line was cut into approximately 50 bands. After trypsin digestion, the resulting peptides from each band were analyzed by RP-LC/ESI-MS/MS using an LTQ spectrometer. The SEQUEST search program, which searched against the IPI database, was used for peptide sequence identification, and peptide sequences were validated by reversed sequence database search and filtered by the Protein Hit Score. Ultimately, 1533 proteins could be detected from the human brain. We classified the identified proteins according to their distribution on cellular components. Among these proteins, 24% were membrane proteins. Our results show that the multiple separation strategy is effective for high-throughput characterization of proteins from complex proteomic mixtures.     

Sequence and hydropathy profile analysis of two classes of secondary transporters

A structural class in the MemGen classification of membrane proteins is a set of evolutionary related proteins sharing a similar global fold. A structural class contains both closely related pairs of proteins for which homology is clear from sequence comparison and very distantly related pairs, for which it is not possible to establish homology based on sequence similarity alone. In the latter case the evolutionary link is based on hydropathy profile analysis. Here, we use these evolutionary related sets of proteins to analyze the relationship between E-values in BLAST searches, sequence similarities in multiple sequence alignments and structural similarities in hydropathy profile analyses. Two structural classes of secondary transporters termed ST[3], which includes the Ion Transporter (IT) superfamily and ST[4], which includes the DAACS family (TC# 2.A.23) were extracted from the NCBI protein database. ST[3] contains 2051 unique sequences distributed over 32 families and 59 subfamilies. ST[4] is a smaller class containing 399 unique sequences distributed over 2 families and 7 subfamilies. One subfamily in ST[4] contains a new class of binding protein dependent secondary transporters. Comparison of the averaged hydropathy profiles of the subfamilies in ST[3] and ST[4] revealed that the two classes represent different folds. Divergence of the sequences in ST[4] is much smaller than observed in ST[3], suggesting different constraints on the proteins during evolution. Analysis of the correlation between the evolutionary relationship of pairs of proteins in a class and the BLAST E-value revealed that: (i) the BLAST algorithm is unable to pick up the majority of the links between proteins in structural class ST[3], (ii) "low complexity filtering" and "composition based statistics" improve the specificity, but strongly reduce the sensitivity of BLAST searches for distantly related proteins, indicating that these filters are too stringent for the proteins analyzed, and (iii) the E-value cut-off, which may be used to evaluate evolutionary significance of a hit in a BLAST search is very different for the two structural classes of membrane proteins.     

Sequence and hydropathy profile analysis of two classes of secondary transporters

A structural class in the MemGen classification of membrane proteins is a set of evolutionary related proteins sharing a similar global fold. A structural class contains both closely related pairs of proteins for which homology is clear from sequence comparison and very distantly related pairs, for which it is not possible to establish homology based on sequence similarity alone. In the latter case the evolutionary link is based on hydropathy profile analysis. Here, we use these evolutionary related sets of proteins to analyze the relationship between E-values in BLAST searches, sequence similarities in multiple sequence alignments and structural similarities in hydropathy profile analyses. Two structural classes of secondary transporters termed ST[3], which includes the Ion Transporter (IT) superfamily and ST[4], which includes the DAACS family (TC# 2.A.23) were extracted from the NCBI protein database. ST[3] contains 2051 unique sequences distributed over 32 families and 59 subfamilies. ST[4] is a smaller class containing 399 unique sequences distributed over 2 families and 7 subfamilies. One subfamily in ST[4] contains a new class of binding protein dependent secondary transporters. Comparison of the averaged hydropathy profiles of the subfamilies in ST[3] and ST[4] revealed that the two classes represent different folds. Divergence of the sequences in ST[4] is much smaller than observed in ST[3], suggesting different constraints on the proteins during evolution. Analysis of the correlation between the evolutionary relationship of pairs of proteins in a class and the BLAST E-value revealed that: (i) the BLAST algorithm is unable to pick up the majority of the links between proteins in structural class ST[3], (ii) ‘low complexity filtering’ and ‘composition based statistics’ improve the specificity, but strongly reduce the sensitivity of BLAST searches for distantly related proteins, indicating that these filters are too stringent for the proteins analyzed, and (iii) the E-value cut-off, which may be used to evaluate evolutionary significance of a hit in a BLAST search is very different for the two structural classes of membrane proteins.     

MPtopo: A database of membrane protein topology

The reliability of the transmembrane (TM) sequence assignments for membrane proteins (MPs) in standard sequence databases is uncertain because the vast majority are based on hydropathy plots. A database of MPs with dependable assignments is necessary for developing new computational tools for the prediction of MP structure. We have therefore created MPtopo, a database of MPs whose topologies have been verified experimentally by means of crystallography, gene fusion, and other methods. Tests using MPtopo strongly validated four existing MP topology-prediction algorithms. MPtopo is freely available over the internet and can be queried by means of an SQL-based search engine.     

The uncoupling protein from brown fat mitochondria is related to the mitochondrial ADP/ATP carrier. Analysis of sequence homologies and of folding of the protein in the membrane.

We report here, for the first time, the primary structure of uncoupling protein as established by amino acid sequencing. Like the ADP/ATP carrier, this protein has a tripartite structure comprising three similar sequences of approximately 100 residues each. These six 'repeats' exhibit striking conservation of several residues, in particular glycine and proline, at possible structurally strategic positions. Although the two proteins differ strongly in their amino acid composition, their sequences are distantly homologous. Three membrane-spanning alpha-helices can be deduced from hydropathy plots. A modified plot accounting for amphiphilic helices indicates 5-6 such alpha-segments. In addition an amphiphilic beta-strand of membrane-spanning length can be discerned. The tripartite sequence structure is also distinctly reflected in the hydropathy distribution. Based on the membrane disposition of the segments of the ADP/ATP carrier, a model for the transmembrane folding path of the polypeptide chain of the uncoupling protein is proposed.     

Using hydropathy features for function prediction of membrane proteins

A novel alignment-free method for computing functional similarity of membrane proteins based on features of hydropathy distribution is presented. The features of hydropathy distribution are used to represent protein families as hydropathy profiles. The profiles statistically summarize the hydropathy distribution of member proteins. The summation is made by using hydropathy features that numerically represent structurally/functionally significant portions of protein sequences. The hydropathy profiles are numerical vectors that are points in a high dimensional ‘hydropathy’ space. Their similarities are identified by projection of the space onto principal axes. Here, the approach is applied to the secondary transporters. The analysis using the presented approach is validated by the standard classification of the secondary transporters. The presented analysis allows for prediction of function attributes for proteins of uncharacterized families of secondary transporters. The results obtained using the presented analysis may help to characterize unknown function attributes of secondary transporters. They also show that analysis of hydropathy distribution can be used for function prediction of membrane proteins.     

An Ochrobactrum anthropi gene conferring paraquat resistance to the heterologous host Escherichia coli

A new gene, pqrA, conferring paraquat resistance to the heterologous host Escherichia coli, from a chromosomal DNA library of Ochrobactrum anthropi JW2, was cloned and analyzed. Cells of E. coli transformed with a plasmid carrying the pqrA gene showed elevated resistance to paraquat, but not to hydrogen peroxide. The predicted amino acid sequence of the PqrA polypeptide showed 71% identity with mll7495 hypothetical membrane protein in Mesorhizobium loti, 49% identity with PA2269 protein in Pseudomonas aeruginosa, and significant identity with other previously reported drug transport proteins. The hydropathy pattern of the PqrA polypeptide showed a significant homology to those of 12-transmembrane-segment (TMS) family export proteins. Immunoblot analysis demonstrated that the PqrA protein found in the membrane protein fraction of O. anthropi JW2 has a molecular mass of 42 kDa. These results suggest that the PqrA protein is a membrane protein that plays an important role in protecting cells against paraquat toxicity.     

Using hydropathy features for function prediction of membrane proteins

A novel alignment-free method for computing functional similarity of membrane proteins based on features of hydropathy distribution is presented. The features of hydropathy distribution are used to represent protein families as hydropathy profiles. The profiles statistically summarize the hydropathy distribution of member proteins. The summation is made by using hydropathy features that numerically represent structurally/functionally significant portions of protein sequences. The hydropathy profiles are numerical vectors that are points in a high dimensional 'hydropathy' space. Their similarities are identified by projection of the space onto principal axes. Here, the approach is applied to the secondary transporters. The analysis using the presented approach is validated by the standard classification of the secondary transporters. The presented analysis allows for prediction of function attributes for proteins of uncharacterized families of secondary transporters. The results obtained using the presented analysis may help to characterize unknown function attributes of secondary transporters. They also show that analysis of hydropathy distribution can be used for function prediction of membrane proteins.     

The LysE superfamily: topology of the lysine exporter LysE of Corynebacterium glutamicum, a paradyme for a novel superfamily of transmembrane solute translocators

In Corynebacterium glutamicum the LysE carrier protein exhibits the unique function of exporting L-lysine. We here analyze the membrane topology of LysE, a protein of 236 amino acyl residues, using PhoA- and LacZ-fusions. The amino-terminal end of LysE is located in the cytoplasm whereas the carboxy-terminal end is found in the periplasm. Although 6 hydrophobic domains were identified based on hydropathy analyses, only five transmembrane spanning helices appear to be present. The additional hydrophobic segment may dip into the membrane or be surface localized. We show that LysE is a member of a family of proteins found, for example, in Escherichia coil, Bacillus subtilis, Mycobacterium tuberculosis and Helicobacter pylori. This family, which we have designated the LysE family, is distantly related to two additional protein families which we have designated the YahN and CadD families. These three families, the members of which exhibit similar sizes, hydropathy profiles, and sequence motifs comprise the LysE superfamily. Functionally characterized members of the LysE superfamily export L-lysine, cadmium and possibly quarternary amines. We suggest that LysE superfamily members will prove to catalyze export of a variety of biologically important solutes.     

Protein biophysical properties that correlate with crystallization success in Thermotoga maritima: maximum clustering strategy for structural genomics

Cost and time reduction are two of the driving forces in the development of new strategies for protein crystallization and subsequent structure determination. Here, we report the analysis of the Thermotoga maritima proteome, in which we compare the proteins that were successfully expressed, purified and crystallized versus the rest of the proteome. This set of almost 500 proteins represents one of the largest, internally consistent, protein expression and crystallization datasets available. The analysis shows that individual parameters, such as isoelectric point, sequence length, average hydropathy, low complexity regions (SEG), and combinations of these biophysical properties for crystallized proteins define a distinct subset of the T. maritima proteome. The distribution profiles of the various biophysical properties in the expression/crystallization set are then used to extract rules to improve target selection and improve the efficiency and output of structural genomics, as well as general structural biology efforts.     

Developing the wool proteome

The wool proteome has been largely uncharted due to a lack of database coverage, poor protein extractability and dynamic range issues. Yet, investigating correlations between wool physical properties and protein content, or characterising UV-, heat- or processing-induced protein damage requires the availability of an identifiable and identified proteome.In this study we have achieved unprecedented wool proteome identification through a strategy of comprehensive data acquisition, iterative protein identification/validation and concurrent augmentation of the sequence database. Data acquisition comprised a range of different hyphenated MS techniques including LC–MS/MS, LC–MALDI, 2D-LC–MS/MS and SDS-PAGE LC–MS. Using iterative searching of databases and search result combination using ProteinScape, a systematic expansion of identifiable proteins in the sequence database was achieved. This was followed by extensive validation and rationalisation of the protein identifications. In total, 72 complete and 30 partial ovine-specific protein sequences were added to the database, and 113 wool proteins were identified.Enhanced access to ovine-specific protein identification and characterisation will facilitate all wool fibre protein chemistry and proteomics research.     

Human plasma proteome analysis by reversed sequence database search and molecular weight correlation based on a bacterial proteome analysis

In shotgun proteomics, proteins can be fractionated by 1-D gel electrophoresis and digested into peptides, followed by liquid chromatography to separate the peptide mixture. Mass spectrometry generates hundreds of thousands of tandem mass spectra from these fractions, and proteins are identified by database searching. However, the search scores are usually not sufficient to distinguish the correct peptides. In this study, we propose a confident protein identification method for high-throughput analysis of human proteome. To build a filtering protocol in database search, we chose Pseudomonas putida KT2440 as a reference because this bacterial proteome contains fewer modifications and is simpler than the human proteome. First, the P. putida KT2440 proteome was filtered by reversed sequence database search and correlated by the molecular weight in 1-D-gel band positions. The characterization protocol was then applied to determine the criteria for clustering of the human plasma proteome into three different groups. This protein filtering method, based on bacterial proteome data analysis, represents a rapid way to generate higher confidence protein list of the human proteome, which includes some of heavily modified and cleaved proteins.     

Proteome Analysis of Light-induced Proteins in Synechocystis sp. PCC 6803: Identification of Proteins Separated by 2D-PAGE Using N-terminal Sequencing and MALDI-TOF MS

The cyanobacterium Synechocystis sp. PCC 6803 is an ideal model organism for the proteome study of light-induced gene expression because the whole genomic sequence has been determined. The soluble proteins extracted from light- and dark-cultured cells were separated by two-dimensional polyacrylamide gel electrophoresis. Light-induced protein spots electroblotted on a polyvinyldiene difluoride membrane were analyzed by N-terminal Edman sequence determination and followed by CyanoBase. The tryptic digests of some proteins were also confirmed by matrix-assisted laser desorption ionization/time-of-flight (MALDI-TOF) and MS-Fit search. Interestingly, eight proteins were related to photosynthesis and respiration (RbcS/L, CbbA, Gap2, AtpB, CpcB, PsbO, and PsbU). Four proteins (SodB, DnaK, GroEL2, and Tig) were involved in cellular processes and the functions of another two proteins (rehydrin and membrane protein) were unknown. The proteome analysis by N-terminal Edman sequencing and MALDI-TOF enabled us to characterize one-shot protein profiles expressed under different physiological conditions.     

Nucleotide sequence of the gene for the ferrienterochelin receptor FepA in Escherichia coli. Homology among outer membrane receptors that interact with TonB

We have determined the nucleotide sequence of the Escherichia coli fepA gene, which codes for the outer membrane receptor for ferrienterochelin and colicins B and D. The predicted FepA polypeptide has a molecular weight of 79,908 and consists of 723 amino acids. A 22-amino acid leader or signal peptide preceded the mature protein. With respect to overall composition, hydropathy, net charge and distribution of nonpolar segments, the FepA polypeptide was typical of other E. coli outer membrane proteins, except that FepA contained 2 cysteine residues. Comparison of the deduced amino acid sequence of FepA with that of three other TonB-dependent receptors (BtuB, FhuA, and IutA) revealed only a few regions of sequence homology; one of these included the amino-termini. An amino acid substitution within the conserved amino-terminal region of BtuB resulted in production of a receptor that had normal binding functions but was incapable of energy-dependent transport of vitamin B12. This result suggests that the amino-terminal end of these four polypeptides is involved in interaction with the TonB protein or another step of energy transduction. Three other regions of homology were shared among the four proteins: one near residues 50 to 70, another at about residue 100 to 140, and the last between 20 and 40 amino acid residues from the carboxyl terminus. The function of these three regions remains speculative.     

Structural features of the glutamate transporter family.

Neuronal and glial glutamate transporters remove the excitatory neurotransmitter glutamate from the synaptic cleft and thus prevent neurotoxicity. The proteins belong to a large and widespread family of secondary transporters, including bacterial glutamate, serine, and C4-dicarboxylate transporters; mammalian neutral-amino-acid transporters; and an increasing number of bacterial, archaeal, and eukaryotic proteins that have not yet been functionally characterized. Sixty members of the glutamate transporter family were found in the databases on the basis of sequence homology. The amino acid sequences of the carriers have diverged enormously. Homology between the members of the family is most apparent in a stretch of approximately 150 residues in the C-terminal part of the proteins. This region contains four reasonably well-conserved sequence motifs, all of which have been suggested to be part of the translocation pore or substrate binding site. Phylogenetic analysis of the C-terminal stretch revealed the presence of five subfamilies with characterized members: (i) the eukaryotic glutamate transporters, (ii) the bacterial glutamate transporters, (iii) the eukaryotic neutral-amino-acid transporters, (iv) the bacterial C4-dicarboxylate transporters, and (v) the bacterial serine transporters. A number of other subfamilies that do not contain characterized members have been defined. In contrast to their amino acid sequences, the hydropathy profiles of the members of the family are extremely well conserved. Analysis of the hydropathy profiles has suggested that the glutamate transporters have a global structure that is unique among secondary transporters. Experimentally, the unique structure of the transporters was recently confirmed by membrane topology studies. Although there is still controversy about part of the topology, the most likely model predicts the presence of eight membrane-spanning alpha-helices and a loop-pore structure which is unique among secondary transporters but may resemble loop-pores found in ion channels. A second distinctive structural feature is the presence of a highly amphipathic membrane-spanning helix that provides a hydrophilic path through the membrane. Recent data from analysis of site-directed mutants and studies on the mechanism and pharmacology of the transporters are discussed in relation to the structural model.     

Structural Features of the Glutamate Transporter Family

Neuronal and glial glutamate transporters remove the excitatory neurotransmitter glutamate from the synaptic cleft and thus prevent neurotoxicity. The proteins belong to a large and widespread family of secondary transporters, including bacterial glutamate, serine, and C₄-dicarboxylate transporters; mammalian neutral-amino-acid transporters; and an increasing number of bacterial, archaeal, and eukaryotic proteins that have not yet been functionally characterized. Sixty members of the glutamate transporter family were found in the databases on the basis of sequence homology. The amino acid sequences of the carriers have diverged enormously. Homology between the members of the family is most apparent in a stretch of approximately 150 residues in the C-terminal part of the proteins. This region contains four reasonably well-conserved sequence motifs, all of which have been suggested to be part of the translocation pore or substrate binding site. Phylogenetic analysis of the C-terminal stretch revealed the presence of five subfamilies with characterized members: (i) the eukaryotic glutamate transporters, (ii) the bacterial glutamate transporters, (iii) the eukaryotic neutral-amino-acid transporters, (iv) the bacterial C₄-dicarboxylate transporters, and (v) the bacterial serine transporters. A number of other subfamilies that do not contain characterized members have been defined. In contrast to their amino acid sequences, the hydropathy profiles of the members of the family are extremely well conserved. Analysis of the hydropathy profiles has suggested that the glutamate transporters have a global structure that is unique among secondary transporters. Experimentally, the unique structure of the transporters was recently confirmed by membrane topology studies. Although there is still controversy about part of the topology, the most likely model predicts the presence of eight membrane-spanning α-helices and a loop-pore structure which is unique among secondary transporters but may resemble loop-pores found in ion channels. A second distinctive structural feature is the presence of a highly amphipathic membrane-spanning helix that provides a hydrophilic path through the membrane. Recent data from analysis of site-directed mutants and studies on the mechanism and pharmacology of the transporters are discussed in relation to the structural model.     

Proteome Analysis of Light-induced Proteins in Synechocystis sp. PCC 6803: Identification of Proteins Separated by 2D-PAGE Using N-terminal Sequencing and MALDI-TOF MS

The cyanobacterium Synechocystis sp. PCC 6803 is an ideal model organism for the proteome study of light-induced gene expression because the whole genomic sequence has been determined. The soluble proteins extracted from light- and dark-cultured cells were separated by two-dimensional polyacrylamide gel electrophoresis. Light-induced protein spots electroblotted on a polyvinyldiene difluoride membrane were analyzed by N-terminal Edman sequence determination and followed by CyanoBase. The tryptic digests of some proteins were also confirmed by matrix-assisted laser desorption ionization/time-of-flight (MALDI-TOF) and MS-Fit search. Interestingly, eight proteins were related to photosynthesis and respiration (RbcS/L, CbbA, Gap2, AtpB, CpcB, PsbO, and PsbU). Four proteins (SodB, DnaK, GroEL2, and Tig) were involved in cellular processes and the functions of another two proteins (rehydrin and membrane protein) were unknown. The proteome analysis by N-terminal Edman sequencing and MALDI-TOF enabled us to characterize one-shot protein profiles expressed under different physiological conditions.     

Prediction and analysis of paralogous proteins in Trichomonas vaginalis genome

Trichomonas vaginalis causes trichomoniasis, second most sexually transmitted disease. The genome sequence draft of T. vaginalis was published by The Institute of Genomic Research reveals an abnormally large genome size of 160 Mb. It was speculated that a significant portion of the proteome contains paralogous proteins. The present study was aimed at identification and analysis of the paralogous proteins. The all against all search approach is used to identify the paralogous proteins. The dataset of proteins was retrieved from TIGR and TrichDB FTP server. The BLAST-P program performed all against all database searches against the protein database of Trichomonas vaginalis available at NCBI genome database. In the present study about 50,000 proteins were searched where 2,700 proteins were found to be paralogous under the rigid selection criteria. The Pfam database search has identified significant number of paralogous proteins which were further categorized among different 1496 paralogous protein in pfam families, 1027 paralogous protein contains domain, 60 proteins were having different repeats and 1092 paralogous protein sequences of clans. Such identification and functional annotation of paralogous proteins will also help in removing paralogous proteins from possible drug targets in future. Presence of huge number of paralogous proteins across wide range of gene families and domains may be one of the possible mechanisms involved in the T. vaginalis genome expansion and evolution.     

Protein sequence modules

Conserved protein sequence segments are commonly believed to correspond to functional sites in the protein sequence. A novel approach is proposed to profile the changing degree of conservation along the protein sequence, by evaluating the occurrence frequencies of all short oligopeptides of the given sequence in a large proteome database. Thus, a protein sequence conservation profile can be plotted for every protein. The profile indicates where along the sequences the potential functional (conserved) sites are located. The corresponding oligopeptides belonging to the sites are very frequent across many prokaryotic species. Analysis of a representative set of such profiles reveals a common feature of all examined proteins: they consist of sequence modules represented by the peaks of conservation. Typical size of the modules (peak-to-peak distance) is 25-30 amino acid residues.