BIOCHEMISTRY: How Active Sites Communicate in Thiamine Enzymes

New Insights into the Cooperativity and Dynamics of Dimeric Enzymes

A survey of protein databases indicates that the majority of enzymes exist in oligomeric forms, with about half of those found in the UniProt database being homodimeric. Understanding why many enzymes are in their dimeric form is imperative. Recent developments in experimental and computational techniques have allowed for a deeper comprehension of the cooperative interactions between the subunits of dimeric enzymes. This review aims to succinctly summarize these recent advancements by providing an overview of experimental and theoretical methods, as well as an understanding of cooperativity in substrate binding and the molecular mechanisms of cooperative catalysis within homodimeric enzymes. Focus is set upon the beneficial effects of dimerization and cooperative catalysis. These advancements not only provide essential case studies and theoretical support for comprehending dimeric enzyme catalysis but also serve as a foundation for designing highly efficient catalysts, such as dimeric organic catalysts. Moreover, these developments have significant implications for drug design, as exemplified by Paxlovid, which was designed for the homodimeric main protease of SARS-CoV-2.

Simulations of Pathogenic E1α Variants: Allostery and Impact on Pyruvate Dehydrogenase Complex-E1 Structure and Function

Pyruvate dehydrogenase complex (PDC) deficiency is a major cause of primary lactic acidemia resulting in high morbidity and mortality, with limited therapeutic options. The E1 component of the mitochondrial multienzyme PDC (PDC-E1) is a symmetric dimer of heterodimers (αβ/α'β') encoded by the PDHA1 and PDHB genes, with two symmetric active sites each consisting of highly conserved phosphorylation loops A and B. PDHA1 mutations are responsible for 82-88% of cases. Greater than 85% of E1α residues with disease-causing missense mutations (DMMs) are solvent-inaccessible, with ∼30% among those involved in subunit-subunit interface contact (SSIC). We performed molecular dynamics simulations of wild-type (WT) PDC-E1 and E1 variants with E1α DMMs at R349 and W185 (residues involved in SSIC), to investigate their impact on human PDC-E1 structure. We evaluated the change in E1 structure and dynamics and examined their implications on E1 function with the specific DMMs. We found that the dynamics of phosphorylation Loop A, which is crucial for E1 biological activity, changes with DMMs that are at least about 15 Å away. Because communication is essential for PDC-E1 activity (with alternating active sites), we also investigated the possible communication network within WT PDC-E1 via centrality analysis. We observed that DMMs altered/disrupted the communication network of PDC-E1. Collectively, these results indicate allosteric effect in PDC-E1, with implications for the development of novel small-molecule therapeutics for specific recurrent E1α DMMs such as replacements of R349 responsible for ∼10% of PDC deficiency due to E1α DMMs.

The mechanism of a one-substrate transketolase reaction. Part II

In a recent paper, we showed the difference between the first stage of the one-substrate and the two-substrate transketolase reactions - the possibility of transfer of glycolaldehyde formed as a result of cleavage of the donor substrate from the thiazole ring of thiamine diphosphate to its aminopyrimidine ring through the tricycle formation stage, which is necessary for binding and splitting the second molecule of donor substrate [O.N. Solovjeva et al., The mechanism of a one-substrate transketolase reaction, Biosci. Rep. 40 (8) (2020) BSR20180246]. Here we show that under the action of the reducing agent a tricycle accumulates in a significant amount. Therefore, a significant decrease in the reaction rate of the one-substrate transketolase reaction compared to the two-substrate reaction is due to the stage of transferring the first glycolaldehyde molecule from the thiazole ring to the aminopyrimidine ring of thiamine diphosphate. Fragmentation of the four-carbon thiamine diphosphate derivatives showed that two glycolaldehyde molecules are bound to both coenzyme rings and the erythrulose molecule is bound to a thiazole ring. It was concluded that in the one-substrate reaction erythrulose is formed on the thiazole ring of thiamine diphosphate from two glycol aldehyde molecules linked to both thiamine diphosphate rings. The kinetic characteristics were determined for the two substrates, fructose 6-phosphate and glycolaldehyde.

The pyruvate dehydrogenase complexes: structure-based function and regulation

The pyruvate dehydrogenase complexes (PDCs) from all known living organisms comprise three principal catalytic components for their mission: E1 and E2 generate acetyl-coenzyme A, whereas the FAD/NAD(+)-dependent E3 performs redox recycling. Here we compare bacterial (Escherichia coli) and human PDCs, as they represent the two major classes of the superfamily of 2-oxo acid dehydrogenase complexes with different assembly of, and interactions among components. The human PDC is subject to inactivation at E1 by serine phosphorylation by four kinases, an inactivation reversed by the action of two phosphatases. Progress in our understanding of these complexes important in metabolism is reviewed.

Decarboxylation mechanisms in biological system

