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.     

Self-assembly of a simple membrane protein: coarse-grained molecular dynamics simulations of the influenza M2 channel

The transmembrane (TM) domain of the M2 channel protein from influenza A is a homotetrameric bundle of α-helices and provides a model system for computational approaches to self-assembly of membrane proteins. Coarse-grained molecular dynamics (CG-MD) simulations have been used to explore partitioning into a membrane of M2 TM helices during bilayer self-assembly from lipids. CG-MD is also used to explore tetramerization of preinserted M2 TM helices. The M2 helix monomer adopts a membrane spanning orientation in a lipid (DPPC) bilayer. Multiple extended CG-MD simulations (5 × 5 μs) were used to study the tetramerization of inserted M2 helices. The resultant tetramers were evaluated in terms of the most populated conformations and the dynamics of their interconversion. This analysis reveals that the M2 tetramer has 2× rotationally symmetrical packing of the helices. The helices form a left-handed bundle, with a helix tilt angle of ∼16°. The M2 helix bundle generated by CG-MD was converted to an atomistic model. Simulations of this model reveal that the bundle's stability depends on the assumed protonation state of the H37 side chains. These simulations alongside comparison with recent x-ray (3BKD) and NMR (2RLF) structures of the M2 bundle suggest that the model yielded by CG-MD may correspond to a closed state of the channel.     

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.     

Molecular Dynamics Simulations of Micelle Formation around Dimeric Glycophorin A Transmembrane Helices

Insertion and formation of membrane proteins involves the interaction of protein helices with one another in lipid environments. Researchers have studied glycophorin A (GpA) transmembrane helices embedded in sodium dodecyl sulfate (SDS) micelles to identify contacts significant for helix dimerization. However, a detailed picture of the conformation and dynamics of the GpA-SDS system cannot be obtained solely through experiment. Molecular dynamics simulations of SDS and a GpA dimer can provide an atomic-level picture of SDS aggregation and helix association. We report 2.5-ns simulations of GpA wild-type and mutants in a preformed micelle as well as a 32-ns simulation showing the formation of a complete micelle around wild-type GpA from an initially random placement of SDS molecules in an aqueous environment. In the latter case, an initial instability of GpA helices in water is reversed after the helices become surrounded by SDS. The properties of the spontaneously formed micelle surrounding the GpA are indistinguishable from those of the preformed micelle surrounding the GpA dimer.     

Transmembrane Helix Association Affinity Can Be Modulated by Flanking and Noninterfacial Residues

The GxxxG sequence motif mediates the association of transmembrane (TM) helices by providing a site of close contact between them. However, it is not sufficient for strong association. For example, both bacteriophage M13 major coat protein (MCP) and human erythrocyte protein glycophorin A (GpA) contain a GxxxG motif in their TM domains and form a homodimer, but the association affinity of MCP, measured by the ToxCAT in vivo assay, is dramatically weaker than that of GpA. Even when all interfacial residues of MCP were substituted for those of GpA (MCP-GpA), association remained significantly weaker than in GpA. Here we provide an explanation for these experimental observations using molecular dynamics simulations in an implicit membrane (IMM1-GC). The association free energies of GpA29 (GpA with 29 residues all from the wild-type sequence), GpA15p11 (GpA with 15 residues from the wild-type sequence plus 11 flanking residues from the ToxCAT construct), MCP, and MCP-GpA TM helices were calculated and compared. MCP and MCP-GpA have the same flanking residues used in the ToxCAT assay as those in GpA15p11, but the position of the flanking residues relative to the GxxxG motif is different. The calculated association free energies follow experimental observations: the association affinity of MCP-GpA falls between those of GpA15p11 and MCP wild-type. MCP exhibits an equally strong interhelical interaction in the TM domain. A major reason for the weaker association of MCP in the calculations was the noninterfacial residue Lys-40, which in the dimer structure is forced to be buried in the membrane interior. To alleviate the desolvation cost, in MCP and MCP-GpA dimers, Lys-40 gets deprotonated. A second factor that modulates association affinity is the flanking residues. Thanks to them, GpA15p11 exhibits a much stronger association affinity than GpA29. The positioning of the flanking residues is also important, as evidenced by the difference in association affinity between MCP and MCP-GpA on one hand and GpA15p11 on the other. Thus, residues outside the contact interface can exert a significant influence on transmembrane helix association affinity.     

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.     

Substitution of the insulin receptor transmembrane domain with that of glycophorin A inhibits insulin action.

