Spectroscopy-Based Modelling of the 3D Structure of the β Subunit of the High Affinity IgE Receptor

Abstract

The high affinity IgE receptor, possesses a tetrameric structure. The 243 residue β subunit is a polytopic protein with four hydrophobic membrane-spanning segments, whereas the individual α and γ subunits are bitopic proteins each containing one transmembrane domain in their monomeric form. In the proposed topographical model (Blank et al., 1989), the four trans-membrane α helices of the β subunit are connected by three loop sequences.

To study the individual subunits and intact receptor, this membrane protein was divided into domains such as its loop peptides, cytoplasmic peptides and transmembrane helices according to Blank et al., 1989. The 3D structure of the synthesized loop peptides and cytoplasmic peptides were calculated; CD and/or NMR data were used as appropriate to generate the resultant structures which were then used as data basis for the higher level calculations.

The four individual transmembrane helices of the β subunit were characterised, first of all, by mapping the relative lipophilicity of their surfaces using lipophilic probes. A second procedure, docking of the individual helices in pairs, was used to predict helix–helix interactions.

The data on the relative lipophilicity of the surfaces as well as the surfaces that favoured helix–helix interactions were used in combination with the spectroscopy-based structures of the loops and cytoplasmic domains to calculate via molecular dynamics, the helix arrangement and 3D structure of the β subunit of the high affinity IgE receptor. In the final analysis, the molecular simulations yielded two structures of the β subunit, which should form a basis for the modelling of the whole high affinity IgE receptor.     

Computational differentiation of N-terminal signal peptides and transmembrane helices

Signal peptides and transmembrane helices both contain a stretch of hydrophobic amino acids. This common feature makes it difficult for signal peptide and transmembrane helix predictors to correctly assign identity to stretches of hydrophobic residues near the N-terminal methionine of a protein sequence. The inability to reliably distinguish between N-terminal transmembrane helix and signal peptide is an error with serious consequences for the prediction of protein secretory status or transmembrane topology. In this study, we report a new method for differentiating protein N-terminal signal peptides and transmembrane helices. Based on the sequence features extracted from hydrophobic regions (amino acid frequency, hydrophobicity, and the start position), we set up discriminant functions and examined them on non-redundant datasets with jackknife tests. This method can incorporate other signal peptide prediction methods and achieve higher prediction accuracy. For Gram-negative bacterial proteins, 95.7% of N-terminal signal peptides and transmembrane helices can be correctly predicted (coefficient 0.90). Given a sensitivity of 90%, transmembrane helices can be identified from signal peptides with a precision of 99% (coefficient 0.92). For eukaryotic proteins, 94.2% of N-terminal signal peptides and transmembrane helices can be correctly predicted with coefficient 0.83. Given a sensitivity of 90%, transmembrane helices can be identified from signal peptides with a precision of 87% (coefficient 0.85). The method can be used to complement current transmembrane protein prediction and signal peptide prediction methods to improve their prediction accuracies.     

Interactions of alpha-helices with lipid bilayers: a review of simulation studies

Membrane proteins, of which the majority seem to contain one or more alpha-helix, constitute approx. 30% of most genomes. A complete understanding of the nature of helix/bilayer interactions is necessary for an understanding of the structural principles underlying membrane proteins. This review describes computer simulation studies of helix/bilayer interactions. Key experimental studies of the interactions of alpha-helices and lipid bilayers are briefly reviewed. Surface associated helices are found in some membrane-bound enzymes (e.g. prostaglandin synthase), and as stages in the mechanisms of antimicrobial peptides and of pore-forming bacterial toxins. Transmembrane alpha-helices are found in most integral membrane proteins, and also in channels formed by amphipathic peptides or by bacterial toxins. Mean field simulations, in which the lipid bilayer is approximated as a hydrophobic continuum, have been used in studies of membrane-active peptides (e.g. alamethicin, melittin, magainin and dermaseptin) and of simple membrane proteins (e.g. phage Pf1 coat protein). All atom molecular dynamics simulations of fully solvated bilayers with transmembrane helices have been applied to: the constituent helices of bacteriorhodopsin; peptide-16 (a simple model TM helix); and a number of pore-lining helices from ion channels. Surface associated helices (e.g. melittin and dermaseptin) have been simulated, as have alpha-helical bundles such as bacteriorhodopsin and alamethicin. From comparison of the results from the two classes of simulation, it emerges that a major theoretical challenge is to exploit the results of all atom simulations in order to improve the mean field approach.     

