Direct Kinetic Evidence for Half-Of-The-Sites Reactivity in the E1 Component of the Human Pyruvate Dehydrogenase Multienzyme Complex through Alternating Sites Cofactor Activation†

Importance of Michaelis Constants for Cancer Cell Redox Balance and Lactate Secretion-Revisiting the Warburg Effect

Cancer cells metabolize a large fraction of glucose to lactate, even under a sufficient oxygen supply. This phenomenon-the "Warburg Effect"-is often regarded as not yet understood. Cancer cells change gene expression to increase the uptake and utilization of glucose for biosynthesis pathways and glycolysis, but they do not adequately up-regulate the tricarboxylic acid (TCA) cycle and oxidative phosphorylation (OXPHOS). Thereby, an increased glycolytic flux causes an increased production of cytosolic NADH. However, since the corresponding gene expression changes are not neatly fine-tuned in the cancer cells, cytosolic NAD+ must often be regenerated by loading excess electrons onto pyruvate and secreting the resulting lactate, even under sufficient oxygen supply. Interestingly, the Michaelis constants (KM values) of the enzymes at the pyruvate junction are sufficient to explain the priorities for pyruvate utilization in cancer cells: 1. mitochondrial OXPHOS for efficient ATP production, 2. electrons that exceed OXPHOS capacity need to be disposed of and secreted as lactate, and 3. biosynthesis reactions for cancer cell growth. In other words, a number of cytosolic electrons need to take the "emergency exit" from the cell by lactate secretion to maintain the cytosolic redox balance.

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.

Frontiers in the enzymology of thiamin diphosphate-dependent enzymes

Enzymes that use thiamin diphosphate (ThDP), the biologically active derivative of vitamin B1, as a cofactor play important roles in cellular metabolism in all domains of life. The analysis of ThDP enzymes in the past decades have provided a general framework for our understanding of enzyme catalysis of this protein family. In this review, we will discuss recent advances in the field that include the observation of "unusual" reactions and reaction intermediates that highlight the chemical versatility of the thiamin cofactor. Further topics cover the structural basis of cooperativity of ThDP enzymes, novel insights into the mechanism and structure of selected enzymes, and the discovery of "superassemblies" as reported, for example, acetohydroxy acid synthase. Finally, we summarize recent findings in the structural organisation and mode of action of 2-keto acid dehydrogenase multienzyme complexes and discuss future directions of this exciting research field.

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.

Structural and functional impact of clinically relevant E1α variants causing pyruvate dehydrogenase complex deficiency

Pyruvate dehydrogenase complex (PDC) catalyzes the oxidative decarboxylation of pyruvate to acetyl-coenzyme A, hinging glycolysis and the tricarboxylic acid cycle. PDC deficiency, an inborn error of metabolism, has a broad phenotypic spectrum. Symptoms range from fatal lactic acidosis or progressive neuromuscular impairment in the neonatal period, to chronic neurodegeneration. Most disease-causing mutations in PDC deficiency affect the PDHA1 gene, encoding the α subunit of the PDC-E1 component. Detailed biophysical analysis of pathogenic protein variants is a challenging approach to support the design of therapies based on improving and correcting protein structure and function. Herein, we report the characterization of clinically relevant PDC-E1α variants identified in Portuguese PDC deficient patients. These variants bear amino acid substitutions in different structural regions of PDC-E1α. The structural and functional analyses of recombinant heterotetrameric (αα'ββ') PDC-E1 variants, combined with molecular dynamics (MD) simulations, show a limited impact of the amino acid changes on the conformational stability, apart from the increased propensity for aggregation of the p.R253G variant as compared to wild-type PDC-E1. However, all variants presented a functional impairment in terms of lower residual PDC-E1 enzymatic activity and ≈3-100 × lower affinity for the thiamine pyrophosphate (TPP) cofactor, in comparison with wild-type PDC-E1. MD simulations neatly showed generally decreased stability (increased flexibility) of all variants with respect to the WT heterotetramer, particularly in the TPP binding region. These results are discussed in light of disease severity of the patients bearing such mutations and highlight the difficulty of developing chaperone-based therapies for PDC deficiency.

Assessing the Thiamine Diphosphate Dependent Pyruvate Dehydrogenase E1 Subunit for Carboligation Reactions with Aliphatic Ketoacids

The synthetic properties of the Thiamine diphosphate (ThDP)-dependent pyruvate dehydrogenase E1 subunit from Escherichia coli (EcPDH E1) was assessed for carboligation reactions with aliphatic ketoacids. Due to its role in metabolism, EcPDH E1 was previously characterised with respect to its biochemical properties, but it was never applied for synthetic purposes. Here, we show that EcPDH E1 is a promising biocatalyst for the production of chiral α-hydroxyketones. WT EcPDH E1 shows a 180-250-fold higher catalytic efficiency towards 2-oxobutyrate or pyruvate, respectively, in comparison to engineered transketolase variants from Geobacillus stearothermophilus (TK_GST). Its broad active site cleft allows for the efficient conversion of both (R)- and (S)-configured α-hydroxyaldehydes, next to linear and branched aliphatic aldehydes as acceptor substrates under kinetically controlled conditions. The alternate, thermodynamically controlled self-reaction of aliphatic aldehydes was shown to be limited to low levels of conversion, which we propose to be due to their large hydration constants. Additionally, the thermodynamically controlled approach was demonstrated to suffer from a loss of stereoselectivity, which makes it unfeasible for aliphatic substrates.

