Helix Dimer Research Articles

Six amino acids define a minimal dimerization sequence and stabilize a transmembrane helix dimer by close packing and hydrogen bonding

Distinct amino acid sequences have been described to mediate oligomerization of transmembrane α-helices. However, as the sequence context is crucial to determine specificity in transmembrane helix–helix interaction, the question arises how small a sequence can be without losing specificity. In the present analysis, six amino acids have been identified in the PsbF transmembrane helix dimer, which form the contact region of two interacting helices and are directly involved in helix–helix interactions. However, individual amino acids within the complex sequence pattern only together ensure sequence specificity of the analyzed transmembrane helix–helix interactions by mediating close packing and inter-helical hydrogen bonding.

The Association of Polar Residues in the DAP12 homodimer: TOXCAT and Molecular Dynamics Simulation Studies

Dimerization of the transmembrane (TM) adaptor protein DAP12 plays a key role in mediating activation signals through TM-TM association with cell-surface receptors. Herein, we apply the TOXCAT assay and molecular dynamics simulation to analyze dynamics and dimerization of the TM helix of DAP12 in the membrane bilayer. In the TOXCAT assay, we performed site-specific mutagenesis of potential dimerization motifs in the DAP12 TM domain. Instead of the common GxxxG dimerization motif, mutating either of the polar residues Asp-50 and Thr-54 significantly decreased the TOXCAT signal for the dimerization of DAP12 TM domain. Furthermore, through the conformational difference between wild-type and mutant DAP12 TM homodimers, a combined coarse-grained and atomistic molecular dynamics simulation has identified both Asp-50 and Thr-54 at the dimerization interface. The experimental and computational results of the DAP12 TM dimer are in excellent agreement with the previously reported NMR structure obtained in detergent micelles. Such a combination of dynamics simulation and cell-based experiments can be applied to produce insights at the molecular level into the TM-TM association of many other transmembrane proteins.

Tryptophan at the transmembrane–cytosolic junction modulates thrombopoietin receptor dimerization and activation

Dimerization of single-pass membrane receptors is essential for activation. In the human thrombopoietin receptor (TpoR), a unique amphipathic RWQFP motif separates the transmembrane (TM) and intracellular domains. Using a combination of mutagenesis, spectroscopy, and biochemical assays, we show that W515 of this motif impairs dimerization of the upstream TpoR TM helix. TpoR is unusual in that a specific residue is required for this inhibitory function, which prevents receptor self-activation. Mutations as diverse as W515K and W515L cause oncogenic activation of TpoR and lead to human myeloproliferative neoplasms. Two lines of evidence support a general mechanism in which W515 at the intracellular juxtamembrane boundary inhibits dimerization of the TpoR TM helix by increasing the helix tilt angle relative to the membrane bilayer normal, which prevents the formation of stabilizing TM dimer contacts. First, measurements using polarized infrared spectroscopy show that the isolated TM domain of the active W515K mutant has a helix tilt angle closer to the bilayer normal than that of the wild-type receptor. Second, we identify second-site R514W and Q516W mutations that reverse dimerization and tilt angle changes induced by the W515K and W515L mutations. The second-site mutations prevent constitutive activation of TpoR W515K/L, while preserving ligand-induced signaling. The ability of tryptophan to influence the angle and dimerization of the TM helix in wild-type TpoR and in the second-site revertants is likely associated with its strong preference to be buried in the headgroup region of membrane bilayers.

Transmembrane helix orientation influences membrane binding of the intracellular juxtamembrane domain in Neu receptor peptides

The transmembrane (TM) and juxtamembrane (JM) regions of the ErbB family receptor tyrosine kinases connect the extracellular ligand-binding domain to the intracellular kinase domain. Evidence for the role of these regions in the mechanism of receptor dimerization and activation is provided by TM-JM peptides corresponding to the Neu (or rat ErbB2) receptor. Solid-state NMR and fluorescence spectroscopy show that there are tight interactions of the JM sequence with negatively charged lipids, including phosphatidylinositol 4,5-bisphosphate, in TM-JM peptides corresponding to the wild-type receptor sequence. We observe a release of the JM sequence from the negatively charged membrane surface using peptides containing an activating V664E mutation within the TM domain or in peptides engineered to form TM helix dimers with Val664 in the interface. These results provide the basis of a mechanism for coupling ligand binding to kinase activation in the full-length receptor.

