The strong dimerization of the transmembrane domain of the fibroblast growth factor receptor (FGFR) is modulated by C‐terminal juxtamembrane residues

The fibroblast growth factor receptor 3 (FGFR3) is a member of the FGFR subfamily of the receptor tyrosine kinases (RTKs) involved in signaling across the plasma membrane. Generally, ligand binding leads to receptor dimerization and activation. Dimerization involves the transmembrane (TM) domain, where mutations can lead to constitutive activation in certain cancer types and also in skeletal malformations. Thus, it has been postulated that FGFR homodimerization must be inherently weak to allow regulation, a feature reminiscent of α and β integrin TM interactions. However, we show herein that in FGFR3‐TM, four C‐terminal residues, CRLR, have a profound destabilizing effect in an otherwise strongly dimerizing TM peptide. In the absence of these four residues, the dimerizing propensity of FGFR3‐TM is comparable to glycophorin, as shown using various detergents. In addition, the expected enhanced dimerization induced by the mutation associated to the Crouzon syndrome A391E, was observed only when these four C‐terminal residues were present. In the absence of these four residues, A391E was dimer‐destabilizing. Finally, using site specific infrared dichroism and convergence with evolutionary conservation data, we have determined the backbone model of the FGFR3‐TM homodimer in model lipid bilayers. This model is consistent with, and correlates with the effects of, most known pathological mutations found in FGFR‐TM.     

FGFR3 dimer stabilization due to a single amino acid pathogenic mutation

Mutations in the transmembrane (TM) domains of receptor tyrosine kinases (RTKs) have been implicated in the induction of pathological phenotypes. These mutations are believed to stabilize the RTK dimers, and thus promote unregulated signaling. However, the energetics behind the pathology induction has not been determined. An example of a TM domain pathogenic mutation is the Ala391-->Glu mutation in fibroblast growth factor receptor 3 (FGFR3), linked to Crouzon syndrome with acanthosis nigricans, as well as to bladder cancer. Here, we determine the free energy of dimerization of wild-type and mutant FGFR3 TM domain in lipid bilayers using F?rster resonance energy transfer, and we show that hydrogen bonding between Glu391 and the adjacent helix in the dimer is a feasible mechanism for dimer stabilization. The measured change in the free energy of dimerization due to the Ala391-->Glu pathogenic mutation is -1.3 kcal/mol, consistent with previous reports of hydrogen bond strengths in proteins. This is the first quantitative measurement of mutant RTK stabilization in a membrane environment. We show that this seemingly modest value can lead to a large increase in dimer fraction and thus profoundly affect RTK-mediated signal transduction.     

Molecular basis for dramatic changes in cannabinoid CB1 G protein‐coupled receptor activation upon single and double point mutations

There is considerable interest in determining the activation mechanism of G protein‐coupled receptors (GPCRs), one of the most important types of proteins for intercellular signaling. Recently, it was demonstrated for the cannabinoid CB1 GPCR, that a single mutation T210A could make CB1 completely inactive whereas T210I makes it essentially constitutively active. To obtain an understanding of this dramatic dependence of activity on mutation, we used first‐principles‐based methods to predict the ensemble of low‐energy seven‐helix conformations for the wild‐type (WT) and mutants (T210A and T210I). We find that the transmembrane (TM) helix packings depend markedly on these mutations, leading for T210A to both TM3+TM6 and TM2+TM6 salt‐bridge couplings in the cytoplasmic face that explains the inactivity of this mutant. In contrast T210I has no such couplings across the receptor explaining the ease in activating this mutant. WT has just the TM3+TM6 coupling, known to be broken upon GPCR activation. To test this hypothesis on activity, we predicted double mutants that would convert the inactive mutant to normal activity and then confirmed this experimentally. This CB1 activation mechanism, or one similar to it, is expected to play a role in other constitutively active GPCRs as well.     

His499 Regulates Dimerization and Prevents Oncogenic Activation by Asparagine Mutations of the Human Thrombopoietin Receptor

