X-ray structures of the three Lactococcus lactis dihydroxyacetone kinase subunits and of a transient intersubunit complex

Homozygous variant in TKFC abolishing triokinase activities is associated with isolated immunodeficiency

BACKGROUND

Triokinase and FMN cyclase (TKFC) is a bifunctional enzyme involved in fructose metabolism. Triokinase catalyses the phosphorylation of fructose-derived glyceraldehyde (GA) and exogenous dihydroxyacetone (DHA), while FMN cyclase generates cyclic FMN. TKFC regulates the antiviral immune response by interacting with IFIH1 (MDA5). Previously reported pathogenic variants in TKFC are associated with either a multisystemic disease or isolated hypotrichosis with loose anagen hairs.

METHODS

Whole-exome sequencing identified a homozygous novel variant in TKFC (c.1624G>A; p.Gly542Arg) in an individual with a complex primary immunodeficiency disorder. The variant was characterised using enzymatic assays and yeast studies of mutant recombinant proteins.

RESULTS

The individual presented with chronic active Epstein-Barr virus disease and multiple bacterial and viral infections. Clinical investigations revealed hypogammaglobulinaemia, near absent natural killer cells and decreased memory B cells. Enzymatic assays showed that this variant displayed defective DHA and GA kinase activity while maintaining FMN cyclase activity. An allogenic bone marrow transplantation corrected the patient's immunodeficiency.

CONCLUSION

Our report suggests that TKFC may have a role in the immunological system. The pathological features associated with this variant are possibly linked with DHA/GA kinase inactivation through a yet an unknown mechanism. This report thus adds a possible new pathway of immunometabolism to explore further.

Crystal structure of a nucleotide-binding domain of fatty acid kinase FakA from Thermus thermophilus HB8

Fatty acid kinase is necessary for the incorporation of exogenous fatty acids into membrane phospholipids. Fatty acid kinase consists of two components: a kinase component, FakA, that phosphorylates a fatty acid bound to a fatty acid-binding component, FakB. However, the molecular details underlying the phosphotransfer reaction remain to be resolved. We determined the crystal structure of the N-terminal domain of FakA bound to ADP from Thermus thermophilus HB8. The overall structure of this domain showed that the helical barrel fold is similar to the nucleotide-binding component of dihydroxyacetone kinase. The structure of the nucleotide-binding site revealed the roles of the conserved residues in recognition of ADP and Mg²⁺, but the N-terminal domain of FakA lacked the ADP-capping loop found in the dihydroxyacetone kinase component. Based on the structural similarity to the two subunits of dihydroxyacetone kinase complex, we constructed a model of the complex of T. thermophilus FakB and the N-terminal domain of FakA. In this model, the invariant Arg residue of FakB occupied a position that was spatially similar to that of the catalytically important Arg residue of dihydroxyacetone kinase, which predicted a composite active site in the Fatty acid kinase complex.

Domain architecture and catalysis of the Staphylococcus aureus fatty acid kinase

Fatty acid kinase (Fak) is a two-component enzyme that generates acyl-phosphate for phospholipid synthesis. Fak consists of a kinase domain protein (FakA) that phosphorylates a fatty acid enveloped by a fatty acid binding protein (FakB). The structural basis for FakB function has been established, but little is known about FakA. Here, we used limited proteolysis to define three separate FakA domains: the amino terminal FakA_N, the central FakA_L, and the carboxy terminal FakA_C. The isolated domains lack kinase activity, but activity is restored when FakA_N and FakA_L are present individually or connected as FakA_NL. The X-ray structure of the monomeric FakA_N captures the product complex with ADP and two Mg2+ ions bound at the nucleotide site. The FakA_L domain encodes the dimerization interface along with conserved catalytic residues Cys240, His282, and His284. AlphaFold analysis of FakA_L predicts the catalytic residues are spatially clustered and pointing away from the dimerization surface. Furthermore, the X-ray structure of FakA_C shows that it consists of two subdomains that are structurally related to FakB. Analytical ultracentrifugation demonstrates that FakA_C binds FakB, and site-directed mutagenesis confirms that a positively charged wedge on FakB meshes with a negatively charged groove on FakA_C. Finally, small angle X-ray scattering analysis is consistent with freely rotating FakA_N and FakA_C domains tethered by flexible linkers to FakA_L. These data reveal specific roles for the three independently folded FakA protein domains in substrate binding and catalysis.

