Molecular packing and packing defects in helical membrane proteins

The packing of helices spanning lipid bilayers is crucial for the stability and function of alpha-helical membrane proteins. Using a modified Voronoi procedure, we calculated packing densities for helix-helix contacts in membrane spanning domains. Our results show that the transmembrane helices of protein channels and transporters are significantly more loosely packed compared with helices in globular proteins. The observed packing deficiencies of these membrane proteins are also reflected by a higher amount of cavities at functionally important sites. The cavities positioned along the gated pores of membrane channels and transporters are noticeably lined by polar amino acids that should be exposed to the aqueous medium when the protein is in the open state. In contrast, nonpolar amino acids surround the cavities in those protein regions where large rearrangements are supposed to take place, as near the hinge regions of transporters or at restriction sites of protein channels. We presume that the observed deficiencies of helix-helix packing are essential for the helical mobility that sustains the function of many membrane protein channels and transporters.     

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.     

Self-association of transmembrane alpha-helices in model membranes: importance of helix orientation and role of hydrophobic mismatch

Interactions between transmembrane helices play a key role in almost all cellular processes involving membrane proteins. We have investigated helix-helix interactions in lipid bilayers with synthetic tryptophan-flanked peptides that mimic the membrane spanning parts of membrane proteins. The peptides were functionalized with pyrene to allow the self-association of the helices to be monitored by pyrene fluorescence and Trp-pyrene fluorescence resonance energy transfer (FRET). Specific labeling of peptides at either their N or C terminus has shown that helix-helix association occurs almost exclusively between antiparallel helices. Furthermore, computer modeling suggested that antiparallel association arises primarily from the electrostatic interactions between alpha-helix backbone atoms. We propose that such interactions may provide a force for the preferentially antiparallel association of helices in polytopic membrane proteins. Helix-helix association was also found to depend on the lipid environment. In bilayers of dioleoylphosphatidylcholine, in which the hydrophobic length of the peptides approximately matched the bilayer thickness, association between the helices was found to require peptide/lipid ratios exceeding 1/25. Self-association of the helices was promoted by either increasing or decreasing the bilayer thickness, and by adding cholesterol. These results indicate that helix-helix association in membrane proteins can be promoted by unfavorable protein-lipid interactions.     

How do helix–helix interactions help determine the folds of membrane proteins? Perspectives from the study of homo‐oligomeric helical bundles

The final, structure-determining step in the folding of membrane proteins involves the coalescence of preformed transmembrane helices to form the native tertiary structure. Here, we review recent studies on small peptide and protein systems that are providing quantitative data on the interactions that drive this process. Gel electrophoresis, analytical ultracentrifugation, and fluorescence resonance energy transfer (FRET) are useful methods for examining the assembly of homo-oligomeric transmembrane helical proteins. These methods have been used to study the assembly of the M2 proton channel from influenza A virus, glycophorin, phospholamban, and several designed membrane proteins-all of which have a single transmembrane helix that is sufficient for association into a transmembrane helical bundle. These systems are being studied to determine the relative thermodynamic contributions of van der Waals interactions, conformational entropy, and polar interactions in the stabilization of membrane proteins. Although the database of thermodynamic information is not yet large, a few generalities are beginning to emerge concerning the energetic differences between membrane and water-soluble proteins: the packing of apolar side chains in the interior of helical membrane proteins plays a smaller, but nevertheless significant, role in stabilizing their structure. Polar, hydrogen-bonded interactions occur less frequently, but, nevertheless, they often provide a strong driving force for folding helix-helix pairs in membrane proteins. These studies are laying the groundwork for the design of sequence motifs that dictate the association of membrane helices.     

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.     

Analysis and prediction of helix-helix interactions in membrane channels and transporters

