Search

For many membrane proteins, self-association serves both structural and functional roles. Studies of such association can be simplified by switching to micelles as the membrane-mimicking environment, but native interaction is not preserved in all detergents. The selection of suitable conditions for biochemical experiments would be greatly facilitated by a quantitative high-throughput assay. Here we showed that the fluorescence polarization reduction, which resulted from homo-Förster resonance energy transfer and was measured in a high-throughput compatible format, can be used to determine both association states and constants for membrane proteins in micelles.

Cytochrome b(559)' is a transmembrane protein formed by homodimerization of the 44-residue PsbF polypeptide and noncovalent binding of a heme cofactor. The PsbF polypeptide can dimerize in the absence and presence of heme. To monitor structural alterations associated with binding of heme to the apo-cytochrome, we analyzed the apo- and holo-cytochrome structure by electron paramagnetic resonance spectroscopy. Spin labeling of amino acids located close to the heme binding domain of the cytochrome revealed that the structure of the heme binding domain is unconstrained in the absence of heme. Heme binding restricts the conformational dynamics of the heme binding domain, resulting in the structurally more constricted holo-cytochrome structure.

Oligomerization of transmembrane (TM) helices is a key stage in the folding of membrane proteins. Glycophorin A (GpA) is a well-documented test system for this process. Coarse-grained molecular dynamics (CG-MD) allows us to simulate the self-assembly of TM helices into dimers, for both wild-type (WT) and mutant GpA sequences. For the WT sequences, dimers formed rapidly and remained stable in all simulations. The resultant dimers exhibited right-handed crossing and the same interhelix contacts as in NMR structures. Simulations of disruptive mutants revealed the dimers were less stable, with values of DeltaDelta G dimerization consistent with experimental data. The dimers of disruptive mutants were distorted relative to the WT and showed left-handed crossing of their helices. CG-MD can therefore be used to explore the interactions of TM helices, an important stage in the folding of membrane proteins. In particular, CG-MD has been shown to be sensitive enough to detect disruptions introduced by mutation. Future refinement of such models via atomistic simulations will enable a multiscale approach to predict the folding of membrane proteins.

The cation-pi interaction is an electrostatic attraction between a positive charge and the conjugated pi electrons of an aromatic ring. These interactions are well documented in soluble proteins and can be both structurally and functionally important. Catalyzed by observations in our laboratory that an Ala- and Ile-rich two-helix transmembrane segment tended to form SDS-resistant dimers upon the incorporation of suitably located Trp residues, here we have constructed a library of related constructs to study systematically the impact of aromatic-aromatic and cation-pi interactions on tertiary structure formation within an Escherichia coli membrane. Using the TOXCAT oligomerization assay with the hydrophobic segment AIAIAIIAZAXAIIAIAIAI, where Z = A, W, Y, or F and X = A, H, R, or K in all possible combinations of cation and/or aromatic pairings, to assess the TM-TM dependent expression of the chloramphenicol acetyltransferase reporter gene, we find that cation-pi interactions, particularly between Lys and Trp, Tyr, or Phe, as well as weakly polar interactions between pairs of aromatic residues, significantly enhance the strength of oligomerization of these hydrophobic helices, in some instances forming oligomers four times stronger than the high-affinity glycophorin A dimer. The contribution of these forces to the tertiary structure formation in designed transmembrane segments suggests that similar forces may also be a significant factor in the folding and stability of native membrane proteins.

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.

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.

