Predicting weakly stable regions, oligomerization state, and protein-protein interfaces in transmembrane domains of outer membrane proteins

New insights into the structure and dynamics of the TOM complex in mitochondria

To date, there is no general physical model of the mechanism by which unfolded polypeptide chains with different properties are imported into the mitochondria. At the molecular level, it is still unclear how transit polypeptides approach, are captured by the protein translocation machinery in the outer mitochondrial membrane, and how they subsequently cross the entropic barrier of a protein translocation pore to enter the intermembrane space. This deficiency has been due to the lack of detailed structural and dynamic information about the membrane pores. In this review, we focus on the recently determined sub-nanometer cryo-EM structures and our current knowledge of the dynamics of the mitochondrial two-pore outer membrane protein translocation machinery (TOM core complex), which provide a starting point for addressing the above questions. Of particular interest are recent discoveries showing that the TOM core complex can act as a mechanosensor, where the pores close as a result of interaction with membrane-proximal structures. We highlight unusual and new correlations between the structural elements of the TOM complexes and their dynamic behavior in the membrane environment.

Challenges of Protein-Protein Docking of the Membrane Proteins

Despite the recent advances in the determination of high-resolution membrane protein (MP) structures, the structural and functional characterization of MPs remains extremely challenging, mainly due to the hydrophobic nature, low abundance, poor expression, purification, and crystallization difficulties associated with MPs. Whereby the major challenges/hurdles for MP structure determination are associated with the expression, purification, and crystallization procedures. Although there have been significant advances in the experimental determination of MP structures, only a limited number of MP structures (approximately less than 1% of all) are available in the Protein Data Bank (PDB). Therefore, the structures of a large number of MPs still remain unresolved, which leads to the availability of widely unplumbed structural and functional information related to MPs. As a result, recent developments in the drug discovery realm and the significant biological contemplation have led to the development of several novel, low-cost, and time-efficient computational methods that overcome the limitations of experimental approaches, supplement experiments, and provide alternatives for the characterization of MPs. Whereby the fine tuning and optimizations of these computational approaches remains an ongoing endeavor.Computational methods offer a potential way for the elucidation of structural features and the augmentation of currently available MP information. However, the use of computational modeling can be extremely challenging for MPs mainly due to insufficient knowledge of (or gaps in) atomic structures of MPs. Despite the availability of numerous in silico methods for 3D structure determination the applicability of these methods to MPs remains relatively low since all methods are not well-suited or adequate for MPs. However, sophisticated methods for MP structure predictions are constantly being developed and updated to integrate the modifications required for MPs. Currently, different computational methods for (1) MP structure prediction, (2) stability analysis of MPs through molecular dynamics simulations, (3) modeling of MP complexes through docking, (4) prediction of interactions between MPs, and (5) MP interactions with its soluble partner are extensively used. Towards this end, MP docking is widely used. It is notable that the MP docking methods yet few in number might show greater potential in terms of filling the knowledge gap. In this chapter, MP docking methods and associated challenges have been reviewed to improve the applicability, accuracy, and the ability to model macromolecular complexes.

Two conformations of the Tom20 preprotein receptor in the TOM holo complex

The TOM complex is the main entry point for precursor proteins (preproteins) into mitochondria. Preproteins containing targeting sequences are recognized by the TOM complex and imported into mitochondria. We have determined the structure of the TOM core complex from Neurospora crassa by single-particle electron cryomicroscopy at 3.3 Å resolution, showing its interaction with a bound preprotein at 4 Å resolution, and of the TOM holo complex including the Tom20 receptor at 6 to 7 Å resolution. TOM is a transmembrane complex consisting of two β-barrels, three receptor subunits, and three short transmembrane subunits. Tom20 has a transmembrane helix and a receptor domain on the cytoplasmic side. We propose that Tom20 acts as a dynamic gatekeeper, guiding preproteins into the pores of the TOM complex. We analyze the interactions of Tom20 with other TOM subunits, present insights into the structure of the TOM holo complex, and suggest a translocation mechanism.

Protein-Protein Interface Topology as a Predictor of Secondary Structure and Molecular Function Using Convolutional Deep Learning

To power the specific recognition and binding of protein partners into functional complexes, a wealth of information about the structure and function of the partners is necessarily encoded into the global shape of protein-protein interfaces and their local topological features. To identify whether this is the case, this study uses convolutional deep learning methods (typically leveraged for 2D image recognition) on 3D voxel representations of protein-protein interfaces colored by burial depth. A novel two-stage network fed with voxelizations of each interface at two distinct resolutions achieves balance between performance and computational cost. From the shape of the interfaces, the network tries to predict the presence of secondary structure motifs at the interface and the molecular function of the corresponding complex. Secondary structure and certain classes of function are found to be very well predicted, validating the hypothesis that interface shape is a conveyor of higher-level information. Interface patterns triggering the recognition of specific classes are also identified and described.

The structure of the TOM core complex in the mitochondrial outer membrane

In the past three decades, significant advances have been made in providing the biochemical background of TOM (translocase of the outer mitochondrial membrane)-mediated protein translocation into mitochondria. In the light of recent cryoelectron microscopy-derived structures of TOM isolated from Neurospora crassa and Saccharomyces cerevisiae, the interpretation of biochemical and biophysical studies of TOM-mediated protein transport into mitochondria now rests on a solid basis. In this review, we compare the subnanometer structure of N. crassa TOM core complex with that of yeast. Both structures reveal remarkably well-conserved symmetrical dimers of 10 membrane protein subunits. The structural data also validate predictions of weakly stable regions in the transmembrane β-barrel domains of the protein-conducting subunit Tom40, which signal the existence of β-strands located in interfaces of protein-protein interactions.

Marker residue types at the structural regions of transmembrane alpha-helical and beta-barrel interfaces