Architecture and Membrane Interactions of the EGF Receptor

Dimerization-driven activation of the intracellular kinase domains of the epidermal growth factor receptor (EGFR) upon extracellular ligand binding is crucial to cellular pathways regulating proliferation, migration, and differentiation. Inactive EGFR can exist as both monomers and dimers, suggesting that the mechanism regulating EGFR activity may be subtle. The membrane itself may play a role but creates substantial difficulties for structural studies. Our molecular dynamics simulations of membrane-embedded EGFR suggest that, in ligand-bound dimers, the extracellular domains assume conformations favoring dimerization of the transmembrane helices near their N termini, dimerization of the juxtamembrane segments, and formation of asymmetric (active) kinase dimers. In ligand-free dimers, by holding apart the N termini of the transmembrane helices, the extracellular domains instead favor C-terminal dimerization of the transmembrane helices, juxtamembrane segment dissociation and membrane burial, and formation of symmetric (inactive) kinase dimers. Electrostatic interactions of EGFR's intracellular module with the membrane are critical in maintaining this coupling.

Transmembrane Helix Association of ErbB Receptors Investigated by Coarse-Grained Metadynamics Simulations

Receptor tyrosine kinases (RTKs) are essential components of signaling at cell membranes. Members of the epidermal growth factor receptor family (ErbB receptors) play critical roles in regulating cell proliferation, differentiation and migration. Their activation involves the dimerization of alpha helices in the transmembrane (TM) domain of the receptors. Single point mutations in this region can modulate efficiency of helix dimerization and abnormal activation leads to human cancers. Therefore, it is important to determine the mode and the strength of the association in functional receptors and oncogenic variants. Various experimental and computational techniques have been recently used in this respect, which have provided substantial results. The development of innovative approaches might contribute to improve our understanding of the mechanisms of association.Here, we present a new computational protocol that combines coarse-grained molecular dynamics and metadynamics, to investigate the dimerization of TM helices of ErbB receptors. In such a way the sampling of the dimerization process is highly enhanced and the free energy profile of the association computed. First, different coarse-grained models are assessed by comparison to experimental structures of homo and heterodimers. Then, metadynamics simulations are tested using different set of parameters to analyze the reliability and the reproducibility of the results. The final results are examined in reference to the literature data. This study provides a promising protocol for systematic investigation of RTK transmembrane domains in general. Possible extensions of the methodology are finally considered to treat larger part of RTK dimers as well as more complex clustering in membrane, representing a tool of general interest in biophysical research.

Temperature Controlled Helix-Helix Interactions in Desk Minimal Sensor

DesK is a temperature sensing transmembrane protein that acts as a molecular switch to regulate membrane fluidity as a function of temperature change. The full function of DesK is modelled by a chimeric construct, denominated minimal sensor (MS), formed by the upper and lower halves of transmembrane helices 1 and 5, respectively. The signalling by MS has been explored experimentally providing data for modelling studies. Resulting from this, the current view is that it forms a dimer, switching its conformation depending on the temperature. To further investigate the molecular details of the switch mechanism, we have developed a new method for exploring the energy landscape of interaction, which allows high throughput screening of transmembrane helix dimers. The results show a clear distinction between helix-helix interactions at high and at low temperatures, providing a molecular basis for the functioning of the minimal sensor. These results form the basis for further experimental exploration, as well as for the rational design of other switching sensors.

Sequence dependent lipid-mediated effects modulate the dimerization of ErbB2 and its associative mutants

The association of transmembrane helices is an important event in several biological processes, but the factors governing association, especially the non-specific environmental effects, have still not been elucidated. Here, we use coarse-grain molecular dynamics simulations to study the association of ErbB2 transmembrane helices and three "oncogenic mutants." Self-assembly simulations and the dimerization free-energy profiles confirm an energetically-favorable dimerized state for both the wildtype and the mutants. The dissociation free energy of all three mutants is calculated to be larger than the wildtype peptide. Along with favourable protein-protein interactions, non-specific environmental effects are observed to contribute to the association. In particular, local bilayer thinning along with membrane perturbations are seen around the mutants. The membrane perturbations are reduced upon helix association, suggesting that lipid chain packing is an important driving force for helix dimerization. Our results highlight the importance of both specific as well as non-specific driving forces in the association of transmembrane helices.