Structural analyses of membrane proteins reveal a large number of C(alpha)-H...O contacts between transmembrane helices, presumed to be hydrogen bonds. Recent experiments produced conflicting results for the contribution of such hydrogen bonds to membrane protein stability. An FTIR study estimated an energy of -0.88 kcal/mol for the G79-C(alpha)-H...I76-O hydrogen bond in glycophorin A, whereas a mutagenesis study showed that the A51-C(alpha)-H...T24-O(gamma) hydrogen bond does not stabilize bacteriorhodopsin. Here, we reconcile these results using molecular mechanics calculations and an implicit membrane model (IMM1). With explicit hydrogen atoms, the potential energy of the G79-C(alpha)-H...I76-O interaction in GpA ranges from -0.54 to -0.9 kcal/mol and its contribution to stability (effective energy) from -0.49 to -0.83 kcal/mol, depending on the structural model used. The average values of these quantities in GpA-like motifs are similar. In bR, the corresponding numbers for the A51-C(alpha)-H...T24-O(gamma) interaction are +0.15 and +0.32 kcal/mol. The difference results from the different arrangement of the interacting groups and specifically the position of the acceptor with respect to the C(alpha) and N atoms. This conclusion likely applies to soluble proteins as well.

To clarify the function of the hydrophilic carboxyl-terminal tail of human erythrocyte membrane band 3 protein (HEM-B3), we purified two peptides, C1 (Ala893-Val911) and KS4 (Gly647-Arg656), from human erythrocyte band 3 protein preparations. Purified C1 peptides at concentrations from 5 to 80 microM were incubated with fresh human erythrocyte white ghosts. The C1 peptide demonstrated a novel protease activity, which cleaved glycophorin A (GPA) at Leu118-Ser119 in a dose-dependent manner. This activity was eliminated by trypsin. In a control experiment, the KS4 peptide did not cleave GPA under the same conditions. To help substantiate that the band 3 C-terminal tail peptide (C1) alone possesses the protease activity, two experiments were performed. First, the plasmids pGBKT(7)-GPA-Ct and pGADT(7)-AE1-Ct were cotransformed into the yeast strain AH109. The pGBKT(7)-GPA-Ct plasmid contains the cDNA of the 33 amino acid residue section of GPA (Tyr93-Asn125) fused with the pGBKT(7) vector. The plasmid pGADT(7)-AE1-Ct contains the cDNA of the C-terminal 33 amino acid residues of HEM-B3 fused with the GAL4 DNA-binding domain in the pGADT(7) vector. The results of the cotransformation experiment indicated that the C-terminal 33 amino acid residues of HEM-B3 interacted directly with the GPA C-terminal segment defined above. Second, we used a mammalian two-hybrid analysis to confirm the interaction relationship between the band 3 C-terminal segment and the GPA C-terminus. The C-terminus of GPA and the C-terminal 33 amino acid residues of HEM-B3 were subcloned into the DNA-binding domain and transcription activation domain vectors of the two-hybrid system, respectively. They were then cotransfected along with a chloramphenicol acetyltransferse (CAT) reporter vector into HeLa cells. The CAT activity measured in this experiment also indicated that there was interaction between the C-terminal 33 amino acid residues of HEM-B3 and the C-terminus of GPA.

Synaptobrevin 2 is thought to facilitate fusion of synaptic vesicles with the presynaptic membrane through formation of a soluble NSF attachment protein receptor complex (SNARE) with syntaxin 1a and a synaptosomal associated protein of 25 kDa (SNAP-25). Previous reports have described a homodimer of synaptobrevin that is dependent on the transmembrane domain. However, these reports disagree about the magnitude of dimerization, which makes it difficult to assess the biological relevance of this interaction. We used SDS-PAGE and the TOXCAT genetic assay to reexamine the homodimerization of the synaptobrevin transmembrane domain in detergents and the Escherichia coli inner membrane, respectively. To gauge the magnitude of synaptobrevin homodimerization, we used the well-characterized glycophorin A homodimer as a positive standard. In contrast to previous studies, we found synaptobrevin homodimerization in E. coli is very weak when compared to glycophorin A. Recombinant synaptobrevin forms a small amount of dimer and higher order oligomers in detergents that are highly dependent on solublization conditions. We estimate a dissociation constant of 10 mM for synaptobrevin dimerization in detergent. Thus, the dimerization of synaptobrevin in membranes is very weak, questioning any possible functional role for this association in vivo.

