Helix Dimer Research Articles

Studies of Receptor Tyrosine Kinase Transmembrane Domain Interactions: The EmEx-FRET Method

The energetics of transmembrane (TM) helix dimerization in membranes and the thermodynamic principles behind receptor tyrosine kinase (RTK) TM domain interactions during signal transduction can be studied using Förster resonance energy transfer (FRET). For instance, FRET studies have yielded the stabilities of wild-type fibroblast growth factor receptor 3 (FGFR3) TM domains and two FGFR3 pathogenic mutants, Ala391Glu and Gly380Arg, in the native bilayer environment. To further our understanding of the molecular mechanisms of deregulated FGFR3 signaling underlying different pathologies, we determined the effect of the Gly382Asp FGFR3 mutation, identified in a multiple myeloma cell line, on the energetics of FGFR3 TM domain dimerization. We measured dimerization energetics using a novel FRET acquisition and processing method, termed "emission-excitation FRET (EmEx-FRET)," which improves the precision of thermodynamic measurements of TM helix association. The EmEx-FRET method, verified here by analyzing previously published data for wild-type FGFR3 TM domain, should have broad utility in studies of protein interactions, particularly in cases when the concentrations of fluorophore-tagged molecules cannot be controlled.

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.

Mutagenesis data in the automated prediction of transmembrane helix dimers

We present a molecular modeling protocol that selects modeled protein structures based on experimental mutagenesis results. The computed effect of a point mutation should be consistent with its experimental effect for correct models; mutations that do not affect protein stability and function should not affect the computed energy of a correct model while destabilizing mutations should have unfavorable computed energies. On the other hand, an incorrect model will likely display computed energies that are inconsistent with experimental results. We added terms to our energy function which penalize models that are inconsistent with experimental results. This creates a selective advantage for models that are consistent with experimental results in the Monte Carlo simulated annealing protocol we use to search conformational space. We calibrated our protocol to predict the structure of transmembrane helix dimers using glycophorin A as a model system. Inclusion of mutational data in this protocol compensates for the limitations of our force field and the limitations of our conformational search. We demonstrate an application of this structure prediction protocol by modeling the transmembrane region of the BNIP3 apoptosis factor.

Design of a Functional Membrane Protein by Engineering a Heme-Binding Site in Glycophorin A

We have designed a functional model membrane protein by engineering a bis-Histidine heme-binding site into a natural membrane protein, glycophorin A (GpA), structurally characterized by the dimerization of a single transmembrane helix. Out of the 32 residues comprising the transmembrane helix of GpA, five amino acids were mutated; the resulting protein, ME1, has been characterized in dodecyl phosphocholin (DPC) micelles by UV-vis, CD spectroscopy, gel electrophoresis, and analytical ultracentrifugation. ME1 binds heme with sub-micromolar affinity and maintains the highly helical secondary structure and dimeric oligomerization state of GpA. The ME1-Heme complex exhibits a redox potential of -128 +/- 2 mV vs SHE, indicating that the heme resides in a hydrophobic environment and is well shielded from the aqueous phase. Moreover, ME1 catalyzes the hydrogen peroxide dependent oxidation of organic substrates such as TMB (2,2',5,5'-tetramethyl-benzidine). This protein may provide a useful framework to investigate how the protein matrix tunes the cofactor properties in membrane proteins.

Helix-dimer–Random-coil Thermal Switching Process of Ethynylhelicene Heptamer Highly Sensitive to Its Environment

Abstract An acyclic ethynylhelicene heptamer reversibly switches its structure in solution between a helix dimer at low temperatures and random coils at high temperatures. The structural change can be repeated by alternating heating and cooling. The process monitored by CD exhibited various patterns of Δε/time profiles depending on changes in temperature, concentration, and solvent type. This single compound provides diverse outputs upon thermal input.

Computational Design of a Model for the Platelet Integrin αIIbβ3.

