Parallel topology of genetically fused EmrE homodimers

EmrE is a small H+-coupled multidrug transporter in Escherichia coli. Claims have been made for an antiparallel topology of this homodimeric protein. However, our own biochemical studies performed with detergent-solubilized purified protein support a parallel topology of the protomers. We developed an alternative approach to constrain the relative topology of the protomers within the dimer so that their activity can be assayed also in vivo before biochemical handling. Tandem EmrE was built with two identical monomers genetically fused tail to head (C-terminus of the first to N-terminus of the second monomer) with hydrophilic linkers of varying length. All the constructs conferred resistance to ethidium by actively removing it from the cytoplasm. The purified proteins bound substrate and transported methyl viologen into proteoliposomes by a proton-dependent mechanism. A tandem where one of the essential glutamates was replaced with glutamine transported only monovalent substrates and displayed a modified stoichiometry. The results support a parallel topology of the protomers in the functional dimer. The implications regarding insertion and evolution of membrane proteins are discussed.     

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

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

Application of split-green fluorescent protein for topology mapping membrane proteins in Escherichia coli

A topology map of a membrane protein defines the location of transmembrane helices and the orientation of soluble domains relative to the membrane. In the absence of a high-resolution structure, a topology map is an essential guide for studying structure-function relationships. Although these maps can be predicted directly from amino acid sequence, the predictions are more accurate if combined with experimental data, which are usually obtained by fusing a reporter protein to the C-terminus of the protein. However, as reporter proteins are large, they cannot be used to report on the cytoplasmic/periplasmic location of the N-terminus of a protein. Here, we show that the bimolecular split-green fluorescent protein complementation system can overcome this limitation and can be used to determine the location of both the N- and C-termini of inner membrane proteins in Escherichia coli.     

Quasi-symmetry in the cryo-EM structure of EmrE provides the key to modeling its transmembrane domain

Small multidrug resistance (SMR) transporters contribute to bacterial resistance by coupling the efflux of a wide range of toxic aromatic cations, some of which are commonly used as antibiotics and antiseptics, to proton influx. EmrE is a prototypical small multidrug resistance transporter comprising four transmembrane segments (M1-M4) that forms dimers. It was suggested recently that EmrE molecules in the dimer have different topologies, i.e. monomers have opposite orientations with respect to the membrane plane. A 3-D structure of EmrE acquired by electron cryo-microscopy (cryo-EM) at 7.5 Angstroms resolution in the membrane plane showed that parts of the structure are related by quasi-symmetry. We used this symmetry relationship, combined with sequence conservation data, to assign the transmembrane segments in EmrE to the densities seen in the cryo-EM structure. A C alpha model of the transmembrane region was constructed by considering the evolutionary conservation pattern of each helix. The model is validated by much of the biochemical data on EmrE with most of the positions that were identified as affecting substrate translocation being located around the substrate-binding cavity. A suggested mechanism for proton-coupled substrate translocation in small multidrug resistance antiporters provides a mechanistic rationale to the experimentally observed inverted topology.     

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.     

Remote homology detection of integral membrane proteins using conserved sequence features

Compared with globular proteins, transmembrane proteins are surrounded by a more intricate environment and, consequently, amino acid composition varies between the different compartments. Existing algorithms for homology detection are generally developed with globular proteins in mind and may not be optimal to detect distant homology between transmembrane proteins. Here, we introduce a new profile-profile based alignment method for remote homology detection of transmembrane proteins in a hidden Markov model framework that takes advantage of the sequence constraints placed by the hydrophobic interior of the membrane. We expect that, for distant membrane protein homologs, even if the sequences have diverged too far to be recognized, the hydrophobicity pattern and the transmembrane topology are better conserved. By using this information in parallel with sequence information, we show that both sensitivity and specificity can be substantially improved for remote homology detection in two independent test sets. In addition, we show that alignment quality can be improved for the most distant homologs in a public dataset of membrane protein structures. Applying the method to the Pfam domain database, we are able to suggest new putative evolutionary relationships for a few relatively uncharacterized protein domain families, of which several are confirmed by other methods. The method is called Searcher for Homology Relationships of Integral Membrane Proteins (SHRIMP) and is available for download at http://www.sbc.su.se/shrimp/.     

Predicting the topology of transmembrane helical proteins using mean burial propensity and a hidden-Markov-model-based method

