Energetics of the native and non-native states of the glycophorin transmembrane helix dimer

Using an implicit membrane model (IMM1), we examine whether the structure of the transmembrane domain of Glycophorin A (GpA) could be predicted based on energetic considerations alone. The energetics of native GpA shows that van der Waals interactions make the largest contribution to stability. Although specific electrostatic interactions are stabilizing, the overall electrostatic contribution is close to zero. The GXXXG motif contributes significantly to stability, but residues outside this motif contribute almost twice as much. To generate non-native states a global conformational search was done on two segments of GpA: an 18-residue peptide (GpA74-91) that is embedded in the membrane and a 29-residue peptide (GpA70-98) that has additional polar residues flanking the transmembrane region. Simulated annealing was done on a large number of conformations generated from parallel, antiparallel, left- and right-handed starting structures by rotating each helix at 20 degrees intervals around its helical axis. Several crossing points along the helix dimer were considered. For 18-residue parallel topology, an ensemble of native-like structures was found at the lowest effective energy region; the effective energy is lowest for a right-handed structure with an RMSD of 1.0 A from the solid-state NMR structure with correct orientation of the helices. For the 29-residue peptide, the effective energies of several left-handed structures were lower than that of the native, right-handed structure. This could be due to deficiencies in modeling the interactions between charged sidechains and/or omission of the sidechain entropy contribution to the free energy. For 18-residue antiparallel topology, both IMM1 and a Generalized Born model give effective energies that are lower than that of the native structure. In contrast, the Poisson-Boltzmann solvation model gives lower effective energy for the parallel topology, largely because the electrostatic solvation energy is more favorable for the parallel structure. IMM1 seems to underestimate the solvation free energy advantage when the CO and NH dipoles just outside the membrane are parallel. This highlights the importance of electrostatic interactions even when these are not obvious by looking at the structures.     

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.     

Complex interactions at the helix-helix interface stabilize the glycophorin A transmembrane dimer

To explore the residue interactions in the glycophorin A dimerization motif, an alanine scan double mutant analysis at the helix-helix interface was carried out. These data reveal a combination of additive and coupled effects. The majority of the double mutants are found to be equally or slightly more stable than would be predicted by the sum of the energetic cost of the single-point mutants. The proximity of the mutated sites is not related to the presence of coupling between those sites. Previous studies reveal that a single face of the glycophorin A monomer contains a specific glycine-containing motif (GxxxG) that is thought to be a driving force for the association of transmembrane helices. Double mutant cycles suggest that the relationship of the GxxxG motif to the remainder of the helix-helix interface is complex. Sequences containing mutations that abolish the GxxxG motif retain an ability to dimerize, while a sequence containing a GxxxG motif appears unable to form dimers. The energetic effects of weakly coupled and additive double mutants can be explained by changes in van der Waals interactions at the dimer interface. These results emphasize the fact that the sequence context of the dimer interface modulates the strength of the glycophorin A GxxxG-mediated transmembrane dimerization reaction.     

Thermodynamics of glycophorin A transmembrane helix dimerization in C14 betaine micelles

We have used sedimentation equilibrium analytical ultracentrifugation to measure the free energy change for the glycophorin A transmembrane helix-helix dimerization in C14 betaine micelles. By varying the amount of micellar C14 betaine, we show that the protein association reaction in the micellar C14 phase behaves as an ideal-dilute solution. In this hydrophobic environment, the mole-fraction standard state free energy change for self-association of the SNGpA99 glycophorin A construct is -5.7 (+/-0.3, N=5) kcal mol(-1) at 25 degrees C. Compared with previous results carried out in C(8)E(5) micellar solutions, the free energy of dimerization is 1.3 kcal mol(-1) less favorable in C14 betaine micelles. In contrast, when considered on a per-interface basis, the formation of the glycophorin A transmembrane dimer in C14 betaine micelles may be more favorable than the association of several designed transmembrane peptides.     

