Prediction of Membrane Protein Topology Utilizing Multiple Sequence Alignments

A technique for prediction of protein membrane toplogy (intra- and extraceullular sidedness) has been developed. Membrane-spanning segments are first predicted using an algorithm based upon multiply aligned amino acid sequences. The compositional differences in the protein segments exposed at each side of the membrane are then investigated. The ratios are calculated for Asn, Asp, Gly, Phe, Pro, Trp, Tyr, and Val, mostly found on the extracellular side, and for Ala, Arg, Cys, and Lys, mostly occurring on the intracellular side. The consensus over these 12 residue distributions is used for sidedness prediction. The method was developed with a set of 42 protein families for which all but one were correctly predicted with the new algorithm. This represents an improvement over previous techniques. The new method, applied to a set of 12 membrane protein families different from the test set and with recently determined topologies, performed well, with 11 of 12 sidedness assignments agreeing with experimental results. The method has also been applied to several membrane protein families for which the topology has yet to be determined. An electronic prediction service is available at the E-mail address tmap@embl-heidelberg.de and on WWW via http://www.emblheidelberg.de.     

Sequence similarity as a predictor of the transmembrane topology of membrane-intrinsic subunits of bacterial respiratory chain enzymes

Integral membrane proteins usually have a predominantly alpha-helical secondary structure in which transmembrane segments are connected by membrane-extrinsic loops. Although a number of membrane protein structures have been reported in recent years, in most cases transmembrane topologies are initially predicted using a variety of theoretical techniques, including hydropathy analyses and the "positive inside" rule. We have explored the use of plots of the distribution of sequence similarity within families of membrane proteins comprising homeomorphic domains as a new method for the prediction/verification of the orientation of transmembrane topology models within certain families of multimeric respiratory chain enzymes. Within such proteins, analyses of sequence similarity can: i) identify heme and/or quinol binding sites; ii) identify potential electron-transfer conduits to/from prosthetic groups; and iii) locate regions defining potential subunit-subunit interactions. We mined emerging bioinformatic data for sequences of 11 families of membrane-intrinsic proteins that are part of multimeric respiratory chain complexes that also have membrane-extrinsic subunits. The sequences of each family were then aligned and the resultant alignments converted into a graphical format recording an empirical measure of the sequence similarity plotted versus residue position. In each case, this plot was compared to the predicted transmembrane topology. With one exception, there is a strong correlation between the existence     

Enhanced membrane protein topology prediction using a hierarchical classification method and a new scoring function

The prediction of transmembrane (TM) helix and topology provides important information about the structure and function of a membrane protein. Due to the experimental difficulties in obtaining a high-resolution model, computational methods are highly desirable. In this paper, we present a hierarchical classification method using support vector machines (SVMs) that integrates selected features by capturing the sequence-to-structure relationship and developing a new scoring function based on membrane protein folding. The proposed approach is evaluated on low- and high-resolution data sets with cross-validation, and the topology (sidedness) prediction accuracy reaches as high as 90%. Our method is also found to correctly predict both the location of TM helices and the topology for 69% of the low-resolution benchmark set. We also test our method for discrimination between soluble and membrane proteins and achieve very low overall false positive (0.5%) and false negative rates (0 to approximately 1.2%). Lastly, the analysis of the scoring function suggests that the topogeneses of single-spanning and multispanning TM proteins have different levels of complexity, and the consideration of interloop topogenic interactions for the latter is the key to achieving better predictions. This method can facilitate the annotation of membrane proteomes to extract useful structural and functional information. It is publicly available at http://bio-cluster.iis.sinica.edu.tw/~bioapp/SVMtop.     

Topology prediction of Brucella abortus Omp2b and Omp2a porins after critical assessment of transmembrane beta strands prediction by several secondary structure prediction methods

