Structural, Kinetic, and Thermodynamic Analysis of Glucoimidazole-Derived Glycosidase Inhibitors†,‡

Binding interactions of hydrophobically-modified flavonols with β-glucosidase: fluorescence spectroscopy and molecular modelling study

Natural flavonoids are capable of inhibiting glucosidase activity, so they can be used for treating diabetes mellitus and hypertension. However, molecular-level details of their interactions with glucosidase enzymes remain poorly understood. This paper describes the synthesis and spectral characterization of a series of fluorescent flavonols and their interaction with the β-glucosidase enzyme. To tune flavonol-enzyme interaction modes and affinity, we introduced different polar halogen-containing groups or bulky aromatic/alkyl substituents in the peripheral 2-aryl ring of a flavonol moiety. Using fluorescence spectroscopy methods in combination with molecular docking and molecular dynamics simulations, we examined the binding affinity and identified probe binding patterns, which are critical for steric blockage of the key catalytic residues of the enzyme. Using a fluorescent assay, we demonstrated that the binding of flavonol 2e to β-glucosidase decreased its enzymatic activity up to 3.5 times. In addition, our molecular docking and all-atom molecular dynamics simulations suggest that the probe binding is driven by hydrophobic interactions with aromatic Trp and Tyr residues within the catalytic glycone binding pockets of β-glucosidase. Our study provides a new insight into structure-property relations for flavonol-protein interactions, which govern their enzyme binding, and outlines a framework for a rational design of new flavonol-based potent inhibitors for β-glucosidases.

Reaction Mechanism of Glycoside Hydrolase Family 116 Utilizes Perpendicular Protonation

Retaining glycoside hydrolases use acid/base catalysis with an enzymatic acid/base protonating the glycosidic bond oxygen to facilitate leaving-group departure alongside attack by a catalytic nucleophile to form a covalent intermediate. Generally, this acid/base protonates the oxygen laterally with respect to the sugar ring, which places the catalytic acid/base and nucleophile carboxylates within about 4.5-6.5 Å of each other. However, in glycoside hydrolase (GH) family 116, including disease-related human acid β-glucosidase 2 (GBA2), the distance between the catalytic acid/base and the nucleophile is around 8 Å (PDB: 5BVU) and the catalytic acid/base appears to be above the plane of the pyranose ring, rather than being lateral to that plane, which could have catalytic consequences. However, no structure of an enzyme-substrate complex is available for this GH family. Here, we report the structures of Thermoanaerobacterium xylanolyticum β-glucosidase (TxGH116) D593N acid/base mutant in complexes with cellobiose and laminaribiose and its catalytic mechanism. We confirm that the amide hydrogen bonding to the glycosidic oxygen is in a perpendicular rather than lateral orientation. Quantum mechanics/molecular mechanics (QM/MM) simulations of the glycosylation half-reaction in wild-type TxGH116 indicate that the substrate binds with the nonreducing glucose residue in an unusual relaxed 4C1 chair at the -1 subsite. Nevertheless, the reaction can still proceed through a 4H3 half-chair transition state, as in classical retaining β-glucosidases, as the catalytic acid D593 protonates the perpendicular electron pair. The glucose C6OH is locked in a gauche, trans orientation with respect to the C5-O5 and C4-C5 bonds to facilitate perpendicular protonation. These data imply a unique protonation trajectory in Clan-O glycoside hydrolases, which has strong implications for the design of inhibitors specific to either lateral protonators, such as human GBA1, or perpendicular protonators, such as human GBA2.

Influence of Side Chain Conformation on the Activity of Glycosidase Inhibitors

Substrate side chain conformation impacts reactivity during glycosylation and glycoside hydrolysis and is restricted by many glycosidases and glycosyltransferases during catalysis. We show that the side chains of gluco and manno iminosugars can be restricted to predominant conformations by strategic installation of a methyl group. Glycosidase inhibition studies reveal that iminosugars with the gauche,gauche side chain conformations are 6- to 10-fold more potent than isosteric compounds with the gauche,trans conformation; a manno-configured iminosugar with the gauche,gauche conformation is a 27-fold better inhibitor than 1-deoxymannojirimycin. The results are discussed in terms of the energetic benefits of preorganization, particularly when in synergy with favorable hydrophobic interactions. The demonstration that inhibitor side chain preorganization can favorably impact glycosidase inhibition paves the way for improved inhibitor design through conformational preorganization.

