The Molecular Modeling of Novel Inhibitors of Protein Tyrosine Phosphatase 1B Based on Catechol by MD and MM-GB (PB)/SA Calculations

Abstract

Binding modes of a series of catechol derivatives such as protein tyrosine phosphatase 1B (PTP1B) inhibitors were identified by molecular modeling techniques. Docking, molecular dynamics simulations and free energy calculations were employed to determine the modes of these new inhibitors. Binding free energies were calculated by involving different energy components using the Molecular Mechanics-Poisson-Boltzmann Surface Area and Generalized Born Surface Area methods. Relatively larger binding energies were obtained for the catechol derivatives compared to one of the PTP1B inhibitors already in use. The Molecular Mechanics/ Generalized Born Surface Area (MM/GBSA) free energy decomposition analysis indicated that the hydroxyl functional groups and biphenyl ring system had favorable interactions with Met258, Tyr46, Gln262 and Phe182 residues of PTP1B. The results of hydrogen bound analysis indicated that catechol derivatives, in addition to hydrogen bonding interactions, Val49, Ile219, Gln266, Asp181 and amino acid residues of PTP1B are responsible for governing the inhibitor potency of the compounds. The information generated from the present study should be useful for the design of more potent PTP1B inhibitors as anti-diabetic agents.