Crystal structures of an ATP-dependent hexokinase with broad substrate specificity from the hyperthermophilic archaeon Sulfolobus tokodaii

The crystal structure and biochemical analyses of Escherichia coli YqgF illuminate its diverse functions

The Escherichia coli yqgF gene product is essential for bacterial growth and to confer resistance to multiple antimicrobial agents. Furthermore, additional evidence suggests that YqgF, a RuvC family protein, is required for DNA damage repair, yet the mechanism underlying its action remains elusive. To address this knowledge gap, we conducted structural and biochemical investigations on E. coli YqgF (EcYqgF). Here, we reveal that EcYqgF binds preferentially to branched DNA structures compared to single-stranded (ssDNA) and double-stranded DNA (dsDNA), and that the EcYqgF:DNA complexes formed with branched DNA species were more stable and resistant against high salt and excess of competitor DNA than those formed with dsDNA. We show that EcYqgF has a strong preference towards cleavage of branched DNA structures than dsDNA, ssDNA, and dsDNA with 5'- or 3'-ssDNA overhangs. Crucially, we found that EcYqgF has a DNA-independent, Mg2+-dependent ATPase activity that is tightly coupled to DNA cleavage. We have determined the crystal structure of EcYqgF, developed a model of ATP binding using AI-based methods, and rationalized the impact of site-directed mutants on ATP binding. Furthermore, we discovered two unusual ATPase-defective EcYqgF variants, proficient in ATP-binding but not hydrolysis, which display a modest increase in the DNA-binding affinity, yet are devoid of endonucleolytic activity, thus revealing a previously unappreciated property of YqgF endonucleases. Collectively, our results suggest that despite its overall structural similarity to the well-studied resolvase, RuvC, EcYqgF is functionally distinct. Importantly, the EcYqgF dual activity that couples ATP hydrolysis to endonuclease activity is absent in RuvC.

MolViewSpec: a Mol* extension for describing and sharing molecular visualizations

Data visualization is a pivotal component of a structural biologist's arsenal. The Mol* Viewer makes molecular visualizations available to broader audiences via most web browsers. While Mol* provides a wide range of functionality, it has a steep learning curve and is only available via a JavaScript interface. To enhance the accessibility and usability of web-based molecular visualization, we introduce MolViewSpec (molstar.org/mol-view-spec), a standardized approach for defining molecular visualizations that decouples the definition of complex molecular scenes from their rendering. Scene definition can include references to commonly used structural, volumetric, and annotation data formats together with a description of how the data should be visualized and paired with optional annotations specifying colors, labels, measurements, and custom 3D geometries. Developed as an open standard, this solution paves the way for broader interoperability and support across different programming languages and molecular viewers, enabling more streamlined, standardized, and reproducible visual molecular analyses. MolViewSpec is freely available as a Mol* extension and a standalone Python package.

Structural and biochemical insights into the molecular mechanism of N-acetylglucosamine/N-Acetylmuramic acid kinase MurK from Clostridium acetobutylicum

MurK is a MurNAc- and GlcNAc-specific amino sugar kinase, phosphorylates MurNAc and GlcNAc at the 6-hydroxyl group in an ATP-dependent manner, and contributes to the recovery of both amino sugars during the cell wall turnover in Clostridium acetobutylicum. Herein, we determined the crystal structures of MurK in complex with MurNAc, GlcNAc, and glucose, respectively. MurK represents the V-shaped fold, which is divided into a small N-terminal domain and a large C-terminal domain. The catalytic pocket is located within the deep cavity between the two domains of the MurK monomer. We mapped the significant enzyme-substrate interactions, identified key residues involved in the catalytic activity of MurK, and found that residues Asp77 and Arg78 from the β4-α2-loop confer structural flexibilities to specifically accommodate GlcNAc and MurNAc, respectively. Moreover, structural comparison revealed that MurK adopts closed-active conformation induced by the N-acetyl moiety from GlcNAc/MurNAc, rather than closed-inactive conformation induced by glucose, to carry out its catalytic reaction. Taken together, our study provides structural and functional insights into the molecular mechanism of MurK for the phosphorylation of both MurNAc and GlcNAc, sugar substrate specificity, and conformational changes upon sugar substrate binding.

Binding Affinity Determination in Drug Design: Insights from Lock and Key, Induced Fit, Conformational Selection, and Inhibitor Trapping Models

Binding affinity is a fundamental parameter in drug design, describing the strength of the interaction between a molecule and its target protein. Accurately predicting binding affinity is crucial for the rapid development of novel therapeutics, the prioritization of promising candidates, and the optimization of their properties through rational design strategies. Binding affinity is determined by the mechanism of recognition between proteins and ligands. Various models, including the lock and key, induced fit, and conformational selection, have been proposed to explain this recognition process. However, current computational strategies to predict binding affinity, which are based on these models, have yet to produce satisfactory results. This article explores the connection between binding affinity and these protein-ligand interaction models, highlighting that they offer an incomplete picture of the mechanism governing binding affinity. Specifically, current models primarily center on the binding of the ligand and do not address its dissociation. In this context, the concept of ligand trapping is introduced, which models the mechanisms of dissociation. When combined with the current models, this concept can provide a unified theoretical framework that may allow for the accurate determination of the ligands' binding affinity.

Identifying protein conformational states in the Protein Data Bank: Toward unlocking the potential of integrative dynamics studies

Studying protein dynamics and conformational heterogeneity is crucial for understanding biomolecular systems and treating disease. Despite the deposition of over 215 000 macromolecular structures in the Protein Data Bank and the advent of AI-based structure prediction tools such as AlphaFold2, RoseTTAFold, and ESMFold, static representations are typically produced, which fail to fully capture macromolecular motion. Here, we discuss the importance of integrating experimental structures with computational clustering to explore the conformational landscapes that manifest protein function. We describe the method developed by the Protein Data Bank in Europe - Knowledge Base to identify distinct conformational states, demonstrate the resource's primary use cases, through examples, and discuss the need for further efforts to annotate protein conformations with functional information. Such initiatives will be crucial in unlocking the potential of protein dynamics data, expediting drug discovery research, and deepening our understanding of macromolecular mechanisms.

Computational insights into dynamics and conformational stability of N-acetylmannosamine kinase mutations

