Oligomerization of the fifth transmembrane domain from the adenosine A2A receptor

The human adenosine A_2A receptor (A_2AR) belongs to one of the largest family of membrane proteins, the G-protein coupled receptors (GPCRs), characterized by seven transmembrane (TM) helices. Little is known about the determinants of their structures, folding, assembly, activation mechanisms, and oligomeric states. Previous studies in our group showed that peptides corresponding to all seven TM domains form stable helical structures in detergent micelles and lipid vesicles. However, the peptides behave differently; TM5 is the only peptide to have a ratio [θ]₂₂₂/[θ]₂₀₈ obtained by circular dichroism (CD) spectroscopy>1. This finding suggested to us that TM5 might self-associate. In the present study, we investigate the unique properties of the TM5 domain. We performed detailed analyses of TM5 peptide behavior in membrane-mimetic environments using CD spectroscopy, fluorescence spectroscopy and Förster resonance energy transfer, and gel electrophoresis. We find that TM5 peptide has the ability to self-associate to form oligomeric structures in various hydrophobic milieus and that these oligomers are highly resistant to temperature and chemical denaturation. We also find that mutation of the full-length A_2AR at position M193, which is located in the fifth TM domain, noticeably alters A_2AR monomer: dimer ratio as observed on SDS-PAGE. Our results suggest that parallel association of TM5 dimers may play a role in the known adenosine A_2A receptor dimerization. This study represents the first evidence of an individual GPCR transmembrane domain self-association.     

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     

Conformation of di-n-propylglycine residues (Dpg) in peptides: crystal structures of a type I' beta-turn forming tetrapeptide and an alpha-helical tetradecapeptide.

The crystal structures of two oligopeptides containing di-n-propylglycine (Dpg) residues, Boc-Gly-Dpg-Gly-Leu-OMe (1) and Boc-Val-Ala-Leu-Dpg-Val-Ala-Leu-Val-Ala-Leu-Dpg-Val-Ala-Leu-OMe (2) are presented. Peptide 1 adopts a type I'beta-turn conformation with Dpg(2)-Gly(3) at the corner positions. The 14-residue peptide 2 crystallizes with two molecules in the asymmetric unit, both of which adopt alpha-helical conformations stabilized by 11 successive 5 --> 1 hydrogen bonds. In addition, a single 4 --> 1 hydrogen bond is also observed at the N-terminus. All five Dpg residues adopt backbone torsion angles (phi, psi) in the helical region of conformational space. Evaluation of the available structural data on Dpg peptides confirm the correlation between backbone bond angle N-C(alpha)-C' (tau) and the observed backbone phi,psi values. For tau > 106 degrees, helices are observed, while fully extended structures are characterized by tau < 106 degrees. The mean tau values for extended and folded conformations for the Dpg residue are 103.6 degrees +/- 1.7 degrees and 109.9 degrees +/- 2.6 degrees, respectively.     

Sequence dependence of BNIP3 transmembrane domain dimerization implicates side-chain hydrogen bonding and a tandem GxxxG motif in specific helix-helix interactions

The transmembrane domain of the pro-apoptotic protein BNIP3 self-associates strongly in membranes and in detergents. We have used site-directed mutagenesis to analyze the sequence dependence of BNIP3 transmembrane domain dimerization, from which we infer the physical basis for strong and specific helix-helix interactions in this system. Hydrophobic substitutions identify six residues as critical to dimerization, and the pattern of sensitive residues suggests that the BNIP3 helices interact at a right-handed crossing angle. Based on the dimerization propensities of single point mutants, we propose that: polar residues His173 and Ser172 make inter-monomer hydrogen bonds to one another through their side-chains; Ala176, Gly180, and Gly184 form a tandem GxxxG motif that allows close approach of the helices; and Ile183 makes inter-monomer van der Waals contacts. Since neither the tandem GxxxG motif nor the hydrogen bonding pair is sufficient to drive dimerization, our results demonstrate the importance of sequence context for either hydrogen bonding or GxxxG motif involvement in BNIP3 transmembrane helix-helix interactions. In this study, hydrophobic substitutions away from the six interfacial positions have almost no effect on dimerization, confirming the expectation that hydrophobic replacements affect helix-helix interactions only if they interfere with packing or hydrogen bonding by interfacial residues. However, changes to slightly polar residues are somewhat disruptive even when located away from the interface, and the degree of disruption correlates with the decrease in hydrophobicity. Changing the hydrophobicity of the BNIP3 transmembrane domain alters its helicity and protection of its backbone amides. We suggest that polar substitutions decrease the fraction of dimer by stabilizing an unfolded monomeric state of the transmembrane span, rather than by affecting helix-helix interactions. This result has broad implications for interpreting the sequence dependence of membrane protein stability in detergents.     

