Saturation mutagenesis of putative catalytic residues of benzoylformate decarboxylase provides a challenge to the accepted mechanism

Progress in research on the biosynthesis of 1,2,4-butanetriol by engineered microbes

1,2,4-butanetriol (BT) is a polyol with unique chemical properties, which has a stereocenter and can be divided into D-BT (the S-enantiomer) and L-BT (the R-enantiomer). BT can be used for the synthesis of 1,2,4-butanetriol trinitrate, 3-hydroxytetrahydrofuran, polyurethane, and other chemicals. It is widely used in the military industry, medicine, tobacco, polymer. At present, the BT is mainly synthesized by chemical methods, which are accompanied by harsh reaction conditions, poor selectivity, many by-products, and environmental pollution. Therefore, BT biosynthesis methods with the advantages of mild reaction conditions and green sustainability have become a current research hotspot. In this paper, the research status of microbial synthesis of BT was summarized from the following three aspects: (1) the biosynthetic pathway establishment for BT from xylose; (2) metabolic engineering strategies employed for improving BT production from xylose; (3) other substrates for BT production. Finally, the challenges and prospects of biosynthetic BT were discussed for future methods to improve competitiveness for industrial production.

From in-silico screening to in-vitro evaluation: Enhancing the detection of Microcystins with engineered PP1 mutant variants

Cyanotoxins produced during harmful algal blooms (CyanoHABs) have become a worldwide issue of concern. Microcystins (MC) are the most ubiquitous group of cyanotoxins and have known carcinogenic and hepatotoxic effects. The protein phosphatase inhibition assays (PPIAs), based on the inhibition of Protein Phosphatase 1/2A (PP1/PP2A) by MC, are one of the most cost-effective options for detecting MC. In this work, we aimed to design in-silico and evaluate in-vitro mutant variants of the PP1 protein, in order to enhance their capabilities as a MC biosensor. To this end, we performed an in-silico active site-saturated mutagenesis screening, followed by stability and docking affinity calculation with the MCLR cyanotoxin. Candidates with improved both affinity and stability were further tested in a fully flexible active-site docking. The best-scored mutations (19) were individually analysed regarding their locations and interactions. Four of them (p.D197F; p.Q249Y; p.S129W; p.D220Q) were selected for in-vitro expression and evaluation. Mutant p.D197F, exhibited a significant increment in inhibition by MCLR with respect to the WT, while showing a non-significant difference in stability nor activity. This successful PP1 inhibition enhancement suggests the potential of the p.D197F variant for practical MC detection applications.

Deep learning-based kcat prediction enables improved enzyme-constrained model reconstruction

Enzyme turnover numbers (kcat) are key to understanding cellular metabolism, proteome allocation and physiological diversity, but experimentally measured kcat data are sparse and noisy. Here we provide a deep learning approach (DLKcat) for high-throughput kcat prediction for metabolic enzymes from any organism merely from substrate structures and protein sequences. DLKcat can capture kcat changes for mutated enzymes and identify amino acid residues with a strong impact on kcat values. We applied this approach to predict genome-scale kcat values for more than 300 yeast species. Additionally, we designed a Bayesian pipeline to parameterize enzyme-constrained genome-scale metabolic models from predicted kcat values. The resulting models outperformed the corresponding original enzyme-constrained genome-scale metabolic models from previous pipelines in predicting phenotypes and proteomes, and enabled us to explain phenotypic differences. DLKcat and the enzyme-constrained genome-scale metabolic model construction pipeline are valuable tools to uncover global trends of enzyme kinetics and physiological diversity, and to further elucidate cellular metabolism on a large scale.

Bacterial mandelic acid degradation pathway and its application in biotechnology