The activity of UDP-GlcNAc 2-epimerase/ManNAc kinase (GNE) is essential for the biosynthesis of sialic acid, which is involved in cellular processes in health and diseases. GNE contains an N-terminal epimerase domain and a C-terminal kinase domain (N-acetylmannosamine kinase, MNK). Mutations of the GNE protein led to hypoactivity of the enzyme and cause sialurea or autosomal recessive inclusion body myopathy/Nonaka myopathy. Here, we used all-atom molecular dynamics (MD) simulations to comprehend the folding, dynamics and conformational stability of MNK variants, including the wild type (WT) and three mutants (H677R, V696M and H677R/V696M). The deleterious and destabilizing nature of MNK mutants were predicted using different prediction tools. Results predicted that mutations modulate the stability, flexibility and function of MNK. The effect of mutations on the conformational stability and dynamics of MNK was next studied through the free-energy landscape (FEL), hydrogen-bonds and secondary structure changes. The FEL results show that the mutations interfere with various conformational transitions in both WT and mutants, exposing the structural underpinnings of protein destabilization and unfolding brought on by mutation. We discover that, when compared to the other two mutations, V696M and H677R/V696M, H677R has the most harmful effects. These findings have a strong correlation with published experimental studies that demonstrate how these mutations disrupt MNK activity. Hence, this computational study describes the structural details to unravel the mutant effects at the atomistic resolution and has implications for understanding the GNE's physiological and pathological role.

N-acetylmuramic acid recognition by MurK kinase from the MurNAc auxotrophic oral pathogen Tannerella forsythia

The bacterial cell wall consists of a three-dimensional peptidoglycan layer, composed of peptides linked to the sugars N-acetylmuramic acid (MurNAc) and GlcNAc. Unlike other bacteria, the pathogenic Tannerella forsythia, a member of the red complex group of bacteria associated with the late stages of periodontitis, lacks biosynthetic pathways for MurNAc production and therefore obtains MurNAc from the environment. Sugar kinases play a crucial role in the MurNAc recycling process, activating the sugar molecules by phosphorylation. In this study, we present the first crystal structures of a MurNAc kinase, called murein sugar kinase (MurK), in its unbound state as well as in complexes with the ATP analog β-γ-methylene adenosine triphosphate (AMP-PCP) and with MurNAc. We also determined the crystal structures of K1058, a paralogous MurNAc kinase of T. forsythia, in its unbound state and in complex with MurNAc. We identified the active site and residues crucial for MurNAc specificity as the less bulky side chains of S133, P134, and L135, which enlarge the binding cavity for the lactyl ether group, unlike the glutamate or histidine residues present in structural homologs. In establishing the apparent kinetic parameters for both enzymes, we showed a comparable affinity for MurNAc (Km 180 μM and 30 μM for MurK and K1058, respectively), with MurK being over two hundred times faster than K1058 (Vmax 80 and 0.34 μmol min-1 mg-1, respectively). These data might support a structure-guided approach to development of inhibitory MurNAc analogs for pathogen MurK enzymes.

Structural Analysis of Plasmodium falciparum Hexokinase Provides Novel Information about Catalysis Due to a Plasmodium-Specific Insertion

The protozoan parasite Plasmodium falciparum is the causative pathogen of the most severe form of malaria, for which novel strategies for treatment are urgently required. The primary energy supply for intraerythrocytic stages of Plasmodium is the production of ATP via glycolysis. Due to the parasite's strong dependence on this pathway and the significant structural differences of its glycolytic enzymes compared to its human counterpart, glycolysis is considered a potential drug target. In this study, we provide the first three-dimensional protein structure of P. falciparum hexokinase (PfHK) containing novel information about the mechanisms of PfHK. We identified for the first time a Plasmodium-specific insertion that lines the active site. Moreover, we propose that this insertion plays a role in ATP binding. Residues of the insertion further seem to affect the tetrameric interface and therefore suggest a special way of communication among the different monomers. In addition, we confirmed that PfHK is targeted and affected by oxidative posttranslational modifications (oxPTMs). Both S-glutathionylation and S-nitrosation revealed an inhibitory effect on the enzymatic activity of PfHK.

The ROK kinase N-acetylglucosamine kinase uses a sequential random enzyme mechanism with successive conformational changes upon each substrate binding

N-acetyl-d-glucosamine (GlcNAc) is a major component of bacterial cell walls. Many organisms recycle GlcNAc from the cell wall or metabolize environmental GlcNAc. The first step in GlcNAc metabolism is phosphorylation to GlcNAc-6-phosphate. In bacteria, the ROK family kinase N-acetylglucosamine kinase (NagK) performs this activity. Although ROK kinases have been studied extensively, no ternary complex showing the two substrates has yet been observed. Here, we solved the structure of NagK from the human pathogen Plesiomonas shigelloides in complex with GlcNAc and the ATP analog AMP-PNP. Surprisingly, PsNagK showed distinct conformational changes associated with the binding of each substrate. Consistent with this, the enzyme showed a sequential random enzyme mechanism. This indicates that the enzyme acts as a coordinated unit responding to each interaction. Our molecular dynamics modeling of catalytic ion binding confirmed the location of the essential catalytic metal. Additionally, site-directed mutagenesis confirmed the catalytic base and that the metal-coordinating residue is essential. Together, this study provides the most comprehensive insight into the activity of a ROK kinase.

Glycation Interferes with the Activity of the Bi-Functional UDP-N-Acetylglucosamine 2-Epimerase/N-Acetyl-mannosamine Kinase (GNE)

Mutations in the gene coding for the bi-functional UDP-N-acetylglucosamine 2-epimerase/N-acetylmannosamine kinase (GNE), the key enzyme of the sialic acid biosynthesis, are responsible for autosomal-recessive GNE myopathy (GNEM). GNEM is an adult-onset disease with a yet unknown exact pathophysiology. Since the protein appears to work adequately for a certain period of time even though the mutation is already present, other effects appear to influence the onset and progression of the disease. In this study, we want to investigate whether the late onset of GNEM is based on an age-related effect, e.g., the accumulation of post-translational modifications (PTMs). Furthermore, we also want to investigate what effect on the enzyme activity such an accumulation would have. We will particularly focus on glycation, which is a PTM through non-enzymatic reactions between the carbonyl groups (e.g., of methylglyoxal (MGO) or glyoxal (GO)) with amino groups of proteins or other biomolecules. It is already known that the levels of both MGO and GO increase with age. For our investigations, we express each domain of the GNE separately, treat them with one of the glycation agents, and determine their activity. We demonstrate that the enzymatic activity of the N-acetylmannosamine kinase (GNE-kinase domain) decreases dramatically after glycation with MGO or GO-with a remaining activity of 13% ± 5% (5 mM MGO) and 22% ± 4% (5 mM GO). Whereas the activity of the UDP-N-acetylglucosamine 2-epimerase (GNE-epimerase domain) is only slightly reduced after glycation-with a remaining activity of 60% ± 8% (5 mM MGO) and 63% ± 5% (5 mM GO).