Membrane proteins play a variety of biological functions to the survival of organisms and functionalities of these proteins are often due to their homo- or hetero-complexation. Encoded by ~30% of the genome in most organisms, they represent the target of over half of nowadays drugs. Spanning the entirety of the cell membrane, transmembrane proteins are the most common type of membrane proteins and can be classified by secondary structures: alpha-helical and beta-barrel structures. Protein-protein interaction (PPI) have been widely studied for globular proteins and many computational tools are available for predicting PPI sites and construct models of complexes. Here, the structural regions of a non-redundant set of 232 alpha-helical and 37 beta-barrel transmembrane complexes and their interfaces are analyzed. Using the residue composition, frequency and propensity, this study brings the light on the marker residue types located at the structural regions of alpha-helical and beta-barrel transmembrane homomeric protein complexes and of their interfaces. This study also shows the necessity to relate the frequency to the composition into a ratio for immediately figuring out residue types presenting high frequencies at the interface and/or at one of its structural regions despite being a minor contributor compared to other residue types to that location's residue composition.

Porin from Marine Bacterium Marinomonas primoryensis KMM 3633T: Isolation, Physico-Chemical Properties, and Functional Activity

Marinomonas primoryensis KMM 3633^T, extreme living marine bacterium was isolated from a sample of coastal sea ice in the Amursky Bay near Vladivostok, Russia. The goal of our investigation is to study outer membrane channels determining cell permeability. Porin from M. primoryensis KMM 3633^T (MpOmp) has been isolated and characterized. Amino acid analysis and whole genome sequencing were the sources of amino acid data of porin, identified as Porin_4 according to the conservative domain searching. The amino acid composition of MpOmp distinguished by high content of acidic amino acids and low content of sulfur-containing amino acids, but there are no tryptophan residues in its molecule. The native MpOmp existed as a trimer. The reconstitution of MpOmp into black lipid membranes demonstrated its ability to form ion channels whose conductivity depends on the electrolyte concentration. The spatial structure of MpOmp had features typical for the classical gram-negative porins. However, the oligomeric structure of isolated MpOmp was distinguished by very low stability: heat-modified monomer was already observed at 30 °C. The data obtained suggest the stabilizing role of lipids in the natural membrane of marine bacteria in the formation of the oligomeric structure of porin.

GeTFEP: A general transfer free energy profile of transmembrane proteins

Free energy of transferring amino acid side-chains from aqueous environment into lipid bilayers, known as transfer free energy (TFE), provides important information on the thermodynamic stability of membrane proteins. In this study, we derived a TFE profile named General Transfer Free Energy Profile (GeTFEP) based on computation of the TFEs of 58 β-barrel membrane proteins (βMPs). The GeTFEP agrees well with experimentally measured and computationally derived TFEs. Analysis based on the GeTFEP shows that residues in different regions of the transmembrane (TM) segments of βMPs have different roles during the membrane insertion process. Results further reveal the importance of the sequence pattern of TM strands in stabilizing βMPs in the membrane environment. In addition, we show that GeTFEP can be used to predict the positioning and the orientation of βMPs in the membrane. We also show that GeTFEP can be used to identify structurally or functionally important amino acid residue sites of βMPs. Furthermore, the TM segments of α-helical membrane proteins can be accurately predicted with GeTFEP, suggesting that the GeTFEP is of general applicability in studying membrane protein.

Biologically Active Ultra-Simple Proteins Reveal Principles of Transmembrane Domain Interactions

Specific interactions between the helical membrane-spanning domains of transmembrane proteins play central roles in the proper folding and oligomerization of these proteins. However, the relationship between the hydrophobic amino acid sequences of transmembrane domains and their functional interactions is in most cases unknown. Here, we use ultra-simple artificial proteins to systematically study the sequence basis for transmembrane domain interactions. We show that most short homopolymeric polyleucine transmembrane proteins containing single amino acid substitutions can activate the platelet-derived growth factor β receptor or the erythropoietin receptor in cultured mouse cells, resulting in cell transformation or proliferation. These proteins displayed complex patterns of activity that were markedly affected by seemingly minor sequence differences in the ultra-simple protein itself or in the transmembrane domain of the target receptor, and the effects of these sequence differences are not additive. In addition, specific leucine residues along the length of these proteins are required for activity, and the positions of these required leucines differ based on the identity and position of the central substituted amino acid. Our results suggest that these ultra-simple proteins use a variety of molecular mechanisms to activate the same target and that diversification of transmembrane domain sequences over the course of evolution minimized off-target interactions.

Quaternary structure of the small amino acid transporter OprG from Pseudomonas aeruginosa

Pseudomonas aeruginosa is an opportunistic human pathogen that causes nosocomial infections. The P. aeruginosa outer membrane contains specific porins that enable substrate uptake, with the outer membrane protein OprG facilitating transport of small, uncharged amino acids. However, the pore size of an eight-stranded β-barrel monomer of OprG is too narrow to accommodate even the smallest transported amino acid, glycine, raising the question of how OprG facilitates amino acid uptake. Pro-92 of OprG is critically important for amino acid transport, with a P92A substitution inhibiting transport and the NMR structure of this variant revealing that this substitution produces structural changes in the barrel rim and restricts loop motions. OprG may assemble into oligomers in the outer membrane (OM) whose subunit interfaces could form a transport channel. Here, we explored the contributions of the oligomeric state and the extracellular loops to OprG's function. Using chemical cross-linking to determine the oligomeric structures of both WT and P92A OprG in native outer membranes and atomic force microscopy, and single-molecule fluorescence of the purified proteins reconstituted into lipid bilayers, we found that both protein variants form oligomers, supporting the notion that subunit interfaces in the oligomer could provide a pathway for amino acid transport. Furthermore, performing transport assays with loop-deleted OprG variants, we found that these variants also can transport small amino acids, indicating that the loops are not solely responsible for substrate transport. We propose that OprG functions as an oligomer and that conformational changes in the barrel-loop region might be crucial for its activity.

Efficient computation of transfer free energies of amino acids in beta-barrel membrane proteins

Motivation

Transmembrane beta-barrel proteins (TMBs) serve a multitude of essential cellular functions in Gram-negative bacteria, mitochondria and chloroplasts. Transfer free energies (TFEs) of residues in the transmembrane (TM) region provides fundamental quantifications of thermodynamic stabilities of TMBs, which are important for the folding and the membrane insertion processes, and may help in understanding the structure-function relationship. However, experimental measurement of TFEs of TMBs is challenging. Although a recent computational method can be used to calculate TFEs, the results of which are in excellent agreement with experimentally measured values, this method does not scale up, and is limited to small TMBs.