Ligand binding to the extracellular domain of the thrombopoietin receptor (TpoR) imparts a specific orientation on the transmembrane (TM) and intracellular domains of the receptors that is required for physiologic activation via receptor dimerization. To map the inactive and active dimeric orientations of the TM helices, we performed asparagine (Asn)-scanning mutagenesis of the TM domains of the murine and human TpoR. Substitution of Asn at only one position (S505N) activated the human receptor, whereas Asn substitutions at several positions activated the murine receptor. Second site mutational studies indicate that His⁴⁹⁹ near the N terminus of the TM domain is responsible for protecting the human receptor from activation by Asn mutations. Structural studies reveal that the sequence preceding His⁴⁹⁹ is helical in the murine receptor but non-helical in peptides corresponding to the TM domain of the inactive human receptor. The activating S505N mutation and the small molecule agonist eltrombopag both induce helix in this region of the TM domain and are associated with dimerization and activation of the human receptor. Thus, His⁴⁹⁹ regulates the activation of human TpoR and provides additional protection against activating mutations, such as oncogenic Asn mutations in the TM domain.     

Sodium dodecyl sulfate-polyacrylamide gel electrophoresis and forster resonance energy transfer suggest weak interactions between fibroblast growth factor receptor 3 (FGFR3) transmembrane domains in the absence of extracellular domains and ligands

Lateral dimerization of membrane proteins has evolved as a means of signal transduction across the plasma membrane for all receptor tyrosine kinases (RTKs). The transmembrane (TM) domains of RTKs are proposed to play an important role in the dimerization process. We have investigated whether the TM domains of one RTK, fibroblast growth factor receptor 3 (FGFR3), dimerize in lipid vesicles in the absence of the extracellular domains and ligands. We have performed sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) with peptides produced via solid-phase peptide synthesis that correspond to the TM domain of FGFR3. We have carried out Forster resonance energy transfer (FRET) measurements using two donor-acceptor pairs, fluorescein/rhodamine and Cy3/Cy5, as a function of peptide concentration and donor-to-acceptor mole ratios. Our results suggest that FGFR3 TM domains form sequence-specific dimers in lipid bilayers. However, the dimerization propensity of FGFR3 TM domain is much weaker than the dimerization propensity of glycophorin A (GpA), the well-characterized "membrane dimer standard". We discuss our findings in the context of cell signaling across the plasma membrane and diseases or disorders that occur due to single amino acid mutations in the TM domain of FGFR3.     

Receptor tyrosine kinase transmembrane domains: Function,dimer structure and dimerization energetics

The transmembrane (TM) domains of receptor tyrosine kinases (RTKs) play an active role in signaling. They contribute to the stability of full-length receptor dimers and to maintaining a signaling-competent dimeric receptor conformation. In an exciting new development, two structures of RTK TM domains have been solved, a break-through achievement in the field. Here we review these structures, and we discuss recent studies of RTK TM domain dimerization energetics, possible synergies between domains, and the effects of pathogenic RTK TM mutations on structure and dimerization.Key words: transmembrane domain, dimerization thermodynamics, receptor tyrosine kinases, pathogenic mutations, dimer structure     

Receptor tyrosine kinase transmembrane domains

The transmembrane (TM) domains of receptor tyrosine kinases (RTKs) play an active role in signaling. They contribute to the stability of full-length receptor dimers and to maintaining a signaling-competent dimeric receptor conformation. In an exciting new development, two structures of RTK TM domains have been solved, a break-through achievement in the field. Here we review these structures, and we discuss recent studies of RTK TM domain dimerization energetics, possible synergies between domains, and the effects of pathogenic RTK TM mutations on structure and dimerization.     

Random mutagenesis of the M3 muscarinic acetylcholine receptor expressed in yeast: identification of second-site mutations that restore function to a coupling-deficient mutant M3 receptor

