Transmembrane helix packing of ErbB/Neu receptor in membrane environment: a molecular dynamics study

Dimerization or oligomerization of the ErbB/Neu receptors are necessary but not sufficient for initiation of receptor signaling. The two intracellular domains must be properly oriented for the juxtaposition of the kinase domains allowing trans-phosphorylation. This suggests that the transmembrane (TM) domain acts as a guide for defining the proper orientation of the intracellular domains. Two structural models, with the two helices either in left-handed or in right-handed coiling have been proposed as the TM domain structure of the active receptor. Because experimental data do not distinguish clearly helix-helix packing, molecular dynamics (MD) simulations are used to investigate the energetic factors that drive Neu TM-TM interactions of the wild and the oncogenic receptor (Val664/Glu mutation) in DMPC or in POPC environments. MD results indicate that helix-lipid interactions in the bilayer core are extremely similar in the two environments and raise the role of the juxtamembrane residues in helix insertion and helix-helix packing. The TM domain shows a greater propensity to adopt a left-handed structure in DMPC, with helices in optimal position for strong inter-helical Hbonds induced by the Glu mutation. In POPC, the right-handed structure is preferentially formed with the participation of water in inter-helical Hbonds. The two structural arrangements of the Neu(TM) helices both with GG4 residue motif in close contact at the interface are permissible in the membrane environment. According to the hypothesis of a monomer-dimer equilibrium of the proteins it is likely that the bilayer imposes structural constraints that favor dimerization-competent structure responsible of the proper topology necessary for receptor activation.     

Molecular dynamics (MD) investigations of preformed structures of the transmembrane domain of the oncogenic Neu receptor dimer in a DMPC bilayer

The critical Val/Glu mutation in the membrane spanning domain of the rat Neu receptor confers the ability for ligand-independent signaling and leads to increased dimerization and transforming ability. There is evidence that the two transmembrane interacting helices play a role in receptor activation by imposing orientation constraints to the intracellular tyrosine kinase domains. By using MD simulations we have attempted to discriminate between correct and improper helix-helix packing by examining the structural and energetic properties of preformed left-handed and right-handed structures in a fully hydrated DMPC bilayer. The best energetic balance between the residues at the helix-helix interface and the residues exposed to the lipids is obtained for helices in symmetrical left-handed interactions packed together via Glu side chain/Ala backbone interhelical hydrogen bonds. Analyses demonstrate the importance of the ATVEG motif in helix-helix packing and point to additional contacting residues necessary for association. Our findings, all consistent with experimental data, suggest that a symmetrical left-handed structure of the helices could be the transmembrane domain configuration that promotes receptor activation and transformation. The present study may provide further insight into signal transduction mechanisms of the ErbB/Neu receptors.     

Molecular dynamics simulations of the transmembrane domain of the oncogenic ErbB2 receptor dimer in a DMPC bilayer

Molecular dynamics simulations of an atomic model of the transmembrane domain of the oncogenic ErbB2 receptor dimer embedded in an explicit dimyristoylphosphatidylcholine (DMPC) bilayer were performed for more than 4 ns. The oncogenic Glu mutation in the membrane spanning segment plays a major role in tyrosine kinase activity and receptor dimerization, and is thought to be partly responsible for the structure of the transmembrane domain of the active receptor. MD results show that the interactions between the two transmembrane helices are characteristic of a left-handed packing as previously demonstrated from in vacuo simulations. Moreover, MD results reveal the absence of persistent hydrogen bonds between the Glu side chains in a membrane environment, which raise the question of the ability for Glu alone to stabilize the TM domain of the ErbB2 receptor. Interestingly the formation of the alpha-pi motif in the two ErbB2 transmembrane helices confirms the concept of intrinsic sequence-induced conformational flexibility. From a careful analysis of our MD results, we suggest that the left-handed helix-helix packing could be the key to correctly orient the intracellular domain of the activated receptor dimer. The prediction of such interactions from computer simulations represents a new step towards the understanding of signaling mechanisms.     

Sequence context strongly modulates association of polar residues in transmembrane helices

