Phenotypic Mutation 'Grahamcracker' (pdf version)
AlleleGrahamcracker
Mutation Type missense
Chromosome11
Coordinate5,852,165 bp (GRCm39)
Base Change A ⇒ T (forward strand)
Gene Gck
Gene Name glucokinase
Synonym(s) Gls006, hexokinase 4, HK4, MODY2, Hlb62
Chromosomal Location 5,850,820-5,900,081 bp (-) (GRCm39)
MGI Phenotype FUNCTION: [Summary is not available for the mouse gene. This summary is for the human ortholog.] Hexokinases phosphorylate glucose to produce glucose-6-phosphate, the first step in most glucose metabolism pathways. Alternative splicing of this gene results in three tissue-specific forms of glucokinase, one found in pancreatic islet beta cells and two found in liver. The protein localizes to the outer membrane of mitochondria. In contrast to other forms of hexokinase, this enzyme is not inhibited by its product glucose-6-phosphate but remains active while glucose is abundant. Mutations in this gene have been associated with non-insulin dependent diabetes mellitus (NIDDM), maturity onset diabetes of the young, type 2 (MODY2) and persistent hyperinsulinemic hypoglycemia of infancy (PHHI). [provided by RefSeq, Apr 2009]
PHENOTYPE: Targeted disruption of this gene causes mild hyperglycemia in heterozygous mice and extreme hyperglycemia and embryonic to postnatal lethality in homozygous mice. Hyperglycemic knock-out or ENU-induced mutants may show reduced body weight and liver glycogen level, hepatic steatosis, and glucosuria. [provided by MGI curators]
Accession Number

NCBI RefSeq: NM_010292, NM_001287386; MGI:1270854

MappedYes 
Amino Acid Change Asparagine changed to Lysine
Institutional SourceBeutler Lab
Gene Model predicted gene model for protein(s): [ENSMUSP00000099984] [ENSMUSP00000105447] [ENSMUSP00000105448]
AlphaFold P52792
SMART Domains Protein: ENSMUSP00000099984
Gene: ENSMUSG00000041798
AA Change: N391K

DomainStartEndE-ValueType
Pfam:Hexokinase_1 10 217 4.3e-80 PFAM
Pfam:Hexokinase_2 219 458 1.3e-100 PFAM
Predicted Effect probably damaging

PolyPhen 2 Score 1.000 (Sensitivity: 0.00; Specificity: 1.00)
(Using ENSMUST00000102920)
SMART Domains Protein: ENSMUSP00000105447
Gene: ENSMUSG00000041798
AA Change: N391K

DomainStartEndE-ValueType
Pfam:Hexokinase_1 10 217 1e-79 PFAM
Pfam:Hexokinase_2 219 458 7.8e-101 PFAM
Predicted Effect probably damaging

PolyPhen 2 Score 1.000 (Sensitivity: 0.00; Specificity: 1.00)
(Using ENSMUST00000109822)
SMART Domains Protein: ENSMUSP00000105448
Gene: ENSMUSG00000041798
AA Change: N391K

DomainStartEndE-ValueType
Pfam:Hexokinase_1 15 216 1.9e-74 PFAM
Pfam:Hexokinase_2 221 455 2.2e-79 PFAM
Predicted Effect probably damaging

PolyPhen 2 Score 1.000 (Sensitivity: 0.00; Specificity: 1.00)
(Using ENSMUST00000109823)
SMART Domains Protein: ENSMUSP00000123016
Gene: ENSMUSG00000041798

DomainStartEndE-ValueType
low complexity region 8 28 N/A INTRINSIC
Pfam:Hexokinase_2 45 87 1.1e-8 PFAM
Predicted Effect probably benign
Meta Mutation Damage Score 0.6467 question?
Is this an essential gene? Essential (E-score: 1.000) question?
Phenotypic Category Autosomal Semidominant
Candidate Explorer Status loading ...
Single pedigree
Linkage Analysis Data
Penetrance  
Alleles Listed at MGI