Membrane proteins span a large variety of different functions such as cell-surface receptors, redox proteins, ion channels, and transporters. Proteins with functional pores show different characteristics of helix-helix packing as other helical membrane proteins. We found that the helix-helix contacts of 13 nonhomologous high-resolution structures of membrane channels and transporters are mainly accomplished by weakly polar amino acids (G > S > T > F) that preferably create contacts every fourth residue, typical for right-handed helix crossings. There is a strong correlation between the now available biological hydrophobicity scale and the propensities of the weakly polar and hydrophobic residues to be buried at helix-helix interfaces or to be exposed to the lipids in membrane channels and transporters. The polar residues, however, make no major contribution towards the packing of their transmembrane helices, and are therefore subsumed to be primarily exposed to the polar milieu during the folding process. The contact formation of membrane channels and transporters is therefore ruled by the solubility of the residues, which we suppose to be the driving force for the assembly of their transmembrane helices. By contrast, in 14 nonhomologous high-resolution structures of other membrane protein coils, also large and polar amino acids (D > S > M > Q) create characteristic contacts every 3.5th residues, which is a signature for left-handed helix crossings. Accordingly, it seems that dependent on the function, different concepts of folding and stabilization are realized for helical membrane proteins. Using a sequence-based matrix prediction method these differences are exploited to improve the prediction of buried and exposed residues of transmembrane helices significantly. When the sequence motifs typical for membrane channels and transporters were applied for the prediction of helix-helix contacts the quality of prediction rises by 16% to an average value of 76%, compared to the same approach when only single amino acid positions are taken into account.     

Understanding peptide interactions with the lipid bilayer: a guide to membrane protein engineering

In recent years, studies of the interactions of designed hydrophobic and amphipathic polypeptides with biological membranes have progressed considerably. In-plane and transmembrane helical domains have been engineered as well as sequences that exhibit dynamic distributions of different topologies. These sequences not only help our understanding of the thermodynamic interaction contributions that determine peptide orientation, but also exhibit interesting biological functions as autonomous units. In addition, helices are considered to be folding intermediates of multi-spanning membrane proteins. The specificity and thermodynamic stability of transmembrane helix-helix interactions or the assembly of beta sheets at the membrane surface have, therefore, become a focus of ongoing research.     

Sequence specificity in the dimerization of transmembrane alpha-helices.

While several reports have suggested a role for helix-helix interactions in membrane protein oligomerization, there are few direct biochemical data bearing on this subject. Here, using mutational analysis, we show that dimerization of the transmembrane alpha-helix of glycophorin A in a detergent environment is spontaneous and highly specific. Very subtle changes in the side-chain structure at certain sensitive positions disrupt the helix-helix association. These sensitive positions occur at approximately every 3.9 residues along the helix, consistent with their comprising the interface of a closely fit transmembranous supercoil of alpha-helices. By contrast with other reported cases of interactions between transmembrane helices, the set of interfacial residues in this case contains no highly polar groups. Amino acids with aliphatic side chains define much of the interface, indicating that precise packing interactions between the helices may provide much of the energy for association. These data highlight the potential general importance of specific interactions between the hydrophobic anchors of integral membrane proteins.     

Simulation studies on bacteriorhodopsin bundle of transmembrane α segments

Bacteriorhodopsin (BR) is a membrane protein which pumps protons through the plasma membrane. Seven transmembrane BR helical segments are subjected to simulation studies in order to investigate the packing process of transmembrane helices. A Monte Carlo simulated annealing protocol is employed to optimize the helix bundle system. Helix packing is optimized according to a semi-empirical potential mainly composed of six components: a bilayer potential, a crossing angle potential, a helix dipole potential, a helix-helix distance potential, a helix orientation potential and a helix-helix distance restraint potential (a loop potential). Necessary parameters are derived from theoretical studies and statistical analysis of experimentally determined protein structures. The structures from the simulations are compared with the experimentally determined structures in terms of geometry. The structures generated show similar shapes to the experimentally suggested structure even without the helix-helix distance restraint potential. However, the relative locations of individual helices were reproduced only when the helix-helix distance restraint potential was used with restraint conditions. Our results suggest that transmembrane helix bundles resembling those observed experimentally may be generated by simulations using simple potentials. Received: 19 April 1999 / Revised version: 6 September 1999 / Accepted: 17 September 1999     

