Drug binding and resistance mechanism of KIT tyrosine kinase revealed by hydrogen/deuterium exchange FTICR mass spectrometry

Mutations of the receptor tyrosine kinase KIT are linked to certain cancers such as gastrointestinal stromal tumors (GISTs). Biophysical, biochemical, and structural studies have provided insight into the molecular basis of resistance to the KIT inhibitors, imatinib and sunitinib. Here, solution‐phase hydrogen/deuterium exchange (HDX) and direct binding mass spectrometry experiments provide a link between static structure models and the dynamic equilibrium of the multiple states of KIT, supporting that sunitinib targets the autoinhibited conformation of WT‐KIT. The D816H mutation shifts the KIT conformational equilibrium toward the activated state. The V560D mutant exhibits two low energy conformations: one is more flexible and resembles the D816H mutant shifted toward the activated conformation, and the other is less flexible and resembles the wild‐type KIT in the autoinhibited conformation. This result correlates with the V560D mutant exhibiting a sensitivity to sunitinib that is less than for WT KIT but greater than for KIT D816H. These findings support the elucidation of the resistance mechanism for the KIT mutants.     

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.     

Differential Effects of CSF-1R D802V and KIT D816V Homologous Mutations on Receptor Tertiary Structure and Allosteric Communication

The colony stimulating factor-1 receptor (CSF-1R) and the stem cell factor receptor KIT, type III receptor tyrosine kinases (RTKs), are important mediators of signal transduction. The normal functions of these receptors can be compromised by gain-of-function mutations associated with different physiopatological impacts. Whereas KIT D816V/H mutation is a well-characterized oncogenic event and principal cause of systemic mastocytosis, the homologous CSF-1R D802V has not been identified in human cancers. The KIT D816V oncogenic mutation triggers resistance to the RTK inhibitor Imatinib used as first line treatment against chronic myeloid leukemia and gastrointestinal tumors. CSF-1R is also sensitive to Imatinib and this sensitivity is altered by mutation D802V. Previous in silico characterization of the D816V mutation in KIT evidenced that the mutation caused a structure reorganization of the juxtamembrane region (JMR) and facilitated its departure from the kinase domain (KD). In this study, we showed that the equivalent CSF-1R D802V mutation does not promote such structural effects on the JMR despite of a reduction on some key H-bonds interactions controlling the JMR binding to the KD. In addition, this mutation disrupts the allosteric communication between two essential regulatory fragments of the receptors, the JMR and the A-loop. Nevertheless, the mutation-induced shift towards an active conformation observed in KIT D816V is not observed in CSF-1R D802V. The distinct impact of equivalent mutation in two homologous RTKs could be associated with the sequence difference between both receptors in the native form, particularly in the JMR region. A local mutation-induced perturbation on the A-loop structure observed in both receptors indicates the stabilization of an inactive non-inhibited form, which Imatinib cannot bind.     

Allosteric Communication across the Native and Mutated KIT Receptor Tyrosine Kinase

A fundamental goal in cellular signaling is to understand allosteric communication, the process by which signals originated at one site in a protein propagate dependably to affect remote functional sites. Here, we describe the allosteric regulation of the receptor tyrosine kinase KIT. Our analysis evidenced that communication routes established between the activation loop (A-loop) and the distant juxtamembrane region (JMR) in the native protein were disrupted by the oncogenic mutation D816V positioned in the A-loop. In silico mutagenesis provided a plausible way of restoring the protein communication detected in the native KIT by introducing a counter-balancing second mutation D792E. The communication patterns observed in the native and mutated KIT correlate perfectly with the structural and dynamical features of these proteins. Particularly, a long-distance effect of the D816V mutation manifested as an important structural re-organization of the JMR in the oncogenic mutant was completely vanished in the double mutant D816V/D792E. This detailed characterization of the allosteric communication in the different forms of KIT, native and mutants, was performed by using a modular network representation composed of communication pathways and independent dynamic segments. Such representation permits to enrich a purely mechanistic interaction-based model of protein communication by the introduction of concerted local atomic fluctuations. This method, validated on KIT receptor, may guide a rational modulation of the physiopathological activities of other receptor tyrosine kinases.     

Neoplasia driven by mutant c-KIT is mediated by intracellular, not plasma membrane, receptor signaling