Results

We have developed an approximation method that calculates TFEs of TM residues in TMBs accurately, with which depth-dependent transfer free energy profiles can be derived. Our results are in excellent agreement with experimental measurements. This method is efficient and applicable to all bacterial TMBs regardless of the size of the protein.

Availability and Implementation

An online webserver is available at http://tanto.bioe.uic.edu/tmb-tfe .

Contact

: jliang@uic.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

High-resolution structure prediction of β-barrel membrane proteins

[Formula: see text]-Barrel membrane proteins ([Formula: see text]MPs) play important roles, but knowledge of their structures is limited. We have developed a method to predict their 3D structures. We predict strand registers and construct transmembrane (TM) domains of [Formula: see text]MPs accurately, including proteins for which no prediction has been attempted before. Our method also accurately predicts structures from protein families with a limited number of sequences and proteins with novel folds. An average main-chain rmsd of 3.48 Å is achieved between predicted and experimentally resolved structures of TM domains, which is a significant improvement ([Formula: see text]3 Å) over a recent study. For [Formula: see text]MPs with NMR structures, the deviation between predictions and experimentally solved structures is similar to the difference among the NMR structures, indicating excellent prediction accuracy. Moreover, we can now accurately model the extended [Formula: see text]-barrels and loops in non-TM domains, increasing the overall coverage of structure prediction by [Formula: see text]%. Our method is general and can be applied to genome-wide structural prediction of [Formula: see text]MPs.

Engineering a Novel Porin OmpGF Via Strand Replacement from Computational Analysis of Sequence Motif

β-Barrelmembrane proteins (βMPs) form barrel-shaped pores in the outer membrane of Gram-negative bacteria, mitochondria, and chloroplasts. Because of the robustness of their barrel structures, βMPs have great potential as nanosensors for single-molecule detection. However, natural βMPs currently employed have inflexible biophysical properties and are limited in their pore geometry, hindering their applications in sensing molecules of different sizes and properties. Computational engineering has the promise to generate βMPs with desired properties. Here we report a method for engineering novel βMPs based on the discovery of sequence motifs that predominantly interact with the cell membrane and appear in more than 75% of transmembrane strands. By replacing β1-β6 strands of the protein OmpF that lack these motifs with β1-β6 strands of OmpG enriched with these motifs and computational verification of increased stability of its transmembrane section, we engineered a novel βMP called OmpGF. OmpGF is predicted to form a monomer with a stable transmembrane region. Experimental validations showed that OmpGF could refold in vitro with a predominant β-sheet structure, as confirmed by circular dichroism. Evidence of OmpGF membrane insertion was provided by intrinsic tryptophan fluorescence spectroscopy, and its pore-forming property was determined by a dye-leakage assay. Furthermore, single-channel conductance measurements confirmed that OmpGF function as a monomer and exhibits increased conductance than OmpG and OmpF. These results demonstrated that a novel and functional βMP can be successfully engineered through strand replacement based on sequence motif analysis and stability calculation.

Peculiarities of thermal denaturation of OmpF porin from Yersinia ruckeri

Irreversible denaturation of membrane proteins in detergent solutions is similar to unfolding of water-soluble multidomain proteins and represents a complex, multistage process.

Outer membrane protein design

Membrane proteins are the gateway to the cell. These proteins are also a control center of the cell, as information from the outside is passed through membrane proteins as signals to the cellular machinery. The design of membrane proteins seeks to harness the power of these gateways and signal carriers. This review will focus on the design of the membrane proteins that are in the outer membrane, a membrane which only exists for gram negative bacteria, mitochondria, and chloroplasts. Unlike other membrane proteins, outer membrane proteins are uniquely shaped as β-barrels. Herein, I describe most known examples of membrane β-barrel design to date, focusing particularly on categorizing designs as: Firstly, structural deconstruction; secondly, structural changes; thirdly, chemical function design; and finally, the creation of new folds.

Protein folding in the cell envelope of Escherichia coli

While the entire proteome is synthesized on cytoplasmic ribosomes, almost half associates with, localizes in or crosses the bacterial cell envelope. In Escherichia coli a variety of mechanisms are important for taking these polypeptides into or across the plasma membrane, maintaining them in soluble form, trafficking them to their correct cell envelope locations and then folding them into the right structures. The fidelity of these processes must be maintained under various environmental conditions including during stress; if this fails, proteases are called in to degrade mislocalized or aggregated proteins. Various soluble, diffusible chaperones (acting as holdases, foldases or pilotins) and folding catalysts are also utilized to restore proteostasis. These responses can be general, dealing with multiple polypeptides, with functional overlaps and operating within redundant networks. Other chaperones are specialized factors, dealing only with a few exported proteins. Several complex machineries have evolved to deal with binding to, integration in and crossing of the outer membrane. This complex protein network is responsible for fundamental cellular processes such as cell wall biogenesis; cell division; the export, uptake and degradation of molecules; and resistance against exogenous toxic factors. The underlying processes, contributing to our fundamental understanding of proteostasis, are a treasure trove for the development of novel antibiotics, biopharmaceuticals and vaccines.

Outer Membrane Protein Folding and Topology from a Computational Transfer Free Energy Scale

Knowledge of the transfer free energy of amino acids from aqueous solution to a lipid bilayer is essential for understanding membrane protein folding and for predicting membrane protein structure. Here we report a computational approach that can calculate the folding free energy of the transmembrane region of outer membrane β-barrel proteins (OMPs) by combining an empirical energy function with a reduced discrete state space model. We quantitatively analyzed the transfer free energies of 20 amino acid residues at the center of the lipid bilayer of OmpLA. Our results are in excellent agreement with the experimentally derived hydrophobicity scales. We further exhaustively calculated the transfer free energies of 20 amino acids at all positions in the TM region of OmpLA. We found that the asymmetry of the Gram-negative bacterial outer membrane as well as the TM residues of an OMP determine its functional fold in vivo. Our results suggest that the folding process of an OMP is driven by the lipid-facing residues in its hydrophobic core, and its NC-IN topology is determined by the differential stabilities of OMPs in the asymmetrical outer membrane. The folding free energy is further reduced by lipid A and assisted by general depth-dependent cooperativities that exist between polar and ionizable residues. Moreover, context-dependency of transfer free energies at specific positions in OmpLA predict regions important for protein function as well as structural anomalies. Our computational approach is fast, efficient and applicable to any OMP.