Mandelic acid and its derivatives are an important class of chemical synthetic blocks, which is widely used in drug synthesis and stereochemistry research. In nature, mandelic acid degradation pathway has been widely identified and analyzed as a representative pathway of aromatic compounds degradation. The most studied mandelic acid degradation pathway from Pseudomonas putida consists of mandelate racemase, S-mandelate dehydrogenase, benzoylformate decarboxylase, benzaldehyde dehydrogenase and downstream benzoic acid degradation pathways. Because of the ability to catalyze various reactions of aromatic substrates, pathway enzymes have been widely used in biocatalysis, kinetic resolution, chiral compounds synthesis or construction of new metabolic pathways. In this paper, the physiological significance and the existing range of the mandelic acid degradation pathway were introduced first. Then each of the enzymes in the pathway is reviewed one by one, including the researches on enzymatic properties and the applications in biotechnology as well as efforts that have been made to modify the substrate specificity or improving catalytic activity by enzyme engineering to adapt different applications. The composition of the important metabolic pathway of bacterial mandelic acid degradation pathway as well as the researches and applications of pathway enzymes is summarized in this review for the first time.

Identification of the Gene Responsible for Lignin-Derived Low-Molecular-Weight Compound Catabolism in Pseudomonas sp. Strain LLC-1

Pseudomonas sp. strain LLC-1 (NBRC 111237) is capable of degrading lignin-derived low-molecular-weight compounds (LLCs). The genes responsible for the catabolism of LLCs were characterized in this study using whole-genome sequencing. Despite the close phylogenetic relationship with Pseudomonas putida, strain LLC-1 lacked the genes usually found in the P. putida genome, which included fer, encoding an enzyme for ferulic acid catabolism, and vdh encoding an NAD⁺-dependent aldehyde dehydrogenase specific for its catabolic intermediate, vanillin. Cloning and expression of the 8.5 kb locus adjacent to the van operon involved in vanillic acid catabolism revealed the bzf gene cluster, which is involved in benzoylformic acid catabolism. One of the structural genes identified, bzfC, expresses the enzyme (BzfC) having the ability to transform vanillin and syringaldehyde to corresponding acids, indicating that BzfC is a multifunctional enzyme that initiates oxidization of LLCs in strain LLC-1. Benzoylformic acid is a catabolic intermediate of (R,S)-mandelic acid in P. putida. Strain LLC-1 did not possess the genes for mandelic acid racemization and oxidation, suggesting that the function of benzoylformic acid catabolic enzymes is different from that in P. putida. Genome-wide characterization identified the bzf gene responsible for benzoylformate and vanillin catabolism in strain LLC-1, exhibiting a unique mode of dissimilation for biomass-derived aromatic compounds by this strain.

Characterization of Thiamine Diphosphate-Dependent 4-Hydroxybenzoylformate Decarboxylase Enzymes from Rhodococcus jostii RHA1 and Pseudomonas fluorescens Pf-5 Involved in Degradation of Aryl C2 Lignin Degradation Fragments

A thiamine diphosphate-dependent enzyme annotated as a benzoylformate decarboxylase is encoded by gene cluster ro02984-ro02986 in Rhodococcus jostii RHA1 previously shown to generate vanillin and 4-hydroxybenzaldehyde from lignin oxidation, and a closely related gene cluster is also found in the genome of Pseudomonas fluorescens Pf-5. Two hypotheses for possible pathways involving a thiamine diphosphate-dependent cleavage, either C-C cleavage of a ketol or diketone aryl C₃ substrate or decarboxylation of an aryl C₂ substrate, were investigated by expression and purification of the recombinant enzymes and expression of dehydrogenase and oxidase enzymes also found in the gene clusters. The ThDP-dependent enzymes showed no activity for cleavage of aryl C₃ ketol or diketone substrates but showed activity for decarboxylation of benzoylformate and 4-hydroxybenzoylformate. A flavin-dependent oxidase encoded by gene ro02984 was found to oxidize either mandelic acid or phenylglyoxal. The crystal structure of the P. fluorescens decarboxylase enzyme was determined at 1.69 Å resolution, showing similarity to structures of known benzoylformate decarboxylase enzymes. The P. fluorescens decarboxylase enzyme showed enhanced carboligase activity between vanillin and acetaldehyde, rationalized by the presence of alanine versus serine at residue 73 in the enzyme active site, which was investigated further by site-directed mutagenesis of this residue. A hypothesis for a pathway for degradation of aryl C₂ fragments arising from oxidative cleavage of phenylcoumaran and diarylpropane structures in lignin is proposed.