The Two-Species Model of transketolase explains donor substrate-binding, inhibition and heat-activation

We recently characterised a low-activity form of E. coli transketolase, TK_low, which also binds the cofactor thiamine pyrophosphate (TPP) with an affinity up to two-orders of magnitude lower than the previously known high TPP-affinity and high-activity form, TK_high, in the presence of Mg²⁺. We observed previously that partial oxidation was responsible for increased TK_high activity, while low-activity TK_low was unmodified. In the present study, the fluorescence-based cofactor-binding assay was adapted to detect binding of the β-hydroxypyruvate (HPA) donor substrate to wild-type transketolase and a variant, S385Y/D469T/R520Q, that is active towards aromatic aldehydes. Transketolase HPA affinity again revealed the two distinct forms of transketolase at a TK_high:TK_low ratio that matched those observed previously via TPP binding to each variant. The HPA dissociation constant of TK_low was comparable to the substrate-inhibition dissociation constant, K_i^HPA, determined previously. We provide evidence that K_i^HPA is a convolution of binding to the low-activity TK_low-TK_low dimer, and the TK_low subunit of the partially-active TK_high-TK_low mixed dimer, where HPA binding to the TK_low subunit of the mixed dimer results in inhibition of the active TK_high subunit. Heat-activation of transketolase was similarly investigated and found to convert the TK_low subunit of the mixed dimer to have TK_high-like properties, but without oxidation.

Hepatitis C Virus Downregulates Core Subunits of Oxidative Phosphorylation, Reminiscent of the Warburg Effect in Cancer Cells

Hepatitis C Virus (HCV) mainly infects liver hepatocytes and replicates its single-stranded plus strand RNA genome exclusively in the cytoplasm. Viral proteins and RNA interfere with the host cell immune response, allowing the virus to continue replication. Therefore, in about 70% of cases, the viral infection cannot be cleared by the immune system, but a chronic infection is established, often resulting in liver fibrosis, cirrhosis and hepatocellular carcinoma (HCC). Induction of cancer in the host cells can be regarded to provide further advantages for ongoing virus replication. One adaptation in cancer cells is the enhancement of cellular carbohydrate flux in glycolysis with a reduction of the activity of the citric acid cycle and aerobic oxidative phosphorylation. To this end, HCV downregulates the expression of mitochondrial oxidative phosphorylation complex core subunits quite early after infection. This so-called aerobic glycolysis is known as the “Warburg Effect” and serves to provide more anabolic metabolites upstream of the citric acid cycle, such as amino acids, pentoses and NADPH for cancer cell growth. In addition, HCV deregulates signaling pathways like those of TNF-β and MAPK by direct and indirect mechanisms, which can lead to fibrosis and HCC.

Novel insights into transketolase activation by cofactor binding identifies two native species subpopulations

Transketolase (TK) cofactor binding has been studied extensively over many years, yet certain mysteries remain, such as a lack of consensus on the cooperativity of thiamine pyrophosphate (TPP) binding into the two active sites, in the presence and absence of the divalent cation, Mg²⁺. Using a novel fluorescence-based assay, we determined directly the dissociation constants and cooperativity of TPP binding and provide the first comprehensive study over a broad range of cofactor concentrations. We confirmed the high-affinity dissociation constants and revealed a dependence of both the affinity and cooperativity of binding on [Mg²⁺], which explained the previous lack of consensus. A second, discrete and previously uncharacterised low-affinity TPP binding-site was also observed, and hence indicated the existence of two forms of TK with high- (TK_high) and low-affinity (TK_low). The relative proportions of each dimer were independent of the monomer-dimer transition, as probed by analytical ultracentrifugation at various [TK]. Mass spectrometry revealed that chemical oxidation of TK_low led to the formation of TK_high, which was 22-fold more active than TK_low. Finally, we propose a two-species model of transketolase activation that describes the interconversions between apo-/holo-TK_high and TK_low, and the potential to significantly improve biocatalytic activity by populating only the most active form.

Low-barrier hydrogen bonds in enzyme cooperativity

The underlying molecular mechanisms of cooperativity and allosteric regulation are well understood for many proteins, with haemoglobin and aspartate transcarbamoylase serving as prototypical examples^1,2. The binding of effectors typically causes a structural transition of the protein that is propagated through signalling pathways to remote sites and involves marked changes on the tertiary and sometimes even the quaternary level^1-5. However, the origin of these signals and the molecular mechanism of long-range signalling at an atomic level remain unclear^5-8. The different spatial scales and timescales in signalling pathways render experimental observation challenging; in particular, the positions and movement of mobile protons cannot be visualized by current methods of structural analysis. Here we report the experimental observation of fluctuating low-barrier hydrogen bonds as switching elements in cooperativity pathways of multimeric enzymes. We have observed these low-barrier hydrogen bonds in ultra-high-resolution X-ray crystallographic structures of two multimeric enzymes, and have validated their assignment using computational calculations. Catalytic events at the active sites switch between low-barrier hydrogen bonds and ordinary hydrogen bonds in a circuit that consists of acidic side chains and water molecules, transmitting a signal through the collective repositioning of protons by behaving as an atomistic Newton's cradle. The resulting communication synchronizes catalysis in the oligomer. Our studies provide several lines of evidence and a working model for not only the existence of low-barrier hydrogen bonds in proteins, but also a connection to enzyme cooperativity. This finding suggests new principles of drug and enzyme design, in which sequences of residues can be purposefully included to enable long-range communication and thus the regulation of engineered biomolecules.