The M(3) muscarinic receptor is a prototypical member of the class A family of G protein-coupled receptors (GPCRs). To gain insight into the structural mechanisms governing agonist-mediated M(3) receptor activation, we recently developed a genetically modified yeast strain (Saccharomyces cerevisiae) which allows the efficient screening of large libraries of mutant M(3) receptors to identify mutant receptors with altered/novel functional properties. Class A GPCRs contain a highly conserved Asp residue located in transmembrane domain II (TM II; corresponding to Asp-113 in the rat M(3) muscarinic receptor) which is of fundamental importance for receptor activation. As observed previously with other GPCRs analyzed in mammalian expression systems, the D113N point mutation abolished agonist-induced receptor/G protein coupling in yeast. We then subjected the D113N mutant M(3) receptor to PCR-based random mutagenesis followed by a yeast genetic screen to recover point mutations that can restore G protein coupling to the D113N mutant receptor. A large scale screening effort led to the identification of three such second-site suppressor mutations, R165W, R165M, and Y250D. When expressed in the wild-type receptor background, these three point mutations did not lead to an increase in basal activity and reduced the efficiency of receptor/G protein coupling. Similar results were obtained when the various mutant receptors were expressed and analyzed in transfected mammalian cells (COS-7 cells). Interestingly, like Asp-113, Arg-165 and Tyr-250, which are located at the cytoplasmic ends of TM III and TM V, respectively, are also highly conserved among class A GPCRs. Our data suggest a conformational link between the highly conserved Asp-113, Arg-165, and Tyr-250 residues which is critical for receptor activation.     

Random mutagenesis of the M3 muscarinic acetylcholine receptor expressed in yeast. Identification of point mutations that "silence" a constitutively active mutant M3 receptor and greatly impair receptor/G protein coupling

The M3 muscarinic receptor is a prototypical member of the class I family of G protein-coupled receptors (GPCRs). To facilitate studies on the structural mechanisms governing M3 receptor activation, we generated an M3 receptor-expressing yeast strain (Saccharomyces cerevisiae) that requires agonist-dependent M3 receptor activation for cell growth. By using receptor random mutagenesis followed by a genetic screen in yeast, we initially identified a point mutation at the cytoplasmic end of transmembrane domain (TM) VI (Q490L) that led to robust agonist-independent M3 receptor signaling in both yeast and mammalian cells. To explore further the molecular mechanisms by which point mutations can render GPCRs constitutively active, we subjected a region of the Q490L mutant M3 receptor that included TM V-VII to random mutagenesis. We then applied a yeast genetic screen to identify second-site mutations that could suppress the activating effects of the Q490L mutation and restore wild-type receptor-like function to the Q490L mutant receptor. This analysis led to the identification of 12 point mutations that allowed the Q490L mutant receptor to function in a fashion similar to the wild-type receptor. These amino acid substitutions mapped to two distinct regions of the M3 receptor, the exofacial segments of TM V and VI and the cytoplasmic ends of TM V-VII. Strikingly, in the absence of the activating Q490L mutation, all recovered point mutations severely reduced the efficiency of receptor/G protein coupling, indicating that the targeted residues play important roles in receptor activation and/or receptor/G protein coupling. This strategy should be generally applicable to identify sites in GPCRs that are critically involved in receptor function.     

Role of receptor tyrosine kinase transmembrane domains in cell signaling and human pathologies

Receptor tyrosine kinases (RTKs) conduct biochemical signals via lateral dimerization in the plasma membrane, and their transmembrane (TM) domains play an important role in the dimerization process. Here we present two models of RTK-mediated signaling, and we discuss the role of the TM domains within the framework of these two models. We summarize findings of single-amino acid mutations in RTK TM domains that induce unregulated signaling and, as a consequence, pathological phenotypes. We review the current knowledge of pathology induction mechanisms due to these mutations, focusing on the structural and thermodynamic basis of pathogenic dimer stabilization.     

High-Throughput Selection of Transmembrane Sequences That Enhance Receptor Tyrosine Kinase Activation

