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

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.

Pyranoid Spirosugars as Enzyme Inhibitors

Background:

Pyranoid spirofused sugar derivatives represent a class of compounds with a significant impact in the literature. From the structural point of view, the rigidity inferred by the spirofused entity has made these compounds object of interest mainly as enzymatic inhibitors, in particular, carbohydrate processing enzymes. Among them glycogen phosphorylase and sodium glucose co-transporter 2 are important target enzymes for diverse pathological states. Most of the developed compounds present the spirofused entity at the C1 position of the sugar moiety; nevertheless, spirofused entities can also be found at other sugar ring positions. The main spirofused entities encountered are spiroacetals/thioacetals, spiro-hydantoin and derivatives, spiro-isoxazolines, spiro-aminals, spiro-lactams, spiro-oxathiazole and spiro-oxazinanone, but also others are present.

Objectives:

The present review focuses on the most explored synthetic strategies for the preparation of this class of compounds, classified according to the position and structure of the spirofused moiety on the pyranoid scaffold. Moreover, the structures are correlated to their main biological activities or to their role as chiral auxiliaries.

Conclusion:

It is clear from the review that, among the different derivatives, the spirofused structures at position C1 of the pyranoid scaffold are the most represented and possess the most relevant enzymatic inhibitor activities. Nevertheless, great efforts have been devoted to the introduction of the spirofused entity also in the other positions, mainly for the preparation of biologically active compounds but also for the synthesis of chiral auxiliaries useful in asymmetric reactions; examples of such auxiliaries are the spirofused chiral 1,3-oxazolidin-2-ones and 1,3-oxazolidine-2-thiones.

Identification of Antifungal Targets Based on Computer Modeling

Aspergillus fumigatus is a saprophytic, cosmopolitan fungus that attacks patients with a weak immune system. A rational solution against fungal infection aims to manipulate fungal metabolism or to block enzymes essential for Aspergillus survival. Here we discuss and compare different bioinformatics approaches to analyze possible targeting strategies on fungal-unique pathways. For instance, phylogenetic analysis reveals fungal targets, while domain analysis allows us to spot minor differences in protein composition between the host and fungi. Moreover, protein networks between host and fungi can be systematically compared by looking at orthologs and exploiting information from host⁻pathogen interaction databases. Further data—such as knowledge of a three-dimensional structure, gene expression data, or information from calculated metabolic fluxes—refine the search and rapidly put a focus on the best targets for antimycotics. We analyzed several of the best targets for application to structure-based drug design. Finally, we discuss general advantages and limitations in identification of unique fungal pathways and protein targets when applying bioinformatics tools.

A large-scale test of free-energy simulation estimates of protein-ligand binding affinities

We have performed a large-scale test of alchemical perturbation calculations with the Bennett acceptance-ratio (BAR) approach to estimate relative affinities for the binding of 107 ligands to 10 different proteins. Employing 20-Å truncated spherical systems and only one intermediate state in the perturbations, we obtain an error of less than 4 kJ/mol for 54% of the studied relative affinities and a precision of 0.5 kJ/mol on average. However, only four of the proteins gave acceptable errors, correlations, and rankings. The results could be improved by using nine intermediate states in the simulations or including the entire protein in the simulations using periodic boundary conditions. However, 27 of the calculated affinities still gave errors of more than 4 kJ/mol, and for three of the proteins the results were not satisfactory. This shows that the performance of BAR calculations depends on the target protein and that several transformations gave poor results owing to limitations in the molecular-mechanics force field or the restricted sampling possible within a reasonable simulation time. Still, the BAR results are better than docking calculations for most of the proteins.

Challenges in the determination of the binding modes of non-standard ligands in X-ray crystal complexes

Despite its central role in structure based drug design the determination of the binding mode (position, orientation and conformation in addition to protonation and tautomeric states) of small heteromolecular ligands in protein:ligand complexes based on medium resolution X-ray diffraction data is highly challenging. In this perspective we demonstrate how a combination of molecular dynamics simulations and free energy (FE) calculations can be used to correct and identify thermodynamically stable binding modes of ligands in X-ray crystal complexes. The consequences of inappropriate ligand structure, force field and the absence of electrostatics during X-ray refinement are highlighted. The implications of such uncertainties and errors for the validation of virtual screening and fragment-based drug design based on high throughput X-ray crystallography are discussed with possible solutions and guidelines.

Structural glycobiology: a game of snakes and ladders

Oligo- and polysaccharides are infamous for being extremely flexible molecules, populating a series of well-defined rotational isomeric states under physiological conditions. Characterization of this heterogeneous conformational ensemble has been a major obstacle impeding high-resolution structure determination of carbohydrates and acting as a bottleneck in the effort to understand the relationship between the carbohydrate structure and function. This challenge has compelled the field to develop and apply theoretical and experimental methods that can explore conformational ensembles by both capturing and deconvoluting the structural and dynamic properties of carbohydrates. This review focuses on computational approaches that have been successfully used in combination with experiment to detail the three-dimensional structure of carbohydrates in a solution and in a complex with proteins. In addition, emerging experimental techniques for three-dimensional structural characterization of carbohydrate-protein complexes and future challenges in the field of structural glycobiology are discussed. The review is divided into five sections: (1) The complexity and plasticity of carbohydrates, (2) Predicting carbohydrate-protein interactions, (3) Calculating relative and absolute binding free energies for carbohydrate-protein complexes, (4) Emerging and evolving techniques for experimental characterization of carbohydrate-protein structures, and (5) Current challenges in structural glycoscience.

Synthesis and glycogen phosphorylase inhibitor activity of 2,3-dihydrobenzo[1,4]dioxin derivatives

Novel 5-benzyl and 5-benzylidene-thiazolidine-2,4-diones carrying 2,3-dihydrobenzo[1,4]dioxin pharmacophore were synthesized and their glycogen phosphorylase inhibitor activity was also studied.

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.