Substrate-Based Allosteric Regulation of a Homodimeric Enzyme

Many enzymes operate through half-of-the sites reactivity wherein a single protomer is catalytically engaged at one time. In the case of the homodimeric enzyme, fluoroacetate dehalogenase, substrate binding triggers closing of a regulatory cap domain in the empty protomer, preventing substrate access to the remaining active site. However, the empty protomer serves a critical role by acquiring more disorder upon substrate binding, thereby entropically favoring the forward reaction. Empty protomer dynamics are also allosterically coupled to the bound protomer, driving conformational exchange at the active site and progress along the reaction coordinate. Here, we show that at high concentrations, a second substrate binds along the substrate-access channel of the occupied protomer, thereby dampening interprotomer dynamics and inhibiting catalysis. While a mutation (K152I) abrogates second site binding and removes inhibitory effects, it also precipitously lowers the maximum catalytic rate, implying a role for the allosteric pocket at low substrate concentrations, where only a single substrate engages the enzyme at one time. We show that this outer pocket first desolvates the substrate, whereupon it is deposited in the active site. Substrate binding to the active site then triggers the empty outer pocket to serve as an interprotomer allosteric conduit, enabling enhanced dynamics and sampling of activation states needed for catalysis. These allosteric networks and the ensuing changes resulting from second substrate binding are delineated using rigidity-based allosteric transmission theory and validated by nuclear magnetic resonance and functional studies. The results illustrate the role of dynamics along allosteric networks in facilitating function.

Structural and Functional Analyses of the Human PDH Complex Suggest a "Division-of-Labor" Mechanism by Local E1 and E3 Clusters

The pseudo-atomic structural model of human pyruvate dehydrogenase complex (PDHc) core composed of full-length E2 and E3BP components, calculated from our cryoelectron microscopy-derived density maps at 6-Å resolution, is similar to those of prokaryotic E2 structures. The spatial organization of human PDHc components as evidenced by negative-staining electron microscopy and native mass spectrometry is not homogeneous, and entails the unanticipated formation of local clusters of E1:E2 and E3BP:E3 complexes. Such uneven, clustered organization translates into specific duties for E1-E2 clusters (oxidative decarboxylation and acetyl transfer) and E3BP-E3 clusters (regeneration of reduced lipoamide) corresponding to half-reactions of the PDHc catalytic cycle. The addition of substrate coenzyme A modulates the conformational landscape of PDHc, in particular of the lipoyl domains, extending the postulated multiple random coupling mechanism. The conformational and associated chemical landscapes of PDHc are thus not determined entirely stochastically, but are restrained and channeled through an asymmetric architecture and further modulated by substrate binding.

Impairments in Oxidative Glucose Metabolism in Epilepsy and Metabolic Treatments Thereof

There is mounting evidence that oxidative glucose metabolism is impaired in epilepsy and recent work has further characterized the metabolic mechanisms involved. In healthy people eating a traditional diet, including carbohydrates, fats and protein, the major energy substrate in brain is glucose. Cytosolic glucose metabolism generates small amounts of energy, but oxidative glucose metabolism in the mitochondria generates most ATP, in addition to biosynthetic precursors in cells. Energy is crucial for the brain to signal "normally," while loss of energy can contribute to seizure generation by destabilizing membrane potentials and signaling in the chronic epileptic brain. Here we summarize the known biochemical mechanisms that contribute to the disturbance in oxidative glucose metabolism in epilepsy, including decreases in glucose transport, reduced activity of particular steps in the oxidative metabolism of glucose such as pyruvate dehydrogenase activity, and increased anaplerotic need. This knowledge justifies the use of alternative brain fuels as sources of energy, such as ketones, TCA cycle intermediates and precursors as well as even medium chain fatty acids and triheptanoin.

Part-of-the-sites binding and reactivity in the homooligomeric enzymes - facts and artifacts

For a number of enzymes composed of several subunits with the same amino acid sequence, it was documented, or suggested, that binding of a ligand, or catalysis, is carried out by a single subunit. This phenomenon may be the result of a pre-existent asymmetry of subunits or a limiting case of the negative cooperativity, and is sometimes called "half-of-the-sites binding (or reactivity)" for dimers and could be called "part-of-the-sites binding (or reactivity)" for higher oligomers. In this article, we discuss molecular mechanisms that may result in "part-of-the-sites binding (and reactivity)", offer possible explanations why it may have a beneficial role in enzyme function, and point to experimental problems in documenting this behaviour. We describe some cases, for which such a mechanism was first reported and later disproved. We also give several examples of enzymes, for which this mechanism seems to be well documented, and profitable. A majority of enzymes identified in this study as half-of-the-sites binding (or reactive) use it in the flip-flop version, in which "half-of-the-sites" refers to a particular moment in time. In general, the various variants of the mechanism seems to be employed often by oligomeric enzymes for allosteric regulation to enhance the efficiency of enzymatic reactions in many key metabolic pathways.

