Phosphotyrosine couples peptide binding and SHP2 activation via a dynamic allosteric network

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.

The pathogenic T42A mutation in SHP2 rewires the interaction specificity of its N-terminal regulatory domain

Mutations in the tyrosine phosphatase Src homology-2 domain-containing protein tyrosine phosphatase-2 (SHP2) are associated with a variety of human diseases. Most mutations in SHP2 increase its basal catalytic activity by disrupting autoinhibitory interactions between its phosphatase domain and N-terminal SH2 (phosphotyrosine recognition) domain. By contrast, some disease-associated mutations located in the ligand-binding pockets of the N- or C-terminal SH2 domains do not increase basal activity and likely exert their pathogenicity through alternative mechanisms. We lack a molecular understanding of how these SH2 mutations impact SHP2 structure, activity, and signaling. Here, we characterize five SHP2 SH2 domain ligand-binding pocket mutants through a combination of high-throughput biochemical screens, biophysical and biochemical measurements, and molecular dynamics simulations. We show that while some of these mutations alter binding affinity to phosphorylation sites, the T42A mutation in the N-SH2 domain is unique in that it also substantially alters ligand-binding specificity, despite being 8 to 10 Å from the specificity-determining region of the SH2 domain. This mutation exerts its effect on sequence specificity by remodeling the phosphotyrosine-binding pocket, altering the mode of engagement of both the phosphotyrosine and surrounding residues on the ligand. The functional consequence of this altered specificity is that the T42A mutant has biased sensitivity toward a subset of activating ligands and enhances downstream signaling. Our study highlights an example of a nuanced mechanism of action for a disease-associated mutation, characterized by a change in protein-protein interaction specificity that alters enzyme activation.

Folding and Binding Kinetics of the Tandem of SH2 Domains from SHP2

The SH2 domains of SHP2 play a crucial role in determining the function of the SHP2 protein. While the folding and binding properties of the isolated NSH2 and CSH2 domains have been extensively studied, there is limited information about the tandem SH2 domains. This study aims to elucidate the folding and binding kinetics of the NSH2-CSH2 tandem domains of SHP2 through rapid kinetic experiments, complementing existing data on the isolated domains. The results indicate that while the domains generally fold and unfold independently, acidic pH conditions induce complex scenarios involving the formation of a misfolded intermediate. Furthermore, a comparison of the binding kinetics of isolated NSH2 and CSH2 domains with the NSH2-CSH2 tandem domains, using peptides that mimic specific portions of Gab2, suggests a dynamic interplay between NSH2 and CSH2 in binding Gab2 that modulate the microscopic association rate constant of the binding reaction. These findings, discussed in the context of previous research on the NSH2 and CSH2 domains, enhance our understanding of the function of the SH2 domain tandem of SHP2.

Learning peptide properties with positive examples only

Deep learning can create accurate predictive models by exploiting existing large-scale experimental data, and guide the design of molecules. However, a major barrier is the requirement of both positive and negative examples in the classical supervised learning frameworks. Notably, most peptide databases come with missing information and low number of observations on negative examples, as such sequences are hard to obtain using high-throughput screening methods. To address this challenge, we solely exploit the limited known positive examples in a semi-supervised setting, and discover peptide sequences that are likely to map to certain antimicrobial properties via positive-unlabeled learning (PU). In particular, we use the two learning strategies of adapting base classifier and reliable negative identification to build deep learning models for inferring solubility, hemolysis, binding against SHP-2, and non-fouling activity of peptides, given their sequence. We evaluate the predictive performance of our PU learning method and show that by only using the positive data, it can achieve competitive performance when compared with the classical positive-negative (PN) classification approach, where there is access to both positive and negative examples.

The pathogenic T42A mutation in SHP2 rewires the interaction specificity of its N-terminal regulatory domain