The membrane environment modulates self-association of the human GpA TM domain—Implications for membrane protein folding and transmembrane signaling

The influence of lipid bilayer properties on a defined and sequence-specific transmembrane helix-helix interaction is not well characterized yet. To study the potential impact of changing bilayer properties on a sequence-specific transmembrane helix-helix interaction, we have traced the association of fluorescent-labeled glycophorin A transmembrane peptides by fluorescence spectroscopy in model membranes with varying lipid compositions. The observed changes of the glycophorin A dimerization propensities in different lipid bilayers suggest that the lipid bilayer thickness severely influences the monomer-dimer equilibrium of this transmembrane domain, and dimerization was most efficient under hydrophobic matching conditions. Moreover, cholesterol considerably promotes self-association of transmembrane helices in model membranes by affecting the lipid acyl chain ordering. In general, the order of the lipid acyl chains appears to be an important factor involved in determining the strength and stability of transmembrane helix-helix interactions. As discussed, the described influences of membrane properties on transmembrane helix-helix interactions are highly important for understanding the mechanism of transmembrane protein folding and functioning as well as for gaining a deeper insight into the regulation of signal transduction via membrane integral proteins by bilayer properties.     

Role of the extracellular loop in the folding of a CFTR transmembrane helical hairpin

The folding of membrane-spanning domains into their native functional forms depends on interactions between transmembrane (TM) helices joined by covalent loops. However, the importance of these covalent linker regions in mediating the strength of helix-helix associations has not been systematically addressed. Here we examine the potential structural impact of cystic fibrosis-phenotypic mutations in the extracellular loop 2 (ECL2) on interactions between the TM3 and TM4 helices of the cystic fibrosis transmembrane conductance regulator (CFTR) in constructs containing CFTR residues 194-241. When the effects of replacements in ECL2 (including the CF-phenotypic mutants E217G and Q220R) were evaluated in a library of wild-type and mutant TM3-ECL2-TM4 hairpin constructs, we found that SDS-PAGE gel migration rates differed over a range of nearly 40% +/- the wild-type position and that decreased migration rates correlate with increasing hairpin alpha-helical content as measured by circular dichroism spectra in sodium dodecyl sulfate micelles. The decreased mobility of TM3/4 constructs by introduction of non-native residues is interpreted in terms of an elongation or "opening" of the helical hairpin and concomitant destabilization of membrane-based helix-helix interactions. Our results support a role for short loop regions in dictating the stability of membrane protein folds and highlight the interplay between membrane-embedded helix-helix interactions and loop conformation in influencing the structure of membrane proteins.     

The position of the Gly-xxx-Gly motif in transmembrane segments modulates dimer affinity.

Although the intrinsic low solubility of membrane proteins presents challenges to their high-resolution structure determination, insight into the amino acid sequence features and forces that stabilize their folds has been provided through study of sequence-dependent helix-helix interactions between single transmembrane (TM) helices. While the stability of helix-helix partnerships mediated by the Gly-xxx-Gly (GG4) motif is known to be generally modulated by distal interfacial residues, it has not been established whether the position of this motif, with respect to the ends of a given TM segment, affects dimer affinity. Here we examine the relationship between motif position and affinity in the homodimers of 2 single-spanning membrane protein TM sequences: glycophorin A (GpA) and bacteriophage M13 coat protein (MCP). Using the TOXCAT assay for dimer affinity on a series of GpA and MCP TM segments that have been modified with either 4 Leu residues at each end or with 8 Leu residues at the N-terminal end, we show that in each protein, centrally located GG4 motifs are capable of stronger helix-helix interactions than those proximal to TM helix ends, even when surrounding interfacial residues are maintained. The relative importance of GG4 motifs in stabilizing helix-helix interactions therefore must be considered not only in its specific residue context but also in terms of the location of the interactive surface relative to the N and C termini of alpha-helical TM segments.     