Helices in membrane spanning regions are more tightly packed than the helices in soluble proteins. Thus, we introduce a method that uses a simple scale of burial propensity and a new algorithm to predict transmembrane helical (TMH) segments and a positive-inside rule to predict amino-terminal orientation. The method (the topology predictor of transmembrane helical proteins using mean burial propensity [THUMBUP]) correctly predicted the topology of 55 of 73 proteins (or 75%) with known three-dimensional structures (the 3D helix database). This level of accuracy can be reached by MEMSAT 1.8 (a 200-parameter model-recognition method) and a new HMM-based method (a 111-parameter hidden Markov model, UMDHMM(TMHP)) if they were retrained with the 73-protein database. Thus, a method based on a physiochemical property can provide topology prediction as accurate as those methods based on more complicated statistical models and learning algorithms for the proteins with accurately known structures. Commonly used HMM-based methods and MEMSAT 1.8 were trained with a combination of the partial 3D helix database and a 1D helix database of TMH proteins in which topology information were obtained by gene fusion and other experimental techniques. These methods provide a significantly poorer prediction for the topology of TMH proteins in the 3D helix database. This suggests that the 1D helix database, because of its inaccuracy, should be avoided as either a training or testing database. A Web server of THUMBUP and UMDHMM(TMHP) is established for academic users at http://www.smbs.buffalo.edu/phys_bio/service.htm. The 3D helix database is also available from the same Web site.     

Bax forms multispanning monomers that oligomerize to permeabilize membranes during apoptosis

Bax promotes cell death by permeabilizing mitochondrial outer membranes by an unresolved mechanism. However, in cells lacking the gene c-myc, membrane permeabilization by Bax is blocked by changes in the mitochondria that prevent Bax oligomerization. Drug-treated c-myc null cells and cells expressing Myc were used to map the topology of Bax in membranes prior to and after mitochondrial permeabilization. Chemical labeling of single cysteine mutants of Bax using a membrane bilayer impermeant cysteine-specific modifying agent revealed that Bax inserted both the 'pore domain' (helices alpha5-alpha6), and the tail-anchor (helix alpha9) into membranes prior to oligomerization and membrane permeabilization. Additional topology changes for Bax were not required in Myc-expressing cells to promote oligomerization and cytochrome c release. Our results suggest that unlike most pore-forming proteins, Bax membrane permeabilization results from oligomerization of transmembrane monomers rather than concerted insertion of the pore domains of a preformed oligomer.     

Structural analysis of synthetic peptide fragments from EmrE,a multidrug resistance protein,in a membrane-mimetic environment

EmrE, a multidrug resistance protein from Escherichia coli, renders the bacterium resistant to a variety of cytotoxic drugs by active translocation out of the cell. The 110-residue sequence of EmrE limits the number of structural possibilities that can be envisioned for this membrane protein. Four helix bundle models have been considered [Yerushalmi, H., Lebendiker, M., and Schuldiner, S. (1996) J. Biol. Chem. 271, 31044-31048]. The validity of EmrE structural models has been probed experimentally by investigations on overlapping peptides (ranging in length from 19 to 27 residues), derived from the sequence of EmrE. The choice of peptides was made to provide sequences of two complete, predicted transmembrane helices (peptides H1 and H3) and two helix-loop-helix motifs (peptides A and B). Peptide (B) also corresponds to a putative hairpin in a speculative beta-barrel model, with the "Pro-Thr-Gly" segment forming a turn. Structure determination in SDS micelles using NMR indicates peptide H1 to be predominantly helical, with helix boundaries in the micellar environment corroborating predicted helical limits. Peptide A adopts a helix-loop-helix structure in SDS micelles, and peptide B was also largely helical in micellar environments. An analogue peptide, C, in which the central "Pro-Thr-Gly" was replaced by "(D)Pro-Gly" displays local turn conformation at the (D)Pro-Gly segment, but neither a continuous helical stretch nor beta-hairpin formation was observed. This study implies that the constraints of membrane and micellar environments largely direct the structure of transmembrane peptides and proteins and study of judiciously selected peptide fragments can prove useful in the structural elucidation of membrane proteins.     

Why is the biological hydrophobicity scale more accurate than earlier experimental hydrophobicity scales?