Crystal structures of archaeal 2-oxoacid:ferredoxin oxidoreductases from Sulfolobus tokodaii

As the first three-dimensional structure of the two-subunit type 2-oxoacid:ferredoxin oxidoreductases (OFOR) from archaea, we solved the crystal structures of STK_23000/STK_22980 (StOFOR1) and STK_24350/STK_24330 (StOFOR2) from Sulfolobus tokodaii. They showed similar overall structures, consisting of two a- and b-subunit heterodimers containing thiamin pyrophosphate (TPP) cofactor and [4Fe-4S] cluster, but lack an intramolecular ferredoxin domain. Unlike other OFORs, StOFORs can utilize both pyruvate and 2-oxoglutarate, playing a key role in the central metabolism. In the structure of StOFOR2 in unreacted pyruvate complex form, carboxylate group of pyruvate is recognized by Arg344 and Thr257 from the a-subunit, which are conserved in pyruvate:ferredoxin oxidoreductase from Desulfovbrio africanus (DaPFOR). In the structure of StOFOR1 co-crystallized with 2-oxobutyrate, electron density corresponding to a 1-hydroxypropyl group (post-decarboxylation state) was observed at the thiazole ring of TPP. The binding pockets of the StOFORs surrounding the methyl or propyl group of the ligands are wider than that of DaPFOR. Mutational analyses indicated that several residues were responsible for the broad 2-oxoacid specificity of StOFORs. We also constructed a possible complex structural model by placing a Zn(2+)-containing dicluster ferredoxin of S. tokodaii into the large pocket of StOFOR2, providing insight into the electron transfer between the two redox proteins.

Formation of reactive oxygen species by human and bacterial pyruvate and 2-oxoglutarate dehydrogenase multienzyme complexes reconstituted from recombinant components

Individual recombinant components of pyruvate and 2-oxoglutarate dehydrogenase multienzyme complexes (PDHc, OGDHc) of human and Escherichia coli (E. coli) origin were expressed and purified from E. coli with optimized protocols. The four multienzyme complexes were each reconstituted under optimal conditions at different stoichiometric ratios. Binding stoichiometries for the highest catalytic efficiency were determined from the rate of NADH generation by the complexes at physiological pH. Since some of these complexes were shown to possess 'moonlighting' activities under pathological conditions often accompanied by acidosis, activities were also determined at pH 6.3. As reactive oxygen species (ROS) generation by the E3 component of hOGDHc is a pathologically relevant feature, superoxide generation by the complexes with optimal stoichiometry was measured by the acetylated cytochrome c reduction method in both the forward and the reverse catalytic directions. Various known affectors of physiological activity and ROS production, including Ca(2+), ADP, lipoylation status or pH, were investigated. The human complexes were also reconstituted with the most prevalent human pathological mutant of the E3 component, G194C and characterized; isolated human E3 with the G194C substitution was previously reported to have an enhanced ROS generating capacity. It is demonstrated that: i. PDHc, similarly to OGDHc, is able to generate ROS and this feature is displayed by both the E. coli and human complexes, ii. Reconstituted hPDHc generates ROS at a significantly higher rate as compared to hOGDHc in both the forward and the reverse reactions when ROS generation is calculated for unit mass of their common E3 component, iii. The E1 component or E1-E2 subcomplex generates significant amount of ROS only in hOGDHc; iv. Incorporation of the G194C variant of hE3, the result of a disease-causing mutation, into reconstituted hOGDHc and hPDHc indeed leads to a decreased activity of both complexes and higher ROS generation by only hOGDHc and only in its reverse reaction.

Progress in the 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/carbene/C2-carbanion of the thiazolium ring and the C2α-carbanion/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 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 several intermediates which are stable under acidic conditions, solid-state NMR and circular dichroism detection of the states of ionization and tautomerization of the 4'-aminopyrimidine moiety of ThDP in some of the intermediates. These methods also enabled in some cases determination of the rate-limiting step in the complex series of steps. This review is an update of a review with the same title published by the authors in 2005 in this Journal. Much progress has been made in the intervening decade in the identification of the intermediates and their application to gain additional mechanistic insight.

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.

Catalysis in Enzymatic Decarboxylations: Comparison of Selected Cofactor-dependent and Cofactor-independent Examples

This review is focused on three types of enzymes decarboxylating very different substrates: (1) Thiamin diphosphate (ThDP)-dependent enzymes reacting with 2-oxo acids; (2) Pyridoxal phosphate (PLP)-dependent enzymes reacting with α-amino acids; and (3) An enzyme with no known co-factors, orotidine 5'-monophosphate decarboxylase (OMPDC). While the first two classes have been much studied for many years, during the past decade studies of both classes have revealed novel mechanistic insight challenging accepted understanding. The enzyme OMPDC has posed a challenge to the enzymologist attempting to explain a 10¹⁷-fold rate acceleration in the absence of cofactors or even metal ions. A comparison of the available evidence on the three types of decarboxylases underlines some common features and more differences. The field of decarboxylases remains an interesting and challenging one for the mechanistic enzymologist notwithstanding the large amount of information already available.

Determination of pre-steady-state rate constants on the Escherichia coli pyruvate dehydrogenase complex reveals that loop movement controls the rate-limiting step

