Molecular bases of hexokinase deficiency

The Hexokinase 1 5'-UTR Mutation in Charcot-Marie-Tooth 4G Disease Alters Hexokinase 1 Binding to Voltage-Dependent Anion Channel-1 and Leads to Dysfunctional Mitochondrial Calcium Buffering

Demyelinating Charcot-Marie-Tooth 4G (CMT4G) results from a recessive mutation in the 5'UTR region of the Hexokinase 1 (HK1) gene. HK participates in mitochondrial calcium homeostasis by binding to the Voltage-Dependent Anion Channel (VDAC), through its N-terminal porin-binding domain. Our hypothesis is that CMT4G mutation results in a broken interaction between mutant HK1 and VDAC, disturbing mitochondrial calcium homeostasis. We studied a cohort of 25 CMT4G patients recruited in the French gypsy population. The disease was characterized by a childhood onset, an intermediate demyelinating pattern, and a significant phenotype leading to becoming wheelchair-bound by the fifth decade of life. Co-IP and PLA studies indicated a strong decreased interaction between VDAC and HK1 in the patients' PBMCs and sural nerve. We observed that either wild-type HK1 expression or a peptide comprising the 15 aa of the N-terminal wild-type HK1 administration decreased mitochondrial calcium release in HEK293 cells. However, mutated CMT4G HK1 or the 15 aa of the mutated HK1 was unable to block mitochondrial calcium release. Taken together, these data show that the CMT4G-induced modification of the HK1 N-terminus disrupts HK1-VDAC interaction. This alters mitochondrial calcium buffering that has been shown to be critical for myelin sheath maintenance.

Diet evolution of carnivorous and herbivorous mammals in Laurasiatheria

Background

Laurasiatheria contains taxa with diverse diets, while the molecular basis and evolutionary history underlying their dietary diversification are less clear.

Results

In this study, we used the recently developed molecular phyloecological approach to examine the adaptive evolution of digestive system-related genes across both carnivorous and herbivorous mammals within Laurasiatheria. Our results show an intensified selection of fat and/or protein utilization across all examined carnivorous lineages, which is consistent with their high-protein and high-fat diets. Intriguingly, for herbivorous lineages (ungulates), which have a high-carbohydrate diet, they show a similar selection pattern as that of carnivorous lineages. Our results suggest that for the ungulates, which have a specialized digestive system, the selection intensity of their digestive system-related genes does not necessarily reflect loads of the nutrient components in their diets but appears to be positively related to the loads of the nutrient components that are capable of being directly utilized by the herbivores themselves. Based on these findings, we reconstructed the dietary evolution within Laurasiatheria, and our results reveal the dominant carnivory during the early diversification of Laurasiatheria. In particular, our results suggest that the ancestral bats and the common ancestor of ruminants and cetaceans may be carnivorous as well. We also found evidence of the convergent evolution of one fat utilization-related gene, APOB, across carnivorous taxa.

Conclusions

Our molecular phyloecological results suggest that digestive system-related genes can be used to determine the molecular basis of diet differentiations and to reconstruct ancestral diets.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02033-6.

Novel pathogenic variant c.2714C>A (p. Thr905Lys) in the HK1 gene causing severe haemolytic anaemia with developmental delay in an Indian family

Hexokinase (EC 2.7.1.1, Adenosine Tri Phosphate (ATP): D-hexose-6-phosphotransferase) is a crucial regulatory enzyme of the glycolytic pathway (Embden-Meyerhof pathway). Hexokinase deficiency is associated with chronic non-spherocytic haemolytic anaemia (HA) with some exceptional cases showing psychomotor/mental retardation and fetus death. The proband is a four-and-half-year-old female child born of a four-degree consanguineous marriage hailing from South India with autosomal recessive congenital HA associated with developmental delay. She was well till 3 months of her age post an episode of diarrhoea when she was noted to be severely anaemic and requiring regular transfusions. The common causes of HA, haemoglobinopathies, red cell membranopathies and common red cell enzymopathies (G6PD, GPI, PK and P5N) were ruled out. Targeted analysis of whole exome sequencing (WES) using an insilico gene panel for hereditary anaemia was performed to identify pathogenic variants in the patient. Next-generation sequencing revealed a novel homozygous variant in hexokinase gene c.2714C>A (p. Thr905Lys) in exon-18. The pathogenic nature of the variant p. Thr905Lys in the HK1 gene was confirmed collectively by biochemical and molecular studies. Insilico analysis (PolyPhen-2, Provean, Mutation Taster) predicted the variant to be severe disease causing. Multiple sequence alignment demonstrated the conservation of p. Thr905 across the species. The impact of the mutation on the protein structure was studied by PyMOL and Swiss Protein databank viewer.