Computer modeling of polyleucine-based coiled coil dimers in a realistic membrane environment: insight into helix-helix interactions in membrane proteins

Simulated annealing was performed to model parallel dimers of alpha-helical transmembrane peptides with the sequence L(11)XL(12), predicting left-handed coiled coil geometry in all cases. Insertion of peptides containing threonine, asparagine, alanine, phenylalanine, and leucine in position 12 into realistic model membranes showed these structures were stable for 20 ns of molecular dynamics simulation time. Threonine could participate in intermolecular hydrogen bonds, but predominantly formed hydrogen bonds to the backbone of the helix it resided on. These hydrogen bonds, although infrequent, appeared to promote closer association of polyleucine helices. Asparagine participated in multiple, rapidly fluctuating intermolecular and intramolecular hydrogen bonds, and may have slightly destabilized optimum van der Waals packing in favor of optimum hydrogen bonding. Coordinated rotations of transmembrane helices about their axes were observed, indicating helices may rotate around one another during the folding of membrane proteins or other processes. These rotations were inhibited by phenylalanine, suggesting a role for bulky residues in modulating membrane protein dynamics.     

Structural analysis of peptide helices containing centrally positioned lactic acid residues

The effect of insertion of lactic acid (Lac) residues into peptide helices has been probed using specifically designed sequences. The crystal structures of 11-residue and 14-residue depsipeptides Boc-Val-Val-Ala-Leu-Val-Lac-Leu-Aib-Val-Ala-Leu-OMe (1) and Boc-Val-Ala-Leu-Aib-Val-Ala-Leu-Val-Lac-Leu-Aib-Val-Ala-Leu-OMe (3), containing centrally positioned Lac residues, have been determined. The structure of an 11-residue peptide Boc-Val-Ala-Leu-Aib-Val-Ala-Leu-Aib-Val-Ala-Leu-OMe (2), analog of a which is an amide previously determined Lac-containing depsipeptide, Boc-Val-Ala-Leu-Aib-Val-Lac-Leu-Aib-Val-Ala-Leu-OMe I. L. Karle, C. Das, and P. Balaram, Biopolymers, Vol. 59, (2001) pp. 276-289], is also reported. Peptide 1 adopts a helical fold, which is stabilized by mixture of 4-->1 and 5-->1 hydrogen bonds. Peptide 2 adopts a completely alpha-helical conformation stabilized by eight successive 5-->1 hydrogen bonds. Peptide 3 appears to be predominately alpha-helical, with seven 5-->1 hydrogen bonds and three 4-->1 interaction interspersed in the sequence. In the structure of peptide 3 in addition to water molecules in the head-to-tail region, hydration at an internal segment of the helix is also observed. A comparison of five related peptide helices, containing a single Lac residue, reveals that the hydroxy acid can be comfortably accommodated at interior positions in the helix, with the closest C=O...O distances lying between 2.8 and 3.3 A.     

Helix packing of leucine-rich peptides: a parallel leucine ladder in the structure of Boc-Aib-Leu-Aib-Aib-Leu-Leu-Leu-Aib-Leu-Aib-OMe.

