Structure of the transmembrane dimer interface of glycophorin A in membrane bilayers

The hydrophobic transmembrane domain of glycophorin A contains a sequence motif that mediates dimerization in membrane environments. Long-range interhelical distance measurements using magic angle spinning NMR spectroscopy provide high-resolution structural constraints on the packing of the dimer interface in membrane bilayers. We show that direct packing contacts occur between glycine residues at positions 79 and 83 in the transmembrane sequence. Additional interhelical constraints between Ile76 and Gly79 and between Val80 and Gly83 restrict the rotational orientation and crossing angle of the interacting helices. These results refine our previously proposed structure of the glycophorin A dimer [Smith, S. O., and Bormann, B. J. (1995) Proc. Natl. Acad. Sci. U.S.A. 92, 488-491] which revealed that the methyl groups of Val80 and Val84 are packed against Gly79 and Gly83, respectively.     

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.     

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.     

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.     

Influence of proline residues in transmembrane helix packing

Integral membrane proteins often contain proline residues in their alpha-helical transmembrane (TM) fragments, which may strongly influence their folding and association. Pro-scanning mutagenesis of the helical domain of glycophorin A (GpA) showed that replacement of the residues located at the center abrogates helix packing while substitution of the residues forming the ending helical turns allows dimer formation. Synthetic TM peptides revealed that a point mutation of one of the residues of the dimerization motif (L75P) located at the N-terminal helical turn of the GpA TM fragment, adopts a secondary structure and oligomeric state similar to the wild-type sequence in detergents. In addition, both glycosylation mapping in biological membranes and molecular dynamics showed that the presence of a proline residue at the lipid/water interface has as an effect the extension of the helical end. Thus, helix packing can be an important factor that determines appearance of proline in TM helices. Membrane proteins might accumulate proline residues at the two ends of their TM segments in order to modulate the exposition of key amino acid residues at the interface for molecular recognition events while allowing stable association and native folding.     

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.     

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.     

Influence of hydrophobic matching on association of model transmembrane fragments containing a minimised glycophorin A dimerisation motif

The principles that govern the folding and packing of membrane proteins are still not completely understood. In the present work, we have revisited the glycophorin A (GpA) dimerisation motif that mediates transmembrane (TM) helix association, one of the best-suited models of membrane protein oligomerisation. By using artificial polyleucine TM segments we have demonstrated in this study that a pattern of only five amino acids (GVxxGVxxT) promotes specific dimerisation. Further, we have used this minimised GpA motif to assess the influence of hydrophobic matching on the TM helix packing process in detergent micelles and found that this factor modulates helix-helix association and/or dissociation between TM fragments.     

Structure of the membrane domain of subunit b of the Escherichia coli F0F1 ATP synthase.

The structure of the N-terminal transmembrane domain (residues 1-34) of subunit b of the Escherichia coli F0F1-ATP synthase has been solved by two-dimensional 1H NMR in a membrane mimetic solvent mixture of chloroform/methanol/H2O (4:4:1). Residues 4-22 form an alpha-helix, which is likely to span the hydrophobic domain of the lipid bilayer to anchor the largely hydrophilic subunit b in the membrane. The helical structure is interrupted by a rigid bend in the region of residues 23-26 with alpha-helical structure resuming at Pro-27 at an angle offset by 20 degrees from the transmembrane helix. In native subunit b, the hinge region and C-terminal alpha-helical segment would connect the transmembrane helix to the cytoplasmic domain. The transmembrane domains of the two subunit b in F0 were shown to be close to each other by cross-linking experiments in which single Cys were substituted for residues 2-21 of the native subunit and b-b dimer formation tested after oxidation with Cu(II)(phenanthroline)2. Cys residues that formed disulfide cross-links were found with a periodicity indicative of one face of an alpha-helix, over the span of residues 2-18, where Cys at positions 2, 6, and 10 formed dimers in highest yield. A model for the dimer is presented based upon the NMR structure and distance constraints from the cross-linking data. The transmembrane alpha-helices are positioned at a 23 degrees angle to each other with the side chains of Thr-6, Gln-10, Phe-14, and Phe-17 at the interface between subunits. The change in direction of helical packing at the hinge region may be important in the functional interaction of the cytoplasmic domains.     

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.     

Influence of the C-terminus of the glycophorin A transmembrane fragment on the dimerization process

The monomer-dimer equilibrium of the glycophorin A (GpA) transmembrane (TM) fragment has been used as a model system to investigate the amino acid sequence requirements that permit an appropriate helix-helix packing in a membrane-mimetic environment. In particular, we have focused on a region of the helix where no crucial residues for packing have been yet reported. Various deletion and replacement mutants in the C-terminal region of the TM fragment showed that the distance between the dimerization motif and the flanking charged residues from the cytoplasmic side of the protein is important for helix packing. Furthermore, selected GpA mutants have been used to illustrate the rearrangement of TM fragments that takes place when leucine repeats are introduced in such protein segments. We also show that secondary structure of GpA derivatives was independent from dimerization, in agreement with the two-stage model for membrane protein folding and oligomerization.     

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.     