Integrins, the principal receptors linking the cytoskeleton to extracellular ligands, are responsible, in large part, for platelet adhesion and aggregation. Like many membrane spanning proteins, the molecular structure of integrins has been difficult to determine, although crystal and NMR structures are available for water soluble portions of αVβ3 and αIIbβ3. In addition, models for the αIIb and β3 transmembrane (TM) heterodimer have been obtained using computational methods, but with conflicting results. Accordingly, we employed two different modeling strategies to distinguish between existing models. First, we used a Monte Carlo-simulated annealing (MCSA) protocol that features a selective advantage for structures that are consistent with previously characterized αIIb and β3 TM domain point mutations. During each step of the protocol, we calculated the effects of selected mutations and assigned a penalty to the step if it produced a structure that was inconsistent with experimental results. Second, we threaded the αIIb and β3 TM sequences onto a database of membrane helix dimers and rank-ordered the structures based on their calculated energies and their agreement with the cysteine crosslinking experiments reported by Luo et al (PLoS Biol 2:776–86, 2004). The MCSA and threaded TM models were similar with a root mean squared deviation (RMSD) of 1.5 Å. Moreover, they differed from a model reported by Kay Gottschalk (Structure 13:703–12, 2005) most notably by an ~50° rotation of the αIIb TM helix about its helical axis. This positions αIIb residue Leu980 in the TM interface and in van der Waals contact with β3 residue Gly708, a result consistent with mutagenesis experiments. Further, the αIIb and β3 TM helices pack along conserved glycine residues (Gly972 and Gly976) in the GXXXG motif of αIIb with a right-handed crossing angle of −18°. To create a model for the full-length integrin, the N-termini of the TM model were connected to the crystal structure of the αVβ3 extracellular domain and its C-termini were connected to the corresponding NMR structures of the αIIb and β3 cytosolic domains. The model predicts an αIIb Arg995-β3 Asp723 salt bridge in the membrane proximal cytosolic domains and a change in the secondary structure of the extracellular domains from an extended to a helical conformation that occurs at Pro965 in αIIb and Pro693 in β3, just proximal to the TM domains. Finally, the full length model provides insight into distal interactions since its extracellular portion is held close to the membrane and is therefore sterically inaccessible to large ligands such as fibrinogen, whereas cytosolic regulatory sites such as the talin binding site in β3 are accessible and not sequestered by α/β interactions. In conclusion, we have modeled the αIIbβ3 TM heterodimer using two fundamentally different computational strategies, both of which converge upon similar structures that were used to construct a model for the full-length integrin.

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.

Amide vibrations are delocalized across the hydrophobic interface of a transmembrane helix dimer.

The tertiary interactions between amide-I vibrators on the separate helices of transmembrane helix dimers were probed by ultrafast 2D vibrational photon echo spectroscopy. The 2D IR approach proves to be a useful structural method for the study of membrane-bound structures. The 27-residue human erythrocyte protein Glycophorin A transmembrane peptide sequence: KKITLIIFG(79)VMAGVIGTILLISWG(94)IKK was labeled at G(79) and G(94) with (13)C=(16)O or (13)C=(18)O. The isotopomers and their 50:50 mixtures formed helical dimers in SDS micelles whose 2D IR spectra showed components from homodimers when both helices had either (13)C=(16)O or (13)C=(18)O substitution and a heterodimer when one had (13)C=(16)O substitution and the other had (13)C=(18)O substitution. The cross-peaks in the pure heterodimer 2D IR difference spectrum and the splitting of the homodimer peaks in the linear IR spectrum show that the amide-I mode is delocalized across a pair of helices. The excitation exchange coupling in the range 4.3-6.3 cm(-1) arises from through-space interactions between amide units on different helices. The angle between the two Gly(79) amide-I transition dipoles, estimated at 103 degrees from linear IR spectroscopy and 110 degrees from 2D IR spectroscopy, combined with the coupling led to a structural picture of the hydrophobic interface that is remarkably consistent with results from NMR on helix dimers. The helix crossing angle in SDS is estimated at 45 degrees. Two-dimensional IR spectroscopy also sets limits on the range of geometrical parameters for the helix dimers from an analysis of the coupling constant distribution.

The Stability of Transmembrane Helix Interactions Measured in a Biological Membrane

