Improved prediction for the structure of the dimeric transmembrane domain of glycophorin A obtained through global searching

A Rigorous Framework for Calculating Protein-Protein Binding Affinities in Membranes

Calculating the binding free energy of integral transmembrane (TM) proteins is crucial for understanding the mechanisms by which they recognize one another and reversibly associate. The glycophorin A (GpA) homodimer, composed of two α-helical segments, has long served as a model system for studying TM protein reversible association. The present work establishes a methodological framework for calculating the binding affinity of the GpA homodimer in the heterogeneous environment of a membrane. Our investigation carefully considered a variety of protocols, including the appropriate choice of the force field, rigorous standardization reflecting the experimental conditions, sampling algorithm, anisotropic environment, and collective variables, to accurately describe GpA dimerization via molecular dynamics-based approaches. Specifically, two strategies were explored: (i) an unrestrained potential mean force (PMF) calculation, which merely enhances sampling along the separation of the two binding partners without any restraint, and (ii) a so-called "geometrical route", whereby the α-helices are progressively separated with imposed restraints on their orientational, positional, and conformational degrees of freedom to accelerate convergence. Our simulations reveal that the simplified, unrestrained PMF approach is inadequate for the description of GpA dimerization. Instead, the geometrical route, tailored specifically to GpA in a membrane environment, yields excellent agreement with experimental data within a reasonable computational time. A dimerization free energy of -10.7 kcal/mol is obtained, in fairly good agreement with available experimental data. The geometrical route further helps elucidate how environmental forces drive association before helical interactions stabilize it. Our simulations also brought to light a distinct, long-lived spatial arrangement that potentially serves as an intermediate state during dimer formation. The methodological advances in the generalized geometrical route provide a powerful tool for accurate and efficient binding-affinity calculations of intricate TM protein complexes in inhomogeneous environments.

Structural Determinant of β-Amyloid Formation: From Transmembrane Protein Dimerization to β-Amyloid Aggregates

Most neurodegenerative diseases have the characteristics of protein folding disorders, i.e., they cause lesions to appear in vulnerable regions of the nervous system, corresponding to protein aggregates that progressively spread through the neuronal network as the symptoms progress. Alzheimer's disease is one of these diseases. It is characterized by two types of lesions: neurofibrillary tangles (NFTs) composed of tau proteins and senile plaques, formed essentially of amyloid peptides (Aβ). A combination of factors ranging from genetic mutations to age-related changes in the cellular context converge in this disease to accelerate Aβ deposition. Over the last two decades, numerous studies have attempted to elucidate how structural determinants of its precursor (APP) modify Aβ production, and to understand the processes leading to the formation of different Aβ aggregates, e.g., fibrils and oligomers. The synthesis proposed in this review indicates that the same motifs can control APP function and Aβ production essentially by regulating membrane protein dimerization, and subsequently Aβ aggregation processes. The distinct properties of these motifs and the cellular context regulate the APP conformation to trigger the transition to the amyloid pathology. This concept is critical to better decipher the patterns switching APP protein conformation from physiological to pathological and improve our understanding of the mechanisms underpinning the formation of amyloid fibrils that devastate neuronal functions.

Characterization of interactions within the Igα/Igβ transmembrane domains of the human B-cell receptor provides insights into receptor assembly

The B-cell receptor (BCR), a complex comprised of a membrane-associated immunoglobulin and the Igα/β heterodimer, is one of the most important immune receptors in humans and controls B-cell development, activity, selection, and death. BCR signaling plays key roles in autoimmune diseases and lymphoproliferative disorders, yet, despite the clinical significance of this protein complex, key regions (i.e., the transmembrane domains) have yet to be structurally characterized. The mechanism for BCR signaling also remains unclear and has been variously described by the mutually exclusive cross-linking and dissociation activation models. Common to these models is the significance of local plasma membrane composition, which implies that interactions between BCR transmembrane domains (TMDs) play a role in receptor functionality. Here we used an in vivo assay of TMD oligomerization called GALLEX alongside spectroscopic and computational methods to characterize the structures and interactions of human Igα and Igβ TMDs in detergent micelles and natural membranes. We observed weak self-association of the Igβ TMD and strong self-association of the Igα TMD, which scanning mutagenesis revealed was entirely stabilized by an E–X₁₀–P motif. We also demonstrated strong heterotypic interactions between the Igα and Igβ TMDs both in vitro and in vivo, which scanning mutagenesis and computational models suggest is multiconfigurational but can accommodate distinct interaction sites for self-interactions and heterotypic interactions of the Igα TMD. Taken together, these results demonstrate that the TMDs of the human BCR are sites of strong protein–protein interactions that may direct BCR assembly, endoplasmic reticulum retention, and immune signaling.

Drugging Membrane Protein Interactions

The majority of therapeutics target membrane proteins, accessible on the surface of cells, to alter cellular signaling. Cells use membrane proteins to transduce signals into cells, transport ions and molecules, bind cells to a surface or substrate, and catalyze reactions. Newly devised technologies allow us to drug conventionally "undruggable" regions of membrane proteins, enabling modulation of protein-protein, protein-lipid, and protein-nucleic acid interactions. In this review, we survey the state of the art of high-throughput screening and rational design in drug discovery, and we evaluate the advances in biological understanding and technological capacity that will drive pharmacotherapy forward against unorthodox membrane protein targets.

Prediction of the structures of helical membrane proteins based on a minimum unfavorable contacts approach

An understanding of structure-function relationships of membrane proteins continues to be a challenging problem, owing to the difficulty in obtaining their structures experimentally. This study suggests a method for modeling membrane protein structures that can be used to generate a reliable initial conformation prior to the use of other approaches for sampling conformations. It involves optimizing the orientation of hydrophilic residues so as to minimize unfavorable contacts with the hydrophobic tails of the lipid bilayer. Starting with the optimized initial conformation for three different proteins modeled based on this method, two independent approaches have been used for sampling the conformational space of the proteins. Both approaches are able to predict structures reasonably close to experimental structures, indicating that the initial structure enables the sampling of conformations that are close to the native structure. Possible improvements in the method for making it broadly applicable to helical membrane proteins are discussed.

Mechanism of activation of protein kinase JAK2 by the growth hormone receptor