Computational redesign of the lipid-facing surface of the outer membrane protein OmpA

Advances in computational design methods have made possible extensive engineering of soluble proteins, but designed β-barrel membrane proteins await improvements in our understanding of the sequence determinants of folding and stability. A subset of the amino acid residues of membrane proteins interact with the cell membrane, and the design rules that govern this lipid-facing surface are poorly understood. We applied a residue-level depth potential for β-barrel membrane proteins to the complete redesign of the lipid-facing surface of Escherichia coli OmpA. Initial designs failed to fold correctly, but reversion of a small number of mutations indicated by backcross experiments yielded designs with substitutions to up to 60% of the surface that did support folding and membrane insertion.

TtOmp85, a β-barrel assembly protein, functions by barrel augmentation

Outer membrane proteins are vital for Gram-negative bacteria and organisms that inherited organelles from them. Proteins from the Omp85/BamA family conduct the insertion of membrane proteins into the outer membrane. We show that an eight-stranded outer membrane β-barrel protein, TtoA, is inserted and folded into liposomes by an Omp85 homologue. Furthermore, we recorded the channel conductance of this Omp85 protein in black lipid membranes, alone and in the presence of peptides comprising the sequence of the two N-terminal and the two C-terminal β-strands of TtoA. Only with the latter could a long-living compound channel that exhibits conductance levels higher than those of the Omp85 protein alone be observed. These data support a model in which unfolded outer membrane protein after docking with its C-terminus penetrates into the transmembrane β-barrel of the Omp85 protein and augments its β-sheet at the first strand. Augmentation with successive β-strands leads to a compound, dilated barrel of both proteins.

Structure and function of POTRA domains of Omp85/TPS superfamily

The Omp85/TPS (outer-membrane protein of 85 kDa/two-partner secretion) superfamily is a ubiquitous and major class of β-barrel proteins. This superfamily is restricted to the outer membranes of gram-negative bacteria, mitochondria, and chloroplasts. The common architecture, with an N-terminus consisting of repeats of soluble polypeptide-transport-associated (POTRA) domains and a C-terminal β-barrel pore is highly conserved. The structures of multiple POTRA domains and one full-length TPS protein have been solved, yet discovering roles of individual POTRA domains has been difficult. This review focuses on similarities and differences between POTRA structures, emphasizing POTRA domains in autotrophic organisms including plants and cyanobacteria. Unique roles, specific for certain POTRA domains, are examined in the context of POTRA location with respect to their attachment to the β-barrel pore, and their degree of biological dispensability. Finally, because many POTRA domains may have the ability to interact with thousands of partner proteins, possible modes of these interactions are also explored.

Allosteric signalling in the outer membrane translocation domain of PapC usher

PapC ushers are outer-membrane proteins enabling assembly and secretion of P pili in uropathogenic E. coli. Their translocation domain is a large β-barrel occluded by a plug domain, which is displaced to allow the translocation of pilus subunits across the membrane. Previous studies suggested that this gating mechanism is controlled by a β-hairpin and an α-helix. To investigate the role of these elements in allosteric signal communication, we developed a method combining evolutionary and molecular dynamics studies of the native translocation domain and mutants lacking the β-hairpin and/or the α-helix. Analysis of a hybrid residue interaction network suggests distinct regions (residue 'communities') within the translocation domain (especially around β12-β14) linking these elements, thereby modulating PapC gating. Antibiotic sensitivity and electrophysiology experiments on a set of alanine-substitution mutants confirmed functional roles for four of these communities. This study illuminates the gating mechanism of PapC ushers and its importance in maintaining outer-membrane permeability.

Toward the functional oligomerization state of tryptophan-rich sensory proteins

A conserved family of tryptophan-rich sensory proteins (TspO) mediates the transport of heme degradation intermediates across membranes. In eukaryotes, the homologous mitochondrial translocator protein (TSPO) binds cholesterol and radioligands as monomer. On the basis of the mammalian TSPO structure, bioinformatic analysis, and a 10 Å resolution electron microscopy map of TspO from Rhodobacter sphaeroides, we developed a model of the tertiary and quaternary structure of TspO that is in agreement with available mutagenesis data. Our study provides insight into the conformational basis for the restricted interaction of bacterial TspO with radioligands and the functional oligomerization state of bacterial TspO proteins.

A 3-dimensional trimeric β-barrel model for Chlamydia MOMP contains conserved and novel elements of Gram-negative bacterial porins

Chlamydia trachomatis is the most prevalent cause of bacterial sexually transmitted diseases and the leading cause of preventable blindness worldwide. Global control of Chlamydia will best be achieved with a vaccine, a primary target for which is the major outer membrane protein, MOMP, which comprises ~60% of the outer membrane protein mass of this bacterium. In the absence of experimental structural information on MOMP, three previously published topology models presumed a16-stranded barrel architecture. Here, we use the latest β-barrel prediction algorithms, previous 2D topology modeling results, and comparative modeling methodology to build a 3D model based on the 16-stranded, trimeric assumption. We find that while a 3D MOMP model captures many structural hallmarks of a trimeric 16-stranded β-barrel porin, and is consistent with most of the experimental evidence for MOMP, MOMP residues 320-334 cannot be modeled as β-strands that span the entire membrane, as is consistently observed in published 16-stranded β-barrel crystal structures. Given the ambiguous results for β-strand delineation found in this study, recent publications of membrane β-barrel structures breaking with the canonical rule for an even number of β-strands, findings of β-barrels with strand-exchanged oligomeric conformations, and alternate folds dependent upon the lifecycle of the bacterium, we suggest that although the MOMP porin structure incorporates canonical 16-stranded conformations, it may have novel oligomeric or dynamic structural changes accounting for the discrepancies observed.

Predictions of Protein-Protein Interfaces within Membrane Protein Complexes

