Metabolites produced via traditional biochemical processes affect intracellular communication, inflammation, and malignancy. Unexpectedly, acetyl-CoA, α-ketoglutarate and palmitic acid, which are chemical species of reactions catalyzed by highly abundant, gigantic enzymatic complexes, dubbed as “metabolons”, have broad “nonmetabolic” signaling functions. Conserved unstructured regions within metabolons determine the yield of these metabolites. Unstructured regions tether functional protein domains, act as spatial constraints to confine constituent enzyme communication, and, in the case of acetyl-CoA production, tend to be regulated by intricate phosphorylation patterns. This review presents the multifaceted roles of these three significant metabolites and describes how their perturbation leads to altered or transformed cellular function. Their dedicated enzymatic systems are then introduced, namely, the pyruvate dehydrogenase (PDH) and oxoglutarate dehydrogenase (OGDH) complexes, and the fatty acid synthase (FAS), with a particular focus on their structural characterization and the localization of unstructured regions. Finally, upstream metabolite regulation, in which spatial occupancy of unstructured regions within dedicated metabolons may affect metabolite availability and subsequently alter cell functions, is discussed.
In the second part of the twentieth century, there have been extensive studies that were mainly focused on discovering the molecular mechanisms behind various diseases, and specifically malignancies. Gene knockout, sequencing and, later, −omics methods, as well as network biology, revealed the different underlying pathways that, when perturbed due to various external or internal factors, are the causes of numerous pathological conditions. This collective work led to the creation of a basic paradigm: a cell will either use an extracellular or intracellular “sensor” to measure changes in the environment and trigger a signaling network. This signal transduction network will propagate via protein-protein interactions, until the signal reaches the nucleus, consequently altering gene expression and triggering a coordinated response to altered conditions, similar to an electric circuit triggered by an on/off switch.
Recent studies, however, have illuminated an aspect of cell signaling that had been mostly overlooked. The effectors of signaling pathways not only can be proteins, but also metabolites, the products or byproducts of cellular metabolism. Earlier observations, e.g., those showing the control of gene expression in bacteria by changes in glucose and lactose concentrations  and the covalent modifications of proteins and their implications on protein function in rats , hinted towards functions for metabolites in signaling. The current view of cell function is that cellular signaling and metabolism are tightly interconnected [23, 79] and that metabolic networks are as complex as their signaling counterparts; that is, they have an equally intricate role, not only by producing and storing energy or metabolizing proteins, nucleic acids, lipids, and sugars (or to catabolizing them), but also by signaling through their substrates, intermediates, and products to alter cellular behavior and fate.
This review focuses on three different metabolites, acetyl-CoA, α-ketoglutarate and, palmitic acid, and examines how their altered intracellular concentrations can affect various signaling pathways and cellular processes. In addition, more details are presented to offer perspective on the large enzymatic complexes, dubbed “metabolons”, which are the main players in the production and regulation of these metabolites, namely the pyruvate dehydrogenase complex (PDHc), the 2-oxoglutarate dehydrogenase complex (OGDHc), and fatty acid synthase (FAS). Finally, an effort is made to describe unstructured regions in the metabolons and to explain how they may have a potential effect on the availability and regulation of these metabolites.
Nonmetabolic functions of acetyl-CoA, α-ketoglutarate and palmitic acid
Acetyl-CoA is involved in cell fate and survival and the regulation of protein-protein interactions via post-translational modifications
Acetyl-coenzyme A (ACoA) is the lynchpin of cellular metabolism. It is a central component of a plethora of cellular processes, a main intermediate of multiple metabolic pathways and reactions, and the critical cosubstrate of protein acetylation, among other modifications [56, 70]. ACoA is the molecule through which pyruvate is derived from glycolysis before entering the tricarboxylic acid (TCA) cycle and, along with malonyl-CoA, is a precursor for lipid synthesis . ACoA is comprised of an acetyl group that is connected via a thioester bond to coenzyme A, a molecule that in turn is constituted by 3′,5′-ADP, pantothenic acid (a derivative of vitamin B5) and β-mercaptoethylamine. ACoA is produced by multiple sources, e.g., the oxidation of fatty acids and degradation of amino acids, but the most important route of its biosynthesis is through pyruvate generated during glycolysis in a decarboxylation process carried out by PDHc.
Intracellular levels of ACoA are highly regulated, since changes in cytosolic concentration can greatly affect cell fate through multiple processes (Fig. 1). A reduction in cytosolic ACoA has been shown to induce autophagy . The depletion of ACoA activates the AMP-activated protein kinase, which in turn activates Unc-51-like autophagy activating kinase 1 (ULK1), an important kinase with many downstream phosphorylation targets that play critical roles in the formation of the autophagosomes . To further increase the intensity of the signal, reduced cytosolic ACoA levels also result in the inhibition of mechanistic target of rapamycin complex 1 (mTORC1), a known inhibitor of ULK1, thus preventing mTOR signaling, which inhibits autophagy under normal conditions . Additionally, the cytosolic concentration ratio of ACoA and CoA has other important effects on cell death regulation. An increased ratio can induce an ACoA-dependent metabolic signal that activates Ca2+/calmodulin-dependent protein kinase II (CaMKII), which in turn phosphorylates caspase-2, through which multiple anti-apoptotic responses are triggered .
The most important role of ACoA is arguably as the donor of the acetyl group needed for protein acetylation, which can either occur simultaneously with translation (N-term acetylation) or as posttranslational modifications (PTMs). In the case of N-terminal acetylation, the acetyl group is transferred to the amino group of the first N-terminal residue, usually Ser, Thr, Val, Cys, or Ala, after the initiating Met has been cleaved. This modification contributes not only to protein stability, but also to protein function and localization .
Protein acetylation has been extensively studied as a PTM of Lys residues. This addition of an acetyl group to the amino group of a lysine residue has a profound effect on a protein’s stability, catalytic activity, interaction with other proteins and localization (Fig. 1) . Protein acetylation is mediated by two main classes of enzymes, lysine acetyl-transferases (KATs), which add the group to proteins, and lysine deacetylases (KDACs), which reverse the modification by removing the acetyl- group. Histone acetyl-transferases/deacetylases (HATs/HDACs) are critical for modifying histones , thus altering the epigenetic profile by either “compacting” or “loosening” different chromatin areas and granting or denying access to the transcriptional machinery of the cell . Other interesting acetylation targets was tubulin subunits that build the microtubule cytoskeleton. Acetylation of Lys40 of α-tubulin is the only PTM that takes part inside the microtubule lumen, in contrast with other PTMs that only affect its outer surface . Even though this kind of modification has been known for some time, its impact is now only being understood, and recent studies have connected it to intracellular trafficking , cell migration and autophagy, among other processes . Finally, acetylation is of extreme importance for the function of p53, widely known as “the guardian of the genome” because of its role in tumor suppression. In response to various oncogenic stress factors, p53 activation and the transcription of its downstream targets is directly dependent on its activation by acetylation in multiple Lys sites located at its N-terminal domain .
α-Ketoglutarate is a promoter of cell survival and regulator of hypoxic signaling pathways
α-Ketoglutaric acid (AKG) (or, as previously known, 2-oxoglutaric acid) is a TCA cycle intermediate that is converted to succinate through decarboxylation by OGDHc. AKG conversion is a rate-determining step of the cycle, as anaplerotic reactions can increase its concentration, using glutamate as a source . AKG is an important substrate for the production of multiple amino acids, such as proline and leucine , and in general plays a critical role in cellular metabolic pathways by being an important source of glutamine and glutamate .
Apart from its direct role in energy production and biosynthetic pathways, AKG has also been shown to affect cell signaling pathways (Fig. 1) related to tumorigenesis and aging. In the first case, AKG acts as a substrate for prolyl-hydroxylases (PHDs), which constitute a class of α-ketoglutarate-dependent dioxygenases (AKGDDs); these enzymes require AKG and oxygen to hydroxylate hypoxia-inducible factor 1 alpha (HIF-1α) , producing succinate in the process. HIF-1α is one of the two subunits (with HIF-1β) of the HIF-1 transcription factor, the main player in cellular oxygen sensing . HIF-1, under hypoxic conditions, activates the transcription of genes related to glycolysis, reducing the dependence on oxidative phosphorylation and enabling the production of ATP in the absence of oxygen. Tumor cells, having to survive in hypoxic conditions, rely heavily on the activation of HIF-1 to satisfy their massive needs for available ATP . In the presence of oxygen, AKG activates prolyl hydroxylase 1–3, which in turn hydroxylates HIF-1α. Through this modification, the von Hippel-Lindau factor recognizes and binds to HIF-1α, marking the protein for proteasome-mediated degradation . Cancer cells take advantage of this process by reducing the levels of available AKG, thus inhibiting the action of AKGDDs, stabilizing HIF-1α and allowing it to migrate to the nucleus, dimerize with HIF-1β and activate the transcription of hypoxia related elements (HREs) , thus promoting cell survival and proliferation.
Recent studies have also highlighted the anti-aging effects of AKG . AKG can considerably increase the lifespan and delay the onset of age-related phenotypes of C. elegans through downstream inactivation of the TOR signaling pathway. More specifically, ATP synthase, which is a ubiquitous membrane-embedded enzymatic complex, related to energy metabolism and consequently of paramount importance for all cells, was found to be a novel binding target of AKG. AKG acts as an inhibitor of ATP synthase, reducing ATP cell levels, which in turn reduces oxygen consumption and increases autophagy . Additionally, a decrease in ATP levels changes the ADP/ATP ratio in cells, activating 5′ AMP-activated protein kinase (AMPK), which in turn initiates the phosphorylation of the TOR suppressor TSC2, further strengthening autophagy signaling .
Palmitic acid is a major player in the inflammatory response and regulator of protein localization through palmitoylation
Palmitic acid (PA) is the most common form of saturated fatty acid found inside the human body. ACoA and malonyl-CoA (MCoA) enter the fatty acid synthesis pathway, and MCoA is elongated through a repetitive synthesis cycle until it reaches a length of 16 carbons (16:0). In this process, fatty acid synthase (FAS) is a key enzyme . PA is consistently under tight homeostatic control, and extracellular intake levels are continually balanced by de novo lipogenesis (DNL), hinting at the critical role of PA in various cellular processes . Disturbances to intracellular concentrations can have adverse effects, as it is evident by the implication that PA is involved in various signaling processes associated with cancer phenotypes (Fig. 1). DNL is closely associated with cancer cells through the Warburg effect , through which cells show a high preference for anaerobic glycolysis, even before they reach the stage of oxygen deprivation due to reduced access to oxygenated blood.
Excess PA in cells increases diacylglycerol levels through its diversion to nonoxidative metabolic pathways. This diversion, in turn, activates protein kinase C (PKC) and upregulates the phosphorylation of insulin receptor substrate-1 (IRS-1) , lowering its activation and subseqently inhibiting the PI3K/AKT signaling pathway, which is a key pathway regulating the cell cycle and is thus directly related to cellular quiescence, proliferation, and cancer. Additionally, PA enhances the p38-mediated activation of PTEN , increasing its effect on AKT signaling. In adipose tissue, increased intracellular PA concentrations have been shown to facilitate inflammatory responses, leading to diseases such as insulin resistance and obesity. PA induces the activation of NF-κΒ-, protein kinase C-, and MAP kinase-mediated pathways, which in turn leads to the production of cytokines, such as TNF and interleukin-10, which indicates that cells are in a constant inflammatory state . This inflammatory state is further amplified by the additional phosphorylation of the JNK/ERK kinases, which increase the inflammatory response through their interaction with the MAPK pathway components and can cause pathological conditions such as metabolic syndrome . Interestingly, in the case of hepatocellular carcinoma, high levels of PA have been observed to have antitumor effects by reducing cell proliferation and invasiveness through the downregulation of the mTOR signaling pathway .
PA affects cellular processes not only through direct/indirect signaling but by also as the main substrate of another type of posttranslational modification, palmitoylation, through which a palmitate group (a saturated chain of 16 carbons) is attached to a cysteine residue of a protein target via a thioester bond in a process that also has the advantage of being completely reversible . This PTM has been shown to affect numerous functions (Fig. 1), including protein trafficking to the membrane and docking . Specifically, a palmitate group added to a soluble protein can serve as a hydrophobic anchor that facilitates protein docking to a membrane. Additionally, many transmembrane proteins, such as GPCRs, have been shown to undergo palmitoylation, thereby increasing their stability, preventing their aggregation, and modulating their target-specific binding to lipids or other proteins in the membrane . Palmitoylation has also been shown to modulate intercompartmental trafficking in the cell, as in the case of Ras proteins , regulating their transfer between Golgi and post-Golgi membrane compartments and further “fine-tune” other PTMs on the same protein, especially those in close proximity to the palmitoylation site .
Metabolons involved in the availability of acetyl-CoA, palmitic acid and α-ketoglutarate
Acetyl-CoA: produced by the giant PDHc with > 100 “flexible arms”
The acetylation of CoA is determined by the availability of carbon sources and can have either extramitochondrial or intramitochondrial origins. At low glucose levels, acetyl-CoA synthetase (ACS), with acetate, and alcohol dehydrogenase (ADH), with ethanol, acetylate CoA, coupled with ATP hydrolysis . The degradation of branched-chain ketogenic amino acids (e.g., valine, leucine, and isoleucine) leads to their conversion to α-ketoacids by transamination and eventually to isovaleryl-CoA through oxidative decarboxylation by the mitochondrial branched-chain α-ketoacid dehydrogenase complex (BCKDHc). Isovaleryl-CoA then undergoes dehydrogenation, carboxylation and hydration to form a CoA-derived intermediate before it is cleaved to ACoA and acetoacetate. In mitochondria, also at low glucose levels, the production of ACoA is linked to the β-oxidation of fatty acids. Fatty acids are converted to acyl-CoA and, subsequently, acyl-CoA is degraded in a four-step cycle by four respective enzymes, namely, acyl-CoA dehydrogenase, enoyl-CoA hydratase, 3-hydroxyacyl-CoA dehydrogenase, and thiolase. This enzymatic cycle leads to a new fatty acid chain with two fewer carbons and ACoA as a byproduct. At high glucose levels, glycolysis takes place rapidly, thus increasing the amount of citrate produced by the tricarboxylic acid cycle. Then, citrate is exported from mitochondria to other organelles to be catabolized into ACoA and oxaloacetate by ATP citrate lyase (ACL) coupled with the hydrolysis of ATP. Conversions between pyruvate and ACoA are possible; for example, pyruvate formate lyase converts pyruvate into ACoA and formic acid; however, while ACoA is mainly produced through glycolysis, pyruvate undergoes oxidative decarboxylation, during which it loses its carboxyl group (as carbon dioxide) to become ACoA. This pyruvate dehydrogenase reaction is catalyzed by the pyruvate dehydrogenase complex (PDHc), utilizing pyruvate to produce ACoA, CO2, and NADH.
The PDHc is a 10 MDa assembly, one of the largest soluble molecular machineries in the cell. It is composed of three basic enzymes, pyruvate decarboxylase [pyruvate dehydrogenase (lipoamide), E1p], dihydrolipoyl acetyltransferase (E2p) and dihydrolipoyl dehydrogenase (E3) (Fig. 2). In humans, the structural core of the PDHc is formed by multiple copies of the relevant E2p polypeptide chain and is arranged in pseudo-icosahedral (60-mer) symmetry similar to a regular dodecahedron (Fig. 1), with each vertex composed of an E2p trimer. The peripheral components E1p and E3 are also present in multiple copies, and bound tightly but noncovalently to the outer region of the core. E2p has a pronounced domain-and-linker structure (Figs. 2 and 3, Additional file 1: Table S1). It is composed of 4 ordered domains connected by flexible linker regions. The 4 ordered domains are directly related to the function of the PDHc, and include, from the N- to the C-terminus, two lipoyl domains (LD1 and LD2; 2× ~ 8 kDa), a structured region critical for binding the peripheral E1p subunits (peripheral subunit-binding domain, PSBD; ∼4 kDa), and the catalytic core module of the E2p (∼26 kDa) (Figs. 2 and 3, Additional file 1: Table S1). Ordered structures of human sequences have been resolved by means of X-ray crystallography [12,13,14], NMR spectroscopy  and cryo-EM . For ACoA production, pyruvate is captured and decarboxylated by E1p and then transferred to a mobile lipoyl domain in E2p. The lipoyl domain transports the acyl group to the E2p catalytic domain, which in turn transfers it to coenzyme A. ACoA is released, and the lipoyl domain of E2p is regenerated (reoxidized) for another cycle by E3. Therefore, the lipoyl domain acts as a “swinging arm” that brings substrates, intermediates and products from one enzymatic center to another. This arm includes the rigid lipoyl domain but is flexible overall, connected by structurally unresolved linker regions . These linker segments in the polypeptide chain are 25–30 amino acids long and act to facilitate the coupling of domain movement with an active site as part of the catalytic mechanism . These linker regions are remarkably predicted to include disorder, at least to a certain extent (Figs. 2 and 3a, Additional file 1: Table S1). Considering fully extended linker regions in the human PDHc (dCα-Cα = 3.8 Å), the maximum distance between the ordered domains is calculated to be long as 122 Å between the E2p core and the PSBD, 194 Å between PSBD and LD1 and 129 Å between LD1 and LD2. These distances seem very long, being the maximum for extended conformations; Experimental data suggest an overall diameter of the PDHc of ~ 500 Å, with the external diameter in a region of ~ 300 Å, including a gap between the core structure and the areas of external density . Flexible regions are suggested to be extended  and highly similar but divergent, with limited immunological cross-reactivity . Due to their high copy number, these regions are likely confined by the spatial proximity of the redundant N-ter domains, to prevent unfavorable van der Waals interactions. The linkers are shown to have disordered structures but are not random coils . The disorder predicted in the flexible parts of the E2p sequence (Fig. 3a, Additional file 1: Table S1) likely contributes to the regulated positioning of the LDs by selective conformational ensembles. For example, mutagenesis or constructs of variable lengths indeed were associated with poor growth rates in bacteria . Considering that 48–60 E2p proteins are present in each PDHc, the extent of the linker regions and the manifested disorder is greatly amplified. The effect of cumulated disordered regions on the resulting conformational ensemble of lipoyl domains is unknown, although the association of molecular redundancy, disordered regions, flexible structural parts, and ordered domains seems to be important. Recently performed molecular dynamics simulations of the full complex suggest significant flexibility changes imposed by neighboring E1p and E3 . Another protein that also contributes to the core structure of the PDHc, the E3 Binding Protein (E3BP), is present in ~ 12 copies, further complicating the ability to predict the overall structure critical for the PDHc reaction. E3BP is very similar to the E2p protein, except that it (a) does not actively acetylate CoA, (b) includes a region that binds specifically the E3 protein but not E1p protein and (c) includes 1 lipoyl domain not 2 lipoyl domains. It is also known to include highly flexible linkers, and those are also predicted to be disordered [7, 24] (Fig. 3a, Additional file 1: Table S1). The overall spatial architecture of the PDHc is highly intricate: 48–60 E1p and ~ 12 E3 proteins surround the E2p-E3BP core (composed of 48–60 E2p and ~ 6–12 E3BP proteins), while > 102 LDs are generally involved. It is reasonable to assume that the disorder manifested in the linker regions promotes the regulation of conformational freedom during the transfer of lipoyl intermediates. Therefore, the generation of ACoA, may depend on the lipoyl domain on and off rates, of its interaction with its surrounding polypeptide chains, and its underlying physical chemical and structural properties.
Another important layer to the regulation of the PDHc is based on posttranslational modifications , especially those conferred by its dedicated kinases and phosphatases, which inactivate the PDHc, a phenomenon observed after mitochondrial function is attenuated in cancer . The production of ACoA may be additionally inactivated by phosphorylation of the E1p at distinct Ser sites (Ser-264-α-P, Ser-271-α-P and Ser-203-α-P) present in the two flexible loops (loop A: from 259-α to 282-α, loop B: from 195-α to 205-α), each with a different affinity and specificity . The phosphorylation of loop A plays an integral role in channeling and anchoring thiamine diphosphate to the active site of E1p, whereas the phosphorylation of loop B coordinates a Mg2+ ion that is chelated by the diphosphate group of the cofactor . The disorder of these phosphorylation loops was well-defined by the authors, who interpreted the absence of electron density in the crystallographic structures as characteristic of disorder [36, 80]. The phosphorylation of the E1p loops causes a steric clash between introduced phosphoryl groups and nearby residues, abolishing the hydrogen-bonding network essential for maintaining ordered loop conformations . Surprisingly, the increased regulatory loop disorder is associated with the loss of E1p catalytic efficiency, a molecular cause of PDHc deficiency (PDCD) . Other modifications are also important, e.g. Tyr301 phosphorylation , which has been implicated as drivers of the Warburg effect and resistance to therapy. Overall, various reproducibly identified PTMs correlate with ordered, but also some with flexible regions, and surprisingly, few regions of predicted disorder in the PDHc are also prone to systematic PTM regulation, with currently unclear physiological effects (Fig. 3a, Additional file 1: Table S1).
Notably, the transamination of the amino acids leucine, isoleucine and valine produces short branched-chain 2-oxo acids that are oxidatively decarboxylated by a single multienzyme complex (BCKDHc) in an intermediate step during ACoA generation. The structure of the BCKDHc is comparable to that of the PDHc and the OGDHc. The OGDHc, a key enzyme in the TCA cycle, converts 2-oxoglutaric acid to acyl-CoA, with the release of CO2 and NADH (described below).
α-Ketoglutaric acid and the regulation of its availability by the flexible OGDHc
AKG is produced by the oxidative deamination of glutamate by glutamate dehydrogenase and the oxidative decarboxylation of isocitrate by isocitrate dehydrogenase. In addition, alanine transaminase converts AKG and L-alanine to L-glutamate and pyruvate, respectively, in a reversible process. Producing key intermediates in the Krebs cycle, anaplerotic reactions can replenish the cycle after the isocitrate step but before succinyl-CoA step by synthesizing AKG through the transamination of glutamate or through the action of glutamate dehydrogenase on glutamate. Another MDa multienzyme complex, the 2-oxoglutarate dehydrogenase complex (OGDHc), with architecture surprisingly similar to that of the PDHc is involved in the processing of AKG, as described above. It consists of 2-oxoglutarate decarboxylase (E1o), dihydrolipoyl succinyltransferase (E2o) and the same dihydrolipoyl dehydrogenase (E3) enzymes as in the PDHc (Figs. 2 and 3b, Additional file 1: Table S1); these subunits are assembled into a giant 4 MDa complex but with distinct structural alterations, the most prominent of which is the formation of an octahedral E2 core composed of 24 E2 enzymes, not 60 as in the eukaryotic PDHc (Fig. 1). In addition, the OGDHc lacks dedicated kinases and phosphatases, and does not have an E3BP protein equivalent to that of the PDHc . The lack of an E3BP and a previously undiscovered mechanism of peripheral subunit recruitment are of particular interest. It is unknown whether a sequence on the outer region of E2o plays a part in the recruitment and spatial confinement of the E1o/E3 proteins. Due to their inherent similarities, one could hypothesize that the mechanism of OGDHc action is similar to that observed for the PDHc and described above. If this mechanism follows the paradigm of the E2p or E3BP proteins, then a disorder-to-order transition likely occurs in the flexible part of the E2o subunit, allowing a similar function and compensating for the absence of a subunit whose role would be homologous to that of the E3BP. If the mechanism is different, then other molecules not yet been discovered may be involved in regulating the proximity of E1o or/and E3 to the E2o, or direct interactions of the peripheral subunits with the E2o might play a role. A recent discovery may shed some light on the topic, as a novel subunit of OGDHc was identified, Kdg4, which has a moonlighting function in mitochondria, being part of both the mitochondrial ribosome and the OGDHc . It binds the E1o-E2o core with the N-ter, and tethers the E3 subunit with its C-terminal domain .
Overall, E2o includes a currently structurally resolved E2 acyltransferase domain, which tethers the lipoyl domain to a flexible region, which is predicted to be disordered. Again, the lipoyl arm is tethered to this region, which confers the flexibility to undergo various conformational states and regulate intracomplex interactions . E1o is twice the size of the E1p in the PDHc, combining both dehydrogenase and transketolase functions through an unknown structure. However, E1o undergoes allosteric interactions that control OGDHc activity [54, 58].
The flexible lipoyl arm that putatively attaches to the dehydrogenase domain of the E1o is in a conformational space that enables its interaction with the E1o. Therefore, the key disordered region in the fully extended distance of 289 Å may regulate further processing of the AKG to acyl-CoA, in a similar fashion to that in the PDHc. The same questions regarding the concentrations of multiple flexible regions on the surface of the E2p are also applicable to the OGDHc, where 24 lipoyl domains must be coordinated to perform the transport of reaction intermediates. PTMs have been identified on both the E1o and E2o proteins, and few have been predicted to be localized on corresponding disordered regions (Fig. 3b, Additional file 1: Table S1). It is currently unclear how these PTMs might affect the function and/or dynamics/conformational variation of the flexible sequence. Interestingly, the OGDHc can also act as a posttranslational modifier: by producing succinyl-CoA, the OGDHc can also succinylate fumarase or the PDHc and thus regulate the abundance of other metabolites, such as those in neurons .
Palmitic acid and FAS form a highly modular enzymatic machine
Palmitic acid (PA) is a component in the body of animals, and in humans, it may constitute as much as 21–30% (molar) of human fat depots. It is also a major lipid component in human breast milk. These roles are based on carbohydrates being converted to PA, the first fatty acid produced during fatty acid synthesis and the precursor to longer fatty acids. As a consequence, PA negatively mediates acetyl-CoA carboxylase (ACC) in a feedback mechanism critical for converting acetyl-CoA to malonyl-CoA, which in turn adds carbons to the growing acyl chain, thus preventing further palmitate generation. The de novo generation of palmitate is controlled by another giant multienzyme machine, the fatty acid synthase complex (FAS, Figs. 1, 2 and 3c), which catalyzes the conversion of acetyl-CoA and malonyl-CoA into long-chain saturated fatty acids in the presence of NADPH . In prokaryotes, plants and mitochondria, fatty acid synthesis consists of 7 structurally independent monofunctional enzymes and an acyl carrier protein. In fungi and animals, however, all of the component enzymes and the acyl carrier protein are fused and organized in large polypeptide chains. FAS is an A6B6 complex in fungi and an A2 complex in humans. In humans, it is a surprisingly long polypeptide chain consisting of ~ 2500 residues with a mass of ~ 273 kDa. The complete structure of human FAS is not yet known, with only fragments of the overall architecture revealed, similar to the previously described PDHc and OGDHc metabolons, despite its overall smaller size (~ 0.5 MDa).
The FAS reaction starts when the acyl moiety of acetyl-CoA is transferred to the acyl carrier protein (ACP) by the malonyl−/acetyl-CoA-ACP-transacylase (MAT). MAT also transacylates the malonyl of the elongation substrate malonyl-CoA to ACP. β-ketoacyl synthase (KS) then causes the decarboxylative condensation of the acyl intermediate with malonyl-ACP converted to β-ketoacyl-ACP (or an acetoacetyl-ACP in the first elongation cycle). The β-carbon is then reduced through β-ketoacyl reductase (KR) in an NADPH-dependent manner, and the product (β-hydroxyacyl-ACP) is dehydrated by a dehydratase (DH) to become a β-enoyl, an intermediate that is then reduced by β-enoyl reductase (ER) in an NADPH-dependent manner to yield a C4 acyl substrate for the subsequent cyclic elongation. Two-carbon units are derived from MCoA until a length of C16 (palmitic acid) to C18 (stearic acid) is reached. The end product is released from the ACP by the thioesterase (TE).
It was previously thought that this cycle was not regulated by FAS itself, and that the binding of ACP to various sites was stochastic and asymmetric, as indicated by molecular dynamics simulation studies . The recent discovery of the γ subunit in yeast FAS showed that FAS enzymatic activity is regulated by the γ subunit in response to an abundance of its cosubstrate NADP and has impacts on the higher-order structure and consequently the conformational space probed by ACP . ACP (residues 2125–2192) is tethered to flexible regions on either side of its termini (residues 2066–2124 at the N-ter, and residues 2125–2241 at the C-ter), which have not yet been captured by crystallography or electron cryomicroscopy at high resolution. Additionally, mobile ACP is not frequently captured in the context of larger FAS structures [35, 38, 72], which points to a significant mobility caused by the connecting flexible regions. These flexible linkers are also predicted to include disordered regions and have only been seen at low resolution with cryo-EM . PTMs are widely present in the FAS structure, but currently their functions are unknown. However, we can still observe various identified PTMs on the ordered, flexible and disordered parts of the polypeptide chain (Fig. 3c).
Considering the similarity in the shuttling of intermediates between FAS, the PDHc and the OGDHc by lipoyl or acyl-carrier arms (Fig. 2), tethered to flexible regions with the aim of bringing the substrate from one enzymatic site to the other, it is reasonable to assume that the common feature of the manifested disorder observed in these complexes may contribute to the conformational space that is being explored by the flexible arms. Human and fungal FAS have two and six lipoyl arms, respectively, which are spatially confined, but can intercommunicate within their respective spatial constraints as defined by the length of the tethered flexible linker regions  in combination with the suppression of futile catalytic cycles by the γ subunit . The mechanisms of the higher-order regulation of multiple carrier proteins is also currently unknown, since they must be spatially confined in order to avoid unintended protein interactions; in the case of FAS in yeast or C. thermophilum, this is succeeded by a chamber-like structure [35, 42], but in the case of its human equivalent, the identity of the regulation mechanisms remains an open question because of the absence of resolved native FAS complexes. In the cases of the PDHc and OGDHc, regulatory confinement may be possible, but a mechanism driving this remains unknown due to the scarcity of structures from native PDHc and OGDHc complexes.
The evidence shows that multiple metabolites, namely, acetyl-CoA, α-ketoglutarate and palmitic acid, apart from their role in the cellular bioenergetic and biosynthetic pathways, can cause profound changes through their participation in multiple cell signaling pathways . They may affect processes, such as autophagy [4, 67] and inflammation [1, 30], and cell survival and proliferation , whereas their perturbation can lead to pathogenic phenotypes, associated with tumor survival and progression . Notably, these specific metabolites can affect the same metabolic pathways, but in different cell types, adding another level of regulation associated with tissue specificity; this is the case of the mTOR signaling pathway [40, 78], although no study has directly investigated this correlation.
All three metabolites investigated in this review are regulated by large and highly intricate enzymatic complexes  or “metabolons”, the pyruvate dehydrogenase complex, the 2-oxoglutarate dehydrogenase complex, and fatty acid synthase. Although these complexes may seem rigid and highly ordered, they hide a considerable amount of disorder, mainly in regions where conformational flexibility, one of the characteristics of intrinsically disordered regions, is of paramount importance for the facilitation of substrate accessibility, substrate channeling to multiple enzymes, enhanced protein-protein interactions and increased overall reaction speed . Disordered regions have displayed highly distinct sequence features throughout evolution and, in multiple cases, have been shown to contribute to the aforementioned characteristics, which are needed by protein complexes that require every structural advantage available to optimally regulate the availability of highly important metabolites . The flexible linker regions in PDHc and OGDHc E2 proteins are not conserved regarding sequence and length but retain a very high alanine/proline/serine content. These amino acids often underlie the residue content of the flexible regions. One may speculate that the disordered regions are not based on a highly conserved residue pattern, such as those indicated by active sites, but that the overall composition is critical for their function. Another intriguing feature of the flexible regions of PDHc and OGDHc is a relatively high Lys content (~ 10%). Lys can be spontaneously modified by acetyl-CoA , which might thus function as an additional regulator. The overproduction of acetyl-CoA modifies the flexible linker and alters its mobility via a concentration-dependent mechanism, which is not catalyzed by enzymes. Additionally, FAS shows a fascinating evolution path, from a multiprotein complex in bacteria to multidomain organization in mammals, which may correlate with a more controlled mechanism associated with intermediate substrate shuttling and regulation of its activity .
Nevertheless, structural studies that delve into disorder of large biomolecular complexes are limited by the very nature of both the disorder and the size of these large enzymatic complexes. Fortunately, advances in the field of structural and computational biology, as demonstrated in the field of electron cryomicroscopy , allow us to study protein complexes in their native state , preserving the characteristics that are not possible to observe with in vitro protein expression and purification. Studying the organization of the ordered parts of native protein complexes will reveal information on how the flexible, intrinsically disordered regions contribute to regulating enzyme kinetics, protein phosphorylation and disease development. To summarize, understanding the native state of a protein allows us to understand its cell function , and when this function is tightly bound to cellular fate, as is the case of ACoA, AKG and PA implicated in multiple cellular processes, insights such as these allow us to elucidate and counteract previously unknown disease mechanisms.
Availability of data and materials
Datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
ATP citrate lyase
Acyl- carrier protein
Protein kinase B
5′ AMP-activated protein kinase
Branched chain α-ketoacid dehydrogenase complex
Ca2+/calmodulin-dependent protein kinase II
Cryogenic electron microscopy
De novo lipogenesis
Pyruvate dehydrogenase (lipoamide)
E3 binding protein
Fatty acid synthase
G protein-coupled receptor
Insulin receptor substrate-1
C-Jun N-terminal kinase
α-ketoglutarate dehydrogenase complex
Mitogen-activated protein kinase
mechanistic target of rapamycin
mechanistic target of rapamycin complex 1
Nicotinamide adenine dinucleotide
Nicotinamide adenine dinucleotide phosphate
Nuclear factor ‘kappa-light-chain-enhancer’ of activated B-cells
Research collaboratory for structural bioinformatics
Tri Carboxylic Acid
Tumor necrosis factor
Target of rapamycin
Unc-51 like autophagy activating kinase 1
Ajuwon KM, Spurlock ME. Palmitate activates the NF-kappaB transcription factor and induces IL-6 and TNFalpha expression in 3T3-L1 adipocytes. J Nutr. 2005;135(8):1841–6. https://doi.org/10.1093/jn/135.8.1841.
Anselmi C, Grininger M, Gipson P, Faraldo-Gomez JD. Mechanism of substrate shuttling by the acyl-carrier protein within the fatty acid mega-synthase. J Am Chem Soc. 2010;132(35):12357–64. https://doi.org/10.1021/ja103354w.
Belew GD, Silva J, Rito J, Tavares L, Viegas I, Teixeira J, et al. Transfer of glucose hydrogens via acetyl-CoA, malonyl-CoA, and NADPH to fatty acids during de novo lipogenesis. J Lipid Res. 2019;60(12):2050–6. https://doi.org/10.1194/jlr.RA119000354.
Brautigam CA, Wynn RM, Chuang JL, Machius M, Tomchick DR, Chuang DT. Structural insight into interactions between dihydrolipoamide dehydrogenase (E3) and E3 binding protein of human pyruvate dehydrogenase complex. Structure. 2006;14(3):611–21. https://doi.org/10.1016/j.str.2006.01.001.
Chin RM, Fu X, Pai MY, Vergnes L, Hwang H, Deng G, et al. The metabolite alpha-ketoglutarate extends lifespan by inhibiting ATP synthase and TOR. Nature. 2014;510(7505):397–401. https://doi.org/10.1038/nature13264.
Choudhary C, Weinert BT, Nishida Y, Verdin E, Mann M. The growing landscape of lysine acetylation links metabolism and cell signalling. Nat Rev Mol Cell Biol. 2014;15(8):536–50. https://doi.org/10.1038/nrm3841.
Ciszak EM, Makal A, Hong YS, Vettaikkorumakankauv AK, Korotchkina LG, Patel MS. How dihydrolipoamide dehydrogenase-binding protein binds dihydrolipoamide dehydrogenase in the human pyruvate dehydrogenase complex. J Biol Chem. 2006;281(1):648–55. https://doi.org/10.1074/jbc.M507850200.
Devedjiev Y, Steussy CN, Vassylyev DG. Crystal structure of an asymmetric complex of pyruvate dehydrogenase kinase 3 with lipoyl domain 2 and its biological implications. J Mol Biol. 2007;370(3):407–16. https://doi.org/10.1016/j.jmb.2007.04.083.
Fan J, Kang HB, Shan C, Elf S, Lin R, Xie J, et al. Tyr-301 phosphorylation inhibits pyruvate dehydrogenase by blocking substrate binding and promotes the Warburg effect. J Biol Chem. 2014;289(38):26533–41. https://doi.org/10.1074/jbc.M114.593970.
Fatima S, Hu X, Gong RH, Huang C, Chen M, Wong HLX, et al. Palmitic acid is an intracellular signaling molecule involved in disease development. Cell Mol Life Sci. 2019;76(13):2547–57. https://doi.org/10.1007/s00018-019-03092-7.
Flaveny CA, Griffett K, El-Gendy Bel D, Kazantzis M, Sengupta M, Amelio AL, et al. Broad anti-tumor activity of a small molecule that selectively targets the Warburg effect and Lipogenesis. Cancer Cell. 2015;28(1):42–56. https://doi.org/10.1016/j.ccell.2015.05.007.
Gipson P, Mills DJ, Wouts R, Grininger M, Vonck J, Kuhlbrandt W. Direct structural insight into the substrate-shuttling mechanism of yeast fatty acid synthase by electron cryomicroscopy. Proc Natl Acad Sci U S A. 2010;107(20):9164–9. https://doi.org/10.1073/pnas.0913547107.
Green JD, Perham RN, Ullrich SJ, Appella E. Conformational studies of the interdomain linker peptides in the dihydrolipoyl acetyltransferase component of the pyruvate dehydrogenase multienzyme complex of Escherichia coli. J Biol Chem. 1992;267(33):23484–8 Retrieved from https://www.ncbi.nlm.nih.gov/pubmed/1429691.
Harris RA, Bowker-Kinley MM, Wu P, Jeng J, Popov KM. Dihydrolipoamide dehydrogenase-binding protein of the human pyruvate dehydrogenase complex. DNA-derived amino acid sequence, expression, and reconstitution of the pyruvate dehydrogenase complex. J Biol Chem. 1997;272(32):19746–51. https://doi.org/10.1074/jbc.272.32.19746.
Harrison AP, Pierzynowski SG. Biological effects of 2-oxoglutarate with particular emphasis on the regulation of protein, mineral and lipid absorption/metabolism, muscle performance, kidney function, bone formation and cancerogenesis, all viewed from a healthy ageing perspective state of the art--review article. J Physiol Pharmacol. 2008;59(Suppl 1):91–106 Retrieved from https://www.ncbi.nlm.nih.gov/pubmed/18802218.
Heublein M, Burguillos MA, Vogtle FN, Teixeira PF, Imhof A, Meisinger C, Ott M. The novel component Kgd4 recruits the E3 subunit to the mitochondrial alpha-ketoglutarate dehydrogenase. Mol Biol Cell. 2014;25(21):3342–9. https://doi.org/10.1091/mbc.E14-07-1178.
Honjo T, Nishizuka Y, Hayaishi O. Diphtheria toxin-dependent adenosine diphosphate ribosylation of aminoacyl transferase II and inhibition of protein synthesis. J Biol Chem. 1968;243(12):3553–5 Retrieved from https://www.ncbi.nlm.nih.gov/pubmed/4297784.
Jiang J, Baiesc FL, Hiromasa Y, Yu X, Hui WH, Dai X, et al. Atomic structure of the E2 inner Core of human pyruvate dehydrogenase complex. Biochemistry. 2018;57(16):2325–34. https://doi.org/10.1021/acs.biochem.8b00357.
Kalia YN, Brocklehurst SM, Hipps DS, Appella E, Sakaguchi K, Perham RN. The high-resolution structure of the peripheral subunit-binding domain of dihydrolipoamide acetyltransferase from the pyruvate dehydrogenase multienzyme complex of Bacillus stearothermophilus. J Mol Biol. 1993;230(1):323–41. https://doi.org/10.1006/jmbi.1993.1145.
Kastritis PL, O'Reilly FJ, Bock T, Li Y, Rogon MZ, Buczak K, et al. Capturing protein communities by structural proteomics in a thermophilic eukaryote. Mol Syst Biol. 2017;13(7):936. https://doi.org/10.15252/msb.20167412.
Kato M, Wynn RM, Chuang JL, Tso SC, Machius M, Li J, Chuang DT. Structural basis for inactivation of the human pyruvate dehydrogenase complex by phosphorylation: role of disordered phosphorylation loops. Structure. 2008;16(12):1849–59. https://doi.org/10.1016/j.str.2008.10.010.
Kyrilis FL, Meister A, Kastritis PL. Integrative biology of native cell extracts: a new era for structural characterization of life processes. Biol Chem. 2019;400(7):831–46. https://doi.org/10.1515/hsz-2018-0445.
Leibundgut M, Jenni S, Frick C, Ban N. Structural basis for substrate delivery by acyl carrier protein in the yeast fatty acid synthase. Science. 2007;316(5822):288–90. https://doi.org/10.1126/science.1138249.
Lengyel JS, Stott KM, Wu X, Brooks BR, Balbo A, Schuck P, et al. Extended polypeptide linkers establish the spatial architecture of a pyruvate dehydrogenase multienzyme complex. Structure. 2008;16(1):93–103. https://doi.org/10.1016/j.str.2007.10.017.
Lin L, Ding Y, Wang Y, Wang Z, Yin X, Yan G, et al. Functional lipidomics: Palmitic acid impairs hepatocellular carcinoma development by modulating membrane fluidity and glucose metabolism. Hepatology. 2017;66(2):432–48. https://doi.org/10.1002/hep.29033.
Lomakin IB, Xiong Y, Steitz TA. The crystal structure of yeast fatty acid synthase, a cellular machine with eight active sites working together. Cell. 2007;129(2):319–32. https://doi.org/10.1016/j.cell.2007.03.013.
Mahamid J, Pfeffer S, Schaffer M, Villa E, Danev R, Cuellar LK, et al. Visualizing the molecular sociology at the HeLa cell nuclear periphery. Science. 2016;351(6276):969–72. https://doi.org/10.1126/science.aad8857.
Miles JS, Guest JR, Radford SE, Perham RN. Investigation of the mechanism of active site coupling in the pyruvate dehydrogenase multienzyme complex of Escherichia coli by protein engineering. J Mol Biol. 1988;202(1):97–106. https://doi.org/10.1016/0022-2836(88)90522-0.
Mole DR, Maxwell PH, Pugh CW, Ratcliffe PJ. Regulation of HIF by the von Hippel-Lindau tumour suppressor: implications for cellular oxygen sensing. IUBMB Life. 2001;52(1–2):43–7. https://doi.org/10.1080/15216540252774757.
Muhling J, Tussing F, Nickolaus KA, Matejec R, Henrich M, Harbach H, et al. Effects of alpha-ketoglutarate on neutrophil intracellular amino and alpha-keto acid profiles and ROS production. Amino Acids. 2010;38(1):167–77. https://doi.org/10.1007/s00726-008-0224-5.
Palomer X, Pizarro-Delgado J, Barroso E, Vazquez-Carrera M. Palmitic and oleic acid: the Yin and Yang of fatty acids in type 2 diabetes mellitus. Trends Endocrinol Metab. 2018;29(3):178–90. https://doi.org/10.1016/j.tem.2017.11.009.
Perham RN. Domains, motifs, and linkers in 2-oxo acid dehydrogenase multienzyme complexes: a paradigm in the design of a multifunctional protein. Biochemistry. 1991;30(35):8501–12. https://doi.org/10.1021/bi00099a001.
Radford SE, Laue ED, Perham RN, Martin SR, Appella E. Conformational flexibility and folding of synthetic peptides representing an interdomain segment of polypeptide chain in the pyruvate dehydrogenase multienzyme complex of Escherichia coli. J Biol Chem. 1989;264(2):767–75 Retrieved from https://www.ncbi.nlm.nih.gov/pubmed/2642904.
Ryan DG, Murphy MP, Frezza C, Prag HA, Chouchani ET, O'Neill LA, Mills EL. Coupling Krebs cycle metabolites to signalling in immunity and cancer. Nat Metab. 2019;1:16–33. https://doi.org/10.1038/s42255-018-0014-7.
Saunier E, Benelli C, Bortoli S. The pyruvate dehydrogenase complex in cancer: an old metabolic gatekeeper regulated by new pathways and pharmacological agents. Int J Cancer. 2016;138(4):809–17. https://doi.org/10.1002/ijc.29564.
Simic Z, Weiwad M, Schierhorn A, Steegborn C, Schutkowski M. The varepsilon-amino Group of Protein Lysine Residues is Highly Susceptible to nonenzymatic acylation by several physiological acyl-CoA Thioesters. Chembiochem. 2015;16(16):2337–47. https://doi.org/10.1002/cbic.201500364.
Singh K, Graf B, Linden A, Sautner V, Urlaub H, Tittmann K, et al. Discovery of a regulatory subunit of the yeast fatty acid synthase. Cell. 2020;180(6):1130–43 e1120. https://doi.org/10.1016/j.cell.2020.02.034.
Sramek J, Nemcova-Furstova V, Kovar J. Kinase signaling in apoptosis induced by saturated fatty acids in pancreatic beta-cells. Int J Mol Sci. 2016;17(9). https://doi.org/10.3390/ijms17091400.
van der Lee R, Buljan M, Lang B, Weatheritt RJ, Daughdrill GW, Dunker AK, et al. Classification of intrinsically disordered regions and proteins. Chem Rev. 2014;114(13):6589–631. https://doi.org/10.1021/cr400525m.
Vatrinet R, Leone G, De Luise M, Girolimetti G, Vidone M, Gasparre G, Porcelli AM. The alpha-ketoglutarate dehydrogenase complex in cancer metabolic plasticity. Cancer Metab. 2017;5:3. https://doi.org/10.1186/s40170-017-0165-0.
Whitley MJ, Arjunan P, Nemeria NS, Korotchkina LG, Park YH, Patel MS, et al. Pyruvate dehydrogenase complex deficiency is linked to regulatory loop disorder in the alphaV138M variant of human pyruvate dehydrogenase. J Biol Chem. 2018;293(34):13204–13. https://doi.org/10.1074/jbc.RA118.003996.
Zhou J, Yang L, Ozohanics O, Zhang X, Wang J, Ambrus A, et al. A multipronged approach unravels unprecedented protein-protein interactions in the human 2-oxoglutarate dehydrogenase multienzyme complex. J Biol Chem. 2018;293(50):19213–27. https://doi.org/10.1074/jbc.RA118.005432.
We thank the members of the Kastritis laboratory for valuable discussions.
This work was supported by the Federal Ministry for Education and Research (BMBF, ZIK program) [grant number 03Z22HN23 to P.L.K.]; the European Regional Development Funds for Saxony-Anhalt [grant number EFRE: ZS/2016/04/78115 to P.L.K] and funding by DFG (project number 391498659, RTG 2467: “Intrinsically Disordered Proteins – Molecular Principles, Cellular Functions, and Diseases”). Open access funding provided by Projekt DEAL.
Authors and Affiliations
Interdisciplinary Research Center HALOmem, Charles Tanford Protein Center, Martin Luther University Halle-Wittenberg, Kurt-Mothes-Straße 3a, Halle/Saale, Germany
Ioannis Skalidis, Christian Tüting & Panagiotis L. Kastritis
Institute of Biochemistry and Biotechnology, Martin Luther University Halle-Wittenberg, Kurt-Mothes-Straße 3, Halle/Saale, Germany
Ioannis Skalidis, Christian Tüting & Panagiotis L. Kastritis
ZIK HALOmem, Martin Luther University Halle-Wittenberg, Biozentrum, Room A.2.14, Weinbergweg 22, 06120, Halle/Saale, Germany
Domain annotation, location of unresolved, resolved and disordered regions for each subunit of the PDH and OGDH complexes (A) and FAS (B). Each entry is characterized by the number of the starting and ending amino-acid. Methods of annotation of the regions is described in Fig. 3 and its legend.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
Skalidis, I., Tüting, C. & Kastritis, P.L. Unstructured regions of large enzymatic complexes control the availability of metabolites with signaling functions.
Cell Commun Signal18, 136 (2020). https://doi.org/10.1186/s12964-020-00631-9