Discriminating between competing models for the allosteric regulation of oncogenic phosphatase SHP2 by characterizing its active state

Phosphopeptide binding to the N-SH2 domain of tyrosine phosphatase SHP2 correlates with the unzipping of its central β-sheet

SHP2 is a tyrosine phosphatase that plays a regulatory role in multiple intracellular signaling cascades and is known to be oncogenic in certain contexts. In the absence of effectors, SHP2 adopts an autoinhibited conformation with its N-SH2 domain blocking the active site. Given the key role of N-SH2 in regulating SHP2, this domain has been extensively studied, often by X-ray crystallography. Using a combination of structural analyses and molecular dynamics (MD) simulations we show that the crystallographic environment can significantly influence the structure of the isolated N-SH2 domain, resulting in misleading interpretations. As an orthogonal method to X-ray crystallography, we use a combination of NMR spectroscopy and MD simulations to accurately determine the conformation of apo N-SH2 in solution. In contrast to earlier reports based on crystallographic data, our results indicate that apo N-SH2 in solution primarily adopts a conformation with a fully zipped central β-sheet, and that partial unzipping of this β-sheet is promoted by binding of either phosphopeptides or even phosphate/sulfate ions.

Molecular recognition of ITIM/ITSM domains with SHP2 and their allosteric effect

Src homology 2-domain-containing tyrosine phosphatase 2 (SHP2) is a non-receptor protein tyrosine phosphatase that is widely expressed in a variety of cells and regulates the immune response of T cells through the PD-1 pathway. However, the activation mechanism and allosteric effects of SHP2 remain unclear, hindering the development of small molecule inhibitors. For the first time, in this study, the complex structure formed by the intact PD-1 tail and SHP2 was modeled. The molecular recognition and conformational changes of inactive/active SHP2 versus ITIM/ITSM were compared based on prolonged MD simulations. The relative flexibility of the two SH2 domains during MD simulations contributes to the recruitment of ITIM/ITSM and supports the subsequent conformational change of SHP2. The binding free energy calculation shows that inactive SHP2 has a higher affinity for ITIM/ITSM than active SHP2, mainly because the former's N-SH2 refers to the α-state. In addition, a significant decrease in the contribution to the binding energy of certain residues (e.g., R32, S34, K35, T42, and K55) of conformationally transformed SHP2 contributes to the above result. These detailed changes during conformational transition will provide theoretical guidance for the molecular design of subsequent novel anticancer drugs.

Single Ion Pair Is Essential for Stabilizing SHP2's Open Conformation

Src-homology-2-domain-containing PTP-2 (SHP2) is a widely expressed signaling enzyme whose misregulation is associated with multiple human pathologies. SHP2's enzymatic activity is controlled by a conformational equilibrium between its autoinhibited ("closed") state and its activated ("open") state. Although SHP2's closed state has been extensively characterized, the putative structure of its open form has only been revealed in the context of a highly activated mutant (E76K), and no systematic studies of the biochemical determinants of SHP2's open-state stabilization have been reported. To identify amino-acid interactions that are critical for stabilizing SHP2's active state, we carried out a mutagenic study of residues that lie at potentially important interdomain interfaces of the open conformation. The open/closed equilibria of the mutants were evaluated, and we identified several interactions that contribute to the stabilization of SHP2's open state. In particular, our findings establish that an ion pair between glutamate 249 on SHP2's PTP domain and arginine 111 on an interdomain loop is the key determinant of SHP2's open-state stabilization. Mutations that disrupt the R111/E249 ion pair substantially shift SHP2's open/closed equilibrium to the closed state, even compared to wild-type SHP2's basal-state equilibrium, which strongly favors the closed state. To the best of our knowledge, the ion-pair variants uncovered in this study are the first known SHP2 mutants in which autoinhibition is augmented with respect to the wild-type protein. Such "hyperinhibited" mutants may provide useful tools for signaling studies that investigate the connections between SHP2 inhibition and the suppression of human disease progression.