Signaling from JAK (Janus kinase) protein kinases to STAT (signal transducers and activators of transcription) transcription factors is key to many aspects of biology and medicine, yet the mechanism by which cytokine receptors initiate signaling is enigmatic. We present a complete mechanistic model for activation of receptor-bound JAK2, based on an archetypal cytokine receptor, the growth hormone receptor. For this, we used fluorescence resonance energy transfer to monitor positioning of the JAK2 binding motif in the receptor dimer, substitution of the receptor extracellular domains with Jun zippers to control the position of its transmembrane (TM) helices, atomistic modeling of TM helix movements, and docking of the crystal structures of the JAK2 kinase and its inhibitory pseudokinase domain with an opposing kinase-pseudokinase domain pair. Activation of the receptor dimer induced a separation of its JAK2 binding motifs, driven by a ligand-induced transition from a parallel TM helix pair to a left-handed crossover arrangement. This separation leads to removal of the pseudokinase domain from the kinase domain of the partner JAK2 and pairing of the two kinase domains, facilitating trans-activation. This model may well generalize to other class I cytokine receptors.

Hepatitis C virus RNA replication and virus particle assembly require specific dimerization of the NS4A protein transmembrane domain

Hepatitis C virus (HCV) NS4A is a single-pass transmembrane (TM) protein essential for viral replication and particle assembly. The sequence of the NS4A TM domain is highly conserved, suggesting that it may be important for protein-protein interactions. To test this hypothesis, we measured the potential dimerization of the NS4A TM domain in a well-characterized two-hybrid TM protein interaction system. The NS4A TM domain exhibited a strong homotypic interaction that was comparable in affinity to glycophorin A, a well-studied human blood group antigen that forms TM homodimers. Several mutations predicted to cluster on a common surface of the NS4A TM helix caused significant reductions in dimerization, suggesting that these residues form an interface for NS4A dimerization. Mutations in the NS4A TM domain were further examined in the JFH-1 genotype 2a replicon system; importantly, all mutations that destabilized NS4A dimers also caused defects in RNA replication and/or virus assembly. Computational modeling of NS4A TM interactions suggests a right-handed dimeric interaction of helices with an interface that is consistent with the mutational effects. Furthermore, defects in NS4A oligomerization and virus particle assembly of two mutants were rescued by NS4A A15S, a TM mutation recently identified through forward genetics as a cell culture-adaptive mutation. Together, these data provide the first example of a functionally important TM dimer interface within an HCV nonstructural protein and reveal a fundamental role of the NS4A TM domain in coordinating HCV RNA replication and virus particle assembly.

Inside-out signaling promotes dynamic changes in the carcinoembryonic antigen-related cellular adhesion molecule 1 (CEACAM1) oligomeric state to control its cell adhesion properties

Cell-cell contacts are fundamental to multicellular organisms and are subject to exquisite levels of control. The carcinoembryonic antigen-related cell adhesion molecule 1 (CEACAM1) can engage in both cis-homophilic (parallel) oligomerization and trans-homophilic (anti-parallel) binding. In this study, we establish that the CEACAM1 transmembrane domain has a propensity to form cis-dimers via the transmembrane-embedded (432)GXXXG(436) motif and that this basal state is overcome when activated calmodulin binds to the CEACAM1 cytoplasmic domain. Although mutation of the (432)GXXXG(436) motif reduced CEACAM1 oligomerization, it did not affect surface localization of the receptor or influence CEACAM1-dependent cellular invasion by the pathogenic Neisseria. The mutation did, however, have a striking effect on CEACAM1-dependent cellular aggregation, increasing both the kinetics of cell-cell association and the size of cellular aggregates formed. CEACAM1 association with tyrosine kinase c-Src and tyrosine phosphatases SHP-1 and SHP-2 was not affected by the (432)GXXXG(436) mutation, consistent with their association with the monomeric form of wild type CEACAM1. Collectively, our results establish that a dynamic oligomer-to-monomer shift in surface-expressed CEACAM1 facilitates trans-homophilic binding and downstream effector signaling.

Structural stability of proteins in aqueous and nonpolar environments

A protein folds into its native structure with the α-helix and∕or β-sheet in aqueous solution under the physiological condition. The relative content of these secondary structures largely varies from protein to protein. However, such structural variability is not exhibited in nonaqueous environment. For example, there is a strong trend that alcohol induces a protein to form α-helices, and many of the membrane proteins within the lipid bilayer consists of α-helices. Here we investigate the structural stability of proteins in aqueous and nonpolar environments using our recently developed free-energy function F = (Λ - TS)∕(k(B)T(0)) = Λ∕(k(B)T(0)) - S∕k(B) (T(0) = 298 K and the absolute temperature T is set at T(0)) which is based on statistical thermodynamics. Λ∕(k(B)T(0)) and S∕k(B) are the energetic and entropic components, respectively, and k(B) is Boltzmann's constant. A smaller value of the positive quantity, -S, represents higher efficiency of the backbone and side-chain packing promoted by the entropic effect arising from the translational displacement of solvent molecules or the CH(2), CH(3), and CH groups which constitute nonpolar chains of lipid molecules. As for Λ, in aqueous solution, a transition to a more compact structure of a protein accompanies the break of protein-solvent hydrogen bonds: As the number of donors and acceptors buried without protein intramolecular hydrogen bonding increases, Λ becomes higher. In nonpolar solvent, lower Λ simply implies more intramolecular hydrogen bonds formed. We find the following. The α-helix and β-sheet are advantageous with respect to -S as well as Λ and to be formed as much as possible. In aqueous solution, the solvent-entropy effect on the structural stability is so strong that the close packing of side chains is dominantly important, and the α-helix and β-sheet contents are judiciously adjusted to accomplish it. In nonpolar solvent, the solvent-entropy effect is substantially weaker than in aqueous solution. Λ is crucial and the α-helix is more stable than the β-sheet in terms of Λ, which develops a tendency that α-helices are exclusively chosen. For a membrane protein, α-helices are stabilized as fundamental structural units for the same reason, but their arrangement is performed through the entropic effect mentioned above.

The backbone dynamics of the amyloid precursor protein transmembrane helix provides a rationale for the sequential cleavage mechanism of γ-secretase