Effect of detergents on the association of the glycophorin a transmembrane helix

We have examined the role of the environment on the interactions between transmembrane helices using, as a model system, the dimerization of the glycophorin A transmembrane helix. In this study we have focused on micellar environments and have examined a series of detergents that include a range of alkyl chain lengths, combined with ionic, zwitterionic, and nonionic headgroups. For each we have measured how the apparent equilibrium constant depends on the detergent concentration. In two detergents we also measured the thermal sensitivity of the equilibrium constant, from which we derive the van't Hoff enthalpy and entropy. We show that several simple models are inadequate for explaining our results; however, models that include the effect of detergent concentration on detergent binding are able to account for our measurements. Our analysis suggests that the effects of detergents on helix association are due to a pair of opposing effects: an enthalpic effect, which drives association as the detergent concentration is increased and which is sensitive to the chemical nature of the detergent headgroup, opposed by an entropic effect, which drives peptide dissociation as the detergent concentration is raised. Our results also indicate that the monomer-monomer interface is relatively hydrophilic and that association within detergent micelles is driven by the enthalpy change. The wide variations in glycophorin a dimmer, stability with the detergent used, together with the realization that this results from the balance between two opposing effects, suggests that detergents might be selected that drive association rather than dissociation of peptide dimers.     

Unfolding a transmembrane helix dimer: A FRET study in mixed micelles

The exact nature of membrane protein folding and assembly is not understood in detail yet. Addition of SDS to a membrane protein dissolved in mild, non-polar detergent results in formation of mixed micelles and in subsequent denaturation of higher ordered membrane protein structures. The exact nature of this denaturation event is, however, enigmatic, and separation of an individual helix pair in mixed micelles has also not been reported yet. Here we followed unfolding of the human glycophorin A transmembrane helix dimer in mixed micelles by fluorescence spectroscopy. Energy transfer between differently labelled glycophorin A transmembrane helices decreased with increasing SDS mole fractions albeit without modifying the helicity of the peptides. The energetics and kinetics of the dimer dissociation in mixed micelles is analyzed and discussed, and the experimental data demonstrate that mixed micelles can be used as a general method to investigate unfolding of α-helical membrane proteins.     

Modulation of glycophorin A transmembrane helix interactions by lipid bilayers: molecular dynamics calculations

Starting from the glycophorin A dimer structure determined by NMR, we performed simulations of both dimer and monomer forms in explicit lipid bilayers with constant normal pressure, lateral area, and temperature using the CHARMM potential. Analysis of the trajectories in four different lipids reveals how lipid chain length and saturation modulate the structural and energetic properties of transmembrane helices. Helix tilt, helix-helix crossing angle, and helix accessible volume depend on lipid type in a manner consistent with hydrophobic matching concepts: the most relevant lipid property appears to be the bilayer thickness. Although the net helix-helix interaction enthalpy is strongly attractive, analysis of residue-residue interactions reveals significant unfavorable electrostatic repulsion between interfacial glycine residues previously shown to be critical for dimerization. Peptide volume is nearly conserved upon dimerization in all lipid types, indicating that the monomeric helices pack equally well with lipid as dimer helices do with one another. Enthalpy calculations indicate that the helix-environment interaction energy is lower in the dimer than in the monomer form, when solvated by unsaturated lipids. In all lipid environments there is a marked preference for lipids to interact with peptide predominantly through one rather than both acyl chains. Although our trajectories are not long enough to allow a full thermodynamic treatment, these results demonstrate that molecular dynamics simulations are a powerful method for investigating the protein-protein, protein-lipid, and lipid-lipid interactions that determine the structure, stability and dynamics of transmembrane alpha-helices in membranes.     

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.     

Nicotinic receptor fourth transmembrane domain: hydrogen bonding by conserved threonine contributes to channel gating kinetics