Closure of the Human TKFC Active Site: Comparison of the Apoenzyme and the Complexes Formed with Either Triokinase or FMN Cyclase Substrates

Human triokinase/flavin mononucleotide (FMN) cyclase (hTKFC) catalyzes the adenosine triphosphate (ATP)-dependent phosphorylation of D-glyceraldehyde and dihydroxyacetone (DHA), and the cyclizing splitting of flavin adenine dinucleotide (FAD). hTKFC structural models are dimers of identical subunits, each with two domains, K and L, with an L2-K1-K2-L1 arrangement. Two active sites lie between L2-K1 and K2-L1, where triose binds K and ATP binds L, although the resulting ATP-to-triose distance is too large (≈14 Å) for phosphoryl transfer. A 75-ns trajectory of molecular dynamics shows considerable, but transient, ATP-to-DHA approximations in the L2-K1 site (4.83 Å or 4.16 Å). To confirm the trend towards site closure, and its relationship to kinase activity, apo-hTKFC, hTKFC:2DHA:2ATP and hTKFC:2FAD models were submitted to normal mode analysis. The trajectory of hTKFC:2DHA:2ATP was extended up to 160 ns, and 120-ns trajectories of apo-hTKFC and hTKFC:2FAD were simulated. The three systems were comparatively analyzed for equal lengths (120 ns) following the principles of essential dynamics, and by estimating site closure by distance measurements. The full trajectory of hTKFC:2DHA:2ATP was searched for in-line orientations and short distances of DHA hydroxymethyl oxygens to ATP γ-phosphorus. Full site closure was reached only in hTKFC:2DHA:2ATP, where conformations compatible with an associative phosphoryl transfer occurred in L2-K1 for significant trajectory time fractions.

Theoretical Study of the Phosphoryl Transfer Reaction from ATP to Dha Catalyzed by DhaK from Escherichia coli

Protein kinases, representing one of the largest protein families involved in almost all aspects of cell life, have become one of the most important targets for the development of new drugs to be used in, for instance, cancer treatments. In this article an exhaustive theoretical study of the phosphoryl transfer reaction from adenosine triphosphate (ATP) to dihydroxyacetone (Dha) catalyzed by DhaK from Escherichia coli (E. coli) is reported. Two different mechanisms, previously proposed for the phosphoryl transfer from ATP to the hydroxyl side chain of specific serine, threonine, or tyrosine residues, have been explored based on the generation of free energy surfaces (FES) computed with hybrid QM/MM potentials. The results suggest that the substrate-assisted phosphoryl and proton-transfer mechanism is kinetically more favorable than the mechanism where an aspartate would be activating the Dha. Although the details of the mechanisms appear to be dramatically dependent on the level of theory employed in the calculations (PM3/MM, B3LYP:PM3/MM, or B3LYP/MM), the transition states (TSs) for the phosphoryl transfer step appear to be described as a concerted step with different degrees of synchronicity in the breaking and forming bonds process in both explored mechanisms. Residues of the active site belonging to different subunits of the protein, such as Gly78B, Thr79A, Ser80A, Arg178B, and one Mg²⁺ cation, would be stabilizing the transferred phosphate in the TS. Asp109A would have a structural role by posing the Dha and other residues of the active site in the proper orientation. The information derived from our calculations not only reveals the role of the enzyme and the particular residues of its active site, but it can assist in the rational design of new more specific inhibitors.

Glycerol metabolism of haloarchaea

Haloarchaea are heterotrophic members of the Archaea that thrive in hypersaline environments, often feeding off the glycerol that is produced as an osmolyte by eucaryotic Dunaliella during primary production. In this study we analyzed glycerol metabolism genes in closed genomes of haloarchaea and examined published data describing the growth properties of haloarchaea and experimental data for the enzymes involved. By integrating the genomic data with knowledge from the literature, we derived an understanding of the ecophysiology and evolutionary properties of glycerol catabolic pathways in haloarchaea.