Dimerization is a critical requirement for the activation of the intracellular kinase domains of receptor tyrosine kinases (RTKs). The single transmembrane (TM) helices of RTKs contribute to dimerization, but the details are not well understood. Work with TM helices in various model systems has revealed a small number of specific dimerization sequence motifs, and it has been suggested that RTK dimerization is modulated by such motifs. Yet questions remain about the universality of these sequence motifs for RTK dimerization and about how TM domain dimerization in model systems relates to RTK activation in mammalian membranes. To investigate these questions, we designed a 3888-member combinatorial peptide library based on the TM domain of Neu (ErbB2) as a model RTK. The library contains many closely related, Neu-like sequences, including thousands of sequences with known dimerization motifs. We used an SDS-PAGE-based screen to select peptides that dimerize better than the native Neu sequence, and we assayed the activation of chimeric Neu receptors in mammalian cells with TM sequences selected in the screen. Despite the very high abundance of known dimerization motifs in the library, only a very few dimerizing sequences were identified by SDS-PAGE. About half of those sequences activated the Neu kinase significantly more than did the wild-type TM sequence. This work furthers our knowledge about the requirements for membrane protein interactions and the requirements for RTK activation in cells.     

A region in the seven-transmembrane domain of the human Ca2+ receptor critical for response to Ca2+

Of 12 naturally occurring, activating mutations in the seven-transmembrane (7TM) domain of the human Ca2+ receptor (CaR) identified previously in subjects with autosomal dominant hypocalcemia (ADH), five appear at the junction of TM helices 6 and 7 between residue Ile819 and Glu837. After identifying a sixth activating mutation in this region, V836L, in an ADH patient, we studied the remaining residues in this region to determine whether they are potential sites for activating mutations. Alanine-scanning mutagenesis revealed five additional residues in this region that when substituted by alanine led to CaR activation. We also found that, whereas E837A did not activate the receptor, E837D and E837K mutations did. Thus, region Ile819-Glu837 of the 7TM domain represents a "hot spot" for naturally occurring, activating mutations of the receptor, and most of the residues in this region apparently maintain the 7TM domain in its inactive configuration. Unique among the residues in this region, Pro823, which is highly conserved in family 3 of the G protein-coupled receptors, when mutated to either alanine or glycine, despite good expression severely impaired CaR activation by Ca2+. Both the P823A mutation and NPS 2143, a negative allosteric modulator that acts on the 7TM through a critical interaction with Glu837, blocked activation of the CaR by various ADH mutations. These results suggest that the 7TM domain region Ile819-Glu837 plays a key role in CaR activation by Ca2+. The implications of our finding that NPS 2143 corrects the molecular defect of ADH mutations for treatment of this disease are also discussed.     

Characterization of the full‐length human Grb7 protein and a phosphorylation representative mutant

It is well known the dimerization state of receptor tyrosine kinases (RTKs), in conjunction with binding partners such as the growth factor receptor bound protein 7 (Grb7) protein, plays an important role in cell signaling regulation. Previously, we proposed, downstream of RTKs, that the phosphorylation state of Grb7SH2 domain tyrosine residues could control Grb7 dimerization, and dimerization may be an important regulatory step in Grb7 binding to RTKs. In this manner, additional dimerization‐dependent regulation could occur downstream of the membrane‐bound kinase in RTK‐mediated signaling pathways. Extrapolation to the full‐length (FL) Grb7 protein, and the ability to test this hypothesis further, has been hampered by the availability of large quantities of pure and stable FL protein. Here, we report the biophysical characterization of the FL Grb7 protein and also a mutant representing a tyrosine‐phosphorylated Grb7 protein form. Through size exclusion chromatography and analytical ultracentrifugation, we show the phosphorylated‐tyrosine‐mimic Y492E‐FL‐Grb7 protein (Y492E‐FL‐Grb7) is essentially monomeric at expected physiological concentrations. It has been shown previously the wild‐type FL Grb7(WT‐FLGrb7) protein is dimeric with a dissociation constant (Kd) of approximately 11μM. Our studies here measure a FL protein dimerization K_d of WT‐FL‐Grb7 within one order of magnitude at approximately 1μM. The approximate size and shape of the WT‐FL‐Grb7 in comparison the tyrosine‐phosphorylation mimic Y492E‐FL‐Grb7 protein was determined by dynamic light scattering methods. In vitro phosphorylation of the Grb7SH2 domain indicates only one of the available tyrosine residues is phosphorylated, suggesting the same phosphorylation pattern could be relevant in the FL protein. The biophysical characterization studies in total are interpreted with a view towards understanding the functionally active Grb7 protein conformation.     

