A QM/MM study on the catalytic mechanism of pyruvate decarboxylase

Improving tyrosol production efficiency through shortening the allosteric signal transmission distance of pyruvate decarboxylase

Tyrosol is an important chemical in medicine and chemical industries, which can be synthesized by a four-enzyme cascade pathway constructed in our previous study. However, the low catalytic efficiency of pyruvate decarboxylase from Candida tropicalis (CtPDC) in this cascade is a rate-limiting step. In this study, we resolved the crystal structure of CtPDC and investigated the mechanism of allosteric substrate activation and decarboxylation of this enzyme toward 4-hydroxyphenylpyruvate (4-HPP). In addition, based on the molecular mechanism and structural dynamic changes, we conducted protein engineering of CtPDC to improve decarboxylation efficiency. The conversion of the best mutant, CtPDC^{Q112G/Q162H/G415S/I417V} (CtPDC^Mu5), had over two-fold improvement compared to the wild-type. Molecular dynamic (MD) simulation revealed that the key catalytic distances and allosteric transmission pathways were shorter in CtPDC^Mu5 than in the wild type. Furthermore, when CtPDC in the tyrosol production cascade was replaced with CtPDC^Mu5, the tyrosol yield reached 38 g·L^-1 with 99.6% conversion and 1.58 g·L^-1·h^-1 space-time yield in 24 h through further optimization of the conditions. Our study demonstrates that protein engineering of the rate-limiting enzyme in the tyrosol synthesis cascade provides an industrial-scale platform for the biocatalytic production of tyrosol. KEY POINTS: • Protein engineering of CtPDC based on allosteric regulation improved the catalytic efficiency of decarboxylation. • The application of the optimum mutant of CtPDC removed the rate-limiting bottleneck in the cascade. • The final titer of tyrosol reached 38 g·L^-1 in 24 h in 3 L bioreactor.

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.

Molecular Dynamics and Free Energy Simulations of Phenylacetate and CO2 Release from AMDase and Its G74C/C188S Mutant: A Possible Rationale for the Reduced Activity of the Latter

Arylmalonate decarboxylase (AMDase) catalyzes the decarboxylation of α-aryl-α-methyl malonates to produce optically pure α-arylpropionates of industrial and medicinal importance. Herein, atomistic molecular dynamics simulations have been carried out to delineate the mechanism of the release of product molecules phenylacetate (PAC) and carbon dioxide (CO₂), from the wild-type (WT) and its G74C/C188S mutant enzymes. Both of the product molecules follow a crystallographically characterized solvent-accessible channel to come out of the protein interior. A higher free energy barrier for the release of PAC from G74C/C188S compared to that in the WT is consistent with the experimentally observed compromised efficiency of the mutant. The release of CO₂ precedes that of PAC; free energy barriers for CO₂ and PAC release in the WT enzyme are calculated to be ∼1-2 and ∼23 kcal/mol, respectively. Postdecarboxylation, CO₂ moves toward a hydrophobic pocket formed by Pro 14, Leu 38, Leu 40, Leu 77, and the side chain of Tyr 48 which serves as its temporary "reservoir". CO₂ releases following a channel mainly decorated by apolar residues, unlike in the case of oxalate decarboxylase where polar residues mediate its transport.

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.