Sourcing the affinity of flavonoids for the glycogen phosphorylase inhibitor site via crystallography, kinetics and QM/MM-PBSA binding studies: Comparison of chrysin and flavopiridol

Flavonoids have been discovered as novel inhibitors of glycogen phosphorylase (GP), a target to control hyperglycemia in type 2 diabetes. To elucidate the mechanism of inhibition, we have determined the crystal structure of the GPb-chrysin complex at 1.9 angstrom resolution. Chrysin is accommodated at the inhibitor site intercalating between the aromatic side chains of Phe285 and Tyr613 through pi-stacking interactions. Chrysin binds to GPb similar to 15 times weaker (K-i = 19.01 mu M) than flavopiridol (K-i = 1.24 mu M), exclusively at the inhibitor site, and both inhibitors display similar behavior with respect to AMP. To identify the source of flavopiridols' stronger affinity, molecular docking with Glide and postdocking binding free energy calculations using QM/MM-PBSA have been performed and compared. Whereas docking failed to correctly rank inhibitor binding conformations, the QM/MM-PBSA method employing M06-2X/6-31+G* to model the pi-stacking interactions correctly reproduced the experimental results. Flavopiridols' greater binding affinity is sourced to favorable interactions of the cationic 4-hydroxypiperidin-1-yl substituent with GPb, with desolvation effects limited by the substituent conformation adopted in the crystallographic complex. Further successful predictions using QM/MM-PBSA for the flavonoid quercetagetin (which binds at the allosteric site) leads us to propose the methodology as a useful and inexpensive tool to predict flavonoid binding. (C) 2012 Elsevier Ltd. All rights reserved.