Computational study on the carboligation reaction of acetohidroxyacid synthase: new approach on the role of the HEThDP- intermediate

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.

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.

Theoretical study of the catalytic mechanism of E1 subunit of pyruvate dehydrogenase multienzyme complex from Bacillus stearothermophilus

Pyruvate dehydrogenase multienzyme complex (PDHc) is a member of a family of 2-oxo acid dehydrogenase (OADH) multienzyme complexes involved in several central points of oxidative metabolism, and the E1 subunit is the most important component in the entire PDHc catalytic system, which catalyzes the reversible transfer of an acetyl group from a pyruvate to the lipoyl group of E2 subunit lipoly domain. In this article, the catalytic mechanism of the E1 subunit has been systematically studied using density functional theory (DFT). Four possible pathways with different general acid/base catalysts in decarboxylation and reductive acylation processes were explored. Our calculation results indicate that the 4'-amino pyrimidine of ThDP and residue His128 are the most likely proton donors in the decarboxylation and reductive acylation processes, respectively. During the reaction, each C-C and C-S bond formation or cleavage process, except for the liberation of CO2, is always accompanied by a proton transfer between the substrates and proton donors. The liberation of CO2 is calculated to be the rate-limiting step for the overall reaction, with an energy barrier of 13.57 kcal/mol. The decarboxylation process is endothermic by 5.32 kcal/mol, whereas the reductive acylation process is exothermic with a value of 5.74 kcal/mol. The assignment of protonation states of the surrounding residues can greatly influence the reaction. Residues His128 and His271 play roles in positioning the first substrate pyruvate and second substrate lipoyl group, respectively.

Expressed sequence tag analysis of khat (Catha edulis) provides a putative molecular biochemical basis for the biosynthesis of phenylpropylamino alkaloids

Khat (Catha edulis Forsk.) is a flowering perennial shrub cultivated for its neurostimulant properties resulting mainly from the occurrence of (S)-cathinone in young leaves. The biosynthesis of (S)-cathinone and the related phenylpropylamino alkaloids (1S,2S)-cathine and (1R,2S)-norephedrine is not well characterized in plants. We prepared a cDNA library from young khat leaves and sequenced 4,896 random clones, generating an expressed sequence tag (EST) library of 3,293 unigenes. Putative functions were assigned to > 98% of the ESTs, providing a key resource for gene discovery. Candidates potentially involved at various stages of phenylpropylamino alkaloid biosynthesis from L-phenylalanine to (1S,2S)-cathine were identified.