BACKGROUND

Prediction of interaction sites within the membrane protein complexes using the sequence data is of a great importance, because it would find applications in modification of molecules transport through membrane, signaling pathways and drug targets of many diseases. Nevertheless, it has gained little attention from the protein structural bioinformatics community.

METHODS

In this study, a wide variety of prediction and classification tools were applied to distinguish the residues at the interfaces of membrane proteins from those not in the interfaces.

RESULTS

The tuned SVM model achieved the high accuracy of 86.95% and the AUC of 0.812 which outperforms the results of the only previous similar study. Nevertheless, prediction performances obtained using most employed models cannot be used in applied fields and needs more effort to improve.

CONCLUSION

Considering the variety of the applied tools in this study, the present investigation could be a good starting point to develop more efficient tools to predict the membrane protein interaction site residues.

Prediction and design of outer membrane protein-protein interactions

Protein-protein interactions (PPI) play central roles in biological processes, motivating us to understand the structural basis underlying affinity and specificity. In this chapter, we focus on biochemical and computational design strategies of assessing and detecting PPIs of β-barrel outer membrane proteins (OMPs). A few case studies are presented highlighting biochemical techniques used to dissect the energetics of oligomerization and determine amino acids forming the key interactions of the PPI sites. Current computational strategies for detecting/predicting PPIs are introduced, and examples of computational and rational engineering strategies applied to OMPs are presented.

Bacterial origin of a mitochondrial outer membrane protein translocase: new perspectives from comparative single channel electrophysiology

Mitochondria are of bacterial ancestry and have to import most of their proteins from the cytosol. This process is mediated by Tom40, an essential protein that forms the protein-translocating pore in the outer mitochondrial membrane. Tom40 is conserved in virtually all eukaryotes, but its evolutionary origin is unclear because bacterial orthologues have not been identified so far. Recently, it was shown that the parasitic protozoon Trypanosoma brucei lacks a conventional Tom40 and instead employs the archaic translocase of the outer mitochondrial membrane (ATOM), a protein that shows similarities to both eukaryotic Tom40 and bacterial protein translocases of the Omp85 family. Here we present electrophysiological single channel data showing that ATOM forms a hydrophilic pore of large conductance and high open probability. Moreover, ATOM channels exhibit a preference for the passage of cationic molecules consistent with the idea that it may translocate unfolded proteins targeted by positively charged N-terminal presequences. This is further supported by the fact that the addition of a presequence peptide induces transient pore closure. An in-depth comparison of these single channel properties with those of other protein translocases reveals that ATOM closely resembles bacterial-type protein export channels rather than eukaryotic Tom40. Our results support the idea that ATOM represents an evolutionary intermediate between a bacterial Omp85-like protein export machinery and the conventional Tom40 that is found in mitochondria of other eukaryotes.

Predicting protein-protein interface residues using local surface structural similarity

BACKGROUND

Identification of the residues in protein-protein interaction sites has a significant impact in problems such as drug discovery. Motivated by the observation that the set of interface residues of a protein tend to be conserved even among remote structural homologs, we introduce PrISE, a family of local structural similarity-based computational methods for predicting protein-protein interface residues.

RESULTS

We present a novel representation of the surface residues of a protein in the form of structural elements. Each structural element consists of a central residue and its surface neighbors. The PrISE family of interface prediction methods uses a representation of structural elements that captures the atomic composition and accessible surface area of the residues that make up each structural element. Each of the members of the PrISE methods identifies for each structural element in the query protein, a collection of similar structural elements in its repository of structural elements and weights them according to their similarity with the structural element of the query protein. PrISEL relies on the similarity between structural elements (i.e. local structural similarity). PrISEG relies on the similarity between protein surfaces (i.e. general structural similarity). PrISEC, combines local structural similarity and general structural similarity to predict interface residues. These predictors label the central residue of a structural element in a query protein as an interface residue if a weighted majority of the structural elements that are similar to it are interface residues, and as a non-interface residue otherwise. The results of our experiments using three representative benchmark datasets show that the PrISEC outperforms PrISEL and PrISEG; and that PrISEC is highly competitive with state-of-the-art structure-based methods for predicting protein-protein interface residues. Our comparison of PrISEC with PredUs, a recently developed method for predicting interface residues of a query protein based on the known interface residues of its (global) structural homologs, shows that performance superior or comparable to that of PredUs can be obtained using only local surface structural similarity. PrISEC is available as a Web server at http://prise.cs.iastate.edu/

CONCLUSIONS

Local surface structural similarity based methods offer a simple, efficient, and effective approach to predict protein-protein interface residues.

Engineered oligomerization state of OmpF protein through computational design decouples oligomer dissociation from unfolding

Biogenesis of β-barrel membrane proteins is a complex, multistep, and as yet incompletely characterized process. The bacterial porin family is perhaps the best-studied protein family among β-barrel membrane proteins that allows diffusion of small solutes across the bacterial outer membrane. In this study, we have identified residues that contribute significantly to the protein-protein interaction (PPI) interface between the chains of outer membrane protein F (OmpF), a trimeric porin, using an empirical energy function in conjunction with an evolutionary analysis. By replacing these residues through site-directed mutagenesis either with energetically favorable residues or substitutions that do not occur in natural bacterial outer membrane proteins, we succeeded in engineering OmpF mutants with dimeric and monomeric oligomerization states instead of a trimeric oligomerization state. Moreover, our results suggest that the oligomerization of OmpF proceeds through a series of interactions involving two distinct regions of the extensive PPI interface: two monomers interact to form a dimer through the PPI interface near G19. This dimer then interacts with another monomer through the PPI interface near G135 to form a trimer. We have found that perturbing the PPI interface near G19 results in the formation of the monomeric OmpF only. Thermal denaturation of the designed dimeric OmpF mutant suggests that oligomer dissociation can be separated from the process of protein unfolding. Furthermore, the conserved site near G57 and G59 is important for the PPI interface and might provide the essential scaffold for PPIs.

Predicting three-dimensional structures of transmembrane domains of β-barrel membrane proteins

