A frequent, GxxxG-mediated, transmembrane association motif is optimized for the formation of interhelical Cα-H hydrogen bonds

Toward high-resolution computational design of the structure and function of helical membrane proteins

The computational design of α-helical membrane proteins is still in its infancy but has already made great progress. De novo design allows stable, specific and active minimal oligomeric systems to be obtained. Computational reengineering can improve the stability and function of naturally occurring membrane proteins. Currently, the major hurdle for the field is the experimental characterization of the designs. The emergence of new structural methods for membrane proteins will accelerate progress.

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.

A cell-free method for expressing and reconstituting membrane proteins enables functional characterization of the plant receptor-like protein kinase FERONIA

There are more than 600 receptor-like kinases (RLKs) in Arabidopsis, but due to challenges associated with the characterization of membrane proteins, only a few have known biological functions. The plant RLK FERONIA is a peptide receptor and has been implicated in plant growth regulation, but little is known about its molecular mechanism of action. To investigate the properties of this enzyme, we used a cell-free wheat germ-based expression system in which mRNA encoding FERONIA was co-expressed with mRNA encoding the membrane scaffold protein variant MSP1D1. With the addition of the lipid cardiolipin, assembly of these proteins into nanodiscs was initiated. FERONIA protein kinase activity in nanodiscs was higher than that of soluble protein and comparable with other heterologously expressed protein kinases. Truncation experiments revealed that the cytoplasmic juxtamembrane domain is necessary for maximal FERONIA activity, whereas the transmembrane domain is inhibitory. An ATP analogue that reacts with lysine residues inhibited catalytic activity and labeled four lysines; mutagenesis demonstrated that two of these, Lys-565 and Lys-663, coordinate ATP in the active site. Mass spectrometric phosphoproteomic measurements further identified phosphorylation sites that were examined using phosphomimetic mutagenesis. The results of these experiments are consistent with a model in which kinase-mediated phosphorylation within the C-terminal region is inhibitory and regulates catalytic activity. These data represent a step further toward understanding the molecular basis for the protein kinase catalytic activity of FERONIA and show promise for future characterization of eukaryotic membrane proteins.

Membranome: a database for proteome-wide analysis of single-pass membrane proteins

The Membranome database was developed to assist analysis and computational modeling of single-pass (bitopic) transmembrane (TM) proteins and their complexes by providing structural information about these proteins on a genomic scale. The database currently collects data on >6000 bitopic proteins from Homo sapiens, Arabidopsis thaliana, Dictyostelium discoideum, Saccharomyces cerevisiae, Escherichia coli and Methanocaldococcus jannaschii It presents the following data: (i) hierarchical classification of bitopic proteins into 15 functional classes, 689 structural superfamilies and 1404 families; (ii) 446 complexes of bitopic proteins with known three-dimensional (3D) structures classified into 129 families; (iii) computationally generated three-dimensional models of TM α-helices positioned in membranes; (iv) amino acid sequences, domain architecture, functional annotation and available experimental structures of bitopic proteins; (v) TM topology and intracellular localization, (vi) physical interactions between proteins from the database along with links to other resources. The database is freely accessible at http://membranome.org There is a variety of options for browsing, sorting, searching and retrieval of the content, including downloadable coordinate files of TM domains with calculated membrane boundaries.

Interplay of Specific Trans- and Juxtamembrane Interfaces in Plexin A3 Dimerization and Signal Transduction

Plexins are transmembrane proteins that serve as guidance receptors during angiogenesis, lymphangiogenesis, neuronal development, and zebrafish fin regeneration, with a putative role in cancer metastasis. Receptor dimerization or clustering, induced by extracellular ligand binding but modulated in part by the plexin transmembrane (TM) and juxtamembrane (JM) domains, is thought to drive plexin activity. Previous studies indicate that isolated plexin TM domains interact through a conserved, small-x3-small packing motif, and the cytosolic JM region interacts through a hydrophobic heptad repeat; however, the roles and interplay of these regions in plexin signal transduction remain unclear. Using an integrated experimental and simulation approach, we find disruption of the small-x3-small motifs in the Danio rerio Plexin A3 TM domain enhances dimerization of the TM-JM domain by enhancing JM-mediated dimerization. Furthermore, mutations of the cytosolic JM heptad repeat that disrupt dimerization do so even in the presence of TM domain mutations. However, mutations to the small-x3-small TM interfaces also disrupt Plexin A3 signaling in a zebrafish axonal guidance assay, indicating the importance of this TM interface in signal transduction. Collectively, our experimental and simulation results demonstrate that multiple TM and JM interfaces exist in the Plexin A3 homodimer, and these interfaces independently regulate dimerization that is important in Plexin A3 signal transduction.