Design and synthesis of novel benzimidazole-iminosugars linked a substituted phenyl group and their inhibitory activities against β-glucosidase

A series of novel benzimidazole-iminosugars linked a (substuituted) phenyl group on benzene ring of benzimidazole 5(a-p) and 6(a-p) have been rationally designed and conveniently synthesized through Suzuki coupling reaction in high yields. All compounds have been evaluated for their inhibitory activities against β-glucosidase (almond). Six compounds 5d, 6d, 6e, 6i, 6n, and 6p showed more significant inhibitory activities with IC₅₀ values in the range of 0.03-0.08 μM, almost 10-fold improved than that of the parent analogue 4, and much higher than that of the positive control castanospermine. The additional phenyl ring and the electron donating groups on it would be beneficial for the activity. Compounds 6d, 6n, and 4 had been chosen to be tested for their inhibition types against β-glucosidase. Interestingly, three compounds have different inhibition types although they had very similar structure. Their K_i values were calculated to be 0.02 ± 0.01 μM, 0.02 ± 0.01 μM, and 0.66 ± 0.14 μM, respectively. The equilibrium dissociation constant (K_D) for 6d, 6n, and 4 and β-glucosidase was 0.04 μM, 0.03 μM and 0.45 μM by the ITC-based assay, respectively. Molecular docking work suggests that such benzimidazole-iminosugars derivatives might bind to the active site of β-glucosidase mainly through hydrogen bonds, the additional phenyl ring towards the solvent-exposed region played an important effect on their inhibitory activity against β-glucosidase.

Substrate specificity of glycoside hydrolase family 1 β-glucosidase AtBGlu42 from Arabidopsis thaliana and its molecular mechanism

Plants possess many glycoside hydrolase family 1 (GH1) β-glucosidases, which physiologically function in cell wall metabolism and activation of bioactive substances, but most remain uncharacterized. One GH1 isoenzyme AtBGlu42 in Arabidopsis thaliana has been identified to hydrolyze scopolin using the gene deficient plants, but no enzymatic properties were obtained. Its sequence similarity to another functionally characterized enzyme Os1BGlu4 in rice suggests that AtBGlu42 also acts on oligosaccharides. Here, we show that the recombinant AtBGlu42 possesses high kcat/Km not only on scopolin, but also on various β-glucosides, cellooligosaccharides, and laminarioligosaccharides. Of the cellooligosaccharides, cellotriose was the most preferred. The crystal structure, determined at 1.7 Å resolution, suggests that Arg342 gives unfavorable binding to cellooligosaccharides at subsite +3. The mutants R342Y and R342A showed the highest preference on cellotetraose or cellopentaose with increased affinities at subsite +3, indicating that the residues at this position have an important role for chain length specificity.

Cloning, biochemical characterization and molecular docking of novel thermostable β-glucosidase BglA9 from Anoxybacillus ayderensis A9 and its application in de-glycosylation of Polydatin

This study reports a novel BglA9 gene of 1345 bp encoding β-glucosidase from Anoxybacillus ayderensis A9, which was amplified and expressed in E. coli BL21 (DE3): pLysS cells, purified with Ni-NTA column having molecular weight of 52.6 kDa and was used in the bioconversion of polydatin to resveratrol. The kinetic parameters values using pNPG as substrate were K_m (0.28 mM), V_max (43.8 μmol/min/mg), k_cat (38.43 s^-1) and k_cat/K_m (135.5 s^-1 mM^-1). The BglA9 was active in a broad pH range and had an activity half-life around 24 h at 50 °C. The de-glycosylation efficiency of BglA9 for polydatin was determined by estimating the amount of glucose released after enzymatic reaction by a dinitrosalicylic acid (DNS) assay. The kinetic parameters of BglA9 for polydatin were 5.5 mM, 20.84 μmol/min/mg, 18.28 s^-1and 3.27 s^-1 mM^-1 for K_m, V_max, k_cat, and k_cat/K_m values, respectively. The K_i value for glucose was determined to be 1.7 M. The residues Gln19, His120, Glu355, Glu409, Glu178, Asn222 may play a crucial role in the deglycosylation as revealed by the 3D structure of enzyme docked with polydatin.

