The glucokinase regulatory protein (GKRP) is a competitive inhibitor of glucokinase (GCK), triggering its localization to the hepatocyte nucleus upon glucose deprivation. Here we establish the kinetic mechanism of GKRP action by analyzing its association with a genetically encoded, fluorescent variant of human GCK. Our results demonstrate that binding of GKRP to GCK involves two steps, formation of an initial encounter complex followed by conformational equilibration between two GKRP–GCK states. Fructose 6-phosphate, a known enhancer of GKRP action, promotes formation of the initial encounter complex via a 2.6-fold increase in kon and stabilizes the complex through a 60-fold decrease in koff.

Cooperativity in human glucokinase (GCK), the body’s primary glucose sensor and a major determinant of glucose homeostatic diseases, is fundamentally different from textbook models of allostery because GCK is monomeric and contains only one glucose-binding site. Prior work has demonstrated that millisecond timescale order-disorder transitions within the enzyme’s small domain govern cooperativity. Here, using limited proteolysis, we map the site of disorder in unliganded GCK to a 30-residue active-site loop that closes upon glucose binding. Positional randomization of the loop, coupled with genetic selection in a glucokinase-deficient bacterium, uncovers a hyperactive GCK variant with substantially reduced cooperativity. Biochemical and structural analysis of this loop variant and GCK variants associated with hyperinsulinemic hypoglycemia reveal two distinct mechanisms of enzyme activation. In α-type activation, glucose affinity is increased, the proteolytic susceptibility of the active site loop is suppressed and the 1H-13C heteronuclear multiple quantum coherence (HMQC) spectrum of 13C-Ile–labeled enzyme resembles the glucose-bound state. In β-type activation, glucose affinity is largely unchanged, proteolytic susceptibility of the loop is enhanced, and the 1H-13C HMQC spectrum reveals no perturbation in ensemble structure. Leveraging both activation mechanisms, we engineer a fully noncooperative GCK variant, whose functional properties are indistinguishable from other hexokinase isozymes, and which displays a 100-fold increase in catalytic efficiency over wild-type GCK. This work elucidates specific structural features responsible for generating allostery in a monomeric enzyme and suggests a general strategy for engineering cooperativity into proteins that lack the structural framework typical of traditional allosteric systems.

Although more than 109 years have passed since the existence of the last universal common ancestor, proteins have yet to reach the limits of divergence. As a result, metabolic complexity is ever expanding. Identifying and understanding the mechanisms that drive and limit the divergence of protein sequence space impact not only evolutionary biologists investigating molecular evolution but also synthetic biologists seeking to design useful catalysts and engineer novel metabolic pathways. Investigations over the past 50 years indicate that the recruitment of enzymes for new functions is a key event in the acquisition of new metabolic capacity. In this review, we outline the genetic mechanisms that enable recruitment and summarize the present state of knowledge regarding the functional characteristics of extant catalysts that facilitate recruitment. We also highlight recent examples of enzyme recruitment, both from the historical record provided by phylogenetics and from enzyme evolution experiments. We conclude with a look to the future, which promises fruitful consequences from the convergence of molecular evolutionary theory, laboratory-directed evolution, and synthetic biology.

Synthetic allosteric activators of human glucokinase are receiving considerable attention as potential diabetes therapeutic agents. Although their mechanism of action is not fully understood, structural studies suggest that activator association requires prior formation of a binary enzyme–glucose complex. Here, we demonstrate that three previously described activators associate with glucokinase in a glucose-independent fashion. Transient-state kinetic assays reveal a lag in enzyme progress curves that is systematically reduced when the enzyme is preincubated with activators. Isothermal titration calorimetry demonstrates that activator binding is enthalpically driven for all three compounds, whereas the entropic changes accompanying activator binding can be favorable or unfavorable. Viscosity variation experiments indicate that the kcat value of glucokinase is almost fully limited by product release, both in the presence and absence of activators, suggesting that activators impact a step preceding product release. The observation of glucose-independent allosteric activation of glucokinase has important implications for the refinement of future diabetes therapeutics and for the mechanism of kinetic cooperativity of mammalian glucokinase.Keywords: glucokinase activators; hysteretic enzyme; kinetic cooperativity; Maturity onset diabetes of the young type II;

