A Theoretical Study of the Benzoylformate Decarboxylase Reaction Mechanism

Structural Analysis of Phosphonopyruvate Decarboxylase RhiEF: First Insights into an Ancestral Heterooligomeric Thiamine Pyrophosphate-Dependent Decarboxylase

The RhiE and RhiF proteins work together as RhiEF and function as a thiamine pyrophosphate (TPP)-dependent phosphonopyruvate decarboxylase to produce phosphonoacetaldehyde in the rhizocticin biosynthesis pathway. In this study, we determined the crystal structure of the RhiEF complexed with TPP and Mg2+. RhiEF forms a dimer of heterodimers, and the cofactor TPP is bound at the heterotetrameric subunit interface. Structural analysis of RhiEF revealed that the RhiE and RhiF moieties correspond to the pyrimidine-binding (PYR) and pyrophosphate-binding (PP) domains commonly found in TPP-dependent enzymes, respectively, as predicted by amino acid sequence alignment analysis. In contrast to other TPP-dependent enzymes with known structures, RhiEF has no domains other than the PYR and PP domains. Furthermore, structure-based evolutionary and sequence-based phylogenetic analyses have suggested that heteromultimeric enzymes such as RhiEF are ancestral types. These results indicate that RhiEF is one of the smallest and most ancient TPP-dependent decarboxylases. Based on the structural comparisons of RhiEF with other TPP-dependent decarboxylases, we identified the amino acid residues responsible for the catalytic mechanism of TPP-dependent decarboxylation in RhiEF.

Intramolecular Cyclization and a Retro-Ene Reaction Enable the Rapid Fragmentation of a Vitamin B1-Derived Breslow Intermediate

In solution, analogues of the Breslow intermediate formed during catalysis by benzoylformate decarboxylase (BFDC) undergo rapid, irreversible fragmentation. The ability of BFDC to prevent this reaction and preserve its cofactor is a striking example of an enzyme 'steering' a reactive intermediate towards a productive pathway. To understand how BFDC suppresses the off-pathway reactivity of this Breslow intermediate, a clear mechanistic understanding of the fragmentation reaction is required. Here, DFT calculations reveal an unexpected mechanism for the solution-phase fragmentation that involves an intramolecular cyclization and a subsequent retro-ene reaction to release the final products. Free energy profiles demonstrate that this pathway is significantly more facile than the previously proposed mechanism that invoked Breslow intermediate enolates as intermediates. Additional computations have been performed to understand why related Breslow intermediates do not undergo analogous fragmentation reactions. Calculations performed with two closely related Breslow intermediates suggest that subtle differences in the relative values of ΔG≠ for protonation and fragmentation dictate whether a given intermediate will fragment or not. These differences and the fragmentation mechanism unveiled in this work may have ramifications for the mechanism of BFDC and other thiamin-dependent enzymes and could provide general lessons related to the control of reactive intermediates by enzymes.

Eliminating Imaginary Vibrational Frequencies in Quantum-Chemical Cluster Models of Enzymatic Active Sites

In constructing finite models of enzyme active sites for quantum-chemical calculations, atoms at the periphery of the model must be constrained to prevent unphysical rearrangements during geometry relaxation. A simple fixed-atom or "coordinate-lock" approach is commonly employed but leads to undesirable artifacts in the form of small imaginary frequencies. These preclude evaluation of finite-temperature free-energy corrections, limiting thermochemical calculations to enthalpies only. Full-dimensional vibrational frequency calculations are possible by replacing the fixed-atom constraints with harmonic confining potentials. Here, we compare that approach to an alternative strategy in which fixed-atom contributions to the Hessian are simply omitted. While the latter strategy does eliminate imaginary frequencies, it tends to underestimate both the zero-point energy and the vibrational entropy while introducing artificial rigidity. Harmonic confining potentials eliminate imaginary frequencies and provide a flexible means to construct active-site models that can be used in unconstrained geometry relaxations, affording better convergence of reaction energies and barrier heights with respect to the model size, as compared to models with fixed-atom constraints.

The importance of thiamine (vitamin B1) in humans

