Density-functional study on the equilibria in the ThDP activation

Mechanistic study of the biosynthesis of R-phenylacetylcarbinol by acetohydroxyacid synthase enzyme using hybrid quantum mechanics/molecular mechanics simulations

The biosynthesis of R-phenylacetylcarbinol (R-PAC) by the acetohydroxy acid synthase, (AHAS) is addressed by molecular dynamics simulations (MD), hybrid quantum mechanics/molecular mechanics (QM/MM), and QM/MM free energy calculations. The results show the reaction starts with the nucleophilic attack of the C2α atom of the HEThDP intermediate on the C_β atom of the carbonyl group of benzaldehyde substrate via the formation of a transition state (TS1) with the HEThDP intermediate under 4'-aminopyrimidium (APH⁺) form. The calculated activation free energy for this step is 17.4 kcal mol^-1 at 27 °C. From this point, the reaction continues with the abstraction of H_β atom of the HEThDP intermediate by the O_β atom of benzaldehyde to form the intermediate I. The reaction is completed with the cleavage of the bond C2α-C2 to form the product R-PAC and to regenerate the ylide intermediate under the APH⁺ form, allowing in this way to reinitiate to the catalytic cycle once more. The calculated activation barrier for this last step is 15.9 kcal mol^-1 at 27 °C.

Mechanistic study of the biosynthesis of R-phenylcarbinol by acetohydroxyacid synthase enzyme using hybrid quantum mechanics/molecular mechanics simulations

The biosynthesis of R-phenylacetylcarbinol (R-PAC) by the acetohydroxy acid synthase, (AHAS) is addressed by molecular dynamics simulations (MD), hybrid quantum mechanics/molecular mechanics (QM/MM), and QM/MM free energy calculations. The results show the reaction starts with the nucleophilic attack of the C2α atom of the HEThDP intermediate on the C_β atom of the carbonyl group of benzaldehyde substrate via the formation of a transition state (TS1) with the HEThDP intermediate under 4'-aminopyrimidium (APH⁺) form. The calculated activation free energy for this step is 17.4kcal mol^-1 at 27 °C. From this point, the reaction continues with the abstraction of H_β atom of the HEThDP intermediate by the O_β atom of benzaldehyde to form the intermediate I. The reaction is completed with the cleavage of the bond C2α-C2 to form the product R-PAC and to regenerate the ylide intermediate under the APH⁺ form, allowing in this way to reinitiate to the catalytic cycle once more. The calculated activation barrier for this last step is 15.9kcal mol^-1 at 27 °C.

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.

Insights into the Thiamine Diphosphate Enzyme Activation Mechanism: Computational Model for Transketolase Using a Quantum Mechanical/Molecular Mechanical Method

We propose the first computational model for transketolase (TK), a thiamine diphosphate (ThDP)-dependent enzyme, using a quantum mechanical/molecular mechanical method on the basis of crystallographic TK structures from yeast and Escherichia coli, together with experimental kinetic data reported in the literature with wild-type and mutant TK. This model allowed us to define a new route for ThDP activation in the enzyme environment. We evidenced a strong interaction between ThDP and Glu418B of the TK active site, itself stabilized by Glu162A. The crucial point highlighted here is that deprotonation of ThDP C2 is not performed by ThDP N4' as reported in the literature, but by His481B, involving a HOH688A molecule bridge. Thus, ThDP N4' is converted from an amino form to an iminium form, ensuring the stabilization of the C2 carbanion or carbene. Finally, ThDP activation proceeds via an intermolecular process and not by an intramolecular one as reported in the literature. More generally, this proposed ThDP activation mechanism can be applied to some other ThDP-dependent enzymes and used to define the entire TK mechanism with donor and acceptor substrates more accurately.