Mutations in the tyrosine phosphatase SHP2 are associated with a variety of human diseases. Most mutations in SHP2 increase its basal catalytic activity by disrupting auto-inhibitory interactions between its phosphatase domain and N-terminal SH2 (phosphotyrosine recognition) domain. By contrast, some disease-associated mutations located in the ligand-binding pockets of the N- or C-terminal SH2 domains do not increase basal activity and likely exert their pathogenicity through alternative mechanisms. We lack a molecular understanding of how these SH2 mutations impact SHP2 structure, activity, and signaling. Here, we characterize five SHP2 SH2 domain ligand-binding pocket mutants through a combination of high-throughput biochemical screens, biophysical and biochemical measurements, and molecular dynamics simulations. We show that, while some of these mutations alter binding affinity to phosphorylation sites, the T42A mutation in the N-SH2 domain is unique in that it also substantially alters ligand-binding specificity, despite being 8-10 Å from the specificity-determining region of the SH2 domain. This mutation exerts its effect on sequence specificity by remodeling the phosphotyrosine binding pocket, altering the mode of engagement of both the phosphotyrosine and surrounding residues on the ligand. The functional consequence of this altered specificity is that the T42A mutant has biased sensitivity toward a subset of activating ligands and enhances downstream signaling. Our study highlights an example of a nuanced mechanism of action for a disease-associated mutation, characterized by a change in protein-protein interaction specificity that alters enzyme activation.

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.

Anti-Tumor Potential of Post-Translational Modifications of PD-1

Programmed cell death protein-1 (PD-1) is a vital immune checkpoint molecule. The location, stability, and protein-protein interaction of PD-1 are significantly influenced by post-translational modification (PTM) of proteins. The biological information of PD-1, including its gene and protein structures and the PD-1/PD-L1 signaling pathway, was briefly reviewed in this review. Additionally, recent research on PD-1 post-translational modification, including the study of ubiquitination, glycosylation, phosphorylation, and palmitoylation, was summarized, and research strategies for PD-1 PTM drugs were concluded. At present, only a part of PD-1/PD-L1 treated patients (35-45%) are benefited from immunotherapies, and novel strategies targeting PTM of PD-1/PD-L1 may be important for anti-PD-1/PD-L1 non-responders (poor responders).

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.

6-Mercaptopurine potently inhibits recruitment of SHP2 by phosphorylated PD-1 to inhibit PD-1 signalling and enhance T cell function

Antibody inhibitors that block PD-1/PD-L1 interaction have been approved for oncological clinics, yielding impressive treatment effects. Small molecules inhibiting PD-1 signalling are at various stages of development, given that small molecular drugs are expected to outperform protein drugs in several ways. Currently, a significant portion of these small molecular inhibitors achieve this purpose by binding to a limited region of the PD-L1 protein, thereby limiting the choice of chemical structures. Alternative strategies for developing small-molecular PD-1 inhibitors are urgently needed to broaden the choice of chemical structures. Here, we report that 6-mercaptopurine (6-MP) inhibits PD-1 signalling, activates T cell function in vitro and in vivo and shrinks tumours by activating cytotoxic T cells. Mechanistically, 6-MP potently inhibited PD-1 signalling by blocking the recruitment of SHP2 by PD-1. Considering that 6-MP is a chemotherapeutic agent already approved by the FDA for childhood leukaemia, our work revealed a novel anti-tumour mechanism for this drug and suggests that 6-MP warrants further clinical evaluation for other tumour types.

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.

Mechanistic Elucidation of Activation/Deactivation Signal Transduction within Neurotensin Receptor 1 Triggered by 'Driver Chemical Groups' of Modulators: A Comparative Molecular Dynamics Simulation

Small-molecule modulators of neurotensin receptor 1 (NTSR1), a class A G-protein-coupled receptor (GPCR), has emerged as promising therapeutic agent for psychiatric disorders and cancer. Interestingly, a chemical group substitution in NTSR1 modulators can launch different types of downstream regulation, highlighting the significance of deciphering the internal fine-tuning mechanism. Here, we conducted a synergistic application of a Gaussian accelerated molecular dynamics simulation, a conventional molecular dynamics simulation, and Markov state models (MSM) to investigate the underlying mechanism of 'driver chemical groups' of modulators triggering inverse signaling. The results indicated that the flexibility of the leucine moiety in NTSR1 agonists contributes to the inward displacement of TM7 through a loosely coupled allosteric pathway, while the rigidity of the adamantane moiety in NTSR1 antagonists leads to unfavorable downward transduction of agonistic signaling. Furthermore, we found that R322^6.54, Y319^6.51, F353^7.42, R148^3.32, S356^7.45, and S357^7.46 may play a key role in inducing the activation of NTSR1. Together, our findings not only highlight the ingenious signal transduction within class A GPCRs but also lay a foundation for the development of targeted drugs harboring different regulatory functions of NTSR1.

Learning Peptide Properties with Positive Examples Only