Molecular docking and molecular dynamics studies of bioactive compounds contained in noni fruit (Morinda citrifolia L.) against human pancreatic α-amylase

Human pancreatic α-amylase inhibition is currently a promising therapeutic target against type 2 diabetes (DMT2) because it can reduce aggressive digestion of carbohydrates into absorbable monosaccharides. In Indonesia, medicinal plants, e.g. Morinda citrifolia fruit, have been empirically utilized as a blood-sugar reducer, however, the inhibitory activity of compounds in this plant against human pancreatic α-amylase is still limited or none. Therefore, this study aimed to test the interaction of 7 compounds (americanin, asperulosidic acid, damnacanthal, quercetin, rutin, scopoletin, and ursolic acid) contained in noni fruit against human pancreatic α-amylase by molecular docking and molecular dynamics and compared their binding modes with that of acarbose. Results of the molecular docking simulation indicated that the ursolic acid compound possesses the best binding energy (-8.58 kcal/mol) and comparable to that of acarbose (-8.59 kcal/mol). The molecular dynamics study at 100 ns simulation, the values of RMSD, RMSF, the radius of gyration (Rg), the solvent-accessible surface area (SASA), principal component analysis (PCA), and MM-PBSA binding free energy were stable and identical to those of acarbose. It could be concluded that ursolic acid might be potential in inhibiting human pancreatic α-amylase, thus, potential to be developed as an anti-DMT2 drug candidate. Communicated by Ramaswamy H. Sarma.

Forty years of study on the thermostable β-glycosidase from S. solfataricus: Production, biochemical characterization and biotechnological applications

The aim of this paper is to make the point on the fortieth years study on the β-glycosidase from Sulfolobus solfataricus. This enzyme represents one of the thermophilic biocatalysts, which is more extensively studied as witnessed by the numerous literature reports available since 1980. Comprehensive biochemical studies highlighted its broad substrate specificity for β-d-galacto-, gluco-, and fuco-sides and also showed its remarkable exo-glucosidase and transglycosidase activities. The enzyme demonstrated to be active and stable over a wide range of temperature and pHs, withstanding to several drastic conditions comprising solvents and detergents. Over the years, a great deal of studies were focused on its homotetrameric tridimensional structure, elucidating several structural features involved in the enzyme stability, such as ion pairs and post-translational modifications. Several β-glycosidase mutants were produced in the years in order to understand its peculiar behavior in extreme conditions and/or to improve its functional properties. The β-glycosidase overproduction was also afforded reporting numerous studies dealing with its production in the mesophilic host Escherichia coli, Saccharomyces cerevisiae, Pichia pastoris, and Lactococcus lactis. Relevant applications in food, beverages, bioenergy, pharmaceuticals, and nutraceutical fields of this enzyme, both in free and immobilized forms, highlighted its biotechnological relevance.

Structural Basis of Specific Glucoimidazole and Mannoimidazole Binding by Os3BGlu7

β-Glucosidases and β-mannosidases hydrolyze substrates that differ only in the epimer of the nonreducing terminal sugar moiety, but most such enzymes show a strong preference for one activity or the other. Rice Os3BGlu7 and Os7BGlu26 β-glycosidases show a less strong preference, but Os3BGlu7 and Os7BGlu26 prefer glucosides and mannosides, respectively. Previous studies of crystal structures with glucoimidazole (GIm) and mannoimidazole (MIm) complexes and metadynamic simulations suggested that Os7BGlu26 hydrolyzes mannosides via the B_2,5 transition state (TS) conformation preferred for mannosides and glucosides via their preferred ⁴H₃/⁴E TS conformation. However, MIm is weakly bound by both enzymes. In the present study, we found that MIm was not bound in the active site of crystallized Os3BGlu7, but GIm was tightly bound in the -1 subsite in a ⁴H₃/⁴E conformation via hydrogen bonds with the surrounding residues. One-microsecond molecular dynamics simulations showed that GIm was stably bound in the Os3BGlu7 active site and the glycone-binding site with little distortion. In contrast, MIm initialized in the B_2,5 conformation rapidly relaxed to a E₃/⁴H₃ conformation and moved out into a position in the entrance of the active site, where it bound more stably despite making fewer interactions. The lack of MIm binding in the glycone site in protein crystals and simulations implies that the energy required to distort MIm to the B_2,5 conformation for optimal active site residue interactions is sufficient to offset the energy of those interactions in Os3BGlu7. This balance between distortion and binding energy may also provide a rationale for glucosidase versus mannosidase specificity in plant β-glycosidases.

