- Research article
- Open Access
Structural insights into Noonan/LEOPARD syndrome-related mutants of protein-tyrosine phosphatase SHP2 (PTPN11)
BMC Structural Biology volume 14, Article number: 10 (2014)
The ubiquitous non-receptor protein tyrosine phosphatase SHP2 (encoded by PTPN11) plays a key role in RAS/ERK signaling downstream of most, if not all growth factors, cytokines and integrins, although its major substrates remain controversial. Mutations in PTPN11 lead to several distinct human diseases. Germ-line PTPN11 mutations cause about 50% of Noonan Syndrome (NS), which is among the most common autosomal dominant disorders. LEOPARD Syndrome (LS) is an acronym for its major syndromic manifestations: multiple Lentigines, Electrocardiographic abnormalities, Ocular hypertelorism, Pulmonary stenosis, Abnormalities of genitalia, Retardation of growth, and sensorineural Deafness. Frequently, LS patients have hypertrophic cardiomyopathy, and they might also have an increased risk of neuroblastoma (NS) and acute myeloid leukemia (AML). Consistent with the distinct pathogenesis of NS and LS, different types of PTPN11 mutations cause these disorders.
Although multiple studies have reported the biochemical and biological consequences of NS- and LS-associated PTPN11 mutations, their structural consequences have not been analyzed fully. Here we report the crystal structures of WT SHP2 and five NS/LS-associated SHP2 mutants. These findings enable direct structural comparisons of the local conformational changes caused by each mutation.
Our structural analysis agrees with, and provides additional mechanistic insight into, the previously reported catalytic properties of these mutants. The results of our research provide new information regarding the structure-function relationship of this medically important target, and should serve as a solid foundation for structure-based drug discovery programs.
The ubiquitous non-receptor protein tyrosine phosphatase SHP2 (encoded by PTPN11) plays a key role in RAS/ERK signaling downstream of most, if not all growth factors, cytokines and integrins, although its major substrates remain controversial [1, 2]. SHP2 contains two N-terminal SH2 domains, a catalytic (PTP) domain, a C-terminal tail with two tyrosine phosphorylation sites and a proline-rich domain [2–5], and is regulated by an elegant molecular switch mechanism that couples appropriate cellular localization to catalytic activation [3, 5]. In the absence of a tyrosine-phosphorylated binding partner for its SH2 domains (basal state), SHP2 assumes a “closed” conformation wherein the N-terminal SH2 (N-SH2) domain is wedged into the PTP domain, blocking substrate access (Figure 1a). Upon agonist stimulation, recruitment of the N-SH2 domain to specific phosphotyrosyl (pTyr-) peptides disrupts this self-locking conformation, freeing the PTP domain for catalysis [3, 5, 6].
Mutations in PTPN11 cause several human diseases. Germ-line PTPN11 mutations cause ~50% of Noonan Syndrome (NS), which is among the most common autosomal dominant disorders [7, 8]. Gain-of-function mutations in other RAS-RAF-MEK-ERK pathway members, including SOS1[9, 10], KRAS, NRAS, SHOC2, and RAF1[14, 15], are responsible for most remaining NS cases. With an estimated incidence of 1/2,000 live births , NS is characterized by facial dysmorphism, proportional short stature, cardiac anomalies, and various less penetrant phenotypes, such as webbed neck, deafness, and motor delay. Many (20-50%) NS patients develop some type of myeloproliferative disorder (MPD), which is typically mild and self-limited . Rare NS patients progress to Juvenile Myelomonocytic Leukemia (JMML), which is fatal if not treated by bone marrow transplantation, somatic PTPN11 mutations are the single most common cause of sporadic JMML [7, 18–20]. LEOPARD Syndrome (LS), a much less common autosomal dominant disorder, is almost always caused by PTPN11 mutations, and is related to, but distinguishable from, NS [7, 16, 21]. LEOPARD is an acronym for its major syndromic manifestations: multiple Lentigines, Electrocardiographic abnormalities, Ocular hypertelorism, Pulmonary stenosis, Abnormalities of genitalia, Retardation of growth, and Deafness . These patients often have hypertrophic cardiomyopathy (HCM), and might also have an increased risk of neuroblastoma (NS) and acute myeloid leukemia (AML) [23, 24]. Knock-in mouse models have been generated for NS and LS alleles of Ptpn11 and generally reproduce the phenotypes seen in the cognate human syndromes [25–27].
Consistent with the distinct pathogenesis of NS and LS, different types of PTPN11 mutations cause these disorders. Most NS-associated PTPN11 mutations alter residues that reside at the interface between the N-SH2 and PTP domains , resulting in elevated enzymatic activity and enhanced RAS/ERK activation [27–31]. These data suggest that NS mutations disrupt the intramolecular interaction between the N-SH2 and PTP domains, shifting the equilibrium between the closed and open conformations and lowering the activation threshold for SHP2. By contrast, LS mutations typically affect PTP domain residues, result in markedly decreased catalytic activity, and lower RAS/ERK activation in transient transfection assays [31–33]. Studies of the LS Y279C mouse model also indicate that LS mutants may have dominant negative effects in at least some tissues in vivo. Whereas NS phenotypes arise from enhanced MEK/ERK activation and can be prevented or reversed by MEK inhibition [34–36], LS-associated HCM is caused by enhanced PI3K/AKT/mTORC1 activity and can be reversed by rapamycin .
Although multiple studies have reported the biochemical and biological consequences of NS- and LS-associated PTPN11 mutations, their structural consequences have not been analyzed. Here, we report the X-ray structures of five NS/LS SHP2 mutants and discuss how these mutations affect the interaction between different SHP2 domains and its catalytic activity.
A wild type (WT) SHP2 expression construct 1–539 (comprising the N + C SH2 and PTP domains) was PCR-amplified from PTPN11 cDNA  with a set of custom-designed primers (see Additional file 1: Description S1). The resultant PCR fragment was cloned into a modified version of the plasmid pET28b (Novagen) that generates a fusion protein with an N-terminal hexahistidine tag. The SHP2 catalytic domain expression construct (a.a. 221–524) was cloned into pGEX4T, which introduces a GST-tag at its N-terminus. Mutations were introduced into these expression constructs by site directed-mutagenesis with specifically designed primers bearing one substitution each (see Additional file 1: Description S1). Pfu Ultra II high fidelity DNA polymerase (Stratagene) was used for PCR, with an extension temperature of 68°C over 10 minutes. To remove any traces of the original cDNA, all reactions were subjected to digestion with DpnI (New England Biolabs) for 1 hour at 37°C. Reaction mixtures were transformed into DH5α cells, and the genetic content of all constructs was verified by Sanger sequencing.
Protein expression & purification
Vectors encoding full length versions of SHP2 mutants were transferred into E. coli BL21(DE3). Cells were grown in Terrific Broth containing kanamycin (50 mg/l) in 1 L Tunair flasks at 37°C to an OD600 of 3–5, after which the temperature was lowered to 16°C, and isopropyl-1-thio-β-D-galactopyranoside (IPTG) was added to 0.2 mM. Expression was allowed to proceed overnight, then cells were harvested by centrifugation, flash-frozen in liquid nitrogen, and stored at −80°C. Due to the low level of expression of the Q506P construct, cells expressing this mutant were washed using the osmotic shock technique  prior to freezing. Unless stated otherwise, all purification procedures were carried out at 4°C. Cells were thawed on ice and resuspended in Binding Buffer (see Additional file 2: Table S1 for detailed buffer components), supplemented with phenylmethylsulfonylfluoride and benzamidine. After disruption by sonication and centrifugation at 60,000 g for 40 min, cell-free extracts were passed through a DE-52 column (2.6 × 7 cm) that had been pre-equilibrated with the same buffer, and loaded by gravity flow onto a Ni-nitrilotriacetic acid (NTA) column (Qiagen, Germantown, MD). The latter column was washed with 20–25 volumes of Wash Buffer A, followed by 20–25 volumes of Wash Buffer B and finally with Elution Buffer. N308D and E139D eluted in Elution Buffer, whereas the other three mutants eluted with the Wash buffers. For the latter proteins, the wash fractions were diluted 15-fold and reloaded on fresh Ni-NTA columns. After washing with 10 column volumes of Binding Buffer, N308D and E139D proteins were eluted in elution buffer. These samples were concentrated using a VIVASpin unit (Sartorius NA, Edgewood, NY), and loaded onto a 2.6 × 60 cm Superdex 200 column (GE Healthcare), equilibrated with Gel Filtration buffer. Elution was performed at a flow rate of 3 ml/min at 8°C, with the SHP2 proteins behaving as apparent monomers. Final protein samples were concentrated to 20–40 mg/ml, divided into 1.5 mg aliquots, flash-frozen and stored at −80°C.
SHP2 catalytic domain mutants were transformed into Escherichia coli strain BL21(DE3). A 25 ml aliquot of an overnight culture from a single colony was added to 500 ml of LB/ampicillin (50 μg/ml) and grown at 37°C to A600 = 0.8. IPTG was added to a final concentration of 0.1 mM, and the bacteria were maintained for 16 h at 25°C with shaking, then centrifuged at 6,000 × g for 10 min. at 4°C. Pellets were resuspended in 12.5 ml of a buffer containing 50 mM Tris–HCl, pH 7.5, 150 mM NaCl, 5 mM MgCl2, 1% Triton X-100, 10% glycerol, 5 mM dithiothreitol, 2 μg/ml aprotinin, 10 μg/ml leupeptin, 1 μg/ml antipain, 1 μg/ml pepstatin A, 0.5 mg/ml lysozyme and 1 mg/mL DNase I. Suspensions were incubated on ice for 30 minutes, and then sonicated for 10 seconds on ice. Lysates were centrifuged at 14,000 × g for 30 min. at 4°C, and supernatants were transferred to a fresh 15-ml polypropylene tube containing 0.5 ml of glutathione-Sepharose 4B (GE Healthcare Life Sciences). This suspension was rotated end-over-end overnight at 4°C, and then centrifuged at 1000 × g for 1 min. at 4°C. The supernatants were discarded, and the beads were washed 3 times for 5 min. each at 4°C with 10 ml of wash buffer (25 mM Tris–HCl, pH 7.5, 150 mM NaCl, 5 mM MgCl2, 1% Triton X-100, 10% glycerol, 5 mM dithiothreitol), and then once with PTP assay buffer (25 mM Hepes, pH 7.5, 100 mM NaCl, 2 mM EDTA and 5 mM dithiothreitol). Bound GST fusion proteins were resuspended 1:1 in PTP assay buffer. A 20 uL aliquot of slurry for each mutant was separated on a 10% SDS-polyacrylamide gel, together with different amounts of BSA. The gel was washed in water for 10 minutes, and stained with Colloidal Coomassie Blue for 1 hour at room temperature. Bands were quantified using a LI-COR Odyssey.
To determine kinetic parameters, fixed amounts of purified GST-WT or -mutant SHP2 catalytic domains (1.6 pmol of WT and N308D, 115pmol of Y279C and 16.3 pmol of Q506P) were incubated with variable concentrations of substrate peptide (R-R-L-I-E-D-A-E-pY-A-A-R-G, Millipore #12-217; Kit #12-217) in PTP assay buffer in a total volume of 50 uL. Reactions were carried out for 10 minutes at 25°C, and phosphate release was quantified by adding Malachite Green (Millipore #17-125) to the supernatants, measuring absorbance at 620 nm, and comparing values to a standard curve generated with varying amounts of KH2PO4. All reactions fell within the linear range. Phosphatase activity is expressed in pmol Pi released/min/pmol enzyme.
Mutant SHP2 proteins were crystallized under conditions similar to those reported previously . In order to obtain the best diffracting crystals, 0.1 M LiCl was added to the literature crystallization conditions for D61G, 5% glycerol for N308D, and 10% glycerol and 0.3 M cycohexyl-methyl-β-D-maltoside for Q506P. The other two mutant proteins and the WT protein were crystallized under literature conditions with optimized precipitant concentrations. Crystals appeared overnight, and reached their full size of about 300 × 300 × 30 microns in one week at room temperature. The stacked plate crystals were separated and flash frozen in liquid nitrogen, using paratone-N oil (Hampton Research Inc.) as a cryo-protectant.
Data were collected at 100 K with a wavelength of 1.0 Å on the Industrial Macromolecular Crystallography Association (IMCA-CAT) beam line at the Advanced Photon Source (Argonne National Laboratory, IL USA). The data were indexed, integrated, and scaled with XDS and XSCALE.
The first mutant structure of N308D was determined by molecular replacement, using the previously solved structure of SHP2 (PDB access code 2SHP)  as the search model. The wild type SHP2 and D61G, E139D, Y279C and Q506P mutant structures were determined by molecular replacement, using the structure of the N308D SHP2 mutant as a search model. Following the initial rigid body refinement, interactive cycles of model building and refinement were performed by using COOT and Buster-TNT software. Special attention was paid to the mutation sites, which were initially replaced with alanine to reduce model bias and later positioned based on the 2mF o -DF c and difference Fourier electron density maps after a few rounds of refinement to confirm that those amino acids had, indeed, been mutated. Data collection and refinement statistics are shown in Table 1. The atomic coordinates have been deposited in the RCSB Protein Data Bank under accession numbers 4NXD, 4H10, 4NWG, 4GWF, 4NWF, and 4H34. All figures were produced using PyMOL (http://www.pymol.org).
Results and discussion
We determined the crystal structures of WT SHP2 (residues 1–539), as well as five mutants (D61G, E139D, Y279C, N308D, and Q506P), chosen to represent the spectrum of disease-associated PTPN11 mutations. Mutants D61G, E139D, and N308D are found in NS, Y279C is a canonical LS mutation [16, 21, 42], and Q506P has been reported in both disorders , although it is unclear whether this reflects misdiagnosis or true bi-potentiality of this allele. The D61G mutation affects the N-SH2 domain, E139D lies within the C-SH2 domain and the other three mutations alter the PTP domain (Figure 1b). The enzymatic properties of the full-length versions of these mutants (including the C-terminal tail, which is missing in our crystallization constructs) were characterized previously by our group [29, 33] (Additional file 3: Figure S1), and range from strongly activated (D61G), to mildly activated (N308D), to catalytically impaired (Y279C). Q506P shows altered specificity for some substrates .
The SHP2 structure published by Hof et al. (PDB accession code: 2SHP; hereafter termed “2SHP”) has three mutations (T2K, F41L and F513S) and is in complex with a detergent molecule (CTAB). We corrected these mutations, and crystallized WT SHP2 under detergent-free conditions. In our WT structure, the SH2 domains and the PTP domain assume a “closed” conformation with the N-SH2 domain locked into the PTP catalytic site, similar to the 2SHP structure (Figure 2). Superimposition of our bona fide WT structure with the earlier “WT” SHP2 structure revealed an overall root mean square deviation (rmsd) of 0.59 Å over 487 aligned residues. The N-SH2 domain had a smaller rmsd of 0.35 Å, whereas the rmsd for the C-SH2 was larger (0.74 Å), implying that the C-SH2 domain might have a higher degree of flexibility than the other two domains in the closed conformation. The three mutations and the CTAB binding site in the 2SHP structure are distant from the C-SH2 domain, so we do not think it likely that they affect C-SH2 domain flexibility directly. In the PTP domain, the major deviations between the two structures were seen in residues 425–436 of the αF helix (average rmsd =1.33 Å), corresponding to the CTAB binding site in the 2SHP structure (Figure 2). Although the overall structural difference between WT and 2SHP was small, we used our WT structure as the reference for comparison with the mutant structures to exclude any potential structural changes induced by the three mutations and the detergent molecule (CTAB) in 2SHP.
In the WT structure, Asp61 was located on the surface of the N-SH2 domain. The side chain of Asp61 formed hydrogen bonds with Ser460 from the catalytic P-loop (residues 458–464), a water-mediated hydrogen bond with the catalytic cysteinyl residue, Cys459, two water-mediated hydrogen bonds with Arg465, and another water-mediated hydrogen bond with Asp425 (Figure 3a). Consequently, Asp61 plays an important role in the N-SH2 and PTP domain interaction. In the D61G mutant structure, these hydrogen bonds were abolished. The change from aspartate to glycine also altered the surface charge from very negative to neutral (Figure 3b). Opposite D61G on the interface surface, the PTP domain presents a predominantly positively charged pocket (Figure 3c). Thus, the D61G mutation greatly loosened the interactions between N-SH2 and PTP domains. These data are consistent with previous publications that observed increased basal activity for this mutant against artificial (e.g., pNPP) and pTyr-peptide substrates [29, 33].
Residue Glu139 was located on the surface of the C-SH2 domain. The overall crystal structure of the E139D mutant was very similar to that of WT SHP2, with an rmsd of 0.4 Å (Figure 4a). Glu139 was about 40 Å away from the catalytic site, with the N-SH2 domain interposed between these domains. E139 is, however, located in the vicinity of the phosphate group of the pTyr peptide-binding site of the C-SH2 domain, and the E139D mutation stabilizes the conformation of the 139–147 loop that plays an essential role in pTyr peptide-binding. The mutant structure has well defined electron density for this loop, whereas this loop is disordered in the WT structure. Compared with WT SHP2, there also were some local structural rearrangements in the E139D mutant, with the most noticeable difference being a conformational change of His116 (Figure 4b). The E139D mutant shows a small increase in basal activity, but when it is assayed with a pTyr peptide (pTyr1172 from IRS-1) that can bind both SH2 domains and be dephosphorylated by the catalytic domain, the activity of this mutant is more than 5 times higher than that of WT SHP2 . Although the side chain charge remained unchanged when glutamic acid was changed to aspartic acid, the size of the side chain was reduced. This change could alter the surface of the adjacent residue, Arg138, a key residue for pTyr binding . As shown in Figure 4b, aspartic acid 139 (located in the βB strand) contributes less hydrogen bonding to residues His114-His116 of the βA strand than does glutamic acid. This subtle change might loosen the connection between the βA and βB strands, helping to expose the side chain of Arg138, and thus enhancing the affinity for pTyr-peptide binding. The binding of pTyr -peptides to the C-SH2 domain also could affect the interaction between the N-SH2 and C-SH2 that is critical for enzymatic activation . The mutant also could facilitate the binding of the C-SH2 domain to certain physiologically important binding partners; e.g., IRS-1, in which pTyr-1222 binds to the C-SH2 domain, while pTyr-1172 binds to the N-SH2 domain. The E139D mutation could indirectly increase the binding affinity of the N-SH2 domain for SHP2 substrates (for the reasons discussed above), and therefore increase catalytic activity. In concert, these effects likely explain why E139D is activated by pTyr-peptide binding more than WT SHP2.
Y279C is a catalytically impaired mutant associated with LS [16, 21, 42]. In the crystal structure of this mutant, Tyr279 was located in the long pTyr loop (residues 277–288) with two small α-helices, αC (residues 265–269) and αD (residues 271–276), at its upstream end and one small β-strand, βB (residues 289–292), at its downstream end. In general, because the intervening region (residues 262–288) lacks a structurally stable long α-helix or β-strand, it is likely to be flexible. It also contained a large number of positively charged side chains pointing toward the surface of the PTP domain that might interact with solvent molecules or other binding partners and thus increase the mobility of the pTyr loop. In the WT structure, the side chain of Tyr279 makes van der Waals contacts with Ser460 and Ala461 of the catalytic P-loop, as well as with Tyr62 and Lys70 from the N-SH2 domain (Figure 5). The –OH group of Tyr62 interacts with the π-electrons of the Tyr279’s aromatic ring. Together with Q506 in the Q-loop (residues 501–507), Tyr279 is believed to play a key role in binding the tyrosine side chains of substrate proteins/peptides during catalysis. In the PTP1B structure , Tyr46 is the residue equivalent to Tyr279 in the SHP2 structure [3, 45]. The stacking interaction of three residues, Tyr46, the substrate pTyr, and Phe182 from the “WPDF” loop, help to properly position the substrate for catalysis. In the Y279C structure, the interactions of Tyr279 with the P-loop and the N-SH2 domain were disrupted, due to the significantly shorter side chain of cysteine compared with that of tyrosine. SHP2 has a WPDH loop that corresponds to WPDF in PTP1B. The Y279C mutation would have less stacking interactions with Tyr279, a bound pTyr substrate, and His426 (the equivalent of Phe182 in PTP1B). The mutation also distorts the pTyr substrate/SHP2 interaction, and thus would be expected to disrupt catalysis significantly. At the same time, the Y279C mutation results in loss of the Tyr279/Tyr62 interaction, diminishing the strength of intramolecular binding between the PTP and N-SH2 domains. This would be expected to facilitate the “open” conformation, and can explain the enhanced interaction of this mutant with binding partners (e.g., GAB1) observed previously [33, 45, 46]. While this manuscript was in preparation, Yu et al.  reported a very similar model of Y279C (PDB accession code: 4DGX), as well as WT SHP2 (PDB accession code: 4DGP). The Cα carbon atom comparison for these structures revealed RMSD of 0.42 and 0.47 Å for the corresponding mutant and wild-type crystal structure pairs correspondingly. Both published crystal structures belong to the P21212 space group, which is different from the space group (P21) for our Y279C mutant and WT structures. Importantly, the residues surrounding the mutant residue are in a similar conformation. Our WT structure has two disordered regions (perhaps due to their flexibility) lacking electron density (88–95 and 156–164). From both Y279C structures, we can see that replacing tyrosine with cysteine at position 279 does not block the accessibility to the PTP catalytic site; instead, it could facilitate local conformational changes that lead to the release of the N-SH2 domain, and thereby open the conformation of the substrate-binding site.
To test the predictions of these structural studies, we assayed the enzymatic activity of isolated SHP2 catalytic domains (present as part of GST-fusion proteins). Compared with WT GST-SHP2, the Y279C mutant had ~80 times lower kcat and 7 times higher Km towards a pY peptide (Figure 5, Table 2). These results are quite similar to those were reported by Yu et al.. Yu et al. also provided hydrogen/deuterium exchange mass spectrometry experiments and molecular dynamics simulations showing that the N-SH2 and PTP domain interaction was decreased in the Y279C mutant. Moreover, they found that the Y279C mutant displayed higher affinity for, and was preferentially activated by, a non-hydrolyzable N-SH2 ligand.
Residue 308 was located in the βC strand of the PTP domain, and was not involved in direct interactions with the N-SH2 domain. However, the Oδ1 atom from the side chain of Asn308 formed a strong hydrogen bond with the side chain of the conserved Arg501. Arg 501 also made direct hydrogen bonds with the main chains of the P-loop residues Ala461 and Gly462 (Figure 6a). The Nδ2 atom of Asn308 formed two hydrogen bonds with the main chain of Phe285 from the pTyr-loop. In the N308D mutant, the charge of the side chain changes from neutral to negative, whereas the side chain polarity changes from polar to acidic polar. Compared with the WT Asn residue, Asp308 formed more hydrogen bonds with pTyr-loop residues, and an especially strong one (2.5 Å) with Thr288 (Figure 6b). Consequently, this mutation could make the pTyr- and P-loops less flexible, locking the enzyme in a more favorable position for catalysis. The greater rigidity of the pTyr and P loops makes it more difficult for the N-SH2 domain to close back on the PTP domain once it opens. Since the open and closed forms are in equilibrium, this could mean that the ability (i.e., rate constant) for closing back is significantly diminished, hence favoring the open form, leading to the increased activity (basally and in response to pTyr peptide) of the N308D mutant as a full-length enzyme. The residue Asn308 is a hot spot for NS mutations, with N308D accounting for 25% of NS cases. Previous studies showed that this mutant (as a full length protein) has a 3-fold higher basal activity than WT . Our PTP assay showed that catalytic domain of this mutation had slightly higher kcat and Km values when compared with the WT PTP domain (Figure 1, Table 2).
Residue 506 was located in the interface between the PTP and N-SH2 domains. In the WT structure, Q506 formed two important hydrogen bonds, with the main chain of Ala72 and with the side chain of Asn58 in the N-SH2 domain. Together with Tyr279, Gln506 also plays an important role in PTP catalysis by binding the tyrosine side chain of the substrate  and by helping to properly position a water molecule for hydrolysis of the thiophosphate intermediate (Figure 7). In the Q506P mutant structure, the proline mutation abolished the two hydrogen bonds between the PTP and N-SH2 domain, which predicts that this mutant should be more “open” than WT SHP2. However, this mutation also disrupts the C459-H2O-Q506 interaction. As a result, the water molecule needed in the second step of catalysis probably cannot be positioned properly. Consistent with this notion, basal PTP activity and PTP activity in the presence of an N-SH2 domain binding pTyr peptide are slightly lower (when measured against the artificial substrate pNPP) than in WT SHP2 . We also monitored the activity of the isolated catalytic domain of the Q506P mutant. In accord with our structural data, the kcat of Q506P was ~15–fold lower and the Km about 2× higher than in WT SHP2 (Table 2).
SHP2 is regulated by a molecular switch mechanism that controls its catalytic activity. Upon binding to a tyrosine-phosphorylated binding partner for its SH2 domains, the N-terminal SH2 domain is released from the PTP domain, activating the enzyme. This elegant mechanism ensures that PTP activity is delivered to the right place in the cell at the right time. Remarkably, germ line mutations that disrupt this regulatory machinery in different ways result in distinct disease syndromes. The crystal structures of “true” WT SHP2 and five NS/LS-associated SHP2 mutants reported herein provide direct comparisons of the local conformational changes caused by each mutation. Our structural observations are in agreement with, and can provide mechanistic insight into, the previously reported catalytic properties of these mutants. For example, mutation of D61G in the N-SH2 domain significantly impacts SHP2 activity because this residue is located at the N-SH2/PTP domain interface and its alteration weakens key interactions between the two domains. On the other hand, our data suggest that the C-SH2 domain mutation E139D might interfere with SHP2 binding to tryrosine-phosphorylated ligands. The other three mutants, Y279C, N308D and Q506P, are located in PTP domain, and the local conformational changes induced by each mutation provide insight into their abnormal catalytic properties. The results of our research provide structural insights into this medically important target and could aid in future structure-based drug discovery programs.
Availability of supporting data
The coordinates and diffraction data for SHP2 wild type and mutant crystal structures are available in Protein Data Bank (http://www.rcsb.org/pdb).
Neel BG, Gu H, Pao L: The ‘Shp’ing news: SH2 domain-containing tyrosine phosphatases in cell signaling. Trends Biochem Sci 2003, 28(6):284–293.
Feng GS: Shp-2 tyrosine phosphatase: signaling one cell or many. Exp Cell Res 1999, 253(1):47–54.
Hof P, Pluskey S, Dhe-Paganon S, Eck MJ, Shoelson SE: Crystal structure of the tyrosine phosphatase SHP-2. Cell 1998, 92(4):441–450.
Tonks NK, Neel BG: Combinatorial control of the specificity of protein tyrosine phosphatases. Curr Opin Cell Biol 2001, 13(2):182–195.
Barford D, Neel BG: Revealing mechanisms for SH2 domain mediated regulation of the protein tyrosine phosphatase SHP-2. Structure 1998, 6(3):249–254.
O’Reilly AM, Neel BG: Structural determinants of SHP-2 function and specificity in Xenopus mesoderm induction. Mol Cell Biol 1998, 18(1):161–177.
Tartaglia M, Gelb BD: Germ-line and somatic PTPN11 mutations in human disease. Eur J Med Genet 2005, 48(2):81–96.
Tartaglia M, Mehler EL, Goldberg R, Zampino G, Brunner HG, Kremer H, van der Burgt I, Crosby AH, Ion A, Jeffery S, Kalidas K, Patton MA, Kucherlapati RS, Gelb BD: Mutations in PTPN11, encoding the protein tyrosine phosphatase SHP-2, cause Noonan syndrome. Nature genetics 2001, 29(4):465–468.
Tartaglia M, Pennacchio LA, Zhao C, Yadav KK, Fodale V, Sarkozy A, Pandit B, Oishi K, Martinelli S, Schackwitz W, Ustaszewska A, Martin J, Bristow J, Carta C, Lepri F, Neri C, Vasta I, Gibson K, Curry CJ, Siguero JP, Digilio MC, Zampino G, Dallapiccola B, Bar-Sagi D, Gelb BD: Gain-of-function SOS1 mutations cause a distinctive form of Noonan syndrome. Nat Genet 2007, 39(1):75–79.
Roberts AE, Araki T, Swanson KD, Montgomery KT, Schiripo TA, Joshi VA, Li L, Yassin Y, Tamburino AM, Neel BG, Kucherlapati RS: Germline gain-of-function mutations in SOS1 cause Noonan syndrome. Nat Genet 2007, 39(1):70–74.
Schubbert S, Zenker M, Rowe SL, Boll S, Klein C, Bollag G, van der Burgt I, Musante L, Kalscheuer V, Wehner LE, Nguyen H, West B, Zhang KY, Sistermans E, Rauch A, Niemeyer CM, Shannon K, Kratz CP: Germline KRAS mutations cause Noonan syndrome. Nat Genet 2006, 38(3):331–336.
Cirstea IC, Kutsche K, Dvorsky R, Gremer L, Carta C, Horn D, Roberts AE, Lepri F, Merbitz-Zahradnik T, Konig R, Kratz CP, Pantaleoni F, Dentici ML, Joshi VA, Kucherlapati RS, Mazzanti L, Mundlos S, Patton MA, Silengo MC, Rossi C, Zampino G, Digilio C, Stuppia L, Seemanova E, Pennacchio LA, Gelb BD, Dallapiccola B, Wittinghofer A, Ahmadian MR, Tartaglia M, Zenker M: A restricted spectrum of NRAS mutations causes Noonan syndrome. Nat Genet 2010, 42(1):27–29.
Cordeddu V, Di Schiavi E, Pennacchio LA, Ma'ayan A, Sarkozy A, Fodale V, Cecchetti S, Cardinale A, Martin J, Schackwitz W, Lipzen A, Zampino G, Mazzanti L, Digilio MC, Martinelli S, Flex E, Lepri F, Bartholdi D, Kutsche K, Ferrero GB, Anichini C, Selicorni A, Rossi C, Tenconi R, Zenker M, Merlo D, Dallapiccola B, Iyengar R, Bazzicalupo P, Gelb BD, Tartaglia M: Mutation of SHOC2 promotes aberrant protein N-myristoylation and causes Noonan-like syndrome with loose anagen hair. Nat Genet 2009, 41(9):1022–1026.
Pandit B, Sarkozy A, Pennacchio LA, Carta C, Oishi K, Martinelli S, Pogna EA, Schackwitz W, Ustaszewska A, Landstrom A, Bos JM, Ommen SR, Esposito G, Lepri F, Faul C, Mundel P, López Siguero JP, Tenconi R, Selicorni A, Rossi C, Mazzanti L, Torrente I, Marino B, Digilio MC, Zampino G, Ackerman MJ, Dallapiccola B, Tartaglia M, Gelb BD: Gain-of-function RAF1 mutations cause Noonan and LEOPARD syndromes with hypertrophic cardiomyopathy. Nat Genet 2007, 39(8):1007–1012.
Razzaque MA, Nishizawa T, Komoike Y, Yagi H, Furutani M, Amo R, Kamisago M, Momma K, Katayama H, Nakagawa M, Fujiwara Y, Matsushima M, Mizuno K, Tokuyama M, Hirota H, Muneuchi J, Higashinakagawa T, Matsuoka R: Germline gain-of-function mutations in RAF1 cause Noonan syndrome. Nat Genet 2007, 39(8):1013–1017.
Tartaglia M, Gelb BD: Noonan syndrome and related disorders: genetics and pathogenesis. Annu Rev Genomics Hum Genet 2005, 6: 45–68.
Bader-Meunier B, Tchernia G, Mielot F, Fontaine JL, Thomas C, Lyonnet S, Lavergne JM, Dommergues JP: Occurrence of myeloproliferative disorder in patients with Noonan syndrome. J Pediatr 1997, 130(6):885–889.
Tartaglia M, Niemeyer CM, Shannon KM, Loh ML: SHP-2 and myeloid malignancies. Curr Opin Hematol 2004, 11(1):44–50.
Loh ML, Vattikuti S, Schubbert S, Reynolds MG, Carlson E, Lieuw KH, Cheng JW, Lee CM, Stokoe D, Bonifas JM, Curtiss NP, Gotlib J, Meshinchi S, Le Beau MM, Emanuel PD, Shannon KM: Mutations in PTPN11 implicate the SHP-2 phosphatase in leukemogenesis. Blood 2004, 103(6):2325–2331.
Chan G, Kalaitzidis D, Neel BG: The tyrosine phosphatase Shp2 (PTPN11) in cancer. Cancer Metastasis Rev 2008, 27(2):179–192.
Legius E, Schrander-Stumpel C, Schollen E, Pulles-Heintzberger C, Gewillig M, Fryns JP: PTPN11 mutations in LEOPARD syndrome. J Med Genet 2002, 39(8):571–574.
Gorlin RJ, Anderson RC, Moller JH: The leopard (multiple lentigines) syndrome revisited. Laryngoscope 1971, 81(10):1674–1681.
Merks JH, Caron HN, Hennekam RC: High incidence of malformation syndromes in a series of 1,073 children with cancer. Am J Med Genet A 2005, 134A(2):132–143.
Ucar C, Calyskan U, Martini S, Heinritz W: Acute myelomonocytic leukemia in a boy with LEOPARD syndrome (PTPN11 gene mutation positive). J Pediatr Hematol Oncol 2006, 28(3):123–125.
Marin TM, Keith K, Davies B, Conner DA, Guha P, Kalaitzidis D, Wu X, Lauriol J, Wang B, Bauer M, Bronson R, Franchini KG, Neel BG, Kontaridis MI: Rapamycin reverses hypertrophic cardiomyopathy in a mouse model of LEOPARD syndrome-associated PTPN11 mutation. The Journal of clinical investigation 2011, 121(3):1026–1043.
Krenz M, Yutzey KE, Robbins J: Noonan syndrome mutation Q79R in Shp2 increases proliferation of valve primordia mesenchymal cells via extracellular signal-regulated kinase 1/2 signaling. Circ Res 2005, 97(8):813–820.
Araki T, Mohi MG, Ismat FA, Bronson RT, Williams IR, Kutok JL, Yang W, Pao LI, Gilliland DG, Epstein JA, Neel BG: Mouse model of Noonan syndrome reveals cell type- and gene dosage-dependent effects of Ptpn11 mutation. Nat Med 2004, 10(8):849–857.
Fragale A, Tartaglia M, Wu J, Gelb BD: Noonan syndrome-associated SHP2/PTPN11 mutants cause EGF-dependent prolonged GAB1 binding and sustained ERK2/MAPK1 activation. Hum Mutat 2004, 23(3):267–277.
Keilhack H, David FS, McGregor M, Cantley LC, Neel BG: Diverse biochemical properties of Shp2 mutants. Implications for disease phenotypes. J Biol Chem 2005, 280(35):30984–30993.
Niihori T, Aoki Y, Ohashi H, Kurosawa K, Kondoh T, Ishikiriyama S, Kawame H, Kamasaki H, Yamanaka T, Takada F, Nishio K, Sakurai M, Tamai H, Nagashima T, Suzuki Y, Kure S, Fujii K, Imaizumi M, Matsubara Y: Functional analysis of PTPN11/SHP-2 mutants identified in Noonan syndrome and childhood leukemia. J Hum Genet 2005, 50(4):192–202.
Tartaglia M, Martinelli S, Stella L, Bocchinfuso G, Flex E, Cordeddu V, Zampino G, Burgt I, Palleschi A, Petrucci TC, Alsorcini M, Schoch C, Foa R, Emanuel PD, Gelb BD: Diversity and functional consequences of germline and somatic PTPN11 mutations in human disease. Am J Hum Genet 2006, 78(2):279–290.
Hanna N, Montagner A, Lee WH, Miteva M, Vidal M, Vidaud M, Parfait B, Raynal P: Reduced phosphatase activity of SHP-2 in LEOPARD syndrome: consequences for PI3K binding on Gab1. FEBS Lett 2006, 580(10):2477–2482.
Kontaridis MI, Swanson KD, David FS, Barford D, Neel BG: PTPN11 (Shp2) mutations in LEOPARD syndrome have dominant negative, not activating, effects. J Biol Chem 2006, 281(10):6785–6792.
Krenz M, Gulick J, Osinska HE, Colbert MC, Molkentin JD, Robbins J: Role of ERK1/2 signaling in congenital valve malformations in Noonan syndrome. Proc Natl Acad Sci U S A 2008, 105(48):18930–18935.
Wu X, Simpson J, Hong JH, Kim KH, Thavarajah NK, Backx PH, Neel BG, Araki T: MEK-ERK pathway modulation ameliorates disease phenotypes in a mouse model of Noonan syndrome associated with the Raf1(L613V) mutation. J Clin Invest 2011, 121(3):1009–1025.
Lauriol J, Kontaridis MI: PTPN11-associated mutations in the heart: has LEOPARD changed Its RASpots? Trends Cardiovasc Med 2011, 21(4):97–104.
Freeman RM Jr, Plutzky J, Neel BG: Identification of a human src homology 2-containing protein-tyrosine-phosphatase: a putative homolog of Drosophila corkscrew. Proc Natl Acad Sci U S A 1992, 89(23):11239–11243.
Magnusdottir A, Johansson I, Dahlgren LG, Nordlund P, Berglund H: Enabling IMAC purification of low abundance recombinant proteins from E. coli lysates. Nature methods 2009, 6(7):477–478.
Kabsch W: Integration, scaling, space-group assignment and post-refinement. Acta crystallographica Section D, Biol Crystallogr 2010, 66(Pt 2):133–144.
Emsley P, Lohkamp B, Scott W, Cowtan K: Features and development of Coot. Acta Crystallogr D Biol Crystallogr 2010, 66(4):486–501.
Bricogne GBE, Brandl M, Flensburg C, Keller P, Paciorek W, Roversi P, Sharff A, Smart OS, Vonrhein C, Womack TO: Buster 2.10.0. In 2.10.0 edn. Cambridge: United Kingdom: Global Phasing Ltd; 2011.
Digilio MC, Conti E, Sarkozy A, Mingarelli R, Dottorini T, Marino B, Pizzuti A, Dallapiccola B: Grouping of multiple-lentigines/LEOPARD and Noonan syndromes on the PTPN11 gene. Am J Hum Genet 2002, 71(2):389–394.
Eck MJ, Pluskey S, Trub T, Harrison SC, Shoelson SE: Spatial constraints on the recognition of phosphoproteins by the tandem SH2 domains of the phosphatase SH-PTP2. Nature 1996, 379(6562):277–280.
Tonks NK: PTP1B: from the sidelines to the front lines! FEBS Lett 2003, 546(1):140–148.
Edouard T, Combier JP, Nedelec A, Bel-Vialar S, Metrich M, Conte-Auriol F, Lyonnet S, Parfait B, Tauber M, Salles JP, Lezoualc'h F, Yart A, Raynal P: Functional effects of PTPN11 (SHP2) mutations causing LEOPARD syndrome on epidermal growth factor-induced phosphoinositide 3-kinase/AKT/glycogen synthase kinase 3beta signaling. Mol Cell Biol 2010, 30(10):2498–2507.
Yu ZH, Xu J, Walls CD, Chen L, Zhang S, Zhang R, Wu L, Wang L, Liu S, Zhang ZY: Structural and mechanistic insights into LEOPARD syndrome-associated SHP2 mutations. J Biol Chem 2013, 288(15):10472–10482.
These studies were supported by the Ontario Research and Development Challenge Fund (99-SEP-0512) and R37 CA49152 (to BGN). BGN and EFP are Canada Research Chairs, Tier 1, and work in their laboratories is partially supported by the Ontario Ministry of Health and Long Term Care and The Princess Margaret Cancer Foundation. Use of the IMCA-CAT beamline 17-ID at the Advanced Photon Source was supported by the companies of the Industrial Macromolecular Crystallography Association through a contract with the Hauptman-Woodward Medical Research Institute. Use of the Advanced Photon Source was supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, under Contract No. DE-AC02-06CH11357. We thank Aiping Dong for providing technical support and the Structural Genomics Consortium, University of Toronto, for the use of their X-ray facilities.
The authors declare that they have no competing interests.
QW carried out crystal structure determination, crystallographic refinement, structure analysis, and drafted the manuscript. WX performed assay studies and participated in the drafting of the manuscript. RV and HA performed protein sample preparation and structure analysis. LA carried out crystallization studies. RM performed construct design and cloning experiments for the mutants and wild type structures. BKP performed X-ray diffraction synchrotron experiments and participated in structure analysis. EFP performed structure analysis and participated in drafting the manuscript. BGN, NYC conceived of the study, participated in its design and coordination, and drafted the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 3: Figure S1: Activities of full-length WT and mutant SHP2 studied in this manuscript (from references [29, 33]). The in vitro catalytic activities of the indicated GST-SHP2-FLAG proteins were measured using the artificial substrate 32P-labeled reduced carboxamido-methylated and –maleylated lysozyme (32P-RCML) in the absence or presence of an insulin receptor substrate-1-derived peptide containing phospho-tyrosine-1172 (pY1172) (100 μM). The pY1172 peptide binds to the N-SH2 domain, which in turn “opens up” the enzyme. (TIFF 1 MB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.