Polar residues are capable of mediating the association of membrane-embedded helices through the formation of side-chain/side-chain inter-helical hydrogen bonds. However, the extent to which native van der Waals packing of the residues surrounding the polar locus can enhance, or interfere with, the interaction of polar residues has not yet been studied. We examined the propensities of four polar residues (aspartic acid, asparagine, glutamic acid, and glutamine) to promote self-association of transmembrane (TM) domains in several biologically derived sequence environments, including (i). four naturally occurring TM domains that contain a Glu or Gln residue (Tnf5/CD40 ligand, C79a/Ig-alpha, C79b/Ig-beta, and Fut3/alpha-fucosyltransferase); and (ii). variants of bacteriophage M13 major coat protein TM segment with Asp and Asn at interfacial and non-interfacial positions. Self-association was quantified by the TOXCAT assay, which measures TM helix self-oligomerization in the Escherichia coli inner membrane. While an appropriately placed polar residue was found in several cases to significantly stabilize TM helix-helix interactions through the formation of an interhelical hydrogen bond, in other cases the strongly polar residues did not enhance the association of the two helices. Overall, these results suggest that an innate structural mechanism may operate to control non-specific association of membrane-embedded polar residues.     

Dimer interface of transmembrane domains for neu/erbB-2 receptor dimerization and transforming activation: a model revealed by molecular dynamics simulations

The specific point mutation Val-->Glu664 within the transmembrane domain of the neu/erbB-2 receptor is associated with increased receptor dimerization and increased receptor tyrosine kinase activity resulting in malignant transformation of cells. It is well established that Glu and residues in proximity are necessary for receptor dimerization but many studies suggest that other intramembrane constraints, not yet elucidated, are determinant for transformation. In this work, we investigated dimer models both to understand the structural role of the Glu mutation in the transmembrane domain association and to determine helix-helix contacts required for oncogenic transformation. Different types of helix-helix association based on data resulting from Cys mutational studies of the full wild receptor and spectroscopic data of transmembrane neu peptides have been explored by molecular dynamics simulations. The study leads to propose a model for the dimeric association of the transmembrane domains of the oncogenic neu receptor showing left-handed interactions of the two helices stabilized by symmetrical hydrogen bonding interactions involving the Glu side chain on one helix and the facing carbonyl of Ala661 on the second helix. Contacting residues observed in the symmetric interface explain the transforming activity or the non transforming activity of many neu mutants. Moreover the left-handed coiled coil structure is fully consistent with recent results proving the role of rotational linkage of the transmembrane domain with the kinase domain. Comparison between the predicted dimer model and those presumed from experiments strongly suggests helix flexibility in the extracellular juxtamembrane region.     

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.     

Comparison of helix interactions in membrane and soluble alpha-bundle proteins

Helix-helix interactions are important for the folding, stability, and function of membrane proteins. Here, two independent and complementary methods are used to investigate the nature and distribution of amino acids that mediate helix-helix interactions in membrane and soluble alpha-bundle proteins. The first method characterizes the packing density of individual amino acids in helical proteins based on the van der Waals surface area occluded by surrounding atoms. We have recently used this method to show that transmembrane helices pack more tightly, on average, than helices in soluble proteins. These studies are extended here to characterize the packing of interfacial and noninterfacial amino acids and the packing of amino acids in the interfaces of helices that have either right- or left-handed crossing angles, and either parallel or antiparallel orientations. We show that the most abundant tightly packed interfacial residues in membrane proteins are Gly, Ala, and Ser, and that helices with left-handed crossing angles are more tightly packed on average than helices with right-handed crossing angles. The second method used to characterize helix-helix interactions involves the use of helix contact plots. We find that helices in membrane proteins exhibit a broader distribution of interhelical contacts than helices in soluble proteins. Both helical membrane and soluble proteins make use of a general motif for helix interactions that relies mainly on four residues (Leu, Ala, Ile, Val) to mediate helix interactions in a fashion characteristic of left-handed helical coiled coils. However, a second motif for mediating helix interactions is revealed by the high occurrence and high average packing values of small and polar residues (Ala, Gly, Ser, Thr) in the helix interfaces of membrane proteins. Finally, we show that there is a strong linear correlation between the occurrence of residues in helix-helix interfaces and their packing values, and discuss these results with respect to membrane protein structure prediction and membrane protein stability.     