In order to propose a reliable model for Brucella porin topology, several structure prediction methods were evaluated in their ability to predict porin topology. Four porins of known structure were selected as test-cases and their secondary structure delineated. The specificity and sensitivity of 11 methods were separately evaluated. Our critical assessment shows that some secondary structure prediction methods (PHD, Dsc, Sopma) originally designed to predict globular protein structure are useful on porin topology prediction. The overall best prediction is obtained by combining these three "generalist" methods with a transmembrane beta strand prediction technique. This "consensus" method was applied to Brucella porins Omp2b and Omp2a, sharing no sequence homology with any other porin. The predicted topology is a 16-stranded antiparallel beta barrel with Omp2a showing a higher number of negatively charged residue in the exposed loops than Omp2b. Experiments are in progress to validate the proposed topology and the functional hypotheses. The ability of the proposed consensus method to predict topology of complex outer membrane protein is briefly discussed.     

WRF-TMH: predicting transmembrane helix by fusing composition index and physicochemical properties of amino acids

Membrane protein is the prime constituent of a cell, which performs a role of mediator between intra and extracellular processes. The prediction of transmembrane (TM) helix and its topology provides essential information regarding the function and structure of membrane proteins. However, prediction of TM helix and its topology is a challenging issue in bioinformatics and computational biology due to experimental complexities and lack of its established structures. Therefore, the location and orientation of TM helix segments are predicted from topogenic sequences. In this regard, we propose WRF-TMH model for effectively predicting TM helix segments. In this model, information is extracted from membrane protein sequences using compositional index and physicochemical properties. The redundant and irrelevant features are eliminated through singular value decomposition. The selected features provided by these feature extraction strategies are then fused to develop a hybrid model. Weighted random forest is adopted as a classification approach. We have used two benchmark datasets including low and high-resolution datasets. tenfold cross validation is employed to assess the performance of WRF-TMH model at different levels including per protein, per segment, and per residue. The success rates of WRF-TMH model are quite promising and are the best reported so far on the same datasets. It is observed that WRF-TMH model might play a substantial role, and will provide essential information for further structural and functional studies on membrane proteins. The accompanied web predictor is accessible at http://111.68.99.218/WRF-TMH/.     

A guideline to proteome-wide α-helical membrane protein topology predictions

For current state-of-the-art methods, the prediction of correct topology of membrane proteins has been reported to be above 80%. However, this performance has only been observed in small and possibly biased data sets obtained from protein structures or biochemical assays. Here, we test a number of topology predictors on an "unseen" set of proteins of known structure and also on four "genome-scale" data sets, including one recent large set of experimentally validated human membrane proteins with glycosylated sites. The set of glycosylated proteins is also used to examine the ability of prediction methods to separate membrane from nonmembrane proteins. The results show that methods utilizing multiple sequence alignments are overall superior to methods that do not. The best performance is obtained by TOPCONS, a consensus method that combines several of the other prediction methods. The best methods to distinguish membrane from nonmembrane proteins belong to the "Phobius" group of predictors. We further observe that the reported high accuracies in the smaller benchmark sets are not quite maintained in larger scale benchmarks. Instead, we estimate the performance of the best prediction methods for eukaryotic membrane proteins to be between 60% and 70%. The low agreement between predictions from different methods questions earlier estimates about the global properties of the membrane proteome. Finally, we suggest a pipeline to estimate these properties using a combination of the best predictors that could be applied in large-scale proteomics studies of membrane proteins.     

The first cytoplasmic loop of the mannitol permease from Escherichia coli is accessible for sulfhydryl reagents from the periplasmic side of the membrane

The mannitol permease (EII(Mtl)) from Escherichia coli couples mannitol transport to phosphorylation of the substrate. Renewed topology prediction of the membrane-embedded C domain suggested that EII(Mtl) contains more membrane-embedded segments than the six proposed previously on the basis of a PhoA fusion study. Cysteine accessibility was used to confirm this notion. Since cysteine 384 in the cytoplasmic B domain is crucial for the phosphorylation activity of EII(Mtl), all cysteine mutants contained this activity-linked cysteine residue in addition to those introduced for probing the membrane topology of the protein. To distinguish between the activity-linked cysteine and the probed cysteine, either trypsin was used to specifically digest the two cytoplasmic domains (A and B), thereby removing Cys384, or Cys384 was protected by phosphorylation from alkylation by N-ethylmaleimide (NEM). Our data show that upon phosphorylation EII(Mtl) undergoes major conformational changes, whereby residues in the putative first cytoplasmic loop become accessible to NEM. Other residues in this loop were accessible to NEM in intact cells and inside-out membrane vesicles, but cysteine residues at these positions only reacted with the membrane-impermeable sulfhydryl reagent from the periplasmic side of the protein. These and other results suggest that the predicted loop between TM2 and TM3 may fold back into the membrane and form part of the translocation path.     