Identification of new autoantigens for primary biliary cirrhosis using human proteome microarrays

Primary biliary cirrhosis (PBC) is a chronic cholestatic liver disease of unknown etiology and is considered to be an autoimmune disease. Autoantibodies are important tools for accurate diagnosis of PBC. Here, we employed serum profiling analysis using a human proteome microarray composed of about 17,000 full-length unique proteins and identified 23 proteins that correlated with PBC. To validate these results, we fabricated a PBC-focused microarray with 21 of these newly identified candidates and nine additional known PBC antigens. By screening the PBC microarrays with additional cohorts of 191 PBC patients and 321 controls (43 autoimmune hepatitis, 55 hepatitis B virus, 31 hepatitis C virus, 48 rheumatoid arthritis, 45 systematic lupus erythematosus, 49 systemic sclerosis, and 50 healthy), six proteins were confirmed as novel PBC autoantigens with high sensitivities and specificities, including hexokinase-1 (isoforms I and II), Kelch-like protein 7, Kelch-like protein 12, zinc finger and BTB domain-containing protein 2, and eukaryotic translation initiation factor 2C, subunit 1. To facilitate clinical diagnosis, we developed ELISA for Kelch-like protein 12 and zinc finger and BTB domain-containing protein 2 and tested large cohorts (297 PBC and 637 control sera) to confirm the sensitivities and specificities observed in the microarray-based assays. In conclusion, our research showed that a strategy using high content protein microarray combined with a smaller but more focused protein microarray can effectively identify and validate novel PBC-specific autoantigens and has the capacity to be translated to clinical diagnosis by means of an ELISA-based method.

Studies of a genetic variant in HK1 in relation to quantitative metabolic traits and to the prevalence of type 2 diabetes

BACKGROUND

Single nucleotide polymorphisms (SNPs) within the gene encoding Hexokinase 1 (HK1) are associated with changes in glycated haemoglobin (HbA1c) levels. Our aim was to investigate the effect of HK1 rs7072268 on measures of glucose- and lipid-metabolism in a Danish non-diabetic population and combine the outcome of these analyses in a meta-analysis with previously published results. Furthermore, our aim was to perform a type 2 diabetes case-control analysis and meta-analysis with two previous case-control studies.

METHODS

SNP rs7072268 was genotyped in 9,724 Danes. The quantitative trait study included 5,604 non-diabetic individuals from the Inter99 cohort. The case-control study included 4,449 glucose tolerant individuals and 3,398 patients with type 2 diabetes. Meta-analyses on quantitative traits included 24,560 Caucasian individuals and 30,802 individuals were included in the combined analysis of present and previous type 2 diabetes case-control studies.

RESULTS

Using an additive model, we confirmed that the T-allele of rs7072268 associates with increased HbA1c of 0.6% (CI: 0.4-0.9), p = 3*10-7 per allele. The same allele associated with an increased area under the curve (AUC) for glucose of 5.0 mmol/l*min (0.1-10.0), p = 0.045 following an oral glucose tolerance test (OGTT) and increased fasting levels of cholesterol of 0.06 mmol/l (0.03-1.0), p = 0.001 and triglycerides of 2.0% (0.2-3.8), p = 0.03 per allele in the same study sample of non-diabetic individuals from the Inter99 cohort. However, the T-allele did not show any association with estimates of insulin release or insulin sensitivity neither in Inter99 nor in combined analyses. The prevalence of type 2 diabetes was increased among carriers of the rs7072268 T-allele both in the Danish study-population with an OR of 1.11 (1.02-1.21) and in a meta-analysis including the two additional sample sets with an OR of 1.06 (1.02-1.11). However, after Bonferroni correction the T-allele only remained associated to HbA1c and fasting cholesterol.

