Role of dimerization efficiency of transmembrane domains in activation of fibroblast growth factor receptor 3

Diversity of Structural, Dynamic, and Environmental Effects Explain a Distinctive Functional Role of Transmembrane Domains in the Insulin Receptor Subfamily

Human InsR, IGF1R, and IRR receptor tyrosine kinases (RTK) of the insulin receptor subfamily play an important role in signaling pathways for a wide range of physiological processes and are directly associated with many pathologies, including neurodegenerative diseases. The disulfide-linked dimeric structure of these receptors is unique among RTKs. Sharing high sequence and structure homology, the receptors differ dramatically in their localization, expression, and functions. In this work, using high-resolution NMR spectroscopy supported by atomistic computer modeling, conformational variability of the transmembrane domains and their interactions with surrounding lipids were found to differ significantly between representatives of the subfamily. Therefore, we suggest that the heterogeneous and highly dynamic membrane environment should be taken into account in the observed diversity of the structural/dynamic organization and mechanisms of activation of InsR, IGF1R, and IRR receptors. This membrane-mediated control of receptor signaling offers an attractive prospect for the development of new targeted therapies for diseases associated with dysfunction of insulin subfamily receptors.

Role of the N-Terminal Transmembrane Helix Contacts in the Activation of FGFR3

Fibroblast growth factor receptor 3 (FGFR3) is a member of receptor tyrosine kinases, which is involved in skeletal cell growth, differentiation, and migration. FGFR3 transduces biochemical signals from the extracellular ligand-binding domain to the intracellular kinase domain through the conformational changes of the transmembrane (TM) helix dimer. Here, we apply generalized replica exchange with solute tempering method to wild type (WT) and G380R mutant (G380R) of FGFR3. The dimer interface in G380R is different from WT and the simulation results are in good agreement with the solid-state nuclear magnetic resonance (NMR) spectroscopy. TM helices in G380R are extended more than WT, and thereby, G375 in G380R contacts near the N-termini of the TM helix dimer. Considering that both G380R and G375C show the constitutive activation, the formation of the N-terminal contacts of the TM helices can be generally important for the activation mechanism. © 2019 Wiley Periodicals, Inc.

Resolving the conformational dynamics of ErbB growth factor receptor dimers

The combinatorial dimerization of the ErbB growth factor receptors (ErbB¹- ErbB⁴) are critical for their function. Here, we have characterized the conformational dynamics of ErbB transmembrane homo-dimers and hetero-dimers by using a coarse-grain simulation framework. All dimers, except ErbB^4-4 and ErbB^1-4, exhibit at least two conformations. The reported NMR structures correspond to one of these conformations, representing the N-terminal active state in ErbB^1-1 (RH2), ErbB^2-2 (RH1) and ErbB^4-4 (RH) homo-dimers and the LH dimer in ErbB^3-3 homo-dimer, validating the computational approach. Further, we predict a right-handed ErbB^3-3 dimer conformer that warrants experimental testing. The five hetero-dimers that have not yet been experimentally resolved display prominent right-handed dimers associating by the SmXXXSm motif. Our results provide insights into the constitutive signaling of ErbB⁴ after cleavage of the extracellular region. The presence of the inactive-like dimer conformers leading to symmetric kinase domains gives clues on the autoinhibition of the receptor dimers. The dimer states characterized here represent an important step towards understanding the combinatorial cross associations in the ErbB family.

Atomistic mechanism of the constitutive activation of PDGFRA via its transmembrane domain

Single-point mutations in the transmembrane (TM) region of receptor tyrosine kinases (RTKs) can lead to abnormal ligand-independent activation. We use a combination of computational modeling, NMR spectroscopy and cell experiments to analyze in detail the mechanism of how TM domains contribute to the activation of wild-type (WT) PDGFRA and its oncogenic V536E mutant. Using a computational framework, we scan all positions in PDGFRA TM helix for identification of potential functional mutations for the WT and the mutant and reveal the relationship between the receptor activity and TM dimerization via different interfaces. This strategy also allows us design a novel activating mutation in the WT (I537D) and a compensatory mutation in the V536E background eliminating its constitutive activity (S541G). We show both computationally and experimentally that single-point mutations in the TM region reshape the TM dimer ensemble and delineate the structural and dynamic determinants of spontaneous activation of PDGFRA via its TM domain. Our atomistic picture of the coupling between TM dimerization and PDGFRA activation corroborates the data obtained for other RTKs and provides a foundation for developing novel modulators of the pathological activity of PDGFRA.

Death Receptor 5 Activation Is Energetically Coupled to Opening of the Transmembrane Domain Dimer