Spectroscopic identification and characterization of covalent and noncovalent intermediates on large enzyme complexes is an exciting and challenging area of modern enzymology. The Escherichia coli pyruvate dehydrogenase multienzyme complex (PDHc), consisting of multiple copies of enzymic components and coenzymes, performs the oxidative decarboxylation of pyruvate to acetyl-CoA and is central to carbon metabolism linking glycolysis to the Krebs cycle. On the basis of earlier studies, we hypothesized that the dynamic regions of the E1p component, which undergo a disorder-order transition upon substrate binding to thiamin diphosphate (ThDP), play a critical role in modulation of the catalytic cycle of PDHc. To test our hypothesis, we kinetically characterized ThDP-bound covalent intermediates on the E1p component, and the lipoamide-bound covalent intermediate on the E2p component in PDHc and in its variants with disrupted active-site loops. Our results suggest that formation of the first covalent predecarboxylation intermediate, C2α-lactylthiamin diphosphate (LThDP), is rate limiting for the series of steps culminating in acetyl-CoA formation. Substitutions in the active center loops produced variants with up to 900-fold lower rates of formation of the LThDP, demonstrating that these perturbations directly affected covalent catalysis. This rate was rescued by up to 5-fold upon assembly to PDHc of the E401K variant. The E1p loop dynamics control covalent catalysis with ThDP and are modulated by PDHc assembly, presumably by selection of catalytically competent loop conformations. This mechanism could be a general feature of 2-oxoacid dehydrogenase complexes because such interfacial dynamic regions are highly conserved.

Determinants of brain cell metabolic phenotypes and energy substrate utilization unraveled with a modeling approach

Although all brain cells bear in principle a comparable potential in terms of energetics, in reality they exhibit different metabolic profiles. The specific biochemical characteristics explaining such disparities and their relative importance are largely unknown. Using a modeling approach, we show that modifying the kinetic parameters of pyruvate dehydrogenase and mitochondrial NADH shuttling within a realistic interval can yield a striking switch in lactate flux direction. In this context, cells having essentially an oxidative profile exhibit pronounced extracellular lactate uptake and consumption. However, they can be turned into cells with prominent aerobic glycolysis by selectively reducing the aforementioned parameters. In the case of primarily oxidative cells, we also examined the role of glycolysis and lactate transport in providing pyruvate to mitochondria in order to sustain oxidative phosphorylation. The results show that changes in lactate transport capacity and extracellular lactate concentration within the range described experimentally can sustain enhanced oxidative metabolism upon activation. Such a demonstration provides key elements to understand why certain brain cell types constitutively adopt a particular metabolic profile and how specific features can be altered under different physiological and pathological conditions in order to face evolving energy demands.

In an environment with appropriate oxygen levels (normoxia), most eukaryotic cells produce energy by oxidizing glucose into carbon dioxide and water. In this process, glucose is transformed into pyruvate, which then fuels oxidative phosphorylation in the mitochondria. Interestingly, Otto Warburg reported back in the 1920's that some eukaryotic cells prominently process glucose-derived pyruvate into lactate, hence “avoiding" the mitochondrial oxidation despite adequate oxygen concentrations. This phenomenon was termed aerobic glycolysis and was first observed in cancer cells. Since then, it has also been described in several normal tissues including the central nervous system. The biochemical basis of aerobic glycolysis has remained elusive until now. Taking advantage of a modeling approach, we unraveled the main metabolic characteristics that determine whether a cell will be strictly oxidative or rather will exhibit aerobic glycolysis. When applied in the context of the central nervous system, our findings not only provide a theoretical demonstration of why neurons and astrocytes differ in terms of metabolic profile, but also suggest that such complementarity forms the basis for metabolic cooperation between the two cell types.

Behavior of the ATP grasp domain of biotin carboxylase monomers and dimers studied using molecular dynamics simulations

The enzyme biotin carboxylase (BC) uses adenosine triphosphate (ATP) to carboxylate biotin and is involved in fatty acid synthesis. Structural evidence suggests that the B domain of BC undergoes a large hinge motion of ∼45° when binding and releasing substrates. Escherichia coli BC can function as a natural homodimer and as a mutant monomer. Using molecular dynamics simulations, we evaluate the free energy profile along a closure angle of the B domain of E. coli BC for three cases: a monomer without bound Mg(2)ATP, a monomer with bound Mg(2)ATP, and a homodimer with bound Mg(2)ATP in one subunit. The simulation results show that a closed state is the most probable for the monomer with or without bound Mg(2)ATP. For the dimer with Mg(2)ATP in one of its subunits, communication between the two subunits was observed. Specifically, in the dimer, the opening of the subunit without Mg(2)ATP caused the other subunit to open, and hysteresis was observed upon reclosing it. The most stable state of the dimer is one in which the B domain of both subunits is closed; however, the open state for the B domain without Mg(2)ATP is only approximately 2k(B)T higher in free energy than the closed state. A simple diffusion model indicates that the mean times for opening and closing of the B domain in the monomer with and without Mg(2)ATP are much smaller than the overall reaction time, which is on the order of seconds.

Quantum-chemical calculations of a long proton wire. Application of a harmonic model to analysis of the structure of an ionic defect in a water chain with an excess proton