Impact of Differential Detergent Interactions on Transmembrane Helix Dimerization Affinities

Interactions between transmembrane (TM) helices play a critical role in the fundamental processes required for cells to communicate and exchange materials with their surroundings. Our understanding of the factors that promote TM helix interactions has greatly benefited from our ability to study these interactions in the solution phase through the use of membrane-mimetic micelles. However, less is known about the potential influence of juxtamembrane regions flanking the interacting TM helices that may modulate dimerization affinities, even when the interacting surface itself is not altered. To investigate this question, we used solution NMR to quantitate the dimerization affinity of the major coat protein from the M13 bacteriophage in sodium dodecyl sulfate (SDS), a well-characterized model of a single-spanning self-associating TM protein. Here, we showed that a shorter construct lacking the N-terminal amphipathic helix has a higher dimerization affinity relative to that of the full-length protein, with no change in the helical structure between the monomeric and dimeric states in both cases. Although this translated into a 0.6 kcal/mol difference in free energy when the SDS solvent was approximated as a continuous phase, there were deviations from this model at high protein to detergent ratios. Instead, the equilibria were better fit to a model that treats the empty micelle as an active participant in the reaction, giving rise to standard free energies of association that were the same for both full-length and TM-segment constructs. According to this model, the higher apparent affinity of the shorter peptide could be completely explained by the enhanced detergent binding by the monomer relative to that bound by the dimer. Therefore, differential detergent binding between the monomeric and dimeric states provides a mechanism by which TM helix interactions can be modulated by noninteracting juxtamembrane regions.

Understanding single-pass transmembrane receptor signaling from a structural viewpoint-what are we missing?

Single-pass transmembrane receptors are involved in essential processes of both physiological and pathological nature and represent more than 1300 proteins in the human genome. Despite the high biological relevance of these receptors, the mechanisms of the signal transductions they facilitate are incompletely understood. One major obstacle is the lack of structures of the transmembrane domains that connect the extracellular ligand-binding domains to the intracellular signaling platforms. Over a period of almost 20 years since the first structure was reported, only 21 of these receptors have become represented by a transmembrane domain structure. This scarceness stands in strong contrast to the significance of these transmembrane α-helices for receptor functionality. In this review, we explore the properties and qualities of the current set of structures, as well as the methodological difficulties associated with their characterization and the challenges left to be overcome. Without an increased and focused effort to bring this class of proteins on par with the remaining membrane protein field, a serious lag in their biological understanding looms. Design of pharmaceutical agents, prediction of mutational affects in relation to disease, and deciphering of functional mechanisms require high-resolution structural information, especially when dealing with a domain carrying so much functionality in so few residues.

Conformational Heterogeneity of Bax Helix 9 Dimer for Apoptotic Pore Formation

Helix α9 of Bax protein can dimerize in the mitochondrial outer membrane (MOM) and lead to apoptotic pores. However, it remains unclear how different conformations of the dimer contribute to the pore formation on the molecular level. Thus we have investigated various conformational states of the α9 dimer in a MOM model — using computer simulations supplemented with site-specific mutagenesis and crosslinking of the α9 helices. Our data not only confirmed the critical membrane environment for the α9 stability and dimerization, but also revealed the distinct lipid-binding preference of the dimer in different conformational states. In our proposed pathway, a crucial iso-parallel dimer that mediates the conformational transition was discovered computationally and validated experimentally. The corroborating evidence from simulations and experiments suggests that, helix α9 assists Bax activation via the dimer heterogeneity and interactions with specific MOM lipids, which eventually facilitate proteolipidic pore formation in apoptosis regulation.