The precise mechanism by which binding of tumor necrosis factor ligands to the extracellular domain of their corresponding receptors transmits signals across the plasma membrane has remained elusive. Recent studies have proposed that activation of several tumor necrosis factor receptors, including Death Receptor 5, involves a scissorlike opening of the disulfide-linked transmembrane (TM) dimer. Using time-resolved fluorescence resonance energy transfer, we provide, to our knowledge, the first direct biophysical evidence that Death Receptor 5 TM-dimers open in response to ligand binding. Then, to probe the importance of the closed-to-open TM domain transition in the overall energetics of receptor activation, we designed point-mutants (alanine to phenylalanine) in the predicted, tightly packed TM domain dimer interface. We hypothesized that the bulky residues should destabilize the closed conformation and eliminate the ∼3 kcal/mol energy barrier to TM domain opening and the ∼2 kcal/mol energy difference between the closed and open states, thus oversensitizing the receptor. To test this, we used all-atom molecular dynamics simulations of the isolated TM domain in explicit lipid bilayers coupled to thermodynamic potential of mean force calculations. We showed that single point mutants at the interface altered the energy landscape as predicted, but were not enough to completely eliminate the barrier to opening. However, the computational model did predict that a double mutation at i, i+4 positions at the center of the TM domain dimer eliminates the barrier and stabilizes the open conformation relative to the closed. We tested these mutants in cells with time-resolved fluorescence resonance energy transfer and death assays, and show remarkable agreement with the calculations. The single mutants had a small effect on TM domain separation and cell death, whereas the double mutant significantly increased the TM domain separation and more than doubled the sensitivity of cells to ligand stimulation.

Identification and in silico characterization of p.G380R substitution in FGFR3, associated with achondroplasia in a non-consanguineous Pakistani family

Background

The dimerization efficiency of FGFR3 transmembrane domain plays a critical role in the formation of a normal skeleton through the negative regulation of bone development. Recently, gain-of-function mutations in the transmembrane domain of FGFR3 has been described associated with an aberrant negative regulation, leading to the development of achondroplasia-group disorders, including achondroplasia (ACH), hypochondroplasia (HCH) and thanatophoric dysplasia (TD). Here, we describe a non-consanguineous Pakistani family with achondroplasia to explain hereditary basis of the disease.

Methods

PCR-based linkage analysis using microsatellite markers was employed to localize the disease gene. Gene specific intronic primers were used to amplify the genomic DNA from all affected as well as phenotypically healthy individuals. Amplified PCR products were then subjected to Sanger sequencing and RFLP analysis to identify a potentially pathogenic mutation. The impact of identified mutation on FGFR3 protein’s structure and stability was highlighted through different bioinformatics tools.

Results

Genetic screening of the family revealed a previously reported heterozygous c.1138 G > A (p.G380R) mutation in the coding exon 8 of FGFR3 gene. Identified genetic variation was confirmed in all affected individuals while healthy individuals and controls were found genotypically normal. The results were further validated by RFLP analysis as c.1138 G > A substitution generates a unique recognition site for SfcI endonuclease. Following SfcI digestion, the electrophoretic pattern of three bands/DNA fragments for each patient is indicative of heterozygous status of the disease allele. In silico studies of the mutant FGFR3 protein predicted to adversely affect the stability of FGFR3 protein.

Conclusions

Mutation in the transmembrane domain may adversely affect the dimerization efficiency and overall stability of the FGFR3, leading to a constitutively active protein. As a result, an uncontrolled intracellular signaling or negative bone growth regulation leads to achondroplasia. Our findings support the fact that p.G380R is a common mutation among diverse population of the world and like other countries, can be used as a molecular diagnosis marker for achondroplasia in Pakistan.

Electronic supplementary material

The online version of this article (doi:10.1186/s13000-017-0642-3) contains supplementary material, which is available to authorized users.

Helix-helix interactions in membrane domains of bitopic proteins: Specificity and role of lipid environment

Interaction between transmembrane helices often determines biological activity of membrane proteins. Bitopic proteins, a broad subclass of membrane proteins, form dimers containing two membrane-spanning helices. Some aspects of their structure-function relationship cannot be fully understood without considering the protein-lipid interaction, which can determine the protein conformational ensemble. Experimental and computer modeling data concerning transmembrane parts of bitopic proteins are reviewed in the present paper. They highlight the importance of lipid-protein interactions and resolve certain paradoxes in the behavior of such proteins. Besides, some properties of membrane organization provided a clue to understanding of allosteric interactions between distant parts of proteins. Interactions of these kinds appear to underlie a signaling mechanism, which could be widely employed in the functioning of many membrane proteins. Treatment of membrane proteins as parts of integrated fine-tuned proteolipid system promises new insights into biological function mechanisms and approaches to drug design. This article is part of a Special Issue entitled: Lipid order/lipid defects and lipid-control of protein activity edited by Dirk Schneider.

Conformational Changes in the Epidermal Growth Factor Receptor: Role of the Transmembrane Domain Investigated by Coarse-Grained MetaDynamics Free Energy Calculations