Detergent Properties Influence the Stability of the Glycophorin A Transmembrane Helix Dimer in Lysophosphatidylcholine Micelles

Detergents might affect membrane protein structures by promoting intramolecular interactions that are different from those found in native membrane bilayers, and fine-tuning detergent properties can be crucial for obtaining structural information of intact and functional transmembrane proteins. To systematically investigate the influence of the detergent concentration and acyl-chain length on the stability of a transmembrane protein structure, the stability of the human glycophorin A transmembrane helix dimer has been analyzed in lyso-phosphatidylcholine micelles of different acyl-chain length. While our results indicate that the transmembrane protein is destabilized in detergents with increasing chain-length, the diameter of the hydrophobic micelle core was found to be less crucial. Thus, hydrophobic mismatch appears to be less important in detergent micelles than in lipid bilayers and individual detergent molecules appear to be able to stretch within a micelle to match the hydrophobic thickness of the peptide. However, the stability of the GpA TM helix dimer linearly depends on the aggregation number of the lyso-PC detergents, indicating that not only is the chemistry of the detergent headgroup and acyl-chain region central for classifying a detergent as harsh or mild, but the detergent aggregation number might also be important.

Prediction, refinement, and persistency of transmembrane helix dimers in lipid bilayers using implicit and explicit solvent/lipid representations: Microsecond molecular dynamics simulations of ErbB1/B2 and EphA1

All-atom simulations are carried out on ErbB1/B2 and EphA1 transmembrane helix dimers in lipid bilayers starting from their solution/DMPC bicelle NMR structures. Over the course of microsecond trajectories, the structures remain in close proximity to the initial configuration and satisfy the majority of experimental tertiary contact restraints. These results further validate CHARMM protein/lipid force fields and simulation protocols on Anton. Separately, dimer conformations are generated using replica exchange in conjunction with an implicit solvent and lipid representation. The implicit model requires further improvement, and this study investigates whether lengthy all-atom molecular dynamics simulations can alleviate the shortcomings of the initial conditions. The simulations correct many of the deficiencies. For example, excessive helix twisting is eliminated over a period of hundreds of nanoseconds. The helix tilt, crossing angles, and dimer contacts approximate those of the NMR-derived structure, although the detailed contact surface remains off-set for one of two helices in both systems. Hence, even microsecond simulations are not long enough for extensive helix rotations. The alternate structures can be rationalized with reference to interaction motifs and may represent still sought after receptor states that are important in ErbB1/B2 and EphA1 signaling.

The DPY-30 Domain and Its Flanking Sequence Mediate the Assembly and Modulation of Flagellar Radial Spoke Complexes

RIIa is known as the dimerization and docking (D/D) domain of the cyclic AMP (cAMP)-dependent protein kinase. However, numerous molecules, including radial spoke protein 2 (RSP2) in Chlamydomonas flagella, also contain an RIIa or a similar DPY-30 domain. To elucidate new roles of D/D domain-containing proteins, we investigated a panel of RSP2 mutants. An RSP2 mutant had paralyzed flagella defective in RSP2 and multiple subunits near the spokehead. New transgenic strains lacking only the DPY-30 domain in RSP2 were also paralyzed. In contrast, motility was restored in strains that lacked only RSP2's calmodulin-binding C-terminal region. These cells swam normally in dim light but could not maintain typical swimming trajectories under bright illumination. In both deletion transgenic strains, the subunits near the spokehead were restored, but their firm attachment to the spokestalk required the DPY-30 domain. We postulate that the DPY-30-helix dimer is a conserved two-prong linker, required for normal motility, organizing duplicated subunits in the radial spoke stalk and formation of a symmetrical spokehead. Further, the dispensable calmodulin-binding region appears to fine-tune the spokehead for regulation of "steering" motility in the green algae. Thus, in general, D/D domains may function to localize molecular modules for both the assembly and modulation of macromolecular complexes.