Probing the role of an invariant active site His in family GH1 β-glycosidases

The reaction mechanism of glycoside hydrolases belonging to family 1 (GH1) of carbohydrate-active enzymes classification, hydrolysing β-O-glycosidic bonds, is well characterised. This family includes several thousands of enzymes with more than 20 different EC numbers depending on the sugar glycone recognised as substrate. Most GH1 β-glycosidases bind their substrates with similar specificity through invariant amino acid residues. Despite extensive studies, the clear identification of the roles played by each of these residues in the recognition of different glycones is not always possible. We demonstrated here that a histidine residue, completely conserved in the active site of the enzymes of this family, interacts with the C2-OH of the substrate in addition to the C3-OH as previously shown by 3 D-structure determination.

Effects of Random Mutagenesis and In Vivo Selection on the Specificity and Stability of a Thermozyme

Factors that give enzymes stability, activity, and substrate recognition result from the combination of few weak molecular interactions, which can be difficult to study through rational protein engineering approaches. We used irrational random mutagenesis and in vivo selection to test if a β-glycosidase from the thermoacidophile Saccharolobus solfataricus (Ssβ-gly) could complement an Escherichia coli strain unable to grow on lactose. The triple mutant of Ssβ-gly (S26L, P171L, and A235V) was more active than the wild type at 85 °C, inactivated at this temperature almost 300-fold quicker, and showed a 2-fold higher kcat on galactosides. The three mutations, which were far from the active site, were analyzed to test their effect at the structural level. Improved activity on galactosides was induced by the mutations. The S26L and P171L mutations destabilized the enzyme through the removal of a hydrogen bond and increased flexibility of the peptide backbone, respectively. However, the flexibility added by S26L mutation improved the activity at T > 60 °C. This study shows that random mutagenesis and biological selection allowed the identification of residues that are critical in determining thermal activity, stability, and substrate recognition.

Gluco-1 H-imidazole: A New Class of Azole-Type β-Glucosidase Inhibitor

Gluco-azoles competitively inhibit glucosidases by transition-state mimicry and their ability to interact with catalytic acid residues in glucosidase active sites. We noted that no azole-type inhibitors described, to date, possess a protic nitrogen characteristic for 1 H-imidazoles. Here, we present gluco-1 H-imidazole, a gluco-azole bearing a 1 H-imidazole fused to a glucopyranose-configured cyclitol core, and three close analogues as new glucosidase inhibitors. All compounds inhibit human retaining β-glucosidase, GBA1, with the most potent ones inhibiting this enzyme (deficient in Gaucher disease) on a par with glucoimidazole. None inhibit glucosylceramide synthase, cytosolic β-glucosidase GBA2 or α-glucosidase GAA. Structural, physical and computational studies provide first insights into the binding mode of this conceptually new class of retaining β-glucosidase inhibitors.

Improved conversion of ginsenoside Rb1 to compound K by semi-rational design of Sulfolobus solfataricus β-glycosidase

Ginsenoside compound K has been used as a key nutritional and cosmetic component because of its anti-fatigue and skin anti-aging effects. β-Glycosidase from Sulfolobus solfataricus (SS-BGL) is known as the most efficient enzyme for compound K production. The hydrolytic pathway from ginsenoside Rb₁ to compound K via Rd and F₂ is the most important because Rb₁ is the most abundant component in ginseng extract. However, the enzymatic conversion of ginsenoside Rd to F₂ is a limiting step in the hydrolytic pathway because of the relatively low activity for Rd. A V209 residue obtained from error-prone PCR was related to Rd-hydrolyzing activity, and a docking pose showing an interaction with Val209 was selected from numerous docking poses. W361F was obtained by rational design using the docking pose that exhibited 4.2-fold higher activity, 3.7-fold higher catalytic efficiency, and 3.1-fold lower binding energy for Rd than the wild-type enzyme, indicating that W361F compensated for the limiting step. W361F completely converted Rb₁ to compound K with a productivity of 843 mg l⁻¹ h⁻¹ in 80 min, and showed also 7.4-fold higher activity for the flavanone, hesperidin, than the wild-type enzyme. Therefore, the W361F variant SS-BGL can be useful for hydrolysis of other glycosides as well as compound K production from Rb₁, and semi-rational design is a useful tool for enhancing hydrolytic activity of β-glycosidase.