To study the role of transmembrane (TM) domains interactions in the activation of the insulin receptor, we have replaced the insulin receptor TM domain with that of glycophorin A (GpA), an erythrocyte protein that spontaneously forms detergent-resistant dimers through TM-TM interactions. Insulin receptor cDNA sequences with the TM domain replaced by that of GpA were constructed and stably transfected in CHO cells. Insulin binding to cells and solubilized receptors was not modified. Electrophoresis after partial reduction of disulfide bonds revealed an altered structure for the soluble chimeric receptors, seen as an altered mobility apparently due to increased interactions between the beta subunits of the receptor. Insulin signaling was markedly decreased for cells transfected with chimeric receptors compared with cells transfected with normal receptors. A decrease in insulin-induced receptor kinase activity was observed for solubilized chimeric receptors. In conclusion, substitution by the native GpA TM domain of the insulin receptor results in structurally modified chimeric receptors that are unable to transmit the insulin signal properly. It is hypothesized that this substitution may impose structural constraints that prevent the proper changes in conformation necessary for activation of the receptor kinase. Other mutants modifying the structure or the membrane orientation of the glycophorin A TM domain are required to better understand these constraints.     

Transmembrane helix dimerization: beyond the search for sequence motifs

Studies of the dimerization of transmembrane (TM) helices have been ongoing for many years now, and have provided clues to the fundamental principles behind membrane protein (MP) folding. Our understanding of TM helix dimerization has been dominated by the idea that sequence motifs, simple recognizable amino acid sequences that drive lateral interaction, can be used to explain and predict the lateral interactions between TM helices in membrane proteins. But as more and more unique interacting helices are characterized, it is becoming clear that the sequence motif paradigm is incomplete. Experimental evidence suggests that the search for sequence motifs, as mediators of TM helix dimerization, cannot solve the membrane protein folding problem alone. Here we review the current understanding in the field, as it has evolved from the paradigm of sequence motifs into a view in which the interactions between TM helices are much more complex. This article is part of a Special Issue entitled: Membrane protein structure and function.     

The affinity of GXXXG motifs in transmembrane helix-helix interactions is modulated by long-range communication

Sequence motifs are responsible for ensuring the proper assembly of transmembrane (TM) helices in the lipid bilayer. To understand the mechanism by which the affinity of a common TM-TM interactive motif is controlled at the sequence level, we compared two well characterized GXXXG motif-containing homodimers, those formed by human erythrocyte protein glycophorin A (GpA, high-affinity dimer) and those formed by bacteriophage M13 major coat protein (MCP, low affinity dimer). In both constructs, the GXXXG motif is necessary for TM-TM association. Although the remaining interfacial residues (underlined) in GpA (LIXXGVXXGVXXT) differ from those in MCP (VVXXGAXXGIXXF), molecular modeling performed here indicated that GpA and MCP dimers possess the same overall fold. Thus, we could introduce GpA interfacial residues, alone and in combination, into the MCP sequence to help decrypt the determinants of dimer affinity. Using both in vivo TOXCAT assays and SDS-PAGE gel migration rates of synthetic peptides derived from TM regions of the proteins, we found that the most distal interfacial sites, 12 residues apart (and approximately 18 A in structural space), work in concert to control TM-TM affinity synergistically.     

Polar residue tagging of transmembrane peptides

Studies that focus on packing interactions between transmembrane (TM) helices in membrane proteins would greatly benefit from the ability to investigate their association and packing interactions in multi-spanning TM domains. However, the production, purification, and characterization of such units have been impeded by their high intrinsic hydrophobicity. We describe the polar tagging approach to biophysical analysis of TM segment peptides, where incorporation of polar residues of suitable type and number at one or both peptide N- and C-termini can serve to counterbalance the apolar nature of a native TM segment, and render it aqueous-soluble. Using the native TM sequences of the human erythrocyte protein glycophorin A (GpA) and bacteriophage M13 major coat protein (MCP), properties of tags such as Lys, His, Asp, sarcosine, and Pro-Gly are evaluated, and general procedures for tagging a given TM segment are presented. Gel-shift assays on sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) establish that various tagged GpA TM segments spontaneously insert into micellar membranes, and exhibit native TM dimeric states. Sedimentation equilibrium analytical centrifugation is used to confirm that Lys-tagged GpA peptides retain the native dimer state. Two-dimensional nuclear magnetic resonance (NMR) spectroscopy studies on Lys-tagged TM MCP peptides selectively enriched with N-15 illustrate the usefulness of this system for evaluating monomer-dimer equilibria in micelle environments. The overall results suggest that polar-tagging of hydrophobic (TM) peptides approach constitutes a valuable tool for the study of protein-protein interactions in membranes.     

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

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

Self-association of models of transmembrane domains of ErbB receptors in a lipid bilayer

