Crystal structure of dihydroneopterin aldolase from Mycobacterium tuberculosis associated with 8-mercaptoguanine, and development of novel S8-functionalized analogues as inhibitors: Synthesis, enzyme inhibition, in vitro toxicity and antitubercular activity.

The crystallographic structure of the FolB enzyme from Mycobacterium tuberculosis (MtFolB), complexed with its inhibitor 8-mercaptoguanine (8-MG), was elucidated at a resolution of 1.95 Å. A novel series of S8-functionalized 8-MG derivatives were synthesised and evaluated as in vitro inhibitors of dihydroneopterin aldolase (DHNA, EC 4.1.2.25) activity of MtFolB. These compounds exhibited IC50 values in the submicromolar range. Evaluation of the activity for five compounds indicated their inhibition mode and inhibition constants. Molecular docking analyses were performed to determine the enzyme-inhibitor intermolecular interactions and ligand conformations upon complex formation. The inhibitory activities of all compounds against the M. tuberculosis H37Rv strain were evaluated. Compound 3e exhibited a minimum inhibitory concentration in the micromolar range. Finally, Compound 3e showed no apparent toxicity in both HepG2 and Vero cells. The findings presented herein will advance the quest for novel, specific inhibitors targeting MtFolB, an attractive molecular target for TB drug development.

Synthesis of Chiral Acyclic Pyrimidine Nucleoside Analogues from DHAP-Dependent Aldolases.

Dihydroxyacetone phosphate (DHAP)-dependent aldolases catalyze the aldol addition of DHAP to a variety of aldehydes and generate compounds with two stereocenters. This reaction is useful to synthesize chiral acyclic nucleosides, which constitute a well-known class of antiviral drugs currently used. In such compounds, the chirality of the aliphatic chain, which mimics the open pentose residue, is crucial for activity. In this work, three DHAP-dependent aldolases: fructose-1,6-biphosphate aldolase from rabbit muscle, rhanmulose-1-phosphate aldolase from Thermotoga maritima, and fuculose-1-phosphate aldolase from Escherichia coli, were used as biocatalysts. Aldehyde derivatives of thymine and cytosine were used as acceptor substrates, generating new acyclic nucleoside analogues containing two new stereocenters with conversion yields between 70% and 90%. Moreover, structural analyses by molecular docking were carried out to gain insights into the diasteromeric excess observed.

8-Mercaptoguanine-based inhibitors of Mycobacterium tuberculosis dihydroneopterin aldolase: synthesis, in vitro inhibition and docking studies.

The dihydroneopterin aldolase (DHNA, EC 4.1.2.25) activity of FolB protein is required for the conversion of 7,8-dihydroneopterin (DHNP) to 6-hydroxymethyl-7,8-dihydropterin (HP) and glycolaldehyde (GA) in the folate pathway. FolB protein from Mycobacterium tuberculosis (MtFolB) is essential for bacilli survival and represents an important molecular target for drug development. S8-functionalized 8-mercaptoguanine derivatives were synthesised and evaluated for inhibitory activity against MtFolB. The compounds showed IC50 values in the submicromolar range. The inhibition mode and inhibition constants were determined for compounds that exhibited the strongest inhibition. Additionally, molecular docking analyses were performed to suggest enzyme-inhibitor interactions and ligand conformations. To the best of our knowledge, this study describes the first class of MtFolB inhibitors.

Computational Investigation of Bisphosphate Inhibitors of 3-Deoxy-d-manno-octulosonate 8-phosphate Synthase.

