Consensus predictions of membrane protein topology

Defining Membrane Protein Topology Using pho-lac Reporter Fusions

Experimental determination of membrane protein topology can be achieved using various techniques. Here, we present the pho-lac dual reporter system, a simple, convenient, and reliable tool to analyze the topology of membrane protein in vivo. The system is based on the use of two topological markers with complementary properties: The Escherichia coli β-galactosidase, LacZ, which is active in the cytoplasm, and the E. coli alkaline phosphatase, PhoA, which is active in the bacterial periplasm. Specifically, in this pho-lac gene system, the reporter molecule is a chimera composed of the mature PhoA that is in-frame with the β-galactosidase α-peptide (LacZα). Hence, when targeted to the periplasm, the PhoA-LacZα dual reporter displays high alkaline phosphatase activity but no β-galactosidase activity. Conversely, when located in the cytoplasm, PhoA-LacZα has no phosphatase activity but exhibits high β-galactosidase activity in E. coli cells expressing the α-fragment of LacZ, LacZα (via the α-complementation phenomenon). The dual nature of the PhoA-LacZα reporter allows a simple way to normalize both enzymatic activities to obtain readily interpretable information about the subcellular location of the fusion site between the membrane protein under study and the reporter. In addition, the PhoA-LacZα reporter permits the utilization of dual indicator agar plates to easily discriminate between colonies bearing cytoplasmic fusions, periplasmic fusions, or out-of-frame fusions. In total, the phoA-lacZα fusion reporter approach is a direct and rather inexpensive method to characterize the topology of membrane proteins in vivo.

Site-Directed Alkylation Detected by In-Gel Fluorescence (SDAF) to Determine the Topology Map and Probe the Solvent Accessibility of Membrane Proteins

The topology of helix-bundle membrane proteins provides low-resolution structural information with regard to the number and orientation of membrane-spanning helices, as well as the sidedness of intra/extra-cellular domains. In the past decades, several strategies have been developed to experimentally determine the topology of membrane proteins. However, generally, these methods are labour-intensive, time-consuming and difficult to implement for quantitative analysis. Here, we report a novel approach, site-directed alkylation detected by in-gel fluorescence (SDAF), which monitors the fluorescent band shift caused by alkylation of the EGFP-fused target membrane protein bearing one single introduced cysteine. In-gel fluorescence provides a unique readout of target membrane proteins with EGFP fusion from non-purified samples, revealing a distinct 5 kDa shift on SDS-PAGE gel due to conjugation with mPEG-MAL-5K. Using the structurally characterised bile acid transporter ASBT_NM as an example, we demonstrate that SDAF generates a topology map consistent with the crystal structure. The efficiency of mPEG-MAL-5K modification at each introduced cysteine can easily be quantified and analysed, providing a useful tool for probing the solvent accessibility at a specific position of the target membrane protein.

Computational Approaches for Revealing the Structure of Membrane Transporters: Case Study on Bilitranslocase

The structural and functional details of transmembrane proteins are vastly underexplored, mostly due to experimental difficulties regarding their solubility and stability. Currently, the majority of transmembrane protein structures are still unknown and this present a huge experimental and computational challenge. Nowadays, thanks to X-ray crystallography or NMR spectroscopy over 3000 structures of membrane proteins have been solved, among them only a few hundred unique ones. Due to the vast biological and pharmaceutical interest in the elucidation of the structure and the functional mechanisms of transmembrane proteins, several computational methods have been developed to overcome the experimental gap. If combined with experimental data the computational information enables rapid, low cost and successful predictions of the molecular structure of unsolved proteins. The reliability of the predictions depends on the availability and accuracy of experimental data associated with structural information. In this review, the following methods are proposed for in silico structure elucidation: sequence-dependent predictions of transmembrane regions, predictions of transmembrane helix–helix interactions, helix arrangements in membrane models, and testing their stability with molecular dynamics simulations. We also demonstrate the usage of the computational methods listed above by proposing a model for the molecular structure of the transmembrane protein bilitranslocase. Bilitranslocase is bilirubin membrane transporter, which shares similar tissue distribution and functional properties with some of the members of the Organic Anion Transporter family and is the only member classified in the Bilirubin Transporter Family. Regarding its unique properties, bilitranslocase is a potentially interesting drug target.

Defining Membrane Protein Topology Using pho-lac Reporter Fusions

Experimental determination of membrane protein topology can be achieved using various techniques. Here we present the pho-lac dual reporter system, a simple, convenient, and reliable tool to analyze the topology of membrane proteins in vivo. The system is based on the use of two topological markers with complementary properties, the Escherichia coli β-galactosidase LacZ, which is active in the cytoplasm, and the E. coli alkaline phosphatase PhoA, which is active in the bacterial periplasm. Specifically, in this pho-lac gene system, the reporter molecule is a chimera composed of the mature PhoA that is in frame with the β-galactosidase α-peptide, LacZα. Hence, when targeted to the periplasm, the PhoA-LacZα dual reporter displays high alkaline phosphatase activity but no β-galactosidase activity. Conversely, when located in the cytoplasm, PhoA-LacZα has no phosphatase activity but exhibits high β-galactosidase activity in E. coli cells expressing the ω fragment of LacZ, LacZω (via the α-complementation phenomenon). The dual nature of the PhoA-LacZα reporter allows a simple way to normalize both enzymatic activities to obtain readily interpretable information about the subcellular location of the fusion site between the membrane protein under study and the reporter. In addition, the PhoA-LacZα reporter permits utilization of dual-indicator agar plates to easily discriminate between colonies bearing cytoplasmic fusions, periplasmic fusions, or out-of-frame fusions. In total, the phoA-lacZα fusion reporter approach is a straightforward and rather inexpensive method of characterizing the topology of membrane proteins in vivo.

Enquiry into the Topology of Plasma Membrane-Localized PIN Auxin Transport Components

Auxin directs plant ontogenesis via differential accumulation within tissues depending largely on the activity of PIN proteins that mediate auxin efflux from cells and its directional cell-to-cell transport. Regardless of the developmental importance of PINs, the structure of these transporters is poorly characterized. Here, we present experimental data concerning protein topology of plasma membrane-localized PINs. Utilizing approaches based on pH-dependent quenching of fluorescent reporters combined with immunolocalization techniques, we mapped the membrane topology of PINs and further cross-validated our results using available topology modeling software. We delineated the topology of PIN1 with two transmembrane (TM) bundles of five α-helices linked by a large intracellular loop and a C-terminus positioned outside the cytoplasm. Using constraints derived from our experimental data, we also provide an updated position of helical regions generating a verisimilitude model of PIN1. Since the canonical long PINs show a high degree of conservation in TM domains and auxin transport capacity has been demonstrated for Arabidopsis representatives of this group, this empirically enhanced topological model of PIN1 will be an important starting point for further studies on PIN structure-function relationships. In addition, we have established protocols that can be used to probe the topology of other plasma membrane proteins in plants.