Biochemical Roles for Conserved Residues in the Bacterial Fatty Acid-binding Protein Family

Fatty acid kinase (Fak) is a ubiquitous Gram-positive bacterial enzyme consisting of an ATP-binding protein (FakA) that phosphorylates the fatty acid bound to FakB. In Staphylococcus aureus, Fak is a global regulator of virulence factor transcription and is essential for the activation of exogenous fatty acids for incorporation into phospholipids. The 1.2-Å x-ray structure of S. aureus FakB2, activity assays, solution studies, site-directed mutagenesis, and in vivo complementation were used to define the functions of the five conserved residues that define the FakB protein family (Pfam02645). The fatty acid tail is buried within the protein, and the exposed carboxyl group is bound by a Ser-93-fatty acid carboxyl-Thr-61-His-266 hydrogen bond network. The guanidinium of the invariant Arg-170 is positioned to potentially interact with a bound acylphosphate. The reduced thermal denaturation temperatures of the T61A, S93A, and H266A FakB2 mutants illustrate the importance of the hydrogen bond network in protein stability. The FakB2 T61A, S93A, and H266A mutants are 1000-fold less active in the Fak assay, and the R170A mutant is completely inactive. All FakB2 mutants form FakA(FakB2)2 complexes except FakB2(R202A), which is deficient in FakA binding. Allelic replacement shows that strains expressing FakB2 mutants are defective in fatty acid incorporation into phospholipids and virulence gene transcription. These conserved residues are likely to perform the same critical functions in all bacterial fatty acid-binding proteins.

Identification of a two-component fatty acid kinase responsible for host fatty acid incorporation by Staphylococcus aureus

Extracellular fatty acid incorporation into the phospholipids of Staphylococcus aureus occurs via fatty acid phosphorylation. We show that fatty acid kinase (Fak) is composed of two dissociable protein subunits encoded by separate genes. FakA provides the ATP binding domain and interacts with two distinct FakB proteins to produce acyl-phosphate. The FakBs are fatty acid binding proteins that exchange bound fatty acid/acyl-phosphate with fatty acid/acyl-phosphate presented in detergent micelles or liposomes. The ΔfakA and ΔfakB1 ΔfakB2 strains were unable to incorporate extracellular fatty acids into phospholipid. FakB1 selectively bound saturated fatty acids whereas FakB2 preferred unsaturated fatty acids. Affymetrix array showed a global perturbation in the expression of virulence genes in the ΔfakA strain. The severe deficiency in α-hemolysin protein secretion in ΔfakA and ΔfakB1 ΔfakB2 mutants coupled with quantitative mRNA measurements showed that fatty acid kinase activity was required to support virulence factor transcription. These data reveal the function of two conserved gene families, their essential role in the incorporation of host fatty acids by Gram-positive pathogens, and connects fatty acid kinase to the regulation of virulence factor transcription in S. aureus.

Bifunctional homodimeric triokinase/FMN cyclase: contribution of protein domains to the activities of the human enzyme and molecular dynamics simulation of domain movements

Mammalian triokinase, which phosphorylates exogenous dihydroxyacetone and fructose-derived glyceraldehyde, is neither molecularly identified nor firmly associated to an encoding gene. Human FMN cyclase, which splits FAD and other ribonucleoside diphosphate-X compounds to ribonucleoside monophosphate and cyclic X-phosphodiester, is identical to a DAK-encoded dihydroxyacetone kinase. This bifunctional protein was identified as triokinase. It was modeled as a homodimer of two-domain (K and L) subunits. Active centers lie between K1 and L2 or K2 and L1: dihydroxyacetone binds K and ATP binds L in different subunits too distant (≈ 14 Å) for phosphoryl transfer. FAD docked to the ATP site with ribityl 4'-OH in a possible near-attack conformation for cyclase activity. Reciprocal inhibition between kinase and cyclase reactants confirmed substrate site locations. The differential roles of protein domains were supported by their individual expression: K was inactive, and L displayed cyclase but not kinase activity. The importance of domain mobility for the kinase activity of dimeric triokinase was highlighted by molecular dynamics simulations: ATP approached dihydroxyacetone at distances below 5 Å in near-attack conformation. Based upon structure, docking, and molecular dynamics simulations, relevant residues were mutated to alanine, and kcat and Km were assayed whenever kinase and/or cyclase activity was conserved. The results supported the roles of Thr(112) (hydrogen bonding of ATP adenine to K in the closed active center), His(221) (covalent anchoring of dihydroxyacetone to K), Asp(401) and Asp(403) (metal coordination to L), and Asp(556) (hydrogen bonding of ATP or FAD ribose to L domain). Interestingly, the His(221) point mutant acted specifically as a cyclase without kinase activity.

