Binding of β-d-glucopyranosyl bismethoxyphosphoramidate to glycogen phosphorylase b: kinetic and crystallographic studies

Outliers in SAR and QSAR: 4. effects of allosteric protein-ligand interactions on the classical quantitative structure-activity relationships

Effects of allosteric interactions on the classical structure-activity relationship (SAR) and quantitative SAR (QSAR) have been investigated. Apprehending the outliers in SAR and QSAR studies can improve the quality, predictability, and use of QSAR in designing unknown compounds in drug discovery research. We explored allosteric protein-ligand interactions as a possible source of outliers in SAR/QSAR. We used glycogen phosphorylase as an example of a protein that has an allosteric site. Examination of the ligand-bound x-ray crystal structures of glycogen phosphorylase revealed that many inhibitors bound at more than one binding site. The results of QSAR analyses of the inhibitors included a QSAR that recognized an outlier bound at a distinctive allosteric binding site. The case provided an example of constructive use of QSAR identifying outliers with alternative binding modes. Other allosteric QSARs that captured our attention were the inverted parabola/bilinear QSARs. The x-ray crystal structures and the QSAR analyses indicated that the inverted parabola QSARs could be associated with the conformational changes in the allosteric interactions. Our results showed that the normal parabola, as well as the inverted parabola QSARs, can describe the allosteric interactions. Examination of the ligand-bound X-ray crystal structures of glycogen phosphorylase revealed that many inhibitors bound at more than one binding site. The results of QSAR analyses of the inhibitors included a QSAR that recognized an outlier bound at a distinctive allosteric binding site.

Iminosugars as Potential Inhibitors of Glycogenolysis: Structural Insights into the Molecular Basis of Glycogen Phosphorylase Inhibition

Iminosugars DAB (5), isofagomine (9), and several N-substituted derivatives have been identified as potent inhibitors of liver glycogen phosphorylase a (IC(50) = 0.4-1.2 microM) and of basal and glucagon-stimulated glycogenolysis (IC(50) = 1-3 microM). The X-ray structures of 5, 9, and its N-3-phenylpropyl analogue 8 in complex with rabbit muscle glycogen phosphorylase (GPb) shows that iminosugars bind tightly at the catalytic site in the presence of the substrate phosphate and induce conformational changes that characterize the R-state conformation of the enzyme. Charged nitrogen N1 is within hydrogen-bonding distance with the carbonyl oxygen of His377 (5) and in ionic contact with the substrate phosphate oxygen (8 and 9). Our findings suggest that the inhibitors function as oxocarbenium ion transition-state analogues. The conformational change to the R state provides an explanation for previous findings that 5, unlike inhibitors that favor the T state, promotes phosphorylation of GPb in hepatocytes with sequential inactivation of glycogen synthase.

Crystallographic studies on N-azidoacetyl-β-d-glucopyranosylamine, an inhibitor of glycogen phosphorylase: Comparison with N-acetyl-β-d-glucopyranosylamine

N-acetyl-beta-D-glucopyranosylamine (NAG) is a potent inhibitor (Ki=32 microM) of glycogen phosphorylase b (GPb), and has been employed as a lead compound for the structure-based design of new analogues, in an effort to utilize its potential as a hypoglycaemic agent. Replacement of the acetamido group by azidoacetamido group resulted in an inhibitor, N-azidoacetyl-beta-D-glucopyranosylamine (azido-NAG), with a Ki value of 48.7 microM, in the direction of glycogen synthesis. In order to elucidate the mechanism of inhibition, we determined the ligand structure in complex with GPb at 2.03 A resolution, and the structure of the fully acetylated derivative in the free form. The molecular packing of the latter is stabilized by a number of bifurcated hydrogen bonds of which the one involving a bifurcated C-H...N...H-C type hydrogen bonding is rather unique in organic azides. Azido-NAG can be accommodated in the catalytic site of T-state GPb at approximately the same position as that of NAG and stabilizes the T-state conformation of the 280 s loop by making several favourable contacts to residues of this loop. The difference observed in the Ki values of the two analogues can be interpreted in terms of desolvation effects, subtle structural changes of protein residues and changes in water structure.

Binding of oxalyl derivatives of β-d-glucopyranosylamine to muscle glycogen phosphorylase b