VirB6 and VirB10 from the Brucella type IV secretion system interact via the N-terminal periplasmic domain of VirB6

Type IV secretion systems are multi-protein complexes that transfer macromolecules across the cell envelope of bacteria. Identifying the sites of interaction between the twelve proteins (VirB1-VirB11 and VirD4) that form these complexes is key to understanding their assembly and function. We have here used phage display, bacterial two-hybrid and fluorescence-based interaction assays to identify an N-terminal domain of the inner membrane protein VirB6 as a site of interaction with the envelope-spanning VirB10 protein. Our results are consistent with the notion that VirB6 acts in concert with VirB10 as well as with VirB8 during secretion system assembly and function.

The human transmembrane proteome

Background

Transmembrane proteins have important roles in cells, as they are involved in energy production, signal transduction, cell-cell interaction, cell-cell communication and more. In human cells, they are frequently targets for pharmaceuticals; therefore, knowledge about their properties and structure is crucial. Topology of transmembrane proteins provide a low resolution structural information, which can be a starting point for either laboratory experiments or modelling their 3D structures.

Results

Here, we present a database of the human α-helical transmembrane proteome, including the predicted and/or experimentally established topology of each transmembrane protein, together with the reliability of the prediction. In order to distinguish transmembrane proteins in the proteome as well as for topology prediction, we used a newly developed consensus method (CCTOP) that incorporates recent state of the art methods, with tested accuracies on a novel human benchmark protein set. CCTOP utilizes all available structure and topology data as well as bioinformatical evidences for topology prediction in a probabilistic framework provided by the hidden Markov model. This method shows the highest accuracy (98.5 % for discrinimating between transmembrane and non-transmembrane proteins and 84 % for per protein topology prediction) among the dozen tested topology prediction methods. Analysis of the human proteome with the CCTOP indicates that it contains 4998 (26 %) transmembrane proteins. Besides predicting topology, reliability of the predictions is estimated as well, and it is demonstrated that the per protein prediction accuracies of more than 60 % of the predictions are over 98 % on the benchmark sets and most probably on the predicted human transmembrane proteome too.

Conclusions

Here, we present the most accurate prediction of the human transmembrane proteome together with the experimental topology data. These data, as well as various statistics about the human transmembrane proteins and their topologies can be downloaded from and can be visualized at the website of the human transmembrane proteome (http://htp.enzim.hu).

Reviewers

This article was reviewed by Dr. Sandor Pongor, Dr. Michael Galperin and Dr. Pascale Gaudet (nominated by Dr Michael Galperin).

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-015-0061-x) contains supplementary material, which is available to authorized users.

Conserved loop cysteines of vitamin K epoxide reductase complex subunit 1-like 1 (VKORC1L1) are involved in its active site regeneration

Vitamin K epoxide reductase complex subunit 1 (VKORC1) reduces vitamin K epoxide in the vitamin K cycle for post-translational modification of proteins that are involved in a variety of biological functions. However, the physiological function of VKORC1-like 1 (VKORC1L1), a paralogous enzyme sharing about 50% protein identity with VKORC1, is unknown. Here we determined the structural and functional differences of these two enzymes using fluorescence protease protection (FPP) assay and an in vivo cell-based activity assay. We show that in vivo VKORC1L1 reduces vitamin K epoxide to support vitamin K-dependent carboxylation as efficiently as does VKORC1. However, FPP assays show that unlike VKORC1, VKORC1L1 is a four-transmembrane domain protein with both its termini located in the cytoplasm. Moreover, the conserved loop cysteines, which are not required for VKORC1 activity, are essential for VKORC1L1's active site regeneration. Results from domain exchanges between VKORC1L1 and VKORC1 suggest that it is VKORC1L1's overall structure that uniquely allows for active site regeneration by the conserved loop cysteines. Intermediate disulfide trapping results confirmed an intra-molecular electron transfer pathway for VKORC1L1's active site reduction. Our results allow us to propose a concerted action of the four conserved cysteines of VKORC1L1 for active site regeneration; the second loop cysteine, Cys-58, attacks the active site disulfide, forming an intermediate disulfide with Cys-139; the first loop cysteine, Cys-50, attacks the intermediate disulfide resulting in active site reduction. The different membrane topologies and reaction mechanisms between VKORC1L1 and VKORC1 suggest that these two proteins might have different physiological functions.

Large-scale proteomic analysis of membrane proteins

Proteomic analysis of membrane proteins is a promising approach for the identification of novel drug targets and/or disease biomarkers. Despite notable technological developments, obstacles related to extraction and solublization of membrane proteins are encountered. A critical discussion of the different preparative methods of membrane proteins is offered in relation to downstream proteomic applications, mainly gel-based analyses and mass spectrometry. Frequently, unknown proteins are identified by high-throughput profiling of membrane proteins. In search for novel membrane proteins, analysis of protein sequences using computational tools is performed to predict the presence of transmembrane domains. This review also presents these bioinformatic tools with the human proteome as a case study. Along with technological innovations, advancements in the areas of sample preparation and computational prediction of membrane proteins will lead to exciting discoveries.

The D3 bacteriophage α-polymerase inhibitor (Iap) peptide disrupts O-antigen biosynthesis through mimicry of the chain length regulator Wzz in Pseudomonas aeruginosa

Lysogenic bacteriophage D3 causes seroconversion of Pseudomonas aeruginosa PAO1 from serotype O5 to O16 by inverting the linkage between O-specific antigen (OSA) repeat units from α to β. The OSA units are polymerized by Wzy to modal lengths regulated by Wzz1 and Wzz2. A key component of the D3 seroconversion machinery is the inhibitor of α-polymerase (Iap) peptide, which is able to solely suppress α-linked long-chain OSA production in P. aeruginosa PAO1. To establish the target specificity of Iap for Wzyα, changes in OSA phenotypes were examined via Western immunoblotting for wzz1 and wzz2 single-knockout strains, as well as a wzz1 wzz2 double knockout, following the expression of iap from a tuneable vector. Increased induction of Iap expression completely abrogated OSA production in the wzz1 wzz2 double mutant, while background levels of OSA production were still observed in either of the single mutants. Therefore, Iap inhibition of OSA biosynthesis was most effective in the absence of both Wzz proteins. Sequence alignment analyses revealed a high degree of similarity between Iap and the first transmembrane segment (TMS) of either Wzz1 or Wzz2. Various topology prediction analyses of the Iap sequence consistently predicted the presence of a single TMS, suggesting a propensity for Iap to insert itself into the inner membrane (IM). The compromised ability of Iap to abrogate Wzyα function in the presence of Wzz1 or Wzz2 provides compelling evidence that inhibition occurs after Wzyα inserts itself into the IM and is achieved through mimicry of the first TMS from the Wzz proteins of P. aeruginosa PAO1.