The GXXXG motif is a frequently occurring sequence of residues that is known to favor helix-helix interactions in membrane proteins. Here we show that the GXXXG motif is also prevalent in soluble proteins whose structures have been determined. Some 152 proteins from a non-redundant PDB set contain at least one alpha-helix with the GXXXG motif, 41 +/- 9% more than expected if glycine residues were uniformly distributed in those alpha-helices. More than 50% of the GXXXG-containing alpha-helices participate in helix-helix interactions. In fact, 26 of those helix-helix interactions are structurally similar to the helix-helix interaction of the glycophorin A dimer, where two transmembrane helices associate to form a dimer stabilized by the GXXXG motif. As for the glycophorin A structure, we find backbone-to-backbone atomic contacts of the C alpha-H...O type in each of these 26 helix-helix interactions that display the stereochemical hallmarks of hydrogen bond formation. These glycophorin A-like helix-helix interactions are enriched in the general set of helix-helix interactions containing the GXXXG motif, suggesting that the inferred C alpha-H...O hydrogen bonds stabilize the helix-helix interactions. In addition to the GXXXG motif, some 808 proteins from the non-redundant PDB set contain at least one alpha-helix with the AXXXA motif (30 +/- 3% greater than expected). Both the GXXXG and AXXXA motifs occur frequently in predicted alpha-helices from 24 fully sequenced genomes. Occurrence of the AXXXA motif is enhanced to a greater extent in thermophiles than in mesophiles, suggesting that helical interaction based on the AXXXA motif may be a common mechanism of thermostability in protein structures. We conclude that the GXXXG sequence motif stabilizes helix-helix interactions in proteins, and that the AXXXA sequence motif also stabilizes the folded state of proteins.

As part of a structural genomics project, we have determined the 2.0 A structure of the E1beta subunit of pyruvate dehydrogenase from Pyrobaculum aerophilum (PA), a thermophilic archaeon. The overall fold of E1beta from PA is closely similar to the previously determined E1beta structures from humans (HU) and P. putida (PP). However, unlike the HU and PP structures, the PA structure was determined in the absence of its partner subunit, E1alpha. Significant structural rearrangements occur in E1beta when its E1alpha partner is absent, including rearrangement of several secondary structure elements such as helix C. Helix C is buried by E1alpha in the HU and PP structures, but makes crystal contacts in the PA structure that lead to an apparent beta(4) tetramer. Static light scattering and sedimentation velocity data are consistent with the formation of PA E1beta tetramers in solution. The interaction of helix C with its symmetry-related counterpart stabilizes the tetrameric interface, where two glycine residues on the same face of one helix create a packing surface for the other helix. This GPhiXXG helix-helix interaction motif has previously been found in interacting transmembrane helices, and is found here at the E1alpha-E1beta interface for both the HU and PP alpha(2)beta(2) tetramers. As a case study in structural genomics, this work illustrates that comparative analysis of protein structures can identify the structural significance of a sequence motif.

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.

Glycophorin A (GPA) has been reconstituted into dimyristoylphosphatidylcholine vesicles and digested with proteinase K to identify the membrane domain and to characterize its structure and orientation. After digestion of the inner and outer domain of GPA by protease action restricted to the aqueous phase, a protected peptide migrates on an electrophoresis gel as a 7.5-kDa dimer (His66-Ile95). The secondary structure and orientation in a lipid bilayer of the 7.5-kDa dimer have been studied by Fourier transform infrared spectroscopy. Our proteolytic and spectroscopic data are in agreement with a topological model in which the His66-Glu72 peptide adopts a beta-sheet conformation and is oriented parallel to the lipid-water interface and the Ile73-Ile95 domain is helical and oriented parallel to the lipid acyl chains, in a transmembrane configuration. Digestion of the domain protruding to the outside of the liposome generates "head-head" and "head-tail" dimers of 16 and 38 kDa, respectively. This observation is discussed in terms of the specificity of the dimer formation process.