Association of transmembrane (TM) helices is facilitated by the close packing of small residues present along the amino-acid sequence. Extensive studies have established the role of such small residue motifs (GxxxG) in the dimerization of Glycophorin A (GpA) and helped to elucidate the association of TM domains in the epidermal growth factor family of receptors (ErbBs). Although membrane-mediated interactions are known to contribute under certain conditions to the dimerization of proteins, their effect is often considered nonspecific, and any potential dependence on protein sequence has not been thoroughly investigated. We recently reported that the association of GpA is significantly assisted by membrane-induced contributions as quantified in different lipid bilayers. Herein we extend our studies to explore the origin of these effects and quantify their magnitude using different amino-acid sequences in the same lipid environment. Using a coarse-grained model that accounts for amino-acid specificity, we perform extensive parallel Monte Carlo simulations of ErbB homodimerization in dipalmitoyl-phosphatidylcholine lipid bilayers. A detailed characterization of dimer formation and estimates of the free energy of association reveal that the TM domains show a significant affinity to self-associate in lipid bilayers, in qualitative agreement with experimental findings. The presence of GxxxG motifs enhances favorable protein-protein interactions at short separations. However, the lipid-induced attraction presents a more complex character than anticipated. Depending on the interfacial residues, lipid-entropic contributions support a decrease of separation or a parallel orientation to the membrane normal, with important implications for protein function.     

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

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

Alpha-hairpin stability and folding of transmembrane segments

Molecular Dynamics (MD) simulations at low dielectric constant have been carried out for peptides matching the double spanning segments of transmembrane proteins. Different folding dynamics have been observed. The peptides folded into the stable helix-turn-helix conformation-alpha-hairpin-with antiparallel-oriented strands or unstable alpha-hairpin conformation that unfolded later into the straight helical structure. The peptide having flexible residues in the TM helices often misfolded into a tangled structure that can be avoided by restricting the flexibility of these residues. General conclusions can be drawn from the observed folding dynamics. The stability and folding of some double spanning transmembrane fragments are self-assembling. The following and/or neighboring peptide chains of the protein may support the stability of the hairpin structure of other fragments. The stability of the TM helices containing flexible residues could be maintained due to contacts with neighboring TM segments.     

Beta-branched residues adjacent to GG4 motifs promote the efficient association of glycophorin A transmembrane helices

Protein transmenembrane (TM) segments participating in helix-helix packing commonly contain small residue patterns (termed GG4 or "small-xxx-small" motifs) at i and i + 4 positions. Within many TM segments - such as the glycophorin A (GpA) sequence L75IxxGVxxGVxxT87- the G17y-xxx-Gly83 motif often occurs in combination with large, usually beta3-branched aliphatic residues at adjacent positions, typified here by Val30 and Val84 residues. To explore the importance of local P-branched character on GpA dimerization, we made systematic replacements to all 16 combinations of single or double Ile, Leu, and AIa residues at GpA TM Val/Val positions 80 and 84. Using the TOXCAT system to assay self-oligomerization in the Escherichia coli inner membrane--we observed that (i) combinations of Val and lie residues maintained, or improved dimerization levels; (ii) single Ala or Leu mutant combinations with Val or Ile maintained near-wild type dimerization affinities; and (iii) in the absence of beta-branching, i.e., Leu/Leu, Ala/Ala and Ala/Leu combinations, GpA dimerization was significantly diminished. An apparent capacity of lle-containing mutants to increase GpA dimerization versus WT likely arises from improved van der Waals packing (vs. Val) within the locus of helix contact, consistent with correlations we noted in lipid accessibility measurements. Examination of several synthetic peptides with sequences corresponding to selected GpA mutants (VV VI, IV II, and LL) confirmed their dimerization on sodium dodecylsulfate polyacrylamide gel electrophoresis (SDS-PAGE). The overall results reinforce the importance of a beta-branch-containing "ridge" residue to complement a "small-xxx-small groove" in promotion of TM-TM interactions.     

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.     

Probing the energy landscape of the membrane protein bacteriorhodopsin

The folding and stability of transmembrane proteins is a fundamental and unsolved biological problem. Here, single bacteriorhodopsin molecules were mechanically unfolded from native purple membranes using atomic force microscopy and force spectroscopy. The energy landscape of individual transmembrane alpha helices and polypeptide loops was mapped by monitoring the pulling speed dependence of the unfolding forces and applying Monte Carlo simulations. Single helices formed independently stable units stabilized by a single potential barrier. Mechanical unfolding of the helices was triggered by 3.9-7.7 A extension, while natural unfolding rates were of the order of 10(-3) s(-1). Besides acting as individually stable units, helices associated pairwise, establishing a collective potential barrier. The unfolding pathways of individual proteins reflect distinct pulling speed-dependent unfolding routes in their energy landscapes. These observations support the two-stage model of membrane protein folding in which alpha helices insert into the membrane as stable units and then assemble into the functional protein.     