This review examines the mechanisms propelling cofactor-independent, organic cofactor-dependent and metal-dependent decarboxylase chemistry. Decarboxylation, the removal of carbon dioxide from organic acids, is a fundamentally important reaction in biology. Numerous decarboxylase enzymes serve as key components of aerobic and anaerobic carbohydrate metabolism and amino acid conversion. In the past decade, our knowledge of the mechanisms enabling these crucial decarboxylase reactions has continued to expand and inspire. This review focuses on the organic cofactors biotin, flavin, NAD, pyridoxal 5'-phosphate, pyruvoyl, and thiamin pyrophosphate as catalytic centers. Significant attention is also placed on the metal-dependent decarboxylase mechanisms.

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.

Communication between thiamin cofactors in the Escherichia coli pyruvate dehydrogenase complex E1 component active centers: evidence for a "direct pathway" between the 4'-aminopyrimidine N1' atoms

Kinetic, spectroscopic, and structural analysis tested the hypothesis that a chain of residues connecting the 4'-aminopyrimidine N1' atoms of thiamin diphosphates (ThDPs) in the two active centers of the Escherichia coli pyruvate dehydrogenase complex E1 component provides a signal transduction pathway. Substitution of the three acidic residues (Glu(571), Glu(235), and Glu(237)) and Arg(606) resulted in impaired binding of the second ThDP, once the first active center was filled, suggesting a pathway for communication between the two ThDPs. 1) Steady-state kinetic and fluorescence quenching studies revealed that upon E571A, E235A, E237A, and R606A substitutions, ThDP binding in the second active center was affected. 2) Analysis of the kinetics of thiazolium C2 hydrogen/deuterium exchange of enzyme-bound ThDP suggests half-of-the-sites reactivity for the E1 component, with fast (activated site) and slow exchanging sites (dormant site). The E235A and E571A variants gave no evidence for the slow exchanging site, indicating that only one of two active sites is filled with ThDP. 3) Titration of the E235A and E237A variants with methyl acetylphosphonate monitored by circular dichroism suggested that only half of the active sites were filled with a covalent predecarboxylation intermediate analog. 4) Crystal structures of E235A and E571A in complex with ThDP revealed the structural basis for the spectroscopic and kinetic observations and showed that either substitution affects cofactor binding, despite the fact that Glu(235) makes no direct contact with the cofactor. The role of the conserved Glu(571) residue in both catalysis and cofactor orientation is revealed by the combined results for the first time.

Multiple roles of mobile active center loops in the E1 component of the Escherichia coli pyruvate dehydrogenase complex - Linkage of protein dynamics to catalysis

The region encompassing residues 401-413 on the E1 component of the pyruvate dehydrogenase multienzyme complex from Escherichia coli comprises a loop (the inner loop) which was not seen in the X-ray structure in the presence of thiamin diphosphate, the required cofactor for the enzyme. This loop is seen in the presence of a stable analogue of the pre-decarboxylation intermediate, the covalent adduct between the substrate analogue methyl acetylphosphonate and thiamin diphosphate, C2α-phosphonolactylthiamin diphosphate. It has been shown that the residue H407 and several other residues on this loop are required to reduce the mobility of the loop so electron density corresponding to it can be seen once the pre-decarboxylation intermediate is formed. Concomitantly, the loop encompassing residues 541-557 (the outer loop) appears to work in tandem with the inner loop and there is a hydrogen bond between the two loops ensuring their correlated motion. The inner loop was shown to: a) sequester the active center from carboligase side reactions; b) assist the interaction between the E1 and the E2 components, thereby affecting the overall reaction rate of the entire multienzyme complex; c) control substrate access to the active center. Using viscosity effects on kinetics it was shown that formation of the pre-decarboxylation intermediate is specifically affected by loop movement. A cysteine-less variant was created for the E1 component, onto which cysteines were substituted at selected loop positions. Introducing an electron spin resonance spin label and an (19)F NMR label onto these engineered cysteines, the loop mobility was examined: a) both methods suggested that in the absence of ligand, the loop exists in two conformations; b) line-shape analysis of the NMR signal at different temperatures, enabled estimation of the rate constant for loop movement, and this rate constant was found to be of the same order of magnitude as the turnover number for the enzyme under the same conditions. Furthermore, this analysis gave important insights into rate-limiting thermal loop dynamics. Overall, the results suggest that the dynamic properties correlate with catalytic events on the E1 component of the pyruvate dehydrogenase complex.

Conformational ensemble modulates cooperativity in the rate-determining catalytic step in the E1 component of the Escherichia coli pyruvate dehydrogenase multienzyme complex