Glucokinase (GCK) is responsible for maintaining glucose homeostasis in the human body. Dysfunction or misregulation of GCK causes hyperinsulinemia, hypertriglyceridemia, and type 2 diabetes. In the liver, GCK is regulated by interaction with the glucokinase regulatory protein (GKRP), a 68 kDa polypeptide that functions as a competitive inhibitor of glucose binding to GCK. Formation of the mammalian GCK–GKRP complex is stimulated by fructose 6-phosphate and antagonized by fructose 1-phosphate. Here we report the crystal structure of the mammalian GCK–GKRP complex in the presence of fructose 6-phosphate at a resolution of 3.50 Å. The interaction interface, which totals 2060 Å2 of buried surface area, is characterized by a small number of polar contacts and substantial hydrophobic interactions. The structure of the complex reveals the molecular basis of disease states associated with impaired regulation of GCK by GKRP. It also offers insight into the modulation of complex stability by sugar phosphates. The atomic description of the mammalian GCK–GKRP complex provides a framework for the development of novel diabetes therapeutic agents that disrupt this critical macromolecular regulatory unit.

Cooperativity is widespread in biology. It empowers a variety of regulatory mechanisms and impacts both the kinetic and thermodynamic properties of macromolecular systems. Traditionally, cooperativity is viewed as requiring the participation of multiple, spatially distinct binding sites that communicate via ligand-induced structural rearrangements; however, cooperativity requires neither multiple ligand binding events nor multimeric assemblies. An underappreciated manifestation of cooperativity has been observed in the non-Michaelis–Menten kinetic response of certain monomeric enzymes that possess only a single ligand-binding site. In this review, we present an overview of kinetic cooperativity in monomeric enzymes. We discuss the primary mechanisms postulated to give rise to monomeric cooperativity and highlight modern experimental methods that could offer new insights into the nature of this phenomenon. We conclude with an updated list of single subunit enzymes that are suspected of displaying cooperativity, and a discussion of the biological significance of this unique kinetic response.Graphical abstractHighlights► Cooperativity does not require multiple ligand-binding sites or oligomerization. ► Kinetic cooperativity in monomeric enzymes relies on slow conformational transitions. ► NMR and single-molecule methods can offer new insights into monomeric cooperativity.

High-resolution nuclear magnetic resonance is used to investigate the conformational dynamics of human glucokinase, a 52 kDa monomeric enzyme that displays kinetic cooperativity. 1H−15N transverse relaxation optimized spectra of uniformly labeled glucokinase, recorded in the absence and presence of glucose, reveal significant cross-peak overlap and heterogeneous peak intensities that persist over a range of temperatures. 15N-specific labeling of isoleucines and tryptophans, reporting on backbone and side chain dynamics, respectively, demonstrates that both unliganded and glucose-bound enzymes sample multiple conformations, although glucose stabilizes certain conformations. These results provide the first direct evidence of glucokinase conformational heterogeneity and hence shed light on the molecular basis of cooperativity.