Deep learning can create accurate predictive models by exploiting existing large-scale experimental data, and guide the design of molecules. However, a major barrier is the requirement of both positive and negative examples in the classical supervised learning frameworks. Notably, most peptide databases come with missing information and low number of observations on negative examples, as such sequences are hard to obtain using high-throughput screening methods. To address this challenge, we solely exploit the limited known positive examples in a semi-supervised setting, and discover peptide sequences that are likely to map to certain antimicrobial properties via positive-unlabeled learning (PU). In particular, we use the two learning strategies of adapting base classifier and reliable negative identification to build deep learning models for inferring solubility, hemolysis, binding against SHP-2, and non-fouling activity of peptides, given their sequence. We evaluate the predictive performance of our PU learning method and show that by only using the positive data, it can achieve competitive performance when compared with the classical positive-negative (PN) classification approach, where there is access to both positive and negative examples.

Experiment-guided molecular simulations define a heterogeneous structural ensemble for the PTPN11 tandem SH2 domains

SHP2 plays an important role in regulating cellular processes, and its pathogenic mutations cause developmental disorders and are linked to cancer. SHP2 is a multidomain protein, comprising two SH2 domains arranged in tandem, a catalytic PTP domain, and a disordered C-terminal tail. SHP2 is activated upon binding two linked phosphopeptides to its SH2 domains, and the peptide orientation and spacing between binding sites are critical for enzymatic activation. For decades, the tandem SH2 has been extensively studied to identify the relative orientation of the two SH2 domains that most effectively binds effectors. So far, neither crystallography nor experiments in solution have provided conclusive results. Using experiment-guided molecular simulations, we determine the heterogeneous structural ensemble of the tandem SH2 in solution in agreement with experimental data from small-angle X-ray scattering and NMR residual dipolar couplings. In the solution ensemble, N-SH2 adopts different orientations and positions relative to C-SH2. We suggest that the intrinsic structural plasticity of the tandem SH2 allows SHP2 to respond to external stimuli and is essential for its functional activity.

Setting sail: Maneuvering SHP2 activity and its effects in cancer

Since the discovery of tyrosine phosphorylation being a critical modulator of cancer signaling, proteins regulating phosphotyrosine levels in cells have fast become targets of therapeutic intervention. The nonreceptor protein tyrosine phosphatase (PTP) coded by the PTPN11 gene "SHP2" integrates phosphotyrosine signaling from growth factor receptors into the RAS/RAF/ERK pathway and is centrally positioned in processes regulating cell development and oncogenic transformation. Dysregulation of SHP2 expression or activity is linked to tumorigenesis and developmental defects. Even as a compelling anti-cancer target, SHP2 was considered "undruggable" for a long time owing to its conserved catalytic PTP domain that evaded drug development. Recently, SHP2 has risen from the "undruggable curse" with the discovery of small molecules that manipulate its intrinsic allostery for effective inhibition. SHP2's unique domain arrangement and conformation(s) allow for a truly novel paradigm of inhibitor development relying on skillful targeting of noncatalytic sites on proteins. In this review we summarize the biological functions, signaling properties, structural attributes, allostery and inhibitors of SHP2.

MUC1-C is necessary for SHP2 activation and BRAF inhibitor resistance in BRAF(V600E) mutant colorectal cancer

Colorectal cancers (CRCs) harboring the BRAF(V600E) mutation are associated with aggressive disease and resistance to BRAF inhibitors by feedback activation of the receptor tyrosine kinase (RTK)→RAS→MAPK pathway. The oncogenic MUC1-C protein promotes progression of colitis to CRC; whereas there is no known involvement of MUC1-C in BRAF(V600E) CRCs. The present work demonstrates that MUC1 expression is significantly upregulated in BRAF(V600E) vs wild-type CRCs. We show that BRAF(V600E) CRC cells are dependent on MUC1-C for proliferation and BRAF inhibitor (BRAFi) resistance. Mechanistically, MUC1-C integrates induction of MYC in driving cell cycle progression with activation of the SHP2 phosphotyrosine phosphatase, which enhances RTK-mediated RAS→ERK signaling. We demonstrate that targeting MUC1-C genetically and pharmacologically suppresses (i) activation of MYC, (ii) induction of the NOTCH1 stemness factor, and (iii) the capacity for self-renewal. We also show that MUC1-C associates with SHP2 and is required for SHP2 activation in driving BRAFi-induced feedback of ERK signaling. In this way, targeting MUC1-C in BRAFi-resistant BRAF(V600E) CRC tumors inhibits growth and sensitizes to BRAF inhibition. These findings demonstrate that MUC1-C is a target for the treatment of BRAF(V600E) CRCs and for reversing their resistance to BRAF inhibitors by suppressing the feedback MAPK pathway.