β-Barrel membrane proteins are found in the outer membrane of gram-negative bacteria, mitochondria, and chloroplasts. They are important for pore formation, membrane anchoring, and enzyme activity. These proteins are also often responsible for bacterial virulence. Due to difficulties in experimental structure determination, they are sparsely represented in the protein structure databank. We have developed a computational method for predicting structures of the transmembrane (TM) domains of β-barrel membrane proteins. Based on physical principles, our method can predict structures of the TM domain of β-barrel membrane proteins of novel topology, including those from eukaryotic mitochondria. Our method is based on a model of physical interactions, a discrete conformational state space, an empirical potential function, as well as a model to account for interstrand loop entropy. We are able to construct three-dimensional atomic structure of the TM domains from sequences for a set of 23 nonhomologous proteins (resolution 1.8-3.0 Å). The median rmsd of TM domains containing 75-222 residues between predicted and measured structures is 3.9 Å for main chain atoms. In addition, stability determinants and protein-protein interaction sites can be predicted. Such predictions on eukaryotic mitochondria outer membrane protein Tom40 and VDAC are confirmed by independent mutagenesis and chemical cross-linking studies. These results suggest that our model captures key components of the organization principles of β-barrel membrane protein assembly.

Identification of two voltage-dependent anion channel-like protein sequences conserved in Kinetoplastida

The eukaryotic porin superfamily consists of two families, voltage-dependent anion channel (VDAC) and Tom40, which are both located in the mitochondrial outer membrane. In Trypanosoma brucei, only a single member of the VDAC family has been described. We report the detection of two additional eukaryotic porin-like sequences in T. brucei. By bioinformatic means, we classify both as putative VDAC isoforms.

Structure-based analysis of VDAC1 protein: defining oligomer contact sites

The outer mitochondrial membrane protein, the voltage-dependent anion channel (VDAC), is increasingly implicated in the control of apoptosis. Oligomeric assembly of VDAC1 was shown to be coupled to apoptosis induction, with oligomerization increasing substantially upon apoptosis induction and inhibited by apoptosis blockers. In this study, structure- and computation-based selection of the predicated VDAC1 dimerization site, in combination with site-directed mutagenesis, cysteine replacement, and chemical cross-linking, were employed to identify contact sites between VDAC1 molecules in dimers and higher oligomers. The predicted weakly stable β-strands were experimentally found to represent the interfaces between VDAC1 monomers composing the oligomer. Replacing hydrophobic amino acids with charged residues in β-strands 1, 2, and 19 interfered with VDAC1 oligomerization. The proximity of β-strands 1, 2, and 19 within the VDAC1 dimer and the existence of other association sites involving β-strand 16 were confirmed when a cysteine was introduced at defined positions in cysteineless VDAC1 mutants, together with the use of cysteine-specific cross-linker bis(maleimido)ethane. Moreover, the results suggest that VDAC1 also exists as a dimer that upon apoptosis induction undergoes conformational changes and that its oligomerization proceeds through a series of interactions involving two distinct interfaces. Dissection of VDAC1 dimerization/oligomerization as presented here provides structural insight into the oligomeric status of cellular VDAC1 under physiological and apoptotic conditions.

A knowledge-based potential highlights unique features of membrane α-helical and β-barrel protein insertion and folding

Outer membrane β-barrel proteins differ from α-helical inner membrane proteins in lipid environment, secondary structure, and the proposed processes of folding and insertion. It is reasonable to expect that outer membrane proteins may contain primary sequence information specific for their folding and insertion behavior. In previous work, a depth-dependent insertion potential, E(z) , was derived for α-helical inner membrane proteins. We have generated an equivalent potential for TM β-barrel proteins. The similarities and differences between these two potentials provide insight into unique aspects of the folding and insertion of β-barrel membrane proteins. This potential can predict orientation within the membrane and identify functional residues involved in intermolecular interactions.

Pattern of amino acid substitutions in transmembrane domains of β-barrel membrane proteins for detecting remote homologs in bacteria and mitochondria

-barrel membrane proteins play an important role in controlling the exchange and transport of ions and organic molecules across bacterial and mitochondrial outer membranes. They are also major regulators of apoptosis and are important determinants of bacterial virulence. In contrast to -helical membrane proteins, their evolutionary pattern of residue substitutions has not been quantified, and there are no scoring matrices appropriate for their detection through sequence alignment. Using a Bayesian Monte Carlo estimator, we have calculated the instantaneous substitution rates of transmembrane domains of bacterial -barrel membrane proteins. The scoring matrices constructed from the estimated rates, called bbTM for -barrel Transmembrane Matrices, improve significantly the sensitivity in detecting homologs of -barrel membrane proteins, while avoiding erroneous selection of both soluble proteins and other membrane proteins of similar composition. The estimated evolutionary patterns are general and can detect -barrel membrane proteins very remote from those used for substitution rate estimation. Furthermore, despite the separation of 2–3 billion years since the proto-mitochondrion entered the proto-eukaryotic cell, mitochondria outer membrane proteins in eukaryotes can also be detected accurately using these scoring matrices derived from bacteria. This is consistent with the suggestion that there is no eukaryote-specific signals for translocation. With these matrices, remote homologs of -barrel membrane proteins with known structures can be reliably detected at genome scale, allowing construction of high quality structural models of their transmembrane domains, at the rate of 131 structures per template protein. The scoring matrices will be useful for identification, classification, and functional inference of membrane proteins from genome and metagenome sequencing projects. The estimated substitution pattern will also help to identify key elements important for the structural and functional integrity of -barrel membrane proteins, and will aid in the design of mutagenesis studies.

Computational studies of membrane proteins: models and predictions for biological understanding

We discuss recent progresses in computational studies of membrane proteins based on physical models with parameters derived from bioinformatics analysis. We describe computational identification of membrane proteins and prediction of their topology from sequence, discovery of sequence and spatial motifs, and implications of these discoveries. The detection of evolutionary signal for understanding the substitution pattern of residues in the TM segments and for sequence alignment is also discussed. We further discuss empirical potential functions for energetics of inserting residues in the TM domain, for interactions between TM helices or strands, and their applications in predicting lipid-facing surfaces of the TM domain. Recent progresses in structure predictions of membrane proteins are also reviewed, with further discussions on calculation of ensemble properties such as melting temperature based on simplified state space model. Additional topics include prediction of oligomerization state of membrane proteins, identification of the interfaces for protein-protein interactions, and design of membrane proteins. This article is part of a Special Issue entitled: Protein Folding in Membranes.

