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

Biochemical, Bioinformatic, and Structural Comparisons of Transketolases and Position of Human Transketolase in the Enzyme Evolution

Transketolases (TKs) are key enzymes of the pentose phosphate pathway, regulating several other critical pathways in cells. Considering their metabolic importance, TKs are expected to be conserved throughout evolution. However, Tittmann et al. (J Biol Chem, 2010, 285(41): 31559-31570) demonstrated that Homo sapiens TK (hsTK) possesses several structural and kinetic differences compared to bacterial TKs. Here, we study 14 TKs from pathogenic bacteria, fungi, and parasites and compare them with hsTK using biochemical, bioinformatic, and structural approaches. For this purpose, six new TK structures are solved by X-ray crystallography, including the TK of Plasmodium falciparum. All of these TKs have the same general fold as bacterial TKs. This comparative study shows that hsTK greatly differs from TKs from pathogens in terms of enzymatic activity, spatial positions of the active site, and monomer-monomer interface residues. An ubiquitous structural pattern is identified in all TKs as a six-residue histidyl crown around the TK cofactor (thiamine pyrophosphate), except for hsTK containing only five residues in the crown. Residue mapping of the monomer-monomer interface and the active site reveals that hsTK contains more unique residues than other TKs. From an evolutionary standpoint, TKs from animals (including H. sapiens) and Schistosoma sp. belong to a distinct structural group from TKs of bacteria, plants, fungi, and parasites, mostly based on a different linker between domains, raising hypotheses regarding evolution and regulation.

In Vitro and In Silico Studies on Cytotoxic Properties of Oxythiamine and 2'-Methylthiamine

It is important to search for cytostatic compounds in order to fight cancer. One of them could be 2'-methylthiamine, which is a thiamine antimetabolite with an additional methyl group at the C-2 carbon of thiazole. So far, the cytostatic potential of 2'-methylthiamine has not been studied. We have come forward with a simplified method of synthesis using commercially available substrates and presented a comparison of its effects, as boosted by oxythiamine, on normal skin fibroblasts and HeLa cancer cells, having adopted in vitro culture techniques. Oxythiamine has been found to inhibit the growth and metabolism of cancer cells significantly better than 2'-methylthiamine (GI50 36 and 107 µM, respectively), while 2'-methylthiamine is more selective for cancer cells than oxythiamine (SI = 180 and 153, respectively). Docking analyses have revealed that 2'-methylthiamine (ΔG -8.2 kcal/mol) demonstrates a better affinity with thiamine pyrophosphokinase than thiamine (ΔG -7.5 kcal/mol ) and oxythiamine (ΔG -7.0 kcal/mol), which includes 2'-methylthiamine as a potential cytostatic. Our results suggest that the limited effect of 2'-methylthiamine on HeLa arises from the related arduous transport as compared to oxythiamine. Given that 2'-methylthiamine may possibly inhibit thiamine pyrophosphokinase, it could once again be considered a potential cytostatic. Thus, research should be carried out in order to find the best way to improve the transport of 2'-methylthiamine into cells, which may trigger its cytostatic properties.

Dynamic Protonation States Underlie Carbene Formation in ThDP-Dependent Enzymes: A Theoretical Study

The activation mechanism of thiamine diphosphate (ThDP) in enzymes has long been the subject of intense research and controversial discussion. Particularly contentious is the formation of a carbene intermediate, the first one observed in an enzyme. For the formation of the carbene to take place, both intramolecular and intermolecular proton transfer pathways have been proposed. However, the physiologically relevant pH of ThDP-dependent enzymes around neutrality does not seem to be suitable for the formation of such reactive chemical species. Herein, we investigate the general mechanism of activation of the ThDP cofactor in human transketolase (TKT), by means of electronic structure methods. We show that in the case of the human TKT, the carbene species is accessible through a pKa shift induced by the electrostatics of a neighboring histidine residue (H110), whose protonation state change modulates the pKa of ThDP and suppresses the latter by more than 6 pH units. Our findings highlight that ThDP enzymes activate the cofactor beyond simple geometric constraints and the canonical glutamate. Such observations in nature can pave the way for the design of biomimetic carbene catalysts and the engineering of tailored enzymatic carbenes.

New Role of Water in Transketolase Catalysis

Transketolase catalyzes the interconversion of keto and aldo sugars. Its coenzyme is thiamine diphosphate. The binding of keto sugar with thiamine diphosphate is possible only after C2 deprotonation of its thiazole ring. It is believed that deprotonation occurs due to the direct transfer of a proton to the amino group of its aminopyrimidine ring. Using mass spectrometry, it is shown that a water molecule is directly involved in the deprotonation process. After the binding of thiamine diphosphate with transketolase and its subsequent cleavage, a thiamine diphosphate molecule is formed with a mass increased by one oxygen molecule. After fragmentation, a thiamine diphosphate molecule is formed with a mass reduced by one and two hydrogen atoms, that is, HO and H₂O are split off. Based on these data, it is assumed that after the formation of holotransketolase, water is covalently bound to thiamine diphosphate, and carbanion is formed as a result of its elimination. This may be a common mechanism for other thiamine enzymes. The participation of a water molecule in the catalysis of the one-substrate transketolase reaction and a possible reason for the effect of the acceptor substrate on the affinity of the donor substrate for active sites are also shown.

QM/MM Study of Human Transketolase: Thiamine Diphosphate Activation Mechanism and Complete Catalytic Cycle

A computational model for human transketolase was proposed, showing that thiamine diphosphate activation was based on His110 in place of His481 reported in yeast transketolase. In addition, a complete catalytic reaction pathway was investigated using d-xylulose-5-phosphate and d-ribose-5-phosphate as substrates, showing at every step a perfect superimposition of our model with high-resolution crystallographic structures 3MOS, 4KXV, and 4KXX. This study shows that H2N4' of the active thiamine diphosphate "V form" no longer has a self-activating role but allows self-stabilization of the cofactor and of the Breslow intermediate. These advances in our knowledge of the human transketolase mechanism offer interesting prospects for the design of new drugs, this enzyme being involved in several diseases, and for a better understanding of the reactions catalyzed by transketolases from other sources.

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.

Constructing a synthetic pathway for acetyl-coenzyme A from one-carbon through enzyme design

Acetyl-CoA is a fundamental metabolite for all life on Earth, and is also a key starting point for the biosynthesis of a variety of industrial chemicals and natural products. Here we design and construct a Synthetic Acetyl-CoA (SACA) pathway by repurposing glycolaldehyde synthase and acetyl-phosphate synthase. First, we design and engineer glycolaldehyde synthase to improve catalytic activity more than 70-fold, to condense two molecules of formaldehyde into one glycolaldehyde. Second, we repurpose a phosphoketolase to convert glycolaldehyde into acetyl-phosphate. We demonstrated the feasibility of the SACA pathway in vitro, achieving a carbon yield ~50%, and confirmed the SACA pathway by ¹³C-labeled metabolites. Finally, the SACA pathway was verified by cell growth using glycolaldehyde, formaldehyde and methanol as supplemental carbon source. The SACA pathway is proved to be the shortest, ATP-independent, carbon-conserving and oxygen-insensitive pathway for acetyl-CoA biosynthesis, opening possibilities for producing acetyl-CoA-derived chemicals from one-carbon resources in the future.

The microbial synthesis of carbon-containing compounds from single carbon precursors is desirable, yet designed pathways to achieve this goal overlap with host metabolism. Here the authors design a de novo metabolic pathway to assimilate formaldehyde into acetyl-CoA that does not overlap with known metabolic networks.

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.