Transmembrane peptides from tyrosine kinase receptor. Mutation-related behavior in a lipid bilayer investigated by molecular dynamics simulations

Polar mutations in transmembrane alpha helices may alter the structural details of the hydrophobic sequences and control intermolecular contacts. We have performed molecular dynamics simulations on the transmembrane domain of the proto-oncogenic and the oncogenic forms of the Neu receptor in a fluid DMPC bilayer to test whether the Glu mutation which replaces the Val residue at position 664 may alter the helical structure and its insertion in the membrane. The simulations show that the wild and the mutant forms of the transmembrane domain have a different behavior in the bilayer. The native transmembrane sequence is found to be more flexible than in the presence of the Glu mutation, characterized by a tendency to pi deformation to accommodate the helix length to the membrane thickness. The mutant form of this domain does not evidence helical deformation in the present simulation. Hydrophobic matching is achieved both by a larger helix tilt and a vertical shift of the helix towards the membrane interface, favoring the accessibility of the Glu side chain to the membrane environment. A rapid exchange of hydrogen bond interactions with the surrounding water molecules and the lipid headgroups is observed. The difference in the behavior between the two peptides in a membrane environment was also observed experimentally. Both simulation and experimental results agree with the hypothesis that water may act as an intermediate for the formation of cross links between the facing Glu side chains stabilizing the dimer.     

Structure prediction of the dimeric neu/ErbB-2 transmembrane domain from multi-nanosecond molecular dynamics simulations

Dimerization of the neu/ErbB-2 receptor tyrosine kinase is a necessary but not a sufficient step for signaling. Despite the efforts expended to identify the molecular interactions responsible for receptor-receptor contacts and particularly those involving the transmembrane domain, structural details are still unknown. In this work, molecular dynamics simulations of the helical transmembrane domain (TM) of neu and ErbB-2 receptors are used to predict their dimer structure both in the wild and oncogenic forms. A global conformational search method, applied to define the best orientations of parallel helices, showed an energetically favorable configuration with the specific mutation site within the interface, common for both the nontransforming and the transforming neu/ErbB-2 TM dimers. Starting from this configuration, a total of 10 simulations, about 1.4 ns each, performed in vacuum, without any constraints, show that the two helices preferentially wrap in left-handed interactions with a packing angle at about 20°. The resulting structures are nonsymmetric and the hydrogen bond network analysis shows that helices experience π local distortions that facilitate inter-helix hydrogen bond interactions and may result in a change in the helix packing, leading to a symmetric interface. For the mutated sequences, we show that the Glu side chain interacts directly with its cognate or with carbonyl groups of the facing backbone. We show that the connectivity between interfacial residues conforms to the knobs-into-holes packing mode of transmembrane helices. The dimeric interface described in our models is discussed with respect to mutagenesis studies. Received: 12 March 1999 / Revised version: 23 August 1999 / Accepted: 23 August 1999     

Structure and dynamics of the second and third transmembrane domains of human glycine receptor

A 61-residue polypeptide resembling the second and third transmembrane domains (TM23) of the alpha-1 subunit of human glycine receptor and its truncated form, both with the wild-type loop linking the two TM domains (the "23" loop), were studied using high-resolution NMR. Well-defined domain structures can be identified for the TM2, 23 loop, and TM3 regions. Contrary to the popular model of a long and straight alpha-helical structure for the pore-lining TM2 domain for the Cys-loop receptor family, the last three residues of the TM2 domain and the first eight residues of the 23 loop (S16-S26) seem to be intrinsically nonhelical and highly flexible even in trifluoroethanol, a solvent known to promote and stabilize alpha-helical structures. The six remaining residues of the 23 loop and most of the TM3 domain exhibit helical structures with a kinked pi-helix (or a pi-turn) from W34 to C38 and a kink angle of 159 +/- 3 degrees . The tertiary fold of TM3 relative to TM2 is defined by several unambiguously identified long-range NOE cross-peaks within the loop region and between TM2 and TM3 domains. The 20 lowest-energy structures show a left-handed tilt of TM3 relative to TM2 with a tilting angle of 44 +/- 2 degrees between TM2 (V1-Q14) and TM3 (L39-E48) helix axes. This left-handed TM2-TM3 arrangement ensures a neatly packed right-handed quaternary structure of five subunits to form an ion-conducting pore. This is the first time that two TM domains of the glycine receptor linked by the important 23 loop have ever been analyzed at atomistic resolution. Many structural characteristics of the receptor can be inferred from the structural and dynamical features identified in this study.     

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.     

