Implications of threonine hydrogen bonding in the glycophorin A transmembrane helix dimer

The transmembrane helix of glycophorin A contains a seven-residue motif, LIxxGVxxGVxxT, that mediates protein dimerization. Threonine is the only polar amino acid in this motif with the potential to stabilize the dimer through hydrogen-bonding interactions. Polarized Fourier transform infrared spectroscopy is used to establish a robust protocol for incorporating glycophorin A transmembrane peptides into membrane bilayers. Analysis of the dichroic ratio of the 1655-cm(-1) amide I vibration indicates that peptides reconstituted by detergent dialysis have a transmembrane orientation with a helix crossing angle of <35 degrees. Solid-state nuclear magnetic resonance spectroscopy is used to establish high resolution structural restraints on the conformation and packing of Thr-87 in the dimer interface. Rotational resonance measurement of a 2.9-A distance between the gamma-methyl and backbone carbonyl carbons of Thr-87 is consistent with a gauche- conformation for the chi1 torsion angle. Rotational-echo double-resonance measurements demonstrate close packing (4.0 +/- 0.2 A) of the Thr-87 gamma-methyl group with the backbone nitrogen of Ile-88 across the dimer interface. The short interhelical distance places the beta-hydroxyl of Thr-87 within hydrogen-bonding range of the backbone carbonyl of Val-84 on the opposing helix. These results refine the structure of the glycophorin A dimer in membrane bilayers and highlight the complementary role of small and polar residues in the tight association of transmembrane helices in membrane proteins.     

Role of Side-Chain Conformational Entropy in Transmembrane Helix Dimerization of Glycophorin A

Dimerization of the transmembrane domain of glycophorin A is mediated by a seven residue motif LIxxGVxxGVxxT through a combination of van der Waals and hydrogen bonding interactions. One of the unusual features of the motif is the large number of β-branched amino acids that may limit the entropic cost of dimerization by restricting side-chain motion in the monomeric transmembrane helix. Deuterium NMR spectroscopy is used to characterize the dynamics of fully deuterated Val80 and Val84, two essential amino acids of the dimerization motif. Deuterium spectra of the glycophorin A transmembrane dimer were obtained using synthetic peptides corresponding to the transmembrane sequence containing either perdeuterated Val80 or Val84. These data were compared with spectra of monomeric glycophorin A peptides deuterated at Val84. In all cases, the deuterium line shapes are characterized by fast methyl group rotation with virtually no motion about the C_α-C_β bond. This is consistent with restriction of the side chain in both the monomer and dimer due to intrahelical packing interactions involving the β-methyl groups, and indicates that there is no energy cost associated with dimerization due to loss of conformational entropy. In contrast, deuterium NMR spectra of Met81 and Val82, in the lipid interface, reflected greater motional averaging and fast exchange between different side-chain conformers.     

Transmembrane interactions in the activation of the Neu receptor tyrosine kinase

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

Computation and mutagenesis suggest a right-handed structure for the synaptobrevin transmembrane dimer.

Biological membrane fusion involves a highly precise and ordered set of protein-protein interactions. Synaptobrevin is a key player in this process. Mutagenesis studies of its single transmembrane segment suggest that it dimerizes in a sequence specific manner. Using the computational methods developed for the successful structure prediction of the glycophorin A transmembrane dimer, we have calculated a structural model for the synaptobrevin dimer. Our computational search yields a well-populated cluster of right-handed structures consistent with the experimentally determined dimerization motif. The three-dimensional structure contains an interface formed primarily by leucine and isoleucine side-chain atoms and has no interhelical hydrogen bonds. The model is the first three-dimensional picture of the synaptobrevin transmembrane dimer and provides a basis for further focused experimentation on its structure and association thermodynamics.     

Sequence context modulates the stability of a GxxxG-mediated transmembrane helix-helix dimer