Dihydroxyacetone metabolism in Haloferax volcanii

Dihydroxyacetone (DHA) is a ketose sugar that can be produced by oxidizing glycerol. DHA in the environment is taken up and phosphorylated to DHA-phosphate by glycerol kinase or DHA kinase. In hypersaline environments, it is hypothesized that DHA is produced as an overflow product from glycerol utilization by organisms such as Salinibacter ruber. Previous research has demonstrated that the halobacterial species Haloquadratum walsbyi can use DHA as a carbon source, and putative DHA kinase genes were hypothesized to be involved in this process. However, DHA metabolism has not been demonstrated in other halobacterial species, and the role of the DHA kinase genes was not confirmed. In this study, we examined the metabolism of DHA in Haloferax volcanii because putative DHA kinase genes were annotated in its genome, and it has an established genetic system to assay growth of mutant knockouts. Experiments in which Hfx. volcanii was grown on DHA as the sole carbon source demonstrated growth, and that it is concentration dependent. Three annotated DHA kinase genes (HVO_1544, HVO_1545, and HVO_1546), which are homologous to the putative DHA kinase genes present in Hqm. walsbyi, as well as the glycerol kinase gene (HVO_1541), were deleted to examine the effect of these genes on the growth of Hfx. volcanii on DHA. Experiments demonstrated that the DHA kinase deletion mutant exhibited diminished, but not absence of growth on DHA compared to the parent strain. Deletion of the glycerol kinase gene also reduced growth on DHA, and did so more than deletion of the DHA kinase. The results indicate that Hfx. volcanii can metabolize DHA and that DHA kinase plays a role in this metabolism. However, the glycerol kinase appears to be the primary enzyme involved in this process. BLASTp analyses demonstrate that the DHA kinase genes are patchily distributed among the Halobacteria, whereas the glycerol kinase gene is widely distributed, suggesting a widespread capability for DHA metabolism.

Channel-mediated lactic acid transport: a novel function for aquaglyceroporins in bacteria

MIPs (major intrinsic proteins), also known as aquaporins, are membrane proteins that channel water and/or uncharged solutes across membranes in all kingdoms of life. Considering the enormous number of different bacteria on earth, functional information on bacterial MIPs is scarce. In the present study, six MIPs [glpF1 (glycerol facilitator 1)–glpF6] were identified in the genome of the Gram-positive lactic acid bacterium Lactobacillus plantarum. Heterologous expression in Xenopus laevis oocytes revealed that GlpF2, GlpF3 and GlpF4 each facilitated the transmembrane diffusion of water, dihydroxyacetone and glycerol. As several lactic acid bacteria have GlpFs in their lactate racemization operon (GlpF1/F4 phylogenetic group), their ability to transport this organic acid was tested. Both GlpF1 and GlpF4 facilitated the diffusion of D/L-lactic acid. Deletion of glpF1 and/or glpF4 in Lb. plantarum showed that both genes were involved in the racemization of lactic acid and, in addition, the double glpF1 glpF4 mutant showed a growth delay under conditions of mild lactic acid stress. This provides further evidence that GlpFs contribute to lactic acid metabolism in this species. This lactic acid transport capacity was shown to be conserved in the GlpF1/F4 group of Lactobacillales. In conclusion, we have functionally analysed the largest set of bacterial MIPs and demonstrated that the lactic acid membrane permeability of bacteria can be regulated by aquaglyceroporins.