Despite some promising progress in the understanding of membrane protein folding and assembly, there is little experimental information regarding the thermodynamic stability of transmembrane helix interactions and even less on the stability of transmembrane helix-helix interactions in a biological membrane. Here we describe an approach that allows quantitative measurement of transmembrane helix interactions in a biological membrane, and calculation of changes in the interaction free energy resulting from substitution of single amino acids. Dimerization of several variants of the glycophorin A transmembrane domain are characterized and compared to the wild-type (wt) glycophorin A transmembrane helix dimerization. The calculated DeltaDeltaG(app) values are further compared with values found in the literature. In addition, we compare interactions between the wt glycophorin A transmembrane domain and helices in which critical glycine residues are replaced by alanine or serine, respectively. The data demonstrate that replacement of the glycine residues by serine is less destabilizing than replacement by alanine with a DeltaDeltaG(app) value of about 0.4 kcal/mol. Our study comprises the first measurement of a transmembrane helix interaction in a biological membrane, and we are optimistic that it can be further developed and applied.

Measurement of Thermodynamic Parameters for Hydrophobic Mismatch 1: Self-Association of a Transmembrane Helix†

Membrane partitioning and self-association of transmembrane helices are crucial thermodynamic steps for membrane protein folding, although experimental difficulties have hampered quantitative estimations of related thermodynamic parameters, especially in lipid bilayer environments. This article reports for the first time, the complete set of thermodynamic parameters (DeltaG, DeltaH, DeltaS, and DeltaC(p)) for the formation of the antiparallel dimer of the inert hydrophobic model transmembrane helix X-(AALALAA)(3)-Y (X = 7-nitro-2-1, 3-benzoxadiazol-4-yl (NBD) and Y = NH(2) (I) or X = Ac and Y = NHCH(2)CH(2)-S-N-[4-[[4-(dimethylamino)phenyl]azo]phenyl]maleimide (DABMI) (II)) in dimonounsaturated phosphocholine lipid bilayers with different hydrophobic thicknesses (C14-C22) at 5-55 degrees C, as evaluated by fluorescence resonance energy transfer from I to II. Stronger dimerization was observed in thicker membranes and at lower temperatures (DeltaG = -9 to -26 kJ mol(-)(1)), driven by large negative DeltaH values (-18 to -80 kJ mol(-)(1)). Fourier transform infrared-polarized spectroscopy revealed that the peptide formed a stable transmembrane helix with an orientation angle of approximately 15 degrees in all bilayers without significant effects on lipid structures, suggesting that the depth to which the helix termini penetrate changes depending on the degree of hydrophobic mismatch. The enthalpy changes for helix-helix interactions can be well explained by the electrostatic interactions between helix macrodipoles in different dielectric environments. The new concept of dipole-dipole interaction as a basic driving force of helix dimerization will become a basis for understanding the structural and functional modifications in response to hydrophobic mismatch.

Insertion and assembly of membrane proteins via simulation.

Interactions of lipids are central to the folding and stability of membrane proteins. Coarse-grained molecular dynamics simulations have been used to reveal the mechanisms of self-assembly of protein/membrane and protein/detergent complexes for representatives of two classes of membrane protein, namely, glycophorin (a simple alpha-helical bundle) and OmpA (a beta-barrel). The accuracy of the coarse-grained simulations is established via comparison with the equivalent atomistic simulations of self-assembly of protein/detergent micelles. The simulation of OmpA/bilayer self-assembly reveals how a folded outer membrane protein can be inserted in a bilayer. The glycophorin/bilayer simulation supports the two-state model of membrane folding, in which transmembrane helix insertion precedes dimer self-assembly within a bilayer. The simulations also suggest that a dynamic equilibrium exists between the glycophorin helix monomer and dimer within a bilayer. The simulated glycophorin helix dimer is remarkably close in structure to that revealed by NMR. Thus, coarse-grained methods may help to define mechanisms of membrane protein (re)folding and will prove suitable for simulation of larger scale dynamic rearrangements of biological membranes.

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.

Förster resonance energy transfer in liposomes: Measurements of transmembrane helix dimerization in the native bilayer environment

The lipid bilayer vesicle is a model of the cellular membrane. Even in this simple system, however, measuring the thermodynamics of membrane protein association is a challenge. Here we discuss Förster resonance energy transfer (FRET) in liposomes as a method to probe the dimerization of transmembrane helices in a membrane environment. Although the measurements are labor intensive, FRET in liposomes can be measured accurately provided that attention is paid to sample homogeneity and sample equilibration. One must also take into account statistical expectations and the FRET that results from random colocalization of donors and acceptors in the bilayer. Without careful attention to these details, misleading results are easy to obtain in membrane FRET experiments. The results that we obtain in model systems are reproducible and depend solely on the concentration of the protein in the bilayer (i.e., on the protein-to-lipid ratio), thereby yielding thermodynamic parameters that are directly relevant to processes in biological membranes.

