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.     

Biochemical evidence for the presence of mixed membrane topologies of the severe acute respiratory syndrome coronavirus envelope protein expressed in mammalian cells

Coronavirus envelope (E) protein is a small integral membrane protein with multi-functions in virion assembly, morphogenesis and virus-host interaction. Different coronavirus E proteins share striking similarities in biochemical properties and biological functions, but seem to adopt distinct membrane topology. In this report, we study the membrane topology of the SARS-CoV E protein by immunofluorescent staining of cells differentially permeabilized with detergents and proteinase K protection assay. It was revealed that both the N- and C-termini of the SARS-CoV E protein are exposed to the cytoplasmic side of the membranes (N(cyto)C(cyto)). In contrast, parallel experiments showed that the E protein from infectious bronchitis virus (IBV) spanned the membranes once, with the N-terminus exposed luminally and the C-terminus exposed cytoplasmically (N(exo(lum)-)C(cyto)). Intriguingly, a minor proportion of the SARS-CoV E protein was found to be modified by N-linked glycosylation on Asn 66 and inserted into the membranes once with the C-terminus exposed to the luminal side. The presence of two distinct membrane topologies of the SARS-CoV E protein may provide a useful clue to the pathogenesis of SARS-CoV.     

Reliability measures for membrane protein topology prediction algorithms