The etiology of Alzheimer's disease depends on the relative abundance of different amyloid-β (Aβ) peptide species. These peptides are produced by sequential proteolytic cleavage within the transmembrane helix of the 99 residue C-terminal fragment of the amyloid precursor protein (C99) by the intramembrane protease γ-secretase. Intramembrane proteolysis is thought to require local unfolding of the substrate helix, which has been proposed to be cleaved as a homodimer. Here, we investigated the backbone dynamics of the substrate helix. Amide exchange experiments of monomeric recombinant C99 and of synthetic transmembrane domain peptides reveal that the N-terminal Gly-rich homodimerization domain exchanges much faster than the C-terminal cleavage region. MD simulations corroborate the differential backbone dynamics, indicate a bending motion at a diglycine motif connecting dimerization and cleavage regions, and detect significantly different H-bond stabilities at the initial cleavage sites. Our results are consistent with the following hypotheses about cleavage of the substrate: First, the GlyGly hinge may precisely position the substrate within γ-secretase such that its catalytic center must start proteolysis at the known initial cleavage sites. Second, the ratio of cleavage products formed by subsequent sequential proteolysis could be influenced by differential extents of solvation and by the stabilities of H-bonds at alternate initial sites. Third, the flexibility of the Gly-rich domain may facilitate substrate movement within the enzyme during sequential proteolysis. Fourth, dimerization may affect substrate processing by decreasing the dynamics of the dimerization region and by increasing that of the C-terminal part of the cleavage region.

Roles of hydrophobicity and charge distribution of cationic antimicrobial peptides in peptide-membrane interactions

Cationic antimicrobial peptides (CAPs) occur as important innate immunity agents in many organisms, including humans, and offer a viable alternative to conventional antibiotics, as they physically disrupt the bacterial membranes, leading to membrane lysis and eventually cell death. In this work, we studied the biophysical and microbiological characteristics of designed CAPs varying in hydrophobicity levels and charge distributions by a variety of biophysical and biochemical approaches, including in-tandem atomic force microscopy, attenuated total reflection-FTIR, CD spectroscopy, and SDS-PAGE. Peptide structural properties were correlated with their membrane-disruptive abilities and antimicrobial activities. In bacterial lipid model membranes, a time-dependent increase in aggregated β-strand-type structure in CAPs with relatively high hydrophobicity (such as KKKKKKALFALWLAFLA-NH(2)) was essentially absent in CAPs with lower hydrophobicity (such as KKKKKKAAFAAWAAFAA-NH(2)). Redistribution of positive charges by placing three Lys residues at both termini while maintaining identical sequences minimized self-aggregation above the dimer level. Peptides containing four Leu residues were destructive to mammalian model membranes, whereas those with corresponding Ala residues were not. This finding was mirrored in hemolysis studies in human erythrocytes, where Ala-only peptides displayed virtually no hemolysis up to 320 μM, but the four-Leu peptides induced 40-80% hemolysis at the same concentration range. All peptides studied displayed strong antimicrobial activity against Pseudomonas aeruginosa (minimum inhibitory concentrations of 4-32 μM). The overall findings suggest optimum routes to balancing peptide hydrophobicity and charge distribution that allow efficient penetration and disruption of the bacterial membranes without damage to mammalian (host) membranes.

STRUCTURE PREDICTION OF HELICAL TRANSMEMBRANE PROTEINS AT TWO LENGTH SCALES

As the first step toward a multi-scale, hierarchical computational approach for membrane protein structure prediction, the packing of transmembrane helices was modeled at the residue and atom levels, respectively. For predictions at the residue level, the helix-helix and helix-membrane interactions were described by a set of knowledge-based energy functions. For predictions at the atom level, CHARMM19 force field was used. To facilitate the system to overcome energy barriers, the Wang–Landau method was employed, where a random walk is performed in the energy space with a uniform probability. Native-like structures were predicted at both levels for two model systems, each of which consists of two transmembrane helices. Interestingly, consistent results were obtained from simulations at the residue and atom levels for the same system, strongly suggesting the feasibility of a hierarchical approach for membrane protein structure predictions.

Orientation-specific signalling by thrombopoietin receptor dimers

Ligand binding to the thrombopoietin receptor is thought to stabilize an active receptor dimer that regulates megakaryocyte differentiation and platelet formation, as well as haematopoietic stem cell renewal. By fusing a dimeric coiled coil in all seven possible orientations to the thrombopoietin receptor transmembrane (TM)-cytoplasmic domains, we show that specific biological effects and in vivo phenotypes are imparted by distinct dimeric orientations, which can be visualized by cysteine mutagenesis and crosslinking. Using functional assays and computational searches, we identify one orientation that represents the inactive dimeric state and another similar to a physiologically activated receptor. Several other dimeric orientations are identified that induce proliferation and in vivo myeloproliferative and myelodysplastic disorders, indicating the receptor can signal from several dimeric interfaces. The set of dimeric thrombopoietin receptors with different TM orientations may offer new insights into the activation of distinct signalling pathways by a single receptor and suggests that subtle differences in cytokine receptor dimerization provide a new layer of signalling regulation that is relevant for disease.

Thrombopoietin receptor activation: transmembrane helix dimerization, rotation, and allosteric modulation

We report how rotational variations in transmembrane (TM) helix interactions participate in the activity states of the thrombopoietin receptor (TpoR), a type 1 cytokine receptor that controls the production of blood platelets. We also explore the mechanism of small-molecule agonists that do not mimic the natural ligand. We show, by a combination of cysteine cross-linking, alanine-scanning mutagenesis, and computational simulations, that the TpoR TM dimerizes strongly and can adopt 3 different stable, rotationally related conformations, which may correspond to specific states of the full-length receptor (active, inactive, and partially active). Thus, our data suggest that signaling and inactive states of the receptor are related by receptor subunit rotations, rather than a simple monomer-dimer transition. Moreover, results from experiments with and without agonists in vitro and in cells allow us to propose a novel allosteric mechanism of action for a class of small molecules, in which they activate TpoR by binding to the TM region and by exploiting the rotational states of the dimeric receptor. Overall, our results support the emerging view of the participation of mutual rotations of the TM domains in cytokine receptor activation.

Molecular determinants and thermodynamics of the amyloid precursor protein transmembrane domain implicated in Alzheimer's disease