Rendezvous in a membrane: close packing, hydrogen bonding, and the formation of transmembrane helix oligomers

The interaction of transmembrane α-helices is promoted by a detailed fit between two helical surfaces, which results in close packing and van der Waals interactions of amino acid side chains between two helices. Recent studies additionally indicate an important role of hydrogen bonding for mediating and stabilizing transmembrane helix–helix interactions. The interplay between close packing and electrostatic interactions in influencing the specificity of helix–helix interactions on the one hand and the stability of an existing interaction on the other hand is still unknown. Here, we suggest that close packing mainly determines the specificity of a helix–helix interaction, whereas hydrogen bonding is important for stabilization of a preformed helix dimer.

Imaging forster resonance energy transfer measurements of transmembrane helix interactions in lipid bilayers on a solid support.

We utilize supported lipid/protein bilayers to probe the dimerization of transmembrane (TM) helices in a membrane environment. The bilayers are formed by incubating substrates with liposomes containing the proteins, and are characterized using fluorescence recovery after photobleaching and imaging Forster resonance energy transfer (FRET). We show that the FRET signal, as a measure of TM helix dimerization, is the same in suspended liposomes and in surface-supported bilayers. This work is the first step toward the development of a new tool for probing the association of TM helices in lipid bilayers.

TEL/AML1 shows dominant-negative effects over TEL as well as AML1

The TEL/AML1 chimeric gene is generated by the t(12;21) translocation in pre-B cell acute lymphoblastic leukemia. TEL/AML1 consists of the helix–loop–helix (HLH) dimerization domain from TEL and almost the entire of AML1, but loses the ETS DNA-binding domain from TEL. Dominant-negative effects of TEL/AML1 over wild-type-AML1 are believed to trigger the development of this type of leukemia. However, it could also be possible that TEL/AML1 affects wild-type-TEL’s molecular and tumor suppressive functions through the HLH domain. To test this possibility, we first confirmed that TEL/AML1 associates with wild-type-TEL. TEL/AML1 neither bound to the ETS-binding consensus site nor repressed transcription through it. Regardless, this prevented wild-type-TEL-induced transcriptional repression. Moreover, TEL/AML1 concomitantly inhibited wild-type-TEL-induced growth suppression and wild-type-AML1-mediated transforming activity in NIH3T3 cells. All these data indicate that TEL/AML1 exerts dominant-interfering effects on both AML1 and TEL, and that expression of TEL/AML1 could result in inactivation of TEL’s tumor suppressive functions in t(12;21)-carrying leukemia.

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.

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 A 2 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.

Helix Packing Moments Reveal Diversity and Conservation in Membrane Protein Structure

Helical membrane proteins are more tightly packed and the packing interactions are more diverse than those found in helical soluble proteins. Based on a linear correlation between amino acid packing values and interhelical propensity, we propose the concept of a helix packing moment to predict the orientation of helices in helical membrane proteins and membrane protein complexes. We show that the helix packing moment correlates with the helix interfaces of helix dimers of single pass membrane proteins of known structure. Helix packing moments are also shown to help identify the packing interfaces in membrane proteins with multiple transmembrane helices, where a single helix can have multiple contact surfaces. Analyses are described on class A G protein-coupled receptors (GPCRs) with seven transmembrane helices. We show that the helix packing moments are conserved across the class A family of GPCRs and correspond to key structural contacts in rhodopsin. These contacts are distinct from the highly conserved signature motifs of GPCRs and have not previously been recognized. The specific amino acid types involved in these contacts, however, are not necessarily conserved between subfamilies of GPCRs, indicating that the same protein architecture can be supported by a diverse set of interactions. In GPCRs, as well as membrane channels and transporters, amino acid residues with small side-chains (Gly, Ala, Ser, Cys) allow tight helix packing by mediating strong van der Waals interactions between helices. Closely packed helices, in turn, facilitate interhelical hydrogen bonding of both weakly polar (Ser, Thr, Cys) and strongly polar (Asn, Gln, Glu, Asp, His, Arg, Lys) amino acid residues. We propose the use of the helix packing moment as a complementary tool to the helical hydrophobic moment in the analysis of transmembrane sequences.

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 °C. Compared with previous results carried out in C8E5 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.