Transmembrane proteins in the Protein Data Bank: identification and classification

MOTIVATION: Integral membrane proteins play important roles in living cells. Although these proteins are estimated to constitute 25% of proteins at a genomic scale, the Protein Data Bank (PDB) contains only a few hundred membrane proteins due to the difficulties with experimental techniques. The presence of transmembrane proteins in the structure data bank, however, is quite invisible, as the annotation of these entries is rather poor. Even if a protein is identified as a transmembrane one, the possible location of the lipid bilayer is not indicated in the PDB because these proteins are crystallized without their natural lipid bilayer, and currently no method is publicly available to detect the possible membrane plane using the atomic coordinates of membrane proteins. RESULTS: Here, we present a new geometrical approach to distinguish between transmembrane and globular proteins using structural information only and to locate the most likely position of the lipid bilayer. An automated algorithm (TMDET) is given to determine the membrane planes relative to the position of atomic coordinates, together with a discrimination function which is able to separate transmembrane and globular proteins even in cases of low resolution or incomplete structures such as fragments or parts of large multi chain complexes. This method can be used for the proper annotation of protein structures containing transmembrane segments and paves the way to an up-to-date database containing the structure of all known transmembrane proteins and fragments (PDB_TM) which can be automatically updated. The algorithm is equally important for the purpose of constructing databases purely of globular proteins.     

Localization of organelle proteins by isotope tagging (LOPIT)

We describe a proteomics method for determining the subcellular localization of membrane proteins. Organelles are partially separated using centrifugation through self-generating density gradients. Proteins from each organelle co-fractionate and therefore exhibit similar distributions in the gradient. Protein distributions can be determined through a series of pair-wise comparisons of gradient fractions, using cleavable ICAT to enable relative quantitation of protein levels by MS. The localization of novel proteins is determined using multivariate data analysis techniques to match their distributions to those of proteins that are known to reside in specific organelles. Using this approach, we have simultaneously demonstrated the localization of membrane proteins in both the endoplasmic reticulum and the Golgi apparatus in Arabidopsis. Localization of organelle proteins by isotope tagging is a new tool for high-throughput protein localization, which is applicable to a wide range of research areas such as the study of organelle function and protein trafficking.     

Recognition of alpha-helix transmembrane domains with an amphipathy scale generated by molecular dynamics using only the primary sequence of proteins

We recently developed an amphipathy scale, elaborated from molecular dynamics data that can be used for the identification of hydrophobic or hydrophilic regions in proteins. This amphipathy scale reflects side chain/water molecule interaction energies. We have now used this amphipathy scale to find candidates for transmembrane segments, by examining a large sample of membrane proteins with alpha-helix segments. The candidates were selected based on an amphipathy coefficient value range and the minimum number of residues in a segment. We compared our results with the transmembrane segments previously identified in the PDB_TM database by the TMDET algorithm. We expected that the hydrophobic segments would be identified using only the primary structures of the proteins and the amphipathy scale. However, some of these hydrophobic segments may pertain to hydrophobic pockets not included in transmembrane regions. We found that our amphipathy scale could identify alpha-helix transmembrane regions with a probability of success of 76% when all segments were included and 90% when all membrane proteins were included.     

Fuzzy KNN for predicting membrane protein types from pseudo-amino acid composition

Cell membranes are vitally important to the life of a cell. Although the basic structure of biological membrane is provided by the lipid bilayer, membrane proteins perform most of the specific functions. Membrane proteins are putatively classified into five different types. Identification of their types is currently an important topic in bioinformatics and proteomics. In this paper, based on the concept of representing protein samples in terms of their pseudo-amino acid composition (Chou, K.C., 2001. Prediction of protein cellular attributes using pseudo amino acid composition. Proteins: Struct. Funct. Genet. 43, 246-255), the fuzzy K-nearest neighbors (KNN) algorithm has been introduced to predict membrane protein types, and high success rates were observed. It is anticipated that, the current approach, which is based on a branch of fuzzy mathematics and represents a new strategy, may play an important complementary role to the existing methods in this area. The novel approach may also have notable impact on prediction of the other attributes, such as protein structural class, protein subcellular localization, and enzyme family class, among many others.     

Tolerance to changes in membrane lipid composition as a selected trait of membrane proteins