Efflux by Small Multidrug Resistance Proteins Is Inhibited by Membrane-interactive Helix-stapled Peptides

Bacterial cell membranes contain several protein pumps that resist the toxic effects of drugs by efficiently extruding them. One family of these pumps, the small multidrug resistance proteins (SMRs), consists of proteins of about 110 residues that need to oligomerize to form a structural pathway for substrate extrusion. As such, SMR oligomerization sites should constitute viable targets for efflux inhibition, by disrupting protein-protein interactions between helical segments. To explore this proposition, we are using Hsmr, an SMR from Halobacter salinarum that dimerizes to extrude toxicants. Our previous work established that (i) Hsmr dimerization is mediated by a helix-helix interface in Hsmr transmembrane (TM) helix 4 (residues ⁹⁰GLALIVAGV⁹⁸); and (ii) a peptide comprised of the full TM4(85–105) sequence inhibits Hsmr-mediated ethidium bromide efflux from bacterial cells. Here we define the minimal linear sequence for inhibitor activity (determined as TM4(88–100), and then “staple” this sequence via Grubbs metathesis to produce peptides typified by acetyl-A-(Sar)₃-⁸⁸VVGLXLIZXGVVV¹⁰⁰-KKK-NH₂ (X = 2-(4′-pentenyl)alanine at positions 92 and 96; Z = Val, Gly, or Asn at position 95)). The Asn⁹⁵ peptide displayed specific efflux inhibition and resensitization of Hsmr-expressing cells to ethidium bromide; and was non-hemolytic to human red blood cells. Stapling essentially prevented peptide degradation in blood plasma and liver homogenates versus an unstapled counterpart. The overall results confirm that the stapled analog of TM4(88–100) retains the structural complementarity required to disrupt the Hsmr TM4-TM4 locus in Hsmr, and portend the general validity of stapled peptides as therapeutics for the disruption of functional protein-protein interactions in membranes.     

Characterization of an archaeal multidrug transporter with a unique amino acid composition

The Smr family of multidrug transporters consists of small membrane proteins that extrude various drugs in exchange with protons rendering cells resistant to these drugs. Smr proteins identified to date have been found only in Eubacteria. In this work we present the cloning and characterization of an Smr protein from the archaeon Halobacterium salinarum, the first Smr in the archaeal kingdom. The protein, named Hsmr, was identified through sequence similarity to the Smr family, and the DNA sequence was cloned into an Escherichia coli expression system. Hsmr is heterologously expressed in a functional form despite the difference in lipid composition of the membrane and the lower salt in the cell and its environment. Cells harboring the Hsmr plasmid transport ethidium bromide in an uncoupler-sensitive process and gain resistance to ethidium bromide and acriflavine. Hsmr binds tetraphenylphosphonium (TPP(+)) with a relatively low affinity (K(D) approximately 200 nm) at low salt concentration that increases (K(D) approximately 40 nm) upon the addition of 2 m of either NaCl or KCl. The Hsmr protein contains many of the signature sequence elements of the Smr family and also a high content of negative residues in the loops, characteristic of extreme halophiles. Strikingly, Hsmr is composed of over 40% valine and alanine residues. These residues are clustered at certain regions of the protein in domains that are not important for activity, as judged from lack of conservation and from previous studies with other Smr proteins. We suggest that this high content of alanine and valine residues is a reflection of a "natural" alanine and valine scanning necessitated by the high GC content of the gene. This phenomenon reveals significant sequence elements in small multidrug transporters.     

Helix packing in polytopic membrane proteins: role of glycine in transmembrane helix association.

The nature and distribution of amino acids in the helix interfaces of four polytopic membrane proteins (cytochrome c oxidase, bacteriorhodopsin, the photosynthetic reaction center of Rhodobacter sphaeroides, and the potassium channel of Streptomyces lividans) are studied to address the role of glycine in transmembrane helix packing. In contrast to soluble proteins where glycine is a noted helix breaker, the backbone dihedral angles of glycine in transmembrane helices largely fall in the standard alpha-helical region of a Ramachandran plot. An analysis of helix packing reveals that glycine residues in the transmembrane region of these proteins are predominantly oriented toward helix-helix interfaces and have a high occurrence at helix crossing points. Moreover, packing voids are generally not formed at the position of glycine in folded protein structures. This suggests that transmembrane glycine residues mediate helix-helix interactions in polytopic membrane proteins in a fashion similar to that seen in oligomers of membrane proteins with single membrane-spanning helices. The picture that emerges is one where glycine residues serve as molecular notches for orienting multiple helices in a folded protein complex.     