Thiamine (thiamin, B1) is a vitamin necessary for proper cell function. It exists in a free form as a thiamine, or as a mono-, di- or triphosphate. Thiamine plays a special role in the body as a coenzyme necessary for the metabolism of carbohydrates, fats and proteins. In addition, it participates in the cellular respiration and oxidation of fatty acids: in malnourished people, high doses of glucose result in acute thiamine deficiency. It also participates in energy production in the mitochondria and protein synthesis. In addition, it is also needed to ensure the proper functioning of the central and peripheral nervous system, where it is involved in neurotransmitter synthesis. Its deficiency leads to mitochondrial dysfunction, lactate and pyruvate accumulation, and consequently to focal thalamic degeneration, manifested as Wernicke's encephalopathy or Wernicke-Korsakoff syndrome. It can also lead to severe or even fatal neurologic and cardiovascular complications, including heart failure, neuropathy leading to ataxia and paralysis, confusion, or delirium. The most common risk factor for thiamine deficiency is alcohol abuse. This paper presents current knowledge of the biological functions of thiamine, its antioxidant properties, and the effects of its deficiency in the body.

The Quantum Chemical Cluster Approach in Biocatalysis

ConspectusThe quantum chemical cluster approach has been used for modeling enzyme active sites and reaction mechanisms for more than two decades. In this methodology, a relatively small part of the enzyme around the active site is selected as a model, and quantum chemical methods, typically density functional theory, are used to calculate energies and other properties. The surrounding enzyme is modeled using implicit solvation and atom fixing techniques. Over the years, a large number of enzyme mechanisms have been solved using this method. The models have gradually become larger as a result of the faster computers, and new kinds of questions have been addressed. In this Account, we review how the cluster approach can be utilized in the field of biocatalysis. Examples from our recent work are chosen to illustrate various aspects of the methodology. The use of the cluster model to explore substrate binding is discussed first. It is emphasized that a comprehensive search is necessary in order to identify the lowest-energy binding mode(s). It is also argued that the best binding mode might not be the productive one, and the full reactions for a number of enzyme-substrate complexes have therefore to be considered to find the lowest-energy reaction pathway. Next, examples are given of how the cluster approach can help in the elucidation of detailed reaction mechanisms of biocatalytically interesting enzymes, and how this knowledge can be exploited to develop enzymes with new functions or to understand the reasons for lack of activity toward non-natural substrates. The enzymes discussed in this context are phenolic acid decarboxylase and metal-dependent decarboxylases from the amidohydrolase superfamily. Next, the application of the cluster approach in the investigation of enzymatic enantioselectivity is discussed. The reaction of strictosidine synthase is selected as a case study, where the cluster calculations could reproduce and rationalize the selectivities of both the natural and non-natural substrates. Finally, we discuss how the cluster approach can be used to guide the rational design of enzyme variants with improved activity and selectivity. Acyl transferase from Mycobacterium smegmatis serves as an instructive example here, for which the calculations could pinpoint the factors controlling the reaction specificity and enantioselectivity. The cases discussed in this Account highlight thus the value of the cluster approach as a tool in biocatalysis. It complements experiments and other computational techniques in this field and provides insights that can be used to understand existing enzymes and to develop new variants with tailored properties.

A Combined Computational–Experimental Study on the Substrate Binding and Reaction Mechanism of Salicylic Acid Decarboxylase

Salicylic acid decarboxylase (SDC) from the amidohydrolase superfamily (AHS) catalyzes the reversible decarboxylation of salicylic acid to form phenol. In this study, the substrate binding mode and reaction mechanism of SDC were investigated using computational and crystallographic methods. Quantum chemical calculations show that the enzyme follows the general mechanism of AHS decarboxylases. Namely, the reaction begins with proton transfer from a metal-coordinated aspartic acid residue (Asp298 in SDC) to the C1 of salicylic acid, which is followed by the C–C bond cleavage, to generate the phenol product and release CO2. Interestingly, the calculations show that SDC is a Mg-dependent enzyme rather than the previously proposed Zn-dependent, and the substrate is shown to be bidentately coordinated to the metal center in the catalysis, which is also different from the previous proposal. These predictions are corroborated by the crystal structure of SDC solved in complex with the substrate analogue 2-nitrophenol. The mechanistic insights into SDC in the present study provide important information for the rational design of the enzyme.

Reaction Mechanism of Human PAICS Elucidated by Quantum Chemical Calculations