To quantify the relationship between sequence and transmembrane dimer stability, a systematic mutagenesis and thermodynamic study of the protein-protein interaction residues in the glycophorin A transmembrane helix-helix dimer was carried out. The results demonstrate that the glycophorin A transmembrane sequence dimerizes when its GxxxG motif is abolished by mutation to large aliphatic residues, suggesting that the sequence encodes an intrinsic propensity to self-associate independent of a GxxxG motif. In the presence of an intact GxxxG motif, the glycophorin A dimer stability can be modulated over a span of -0.5 kcal mol(-1) to +3.2 kcal mol(-1) by mutating the surrounding sequence context. Thus, these flanking residues play an active role in determining the transmembrane dimer stability. To assess the structural consequences of the thermodynamic effects of mutations, molecular models of mutant transmembrane domains were constructed, and a structure-based parameterization of the free energy change due to mutation was carried out. The changes in association free energy for glycophorin A mutants can be explained primarily by changes in packing interactions at the protein-protein interface. The energy cost of removing favorable van der Waals interactions was found to be 0.039 kcal mol(-1) per A2 of favorable occluded surface area. The value corresponds well with estimates for mutations in bacteriorhodopsin as well as for those mutations in the interiors of soluble proteins that create packing defects.     

The GxxxG motif: a framework for transmembrane helix-helix association

In order to identify strong transmembrane helix packing motifs, we have selected transmembrane domains exhibiting high-affinity homo-oligomerization from a randomized sequence library based on the right-handed dimerization motif of glycophorin A. Sequences were isolated using the TOXCAT system, which measures transmembrane helix-helix association in the Escherichia coli inner membrane. Strong selection was applied to a large range of sequences ( approximately 10(7) possibilities) and resulted in the identification of sequence patterns that mediate high-affinity helix-helix association. The most frequent motif isolated, GxxxG, occurs in over 80% of the isolates. Additional correlations suggest that flanking residues act in concert with the GxxxG motif, and that size complementarity is maintained at the interface, consistent with the idea that the identified sequence patterns represent packing motifs. The convergent identification of similar sequence patterns from an analysis of the transmembrane domains in the SwissProt sequence database suggests that these packing motifs are frequently utilized in naturally occurring helical 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.     

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.     

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.     

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.     

Transmembrane helix-helix interactions: comparative simulations of the glycophorin a dimer

The glycophorin helix dimer is a paradigm for the exploration of helix-helix interactions in integral membrane proteins. Two NMR structures of the dimer are known, one in a detergent micelle and one in a lipid bilayer. Multiple (4 x 50 ns) molecular dynamics simulations starting from each of the two NMR structures, with each structure in either a dodecyl phosphocholine (DPC) micelle or a dimyristoyl phosphatidylcholine (DMPC) bilayer, have been used to explore the conformational dynamics of the helix dimer. Analysis of the helix-helix interaction, mediated by the GxxxG sequence motif, suggests convergence of the simulations to a common model. This is closer to the NMR structure determined in a bilayer than to micelle structure. The stable dimer interface in the final simulation model is characterized by (i) Gly/Gly packing and (ii) Thr/Thr interhelix H-bonds. These results demonstrate the ability of extended molecular dynamics simulations in a lipid bilayer environment to refine membrane protein structures or models derived from experimental data obtained in protein/detergent micelles.     

Lipid solvation effects contribute to the affinity of Gly-xxx-Gly motif-mediated helix-helix interactions

Interactions between transmembrane helices are mediated by the concave Gly-xxx-Gly motif surface. Whether Gly residues per se are sufficient for selection of this motif has not been established. Here, we used the in vivo TOXCAT assay to measure the relative affinities of all 18 combinations of Gly, Ala, and Ser "small-xxx-small" mutations in glycophorin A (GpA) and bacteriophage M13 major coat protein (MCP) homodimers. Affinity values were compared with the accessibility to a methylene-sized probe of the total surface area of each helix monomer as a measure of solvation by membrane components. A strong inverse correlation was found between nonpolar-group lipid accessibility and dimer affinity (R = 0.75 for GpA, p = 0.013, and R = 0.81 for MCP, p = 0.004), suggesting that lipid as a poor membrane protein solvent, conceptually analogous to water in soluble protein folding, can contribute to dimer stability and help to define helix-helix interfaces.     