Membrane lipid composition can vary dramatically across the three domains of life and even within single organisms. Here we review evidence that the lipid-exposed surfaces of membrane proteins have generally evolved to maintain correct structure and function in the face of major changes in lipid composition. Such tolerance has allowed evolution to extensively remodel membrane lipid compositions during the emergence of new species without having to extensively remodel the associated membrane proteins. The tolerance of membrane proteins also permits single-cell organisms to vary their membrane lipid composition in response to their changing environments and allows dynamic and organelle-specific variations in the lipid compositions of eukaryotic cells. Membrane protein structural biology has greatly benefited from this seemingly intrinsic property of membrane proteins: the majority of structures determined to date have been characterized under model membrane conditions that little resemble those of native membranes. Nevertheless, with a few notable exceptions, most experimentally determined membrane protein structures appear, to a good approximation, to faithfully report on native structure.     

Visualization of HRas Domains in the Plasma Membrane of Fibroblasts

The plasma membrane is a highly complex, organized structure where the lateral organization of signaling proteins is tightly regulated. In the case of Ras proteins, it has been suggested that the differential activity of the various isoforms is due to protein localization in separate membrane compartments. To date, direct visualization of such compartmentalization has been achieved only by electron microscopy on membrane sheets. Here, we combine photoactivated light microscopy with quantitative statistical analysis to visualize protein distribution in intact cells. In particular, we focus on the localization of HRas and its minimal anchoring domain, CAAX. We demonstrate the existence of a complex partitioning behavior, where small domains coexist with larger ones. The protein content in these domains varied from two molecules to tens of molecules. We found that 40% of CAAX and 60% of HRas were localized in domains. Subsequently, we were able to manipulate protein distributions by inducing coalescence of supposedly cholesterol-enriched domains. Clustering resulted in an increase of the localized fraction by 15%.     

Characterization of the preprotein and amino acid transporter gene family in Arabidopsis

Seventeen loci encode proteins of the preprotein and amino acid transporter family in Arabidopsis (Arabidopsis thaliana). Some of these genes have arisen from recent duplications and are not in annotated duplicated regions of the Arabidopsis genome. In comparison to a number of other eukaryotic organisms, this family of proteins has greatly expanded in plants, with 24 loci in rice (Oryza sativa). Most of the Arabidopsis and rice genes are orthologous, indicating expansion of this family before monocot and dicot divergence. In vitro protein uptake assays, in vivo green fluorescent protein tagging, and immunological analyses of selected proteins determined either mitochondrial or plastidic localization for 10 and six proteins, respectively. The protein encoded by At5g24650 is targeted to both mitochondria and chloroplasts and, to our knowledge, is the first membrane protein reported to be targeted to mitochondria and chloroplasts. Three genes encoded translocase of the inner mitochondrial membrane (TIM)17-like proteins, three TIM23-like proteins, and three outer envelope protein16-like proteins in Arabidopsis. The identity of Arabidopsis TIM22-like proteins is most likely a protein encoded by At3g10110/At1g18320, based on phylogenetic analysis, subcellular localization, and complementation of a yeast (Saccharomyces cerevisiae) mutant and coexpression analysis. The lack of a preprotein and amino acid transporter domain in some proteins, localization in mitochondria, plastids, or both, variation in gene structure, and the differences in expression profiles indicate that the function of this family has diverged in plants beyond roles in protein translocation.     

A novel cholesterol-producing Pichia pastoris strain is an ideal host for functional expression of human Na,K-ATPase α3β1 isoform

The heterologous expression of mammalian membrane proteins in lower eukaryotes is often hampered by aberrant protein localization, structure, and function, leading to enhanced degradation and, thus, low expression levels. Substantial quantities of functional membrane proteins are necessary to elucidate their structure–function relationships. Na,K-ATPases are integral, human membrane proteins that specifically interact with cholesterol and phospholipids, ensuring protein stability and enhancing ion transport activity. In this study, we present a Pichia pastoris strain which was engineered in its sterol pathway towards the synthesis of cholesterol instead of ergosterol to foster the functional expression of human membrane proteins. Western blot analyses revealed that cholesterol-producing yeast formed enhanced and stable levels of human Na,K-ATPase α3β1 isoform. ATPase activity assays suggested that this Na,K-ATPase isoform was functionally expressed in the plasma membrane. Moreover, [³H]-ouabain cell surface-binding studies underscored that the Na,K-ATPase was present in high numbers at the cell surface, surpassing reported expression strains severalfold. This provides evidence that the humanized sterol composition positively influenced Na,K-ATPase α3β1 stability, activity, and localization to the yeast plasma membrane. Prospectively, cholesterol-producing yeast will have high potential for functional expression of many mammalian membrane proteins.     