Cofactor Specificity Switch on Peach Glucitol Dehydrogenase

Most oxidoreductases that use NAD⁺ or NADP⁺ to transfer electrons in redox reactions display a strong preference for the cofactor. The catalytic efficiency of peach glucitol dehydrogenase (GolDHase) for NAD⁺ is 1800-fold higher than that for NADP⁺. Herein, we combined structural and kinetic data to reverse the cofactor specificity of this enzyme. Using site-saturation mutagenesis, we obtained the D216A mutant, which uses both NAD⁺ and NADP⁺, although with different catalytic efficiencies (1000 ± 200 and 170 ± 30 M^-1 s^-1, respectively). This mutant was used as a template to introduce further mutations by site-directed mutagenesis, using information from the fruit fly NADP-dependent GolDHase. The D216A/V217R/D218S triple mutant displayed a 2-fold higher catalytic efficiency with NADP⁺ than with NAD⁺. Overall, our results indicate that the triple mutant has the potential to be used for metabolic and cellular engineering and for cofactor recycling in industrial processes.

Metabolic Engineering of Escherichia coli for para-Amino-Phenylethanol and para-Amino-Phenylacetic Acid Biosynthesis

Aromatic amines are an important class of chemicals which are used as building blocks for the synthesis of polymers and pharmaceuticals. In this study we establish a de novo pathway for the biosynthesis of the aromatic amines para-amino-phenylethanol (PAPE) and para-amino-phenylacetic acid (4-APA) in Escherichia coli. We combined a synthetic para-amino-l-phenylalanine pathway with the fungal Ehrlich pathway. Therefore, we overexpressed the heterologous genes encoding 4-amino-4-deoxychorismate synthase (pabAB from Corynebacterium glutamicum), 4-amino-4-deoxychorismate mutase and 4-amino-4-deoxyprephenate dehydrogenase (papB and papC from Streptomyces venezuelae) and ThDP-dependent keto-acid decarboxylase (aro10 from Saccharomyces cerevisiae) in E. coli. The resulting para-amino-phenylacetaldehyde either was reduced to PAPE or oxidized to 4-APA. The wild type strain E. coli LJ110 with a plasmid carrying these four genes produced (in shake flask cultures) 11 ± 1.5 mg l⁻¹ of PAPE from glucose (4.5 g l⁻¹). By the additional cloning and expression of feaB (phenylacetaldehyde dehydrogenase from E. coli) 36 ± 5 mg l⁻¹ of 4-APA were obtained from 4.5 g l⁻¹ glucose. Competing reactions, such as the genes for aminotransferases (aspC and tyrB) or for biosynthesis of L-phenylalanine and L-tyrosine (pheA, tyrA) and for the regulator TyrR were removed. Additionally, the E. coli genes aroFBL were cloned and expressed from a second plasmid. The best producer strains of E. coli showed improved formation of PAPE and 4-APA, respectively. Plasmid-borne expression of an aldehyde reductase (yahK from E. coli) gave best values for PAPE production, whereas feaB-overexpression led to best values for 4-APA. In fed-batch cultivation, the best producer strains achieved 2.5 ± 0.15 g l⁻¹ of PAPE from glucose (11% C mol mol-1 glucose) and 3.4 ± 0.3 g l⁻¹ of 4-APA (17% C mol mol⁻¹ glucose), respectively which are the highest values for recombinant strains reported so far.

Exposing catalytic versatility of GTPases: taking reaction detours in mutants of hGBP1 enzyme without additional energetic cost

Our study reveals that the replacement of catalytically competent residues by the inert amino acid alanine, S73A and E99A, in hGBP1 opens a plethora of molecularly different reaction pathways featuring very similar energy barriers as the wild type.

A Theoretical Study of the Benzoylformate Decarboxylase Reaction Mechanism