The fourth transmembrane domain (M4) of the nicotinic acetylcholine receptor (AChR) contributes to the kinetics of activation, yet its close association with the lipid bilayer makes it the outermost of the transmembrane domains. To investigate mechanistic and structural contributions of M4 to AChR activation, we systematically mutated alphaT422, a conserved residue that has been labeled by hydrophobic probes, and evaluated changes in rate constants underlying ACh binding and channel gating steps. Aromatic and nonpolar mutations of alphaT422 selectively affect the channel gating step, slowing the rate of opening two- to sevenfold, and speeding the rate of closing four- to ninefold. Additionally, kinetic modeling shows a second doubly liganded open state for aromatic and nonpolar mutations. In contrast, serine and asparagine mutations of alphaT422 largely preserve the kinetics of the wild-type AChR. Thus, rapid and efficient gating of the AChR channel depends on a hydrogen bond involving the side chain at position 422 of the M4 transmembrane domain.     

Subunit positioning and transmembrane helix organisation in the core dimer of photosystem II

Recently 3D structural models of the photosystem II (PSII) core dimer complexes of higher plants (spinach) and cyanobacteria (Synechococcus elongatus) have been derived by electron [Rhee et al. (1998) Nature 396, 283-286; Hankamer et al. (2001) J. Struct. Biol., in press] and X-ray [Zouni et al. (2001) Nature 409, 739-743] crystallography respectively. The intermediate resolutions of these structures do not allow direct identification of side chains and therefore many of the individual subunits within the structure are unassigned. Here we review the structure of the higher plant PSII core dimer and provide evidence for the tentative assignment of the low molecular weight subunits. In so doing we highlight the similarities and differences between the higher plant and cyanobacterial structures.     

The optimisation of the helix/helix interaction of a transmembrane dimer is improved by the IMPALA restraint field

A continuous membrane model (IMPALA) was previously developed to predict how hydrophobic spans of proteins insert in membranes (Mol. Mod. 2 (1996) 27). Using that membrane model, we looked for the interactions between several hydrophobic spans. We used the glycophorin A dimer as an archetype of polytopic protein to validate the approach. We find that the native complex do not dislocate when it is submitted to a 10(5) steps optimisation whereas separated spans converge back to a native-like complex in the same conditions. We also observe that IMPALA restraints are not strictly mandatory but do increase the efficiency of the procedure.     

Ala-insertion scanning mutagenesis of the glycophorin A transmembrane helix: a rapid way to map helix-helix interactions in integral membrane proteins.

Alanine insertions into the glycophorin A transmembrane helix are found to disrupt helix-helix dimerization in a way that is fully consistent with earlier saturation mutagenesis data, suggesting that Ala-insertion scanning can be used to rapidly map the approximate location of structurally and/or functionally important segments in transmembrane helices.     

Probing the structure of the Ff bacteriophage major coat protein transmembrane helix dimer by solution NMR

The transmembrane (TM) segment of the major coat protein from Ff bacteriophage has been extensively studied as an example of dimerization in detergent and lipid bilayer systems. However, almost all the information regarding this interaction has been gained through mutagenesis studies, with little direct structural information being available. To this end solution NMR has the potential to provide new insights into structure of the dimer. In order to evaluate the utility of this approach we have studied a selectively 15N-labeled peptide containing the TM segment of MCP (MCPTM) by solution NMR. This peptide was found to give rise to detergent concentration-dependent spectra that were assigned to monomeric and dimeric forms. The standard free energy of this interaction in SDS was estimated from these spectra and found to be consistent with weak but specific dimerization. In addition, similar spectra could be obtained in beta-octyl glucoside with intermolecular paramagnetic relaxation experiments demonstrating a parallel arrangement of TM helices in the dimer. In both detergents backbone chemical shift differences between monomeric and dimeric forms of MCPTM showed that the largest changes occur around its GXXXG motif. The resulting structural model is consistent with observations made for MCP mutants previously characterized in biological membranes, opening the door to detailed structural characterization of this form of MCP. These results also have general implications for the study of weakly interacting TM segments by solution NMR since the use of similar sample conditions should allow structural data to be accessed for oligomeric states from a wide range systems that undergo biologically relevant but weak associations in the membrane.     