Alpha-periodicity analysis of small multidrug resistance (SMR) efflux transporters.

Proteins in the small multidrug resistance (SMR) family of transport proteins are about 110 amino acids in length and are predicted to have four transmembrane helices. This family is divided into a two groups, one of which we have referred to as small multidrug pumps (Smp) and confer resistance to a wide variety of quaternary ammonium compounds through a proton-drug efflux antiport mechanism. Members of the second group within this family have, as yet, not had their substrate profile characterized and are referred to as Sug proteins. Alpha-periodicity analysis was conducted on a set of six homologous proteins of the SMR family consisting of three established Smp and three Sug proteins. Several amino acid properties were used in the analysis including hydropathy, variability, and a substitution matrix for lipid exposed amino acids. The scanning window was varied between 8 and 14 residues and the alpha-periodicity was calculated from the peaks in the Fourier transform power spectra in the region between 3.0 and 4.3 residues/turn. This analysis adds to the hydropathy analysis to give a more confident prediction of which residues are within the lipid bilayer for each of the four transmembrane helices. Information was also obtained that allowed for the identification of zones within each transmembrane helix that face the interior of the helical bundle on one side and are lipid exposed on the other face.     

Hydrophobic organization of α-helix membrane bundle in bacteriorhodopsin

The hydrophobic organization of the intramembraneα-helical bundle in bacteriorhodopsin (BRh) was assessed based on a new approach to characterization of spatial hydrophobic properties of transmembrane (TM)α-helical peptides. The method employs two independent techniques: Monte Carlo simulations of nonpolar solvent around TM peptides and analysis of molecular hydrophobicity potential on their surfaces. The results obtained by the two methods agree with each other and permit precise hydrophobicity mapping of TM peptides. Superimposition of such data on the experimentally derived spatial model of the membrane moiety together with 2D maps of hydrophobic hydrophilic contacts provide considerable insight into the hydrophobic organization of BRh. The helix bundle is stabilized to a large extent by hydrophobic interactions between helices—neighbors in the sequence of BRh, by long-range interactions in helix pairs C-E, C-F, and C-G, and by nonpolar contracts between retinal and helices C, D, E, F. Unlike globular proteins, no polar contacts between residues distantly separated in the sequence of BRh were found in the bundle. One of the most striking results of this study is the finding that the hydrophobic organization of BRh is significantly different from those in bacterial photoreaction centers. Thus, TMα-helices in BRh expose their most nonpolar sides to the bilayer as well as to the neighboring helices and to the interior of the bundle. Some of them contact lipids with their relatively hydrophilic surfaces. No correlation was found between disposition of the most hydrophobic and the most variable sides of the TM helices.     

Structural analysis of synthetic peptide fragments from EmrE,a multidrug resistance protein,in a membrane-mimetic environment

EmrE, a multidrug resistance protein from Escherichia coli, renders the bacterium resistant to a variety of cytotoxic drugs by active translocation out of the cell. The 110-residue sequence of EmrE limits the number of structural possibilities that can be envisioned for this membrane protein. Four helix bundle models have been considered [Yerushalmi, H., Lebendiker, M., and Schuldiner, S. (1996) J. Biol. Chem. 271, 31044-31048]. The validity of EmrE structural models has been probed experimentally by investigations on overlapping peptides (ranging in length from 19 to 27 residues), derived from the sequence of EmrE. The choice of peptides was made to provide sequences of two complete, predicted transmembrane helices (peptides H1 and H3) and two helix-loop-helix motifs (peptides A and B). Peptide (B) also corresponds to a putative hairpin in a speculative beta-barrel model, with the "Pro-Thr-Gly" segment forming a turn. Structure determination in SDS micelles using NMR indicates peptide H1 to be predominantly helical, with helix boundaries in the micellar environment corroborating predicted helical limits. Peptide A adopts a helix-loop-helix structure in SDS micelles, and peptide B was also largely helical in micellar environments. An analogue peptide, C, in which the central "Pro-Thr-Gly" was replaced by "(D)Pro-Gly" displays local turn conformation at the (D)Pro-Gly segment, but neither a continuous helical stretch nor beta-hairpin formation was observed. This study implies that the constraints of membrane and micellar environments largely direct the structure of transmembrane peptides and proteins and study of judiciously selected peptide fragments can prove useful in the structural elucidation of membrane proteins.     