A specific, transmembrane interface regulates fibroblast activation protein (FAP) homodimerization, trafficking and exopeptidase activity

Fibroblast activation protein (FAP) is a cell-surface serine protease which promotes invasiveness of certain epithelial cancers and is therefore a potential target for cancer drug development and delivery. Unlike dipeptidyl peptidase IV (DPPIV), FAP exhibits prolyl endopeptidase activity and is active as a homodimer with specificity for type I collagen. The mechanism that regulates FAP homodimerization and its relation to prolyl endopeptidase activity is not completely understood. Here, we investigate key residues in the FAP TM domain that may be significant for FAP homodimerization. Mutations to predicted TM interfacial residues (G10L, S14L, and A18L) comprising a small-X3-small motif reduced FAP TM-CYTO dimerization relative to wild type as measured using the AraTM assay, whereas predicted off-interface residues showed no significant change from wild type. The results implied that the predicted small-X3-small dimer interface affect stabilization of FAP TM-CYTO homodimerization. Compared with FAPwild-type, the interfacial TM residue G10L significantly decreased FAP endopeptidase activity more than 25%, and also reduced cell-surface versus intracellular expression relative to other interfacial residues S14L and A18L. Thus, our results suggest FAP dimerization is important for both trafficking and protease activity, and is dependent on a specific TM interface.

Identification of a novel cell death-inducing domain reveals that fungal amyloid-controlled programmed cell death is related to necroptosis

Recent findings have revealed the role of prion-like mechanisms in the control of host defense and programmed cell death cascades. In fungi, HET-S, a cell death-inducing protein containing a HeLo pore-forming domain, is activated through amyloid templating by a Nod-like receptor (NLR). Here we characterize the HELLP protein behaving analogously to HET-S and bearing a new type of N-terminal cell death-inducing domain termed HeLo-like (HELL) and a C-terminal regulatory amyloid motif known as PP. The gene encoding HELLP is part of a three-gene cluster also encoding a lipase (SBP) and a Nod-like receptor, both of which display the PP motif. The PP motif is similar to the RHIM amyloid motif directing formation of the RIP1/RIP3 necrosome in humans. The C-terminal region of HELLP, HELLP(215-278), encompassing the motif, allows prion propagation and assembles into amyloid fibrils, as demonstrated by X-ray diffraction and FTIR analyses. Solid-state NMR studies reveal a well-ordered local structure of the amyloid core residues and a primary sequence that is almost entirely arranged in a rigid conformation, and confirm a β-sheet structure in an assigned stretch of three amino acids. HELLP is activated by amyloid templating and displays membrane-targeting and cell death-inducing activity. HELLP targets the SBP lipase to the membrane, suggesting a synergy between HELLP and SBP in membrane dismantling. Remarkably, the HeLo-like domain of HELLP is homologous to the pore-forming domain of MLKL, the cell death-execution protein in necroptosis, revealing a transkingdom evolutionary relationship between amyloid-controlled fungal programmed cell death and mammalian necroptosis.

A Review of Functional Motifs Utilized by Viruses

Short linear motifs (SLiM) are short peptides that facilitate protein function and protein-protein interactions. Viruses utilize these motifs to enter into the host, interact with cellular proteins, or egress from host cells. Studying functional motifs may help to predict protein characteristics, interactions, or the putative cellular role of a protein. In virology, it may reveal aspects of the virus tropism and help find antiviral therapeutics. This review highlights the recent understanding of functional motifs utilized by viruses. Special attention was paid to the function of proteins harboring these motifs, and viruses encoding these proteins. The review highlights motifs involved in (i) immune response and post-translational modifications (e.g., ubiquitylation, SUMOylation or ISGylation); (ii) virus-host cell interactions, including virus attachment, entry, fusion, egress and nuclear trafficking; (iii) virulence and antiviral activities; (iv) virion structure; and (v) low-complexity regions (LCRs) or motifs enriched with residues (Xaa-rich motifs).