Probing the structure of the Ff bacteriophage major coat protein transmembrane helix dimer by solution NMR

The transmembrane (TM) segment of the major coat protein from Ff bacteriophage has been extensively studied as an example of dimerization in detergent and lipid bilayer systems. However, almost all the information regarding this interaction has been gained through mutagenesis studies, with little direct structural information being available. To this end solution NMR has the potential to provide new insights into structure of the dimer. In order to evaluate the utility of this approach we have studied a selectively ¹⁵N-labeled peptide containing the TM segment of MCP (MCP_TM) by solution NMR. This peptide was found to give rise to detergent concentration-dependent spectra that were assigned to monomeric and dimeric forms. The standard free energy of this interaction in SDS was estimated from these spectra and found to be consistent with weak but specific dimerization. In addition, similar spectra could be obtained in β-octyl glucoside with intermolecular paramagnetic relaxation experiments demonstrating a parallel arrangement of TM helices in the dimer. In both detergents backbone chemical shift differences between monomeric and dimeric forms of MCP_TM showed that the largest changes occur around its GXXXG motif. The resulting structural model is consistent with observations made for MCP mutants previously characterized in biological membranes, opening the door to detailed structural characterization of this form of MCP. These results also have general implications for the study of weakly interacting TM segments by solution NMR since the use of similar sample conditions should allow structural data to be accessed for oligomeric states from a wide range systems that undergo biologically relevant but weak associations in the membrane.     

The role of a hydrogen bonding network in the transmembrane beta-barrel OMPLA

The hydrogen bonding of polar side-chains has emerged as an important theme for membrane protein interactions. The crystal structure of the dimeric state of the transmembrane beta-barrel protein outer membrane phospholipase A (OMPLA) revealed an intermolecular hydrogen bond mediated by a highly conserved glutamine side-chain (Q94). It has been shown that the introduction of a polar residue can drive the association of model helices, and by extension it was presumed that the glutamine hydrogen bond played a key role in stabilizing the OMPLA dimer. However, a thermodynamic investigation using sedimentation equilibrium ultracentrifugation in detergent micelles reveals that the hydrogen bond plays only a very modest role in stabilizing the dimer. The Q94 side-chain is hydrogen bonded intramolecularly to residues Y92 and S96, but amino acid substitutions at these positions suggest these intramolecular interactions are not responsible for attenuating the strength of the intermolecular Q94 hydrogen bond. Other substitutions suggested that hydration of the local environment around Q94 may be responsible for the modest strength of the hydrogen bond. Heat inactivation experiments with the variants suggest that the Y92-Q94-S96 network may instead be important for thermal stability of the monomer. These results highlight the context dependence and broad range of interactions that can be mediated by polar residues in membrane proteins.     

Determinants of the substrate specificity of multidrug resistance protein 1: role of amino acid residues with hydrogen bonding potential in predicted transmembrane helix 17

Human multidrug resistance protein 1 (MRP1) confers resistance to many natural product chemotherapeutic agents and actively transports structurally diverse organic anion conjugates. We previously demonstrated that two hydrogen-bonding amino acid residues in the predicted transmembrane 17 (TM17) of MRP1, Thr(1242) and Trp(1246), were important for drug resistance and 17beta-estradiol 17-(beta-d-glucuronide) (E(2)17betaG) transport. To determine whether other residues with hydrogen bonding potential within TM17 influence substrate specificity, we replaced Ser(1233), Ser(1235), Ser(1237), Gln(1239), Thr(1241), and Asn(1245) with Ala and Tyr(1236) and Tyr(1243) with Phe. Mutations S1233A, S1235A, S1237A, and Q1239A had no effect on any substrate tested. In contrast, mutations Y1236F and T1241A decreased resistance to vincristine but not to VP-16, doxorubicin, and epirubicin. Mutation Y1243F reduced resistance to all drugs tested by 2-3-fold. Replacement of Asn(1245) with Ala also decreased resistance to VP-16, doxorubicin, and epirubicin but increased resistance to vincristine. This mutation also decreased E(2)17betaG transport approximately 5-fold. Only mutation Y1243F altered the ability of MRP1 to transport both leukotriene 4 and E(2)17betaG. Together with our previous results, these findings suggest that residues with side chain hydrogen bonding potential, clustered in the cytoplasmic half of TM17, participate in the formation of a substrate binding site.     