Lattice model simulations of the effects of the position of a peptide trigger segment on helix folding and dimerization

The folding and dimerization of proteins is greatly facilitated by the presence of a trigger site, a segment of amino acids that has a higher propensity for forming α-helix structure as compared to the rest of the chain. In addition to the helical propensity of each chain, dimerization can also be facilitated by interhelical interactions such as saltbridges, and interfacial contacts of different strengths. In this work, we are interested in understanding the interplay of these interactions in a model peptide system. We investigate how these different interactions influence the kinetics of dimer formation and the stability of the fully formed dimer. We use lattice model computer simulations to investigate how the effectiveness of the trigger segment and its saltbridges depends on the location along the protein primary sequence. For different positions of the trigger segment, heat capacity and free energy of unfolded and folded configurations are calculated to study the thermodynamics of folding and dimerization. The kinetics of the process is investigated by calculating characteristic folding times. The thermodynamic and kinetic data from the simulations combine to show that the dimerization process of the model system is faster when the segment with high helical propensity is located near either end of the peptide, as compared to the middle of the chain. The dependence of the stability of the dimer on the trigger segment's position is also studied. The stability can play a role in the ability of the dimer to perform a biological function that involves partial unzipping. The results on folding and dimer stability provide important insights for designing proteins that involve trigger sites.

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 b559′ 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.

Vibrational Coupling between Helices Influences the Amide I Infrared Absorption of Proteins: Application to Bacteriorhodopsin and Rhodopsin

The amide I spectrum of multimers of helical protein segments was simulated using transition dipole coupling (TDC) for long-range interactions between individual amide oscillators and DFT data from dipeptides (la Cour Jansen et al. J. Chem. Phys.2006, 125, 44312) for nearest neighbor interactions. Vibrational coupling between amide groups on different helices shift the helix absorption to higher wavenumbers. This effect is small for helix dimers (1 cm(-1)) at 10 Å distance and only moderately affected by changes in the relative orientation between the helices. However, the effect becomes considerable when several helices are bundled in membrane proteins. Particular examples are the 7-helix membrane proteins bacteriorhodopsin (BR) and rhodopsin, where the upshift is 4.3 and 5.3 cm(-1), respectively, due to interhelical coupling within a BR monomer. A further upshift of 4.0 cm(-1) occurs when BR monomers associate to trimers. We propose that interhelical vibrational coupling explains the experimentally observed unusually high wavenumber of the amide I band of BR.

Molecular Dynamics Prediction and Refinement of Transmembrane Helix Dimers

Membrane proteins represent around one fourth of all proteins encoded in genomes. Amongst these transmembrane (TM) proteins span the entire lipid bilayer, and -by allowing information to be transmitted across the membrane- play important roles in cell biology and physiology. Several mechanisms for such signal transduction have been revealed, ranging from conformational changes in protein channels and pumps to alterations in helix structures that regulate the activity of intra- and extracellular domains in TM receptors. We have begun to investigate the role of TM helix association in single pass receptors, such as the Eph and plexin receptor. Often the structures of the TM helix associated states are not known and a number of laboratories have been trying to fill this void by use of molecular modeling and dynamics for structure prediction.This work explores a plausible route to predicting structures of transmembrane helix dimers that combines two layers of modeling. The method is applied to known helix dimer structures as a validation study. Starting from a pair of unfolded peptides, thorough sampling is performed using an implicit bilayer potential to generate starting configurations that are then relaxed in an explicit bilayer. Helix crossing angles for EphA1 and ErbB1/2 dimers obtained from the explicit bilayer ensembles are in very good agreement with experiment. (Crossing angles from the implicit simulations are incorrect.) The predicted helix packing residues differ from experiment. It is possible the 100 ns simulation times are insufficient for a sampling of helix rotation and thereby a satisfactory determination of the residues involved in helix-helix packing.

Glycophorin a Helix Insertion, Positioning, and Dimerization in Model Membranes