All mutations/alleles(34) : Chemically induced (ENU)(18) Gene trapped(5) Targeted(11)

Lab Alleles
AlleleSourceChrCoordTypePredicted EffectPPH Score
IGL01624:Gck APN 11 5853106 missense possibly damaging 0.67
IGL01647:Gck APN 11 5854472 missense probably damaging 0.97
IGL03145:Gck APN 11 5859093 missense probably damaging 0.99
Tootsie UTSW 11 5859150 missense possibly damaging 0.63
R0139:Gck UTSW 11 5860370 missense probably damaging 1.00
R0139:Gck UTSW 11 5859139 nonsense probably null
R0691:Gck UTSW 11 5856691 missense probably damaging 1.00
R1829:Gck UTSW 11 5860984 missense probably damaging 0.97
R1866:Gck UTSW 11 5853253 missense probably benign 0.02
R1868:Gck UTSW 11 5852165 missense probably damaging 1.00
R1992:Gck UTSW 11 5856515 missense probably damaging 1.00
R3885:Gck UTSW 11 5860318 missense probably damaging 1.00
R4179:Gck UTSW 11 5860295 missense probably benign 0.43
R4888:Gck UTSW 11 5859150 missense possibly damaging 0.63
R7034:Gck UTSW 11 5851747 missense probably damaging 1.00
R7155:Gck UTSW 11 5899705 start gained probably benign
R7548:Gck UTSW 11 5852040 missense
R8039:Gck UTSW 11 5860301 missense probably benign 0.12
R8891:Gck UTSW 11 5851733 missense probably damaging 1.00
R9100:Gck UTSW 11 5856516 missense probably damaging 1.00
R9101:Gck UTSW 11 5856516 missense probably damaging 1.00
R9102:Gck UTSW 11 5856516 missense probably damaging 1.00
R9116:Gck UTSW 11 5854377 missense possibly damaging 0.71
R9370:Gck UTSW 11 5852244 missense possibly damaging 0.78
R9420:Gck UTSW 11 5899553 critical splice donor site probably null
R9536:Gck UTSW 11 5852307 missense possibly damaging 0.92
Z1176:Gck UTSW 11 5856526 missense probably damaging 1.00
Z1177:Gck UTSW 11 5860958 missense possibly damaging 0.81
Mode of Inheritance Autosomal Semidominant
Local Stock Live Mice, Sperm, gDNA
MMRRC Submission 038191-MU
Last Updated 2019-09-04 9:46 PM by Katherine Timer
Record Created 2015-01-16 9:53 AM by Jeff SoRelle
Record Posted 2015-05-04
Phenotypic Description

Figure 1. Grahamcracker mice exhibit impaired glucose tolerance 30 minutes after glucose stimulation. Normalized data are shown for pedigrees R0691 (orange) and R1868 (cyan). Abbreviations: WT, wild-type; REF, homozygous reference mice; HET, heterozygous variant mice; VAR, homozygous variant mice. Mean (μ) and standard deviation (σ) are indicated.

The Grahamcracker phenotype was identified among N-ethyl-N-nitrosourea (ENU)-mutagenized G3 mice of the pedigree R1868, some of which exhibited impaired glucose tolerance (i.e., hyperglycemia) 30 minutes after glucose administration (Figure 1).

Nature of Mutation
Figure 2. Linkage mapping of the impaired glucose tolerance 30 minutes after glucose stimulation using an additive model of inheritance. Manhattan plot shows -log10 P values (Y-axis) plotted against the chromosome positions of 170 mutations (X-axis) identified in the G1 males of pedigrees R0691 and R1868 using Superpedigree Analysis. Normalized phenotype data are shown for single locus linkage analysis without consideration of G2 dam identity.  Horizontal pink and red lines represent thresholds of P = 0.05, and the threshold for P = 0.05 after applying Bonferroni correction, respectively.