The deposition of toxic amyloid-β (Aβ) peptide aggregates in the brain is a hallmark of Alzheimer's disease. The intramembrane proteolysis by γ-secretase of the amyloid precursor protein β-carboxy-terminal fragment (APP-βCTF) constitutes the final step in the production of Aβ peptides. Mounting evidence suggests that APP-βCTF is a transmembrane domain (TMD) dimer, and that dimerization might modulate the production of Aβ species that are prone to aggregation and are therefore most toxic. We combined experimental and computational approaches to study the molecular determinants and thermodynamics of APP-βCTF dimerization, and we produced a unifying structural model that reconciles much of the published data. Using a cell assay that exploits a dimerization-dependent activator of transcription, we identified specific dimerization-affecting mutations located mostly at the N-terminus of the TMD of APP-βCTF. The ability of selected mutants to affect the dimerization of full-length APP-βCTF was confirmed by fluorescence resonance energy transfer experiments. Free-energy estimates of the wild type and mutants of the TMD of APP-βCTF derived from enhanced molecular dynamics simulations showed that the dimeric state is composed of different arrangements, in which either (709)GXXXA(713) or (700)GXXXG(704)GXXXG(708) interaction motifs can engage in symmetric or asymmetric associations. Mutations along the TMD of APP-βCTF were found to modulate the relative free energy of the dimeric configurations and to differently affect the distribution of interfaces within the dimeric state. This observation might have important biological implications, since dimers with a different arrangement of the transmembrane helices are likely to be recognized differently by γ-secretase and to lead to a variation in Aβ levels.

The dimeric transmembrane domain of prolyl dipeptidase DPP-IV contributes to its quaternary structure and enzymatic activities

Dipeptidyl peptidase IV (DPP-IV) is a drug target in the treatment of human type II diabetes. It is a type II membrane protein with a single transmembrane domain (TMD) anchoring the extracellular catalytic domain to the membrane. DPP-IV is active as a dimer, with two dimer interacting surfaces located extracellularly. In this study, we demonstrate that the TM of DPP-IV promotes DPP-IV dimerization and rescues monomeric DPP-IV mutants into partial dimers, which is specific and irreplaceable by TMs of other type II membrane proteins. By bioluminescence resonance energy transfer (BRET) and peptide electrophoresis, we found that the TM domain of DPP-IV is dimerized in mammalian cells and in vitro. The TM dimer interaction is very stable, based on our results with TM site-directed mutagenesis. None of the mutations, including the introduction of two prolines, resulted in their complete disruption to monomers. However, these TM proline mutations result in a significant reduction of DPP-IV enzymatic activity, comparable to what is found with mutations near the active site. A systematic analysis of TM structures deposited in the Protein Data Bank showed that prolines in the TM generally produce much bigger kinking angles than occur in nonproline-containing TMs. Thus, the proline-dependent reduction in enzyme activity may result from propagated conformational changes from the TM to the extracellular active site. Our results demonstrate that TM dimerization and conformation contribute significantly to the structure and activity of DPP-IV. Optimal enzymatic activity of DPP-IV requires an optimal interaction of all three dimer interfaces, including its TM.

Lipid-modulated sequence-specific association of glycophorin A in membranes

Protein association in lipid membranes is a complex process with thermodynamics directed by a multitude of different factors. Amino-acid sequence is a molecular parameter that affects dimerization as shown by limited directed mutations along the transmembrane domains. Membrane-mediated interactions are also important although details of such contributions remain largely unclear. In this study, we probe directly the free energy of association of Glycophorin A by means of extensive parallel Monte Carlo simulations with recently developed methods and a model that accounts for sequence-specificity while representing lipid membranes faithfully. We find that lipid-induced interactions are significant both at short and intermediate separations. The ability of molecules to tilt in a specific hydrophobic environment extends their accessible interfaces, leading to intermittent contacts during protein recognition. The dimer with the lowest free energy is largely determined by the favorable lipid-induced attractive interactions at the closest distance. Finally, the coarse-grained model employed herein, together with the extensive sampling performed, provides estimates of the free energy of association that are in excellent agreement with existing data.

Physical basis behind achondroplasia, the most common form of human dwarfism

Fibroblast growth factor receptor 3 (FGFR3) is a receptor tyrosine kinase that plays an important role in long bone development. The G380R mutation in FGFR3 transmembrane domain is known as the genetic cause for achondroplasia, the most common form of human dwarfism. Despite many studies, there is no consensus about the exact mechanism underlying the pathology. To gain further understanding into the physical basis behind the disorder, here we measure the activation of wild-type and mutant FGFR3 in mammalian cells using Western blots, and we analyze the activation within the frame of a physical-chemical model describing dimerization, ligand binding, and phosphorylation probabilities within the dimers. The data analysis presented here suggests that the mutation does not increase FGFR3 dimerization, as proposed previously. Instead, FGFR3 activity in achondroplasia is increased due to increased probability for phosphorylation of the unliganded mutant dimers. This finding has implications for the design of targeted molecular treatments for achondroplasia.

Structure elucidation of dimeric transmembrane domains of bitopic proteins

The interaction between transmembrane helices is of great interest because it directly determines biological activity of a membrane protein. Either destroying or enhancing such interactions can result in many diseases related to dysfunction of different tissues in human body. One much studied form of membrane proteins known as bitopic protein is a dimer containing two membrane-spanning helices associating laterally. Establishing structure-function relationship as well as rational design of new types of drugs targeting membrane proteins requires precise structural information about this class of objects. At present time, to investigate spatial structure and internal dynamics of such transmembrane helical dimers, several strategies were developed based mainly on a combination of NMR spectroscopy, optical spectroscopy, protein engineering and molecular modeling. These approaches were successfully applied to homo- and heterodimeric transmembrane fragments of several bitopic proteins, which play important roles in normal and in pathological conditions of human organism.

Helical membrane peptides to modulate cell function

In recent years there has been an abundance of research into the potential of helical peptides to influence cell function. These peptides have been used to achieve a variety of different outcomes from cell repair to cell death, depending upon the peptide sequence and the nature of its interactions with cell membranes and membrane proteins. In this critical review, we summarise several mechanisms by which helical peptides, acting as either transporters, inhibitors, agonists or antibiotics, can have significant effects on cell membranes and can radically affect the internal mechanisms of the cell. The various approaches to peptide design are discussed, including the role of naturally-occurring proteins in the design of these helical peptides and current breakthroughs in the use of non-natural (and therefore more stable) peptide scaffolds. Most importantly, the current successful applications of these peptides, and their potential uses in the field of medicine, are reviewed (131 references).

Contribution of charged and polar residues for the formation of the E1-E2 heterodimer from Hepatitis C Virus