The recognition of transmembrane helices by the translocon is primarily guided by the average hydrophobicity of the potential transmembrane helix. However, the exact hydrophobicity of each amino acid can be identified in several different ways. The free energy of transfer for amino acid analogues between a hydrophobic media, for example, octanol and water can be measured or obtained from simulations, the hydrophobicity can also be estimated by statistical properties from known transmembrane segments and finally the contribution of each amino acid type for the probability of translocon recognition has recently been measured directly. Although these scales correlate quite well, there are clear differences between them and it is not well understood which scale represents neither the biology best nor what the differences are. Here, we try to provide some answers to this by studying the ability of different scales to recognize transmembrane helices and predict the topology of transmembrane proteins. From this analysis it is clear that the biological hydrophobicity scale as well scales created from statistical analysis of membrane helices perform better than earlier experimental scales that are mainly based on measurements of amino acid analogs and not directly on transmembrane helix recognition. Using these results we identified the properties of the scales that perform better than other scales. We find, for instance, that the better performing scales consider proline more hydrophilic. This shows that transmembrane recognition is not only governed by pure hydrophobicity but also by the helix preferences for amino acids, as proline is a strong helix breaker. Proteins 2014; 82:2190–2198. © 2014 Wiley Periodicals, Inc.     

Topology and immersion depth of an integral membrane protein by paramagnetic rates from dissolved oxygen

In studies of membrane proteins, knowledge of protein topology can provide useful insight into both structure and function. In this work, we present a solution NMR method for the measurement the tilt angle and average immersion depth of alpha helices in membrane proteins, from analysis of the paramagnetic relaxation rate enhancements arising from dissolved oxygen. No modification to the micelle or protein is necessary, and the topology of both transmembrane and amphipathic helices are readily determined. We apply this method to the measure the topology of a monomeric mutant of phospholamban (AFA-PLN), a 52-residue membrane protein containing both an amphipathic and a transmembrane alpha helix. In dodecylphosphocholine micelles, the amphipathic helix of AFA-PLN was found to have a tilt angle of 87° ± 1° and an average immersion depth of 13.2 ?. The transmembrane helix was found to have an average immersion depth of 5.4 ?, indicating residues 41 and 42 are closest to the micelle centre. The resolution of paramagnetic relaxation rate enhancements from dissolved oxygen compares favourably to those from Ni (II), a hydrophilic paramagnetic species.     

Membrane topology of transmembrane proteins: determinants and experimental tools

Membrane topology refers to the two-dimensional structural information of a membrane protein that indicates the number of transmembrane (TM) segments and the orientation of soluble domains relative to the plane of the membrane. Since membrane proteins are co-translationally translocated across and inserted into the membrane, the TM segments orient themselves properly in an early stage of membrane protein biogenesis. Each membrane protein must contain some topogenic signals, but the translocation components and the membrane environment also influence the membrane topology of proteins. We discuss the factors that affect membrane protein orientation and have listed available experimental tools that can be used in determining membrane protein topology.     

Marginally hydrophobic transmembrane α‐helices shaping membrane protein folding

Cells have developed an incredible machinery to facilitate the insertion of membrane proteins into the membrane. While we have a fairly good understanding of the mechanism and determinants of membrane integration, more data is needed to understand the insertion of membrane proteins with more complex insertion and folding pathways. This review will focus on marginally hydrophobic transmembrane helices and their influence on membrane protein folding. These weakly hydrophobic transmembrane segments are by themselves not recognized by the translocon and therefore rely on local sequence context for membrane integration. How can such segments reside within the membrane? We will discuss this in the light of features found in the protein itself as well as the environment it resides in. Several characteristics in proteins have been described to influence the insertion of marginally hydrophobic helices. Additionally, the influence of biological membranes is significant. To begin with, the actual cost for having polar groups within the membrane may not be as high as expected; the presence of proteins in the membrane as well as characteristics of some amino acids may enable a transmembrane helix to harbor a charged residue. The lipid environment has also been shown to directly influence the topology as well as membrane boundaries of transmembrane helices—implying a dynamic relationship between membrane proteins and their environment.     

Helical integrity and microsolvation of transmembrane domains from Flaviviridae envelope glycoproteins