A topological and conformational stability alphabet for multipass membrane proteins

Multipass membrane proteins perform critical signal transduction and transport across membranes. How transmembrane helix (TMH) sequences encode the topology and conformational flexibility regulating these functions remains poorly understood. Here we describe a comprehensive analysis of the sequence-structure relationships at multiple interacting TMHs from all membrane proteins with structures in the Protein Data Bank (PDB). We found that membrane proteins can be deconstructed in interacting TMH trimer units, which mostly fold into six distinct structural classes of topologies and conformations. Each class is enriched in recurrent sequence motifs from functionally unrelated proteins, revealing unforeseen consensus and evolutionary conserved networks of stabilizing interhelical contacts. Interacting TMHs' topology and local protein conformational flexibility were remarkably well predicted in a blinded fashion from the identified binding-hotspot motifs. Our results reveal universal sequence-structure principles governing the complex anatomy and plasticity of multipass membrane proteins that may guide de novo structure prediction, design, and studies of folding and dynamics.

BH3-in-groove dimerization initiates and helix 9 dimerization expands Bax pore assembly in membranes

Pro-apoptotic Bax induces mitochondrial outer membrane permeabilization (MOMP) by forming oligomers through a largely undefined process. Using site-specific disulfide crosslinking, compartment-specific chemical labeling, and mutational analysis, we found that activated integral membrane Bax proteins form a BH3-in-groove dimer interface on the MOM surface similar to that observed in crystals. However, after the α5 helix was released into the MOM, the remaining interface with α2, α3, and α4 helices was rearranged. Another dimer interface was formed inside the MOM by two intersected or parallel α9 helices. Combinations of these interfaces generated oligomers in the MOM. Oligomerization was initiated by BH3-in-groove dimerization, without which neither the other dimerizations nor MOMP occurred. In contrast, α9 dimerization occurred downstream and was required for release of large but not small proteins from mitochondria. Moreover, the release of large proteins was facilitated by α9 insertion into the MOM and localization to the pore rim. Therefore, the BH3-in-groove dimerization on the MOM nucleates the assembly of an oligomeric Bax pore that is enlarged by α9 dimerization at the rim.

Role of the Lipid Environment in the Dimerization of Transmembrane Domains of Glycophorin A

An efficient computational approach is developed to quantify the free energy of a spontaneous association of the α-helices of proteins in the membrane environment. The approach is based on the numerical decomposition of the free energy profiles of the transmembrane (TM) helices into components corresponding to protein-protein, protein-lipid, and protein-water interactions. The method was tested for the TM segments of human glycophorin A (GpA) and two mutant forms, Gly83Ala and Thr87Val. It was shown that lipids make a significant negative contribution to the free energy of dimerization, while amino acid residues forming the interface of the helix-helix contact may be unfavorable in terms of free energy. The detailed balance between different energy contributions is highly dependent on the amino acid sequence of the TM protein segment. The results show the dominant role of the environment in the interaction of membrane proteins that is changing our notion of the driving force behind the spontaneous association of TM α-helices. Adequate estimation of the contribution of the water-lipid environment thus becomes an extremely urgent task for a rational design of new molecules targeting bitopic membrane proteins, including receptor tyrosine kinases.

Role of GxxxG Motifs in Transmembrane Domain Interactions

Transmembrane (TM) helices of integral membrane proteins can facilitate strong and specific noncovalent protein-protein interactions. Mutagenesis and structural analyses have revealed numerous examples in which the interaction between TM helices of single-pass membrane proteins is dependent on a GxxxG or (small)xxx(small) motif. It is therefore tempting to use the presence of these simple motifs as an indicator of TM helix interactions. In this Current Topic review, we point out that these motifs are quite common, with more than 50% of single-pass TM domains containing a (small)xxx(small) motif. However, the actual interaction strength of motif-containing helices depends strongly on sequence context and membrane properties. In addition, recent studies have revealed several GxxxG-containing TM domains that interact via alternative interfaces involving hydrophobic, polar, aromatic, or even ionizable residues that do not form recognizable motifs. In multipass membrane proteins, GxxxG motifs can be important for protein folding, and not just oligomerization. Our current knowledge thus suggests that the presence of a GxxxG motif alone is a weak predictor of protein dimerization in the membrane.