The transmembrane domains of the envelope glycoprotein E1 and E2 have crucial multifunctional roles in the biogenesis of hepatitis C virus. We have performed molecular dynamics simulations to investigate a structural model of the transmembrane segments of the E1-E2 heterodimer. The simulations support the key role of the Lys370-Asp728 ion pair for mediating the E1-E2 heterodimerization. In comparison to these two residues, the simulation results also reveal the differential effect of the conserved Arg730 residue that has been observed in experimental studies. Furthermore, we discovered the formation of inter-helical hydrogen bonds via Asn367 that stabilize dimer formation. Simulations of single and double mutants further demonstrate the importance of the ion-pair and polar interactions between the interacting helix monomers. The conformation of the E1 fragment in the simulation of the E1-E2 heterodimer is in close agreement with an NMR structure of the E1 transmembrane segment. The proposed model of the E1-E2 heterodimer supports the postulated cooperative insertion of both helices by the translocon complex into the bilayer.

Helix-helix interaction patterns in membrane proteins

Membrane-spanning α-helices represent major sites of protein-protein interaction in membrane protein oligomerization and folding. As such, these interactions may be of exquisite specificity. Specificity often rests on a complex interplay of different types of residues forming the helix-helix interfaces via dense packing and different non-covalent forces, including van der Waal’s forces, hydrogen bonding, charge-charge interactions, and aromatic interactions. These interfaces often contain complex residue motifs where the contribution of constituent amino acids depends on the context of the surrounding sequence. Moreover, transmembrane helix-helix interactions are increasingly recognized as being dynamic and dependent on the functional state of a given protein.

Structure of the lethal phage pinhole

Perhaps the simplest of biological timing systems, bacteriophage holins accumulate during the phage morphogenesis period and then trigger to permeabilize the cytoplasmic membrane with lethal holes; thus, terminating the infection cycle. Canonical holins form very large holes that allow nonspecific release of fully-folded proteins, but a recently discovered class of holins, the pinholins, make much smaller holes, or pinholes, that serve only to depolarize the membrane. Here, we interrogate the structure of the prototype pinholin by negative-stain transmission electron-microscopy, cysteine-accessibility, and chemical cross-linking, as well as by computational approaches. Together, the results suggest that the pinholin forms symmetric heptameric structures with the hydrophilic surface of one transmembrane domain lining the surface of a central channel approximately 15 A in diameter. The structural model also suggests a rationale for the prehole state of the pinholin, the persistence of which defines the duration of the viral latent period, and for the sensitivity of the holin timing system to the energized state of the membrane.

Model generation of viral channel forming 2B protein bundles from polio and coxsackie viruses

2B is a 99 amino acid membrane protein encoded by enteroviruses such as polio and coxsackie viruses with two transmembrane domains. The protein is found to make membranes of infected cells permeable. Using a computational approach which positions the models and assesses stability by molecular dynamics (MD) simulations a putative tetrameric bundle model of 2B is generated. The bundles show a pore lining motif of three lysines followed by a serine. The bundle is discussed in terms of different possible orientations of the helices in the membrane and the consequences this has on the in vivo activity of 2B.

Specificity of helix packing in transmembrane dimer of the cell death factor BNIP3: a molecular modeling study

BNIP3 is a mitochondrial 19-kDa proapoptotic protein, a member of the Bcl-2 family. It has a single COOH-terminal transmembrane (TM) alpha-helical domain, which is required for membrane targeting, proapoptotic activity, hetero- and homo-dimerization in membrane. The role and the molecular details of association of TM helices of BNIP3 are yet to be established. Here, we present a molecular modeling study of helix interactions in its membrane domain. The approach combines Monte Carlo conformational search in an implicit hydrophobic slab followed by molecular dynamics simulations in a hydrated full-atom lipid bilayer. The former technique was used for exhaustive sampling of the peptides' conformational space and for generation of putative "native-like" structures of the dimer. The latter ones were taken as realistic starting points to assess stability and dynamic behavior of the complex in explicit lipid-water surrounding. As a result, several groups of tightly packed right-handed structures of the dimer were proposed. They have almost similar helix-helix interface, which includes the motif A(176)xxxG(180)xxxG(184) and agrees well with previous mutagenesis data and preliminary NMR analysis. Molecular dynamics simulations of these structures reveal perfect adaptation of most of them to heterogeneous membrane environment. A remarkable feature of the predicted dimeric structures is the occurrence of a cluster of H-bonded histidine 173 and serines 168 and 172 on the helix interface, near the N-terminus. Because of specific polar interactions between the monomers, this part of the dimer has no such dense packing as the C-terminal one, thus allowing penetration of water from the extramembrane side into the membrane interior. We propose that the ionization state of His(173) can mediate structural and dynamic properties of the dimer. This, in turn, may be related to pH-dependent proapoptotic activity of BNIP3, which is triggering on by acidosis appearing under hypoxic conditions.

An activating mutation in the CSF3R gene induces a hereditary chronic neutrophilia

We identify an autosomal mutation in the CSF3R gene in a family with a chronic neutrophilia. This T617N mutation energetically favors dimerization of the granulocyte colony-stimulating factor (G-CSF) receptor transmembrane domain, and thus, strongly promotes constitutive activation of the receptor and hypersensitivity to G-CSF for proliferation and differentiation, which ultimately leads to chronic neutrophilia. Mutant hematopoietic stem cells yield a myeloproliferative-like disorder in xenotransplantation and syngenic mouse bone marrow engraftment assays. The survey of 12 affected individuals during three generations indicates that only one patient had a myelodysplastic syndrome. Our data thus indicate that mutations in the CSF3R gene can be responsible for hereditary neutrophilia mimicking a myeloproliferative disorder.

Artificial transmembrane oncoproteins smaller than the bovine papillomavirus E5 protein redefine sequence requirements for activation of the platelet-derived growth factor beta receptor

The bovine papillomavirus E5 protein (BPV E5) is a 44-amino-acid homodimeric transmembrane protein that binds directly to the transmembrane domain of the platelet-derived growth factor (PDGF) beta receptor and induces ligand-independent receptor activation. Three specific features of BPV E5 are considered important for its ability to activate the PDGF beta receptor and transform mouse fibroblasts: a pair of C-terminal cysteines, a transmembrane glutamine, and a juxtamembrane aspartic acid. By using a new genetic technique to screen libraries expressing artificial transmembrane proteins for activators of the PDGF beta receptor, we isolated much smaller proteins, from 32 to 36 residues, that lack all three of these features yet still dimerize noncovalently, specifically activate the PDGF beta receptor via its transmembrane domain, and transform cells efficiently. The primary amino acid sequence of BPV E5 is virtually unrecognizable in some of these proteins, which share as few as seven consecutive amino acids with the viral protein. Thus, small artificial proteins that bear little resemblance to a viral oncoprotein can nevertheless productively interact with the same cellular target. We speculate that similar cellular proteins may exist but have been overlooked due to their small size and hydrophobicity.