The packing of peptide helices in crystals of the leucine-rich decapeptide Boc-Aib-Leu-Aib-Aib-Leu-Leu-Leu-Aib-Leu-Aib-OMe provides an example of ladder-like leucylleucyl interactions between neighboring molecules. The peptide molecule forms a helix with five 5----1 hydrogen bonds and two 4----1 hydrogen bonds near the C terminus. Three head-to-tail NH ... O = C hydrogen bonds between helices form continuous columns of helices in the crystal. The helicial columns associate in an antiparallel fashion, except for the association of Leu ... Leu side chains, which occurs along the diagonal of the cell where the peptide helices are parallel. The peptide, with formula C56H102N10O13, crystallizes in space group P2(1)2(1)2(1) with Z = 4 and cell parameters a = 16.774(3) A, b = 20.032(3) A and c = 20.117(3) A; overall agreement factor R = 10.7% for 2014 data with magnitude of F(obs) greater than 3 sigma (F); resolution 1.0 A.     

Specific heteromeric association of four transmembrane peptides derived from platelet glycoprotein Ib-IX complex

As the receptor on the platelet surface for von Willebrand factor, glycoprotein (GP) Ib-IX complex is critically involved in hemostasis and thrombosis. How the complex is assembled from GP Ibα, GP Ibβ and GP IX subunits, all of which are type I transmembrane proteins, is not entirely clear. Genetic and mutational analyses have identified the transmembrane (TM) domains of these subunits as active participants in assembly of the complex. In this study, peptides containing the transmembrane domain of each subunit have been produced and their interaction with one another characterized. Only the Ibβ TM sequence, but not the Ibα and IX counterparts, can form homo-oligomers in SDS-PAGE and TOXCAT assays. Following up on our earlier observation that a Ibβ-Ibα-Ibβ peptide complex (αβ₂) linked through native juxtamembrane disulfide bonds could be produced from isolated Ibα and Ibβ TM peptides in detergent micelles, we show here that addition of the IX TM peptide facilitates formation of the native αβ₂ complex, reproducing the same effect by the IX subunit in cells expressing the GP Ib-IX complex. Specific fluorescence resonance energy transfer was observed between donor-labeled αβ₂ peptide complex and acceptor-conjugated IX TM peptide in micelles. Finally, the mutation D135K in the IX TM peptide could hamper both the formation of the αβ₂ complex and the energy transfer, consistent with its reported effect in the full-length complex. Overall, our results have demonstrated directly the native-like heteromeric interaction among the isolated Ibα, Ibβ and IX TM peptides, which provides support for the four-helix bundle model of the TM domains in the GP Ib-IX complex and paves the way for further structural analysis. The methods developed in this study may be applicable to other studies of heteromeric interaction among multiple TM helices.     

The determinants of hydrophobic mismatch response for transmembrane helices

Hydrophobic mismatch arises from a difference in the hydrophobic thickness of a lipid membrane and a transmembrane protein segment, and is thought to play an important role in the folding, stability and function of membrane proteins. We have investigated the possible adaptations that lipid bilayers and transmembrane α-helices undergo in response to mismatch, using fully-atomistic molecular dynamics simulations totaling 1.4 μs. We have created 25 different tryptophan-alanine-leucine transmembrane α-helical peptide systems, each composed of a hydrophobic alanine–leucine stretch, flanked by 1–4 tryptophan side chains, as well as the β-helical peptide dimer, gramicidin A. Membrane responses to mismatch include changes in local bilayer thickness and lipid order, varying systematically with peptide length. Adding more flanking tryptophan side chains led to an increase in bilayer thinning for negatively mismatched peptides, though it was also associated with a spreading of the bilayer interface. Peptide tilting, bending and stretching were systematic, with tilting dominating the responses, with values of up to ~ 45° for the most positively mismatched peptides. Peptide responses were modulated by the number of tryptophan side chains due to their anchoring roles and distributions around the helices. Potential of mean force calculations for local membrane thickness changes, helix tilting, bending and stretching revealed that membrane deformation is the least energetically costly of all mismatch responses, except for positively mismatched peptides where helix tilting also contributes substantially. This comparison of energetic driving forces of mismatch responses allows for deeper insight into protein stability and conformational changes in lipid membranes.     

Intermonomer Hydrogen Bonds Enhance GxxxG-Driven Dimerization of the BNIP3 Transmembrane Domain: Roles for Sequence Context in Helix-Helix Association in Membranes