Five oxalyl derivatives of beta-d-glucopyranosylamine were synthesized as potential inhibitors of glycogen phosphorylase (GP). The compounds 1-4 were competitive inhibitors of rabbit muscle GPb (with respect to alpha-d-glucose-1-phosphate) with K(i) values of 0.2-1.4 mM, while compound 5 was not effective up to a concentration of 10 mM. In order to elucidate the structural basis of their inhibition, we analysed the structures of compounds 1-4 in complex with GPb at 1.93-1.96 Angstrom resolution. The complex structures reveal that the inhibitors can be accommodated at the catalytic site at approximately the same position as alpha-d-glucose and stabilize the T-state conformation of the 280 s loop by making several favourable contacts to Asp283 and Asn284 of this loop. Comparison with the lead compound N-acetyl-beta-d-glucopyranosylamine (6) shows that the hydrogen bonding interaction of the amide nitrogen with the main-chain carbonyl oxygen of His377 is not present in these complexes. The differences observed in the K(i) values of the four analogues can be interpreted in terms of subtle conformational changes of protein residues and shifts of water molecules in the vicinity of the catalytic site, variations in van der Waals interactions, conformational entropy and desolvation effects.

Synthesis of N-(β-d-glucopyranosyl) monoamides of dicarboxylic acids as potential inhibitors of glycogen phosphorylase

O-peracetylated N-(beta-D-glucopyranosyl)imino trimethylphosphorane obtained in situ from 2,3,4,6-tetra-O-acetyl-beta-D-glucopyranosyl azide and PMe3 was reacted with saturated and unsaturated aliphatic and aromatic dicarboxylic acids, or their anhydrides, or monoesters to give the corresponding N-(beta-D-glucopyranosyl) monoamides of dicarboxylic acids or derivatives. The acetyl protecting groups were removed according to the Zemplén protocol to give a series of compounds which showed moderate inhibitory effects against rabbit muscle glycogen phosphorylase b. The best inhibitor was 3-(N-beta-D-glucopyranosyl-carbamoyl)propanoic acid (7) with Ki = 20 microM.

Glycogen phosphorylase inhibitors: a free energy perturbation analysis of glucopyranose spirohydantoin analogues

GP catalyzes the phosphorylation of glycogen to Glc-1-P. Because of its fundamental role in the metabolism of glycogen, GP has been the target for a systematic structure-assisted design of inhibitory compounds, which could be of value in the therapeutic treatment of type 2 diabetes mellitus. The most potent catalytic-site inhibitor of GP identified to date is spirohydantoin of glucopyranose (hydan). In this work, we employ MD free energy simulations to calculate the relative binding affinities for GP of hydan and two spirohydantoin analogues, methyl-hydan and n-hydan, in which a hydrogen atom is replaced by a methyl- or amino group, respectively. The results are compared with the experimental relative affinities of these ligands, estimated by kinetic measurements of the ligand inhibition constants. The calculated binding affinity for methyl-hydan (relative to hydan) is 3.75 +/- 1.4 kcal/mol, in excellent agreement with the experimental value (3.6 +/- 0.2 kcal/mol). For n-hydan, the calculated value is 1.0 +/- 1.1 kcal/mol, somewhat smaller than the experimental result (2.3 +/- 0.1 kcal/mol). A free energy decomposition analysis shows that hydan makes optimum interactions with protein residues and specific water molecules in the catalytic site. In the other two ligands, structural perturbations of the active site by the additional methyl- or amino group reduce the corresponding binding affinities. The computed binding free energies are sensitive to the preference of a specific water molecule for two well-defined positions in the catalytic site. The behavior of this water is analyzed in detail, and the free energy profile for the translocation of the water between the two positions is evaluated. The results provide insights into the role of water molecules in modulating ligand binding affinities. A comparison of the interactions between a set of ligands and their surrounding groups in X-ray structures is often used in the interpretation of binding free energy differences and in guiding the design of new ligands. For the systems in this work, such an approach fails to estimate the order of relative binding strengths, in contrast to the rigorous free energy treatment.

Kinetic and crystallographic studies of glucopyranose spirohydantoin and glucopyranosylamine analogs inhibitors of glycogen phosphorylase