Peptide mimics of the M13 coat protein transmembrane segment. Retention of helix-helix interaction motifs

Sequence-specific noncovalent helix-helix interactions between transmembrane (TM) segments in proteins are investigated by incorporating selected TM sequences into synthetic peptides using the construct CKKK-TM-KKK. The peptides are of suitable hydrophobicity for spontaneous membrane insertion, whereas formation of an N-terminal S-S bond can bring pairs of TM helices into proximity and promote their parallel orientation. Using the propensity of the protein to undergo thermally induced alpha-helix --> beta-sheet transitions as a parameter for helix stability, we compared the wild type and mutant (V29A and V31A) bacteriophage M13 coat proteins with their corresponding TM peptide constructs (M13 residues 24-42). Our results demonstrated that the relevant helix-helix tertiary contacts found in the intact proteins persist in the peptide mimics. Molecular dynamics simulations support the tight "two in-two out" dimerization motif for V31A consistent with mutagenesis data. The overall results reinforce the notion of TM segments as autonomous folding domains and suggest that the generic peptide construct provides a viable reductionist system for membrane protein structural and computational analysis.     

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.     

Defining the structural basis for assembly of a transmembrane cytochrome

To define the structural basis for cofactor binding to membrane proteins, we introduce a manageable model system, which allows us, for the first time, to study the influence of individual transmembrane helices and of single amino acid residues on the assembly of a transmembrane cytochrome. In vivo as well as in vitro analyses indicate central roles of single amino acid residues for either interaction of the transmembrane helices or for binding of the cofactor. The results clearly show that interaction of the PsbF transmembrane helix is independent from binding of the heme cofactor. On the other hand, binding of the cofactor highly depends on helix-helix interactions. By site-directed mutagenesis critical amino acid residues were identified, which are involved in the assembly of a functional transmembrane cytochrome. Especially, a highly conserved glycine residue is critical for interaction of the transmembrane helices and assembly of the cytochrome. Based on the two-stage-model of alpha-helical membrane protein folding, the presented results clearly indicate a third stage of membrane protein folding, in which a cofactor binds to a pre-assembled transmembrane protein.     

Discovery of a significant,nontopological preference for antiparallel alignment of helices with parallel regions in sheets

To help elucidate the role of secondary structure packing preferences in protein folding, here we present an analysis of the packing geometry observed between alpha-helices and between alpha-helices and beta-sheets in 1316 diverse, nonredundant protein structures. Finite-length vectors were fit to the alpha-carbon atoms in each of the helices and strands, and the packing angle between the vectors, Omega, was determined at the closest point of approach within each helix-helix or helix-sheet pair. Helix-sheet interactions were found in 391 of the proteins, and the distributions of Omega values were calculated for all the helix-sheet and helix-helix interactions. The packing angle preferences for helix-helix interactions are similar to those previously observed. However, analysis of helix-strand packing preferences uncovered a remarkable tendency for helices to align antiparallel to parallel regions of beta-sheets, independent of the topological constraints or prevalence of beta-alpha-beta motifs in the proteins. This packing angle preference is significantly diminished in helix interactions involving mixed and antiparallel beta-sheets, suggesting a role for helix-sheet dipole alignment in guiding supersecondary structure formation in protein folding. This knowledge of preferred packing angles can be used to guide the engineering of stable protein modules.     

Energetics and stability of transmembrane helix packing: a replica-exchange simulation with a knowledge-based membrane potential