We determined the sequence dependence of human BNIP3 transmembrane domain dimerization using the biological assay TOXCAT. Mutants in which intermonomer hydrogen bonds between Ser172 and His173 are abolished show moderate interaction, indicating that side-chain hydrogen bonds contribute to dimer stability but are not essential to dimerization. Mutants in which a GxxxG motif composed of Gly180 and Gly184 has been abolished show little or no interaction, demonstrating the critical nature of the GxxxG motif to BNIP3 dimerization. These findings show that side-chain hydrogen bonds can enhance the intrinsic dimerization of a GxxxG motif and that sequence context can control how hydrogen bonds influence helix-helix interactions in membranes. The dimer interface mapped by TOXCAT mutagenesis agrees closely with the interfaces observed in the NMR structure and inferred from mutational analysis of dimerization on SDS-PAGE, showing that the native dimer structure is retained in detergents. We show that TOXCAT and SDS-PAGE give complementary and consistent information about BNIP3 transmembrane domain dimerization: TOXCAT is insensitive to mutations that have modest effects on self-association in detergents but readily discriminates among mutations that completely disrupt detergent-resistant dimerization. The close agreement between conclusions reached from TOXCAT and SDS-PAGE data for BNIP3 suggests that accurate estimates of the relative effects of mutations on native-state protein-protein interactions can be obtained even when the detergent environment is strongly disruptive.     

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     

The structure of holo and metal-deficient wild-type human Cu,Zn superoxide dismutase and its relevance to familial amyotrophic lateral sclerosis

Cu, Zn superoxide dismutase (SOD1) forms a crucial component of the cellular defence against oxidative stress. Zn-deficient wild-type and mutant human SOD1 have been implicated in the disease familial amyotrophic lateral sclerosis (FALS). We present here the crystal structures of holo and metal-deficient (apo) wild-type protein at 1.8A resolution. The P21 wild-type holo enzyme structure has nine independently refined dimers and these combine to form a "trimer of dimers" packing motif in each asymmetric unit. There is no significant asymmetry between the monomers in these dimers, in contrast to the subunit structures of the FALS G37R mutant of human SOD1 and in bovine Cu,Zn SOD. Metal-deficient apo SOD1 crystallizes with two dimers in the asymmetric unit and shows changes in the metal-binding sites and disorder in the Zn binding and electrostatic loops of one dimer, which is devoid of metals. The second dimer lacks Cu but has approximately 20% occupancy of the Zn site and remains structurally similar to wild-type SOD1. The apo protein forms a continuous, extended arrangement of beta-barrels stacked up along the short crystallographic b-axis, while perpendicular to this axis, the constituent beta-strands form a zig-zag array of filaments, the overall arrangement of which has a similarity to the common structure associated with amyloid-like fibrils.     

Synthesis and characterization of β‐peptide helices as transmembrane domains in lipid model membranes

Aggregation, orientation and dynamics of transmembrane helices are of relevance for protein function and transmembrane signaling. To explore the interactions of transmembrane helices and the interdependence of peptide structure and lipid composition of the membranes, β‐peptides were explored as model transmembrane domains. Various hydrophobic β‐peptide sequences were synthesized by solid phase peptide synthesis. Conformational analyses of β‐peptide helices were performed in organic solvents (methanol and 2,2,2‐trifluoroethanol) and in large unilamellar liposomes (dimyristoylphosphatidylcholine, dipalmitoylphosphatidylcholine and dioleoylphosphatidylcholine) indicating 12‐ and 14‐helix conformations, depending on β³‐amino acid sequences. The intrinsic tryptophan fluorescence of β³‐homotryptophan units inserted in the center or near the end of the sequence was used to verify the membrane insertion of the β‐peptides. A characteristic blue shift with peripheral β³‐homotryptophan compared with β‐peptides with central tryptophan served as indication for a transmembrane orientation of the β‐peptides within the lipid bilayer. Copyright © 2016 European Peptide Society and John Wiley & Sons, Ltd.     