Density functional theory calculations are used to investigate the detailed reaction mechanism of benzoylformate decarboxylase, a thiamin diphosphate (ThDP)-dependent enzyme that catalyzes the nonoxidative decarboxylation of benzoylformate yielding benzaldehyde and carbon dioxide. A large model of the active site is constructed on the basis of the X-ray structure, and it is used to characterize the involved intermediates and transition states and evaluate their energies. There is generally good agreement between the calculations and available experimental data. The roles of the various active site residues are discussed and the results are compared to mutagenesis experiments. Importantly, the calculations identify off-cycle intermediate species of the ThDP cofactor that can have implications on the kinetics of the reaction.

Progress in the experimental observation of thiamin diphosphate-bound intermediates on enzymes and mechanistic information derived from these observations

Thiamin diphosphate (ThDP), the vitamin B1 coenzyme is an excellent representative of coenzymes, which carry out electrophilic catalysis by forming a covalent complex with their substrates. The function of ThDP is to greatly increase the acidity of two carbon acids by stabilizing their conjugate bases, the ylide/carbene/C2-carbanion of the thiazolium ring and the C2α-carbanion/enamine, once the substrate binds to ThDP. In recent years, several ThDP-bound intermediates on such pathways have been characterized by both solution and solid-state methods. Prominent among these advances are X-ray crystallographic results identifying both oxidative and non-oxidative intermediates, rapid chemical quench followed by NMR detection of several intermediates which are stable under acidic conditions, solid-state NMR and circular dichroism detection of the states of ionization and tautomerization of the 4'-aminopyrimidine moiety of ThDP in some of the intermediates. These methods also enabled in some cases determination of the rate-limiting step in the complex series of steps. This review is an update of a review with the same title published by the authors in 2005 in this Journal. Much progress has been made in the intervening decade in the identification of the intermediates and their application to gain additional mechanistic insight.

Perturbation of the monomer-monomer interfaces of the benzoylformate decarboxylase tetramer

The X-ray structure of benzoylformate decarboxylase (BFDC) from Pseudomonas putida ATCC 12633 shows it to be a tetramer. This was believed to be typical of all thiamin diphosphate-dependent decarboxylases until recently when the structure of KdcA, a branched-chain 2-keto acid decarboxylase from Lactococcus lactis, showed it to be a homodimer. This lent credence to earlier unfolding experiments on pyruvate decarboxylase from Saccharomyces cerevisiae that indicated that it might be active as a dimer. To investigate this possibility in BFDC, we sought to shift the equilibrium toward dimer formation. Point mutations were made in the noncatalytic monomer–monomer interfaces, but these had a minimal effect on both tetramer formation and catalytic activity. Subsequently, the R141E/Y288A/A306F variant was shown by analytical ultracentrifugation to be partially dimeric. It was also found to be catalytically inactive. Further experiments revealed that just two mutations, R141E and A306F, were sufficient to markedly alter the dimer–tetramer equilibrium and to provide an ∼450-fold decrease in k_cat. Equilibrium denaturation studies suggested that the residual activity was possibly due to the presence of residual tetramer. The structures of the R141E and A306F variants, determined to <1.5 Å resolution, hinted that disruption of the monomer interfaces will be accompanied by movement of a loop containing Leu109 and Leu110. As these residues contribute to the hydrophobicity of the active site and the correct positioning of the substrate, it seems that tetramer formation may well be critical to the catalytic activity of BFDC.

Catalysis in Enzymatic Decarboxylations: Comparison of Selected Cofactor-dependent and Cofactor-independent Examples

This review is focused on three types of enzymes decarboxylating very different substrates: (1) Thiamin diphosphate (ThDP)-dependent enzymes reacting with 2-oxo acids; (2) Pyridoxal phosphate (PLP)-dependent enzymes reacting with α-amino acids; and (3) An enzyme with no known co-factors, orotidine 5'-monophosphate decarboxylase (OMPDC). While the first two classes have been much studied for many years, during the past decade studies of both classes have revealed novel mechanistic insight challenging accepted understanding. The enzyme OMPDC has posed a challenge to the enzymologist attempting to explain a 10¹⁷-fold rate acceleration in the absence of cofactors or even metal ions. A comparison of the available evidence on the three types of decarboxylases underlines some common features and more differences. The field of decarboxylases remains an interesting and challenging one for the mechanistic enzymologist notwithstanding the large amount of information already available.