Recognition and binding of apocytochrome c to P. aeruginosa CcmI, a component of cytochrome c maturation machinery

The biogenesis of c-type cytochromes (Cytc) is a process that in Gram-negative bacteria demands the coordinated action of different periplasmic proteins (CcmA-I), whose specific roles are still being investigated. Activities of Ccm proteins span from the chaperoning of heme b in the periplasm to the specific reduction of oxidized apocytochrome (apoCyt) cysteine residues and to chaperoning and recognition of the unfolded apoCyt before covalent attachment of the heme to the cysteine thiols can occur. We present here the functional characterization of the periplasmic domain of CcmI from the pathogen Pseudomonas aeruginosa (Pa-CcmI*). Pa-CcmI* is composed of a TPR domain and a peculiar C-terminal domain. Pa-CcmI* fulfills both the ability to recognize and bind to P. aeruginosa apo-cytochrome c551 (Pa-apoCyt) and a chaperoning activity towards unfolded proteins, as it prevents citrate synthase aggregation in a concentration-dependent manner. Equilibrium and kinetic experiments with Pa-CcmI*, or its isolated domains, with peptides mimicking portions of Pa-apoCyt sequence allow us to quantify the molecular details of the interaction between Pa-apoCyt and Pa-CcmI*. Binding experiments show that the interaction occurs at the level of the TPR domain and that the recognition is mediated mainly by the C-terminal sequence of Pa-apoCyt. The affinity of Pa-CcmI* to full-length Pa-apoCyt or to its C-terminal sequence is in the range expected for a component of a multi-protein complex, whose task is to receive the apoCyt and to deliver it to other components of the apoCyt:heme b ligation protein machinery.

Protein Machineries Involved in the Attachment of Heme to Cytochrome c: Protein Structures and Molecular Mechanisms

Cytochromes c (Cyt c) are ubiquitous heme-containing proteins, mainly involved in electron transfer processes, whose structure and functions have been and still are intensely studied. Surprisingly, our understanding of the molecular mechanism whereby the heme group is covalently attached to the apoprotein (apoCyt) in the cell is still largely unknown. This posttranslational process, known as Cyt c biogenesis or Cyt c maturation, ensures the stereospecific formation of the thioether bonds between the heme vinyl groups and the cysteine thiols of the apoCyt heme binding motif. To accomplish this task, prokaryotic and eukaryotic cells have evolved distinctive protein machineries composed of different proteins. In this review, the structural and functional properties of the main maturation apparatuses found in gram-negative and gram-positive bacteria and in the mitochondria of eukaryotic cells will be presented, dissecting the Cyt c maturation process into three functional steps: (i) heme translocation and delivery, (ii) apoCyt thioreductive pathway, and (iii) apoCyt chaperoning and heme ligation. Moreover, current hypotheses and open questions about the molecular mechanisms of each of the three steps will be discussed, with special attention to System I, the maturation apparatus found in gram-negative bacteria.

Improving transmembrane protein consensus topology prediction using inter-helical interaction

Alpha helix transmembrane proteins (αTMPs) represent roughly 30% of all open reading frames (ORFs) in a typical genome and are involved in many critical biological processes. Due to the special physicochemical properties, it is hard to crystallize and obtain high resolution structures experimentally, thus, sequence-based topology prediction is highly desirable for the study of transmembrane proteins (TMPs), both in structure prediction and function prediction. Various model-based topology prediction methods have been developed, but the accuracy of those individual predictors remain poor due to the limitation of the methods or the features they used. Thus, the consensus topology prediction method becomes practical for high accuracy applications by combining the advances of the individual predictors. Here, based on the observation that inter-helical interactions are commonly found within the transmembrane helixes (TMHs) and strongly indicate the existence of them, we present a novel consensus topology prediction method for αTMPs, CNTOP, which incorporates four top leading individual topology predictors, and further improves the prediction accuracy by using the predicted inter-helical interactions. The method achieved 87% prediction accuracy based on a benchmark dataset and 78% accuracy based on a non-redundant dataset which is composed of polytopic αTMPs. Our method derives the highest topology accuracy than any other individual predictors and consensus predictors, at the same time, the TMHs are more accurately predicted in their length and locations, where both the false positives (FPs) and the false negatives (FNs) decreased dramatically. The CNTOP is available at: http://ccst.jlu.edu.cn/JCSB/cntop/CNTOP.html.

Phosphatidylglycerol::prolipoprotein diacylglyceryl transferase (Lgt) of Escherichia coli has seven transmembrane segments, and its essential residues are embedded in the membrane

Lgt of Escherichia coli catalyzes the transfer of an sn-1,2-diacylglyceryl group from phosphatidylglycerol to prolipoproteins. The enzyme is essential for growth, as demonstrated here by the analysis of an lgt depletion strain. Cell fractionation demonstrated that Lgt is an inner membrane protein. Its membrane topology was determined by fusing Lgt to β-galactosidase and alkaline phosphatase and by substituted cysteine accessibility method (SCAM) studies. The data show that Lgt is embedded in the membrane by seven transmembrane segments, that its N terminus faces the periplasm, and that its C terminus faces the cytoplasm. Highly conserved amino acids in Lgt of both Gram-negative and Gram-positive bacteria were identified. Lgt enzymes are characterized by a so-called Lgt signature motif in which four residues are invariant. Ten conserved residues were replaced with alanine, and the activity of these Lgt variants was analyzed by their ability to complement the lgt depletion strain. Residues Y26, N146, and G154 are absolutely required for Lgt function, and R143, E151, R239, and E243 are important. The results demonstrate that the majority of the essential residues of Lgt are located in the membrane and that the Lgt signature motif faces the periplasm.

RECOGNITION OF TRANSMEMBRANE SEGMENTS IN PROTEINS: REVIEW AND CONSISTENCY-BASED BENCHMARKING OF INTERNET SERVERS

Membrane proteins perform a number of crucial functions as transporters, receptors, and components of enzyme complexes. Identification of membrane proteins and prediction of their topology is thus an important part of genome annotation. We present here an overview of transmembrane segments in protein sequences, summarize data from large-scale genome studies, and report results of benchmarking of several popular internet servers.

Revised transmembrane orientation of the NADH:quinone oxidoreductase subunit NuoA

NuoA is a small membrane spanning subunit of respiratory chain NADH:quinone oxidoreductase (complex I). Unlike the other complex I core protein subunits, the NuoA protein has no known homologue in other enzyme systems. The transmembrane orientation of NuoA cannot be unambiguously predicted, due to the small size of the polypeptide and the varying distribution of charged amino acid residues in NuoA from different organisms. Novel analyses of NuoA from Escherichia coli complex I expressed as fusion proteins to cytochrome c and to alkaline phosphatase demonstrated that the c-terminal end of the polypeptide is localized in the bacterial cytoplasm, in contrast to what was previously reported for the homologous NQO7 subunit from Paracoccus denitrificans complex I.