Rotational resonance (RR) NMR, circular dichroism (CD), and attenuated total reflection Fourier transform infrared (ATR-FTIR) spectroscopy are used to establish the secondary structure and orientation of peptides corresponding to the transmembrane domain of human glycophorin A in dimyristoylphosphatidylcholine bilayers. An amide I vibrational frequency of 1650 cm-1 and negative CD absorption bands at 208 and 222 nm indicate that the peptide is largely alpha-helical, while an order parameter of 0.35-0.50 in the ATR-FTIR measurements indicates that the peptide orientation is generally perpendicular to the bilayer plane. High-resolution structural data on the glycophorin A transmembrane (GPA-TM) peptides were obtained by measuring the rate of magnetization exchange between pairs of specific 13C labels using RR NMR. The exchange rates are translated into internuclear distances with a resolution on the order of 0.3 A. These experiments are similar in design to previous experiments on crystalline peptides where the 13C labels were incorporated into amino acids separated by 2-3 residues in the peptide sequence but close together in space due to a helical peptide geometry [Peersen, O.B., Yoshimura, S., Hojo, H., Aimoto, S., & Smith, S.O. (1992) J. Am. Chem. Soc. 114, 4332-4335]. In the GPA-TM peptides, magnetization exchange rates measured between [1-13C]V80 and [2-13C] G83 between [1-13C]M81 and [2-13C]G83 in the middle of the transmembrane sequence correspond to internuclear distances of approximately 4.5 A and are consistent with a helical peptide structure.(ABSTRACT TRUNCATED AT 250 WORDS)

Rotational resonance (RR) NMR methods are explored for determining intramolecular and intermolecular distances between 13C-labeled sites in membrane bilayers. Specific 13C labels have been incorporated into dipalmitoylphosphatidylcholine (DPPC) and a hydrophobic peptide corresponding to the transmembrane domain of human glycophorin A. The exchange of magnetization between these labels in the RR NMR experiment can be related to their internuclear separation. Intramolecular magnetization exchange rates have been measured between 13C labels incorporated at the 1-position of the sn-1 acyl chain (1-[1-13C]) and the 2-position of the sn-2 acyl chain (2-[2-13C]) of DPPC. These positions are 5.3-5.5 A apart in the crystal structure of dimyristoylphosphatidylcholine (DMPC), but are estimated to be 4.5-5.0 A apart in DPPC below the subphase transition temperature. These results are consistent with a smaller axial displacement between the sn-1 and sn-2 acyl chains than the approximately 3.6-A displacement observed in the DMPC crystal structure. Intermolecular magnetization transfer rates have been measured between 1,2-[2-13C]DPPC and 2-[1,13C]DPPC and between 1,2-[1-13C]DPPC and 2-[2-13C]DPPC. In addition, intermolecular magnetization exchange rates have been measured between 1,2-[2-13C]DPPC and the 13C-OH position of Y93 in the glycophorin transmembrane domain. These intermolecular distance measurements demonstrate that the relative orientation and location of membrane lipids and peptides can be established using RR NMR methods.

Glycophorin A and its isolated transmembrane region (GpATM) are each known to form sequence-specific dimers in SDS micelles. Whether this behavior accurately reflects behavior in red cell membranes or lipid bilayers, however, has remained unclear. Resonance energy transfer between labeled GpATM peptides has been used to observe dimerization of GpATM in bilayers. Separate populations of GpATM peptides were labeled with 2,6-dansyl chloride as the donor chromophore and dabsyl chloride as the acceptor. Quenching of the 2,6-dansyl chloride by the dabsyl group demonstrated an association of the labeled peptides. The quenching was not affected by increases in the amount of lipid present or by unlabeled heterologous peptides but was decreased by the addition of unlabeled GpATM. GpATM was determined to form dimers by fitting the observed energy transfer for a number of donor to acceptor ratios and fitting to the expected number of donor labeled peptides in an oligomer with an acceptor as a function of oligomer number. The finding that the GpATM peptide forms helical dimers in lipid bilayers supports the idea that GpA is a dimer in the erythrocyte membrane. The resonance energy transfer approach may extend to the study of other oligomeric complexes.

