A DFT study of solvation effects on the tautomeric equilibrium and catalytic ylide generation of thiamin models

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.

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.

Thiamine Deficiency and Neurodegeneration: the Interplay Among Oxidative Stress, Endoplasmic Reticulum Stress, and Autophagy

Thiamine (vitamin B1) is an essential nutrient and indispensable for normal growth and development of the organism due to its multilateral participation in key biochemical and physiological processes. Humans must obtain thiamine from their diet since it is synthesized only in bacteria, fungi, and plants. Thiamine deficiency (TD) can result from inadequate intake, increased requirement, excessive deletion, and chronic alcohol consumption. TD affects multiple organ systems, including the cardiovascular, muscular, gastrointestinal, and central and peripheral nervous systems. In the brain, TD causes a cascade of events including mild impairment of oxidative metabolism, neuroinflammation, and neurodegeneration, which are commonly observed in neurodegenerative diseases, such as Alzheimer's disease (AD), Parkinson's disease (PD), and Huntington's disease (HD). Thiamine metabolites may serve as promising biomarkers for neurodegenerative diseases, and thiamine supplementations exhibit therapeutic potential for patients of some neurodegenerative diseases. Experimental TD has been used to model aging-related neurodegenerative diseases. However, to date, the cellular and molecular mechanisms underlying TD-induced neurodegeneration are not clear. Recent research evidence indicates that TD causes oxidative stress, endoplasmic reticulum (ER) stress, and autophagy in the brain, which are known to contribute to the pathogenesis of various neurodegenerative diseases. In this review, we discuss the role of oxidative stress, ER stress, and autophagy in TD-mediated neurodegeneration. We propose that it is the interplay of oxidative stress, ER stress, and autophagy that contributes to TD-mediated neurodegeneration.

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.

Electron density reactivity indexes of the tautomeric/ionization forms of thiamin diphosphate

The generation of the highly reactive ylide in thiamin diphosphate catalysis is analyzed in terms of the nucleophilicity of key atoms, by means of density functional calculations at X3LYP/6-31++G(d,p) level of theory. The Fukui functions of all tautomeric/ionization forms are calculated in order to assess their reactivity. The results allow to conclude that the highly conserved glutamic residue does not protonate the N1' atom of the pyrimidyl ring, but it participates in a strong hydrogen bonding, stabilizing the eventual negative charge on the nitrogen, in all forms involved in the ylide generation. This condition provides the necessary reactivity on key atoms, N4' and C2, to carry out the formation of the ylide required to initiate the catalytic cycle of ThDP-dependent enzymes. This study represents a new approach for the ylide formation in ThDP catalysis.

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

Acetohydroxyacid synthase (AHAS) is a thiamin diphosphate dependent enzyme that catalyses the decarboxylation of pyruvate to yield the hydroxyethyl-thiamin diphosphate (ThDP) anion/enamine intermediate (HEThDP(-)). This intermediate reacts with a second ketoacid to form acetolactate or acetohydroxybutyrate as products. Whereas the mechanism involved in the formation of HEThDP(-) from pyruvate is well understood, the role of the enzyme in controlling the carboligation reaction of HEThDP(-) has not been determined yet. In this work, molecular dynamics (MD) simulations were employed to identify the aminoacids involved in the carboligation stage. These MD studies were carried out over the catalytic subunit of yeast AHAS containing the reaction intermediate (HEThDP(-)) and a second pyruvate molecule. Our results suggest that additional acid-base ionizable groups are not required to promote the catalytic cycle, in contrast with earlier proposals. This finding leads us to postulate that the formation of acetolactate relies on the acid-base properties of the HEThDP(-) intermediate itself. PM3 semiempirical calculations were employed to obtain the energy profile of the proposed mechanism on a reduced model of the active site. These calculations confirm the role of HEThDP(-) intermediate as the ionizable group that promotes the carboligation and product formation steps of the catalytic cycle.