Functional analyses of a highly thermostable hexokinase from Pyrobaculum calidifontis

The gene encoding a repressor open reading frame sugar kinase (ROK) family protein from hyperthermophilic crenarchaeon Pyrobaculum calidifontis, Pcal-HK, was cloned and expressed in Escherichia coli. The recombinant protein was produced in soluble and highly active form. Purified Pcal-HK was highly thermostable and existed in a monomeric form in solution. The enzyme was specific to ATP as phosphoryl donor but showed broad specificity to phosphoryl acceptors. It catalyzed the phosphorylation of a number of hexoses, including glucose, glucosamine, N-acetyl glucosamine, fructose and mannose, at nearly the same rate and similar affinity. The enzyme was metal ion dependent exhibiting highest activity at 90-95 °C and pH 8.5. Mg²⁺ was most effective metal ion, which could be partially replaced by Mn²⁺, Ni²⁺ or Zn²⁺. Kinetic parameters were determined at 90 °C and the enzyme showed almost similar catalytic efficiency (k_cat/K_m) towards the above mentioned hexoses. To the best of our knowledge, Pcal-HK is the most active thermostable ROK family hexokinase characterized to date which catalyzes the phosphorylation of various hexoses with nearly similar affinity.

Novel mutation in hexokinase 2 confers resistance to 2-deoxyglucose by altering protein dynamics

Glucose is central to many biological processes, serving as an energy source and a building block for biosynthesis. After glucose enters the cell, hexokinases convert it to glucose-6-phosphate (Glc-6P) for use in anaerobic fermentation, aerobic oxidative phosphorylation, and the pentose-phosphate pathway. We here describe a genetic screen in Saccharomyces cerevisiae that generated a novel spontaneous mutation in hexokinase-2, hxk2^G238V, that confers resistance to the toxic glucose analog 2-deoxyglucose (2DG). Wild-type hexokinases convert 2DG to 2-deoxyglucose-6-phosphate (2DG-6P), but 2DG-6P cannot support downstream glycolysis, resulting in a cellular starvation-like response. Curiously, though the hxk2^G238V mutation encodes a loss-of-function allele, the affected amino acid does not interact directly with bound glucose, 2DG, or ATP. Molecular dynamics simulations suggest that Hxk2^G238V impedes sugar binding by altering the protein dynamics of the glucose-binding cleft, as well as the large-scale domain-closure motions required for catalysis. These findings shed new light on Hxk2 dynamics and highlight how allosteric changes can influence catalysis, providing new structural insights into this critical regulator of carbohydrate metabolism. Given that hexokinases are upregulated in some cancers and that 2DG and its derivatives have been studied in anti-cancer trials, the present work also provides insights that may apply to cancer biology and drug resistance.

Glucose fuels many of the energy-production processes required for normal cell growth. Before glucose can participate in these processes, it must first be chemically modified by proteins called hexokinases. To better understand how hexokinases modify glucose—and how mutations in hexokinase genes might confer drug resistance—we evolved resistance in yeast to a toxic hexokinase-binding molecule called 2DG. We discovered a mutation in the hexokinase gene that confers 2DG resistance and reduces the protein’s ability to modify glucose. Biochemical analyses and computer simulations of the hexokinase protein suggest that the mutation diminishes glucose binding by altering enzyme flexibility. This work shows how cells can evolve resistance to toxins via only modest changes to protein structures. Furthermore, because cancer-cell hexokinases are particularly active, 2DG has been studied as cancer chemotherapy. Thus, the insights this work provides might also apply to cancer biology.

Molecular mechanism of selective substrate engagement and inhibitor disengagement of cysteine synthase

O-acetyl serine sulfhydrylase (OASS), referred to as cysteine synthase (CS), synthesizes cysteine from O-acetyl serine (OAS) and sulfur in bacteria and plants. The inherent challenge for CS is to overcome 4 to 6 log-folds stronger affinity for its natural inhibitor, serine acetyltransferase (SAT), as compared with its affinity for substrate, OAS. Our recent study showed that CS employs a novel competitive-allosteric mechanism to selectively recruit its substrate in the presence of natural inhibitor. In this study, we trace the molecular features that control selective substrate recruitment. To generalize our findings, we used CS from three different bacteria (Haemophilus, Salmonella, and Mycobacterium) as our model systems and analyzed structural and substrate-binding features of wild-type CS and its ∼13 mutants. Results show that CS uses a noncatalytic residue, M120, located 20 Å away from the reaction center, to discriminate in favor of substrate. M120A and background mutants display significantly reduced substrate binding, catalytic efficiency, and inhibitor binding. Results shows that M120 favors the substrate binding by selectively enhancing the affinity for the substrate and disengaging the inhibitor by 20 to 286 and 5- to 3-folds, respectively. Together, M120 confers a net discriminative force in favor of substrate by 100- to 858-folds.

An intrinsic FRET sensor of protein-ligand interactions

We describe an approach for the development of fluorescent sensors of metabolite binding in which a genetically encoded fluorescent non-canonical amino acid (fNCAA) containing a 7-hydroxycoumarin moiety (7-HCAA) forms a FRET pair with native tryptophan residues. Although previous studies demonstrated the potential for using 7-HCAA as an acceptor for tryptophan, this approach has not yet been explored within a single protein containing multiple tryptophan residues. A structure-based analysis of a hexokinase enzyme with multiple native tryptophan residues identified glutamate 50 as a potential site of 7-HCAA incorporation; Glu50 moves closer to the native tryptophans upon substrate binding. Substitution of 7-HCAA at residue 50 led to an increase in FRET efficiency in the presence of the substrate; this effect was not observed in a control protein where no change in distance between 7-HCAA and the native tryptophans occurs on substrate binding. This system was then used to directly observe differences in binding affinity of the hexokinase that occur at a number of pH values. Our approach builds on previous research in that it eliminates the need for the incorporation of multiple fNCAAs or fluorescent labels within a target protein and can be used to study substrate binding with native ligands. As such, it serves to expand the versatility of FRET-based techniques.

Plasmodium vivax and human hexokinases share similar active sites but display distinct quaternary architectures

Malaria is a devastating disease caused by a protozoan parasite. It affects over 300 million individuals and results in over 400 000 deaths annually, most of whom are young children under the age of five. Hexokinase, the first enzyme in glucose metabolism, plays an important role in the infection process and represents a promising target for therapeutic intervention. Here, cryo-EM structures of two conformational states of Plasmodium vivax hexokinase (PvHK) are reported at resolutions of ∼3 Å. It is shown that unlike other known hexokinase structures, PvHK displays a unique tetrameric organization (∼220 kDa) that can exist in either open or closed quaternary conformational states. Despite the resemblance of the active site of PvHK to its mammalian counterparts, this tetrameric organization is distinct from that of human hexokinases, providing a foundation for the structure-guided design of parasite-selective antimalarial drugs.