Estimation of membrane proteins in the human proteome

Genomics and proteomics have added valuable information to our knowledgebase of the human biological system including the discovery of therapeutic targets and disease biomarkers. However, molecular profiling studies commonly result in the identification of novel proteins of unknown localization. A class of proteins of special interest is membrane proteins, in particular plasma membrane proteins. Despite their biological and medical significance, the 3-dimensional structures of less than 1% of plasma membrane proteins have been determined. In order to aid in identification of membrane proteins, a number of computational methods have been developed. These tools operate by predicting the presence of transmembrane segments. Here, we utilized five topology prediction methods (TMHMM, SOSUI, waveTM, HMMTOP, and TopPred II) in order to estimate the ratio of integral membrane proteins in the human proteome. These methods employ different algorithms and include a newly-developed method (waveTM) that has yet to be tested on a large proteome database. Since these tools are prone for error mainly as a result of falsely predicting signal peptides as transmembrane segments, we have utilized an additional method, SignalP. Based on our analyses, the ratio of human proteins with transmembrane segments is estimated to fall between 15% and 39% with a consensus of 13%. Agreement among the programs is reduced further when both a positive identification of a membrane protein and the number of transmembrane segments per protein are considered. Such a broad range of prediction depends on the selectivity of the individual method in predicting integral membrane proteins. These methods can play a critical role in determining protein structure and, hence, identifying suitable drug targets in humans.     

Subcellular localization of a sporulation membrane protein is achieved through a network of interactions along and across the septum

During the process of spore formation in Bacillus subtilis many membrane proteins localize to the sporulation septum where they play key roles in morphogenesis and cell-cell signalling. However, the mechanism by which these proteins are anchored at this site is not understood. In this report we have defined the localization requirements for the mother-cell membrane protein SpoIVFA, which anchors a signalling complex in the septal membrane on the mother cell side. We have identified five proteins (SpoIID, SpoIIP, SpoIIM, BofA and SpoIIIAH) synthesized in the mother cell under the control of sigma(E) and one protein (SpoIIQ) synthesized in the forespore under the control of sigma(F) that are all required for the proper localization of SpoIVFA. Surprisingly, these proteins appear to have complementary and overlapping anchoring roles suggesting that SpoIVFA is localized in the septal membrane through a web of protein interactions. Furthermore, we demonstrate a direct biochemical interaction between the extracellular domains of two of the proteins required to anchor SpoIVFA: the forespore protein SpoIIQ and the mother-cell protein SpoIIIAH. This result supports the idea that the web of interactions that anchors SpoIVFA is itself held in the septal membrane through a zipper-like interaction across the sporulation septum. Importantly, our results suggest that a second mechanism independent of forespore proteins participates in anchoring SpoIVFA. Finally, we show that the dynamic localization of SpoIIQ in the forespore is impaired in the absence of SpoIVFA but not SpoIIIAH. Thus, a complex web of interactions among mother cell and forespore proteins is responsible for static and dynamic protein localization in both compartments of the sporangium. We envision that this proposed network is involved in anchoring other sporulation proteins in the septum and that protein networks with overlapping anchoring capacity is a feature of protein localization in all bacteria.     

MetaLocGramN: A meta-predictor of protein subcellular localization for Gram-negative bacteria