Mechanistic insights into the bifunctional non-heme iron oxygenase carbapenem synthase by active site saturation mutagenesis

The carbapenem class of β-lactam antibiotics is known for its remarkable potency, antibacterial spectrum, and resistance to β-lactamase-mediated inactivation. While the biosynthesis of structurally "complex" carbapenems, such as thienamycin, share initial biochemical steps with carbapenem-3-carboxylate ("simple" carbapenem), the requisite inversion at C5 and formation of the characteristic α,β-unsaturated carboxylate are different in origin between the two groups. Here, we consider carbapenem synthase, a mechanistically distinct bifunctional non-heme iron α-ketoglutarate-dependent enzyme responsible for the terminal reactions, C5 epimerization and desaturation, in simple carbapenem production. Interestingly, this enzyme accepts two stereoisomeric substrates and transforms each to a common active antibiotic. Owing both to enzyme and product instability, resorting to saturation mutagenesis of active site and selected second-sphere residues gave clearly differing profiles of CarC tolerance to structural modification. Guided by a crystal structure and the mutational data, in silico docking was used to suggest the positioning of each disastereomeric substrate in the active site. The two orientations relative to the reactive iron-oxo center are manifest in the two distinct reactions, C5-epimerization and C2/3-desaturation. These observations favor a two-step reaction scheme involving two complete oxidative cycles as opposed to a single catalytic cycle in which an active site tyrosine, Tyr67, after hydrogen donation to achieve bicyclic ring inversion, is further hypothesized to serve as a radical carrier.

Catalytic effects of mutations of distant protein residues in human DNA polymerase β: theory and experiment

We carried out free-energy calculations and transient kinetic experiments for the insertion of the right (dC) and wrong (dA) nucleotides by wild-type (WT) and six mutant variants of human DNA polymerase β (Pol β). Since the mutated residues in the point mutants, I174S, I260Q, M282L, H285D, E288K, and K289M, were not located in the Pol β catalytic site, we assumed that the WT and its point mutants share the same dianionic phosphorane transition-state structure of the triphosphate moiety of deoxyribonucleotide 5'-triphosphate (dNTP) substrate. On the basis of this assumption, we have formulated a thermodynamic cycle for calculating relative dNTP insertion efficiencies, Ω = (k(pol)/K(D))(mut)/(k(pol)/K(D))(WT) using free-energy perturbation (FEP) and linear interaction energy (LIE) methods. Kinetic studies on five of the mutants have been published previously using different experimental conditions, e.g., primer-template sequences. We have performed a presteady kinetic analysis for the six mutants for comparison with wild-type Pol β using the same conditions, including the same primer/template DNA sequence proximal to the dNTP insertion site used for X-ray crystallographic studies. This consistent set of kinetic and structural data allowed us to eliminate the DNA sequence from the list of factors that can adversely affect calculated Ω values. The calculations using the FEP free energies scaled by 0.5 yielded 0.9 and 1.1 standard deviations from the experimental log Ω values for the insertion of the right and wrong dNTP, respectively. We examined a hybrid FEP/LIE method in which the FEP van der Waals term for the interaction of the mutated amino acid residue with its surrounding environment was replaced by the corresponding van der Waals term calculated using the LIE method, resulting in improved 0.4 and 1.0 standard deviations from the experimental log Ω values. These scaled FEP and FEP/LIE methods were also used to predict log Ω for R283A and R283L Pol β mutants.

Random mutagenesis MAPPIT analysis identifies binding sites for Vif and Gag in both cytidine deaminase domains of Apobec3G