Crystal structures of rice hexokinase 6 with a series of substrates shed light on its enzymatic mechanism

Hexokinases (HXKs) have determined to be multifaceted proteins, and they are the only ones able to phosphorylate glucose in plants. However, the binding mode for ATP to plant HXKs remains unclear. Here, we report the crystal structures of rice hexokinase 6 (OsHXK6) in four different forms: (i) apo-form, (ii) binary complex with D-Glc, (iii) quaternary complex with ADP, PO₄ and Mg²⁺, and (iv) pentanary complex with D-Glc, ADP, PO₄, and Mg²⁺. The apo form is in the open state conformation, and the three others are in the closed state, indicating that glucose and ADP-PO₄ binding induces a large conformational change by domain rearrangement. The quaternary complex is a novel intermediate during the catalytic reaction we trapped for the first time, which provides a new evidence for the enzymatic mechanism of HXKs. In addition, the latter two complexes reveal the binding mode for ADP-PO₄ to plant HXKs, which provide the structural explanation for the dual-function of OsHXK6. In addition, we identified that residues Gly112, Thr261, Gly262, and Gly450 are essential to the binding between ADP-PO₄ and OsHXK6 by a series of single mutations and enzymatic assays. Our study provide structural basis for the other functional studies of OsHXK6 in rice.

Carbohydrate Kinases: A Conserved Mechanism Across Differing Folds

Carbohydrate kinases activate a wide variety of monosaccharides by adding a phosphate group, usually from ATP. This modification is fundamental to saccharide utilization, and it is likely a very ancient reaction. Modern organisms contain carbohydrate kinases from at least five main protein families. These range from the highly specialized inositol kinases, to the ribokinases and galactokinases, which belong to families that phosphorylate a wide range of substrates. The carbohydrate kinases utilize a common strategy to drive the reaction between the sugar hydroxyl and the donor phosphate. Each sugar is held in position by a network of hydrogen bonds to the non-reactive hydroxyls (and other functional groups). The reactive hydroxyl is deprotonated, usually by an aspartic acid side chain acting as a catalytic base. The deprotonated hydroxyl then attacks the donor phosphate. The resulting pentacoordinate transition state is stabilized by an adjacent divalent cation, and sometimes by a positively charged protein side chain or the presence of an anion hole. Many carbohydrate kinases are allosterically regulated using a wide variety of strategies, due to their roles at critical control points in carbohydrate metabolism. The evolution of a similar mechanism in several folds highlights the elegance and simplicity of the catalytic scheme.

Producing glucose 6-phosphate from cellulosic biomass: structural insights into levoglucosan bioconversion

The most abundant carbohydrate product of cellulosic biomass pyrolysis is the anhydrosugar levoglucosan (1,6-anhydro-β-d-glucopyranose), which can be converted to glucose 6-phosphate by levoglucosan kinase (LGK). In addition to the canonical kinase phosphotransfer reaction, the conversion requires cleavage of the 1,6-anhydro ring to allow ATP-dependent phosphorylation of the sugar O6 atom. Using x-ray crystallography, we show that LGK binds two magnesium ions in the active site that are additionally coordinated with the nucleotide and water molecules to result in ideal octahedral coordination. To further verify the metal binding sites, we co-crystallized LGK in the presence of manganese instead of magnesium and solved the structure de novo using the anomalous signal from four manganese atoms in the dimeric structure. The first metal is required for catalysis, whereas our work suggests that the second is either required or significantly promotes the catalytic rate. Although the enzyme binds its sugar substrate in a similar orientation to the structurally related 1,6-anhydro-N-acetylmuramic acid kinase (AnmK), it forms markedly fewer bonding interactions with the substrate. In this orientation, the sugar is in an optimal position to couple phosphorylation with ring cleavage. We also observed a second alternate binding orientation for levoglucosan, and in these structures, ADP was found to bind with lower affinity. These combined observations provide an explanation for the high Km of LGK for levoglucosan. Greater knowledge of the factors that contribute to the catalytic efficiency of LGK can be used to improve applications of this enzyme for levoglucosan-derived biofuel production.

Transcriptome analysis of mRNA and miRNA in skeletal muscle indicates an important network for differential Residual Feed Intake in pigs

Feed efficiency (FE) can be measured by feed conversion ratio (FCR) or residual feed intake (RFI). In this study, we measured the FE related phenotypes of 236 castrated purebred Yorkshire boars, and selected 10 extreme individuals with high and low RFI for transcriptome analysis. We used RNA-seq analyses to determine the differential expression of genes and miRNAs in skeletal muscle. There were 99 differentially expressed genes identified (q ≤ 0.05). The down-regulated genes were mainly involved in mitochondrial energy metabolism, including FABP3, RCAN, PPARGC1 (PGC-1A), HK2 and PRKAG2. The up-regulated genes were mainly involved in skeletal muscle differentiation and proliferation, including IGF2, PDE7A, CEBPD, PIK3R1 and MYH6. Moreover, 15 differentially expressed miRNAs (|log2FC| ≥ 1, total reads count ≥ 20, p ≤ 0.05) were identified. Among them, miR-136, miR-30e-5p, miR-1, miR-208b, miR-199a, miR-101 and miR-29c were up-regulated, while miR-215, miR-365-5p, miR-486, miR-1271, miR-145, miR-99b, miR-191 and miR-10b were down-regulated in low RFI pigs. We conclude that decreasing mitochondrial energy metabolism, possibly through AMPK - PGC-1A pathways, and increasing muscle growth, through IGF-1/2 and TGF-β signaling pathways, are potential strategies for the improvement of FE in pigs (and possibly other livestock). This study provides new insights into the molecular mechanisms that determine RFI and FE in pigs.

Open-close structural change upon ligand binding and two magnesium ions required for the catalysis of N-acetylhexosamine 1-kinase