Quantum-chemical calculations of molecular complexes (NH(3))(3)Zn(2+)...(H(2)O)(n)...NH(3) (C(n), n = 11, 16, 21, and 30) simulating a proton wire donor-water chain-acceptor were carried out. Earlier found periodicity in the length of the O-H bonds in water chain is explained within the framework of a one-component harmonic model. In complexes C(n), the geometry and electronic structure of ionic defect in water chain with an excess proton were studied. Calculations carried out at ab initio (B3LYP/6-31+G**) and semiempirical (PM3) levels of theory predict different patterns of distribution of the O-H bonds lengths and positive charge on the H-bond hydrogen atoms in the region of ionic defect. The obtained data show how a length of water chain and position of a protonated water link in the chain influence the ionic defect structure. To describe the observed structures of ionic defect, the harmonic model was used and the role of parameters of the H-bonded chain was investigated. The performed analysis explains different mechanisms (concerted and stepwise) of proton transfer along the H-bonded chain derived from ab initio and semiempirical calculation schemes.

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.

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.

Snapshot of a reaction intermediate: analysis of benzoylformate decarboxylase in complex with a benzoylphosphonate inhibitor

Benzoylformate decarboxylase (BFDC) is a thiamin diphosphate- (ThDP-) dependent enzyme acting on aromatic substrates. In addition to its metabolic role in the mandelate pathway, BFDC shows broad substrate specificity coupled with tight stereo control in the carbon-carbon bond-forming reverse reaction, making it a useful biocatalyst for the production of chiral alpha-hydroxy ketones. The reaction of methyl benzoylphosphonate (MBP), an analogue of the natural substrate benzoylformate, with BFDC results in the formation of a stable analogue (C2alpha-phosphonomandelyl-ThDP) of the covalent ThDP-substrate adduct C2alpha-mandelyl-ThDP. Formation of the stable adduct is confirmed both by formation of a circular dichroism band characteristic of the 1',4'-iminopyrimidine tautomeric form of ThDP (commonly observed when ThDP forms tetrahedral complexes with its substrates) and by high-resolution mass spectrometry of the reaction mixture. In addition, the structure of BFDC with the MBP inhibitor was solved by X-ray crystallography to a spatial resolution of 1.37 A (PDB ID 3FSJ). The electron density clearly shows formation of a tetrahedral adduct between the C2 atom of ThDP and the carbonyl carbon atom of the MBP. This adduct resembles the intermediate from the penultimate step of the carboligation reaction between benzaldehyde and acetaldehyde. The combination of real-time kinetic information via stopped-flow circular dichroism with steady-state data from equilibrium circular dichroism measurements and X-ray crystallography reveals details of the first step of the reaction catalyzed by BFDC. The MBP-ThDP adduct on BFDC is compared to the recently solved structure of the same adduct on benzaldehyde lyase, another ThDP-dependent enzyme capable of catalyzing aldehyde condensation with high stereospecificity.

Reaction mechanisms of thiamin diphosphate enzymes: defining states of ionization and tautomerization of the cofactor at individual steps

We summarize the currently available information regarding the state of ionization and tautomerization of the 4'-aminopyrimidine ring of the thiamine diphosphate on enzymes requiring this coenzyme. This coenzyme forms a series of covalent intermediates with its substrates as an electrophilic catalyst, and the coenzyme itself also carries out intramolecular proton transfers, which is virtually unprecedented in coenzyme chemistry. An understanding of the state of ionization and tautomerization of the 4'-aminopyrimidine ring in each of these intermediates provides important details about proton movements during catalysis. CD spectroscopy, both steady-state and time-resolved, has proved crucial for obtaining this information because no other experimental method has provided such atomic detail so far.

Structural basis for inactivation of the human pyruvate dehydrogenase complex by phosphorylation: role of disordered phosphorylation loops

We report the crystal structures of the phosporylated pyruvate dehydrogenase (E1p) component of the human pyruvate dehydrogenase complex (PDC). The complete phosphorylation at Ser264-alpha (site 1) of a variant E1p protein was achieved using robust pyruvate dehydrogenase kinase 4 free of the PDC core. We show that unlike its unmodified counterpart, the presence of a phosphoryl group at Ser264-alpha prevents the cofactor thiamine diphosphate-induced ordering of the two loops carrying the three phosphorylation sites. The disordering of these phosphorylation loops is caused by a previously unrecognized steric clash between the phosphoryl group at site 1 and a nearby Ser266-alpha, which nullifies a hydrogen-bonding network essential for maintaining the loop conformations. The disordered phosphorylation loops impede the binding of lipoyl domains of the PDC core to E1p, negating the reductive acetylation step. This results in the disruption of the substrate channeling in the PDC, leading to the inactivation of this catalytic machine.

Snapshots of catalysis in the E1 subunit of the pyruvate dehydrogenase multienzyme complex

The pyruvate dehydrogenase multienzyme assembly (PDH) generates acetyl coenzyme A and reducing equivalents from pyruvate in a multiple-step process that is a nexus of central metabolism. We report crystal structures of the Geobacillus stearothermophilus PDH E1p subunit with ligands that mimic the prereaction complex and the postdecarboxylation product. The structures implicate residues that help to orient substrates, nurture intermediates, and organize surface loops so that they can engage a mobile lipoyl domain that receives the acetyl group and shuttles it to the next active site. The structural and enzymatic data suggest that H128beta performs a dual role: first, as electrostatic catalyst of the reaction of pyruvate with the thiamine cofactor; and second, as a proton donor in the second reaction of acetyl group with the lipoate. We also identify I206alpha as a key residue in mediating the conformation of active-site loops. We propose that a simple conformational flip of the H271alpha side chain assists transfer of the acetyl group from thiamine cofactor to lipoyl domain in synchrony with reduction of the dithiolane ring.