A Ser residue influences the structure and stability of a Pro-kinked transmembrane helix dimer

When localized adjacent to a Pro-kink, Thr and Ser residues can form hydrogen bonds between their polar hydroxyl group and a backbone carbonyl oxygen and thereby modulate the actual bending angle of a distorted transmembrane α-helix. We have used the homo-dimeric transmembrane cytochrome b(559)' to analyze the potential role of a highly conserved Ser residue for assembly and stabilization of transmembrane proteins. Mutation of the conserved Ser residue to Ala resulted in altered heme binding properties and in increased stability of the holo-protein, most likely by tolerating subtle structural rearrangements upon heme binding. The results suggest a crucial impact of an intrahelical Ser hydrogen bond in defining the structure of a Pro-kinked transmembrane helix dimer.     

Structure of the 1-36 N-terminal fragment of human phospholamban phosphorylated at Ser-16 and Thr-17

The structure of a 36-amino-acid-long N-terminal fragment of human phospholamban phosphorylated at Ser-16 and Thr-17 and Cys-36-->Ser mutated was determined from nuclear magnetic resonance data in aqueous solution containing 30% trifluoroethanol. The peptide assumes a conformation characterized by two alpha-helices connected by an irregular strand, which comprises the amino acids from Arg-13 to Pro-21. The proline is in a trans conformation. The two phosphate groups on Ser-16 and Thr-17 are shown to interact preferably with the side chains of Arg-14 and Arg-13, respectively. The helix comprising amino acids 22 to 35 is well determined (the rmsd for the backbone atoms, calculated for a family of 24 nuclear magnetic resonance structures is 0.69 +/- 0.28 A). The structures of phosphorylated and unphosphorylated phospholamban are compared, and the effect of the two phosphate groups on the relative spatial position of the two helices is examined. The packing parameters Omega (interhelical angle) and d (minimal interhelical distance) are calculated: in the case of the phosphorylated phospholamban, Omega = 100 +/- 35 degrees and d = 7.9 +/- 4.6 A, whereas for the unphosphorylated peptide the values are Omega = 80 +/- 20 degrees and d = 7.0 +/- 4.0 A. We conclude that 1) the phosphorylation does not affect the structure of the C terminus between residues 21 and 36 and 2) the phosphorylated phospholamban has more loose helical packing than the nonphosphorylated.     

The role of a conserved proline residue in mediating conformational changes associated with voltage gating of Cx32 gap junctions.

We have explored the role of a proline residue located at position 87 in the second transmembrane segment (TM2) of gap junctions in the mechanism of voltage-dependent gating of connexin32 (Cx32). Substitution of this proline (denoted Cx32P87) with residues G, A, or V affects channel function in a progressive manner consistent with the expectation that a proline kink (PK) motif exists in the second transmembrane segment (TM2) of this connexin. Mutations of the preceding threonine residue T86 to S, A, C, V, N, or L shift the conductance-voltage relation of wild-type Cx32, such that the mutated channels close at smaller transjunctional voltages. The observed shift in voltage dependence is consistent with a reduction in the open probability of the mutant hemichannels at a transjunctional voltage (Vj) of 0 mV. In both cases in which kinetics were examined, the time constants for reaching steady state were faster for T86N and T86A than for wild type at comparable voltages, suggesting that the T86 mutations cause the energetic destabilization of the open state relative to the other states of the channel protein. The structural underpinnings of the observed effects were explored with Monte Carlo simulations. The conformational space of TM2 helices was found to differ for the T86A, V, N, and L mutants, which produce a less bent helix ( approximately 20 degrees bend angle) compared to the wild type, which has a approximately 37 degrees bend angle. The greater bend angle of the wild-type helix reflects the propensity of the T86 residue to hydrogen bond with the backbone carbonyl of amino acid residue I82. The relative differences in propensity for hydrogen bonding of the mutants relative to the wild-type threonine residue in the constructs we studied (T86A, V, N, L, S, and C) correlate with the shift in the conductance-voltage relation observed for T86 mutations. The data are consistent with a structural model in which the open conformation of the Cx32 channel corresponds to a more bent TM2 helix, and the closed conformation corresponds to a less bent helix. We propose that the modulation of the hydrogen-bonding potential of the T86 residue alters the bend angle of the PK motif and mediates conformational changes between open and closed channel states.     

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.     

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

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

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.