Prediction of SHP2-E76K binding sites based on molecular dynamics simulation and Markov algorithm

J Biomol Struct Dyn. 2024 Nov 19:1-12. doi: 10.1080/07391102.2024.2431193. Online ahead of print.

Abstract

SHP2-E76K, a mutant encoded by the PTPN11 gene, was associated with various solid tumors, such as lung cancer, glioblastoma, and intellectual disability. SHP2-E76K has become potential drug targets, while there was no effective inhibitor against the mutant currently. At present, the crystal complex structure of SHP099 with SHP2-E76K has been reported in the RCSB PDB protein data bank, however, the dynamic structure of SHP099 binding to the active center of SHP2-E76K protein was still lacking. Therefore, this study used molecular dynamics simulation and Markov model to characterize the kinetics of the inhibitor SHP099 with SHP2-E76K enzyme and to determine the active binding site, which would give a hint of a vital enzyme-substrate interaction in atomistic detail that proposed the potential to be applied for the discovery of more effective SHP2-E76K inhibitors and, in broader terms, dynamic protein-drug interactions.

Keywords: MM/PBSA; Markov model; SHP2-E76K; transition path theory; umbrella sampling.