Human pancreatic glucokinase is a monomeric enzyme that displays kinetic cooperativity, a feature that facilitates enzyme-mediated regulation of blood glucose levels in the body. Two theoretical models have been proposed to describe the non-Michaelis−Menten behavior of human glucokinase. The mnemonic mechanism postulates the existence of one thermodynamically favored enzyme conformation in the absence of glucose, whereas the ligand-induced slow transition model (LIST) requires a preexisting equilibrium between two enzyme species that interconvert with a rate constant slower than turnover. To investigate whether either of these mechanisms is sufficient to describe glucokinase cooperativity, a transient-state kinetic analysis of glucose binding to the enzyme was undertaken. A complex, time-dependent change in enzyme intrinsic fluorescence was observed upon exposure to glucose, which is best described by an analytical solution comprised of the sum of four exponential terms. Transient-state glucose binding experiments conducted in the presence of increasing glycerol concentrations demonstrate that three of the observed rate constants decrease with increasing viscosity. Global fit analyses of experimental glucose binding curves are consistent with a kinetic model that is an extension of the LIST mechanism with a total of four glucose-bound binary complexes. The kinetic model presented herein suggests that glucokinase samples multiple conformations in the absence of ligand and that this conformational heterogeneity persists even after the enzyme associates with glucose.

Recently, we reported that YghZ from Escherichia coli functions as an efficient l-glyceraldehyde 3-phosphate reductase (Gpr). Here we show that Gpr co-purifies with a b-type heme cofactor. Gpr associates with heme in a 1:1 stoichiometry to form a complex that is characterized by a Kd value of 5.8 ± 0.2 μM in the absence of NADPH and a Kd value of 11 ± 1.3 μM in the presence of saturating NADPH. The absorbance spectrum of reconstituted Gpr indicates that heme is bound in a hexacoordinate low-spin state under both oxidizing and reducing conditions. The physiological function of heme association with Gpr is unclear, as the l-glyceraldehyde 3-phosphate reductase activity of Gpr does not require the presence of the cofactor. Bioinformatics analysis reveals that Gpr clusters with a family of putative monooxygenases in several organisms, suggesting that Gpr may act as a heme-dependent monooxygenase. The discovery that Gpr associates with heme is interesting because Gpr shares 35% amino acid identity with the mammalian voltage-gated K+ channel β-subunit, an NADPH-dependent oxidoreductase that endows certain voltage-gated K+ channels with hemoprotein-like, O2-sensing properties. To date the molecular origin of O2 sensing by voltage-gated K+ channels is unknown and the results presented herein suggest a role for heme in this process.l-Glyceraldehyde 3-phosphate reductase (Gpr) is capable of binding heme in a hexacoordinate low-spin state.

Microbial niches contain toxic chemicals capable of forcing organisms into periods of intense natural selection to afford survival. Elucidating the mechanisms by which microbes evade environmental threats has direct relevance for understanding and combating the rise of antibiotic resistance. In this study we used a toxic small-molecule, bromoacetate, to model the selective pressures imposed by antibiotics and anthropogenic toxins. We report the results of genetic selection experiments that identify nine genes from Escherichia coli whose overexpression affords survival in the presence of a normally lethal concentration of bromoacetate. Eight of these genes encode putative transporters or transmembrane proteins, while one encodes the essential peptidoglycan biosynthetic enzyme, UDP-N-acetylglucosamine enolpyruvoyl transferase (MurA). Biochemical studies demonstrate that the primary physiological target of bromoacetate is MurA, which becomes irreversibly inactivated via alkylation of a critical active-site cysteine. We also screened a comprehensive library of E. coli single-gene deletion mutants and identified 63 strains displaying increased susceptibility to bromoacetate. One hypersensitive bacterium lacks yliJ, a gene encoding a predicted glutathione transferase. Herein, YliJ is shown to catalyze the glutathione-dependent dehalogenation of bromoacetate with a kcat/Km value of 5.4 × 103 M-1 s-1. YliJ displays exceptional substrate specificity and produces a rate enhancement exceeding 5 orders of magnitude, remarkable characteristics for reactivity with a nonnatural molecule. This study illustrates the wealth of intrinsic survival mechanisms that can be exploited by bacteria when they are challenged with toxins.