Viroporins, Examples of the Two-Stage Membrane Protein Folding Model

Viroporins are small, α-helical, hydrophobic virus encoded proteins, engineered to form homo-oligomeric hydrophilic pores in the host membrane. Viroporins participate in multiple steps of the viral life cycle, from entry to budding. As any other membrane protein, viroporins have to find the way to bury their hydrophobic regions into the lipid bilayer. Once within the membrane, the hydrophobic helices of viroporins interact with each other to form higher ordered structures required to correctly perform their porating activities. This two-step process resembles the two-stage model proposed for membrane protein folding by Engelman and Poppot. In this review we use the membrane protein folding model as a leading thread to analyze the mechanism and forces behind the membrane insertion and folding of viroporins. We start by describing the transmembrane segment architecture of viroporins, including the number and sequence characteristics of their membrane-spanning domains. Next, we connect the differences found among viroporin families to their viral genome organization, and finalize focusing on the pathways used by viroporins in their way to the membrane and on the transmembrane helix-helix interactions required to achieve proper folding and assembly.

Filopodia formation driven by membrane glycoprotein M6a depends on the interaction of its transmembrane domains

Membrane glycoprotein M6a, which belongs to the tetraspan proteolipid protein family, promotes structural plasticity in neurons and cell lines by unknown mechanisms. This glycoprotein is encoded by Gpm6a, a stress-regulated gene. The hippocampus of animals chronically stressed by either psychosocial or physical stressors shows decreased M6a expression. Stressed Gpm6a-null mice develop a claustrophobia-like phenotype. In humans, de novo duplication of GPM6A results in learning/behavioral abnormalities, and two single-nucleotide polymorphisms (SNPs) in the non-coding region are linked to mood disorders. Here, we studied M6a dimerization in neuronal membranes and its functional relevance. We showed that the self-interaction of M6a transmembrane domains (TMDs) might be driving M6a dimerization, which is required to induce filopodia formation. Glycine mutants located in TMD2 and TMD4 of M6a affected its dimerization, thus preventing M6a-induced filopodia formation in neurons. In silico analysis of three non-synonymous SNPs located in the coding region of TMDs suggested that these mutations induce protein instability. Indeed, these SNPs prevented M6a from being functional in neurons, owing to decreased stability, dimerization or improper folding. Interestingly, SNP3 (W141R), which caused endoplasmic reticulum retention, is equivalent to that mutated in PLP1, W161L, which causes demyelinating Pelizaeus-Merzbacher disease. In this work we analyzed the functional contribution of transmembrane domains (TMDs) of the neuronal membrane glycoprotein M6a. We determined that certain glycines present in TMD2 and TMD4 are critical for filopodia induction in neurons. In addition, three nsSNPs located in the coding region of TMD2 and TMD3 of GPM6A impair M6a function by affecting its stability, folding and dimer formation.

Evolutionary-guided de novo structure prediction of self-associated transmembrane helical proteins with near-atomic accuracy

How specific protein associations regulate the function of membrane receptors remains poorly understood. Conformational flexibility currently hinders the structure determination of several classes of membrane receptors and associated oligomers. Here we develop EFDOCK-TM, a general method to predict self-associated transmembrane protein helical (TMH) structures from sequence guided by co-evolutionary information. We show that accurate intermolecular contacts can be identified using a combination of protein sequence covariation and TMH binding surfaces predicted from sequence. When applied to diverse TMH oligomers, including receptors characterized in multiple conformational and functional states, the method reaches unprecedented near-atomic accuracy for most targets. Blind predictions of structurally uncharacterized receptor tyrosine kinase TMH oligomers provide a plausible hypothesis on the molecular mechanisms of disease-associated point mutations and binding surfaces for the rational design of selective inhibitors. The method sets the stage for uncovering novel determinants of molecular recognition and signalling in single-spanning eukaryotic membrane receptors.