Activating mutations in c-KIT are associated with gastrointestinal stromal tumors, mastocytosis, and acute myeloid leukemia. In attempting to establish a murine model of human KIT(D816V) (hKIT(D816V))-mediated leukemia, we uncovered an unexpected relationship between cellular transformation and intracellular trafficking. We found that transport of hKIT(D816V) protein was blocked at the endoplasmic reticulum in a species-specific fashion. We exploited these species-specific trafficking differences and a set of localization domain-tagged KIT mutants to explore the relationship between subcellular localization of mutant KIT and cellular transformation. The protein products of fully transforming KIT mutants localized to the Golgi apparatus and to a lesser extent the plasma membrane. Domain-tagged KIT(D816V) targeted to the Golgi apparatus remained constitutively active and transforming. Chemical inhibition of intracellular transport demonstrated that Golgi localization is sufficient, but plasma membrane localization is dispensable, for downstream signaling mediated by KIT mutation. When expressed in murine bone marrow, endoplasmic reticulum-localized hKIT(D816V) failed to induce disease in mice, while expression of either Golgi-localized HyKIT(D816V) or cytosol-localized, ectodomain-deleted KIT(D816V) uniformly caused fatal myeloproliferative diseases. Taken together, these data demonstrate that intracellular, non-plasma membrane receptor signaling is sufficient to drive neoplasia caused by mutant c-KIT and provide the first animal model of myelomonocytic neoplasia initiated by human KIT(D816V).     

Conformation analysis of a surface loop that controls active site access in the GH11 xylanase A from Bacillus subtilis

Xylanases (EC 3.2.1.8 endo-1,4-glycosyl hydrolase) catalyze the hydrolysis of xylan, an abundant hemicellulose of plant cell walls. Access to the catalytic site of GH11 xylanases is regulated by movement of a short β-hairpin, the so-called thumb region, which can adopt open or closed conformations. A crystallographic study has shown that the D11F/R122D mutant of the GH11 xylanase A from Bacillus subtilis (BsXA) displays a stable “open” conformation, and here we report a molecular dynamics simulation study comparing this mutant with the native enzyme over a range of temperatures. The mutant open conformation was stable at 300 and 328 K, however it showed a transition to the closed state at 338 K. Analysis of dihedral angles identified thumb region residues Y113 and T123 as key hinge points which determine the open-closed transition at 338 K. Although the D11F/R122D mutations result in a reduction in local inter-intramolecular hydrogen bonding, the global energies of the open and closed conformations in the native enzyme are equivalent, suggesting that the two conformations are equally accessible. These results indicate that the thumb region shows a broader degree of energetically permissible conformations which regulate the access to the active site region. The R122D mutation contributes to the stability of the open conformation, but is not essential for thumb dynamics, i.e., the wild type enzyme can also adapt to the open conformation.     

A shift of dynamic equilibrium between the KIT active and inactive states causes drug resistance

Tyrosine phosphorylation, a highly regulated post-translational modification, is carried out by the enzyme tyrosine kinase (TK). TKs are important mediators in signaling cascades, facilitating diverse biological processes in response to stimuli. TKs may acquire mutations leading to malignancy and are viable targets for anti-cancer drugs. Mast/stem cell growth factor receptor KIT is a TK involved in cell differentiation, whose dysregulation leads to various types of cancer, including gastrointestinal stromal tumors, leukemia, and melanoma. KIT can be targeted by a range of inhibitors that predominantly bind to the inactive state of the enzyme. A mutation Y823D in the activation loop of KIT is known to be responsible for the loss of sensitivity to some drugs in metastatic tumors. We used all-atom molecular dynamics simulations to study the impact of Y823D on the KIT conformation and dynamics and compared it to the effect of phosphorylation of Y823. We simulated in total 6.4 μs of wild-type, mutant and phosphorylated KIT in the active- and inactive-state conformations. We found that Y823D affects the protein dynamics differently: in the active state, the mutation increases the protein stability, whereas in the inactive state it induces local destabilization, thus shifting the dynamic equilibrium towards the active state, altering the communication between distant regulatory regions. The observed dynamics of the Y823D mutant is similar to the dynamics of KIT phosphorylated at position Y823, thus we hypothesize that this mutation mimics a constitutively active kinase, which is not responsive to inhibitors that bind its inactive conformation.     

Anti-insulin antibody structure and conformation. II. Molecular dynamics with explicit solvent.