Development of inhibitors as research tools for carbohydrate-processing enzymes

Carbohydrates, which are present in all domains of life, play important roles in a host of cellular processes. These ubiquitous biomolecules form highly diverse and often complex glycan structures without the aid of a template. The carbohydrate structures are regulated solely by the location and specificity of the enzymes responsible for their synthesis and degradation. These enzymes, glycosyltransferases and glycoside hydrolases, need to be functionally well characterized in order to investigate the structure and function of glycans. The use of enzyme inhibitors, which target a particular enzyme, can significantly aid this understanding, and may also provide insights into therapeutic applications. The present article describes some of the approaches used to design and develop enzyme inhibitors as tools for investigating carbohydrate-processing enzymes.

Designing allosteric control into enzymes by chemical rescue of structure

Ligand-dependent activity has been engineered into enzymes for purposes ranging from controlling cell morphology to reprogramming cellular signaling pathways. Where these successes have typically fused a naturally allosteric domain to the enzyme of interest, here we instead demonstrate an approach for designing a de novo allosteric effector site directly into the catalytic domain of an enzyme. This approach is distinct from traditional chemical rescue of enzymes in that it relies on disruption and restoration of structure, rather than active site chemistry, as a means to achieve modulate function. We present two examples, W33G in a β-glycosidase enzyme (β-gly) and W492G in a β-glucuronidase enzyme (β-gluc), in which we engineer indole-dependent activity into enzymes by removing a buried tryptophan side chain that serves as a buttress for the active site architecture. In both cases, we observe a loss of function, and in both cases we find that the subsequent addition of indole can be used to restore activity. Through a detailed analysis of β-gly W33G kinetics, we demonstrate that this rescued enzyme is fully functionally equivalent to the corresponding wild-type enzyme. We then present the apo and indole-bound crystal structures of β-gly W33G, which together establish the structural basis for enzyme inactivation and rescue. Finally, we use this designed switch to modulate β-glycosidase activity in living cells using indole. Disruption and recovery of protein structure may represent a general technique for introducing allosteric control into enzymes, and thus may serve as a starting point for building a variety of bioswitches and sensors.

Aglycone specificity of Thermotoga neapolitana β-glucosidase 1A modified by mutagenesis, leading to increased catalytic efficiency in quercetin-3-glucoside hydrolysis

BACKGROUND

The thermostable β-glucosidase (TnBgl1A) from Thermotoga neapolitana is a promising biocatalyst for hydrolysis of glucosylated flavonoids and can be coupled to extraction methods using pressurized hot water. Hydrolysis has however been shown to be dependent on the position of the glucosylation on the flavonoid, and e.g. quercetin-3-glucoside (Q3) was hydrolysed slowly. A set of mutants of TnBgl1A were thus created to analyse the influence on the kinetic parameters using the model substrate para-nitrophenyl-β-D-glucopyranoside (pNPGlc), and screened for hydrolysis of Q3.

RESULTS