Polar/Ionizable Residues in Transmembrane Segments: Effects on Helix-Helix Packing

The vast majority of membrane proteins are anchored to biological membranes through hydrophobic α-helices. Sequence analysis of high-resolution membrane protein structures show that ionizable amino acid residues are present in transmembrane (TM) helices, often with a functional and/or structural role. Here, using as scaffold the hydrophobic TM domain of the model membrane protein glycophorin A (GpA), we address the consequences of replacing specific residues by ionizable amino acids on TM helix insertion and packing, both in detergent micelles and in biological membranes. Our findings demonstrate that ionizable residues are stably inserted in hydrophobic environments, and tolerated in the dimerization process when oriented toward the lipid face, emphasizing the complexity of protein-lipid interactions in biological membranes.     

Assessing the relative stability of dimer interfaces in g protein-coupled receptors

Considerable evidence has accumulated in recent years suggesting that G protein-coupled receptors (GPCRs) associate in the plasma membrane to form homo- and/or heteromers. Nevertheless, the stoichiometry, fraction and lifetime of such receptor complexes in living cells remain topics of intense debate. Motivated by experimental data suggesting differing stabilities for homomers of the cognate human β1- and β2-adrenergic receptors, we have carried out approximately 160 microseconds of biased molecular dynamics simulations to calculate the dimerization free energy of crystal structure-based models of these receptors, interacting at two interfaces that have often been implicated in GPCR association under physiological conditions. Specifically, results are presented for simulations of coarse-grained (MARTINI-based) and atomistic representations of each receptor, in homodimeric configurations with either transmembrane helices TM1/H8 or TM4/3 at the interface, in an explicit lipid bilayer. Our results support a definite contribution to the relative stability of GPCR dimers from both interface sequence and configuration. We conclude that β1- and β2-adrenergic receptor homodimers with TM1/H8 at the interface are more stable than those involving TM4/3, and that this might be reconciled with experimental studies by considering a model of oligomerization in which more stable TM1 homodimers diffuse through the membrane, transiently interacting with other protomers at interfaces involving other TM helices.     

Functional analysis of the transmembrane (TM) domain of the Autographa californica multicapsid nucleopolyhedrovirus GP64 protein: substitution of heterologous TM domains

GP64, the major envelope glycoprotein of the Autographa californica multicapsid nucleopolyhedrovirus (AcMNPV) budded virion, is important for host cell receptor binding and mediates low-pH-triggered membrane fusion during entry by endocytosis. In the current study, we examined the functional role of the AcMNPV GP64 transmembrane (TM) domain by replacing the 23-amino-acid GP64 TM domain with corresponding TM domain sequences from a range of viral and cellular type I membrane proteins, including Orgyia pseudotsugata MNPV (OpMNPV) GP64 and F, thogotovirus GP75, Lymantria dispar MNPV (LdMNPV) F, human immunodeficiency virus type 1 (HIV-1) GP41, human CD4 and glycophorin A (GpA), and influenza virus hemagglutinin (HA), and with a glycosylphosphatidylinositol (GPI) anchor addition sequence. In transient expression experiments with Sf9 cells, chimeric GP64 proteins containing either a GPI anchor or TM domains from LdMNPV F or HIV-1 GP41 failed to localize to the cell surface and thus appear to be incompatible with either GP64 structure or cell transport. All of the mutant constructs detected at the cell surface mediated hemifusion (outer leaflet merger) upon low-pH treatment, but only those containing TM domains from CD4, GpA, OpMNPV GP64, and thogotovirus GP75 mediated pore formation and complete membrane fusion activity. This supports a model in which partial fusion (hemifusion) proceeds by a mechanism that is independent of the TM domain and the TM domain participates in the enlargement or expansion of fusion pores after hemifusion. GP64 proteins containing heterologous TM domains mediated virion budding with dramatically differing levels of efficiency. In addition, chimeric GP64 proteins containing TM domains from CD4, GpA, HA, and OpMNPV F were incorporated into budded virions but were unable to rescue the infectivity of a gp64 null virus, whereas those with TM domains from OpMNPV GP64 and thogotovirus GP75 rescued infectivity. These results show that in addition to its basic role in membrane anchoring, the GP64 TM domain is critically important for GP64 trafficking, membrane fusion, virion budding, and virus infectivity. These critical functions were replaced only by TM domains from related viral membrane proteins.