Molecular dynamics at 300 K was used as a conformation searching tool to analyze a knowledge-based structure prediction of an anti-insulin antibody. Solvation effects were modeled by packing water molecules around the antigen binding loops. Some loops underwent backbone and side-chain conformational changes during the 95-ps equilibration, and most of these new, lower potential energy conformations were stable during the subsequent 200-ps simulation. Alterations to the model include changes in the intraloop, main-chain hydrogen bonding network of loop H3, and adjustments of Tyr and Lys side chains of H3 induced by hydrogen bonding to water molecules. The structures observed during molecular dynamics support the conclusion of the previous paper that hydrogen bonding will play the dominant role in antibody-insulin recognition. Determination of the structure of the antibody by x-ray crystallography is currently being pursued to provide an experimental test of these results. The simulation appears to improve the model, but longer simulations at higher temperatures should be performed.     

Comparative molecular dynamics simulation analysis of G20 and C92 mutations in c-di-GMP I riboswitch and the wild type with docked c-di-GMP ligand

Riboswitch, a bacterial regulatory RNA consists of an aptamer (specific ligand binding unit) and an expression platform (gene expression modulation unit), which act as a potential drug target as it regulates critical genes. Therefore, it is of interest to glean information on the binding of c-di-GMP ligand to mutated conserved G20 and C92 residues of cyclic diguanosine monophosphate I (c-di-GMP I) riboswitch using molecular dynamics simulation. The result shows that the binding energy of wild/native type riboswitch-ligand complex (3IRW) is lower than the mutant complexes suggesting that the binding affinity for c-di-GMP ligand decreases in case of mutant riboswitches. The hydrogen bonding interactions analysis also showed a high number of hydrogen bonds formation in the wild type riboswitch-ligand complex as compared to the mutant complexes illustrating stronger interaction of ligand to wild type riboswitch than the mutants. The simulation result shows that the mutations affected riboswitch-ligand interactions. The residues G14, G21, C46, A47, and U92 were identified as the key residues which contributed effectively to the binding of c-di-GMP I riboswitch with the natural ligand.     

Molecular modeling of Gly80 and Ser80 variants of human group IID phospholipase A2 and their receptor complexes: potential basis for weight loss in chronic obstructive pulmonary disease

Weight loss is a well known systemic manifestation of chronic obstructive pulmonary disease (COPD). A Gly80Ser mutation on human group IID secretory phospholipase A2 (sPLA₂) enhances expression of the cytokines that are responsible for weight loss. In this study, we seek to establish a structural correlation of wild type sPLA₂ and the Gly80Ser mutation with function. sPLA₂ with glycine and serine at the 80th positions and the M-type receptor were modelled. The enzymes were docked to the receptor and molecular dynamics was carried out to 70 ns. Structural analysis revealed the enzymes to comprise three helices (H1–H3), two short helices (SH1 and SH2), and five loops including a calcium binding loop (L1–L5), and to be stabilized by seven disulfide bonds. The overall backbone folds of the two models are very similar, with main chain RMSD of less than 1 Å. The active site within the substrate binding channel shows a catalytic triad of water–His67–Asp112, showing a hydrogen bonded network. Major structural differences between wild type and mutant enzymes were observed locally at the site of the mutation and in their global conformations. These differences include: (1) loop-L3 between H2 and H3, which bears residue Gly80 in the wild type, is in a closed conformation with respect to the channel opening, while in the mutant enzyme it adopts a relatively open conformation; (2) the mutant enzyme is less compact and has higher solvent accessible surface area; and (3) interfacial binding contact surface area is greater, and the quality of interactions with the receptor is better in the mutant enzyme as compared to the wild type. Therefore, the structural differences delineated in this study are potential biophysical factors that could determine the increased potency of the mutant enzyme with macrophage receptor for cytokine secreting function, resulting in exacerbation of cachexia in COPD.     

Theoretical investigation on the glycan‐binding specificity of Agrocybe cylindracea galectin using molecular modeling and molecular dynamics simulation studies