Cooperativity is extensively used by enzymes, particularly those acting at key metabolic branch points, to "fine tune" catalysis. Thus, cooperativity and enzyme catalysis are intimately linked, yet their linkage is poorly understood. Here we show that negative cooperativity in the rate-determining step in the E1 component of the Escherichia coli pyruvate dehydrogenase multienzyme complex is an outcome of redistribution of a "rate-promoting" conformational pre-equilibrium. An array of biophysical and biochemical studies indicates that non-catalytic but conserved residues directly regulate the redistribution. Furthermore, factors such as ligands and temperature, individually or in concert, also strongly influence the redistribution. As a consequence, these factors also exert their influence on catalysis by profoundly influencing the pre-equilibrium facilitated dynamics of communication between multienzyme components. Our observations suggest a mode of cooperativity in the E1 component that is consistent with the dynamical hypothesis shown to satisfactorily explain cooperativity in many well studied enzymes. The results point to the likely existence of multiple modes of communication between subunits when the entire class of thiamin diphosphate-dependent enzymes is considered.

Cofactor chemogenomics

Cofactors are organic molecules, most of them originating from vitamins, that bind to enzymes making them able to catalyze defined reactions. A cofactor-based chemogenomics approach exploits the presence of a cofactor-binding domain to develop compound scaffolds tailored to mimic the cofactor and to replace it within target enzyme classes. As a result, a loss of function is observed. An expansion of the cofactor scaffold to include structural/chemical features derived from the substrate, that usually binds at cofactor adjacent sites, increases the specificity of the enzyme fishing. This approach has been so far applied only to NAD(P)(+)-dependent enzymes. However, it is suitable for all other cofactors, with difficulties, for some of them, originated by very tight binding. In the case of cofactors covalently bound to the enzyme, the competition between the natural cofactor and the cofactor scaffold mimic can only occur during enzyme folding.

A multistate empirical valence bond description of protonatable amino acids

The multistate empirical valence bond (MS-EVB) model, which was developed for molecular dynamics simulations of proton transport in water and biomolecular systems, is extended for the modeling of protonatable amino acid residues in aqueous environments, specifically histidine and glutamic acid. The parameters of the MS-EVB force field are first determined to reproduce the geometries and energetics of the gas phase amino acid-water clusters. These parameters are then optimized to reproduce experimental pK(a) values. The free energy profiles for acid ionization and the corresponding pK(a) values are calculated by MS-EVB molecular dynamics simulations utilizing the umbrella sampling technique, with the center of excess charge coordinate chosen as the dissociation reaction coordinate. A general procedure for fitting the MS-EVB parameters is formulated, which allows for the parametrization of other amino acid residues with protonatable groups and the subsequent use of the MS-EVB approach for molecular dynamics simulations of proton transfer processes in proteins involving protonation/deprotonation of the protonatable amino acid groups.

Isotopic and other studies on the molecular origins of substrate regulation of some pyruvate decarboxylases: a reconsideration

Solvent isotope effect and beta-secondary isotope effect studies of the steady-state reaction of the pyruvate decarboxylase of Saccharomyces cerevisiae (ScPDC), which gains its full activity in a process requiring some seconds, suggested a model that ascribed the isotope effects to a reversible carbonyl-addition reaction between the 'regulatory pyruvate' molecule and the sulfhydryl group of Cys 221, occurring at the beginning and end of each catalytic cycle. The slow ('hysteretic') regulatory activation was thought to consist only of binding the regulatory molecule properly. Since the proposal of this model, much important progress in understanding both the structural and dynamic chemistry of ScPDC has been made in several laboratories. The purpose of this article is to review this new information and to evaluate the above model and others in the light of currently available evidence.

Inhibition of pyruvate decarboxylase from Z. mobilis by novel analogues of thiamine pyrophosphate: investigating pyrophosphate mimics

Replacement of the thiazolium ring of thiamine pyrophosphate with a triazole gives extremely potent inhibitors of pyruvate decarboxylase from Z. mobilis, with K(I) values down to 20 pM; this system was used to explore pyrophosphate mimics and several effective analogues were discovered.

Experimental observation of thiamin diphosphate-bound intermediates on enzymes and mechanistic information derived from these observations

Thiamin diphosphate (ThDP), the vitamin B1 coenzyme, is an excellent representative of coenzymes, which carry out electrophilic catalysis by forming a covalent complex with their substrates. The function of ThDP is to greatly increase the acidity of two carbon acids by stabilizing their conjugate bases, the ylide/C2-carbanion of the thiazolium ring and the C2alpha-carbanion (or enamine) once the substrate binds to ThDP. In recent years, several ThDP-bound intermediates on such pathways have been characterized by both solution and solid-state (X-ray) methods. Prominent among these advances are X-ray crystallographic results identifying both oxidative and non-oxidative intermediates, rapid chemical quench followed by NMR detection of a several intermediates which are stable under acidic conditions, and circular dichroism detection of the 1',4'-imino tautomer of ThDP in some of the intermediates. Some of these methods also enable the investigator to determine the rate-limiting step in the complex series of steps.