Effects of branched beta-carbon dehydro-residues on peptide conformations: syntheses, crystal structures and molecular conformations of two tetrapeptides: (a) N-(benzyloxycarbonyl)-DeltaVal-Leu-DeltaPhe-Leu-OCH3 and (b) N-(benzyloxycarbonyl)-DeltaIle-Ala-DeltaPhe-Ala-OCH3.

The roles of branched beta-carbon dehydro-residues in the design of peptide conformations have not been systematically explored so far. In order to determine the effects of branched beta-carbon dehydro-residues on the peptide conformations, two N-protected tetrapeptides containing new combinations of DeltaVal and DeltaPhe in (a) N-(benzyloxycarbonyl)-DeltaVal-Leu-DeltaPhe-Leu-OCH(3) and DeltaIle and DeltaPhe in (b) N-(benzyloxycarbonyl)-DeltaIle-Ala-DeltaPhe-Ala-OCH(3) were synthesized by solution procedure. The crystal structures of these peptides were determined by X-ray diffraction methods. Single crystals of both peptides were grown by slow evaporation method from their solutions in acetone-water mixtures (80 : 20) at 25 degrees C. The crystals of these peptides belong to the orthorhombic space group P2(1)2(1)2(1) with cell dimensions of a = 12.342(1) A, b = 15.659(1) A, c = 18.970(1) A for peptide (a) and a = 8.093(1) A, b = 15.791(1) A, c = 23.816(1) A for peptide (b) having Z = 4 in the unit cells of both peptides. The structures were refined by full-matrix least-squares procedure to R-factors of 0.076 and 0.052 respectively. Both peptides adopt the right-handed 3(10)-helical conformations stabilized by two intramolecular (i + 3-->i) hydrogen bonds between the CO of N-terminal benzyloxycarbonyl (Cbz) group and the NH of residue at position 3, and between the CO of residue at position 1 and NH of the residue at position 4. The two consecutive 10-membered rings formed by the hydrogen bonds have dihedral angles corresponding to the standard values for type III beta-turns. DeltaVal and DeltaIle in peptides (a) and (b) respectively are located at the (i + 1) position of the first beta-turn while DeltaPhe is located at the (i + 2) position of the second beta-turn. In the crystals, the molecules are linked head to tail by intermolecular hydrogen bonds to form long helical chains. The axes of helices are parallel to the b-axes while the neighbouring helices run in the opposite directions. The crystal packings are further stabilized by van der Waals forces between the columns of molecular packings.     

Rotational orientation of monomers within a designed homo-oligomer transmembrane helical bundle

A peptide designed to form a homo-oligomeric transmembrane helical bundle was reconstituted into lipid bilayers and studied by using (2)H NMR (nuclear magnetic resonance) with magic angle spinning to confirm that the helical interface corresponds to the interface intended in the design. The peptide belongs to a family of model peptides derived from a membrane-solubilized version of the water-soluble coiled-coil GCN4-P1. The variant studied here contains two asparagines thought to engage in interhelical hydrogen bonding critical to the formation of a stable trimer. For the NMR studies, three different peptides were synthesized, each with one of three consecutive leucines in the transmembrane region deuterium labeled. Prior to NMR data collection, polarized infrared spectroscopy was used to establish that the peptides were reconstituted in lipid bilayers in a transmembrane helical conformation. The (2)H NMR line shapes of the three different peptides are consistent with a trimer structure formed by the designed peptide that is stabilized by inter-helical hydrogen bonding of asparagines at positions 7 and 14.     

Structures of the intradiskal loops and amino terminus of the G-protein receptor, rhodopsin.

The intradiskal surface of the transmembrane protein, rhodopsin, consists of the amino terminal domain and three loops connecting six of the seven transmembrane helices. This surface corresponds to the extracellular surface of other G-protein receptors. Peptides that represent each of the extramembraneous domains on this surface (three loops and the amino terminus) were synthesized. These peptides also included residues which, based on a hydrophobic plot, could be expected to be part of the transmembrane helix. The structure of each of these peptides in solution was then determined using two-dimensional 1H nuclear magnetic resonance. All peptide domains showed ordered structures in solution. The structures of each of the peptides from intradiskal loops of rhodopsin exhibited a turn in the central region of the peptide. The ends of the peptides show an unwinding of the transmembrane helices to form this turn. The amino terminal domain peptide exhibited alpha-helical regions with breaks and bends at proline residues. This region forms a compact domain. Together, the structures for the loop and amino terminus domains indicate that the intradiskal surface of rhodopsin is ordered. These data further suggest a structural motif for short loops in transmembrane proteins. The ordered structures of these loops, in the absence of the transmembrane helices, indicate that the primary sequences of these loops are sufficient to code for the turn.     