Structural basis for membrane binding and catalytic activation of the peripheral membrane enzyme pyruvate oxidase from Escherichia coli

The thiamin- and flavin-dependent peripheral membrane enzyme pyruvate oxidase from E. coli catalyzes the oxidative decarboxylation of the central metabolite pyruvate to CO(2) and acetate. Concomitant reduction of the enzyme-bound flavin triggers membrane binding of the C terminus and shuttling of 2 electrons to ubiquinone 8, a membrane-bound mobile carrier of the electron transport chain. Binding to the membrane in vivo or limited proteolysis in vitro stimulate the catalytic proficiency by 2 orders of magnitude. The molecular mechanisms by which membrane binding and activation are governed have remained enigmatic. Here, we present the X-ray crystal structures of the full-length enzyme and a proteolytically activated truncation variant lacking the last 23 C-terminal residues inferred as important in membrane binding. In conjunction with spectroscopic results, the structural data pinpoint a conformational rearrangement upon activation that exposes the autoinhibitory C terminus, thereby freeing the active site. In the activated enzyme, Phe-465 swings into the active site and wires both cofactors for efficient electron transfer. The isolated C terminus, which has no intrinsic helix propensity, folds into a helical structure in the presence of micelles.

Structural evidence for substrate-induced synergism and half-sites reactivity in biotin carboxylase

Bacterial acetyl-CoA carboxylase is a multifunctional biotin-dependent enzyme that consists of three separate proteins: biotin carboxylase (BC), biotin carboxyl carrier protein (BCCP), and carboxyltransferase (CT). Acetyl-CoA carboxylase is a potentially attractive target for novel antibiotics because it catalyzes the first committed step in fatty acid biosynthesis. In the first half-reaction, BC catalyzes the ATP-dependent carboxylation of BCCP. In the second half-reaction, the carboxyl group is transferred from carboxybiotinylated BCCP to acetyl-CoA to produce malonyl-CoA. A series of structures of BC from several bacteria crystallized in the presence of various ATP analogs is described that addresses three major questions concerning the catalytic mechanism. The structure of BC bound to AMPPNP and the two catalytically essential magnesium ions resolves inconsistencies between the kinetics of active-site BC mutants and previously reported BC structures. Another structure of AMPPNP bound to BC shows the polyphosphate chain folded back on itself, and not in the correct (i.e., extended) conformation for catalysis. This provides the first structural evidence for the hypothesis of substrate-induced synergism, which posits that ATP binds nonproductively to BC in the absence of biotin. The BC homodimer has been proposed to exhibit half-sites reactivity where the active sites alternate or "flip-flop" their catalytic cycles. A crystal structure of BC showed the ATP analog AMPPCF(2)P bound to one subunit while the other subunit was unliganded. The liganded subunit was in the closed or catalytic conformation while the unliganded subunit was in the open conformation. This provides the first structural evidence for half-sites reactivity in BC.

Extended polypeptide linkers establish the spatial architecture of a pyruvate dehydrogenase multienzyme complex

Icosahedral pyruvate dehydrogenase (PDH) enzyme complexes are molecular machines consisting of a central E2 core decorated by a shell of peripheral enzymes (E1 and E3) found localized at a distance of approximately 75-90 A from the core. Using a combination of biochemical, biophysical, and cryo-electron microscopic techniques, we show here that the gap between the E2 core and the shell of peripheral enzymes is maintained by the flexible but extended conformation adopted by 60 linker polypeptides that radiate outwards from the inner E2 core, irrespective of the E1 or E3 occupancy. The constancy of the gap is thus not due to protein-protein interactions in the outer protein shell. The extended nature of the E2 inner-linker regions thereby creates the restricted annular space in which the lipoyl domains of E2 that carry catalytic intermediates shuttle between E1, E2, and E3 active sites, while their conformational flexibility facilitates productive encounters.

The crystal structure of phenylpyruvate decarboxylase from Azospirillum brasilense at 1.5 Å resolution

Phenylpyruvate decarboxylase (PPDC) of Azospirillum brasilense, involved in the biosynthesis of the plant hormone indole-3-acetic acid and the antimicrobial compound phenylacetic acid, is a thiamine diphosphate-dependent enzyme that catalyses the nonoxidative decarboxylation of indole- and phenylpyruvate. Analogous to yeast pyruvate decarboxylases, PPDC is subject to allosteric substrate activation, showing sigmoidal v versus [S] plots. The present paper reports the crystal structure of this enzyme determined at 1.5 A resolution. The subunit architecture of PPDC is characteristic for other members of the pyruvate oxidase family, with each subunit consisting of three domains with an open alpha/beta topology. An active site loop, bearing the catalytic residues His112 and His113, could not be modelled due to flexibility. The biological tetramer is best described as an asymmetric dimer of dimers. A cysteine residue that has been suggested as the site for regulatory substrate binding in yeast pyruvate decarboxylase is not conserved, requiring a different mechanism for allosteric substrate activation in PPDC. Only minor changes occur in the interactions with the cofactors, thiamine diphosphate and Mg2+, compared to pyruvate decarboxylase. A greater diversity is observed in the substrate binding pocket accounting for the difference in substrate specificity. Moreover, a catalytically important glutamate residue conserved in nearly all decarboxylases is replaced by a leucine in PPDC. The consequences of these differences in terms of the catalytic and regulatory mechanism of PPDC are discussed.