Human PAICS is a bifunctional enzyme that is involved in the de novo purine biosynthesis, catalyzing the conversion of aminoimidazole ribonucleotide (AIR) into N-succinylcarboxamide-5-aminoimidazole ribonucleotide (SAICAR). It comprises two distinct active sites, AIR carboxylase (AIRc) where the AIR is initially converted to carboxyaminoimidazole ribonucleotide (CAIR) by reaction with CO₂ and SAICAR synthetase (SAICARs) in which CAIR then reacts with an aspartate to form SAICAR, in an ATP-dependent reaction. Human PAICS is a promising target for the treatment of various types of cancer, and it is therefore of high interest to develop a detailed understanding of its reaction mechanism. In the present work, density functional theory calculations are employed to investigate the PAICS reaction mechanism. Starting from the available crystal structures, two large models of the AIRc and SAICARs active sites are built and different mechanistic proposals for the carboxylation and phosphorylation–condensation mechanisms are examined. For the carboxylation reaction, it is demonstrated that it takes place in a two-step mechanism, involving a C–C bond formation followed by a deprotonation of the formed tetrahedral intermediate (known as isoCAIR) assisted by an active site histidine residue. For the phosphorylation–condensation reaction, it is shown that the phosphorylation of CAIR takes place before the condensation reaction with the aspartate. It is further demonstrated that the three active site magnesium ions are involved in binding the substrates and stabilizing the transition states and intermediates of the reaction. The calculated barriers are in good agreement with available experimental data.

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.

Using Atomic Confining Potentials for Geometry Optimization and Vibrational Frequency Calculations in Quantum-Chemical Models of Enzyme Active Sites

Quantum-chemical studies of enzymatic reaction mechanisms sometimes use truncated active-site models as simplified alternatives to mixed quantum mechanics molecular mechanics (QM/MM) procedures. Eliminating the MM degrees of freedom reduces the complexity of the sampling problem, but the trade-off is the need to introduce geometric constraints in order to prevent structural collapse of the model system during geometry optimizations that do not contain a full protein backbone. These constraints may impair the efficiency of the optimization, and care must be taken to avoid artifacts such as imaginary vibrational frequencies. We introduce a simple alternative in which terminal atoms of the model system are placed in soft harmonic confining potentials rather than being rigidly constrained. This modification is simple to implement and straightforward to use in vibrational frequency calculations, unlike iterative constraint-satisfaction algorithms, and allows the optimization to proceed without constraint even though the practical result is to fix the anchor atoms in space. The new approach is more efficient for optimizing minima and transition states, as compared to the use of fixed-atom constraints, and also more robust against unwanted imaginary frequencies. We illustrate the method by application to several enzymatic reaction pathways where entropy makes a significant contribution to the relevant reaction barriers. The use of confining potentials correctly describes reaction paths and facilitates calculation of both vibrational zero-point and finite-temperature entropic corrections to barrier heights.

The Catalytic Mechanism of Human Transketolase

We have computationally determined the catalytic mechanism of human transketolase (hTK) using a cluster model approach and density functional theory calculations. We were able to determine all the relevant structures, bringing solid evidences to the proposed experimental mechanism, and to add important detail to the structure of the transition states and the energy profile associated with catalysis. Furthermore, we have established the existence of a crucial intermediate of the catalytic cycle, in agreement with experiments. The calculated data brought new insights to hTK's catalytic mechanism, providing free-energy values for the chemical reaction, as well as adding atomistic detail to the experimental mechanism.

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.

Computational characterization of enzyme-bound thiamin diphosphate reveals a surprisingly stable tricyclic state: implications for catalysis

Thiamin diphosphate (ThDP)-dependent enzymes constitute a large class of enzymes that catalyze a diverse range of reactions. Many are involved in stereospecific carbon–carbon bond formation and, consequently, have found increasing interest and utility as chiral catalysts in various biocatalytic applications. All ThDP-catalyzed reactions require the reaction of the ThDP ylide (the activated state of the cofactor) with the substrate. Given that the cofactor can adopt up to seven states on an enzyme, identifying the factors affecting the stability of the pre-reactant states is important for the overall understanding of the kinetics and mechanism of the individual reactions.

In this paper we use density functional theory calculations to systematically study the different cofactor states in terms of energies and geometries. Benzoylformate decarboxylase (BFDC), which is a well characterized chiral catalyst, serves as the prototypical ThDP-dependent enzyme. A model of the active site was constructed on the basis of available crystal structures, and the cofactor states were characterized in the presence of three different ligands (crystallographic water, benzoylformate as substrate, and (R)-mandelate as inhibitor). Overall, the calculations reveal that the relative stabilities of the cofactor states are greatly affected by the presence and identity of the bound ligands. A surprising finding is that benzoylformate binding, while favoring ylide formation, provided even greater stabilization to a catalytically inactive tricyclic state. Conversely, the inhibitor binding greatly destabilized the ylide formation. Together, these observations have significant implications for the reaction kinetics of the ThDP-dependent enzymes, and, potentially, for the use of unnatural substrates in such reactions.