Charged and polar amino acids in the transmembrane domains of integral membrane proteins can be crucial for protein function and also promote helix-helix association or protein oligomerization. Yet, our current understanding is still limited on how these hydrophilic amino acids are efficiently translocated from the Sec61/SecY translocon into the cell membrane during the biogenesis of membrane proteins. In hepatitis C virus, the putative transmembrane segments of envelope glycoproteins E1 and E2 were suggested to heterodimerize via a Lys-Asp ion-pair in the host endoplasmic reticulum. Therefore in this work, we carried out molecular dynamic simulations in explicit lipid bilayer and solvent environment to explore the stability of all possible bridging ion-pairs using the model of H-segment helix dimers. We observed that, frequently, several water molecules penetrated from the interface into the membrane core to stabilize the charged and polar pairs. The hydration time and amount of water molecules in the membrane core depended on the position of the charged residues as well as on the type of ion-pairs. Similar microsolvation events were observed in simulations of the putative E1-E2 transmembrane helix dimers. Simulations of helix monomers from other members of the Flaviviridae family suggest that these systems show similar behaviors. Thus this study illustrates the important contribution of water microsolvation to overcome the unfavorable energetic cost of burying charged and polar amino acids in membrane lipid bilayers. Also, it emphasizes the novel role of bridging charged or polar interactions stabilized by water molecules in the hydrophobic lipid bilayer core that has an important biological function for helix dimerization in several envelope glycoproteins from the family of Flaviviridae viruses.     

How important are transmembrane helices of bitopic membrane proteins?

The topology of a bitopic membrane protein consists of a single transmembrane helix connecting two extra-membranous domains. As opposed to helices from polytopic proteins, the transmembrane helices of bitopic proteins were initially considered as merely hydrophobic anchors, while more recent studies have begun to shed light on their role in the protein's function. Herein the overall importance of transmembrane helices from bitopic membrane proteins was analyzed using a relative conservation analysis. Interestingly, the transmembrane domains of bitopic proteins are on average, significantly more conserved than the remainder of the protein, even when taking into account their smaller amino acid repertoire. Analysis of highly conserved transmembrane domains did not reveal any unifying consensus, pointing to a great diversity in their conservation patterns. However, Fourier power spectrum analysis was able to show that regardless of the conservation motif, in most sequences a significant conservation moment was observed, in that one side of the helix was conserved while the other was not. Taken together, it may be possible to conclude that a significant proportion of transmembrane helices from bitopic membrane proteins participate in specific interactions, in a variety of modes in the plane of the lipid bilayer.     

How important are transmembrane helices of bitopic membrane proteins?

The topology of a bitopic membrane protein consists of a single transmembrane helix connecting two extra-membranous domains. As opposed to helices from polytopic proteins, the transmembrane helices of bitopic proteins were initially considered as merely hydrophobic anchors, while more recent studies have begun to shed light on their role in the protein's function. Herein the overall importance of transmembrane helices from bitopic membrane proteins was analyzed using a relative conservation analysis. Interestingly, the transmembrane domains of bitopic proteins are on average, significantly more conserved than the remainder of the protein, even when taking into account their smaller amino acid repertoire. Analysis of highly conserved transmembrane domains did not reveal any unifying consensus, pointing to a great diversity in their conservation patterns. However, Fourier power spectrum analysis was able to show that regardless of the conservation motif, in most sequences a significant conservation moment was observed, in that one side of the helix was conserved while the other was not. Taken together, it may be possible to conclude that a significant proportion of transmembrane helices from bitopic membrane proteins participate in specific interactions, in a variety of modes in the plane of the lipid bilayer.     

Tsc10p and FVT1: topologically distinct short-chain reductases required for long-chain base synthesis in yeast and mammals

In yeast, Tsc10p catalyzes reduction of 3-ketosphinganine to dihydrosphingosine. In mammals, it has been proposed that this reaction is catalyzed by FVT1, which despite limited homology and a different predicted topology, can replace Tsc10p in yeast. Silencing of FVT1 revealed a direct correlation between FVT1 levels and reductase activity, showing that FVT1 is the principal 3-ketosphinganine reductase in mammalian cells. Localization and topology studies identified an N-terminal membrane-spanning domain in FVT1 (absent in Tsc10p) oriented to place it in the endoplasmic reticulum (ER) lumen. In contrast, protease digestion studies showed that the N terminus of Tsc10p is cytoplasmic. Fusion of the N-terminal domain of FVT1 to green fluorescent protein directed the fusion protein to the ER, demonstrating that it is sufficient for targeting. Although both proteins have two predicted transmembrane domains C-terminal to a cytoplasmic catalytic domain, neither had an identifiable lumenal loop. Nevertheless, both Tsc10p and the residual fragment of FVT1 produced by removal of the N-terminal domain with factor Xa protease behave as integral membrane proteins. In addition to their topological differences, mutation of conserved catalytic residues had different effects on the activities of the two enzymes. Thus, while FVT1 can replace Tsc10p in yeast, there are substantial differences between the two enzymes that may be important for regulation of sphingolipid biosynthesis in higher eukaryotes.     