NMR Methods to Study the Dynamics of SH2 Domain-Phosphopeptide Complexes

Nuclear magnetic resonance (NMR) spectroscopy is the method of choice for studying the dynamics of biological macromolecules in solution. By exploiting the intricate interplay between the effects of protein motion (both overall rotational diffusion and internal mobility) and nuclear spin relaxation, NMR allows molecular motion to be probed at atomic resolution over a wide range of timescales, including picosecond (bond vibrations and methyl-group rotations), nanosecond (loop motions and rotational diffusion), and microsecond-millisecond (ligand binding, allostery). In this chapter, we describe different NMR pulse schemes (R1, R1ρ, heteronuclear NOE, and CPMG relaxation dispersion) to characterize the dynamics of SH2 domains. As an example, we use the N-SH2 domain of protein tyrosine phosphatase SHP2 in complex with two phosphopeptides derived from immune checkpoint receptor PD-1 (ITIM and ITSM).

Modeling (not so) rare developmental disorders associated with mutations in the protein-tyrosine phosphatase SHP2

Src homology region 2 (SH2)-containing protein tyrosine phosphatase 2 (SHP2) is a highly conserved protein tyrosine phosphatase (PTP), which is encoded by PTPN11 and is indispensable during embryonic development. Mutations in PTPN11 in human patients cause aberrant signaling of SHP2, resulting in multiple rare hereditary diseases, including Noonan Syndrome (NS), Noonan Syndrome with Multiple Lentigines (NSML), Juvenile Myelomonocytic Leukemia (JMML) and Metachondromatosis (MC). Somatic mutations in PTPN11 have been found to cause cancer. Here, we focus on the role of SHP2 variants in rare diseases and advances in the understanding of its pathogenesis using model systems.

A comprehensive review of SHP2 and its role in cancer

Src homology 2-containing protein tyrosine phosphatase 2 (SHP2) is a non-receptor protein tyrosine phosphatase ubiquitously expressed mainly in the cytoplasm of several tissues. SHP2 modulates diverse cell signaling events that control metabolism, cell growth, differentiation, cell migration, transcription and oncogenic transformation. It interacts with diverse molecules in the cell, and regulates key signaling events including RAS/ERK, PI3K/AKT, JAK/STAT and PD-1 pathways downstream of several receptor tyrosine kinases (RTKs) upon stimulation by growth factors and cytokines. SHP2 acts as both a phosphatase and a scaffold, and plays prominently oncogenic functions but can be tumor suppressor in a context-dependent manner. It typically acts as a positive regulator of RTKs signaling with some inhibitory functions reported as well. SHP2 expression and activity is regulated by such factors as allosteric autoinhibition, microRNAs, ubiquitination and SUMOylation. Dysregulation of SHP2 expression or activity causes many developmental diseases, and hematological and solid tumors. Moreover, upregulated SHP2 expression or activity also decreases sensitivity of cancer cells to anticancer drugs. SHP2 is now considered as a compelling anticancer drug target and several classes of SHP2 inhibitors with different mode of action are developed with some already in clinical trial phases. Moreover, novel SHP2 substrates and functions are rapidly growing both in cell and cancer. In view of this, we comprehensively and thoroughly reviewed literatures about SHP2 regulatory mechanisms, substrates and binding partners, biological functions, roles in human cancers, and different classes of small molecule inhibitors target this oncoprotein in cancer.

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.

Mechanistic Insights into the Long-range Allosteric Regulation of KRAS Via Neurofibromatosis Type 1 (NF1) Scaffold Upon SPRED1 Loading

