Lead optimization of Allium sativum L. compounds for PTP1B inhibition in diabetes treatment: in silico molecular docking and dynamics simulation

Abstract

Protein tyrosine phosphatase 1B (PTP1B) has been identified as a promising drug target for the development of diabetes medications via an inhibition mechanism. Using a computational approach, this study investigates the binding mechanism of lead optimized natural compounds from Allium sativum against the human PTP1B. The molecular docking, induced-fit docking, and binding free energy calculations were analyzed using Schrödinger Suite 2021-2. MD simulation, and gene enrichment analysis was achieved via the Desmond module of Schrödinger to identify best compounds as inhibitors against PTP1B in diabetes management. The docking scores of the lead optimized compounds were good; 5280443_121 from apigenin had the best binding score of −9.345 kcal/mol, followed by 5280443_129 with a binding score of −9.200 kcal/mol, and 5280863_177 from kaempferol had a binding score of −8.528 kcal/mol, followed by 5280863_462 with a binding score of −8.338 kcal/mol. The top two lead optimized compounds, docked better than the standard PTP1B inhibitor (−7.155 kcal/mol), suggesting them as potent inhibitors than the standard PTP1B inhibitor. The outcomes of the induced-fit docking were consistent with the increased binding affinity used in the Glide computation of the five conformed poses between the derivatives (5280443_121, 5280443_129, 5280863_177, and 5280863_462) and the protein (PTP1B). Based on the binding fee energies (MM-GBSA), the lead optimized compounds from kaempferol exhibited more stability than those from apigenin. In the pharmacophore development, all the models exhibit good results across the different metrics. The best performing model with five of five matches on a 1.34 and 1.33 phase score was DDRRR_1, DDRRR_2, and DDDRR_1. The average BEDROC value (= 160.9) was 1, while the average EF 1% value across all models was 101. There were no substantial conformational modifications during the MD simulation process, indicating that the apigenin derivatives (5280443_121) was stable in the protein’s active site in 100 ns. IGF1R, EGFR, INSR, PTPN1, SRC, JAK2, GRB2, BCAR1, and IRS1 are among the 11 potential targets found in the protein–protein interaction (PPI) of A. sativum against PTP1B that may be important in A. sativum’s defense against PTP1B. Sixty-four (64) pathways were found by KEGG pathway enrichment analysis to be potentially involved in the anti-PTP1B of A. sativum. Consequently, data obtained indicates the effectiveness of the in silico studies in identifying potential lead compounds in A. sativum against PTP1B target. Communicated by Ramaswamy H. Sarma