Abstract
Constitutive activation of the RET receptor tyrosine kinase underlies the genesis and progression of multiple endocrine neoplasia type 2 (MEN 2), a dominantly inherited cancer predisposition. Importantly, although kinase activation represents a common theme in neoplasias, not all activating mutations are functionally equivalent. Consistent with this, we ascertained a patient with classical features of MEN 2B, but lacking either of the classical mutations in RET (M918T or A883F). Instead, the patient harbors a novel pair of germ line missense mutations in cis at codons 804 and 805. We evaluated the potential physiochemical effects of these substitutions in silico, predicting both to be moderately deleterious in isolation, but severely deleterious in combination. Consistent with this postulate, we show that the identified tandem mutations (V804M/E805K) are biologically active, transforming cells in culture and that their transforming capacity in combination is distinctly synergistic. Furthermore, the V804M/E805K tandem lesion confers resistance to the small molecule receptor tyrosine kinase inhibitor, PP1, suggesting a mode of action distinct from that known for classical MEN 2B mutations. To address this question, we used homology molecular modeling in silico to model the active site of RET. We predict that RET804 constitutes a critical gatekeeper residue that, when mutated in combination with RET805, induces a conformational change in the hinge region that locks the active site in a position permissive for ATP hydrolysis. Our findings have implications both in the clinic and in the successful development of novel kinase-targeted anticancer drugs. (Cancer Res 2006; 66(20): 10179-87)
Introduction
Inappropriate activation of receptor tyrosine kinases (RTK) is central to the development and progression of many human neoplastic diseases. One RTK implicated in a series of cancers is RET, a 120-kDa membrane-spanning receptor that transduces signals through multiple intracellular pathways upon activation by the glial cell line–derived neurotrophic factor (GDNF) family of ligands, and in complex with the glycosyl phosphatidyl inositol–anchored glial cell–derived neurotrophic factor family receptor-α coreceptors. RET is expressed at various stages of development and, in the adult, in cell populations derived from the branchial arches, including the parathyroid; from the neural crest, including the brain, parasympathetic and sympathetic ganglia, thyroid C-cells, adrenal medulla, and enteric ganglia; and in the urogenital system (1), and mediates signals proposed to influence cell proliferation, differentiation, migration, and apoptosis (2–6).
RET (re-arranged during transfection), was originally identified as a proto-oncogene in a classic cell transformation study (7, 8), and then later in a proportion of papillary thyroid carcinomas (9). The gene comprises 21 exons encompassing approximately 52 kb of genomic DNA on human chromosome 10q11.2 and generates a transcript that is subject to alternative splicing at both the amino and carboxyl termini, resulting in multiple isoforms. The largest of these encodes a 1,114-residue protein possessing 51 COOH-terminal residues (RET51), whereas the COOH-terminal of the short isoform comprises nine amino acids (RET9; 1,072 aa). Constitutive activating mutations in the RET gene were subsequently identified in patients with a dominantly inherited cancer syndrome called multiple endocrine neoplasia type 2 (MEN 2; OMIM 162300 and 171400; refs. 10–12).
MEN 2 is characterized by medullary thyroid carcinoma (MTC) with or without pheochromocytoma and hyperparathyroidism and is subclassified as either (a) familial medullary thyroid carcinoma (FMTC), where MTC is the sole indicator of disease; (b) type 2A (MEN 2A), where, in addition to MTC, pheochromocytomas and/or hyperparathyroidism associate, albeit with incomplete penetrance; or (c) type 2B (MEN 2B), in which MTC is accompanied by developmental abnormalities such as a marfanoid habitus, intestinal dysfunction and ganglioneuromas, thickening of the corneal nerves, neuromas of the lips, tongue, and conjunctiva, and skeletal abnormalities typified by a tall slender physique (5, 13, 14).
Almost all of the MEN 2-causative germ line mutations identified in RET are missense mutations concentrated in a small fraction of the open reading frame. Mutations in patients with MEN 2A exclusively affect just five cysteine residues in the extracellular domain (15). Mutation of any one of these cysteines creates an unpaired residue, which can dimerize with another mutant RET molecule via an illegitimate disulfide bond resulting in constitutive activation of the signaling pathway.
However, the more aggressive MEN 2B subtype is caused exclusively by germ line mutations in the intracellular tyrosine kinase domain of RET (13–17). Almost all (∼95%) of the RET mutations in MEN 2B are confined to a single amino acid in the intracellular region, Met918Thr (12, 18–20). This residue lies within the substrate-binding pocket of the RET tyrosine kinase and substitution by threonine is thought to alter the substrate specificity of the kinase (21, 22). We (23) and others (24) have previously identified a germ line dinucleotide mutation at codon 883 (A883F) in patients with MEN 2B lacking the M918T mutation, suggesting that MEN 2B could be caused by RET mutations other than M918T. This possibility was further strengthened when a combination of two germ line missense mutations in cis at valine 804 and tyrosine 806 was reported in a single patient with MEN 2B-like clinical presentation but without either the M918T or A883F mutations (25). These mutations were then shown to be transforming (26). Subsequently, germ line missense mutations at valine 804 and serine 904 were shown to cosegregate with MTC and mucosal neurilemomas in a family with MEN 2B (27). Interestingly, mutation of codon 805 either singly or in combination with codon 804 has not been reported. Furthermore, mutation of valine 804 also characterizes FMTC cases (28) but, until now, the mechanism by which this mutation activates RET was unknown.
Here, we propose a new mechanism by which novel tandem mutations in the tyrosine kinase domain of RET activate this classical RTK. Mutation of valine 804 and glutamate 805 in tandem affects the hinging motion of the kinase lobes, thereby altering the active site and triggering unsolicited activation of a signaling pathway causing the cancer syndrome MEN 2B. This mechanism is distinct from previously reported activating mutations in RET and MEN 2.
Materials and Methods
Mutation detection. Genomic DNA extracted from the proband's peripheral blood was analyzed for mutations in exons 10 to 16 of RET by DNA sequencing. Blood from her children, sister, and mother were also analyzed; paternal tissue was not available. Sequencing was done in the forward and reverse direction using the BigDye Terminator cycle sequencing kit (Applied Biosystems, Foster City, CA).
Physicochemical conservation and mutational tolerance. We quantified the physicochemical changes in mutant RET proteins against evolutionary variation using multivariate analysis of protein polymorphism (MAPP, see ref. 29). MAPP uses a multiple sequence alignment of protein orthologs to estimate the physicochemical constraints acting on the amino acids observed at each homologous position. The fitness of a substituted protein can then be predicted by a score quantifying the degree to which its amino acid sequence violates these estimated constraints. For this study, the MAPP scores of each substituted variant were calculated as in ref. 30, using the published alignment of human RET with its 11 orthologs (30). Following ref. 29, we obtained conservative P values for each score by comparing one-thirteenth its square to an F6,6 distribution (29, 30). To obtain a composite score for tandem substitutions in cis, we summed the squared MAPP scores of each substitution and took the square root of the result. We obtained a conservative P value for each composite score by comparing one-thirteenth of its square to the convolution of two F6,6 distributions.
Site-directed mutagenesis. Mutations were introduced into both the short (RET9) and long (RET51) 3′-splice isoforms of the human RET cDNA in the plasmid, pRCCMV, using the QuikChange site-directed mutagenesis kit (Stratagene, La Jolla, CA).
Cell lines. Cell lines were maintained in DMEM supplemented with 10% Colorado calf serum (Celbio, Italy). Cells were grown at 37°C in 10% CO2/90% air.
In vitro focus formation, kinase, and PP1 inhibitor assays. The expression constructs were transfected into mouse fibroblast NIH-3T3 cells and their ability to transform assayed as previously described (31). Growth factor stimulation was achieved by adding either GDNF or neurturin (10 ng/μL) to the growth medium. Transformed foci were counted after 2 to 3 weeks and RET expression and phosphorylation checked for each cell line. PP1 inhibitor [4-amino-5-(4-methylphenyl)-7-(-t-butyl)-pyrazolo-d-3,4-pyrimidine; BioMol Research Laboratories, Plymouth Meeting, PA], was dissolved in DMSO and diluted into the culture medium to 1 μmol/L for 12 hours. Cells were lysed as reported (32), immunoprecipitated with anti-RET polyclonal antibody (C-19; Santa Cruz Biotechnology, Santa Cruz, CA), and protein concentration estimated by a modified Bradford assay (Bio-Rad, Munich, Germany). Immune complexes were separated by SDS-PAGE and immunoblotted using a mouse antiphosphotyrosine antibody (4G10; Upstate Biotechnology Inc., Lake Placid, NY).
In silico homology molecular modeling. Based on the coordinates of the active insulin receptor kinase (IRK; ref. 33), a homology model of RET was constructed using the Homology module in Insight II, a homology molecular modeling software program (Accelrys, Inc., San Diego, CA). In this model (residues 713-1018), Val804 and Glu805 correspond to residues Met1076 and Glu1077 in IRK, and Met120 and Glu121 in protein kinase A (PKA).
Results
Clinical history of the patient. The female proband (Fig. 1A) presented at 50 years of age with diarrhea and was diagnosed with endometrial adenocarcinoma. She was also found to have a clear cell carcinoma of one ovary, and endometrial cancer that had spread to the contralateral ovary. She was treated by hysterectomy and bilateral oophorectomy followed by radiotherapy. On follow-up 10 months later, she had a persistently increased CEA level and an 18FDG PET scan showed uptake in the thyroid and adrenals. Subsequently, she was also found to have a thyroid mass, with an increased blood calcitonin level. A fine-needle aspirate showed medullary carcinoma of the thyroid, and a methyliodobenzylguanidine scan preoperatively showed bilateral pheochromocytomas, confirmed biochemically by increased urinary catecholamines and by histology. She was treated by bilateral laparoscopic adrenalectomies, during which a liver biopsy showed metastatic medullary cell thyroid carcinoma. Total thyroidectomy was done and histology confirmed multifocal medullary thyroid carcinoma invading into local nerve and muscle. Consistent with these observations, the patient possessed a marfanoid habitus with a typical MEN 2B facial appearance. She had prominent nodules on her lips and tongue (Fig. 1A), a small pigmented lesion on the skin of her right calf and ophthalmologic examination revealed thickening of the corneal nerves, all features of the MEN 2B syndrome. Apart from diverticulitis, she has remained well 4 years later. There was no family history of cancer except for her father, who died of carcinoma of the bladder at 60 years of age.
Clinical molecular genetics of the patient. A, pedigree of the MEN 2B family (arrow, proband); clinical features of MEN 2B (bumpy lips and tongue) in the proband. B, the location of the tandem mutations in the RET intracellular TK domain (circled). C, DNA sequencing of one allele of the proband's subcloned DNA revealed the presence of tandem missense mutations at codons 804 and 805 in exon 14 of the RET gene. D, MAPP-based P values of substituted RET variants quantify their predicted phenotypic effect. Tandem mutations identified in FMTC are less radical than those identified in MEN 2B. Columns, P values of five RET tandem mutants, from left to right: RETV804M/V778I and RETV804M/R844L, identified in FMTC, and RETV804M/E805K, RETV804M/S904C, and RETV804M/Y806C, identified in MEN 2B.
Clinical molecular genetics of the patient. A, pedigree of the MEN 2B family (arrow, proband); clinical features of MEN 2B (bumpy lips and tongue) in the proband. B, the location of the tandem mutations in the RET intracellular TK domain (circled). C, DNA sequencing of one allele of the proband's subcloned DNA revealed the presence of tandem missense mutations at codons 804 and 805 in exon 14 of the RET gene. D, MAPP-based P values of substituted RET variants quantify their predicted phenotypic effect. Tandem mutations identified in FMTC are less radical than those identified in MEN 2B. Columns, P values of five RET tandem mutants, from left to right: RETV804M/V778I and RETV804M/R844L, identified in FMTC, and RETV804M/E805K, RETV804M/S904C, and RETV804M/Y806C, identified in MEN 2B.
The V804M/E805K tandem mutations are present in cis and cause MEN 2B. DNA sequencing revealed the presence of two missense mutations in exon 14 of RET. At codon 804, methionine was substituted for valine (GTG → ATG), whereas at codon 805, lysine was substituted for glutamate (CAG → AAG). The tandem V804M/E805K mutations lie within the RET tyrosine kinase domain (Fig. 1B). No further mutations were detected in exons 10 to 16. Although DNA from the proband's father was not available for analysis, neither mutation was found in the proband's children or her mother, suggesting that the tandem mutations were present on the same allele.
To confirm this, a region spanning exon 14 was amplified from genomic DNA by PCR, subcloned and analyzed by sequencing and digestion with RcaI (BspHI). DNA sequencing of the amplicon revealed the presence of both mutations in cis (Fig. 1C). RcaI restriction analysis of the amplicon produced products of 230 and 73 bp (data not shown), further confirming that both mutations were present in the cis configuration; mutation of codons 804 or 805 alone is not sufficient to create the RcaI recognition site.
Correlating tandem mutations with observed phenotypes. The phenotypic effect of disease-causing changes in the RET protein sequence varies inversely with the degree to which the physicochemical properties of the substituted amino acid are conserved throughout vertebrate evolution (30). We reasoned that the effect of tandem substitutions in RET, including RETV804M/E805K, should similarly correlate with aggregate physicochemical change. To substantiate this, we identified, from the literature, four additional reports of tandem mutations in RET: (a) RETV804M/Y806C (25) and (b) RETV804M/S904C (27) in MEN 2B, and (c) RETV804M/V778I (34) and (d) RETV804M/R844L (35) in FMTC. We then extended the physicochemical conservation metric used in the previous study (MAPP, see Materials and Methods) to accommodate multiple substitutions in cis. Insofar as MEN 2B is a more severe phenotype than FMTC, we expected the tandem mutations identified in FMTC to be less radical than their MEN 2B counterparts in terms of aggregate physicochemical change from wild-type. This expectation was strongly confirmed (Fig. 1D), and we proceeded to predict the effect of RETV804M/E805K using the same physicochemical criteria. Consistent with our hypothesis, the physicochemical change in RETV804M/E805K is consistent with the other MEN 2B cases and is greater than that observed in FMTC (Fig. 1D).
We next sought physicochemical context for the synergy between V804M and E805K in transforming potential. To accomplish this, we scored RETV804M/E805K against its constituent single substitution TK mutants, V804M and E805K, normalizing the tandem mutant appropriately for comparison (see Materials and Methods). The aggregate physicochemical change in the tandem mutant, as quantified by its MAPP score (31.56) and associated P values (6.5 × 10−4), suggested a phenotype more severe than that of either RETV804M (score, 24.43; P = 9.4 × 10−5) or RETE805K (score, 19.98; P = 3.0 × 10−4). Equating this severity to transforming potential, we then tested these predictions directly on cells in culture.
The V804M/E805K tandem mutations are biologically active and transform cells. The transforming activity of the RETV804M/E805K mutations in both short (RET9) and long (RET51) 3′-splice isoforms (Fig. 1B) was assayed by transfecting NIH-3T3 cells in parallel with the wild-type cDNA and with other known activating mutations (i.e., E768D, V804M, A883F, and M918T). In order to quantify the transforming potential of these mutations, we assayed the number of focus-forming units per microgram of DNA transfected, either in the absence or presence of chronic growth factor stimulation. We used the ligands, GDNF and neurturin, because experiments in our lab indicated that the receptors for these particular growth factors were expressed in NIH-3T3 cells.10
Unpublished results.
The transforming potential of the kinase mutants in the RET 3′-short isoform is very low in the absence of either ligand (Fig. 2A). Upon chronic stimulation with either GDNF or neurturin, each of the mutations, with the exception of E805K, displayed increased transforming activity (Fig. 2A). In contrast, the E805K mutation in the short isoform is barely stimulated above a basal level of activity by the addition of either ligand (Fig. 2A).
Transforming activity of the RET mutants. Columns, transformation potential of V804M, E805K, or the V804M/E805K mutations in a RET cDNA encoding the short isoform (A), or the long isoform (B). Other known RET gain-of-function mutations were assayed for comparison. White columns, no growth factor ligand present; black columns, GDNF added to the culture medium at 10 ng/μL; gray columns, neurturin added to the culture medium at 10 ng/μL.
Transforming activity of the RET mutants. Columns, transformation potential of V804M, E805K, or the V804M/E805K mutations in a RET cDNA encoding the short isoform (A), or the long isoform (B). Other known RET gain-of-function mutations were assayed for comparison. White columns, no growth factor ligand present; black columns, GDNF added to the culture medium at 10 ng/μL; gray columns, neurturin added to the culture medium at 10 ng/μL.
The tandem V804M/E805K mutations in the long isoform, examined in the absence of ligand stimulation, induced 42 transformed foci per microgram of DNA transfected. This is about half the activity of the classical MEN 2B mutations, A883F and M918T (Fig. 2B). Moreover, the transforming potential of the tandem mutations is higher than the sum of the individual mutations alone (V804M/E805K, V804M, and E805K; 42 foci versus 9 versus 8, respectively), indicating that there is synergy when the mutations are in combination. Importantly, the transforming potential of the other TK mutants tested (E768D, V804M, and E805K) is between 8- and 16-fold lower than that of the RETA883F mutation and between 10- and 19-fold lower than that of the RETM918T mutation on the same isoform, confirming that these kinase mutants are only capable of conferring weak oncogenic potential to RET (Fig. 2B). With the exception of the A883F, M918T, and tandem V804M/E805K mutations, the transforming activity of the kinase mutations is further increased by the addition of either GDNF or neurturin. Thus, the transforming activity of the tandem V804M/E805K mutations or classical MEN 2B mutations (A883F and M918T) is growth factor stimulation–independent.
The PP1 inhibitor has specific effects on mutant RET oncoprotein transformation. The small molecular weight pyrazolo-pyrimidine inhibitor, PP1, abrogates the morphologic changes induced by the expression of RETMEN 2A and RETMEN 2B oncoproteins in NIH-3T3 cells and induced dephosphorylation of these mutant proteins (36). To determine whether the V804M, E805K, and V804M/E805K-transformed cells would respond similarly, we cultured each in the presence of PP1.
Consistent with the finding for RETMEN 2A and RETMEN 2B oncoproteins, when examined by phase-contrast microscopy, the morphology of the cells harboring the E805K mutation reverted upon PP1 treatment. In contrast, no morphologic reversion was observed for the cells harboring either the V804M or the tandem mutations (Fig. 3A).
Effects of PP1 on RETV804M-, RETE805K-, and RETV804M/E805K-transformed cells. A, PP1 reverts the transformed morphology of RETE805K, but not that of RETV804M and RETV804M/E805K mutant NIH-3T3 fibroblasts. B, PP1 induces the dephosphorylation of RETE805K but not RETV804M or RETV804M/E805K. Anti-RET immunocomplexes (IP) for RETV804M-, RETE805K-, and RETV804M/E805K-expressing fibroblasts.
Effects of PP1 on RETV804M-, RETE805K-, and RETV804M/E805K-transformed cells. A, PP1 reverts the transformed morphology of RETE805K, but not that of RETV804M and RETV804M/E805K mutant NIH-3T3 fibroblasts. B, PP1 induces the dephosphorylation of RETE805K but not RETV804M or RETV804M/E805K. Anti-RET immunocomplexes (IP) for RETV804M-, RETE805K-, and RETV804M/E805K-expressing fibroblasts.
To assess the effects of PP1 on RET phosphorylation, cells were incubated in the absence or presence of PP1 prior to lysis. RET tyrosine-phosphorylated proteins were analyzed by Western blotting with anti-pTyr antibodies. Consistent with cellular morphology data, PP1 dramatically reduced the amount of phosphorylated RET in cultures harboring the E805K mutation; PP1 treatment was ineffective in reducing phosphorylation in RET proteins harboring the V804M or tandem V804M/E805K mutations (Fig. 3B).
In silico structural modeling predicts the V804M/E805K tandem mutations alter the RET kinase hinge region and reposition the active site cleft such that it is constitutively active. In our homology molecular model of the RET kinase, based on the structure of active IRK (33), and a sequence alignment with IRK (37) and PKA (Fig. 4), the kinase domain starts at position 713 and ends at position 1018. Residues Val804 and Glu805 correspond to Met1076 and Glu1077 in IRK, and Met120 and Glu121 in PKA, respectively.
Homology alignment of the RET tyrosine kinase domain. Our model is based on the active IRK and the RET kinase starts at position 713 and ends at position 1018. Tandem V804M/E805K mutations (encircled); conserved residues (•); mutated residues (▴). The alignment for RET and IRK kinases was taken from ref. 33. PKA was aligned to RET based on the alignment from PileUp multiple sequence alignment program from Genetic Computer Group Wisconsin Package Programs.
Homology alignment of the RET tyrosine kinase domain. Our model is based on the active IRK and the RET kinase starts at position 713 and ends at position 1018. Tandem V804M/E805K mutations (encircled); conserved residues (•); mutated residues (▴). The alignment for RET and IRK kinases was taken from ref. 33. PKA was aligned to RET based on the alignment from PileUp multiple sequence alignment program from Genetic Computer Group Wisconsin Package Programs.
The model reveals a two-lobed structure typical of tyrosine kinase catalytic domains, in which Val804 and Glu805 are located at the junction of the fifth β-strand and the linker segment, which joins the small and large lobes of the kinase core (Fig. 5A and B). This linker is the major component of the hinge mechanism that facilitates the opening and closing of the catalytic cleft (38, 39), hence, mutation of this region could affect the kinase hinging mechanism. It also forms a critical part of the adenine binding pocket, in which the backbone carbonyl hydrogen of Glu805 bonds to ATP (via the N-hydrogen of the 6-amino adenine ring in ATP). In our structural model, Glu805 lies in close proximity to Lys889 in the eighth β-strand at distances of R = 2.9 Å and R = 3.3 Å. This interaction is conserved in the IRK in which Glu1077 (corresponding to Glu805 in the RET model) interacts with Lys1147 (corresponding to Lys889 in RET). There is a similar interaction in PKA albeit at longer distances: Glu121 interacts with Gln181 (corresponding to Lys889 in RET) at R = 3.6 Å. In PKA, the closest positively charged residue is Lys105 at the beginning of the fourth β-strand at R = 3.9 Å. The E805K mutation changes the character of the amino acid at this position, conferring a positive charge and thus creating charge repulsion with Lys889. This newly created Lys805 has the potential to interact with another nearby negatively charged residue, the closest candidate being Glu884, positioned between the seventh and eighth β-strands. The reduction or inhibition of activity conferred by this mutation in the shorter RET isoform could thus be due to charge repulsion between Lys805 and Lys889 and/or formation of an inappropriate ion pair that locks the enzyme into an inactive conformation. As this site is close to the hinge and also contributes to the adenine binding pocket, it is likely to be critical for generating an active enzyme.
Homology molecular modeling of the RET tyrosine kinase domain. A, ribbon diagram of the wild-type RET tyrosine kinase domain (residues 713-1018). Small lobe (white); large lobe (bronze); the linker, or hinge, region (residues 806-811; red); and the insert region (green). The conserved kinase residues (Gly731, Gly733, Gly736, Lys758, Glu775, Asp874, Asn879, Asp892, Asp933, Glu921, and Arg995; gray balls representing the α-carbon atoms). The mutated residues E768D, V804M, E805K, A883F, and M918T (red balls). The Tyr905 residue, which is to be phosphorylated (turquoise ball). B, two ribbon structure models are compared in the tandem mutation region: the wild-type model (gray) and the tandem mutated model (bronze). Val804 and Glu805 (gray), the mutated residues Met804 and Lys805 (green). Lys889, which potentially interacts with Glu805 in the wild-type model, and Glu884 in the tandem mutated structure, which possibly interacts with Lys805 (yellow). The ATP in both structures (thick yellow line, wild-type structure; thin yellow line, tandem mutated structure). Red, the linker, or hinge, region (residues 806-811); white, the DFG region (residues 892-894). C, superimposition of the α-carbon backbones of the inactive (pink) and active (white) IRK. Ribbons, the α-helix C (residues 1038-1052) and hinge regions (residues 1078-1083). The side chains of Met1051 and Met1076 in the inactive IRK (turquoise) and in the active IRK (yellow). Dashed line, the distances between Met1051 and Met1076. D, in the enlargement of (C), the hydrophobic residues Val1074, Leu1062, Val1010, and Val1060 are added to both structures.
Homology molecular modeling of the RET tyrosine kinase domain. A, ribbon diagram of the wild-type RET tyrosine kinase domain (residues 713-1018). Small lobe (white); large lobe (bronze); the linker, or hinge, region (residues 806-811; red); and the insert region (green). The conserved kinase residues (Gly731, Gly733, Gly736, Lys758, Glu775, Asp874, Asn879, Asp892, Asp933, Glu921, and Arg995; gray balls representing the α-carbon atoms). The mutated residues E768D, V804M, E805K, A883F, and M918T (red balls). The Tyr905 residue, which is to be phosphorylated (turquoise ball). B, two ribbon structure models are compared in the tandem mutation region: the wild-type model (gray) and the tandem mutated model (bronze). Val804 and Glu805 (gray), the mutated residues Met804 and Lys805 (green). Lys889, which potentially interacts with Glu805 in the wild-type model, and Glu884 in the tandem mutated structure, which possibly interacts with Lys805 (yellow). The ATP in both structures (thick yellow line, wild-type structure; thin yellow line, tandem mutated structure). Red, the linker, or hinge, region (residues 806-811); white, the DFG region (residues 892-894). C, superimposition of the α-carbon backbones of the inactive (pink) and active (white) IRK. Ribbons, the α-helix C (residues 1038-1052) and hinge regions (residues 1078-1083). The side chains of Met1051 and Met1076 in the inactive IRK (turquoise) and in the active IRK (yellow). Dashed line, the distances between Met1051 and Met1076. D, in the enlargement of (C), the hydrophobic residues Val1074, Leu1062, Val1010, and Val1060 are added to both structures.
In contrast to Glu805, Val804 is hydrophobic and is part of a hydrophobic nodule that is comprised of residues from different parts of the protein: Ile788 (loop between C helix and fourth β-strand), Ala756 (third β-strand), Leu779 (C helix), Leu790 (fourth β-strand), and Leu802 (fifth β-strand). This hydrophobic pocket is conserved in IRK (Val1060, Val1010, Met1051, Leu1062, and Val1074) and in PKA (Val104, Val57, Leu95, Leu106, and Met118) and may be important for the hinging motions discussed above. In fact, this hydrophobic pocket appears to be a conserved feature of most, if not all, protein kinases. In the IRK structure, the residue corresponding to Val804 is actually a methionine (Met1076), and this methionine interacts with the hydrophobic pocket described above. Following mutation, the interactions of V804M are similar with residues Ile788, Ala756, Leu779, Leu790, and Leu802. Even though the methionine residue is larger than the valine, it seems to fit well in this hydrophobic environment. However, the longer methionine does reach better to the COOH terminus of the C helix. Indeed, both residues, 804 and 805, are close to the COOH terminus of the C helix, and it is this helix that must be aligned properly in the active enzyme. Correct positioning of the C helix is an essential requirement of the activation process for most protein kinases. The evolutionary physicochemical and steric constraints evident within this region strongly suggest a critical role for this configuration and a significantly deleterious outcome to mutations therein.
Discussion
Clearly, the ability to discriminate points along the spectrum between neutral and deleterious amino acid substitutions in the genes of patients remains a significant challenge to our understanding of disease mechanisms and consequently to patient care. Importantly, although so-called classical mutations can teach us much about the pathogenesis of disease, it is the exceptions to a rule that are frequently responsible for expanding the boundaries of our understanding. We have identified novel missense mutations affecting codons 804 and 805 on the same allele of the RET proto-oncogene. We report detailed evaluation of these tandem mutations, elucidating their biological characteristics both in isolation and in combination and propose a novel mechanism by which they trigger the activation of RET resulting in MEN 2B.
To date, germ line point mutations in exon 14 leading to substitution of RET Val804 by either leucine or methionine have been reported in patients with sporadic MTC or FMTC (28, 40–44). However, the V804M mutation has only been reported in patients with MEN 2B when it is present on the same RET allele with a missense mutation at either codon 806 (25) or codon 904 (27). Neither germ line nor somatic mutations of codon 805 have been reported previously in MEN 2 or FMTC. Notwithstanding the possibility that the mutations could have arisen de novo in previous generations, we proposed that the pairing of V804M and E805K may exemplify a novel class of RET mutations, whose tumorigenic activity in combination is uncharacteristically severe. We therefore set out to evaluate the effect of these novel tandem mutations on the MEN 2 disease phenotype.
Given that the phenotypic effect of disease-causing substitutions in the RET protein sequence varies inversely with physicochemical and evolutionary constraint (30), we posited, that the effect of the tandem mutations in RET correlates strongly with the aggregate physicochemical change. Furthermore, insofar as MEN 2B is a more severe phenotype than FMTC, we showed that the physicochemical consequences of RETV804M/E805K more closely resemble the other MEN 2B cases than those observed in FMTC.
To test our hypotheses directly, we did classic transformation studies introducing the mutations into the short (RET9) and long (RET51) splice isoforms, producing data that strongly supported our expectations. For example, in the long isoform, the tandem RETV804M/E805K mutations displayed a level of transforming activity greater than that of either mutation alone, and greater than the sum of its parts, indicating that this combination of mutations is highly synergistic. This transforming activity is also stronger than that of two weakly activating RET mutations (E768D and V804M), but not quite as potent as the highly activating mutations, A883F and M918T, which are historically associated with MEN 2B; this may also explain the relatively late age-of-onset of MEN 2B in this patient. However, like the A883F and M918T mutations, the level of activation of V804M/E805K is not enhanced greatly in the presence of ligand. This suggests that, at least in the long isoform, these mutant receptors do not need to dimerize to achieve activation.
To address the mechanism by which the tandem V804M/E805K mutations activate RET, we did comparative homology molecular modeling experiments in which we generated a structural model of the RET tyrosine kinase domain (Fig. 5A). The model clearly predicts a classical two-lobed structure with a relatively open ATP binding site (Fig. 5A). The tandem V804M/E805K mutations are located at a critical “hub” site that lies at the major hinge point between the small and large kinase lobes, which form the cytoplasmic TK catalytic domains (Fig. 5A and B). Our results show that mutations in this region have the potential to be activating or inhibiting. Several of the major elements of secondary structure converge at this point, therefore, small changes could be propagated over long distances. Val804 and Glu805 contribute in quite different ways to this hinge region.
Val804 contributes to a hydrophobic region, and its replacement by methionine introduces a bulkier group but does not change the properties of the side chain. Most likely, it helps to correctly position the α-helix C for catalysis. One can appreciate the crucial position of this residue and its relationship to the C helix when inactive IRK is compared with activated IRK (Fig. 5C and D). Additionally, the V804M mutation is activating in both isoforms.
Not only is Glu805 located at the junction of the fifth β-strand and the linker segment (residues 806-811), which joins the small and large kinase lobes, but it also contributes to the ATP binding pocket. Specifically, Glu805 hydrogen bonds to the adenine ring in ATP. A similar interaction exists when the active IRK and PKA structures are compared (33). In fact, this interaction is conserved in most protein kinases. Mutation of this Glu805 to lysine changes the electrostatic properties of the region. Lys889, which seems to interact with Glu805 in our wild-type RET model, an interaction that is also conserved in IRK, lies in the eighth β-strand preceding the conserved DFG loop (also termed the magnesium positioning loop). When the ionic interaction between Lys889 and Glu805 is broken following mutation, not only would there be charge repulsion between Lys805 and Lys889, but Lys889 does not have another negatively charged residue with which to bind. This could distort the position of the DFG loop. Clearly then, this region is important for activation, most likely for correctly orienting the C helix.
The underlying mechanism seems to be at least partly dependent on the length of the carboxyl-terminal tail because transforming activity is much higher in the long isoform than in the short isoform, suggesting that the longer carboxyl-terminal tail might directly or indirectly interface with this region. Furthermore, consistent with a mechanism in which the receptor does not need to dimerize to become constitutively activated, the tandem mutations are highly activating even in the absence of exogenous growth factor stimulation. This is consistent with a prediction that the primary role of dimerization is to stabilize an inhibited form of RET rather than to specifically stabilize an active conformation.
We synthesized these results into a working model illustrating the novel hinging mechanism by which the newly identified tandem mutations could constitutively activate the RET kinase domain in the absence of growth factor stimulation (Fig. 6). Mutation of V804 and E805 in the critical kinase hinge region has a 2-fold effect: first, it results in a conformational change, which closes the catalytic cleft; and second, this in turn correctly repositions the C helix and DFG residues, thereby facilitating catalysis. The outcome, release of the kinase from its normal autoinhibited state (45), which can occur in the absence of a growth factor signal, is orientation of the active site cleft into a more optimal position for the binding of ATP and subsequent autophosphorylation of the catalytic tyrosine (Y905). Consequently, the mutant monomeric RETV804M/E805K kinase is “switched on” in the absence of exogenous growth factor signals.
A mechanism for the activation of the RET tyrosine kinase domain by the tandem V804M/E805K mutations. The structure of the kinase is drawn as an interpretation of our homology model as two lobes linked by a hinge region (red). The α-helix C (pink, inactive kinase; white, activated kinase) and DFG conserved sequence (gray) are also shown. Val804 and Glu805 (black “sticks”) are replaced by Met804 and Lys805 (blue “sticks”) in the mutated kinase; the change in conformation of the critical hinge region induced by the tandem mutations, such that the orientation of the kinase lobes are made permissive for ATP (orange) binding, is represented by the direction of the purple arrow in (iii). The catalytic tyrosine (Y), phosphorylation events (P) and carboxyl-terminal tail are also shown. The inactive RET kinase (i) is ordinarily activated by ligand-induced dimerization (ii), which induces a conformational change in the NH2-terminal and COOH-terminal lobes of the tyrosine kinase domain such that the active site cleft is permissive for ATP binding. Note that the orientation of the C helix changes in the activated kinase (pink to white) and that the active site cleft is now optimally positioned for binding ATP. However, when present, (iii) the tandem V804M/E805K mutations activate the kinase in the absence of ligand by inducing a conformational change in the hinge. With the kinase lobes “locked” in this position, the C helix is favorably positioned for binding ATP in the active site cleft. This facilitates kinase activation and subsequent autophosphorylation of the catalytic tyrosine. The result of the tandem mutations is to switch the kinase into the “on” position despite the absence of ligand. For simplicity, the glycosyl phosphatidyl inositol–linked coreceptors are not shown. Pink, cysteine-rich extracellular domain; blue, the ligand homodimer; yellow spheres, the plasma membrane; green, the intracellular kinase domain.
A mechanism for the activation of the RET tyrosine kinase domain by the tandem V804M/E805K mutations. The structure of the kinase is drawn as an interpretation of our homology model as two lobes linked by a hinge region (red). The α-helix C (pink, inactive kinase; white, activated kinase) and DFG conserved sequence (gray) are also shown. Val804 and Glu805 (black “sticks”) are replaced by Met804 and Lys805 (blue “sticks”) in the mutated kinase; the change in conformation of the critical hinge region induced by the tandem mutations, such that the orientation of the kinase lobes are made permissive for ATP (orange) binding, is represented by the direction of the purple arrow in (iii). The catalytic tyrosine (Y), phosphorylation events (P) and carboxyl-terminal tail are also shown. The inactive RET kinase (i) is ordinarily activated by ligand-induced dimerization (ii), which induces a conformational change in the NH2-terminal and COOH-terminal lobes of the tyrosine kinase domain such that the active site cleft is permissive for ATP binding. Note that the orientation of the C helix changes in the activated kinase (pink to white) and that the active site cleft is now optimally positioned for binding ATP. However, when present, (iii) the tandem V804M/E805K mutations activate the kinase in the absence of ligand by inducing a conformational change in the hinge. With the kinase lobes “locked” in this position, the C helix is favorably positioned for binding ATP in the active site cleft. This facilitates kinase activation and subsequent autophosphorylation of the catalytic tyrosine. The result of the tandem mutations is to switch the kinase into the “on” position despite the absence of ligand. For simplicity, the glycosyl phosphatidyl inositol–linked coreceptors are not shown. Pink, cysteine-rich extracellular domain; blue, the ligand homodimer; yellow spheres, the plasma membrane; green, the intracellular kinase domain.
Our data suggest that the RET kinase domain harbors at least two regions in which substitutions might maximally affect protein function. The common M918T mutation is clearly in the first “hot-spot” region, the peptide binding site, but this residue is also very important for activity. Indeed, in the activation loop of Ser/Thr kinases, this threonine is absolutely essential for activity (46), but it does not determine substrate specificity. Equally, RETV804M/E805K and the other MEN 2B mutation, RETA833F, reside around the hinge region, and mutation of this second mutation hot-spot region most likely influences conformational states. The hinge region forms part of the ATP binding pocket for the adenine ring and these mutations seem to shift the enzyme into an active conformation, which is then independent of GDNF. Thus, they release an inhibition that is imposed by the longer COOH-terminal tail in the absence of ligand.
The results obtained with the small molecule TK inhibitor, PP1, underscore the importance of understanding the biochemical events that underlie the activation of protein kinases in cancer. Although effective against the E805K mutation, RET/PTC3—a fusion protein of the RET TK domain with the RFG gene (47, 48), and the Src and PDGFR kinases (49), the kinase inhibitor PP1 did not inhibit RETV804M/E805K or RETV804M. Recently, disease-associated mutations at valine 804 conferred resistance to a series of small molecule kinase inhibitors (50). Our homology model data suggests that, in these mutants, the methionine lies closer to PP1 than did the deposed valine. Specifically, we predict that the methionine of codon 804 lies over the methyl-phenyl ring of PP1 in the ATP-binding pocket and that the close proximity of the mutant residue to the inhibitor is sufficient to prevent competitive inhibition of the active site. Consistent with this postulate, we show that differences exist between the mutated (V804M, E805K, and V804M/E805K) and wild-type RET when they are modeled in combination with ATP (data not shown). Overall, ATP is shifted, especially in the V804M and V804M/E805K mutated kinases, when compared with PP1. This indicates that the difference in response to PP1 is due to a mutation-induced difference in the conformation of the ATP-binding site and suggests that PP1 is noncompetitive with ATP for the inhibition of RETV804M and RETV804M/E805K, but that it is competitive with ATP for the inhibition of RETE805K. Therefore, PP1 is a mixed competitive inhibitor vis-à-vis the substrate. Although it is possible that residue 768 could hydrolyze ATP, mutation of 768 does not lock the hinge region in an open aspect (data not shown). This could explain why mutation of this particular residue is milder and only found in cases of FMTC.
Crucially, in our model of the RET catalytic subunit, Val804 is the “gatekeeper”11
Gatekeeper—a single active site residue that governs nucleotide substrate specificity.
Unpublished data.
Our data also illuminate the differences between the RET V804M, E805K, and V804M/E805K mutations in response to drug treatment. Therefore, successful extension of kinases as druggable targets for other cancers (52) raises many challenges, including the need for molecular analysis of tumors from patients and the development of specific kinase inhibitors with requisite target specificity. This is further illustrated by Gleevec (Glivec/STI-571/imatinib mesylate), a small molecule inhibitor, which initially showed dramatic clinical value in the treatment of patients with BCR-ABL-positive chronic myeloid leukemia and c-KIT-positive gastrointestinal stromal tumors (53), but in which some tumors were clinically resistant to treatment. Investigation into the mechanism of resistance has shown that this occurs through point mutations in the catalytic sites of the BCR-ABL and c-KIT kinase domains (54). Significantly, the tandem V804M/E805K mutations in RET correspond to the key site that most likely accounts for the specificity of Gleevec. In ABL, the gatekeeper residue is a threonine, which allows Gleevec to bind, whereas the most common mutation found in patients resistant to treatment is a threonine to tyrosine mutation, which sterically blocks Gleevec binding and which is functionally equivalent to the Val804Met mutation in RET.
A fuller understanding of the functional consequence of such RET mutations will undoubtedly promote important clinical interventions and should lead to the development of tailored drug therapies for MEN 2-related cancers and its associated sporadic cancers. More broadly, the emerging power of chemical genetics (55), when combined with classical genetic and biochemical data, facilitates a better understanding of protein function and dysfunction.
Acknowledgments
Grant support: Cancer Research UK (B.A.J. Ponder), Télethon (I. Bongarzone), U.S. Army Grant AIBS1762 (S.S. Taylor), and March of Dimes (A.S. McCallion).
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
The authors thank the patient and her family for their continued cooperation. B.A.J. Ponder is a Gibb fellow of Cancer Research UK.