The epidermal growth factor receptor (EGFR) is a dimeric membrane protein that regulates key aspects of cellular function. Activation of the EGFR is linked to changes in the conformation of the transmembrane (TM) domain, brought about by changes in interactions of the TM helices of the membrane lipid bilayer. Using an advanced computational approach that combines Coarse-Grained molecular dynamics and well-tempered MetaDynamics (CG-MetaD), we characterize the large-scale motions of the TM helices, simulating multiple association and dissociation events between the helices in membrane, thus leading to a free energy landscape of the dimerization process. The lowest energy state of the TM domain is a right-handed dimer structure in which the TM helices interact through the N-terminal small-X₃-small sequence motif. In addition to this state, which is thought to correspond to the active form of the receptor, we have identified further low-energy states that allow us to integrate with a high level of detail a range of previous experimental observations. These conformations may lead to the active state via two possible activation pathways, which involve pivoting and rotational motions of the helices, respectively. Molecular dynamics also reveals correlation between the conformational changes of the TM domains and of the intracellular juxtamembrane domains, paving the way for a comprehensive understanding of EGFR signaling at the cell membrane.

Mutational mapping of the transmembrane binding site of the G-protein coupled receptor TGR5 and binding mode prediction of TGR5 agonists

TGR5 (Gpbar-1, M-Bar) is a class A G-protein coupled bile acid-sensing receptor predominately expressed in brain, liver and gastrointestinal tract, and a promising drug target for the treatment of metabolic disorders. Due to the lack of a crystal structure of TGR5, the development of TGR5 agonists has been guided by ligand-based approaches so far. Three binding mode models of bile acid derivatives have been presented recently. However, they differ from one another in terms of overall orientation or with respect to the location and interactions of the cholane scaffold, or cannot explain all results from mutagenesis experiments. Here, we present an extended binding mode model based on an iterative and integrated computational and biological approach. An alignment of 68 TGR5 agonists based on this binding mode leads to a significant and good structure-based 3D QSAR model, which constitutes the most comprehensive structure-based 3D-QSAR study of TGR5 agonists undertaken so far and suggests that the binding mode model is a close representation of the "true" binding mode. The binding mode model is further substantiated in that effects predicted for eight mutations in the binding site agree with experimental analyses on the impact of these TGR5 variants on receptor activity. In the binding mode, the hydrophobic cholane scaffold of taurolithocholate orients towards the interior of the orthosteric binding site such that rings A and B are in contact with TM5 and TM6, the taurine side chain orients towards the extracellular opening of the binding site and forms a salt bridge with R79(EL1), and the 3-hydroxyl group forms hydrogen bonds with E169(5.44) and Y240(6.51). The binding mode thus differs in important aspects from the ones recently presented. These results are highly relevant for the development of novel, more potent agonists of TGR5 and should be a valuable starting point for the development of TGR5 antagonists, which could show antiproliferative effects in tumor cells.

Functional competition within a membrane: Lipid recognition vs. transmembrane helix oligomerization

Binding of specific lipids to large, polytopic membrane proteins is well described, and it is clear that such lipids are crucial for protein stability and activity. In contrast, binding of defined lipid species to individual transmembrane helices and regulation of transmembrane helix monomer-oligomer equilibria by binding of distinct lipids is a concept, which has emerged only lately. Lipids bind to single-span membrane proteins, both in the juxta-membrane region as well as in the hydrophobic membrane core. While some interactions counteract transmembrane helix oligomerization, in other cases lipid binding appears to enhance oligomerization. As reversible oligomerization is involved in activation of many membrane proteins, binding of defined lipids to single-span transmembrane proteins might be a mechanism to regulate and/or fine-tune the protein activity. But how could lipid binding trigger the activity of a protein? How can binding of a single lipid molecule to a transmembrane helix affect the structure of a transmembrane helix oligomer, and consequently its signaling state? These questions are discussed in the present article based on recent results obtained with simple, single-span transmembrane proteins. This article is part of a Special Issue entitled: Lipid-protein interactions.

Dimerization of the EphA1 receptor tyrosine kinase transmembrane domain: Insights into the mechanism of receptor activation