Membrane proteins play vital and diverse roles in the cell, e.g., in membrane transport, signal transduction and molecular recognition, and inter-cellular communication in higher organisms. The membrane's lipid composition affects membrane protein helix insertion, positioning and oligomerization, but molecular mechanisms of these effects are largely unknown. We employ our novel highly mobile membrane mimetic (HMMM) model combined with molecular dynamics simulations to investigate structural and dynamic properties determining these lipid-protein interactions during insertion, hydrophobically driven positioning and dimer formation of transmembrane helices of glycophorin A (GpA).

Transmembrane helix–helix interactions are modulated by the sequence context and by lipid bilayer properties

Folding of polytopic transmembrane proteins involves interactions of individual transmembrane helices, and multiple TM helix–helix interactions need to be controlled and aligned to result in the final TM protein structure. While defined interaction motifs, such as the GxxxG motif, might be critically involved in transmembrane helix–helix interactions, the sequence context as well as lipid bilayer properties significantly modulate the strength of a sequence specific transmembrane helix–helix interaction. Structures of 11 transmembrane helix dimers have been described today, and the influence of the sequence context as well as of the detergent and lipid environment on a sequence specific dimerization is discussed in light of the available structural information. This article is part of a Special Issue entitled: Protein Folding in Membranes.

Dramatic Destabilization of Transmembrane Helix Interactions by Features of Natural Membrane Environments

Membrane proteins have evolved to fold and function in a lipid bilayer, so it is generally assumed that their stability should be optimized in a natural membrane environment. Yet optimal stability is not always in accord with optimization of function, so evolutionary pressure, occurring in a complex membrane environment, may favor marginal stability. Here, we find that the transmembrane helix dimer, glycophorin A (GpATM), is actually much less stable in the heterogeneous environment of a natural membrane than it is in model membranes and even common detergents. The primary destabilizing factors are electrostatic interactions between charged lipids and charged GpATM side chains, and nonspecific competition from other membrane proteins. These effects overwhelm stabilizing contributions from lateral packing pressure and excluded volume. Our work illustrates how evolution can employ membrane composition to modulate protein stability.

Residue-specific vibrational echoes yield 3D structures of a transmembrane helix dimer.

Two-dimensional (2D) vibrational echo spectroscopy has previously been applied to structural determination of small peptides. Here we extend the technique to a more complex, biologically important system: the homodimeric transmembrane dimer from the α chain of the integrin α(IIb)β(3). We prepared micelle suspensions of the pair of 30-residue chains that span the membrane in the native structure, with varying levels of heavy ((13)C=(18)O) isotopes substituted in the backbone of the central 10th through 20th positions. The constraints derived from vibrational coupling of the precisely spaced heavy residues led to determination of an optimized structure from a range of model candidates: Glycine residues at the 12th, 15th, and 16th positions form a tertiary contact in parallel right-handed helix dimers with crossing angles of -58° ± 9° and interhelical distances of 7.7 ± 0.5 angstroms. The frequency correlation established the dynamical model used in the analysis, and it indicated the absence of mobile water associated with labeled residues. Delocalization of vibrational excitations between the helices was also quantitatively established.

The Amyloid Formation Mechanism in Human IAPP: Dimers Have β-Strand Monomer−Monomer Interfaces

Early oligomerization of human IAPP (hIAPP) is responsible for β-cell death in the pancreas and is increasingly considered a primary pathological process linked to Type II Diabetes (T2D). Yet, the assembly mechanism remains poorly understood, largely due to the inability of conventional techniques to probe distributions or detailed structures of early oligomeric species. Here, we describe the first experimental data on the isolated and unmodified dimers of human (hIAPP) and nonamyloidogenic rat IAPP (rIAPP). The experiments reveal that the human IAPP dimers are more extended than those formed by rat IAPP and likely descend from extended monomers. Independent all-atom molecular dynamics simulations show that rIAPP forms compact helix and coil rich dimers, whereas hIAPP forms β-strand rich dimers that are generally more extended. Also, the simulations reveal that the monomer-monomer interfaces of the hIAPP dimers are dominated by β-strands and that β-strands can recruit coil or helix structured regions during the dimerization process. Our β-rich interface contrasts with an N-terminal helix-to-helix interface proposed in the literature but is consistent with existing experimental data on the self-interaction pattern of hIAPP, mutation effects, and inhibition effects of the N-methylation in the mutation region.