The synthase, 3-deoxy-d-manno-octulosonate 8-phosphate (KDO8P), is a key enzyme for the lipopolysaccharide (LPS) biosynthesis of gram-negative bacteria and a potential target for developing new antimicrobial agents. In this study, computational molecular modeling methods were used to determine the complete structure of the KDO8P synthase from Neisseria meningitidis and to investigate the molecular mechanism of its inhibition by three bisphosphate inhibitors: BPH1, BPH2, and BPH3. Our results showed that BPH1 presented a protein-ligand complex with the highest affinity, which is in agreement with experimental data. Furthermore, molecular dynamics (MD) simulations showed that BPH1 is more active due to the many effective interactions, most of which are derived from its phosphoenolpyruvate moiety. Conversely, BPH2 exhibited few hydrogen interactions during the MD simulations with key residues located at the active sites of the KDO8P synthase. In addition, we hydroxylated BPH2 to create the hypothetical molecule named BPH3, to investigate the influence of the hydroxyl groups on the affinity of the bisphosphate inhibitors toward the KDO8P synthase. Overall, we discuss the main interactions between the KDO8P synthase and the bisphosphate inhibitors that are potential starting points for the design of new molecules with significant antibiotic activities.

Structure-Based Virtual Screening and Biochemical Evaluation for the Identification of Novel Trypanosoma brucei Aldolase Inhibitors.

INTRODUCTION: The glycolytic enzyme fructose-1,6-bisphosphate aldolase is a validated molecular target in human African trypanosomiasis (HAT) drug discovery, a neglected tropical disease (NTD) caused by the protozoan Trypanosoma brucei. Herein, a structure-based virtual screening (SBVS) approach to the identification of novel T. brucei aldolase inhibitors is described. Distinct molecular docking algorithms were used to screen more than 500,000 compounds against the X-ray structure of the enzyme. This SBVS strategy led to the selection of a series of molecules which were evaluated for their activity on recombinant T. brucei aldolase. The effort led to the discovery of structurally new ligands able to inhibit the catalytic activity of the enzyme. RESULTS: The predicted binding conformations were additionally investigated in molecular dynamics simulations, which provided useful insights into the enzyme-inhibitor intermolecular interactions. CONCLUSION: The molecular modeling results along with the enzyme inhibition data generated practical knowledge to be explored in further structure-based drug design efforts in HAT drug discovery.

Validation of Mycobacterium tuberculosis dihydroneopterin aldolase as a molecular target for anti-tuberculosis drug development.

An early step of target validation in antimicrobial drug discovery is to prove that a gene coding for a putative target is essential for pathogen's viability. However, little attention has been paid to demonstrate the causal links between gene essentiality and a particular protein function that will be the focus of a drug discovery effort. This should be considered an important step in target validation since a growing number of proteins are found to exhibit multiple and unrelated tasks. Here, we show that the Mycobacterium tuberculosis (Mtb) folB gene is essential and that this essentiality depends on the dihydroneopterin aldolase/epimerase activities of its protein product, the FolB protein from the folate biosynthesis pathway. The wild-type (WT) MtFolB and point mutants K99A and Y54F were cloned, expressed, purified and monitored for the aldolase, epimerase and oxygenase activities using HPLC. In contrast to the WT MtFolB, both mutants have neither aldolase nor epimerase activities in the conditions assayed. We then performed gene knockout experiments and showed that folB gene is essential for Mtb survival under the conditions tested. Moreover, only the WT folB sequence could be used as a rescue copy in gene complementation studies. When the sequences of mutants K99A or Y54F were used for complementation, no viable colonies were obtained, indicating that aldolase and/or epimerase activities are crucial for Mtb survival. These results provide a solid basis for further work aiming to develop new anti-TB agents acting as inhibitors of the aldolase/epimerase activities of MtFolB.

Expression of the 2-dehydro-3-deoxyphosphooc-tonate aldolase (KdsA) gene in mulberry leaves (Morus alba L.) is down-regulated under high salt and drought stress.