Infant gut-associated bifidobacteria possess a metabolic pathway to utilize lacto-N-biose (Gal-β1,3-GlcNAc) and galacto-N-biose (Gal-β1,3-GalNAc) from human milk and glycoconjugates specifically. In this pathway, N-acetylhexosamine 1-kinase (NahK) catalyzes the phosphorylation of GlcNAc or GalNAc at the anomeric C1 position with ATP. Crystal structures of NahK have only been determined in the closed state. In this study, we determined open state structures of NahK in three different forms (apo, ADP complex, and ATP complex). A comparison of the open and closed state structures revealed an induced fit structural change defined by two rigid domains. ATP binds to the small N-terminal domain, and binding of the N-acetylhexosamine substrate to the large C-terminal domain induces a closing conformational change with a rotation angle of 16°. In the nucleotide binding site, two magnesium ions bridging the α-γ and β-γ phosphates were identified. A mutational analysis indicated that a residue coordinating both of the two magnesium ions (Asp228) is essential for catalysis. The involvement of two magnesium ions in the catalytic machinery is structurally similar to the catalytic structures of protein kinases and aminoglycoside phosphotransferases, but distinct from the structures of other anomeric kinases or sugar 6-kinases. These findings help to elucidate the possible evolutionary adaptation of substrate specificities and induced fit mechanism.

Theoretical identification on the role of Lys15 for Sulfolobus tokodaii hexokinase

QM/MM studies indicate that the phosphate transfer process is rate-determining, while the existence of K15 facilitates the reaction to proceed.

Interplay of Mg2+, ADP, and ATP in the cytosol and mitochondria: unravelling the role of Mg2+ in cell respiration

In animal and plant cells, the ATP/ADP ratio and/or energy charge are generally considered key parameters regulating metabolism and respiration. The major alternative issue of whether the cytosolic and mitochondrial concentrations of ADP and ATP directly mediate cell respiration remains unclear, however. In addition, because only free nucleotides are exchanged by the mitochondrial ADP/ATP carrier, whereas MgADP is the substrate of ATP synthase (EC 3.6.3.14), the cytosolic and mitochondrial Mg(2+) concentrations must be considered as well. Here we developed in vivo/in vitro techniques using (31)P-NMR spectroscopy to simultaneously measure these key components in subcellular compartments. We show that heterotrophic sycamore (Acer pseudoplatanus L.) cells incubated in various nutrient media contain low, stable cytosolic ADP and Mg(2+) concentrations, unlike ATP. ADP is mainly free in the cytosol, but complexed by Mg(2+) in the mitochondrial matrix, where [Mg(2+)] is tenfold higher. In contrast, owing to a much higher affinity for Mg(2+), ATP is mostly complexed by Mg(2+) in both compartments. Mg(2+) starvation used to alter cytosolic and mitochondrial [Mg(2+)] reversibly increases free nucleotide concentration in the cytosol and matrix, enhances ADP at the expense of ATP, decreases coupled respiration, and stops cell growth. We conclude that the cytosolic ADP concentration, and not ATP, ATP/ADP ratio, or energy charge, controls the respiration of plant cells. The Mg(2+) concentration, remarkably constant and low in the cytosol and tenfold higher in the matrix, mediates ADP/ATP exchange between the cytosol and matrix, [MgADP]-dependent mitochondrial ATP synthase activity, and cytosolic free ADP homeostasis.

Differential effects of a mutation on the normal and promiscuous activities of orthologs: implications for natural and directed evolution

Neutral drift occurring over millions or billions of years results in substantial sequence divergence among enzymes that catalyze the same reaction. Although natural selection maintains the primary activity of orthologous enzymes, there is, by definition, no selective pressure to maintain physiologically irrelevant promiscuous activities. Thus, the levels and the evolvabilities of promiscuous activities may vary among orthologous enzymes. Consistent with this expectation, we have found that the levels of a promiscuous activity in nine gamma-glutamyl phosphate reductase (ProA) orthologs vary by about 50-fold. Remarkably, a single amino acid change from Glu to Ala near the active site appeared to be critical for improvement of the promiscuous activity in every ortholog. The effects of this change varied dramatically. The improvement in the promiscuous activity varied from 50- to 770-fold, and, importantly, was not correlated with the initial level of the promiscuous activity. The decrease in the original activity varied from 190- to 2,100-fold. These results suggest that evolution of a novel enzyme may be possible in some microbes, but not in others. Further, these results underscore the importance of using multiple orthologs as starting points for directed evolution of novel enzyme activities.

Carbohydrate metabolism in Archaea: current insights into unusual enzymes and pathways and their regulation

The metabolism of Archaea, the third domain of life, resembles in its complexity those of Bacteria and lower Eukarya. However, this metabolic complexity in Archaea is accompanied by the absence of many "classical" pathways, particularly in central carbohydrate metabolism. Instead, Archaea are characterized by the presence of unique, modified variants of classical pathways such as the Embden-Meyerhof-Parnas (EMP) pathway and the Entner-Doudoroff (ED) pathway. The pentose phosphate pathway is only partly present (if at all), and pentose degradation also significantly differs from that known for bacterial model organisms. These modifications are accompanied by the invention of "new," unusual enzymes which cause fundamental consequences for the underlying regulatory principles, and classical allosteric regulation sites well established in Bacteria and Eukarya are lost. The aim of this review is to present the current understanding of central carbohydrate metabolic pathways and their regulation in Archaea. In order to give an overview of their complexity, pathway modifications are discussed with respect to unusual archaeal biocatalysts, their structural and mechanistic characteristics, and their regulatory properties in comparison to their classic counterparts from Bacteria and Eukarya. Furthermore, an overview focusing on hexose metabolic, i.e., glycolytic as well as gluconeogenic, pathways identified in archaeal model organisms is given. Their energy gain is discussed, and new insights into different levels of regulation that have been observed so far, including the transcript and protein levels (e.g., gene regulation, known transcription regulators, and posttranslational modification via reversible protein phosphorylation), are presented.

Conformational itinerary of Pseudomonas aeruginosa 1,6-anhydro-N-acetylmuramic acid kinase during its catalytic cycle

Anhydro-sugar kinases are unique from other sugar kinases in that they must cleave the 1,6-anhydro ring of their sugar substrate to phosphorylate it using ATP. Here we show that the peptidoglycan recycling enzyme 1,6-anhydro-N-acetylmuramic acid kinase (AnmK) from Pseudomonas aeruginosa undergoes large conformational changes during its catalytic cycle, with its two domains rotating apart by up to 32° around two hinge regions to expose an active site cleft into which the substrates 1,6-anhydroMurNAc and ATP can bind. X-ray structures of the open state bound to a nonhydrolyzable ATP analog (AMPPCP) and 1,6-anhydroMurNAc provide detailed insight into a ternary complex that forms preceding an operative Michaelis complex. Structural analysis of the hinge regions demonstrates a role for nucleotide binding and possible cross-talk between the bound ligands to modulate the opening and closing of AnmK. Although AnmK was found to exhibit similar binding affinities for ATP, ADP, and AMPPCP according to fluorescence spectroscopy, small angle x-ray scattering analyses revealed that AnmK adopts an open conformation in solution in the absence of ligand and that it remains in this open state after binding AMPPCP, as we had observed for our crystal structure of this complex. In contrast, the enzyme favored a closed conformation when bound to ADP in solution, consistent with a previous crystal structure of this complex. Together, our findings show that the open conformation of AnmK facilitates binding of both the sugar and nucleotide substrates and that large structural rearrangements must occur upon closure of the enzyme to correctly align the substrates and residues of the enzyme for catalysis.