EphA1 is a receptor tyrosine kinase (RTK) that plays a key role in developmental processes, including guidance of the migration of axons and cells in the nervous system. EphA1, in common with other RTKs, contains an N-terminal extracellular domain, a single transmembrane (TM) α-helix, and a C-terminal intracellular kinase domain. The TM helix forms a dimer, as seen in recent NMR studies. We have modeled the EphA1 TM dimer using a multiscale approach combining coarse-grain (CG) and atomistic molecular dynamics (MD) simulations. The one-dimensional potential of mean force (PMF) for this system, based on interhelix separation, has been calculated using CG MD simulations. This provides a view of the free energy landscape for helix–helix interactions of the TM dimer in a lipid bilayer. The resulting PMF profiles suggest two states, consistent with a rotation-coupled activation mechanism. The more stable state corresponds to a right-handed helix dimer interacting via an N-terminal glycine zipper motif, consistent with a recent NMR structure (2K1K). A second metastable state corresponds to a structure in which the glycine zipper motif is not involved. Analysis of unrestrained CG MD simulations based on representative models from the PMF calculations or on the NMR structure reveals possible pathways of interconversion between these two states, involving helix rotations about their long axes. This suggests that the interaction of TM helices in EphA1 dimers may be intrinsically dynamic. This provides a potential mechanism for signaling whereby extracellular events drive a shift in the repopulation of the underlying TM helix dimer energy landscape.

Point mutations in dimerization motifs of the transmembrane domain stabilize active or inactive state of the EphA2 receptor tyrosine kinase

The EphA2 receptor tyrosine kinase plays a central role in the regulation of cell adhesion and guidance in many human tissues. The activation of EphA2 occurs after proper dimerization/oligomerization in the plasma membrane, which occurs with the participation of extracellular and cytoplasmic domains. Our study revealed that the isolated transmembrane domain (TMD) of EphA2 embedded into the lipid bicelle dimerized via the heptad repeat motif L(535)X3G(539)X2A(542)X3V(546)X2L(549) rather than through the alternative glycine zipper motif A(536)X3G(540)X3G(544) (typical for TMD dimerization in many proteins). To evaluate the significance of TMD interactions for full-length EphA2, we substituted key residues in the heptad repeat motif (HR variant: G539I, A542I, G553I) or in the glycine zipper motif (GZ variant: G540I, G544I) and expressed YFP-tagged EphA2 (WT, HR, and GZ variants) in HEK293T cells. Confocal microscopy revealed a similar distribution of all EphA2-YFP variants in cells. The expression of EphA2-YFP variants and their kinase activity (phosphorylation of Tyr(588) and/or Tyr(594)) and ephrin-A3 binding were analyzed with flow cytometry on a single cell basis. Activation of any EphA2 variant is found to occur even without ephrin stimulation when the EphA2 content in cells is sufficiently high. Ephrin-A3 binding is not affected in mutant variants. Mutations in the TMD have a significant effect on EphA2 activity. Both ligand-dependent and ligand-independent activities are enhanced for the HR variant and reduced for the GZ variant compared with the WT. These findings allow us to suggest TMD dimerization switching between the heptad repeat and glycine zipper motifs, corresponding to inactive and active receptor states, respectively, as a mechanism underlying EphA2 signal transduction.

Primary and secondary dimer interfaces of the fibroblast growth factor receptor 3 transmembrane domain: characterization via multiscale molecular dynamics simulations

Receptor tyrosine kinases are single-pass membrane proteins that form dimers within the membrane. The interactions of their transmembrane domains (TMDs) play a key role in dimerization and signaling. Fibroblast growth factor receptor 3 (FGFR3) is of interest as a G380R mutation in its TMD is the underlying cause of ~99% of the cases of achondroplasia, the most common form of human dwarfism. The structural consequences of this mutation remain uncertain: the mutation shifts the position of the TMD relative to the lipid bilayer but does not alter the association free energy. We have combined coarse-grained and all-atom molecular dynamics simulations to study the dimerization of wild-type, heterodimer, and mutant FGFR3 TMDs. The simulations reveal that the helices pack together in the dimer to form a flexible interface. The primary packing mode is mediated by a Gx3G motif. There is also a secondary dimer interface that is more highly populated in heterodimer and mutant configurations that may feature in the molecular mechanism of pathology. Both coarse-grained and atomistic simulations reveal a significant shift of the G380R mutant dimer TMD relative to the bilayer to allow interactions of the arginine side chain with lipid headgroup phosphates.

Structure of FGFR3 transmembrane domain dimer: implications for signaling and human pathologies

Fibroblast growth factor receptor 3 (FGFR3) transduces biochemical signals via lateral dimerization in the plasma membrane, and plays an important role in human development and disease. Eight different pathogenic mutations, implicated in cancers and growth disorders, have been identified in the FGFR3 transmembrane segment. Here, we describe the dimerization of the FGFR3 transmembrane domain in membrane-mimicking DPC/SDS (9/1) micelles. In the solved NMR structure, the two transmembrane helices pack into a symmetric left-handed dimer, with intermolecular stacking interactions occurring in the dimer central region. Some pathogenic mutations fall within the helix-helix interface, whereas others are located within a putative alternative interface. This implies that although the observed dimer structure is important for FGFR3 signaling, the mechanism of FGFR3-mediated transduction across the membrane is complex. We propose an FGFR3 signaling mechanism that is based on the solved structure, available structures of isolated soluble FGFR domains, and published biochemical and biophysical data.