Structural determinants of human proton-coupled folate transporter oligomerization: role of GXXXG motifs and identification of oligomeric interfaces at transmembrane domains 3 and 6

The human proton-coupled folate transporter (hPCFT) is expressed in solid tumours and is active at pHs characterizing the tumour microenvironment. Recent attention focused on exploiting hPCFT for targeting solid tumours with novel cytotoxic anti-folates. hPCFT has 12 transmembrane domains (TMDs) and forms homo-oligomers with functional significance. The hPCFT primary sequence includes GXXXG motifs in TMD2 (G(93)XXXG(97)) and TMD4 (G(155)XXXG(159)). To investigate roles of these motifs in hPCFT function, stability and surface expression, we mutated glycine to leucine to generate single or multiple substitution mutants. Only the G93L and G159L mutants preserved substantial [(3)H]methotrexate (Mtx) transport when expressed in hPCFT-null (R1-11) HeLa cells. Transport activity of the glycine-to-leucine mutants correlated with surface hPCFT by surface biotinylation and confocal microscopy with ECFP*-tagged hPCFTs, suggesting a role for GXXXG in hPCFT stability and intracellular trafficking. When co-expressed in R1-11 cells, haemagglutinin-tagged glycine-to-leucine mutants and His10-tagged wild-type (WT) hPCFT co-associated on nickel affinity columns, suggesting that the GXXXG motifs are not directly involved in hPCFT oligomerization. This was substantiated by in situ FRET experiments with co-expressed ECFP*- and YFP-tagged hPCFT. Molecular modelling of dimeric hPCFT structures showed juxtaposed TMDs 2, 3, 4 and 6 as potential structural interfaces between monomers. hPCFT cysteine insertion mutants in TMD3 (Q136C and L137C) and TMD6 (W213C, L214C, L224C, A227C, F228C, F230C and G231C) were expressed in R1-11 cells and cross-linked with 1,6-hexanediyl bismethanethiosulfonate, confirming TMD juxtapositions. Altogether, our results imply that TMDs 3 and 6 provide critical interfaces for formation of hPCFT oligomers, which might be facilitated by the GXXXG motifs in TMD2 and TMD4.

A structural role for the synaptobrevin 2 transmembrane domain in dense-core vesicle fusion pores

Ca(2+)-triggered release of neurotransmitters and hormones depends on soluble N-ethylmaleimide-sensitive factor attachment protein receptors (SNAREs) to drive the fusion of the vesicle and plasma membranes. The formation of the SNARE complex by the vesicle SNARE synaptobrevin 2 (syb2) and the two plasma membrane SNAREs syntaxin (syx) and SNAP-25 draws the two membranes together, but the events that follow membrane juxtaposition, and the ways that SNAREs remodel lipid membranes remain poorly understood. The SNAREs syx and syb2 have transmembrane domains (TMDs) that can exert force directly on the lipid bilayers. The TMD of syx influences fusion pore flux in a manner that suggests it lines the nascent fusion pore through the plasma membrane. The TMD of syb2 traverses the vesicle membrane and is the most likely partner to syx in completing a proteinaceous fusion pore through the vesicle membrane, but the role of this vesicle SNARE in fusion pores has yet to be tested. Here amperometry and conductance measurements were performed to probe the function of the syb2 TMD in fusion pores formed during catecholamine exocytosis in mouse chromaffin cells. Fusion pore flux was sensitive to the size and charge of TMD residues near the N terminus; fusion pore conductance was altered by substitutions at these sites. Unlike syx, the syb2 residues that influence fusion pore permeation fell along two α-helical faces of its TMD, rather than one. These results indicate a role for the syb2 TMD in nascent fusion pores, but in a very different structural arrangement from that of the syx TMD.

The membrane- and soluble-protein helix-helix interactome: similar geometry via different interactions