We have developed reliability scores for five widely used membrane protein topology prediction methods, and have applied them both on a test set of 92 bacterial plasma membrane proteins with experimentally determined topologies and on all predicted helix bundle membrane proteins in three fully sequenced genomes: Escherichia coli, Saccharomyces cerevisiae and Caenorhabditis elegans. We show that the reliability scores work well for the TMHMM and MEMSAT methods, and that they allow the probability that the predicted topology is correct to be estimated for any protein. We further show that the available test set is biased towards high-scoring proteins when compared to the genome-wide data sets, and provide estimates for the expected prediction accuracy of TMHMM across the three genomes. Finally, we show that the performance of TMHMM is considerably better when limited experimental information (such as the in/out location of a protein's C terminus) is available, and estimate that at least ten percentage points in overall accuracy in whole-genome predictions can be gained in this way.     

Topology-informed strategies for the overexpression and purification of membrane proteins

Membrane proteins represent a significant fraction of all genomes and play key roles in many aspects of biology, but their structural analysis has been hampered by difficulties in large-scale production and crystallisation. To overcome the first of these hurdles, we present here a systematic approach for expression and affinity-tagging which takes into account transmembrane topology. Using a set of bacterial transporters with known topologies, we tested the efficacy of a panel of conventional and Gateway? recombinational cloning vectors designed for protein expression under the control of the tac promoter, and for the addition of differing N- and C-terminal affinity tags. For transporters in which both termini are cytoplasmic, C-terminal oligohistidine tagging by recombinational cloning typically yielded functional protein at levels equivalent to or greater than those achieved by conventional cloning. In contrast, it was not effective for examples of the substantial minority of proteins that have one or both termini located on the periplasmic side of the membrane, possibly because of impairment of membrane insertion by the tag and/or att-site-encoded sequences. However, fusion either of an oligohistidine tag to cytoplasmic (but not periplasmic) termini, or of a Strep-tag II peptide to periplasmic termini using conventional cloning vectors did not interfere with membrane insertion, enabling high-level expression of such proteins. In conjunction with use of a C-terminal Lumio? fluorescence tag, which we found to be compatible with both periplasmic and cytoplasmic locations, these findings offer a system for strategic planning of construct design for high throughput expression of membrane proteins for structural genomics projects.     

Scanning N-glycosylation mutagenesis of membrane proteins

N-Glycosylation of eukaryotic membrane proteins is a co-translational event that occurs in the lumen of the endoplasmic reticulum (ER). This process is catalyzed by a membrane-associated oligosaccharyl transferase (OST) complex that transfers a preformed oligosaccharide (Glc₃Man₉GlcNAc₂-) to an asparagine (Asn) side-chain acceptor located within the sequon (-Asn-X-Ser/Thr-). Scanning N-glycosylation mutagenesis experiments, where novel acceptor sites are introduced at unique sites within membrane proteins, have shown that the acceptor sites must be located a minimum distance (12–14 amino acids) away from the luminal membrane surface of the ER in order to be efficiently N-glycosylated. Scanning N-glycosylation mutagenesis can therefore be used to determine membrane protein topology and it can also serve as a molecular ruler to define the ends of transmembrane (TM) segments. Furthermore, since N-glycosylation is a co-translational event, N-glycosylation mutagenesis can be used to identify folding intermediates in membrane proteins that may expose segments to the ER lumen transiently during biosynthesis.     

Mode of membrane insertion and sequence of a 32-amino acid peptide stretch of the penicillin-binding protein 4 of Enterococcus hirae

Analysis of water-soluble derivatives of the Enterococcus hirae 75-kDa membrane-bound penicillin-binding protein 4 (PBP4) has yielded the amino acid sequence of a 32-amino acid polypeptide stretch. This peptide is similar to peptide segments known to occur in the N-terminal domain of high-Mr PBPs of class B. The E. hirae PBP4 probably belongs to the same class. It is anchored in the membrane at the N-terminus of the polypeptide chain.     

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

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

Predicting secondary structures of membrane proteins with neural networks

Back-propagation, feed-forward neural networks are used to predict the secondary structures of membrane proteins whose structures are known to atomic resolution. These networks are trained on globular proteins and can predict globular protein structures having no homology to those of the training set with correlation coefficients (C) of 0.45, 0.32 and 0.43 for a-helix, -strand and random coil structures, respectively. When tested on membrane proteins, neural networks trained on globular proteins do, on average, correctly predict (Q_i) 62%, 38% and 69% of the residues in the -helix, -strand and random coil structures. These scores rank higher than those obtained with the currently used statistical methods and are comparable to those obtained with the joint approaches tested so far on membrane proteins. The lower success score for -strand as compared to the other structures suggests that the sample of -strand patterns contained in the training set is less representative than those of a-helix and random coil. Our analysis, which includes the effects of the network parameters and of the structural composition of the training set on the prediction, shows that regular patterns of secondary structures can be successfully extrapolated from globular to membrane proteins.Correspondence to: R. Casadio     

Topological analysis of DctQ, the small integral membrane protein of the C4-dicarboxylate TRAP transporter of Rhodobacter capsulatus

Tripartite ATP-independent periplasmic ('TRAP') transporters are a novel group of bacterial and archaeal secondary solute uptake systems which possess a periplasmic binding protein, but which are unrelated to ATP-binding cassette (ABC) systems. In addition to the binding protein, TRAP transporters contain two integral membrane proteins or domains, one of which is 40-50 kDa with 12 predicted transmembrane (TM) helices, thought to be the solute import protein, while the other is 20-30 kDa and of unknown function. Using a series of plasmid-encoded beta-lactamase fusions, we have determined the topology of DctQ, the smaller integral membrane protein from the high-affinity C4-dicarboxylate transporter of Rhodobacter capsulatus, which to date is the most extensively characterised TRAP transporter. DctQ was predicted by several topology prediction programmes to have four TM helices with the N- and C-termini located in the cytoplasm. The levels of ampicillin resistance conferred by the fusions when expressed in Escherichia coli were found to correlate with this predicted topology. The data have provided a topological model which can be used to test hypotheses concerning the function of the different regions of DctQ and which can be applied to other members of the DctQ family.     

SimFold energy function for de novo protein structure prediction: consensus with Rosetta

Predicting protein tertiary structures by in silico folding is still very difficult for proteins that have new folds. Here, we developed a coarse-grained energy function, SimFold, for de novo structure prediction, performed a benchmark test of prediction with fragment assembly simulations for 38 test proteins, and proposed consensus prediction with Rosetta. The SimFold energy consists of many terms that take into account solvent-induced effects on the basis of physicochemical consideration. In the benchmark test, SimFold succeeded in predicting native structures within 6.5 A for 12 of 38 proteins; this success rate was the same as that by the publicly available version of Rosetta (ab initio version 1.2) run with default parameters. We investigated which energy terms in SimFold contribute to structure prediction performance, finding that the hydrophobic interaction is the most crucial for the prediction, whereas other sequence-specific terms have weak but positive roles. In the benchmark, well-predicted proteins by SimFold and by Rosetta were not the same for 5 of 12 proteins, which led us to introduce consensus prediction. With combined decoys, we succeeded in prediction for 16 proteins, four more than SimFold or Rosetta separately. For each of 38 proteins, structural ensembles generated by SimFold and by Rosetta were qualitatively compared by mapping sampled structural space onto two dimensions. For proteins of which one of the two methods succeeded and the other failed in prediction, the former had a less scattered ensemble located around the native. For proteins of which both methods succeeded in prediction, often two ensembles were mixed up.     

Analysis of structure-function relationships in Escherichia coli K12 outer membrane porins with the aid of ompC-phoE and phoE-ompC hybrid genes

Summary To study structure-function relationships in the outer membrane pore proteins OmpC and PhoE of Escherichia coli K12, we have constructed a series of phoE-ompC hybrid genes in which DNA encoding part of one protein is replaced by the homologous part of the other gene. The hybrid gene products were incorporated normally into the outer membrane, allowing their functional characterization. Combined with previous studies, the present results permit the identification of regions involved in determining functions and properties in which the native PhoE and OmpC proteins differ, such as pore characteristics, receptor activity for phages and binding of monoclonal antibodies. Most of these properties were found to be determined by multiple regions clearly separated in the primary structure. The combined phage and antibody binding data have demonstrated that at least five distinct regions in PhoE and OmpC are exposed at the cell surface. The locations of these regions are in agreement with a previously proposed model for porin topology.