Membrane topological structure of neutral system N/A amino acid transporter 4 (SNAT4) protein

Members of system N/A amino acid transporter (SNAT) family mediate transport of neutral amino acids, including l-alanine, l-glutamine, and l-histidine, across the plasma membrane and are involved in a variety of cellular functions. By using chemical labeling, glycosylation, immunofluorescence combined with molecular modeling approaches, we resolved the membrane topological structure of SNAT4, a transporter expressed predominantly in liver. To analyze the orientation using the chemical labeling and biotinylation approach, the "Cys-null" mutant of SNAT4 was first generated by mutating all five endogenous cysteine residues. Based on predicted topological structures, a single cysteine residue was introduced individually into all possible nontransmembrane domains of the Cys-null mutant. The cells expressing these mutants were labeled with N-biotinylaminoethyl methanethiosulfonate, a membrane-impermeable cysteine-directed reagent. We mapped the orientations of N- and C-terminal domains. There are three extracellular loop domains, and among them, the second loop domain is the largest that spans from amino acid residue ∼242 to ∼335. The orientation of this domain was further confirmed by the identification of two N-glycosylated residues, Asn-260 and Asn-264. Together, we showed that SNAT4 contains 10 transmembrane domains with extracellular N and C termini and a large N-glycosylated, extracellular loop domain. This is the first report concerning membrane topological structure of mammalian SNAT transporters, which will provide important implications for our understanding of structure-function of the members in this amino acid transporter family.

Genome-wide inventory of metal homeostasis-related gene products including a functional phytochelatin synthase in the hypogeous mycorrhizal fungus Tuber melanosporum

Ectomycorrhizal fungi are thought to enhance mineral nutrition of their host plants and to confer increased tolerance toward toxic metals. However, a global view of metal homeostasis-related genes and pathways in these organisms is still lacking. Building upon the genome sequence of Tuber melanosporum and on transcriptome analyses, we set out to systematically identify metal homeostasis-related genes in this plant-symbiotic ascomycete. Candidate gene products (101) were subdivided into three major functional classes: (i) metal transport (58); (ii) oxidative stress defence (32); (iii) metal detoxification (11). The latter class includes a small-size metallothionein (TmelMT) that was functionally validated in yeast, and phytochelatin synthase (TmelPCS), the first enzyme of this kind to be described in filamentous ascomycetes. Recombinant TmelPCS was shown to support GSH-dependent, metal-activated phytochelatin synthesis in vitro and to afford increased Cd/Cu tolerance to metal hypersensitive yeast strains. Metal transporters, especially those related to Cu and Zn trafficking, displayed the highest expression levels in mycorrhizae, suggesting extensive translocation of both metals to root cells as well as to fungal metalloenzymes (e.g., laccase) that are strongly upregulated in symbiotic hyphae.

Spin labeling of the Escherichia coli NADH ubiquinone oxidoreductase (complex I)

The proton-pumping NADH:ubiquinone oxidoreductase, the respiratory complex I, couples the transfer of electrons from NADH to ubiquinone with the translocation of protons across the membrane. Electron microscopy revealed the two-part structure of the complex with a peripheral arm involved in electron transfer and a membrane arm most likely involved in proton translocation. It was proposed that the quinone binding site is located at the joint of the two arms. Most likely, proton translocation in the membrane arm is enabled by the energy of the electron transfer reaction in the peripheral arm transmitted by conformational changes. For the detection of the conformational changes and the localization of the quinone binding site, we set up a combination of site-directed spin labeling and EPR spectroscopy. Cysteine residues were introduced to the surface of the Escherichia coli complex I. The spin label (1-oxyl-2,2,5,5-tetramethyl-Δ3-pyrroline-3-methyl)-methanethiosulfonate (MTSL) was exclusively bound to the engineered positions. Neither the mutation nor the labeling had an effect on the NADH:decyl-ubiquinone oxidoreductase activity. The characteristic signals of the spin label were detected by EPR spectroscopy, which did not change by reducing the preparation with NADH. A decyl-ubiquinone derivative with the spin label covalently attached to the alkyl chain was synthesized in order to localize the quinone binding site. The distance between a MTSL labeled complex I variant and the bound quinone was determined by continuous-wave (cw) EPR allowing an inference on the location of the quinone binding site. The distances between the labeled quinone and other complex I variants will be determined in future experiments to receive further geometry information by triangulation.

Membrane topology and identification of critical amino acid residues in the Wzx O-antigen translocase from Escherichia coli O157:H4

Wzx belongs to a family of membrane proteins involved in the translocation of isoprenoid lipid-linked glycans, which is loosely related to members of the major facilitator superfamily. Despite Wzx homologs performing a conserved function, it has been difficult to pinpoint specific motifs of functional significance in their amino acid sequences. Here, we elucidate the topology of the Escherichia coli O157 Wzx (Wzx(EcO157)) by a combination of bioinformatics and substituted cysteine scanning mutagenesis, as well as targeted deletion-fusions to green fluorescent protein and alkaline phosphatase. We conclude that Wzx(EcO157) consists of 12 transmembrane (TM) helices and six periplasmic and five cytosolic loops, with N and C termini facing the cytoplasm. Four TM helices (II, IV, X, and XI) contain polar residues (aspartic acid or lysine), and they may form part of a relatively hydrophilic core. Thirty-five amino acid replacements to alanine or serine were targeted to five native cysteines and most of the aspartic acid, arginine, and lysine residues. From these, only replacements of aspartic acid-85, aspartic acid-326, arginine-298, and lysine-419 resulted in a protein unable to support O-antigen production. Aspartic acid-85 and lysine-419 are located in TM helices II and XI, while arginine-298 and aspartic acid-326 are located in periplasmic and cytosolic loops 4, respectively. Further analysis revealed that the charge at these positions is required for Wzx function since conservative substitutions maintaining the same charge polarity resulted in a functional protein, whereas those reversing or eliminating polarity abolished function. We propose that the functional requirement of charged residues at both sides of the membrane and in two TM helices could be important to allow the passage of the Und-PP-linked saccharide substrate across the membrane.

Membrane topology of conserved components of the type III secretion system from the plant pathogen Xanthomonas campestris pv. vesicatoria

Type III secretion (T3S) systems play key roles in the assembly of flagella and the translocation of bacterial effector proteins into eukaryotic host cells. Eleven proteins which are conserved among Gram-negative plant and animal pathogenic bacteria have been proposed to build up the basal structure of the T3S system, which spans both inner and outer bacterial membranes. We studied six conserved proteins, termed Hrc, predicted to reside in the inner membrane of the plant pathogen Xanthomonas campestris pv. vesicatoria. The membrane topology of HrcD, HrcR, HrcS, HrcT, HrcU and HrcV was studied by translational fusions to a dual alkaline phosphatase–β-galactosidase reporter protein. Two proteins, HrcU and HrcV, were found to have the same membrane topology as the Yersinia homologues YscU and YscV. For HrcR, the membrane topology differed from the model for the homologue from Yersinia, YscR. For our data on three other protein families, exemplified by HrcD, HrcS and HrcT, we derived the first topology models. Our results provide what is believed to be the first complete model of the inner membrane topology of any bacterial T3S system and will aid in elucidating the architecture of T3S systems by ultrastructural analysis.