High-sensitivity differential scanning calorimetry (DSC) was used to study the mutual interactions between a glycoprotein (human glycophorin, GPA) and a glycosphingolipid (GM1 ganglioside) embedded in large unilamellar vesicles composed of dimyristoylphosphatidylcholine (DMPC). The DSC thermograms exhibited by DMPC/GM1 vesicles, either in the presence or in the absence of GPA, are resolvable into two components. The relative contribution of the minor component, centered at higher temperature, to the total enthalpy and its transition temperature increase with the concentration of the glycolipid embedded in the vesicles. This minor peak, undetectable in the absence of ganglioside, is indicative of the occurrence of lateral phase separation and suggests that GM1 ganglioside-enriched domains are present within the bilayer. At a given concentration of GM1 embedded in the vesicles, the proportion of the phase-separated peak is higher in the presence of GPA, suggesting that the glycoprotein enhances the tendency of GM1 to segregate. Experiments investigating the thermotropic behavior of GPA show that the temperature of irreversible thermal unfolding of the glycoprotein inserted in DMPC vesicles, centered at 65.9 degrees C in the absence of GM1, is shifted to 57.6 degrees C when GM1 is present in the bilayer. These results indicate that, at least in this experimental system, on the one hand, GPA enhances the tendency of the glycolipid to segregate within the membrane, and on the other hand, the glycolipid clusters affect the protein conformation and oligomerization in the membrane.

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.

Recent studies suggest specific roles for transmembrane helix association in a range of functions, but understanding of the conformation and energetics of these interactions has been elusive. We have studied the specific dimerization of the transmembrane helix of glycophorin A by calculating the minimized interaction energies of a large number of conformations using simulated annealing techniques and tested the models against mutational analysis data. We find that the dimer is best modeled as a right-handed supercoil with an extensive region of close packing along the dimer interface. Furthermore, we observe a sequence-specific propensity for a right-handed supercoil to form when starting the simulated annealing modeling from a dimer of helices with parallel axes, in contrast with the dimerization region of the transcription factor GCN4 which shows a high propensity for the more prevalent left-handed supercoiling.

Glycophorin has been incorporated into unilamellar cholesterol-containing dimyristoylphosphatidylcholine vesicles that were reconstituted by the freeze and thaw technique. Evidence was obtained for a protein-induced structural reorganization of these mixed membranes. By differential scanning calorimetry, we were able to construct a phase diagram for the phospholipid/cholesterol mixture consisting of a liquid-ordered, a solid-ordered, and a liquid-disordered phase. Glycophorin at low molar fractions (XG less than 3 X 10(-3)) increases the relative amount of lipid in the liquid-ordered phase, which is interpreted as an enrichment of cholesterol in the vicinity of the protein. The formation of such steroid-enriched domains could be demonstrated directly by electron paramagnetic resonance using a spin-labeled cholesterol analogue. A drastic increase of the spin-spin interaction of the labeled steroid was observed in the presence of glycophorin.

The asparagine-linked sugar chain of glucose transporter from human erythrocytes was quantitatively released as oligosaccharides from the polypeptide backbone by hydrazinolysis. They were converted to radioactive oligosaccharides by NaB3H4 reduction after N-acetylation and fractionated by anion-exchange column chromatography and Bio-Gel P-4 column chromatography after sialidase treatment. Structural study of each oligosaccharide by exo- and endoglycosidase digestion and methylation analysis indicated that the glycoprotein contains a high-mannose-type oligosaccharide, Man9.GlcNAc.GlcNAc, and biantennary complex-type oligosaccharides with Man alpha 1----6(+/- GlcNAc beta 1----4)(Man alpha 1----3) Man beta beta 1----4GlcNAc beta 1----4(+/- Fuc alpha 1----6)GlcNAc as their cores and the poly-N-acetyllactosamine composed of about 16 N-acetyllactosaminyl units as their outer chains. These structural features of the sugar moiety of glucose transporter are quite different from those of two major intrinsic glycoproteins of human erythrocytes, glycophorin A and band 3.