The oligosaccharide 3-deoxy-d-manno-oct-2-ulosonic acid (Kdo) is a key component of lipopolysaccharide in Gram-negative bacteria, and is also part of the pectic polysaccharide rhamnogalacturonan (RG-II) of the plant cell wall. The enzyme KDO-8-phosphate synthase (KDO8Ps), encoded by the 2-dehydro-3-deoxyphosphooctonate aldolase (KdsA) gene, catalyzes the first step in the synthesis of Kdo. In this study, the complete coding sequence of the KdsA gene from mulberry leaves was cloned and the primary structure of KDO8Ps was deduced. Alignment of the amino acid sequence of KDO8Ps from mulberry with those of five other plant species revealed a high level of evolutionary conservation. A phylogenetic tree analysis demonstrated a short genetic distance among KDO8Ps proteins of different species. Expression of the KdsA gene was higher in the second leaves than in the eighth leaves of mulberry, and was down-regulated under conditions of high salt or drought stress. Our results suggest that KdsA expression is important for the growth of new plant tissues, and is sensitive to harsh environments.

A full picture of enzymatic catalysis by hydroxynitrile lyases from Hevea brasiliensis: protonation dependent reaction steps and residue-gated movement of the substrate and the product.

Hydroxynitrile lyases (HNLs) defend plants from herbivores and microbial attack by releasing cyanide from hydroxynitriles. The reverse process has been productively applied to bioorganic syntheses of pharmaceuticals and agrochemicals. To improve our understanding of the catalytic mechanism of HNLs, extensive ab initio QM/MM and classical MM molecular dynamics simulations have been performed to explore the catalytic conversion of cyanohydrins into aldehyde (or ketone) and HCN by hydroxynitrile lyases from Hevea brasiliensis (HbHNLs). It was found that the catalytic reaction approximately follows a two-stage mechanism. The first stage involves two fast processes including the proton abstraction of the substrate through a double-proton transfer and the C-CN bond cleavage, while the second stage concerns HCN formation and is rate-determining. The complete free energy profile exhibits a peak of ∼18 kcal mol(-1). Interestingly, the protonation state of Lys236 influences the efficiency of the enzyme only to some extent, but it changes the entire catalytic mechanism. The dynamical behaviors of substrate delivery and HCN release are basically modulated by the gate movement of Trp128. The remarkable exothermicity of substrate binding and the facile release of HCN may drive the enzyme-catalyzed reaction to proceed along the substrate decomposition efficiently. Computational mutagenesis reveals the key residues which play an important role in the catalytic process.

Fructose 6-phosphate phosphoketolase activity in wild-type strains of Lactobacillus, isolated from the intestinal tract of pigs.

Phosphoketolases are key enzymes of the phosphoketolase pathway of heterofermentative lactic acid bacteria, which include lactobacilli. In heterofermentative lactobacilli xylulose 5-phosphate phosphoketolase (X5PPK) is the main enzyme of the phosphoketolase pathway. However, activity of fructose 6-phosphate phosphoketolase (F6PPK) has always been considered absent in lactic acid bacteria. In this study, the F6PPK activity was detected in 24 porcine wild-type strains of Lactobacillus reuteri and Lactobacillus mucosae, but not in the Lactobacillus salivarius or in L. reuteri ATCC strains. The activity of F6PPK increased after treatment of the culture at low-pH and diminished after porcine bile-salts stress conditions in wild-type strains of L. reuteri. Colorimetric quantification at 505 nm allowed to differentiate between microbial strains with low activity and without the activity of F6PPK. Additionally, activity of F6PPK and the X5PPK gene expression levels were evaluated by real time PCR, under stress and nonstress conditions, in 3 L. reuteri strains. Although an exact correlation, between enzyme activity and gene expression was not obtained, it remains possible that the xpk gene codes for a phosphoketolase with dual substrate, at least in the analyzed strains of L. reuteri.

The cyanogenic syndrome in rubber tree Hevea brasiliensis: tissue-damage-dependent activation of linamarase and hydroxynitrile lyase accelerates hydrogen cyanide release.