Structural analysis pinpointed an area in the active site pocket with non-conserved residues between specificity groups in glycoside hydrolase family 1 (GH1). Three residues in this area located on β-strand 5 (F219, N221, and G222) close to sugar binding sub-site +2 were selected for mutagenesis and amplified in a protocol that introduced a few spontaneous mutations. Eight mutants (four triple: F219L/P165L/M278I, N221S/P165L/M278I, G222Q/P165L/M278I, G222Q/V203M/K214R, two double: F219L/K214R, N221S/P342L and two single: G222M and N221S) were produced in E. coli, and purified to apparent homogeneity. Thermostability, measured as Tm by differential scanning calorimetry (101.9°C for wt), was kept in the mutated variants and significant decrease (ΔT of 5-10°C) was only observed for the triple mutants. The exchanged residue(s) in the respective mutant resulted in variations in KM and turnover. The KM-value was only changed in variants mutated at position 221 (N221S) and was in all cases monitored as a 2-3 × increase for pNPGlc, while the KM decreased a corresponding extent for Q3.Turnover was only significantly changed using pNPGlc, and was decreased 2-3 × in variants mutated at position 222, while the single, double and triple mutated variants carrying a mutation at position 221 (N221S) increased turnover up to 3.5 × compared to the wild type. Modelling showed that the mutation at position 221, may alter the position of N291 resulting in increased hydrogen bonding of Q3 (at a position corresponding to the +1 subsite) which may explain the decrease in KM for this substrate.

CONCLUSION

These results show that residues at the +2 subsite are interesting targets for mutagenesis and mutations at these positions can directly or indirectly affect both KM and turnover. An affinity change, leading to a decreased KM, can be explained by an altered position of N291, while the changes in turnover are more difficult to explain and may be the result of smaller conformational changes in the active site.

Imino sugars and glycosyl hydrolases: historical context, current aspects, emerging trends

Forty years of discoveries and research on imino sugars, which are carbohydrate analogues having a basic nitrogen atom instead of oxygen in the sugar ring and, acting as potent glycosidase inhibitors, have made considerable impact on our contemporary understanding of glycosidases. Imino sugars have helped to elucidate the catalytic machinery of glycosidases and have refined our methods and concepts of utilizing them. A number of new aspects have emerged for employing imino sugars as pharmaceutical compounds, based on their profound effects on metabolic activities in which glycosidases are involved. From the digestion of starch to the fight against viral infections, from research into malignant diseases to potential improvements in hereditary storage disorders, glycosidase action and inhibition are essential issues. This account aims at combining general developments with a focus on some niches where imino sugars have become useful tools for glycochemistry and glycobiology.

β-Glucosidases

β-Glucosidases (3.2.1.21) are found in all domains of living organisms, where they play essential roles in the removal of nonreducing terminal glucosyl residues from saccharides and glycosides. β-Glucosidases function in glycolipid and exogenous glycoside metabolism in animals, defense, cell wall lignification, cell wall β-glucan turnover, phytohormone activation, and release of aromatic compounds in plants, and biomass conversion in microorganisms. These functions lead to many agricultural and industrial applications. β-Glucosidases have been classified into glycoside hydrolase (GH) families GH1, GH3, GH5, GH9, and GH30, based on their amino acid sequences, while other β-glucosidases remain to be classified. The GH1, GH5, and GH30 β-glucosidases fall in GH Clan A, which consists of proteins with (β/α)(8)-barrel structures. In contrast, the active site of GH3 enzymes comprises two domains, while GH9 enzymes have (α/α)(6) barrel structures. The mechanism by which GH1 enzymes recognize and hydrolyze substrates with different specificities remains an area of intense study.

Structural studies of phosphoglucose isomerase from Mycobacterium tuberculosis H37Rv

Phosphoglucose isomerase (PGI) plays a key role in both glycolysis and gluconeogenesis inside the cell, whereas outside the cell it exhibits cytokine properties. PGI is also known to act as an autocrine motility factor, a neuroleukin agent and a differentiation and maturation mediator. Here, the first crystal structure of PGI from Mycobacterium tuberculosis H37Rv (Mtb) is reported. The structure was refined at 2.25 A resolution and revealed the presence of one molecule in the asymmetric unit with two globular domains. As known previously, the active site of Mtb PGI contains conserved residues including Glu356, Glu216 and His387 (where His387 is from the neighbouring molecule). The crystal structure of Mtb PGI was observed to be rather more similar to human PGI than other nonbacterial PGIs, with only a few differences being detected in the loops, arm and hook regions of the human and Mtb PGIs, suggesting that the M. tuberculosis enzyme uses the same enzyme mechanism.

Glycosidase inhibition: assessing mimicry of the transition state