Galectins are β‐galactoside binding proteins which have the ability to serve as potent antitumor, cancer biomarker, and induce tumor cell apoptosis. Agrocybe cylindracea galectin (ACG) is a fungal galectin which specifically recognizes α(2,3)‐linked sialyllactose at the cell surface that plays extensive roles in the biological recognition processes. To investigate the change in glycan‐binding specificity upon mutations, single point and double point site‐directed in silico mutations are performed at the binding pocket of ACG. Molecular dynamics (MD) simulation studies are carried out for the wild‐type (ACG) and single point (ACG1) and double point (ACG2) mutated ACGs to investigate the dynamics of substituted mutants and their interactions with the receptor sialyllactose. Plausible binding modes are proposed for galectin–sialylglycan complexes based on the analysis of hydrogen bonding interactions, total pair‐wise interaction energy between the interacting binding site residues and sialyllactose and binding free energy of the complexes using molecular mechanics–Poisson–Boltzmann surface area. Our result shows that high contribution to the binding in different modes is due to the direct and water‐mediated hydrogen bonds. The binding specificity of double point mutant Y59R/N140Q of ACG2 is found to be high, and it has 26 direct and water‐mediated hydrogen bonds with a relatively low‐binding free energy of −47.52 ± 5.2 kcal/mol. We also observe that the substituted mutant Arg59 is crucial for glycan‐binding and for the preference of α(2,3)‐linked sialyllactose at the binding pocket of ACG2 galectin. When compared with the wild‐type and single point mutant, the double point mutant exhibits enhanced affinity towards α(2,3)‐linked sialyllactose, which can be effectively used as a model for biological cell marker in cancer therapeutics. Copyright © 2015 John Wiley & Sons, Ltd.     

Protein tyrosine phosphatase receptor type E (PTPRE) regulates the activation of wild-type KIT and KIT mutants differently

Activation of receptor tyrosine kinases needs tight control by tyrosine phosphatases to keep their normal function. In this study, we investigated the regulation of activation of the type III receptor tyrosine kinase KIT by protein tyrosine phosphatase receptor type E (PTPRE). We found that PTPRE can associate with wild-type KIT and inhibit KIT activation in a dose-dependent manner, although the activation of wild-type KIT is dramatically inhibited even when PTPRE is expressed at low level. The D816V mutation of KIT is the most frequently found oncogenic mutation in mastocytosis, and we found that PTPRE can associate and inhibit the activation of KIT/D816V in a dose dependent manner, but the inhibition is much weaker compared with wild-type KIT. Similar to mastocytosis, KIT mutations are the main oncogenic mutations in gastrointestinal stromal tumors (GISTs) although GISTs carry different types of KIT mutations. We further studied the regulation of the activation of GISTs-type KIT mutants and other mastocytosis-type KIT mutants by PTPRE. Indeed, PTPRE can almost block the activation of GISTs-type KIT mutants, while the activation of mastocytosis-type KIT mutants is more resistant to the inhibition of PTPRE. Taken together, our results suggest that PTPRE can associate with KIT, and inhibit the activation of both wild-type KIT and GISTs-type KIT mutants, while the activation of mastocytosis-type KIT mutants is more resistant to PTPRE.     

Characterization of mutants that change the hydrogen bonding of the semiquinone radical at the QH site of the cytochrome bo3 from Escherichia coli

The cytochrome bo3 ubiquinol oxidase catalyzes the two-electron oxidation of ubiquinol in the cytoplasmic membrane of Escherichia coli, and reduces O2 to water. This enzyme has a high affinity quinone binding site (QH), and the quinone bound to this site acts as a cofactor, necessary for rapid electron transfer from substrate ubiquinol, which binds at a separate site (QL), to heme b. Previous pulsed EPR studies have shown that a semiquinone at the QH site formed during the catalytic cycle is a neutral species, with two strong hydrogen bonds to Asp-75 and either Arg-71 or Gln-101. In the current work, pulsed EPR studies have been extended to two mutants at the QH site. The D75E mutation has little influence on the catalytic activity, and the pattern of hydrogen bonding is similar to the wild type. In contrast, the D75H mutant is virtually inactive. Pulsed EPR revealed significant structural changes in this mutant. The hydrogen bond to Arg-71 or Gln-101 that is present in both the wild type and D75E mutant oxidases is missing in the D75H mutant. Instead, the D75H has a single, strong hydrogen bond to a histidine, likely His-75. The D75H mutant stabilizes an anionic form of the semiquinone as a result of the altered hydrogen bond network. Either the redistribution of charge density in the semiquinone species, or the altered hydrogen bonding network is responsible for the loss of catalytic function.     

Exploring single-domain antibody thermostability by molecular dynamics simulation

Abstract