The energetics and stability of the packing of transmembrane helices were investigated by Monte Carlo simulations with the replica-exchange method. The helices were modeled with a united atom representation, and the CHARMM19 force field was employed. Based on known experimental structures of membrane proteins, an implicit knowledge-based potential was developed to describe the helix-membrane interactions at the residue level, whose validity was tested through prediction of the orientations when single helices were inserted into a membrane. Two systems were studied in this article, namely the glycophorin A dimer, and helices A and B of Bacteriorhodopsin. For the glycophorin A dimer, the most stable structure (0.5 A away from the experimental structure) is mainly stabilized by the favorable helix-helix interactions, and has the most population regardless of the helix-membrane interaction. However, for helices A and B of Bacteriorhodopsin, it was found that the packing determined by helix-helix interactions is nonspecific, and a native-like structure (0.2 A from the experimental one) can be identified from several structural analogs as the most stable one only after applying the membrane potential. Our results suggest that the contribution from the helix-membrane interaction could be critical in the correct packing of transmembrane helices in the membrane.     

Shorter side chains optimize helix-helix packing

A systematic study of helix-helix packing in a comprehensive database of protein structures revealed that the side chains inside helix-helix interfaces on average are shorter than those in the noninterface parts of the helices. The study follows our earlier study of this effect in transmembrane helices. The results obtained on the entire database of protein structures are consistent with those obtained on the transmembrane helices. The difference in the length of interface and noninterface side chains is small but statistically significant. It indicates that helices, if viewed along their main axis, statistically are not circular, but have a flattened interface. This effect brings the helices closer to each other and creates a tighter structural packing. The results provide an interesting insight into the aspects of protein structure and folding.     

Distinct protein interfaces in transmembrane domains suggest an in vivo folding model

Given the known high-resolution structures of alpha-helical transmembrane domains, we show that there are statistically distinct classes of transmembrane interfaces which relate to the folding and oligomerization of transmembrane domains. Distinct types of interfaces have been categorized and refer to those between: the same polypeptide chain, different polypeptide chains, helices that are sequential neighbors, and those that are nonsequential. These different interfaces may reflect different phases in the mechanism of transmembrane domain folding and are consistent with the current experimental evidence pertaining to the folding and oligomerization of transmembrane domains. The classes of helix-helix interfaces have been identified in terms of the numbers and different types of pairwise amino acid interactions. The specific measures used are interaction entropy, the information content of interacting partners compared to a random set of contacts, the amino acid composition of the classes and the abundances of specific amino acid pairs in close contact. Knowledge of the clear differences in the types of helix-helix contacts helps with the derivation of knowledge-based constraints which until now have focused on only the interiors of transmembrane domains as compared to the exterior. Taken together, an in vivo model for membrane protein folding is presented, which is distinct from the familiar two-stage model. The model takes into account the different interfaces of membrane helices defined herein, and the available data regarding folding in the translocation channel.     

The stability of transmembrane helix interactions measured in a biological membrane

Despite some promising progress in the understanding of membrane protein folding and assembly, there is little experimental information regarding the thermodynamic stability of transmembrane helix interactions and even less on the stability of transmembrane helix-helix interactions in a biological membrane. Here we describe an approach that allows quantitative measurement of transmembrane helix interactions in a biological membrane, and calculation of changes in the interaction free energy resulting from substitution of single amino acids. Dimerization of several variants of the glycophorin A transmembrane domain are characterized and compared to the wild-type (wt) glycophorin A transmembrane helix dimerization. The calculated DeltaDeltaG(app) values are further compared with values found in the literature. In addition, we compare interactions between the wt glycophorin A transmembrane domain and helices in which critical glycine residues are replaced by alanine or serine, respectively. The data demonstrate that replacement of the glycine residues by serine is less destabilizing than replacement by alanine with a DeltaDeltaG(app) value of about 0.4 kcal/mol. Our study comprises the first measurement of a transmembrane helix interaction in a biological membrane, and we are optimistic that it can be further developed and applied.