Anti-apoptotic Bcl-XL protein in complex with BH3 peptides of pro-apoptotic Bak, Bad, and Bim proteins: comparative molecular dynamics simulations

The Bcl-2 family of proteins plays a central role in the regulation of mitochondrial outer-membrane permeabilization, a critical step in apoptosis. Heterodimerization between the pro- and anti-apoptotic members of Bcl-2 family is a key event in this process. Anti-apoptotic proteins have high levels of expression in many cancers and they have different affinities for different pro-apoptotic proteins. Experimentally determined structures of all members of Bcl-2 proteins have remarkably similar helical fold despite poor amino acid sequence identity. Peptides representing BH3 region of pro-apoptotic proteins have been shown to bind the hydrophobic cleft of anti-apoptotic proteins and this segment is responsible in modulating the apoptotic pathways in living cells. Understanding the molecular basis of protein-protein recognition is required to develop inhibitors specific to a particular anti-apoptotic protein. We have carried out molecular dynamics simulations on the anti-apoptotic Bcl-X(L) protein in complex with three different BH3 peptides derived from pro-apoptotic Bak, Bad and Bim proteins. Each complex structure was simulated for a period of 50 ns after 2.5 ns equilibration. Analysis of the simulation results showed that in the Bcl-X(L) protein, the helix containing the BH3 region is more flexible than other helices in all three simulations. A network of strong hydrophobic interactions exists between four of the six helices and they contribute significantly to the stability of this helix bundle protein. Analysis of Bcl-X(L)-BH3 peptide interactions reveals the role of loop residues in the protein-peptide interactions in all three simulations. Bad and Bim peptides maintain strong hydrophobic and hydrophilic interactions with the helix preceding the central hydrophobic helix. Residues from this helix interact with an Arg residue in Bad and Bim peptides. This Arg residue is next to the conserved Leu residue and is replaced by Ala in Bak. Absence of these interactions and the helix propensity are likely to be the cause for Bak peptide's weaker binding affinity with the Bcl-X(L) protein. The results of this study have implications in the design of Bcl-X(L)-specific inhibitors.     

Modulation of substrate efflux in bacterial small multidrug resistance proteins by mutations at the dimer interface

Bacteria evade the effects of cytotoxic compounds through the efflux activity of membrane-bound transporters such as the small multidrug resistance (SMR) proteins. Consisting typically of ca. 110 residues with four transmembrane (TM) α-helices, crystallographic studies have shown that TM helix 1 (TM1) through TM helix 3 (TM3) of each monomer create a substrate binding "pocket" within the membrane bilayer, while a TM4-TM4 interaction accounts for the primary dimer formation. Previous work from our lab has characterized a highly conserved small-residue heptad motif in the Halobacterium salinarum transporter Hsmr as (90)GLXLIXXGV(98) that lies along the TM4-TM4 dimer interface of SMR proteins as required for function. Focusing on conserved positions 91, 93, 94, and 98, we substituted the naturally occurring Hsmr residue for Ala, Phe, Ile, Leu, Met, and Val at each position in the Hsmr TM4-TM4 interface. Large-residue replacements were studied for their ability to dimerize on SDS-polyacrylamide gels, to bind the cytotoxic compound ethidium bromide, and to confer resistance by efflux. Although the relative activity of mutants did not correlate with dimer strength for all mutants, all functional mutants lay within 10% of dimerization relative to the wild type (WT), suggesting that the optimal dimer strength at TM4 is required for proper efflux. Furthermore, nonfunctional substitutions at the center of the dimerization interface that do not alter dimer strength suggest a dynamic TM4-TM4 "pivot point" that responds to the efflux requirements of different substrates. This functionally critical region represents a potential target for inhibiting the ability of bacteria to evade the effects of cytotoxic compounds.     