Prediction of the human membrane proteome

Membrane proteins are key molecules in the cell, and are important targets for pharmaceutical drugs. Few three-dimensional structures of membrane proteins have been obtained, which makes computational prediction of membrane proteins crucial for studies of these key molecules. Here, seven membrane protein topology prediction methods based on different underlying algorithms, such as hidden Markov models, neural networks and support vector machines, have been used for analysis of the protein sequences from the 21,416 annotated genes in the human genome. The number of genes coding for a protein with predicted alpha-helical transmembrane region(s) ranged from 5508 to 7651, depending on the method used. Based on a majority decision method, we estimate 5539 human genes to code for membrane proteins, corresponding to approximately 26% of the human protein-coding genes. The largest fraction of these proteins has only one predicted transmembrane region, but there are also many proteins with seven predicted transmembrane regions, including the G-protein coupled receptors. A visualization tool displaying the topologies suggested by the eight prediction methods, for all predicted membrane proteins, is available on the public Human Protein Atlas portal (www.proteinatlas.org).

Bioinformatic approaches for the structure and function of membrane proteins

Membrane proteins play important roles in the biology of the cell, including intercellular communication and molecular transport. Their well-established importance notwithstanding, the high-resolution structures of membrane proteins remain elusive due to difficulties in protein expression, purification and crystallization. Thus, accurate prediction of membrane protein topology can increase the understanding of membrane protein function. Here, we provide a brief review of the diverse computational methods for predicting membrane protein structure and function, including recent progress and essential bioinformatics tools. Our hope is that this review will be instructive to users studying membrane protein biology in their choice of appropriate bioinformatics methods.

MetaTM - a consensus method for transmembrane protein topology prediction

Background

Transmembrane (TM) proteins are proteins that span a biological membrane one or more times. As their 3-D structures are hard to determine, experiments focus on identifying their topology (i. e. which parts of the amino acid sequence are buried in the membrane and which are located on either side of the membrane), but only a few topologies are known. Consequently, various computational TM topology predictors have been developed, but their accuracies are far from perfect. The prediction quality can be improved by applying a consensus approach, which combines results of several predictors to yield a more reliable result.

Results

A novel TM consensus method, named MetaTM, is proposed in this work. MetaTM is based on support vector machine models and combines the results of six TM topology predictors and two signal peptide predictors. On a large data set comprising 1460 sequences of TM proteins with known topologies and 2362 globular protein sequences it correctly predicts 86.7% of all topologies.

Conclusion

Combining several TM predictors in a consensus prediction framework improves overall accuracy compared to any of the individual methods. Our proposed SVM-based system also has higher accuracy than a previous consensus predictor. MetaTM is made available both as downloadable source code and as DAS server at

TOPCONS: consensus prediction of membrane protein topology