Molecular dynamics simulations of the E1/E2 transmembrane domain of the Semliki Forest virus

Transmembrane (TM) helix-helix interactions are important for virus budding and fusion. We have developed a simulation strategy that reveals the main features of the helical packing between the TM domains of the two glycoproteins E1 and E2 of the alpha-virus Semliki Forest virus and that can be extrapolated to sketch TM helical packing in other alpha-viruses. Molecular dynamics simulations were performed in wild-type and mutant peptides, both isolated and forming E1/E2 complexes. The simulations revealed that the isolated wild-type E1 peptide formed a more flexible helix than the rest of peptides and that the wild-type E1/E2 complex consists of two helices that intimately pack their N-terminals. The residues located at the interhelical interface displayed the typical motif of the left-handed coiled-coils. These were small and medium residues as Gly, Ala, Ser, and Leu, which also had the possibility to form interhelical Calpha-H...O hydrogen bonds. Results from the mutant complexes suggested that correct packing is a compromise between these residues at both E1 and E2 interhelical interfaces. This compromise allowed prediction of E1-E2 contact residues in the TM spanning domain of other alphaviruses even though the sequence identity of E2 peptides is low in this domain.     

Transmembrane interactions in the activation of the Neu receptor tyrosine kinase

The Neu receptor tyrosine kinase is constitutively activated by a single amino acid change in the transmembrane domain of the receptor. The mutation of Val664 to glutamate or glutamine induces receptor dimerization and autophosphorylation of the receptor's intracellular kinase domain. The ability of this single mutation to activate the receptor is sequence-dependent, suggesting that specific helix-helix interactions stabilize the transmembrane dimer. We have determined the local secondary structure and interhelical contacts in the region of position 664 in peptide models of the activated receptor using solid-state rotational resonance and rotational echo double-resonance (REDOR) NMR methods. Intrahelical (13)C rotational resonance distance measurements were made between 1-(13)C-Thr662 and 2-(13)C-Gly665 on peptides corresponding to the wild-type Neu and activated Neu transmembrane sequences containing valine and glutamate at position 664, respectively. We observed similar internuclear distances (4.5 +/- 0.2 A) in both Neu and Neu*, indicating that the region near residue 664 is helical and is not influenced by mutation. Interhelical (15)N...(13)C REDOR measurements between Gln664 side chains on opposing helices were not consistent with hydrogen bonding between the side chain functional groups. However, interhelical rotational resonance measurements between 1-(13)C-Glu664 and 2-(13)C-Gly665 and between 1-(13)C-Gly665 and 2-(13)C-Gly665 demonstrated close contacts (4.3-4.5 A) consistent with the packing of Gly665 in the Neu* dimer interface. These measurements provide structural constraints for modeling the transmembrane dimer and define the rotational orientation of the transmembrane helices in the activated receptor.     

Transmembrane Recognition of the Semaphorin Co-Receptors Neuropilin 1 and Plexin A1: Coarse-Grained Simulations

The cancer associated class 3 semaphorins require direct binding to neuropilins and association to plexins to trigger cell signaling. Here, we address the role of the transmembrane domains of neuropilin 1 and plexin A1 for the dimerization of the two receptors by characterizing the assembly in lipid bilayers using coarse-grained molecular dynamics simulations. From experimental evidence using a two-hybrid system showing the biochemical association of the two receptors transmembrane domains, we performed molecular simulations in DOPC and POPC demonstrating spontaneously assembly to form homodimers and heterodimers with a very high propensity for right-handed packing of the helices. Inversely, left-handed packing was observed with a very low propensity. This mode of packing was observed uniquely when the plexin A1 transmembrane domain was involved in association. Potential of mean force calculations were used to predict a hierarchy of self-association for the monomers: the two neuropilin 1 transmembrane domains strongly associated, neuropilin 1 and plexin A1 transmembrane domains associated less and the two plexin A1 transmembrane domains weakly but significantly associated. We demonstrated that homodimerization and heterodimerization are driven by GxxxG motifs, and that the sequence context modulates the packing mode of the plexin A1 transmembrane domains. This work presents major advances towards our understanding of membrane signaling platforms assembly through membrane domains and provides exquisite information for the design of antagonist drugs defining a novel class of therapeutic agents.     