Dynamic regulation of lipid–protein interactions

We review the importance of helix motions for the function of several important categories of membrane proteins and for the properties of several model molecular systems. For voltage-gated potassium or sodium channels, sliding, tilting and/or rotational movements of the S4 helix accompanied by a swapping of cognate side-chain ion-pair interactions regulate the channel gating. In the seven-helix G protein-coupled receptors, exemplified by the rhodopsins, collective helix motions serve to activate the functional signaling. Peptides which initially associate with lipid-bilayer membrane surfaces may undergo dynamic transitions from surface-bound to tilted-transmembrane orientations, sometimes accompanied by changes in the molecularity, formation of a pore or, more generally, the activation of biological function. For single-span membrane proteins, such as the tyrosine kinases, an interplay between juxtamembrane and transmembrane domains is likely to be crucial for the regulation of dimer assembly that in turn is associated with the functional responses to external signals. Additionally, we note that experiments with designed single-span transmembrane helices offer fundamental insights into the molecular features that govern protein–lipid interactions.This article is part of a Special Issue entitled: Lipid–protein interactions.     

Structures of the transmembrane helices of the G-protein coupled receptor, rhodopsin.

An hypothesis is tested that individual peptides corresponding to the transmembrane helices of the membrane protein, rhodopsin, would form helices in solution similar to those in the native protein. Peptides containing the sequences of helices 1, 4 and 5 of rhodopsin were synthesized. Two peptides, with overlapping sequences at their termini, were synthesized to cover each of the helices. The peptides from helix 1 and helix 4 were helical throughout most of their length. The N- and C-termini of all the peptides were disordered and proline caused opening of the helical structure in both helix 1 and helix 4. The peptides from helix 5 were helical in the middle segment of each peptide, with larger disordered regions in the N- and C-termini than for helices 1 and 4. These observations show that there is a strong helical propensity in the amino acid sequences corresponding to the transmembrane domain of this G-protein coupled receptor. In the case of the peptides from helix 4, it was possible to superimpose the structures of the overlapping sequences to produce a construct covering the whole of the sequence of helix 4 of rhodopsin. As similar superposition for the peptides from helix 1 also produced a construct, but somewhat less successfully because of the disordering in the region of sequence overlap. This latter problem was more severe for helix 5 and therefore a single peptide was synthesized for the entire sequence of this helix, and its structure determined. It proved to be helical throughout. Comparison of all these structures with the recent crystal structure of rhodopsin revealed that the peptide structures mimicked the structures seen in the whole protein. Thus similar studies of peptides may provide useful information on the secondary structure of other transmembrane proteins built around helical bundles.     

Side chains at the membrane-water interface modulate the signaling state of a transmembrane receptor

Previous model studies of peptides and proteins have shown that protein-lipid interactions, primarily involving amino acid side chains near the membrane-water interface, modulate the position of transmembrane helices in bilayers. The present study examines whether such interfacial side chains stabilize the signaling states of a transmembrane signaling helix in a representative receptor, the aspartate receptor of bacterial chemotaxis. To examine the functional roles of signaling helix side chains at the periplasmic and cytoplasmic membrane-water interfaces, arginine and cysteine substitutions were scanned through these two interfacial regions. The chemical reactivities of the cysteine residues were first measured to determine the positions at which the helix crosses the membrane-water interface in both the periplasmic and cytoplasmic compartments. Subsequently, two antisymmetric in vitro activity measurements were carried out to determine the effect of each interfacial arginine or cysteine substitution on receptor signaling. Substitutions that stabilize the receptor on-state cause upregulation of receptor-coupled kinase activity and inhibition of methylation at receptor adaptation sites, while substitutions that stabilize the off-state have the opposite effects on these two activities. Notably, four substitutions at aromatic tryptophan and phenylalanine positions buried in the membrane near the membrane-water interface were found to stabilize the native on- or off-signaling state. The striking ability of these substitutions to drive the receptor toward a specific signaling state indicates that interfacial side chains are highly optimized to correctly position the native signaling helix in the membrane and to allow normal switching between the on- and off-signaling states. The analogous substitutions in model transmembrane helices are known to drive small piston-type displacements of the helix normal to the membrane. Thus, the simplest molecular interpretation of the present findings is that the signal-stabilizing substitutions drive piston displacements of the signaling helix, providing further support for the piston model for transmembrane signaling in bacterial chemoreceptors. More generally, the findings indicate that the interfacial phenylalanine, tryptophan, and arginine side chains widespread in the transmembrane alpha-helices of receptors, channels, and transporters can play important roles in modulating transitions between signaling and conformational states.     