The ROK (repressor, open reading frame, kinase) protein family (Pfam 00480) is a large collection of bacterial polypeptides that includes sugar kinases, carbohydrate responsive transcriptional repressors, and many functionally uncharacterized gene products. ROK family sugar kinases phosphorylate a range of structurally distinct hexoses including the key carbon source d-glucose, various glucose epimers, and several acetylated hexosamines. The primary sequence elements responsible for carbohydrate recognition within different functional categories of ROK polypeptides are largely unknown due to a limited structural characterization of this protein family. In order to identify the structural bases for substrate discrimination in individual ROK proteins, and to better understand the evolutionary processes that led to the divergent evolution of function in this family, we constructed an inclusive alignment of 227 representative ROK polypeptides. Phylogenetic analyses and ancestral sequence reconstructions of the resulting tree reveal a discrete collection of active site residues that dictate substrate specificity. The results also suggest a series of mutational events within the carbohydrate-binding sites of ROK proteins that facilitated the expansion of substrate specificity within this family. This study provides new insight into the evolutionary relationship of ROK glucokinases and non-ROK glucokinases (Pfam 02685), revealing the primary sequence elements shared between these two protein families, which diverged from a common ancestor in ancient times.

We describe the discovery of 11 new activating mutations in the human glk gene associated with the disease persistent hyperinsulinemic hypoglycemia of infancy (PHHI). Three of the newly identified substitutions colocalize to a region of the glucokinase polypeptide where a synthetic allosteric activator binds. Of these substitutions, I211F is the most active variant identified to date, with a kcat/K0.5,glucose value (6.6 × 104 M−1 s−1) that is 12-fold higher than that of wild-type glucokinase. The stimulatory mutations described herein represent surreptitious genetic determinants of PHHI. They also identify novel features of the glucokinase scaffold that could be targeted during the development of diabetes therapeutics.

Human glucokinase is a monomeric enzyme that displays a sigmoidal steady-state kinetic response toward increasing glucose concentrations. The allosteric regulation produced by glucose is postulated to arise from the slow interconversion of multiple enzyme conformations during the course of catalysis. Crystallographic data suggest that structural rearrangements linked to glucokinase cooperativity involve a substrate-induced repositioning of an α-helix (α13) located at the C-terminus of the polypeptide. Here, we show that removal of helix α13 abolishes cooperativity and restores Michaelis−Menten kinetics, while reducing the kcat value of the wild-type enzyme by 160-fold. The impaired catalytic activity of the truncated enzyme is not rescued by the trans addition of a synthetic α13 peptide. Unexpectedly, the Km glucose value of a glucokinase variant lacking α13 is equivalent to the K0.5 glucose value of the full-length enzyme. Glucokinase steady-state kinetics is unaffected by the elongation of α13 via the addition of a C-terminal polyalanine tail. To explore the link between cooperativity and the primary sequence of α13, we randomized seven residues within the helix core. Genetic selection experiments in a glucokinase-deficient bacterium identified a variety of hyperactive α13 variants that display lower K0.5 glucose values, Hill coefficients near unity, and enhanced equilibrium binding affinities for glucose. The present results demonstrate that α13 plays an essential role in facilitating cooperativity. Our findings also establish a link between the primary amino acid sequence of helix α13 and the functional dynamics of the glucokinase scaffold that are required for allostery.

Triosephosphate isomerase (TIM) catalyzes the interconversion of d-glyceraldehyde 3-phosphate and dihydroxyacetone phosphate, an essential step in glycolytic and gluconeogenic metabolism. To uncover promiscuous isomerases embedded within the Escherichia coli genome, we searched for genes capable of restoring growth of a TIM-deficient bacterium under gluconeogenic conditions. Rather than discovering an isomerase, we selected yghZ, a gene encoding a member of the aldo-keto reductase superfamily. Here we show that YghZ catalyzes the stereospecific, NADPH-dependent reduction of l-glyceraldehyde 3-phosphate, the enantiomer of the TIM substrate. This transformation provides an alternate pathway to the formation of dihydroxyacetone phosphate.