TOPCONS (http://topcons.net/) is a web server for consensus prediction of membrane protein topology. The underlying algorithm combines an arbitrary number of topology predictions into one consensus prediction and quantifies the reliability of the prediction based on the level of agreement between the underlying methods, both on the protein level and on the level of individual TM regions. Benchmarking the method shows that overall performance levels match the best available topology prediction methods, and for sequences with high reliability scores, performance is increased by ∼10 percentage points. The web interface allows for constraining parts of the sequence to a known inside/outside location, and detailed results are displayed both graphically and in text format.

The adaptive evolution of the mammalian mitochondrial genome

Background

The mitochondria produce up to 95% of a eukaryotic cell's energy through oxidative phosphorylation. The proteins involved in this vital process are under high functional constraints. However, metabolic requirements vary across species, potentially modifying selective pressures. We evaluate the adaptive evolution of 12 protein-coding mitochondrial genes in 41 placental mammalian species by assessing amino acid sequence variation and exploring the functional implications of observed variation in secondary and tertiary protein structures.

Results

Wide variation in the properties of amino acids were observed at functionally important regions of cytochrome b in species with more-specialized metabolic requirements (such as adaptation to low energy diet or large body size, such as in elephant, dugong, sloth, and pangolin, and adaptation to unusual oxygen requirements, for example diving in cetaceans, flying in bats, and living at high altitudes in alpacas). Signatures of adaptive variation in the NADH dehydrogenase complex were restricted to the loop regions of the transmembrane units which likely function as protons pumps. Evidence of adaptive variation in the cytochrome c oxidase complex was observed mostly at the interface between the mitochondrial and nuclear-encoded subunits, perhaps evidence of co-evolution. The ATP8 subunit, which has an important role in the assembly of F₀, exhibited the highest signal of adaptive variation. ATP6, which has an essential role in rotor performance, showed a high adaptive variation in predicted loop areas.

Conclusion

Our study provides insight into the adaptive evolution of the mtDNA genome in mammals and its implications for the molecular mechanism of oxidative phosphorylation. We present a framework for future experimental characterization of the impact of specific mutations in the function, physiology, and interactions of the mtDNA encoded proteins involved in oxidative phosphorylation.

The Arabidopsis vacuolar malate channel is a member of the ALMT family

In plants, malate is a central metabolite and fulfills a large number of functions. Vacuolar malate may reach very high concentrations and fluctuate rapidly, whereas cytosolic malate is kept at a constant level allowing optimal metabolism. Recently, a vacuolar malate transporter (Arabidopsis thaliana tonoplast dicarboxylate transporter, AttDT) was identified that did not correspond to the well-characterized vacuolar malate channel. We therefore hypothesized that a member of the aluminum-activated malate transporter (ALMT) gene family could code for a vacuolar malate channel. Using GFP fusion constructs, we could show that AtALMT9 (A. thaliana ALMT9) is targeted to the vacuole. Promoter-GUS fusion constructs demonstrated that this gene is expressed in all organs, but is cell-type specific as GUS activity in leaves was detected nearly exclusively in mesophyll cells. Patch-clamp analysis of an Atalmt9 T-DNA insertion mutant exhibited strongly reduced vacuolar malate channel activity. In order to functionally characterize AtALMT9 as a malate channel, we heterologously expressed this gene in tobacco and in oocytes. Overexpression of AtALMT9-GFP in Nicotiana benthamiana leaves strongly enhanced the malate current densities across the mesophyll tonoplasts. Functional expression of AtALMT9 in Xenopus oocytes induced anion currents, which were clearly distinguishable from endogenous oocyte currents. Our results demonstrate that AtALMT9 is a vacuolar malate channel. Deletion mutants for AtALMT9 exhibit only slightly reduced malate content in mesophyll protoplasts and no visible phenotype, indicating that AttDT and the residual malate channel activity are sufficient to sustain the transport activity necessary to regulate the cytosolic malate homeostasis.

Functional and evolutionary implications of enhanced genomic analysis of rhomboid intramembrane proteases

Rhomboids are a recently discovered family of widely distributed intramembrane serine proteases. They have diverse biological functions, including the regulation of growth factor signaling, mitochondrial fusion, and parasite invasion. Despite their existence in all branches of life, the sequence identity between rhomboids is low. We have combined BLAST-based database mining with functional and structural data to generate a comprehensive genomic analysis of eukaryotic rhomboid-like proteins. We show that robust membrane topology models are necessary to classify active rhomboid proteases unambiguously, and we define rules for distinguishing predicted active proteases from the larger evolutionary group of rhomboid-like proteins. This leads to a revision of estimates of numbers of proteolytically active rhomboids. We identify three groups of eukaryotic rhomboid-like proteins: true active rhomboids, a tightly clustered group of novel inactive rhomboids that we name the iRhoms, and a small number of other inactive rhomboid-like proteins. The active proteases are themselves subdivided into secretase and PARL-type (mitochondrial) subfamilies; these have distinct transmembrane topologies. This enhanced genomic analysis leads to conclusions about rhomboid enzyme function. It suggests that a given rhomboid can only cleave a single orientation of substrate, and that both products of rhomboid catalyzed intramembrane cleavage can be released from the membrane. Our phylogeny predictions also have evolutionary implications: Despite the complex classification of rhomboids, our data suggest that a rhomboid-type intramembrane protease may have been present in the last eukaryotic common ancestor.

Membrane insertion and topology of the p7B movement protein of Melon Necrotic Spot Virus (MNSV)

Cell-to-cell movement of the Melon Necrotic Spot Virus (MNSV) is controlled by two small proteins working in trans, an RNA-binding protein (p7A) and an integral membrane protein (p7B) separated by an amber stop codon. p7B contains a single hydrophobic region. Membrane integration of this region was observed when inserted into model proteins in the presence of microsomal membranes. Furthermore, we explored the topology and targeting mechanisms of full-length p7B. Here we present evidence that p7B integrates in vitro into the ER membrane cotranslationally and with an Nt-cytoplasmic/Ct-luminal orientation. The observed topology was monitored in vivo by fusing GFP to the Ct of p7B, enabling the overexpression in Escherichia coli cultures. Finally, the topology of a putative p14 movement protein was established by replacing the amber stop codon located between p7A and p7B.

Membrane Protein Structure: Prediction versus Reality

Since high-resolution structural data are still scarce, different kinds of theoretical structure prediction algorithms are of major importance in membrane protein biochemistry. But how well do the current prediction methods perform? Which structural features can be predicted and which cannot? And what can we expect in the next few years?

Membrane protein prediction methods

We survey computational approaches that tackle membrane protein structure and function prediction. While describing the main ideas that have led to the development of the most relevant and novel methods, we also discuss pitfalls, provide practical hints and highlight the challenges that remain. The methods covered include: sequence alignment, motif search, functional residue identification, transmembrane segment and protein topology predictions, homology and ab initio modeling. In general, predictions of functional and structural features of membrane proteins are improving, although progress is hampered by the limited amount of high-resolution experimental information available. While predictions of transmembrane segments and protein topology rank among the most accurate methods in computational biology, more attention and effort will be required in the future to ameliorate database search, homology and ab initio modeling.

A phylogenomic analysis of the Actinomycetales mce operons

Background

The genome of Mycobacterium tuberculosis harbors four copies of a cluster of genes termed mce operons. Despite extensive research that has demonstrated the importance of these operons on infection outcome, their physiological function remains obscure. Expanding databases of complete microbial genome sequences facilitate a comparative genomic approach that can provide valuable insight into the role of uncharacterized proteins.

Results

The M. tuberculosis mce loci each include two yrbE and six mce genes, which have homology to ABC transporter permeases and substrate-binding proteins, respectively. Operons with an identical structure were identified in all Mycobacterium species examined, as well as in five other Actinomycetales genera. Some of the Actinomycetales mce operons include an mkl gene, which encodes an ATPase resembling those of ABC uptake transporters. The phylogenetic profile of Mkl orthologs exactly matched that of the Mce and YrbE proteins. Through topology and motif analyses of YrbE homologs, we identified a region within the penultimate cytoplasmic loop that may serve as the site of interaction with the putative cognate Mkl ATPase. Homologs of the exported proteins encoded adjacent to the M. tuberculosis mce operons were detected in a conserved chromosomal location downstream of the majority of Actinomycetales operons. Operons containing linked mkl, yrbE and mce genes, resembling the classic organization of an ABC importer, were found to be common in Gram-negative bacteria and appear to be associated with changes in properties of the cell surface.

Conclusion

Evidence presented suggests that the mce operons of Actinomycetales species and related operons in Gram-negative bacteria encode a subfamily of ABC uptake transporters with a possible role in remodeling the cell envelope.

The use of functional domains to improve transmembrane protein topology prediction

Transmembrane proteins affect vital cellular functions and pathogenesis, and are a focus of drug design. It is difficult to obtain diffraction quality crystals to study transmembrane protein structure. Computational tools for transmembrane protein topology prediction fill in the gap between the abundance of transmembrane proteins and the scarcity of known membrane protein structures. Their prediction accuracy is still inadequate: TMHMM, the current state-of-the-art method, has less than 52% accuracy in topology prediction on one set of transmembrane proteins of known topology. Based on the observation that there are functional domains that occur preferentially internal or external to the membrane, we have extended the model of TMHMM to incorporate functional domains, using a probabilistic approach originally developed for computational gene finding. Our extension is better than TMHMM in predicting the topology of transmembrane proteins. As prediction of functional domain improves, our system's prediction accuracy will likely improve as well.

Membrane-associated transcripts in Arabidopsis; their isolation and characterization by DNA microarray analysis and bioinformatics

Membrane-associated, integral membrane and secreted proteins are of key importance in many cellular processes. For most of the 28,952 predicted proteins in Arabidopsis, the actual subcellular localization has not been demonstrated experimentally. So far, their potential membrane-association has been deduced from algorithms that predict transmembrane domains and signal peptides. However, the comprehensiveness and accuracy of these algorithms is still limited. The majority of membrane-associated and secreted proteins is synthesized on membrane-bound polysomes. Therefore, the isolation and characterization of mRNA associated with membrane-bound polysomes offers an experimental tool for the genome-wide identification of these proteins. Here we describe an efficient method to isolate mRNA from membrane-bound polysomes and report on the validation of the method to enrich for transcripts encoding membrane-associated and secreted proteins. The sensitivity and reproducibility of the isolation method was investigated by DNA microarray analysis. Pearson correlations between transcript levels obtained from three replicate isolations showed that the method is highly reproducible. A significant enrichment for mRNAs encoding proteins containing predicted transmembrane domains and signal peptides was observed in the membrane-bound polysomal fraction. In this fraction, 301 transcripts were classified by gene ontologies as 'cellular component unknown', and potentially encode previously unrecognized secreted or membrane-associated proteins.

Algorithms for incorporating prior topological information in HMMs: application to transmembrane proteins

Background

Hidden Markov Models (HMMs) have been extensively used in computational molecular biology, for modelling protein and nucleic acid sequences. In many applications, such as transmembrane protein topology prediction, the incorporation of limited amount of information regarding the topology, arising from biochemical experiments, has been proved a very useful strategy that increased remarkably the performance of even the top-scoring methods. However, no clear and formal explanation of the algorithms that retains the probabilistic interpretation of the models has been presented so far in the literature.

Results

We present here, a simple method that allows incorporation of prior topological information concerning the sequences at hand, while at the same time the HMMs retain their full probabilistic interpretation in terms of conditional probabilities. We present modifications to the standard Forward and Backward algorithms of HMMs and we also show explicitly, how reliable predictions may arise by these modifications, using all the algorithms currently available for decoding HMMs. A similar procedure may be used in the training procedure, aiming at optimizing the labels of the HMM's classes, especially in cases such as transmembrane proteins where the labels of the membrane-spanning segments are inherently misplaced. We present an application of this approach developing a method to predict the transmembrane regions of alpha-helical membrane proteins, trained on crystallographically solved data. We show that this method compares well against already established algorithms presented in the literature, and it is extremely useful in practical applications.

Conclusion

The algorithms presented here, are easily implemented in any kind of a Hidden Markov Model, whereas the prediction method (HMM-TM) is freely available for academic users at , offering the most advanced decoding options currently available.

Proteomics-based consensus prediction of protein retention in a bacterial membrane

The availability of complete bacterial genome sequences allows proteome-wide predictions of exported proteins that are potentially retained in the cytoplasmic membranes of the corresponding organisms. In practice, however, major problems are encountered with the computer-assisted distinction between (Sec-type) signal peptides that direct exported proteins into the growth medium and lipoprotein signal peptides or amino-terminal membrane anchors that cause protein retention in the membrane. In the present studies, which were aimed at improving methods to predict protein retention in the bacterial cytoplasmic membrane, we have compared sets of membrane-attached and extracellular proteins of Bacillus subtilis that were recently identified through proteomics approaches. The results showed that three classes of membrane-attached proteins can be distinguished. Two classes include 43 lipoproteins and 48 proteins with an amino-terminal transmembrane segment, respectively. Remarkably, a third class includes 31 proteins that remain membrane-retained despite the presence of typical Sec-type signal peptides with consensus signal peptidase recognition sites. This unprecedented finding indicates that unknown mechanisms are involved in membrane retention of this class of proteins. A further novelty is a consensus sequence indicative for release of certain lipoproteins from the membrane by proteolytic shaving. Finally, using non-overlapping sets of secreted and membrane-retained proteins, the accuracy of different signal peptide prediction algorithms was assessed. Accuracy for the prediction of protein retention in the membrane was increased to 82% using a majority-vote approach. Our findings provide important leads for future identification of surface proteins from pathogenic bacteria, which are attractive candidate infection markers and potential targets for drugs or vaccines.

Functional properties and differential mode of regulation of the nitrate transporter from a plant symbiotic ascomycete

Nitrogen assimilation by plant symbiotic fungi plays a central role in the mutualistic interaction established by these organisms, as well as in nitrogen flux in a variety of soils. In the present study, we report on the functional properties, structural organization and distinctive mode of regulation of TbNrt2 (Tuber borchii NRT2 family transporter), the nitrate transporter of the mycorrhizal ascomycete T. borchii. As revealed by experiments conducted in a nitrate-uptake-defective mutant of the yeast Hansenula polymorpha, TbNrt2 is a high-affinity transporter (K(m)=4.7 microM nitrate) that is bispecific for nitrate and nitrite. It is expressed in free-living mycelia and in mycorrhizae, where it preferentially accumulates in the plasma membrane of root-contacting hyphae. The TbNrt2 mRNA, which is transcribed from a single-copy gene clustered with the nitrate reductase gene in the T. borchii genome, was specifically up-regulated following transfer of mycelia to nitrate- (or nitrite)-containing medium. However, at variance with the strict nitrate-dependent induction commonly observed in other organisms, TbNrt2 was also up-regulated (at both the mRNA and the protein level) following transfer to a nitrogen-free medium. This unusual mode of regulation differs from that of the adjacent nitrate reductase gene, which was expressed at basal levels under nitrogen deprivation conditions and required nitrate for induction. The functional and expression properties, described in the present study, delineate TbNrt2 as a versatile transporter that may be especially suited to cope with the fluctuating (and often low) mineral nitrogen concentrations found in most natural, especially forest, soils.

Transmembrane protein topology mapping by the substituted cysteine accessibility method (SCAM(TM)): application to lipid-specific membrane protein topogenesis

We provide an overview of lipid-dependent polytopic membrane protein topogenesis, with particular emphasis on Escherichia coli strains genetically altered in their lipid composition and strategies for experimentally determining the transmembrane organization of proteins. A variety of reagents and experimental strategies are described including the use of lipid mutants and thiol-specific chemical reagents to study lipid-dependent and host-specific membrane protein topogenesis by substituted cysteine site-directed chemical labeling. Employing strains in which lipid composition can be controlled temporally during membrane protein synthesis and assembly provides a means to observe dynamic changes in protein topology as a function of membrane lipid composition.

Dynamic Membrane Topology of the Escherichia coli β-Glucoside Transporter BglF

The Escherichia coli BglF protein, a permease of the phosphoenolpyruvate-dependent phosphotransferase system, catalyzes transport and phosphorylation of beta-glucosides. In addition, BglF regulates bgl operon expression by controlling the activity of the transcriptional regulator BglG via reversible phosphorylation. BglF is composed of three domains; one is hydrophobic, which presumably forms the sugar translocation channel. We studied the topology of this domain by Cys-replacement mutagenesis and chemical modification by thiol reagents. Most Cys substitutions were well tolerated, as demonstrated by the ability of the mutant proteins to catalyze BglF activities. Our results suggest that the membrane domain contains eight transmembrane helices and an alleged cytoplasmic loop that contains two additional helices. The latter region forms a dynamic structure, as evidenced by the alternation of residues near its ends between faced-in and faced-out states. We suggest that this region, together with the two transmembrane helices encompassing it, forms the sugar translocation channel. BglF periplasmic loops are close to the membrane, the first being a reentrant loop. This is the first systematic topological study carried out with an intact phosphotransferase system permease and the first demonstration of a reentrant loop in this group of proteins.

Modeling membrane proteins based on low-resolution electron microscopy maps: a template for the TM domains of the oxalate transporter OxlT

The availability of both EM and high-resolution crystallographic data for several membrane proteins (MPs) permits a detailed evaluation of the ability of molecular modeling techniques to complement EM data in the development of models of MPs. A protocol for this purpose is presented, consisting of (1) identifying transmembrane (TM) domains from sequence; (2) assigning buried and lipid-exposed faces of the TM domains; and (3) assembling the TM domains into a bundle, based on geometric restraints obtained from the EM data. The protocol is validated by predicting the structures of several 7- and 12-TM MPs to within 3-5 A r.m.s.d. from their crystal structures. The protocol is applied to generate a model of the oxalate transporter OxlT, for which a high-resolution structure is not yet available.

Topological analysis of GtrA and GtrB proteins encoded by the serotype-converting cassette of Shigella flexneri

Serotype conversion (O-antigen glucosylation) in Shigella flexneri is mediated by temperate bacteriophages, which encode a three-gene cluster that contains gtrA, gtrB, and gtr([type]). Sequence analysis has revealed that gtrA and gtrB are conserved and readily interchangeable between serotypes. The gtr([type]) is unique in each serotype and responsible for specifically mediating conversion by the addition of a glucosyl group to the O-antigen units. Analysis of the GtrA and GtrB amino acid sequence using computer prediction programs indicated that GtrA and GtrB have four and two transmembrane segments, respectively. The topology model of GtrA was analyzed by constructing consecutive sandwich fusions using a dual reporter PhoA/LacZ at predetermined positions targeting each of the 3 cytoplasmic and 2 periplasmic hypothetical loops. The topology of GtrB was determined by constructing C-terminal truncated fusions of GtrB to full-length PhoA and LacZ by a PCR-mediated method. These approaches revealed that GtrA consists of four transmembrane segments with both the N-terminal and C-terminal ends in the cytoplasm. Accordingly, GtrB consists of two transmembrane segments with both ends also in the cytoplasm. Furthermore, membrane anchorage of the extended N-terminal end of GtrB was found to be important in catalysis. This study completes the topology of all three proteins (GtrA, GtrB, and the gtr([type]): GtrV) involved in the glucosyltransferase activity that results in serotype conversion of S. flexneri. A model is proposed showing how both O-antigen synthesis and modification take place in S. flexneri.

Membrane Topology Mapping of Vitamin K Epoxide Reductase by in Vitro Translation/Cotranslocation

Vitamin K epoxide reductase (VKOR) catalyzes the conversion of vitamin K 2,3-epoxide into vitamin K in the vitamin K redox cycle. Recently, the gene encoding the catalytic subunit of VKOR was identified as a 163-amino acid integral membrane protein. In this study we report the experimentally derived membrane topology of VKOR. Our results show that four hydrophobic regions predicted as the potential transmembrane domains in VKOR can individually insert across the endoplasmic reticulum membrane in vitro. However, in the intact enzyme there are only three transmembrane domains, residues 10-29, 101-123, and 127-149, and membrane-integration of residues 75-97 appears to be suppressed by the surrounding sequence. Results of N-linked glycosylation-tagged full-length VKOR shows that the N terminus of VKOR is located in the endoplasmic reticulum lumen, and the C terminus is located in the cytoplasm. Further evidence for this topological model of VKOR was obtained with freshly prepared intact microsomes from insect cells expressing HPC4-tagged full-length VKOR. In these experiments an HPC4 tag at the N terminus was protected from proteinase K digestion, whereas an HPC4 tag at the C terminus was susceptible. Altogether, our results suggest that VKOR is a type III membrane protein with three transmembrane domains, which agrees well with the prediction by the topology prediction program TMHMM.

Functional properties of the alternative NADH:ubiquinone oxidoreductase from E. coli through comparative 3-D modelling

The alternative NADH:ubiquinone oxidoreductase (NDH-2) from Escherichia coli is a membrane protein playing a prominent role in respiration by linking the reduction of NADH to the quinone pool. Remote sequence similarity reveals an evolutionary relation between alternative NADH:quinone oxidoreductases and the SCOP-family "FAD/NAD-linked reductases". We have created a structural model for NDH-2 from E. coli through comparative modelling onto a template from this family. Combined analysis of our model and sequence conservation allowed us to include the cofactor FAD and the substrate NADH in atomic detail. Furthermore, we propose the most plausible orientation of NDH-2 relative to the membrane and specify a region of the protein potentially involved in ubiquinone binding.

A hidden Markov model with molecular mechanics energy-scoring function for transmembrane helix prediction

A range of methods has been developed to predict transmembrane helices and their topologies. Although most of these algorithms give good predictions, no single method consistently outperforms the others. However, combining different algorithms is one approach that can potentially improve the accuracy of the prediction. We developed a new method that initially uses a hidden Markov model to predict alternative models for membrane spanning helices in proteins. The algorithm subsequently identifies the best among models by ranking them using a novel scoring function based on the folding energy of transmembrane helical fragments. This folding of helical fragments and the incorporation into membrane is modeled using CHARMm, extended with the Generalized Born surface area solvent model (GBSA/IM) with implicit membrane. The combined method reported here, TMHGB significantly increases the accuracy of the original hidden Markov model-based algorithm.

Evaluation of methods for predicting the topology of beta-barrel outer membrane proteins and a consensus prediction method

Background

Prediction of the transmembrane strands and topology of β-barrel outer membrane proteins is of interest in current bioinformatics research. Several methods have been applied so far for this task, utilizing different algorithmic techniques and a number of freely available predictors exist. The methods can be grossly divided to those based on Hidden Markov Models (HMMs), on Neural Networks (NNs) and on Support Vector Machines (SVMs). In this work, we compare the different available methods for topology prediction of β-barrel outer membrane proteins. We evaluate their performance on a non-redundant dataset of 20 β-barrel outer membrane proteins of gram-negative bacteria, with structures known at atomic resolution. Also, we describe, for the first time, an effective way to combine the individual predictors, at will, to a single consensus prediction method.

Results

We assess the statistical significance of the performance of each prediction scheme and conclude that Hidden Markov Model based methods, HMM-B2TMR, ProfTMB and PRED-TMBB, are currently the best predictors, according to either the per-residue accuracy, the segments overlap measure (SOV) or the total number of proteins with correctly predicted topologies in the test set. Furthermore, we show that the available predictors perform better when only transmembrane β-barrel domains are used for prediction, rather than the precursor full-length sequences, even though the HMM-based predictors are not influenced significantly. The consensus prediction method performs significantly better than each individual available predictor, since it increases the accuracy up to 4% regarding SOV and up to 15% in correctly predicted topologies.

Conclusions

The consensus prediction method described in this work, optimizes the predicted topology with a dynamic programming algorithm and is implemented in a web-based application freely available to non-commercial users at .