Association of a model transmembrane peptide containing gly in a heptad sequence motif

A peptide containing glycine at a and d positions of a heptad motif was synthesized to investigate the possibility that membrane-soluble peptides with a Gly-based, left-handed helical packing motif would associate. Based on analytical ultracentrifugation in C14-betaine detergent micelles, the peptide did associate in a monomer-dimer equilibrium, although the association constant was significantly less than that reported for the right-handed dimer of the glycophorin A transmembrane peptide in similar detergents. Fluorescence resonance energy transfer (FRET) experiments conducted on peptides labeled at their N-termini with either tetramethylrhodamine (TMR) or 7-nitrobenz-2-oxa-1,3-diazole (NBD) also indicated association. However, analysis of the FRET data using the usual assumption of complete quenching for NBD-TMR pairs in the dimer could not be quantitatively reconciled with the analytical ultracentrifugation-measured dimerization constant. This led us to develop a general treatment for the association of helices to either parallel or antiparallel structures of any aggregation state. Applying this treatment to the FRET data, constraining the dimerization constant to be within experimental uncertainty of that measured by analytical ultracentrifugation, we found the data could be well described by a monomer-dimer equilibrium with only partial quenching of the dimer, suggesting that the helices are most probably antiparallel. These results also suggest that a left-handed Gly heptad repeat motif can drive membrane helix association, but the affinity is likely to be less strong than the previously reported right-handed motif described for glycophorin A.     

Statistical analysis of amino acid patterns in transmembrane helices: the GxxxG motif occurs frequently and in association with beta-branched residues at neighboring positions

To find motifs that mediate helix-helix interactions in membrane proteins, we have analyzed frequently occurring combinations of residues in a database of transmembrane domains. Our analysis was performed with a novel formalism, which we call TMSTAT, for exactly calculating the expectancies of all pairs and triplets of residues in individual sequences, taking into account differential sequence composition and the substantial effect of finite length in short segments. We found that the number of significantly over and under-represented pairs and triplets was much greater than the random expectation. Isoleucine, glycine and valine were the most common residues in these extreme cases. The main theme observed is patterns of small residues (Gly, Ala and Ser) at i and i+4 found in association with large aliphatic residues (Ile, Val and Leu) at neighboring positions (i.e. i+/-1 and i+/-2). The most over-represented pair is formed by two glycine residues at i and i+4 (GxxxG, 31.6 % above expectation, p<1x10(-33)) and it is strongly associated with the neighboring beta-branched residues Ile and Val. In fact, the GxxxG pair has been described as part of the strong interaction motif in the glycophorin A transmembrane dimer, in which the pair is associated with two Val residues (GVxxGV). GxxxG is also the major motif identified using TOXCAT, an in vivo selection system for transmembrane oligomerization motifs. In conjunction with these experimental observations, our results highlight the importance of the GxxxG+beta-branched motif in transmembrane helix-helix interactions. In addition, the special role for the beta-branched residues Ile and Val suggested here is consistent with the hypothesis that residues with constrained rotameric freedom in helical conformation might reduce the entropic cost of folding in transmembrane proteins. Additional material is available at http://engelman.csb.yale. edu/tmstat and http://bioinfo.mbb.yale. edu/tmstat.     

Changes in apparent free energy of helix-helix dimerization in a biological membrane due to point mutations

We present an implementation of the TOXCAT membrane protein self-association assay that measures the change in apparent free energy of transmembrane helix dimerization caused by point mutations. Quantifying the reporter gene expression from cells carrying wild-type and mutant constructs shows that single point mutations that disrupt dimerization of the transmembrane domain of glycophorin A reproducibly lower the TOXCAT signal more than 100-fold. Replicate cultures can show up to threefold changes in the level of expression of the membrane bound fusion construct, and correcting for these variations improves the precision of the calculated apparent free energy change. The remarkably good agreement between our TOXCAT apparent free energy scale and free energy differences from sedimentation equilibrium studies for point mutants of the glycophorin A transmembrane domain dimer indicate that sequence changes usually affect membrane helix-helix interactions quite similarly in these two very different environments. However, the effects of point mutations at threonine 87 suggest that intermonomer polar contacts by this side-chain contribute significantly to dimer stability in membranes but not in detergents. Our findings demonstrate that a comparison of quantitative measurements of helix-helix interactions in biological membranes and genuine thermodynamic data from biophysical measurements on purified proteins can elucidate how changes in the lipidic environment modulate membrane protein stability.     