BACKGROUND AND AIMS: The release of hydrogen cyanide (HCN) from injured plant tissue affects multiple ecological interactions. Plant-derived HCN can act as a defence against herbivores and also plays an important role in plant-pathogen interactions. Crucial for activity as a feeding deterrent is the amount of HCN generated per unit time, referred to as cyanogenic capacity (HCNc). Strong intraspecific variation in HCNc has been observed among cyanogenic plants. This variation, in addition to genotypic variability (e.g. in Trifolium repens), can result from modifications in the expression level of the enzymes involved in either cyanogenic precursor formation or HCN release (as seen in Sorghum bicolor and Phaseolus lunatus). Thus, a modification or modulation of HCNc in reaction to the environment can only be achieved from one to the next generation when under genetic control and within days or hours when transcriptional regulations are involved. In the present study, it is shown that in rubber tree (Hevea brasiliensis) HCNc is modulated by post-translational activity regulation of the key enzymes for cyanide release. METHODS: Linamarase (LIN) and hydroxynitrile lyase (HNL) activity was determined by colorimetric assays utilizing dissociation of the substrates p-nitrophenyl-ß-d-glucopyranoside and acetone cyanohydrin, respectively. KEY RESULTS: In rubber tree leaves, LIN and HNL show up to ten-fold increased activity in response to tissue damage. This enzyme activation occurs within seconds and results in accelerated HCN formation. It is restricted to the damaged leaf area and depends on the severity of tissue damage. CONCLUSIONS: LIN and HNL activation (in contrast to genetic and transcriptional regulations) allows an immediate, local and damage type-dependent modulation of the cyanogenic response. Accordingly, this post-translational activation plays a decisive role in the defence of H. brasiliensis against herbivores as well as pathogens and may allow more flexible reactions in response to these different antagonists.

Selection of a new whole cell biocatalyst for the synthesis of 2-deoxyribose 5-phosphate.

2-deoxyribose 5-phosphate (DR5P) is a key intermediate in the biocatalyzed preparation of deoxyribonucleosides. Therefore, DR5P production by means of simpler, cleaner, and economic pathways becomes highly interesting. One strategy involves the use of bacterial whole cells containing DR5P aldolase as biocatalyst for the aldol addition between acetaldehyde and D: -glyceraldehyde 3-phosphate or glycolytic intermediates that in situ generate the acceptor substrate. In this work, diverse microorganisms capable of synthesizing DR5P were selected by screening several bacteria genera. In particular, Erwinia carotovora ATCC 33260 was identified as a new biocatalyst that afforded 14.1-mM DR5P starting from a cheap raw material like glucose.

Hydroxynitrile lyases from prunus seeds in the preparation of cyanohydrins.

The hydroxynitrile lyase (HNL) activity of nine defatted Prunus seeds was compared for catalyzing the addition of HCN to aromatic, heteroaromatic and α,ß-unsaturated aldehydes. Although the conversion and enantiomeric excess (ee) of the corresponding cyanohydrins were both influenced by the HNL source and the chemical structure of the aldehyde, Prunus HNLs were all suitable for the enantioselective preparation of cyanohydrins.

Catalytic mechanism of hydroxynitrile lyase from Hevea brasiliensis: a theoretical investigation.

Density functional theory (DFT) calculations using the hybrid functional B3LYP have been performed to investigate the catalytic mechanism of hydroxynitrile lyase from Hevea brasiliensis (Hb-HNL). This enzyme catalyzes the cleavage of acetone cyanohydrin to hydrocyanic acid plus acetone. Two models (A and B) of the active site consisting of 105 and 155 atoms, respectively, were constructed on the basis of the crystal structure. Good consistency between the two models provides a verification of the proposed mechanism. Our calculations show that the catalytic reaction proceeds via three elementary steps: (1) deprotonation of the OH-Ser80 by His235 and concomitant abstraction of a proton from the substrate hydroxyl by Ser80; (2) the C-C bond cleavage of the acetone cyanohydrin; and (3) protonation of the cleaved cyanide by His235. The cleavage of the C-C bond is the rate-limiting step with the overall free energy barrier of 13.5 kcal/mol for relatively smaller model A (14.9 kcal/mol for a larger model B) in the protein environment, which is in good agreement with experimental rate. The present results give support to the previously proposed general acid/base catalytic mechanism, in which the catalytic triad acts as a general acid/base. Moreover, the calculated results for model C, with the positive charge of Lys236 removed from model A, show that Lys236 with the positive charge plays a vital role in lowering the reaction barrier of the rate-determining and helps in stabilizing the negatively charged CN(-) by forming a hydrogen bond with the substrate, consistent with the experimental analysis.