Atomistic ensemble of active SHP2 phosphatase

SHP2 phosphatase plays an important role in regulating several intracellular signaling pathways. Pathogenic mutations of SHP2 cause developmental disorders and are linked to hematological malignancies and cancer. SHP2 comprises two tandemly-arranged SH2 domains, a catalytic PTP domain, and a disordered C-terminal tail. Under physiological, non-stimulating conditions, the catalytic site of PTP is occluded by the N-SH2 domain, so that the basal activity of SHP2 is low. Whereas the autoinhibited structure of SHP2 has been known for two decades, its active, open structure still represents a conundrum. Since the oncogenic mutant SHP2E76K almost completely populates the active, open state, this mutant has been extensively studied as a model for activated SHP2. By molecular dynamics simulations and accurate explicit-solvent SAXS curve predictions, we present the heterogeneous atomistic ensemble of constitutively active SHP2E76K in solution, encompassing a set of conformational arrangements and radii of gyration in agreement with experimental SAXS data.

SHP2 clinical phenotype, cancer, or RASopathies, can be predicted by mutant conformational propensities

SHP2 phosphatase promotes full activation of the RTK-dependent Ras/MAPK pathway. Its mutations can drive cancer and RASopathies, a group of neurodevelopmental disorders (NDDs). Here we ask how same residue mutations in SHP2 can lead to both cancer and NDD phenotypes, and whether we can predict what the outcome will be. We collected and analyzed mutation data from the literature and cancer databases and performed molecular dynamics simulations of SHP2 mutants. We show that both cancer and Noonan syndrome (NS, a RASopathy) mutations favor catalysis-prone conformations. As to cancer versus RASopathies, we demonstrate that cancer mutations are more likely to accelerate SHP2 activation than the NS mutations at the same genomic loci, in line with NMR data for K-Ras4B more aggressive mutations. The compiled experimental data and dynamic features of SHP2 mutants lead us to propose that different from strong oncogenic mutations, SHP2 activation by NS mutations is less likely to induce a transition of the ensemble from the SHP2 inactive state to the active state. Strong signaling promotes cell proliferation, a hallmark of cancer. Weak, or moderate signals are associated with differentiation. In embryonic neural cells, dysregulated differentiation is connected to NDDs. Our innovative work offers structural guidelines for identifying and correlating mutations with clinical outcomes, and an explanation for why bearers of RASopathy mutations may have a higher probability of cancer. Finally, we propose a drug strategy against SHP2 variants-promoting cancer and RASopathies.

Quantitation of autoinhibitory defects in pathogenic SHP2 mutants by differential scanning fluorimetry

Src-homology-2-domain-containing protein tyrosine phosphatase-2 (SHP2) is a signaling enzyme whose activity is governed by an equilibrium between autoinhibited and activated states. Regulation of SHP2 activity is critical for cellular homeostasis, and mutations that alter its autoregulatory equilibrium cause cancers and developmental disorders. Several methods for assessing the strength of autoinhibitory interactions in SHP2 mutants have been previously reported, but each has limitations. We show that differential scanning fluorimetry provides a rapid, quantitative measure of SHP2 autoinhibition that is independent of the intrinsic activity of the SHP2 mutant being analyzed, does not involve protein labeling, and does not require specialized instrumentation.

Study on the allosteric activation mechanism of SHP2 via elastic network models and neural relational inference molecular dynamics simulation