CONCLUSIONS

The present study provides suggestive evidence of an association of the rs7072268 T-allele in HK1 with increased AUC glucose following an OGTT in non-diabetic individuals and a nominal association with type 2 diabetes prior to Bonferroni correction. The latter was confirmed in combined analyses involving 16,445 cases and 14,357 control subjects.

Bovine hexokinase type I: full-length cDNA sequence and characterisation of the recombinant enzyme

This study reports the revised and full-length cDNA sequence of bovine hexokinase type I obtained from bovine brain. Since dissimilarities have been observed between the published bovine hexokinase type I coding sequence (GenBank accession no. M65140) (Genomics 11: 1014-1024, 1991) and an analysed portion of bovine hexokinase type I gene, the entire open reading frame was re-sequenced and the ends of cDNA isolated by rapid amplification of cDNA ends. The coding sequences, when compared with the published bovine hexokinase type I, contained a large number of mismatches that lead to changes in the resulting amino acid sequence. The revisions result in a hexokinase type I cDNA of 3619 bp that encodes a protein of 917 amino acids highly homologous to human hexokinase type I. The expression of the recombinant full-length enzyme demonstrated that it was a catalytically active hexokinase. When characterised for its kinetic and regulatory properties, it displayed the same affinity for glucose and MgATP as the human hexokinase type I and was inhibited by glucose 6-phosphate competitively versus MgATP. The production of the N- and C-terminal recombinant halves of the enzyme followed by comparison with the full-length hexokinase indicated that the catalytic activity is located in the C-terminal domain.

Hexokinase: gene structure and mutations

Hexokinase (HK) deficiency is a rare red cell enzyme deficiency associated with hereditary non-spherocytic haemolytic anaemia; to date, only 17 affected families have been reported. Human HK has four major isozymes, each of which is encoded by a separate gene. Recent studies have shown that both ubiquitously expressed type I HK (HK-I) and erythroid-specific HK-R are expressed in erythrocytes, and that these isozymes are encoded by the single HK-I gene. The human HK-I gene has 19 exons, the HK-I and HK-R transcripts being produced by using two distinct promoters. Thus, the first and second exons are specifically utilized for the erythroid-specific HK-R and ubiquitously expressed HK-I isozymes respectively. So far, only two HK variants have been analysed at the molecular level. Since the human HK-I crystal structure has recently been elucidated, the molecular analysis of the HK variants will be useful for discussing the structure-function relationship of the enzyme.

Binding of non-catalytic ATP to human hexokinase I highlights the structural components for enzyme-membrane association control

BACKGROUND

Hexokinase I sets the pace of glycolysis in the brain, catalyzing the ATP-dependent phosphorylation of glucose. The catalytic properties of hexokinase I are dependent on product inhibition as well as on the action of phosphate. In vivo, a large fraction of hexokinase I is bound to the mitochondrial outer membrane, where the enzyme adopts a tetrameric assembly. The mitochondrion-bound hexokinase I is believed to optimize the ATP/ADP exchange between glucose phosphorylation and the mitochondrial oxidative phosphorylation reactions.

RESULTS

The crystal structure of human hexokinase I has been determined at 2.25 A resolution. The overall structure of the enzyme is in keeping with the closed conformation previously observed in yeast hexokinase. One molecule of the ATP analogue AMP-PNP is bound to each N-terminal domain of the dimeric enzyme in a surface cleft, showing specific interactions with the nucleotide, and localized positive electrostatic potential. The molecular symmetry brings the two bound AMP-PNP molecules, at the centre of two extended surface regions, to a common side of the dimeric hexokinase I molecule.

CONCLUSIONS

The binding of AMP-PNP to a protein site separated from the catalytic centre of human hexokinase I can be related to the role played by some nucleotides in dissociating the enzyme from the mitochondrial membrane, and helps in defining the molecular regions of hexokinase I that are expected to be in contact with the mitochondrion. The structural information presented here is in keeping with monoclonal antibody mapping of the free and mitochondrion-bound forms of the enzyme, and with sequence analysis of hexokinases that differ in their mitochondria binding properties.