Glycoside hydrolases, the enzymes responsible for hydrolysis of the glycosidic bond in di-, oligo- and polysaccharides, and glycoconjugates, are ubiquitous in Nature and fundamental to existence. The extreme stability of the glycosidic bond has meant these enzymes have evolved into highly proficient catalysts, with an estimated 10(17) fold rate enhancement over the uncatalysed reaction. Such rate enhancements mean that enzymes bind the substrate at the transition state with extraordinary affinity; the dissociation constant for the transition state is predicted to be 10(-22) M. Inhibition of glycoside hydrolases has widespread application in the treatment of viral infections, such as influenza and HIV, lysosomal storage disorders, cancer and diabetes. If inhibitors are designed to mimic the transition state, it should be possible to harness some of the transition state affinity, resulting in highly potent and specific drugs. Here we examine a number of glycosidase inhibitors which have been developed over the past half century, either by Nature or synthetically by man. A number of criteria have been proposed to ascertain which of these inhibitors are true transition state mimics, but these features have only be critically investigated in a very few cases.

Insight into a strategy for attenuating AmpC-mediated beta-lactam resistance: structural basis for selective inhibition of the glycoside hydrolase NagZ

NagZ is an exo-N-acetyl-beta-glucosaminidase, found within Gram-negative bacteria, that acts in the peptidoglycan recycling pathway to cleave N-acetylglucosamine residues off peptidoglycan fragments. This activity is required for resistance to cephalosporins mediated by inducible AmpC beta-lactamase. NagZ uses a catalytic mechanism involving a covalent glycosyl enzyme intermediate, unlike that of the human exo-N-acetyl-beta-glucosaminidases: O-GlcNAcase and the beta-hexosaminidase isoenzymes. These latter enzymes, which remove GlcNAc from glycoconjugates, use a neighboring-group catalytic mechanism that proceeds through an oxazoline intermediate. Exploiting these mechanistic differences we previously developed 2-N-acyl derivatives of O-(2-acetamido-2-deoxy-D-glucopyranosylidene)amino-N-phenylcarbamate (PUGNAc), which selectively inhibits NagZ over the functionally related human enzymes and attenuate antibiotic resistance in Gram-negatives that harbor inducible AmpC. To understand the structural basis for the selectivity of these inhibitors for NagZ, we have determined its crystallographic structure in complex with N-valeryl-PUGNAc, the most selective known inhibitor of NagZ over both the human beta-hexosaminidases and O-GlcNAcase. The selectivity stems from the five-carbon acyl chain of N-valeryl-PUGNAc, which we found ordered within the enzyme active site. In contrast, a structure determination of a human O-GlcNAcase homologue bound to a related inhibitor N-butyryl-PUGNAc, which bears a four-carbon chain and is selective for both NagZ and O-GlcNAcase over the human beta-hexosamnidases, reveals that this inhibitor induces several conformational changes in the active site of this O-GlcNAcase homologue. A comparison of these complexes, and with the human beta-hexosaminidases, reveals how selectivity for NagZ can be engineered by altering the 2-N-acyl substituent of PUGNAc to develop inhibitors that repress AmpC mediated beta-lactam resistance.

Application of the PM6 semi-empirical method to modeling proteins enhances docking accuracy of AutoDock

BACKGROUND

Molecular docking methods are commonly used for predicting binding modes and energies of ligands to proteins. For accurate complex geometry and binding energy estimation, an appropriate method for calculating partial charges is essential. AutoDockTools software, the interface for preparing input files for one of the most widely used docking programs AutoDock 4, utilizes the Gasteiger partial charge calculation method for both protein and ligand charge calculation. However, it has already been shown that more accurate partial charge calculation - and as a consequence, more accurate docking- can be achieved by using quantum chemical methods. For docking calculations quantum chemical partial charge calculation as a routine was only used for ligands so far. The newly developed Mozyme function of MOPAC2009 allows fast partial charge calculation of proteins by quantum mechanical semi-empirical methods. Thus, in the current study, the effect of semi-empirical quantum-mechanical partial charge calculation on docking accuracy could be investigated.

RESULTS