Allosteric regulation is the most direct and efficient way of regulating protein function, wherein proteins transmit the perturbations at one site to another distinct functional site. Deciphering the mechanism of allosteric regulation is of vital importance for the comprehension of both physiological and pathological events in vivo as well as the rational allosteric drug design. However, it remains challenging to elucidate dominant allosteric signal transduction pathways, especially for large and multi-component protein machineries where long-range allosteric regulation exits. One of the quintessential examples having long-range allosteric regulation is the ternary complex, SPRED1-RAS-neurofibromin type 1 (NF1, a RAS GTPase-activating protein), in which SPRED1 facilitates RAS-GTP hydrolysis by interacting with NF1 at a distal, allosteric site from the RAS binding site. To address the underlying mechanism, we performed extensive Gaussian accelerated molecular dynamics simulations and Markov state model analysis of KRAS-NF1 complex in the presence and absence of SPRED1. Our findings suggested that SPRED1 loading allosterically enhanced KRAS-NF1 binding, but hindered conformational transformation of the NF1 catalytic center for RAS hydrolysis. Moreover, we unveiled the possible allosteric pathways upon SPRED1 binding through difference contact network analysis. This study not only provided an in-depth mechanistic insight into the allosteric regulation of KRAS by SPRED1, but also shed light on the investigation of long-range allosteric regulation among complex macromolecular systems.

Understanding the P-Loop Conformation in the Determination of Inhibitor Selectivity Toward the Hepatocellular Carcinoma-Associated Dark Kinase STK17B

As a member of the death-associated protein kinase family of serine/threonine kinases, the STK17B has been associated with diverse diseases such as hepatocellular carcinoma. However, the conformational dynamics of the phosphate-binding loop (P-loop) in the determination of inhibitor selectivity profile to the STK17B are less understood. Here, a multi-microsecond length molecular dynamics (MD) simulation of STK17B in the three different states (ligand-free, ADP-bound, and ligand-bound states) was carried out to uncover the conformational plasticity of the P-loop. Together with the analyses of principal component analysis, cross-correlation and generalized correlation motions, secondary structural analysis, and community network analysis, the conformational dynamics of the P-loop in the different states were revealed, in which the P-loop flipped into the ADP-binding site upon the inhibitor binding and interacted with the inhibitor and the C-lobe, strengthened the communication between the N- and C-lobes. These resulting interactions contributed to inhibitor selectivity profile to the STK17B. Our results may advance our understanding of kinase inhibitor selectivity and offer possible implications for the design of highly selective inhibitors for other protein kinases.

Autopromotion of K-Ras4B Feedback Activation Through an SOS-Mediated Long-Range Allosteric Effect

The Ras-specific guanine nucleotide exchange factors Son of Sevenless (SOS) regulates Ras activation by converting inactive GDP-bound to active GTP-bound states. The catalytic activity of Ras is further allosterically regulated by GTP-Ras bound to a distal site through a positive feedback loop. To address the mechanism underlying the long-range allosteric activation of the catalytic K-Ras4B by an additional allosteric GTP-Ras through SOS, we employed molecular dynamics simulation of the K-Ras4BG13D•SOScat complex with and without an allosteric GTP-bound K-Ras4BG13D. We found that the binding of an allosteric GTP-K-Ras4BG13D enhanced the affinity between the catalytic K-Ras4BG13D and SOScat, forming a more stable conformational state. The peeling away of the switch I from the nucleotide binding site facilitated the dissociation of GDP, thereby contributing to the increased nucleotide exchange rate. The community networks further showed stronger edge connection upon allosteric GTP-K-Ras4BG13D binding, which represented an increased interaction between catalytic K-Ras4BG13D and SOScat. Moreover, GTP-K-Ras4BG13D binding transmitted allosteric signaling pathways though the Cdc25 domain of SOS that enhanced the allosteric regulatory from the K-Ras4BG13D allosteric site to the catalytic site. This study may provide an in-depth mechanism for abnormal activation and allosteric regulation of K-Ras4BG13D.

Delineating the activation mechanism and conformational landscape of a class B G protein-coupled receptor glucagon receptor

Class B G protein-coupled receptors (GPCRs) are important targets in the treatment of metabolic syndrome and diabetes. Although multiple structures of class B GPCRs-G protein complexes have been elucidated, the detailed activation mechanism of the receptors remains unclear. Here, we combine Gaussian accelerated molecular dynamics simulations and Markov state models (MSM) to investigate the activation mechanism of a canonical class B GPCR, human glucagon receptor-GCGR, including the negative allosteric modulator-bound inactive state, the agonist glucagon-bound active state, and both glucagon- and Gs-bound fully active state. The free-energy landscapes of GCGR show the conformational ensemble consisting of three activation-associated states: inactive, active, and fully active. The structural analysis indicates the high dynamics of GCGR upon glucagon binding with both active and inactive conformations in the ensemble. Significantly, the H8 and TM6 exhibits distinct features from the inactive to the active states. The additional simulations demonstrate the role of H8 in the recruitment of Gs. Gs binding presents a crucial function of stabilizing the glucagon binding site and MSM highlights the absolute requirement of Gs to help the GCGR reach the fully active state. Together, our results reveal the detailed activation mechanism of GCGR from the view of conformational dynamics.