Glycogen phosphorylase (GP) is currently exploited as a target for inhibition of hepatic glycogenolysis under high glucose conditions. Spirohydantoin of glucopyranose and N-acetyl-beta-D-glucopyranosylamine have been identified as the most potent inhibitors of GP that bind at the catalytic site. Four spirohydantoin and three beta-D-glucopyranosylamine analogs have been designed, synthesized and tested for inhibition of GP in kinetic experiments. Depending on the functional group introduced, the K(i) values varied from 16.5 microM to 1200 microM. In order to rationalize the kinetic results, we determined the crystal structures of the analogs in complex with GP. All the inhibitors bound at the catalytic site of the enzyme, by making direct and water-mediated hydrogen bonds with the protein and by inducing minor movements of the side chains of Asp283 and Asn284, of the 280s loop that blocks access of the substrate glycogen to the catalytic site, and changes in the water structure in the vicinity of the site. The differences observed in the Ki values of the analogs can be interpreted in terms of variations in hydrogen bonding and van der Waals interactions, desolvation effects, ligand conformational entropy, and displacement of water molecules on ligand binding to the catalytic site.

Crystallographic studies on two bioisosteric analogues, N-acetyl-beta-D-glucopyranosylamine and N-trifluoroacetyl-beta-D-glucopyranosylamine, potent inhibitors of muscle glycogen phosphorylase

Structure-based inhibitor design has led to the discovery of a number of potent inhibitors of glycogen phosphorylase b (GPb), N-acyl derivatives of beta-D-glucopyranosylamine, that bind at the catalytic site of the enzyme. The first good inhibitor in this class of compounds, N-acetyl-beta-D-glucopyranosylamine (NAG) (K(i) = 32 microM), has been previously characterized by biochemical, biological and crystallographic experiments at 2.3 angstroms resolution. Bioisosteric replacement of the acetyl group by trifluoroacetyl group resulted in an inhibitor, N-trifluoroacetyl-beta-D-glucopyranosylamine (NFAG), with a K(i) = 75 microM. To elucidate the structural basis of its reduced potency, we determined the ligand structure in complex with GPb at 1.8 angstroms resolution. To compare the binding mode of N-trifluoroacetyl derivative with that of the lead molecule, we also determined the structure of GPb-NAG complex at a higher resolution (1.9 angstroms). NFAG can be accommodated in the catalytic site of T-state GPb at approximately the same position as that of NAG and stabilize the T-state conformation of the 280 s loop by making several favourable contacts to Asn284 of this loop. The difference observed in the K(i) values of the two analogues can be interpreted in terms of subtle conformational changes of protein residues and shifts of water molecules in the vicinity of the catalytic site, variations in van der Waals interaction, and desolvation effects.

Kinetic and crystallographic studies on 2-(beta-D-glucopyranosyl)-5-methyl-1, 3, 4-oxadiazole, -benzothiazole, and -benzimidazole, inhibitors of muscle glycogen phosphorylase b. Evidence for a new binding site

In an attempt to identify leads that would enable the design of inhibitors with enhanced affinity for glycogen phosphorylase (GP), that might control hyperglycaemia in type 2 diabetes, three new analogs of beta-D-glucopyranose, 2-(beta-D-glucopyranosyl)-5-methyl-1, 3, 4-oxadiazole, -benzothiazole, and -benzimidazole were assessed for their potency to inhibit GPb activity. The compounds showed competitive inhibition (with respect to substrate Glc-1-P) with K(i) values of 145.2 (+/-11.6), 76 (+/-4.8), and 8.6 (+/-0.7) muM, respectively. In order to establish the mechanism of this inhibition, crystallographic studies were carried out and the structures of GPb in complex with the three analogs were determined at high resolution (GPb-methyl-oxadiazole complex, 1.92 A; GPb-benzothiazole, 2.10 A; GPb-benzimidazole, 1.93 A). The complex structures revealed that the inhibitors can be accommodated in the catalytic site of T-state GPb with very little change of the tertiary structure, and provide a rationalization for understanding variations in potency of the inhibitors. In addition, benzimidazole bound at the new allosteric inhibitor or indole binding site, located at the subunit interface, in the region of the central cavity, and also at a novel binding site, located at the protein surface, far removed (approximately 32 A) from the other binding sites, that is mostly dominated by the nonpolar groups of Phe202, Tyr203, Val221, and Phe252.