Interaction and conformational dynamics of membrane-spanning protein helices

Within 1 or 2 decades, the reputation of membrane-spanning alpha-helices has changed dramatically. Once mostly regarded as dull membrane anchors, transmembrane domains are now recognized as major instigators of protein-protein interaction. These interactions may be of exquisite specificity in mediating assembly of stable membrane protein complexes from cognate subunits. Further, they can be reversible and regulatable by external factors to allow for dynamic changes of protein conformation in biological function. Finally, these helices are increasingly regarded as dynamic domains. These domains can move relative to each other in different functional protein conformations. In addition, small-scale backbone fluctuations may affect their function and their impact on surrounding lipid shells. Elucidating the ways by which these intricate structural features are encoded by the amino acid sequences will be a fascinating subject of research for years to come.

Molecular dynamics studies of the transmembrane domain of gp41 from HIV-1

Helix-helix interactions in the putative three-helix bundle formation of the gp41 transmembrane (TM) domain may contribute to the process of virus-cell membrane fusion in HIV-1 infection. In this study, molecular dynamics is used to analyze and compare the conformations of monomeric and trimeric forms of the TM domain in various solvent systems over the course of 4 to 23-ns simulations. The trimeric bundles of the TM domain were stable as helices and remained associated in a hydrated POPE lipid bilayer for the duration of the 23-ns simulation. Several stable inter-chain hydrogen bonds, mostly among the three deprotonated arginine residues located at the center of each of the three TM domains, formed in a right-handed bundle embedded in the lipid bilayer. No such bonds were observed when the bundle was left-handed or when the central arginine residue in each of the three TM helices was replaced with isoleucine (R_I mutant), suggesting that the central arginine residues may play an essential role in maintaining the integrity of the three-helix bundle. These observations suggest that formation of the three-helix bundle of the TM domain may play a role in the trimerization of gp41, thought to occur during the virus-cell membrane fusion process.

Multiple approaches converge on the structure of the integrin alphaIIb/beta3 transmembrane heterodimer

Integrins link the cytoskeleton to the extracellular matrix and regulate key signaling events that coordinate cellular processes such as secretion, migration, and proliferation. A single integrin molecule can exist in a resting state that does not bind extracellular ligands or in an active state that can engage ligands and form large signaling complexes. Activation signals are transduced between the cytosolic region and the extracellular region by a binary on/off switch in the integrin's transmembrane (TM) domain; the integrin's alpha and beta subunits each have a single TM helix that forms an alpha/beta heterodimer in the resting state, and the TM heterodimer separates to transduce an activation signal across the membrane. In this article, two methods used to generate models of the TM heterodimer, both converging on the same structure, are described. The first model was generated by a Monte Carlo algorithm that selected conformations based on their agreement with published experimental mutagenesis results. The second model was generated by threading the integrin's sequence onto TM helix dimers parsed from the Protein Data Bank and by selecting conformations based on their agreement with published experimental cysteine crosslinking results. The two models have similar structures; however, they differ markedly from some previously published models. To distinguish conformations that reflect the native integrin, we compared the Monte Carlo model, the threaded model, and four published models with experimental mutagenesis and cysteine crosslinking results. The models presented here had high correlation coefficients when compared with experimental findings, and they are in excellent agreement, both in terms of accuracy and in terms of precision, with a recent NMR structure. These results demonstrate that multiple approaches converged on the same structure of the resting integrin's TM heterodimer, and this conformation likely reflects the integrin's native structure.

The assembly motif of a bacterial small multidrug resistance protein

Multidrug transporters such as the small multidrug resistance (SMR) family of bacterial integral membrane proteins are capable of conferring clinically significant resistance to a variety of common therapeutics. As antiporter proteins of approximately 100 amino acids, SMRs must self-assemble into homo-oligomeric structures for efflux of drug molecules. Oligomerization centered at transmembrane helix four (TM4) has been implicated in SMR assembly, but the full complement of residues required to mediate its self-interaction remains to be characterized. Here, we use Hsmr, the 110-residue SMR family member of the archaebacterium Halobacterium salinarum, to determine the TM4 residue motif required to mediate drug resistance and SMR self-association. Twelve single point mutants that scan the central portion of the TM4 helix (residues 85-104) were constructed and were tested for their ability to confer resistance to the cytotoxic compound ethidium bromide. Six residues were found to be individually essential for drug resistance activity (Gly(90), Leu(91), Leu(93), Ile(94), Gly(97), and Val(98)), defining a minimum activity motif of (90)GLXLIXXGV(98) within TM4. When the propensity of these mutants to dimerize on SDS-PAGE was examined, replacements of all but Ile resulted in approximately 2-fold reduction of dimerization versus the wild-type antiporter. Our work defines a minimum activity motif of (90)GLXLIXXGV(98) within TM4 and suggests that this sequence mediates TM4-based SMR dimerization along a single helix surface, stabilized by a small residue heptad repeat sequence. These TM4-TM4 interactions likely constitute the highest affinity locus for disruption of SMR function by directly targeting its self-assembly mechanism.

Structural and functional importance of transmembrane domain 3 (TM3) in the aspartate:alanine antiporter AspT: topology and function of the residues of TM3 and oligomerization of AspT

AspT, the aspartate:alanine antiporter of Tetragenococcus halophilus, a membrane protein of 543 amino acids with 10 putative transmembrane (TM) helices, is the prototype of the aspartate:alanine exchanger (AAE) family of transporters. Because TM3 (isoleucine 64 to methionine 85) has many amino acid residues that are conserved among members of the AAE family and because TM3 contains two charged residues and four polar residues, it is thought to be located near (or to form part of) the substrate translocation pathway that includes the binding site for the substrates. To elucidate the role of TM3 in the transport process, we carried out cysteine-scanning mutagenesis. The substitutions of tyrosine 75 and serine 84 had the strongest inhibitory effects on transport (initial rates of l-aspartate transport were below 15% of the rate for cysteine-less AspT). Considerable but less-marked effects were observed upon the replacement of methionine 70, phenylalanine 71, glycine 74, arginine 76, serine 83, and methionine 85 (initial rates between 15% and 30% of the rate for cysteine-less AspT). Introduced cysteine residues at the cytoplasmic half of TM3 could be labeled with Oregon green maleimide (OGM), whereas cysteines close to the periplasmic half (residues 64 to 75) were not labeled. These results suggest that TM3 has a hydrophobic core on the periplasmic half and that hydrophilic residues on the cytoplasmic half of TM3 participate in the formation of an aqueous cavity in membranes. Furthermore, the presence of l-aspartate protected the cysteine introduced at glycine 62 against a reaction with OGM. In contrast, l-aspartate stimulated the reactivity of the cysteine introduced at proline 79 with OGM. These results demonstrate that TM3 undergoes l-aspartate-induced conformational alterations. In addition, nonreducing sodium dodecyl sulfate-polyacrylamide gel electrophoresis analyses and a glutaraldehyde cross-linking assay suggest that functional AspT forms homo-oligomers as a functional unit.