Whole exome HiSeq sequencing of the G1 grandsire identified 120 mutations. The hyperglycemia phenotype was linked by continuous variable mapping to a mutation in Gck by gene-based Superpedigree Analysis of pedigrees R0691 and R1868. The Gck mutation in pedigree R1868 is a T to A transversion at base pair 5,902,165 (v38) on chromosome 11, or base pair 47,917 in the GenBank genomic region NC_000077 encoding Gck. Linkage was found with an additive model of inheritance in the Superpedigree Analysis, wherein one variant homozygote and 14 heterozygotes departed phenotypically from 7 homozygous reference mice with a P value of 4.175 x 10-5 (Figure 2).  

The mutation corresponds to residue 1,655 in the mRNA sequence NM_010292 (variant 1) within exon 9 of 10 total exons and to residue 1,297 in the mRNA sequence NM_001287386 (variant 2) within exon 9 of 10 total exons.

47900 CTAGCGGGGGTCATAAATCGCATGCGCGAAAGC

386   -L--A--G--V--I--N--R--M--R--E--S-

Genomic numbering corresponds to NC_000077. The mutated nucleotide is indicated in red. The mutation results in an asparagine (N) to lysine (K) substitution at position 391 (N391K) in both isoforms of the glucokinase (GCK) protein, and is strongly predicted by PolyPhen-2 to cause loss of function (score = 1.00) (1).

Illustration of Mutations in
Gene & Protein
Protein Prediction

Figure 3. Domain stucture of GCK. The residues involved in ATP and glucose binding are indicated. The small and large domains refer to the corresponding domains in the crystal structures shown in Figure 4. The Grahamcracker mutation results in an asparagine (N) to lysine (K) substitution at position 391.