alpha-Hydroxynitrile lyase protein from Xylella fastidiosa: Cloning, expression, and characterization.

Xylella fastidiosa is a xylem-restricted plant pathogen that causes a range of diseases in several and important crops. Through comparative genomic sequence analysis many genes were identified and, among them, several potentially involved in plant-pathogen interaction. The experimental determination of the primary sequence of some markedly expressed proteins for X. fastidiosa and the comparison with the nucleic acids sequence of genome identified one of them as being SCJ21.16 (XFa0032) gene product. The comparative analysis of this protein against SWISSPROT database, in special, resulted in similarity with alpha-hydroxynitrile lyase enzyme (HNL) from Arabidopsis thaliana, causing interest for being one of the most abundant proteins both in the whole cell extract as well as in the extracellular protein fraction. It is known that HNL enzyme are involved in a process termed "cyanogenesis", which catalyzes the dissociation of alpha-hydroxinitrile into carbonyle and HCN when plant tissue is damaged. Although the complete genome sequences of X. fastidiosa are available and the cyanogenesis process is well known, the biological role of this protein in this organism is not yet functionally characterized. In this study we presented the cloning, expression, characterization of recombinant HNL from X. fastidiosa, and its probable function in the cellular metabolism. The successful cloning and heterologous expression in Escherichia coli resulted in a satisfactory amount of the recombinant HNL expressed in a soluble, and active form giving convenient access to pure enzyme for biochemical and structural studies. Finally, our results confirmed that the product of the gene XFa0032 can be positively assigned as FAD-independent HNLs.

Observation of a short, strong hydrogen bond in the active site of hydroxynitrile lyase from Hevea brasiliensis explains a large pKa shift of the catalytic base induced by the reaction intermediate.

The hydroxynitrile lyase from Hevea brasiliensis (HbHNL) uses a catalytic triad consisting of Ser(80)-His(235)-Asp(207) to enhance the basicity of Ser(80)-O gamma for abstracting a proton from the OH group of the substrate cyanohydrin. Following the observation of a relatively short distance between a carboxyl oxygen of Asp(207) and the N delta(1)(His(235)) in a 1.1 A crystal structure of HbHNL, we here show by (1)H and (15)N-NMR spectroscopy that a short, strong hydrogen bond (SSHB) is formed between the two residues upon binding of the competitive inhibitor thiocyanate to HbHNL: the proton resonance of H-N delta 1(His(235)) moves from 15.41 ppm in the free enzyme to 19.35 ppm in the complex, the largest downfield shift observed so far upon inhibitor binding. Simultaneously, the D/H fractionation factor decreases from 0.98 to 0.35. In the observable pH range, i.e. between pH 4 and 10, no significant changes in chemical shifts (and therefore hydrogen bond strength) were observed for free HbHNL. For the complex with thiocyanate, the 19.35 ppm signal returned to 15.41 ppm at approximately pH 8, which indicates a pK(a) near this value for the H-N epsilon(2)(His(235)). These NMR results were analyzed on the basis of finite difference Poisson-Boltzmann calculations, which yielded the relative free energies of four protonation states of the His(235)-Asp(207) pair in solution as well as in the protein environment with and without bound inhibitor. The calculations explain all the NMR features, i.e. they suggest why a short, strong hydrogen bond is formed upon inhibitor binding and why this short, strong hydrogen bond reverts back to a normal one at approximately pH 8. Importantly, the computations also yield a shift of the free energy of the anionic state relative to the zwitterionic reference state by about 10.6 kcal/mol, equivalent to a shift in the apparent pK(a) of His(235) from 2.5 to 10. This huge inhibitor-induced increase in basicity is a prerequisite for His(235) to act as general base in the HbHNL-catalyzed cyanohydrin reaction.