Fluorescence resonance energy transfer shows a close helix-helix distance in the transmembrane M13 procoat protein

During the membrane insertion process the major coat protein of bacteriophage M13 assumes a conformation in which two transmembrane helices corresponding to the leader sequence and the anchor region in the mature part of the protein coming into close contact with each other. Previous studies on the molecular mechanism of membrane insertion of M13 procoat protein have shown that this interaction between the two helices might drive the actual translocation process. We investigated the intramolecular distance between the two helices of the transmembrane procoat protein by measuring fluorescence resonance energy transfer (FRET) between the donor (Tyr) placed in one helix and the acceptor (Trp) placed in the other helix. Various mutant procoat proteins with differently positioned donor-acceptor pairs were generated, purified, and reconstituted into artificial lipid bilayers. The results obtained from the FRET measurements, combined with molecular modeling, show that the transmembrane helices are in close contact on the order of 1-1.5 nm. The present approach might be of general interest for determining the topology and the folding of membrane proteins.     

Helix packing moments reveal diversity and conservation in membrane protein structure

Helical membrane proteins are more tightly packed and the packing interactions are more diverse than those found in helical soluble proteins. Based on a linear correlation between amino acid packing values and interhelical propensity, we propose the concept of a helix packing moment to predict the orientation of helices in helical membrane proteins and membrane protein complexes. We show that the helix packing moment correlates with the helix interfaces of helix dimers of single pass membrane proteins of known structure. Helix packing moments are also shown to help identify the packing interfaces in membrane proteins with multiple transmembrane helices, where a single helix can have multiple contact surfaces. Analyses are described on class A G protein-coupled receptors (GPCRs) with seven transmembrane helices. We show that the helix packing moments are conserved across the class A family of GPCRs and correspond to key structural contacts in rhodopsin. These contacts are distinct from the highly conserved signature motifs of GPCRs and have not previously been recognized. The specific amino acid types involved in these contacts, however, are not necessarily conserved between subfamilies of GPCRs, indicating that the same protein architecture can be supported by a diverse set of interactions. In GPCRs, as well as membrane channels and transporters, amino acid residues with small side-chains (Gly, Ala, Ser, Cys) allow tight helix packing by mediating strong van der Waals interactions between helices. Closely packed helices, in turn, facilitate interhelical hydrogen bonding of both weakly polar (Ser, Thr, Cys) and strongly polar (Asn, Gln, Glu, Asp, His, Arg, Lys) amino acid residues. We propose the use of the helix packing moment as a complementary tool to the helical hydrophobic moment in the analysis of transmembrane sequences.     

Quasi-symmetry in the cryo-EM structure of EmrE provides the key to modeling its transmembrane domain

Small multidrug resistance (SMR) transporters contribute to bacterial resistance by coupling the efflux of a wide range of toxic aromatic cations, some of which are commonly used as antibiotics and antiseptics, to proton influx. EmrE is a prototypical small multidrug resistance transporter comprising four transmembrane segments (M1-M4) that forms dimers. It was suggested recently that EmrE molecules in the dimer have different topologies, i.e. monomers have opposite orientations with respect to the membrane plane. A 3-D structure of EmrE acquired by electron cryo-microscopy (cryo-EM) at 7.5 Angstroms resolution in the membrane plane showed that parts of the structure are related by quasi-symmetry. We used this symmetry relationship, combined with sequence conservation data, to assign the transmembrane segments in EmrE to the densities seen in the cryo-EM structure. A C alpha model of the transmembrane region was constructed by considering the evolutionary conservation pattern of each helix. The model is validated by much of the biochemical data on EmrE with most of the positions that were identified as affecting substrate translocation being located around the substrate-binding cavity. A suggested mechanism for proton-coupled substrate translocation in small multidrug resistance antiporters provides a mechanistic rationale to the experimentally observed inverted topology.     

Molecular Dynamics Simulations of Homo-oligomeric Bundles Embedded Within a Lipid Bilayer

