Search

Glycogen phosphorylase (GP), a validated target for the development of anti-hyperglycaemic agents, has been targeted for the design of novel glycopyranosylamine inhibitors. Exploiting the two most potent inhibitors from our previous study of N-acyl-β-D-glucopyranosylamines (Parmenopoulou et al., Bioorg. Med. Chem. 2014, 22, 4810), we have extended the linking group to -NHCONHCO- between the glucose moiety and the aliphatic/aromatic substituent in the GP catalytic site β-cavity. The N-acyl-N´-(β-D-glucopyranosyl) urea inhibitors were synthesized and their efficiency assessed by biochemical methods, revealing inhibition constant values of 4.95 µM and 2.53 µM. Crystal structures of GP in complex with these inhibitors were determined and analyzed, providing data for further structure based design efforts. A novel Linear Response - Molecular Mechanics Coulomb Surface Area (LR-MM-CBSA) method has been developed which relates predicted and experimental binding free energies for a training set of N-acyl-N´-(β-D-glucopyranosyl) urea ligands with a correlation coefficient R(2) of 0.89 and leave-one-out cross-validation (LOO-cv) Q(2) statistic of 0.79. The method has significant applications to direct future lead optimization studies, where ligand entropy loss on binding is revealed as a key factor to be considered. ADMET property predictions revealed that apart from potential permeability issues, the synthesized N-acyl-N´-(β-D-glucopyranosyl) urea inhibitors have drug-like potential without any toxicity warnings.

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 enzyme's cofactor, pyridoxal 5'-phosphate (PLP), thus enhancing the inhibitory potency. Hence, we have investigated the binding of two ligands, 1-(β-d-glucopyranosyl)5-fluorouracil (GlcFU) and its 3'-CH(2) OH glucopyranose derivative. Both ligands were found to be low micromolar inhibitors with K(i) values of 7.9 and 27.1 μm, respectively. X-ray crystallography revealed that the 3'-CH(2) OH glucopyranose substituent is indeed involved in additional molecular interactions with the PLP γ-phosphate compared with GlcFU. However, it is 3.4 times less potent. To elucidate this discovery, docking followed by postdocking Quantum Mechanics/Molecular Mechanics - Poisson-Boltzmann 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'-CH(2) OH-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., (© 2012 John Wiley & Sons A/S.)

A series of glucopyranosylidene-spiro-isoxazolines was prepared through regio- and stereoselective [3+2]-cycloaddition between the methylene acetylated exo-glucal and aromatic nitrile oxides. The deprotected cycloadducts were evaluated as inhibitors of muscle glycogen phosphorylase b. The carbohydrate-based family of five inhibitors displays K(i) values ranging from 0.63 to 92.5 microM. The X-ray structures of the enzyme-ligand complexes show that the inhibitors bind preferentially at the catalytic site of the enzyme retaining the less active T-state conformation. Docking calculations with GLIDE in extra-precision (XP) mode yielded excellent agreement with experiment, as judged by comparison of the predicted binding modes of the five ligands with the crystallographic conformations and the good correlation between the docking scores and the experimental free binding energies. Use of docking constraints on the well-defined positions of the glucopyranose moiety in the catalytic site and redocking of GLIDE-XP poses using electrostatic potential fit-determined ligand partial charges in quantum polarized ligand docking (QPLD) produced the best results in this regard.

4-Phenyl-N-(beta-D-glucopyranosyl)-1H-1,2,3-triazole-1-acetamide (glucosyltriazolylacetamide) has been studied in kinetic and crystallographic experiments with glycogen phosphorylase b (GPb), in an effort to utilize its potential as a lead for the design of potent antihyperglycaemic agents. Docking and molecular dynamics (MD) calculations have been used to monitor more closely the binding modes in operation and compare the results with experiment. Kinetic experiments in the direction of glycogen synthesis showed that glucosyltriazolylacetamide is a better inhibitor (K(i) = 0.18 mM) than the parent compound alpha-D-glucose (K(i) = 1.7 mM) or beta-D-glucose (K(i) = 7.4 mM) but less potent inhibitor than the lead compound N-acetyl-beta-D-glucopyranosylamine (K(i) = 32 microM). To elucidate the molecular basis underlying the inhibition of the newly identified compound, we determined the structure of GPb in complex with glucosyltriazolylacetamide at 100 K to 1.88 A resolution, and the structure of the compound in the free form. Glucosyltriazolylacetamide is accommodated in the catalytic site of the enzyme and the glucopyranose interacts in a manner similar to that observed in the GPb-alpha-D-glucose complex, while the substituent group in the beta-position of the C1 atom makes additional hydrogen bonding and van der Waals interactions to the protein. A bifurcated donor type hydrogen bonding involving O3H, N3, and N4 is seen as an important structural motif strengthening the binding of glucosyltriazolylacetamide with GP which necessitated change in the torsion about C8-N2 bond by about 62 degrees going from its free to the complex form with GPb. On binding to GP, glucosyltriazolylacetamide induces significant conformational changes in the vicinity of this site. Specifically, the 280s loop (residues 282-288) shifts 0.7 to 3.1 A (CA atoms) to accommodate glucosyltriazolylacetamide. These conformational changes do not lead to increased contacts between the inhibitor and the protein that would improve ligand binding compared with the lead compound. In the molecular modeling calculations, the GOLD docking runs with and without the crystallographic ordered cavity waters using the GoldScore scoring function, and without cavity waters using the ChemScore scoring function successfully reproduced the crystallographic binding conformation. However, the GLIDE docking calculations both with (GLIDE XP) and without (GLIDE SP and XP) the cavity water molecules were, impressively, further able to accurately reproduce the finer details of the GPb-glucosyltriazolylacetamide complex structure. The importance of cavity waters in flexible receptor MD calculations compared to "rigid" (docking) is analyzed and highlighted, while in the MD itself very little conformational flexibility of the glucosyltriazolylacetamide ligand was observed over the time scale of the simulations., (2007 Wiley-Liss, Inc.)