The mammalian two-hybrid system MAPPIT allows the detection of protein-protein interactions in intact human cells. We developed a random mutagenesis screening strategy based on MAPPIT to detect mutations that disrupt the interaction of one protein with multiple protein interactors simultaneously. The strategy was used to detect residues of the human cytidine deaminase Apobec3G that are important for its homodimerization and its interaction with the HIV-1 Gag and Vif proteins. The strategy is able to identify the previously described head-to-head homodimerization interface in the N-terminal domain of Apobec3G. Our analysis further detects two new potential interaction surfaces in the N-and C-terminal domain of Apobec3G for interaction with Vif and Gag or for Apobec3G dimerization.

Assignment of function to histidines 260 and 298 by engineering the E1 component of the Escherichia coli 2-oxoglutarate dehydrogenase complex; substitutions that lead to acceptance of substrates lacking the 5-carboxyl group

The first component (E1o) of the Escherichia coli 2-oxoglutarate dehydrogenase complex (OGDHc) was engineered to accept substrates lacking the 5-carboxylate group by subjecting H260 and H298 to saturation mutagenesis. Apparently, H260 is required for substrate recognition, but H298 could be replaced with hydrophobic residues of similar molecular volume. To interrogate whether the second component would allow synthesis of acyl-coenzyme A derivatives, hybrid complexes consisting of recombinant components of OGDHc (o) and pyruvate dehydrogenase (p) enzymes were constructed, suggesting that a different component is the "gatekeeper" for specificity for these two multienzyme complexes in bacteria, E1p for pyruvate but E2o for 2-oxoglutarate.

In silico saturation mutagenesis and docking screening for the analysis of protein-ligand interaction: the Endothelial Protein C Receptor case study

Background

The design of mutants in protein functional regions, such as the ligand binding sites, is a powerful approach to recognize the determinants of specific protein activities in cellular pathways. For an exhaustive analysis of selected positions of protein structure large scale mutagenesis techniques are often employed, with laborious and time consuming experimental set-up. 'In silico' mutagenesis and screening simulation represents a valid alternative to laboratory methods to drive the 'in vivo' testing toward more focused objectives.

Results

We present here a high performance computational procedure for large-scale mutant modelling and subsequent evaluation of the effect on ligand binding affinity. The mutagenesis was performed with a 'saturation' approach, where all 20 natural amino acids were tested in positions involved in ligand binding sites. Each modelled mutant was subjected to molecular docking simulation and stability evaluation. The simulated protein-ligand complexes were screened for their impairment of binding ability based on change of calculated Ki compared to the wild-type.

An example of application to the Endothelial Protein C Receptor residues involved in lipid binding is reported.

Conclusion

The computational pipeline presented in this work is a useful tool for the design of structurally stable mutants with altered affinity for ligand binding, considerably reducing the number of mutants to be experimentally tested. The saturation mutagenesis procedure does not require previous knowledge of functional role of the residues involved and allows extensive exploration of all possible substitutions and their pairwise combinations. Mutants are screened by docking simulation and stability evaluation followed by a rationally driven selection of those presenting the required characteristics. The method can be employed in molecular recognition studies and as a preliminary approach to select models for experimental testing.

Understanding the mode of action of ThDP in benzoylformate decarboxylase

The mechanism of all elementary steps involved in the catalytic cycle of benzoylformate decarboxylase (BFD, E.C. 4.1.1.7) to generate the acyloin linkage is investigated by extensive molecular dynamics simulations. Models involving different charge states of amino acids and/or mutants of critical residues were constructed to understand the involvement of the catalytically active residues and the reactivity differences between different substrates in this reaction. Our calculations confirm that H70, S26, and H281 are catalytically active amino acids. H281 functions as a base to accept H(donor) in the first nucleophilic attack and as an acid in the second, to donate the proton back to O(acceptor). S26 assists H281 in deprotonation of the donor aldehyde and protonation of the acceptor aldehyde. In both the first and second nucleophilic attacks, H70 interacts with O(aldehyde) and aligns it toward the nucleophilic center. H70 has been found to have an electrostatic effect on the approaching aldehyde whose absence would block the initiation of the reaction. The reactivity difference between benzaldehyde (BA) and acetaldehyde (AA) is mainly explained by the steric interactions of the acceptor aldehyde with the surrounding amino acids in the active center of the enzyme.