Nucleotide binding, ATP hydrolysis, and mutation of the catalytic carboxylates of human P-glycoprotein cause distinct conformational changes in the transmembrane segments

P-Glycoprotein (P-gp, ABCB1) transports a variety of structurally unrelated cytotoxic compounds out of the cell. Each homologous half of P-gp has a transmembrane (TM) domain containing six TM segments and a nucleotide-binding domain (NBD) and is joined by a linker region. It has been postulated that binding of two ATP molecules at the NBD interface to form a "nucleotide sandwich" induces drug efflux by altering packing of the TM segments that make up the drug-binding pocket. To test if ATP binding alone could alter packing of the TM segments, we introduced catalytic carboxylate mutations (E556Q in NBD1 and E1201Q in NBD2) into double-cysteine mutants that exhibited ATP-dependent cross-linking so that the mutants could bind but not hydrolyze ATP. It was found that ATP binding alone could alter disulfide cross-linking between the TM segments. For example, ATP inhibited cross-linking of mutant L339C(TM6)/V982C(TM12)/E556Q(NBD1)/E1201Q(NBD2) but promoted cross-linking of mutant F343C(TM6)/V982C(TM12)/E556Q(NBD1)/E1201Q(NBD2). Cross-linking of some mutants, however, appeared to require ATP hydrolysis as introduction of the catalytic carboxylate mutations into mutant L332C(TM6)/L975C(TM12) inhibited ATP-dependent cross-linking. Cross-linking between cysteines in the TM segments also could be altered via introduction of a single catalytic carboxylate mutation into mutant L332C(TM6)/L975C(TM12) or by using the nonhydrolyzable ATP analogue, AMP.PNP. The results show that the TM segments are quite sensitive to changes within the ATP-binding sites because different conformations could be detected in the presence of ATP, AMP.PNP, during ATP hydrolysis or through mutation of the catalytic carboxylates.     

Soluble Epoxide Hydrolase Dimerization Is Required for Hydrolase Activity

Soluble epoxide hydrolase (sEH) plays a key role in the metabolic conversion of the protective eicosanoid 14,15-epoxyeicosatrienoic acid to 14,15-dihydroxyeicosatrienoic acid. Accordingly, inhibition of sEH hydrolase activity has been shown to be beneficial in multiple models of cardiovascular diseases, thus identifying sEH as a valuable therapeutic target. Recently, a common human polymorphism (R287Q) was identified that reduces sEH hydrolase activity and is localized to the dimerization interface of the protein, suggesting a relationship between sEH dimerization and activity. To directly test the hypothesis that dimerization is essential for the proper function of sEH, we generated mutations within the sEH protein that would either disrupt or stabilize dimerization. We quantified the dimerization state of each mutant using a split firefly luciferase protein fragment-assisted complementation system. The hydrolase activity of each mutant was determined using a fluorescence-based substrate conversion assay. We found that mutations that disrupted dimerization also eliminated hydrolase enzymatic activity. In contrast, a mutation that stabilized dimerization restored hydrolase activity. Finally, we investigated the kinetics of sEH dimerization and found that the human R287Q polymorphism was metastable and capable of swapping dimer partners faster than the WT enzyme. These results indicate that dimerization is required for sEH hydrolase activity. Disrupting sEH dimerization may therefore serve as a novel therapeutic strategy for reducing sEH hydrolase activity.     

Dominant negative dimerization of a mutant homeodomain protein in Axenfeld-Rieger syndrome