Electrospray ionization mass spectrometry, circular dichroism and SAXS studies of the (S)-hydroxynitrile lyase from Hevea brasiliensis.

We report on experiments pertaining to solution properties of the (S)-hydroxynitrile lyase from Hevea brasiliensis (HbHNL). Small angle X-ray scattering unequivocally established the enzyme to occur in solution as a dimer, presumably of the same structure as in the crystal. The acid induced, irreversible deactivation of HbHNL was examined by electrospray ionization mass spectrometry (ESI-MS), circular dichroism (CD) and by measuring the enzyme activity. The deactivation is paralleled by an unfolding of the enzyme. ESI-MS of this 30000 Da per monomer heavy protein demonstrated that unfolding took place in several stages which are paralleled by a decrease in enzyme activity. Unfolding can also be observed by CD spectroscopy, and there is a clear correlation between enzyme activity and unfolding as detected by ESI-MS and CD.

Elucidation of the mode of substrate binding to hydroxynitrile lyase from Hevea brasiliensis.

The hydroxynitrile lyase from Hevea brasiliensis (Hb-HNL) is used as a catalyst in enantiospecific syntheses of alpha-hydroxynitriles from aldehydes and methyl-ketones. The catalyzed reaction represents one of the few industrially relevant examples of enzyme mediated C-C coupling reactions. In this work, we modeled Hb-HNL substrate complexes that have as yet proven inaccessible to experimental structure determination and were able to identify two binding modes for the natural substrate acetone cyanohydrin in docking simulations. Discrimination of the two alternatives was achieved by modeling complexes with two different chiral cyanohydrins followed by an analysis of the respective relative binding energies from molecular mechanics and thermodynamic integration. Only for one of the alternative binding modes the experimentally established S-selectivity of the enzyme was correctly predicted. Our results yielded further support for an enzymatic mechanism involving the catalytic triad Ser80, His235, and Asp207 as a general acid/base. A pivotal role was ascribed to Lys236, which seems to be crucial for enzymatic activity at low pH values. In addition, the modeling calculations provided possible explanations for the observed substrate and enantioselectivity of the enzyme that rationalize available mutational data and will be the basis for future protein engineering efforts.

Enzyme-based most probable number method for the enumeration of Bifidobacterium in dairy products.

An enzyme-based assay in combination with the most probable number (MPN) technique was developed for the enumeration of bifidobacteria. The assay employs the detection of fructose-6-phosphate phosphoketolase (F6PPK) activity as an indicator of the presence of bifidobacteria. The method was validated against viable counts and optimized with respect to selective media in order to quantitatively assess bifidobacteria in dairy products and other probiotic preparations. Several commercial products and homemade fermented milks were analyzed. Counts of bifidobacteria ranged from 10(7) to 10(8) MPN/ml in commercial products and homemade fermented milks. Commercial starters provided by Argentinean industries had between 10(7) and 10(11) MPN/ml. The results obtained in this study suggest that the combination of F6PPK activity determination and the MPN methodology allows an accurate determination of Bifidobacterium in pure cultures, dairy products, and other probiotic preparations.

An enzymatic-colorimetric assay for the quantification of Bifidobacterium.

An enzymatic-colorimetric assay for the quantification of Bifidobacterium was developed. The method, based upon the standard detection of fructose-6-phosphate phosphoketolase activity, was optimized with respect to bacterial cell pretreatment, time of incubation, and substrate concentration. The relationship between bacterial biomass and phosphoketolase activity was linear in a wide spectrum of bacterial densities. Higher sensitivity over the standard method was achieved by using 0.25% Triton X-100 in the reaction mixture to pretreat the bacterial cells. Because autoaggregation is a frequent feature among Bifidobacterium strains, this simple and reproducible method offers good advantage over viable plate count and turbidimetric techniques. The methodology can also be applied to the assessment of adherent Bifidobacterium strains to human epithelial cells.