Aerobic glycerol dissimilation via the Enterococcus faecalis DhaK pathway depends on NADH oxidase and a phosphotransfer reaction from PEP to DhaK via EIIADha

Two pathways for glycerol dissimilation are present in Enterococcus faecalis. Either glycerol is first phosphorylated by glycerol kinase and then oxidized by glycerol-3-phosphate oxidase with molecular oxygen as the electron acceptor (GlpO/GlpK pathway), or it is first oxidized by glycerol dehydrogenase with NAD(+) as the acceptor of the reduction equivalents and then phosphorylated by dihydroxyacetone kinase (GldA/DhaK pathway). The final end product in both cases is dihydroxyacetone phosphate (DHAP). The genes of the GldA/DhaK pathway are present in a four-gene operon structure encoding GldA, a small hypothetical protein (EF1359), and two subunits of dihydroxyacetone kinase (DhaK and DhaL). We demonstrate in this study that protein EF1359 is part of a phosphorylation cascade which phosphorylates dihydroxyacetone in a phosphoenolpyruvate (PEP)-dependent reaction via EI, HPr, EF1359 and DhaLK. Furthermore we show that aerobic dissimilation of glycerol via the GldA/DhaK pathway is dependent on active NADH oxidase to regenerate NADH in Ent. faecalis. A refined model of the aerobic metabolism of glycerol via the GldA/DhaK pathway is presented.

Novel listerial glycerol dehydrogenase- and phosphoenolpyruvate-dependent dihydroxyacetone kinase system connected to the pentose phosphate pathway

Several bacteria use glycerol dehydrogenase to transform glycerol into dihydroxyacetone (Dha). Dha is subsequently converted into Dha phosphate (Dha-P) by an ATP- or phosphoenolpyruvate (PEP)-dependent Dha kinase. Listeria innocua possesses two potential PEP-dependent Dha kinases. One is encoded by 3 of the 11 genes forming the glycerol (gol) operon. This operon also contains golD (lin0362), which codes for a new type of Dha-forming NAD(+)-dependent glycerol dehydrogenase. The subsequent metabolism of Dha requires its phosphorylation via the PEP:sugar phosphotransferase system components enzyme I, HPr, and EIIA(Dha)-2 (Lin0369). P∼EIIA(Dha)-2 transfers its phosphoryl group to DhaL-2, which phosphorylates Dha bound to DhaK-2. The resulting Dha-P is probably metabolized mainly via the pentose phosphate pathway, because two genes of the gol operon encode proteins resembling transketolases and transaldolases. In addition, purified Lin0363 and Lin0364 exhibit ribose-5-P isomerase (RipB) and triosephosphate isomerase activities, respectively. The latter enzyme converts part of the Dha-P into glyceraldehyde-3-P, which, together with Dha-P, is metabolized via gluconeogenesis to form fructose-6-P. Together with another glyceraldehyde-3-P molecule, the transketolase transforms fructose-6-P into intermediates of the pentose phosphate pathway. The gol operon is preceded by golR, transcribed in the opposite orientation and encoding a DeoR-type repressor. Its inactivation causes the constitutive but glucose-repressible expression of the entire gol operon, including the last gene, encoding a pediocin immunity-like (PedB-like) protein. Its elevated level of synthesis in the golR mutant causes slightly increased immunity against pediocin PA-1 compared to the wild-type strain or a pedB-like deletion mutant.

Molecular characterization of the glycerol-oxidative pathway of Clostridium butyricum VPI 1718