Amino acid composition of parallel helix-helix interfaces

Amino acids at helix-helix parallel interfaces influence arrangement of helices and interhelical angles. Parallel interfaces in 79 proteins were considered. Location of amino acids at the positions analogous to a and d in GCN4 leucine zipper nomenclature shows that certain combinations of amino acids characteristic for parallel packing occur more often than could be expected by chance. Repeating sequence combinations occur at a and d positions of parallel helix-helix interfaces with similar values of interhelical angles not only in homologous proteins but also within the same protein and in nonhomologous proteins. Within each group of observed combinations correlation exists between the size of amino acid and magnitude of the interhelical angle.     

Structural role of glycine in amyloid fibrils formed from transmembrane alpha-helices

Amyloid fibrils associated with diseases such as Alzheimer's are often derived from the transmembrane helices of membrane proteins. It is known that the fibrils have a cross-beta-sheet structure where main chain hydrogen bonding occurs between beta-strands in the direction of the fibril axis. However, the structural basis for how the membrane-spanning helix is converted into a beta-sheet or how protofibrils associate into fibrils is not known. Here, we use a model peptide corresponding to a portion of the single transmembrane helix of glycophorin A to investigate the structural role of glycine in amyloid-like fibrils formed from transmembrane helices. Glycophorin A contains a GxxxG motif that is found in many transmembrane sequences including that of the amyloid precursor protein and prion protein. We propose that glycine, which mediates helix interactions in membrane proteins, also provides key packing motifs when it occurs in beta-sheets. We show that glycines in the glycophorin A transmembrane helix promote extended beta-strand formation when the helix partitions into aqueous environments and stabilize the packing of beta-sheets in the formation of amyloid-like fibrils. We demonstrate that fibrillization can be disrupted with a new class of inhibitors that target the molecular grooves created by glycine.     

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.     

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.     

Helix-helix packing and interfacial pairwise interactions of residues in membrane proteins

Helix-helix packing plays a critical role in maintaining the tertiary structures of helical membrane proteins. By examining the overall distribution of voids and pockets in the transmembrane (TM) regions of helical membrane proteins, we found that bacteriorhodopsin and halorhodopsin are the most tightly packed, whereas mechanosensitive channel is the least tightly packed. Large residues F, W, and H have the highest propensity to be in a TM void or a pocket, whereas small residues such as S, G, A, and T are least likely to be found in a void or a pocket. The coordination number for non-bonded interactions for each of the residue types is found to correlate with the size of the residue. To assess specific interhelical interactions between residues, we have developed a new computational method to characterize nearest neighboring atoms that are in physical contact. Using an atom-based probabilistic model, we estimate the membrane helical interfacial pairwise (MHIP) propensity. We found that there are many residue pairs that have high propensity for interhelical interactions, but disulfide bonds are rarely found in the TM regions. The high propensity pairs include residue pairs between an aromatic residue and a basic residue (W-R, W-H, and Y-K). In addition, many residue pairs have high propensity to form interhelical polar-polar atomic contacts, for example, residue pairs between two ionizable residues, between one ionizable residue and one N or Q. Soluble proteins do not share this pattern of diverse polar-polar interhelical interaction. Exploratory analysis by clustering of the MHIP values suggests that residues similar in side-chain branchness, cyclic structures, and size tend to have correlated behavior in participating interhelical interactions. A chi-square test rejects the null hypothesis that membrane protein and soluble protein have the same distribution of interhelical pairwise propensity. This observation may help us to understand the folding mechanism of membrane proteins.