Enzymatic properties of overexpressed human hexokinase fragments

Full-length hexokinase (HK; ATP: D-hexose 6-phosphotransferase, EC 2.7.1.1), a truncate form of the enzyme lacking the first 11 amino acids (HK-11aa) and the 50 kDa C-terminal half ('mini'-HK) containing the catalytic domain, were overexpressed and purified to homogeneity to investigate the influence of the N-terminal region of human hexokinase type I (HK) on its regulatory properties. All forms of the enzyme are catalytically active with the HK-11aa being the most active. All the forms of HK showed the same affinity for glucose and MgATP and were also inhibited by glucose 6-phosphate (Glc 6-P) competitively vs. MgATP with similar Kis (28.5-37 microM). Glucose 1,6-bisphosphate (Glc 1,6-P2) was also a strong inhibitor of all HKs without significant differences among the different truncate forms of the enzyme (Kis 49.5-59 microM). At low concentrations (0-3 mM), Pi was able to reverse the sugar phosphate inhibition of the full-length HK and HK-11aa but not of the 'mini'-HK. In contrast, at high concentrations Pi was an inhibitor of all the hexokinases investigated. These findings confirm that Pi has a low affinity binding site on the C-terminal of HK while counteracts glucose 6-phosphate inhibition by binding to or requiring the N-terminal half of the enzyme. The first 11 N-terminal amino acids influence the specific activity of HK but are unable to affect the kinetic properties investigated.

Expression, purification, and characterization of a recombinant erythroid-specific hexokinase isozyme

Hexokinase type I (HK I; ATP: D-hexose 6-phosphotransferase, EC 2.7.1.1), the predominant glucose-phosphorylating enzyme in red blood cells, exists in human erythrocytes in multiple molecular forms that differ in isoelectric point and are separable by ion-exchange chromatography. The major forms, designated HK Ia, Ib and Ic, have similar kinetic properties but are characterized by different age-dependent decay and different intracellular distribution in reticulocytes. HK Ib, which elutes between HK I and HK II in the DEAE ion-exchange chromatography, appears to be unique to RBCs and different from any other hexokinase isozyme previously described. Indeed, Murakami and Piomelli recently reported the presence of a specific HK isozyme (named HKr) expressed in K562 cells and in human reticulocytes and, moreover, the resolution of the human HK I gene structure provided the direct evidence of an erythroid-specific exon 1. To further investigate the microheterogeneity of HK I in human RBCs we established a prokaryotic expression system for the HKr isozyme, using the pET plasmid, inducible with IPTG. The recombinant HKr, expressed in bacterial cells as a catalytically active enzyme, was purified to homogeneity by a combination of DEAE ionexchange chromatography followed by hydrophobic interaction chromatography and dye-ligand affinity chromatography. The kinetic and chromatographic properties of the homogeneous recombinant HKr suggest that this erythroid-specific HK isozyme in fact corresponds to the HK isoform previously described in human RBCs and referred to as HK Ib.

Mouse spermatogenic cell-specific type 1 hexokinase (mHk1-s) transcripts are expressed by alternative splicing from the mHk1 gene and the HK1-S protein is localized mainly in the sperm tail

Unique type 1 hexokinase (HK1) mRNAs are present in mouse spermatogenic cells (mHk1-s). They encode a spermatogenic cell-specific sequence region (SSR) but not the porin-binding domain (PBD) necessary for HK1 binding to porin on the outer mitochondrial membrane. This study determined the origin of the multiple Hk1-s transcripts in mouse spermatogenic cells and verified that they are translated in mouse spermatogenic cells. It also showed that a single mHk1 gene encodes the mHk1 transcripts of somatic cells and the mHk1-sa and mHk1-sb transcripts of spermatogenic cells, that alternative exons are used during mHk1 gene expression in mouse spermatogenic cells, and that mHK1-S is translated in mouse spermatogenic cells and is localized mainly with the fibrous sheath in the tail region, not with the mitochondria in the midpiece of mouse sperm.