Figure 4. Crystal structures of human GCK in an active/closed form (GCK(Δ1-11)) and in an inactive super-open form (GK(Δ1-15). The large and small subdomains as well as the location of the residue affected by the Grahamcracker mutation. The figures were generated by Chimera and are based on PDB:1V4S (active), PDB:1V4T (inactive), and Kamata et al. 2004. PDB:1V4S comprises amino acids 11-465 and PDB:1V4T comprises amino acids 15-465. The image is interactive; click to rotate.

Glucokinase [GCK; alternatively, hexokinase 4 (HK4)] is a member of the hexokinase family (Figure 3). Alternative splicing of Gck generates three tissue-specific 465-amino acid GCK isoforms: GCK-B1, GCK-B2, and GCK-L1 (2-5). GCK-B1 and GCK-B2 are specific to pancreatic β-cells, while GCK-L1 is specific to the liver. At present, it is unknown whether it is the B1 or B2 isoform (or both) that generate the pancreatic GCK enzymatic activity. The GCK-L1 isoform differs from the GCK-B isoforms within the 16 N-terminal amino acids; the functional significance of the N-terminus of all of the isoforms is unknown. The crystal structure of amino acids 11-465 of the active/closed form (GK(Δ1-11); PDB:1V4S) and amino acids 15-465 of the inactive/super-open form (GK(Δ1-15); PDB:1V4T) of human GCK-L1 have been solved [Figure 4; (6)]. The crystal of GK(Δ1-11) in complex with glucose and compound A, an activator, folded into a large and small subdomain corresponding to amino acids 67-203 and 204-443 in mouse GCK, respectively (6). The large and small subdomains are separated by a deep cleft which forms the active site (6). The glucose binding site within the cleft is composed of residues Glu256 and Glu290 of the large subdomain and Thr168 and Lys169 of the small subdomain (6). Ser65 through Gly72 in GK(Δ1-11) comprises connecting region I, which is a flexible structure exposed to the solvent (6). The allosteric site is surrounded by connecting region I, the large subdomain (β1 strand and α5 helix), and the small subdomain (α13 helix) (6). The GK(Δ1-15) crystal was similar to GK(Δ1-11) in overall conformation. However, the spatial relationship of the small subdomains was significantly different between the proteins. Upon glucose binding, the inactive form undergoes a conformational change to the active form. After transition to the open form, the enzymatic reaction is carried out by changing to the closed form in the presence of ATP (6). After the reaction is complete, GCK returns to the open form in order to release glucose-6-phosphate and ADP (6).

The Grahamcracker mutation results in an asparagine (N) to lysine (K) substitution at position 391. N391 is within the GK large subdomain, but is not predicted to be involved in either ATP or glucose binding.

Expression/Localization

GCK is expressed in the liver (the L1 isoform), pancreatic α-, β-, and δ-cells (the B1 and B2 isoforms), the intestine, pituitary, and the hypothalamus (2;5;7-9).

GCK expression and activity are regulated by several factors in both the liver and pancreas. GCK enzyme activity and expression is reduced after 48 hours of starvation in both hepatocytes and pancreatic beta cells (10;11). After refeeding, the GCK enzyme activity and gene expression was restored. In both the liver and pancreas, GK is activated by the enzyme PFK2/fructose-2,6-bisphosphatase-2 (FBPase-2) (12;13). In hepatocytes, Gck expression is induced in the liver by insulin and glucose (3;4;14;15). Insulin stimulates NO production, subsequently leading to S-nitrosylation of GCK (16). Glucagon-like peptide 1 (GLP-1), a factor that potentiates glucose-stimulated insulin secretion in beta cells, increased GCK activity and enhanced GCK S-nitrosylation in βTC3 insulinoma cells (17) .In addition, GCK-B1/B2 is SUMOylated, which subsequently increases its enzymatic activity (18). Insulin receptor substrate-2 (IRS-2)-associated signaling is essential for maintaining the GCK-L1 activity (19). Gck expression is inhibited by glucagon (20). GCK activity in the liver is inhibited by the nuclear protein glucokinase regulatory protein (GRP), which binds GCK to inactivate GCK during starvation (21). After refeeding, GCK is released from GRP in the nucleus and is returned to an active form in the cytosol (21). GCK is also repressed by cAMP in the liver (3;20). In pancreatic beta cells, GCK is regulated by ubiquitination and proteins that have ubiquitin-like domains (e.g., midnolin and parkin) (22;23). Association of GCK with secretory granules in beta cells is regulated by nitric oxide (NO) (16).

Background
Figure 5. The steps in glycolysis. Glycolysis converts glucose to pyruvate in a 10-step process. Each of the steps are numbered. Enzymes and products from each step within the glycolysis cascade are indicated. Step 1, glucose is broken down by hexokinase and is phosphorylated by ATP to form glucose-6-phosphate (G6P). Step 2, isomerization of the G6P by phosphylglucose isomerase to form fructose-6-phosphosphate (Fru-6-P). Step 3, Fru-6-P is phosphorylated by ATP to form fructose 1,6-bisphosphate (Fru-1,6-P2); the enzyme phosphofructokinase controls the entry of sugars into glycolysis. Step 4, Fru-1,6-P2 is cleaved by aldolase to produce dihydroxyacetone phosphate (DHAP) and glyceraldehyde 3-phosphate (G3P). Step 5, the DHAP is isomerized to form G3P. Step 6, the G3P generated from steps 4 and 5 is oxidized by G3P dehydrogenase (GPDH) to form 1,3-biphosphoglycerate (1,3-BPG). Step 7, the phosphate group on 1,3-BPG is transferred to ADP by phosphoglycerate kinase (PGK) to form 3-phosphoglycerate. Step 8, the phosphate ester linkage in 3-phosphoglycerate is moved from carbon 3 to carbon 2 by phosphoglycerate mutase to form 2-phosphoglycerate. Step 9, water is removed from 2-phosphoglycerate by enolase to form phosphoenolpyruvate. Step 10, the phosphate group generated in step 9 is transferred to ADP by pyruvate kinase to form pyruvate and ATP.  See the text for more details.

Glycolysis is a 10-step process that occurs in all life forms in which glucose (C6H12O6) or glycogen (a glucose polymer) is broken down into pyruvate (CH3COCOO + H+) to generate energy for the cells (Figure 5). Each enzymatic step of glycolysis produces a different sugar intermediate. In the first step of glycolysis, glucose is broken down by hexokinase and is phosphorylated by two ATP (adenosine triphosphate) molecules to form glucose-6-phosphate (G6P). In the second step, isomerization of the G6P moves the carbonyl oxygen from carbon 1 to carbon 2, forming fructose-6-phosphosphate (Fru-6-P). In step 3, the hydroxyl group on carbon 1 of Fru-6-P is phosphorylated by ATP to form fructose 1,6-bisphosphate (Fru-1,6-P2); the enzyme phosphofructokinase regulates sugar entry into glycolysis at this step. In step 4, Fru-1,6-P2 is cleaved by aldolase to produce two three-carbon molecules, dihydroxyacetone phosphate (DHAP) and glyceraldehyde 3-phosphate (G3P). In step 5, the DHAP is isomerized to form G3P. In step 6, the G3P generated from steps 4 and 5 is oxidized by G3P dehydrogenase (GPDH) to form 1,3-biphosphoglycerate (1,3-BPG), producing two NADH (nicotinamide adenine dinucleotide) molecules. In step 7, the high-energy phosphate group on 1,3-BPG is transferred to ADP by phosphoglycerate kinase (PGK) to form 3-phosphoglycerate. In step 8, the phosphate ester linkage in 3-phosphoglycerate is moved from carbon 3 to carbon 2 by phosphoglycerate mutase to form 2-phosphoglycerate. In step 9, water is removed from 2-phosphoglycerate by enolase to form phosphoenolpyruvate. In the final step, the phosphate group generated in step 9 is transferred to ADP by pyruvate kinase to form pyruvate and ATP. In total, two ATP and two NADH molecules are released as energy during glycolysis. The pyruvate generated from glycolysis can subsequently enter the Krebs (alternatively, citric acid or tricaboxylic acid) cycle to generate energy in the form of ATP, amino acid precursors, and the reducing agent NADH that is used in other biochemical pathways (e.g., electron transport pathway).

GCK is one of four members of the hexokinase family. The functional significance of GCK in the gut and hypothalamic nuclei in the brain is unknown. In the gut, GCK may promote the secretion of enteroincretins, which are hormones that stimulate a decrease in blood glucose levels (e.g., GLP-1) (24). In the brain, GCK may assist in the regulation of feeding behavior and counter regulatory responses by functioning as a glucose sensor (8;24). In the liver and pancreas, GCK is a component of the ‘glucose sensor’ that regulates plasma glucose levels. Glucose metabolism results in the closure of ATP-sensitive potassium (KATP) channels in the beta-cell plasma membrane and subsequent beta-cell depolarization, activation of voltage-dependent calcium channels, increased calcium influx, and induction of insulin secretion (25-28). Beta cell-specific Gck heterozygous (β-Gck+/-) mice exhibit moderate hyperglycemia and defective insulin secretion in response to glucose (27;28). Pancreatic beta cells from the β-Gck+/- mice exhibited impaired glucose sensitivity that worsened with age (29). β-Gck+/- mice fed a high-fat diet exhibited reduced beta cell replication and beta cell hyperplasia compared to wild-type mice (30). In addition, islets from the β-Gck+/- mice had diminished expression of IRS-2 compared to that in wild-type mice (30). In the liver, GCK function mediates hepatic glucose uptake, the synthesis and subsequent storage of glycogen in the liver, and the regulation of glucose-responsive genes (27;31). Hepatocyte-specific conditional Gck knockout (Liver-Gck-/- ) mice were mildly hyperglycemic, but exhibited defects in glycogen synthesis, glucose turnover rates during hyperglycemic clamp, and impaired insulin secretion in response to glucose (27;32;33). Hepatocyte-specific Gck knockdown in obesity-prone mice attenuated weight gain with a concomitant increase in adaptive thermogenesis (34).

Mutations in GCK are linked to gestational diabetes mellitus [OMIM: #125851; (35)], reduced birth rate (36), late onset noninsulin-dependent diabetes mellitus [OMIM: #125853; (37)], permanent neonatal diabetes mellitus [OMIM: #606176; (38-40)], familial hyperinsulinemic hypoglycemia 3 [OMIM: #602485; (41;42)], and type II maturity-onset diabetes of the young [MODY 2 (alternatively, GCK-MODY); OMIM: #125851; (39;43-45)]. All of the above-mentioned conditions are the result of variable degrees of glucose intolerance due to impaired glucose-responsive insulin secretion.

Putative Mechanism

Gck knockout (Gck-/-) mice are not viable and die from approximately embryonic day (E) 9 until shortly before birth from severe hyperglycemia (27;46;47). Gck heterozygous (Gck+/-) mice survive, but exhibit reduced islet GCK activity and subsequent elevation in fasting blood glucose levels (46-48). After hyperglycemic clamp, the Gck+/- mice exhibited reduced glucose tolerance and defective liver glucose metabolism (46;48). In addition, Gck+/- mice exhibited reduced fertility, increased levels of plasma corticosterone, increased food intake, and hypothalamic gene expression (e.g., increased hypothalamic neuropeptide Y mRNA and reduced hypothalamic proopiomelanocortin mRNA) (31). ENU-induced mutations in Gck resulted in variable reduced viability, reduced GCK activity, and subsequent impaired glucose-responsive insulin secretion (7;49-51). Several homozygote ENU-induced models (GckGENA348/GENA348 and GckD217VD217V/D217V) exhibited increased viability (e.g., 5-12 weeks of age) compared to the Gck-/- mice (49;52). However, the ENU-induced mutant Munich GckM210R/M210R mice exhibited growth retardation and perinatal lethality (51). All heterozygote ENU-induced models exhibited elevated plasma glucose levels, impaired glucose tolerance, and reduced glucose-induced insulin secretion compared to wild-type mice (49;51;52). Similar to other ENU-induced mutants, the Grahamcracker mice exhibit impaired glucose tolerance. In addition, the low numbers of homozygous Grahamcracker mice indicates that the mutation may result in almost complete loss of GCKGrahamcracker function.

Primers PCR Primer
Grahamcracker_pcr_F: GCTCACAAAAGCCTGGAGTTG
Grahamcracker_pcr_R: TGGCATTGTCATCCTCACG

Sequencing Primer
Grahamcracker_seq_F: CTCACAAAAGCCTGGAGTTGAAAGC
Grahamcracker_seq_R: GCATTGTCATCCTCACGACAGAG
Genotyping

PCR program

1) 94°C 2:00
2) 94°C 0:30
3) 55°C 0:30
4) 72°C 1:00
5) repeat steps (2-4) 40x
6) 72°C 10:00
7) 4°C hold


The following sequence of 444 nucleotides is amplified (chromosome 11, - strand):


1   tggcattgtc atcctcacga cagagaaaac ccgggggctg accaaaggac tgctgggcag
61  ggcagggcat tcccattccc ccaggaccct cagtgacttc cctccctaat accgcccgca
121 gcgactctgg ggaccgaagg cagatcctta acatcctgag cactctgggc cttcgaccct
181 ctgtcgccga ctgcgacatt gtgcgccgtg cctgtgaaag cgtgtccact cgcgccgccc
241 acatgtgctc agcaggacta gcgggggtca taaatcgcat gcgcgaaagc cgcagtgagg
301 acgtgatgcg catcacggtg ggcgtggatg gctccgtgta caagctgcac ccgaggtcag
361 cctccacctt cgctcgcgct ccaacctccc gggctccaag ggaggctcca agggcgcggc
421 tttcaactcc aggcttttgt gagc


Primer binding sites are underlined and the sequencing primers are highlighted; the mutated nucleotide is shown in red.

References
Science Writers Anne Murray
Illustrators Peter Jurek
AuthorsZhe Chen, Jeff SoRelle, Noelle Hutchins, Jianhui Wang