3 '-Axial CH2OH Substitution on Glucopyranose does not Increase Glycogen Phosphorylase Inhibitory Potency. QM/MM-PBSA Calculations Suggest Why
AuthorManta, S.; Xipnitou, A.; Kiritsis, C.; Kantsadi, A. L.; Hayes, J. M.; Skamnaki, V. T.; Lamprakis, C.; Kontou, M.; Zoumpoulakis, P.; Zographos, S. E.; Leonidas, D. D.; Komiotis, D.
Glycogen phosphorylase is a molecular target for the design of potential hypoglycemic agents. Structure-based design pinpointed that the 3'-position of glucopyranose equipped with a suitable group has the potential to form interactions with enzymes cofactor, pyridoxal 5'-phosphate (PLP), thus enhancing the inhibitory potency. Hence, we have investigated the binding of two ligands, 1-(beta-d-glucopyranosyl)5-fluorouracil (GlcFU) and its 3'-CH2OH glucopyranose derivative. Both ligands were found to be low micromolar inhibitors with K-i values of 7.9 and 27.1 mu m, respectively. X-ray crystallography revealed that the 3'-CH2OH glucopyranose substituent is indeed involved in additional molecular interactions with the PLP gamma-phosphate compared with GlcFU. However, it is 3.4 times less potent. To elucidate this discovery, docking followed by postdocking Quantum Mechanics/Molecular Mechanics PoissonBoltzmann Surface Area (QM/MM-PBSA) binding affinity calculations were performed. While the docking predictions failed to reflect the kinetic results, the QM/MM-PBSA revealed that the desolvation energy cost for binding of the 3'-CH2OH-substituted glucopyranose derivative out-weigh the enthalpy gains from the extra contacts formed. The benefits of performing postdocking calculations employing a more accurate solvation model and the QM/MM-PBSA methodology in lead optimization are therefore highlighted, specifically when the role of a highly polar/charged binding interface is significant.