The docking accuracy of AutoDock 4 using the original AutoDock scoring function was investigated on a set of 53 protein ligand complexes using Gasteiger and PM6 partial charge calculation methods. This has enabled us to compare the effect of the partial charge calculation method on docking accuracy utilizing AutoDock 4 software. Our results showed that the docking accuracy in regard to complex geometry (docking result defined as accurate when the RMSD of the first rank docking result complex is within 2 A of the experimentally determined X-ray structure) significantly increased when partial charges of the ligands and proteins were calculated with the semi-empirical PM6 method. Out of the 53 complexes analyzed in the course of our study, the geometry of 42 complexes were accurately calculated using PM6 partial charges, while the use of Gasteiger charges resulted in only 28 accurate geometries. The binding affinity estimation was not influenced by the partial charge calculation method - for more accurate binding affinity prediction development of a new scoring function for AutoDock is needed.

CONCLUSION

Our results demonstrate that the accuracy of determination of complex geometry using AutoDock 4 for docking calculation greatly increases with the use of quantum chemical partial charge calculation on both the ligands and proteins.

Mechanistic insights into glycosidase chemistry

The enzymatic hydrolysis of the glycosidic bond continues to gain importance, reflecting the critically important roles complex glycans play in health and disease as well as the rekindled interest in enzymatic biomass conversion. Recent advances include the broadening of our understanding of enzyme reaction coordinates, through both computational and structural studies, improved understanding of enzyme inhibition through transition state mimicry and fascinating insights into mechanism yielded by physical organic chemistry approaches.

GlcNAcstatin: a picomolar, selective O-GlcNAcase inhibitor that modulates intracellular O-glcNAcylation levels

Many phosphorylation signal transduction pathways in the eukaryotic cell are modulated by posttranslational modification of specific serines/threonines with N-acetylglucosamine (O-GlcNAc). Levels of O-GlcNAc on key proteins regulate biological processes as diverse as the cell cycle, insulin signaling, and protein degradation. The two enzymes involved in this dynamic and abundant modification are the O-GlcNAc transferase and O-GlcNAcase. Structural data have recently revealed that the O-GlcNAcase possesses an active site with significant structural similarity to that of the human lysosomal hexosaminidases HexA/HexB. PUGNAc, an O-GlcNAcase inhibitor widely used to raise levels of O-GlcNAc in human cell lines, also inhibits these hexosaminidases. Here, we have exploited recent structural information of an O-GlcNAcase-PUGNAc complex to design and synthesize a glucoimidazole-based inhibitor, GlcNAcstatin, which is a 5 pM competitive inhibitor of enzymes of the O-GlcNAcase family, shows 100000-fold selectivity over HexA/B, and binds to the O-GlcNAcase active site by mimicking the transition state as revealed by X-ray crystallography. This compound is able to raise O-GlcNAc levels in human HEK 293 and SH-SY5Y neuroblastoma cell lines and thus provides a novel, potent tool for the study of the role of O-GlcNAc in intracellular signal transduction pathways.

Glycosidase Inhibition: An Assessment of the Binding of 18 Putative Transition-State Mimics

The inhibition of glycoside hydrolases, through transition-state mimicry, is important both as a probe of enzyme mechanism and in the continuing quest for new drugs, notably in the treatment of cancer, HIV, influenza, and diabetes. The high affinity with which these enzymes are known to bind the transition state provides a framework upon which to design potent inhibitors. Recent work [for example, Bülow, A. et al. J. Am. Chem. Soc. 2000, 122, 8567-8568; Zechel, D. L. et al. J. Am. Chem. Soc. 2003, 125, 14313-14323] has revealed quite confusing and counter-intuitive patterns of inhibition for a number of glycosidase inhibitors. Here we describe a synergistic approach for analysis of inhibitors with a single enzyme 'model system', the Thermotoga maritima family 1 beta-glucosidase, TmGH1. The pH dependence of enzyme activity and inhibition has been determined, structures of inhibitor complexes have been solved by X-ray crystallography, with data up to 1.65 A resolution, and isothermal titration calorimetry was used to establish the thermodynamic signature. This has allowed the characterization of 18 compounds, all putative transition-state mimics, in order to build an 'inhibition profile' that provides an insight into what governs binding. In contrast to our preconceptions, there is little correlation of inhibitor chemistry with the calorimetric dissection of thermodynamics. The ensemble of inhibitors shows strong enthalpy-entropy compensation, and the random distribution of similar inhibitors across the plot of DeltaH degrees a vs TDeltaS degrees a likely reflects the enormous contribution of solvation and desolvation effects on ligand binding.