Structural and kinetic studies on native intermediates and an intermediate analogue in benzoylformate decarboxylase reveal a least motion mechanism with an unprecedented short-lived predecarboxylation intermediate

The thiamin diphosphate- (ThDP-) dependent enzyme benzoylformate decarboxylase (BFDC) catalyzes the nonoxidative decarboxylation of benzoylformic acid to benzaldehyde and carbon dioxide. To date, no structural information for a cofactor-bound reaction intermediate in BFDC is available. For kinetic analysis, a chromophoric substrate analogue was employed that produces various absorbing intermediates during turnover but is a poor substrate with a 10(4)-fold compromised kcat. Here, we have analyzed the steady-state distribution of native intermediates by a combined chemical quench/1H NMR spectroscopic approach and estimated the net rate constants of elementary catalytic steps. At substrate saturation, carbonyl addition of the substrate to the cofactor (k' approximately 500 s-1 at 30 degrees C) and elimination of benzaldehyde (k' approximately 2.400 s-1) were found to be partially rate-determining for catalysis, whereas decarboxylation of the transient 2-mandelyl-ThDP intermediate is 1 order of magnitude faster with k' approximately 16.000 s-1, the largest rate constant of decarboxylation in any thiamin enzyme characterized so far. The X-ray structure of a predecarboxylation intermediate analogue was determined to 1.6 A after cocrystallization of BFDC from Pseudomonas putida with benzoylphosphonic acid methyl ester. In contrast to the free acid, for which irreversible phosphorylation of active center Ser26 was reported, the methyl ester forms a covalent adduct with ThDP with a similar configuration at C2alpha as observed for other thiamin enzymes. The C2-C2alpha bond of the intermediate analogue is out of plane by 7degrees, indicating strain. The phosphonate part of the adduct forms hydrogen bonds with Ser26 and His281, and the 1-OH group is held in place by interactions with His70 and the 4'-amino group of ThDP. The phenyl ring accommodates in a hydrophobic pocket formed by Phe464, Phe397, Leu109, and Leu403. A comparison with the previously determined structure of BFDC in noncovalent complex with the inhibitor (R)-mandelate suggests a least motion mechanism. Binding of benzoylphosphonic acid methyl ester to BFDC was further characterized by CD spectroscopy and stopped-flow kinetics, indicating a two-step binding mechanism with a 200-fold slower carbonyl addition to ThDP than determined for benzoylformic acid, in line with the observed slight structural reorganization of Phe464 due to steric clashes with the phosphonate moiety.

Snapshot of a reaction intermediate: analysis of benzoylformate decarboxylase in complex with a benzoylphosphonate inhibitor

Benzoylformate decarboxylase (BFDC) is a thiamin diphosphate- (ThDP-) dependent enzyme acting on aromatic substrates. In addition to its metabolic role in the mandelate pathway, BFDC shows broad substrate specificity coupled with tight stereo control in the carbon-carbon bond-forming reverse reaction, making it a useful biocatalyst for the production of chiral alpha-hydroxy ketones. The reaction of methyl benzoylphosphonate (MBP), an analogue of the natural substrate benzoylformate, with BFDC results in the formation of a stable analogue (C2alpha-phosphonomandelyl-ThDP) of the covalent ThDP-substrate adduct C2alpha-mandelyl-ThDP. Formation of the stable adduct is confirmed both by formation of a circular dichroism band characteristic of the 1',4'-iminopyrimidine tautomeric form of ThDP (commonly observed when ThDP forms tetrahedral complexes with its substrates) and by high-resolution mass spectrometry of the reaction mixture. In addition, the structure of BFDC with the MBP inhibitor was solved by X-ray crystallography to a spatial resolution of 1.37 A (PDB ID 3FSJ). The electron density clearly shows formation of a tetrahedral adduct between the C2 atom of ThDP and the carbonyl carbon atom of the MBP. This adduct resembles the intermediate from the penultimate step of the carboligation reaction between benzaldehyde and acetaldehyde. The combination of real-time kinetic information via stopped-flow circular dichroism with steady-state data from equilibrium circular dichroism measurements and X-ray crystallography reveals details of the first step of the reaction catalyzed by BFDC. The MBP-ThDP adduct on BFDC is compared to the recently solved structure of the same adduct on benzaldehyde lyase, another ThDP-dependent enzyme capable of catalyzing aldehyde condensation with high stereospecificity.

