Deciphering SHP-2 conformational switching via elastic network models and transfer entropy analysis.

contextSrc homology 2 domain-containing protein tyrosine phosphatase-2 (SHP-2), a non-receptor protein tyrosine phosphatase encoded by PTPN11, functions as a central regulator of the Ras-MAPK, JAK-STAT, and PI3K-AKT signaling pathways. Its gain-of-function mutations, such as E76K, drive aberrant activation in multiple cancers, making it an attractive therapeutic target. While molecular dynamics simulations have provided atomic-level insights into SHP-2 conformational dynamics, they are computationally prohibitive for systematically characterizing large-scale domain rearrangements and allosteric communication networks. Here, we employed an equal-weighted multi-parameter Gaussian network model (ewmGNM) to investigate the two lowest-frequency modes governing the functional dynamics of SHP-2 in three key systems: the wild-type closed (autoinhibited) conformation, the oncogenic E76K mutant in its open (active) conformation, and its complex with the allosteric inhibitor SHP099. By focusing on low-frequency modes, we computed residue fluctuations, dynamic cross-correlations, and transfer entropy to map directional information flow and allosteric communication landscapes. Our results demonstrate that the E76K mutation globally increases protein flexibility and disrupts inter-domain coupling, particularly at the N-SH2/PTP interface, whereas SHP099 binding largely restores the closed state dynamics but retains mutation-induced residual flexibility. Transfer entropy analysis reveals a conformation-dependent reorganization of information flow, with the open state exhibiting that the N-SH2 domain possibly acts as a new entropy source, tentatively suggesting its potential as a novel allosteric site. These findings provide a mechanistic understanding of SHP-2 conformational regulation and offer a theoretical basis for the design of allosteric inhibitors targeting this open state site.

methodsProtein crystal structures were downloaded from the Protein Data Bank. The ewmGNM and the conventional GNM were implemented using custom scripts written in MATLAB R2024b. Protein dynamics data were mapped using custom Python scripts in PyMOL.