Subcellular localization is a key functional characteristic of proteins. It is determined by signals encoded in the protein sequence. The experimental determination of subcellular localization is laborious. Thus, a number of computational methods have been developed to predict the protein location from sequence. However predictions made by different methods often disagree with each other and it is not always clear which algorithm performs best for the given cellular compartment. We benchmarked primary subcellular localization predictors for proteins from Gram-negative bacteria, PSORTb3, PSLpred, CELLO, and SOSUI-GramN, on a common dataset that included 1056 proteins. We found that PSORTb3 performs best on the average, but is outperformed by other methods in predictions of extracellular proteins. This motivated us to develop a meta-predictor, which combines the primary methods by using the logistic regression models, to take advantage of their combined strengths, and to eliminate their individual weaknesses. MetaLocGramN runs the primary methods, and based on their output classifies protein sequences into one of five major localizations of the Gram-negative bacterial cell: cytoplasm, plasma membrane, periplasm, outer membrane, and extracellular space. MetaLocGramN achieves the average Matthews correlation coefficient of 0.806, i.e. 12% better than the best individual primary method. MetaLocGramN is a meta-predictor specialized in predicting subcellular localization for proteins from Gram-negative bacteria. According to our benchmark, it performs better than all other tools run independently. MetaLocGramN is a web and SOAP server available for free use by all academic users at the URL http://iimcb.genesilico.pl/MetaLocGramN. This article is part of a Special Issue entitled: Computational Methods for Protein Interaction and Structural Prediction.     

Cellular locations of Pseudomonas syringae pv. syringae HrcC and HrcJ proteins, required for harpin secretion via the type III pathway

The complete hrp-hrc-hrmA cluster of Pseudomonas syringae pv. syringae 61 encodes 28 polypeptides. A saprophytic bacterium carrying this cluster is capable of secreting HrpZ-a harpin encoded by hrpZ-in an hrp-dependent manner, which suggests that this cluster contains sufficient components to assemble functional type III secretion machinery. Sequence data show that HrcJ and HrcC are putative outer membrane proteins, and nonpolar mutagenesis demonstrates they are all required for HrpZ secretion. In this study, we investigated the cellular localization of the HrcC and HrcJ proteins by Triton solubilization, sucrose-gradient isopycnic centrifugation, and immunogold labeling of the bacterial cell surface. Our results indicate that HrcC is indeed an outer membrane protein and that HrcJ is located between both membranes. Their membrane localization suggests that they might be involved in the formation of a supramolecular structure for protein secretion.     

A folding pathway-dependent score to recognize membrane proteins

While various approaches exist to study protein localization, it is still a challenge to predict where proteins localize. Here, we consider a mechanistic viewpoint for membrane localization. Taking into account the steps for the folding pathway of α-helical membrane proteins and relating biophysical parameters to each of these steps, we create a score capable of predicting the propensity for membrane localization and call it FP(3)mem. This score is driven from the principal component analysis (PCA) of the biophysical parameters related to membrane localization. FP(3)mem allows us to rationalize the colocalization of a number of channel proteins with the Cav1.2 channel by their fewer propensities for membrane localization.     

Identification of new Golgi complex specific proteins by direct organelle proteomic analysis

The Golgi complex is in the crossroad of the endocytic and secretory pathways. Its function is to post-translationally modify and sort proteins and lipids, and regulate the membrane balance in the cell. To understand the structure-function relationship of the Golgi complex the Golgi proteome has to be identified first. We have used a direct organelle proteomic analysis to identify new Golgi complex proteins. Enriched stacked Golgi membrane fractions from rat livers were isolated, and the proteins from these membranes were subsequently digested into peptides. The peptides were fractionated by cation-exchange chromatography followed by protein identification by automated capillary-LC/ESI-MS/MS analysis and database searches. Two different search programs, ProID and MASCOT were used. This resulted in a total of 1125 protein identifications in two experiments. In addition to the known Golgi resident proteins, a significant number of unknown proteins were identified. Some of these were further characterized in silico using different programs to provide insight into their structure, intracellular localization and biological functions. The Golgi localization of two of these newly identified proteins was also confirmed by indirect immunofluorescence.     

Protein–protein association and cellular localization of four essential gene products encoded by tellurite resistance-conferring cluster “ter” from pathogenic Escherichia coli

Gene cluster “ter” conferring high tellurite resistance has been identified in various pathogenic bacteria including Escherichia coli O157:H7. However, the precise mechanism as well as the molecular function of the respective gene products is unclear. Here we describe protein–protein association and localization analyses of four essential Ter proteins encoded by minimal resistance-conferring fragment (terBCDE) by means of recombinant expression. By using a two-plasmid complementation system we show that the overproduced single Ter proteins are not able to mediate tellurite resistance, but all Ter members play an irreplaceable role within the cluster. We identified several types of homotypic and heterotypic protein–protein associations among the Ter proteins by in vitro and in vivo pull-down assays and determined their cellular localization by cytosol/membrane fractionation. Our results strongly suggest that Ter proteins function involves their mutual association, which probably happens at the interface of the inner plasma membrane and the cytosol.     