Reaction mechanisms of thiamin diphosphate enzymes: defining states of ionization and tautomerization of the cofactor at individual steps

We summarize the currently available information regarding the state of ionization and tautomerization of the 4'-aminopyrimidine ring of the thiamine diphosphate on enzymes requiring this coenzyme. This coenzyme forms a series of covalent intermediates with its substrates as an electrophilic catalyst, and the coenzyme itself also carries out intramolecular proton transfers, which is virtually unprecedented in coenzyme chemistry. An understanding of the state of ionization and tautomerization of the 4'-aminopyrimidine ring in each of these intermediates provides important details about proton movements during catalysis. CD spectroscopy, both steady-state and time-resolved, has proved crucial for obtaining this information because no other experimental method has provided such atomic detail so far.

Detection and time course of formation of major thiamin diphosphate-bound covalent intermediates derived from a chromophoric substrate analogue on benzoylformate decarboxylase

The mechanism of the enzyme benzoylformate decarboxylase (BFDC), which carries out a typical thiamin diphosphate (ThDP)-dependent nonoxidative decarboxylation reaction, was studied with the chromophoric alternate substrate (E)-2-oxo-4(pyridin-3-yl)-3-butenoic acid (3-PKB). Addition of 3-PKB resulted in the appearance of two transient intermediates formed consecutively, the first one to be formed a predecarboxylation ThDP-bound intermediate with lambda(max) at 477 nm, and the second one corresponding to the first postdecarboxylation intermediate the enamine with lambda(max) at 437 nm. The time course of formation/depletion of the PKB-ThDP covalent complex and of the enamine showed that decarboxylation was slower than formation of the PKB-ThDP covalent adduct. When the product of decarboxylation 3-(pyridin-3-yl)acrylaldehyde (PAA) was added to BFDC, again an absorbance with lambda(max) at 473 nm was formed, corresponding to the tetrahedral adduct of PAA with ThDP. Addition of well-formed crystals of BFDC to a solution of PAA resulted in a high resolution (1.34 A) structure of the BFDC-bound adduct of ThDP with PAA confirming the tetrahedral nature at the C2alpha atom, rather than of the enamine, and supporting the assignment of the lambda(max) at 473 nm to the PAA-ThDP adduct. The structure of the PAA-ThDP covalent complex is the first example of a product-ThDP adduct on BFDC. Similar studies with 3-PKB indicated that decarboxylation had taken place. Evidence was also obtained for the slow formation of the enamine intermediate when BFDC was incubated with benzaldehyde, the product of the decarboxylation reaction thus confirming its presence on the reaction pathway.

Experimental observation of thiamin diphosphate-bound intermediates on enzymes and mechanistic information derived from these observations

Thiamin diphosphate (ThDP), the vitamin B1 coenzyme, is an excellent representative of coenzymes, which carry out electrophilic catalysis by forming a covalent complex with their substrates. The function of ThDP is to greatly increase the acidity of two carbon acids by stabilizing their conjugate bases, the ylide/C2-carbanion of the thiazolium ring and the C2alpha-carbanion (or enamine) once the substrate binds to ThDP. In recent years, several ThDP-bound intermediates on such pathways have been characterized by both solution and solid-state (X-ray) methods. Prominent among these advances are X-ray crystallographic results identifying both oxidative and non-oxidative intermediates, rapid chemical quench followed by NMR detection of a several intermediates which are stable under acidic conditions, and circular dichroism detection of the 1',4'-imino tautomer of ThDP in some of the intermediates. Some of these methods also enable the investigator to determine the rate-limiting step in the complex series of steps.