Self-association of transmembrane domain 2 (TM2), but not TM1, in carnitine palmitoyltransferase 1A: role of GXXXG(A) motifs

Carnitine palmitoyltransferase 1 (CPT1) controls the rate of entry of long-chain fatty acids into the mitochondrial matrix for beta-oxidation and has been reported to exist as an oligomer. We have investigated the in vivo oligomerization of full-length rat CPT1A (rCPT1A) along with those of the N-terminal truncation/deletion mutants Delta(1-82), Delta(1-18), and Delta(19-30) expressed in yeast mitochondria. The data indicate that in liver mitochondria in vivo CPT1A exists as a hexamer but that during preparation and storage of mitochondria the order of oligomerization is rapidly reduced to the trimer, such that a mixture of hexamer and trimer is observed in isolated mitochondria in vitro. Mutants bearing deletions of different segments of the N terminus (including the more N-terminal of the two transmembrane domains) have the same pattern of oligomerization when expressed in yeast mitochondria. The self-association of the individual rCPT1A transmembrane (TM) domains (TM1, TM2) was also studied using the TOXCAT assay (which measures TM self-association in the Escherichia coli inner membrane). There was minimal self-association of the sequence corresponding to TM1 but significant self-association of TM2 in TOXCAT. Chemical cross-linking and analytical ultracentrifugation of a TM2-derived synthetic peptide showed oligomerization with a similar trimer/hexamer equilibrium to that observed for native rCPT1A in isolated mitochondria. Therefore, there was a correlation between the oligomerization behavior of TM2 peptide and that of the full-length protein. In silico molecular modeling of rCPT1A TM2 highlighted the favorable orientation of GXXXG and GXXXA motifs in the formation of the TM2 hexamer.

Amyloidogenic processing but not amyloid precursor protein (APP) intracellular C-terminal domain production requires a precisely oriented APP dimer assembled by transmembrane GXXXG motifs

The beta-amyloid peptide (Abeta) is the major constituent of the amyloid core of senile plaques found in the brain of patients with Alzheimer disease. Abeta is produced by the sequential cleavage of the amyloid precursor protein (APP) by beta- and gamma-secretases. Cleavage of APP by gamma-secretase also generates the APP intracellular C-terminal domain (AICD) peptide, which might be involved in regulation of gene transcription. APP contains three Gly-XXX-Gly (GXXXG) motifs in its juxtamembrane and transmembrane (TM) regions. Such motifs are known to promote dimerization via close apposition of TM sequences. We demonstrate that pairwise replacement of glycines by leucines or isoleucines, but not alanines, in a GXXXG motif led to a drastic reduction of Abeta40 and Abeta42 secretion. beta-Cleavage of mutant APP was not inhibited, and reduction of Abeta secretion resulted from inhibition of gamma-cleavage. It was anticipated that decreased gamma-cleavage of mutant APP would result from inhibition of its dimerization. Surprisingly, mutations of the GXXXG motif actually enhanced dimerization of the APP C-terminal fragments, possibly via a different TM alpha-helical interface. Increased dimerization of the TM APP C-terminal domain did not affect AICD production.

Transmembrane protein models based on high-throughput molecular dynamics simulations with experimental constraints

Elucidating the structure of transmembrane proteins domains with high-resolution methods is a difficult and sometimes impossible task. Here, we explain the method of combining a limited amount of experimental data with automated high-throughput molecular dynamics (MD) simulations of alpha-helical transmembrane bundles in an explicit lipid bilayer/water environment. The procedure uses a systematic conformational search of the helix rotation with experimentally constrained MDs simulations. The experimentally determined helix tilt and rotational angle of a labeled residue with site-specific infrared dichroism allows us to select a unique high-resolution model from a number of possible energy minima encountered in the systematic conformational search.

Phenylalanine promotes interaction of transmembrane domains via GxxxG motifs

Interactions of transmembrane helices play a crucial role in the folding and oligomerisation of integral membrane proteins. In order to uncover novel sequence motifs mediating these interactions, we randomised one face of a transmembrane helix with a set of non-polar or moderately polar amino acids. Those sequences capable of self-interaction upon integration into bacterial inner membranes were selected by means of the ToxR/POSSYCCAT system. A comparison between low/medium-affinity and high-affinity sequences reveals that high-affinity sequences are strongly enriched in phenylalanine residues that are frequently observed at the -3 position of GxxxG motifs, thus yielding FxxGxxxG motifs. Mutation of Phe or GxxxG in selected sequences significantly reduces self-interaction of the transmembrane domains without affecting their efficiency of membrane integration. Conversely, grafting FxxGxxxG onto unrelated transmembrane domains strongly enhances their interaction. Further, we find that FxxGxxxG is significantly over-represented in transmembrane domains of bitopic membrane proteins. The same motif contributes to self-interaction of the vesicular stomatitis virus G protein transmembrane domain. We conclude that Phe stabilises membrane-spanning GxxxG motifs. This is one example of how the role of certain side-chains in helix-helix interfaces is modulated by sequence context.

Molecular dynamics simulation approach for the prediction of transmembrane helix–helix heterodimers assembly