Structural characterization of the ribonuclease H-like type ASKHA superfamily kinase MK0840 from Methanopyrus kandleri

Murein recycling is a process in which microorganisms recover peptidoglycan-degradation products in order to utilize them in cell wall biosynthesis or basic metabolic pathways. Methanogens such as Methanopyrus kandleri contain pseudomurein, which differs from bacterial murein in its composition and branching. Here, four crystal structures of the putative sugar kinase MK0840 from M. kandleri in apo and nucleotide-bound states are reported. MK0840 shows high similarity to bacterial anhydro-N-acetylmuramic acid kinase, which is involved in murein recycling. The structure shares a common fold with panthothenate kinase and the 2-hydroxyglutaryl-CoA dehydratase component A, both of which are members of the ASKHA (acetate and sugar kinases/Hsc70/actin) superfamily of phosphotransferases. Local conformational changes in the nucleotide-binding site between the apo and holo forms are observed upon nucleotide binding. Further insight is given into domain movements and putative active-site residues are identified.

ViewMotions Rainbow: a new method to illustrate molecular motions in proteins

The biological functions of many enzymes are often coupled with significant conformational changes. The end states of these conformational changes can often be determined by X-ray crystallography. These X-ray structures are snapshots of the two extreme conformations in which the macromolecule exists, but the dynamic movements between the states are not easily visualized in a two-dimensional illustration. Here we have developed a new method to visualize macromolecular motions called a ViewMotions Rainbow diagram. These diagrams show the initial and final states overlaid along with approximately 30 intermediate structures calculated by linear interpolation of the backbone coordinates of the initial and final states. This group of structures is then spectrally colored from the initial structure in blue to the final structure in red. ViewMotions Rainbow diagrams provide the reader with a much easier way to understand the macromolecular motions using a single two-dimensional illustration. Since producing these diagrams requires a number of different software packages, we have setup the ViewMotions Web Server (http://viewmotions.bc.edu) to automatically generate these diagrams from two Protein Data Bank files or from the Database of Macromolecular Movements (http://molmovdb.org).

Unusual arginine formations in protein function and assembly: rings, strings, and stacks

Protein-protein interfaces are often stabilized by a small number of dominant contacts, exemplified by the overrepresentation of arginine residues at oligomerization interfaces. Positively charged arginines are most commonly involved in ion pairs of opposite charge; however, previous work of Scheraga and co-workers described the stable, close range interaction between guanidinium pairs in a solvated environment. To extend this work, we searched over 70 thousand protein structures and complexes for unusual formations of arginine residues supported by the electron density. Symmetry transformations were used to generate full assemblies. Clusters of four to eight arginine residues with C(ζ)-C(ζ) distances <5 Å, organized as rings with four to eight members, stacks of two arginines, and strings of stacked arginines, are commonly located at the interfaces of oligomeric proteins. The positive charge is properly balanced by negatively charged counterions in about 90% of the cases. We also observed planar stacking of guanidinium groups, bridged by hydrogen bonds and interactions with water molecules. The guanidinium groups are commonly involved in five hydrogen bonds with water molecules and acceptor groups from surrounding amino acids. Water molecules have a bridging effect on the arginine pairs, but in some cases, small molecular weight chemicals in the crystallization buffer may be misinterpreted as water molecules. In summary, despite electrostatic repulsion, arginines do form various clusters that are exposed to interact with and potentially be controlled or switched by charged metabolites, membrane lipids, nucleic acids, or side chains of other proteins. Control of the stability of arginine clusters may play an important role in protein-protein oligomerization, molecular recognition, and ligand binding.

Crystal structures of N-acetylmannosamine kinase provide insights into enzyme activity and inhibition

Sialic acids are essential components of membrane glycoconjugates. They are responsible for the interaction, structure, and functionality of all deuterostome cells and have major functions in cellular processes in health and diseases. The key enzyme of the biosynthesis of sialic acid is the bifunctional UDP-N-acetylglucosamine-2-epimerase/N-acetylmannosamine kinase that transforms UDP-N-acetylglucosamine to N-acetylmannosamine (ManNAc) followed by its phosphorylation to ManNAc 6-phosphate and has a direct impact on the sialylation of cell surface components. Here, we present the crystal structures of the human N-acetylmannosamine kinase (MNK) domain of UDP-N-acetylglucosamine-2-epimerase/N-acetylmannosamine kinase in complexes with ManNAc at 1.64 Å resolution, MNK·ManNAc·ADP (1.82 Å) and MNK·ManNAc 6-phosphate · ADP (2.10 Å). Our findings offer detailed insights in the active center of MNK and serve as a structural basis to design inhibitors. We synthesized a novel inhibitor, 6-O-acetyl-ManNAc, which is more potent than those previously tested. Specific inhibitors of sialic acid biosynthesis may serve to further study biological functions of sialic acid.

Identification, tissue distribution, and molecular modeling of novel human isoforms of the key enzyme in sialic acid synthesis, UDP-GlcNAc 2-epimerase/ManNAc kinase

UDP-GlcNAc 2-epimerase/ManNAc kinase (GNE) catalyzes the first two committed steps in sialic acid synthesis. In addition to the three previously described human GNE isoforms (hGNE1-hGNE3), our database and polymerase chain reaction analysis yielded five additional human isoforms (hGNE4-hGNE8). hGNE1 is the ubiquitously expressed major isoform, while the hGNE2-hGNE8 isoforms are differentially expressed and may act as tissue-specific regulators of sialylation. hGNE2 and hGNE7 display a 31-residue N-terminal extension compared to hGNE1. On the basis of similarities to kinases and helicases, this extension does not seem to hinder the epimerase enzymatic active site. hGNE3 and hGNE8 contain a 55-residue N-terminal deletion and a 50-residue N-terminal extension compared to hGNE1. The size and secondary structures of these fragments are similar, and modeling predicted that these modifications do not affect the overall fold compared to that of hGNE1. However, the epimerase enzymatic activity of GNE3 and GNE8 is likely absent, because the deleted fragment contains important substrate binding residues in homologous bacterial epimerases. hGNE5-hGNE8 have a 53-residue deletion, which was assigned a role in substrate (UDP-GlcNAc) binding. Deletion of this fragment likely eliminates epimerase enzymatic activity. Our findings imply that GNE is subject to evolutionary mechanisms to improve cellular functions, without increasing the number of genes. Our expression and modeling data contribute to elucidation of the complex functional and regulatory mechanisms of human GNE and may contribute to further elucidating the pathology and treatment strategies of the human GNE-opathies sialuria and hereditary inclusion body myopathy.

Structure of 2-oxo-3-deoxygalactonate kinase from Klebsiella pneumoniae

In most organisms, efficient D-galactose utilization requires the highly conserved Leloir pathway that converts D-galactose to D-glucose 1-phosphate. However, in some bacterial and fungal species alternative routes of D-galactose assimilation have been identified. In the so-called De Ley-Doudoroff pathway, D-galactose is metabolized into pyruvate and D-glyceraldehyde 3-phosphate in five consecutive reactions carried out by specific enzymes. The penultimate step in this pathway involves the phosphorylation of 2-oxo-3-deoxygalactonate to 2-oxo-3-deoxygalactonate 6-phosphate catalyzed by 2-oxo-3-deoxygalactonate kinase, with ATP serving as a phosphoryl-group donor. Here, a crystal structure of 2-oxo-3-deoxygalactonate kinase from Klebsiella pneumoniae determined at 2.1 Å resolution is reported, the first structure of an enzyme from the De Ley-Doudoroff pathway. Structural comparison indicates that the enzyme belongs to the ASKHA (acetate and sugar kinases/hsc70/actin) family of phosphotransferases. The protein is composed of two α/β domains, each of which contains a core common to all family members. Additional elements introduced between conserved structural motifs define the unique features of 2-oxo-3-deoxygalactonate kinase and possibly determine the biological function of the protein.

Characterization of an N-acetylmuramic acid/N-acetylglucosamine kinase of Clostridium acetobutylicum

We report here the cloning and characterization of a cytoplasmic kinase of Clostridium acetobutylicum, named MurK (for murein sugar kinase). The enzyme has a unique specificity for both amino sugars of the bacterial cell wall, N-acetylmuramic acid (MurNAc) and N-acetylglucosamine (GlcNAc), which are phosphorylated at the 6-hydroxyl group. Kinetic analyses revealed Km values of 190 and 127 μM for MurNAc and GlcNAc, respectively, and a kcat value (65.0 s(-1)) that was 1.5-fold higher for the latter substrate. Neither the non-N-acetylated forms of the cell wall sugars, i.e., glucosamine and/or muramic acid, nor epimeric hexoses or 1,6-anhydro-MurNAc were substrates for the enzyme. MurK displays low overall amino acid sequence identity (24%) with human GlcNAc kinase and is the first characterized bacterial representative of the BcrAD/BadFG-like ATPase family. We propose a role of MurK in the recovery of muropeptides during cell wall rescue in C. acetobutylicum. The kinase was applied for high-sensitive detection of the amino sugars in cell wall preparations by radioactive phosphorylation.

Molecular basis of 1,6-anhydro bond cleavage and phosphoryl transfer by Pseudomonas aeruginosa 1,6-anhydro-N-acetylmuramic acid kinase

Anhydro-N-acetylmuramic acid kinase (AnmK) catalyzes the ATP-dependent conversion of the Gram-negative peptidoglycan (PG) recycling intermediate 1,6-anhydro-N-acetylmuramic acid (anhMurNAc) to N-acetylmuramic acid-6-phosphate (MurNAc-6-P). Here we present crystal structures of Pseudomonas aeruginosa AnmK in complex with its natural substrate, anhMurNAc, and a product of the reaction, ADP. AnmK is homodimeric, with each subunit comprised of two subdomains that are separated by a deep active site cleft, which bears similarity to the ATPase core of proteins belonging to the hexokinase-hsp70-actin superfamily of proteins. The conversion of anhMurNAc to MurNAc-6-P involves both cleavage of the 1,6-anhydro ring of anhMurNAc along with addition of a phosphoryl group to O6 of the sugar, and thus represents an unusual enzymatic mechanism involving the formal addition of H3PO4 to anhMurNAc. The structural complexes and NMR analysis of the reaction suggest that a water molecule, activated by Asp-182, attacks the anomeric carbon of anhMurNAc, aiding cleavage of the 1,6-anhydro bond and facilitating the capture of the γ phosphate of ATP by O6 via an in-line phosphoryl transfer. AnmK is active only against anhMurNAc and not the metabolically related 1,6-anhydro-N-acetylmuramyl peptides, suggesting that the cytosolic N-acetyl-anhydromuramyl-l-alanine amidase AmpD must first remove the stem peptide from these PG muropeptide catabolites before anhMurNAc can be acted upon by AnmK. Our studies provide the foundation for a mechanistic model for the dual activities of AnmK as a hydrolase and a kinase of an unusual heterocyclic monosaccharide.

Structural studies of ROK fructokinase YdhR from Bacillus subtilis: insights into substrate binding and fructose specificity

The main pathway of bacterial sugar phosphorylation utilizes specific phosphoenolpyruvate phosphotransferase system (PTS) enzymes. In addition to the classic PTS system, a PTS-independent secondary system has been described in which nucleotide-dependent sugar kinases are used for monosaccharide phosphorylation. Fructokinase (FK), which phosphorylates d-fructose with ATP as a cofactor, has been shown to be a member of this secondary system. Bioinformatic analysis has shown that FK is a member of the "ROK" (bacterial Repressors, uncharacterized Open reading frames, and sugar Kinases) sequence family. In this study, we report the crystal structures of ROK FK from Bacillus subtilis (YdhR) (a) apo and in the presence of (b) ADP and (c) ADP/d-fructose. All structures show that YdhR is a homodimer with a monomer composed of two similar α/β domains forming a large cleft between domains that bind ADP and D-fructose. Enzymatic activity assays support YdhR function as an ATP-dependent fructose kinase.

Molecular modeling of the bifunctional enzyme UDP-GlcNAc 2-epimerase/ManNAc kinase and predictions of structural effects of mutations associated with HIBM and sialuria

The bifunctional enzyme UDP-GlcNAc 2-epimerase/ ManNAc kinase (GNE/MNK), encoded by the GNE gene, catalyzes the first two committed, rate-limiting steps in the biosynthesis of N-acetylneuraminic acid (sialic acid). GNE/MNK is feedback inhibited by binding of the downstream product, CMP-sialic acid in its allosteric site. GNE mutations can result in two human disorders, hereditary inclusion body myopathy (HIBM) or sialuria. So far, no active site geometry predictions or conformational transitions involved with function are available for mammalian GNE/MNK. The N-terminal GNE domain is homologous to various prokaryotic 2-epimerases, some of which have solved crystallographic structures. The C-terminal MNK domain belongs to the sugar kinases superfamily; its crystallographic structure is solved at 2.84 A and three-dimensional structures have also been reported for several other kinases. In this work, we employed available structural data of GNE/MNK homologs to model the active sites of human GNE/MNK and identify critical amino acid residues responsible for interactions with substrates. In addition, we modeled effects of GNE/MNK missense mutations associated with HIBM or sialuria on helix arrangement, substrate binding, and enzyme action. We found that all reported mutations are associated with the active sites or secondary structure interfaces of GNE/MNK. The Persian-Jewish HIBM founder mutation p.M712T is located at the interface alpha4alpha10 and likely affects GlcNAc, Mg2+, and ATP binding. This work contributes to further understanding of GNE/MNK function and ligand binding, which may assist future studies for therapeutic options that target misfolded GNE/MNK in HIBM and/or sialuria.

Human sirt-1: molecular modeling and structure-function relationships of an unordered protein

Background

Sirt-1 is a NAD+-dependent nuclear deacetylase of 747 residues that in mammals is involved in various important metabolic pathways, such as glucose metabolism and insulin secretion, and often works on many different metabolic substrates as a multifunctional protein. Sirt-1 down-regulates p53 activity, rising lifespan, and cell survival; it also deacetylases peroxisome proliferator-activated receptor-gamma (PPAR-γ) and its coactivator 1 alpha (PGC-1α), promoting lipid mobilization, positively regulating insulin secretion, and increasing mitochondrial dimension and number. Therefore, it has been implicated in diseases such as diabetes and the metabolic syndrome and, also, in the mechanisms of longevity induced by calorie restriction. Its whole structure is not yet experimentally determined and the structural features of its allosteric site are unknown, and no information is known about the structural changes determined by the binding of its allosteric effectors.

Methodology

In this study, we modelled the whole three-dimensional structure of Sirt-1 and that of its endogenous activator, the nuclear protein AROS. Moreover, we modelled the Sirt-1/AROS complex in order to study the structural basis of its activation and regulation.

Conclusions

Amazingly, the structural data show that Sirt-1 is an unordered protein with a globular core and two large unordered structural regions at both termini, which play an important role in the protein-protein interaction. Moreover, we have found on Sirt-1 a conserved pharmacophore pocket of which we have discussed the implication.

Cohesion group approach for evolutionary analysis of TyrA, a protein family with wide-ranging substrate specificities

Many enzymes and other proteins are difficult subjects for bioinformatic analysis because they exhibit variant catalytic, structural, regulatory, and fusion mode features within a protein family whose sequences are not highly conserved. However, such features reflect dynamic and interesting scenarios of evolutionary importance. The value of experimental data obtained from individual organisms is instantly magnified to the extent that given features of the experimental organism can be projected upon related organisms. But how can one decide how far along the similarity scale it is reasonable to go before such inferences become doubtful? How can a credible picture of evolutionary events be deduced within the vertical trace of inheritance in combination with intervening events of lateral gene transfer (LGT)? We present a comprehensive analysis of a dehydrogenase protein family (TyrA) as a prototype example of how these goals can be accomplished through the use of cohesion group analysis. With this approach, the full collection of homologs is sorted into groups by a method that eliminates bias caused by an uneven representation of sequences from organisms whose phylogenetic spacing is not optimal. Each sufficiently populated cohesion group is phylogenetically coherent and defined by an overall congruence with a distinct section of the 16S rRNA gene tree. Exceptions that occasionally are found implicate a clearly defined LGT scenario whereby the recipient lineage is apparent and the donor lineage of the gene transferred is localized to those organisms that define the cohesion group. Systematic procedures to manage and organize otherwise overwhelming amounts of data are demonstrated.

The structure of bacterial ParM filaments

Bacterial ParM is a homolog of eukaryotic actin and is involved in moving plasmids so that they segregate properly during cell division. Using cryo-EM and three-dimensional reconstruction, we show that ParM filaments have a different structure from F-actin, with very different subunit-subunit interfaces. These interfaces result in the helical handedness of the ParM filament being opposite to that of F-actin. Like F-actin, ParM filaments have a variable twist, and we show that this involves domain-domain rotations within the ParM subunit. The present results yield new insights into polymorphisms within F-actin, as well as the evolution of polymer families.

The crystal structure of Trypanosoma cruzi glucokinase reveals features determining oligomerization and anomer specificity of hexose-phosphorylating enzymes

Glucose is an essential substrate for Trypanosoma cruzi, the protozoan organism responsible for Chagas' disease. The glucose is intracellularly phosphorylated to glucose 6-phosphate. Previously, a hexokinase responsible for this phosphorylation has been characterized. Recently, we identified an ATP-dependent glucokinase in T. cruzi exhibiting a tenfold lower substrate affinity compared to the hexokinase. Both enzymes, which belong to very different groups of the same family, are located inside glycosomes, the peroxisome-like organelles of Kinetoplastida that are known to contain the first seven glycolytic steps as well as enzymes of the oxidative branch of the pentose phosphate pathway. Here, we present the crystallographic structure of T. cruzi glucokinase, in complex with glucose and ADP. The structure suggests a loose tetrameric assembly formed by the association of two tight dimers. TcGlcK was previously reported to exist in a concentration-dependent equilibrium of monomeric and dimeric states. Here, we used mass spectrometry analysis to confirm the existence of TcGlcK monomeric and dimeric states. The analysis of subunit interactions and comparison with the bacterial glucokinases give insights into the forces promoting the stability of the different oligomeric states. Each T. cruzi glucokinase monomer contains one glucose and one ADP molecule. In contrast to hexokinases, which show a moderate preference for the alpha anomer of glucose, the electron density clearly shows the d-glucose bound in the beta configuration in the T.cruzi glucokinase. Kinetic assays with alpha and beta-d-glucose further confirm a moderate preference of the T. cruzi glucokinase for the beta anomer. Structural comparison of the glucokinase and hexokinases permits the identification of a possible mechanism for anomer selectivity in these hexose-phosphorylating enzymes. The preference for distinct anomers suggests that in T. cruzi hexokinase and glucokinase are not directly competing for the same substrate and are probably both present because they exert distinct physiological functions.