Axenfeld-Rieger syndrome is an autosomal-dominant disorder caused by mutations in the PITX2 homeodomain protein. We have studied the mechanism underlying the dominant negative K88E mutation, which occurs at position 50 of the homeodomain. By using yeast two-hybrid and in vitro pulldown assays, we have documented that PITX2a can form homodimers in the absence of DNA. Moreover, the K88E mutant had even stronger dimerization ability, primarily due to interactions involving the C-terminal region. Dimerization allowed cooperative binding of wild-type (WT) PITX2a to DNA containing tandem bicoid sites in a head-to-tail orientation (Hill coefficient, 1.73). In contrast, the WT-K88E heterodimer bound the tandem sites with greatly reduced cooperativity and decreased transactivation activity. To further explore the role of position 50 in PITX2a dimerization, we introduced a charge-conservative mutation of lysine to arginine (K88R). The K88R protein had greatly reduced binding to a TAATCC element and did not specifically bind any other TAATNN motif. Like K88E, K88R formed relatively stronger dimers with WT. As predicted by our model, the K88R protein acted in a dominant negative manner to suppress WT PITX2a activity. These results suggest that the position 50 residue in the PITX2 homeodomain plays an important role in both DNA binding and dimerization activities.     

Effects of the oncogenic V(664)E mutation on membrane insertion, structure, and sequence-dependent interactions of the Neu transmembrane domain in micelles and model membranes: an integrated biophysical and simulation study

Receptor tyrosine kinases bind ligands such as cytokines, hormones, and growth factors and regulate key cellular processes, including cell division. They are also implicated in the development of many types of cancer. One such example is the Neu receptor tyrosine kinase found in rats (homologous to the human ErbB2 protein), which can undergo a valine to glutamic acid (V(664)E) mutation at the center of its α-helical transmembrane domain. This substitution results in receptor activation and oncogenesis. The molecular basis of this dramatic change in behavior upon introduction of the V(664)E mutation has been difficult to pin down, with conflicting results reported in the literature. Here we report the first quantitative, thermodynamic analysis of dimerization and biophysical characterization of the rat Neu transmembrane domain and several mutants in a range of chemical environments. These data have allowed us to identify the effects of the V(664)E mutation in the isolated TM domain with respect to protein-protein and protein-lipid interactions, membrane insertion, and secondary structure. We also report the results from a 100 ns atomistic molecular dynamics simulation of the Neu transmembrane domain in a model membrane bilayer (dipalmitoylphosphatidylcholine). The results from simulation and experiment are in close agreement and suggest that, in the model systems investigated, the V(664)E mutation leads to a weakening of the TM dimer and a change in sequence-dependent interactions. These results are contrary to recent results obtained in mammalian membranes, and the implications of this are discussed.     

Mutation D816V alters the internal structure and dynamics of c-KIT receptor cytoplasmic region: implications for dimerization and activation mechanisms

The type III receptor tyrosine kinase (RTK) KIT plays a crucial role in the transmission of cellular signals through phosphorylation events that are associated with a switching of the protein conformation between inactive and active states. D816V KIT mutation is associated with various pathologies including mastocytosis and cancers. D816V-mutated KIT is constitutively active, and resistant to treatment with the anti-cancer drug Imatinib. To elucidate the activating molecular mechanism of this mutation, we applied a multi-approach procedure combining molecular dynamics (MD) simulations, normal modes analysis (NMA) and binding site prediction. Multiple 50-ns MD simulations of wild-type KIT and its mutant D816V were recorded using the inactive auto-inhibited structure of the protein, characteristic of type III RTKs. Computed free energy differences enabled us to quantify the impact of D816V on protein stability in the inactive state. We evidenced a local structural alteration of the activation loop (A-loop) upon mutation, and a long-range structural re-organization of the juxta-membrane region (JMR) followed by a weakening of the interaction network with the kinase domain. A thorough normal mode analysis of several MD conformations led to a plausible molecular rationale to propose that JMR is able to depart its auto-inhibitory position more easily in the mutant than in wild-type KIT and is thus able to promote kinase mutant dimerization without the need for extra-cellular ligand binding. Pocket detection at the surface of NMA-displaced conformations finally revealed that detachment of JMR from the kinase domain in the mutant was sufficient to open an access to the catalytic and substrate binding sites.