Single-domain antibodies also known as nanobodies are recombinant antigen-binding domains that correspond to the heavy-chain variable region of camelid antibodies. Previous experimental studies showed that the nanobodies have stable and active structures at high temperatures. In this study, the thermal stability and dynamics of nanobodies have been studied by employing molecular dynamics simulation at different temperatures. Variations in root mean square deviation, native contacts, and solvent-accessible surface area of the nanobodies during the simulation were calculated to analyze the effect of different temperatures on the overall conformation of the nanobody. Then, the thermostability mechanism of this protein was studied through calculation of dynamic cross-correlation matrix, principal component analyses, native contact analyses, and root mean square fluctuation. Our results manifest that the side chain conformation of some residues in the complementarity-determining region 3 (CDR3) and also the interaction between α-helix region of CDR3 and framework2 play a critical role to stabilize the protein at a high temperature.

Communicated by Ramaswamy H. Sarma     

Probing subtle conformational changes induced by phosphorylation and point mutations in the TIR domains of TLR2 and TLR3

Extensive research performed on Toll‐like receptor (TLR) signaling has identified residues in the Toll/interleukin‐1 receptor (TIR) domains that are essential for its proper functioning. Among these residues, those in BB loop are particularly significant as single amino acid mutations in this region can cause drastic changes in downstream signaling. However, while the effect of these mutations on the function is well studied (like the P681H mutation in TLR2, the A795P mutation in TLR3, and the P714H mutation in TLR4), their influence on the dynamics and inter‐residue networks is not well understood. The effects of local perturbations induced by these mutations could propagate throughout the TIR domain, influencing interactions with other TIR domain‐containing proteins. The identification of these subtle changes in inter‐residue interactions can provide new insights and structural rationale for how single‐point mutations cause drastic changes in TIR–TIR interactions. We employed molecular dynamics simulations and protein structure network (PSN) analyses to investigate the structural transitions with special emphasis on TLR2 and TLR3. Our results reveal that phosphorylation of the Tyr 759 residue in the TIR domain of TLR3 introduces rigidity to its BB loop. Subtle differences in the intra BB loop hydrogen bonding network between TLR3 and TLR2 are also observed. The PSN analyses indicate that the TIR domain is highly connected and pinpoints key differences in the inter‐residue interactions between the wild‐type and mutant TIR domains, suggesting that TIR domain structure is prone to allosteric effects, consistent with the current view of the influence of allostery on TLR signaling.     

Helix 8 Leu in the CB1 cannabinoid receptor contributes to selective signal transduction mechanisms

The intracellular C-terminal helix 8 (H8) of the CB(1) cannabinoid receptor deviates from the highly conserved NPXXY(X)(5,6)F G-protein-coupled receptor motif, possessing a Leu instead of a Phe. We compared the signal transduction capabilities of CB(1) with those of an L7.60F mutation and an L7.60I mutation that mimics the CB(2) sequence. The two mutant receptors differed from wild type (WT) in their ability to regulate G-proteins in the [(35)S]guanosine 5'-3-O-(thio)triphosphate binding assay. The L7.60F receptor exhibited attenuated stimulation by agonists WIN-55,212-2 and CP-55,940 but not HU-210, whereas the L7.60I receptor exhibited impaired stimulation by all agonists tested as well as by the inverse agonist rimonabant. The mutants internalized more rapidly than WT receptors but could equally sequester G-proteins from the somatostatin receptor. Both the time course and maximal N-type Ca(2+) current inhibition by WIN-55,212-2 were reduced in the mutants. Reconstitution experiments with pertussis toxin-insensitive G-proteins revealed loss of coupling to Galpha(i3) but not Galpha(0A) in the L7.60I mutant, whereas the reduction in the time course for the L7.60F mutant was governed by Galpha(i3). Furthermore, Galpha(i3) but not Galpha(0A) enhanced basal facilitation ratio, suggesting that Galpha(i3) is responsible for CB(1) tonic activity. Co-immunoprecipitation studies revealed that both mutant receptors were associated with Galpha(i1) or Galpha(i2) but not with Galpha(i3). Molecular dynamics simulations of WT CB(1) receptor and each mutant in a 1-palmitoyl-2-oleoylphosphatidylcholine bilayer suggested that the packing of H8 is different in each. The hydrogen bonding patterns along the helix backbones of each H8 also are different, as are the geometries of the elbow region of H8 (R7.56(400)-K7.58(402)). This study demonstrates that the evolutionary modification to NPXXY(X)(5,6)L contributes to maximal activity of the CB(1) receptor and provides a molecular basis for the differential coupling observed with chemically different agonists.