Dislocation by the m-AAA Protease Increases the Threshold Hydrophobicity for Retention of Transmembrane Helices in the Inner Membrane of Yeast Mitochondria

Sorting of mitochondrial inner membrane proteins is a complex process in which translocons and proteases function in a concerted way. Many inner membrane proteins insert into the membrane via the TIM23 translocon, and some are then further acted upon by the mitochondrial m-AAA protease, a molecular motor capable of dislocating proteins from the inner membrane. This raises the possibility that the threshold hydrophobicity for the retention of transmembrane segments in the inner membrane is different depending on whether they belong to membrane proteins that are m-AAA protease substrates or not. Here, using model transmembrane segments engineered into m-AAA protease-dependent proteins, we show that the threshold hydrophobicity for membrane retention measured in yeast cells in the absence of a functional m-AAA protease is markedly lower than that measured in its presence. Whether a given hydrophobic segment in a mitochondrial inner membrane protein will ultimately form a transmembrane helix may therefore depend on whether or not it will be exposed to the pulling force exerted by the m-AAA protease during biogenesis.     

Rapid exploration of the folding topology of helical membrane proteins using paramagnetic perturbation

An understanding of the folding states of α-helical membrane proteins in detergent systems is important for functional and structural studies of these proteins. Here, we present a rapid and simple method for identification of the folding topology and assembly of transmembrane helices using paramagnetic perturbation in nuclear magnetic resonance spectroscopy. By monitoring the perturbation of signals from glycine residues located at specific sites, the folding topology and the assembly of transmembrane helices of membrane proteins were easily identified without time-consuming backbone assignment. This method is validated with Mistic (membrane-integrating sequence for translation of integral membrane protein constructs) of known structure as a reference protein. The folding topologies of two bacterial histidine kinase membrane proteins (SCO3062 and YbdK) were investigated by this method in dodecyl phosphocholine (DPC) micelles. Combing with analytical ultracentrifugation, we identified that the transmembrane domain of YbdK is present as a parallel dimer in DPC micelle. In contrast, the interaction of transmembrane domain of SCO3062 is not maintained in DPC micelle due to disruption of native structure of the periplasmic domain by DPC micelle.     

Helix packing moments reveal diversity and conservation in membrane protein structure

Helical membrane proteins are more tightly packed and the packing interactions are more diverse than those found in helical soluble proteins. Based on a linear correlation between amino acid packing values and interhelical propensity, we propose the concept of a helix packing moment to predict the orientation of helices in helical membrane proteins and membrane protein complexes. We show that the helix packing moment correlates with the helix interfaces of helix dimers of single pass membrane proteins of known structure. Helix packing moments are also shown to help identify the packing interfaces in membrane proteins with multiple transmembrane helices, where a single helix can have multiple contact surfaces. Analyses are described on class A G protein-coupled receptors (GPCRs) with seven transmembrane helices. We show that the helix packing moments are conserved across the class A family of GPCRs and correspond to key structural contacts in rhodopsin. These contacts are distinct from the highly conserved signature motifs of GPCRs and have not previously been recognized. The specific amino acid types involved in these contacts, however, are not necessarily conserved between subfamilies of GPCRs, indicating that the same protein architecture can be supported by a diverse set of interactions. In GPCRs, as well as membrane channels and transporters, amino acid residues with small side-chains (Gly, Ala, Ser, Cys) allow tight helix packing by mediating strong van der Waals interactions between helices. Closely packed helices, in turn, facilitate interhelical hydrogen bonding of both weakly polar (Ser, Thr, Cys) and strongly polar (Asn, Gln, Glu, Asp, His, Arg, Lys) amino acid residues. We propose the use of the helix packing moment as a complementary tool to the helical hydrophobic moment in the analysis of transmembrane sequences.