The Bam machine: a molecular cooper

The bacterial outer membrane (OM) is an exceptional biological structure with a unique composition that contributes significantly to the resiliency of Gram-negative bacteria. Since all OM components are synthesized in the cytosol, the cell must efficiently transport OM-specific lipids and proteins across the cell envelope and stably integrate them into a growing membrane. In this review, we discuss the challenges associated with these processes and detail the elegant solutions that cells have evolved to address the topological problem of OM biogenesis. Special attention will be paid to the Bam machine, a highly conserved multiprotein complex that facilitates OM β-barrel folding. This article is part of a Special Issue entitled: Protein Folding in Membranes.

Improving the resistance of a eukaryotic β-barrel protein to thermal and chemical perturbations

β-Barrel membrane proteins have regular structures with extensive hydrogen-bond networks between their transmembrane (TM) β-strands, which stabilize their protein fold. Nevertheless, weakly stable TM regions, which are important for the protein function and interaction with other proteins, exist. Here, we report on the apparent stability of human Tom40A, a member of the "mitochondrial porin family" and main constituent of the mitochondrial protein-conducting channel TOM (translocase of the outer membrane). Using a physical interaction model, TmSIP, for β-barrel membrane proteins, we have identified three unfavorable β-strands in the TM domain of the protein. Substitution of key residues inside these strands with hydrophobic amino acids results in a decreased sensitivity of the protein to chemical and/or thermal denaturation. The apparent melting temperature observed when denatured at a rate of 1 °C per minute is shifted from 73 to 84 °C. Moreover, the sensitivity of the protein to denaturant agents is significantly lowered. Further, we find a reduced tendency for the mutated protein to form dimers. We propose that the identified weakly stable β-strands 1, 2 and 9 of human Tom40A play an important role in quaternary protein-protein interactions within the mammalian TOM machinery. Our results show that the use of empirical energy functions to model the apparent stability of β-barrel membrane proteins may be a useful tool in the field of nanopore bioengineering.

Prediction of the exposure status of transmembrane beta barrel residues from protein sequence

We present BTMX (Beta barrel TransMembrane eXposure), a computational method to predict the exposure status (i.e. exposed to the bilayer or hidden in the protein structure) of transmembrane residues in transmembrane beta barrel proteins (TMBs). BTMX predicts the exposure status of known TM residues with an accuracy of 84.2% over 2,225 residues and provides a confidence score for all predictions. Predictions made are in concert with the fact that hydrophobic residues tend to be more exposed to the bilayer. The biological relevance of the input parameters is also discussed. The highest prediction accuracy is obtained when a sliding window comprising three residues with similar C(α)-C(β) vector orientations is employed. The prediction accuracy of the BTMX method on a separate unseen non-redundant test dataset is 78.1%. By employing out-pointing residues that are exposed to the bilayer, we have identified various physico-chemical properties that show statistically significant differences between the beta strands located at the oligomeric interfaces compared to the non-oligomeric strands. The BTMX web server generates colored, annotated snake-plots as part of the prediction results and is available under the BTMX tab at http://service.bioinformatik.uni-saarland.de/tmx-site/. Exposure status prediction of TMB residues may be useful in 3D structure prediction of TMBs.

PagP crystallized from SDS/cosolvent reveals the route for phospholipid access to the hydrocarbon ruler

Enzymatic reactions involving bilayer lipids occur in an environment with strict physical and topological constraints. The integral membrane enzyme PagP transfers a palmitoyl group from a phospholipid to lipid A in order to assist Escherichia coli in evading host immune defenses during infection. PagP measures the palmitoyl group with an internal hydrocarbon ruler that is formed in the interior of the eight-stranded antiparallel β barrel. The access and egress of the palmitoyl group is thought to take a lateral route from the bilayer phase to the barrel interior. Molecular dynamics, mutagenesis, and a 1.4 A crystal structure of PagP in an SDS / 2-methyl-2,4-pentanediol (MPD) cosolvent system reveal that phospholipid access occurs at the crenel present between strands F and G of PagP. In this way, the phospholipid head group can remain exposed to the cell exterior while the lipid acyl chain remains in a predominantly hydrophobic environment as it translocates to the protein interior.

Lipid-binding surfaces of membrane proteins: evidence from evolutionary and structural analysis

Membrane proteins function in the diverse environment of the lipid bilayer. Experimental evidence suggests that some lipid molecules bind tightly to specific sites on the membrane protein surface. These lipid molecules often act as co-factors and play important functional roles. In this study, we have assessed the evolutionary selection pressure experienced at lipid-binding sites in a set of α-helical and β-barrel membrane proteins using posterior probability analysis of the ratio of synonymous vs. nonsynonymous substitutions (ω-ratio). We have also carried out a geometric analysis of the membrane protein structures to identify residues in close contact with co-crystallized lipids. We found that residues forming cholesterol-binding sites in both β(2)-adrenergic receptor and Na(+)-K(+)-ATPase exhibit strong conservation, which can be characterized by an expanded cholesterol consensus motif for GPCRs. Our results suggest the functional importance of aromatic stacking interactions and interhelical hydrogen bonds in facilitating protein-cholesterol interactions, which is now reflected in the expanded motif. We also find that residues forming the cardiolipin-binding site in formate dehydrogenase-N γ-subunit and the phosphatidylglycerol binding site in KcsA are under strong purifying selection pressure. Although the lipopolysaccharide (LPS)-binding site in ferric hydroxamate uptake receptor (FhuA) is only weakly conserved, we show using a statistical mechanical model that LPS binds to the least stable FhuA β-strand and protects it from the bulk lipid. Our results suggest that specific lipid binding may be a general mechanism employed by β-barrel membrane proteins to stabilize weakly stable regions. Overall, we find that the residues forming specific lipid binding sites on the surfaces of membrane proteins often experience strong purifying selection pressure.