Structural adaptation of the glycophorin A transmembrane homodimer to D-amino acid modifications

Protein-protein recognition is an essential process in life. The chemistry of these kind of interactions is predominantly stereospecific (i.e. receptor-ligand, antibody-hapten binding). Here, we investigated whether the hydrophobic nature of the membrane affects this stereospecificity. To this end, we synthesized a diastereomer analogue (2D-GPA) of the glycophorin A transmembrane helix, with two l-valine residues replaced by their d-enantiomer. This ensures a disruption of the secondary structure. We investigated the ability of the diastereomer peptide to recognize the GPA chimera in the ToxR homodimer reporting system, in vivo. The peptide demonstrated a dose-dependent dominant negative effect on the GPA transmembrane in the bacterial ToxR system, suggesting a wild-type like interaction. This result was corroborated in vitro by fluorescence energy transfer between 2D-GPA and all-l GPA. Peptide binding to the bacteria was confirmed through confocal imaging, and Western blot confirmed that ToxR GPA receptor levels are not affected by the presence of the exogenous peptide. In order to understand the structural basis for heterodimer formation, homodimer and heterodimer structures, based on the NMR 3D structure of GPA, were subjected to a molecular dynamics simulation. The resulting heterodimer structure maintained most of the original inter-helical interactions, and its structure is similar to that of the homodimer. We postulate that the need to satisfy all H-bonds can compensate for the structural strain induced by the presence of the d-amino acid residues.     

Specificity of helix packing in transmembrane dimer of the cell death factor BNIP3: a molecular modeling study

BNIP3 is a mitochondrial 19-kDa proapoptotic protein, a member of the Bcl-2 family. It has a single COOH-terminal transmembrane (TM) alpha-helical domain, which is required for membrane targeting, proapoptotic activity, hetero- and homo-dimerization in membrane. The role and the molecular details of association of TM helices of BNIP3 are yet to be established. Here, we present a molecular modeling study of helix interactions in its membrane domain. The approach combines Monte Carlo conformational search in an implicit hydrophobic slab followed by molecular dynamics simulations in a hydrated full-atom lipid bilayer. The former technique was used for exhaustive sampling of the peptides' conformational space and for generation of putative "native-like" structures of the dimer. The latter ones were taken as realistic starting points to assess stability and dynamic behavior of the complex in explicit lipid-water surrounding. As a result, several groups of tightly packed right-handed structures of the dimer were proposed. They have almost similar helix-helix interface, which includes the motif A(176)xxxG(180)xxxG(184) and agrees well with previous mutagenesis data and preliminary NMR analysis. Molecular dynamics simulations of these structures reveal perfect adaptation of most of them to heterogeneous membrane environment. A remarkable feature of the predicted dimeric structures is the occurrence of a cluster of H-bonded histidine 173 and serines 168 and 172 on the helix interface, near the N-terminus. Because of specific polar interactions between the monomers, this part of the dimer has no such dense packing as the C-terminal one, thus allowing penetration of water from the extramembrane side into the membrane interior. We propose that the ionization state of His(173) can mediate structural and dynamic properties of the dimer. This, in turn, may be related to pH-dependent proapoptotic activity of BNIP3, which is triggering on by acidosis appearing under hypoxic conditions.     

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.