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

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

An integrated drug repurposing strategy for the rapid identification of potential SARS-CoV-2 viral inhibitors

The Coronavirus disease 2019 (COVID-19) is an infectious disease caused by the severe acute respiratory syndrome-coronavirus 2 (SARS-CoV-2). The virus has rapidly spread in humans, causing the ongoing Coronavirus pandemic. Recent studies have shown that, similarly to SARS-CoV, SARS-CoV-2 utilises the Spike glycoprotein on the envelope to recognise and bind the human receptor ACE2. This event initiates the fusion of viral and host cell membranes and then the viral entry into the host cell. Despite several ongoing clinical studies, there are currently no approved vaccines or drugs that specifically target SARS-CoV-2. Until an effective vaccine is available, repurposing FDA approved drugs could significantly shorten the time and reduce the cost compared to de novo drug discovery. In this study we attempted to overcome the limitation of in silico virtual screening by applying a robust in silico drug repurposing strategy. We combined and integrated docking simulations, with molecular dynamics (MD), Supervised MD (SuMD) and Steered MD (SMD) simulations to identify a Spike protein – ACE2 interaction inhibitor. Our data showed that Simeprevir and Lumacaftor bind the receptor-binding domain of the Spike protein with high affinity and prevent ACE2 interaction

nMoldyn - Interfacing spectroscopic experiments, molecular dynamics simulations and models for time correlation functions

This article gives an introduction into the program nMoldyn, which has been originally conceived to support the interpretation of neutron scattering experiments on complex molecular systems by the calculation of appropriate time correlation functions from classical and quantum molecular dynamics simulations of corresponding model systems. Later the functionality has been extended to include more advanced time series analyses of molecular dynamics trajectories, in particular the calculation of memory functions, which play an essential role in the theory of time correlation functions. Here we present a synoptic view of the range of applications of the latest version of nMoldyn, which includes new modules for a simulation-based interpretation of data from nuclear magnetic resonance spectroscopy, far infrared spectroscopy and for protein secondary structure analysis

Structural determinants of phosphopeptide binding to the N-terminal Src homology 2 domain of the SHP2 phosphatase

SH2 domain-containing tyrosine phosphatase 2 (SHP2), encoded by PTPN11, plays a fundamental role in the modulation of several signaling pathways. Germline and somatic mutations in PTPN11 are associated with different rare diseases and hematologic malignancies, and recent studies have individuated SHP2 as a central node in oncogenesis and cancer drug resistance. The SHP2 structure includes two Src homology 2 domains (N-SH2 and C-SH2) followed by a catalytic protein tyrosine phosphatase (PTP) domain. Under basal conditions, the N-SH2 domain blocks the active site, inhibiting phosphatase activity. Association of the N-SH2 domain with binding partners containing short amino acid motifs comprising a phosphotyrosine residue (pY) leads to N-SH2/PTP dissociation and SHP2 activation. Considering the relevance of SHP2 in signaling and disease and the central role of the N-SH2 domain in its allosteric regulation mechanism, we performed microsecond-long molecular dynamics (MD) simulations of the N-SH2 domain complexed to 12 different peptides to define the structural and dynamical features determining the binding affinity and specificity of the domain. Phosphopeptide residues at position -2 to +5, with respect to pY, have significant interactions with the SH2 domain. In addition to the strong interaction of the pY residue with its conserved binding pocket, the complex is stabilized hydrophobically by insertion of residues +1, +3, and +5 in an apolar groove of the domain and interaction of residue -2 with both the pY and a protein surface residue. Additional interactions are provided by hydrogen bonds formed by the backbone of residues -1, +1, +2, and +4. Finally, negatively charged residues at positions +2 and +4 are involved in electrostatic interactions with two lysines (Lys89 and Lys91) specific for the SHP2 N-SH2 domain. Interestingly, the MD simulations illustrated a previously undescribed conformational flexibility of the domain, involving the core beta sheet and the loop that closes the pY binding pocket