Structure, topology, and tilt of cell-signaling peptides containing nuclear localization sequences in membrane bilayers determined by solid-state NMR and molecular dynamics simulation studies

Cell-signaling peptides have been extensively used to transport functional molecules across the plasma membrane into living cells. These peptides consist of a hydrophobic sequence and a cationic nuclear localization sequence (NLS). It has been assumed that the hydrophobic region penetrates the hydrophobic lipid bilayer and delivers the NLS inside the cell. To better understand the transport mechanism of these peptides, in this study, we investigated the structure, orientation, tilt of the peptide relative to the bilayer normal, and the membrane interaction of two cell-signaling peptides, SA and SKP. Results from CD and solid-state NMR experiments combined with molecular dynamics simulations suggest that the hydrophobic region is helical and has a transmembrane orientation with the helical axis tilted away from the bilayer normal. The influence of the hydrophobic mismatch, between the hydrophobic length of the peptide and the hydrophobic thickness of the bilayer, on the tilt angle of the peptides was investigated using thicker POPC and thinner DMPC bilayers. NMR experiments showed that the hydrophobic domain of each peptide has a tilt angle of 15 +/- 3 degrees in POPC, whereas in DMPC, 25 +/- 3 degree and 30 +/- 3 degree tilts were observed for SA and SKP peptides, respectively. These results are in good agreement with molecular dynamics simulations, which predict a tilt angle of 13.3 degrees (SA in POPC), 16.4 degrees (SKP in POPC), 22.3 degrees (SA in DMPC), and 31.7 degrees (SKP in DMPC). These results and simulations on the hydrophobic fragment of SA or SKP suggest that the tilt of helices increases with a decrease in bilayer thickness without changing the phase, order, and structure of the lipid bilayers.     

Hydrophobic matching controls the tilt and stability of the dimeric platelet-derived growth factor receptor (PDGFR) β transmembrane segment

The platelet-derived growth factor receptor β is a member of the cell surface receptor tyrosine kinase family and dimerizes upon activation. We determined the structure of the transmembrane segment in dodecylphosphocholine micelles by liquid-state NMR and found that it forms a stable left-handed helical dimer. Solid-state NMR and oriented circular dichroism were used to measure the tilt angle of the helical segments in macroscopically aligned model membranes with different acyl chain lengths. Both methods showed that decreasing bilayer thickness (DEPC-POPC-DMPC) led to an increase in the helix tilt angle from 10° to 30° with respect to the bilayer normal. At the same time, reconstitution of the comparatively long hydrophobic segment became less effective, eventually resulting in complete protein aggregation in the short-chain lipid DLPC. Unrestrained molecular dynamics simulations of the dimer were carried out in explicit lipid bilayers (DEPC, POPC, DMPC, sphingomyelin), confirming the observed dependence of the helix tilt angle on bilayer thickness. Notably, molecular dynamics revealed that the left-handed dimer gets tilted en bloc, whereas conformational transitions to alternative (e.g. right-handed dimeric) states were not supported. The experimental data along with the simulation results demonstrate a pronounced interplay between the platelet-directed growth factor receptor β transmembrane segment and the bilayer thickness. The effect of hydrophobic mismatch might play a key role in the redistribution and activation of the receptor within different lipid microdomains of the plasma membrane in vivo.     

Identifying interactions between transmembrane helices from the adenosine A2A receptor