The glycerol oxidative pathway of Clostridium butyricum VPI 1718 plays an important role in glycerol dissimilation. We isolated, sequenced, and characterized the region coding for the glycerol oxidation pathway. Five open reading frames (ORFs) were identified: dhaR, encoding a putative transcriptional regulator; dhaD (1,142 bp), encoding a glycerol dehydrogenase; and dhaK (995 bp), dhaL (629 bp), and dhaM (386 bp), encoding a phosphoenolpyruvate (PEP)-dependent dihydroxyacetone (DHA) kinase enzyme complex. Northern blot analysis demonstrated that the last four genes are transcribed as a 3.2-kb polycistronic operon only in glycerol-metabolizing cultures, indicating that the expression of this operon is regulated at the transcriptional level. The transcriptional start site of the operon was determined by primer extension, and the promoter region was deduced. The glycerol dehydrogenase activity of DhaD and the PEP-dependent DHA kinase activity of DhaKLM were demonstrated by heterologous expression in different Escherichia coli mutants. Based on our complementation experiments, we proposed that the HPr phosphoryl carrier protein and His9 residue of the DhaM subunit are involved in the phosphoryl transfer to dihydroxyacetone-phosphate. DhaR, a potential regulator of this operon, was found to contain conserved transmitter and receiver domains that are characteristic of two-component systems present in the AraC family. To the best of our knowledge, this is the first molecular characterization of a glycerol oxidation pathway in a Gram-positive bacterium.

Structural and mechanistic insight into covalent substrate binding by Escherichia coli dihydroxyacetone kinase

The Escherichia coli dihydroxyacetone (Dha) kinase is an unusual kinase because (i) it uses the phosphoenolpyruvate carbohydrate: phosphotransferase system (PTS) as the source of high-energy phosphate, (ii) the active site is formed by two subunits, and (iii) the substrate is covalently bound to His218(K)* of the DhaK subunit. The PTS transfers phosphate to DhaM, which in turn phosphorylates the permanently bound ADP coenzyme of DhaL. This phosphoryl group is subsequently transferred to the Dha substrate bound to DhaK. Here we report the crystal structure of the E. coli Dha kinase complex, DhaK-DhaL. The structure of the complex reveals that DhaK undergoes significant conformational changes to accommodate binding of DhaL. Combined mutagenesis and enzymatic activity studies of kinase mutants allow us to propose a catalytic mechanism for covalent Dha binding, phosphorylation, and release of the Dha-phosphate product. Our results show that His56(K) is involved in formation of the covalent hemiaminal bond with Dha. The structure of H56N(K) with noncovalently bound substrate reveals a somewhat different positioning of Dha in the binding pocket as compared to covalently bound Dha, showing that the covalent attachment to His218(K) orients the substrate optimally for phosphoryl transfer. Asp109(K) is critical for activity, likely acting as a general base activating the γ-OH of Dha. Our results provide a comprehensive picture of the roles of the highly conserved active site residues of dihydroxyacetone kinases.

Crystal structure of enzyme I of the phosphoenolpyruvate sugar phosphotransferase system in the dephosphorylated state

The bacterial phosphoenolpyruvate (PEP) sugar phosphotransferase system mediates sugar uptake and controls the carbon metabolism in response to carbohydrate availability. Enzyme I (EI), the first component of the phosphotransferase system, consists of an N-terminal protein binding domain (EIN) and a C-terminal PEP binding domain (EIC). EI transfers phosphate from PEP by double displacement via a histidine residue on EIN to the general phosphoryl carrier protein HPr. Here we report the 2.4 A crystal structure of the homodimeric EI from Staphylococcus aureus. EIN consists of the helical hairpin HPr binding subdomain and the phosphorylatable betaalpha phospho-histidine (P-His) domain. EIC folds into an (betaalpha)(8) barrel. The dimer interface of EIC buries 1833 A(2) of accessible surface per monomer and contains two Ca(2+) binding sites per dimer. The structures of the S. aureus and Escherichia coli EI domains (Teplyakov, A., Lim, K., Zhu, P. P., Kapadia, G., Chen, C. C., Schwartz, J., Howard, A., Reddy, P. T., Peterkofsky, A., and Herzberg, O. (2006) Proc. Natl. Acad. Sci. U.S.A. 103, 16218-16223) are very similar. The orientation of the domains relative to each other, however, is different. In the present structure the P-His domain is docked to the HPr binding domain in an orientation appropriate for in-line transfer of the phosphate to the active site histidine of the acceptor HPr. In the E. coli structure the phospho-His of the P-His domain projects into the PEP binding site of EIC. In the S. aureus structure the crystallographic temperature factors are lower for the HPr binding domain in contact with the P-His domain and higher for EIC. In the E. coli structure it is the reverse.