Influence of the extracellular domain size on the dynamic behavior of membrane proteins

The dynamic behavior of plasma membrane proteins mediates various cellular processes such as cellular motility, communication, and signaling. It is widely accepted that the dynamics of the membrane proteins is determined either by the interactions of the transmembrane domain with the surrounding lipids or by the interactions of the intracellular domain with cytosolic components such as cortical actin. Although initiation of different cellular signaling events at the plasma membrane has been attributed to the extracellular domain (ECD) properties recently, the impact of ECDs on the dynamic behavior of membrane proteins is rather unexplored. Here, we investigate how ECD properties influence protein dynamics in the lipid bilayer by reconstituting ECDs of different sizes or glycosylation in model membrane systems and analyzing ECD-driven protein sorting in lipid domains as well as protein mobility. Our data show that increasing the ECD mass or glycosylation leads to a decrease in ordered domain partitioning and diffusivity. Our data reconcile different mechanisms proposed for the initiation of cellular signaling by linking the ECD size of membrane proteins with their localization and diffusion dynamics in the plasma membrane.     

Assembly of a polytopic membrane protein structure from the solution structures of overlapping peptide fragments of bacteriorhodopsin.

Three-dimensional structures of only a handful of membrane proteins have been solved, in contrast to the thousands of structures of water-soluble proteins. Difficulties in crystallization have inhibited the determination of the three-dimensional structure of membrane proteins by x-ray crystallography and have spotlighted the critical need for alternative approaches to membrane protein structure. A new approach to the three-dimensional structure of membrane proteins has been developed and tested on the integral membrane protein, bacteriorhodopsin, the crystal structure of which had previously been determined. An overlapping series of 13 peptides, spanning the entire sequence of bacteriorhodopsin, was synthesized, and the structures of these peptides were determined by NMR in dimethylsulfoxide solution. These structures were assembled into a three-dimensional construct by superimposing the overlapping sequences at the ends of each peptide. Onto this construct were written all the distance and angle constraints obtained from the individual solution structures along with a limited number of experimental inter-helical distance constraints, and the construct was subjected to simulated annealing. A three-dimensional structure, determined exclusively by the experimental constraints, emerged that was similar to the crystal structure of this protein. This result suggests an alternative approach to the acquisition of structural information for membrane proteins consisting of helical bundles.     

Three-dimensional reconstruction of dynamin in the constricted state

Members of the dynamin family of GTPases have unique structural properties that might reveal a general mechanochemical basis for membrane constriction. Receptor-mediated endocytosis, caveolae internalization and certain trafficking events in the Golgi all require dynamin for vesiculation. The dynamin-related protein Drp1 (Dlp1) has been implicated in mitochondria fission and a plant dynamin-like protein phragmoplastin is involved in the vesicular events leading to cell wall formation. A common theme among these proteins is their ability to self-assemble into spirals and their localization to areas of membrane fission. Here we present the first three-dimensional structure of dynamin at a resolution of approximately 20 A, determined from cryo-electron micrographs of tubular crystals in the constricted state. The map reveals a T-shaped dimer consisting of three prominent densities: leg, stalk and head. The structure suggests that the dense stalk and head regions rearrange when GTP is added, a rearrangement that generates a force on the underlying lipid bilayer and thereby leads to membrane constriction. These results indicate that dynamin is a force-generating 'contrictase'.     

Targeting Rab GTPases to distinct membrane compartments

Rab GTPases are key to membrane-trafficking events in eukaryotic cells, and human cells contain more than 60 Rab proteins that are localized to distinct compartments. The recent determination of the structure of a monoprenylated Rab GTPase bound to GDP-dissociation inhibitor provides new molecular details that are relevant to models of Rab delivery. The further discovery of an integral membrane protein that can dissociate prenylated Rab proteins from GDP-dissociation inhibitor gives new insights into the mechanisms of Rab localization.