The human adenosine A(2A) receptor (A(2A)R) is an integral membrane protein and a member of the G-protein-coupled receptor (GPCR) superfamily, characterized by seven transmembrane (TM) helices. Although helix-helix association in the lipid bilayer is known to be an essential step in the folding of GPCRs, the determinants of their structures, folding, and assembly in the cell membrane are poorly understood. Previous studies in our group showed that while peptides corresponding to all seven TM domains of A(2A)R form stable helical structures in detergent micelles and lipid vesicles, they display significant variability in their helical propensity. This finding suggested to us that some TM domains might need to interact with other domains to properly insert and fold in hydrophobic environments. In this study, we assessed the ability of TM peptides to interact in pairwise combinations. We analyzed peptide interactions in hydrophobic milieus using circular dichroism spectroscopy and F?rster resonance energy transfer. We find that specific interactions between TM helices occur, leading to additional helical content, especially in weakly helical TM domains, suggesting that some TM domains need a partner for proper folding in the membrane. The approach developed in this study will enable complete analysis of the TM domain interactions and the modeling of a folding pathway for A(2A)R.     

A helix heterodimer in a lipid bilayer: prediction of the structure of an integrin transmembrane domain via multiscale simulations

Dimerization of transmembrane (TM) α helices of membrane receptors plays a key role in signaling. We show that molecular dynamics simulations yield models of integrin TM helix heterodimers, which agree well with available NMR structures. We use?a multiscale simulation approach, combining coarse-grained and subsequent atomistic simulation, to model the dimerization of wild-type (WT) and mutated sequences of the αIIb and β3 integrin TM helices. The WT helices formed a stable, right-handed dimer with the same helix-helix interface as in the published NMR structure (PDB: 2K9J). In contrast, the presence of disruptive mutations perturbed the interface between the helices, altering the conformational stability of the dimer. The αIIb/β3 interface was more flexible than that of, e.g., glycophorin A. This is suggestive of a role for alternative packing modes of the TM helices in transbilayer signaling.     

Parallel helix bundles and ion channels: molecular modeling via simulated annealing and restrained molecular dynamics.

A parallel bundle of transmembrane (TM) alpha-helices surrounding a central pore is present in several classes of ion channel, including the nicotinic acetylcholine receptor (nAChR). We have modeled bundles of hydrophobic and of amphipathic helices using simulated annealing via restrained molecular dynamics. Bundles of Ala20 helices, with N = 4, 5, or 6 helices/bundle were generated. For all three N values the helices formed left-handed coiled coils, with pitches ranging from 160 A (N = 4) to 240 A (N = 6). Pore radius profiles revealed constrictions at residues 3, 6, 10, 13, and 17. A left-handed coiled coil and a similar pattern of pore constrictions were observed for N = 5 bundles of Leu20. In contrast, N = 5 bundles of Ile20 formed right-handed coiled coils, reflecting loosened packing of helices containing beta-branched side chains. Bundles formed by each of two classes of amphipathic helices were examined: (a) M2a, M2b, and M2c derived from sequences of M2 helices of nAChR; and (b) (LSSLLSL)3, a synthetic channel-forming peptide. Both classes of amphipathic helix formed left-handed coiled coils. For (LSSLLSL)3 the pitch of the coil increased as N increased from 4 to 6. The M2c N = 5 helix bundle is discussed in the context of possible models of the pore domain of nAChR.     

Molecular modeling of c-erbB2 receptor dimerization: Coiled-coil structure of wild and oncogenic transmembrane domains—Stabilization by interhelical hydrogen bonds in the oncogenic form

Dimerization models of c-erbB2 transmembrane domains (Leu651-Ile675) are studied by molecular mechanics and molecular dynamics simulations. Both wild and Glu mutated transmembrane helices exhibit the same relative orientation for favorable associations and dimerize preferentially in left-handed coiled-coil structures. The mutation point 659 belongs to the interfacing residues, and in the transforming domain, symmetric hydrogen bonds between Glu carboxylic groups stabilize the dimeric structure. The same helix packing found for the wild dimers, except side-chain—side-chain hydrogen bonds, suggests that the transmembrane domains dimerize according to similar process. Structural and energetical characterization of the models are presented. © 1997 John Wiley & Sons, Inc. Biopoly 42: 157–168, 1997