A procedure combining molecular docking and semiempirical method PM7 for identification of selective Shp2 inhibitors

Shp2 and Shp1 make up a small family of protein tyrosine phosphatases. Finding selective inhibitors for Shp2 is useful because although its inhibition is advantageous for the treatment of some types of cancer, inhibition of Shp1 may have the opposite effect, since it acts as a suppressor of tumors. We combined molecular docking and semiempirical molecular orbital-based calculations to produce data that were effective for the identification of selective inhibitors for Shp2. After definition of the interaction modes of the inhibitors with Shp2 and Shp1 by molecular docking, the resulting interaction enthalpy values were calculated following refinement of the enzyme/inhibitor complexes' geometries with the PM7 semiempirical molecular orbital method. Despite the complexity of the thermodynamics involved in the enzyme/inhibitor interaction, the selectivity for Shp2 of a series of 76 inhibitors, divided in two groups, could be effectively correlated with the difference in their interaction enthalpy values with both enzymes. For the first group, composed by 52 Shp2-selective indoline inhibitors for which only Shp2 inhibition activity data are available, we demonstrated that the interaction enthalpy can be used as a criterion for identification of selective Shp2 inhibitors, as it is significantly more favorable for Shp2 than Shp1 at a 99% confidence level. For the second group, composed by 24 oxindole derivatives with available Shp2 and Shp1 inhibition activity data, a satisfactory correlation (R = 0.70) could be obtained between the selectivity, based on the IC50 data, and the relative percentage difference of the calculated interaction enthalpies with the two enzymes.