α Helices are a basic unit of protein secondary structure and therefore the interaction between helices is crucial to understanding tertiary and higher-order folds. Comparing subtle variations in the structural and sequence motifs between membrane and soluble proteins sheds light on the different constraints faced by each environment and elucidates the complex puzzle of membrane protein folding. Here, we demonstrate that membrane and water-soluble helix pairs share a small number of similar folds with various interhelical distances. The composition of the residues that pack at the interface between corresponding motifs shows that hydrophobic residues tend to be more enriched in the water-soluble class of structures and small residues in the transmembrane class. The latter group facilitates packing via sidechain- and backbone-mediated hydrogen bonds within the low-dielectric membrane milieu. The helix-helix interactome space, with its associated sequence preferences and accompanying hydrogen-bonding patterns, should be useful for engineering, prediction, and design of protein structure.

A designed buried salt bridge modulates heterodimerization of a membrane peptide

Specific helix-helix interactions underpin the correct assembly of multipass membrane proteins. Here, we show that a designed buried salt bridge mediates heterodimer formation of model transmembrane helical peptides in a pH-dependent manner. The model peptides bear side chains functionalized with either a carboxylic acid or a primary amine within a hydrophobic segment. The association behavior was monitored by Förster resonance energy transfer, revealing that heterodimer formation is maximized at a pH close to neutrality (pH 6.5), at which each peptide is found in a charged state. In contrast, heterodimerization is disfavored at low and high values of pH, because either the carboxylic acid or the primary amine is present in its neutral state, thus preventing the formation of a salt bridge. These findings provide a blueprint for the design and modulation of protein-protein interactions in membrane proteins.

Backbone dependency further improves side chain prediction efficiency in the Energy-based Conformer Library (bEBL)

Side chain optimization is an integral component of many protein modeling applications. In these applications, the conformational freedom of the side chains is often explored using libraries of discrete, frequently occurring conformations. Because side chain optimization can pose a computationally intensive combinatorial problem, the nature of these conformer libraries is important for ensuring efficiency and accuracy in side chain prediction. We have previously developed an innovative method to create a conformer library with enhanced performance. The Energy-based Library (EBL) was obtained by analyzing the energetic interactions between conformers and a large number of natural protein environments from crystal structures. This process guided the selection of conformers with the highest propensity to fit into spaces that should accommodate a side chain. Because the method requires a large crystallographic data-set, the EBL was created in a backbone-independent fashion. However, it is well established that side chain conformation is strongly dependent on the local backbone geometry, and that backbone-dependent libraries are more efficient in side chain optimization. Here we present the backbone-dependent EBL (bEBL), whose conformers are independently sorted for each populated region of Ramachandran space. The resulting library closely mirrors the local backbone-dependent distribution of side chain conformation. Compared to the EBL, we demonstrate that the bEBL uses fewer conformers to produce similar side chain prediction outcomes, thus further improving performance with respect to the already efficient backbone-independent version of the library.

A Gly-zipper motif mediates homodimerization of the transmembrane domain of the mitochondrial kinase ADCK3

Interactions between α-helices within the hydrophobic environment of lipid bilayers are integral to the folding and function of transmembrane proteins; however, the major forces that mediate these interactions remain debated, and our ability to predict these interactions is still largely untested. We recently demonstrated that the frequent transmembrane association motif GASright, the GxxxG-containing fold of the glycophorin A dimer, is optimal for the formation of extended networks of Cα-H hydrogen bonds, supporting the hypothesis that these bonds are major contributors to association. We also found that optimization of Cα-H hydrogen bonding and interhelical packing is sufficient to computationally predict the structure of known GASright dimers at near atomic level. Here, we demonstrate that this computational method can be used to characterize the structure of a protein not previously known to dimerize, by predicting and validating the transmembrane dimer of ADCK3, a mitochondrial kinase. ADCK3 is involved in the biosynthesis of the redox active lipid, ubiquinone, and human ADCK3 mutations cause a cerebellar ataxia associated with ubiquinone deficiency, but the biochemical functions of ADCK3 remain largely undefined. Our experimental analyses show that the transmembrane helix of ADCK3 oligomerizes, with an interface based on an extended Gly-zipper motif, as predicted by our models. The data provide strong evidence for the hypothesis that optimization of Cα-H hydrogen bonding is an important factor in the association of transmembrane helices. This work also provides a structural foundation for investigating the role of transmembrane association in regulating the biological activity of ADCK3.