As a ubiquitous protein tyrosine phosphatase, SHP2 is involved in PD-1/PD-L1 mediated tumor immune escape and undergoes substantial conformational changes. Therefore, it is considered an ideal target for tumor intervention. However, the allosteric mechanisms of SHP2 binding PD-1 intracellular ITIM/ITSM phosphopeptides remain unclear, which greatly hinders the development of novel structure-based anticancer allosteric inhibitors. In this work, the open and closed structural models of SHP2 are first constructed based on this knowledge; next their motion modes are investigated via elastic network models such as the Gaussian network model (GNM), anisotropic network model (ANM) and adaptive anisotropic network model (aANM); and finally, a possible allosteric signaling pathway is proposed using a neural relational inference molecular dynamics (NRI-MD) simulation embedded with an artificial intelligence (AI) strategy. In GNM and ANM, the N-SH2, C-SH2 and PTP domains all exhibit distinct dynamics partitions, and the N-SH2/C-SH2 regions show a rigid rotation relative to PTP. According to a series of intermediate snapshots given by aANM, N-SH2 is first identified with pY223 specifically, inducing a D'E-loop to change from β-sheets to random coils, and then, C-SH2 serves as a fulcrum to drive N-SH2 to rotate 110° completely away from the original active sites of PTP. Finally, a possible allosteric signaling-transfer path for SHP2, namely R220-R138-T108-R32, is proposed based on NRI-MD sampling. This work provides a possible allosteric mechanism of SHP2, which is helpful for the following design of novel allosteric inhibitors and is expected to be used in clinical synergies with PD-1 monoclonal antibody.

SH2 Domains: Folding, Binding and Therapeutical Approaches

SH2 (Src Homology 2) domains are among the best characterized and most studied protein-protein interaction (PPIs) modules able to bind and recognize sequences presenting a phosphorylated tyrosine. This post-translational modification is a key regulator of a plethora of physiological and molecular pathways in the eukaryotic cell, so SH2 domains possess a fundamental role in cell signaling. Consequently, several pathologies arise from the dysregulation of such SH2-domains mediated PPIs. In this review, we recapitulate the current knowledge about the structural, folding stability, and binding properties of SH2 domains and their roles in molecular pathways and pathogenesis. Moreover, we focus attention on the different strategies employed to modulate/inhibit SH2 domains binding. Altogether, the information gathered points to evidence that pharmacological interest in SH2 domains is highly strategic to developing new therapeutics. Moreover, a deeper understanding of the molecular determinants of the thermodynamic stability as well as of the binding properties of SH2 domains appears to be fundamental in order to improve the possibility of preventing their dysregulated interactions.

Allosteric regulation of autoinhibition and activation of c-Abl

c-Abl, a non-receptor tyrosine kinase, regulates cell growth and survival in healthy cells and causes chronic myeloid leukemia (CML) when fused by Bcr. Its activity is blocked in the assembled inactive state, where the SH3 and SH2 domains dock into the kinase domain, reducing its conformational flexibility, resulting in the autoinhibited state. It is active in an extended ‘open’ conformation. Allostery governs the transitions between the autoinhibited and active states. Even though experiments revealed the structural hallmarks of the two states, a detailed grasp of the determinants of c-Abl autoinhibition and activation at the atomic level, which may help innovative drug discovery, is still lacking. Here, using extensive molecular dynamics simulations, we decipher exactly how these determinants regulate it. Our simulations confirm and extend experimental data that the myristoyl group serves as the switch for c-Abl inhibition/activation. Its dissociation from the kinase domain promotes the SH2-SH3 release, initiating c-Abl activation. We show that the precise SH2/N-lobe interaction is required for full activation of c-Abl. It stabilizes a catalysis-favored conformation, priming it for catalytic action. Bcr-Abl allosteric drugs elegantly mimic the endogenous myristoyl-mediated autoinhibition state of c-Abl 1b. Allosteric activating mutations shift the ensemble to the active state, blocking ATP-competitive drugs. Allosteric drugs alter the active-site conformation, shifting the ensemble to re-favor ATP-competitive drugs. Our work provides a complete mechanism of c-Abl activation and insights into critical parameters controlling at the atomic level c-Abl inactivation, leading us to propose possible strategies to counter reemergence of drug resistance.