Using molecular dynamics simulations, we studied the structure, interhelix interactions, and dynamics of transmembrane proteins. Specifically, we investigated homooligomeric helical bundle systems consisting of synthetic α-helices with either the sequence Ac-(LSLLLSL)₃-NH₂ (LS2) or Ac-(LSSLLSL)₃-NH₂ (LS3). The LS2 and LS3 helical peptides are designed to have amphipathic characteristics that form ion channels in membrane. We simulated bundles containing one to six peptides that were embedded in palmitoyl-oleoyl-phosphatidylcholine (POPC) lipid bilayer and placed between two lamellae of water. We aim to provide a fundamental understanding of how amphipathic helical peptides interact with each other and their dynamical behaviors in different homooligomeric states. To understand structural properties, we examined the helix lengths, tilt angles of individual helices and the entire bundle, interhelix distances, interhelix cross-angles, helix hydrophobic-to-hydrophilic vector projections, and the average number of interhelix hydrophilic (serine–serine) contacts lining the pore of the transmembrane channel. To analyze dynamical properties, we calculated the rotational autocorrelation function of each helix and the cross-correlation of the rotational velocity between adjacent helices. The observed structural and dynamical characteristics show that higher order bundles containing four to six peptides are composed of multiple lower order bundles of one to three peptides. For example, the LS2 channel was found to be stable in a tetrameric bundle composed of a “dimer of dimers.” In addition, we observed that there is a minimum of two strong hydrophilic contacts between a pair of adjacent helices in the dimer to tetramer systems and only one strong hydrophilic interhelix contact in helix pairs of the pentamer and hexamer systems. We believe these results are general and can be applied to more complex ion channels, providing insight into ion channel stability and assembly.     

Molecular Dynamics Simulations of Homo-oligomeric Bundles Embedded Within a Lipid Bilayer

Using molecular dynamics simulations, we studied the structure, interhelix interactions, and dynamics of transmembrane proteins. Specifically, we investigated homooligomeric helical bundle systems consisting of synthetic α-helices with either the sequence Ac-(LSLLLSL)₃-NH₂ (LS2) or Ac-(LSSLLSL)₃-NH₂ (LS3). The LS2 and LS3 helical peptides are designed to have amphipathic characteristics that form ion channels in membrane. We simulated bundles containing one to six peptides that were embedded in palmitoyl-oleoyl-phosphatidylcholine (POPC) lipid bilayer and placed between two lamellae of water. We aim to provide a fundamental understanding of how amphipathic helical peptides interact with each other and their dynamical behaviors in different homooligomeric states. To understand structural properties, we examined the helix lengths, tilt angles of individual helices and the entire bundle, interhelix distances, interhelix cross-angles, helix hydrophobic-to-hydrophilic vector projections, and the average number of interhelix hydrophilic (serine–serine) contacts lining the pore of the transmembrane channel. To analyze dynamical properties, we calculated the rotational autocorrelation function of each helix and the cross-correlation of the rotational velocity between adjacent helices. The observed structural and dynamical characteristics show that higher order bundles containing four to six peptides are composed of multiple lower order bundles of one to three peptides. For example, the LS2 channel was found to be stable in a tetrameric bundle composed of a “dimer of dimers.” In addition, we observed that there is a minimum of two strong hydrophilic contacts between a pair of adjacent helices in the dimer to tetramer systems and only one strong hydrophilic interhelix contact in helix pairs of the pentamer and hexamer systems. We believe these results are general and can be applied to more complex ion channels, providing insight into ion channel stability and assembly.     

Interhelical hydrogen bonding drives strong interactions in membrane proteins

Polar residues in transmembrane alpha-helices may strongly influence the folding or association of integral membrane proteins. To test whether a motif that promotes helix association in a soluble protein could do the same within a membrane, we designed a model transmembrane helix based on the GCN4 leucine zipper. We found in both detergent micelles and biological membranes that helix association is driven strongly by asparagine, independent of the rest of the hydrophobic leucine and/or valine sequence. Hydrogen bonding between membrane helices gives stronger associations than the packing of surfaces in glycophorin A helices, creating an opportunity to stabilize structures, but also implying a danger that non-specific interactions might occur. Thus, membrane proteins may fold to avoid exposure of strongly hydrogen bonding groups at their lipid exposed surfaces.