Atomic-scale insights into allosteric inhibition and evolutional rescue mechanism of Streptococcus thermophilus Cas9 by the anti-CRISPR protein AcrIIA6

CRISPR-Cas systems are prokaryotic adaptive immunity against invading phages and plasmids. Phages have evolved diverse protein inhibitors of CRISPR-Cas systems, called anti-CRISPR (Acr) proteins, to neutralize this CRISPR machinery. In response, bacteria have co-evolved Cas variants to escape phage's anti-CRISPR strategies, called anti-anti-CRISPR systems. Here we explore the anti-CRISPR allosteric inhibition and anti-anti-CRISPR rescue mechanisms between Streptococcus thermophilus Cas9 (St1Cas9) and the anti-CRISPR protein AcrIIA6 at the atomic level, by generating mutants of key residues in St1Cas9. Extensive unbiased molecular dynamics simulations show that the functional motions of St1Cas9 in the presence of AcrIIA6 differ substantially from those of St1Cas9 alone. AcrIIA6 binding triggers a shift of St1Cas9 conformational ensemble towards a less catalytically competent state; this state significantly compromises protospacer adjacent motif (PAM) recognition and nuclease activity by altering interdependently conformational dynamics and allosteric signals among nuclease domains, PAM-interacting (PI) regions, and AcrIIA6 binding motifs. Via in vitro DNA cleavage assays, we further elucidate the rescue mechanism of efficiently escaping AcrIIA6 inhibition harboring St1Cas9 triple mutations (G993K/K1008M/K1010E) in the PI domain and identify the evolutionary landscape of such mutational escape within species. Our results provide mechanistic insights into Acr proteins as natural brakes for the CRISPR-Cas systems and a promising potential for the design of allosteric Acr peptidomimetics.

Targeting Oncogenic Src Homology 2 Domain-Containing Phosphatase 2 (SHP2) by Inhibiting Its Protein-Protein Interactions

We developed a new class of inhibitors of protein-protein interactions of the SHP2 phosphatase, which is pivotal in cell signaling and represents a central target in the therapy of cancer and rare diseases. Currently available SHP2 inhibitors target the catalytic site or an allosteric pocket but lack specificity or are ineffective for disease-associated SHP2 mutants. Considering that pathogenic lesions cause signaling hyperactivation due to increased levels of SHP2 association with cognate proteins, we developed peptide-based molecules with nanomolar affinity for the N-terminal Src homology domain of SHP2, good selectivity, stability to degradation, and an affinity for pathogenic variants of SHP2 that is 2-20 times higher than for the wild-type protein. The best peptide reverted the effects of a pathogenic variant (D61G) in zebrafish embryos. Our results provide a novel route for SHP2-targeted therapies and a tool for investigating the role of protein-protein interactions in the function of SHP2.

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

The Src-homology 2 domain containing phosphatase 2 (SHP2) plays a critical role in crucial signaling pathways and is involved in oncogenesis and in developmental disorders. Its structure includes two SH2 domains (N-SH2 and C-SH2), and a protein tyrosine phosphatase (PTP) domain. Under basal conditions, SHP2 is auto-inhibited, with the N-SH2 domain blocking the PTP active site. Activation involves a rearrangement of the domains that makes the catalytic site accessible, coupled to the association between the SH2 domains and cognate proteins containing phosphotyrosines. Several aspects of this transition are debated and competing mechanistic models have been proposed. A crystallographic structure of SHP2 in an active state has been reported (PDB code 6crf), but several lines of evidence suggests that it is not fully representative of the conformations populated in solution. To clarify the structural rearrangements involved in SHP2 activation, enhanced sampling simulations of the autoinhibited and active states have been performed, for wild type SHP2 and its pathogenic E76K variant. Our results demonstrate that the crystallographic conformation of the active state is unstable in solution, and multiple interdomain arrangements are populated, thus allowing association to bisphosphorylated sequences. Contrary to a recent proposal, activation is coupled to the conformational changes of the N-SH2 binding site, which is significantly more accessible in the active sate, rather than to the structure of the central β-sheet of the domain. In this coupling, a previously undescribed role for the N-SH2 BG loop emerged.