Computational methods are useful to identify favorable structures of transmembrane (TM) helix oligomers when experimental data are not available or when they cannot help to interpret helix-helix association. We report here a global search method using molecular dynamics (MD) simulations to predict the structures of transmembrane homo and heterodimers. The present approach is based only on sequence information without any experimental data and is first applied to glycophorin A to validate the protocol and to the HER2-HER3 heterodimer receptor. The method successfully reproduces the experimental structures of the TM domain of glycophorin A (GpA(TM)) with a root mean square deviation of 1.5 A. The search protocol identifies three energetically stable models of the TM domain of HER2-HER3 receptor with favorable helix-helix arrangement, including right-handed and left-handed coiled-coils. The predicted TM structures exhibit the GxxxG-like motif at the dimer interface which is presumed to drive receptor oligomerization. We demonstrate that native structures of TM domain can be predicted without quantitative experimental data. This search protocol could help to predict structures of the TM domain of HER heterodimer family.

Membrane interactions of the hydrophobic segment of diacylglycerol kinase epsilon

Diacylglycerol kinase epsilon (DGKepsilon) is unique among mammalian DGK isoforms in having a segment of hydrophobic amino acids as a putative membrane anchor. To model the conformation, and stoichiometry of this segment in membrane-mimetic environments, we have prepared a peptide corresponding to this hydrophobic segment of DGKepsilon of sequence KKKKLILWTLCSVLLPVFITFWKKKKK-NH(2). Flanking Lys residues mimic the natural setting of this peptide in DGKepsilon, while facilitating peptide synthesis and characterization. Circular dichroism and fluorescence spectroscopic analysis demonstrated that the peptide has increased helical content and significant blue shifts in the presence of anionic--but not zwitterionic--bilayer membranes. When labeled with fluorophores that can undergo fluorescence resonance energy transfer, the peptide was found to dimerize--a result also observed from migration rates on SDS-PAGE gels under both reducing and non-reducing disulfide bridge conditions. The peptide was shown to preferentially interact with cholesterol in lipid films comprised of homogeneous mixtures of cholesterol and phosphatidylcholine, yet the presence of cholesterol in hydrated vesicle bilayers decreases its helical content. The peptide was also able to inhibit the activity of DGKepsilon protein in vitro. Our overall findings suggest that the peptide ultimately cannot leave the bulk water for attachment/insertion into the outer leaflet of an erythrocyte-like bilayer, yet its core sequence is sufficiently hydrophobic to insert into membrane core regions when membrane attachment is promoted by electrostatic attraction to anionic lipid head groups of the inner leaflet of an erythrocyte-like bilayer.

Mutagenesis data in the automated prediction of transmembrane helix dimers

We present a molecular modeling protocol that selects modeled protein structures based on experimental mutagenesis results. The computed effect of a point mutation should be consistent with its experimental effect for correct models; mutations that do not affect protein stability and function should not affect the computed energy of a correct model while destabilizing mutations should have unfavorable computed energies. On the other hand, an incorrect model will likely display computed energies that are inconsistent with experimental results. We added terms to our energy function which penalize models that are inconsistent with experimental results. This creates a selective advantage for models that are consistent with experimental results in the Monte Carlo simulated annealing protocol we use to search conformational space. We calibrated our protocol to predict the structure of transmembrane helix dimers using glycophorin A as a model system. Inclusion of mutational data in this protocol compensates for the limitations of our force field and the limitations of our conformational search. We demonstrate an application of this structure prediction protocol by modeling the transmembrane region of the BNIP3 apoptosis factor.

Ab initio molecular-replacement phasing for symmetric helical membrane proteins

Obtaining phases for X-ray diffraction data can be a rate-limiting step in structure determination. Taking advantage of constraints specific to membrane proteins, an ab initio molecular-replacement method has been developed for phasing X-ray diffraction data for symmetric helical membrane proteins without prior knowledge of their structure or heavy-atom derivatives. The described method is based on generating all possible orientations of idealized transmembrane helices and using each model in a molecular-replacement search. The number of models is significantly reduced by taking advantage of geometrical and structural restraints specific to membrane proteins. The top molecular-replacement results are evaluated based on noncrystallographic symmetry (NCS) map correlation, OMIT map correlation and R(free) value after refinement of a polyalanine model. The feasibility of this approach is illustrated by phasing the mechanosensitive channel of large conductance (MscL) with only 4 A diffraction data. No prior structural knowledge was used other than the number of transmembrane helices. The search produced the correct spatial organization and the position in the asymmetric unit of all transmembrane helices of MscL. The resulting electron-density maps were of sufficient quality to automatically build all helical segments of MscL including the cytoplasmic domain. The method does not require high-resolution diffraction data and can be used to obtain phases for symmetrical helical membrane proteins with one or two helices per monomer.

Membrane assembly of simple helix homo-oligomers studied via molecular dynamics simulations

The assembly of simple transmembrane helix homo-oligomers is studied by combining a generalized Born implicit membrane model with replica exchange molecular dynamics simulations to sample the conformational space of various oligomerization states and the native oligomeric conformation. Our approach is applied to predict the structures of transmembrane helices of three proteins--glycophorin A, the M2 proton channel, and phospholamban--using only peptide sequence and the native oligomerization state information. In every case, the methodology reproduces native conformations that are in good agreement with available experimental structural data. Thus, our method should be useful in the prediction of native structures of transmembrane domains of other peptides. When we ignore the experimental constraint on the native oligomerization state and attempt de novo prediction of the structure and oligomerization state based only on sequence and simple energetic considerations, we identify the pentamer as the most stable oligomer for phospholamban. However, for the glycophorin A and the M2 proton channels, we tend to predict higher oligomers as more stable. Our studies demonstrate that reliable predictions of the structure of transmembrane helical oligomers can be achieved when the observed oligomerization state is imposed as a constraint, but that further efforts are needed for the de novo prediction of both structure and oligomeric state.

Systematic molecular dynamics searching in a lipid bilayer: Application to the glycophorin A and oncogenic ErbB-2 transmembrane domains

Molecular dynamics (MD) simulations of proteins in a lipid bilayer environment are usually undertaken with one or a few starting structures. Here we report a search protocol for systematically exploring the possible interactions in helical bundle transmembrane proteins, a frequently occurring structural motif. The search protocol correctly identifies the experimentally known structure of the dimeric human glycophorin A transmembrane domain as the lowest energy structure among five different models without any prior assumptions, whilst an identical in vacuo search fails to identify the correct structure. The lowest energy structure from the search in a lipid bilayer has a root mean square deviation of 1.1A to the experimental structure. We have applied the same search protocol to the unknown transmembrane structure of the oncogenic mutant ErbB-2 protein, a member of the family of epidermal growth factor receptors. Resulting structures show the role of glutamic acid hydrogen bonding and close helical packing. Water molecules may also play a key role in stabilisation of the transmembrane helix association.