Functional dynamics in the voltage-dependent anion channel

The voltage-dependent anion channel (VDAC), located in the outer mitochondrial membrane, acts as a gatekeeper for the entry and exit of mitochondrial metabolites. Here we reveal functional dynamics of isoform one of VDAC (VDAC1) by a combination of solution NMR spectroscopy, Gaussian network model analysis, and molecular dynamics simulation. Micro- to millisecond dynamics are significantly increased for the N-terminal six β-strands of VDAC1 in micellar solution, in agreement with increased B-factors observed in the same region in the bicellar crystal structure of VDAC1. Molecular dynamics simulations reveal that a charge on the membrane-facing glutamic acid 73 (E73) accounts for the elevation of N-terminal protein dynamics as well as a thinning of the nearby membrane. Mutation or chemical modification of E73 strongly reduces the micro- to millisecond dynamics in solution. Because E73 is necessary for hexokinase-I-induced VDAC channel closure and inhibition of apoptosis, our results imply that micro- to millisecond dynamics in the N-terminal part of the barrel are essential for VDAC interaction and gating.

Mechanisms of protein oligomerization, the critical role of insertions and deletions in maintaining different oligomeric states

The main principles of protein-protein recognition are elucidated by the studies of homooligomers which in turn mediate and regulate gene expression, activity of enzymes, ion channels, receptors, and cell-cell adhesion processes. Here we explore oligomeric states of homologous proteins in various organisms to better understand the functional roles and evolutionary mechanisms of homooligomerization. We observe a great diversity in mechanisms controlling oligomerization and focus in our study on insertions and deletions in homologous proteins and how they enable or disable complex formation. We show that insertions and deletions which differentiate monomers and dimers have a significant tendency to be located on the interaction interfaces and about a quarter of all proteins studied and forty percent of enzymes have regions which mediate or disrupt the formation of oligomers. We suggest that relatively small insertions or deletions may have a profound effect on complex stability and/or specificity. Indeed removal of complex enabling regions from protein structures in many cases resulted in the complete or partial loss of stability. Moreover, we find that insertions and deletions modulating oligomerization have a lower aggregation propensity and contain a larger fraction of polar, charged residues, glycine and proline compared to conventional interfaces and protein surface. Most likely, these regions may mediate specific interactions, prevent nonspecific dysfunctional aggregation and preclude undesired interactions between close paralogs therefore separating their functional pathways. Last, we show how the presence or absence of insertions and deletions on interfaces might be of practical value in annotating protein oligomeric states.

The prediction and characterization of YshA, an unknown outer-membrane protein from Salmonella typhimurium

We have developed an effective pathway for the prediction and characterization of novel transmembrane β-barrel proteins. The Freeman-Wimley algorithm, which is a highly accurate prediction method based on the physicochemical properties of experimentally characterized transmembrane β barrel (TMBB) structures, was used to predict TMBBs in the genome of Salmonella typhimurium LT2. The previously uncharacterized product of gene yshA was tested as a model for validating the algorithm. YshA is a highly conserved 230-residue protein that is predicted to have 10 transmembrane β-strands and an N-terminal signal sequence. All of the physicochemical and spectroscopic properties exhibited by YshA are consistent with the prediction that it is a TMBB. Specifically, recombinant YshA localizes to the outer membrane when expressed in Escherichia coli; YshA has a β-sheet-rich secondary structure with stable tertiary contacts in the presence of detergent micelles or when reconstituted into a lipid bilayer. When in a lipid bilayer, YshA forms a membrane-spanning pore with an effective radius of ~0.7nm. Taken together, these data substantiate the predictions made by the Freeman-Wimley algorithm by showing that YshA is a TMBB protein.

Porins in prokaryotes and eukaryotes: common themes and variations

Gram-negative bacteria and mitochondria are both covered by two distinct biological membranes. These membrane systems have been maintained during the course of evolution from an early evolutionary precursor. Both outer membranes accommodate channels of the porin family, which are designed for the uptake and exchange of metabolites, including ions and small molecules, such as nucleosides or sugars. In bacteria, the structure of the outer membrane porin protein family of β-barrels is generally characterized by an even number of β-strands; usually 14, 16 or 18 strands are observed forming the bacterial porin barrel wall. In contrast, the recent structures of the mitochondrial porin, also known as VDAC (voltage-dependent anion channel), show an uneven number of 19 β-strands, but a similar molecular architecture. Despite the lack of a clear evolutionary link between these protein families, their common principles and differences in assembly, architecture and function are summarized in the present review.

VDAC, a multi-functional mitochondrial protein regulating cell life and death

Research over the past decade has extended the prevailing view of the mitochondrion to include functions well beyond the generation of cellular energy. It is now recognized that mitochondria play a crucial role in cell signaling events, inter-organellar communication, aging, cell proliferation, diseases and cell death. Thus, mitochondria play a central role in the regulation of apoptosis (programmed cell death) and serve as the venue for cellular decisions leading to cell life or death. One of the mitochondrial proteins controlling cell life and death is the voltage-dependent anion channel (VDAC), also known as mitochondrial porin. VDAC, located in the mitochondrial outer membrane, functions as gatekeeper for the entry and exit of mitochondrial metabolites, thereby controlling cross-talk between mitochondria and the rest of the cell. VDAC is also a key player in mitochondria-mediated apoptosis. Thus, in addition to regulating the metabolic and energetic functions of mitochondria, VDAC appears to be a convergence point for a variety of cell survival and cell death signals mediated by its association with various ligands and proteins. In this article, we review what is known about the VDAC channel in terms of its structure, relevance to ATP rationing, Ca(2+) homeostasis, protection against oxidative stress, regulation of apoptosis, involvement in several diseases and its role in the action of different drugs. In light of our recent findings and the recently solved NMR- and crystallography-based 3D structures of VDAC1, the focus of this review will be on the central role of VDAC in cell life and death, addressing VDAC function in the regulation of mitochondria-mediated apoptosis with an emphasis on structure-function relations. Understanding structure-function relationships of VDAC is critical for deciphering how this channel can perform such a variety of functions, all important for cell life and death. This review also provides insight into the potential of VDAC1 as a rational target for new therapeutics.