kPROT: a knowledge-based scale for the propensity of residue orientation in transmembrane segments. Application to membrane protein structure prediction

Modeling of integral membrane proteins and the prediction of their functional sites requires the identification of transmembrane (TM) segments and the determination of their angular orientations. Hydrophobicity scales predict accurately the location of TM helices, but are less accurate in computing angular disposition. Estimating lipid-exposure propensities of the residues from statistics of solved membrane protein structures has the disadvantage of relying on relatively few proteins. As an alternative, we propose here a scale of knowledge-based Propensities for Residue Orientation in Transmembrane segments (kPROT), derived from the analysis of more than 5000 non-redundant protein sequences. We assume that residues that tend to be exposed to the membrane are more frequent in TM segments of single-span proteins, while residues that prefer to be buried in the transmembrane bundle interior are present mainly in multi-span TMs. The kPROT value for each residue is thus defined as the logarithm of the ratio of its proportions in single and multiple TM spans. The scale is refined further by defining it for three discrete sections of the TM segment; namely, extracellular, central, and intracellular. The capacity of the kPROT scale to predict angular helical orientation was compared to that of alternative methods in a benchmark test, using a diversity of multi-span alpha-helical transmembrane proteins with a solved 3D structure. kPROT yielded an average angular error of 41 degrees, significantly lower than that of alternative scales (62 degrees -68 degrees ). The new scale thus provides a useful general tool for modeling and prediction of functional residues in membrane proteins. A WWW server (http://bioinfo.weizmann.ac.il/kPROT) is available for automatic helix orientation prediction with kPROT.     

Prediction of partial membrane protein topologies using a consensus approach

We have developed a method to reliably identify partial membrane protein topologies using the consensus of five topology prediction methods. When evaluated on a test set of experimentally characterized proteins, we find that approximately 90% of the partial consensus topologies are correctly predicted in membrane proteins from prokaryotic as well as eukaryotic organisms. Whole-genome analysis reveals that a reliable partial consensus topology can be predicted for approximately 70% of all membrane proteins in a typical bacterial genome and for approximately 55% of all membrane proteins in a typical eukaryotic genome. The average fraction of sequence length covered by a partial consensus topology is 44% for the prokaryotic proteins and 17% for the eukaryotic proteins in our test set, and similar numbers are found when the algorithm is applied to whole genomes. Reliably predicted partial topologies may simplify experimental determinations of membrane protein topology.     

Membrane topology of human ABC proteins

In this review, we summarize the currently available information on the membrane topology of some key members of the human ABC protein subfamilies, and present the predicted domain arrangements. In the lack of high-resolution structures for eukaryotic ABC transporters this topology is based only on prediction algorithms and biochemical data for the location of various segments of the polypeptide chain, relative to the membrane. We suggest that topology models generated by the available prediction methods should only be used as guidelines to provide a basis of experimental strategies for the elucidation of the membrane topology.     

α-螺旋跨膜蛋白拓扑结构预测方法的评价

在基因组数据中,有20%～30%的产物被预测为跨膜蛋白,本文通过对膜蛋白拓扑结构预测方法进行分析,并评价其结果,为选择更合适的拓扑结构预测方法预测膜蛋白结构。通过对目前已有的拓扑结构预测方法的评价分析,可以为我们在实际工作中提供重要的参考。比如对一个未知拓扑结构的跨膜蛋白序列,我们可以先进行是否含有信号肽的预测,参考Polyphobius和SignalP两种方法,若两种方法预测结果不一致,综合上述对两种方法的评价,Polyphobius预测的综合能力较好,可取其预测的结果,一旦确定含有信号肽,则N端必然位于膜外侧。然后结合序列的长度,判断蛋白是单跨膜还是多重跨膜,即可参照上述评价结果,选择合适的拓扑结构预测方法进行预测。     

Membrane topology of the Streptococcus pneumoniae FtsW division protein

The topology of FtsW from Streptococcus pneumoniae, an essential membrane protein involved in bacterial cell division, was predicted by computational methods and probed by the alkaline phosphatase fusion and cysteine accessibility techniques. Consistent results were obtained for the seven N-terminal membrane-spanning segments. However, the results from alkaline phosphatase fusions did not confirm the hydropathy analysis of the C-terminal part of FtsW, whereas the accessibility of introduced cysteine residues was in agreement with the theoretical prediction. Based on the combined results, we propose the first topological model of FtsW, featuring 10 membrane-spanning segments, a large extracytoplasmic loop, and both N and C termini located in the cytoplasm.     

Orientation of Small Multidrug Resistance Transporter Subunits in the Membrane: Correlation with the Positive-Inside Rule

Small multidrug resistance (SMR) transport proteins provide a model for the evolution of larger two-domain transport proteins. The orientation in the membrane of 27 proteins from the SMR family was determined using the reporter fusion technique. Nine members were encoded monocistronically (singles) and shown to insert in both orientations (dual topology). Eighteen members were encoded in pairs on the chromosome and shown to insert in fixed orientations; the two proteins in each pair invariably had opposite orientations in the membrane. Interaction between the two proteins in pairs was demonstrated by copurification. The orientation in the membrane of either protein in the pair was affected only marginally by the presence of the other protein.For the proteins in pairs, the orientation in the membrane correlated well with the distribution of positively charges residues (R + K) over the cytoplasmic and extracellular loops (positive-inside rule). In contrast, dual-topology insertion of the singles was predicted less well by the positive-inside rule. Three singles were predicted to insert in a single orientation with the N-terminus and the C-terminus at the extracellular side of the membrane. Analysis of charge distributions suggests the requirement of a threshold number of charges in the cytoplasmic loops for the positive-inside rule to be of predictive value. It is concluded that a combined analysis of gene organization on the chromosome and phylogeny is sufficient to distinguish between fixed or dual topology of SMR members and, probably, similar types of membrane proteins. The positive-inside rule can be used to predict the orientation of members in pairs, but is not suitable as a sole predictor of dual topology.     

Transmembrane segment (TMS) VIII of the Na(+)/Citrate transporter CitS requires downstream TMS IX for insertion in the Escherichia coli membrane.

The amino acid sequence of the sodium ion-dependent citrate transporter CitS of K. pneumoniae contains 12 hydrophobic stretches that could form membrane-spanning segments. A previous analysis of the membrane topology in Escherichia coli using the PhoA gene fusion technique indicated that only nine of these hydrophobic segments span the membrane, while three segments, Vb, VIII and IX, were predicted to have a periplasmic location (Van Geest, M., and Lolkema, J. S. (1996) J. Biol. Chem. 271, 25582-25589). A topology study of C-terminally truncated CitS molecules in dog pancreas microsomes revealed that the protein traverses the endoplasmic reticulum membrane 11 times. In agreement with the PhoA fusion data, segment Vb was predicted to have a periplasmic location, but, in contrast, segments VIII and IX were found to be membrane-spanning (Van Geest, M., Nilsson, I., von Heijne, G., and Lolkema, J. S. (1999) J. Biol. Chem. 274, 2816-2823). In the present study, using site-directed Cys labeling, the topology of segments VIII and IX in the full-length CitS protein was determined in the E. coli membrane. Engineered cysteine residues in the loop between the two segments were accessible to a membrane-impermeable thiol reagent exclusively from the cytoplasmic side of the membrane, demonstrating that transmembrane segments (TMSs) VIII and IX are both membrane-spanning. It follows that the folding of CitS in the E. coli and endoplasmic reticulum membrane is the same. Cysteine accessibility studies of CitS-PhoA fusion molecules demonstrated that in the E. coli membrane segment VIII is exported to the periplasm in the absence of the C-terminal CitS sequences, thus explaining why the PhoA fusions do not correctly predict the topology. An engineered cysteine residue downstream of TMS VIII moved from a periplasmic to a cytoplasmic location when the fusion protein containing TMSs I-VIII was extended with segment IX. Thus, downstream segment IX is both essential and sufficient for the insertion of segment VIII of CitS in the E. coli membrane.     

Improved membrane protein topology prediction by domain assignments

Topology predictions for integral membrane proteins can be substantially improved if parts of the protein can be constrained to a given in/out location relative to the membrane using experimental data or other information. Here, we have identified a set of 367 domains in the SMART database that, when found in soluble proteins, have compartment-specific localization of a kind relevant for membrane protein topology prediction. Using these domains as prediction constraints, we are able to provide high-quality topology models for 11% of the membrane proteins extracted from 38 eukaryotic genomes. Two-thirds of these proteins are single spanning, a group of proteins for which current topology prediction methods perform particularly poorly.     

Consensus predictions of membrane protein topology

We have explored the possibility that consensus predictions of membrane protein topology might provide a means to estimate the reliability of a predicted topology. Using five current topology prediction methods and a test set of 60 Escherichia coli inner membrane proteins with experimentally determined topologies, we find that prediction performance varies strongly with the number of methods that agree, and that the topology of nearly half of all E. coli inner membrane proteins can be predicted with high reliability (>90% correct predictions) by a simple majority-vote approach.     

Distant downstream sequence determinants can control N-tail translocation during protein insertion into the endoplasmic reticulum membrane

We have studied the membrane insertion of ProW, an Escherichia coli inner membrane protein with seven transmembrane segments and a large periplasmic N-terminal tail, into endoplasmic reticulum (ER)-derived dog pancreas microsomes. Strikingly, significant levels of N-tail translocation is seen only when a minimum of four of the transmembrane segments are present; for constructs with fewer transmembrane segments, the N-tail remains mostly nontranslocated and the majority of the molecules adopt an "inverted" topology where normally nontranslocated parts are translocated and vice versa. N-tail translocation can also be promoted by shortening of the N-tail and by the addition of positively charged residues immediately downstream of the first trasnmembrane segment. We conclude that as many as four consecutive transmembrane segments may be collectively involved in determining membrane protein topology in the ER and that the effects of downstream sequence determinants may vary depending on the size and charge of the N-tail. We also provide evidence to suggest that the ProW N-tail is translocated across the ER membrane in a C-to-N-terminal direction.     

Improving the accuracy of transmembrane protein topology prediction using evolutionary information

MOTIVATION: Many important biological processes such as cell signaling, transport of membrane-impermeable molecules, cell-cell communication, cell recognition and cell adhesion are mediated by membrane proteins. Unfortunately, as these proteins are not water soluble, it is extremely hard to experimentally determine their structure. Therefore, improved methods for predicting the structure of these proteins are vital in biological research. In order to improve transmembrane topology prediction, we evaluate the combined use of both integrated signal peptide prediction and evolutionary information in a single algorithm. RESULTS: A new method (MEMSAT3) for predicting transmembrane protein topology from sequence profiles is described and benchmarked with full cross-validation on a standard data set of 184 transmembrane proteins. The method is found to predict both the correct topology and the locations of transmembrane segments for 80% of the test set. This compares with accuracies of 62-72% for other popular methods on the same benchmark. By using a second neural network specifically to discriminate transmembrane from globular proteins, a very low overall false positive rate (0.5%) can also be achieved in detecting transmembrane proteins. AVAILABILITY: An implementation of the described method is available both as a web server (http://www.psipred.net) and as downloadable source code from http://bioinf.cs.ucl.ac.uk/memsat. Both the server and source code files are free to non-commercial users. Benchmark and training data are also available from http://bioinf.cs.ucl.ac.uk/memsat.