The two active sites in human branched-chain alpha-keto acid dehydrogenase operate independently without an obligatory alternating-site mechanism

A long standing controversy is whether an alternating activesite mechanism occurs during catalysis in thiamine diphosphate (ThDP)-dependent enzymes. We address this question by investigating the ThDP-dependent decarboxylase/dehydrogenase (E1b) component of the mitochondrial branched-chain alpha-keto acid dehydrogenase complex (BCKDC). Our crystal structure reveals that conformations of the two active sites in the human E1b heterotetramer harboring the reaction intermediate are identical. Acidic residues in the core of the E1b heterotetramer, which align with the proton-wire residues proposed to participate in active-site communication in the related pyruvate dehydrogenase from Bacillus stearothermophilus, are mutated. Enzyme kinetic data show that, except in a few cases because of protein misfolding, these alterations are largely without effect on overall activity of BCKDC, ruling out the requirement of a proton-relay mechanism in E1b. BCKDC overall activity is nullified at 50% phosphorylation of E1b, but it is restored to nearly half of the pre-phosphorylation level after dissociation and reconstitution of BCKDC with the same phosphorylated E1b. The results suggest that the abolition of overall activity likely results from the specific geometry of the half-phosphorylated E1b in the BCKDC assembly and not due to a disruption of the alternating active-site mechanism. Finally, we show that a mutant E1b containing only one functional active site exhibits half of the wild-type BCKDC activity, which directly argues against the obligatory communication between active sites. The above results provide evidence that the two active sites in the E1b heterotetramer operate independently during the ThDP-dependent decarboxylation reaction.

A novel branched-chain amino acid metabolon. Protein-protein interactions in a supramolecular complex

The catabolic pathways of branched-chain amino acids have two common steps. The first step is deamination catalyzed by the vitamin B(6)-dependent branched-chain aminotransferase isozymes (BCATs) to produce branched-chain alpha-keto acids (BCKAs). The second step is oxidative decarboxylation of the BCKAs mediated by the branched-chain alpha-keto acid dehydrogenase enzyme complex (BCKD complex). The BCKD complex is organized around a cubic core consisting of 24 lipoate-bearing dihydrolipoyl transacylase (E2) subunits, associated with the branched-chain alpha-keto acid decarboxylase/dehydrogenase (E1), dihydrolipoamide dehydrogenase (E3), BCKD kinase, and BCKD phosphatase. In this study, we provide evidence that human mitochondrial BCAT (hBCATm) associates with the E1 decarboxylase component of the rat or human BCKD complex with a K(D) of 2.8 microM. NADH dissociates the complex. The E2 and E3 components do not interact with hBCATm. In the presence of hBCATm, k(cat) values for E1-catalyzed decarboxylation of the BCKAs are enhanced 12-fold. Mutations of hBCATm proteins in the catalytically important CXXC center or E1 proteins in the phosphorylation loop residues prevent complex formation, indicating that these regions are important for the interaction between hBCATm and E1. Our results provide evidence for substrate channeling between hBCATm and BCKD complex and formation of a metabolic unit (termed branched-chain amino acid metabolon) that can be influenced by the redox state in mitochondria.

The 1',4'-iminopyrimidine tautomer of thiamin diphosphate is poised for catalysis in asymmetric active centers on enzymes

Thiamin diphosphate, a key coenzyme in sugar metabolism, is comprised of the thiazolium and 4'-aminopyrimidine aromatic rings, but only recently has participation of the 4'-aminopyrimidine moiety in catalysis gained wider acceptance. We report the use of electronic spectroscopy to identify the various tautomeric forms of the 4'-aminopyrimidine ring on four thiamin diphosphate enzymes, all of which decarboxylate pyruvate: the E1 component of human pyruvate dehydrogenase complex, the E1 subunit of Escherichia coli pyruvate dehydrogenase complex, yeast pyruvate decarboxylase, and pyruvate oxidase from Lactobacillus plantarum. It is shown that, according to circular dichroism spectroscopy, both the 1',4'-iminopyrimidine and the 4'-aminopyrimidine tautomers coexist on the E1 component of human pyruvate dehydrogenase complex and pyruvate oxidase. Because both tautomers are seen simultaneously, these two enzymes provide excellent evidence for nonidentical active centers (asymmetry) in solution in these multimeric enzymes. Asymmetry of active centers can also be induced upon addition of acetylphosphinate, an excellent electrostatic pyruvate mimic, which participates in an enzyme-catalyzed addition to form a stable adduct, resembling the common predecarboxylation thiamin-bound intermediate, which exists in its 1',4'-iminopyrimidine form. The identification of the 1',4'-iminopyrimidine tautomer on four enzymes is almost certainly applicable to all thiamin diphosphate enzymes: this tautomer is the intramolecular trigger to generate the reactive ylide/carbene at the thiazolium C2 position in the first fundamental step of thiamin catalysis.

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.