Study on SHP2 Conformational Transition and Structural Characterization of Its High-Potency Allosteric Inhibitors by Molecular Dynamics Simulations Combined with Machine Learning

Abstract

The src-homology 2 domain-containing phosphatase 2 (SHP2) is a human cytoplasmic protein tyrosine phosphatase that plays a crucial role in cellular signal transduction. Aberrant activation and mutations of SHP2 are associated with tumor growth and immune suppression, thus making it a potential target for cancer therapy. Initially, researchers sought to develop inhibitors targeting SHP2’s catalytic site (protein tyrosine phosphatase domain, PTP). Due to limitations such as conservativeness and poor membrane permeability, SHP2 was once considered a challenging drug target. Nevertheless, with the in-depth investigations into the conformational switch mechanism from SHP2’s inactive to active state and the emergence of various SHP2 allosteric inhibitors, new hope has been brought to this target. In this study, we investigated the interaction models of various allosteric inhibitors with SHP2 using molecular dynamics simulations. Meanwhile, we explored the free energy landscape of SHP2 activation using enhanced sampling technique (meta-dynamics simulations), which provides insights into its conformational changes and activation mechanism. Furthermore, to biophysically interpret high-dimensional simulation trajectories, we employed interpretable machine learning methods, specifically extreme gradient boosting (XGBoost) with Shapley additive explanations (SHAP), to comprehensively analyze the simulation data. This approach allowed us to identify and highlight key structural features driving SHP2 conformational dynamics and regulating the activity of the allosteric inhibitor. These studies not only enhance our understanding of SHP2’s conformational switch mechanism but also offer crucial insights for designing potent allosteric SHP2 inhibitors and addressing drug resistance issues.