Characterization of a membrane protein folding motif, the Ser zipper, using designed peptides

Polar residues play important roles in the association of transmembrane helices and the stabilities of membrane proteins. Although a single Ser residue in a transmembrane helix is unable to mediate a strong association of the helices, the cooperative interactions of two or more appropriately placed serine hydroxyl groups per helix has been hypothesized to allow formation of a "serine zipper" that can stabilize transmembrane helix association. In particular, a heptad repeat Sera Xxx Xxx Leud Xxx Xxx Xxx (Xxx is a hydrophobic amino acid) appears in both antiparallel helical pairs of polytopic membrane proteins as well as the parallel helical dimerization motif found in the murine erythropoietin receptor. To examine the intrinsic conformational preferences of this motif independent of its context within a larger protein, we synthesized a peptide containing three copies of a SeraLeud heptad motif. Computational results are consistent with the designed peptide adopting either a parallel or antiparallel structure, and conformational search calculations yield the parallel dimer as the lowest energy configuration, which is also significantly more stable than the parallel trimer. Analytical ultracentrifugation indicated that the peptide exists in a monomer-dimer equilibrium in dodecylphosphocholine micelles. Thiol disulfide interchange studies showed a preference for forming parallel dimers in micelles. In phospholipid vesicles, only the parallel dimer was formed. The stability of the SerZip peptide was studied in vesicles prepared from phosphatidylcholine (PC) lipids of different chain length: POPC (C16:0C18:1 PC) and DLPC (C12:0PC). The stability was greater in POPC, which has a good match between the length of the hydrophobic region of the peptide and the bilayer length. Finally, mutation to Ala of the Ser residues in the SerZip motif gave rise to a relatively small decrease in the stability of the dimer, indicating that packing interactions rather than hydrogen-bonding provided the primary driving force for association.     

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.     

Helix packing motif common to the crystal structures of two undecapeptides containing dehydrophenylalanine residues: implications for the de novo design of helical bundle super secondary structural modules

De novo designed peptide based super secondary structures are expected to provide scaffolds for the incorporation of functional sites as in proteins. Self-association of peptide helices of similar screw sense, mediated by weak interactions, has been probed by the crystal structure determination of two closely related peptides: Ac-Gly1-Ala2-Delta Phe3-Leu4-Val5-DeltaPhe6-Leu7-Val8-DeltaPhe9-Ala10-Gly11-NH2 (I) and Ac-Gly1-Ala2-DeltaPhe3-Leu4-Ala5-DeltaPhe6-Leu7-Ala8-DeltaPhe9-Ala10-Gly11-NH2 (II). The crystal structures determined to atomic resolution and refined to R factors 8.12 and 4.01%, respectively, reveal right-handed 3(10)-helical conformations for both peptides. CD has also revealed the preferential formation of right-handed 3(10)-helical conformations for both molecules. Our aim was to critically analyze the packing of the helices in the solid state with a view to elicit clues for the design of super secondary structural motifs such as two, three, and four helical bundles based on helix-helix interactions. An important finding is that a packing motif could be identified common to both the structures, in which a given peptide helix is surrounded by six other helices reminiscent of transmembrane seven helical bundles. The outer helices are oriented either parallel or antiparallel to the central helix. The helices interact laterally through a combination of N--H...O, C--H...O, and C--H...pi hydrogen bonds. Layers of interacting leucine residues are seen in both peptide crystal structures. The packing of the peptide helices in the solid state appears to provide valuable leads for the design of super secondary structural modules such as two, three, or four helix bundles by connecting adjacent antiparallel helices through suitable linkers such as tetraglycine segments.