Conformational Transition Pathway in the Activation Process of Allosteric Glucokinase

Glucokinase (GK) is a glycolytic enzyme that plays an important role in regulating blood glucose level, thus acting as a potentially attractive target for drug discovery in the treatment of diabetes of the young type 2 and persistent hyperinsulinemic hypoglycemia of infancy. To characterize the activation mechanism of GK from the super-open state (inactive state) to the closed state (active state), a series of conventional molecular dynamics (MD) and targeted MD (TMD) simulations were performed on this enzyme. Conventional MD simulation showed a specific conformational ensemble of GK when the enzyme is inactive. Seven TMD simulations depicted a reliably conformational transition pathway of GK from the inactive state to the active state, and the components important to the conformational change of GK were identified by analyzing the detailed structures of the TMD trajectories. In combination with the inactivation process, our findings showed that the whole conformational pathway for the activation-inactivation-activation of GK is a one-direction circulation, and the active state is less stable than the inactive state in the circulation. Additionally, glucose was demonstrated to gradually modulate its binding pose with the help of residues in the large domain and connecting region of GK during the activation process. Furthermore, the obtained energy barriers were used to explain the preexisting equilibrium and the slow binding kinetic process of the substrate by GK. The simulated results are in accordance with the recent findings from the mutagenesis experiments and kinetic analyses. Our observations reveal a complicated conformational process in the allosteric protein, resulting in new knowledge about the delicate mechanisms for allosteric biological macromolecules that will be useful in drug design for targeting allosteric proteins.

[1]  K. Neet,et al.  Transients and cooperativity. A slow transition model for relating transients and cooperative kinetics of enzymes. , 1972, The Journal of biological chemistry.

[2]  J. Ricard,et al.  Regulatory Behavior of Monomeric Enzymes , 1974 .

[3]  J. Ricard,et al.  Regulatory behavior of monomeric enzymes. 2. A wheat-germ hexokinase as a mnemonical enzyme. , 1974, European journal of biochemistry.

[4]  J. Ricard,et al.  Regulatory behavior of monomeric enzymes. 1. The mnemonical enzyme concept. , 1974, European journal of biochemistry.

[5]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[6]  K. Neet,et al.  [8] Hysteretic enzymes , 1980 .

[7]  K. Neet,et al.  Hysteretic enzymes. , 1980, Methods in enzymology.

[8]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[9]  U. Singh,et al.  A NEW FORCE FIELD FOR MOLECULAR MECHANICAL SIMULATION OF NUCLEIC ACIDS AND PROTEINS , 1984 .

[10]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[11]  K. Neet,et al.  Demonstration of a slow conformational change in liver glucokinase by fluorescence spectroscopy. , 1990, The Journal of biological chemistry.

[12]  P. Tippett,et al.  Observation of a kinetic slow transition in monomeric glucokinase. , 1990, Biochemistry.

[13]  M. Stoffel,et al.  Nonsense mutation in the glucokinase gene causes early-onset non-insulin-dependent diabetes mellitus , 1992, Nature.

[14]  Jürgen Schlitter,et al.  Targeted Molecular Dynamics Simulation of Conformational Change-Application to the T ↔ R Transition in Insulin , 1993 .

[15]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[16]  J. Beckmann,et al.  Familial hyperglycemia due to mutations in glucokinase. Definition of a subtype of diabetes mellitus. , 1993, The New England journal of medicine.

[17]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[18]  P Willett,et al.  Development and validation of a genetic algorithm for flexible docking. , 1997, Journal of molecular biology.

[19]  F M Matschinsky,et al.  Familial hyperinsulinism caused by an activating glucokinase mutation. , 1998, The New England journal of medicine.

[20]  H. Fromm,et al.  Regulation of hexokinase I: crystal structure of recombinant human brain hexokinase complexed with glucose and phosphate. , 1998, Journal of molecular biology.

[21]  B. Glaser,et al.  Pancreatic beta-cell glucokinase: closing the gap between theoretical concepts and experimental realities. , 1998, Diabetes.

[22]  D. Eisenberg,et al.  Protein function in the post-genomic era , 2000, Nature.

[23]  G. Henkelman,et al.  Improved tangent estimate in the nudged elastic band method for finding minimum energy paths and saddle points , 2000 .

[24]  G. Henkelman,et al.  A climbing image nudged elastic band method for finding saddle points and minimum energy paths , 2000 .

[25]  Todd J. A. Ewing,et al.  DOCK 4.0: Search strategies for automated molecular docking of flexible molecule databases , 2001, J. Comput. Aided Mol. Des..

[26]  Klaus Brusgaard,et al.  The second activating glucokinase mutation (A456V): implications for glucose homeostasis and diabetes therapy. , 2002, Diabetes.

[27]  Luhua Lai,et al.  Further development and validation of empirical scoring functions for structure-based binding affinity prediction , 2002, J. Comput. Aided Mol. Des..

[28]  Yufeng Shen,et al.  Conformational pathways in the gating of Escherichia coli mechanosensitive channel , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[29]  M. Tschöp,et al.  Two birds with one stone: novel glucokinase activator stimulates glucose-induced pancreatic insulin secretion and augments hepatic glucose metabolism. , 2003, Molecular interventions.

[30]  E. Schaftingen Short-term regulation of glucokinase , 1994, Diabetologia.

[31]  F. Gago,et al.  Understanding the basis of resistance in the irksome Lys103Asn HIV-1 reverse transcriptase mutant through targeted molecular dynamics simulations. , 2004, Journal of the American Chemical Society.

[32]  Teruyuki Nishimura,et al.  Structural basis for allosteric regulation of the monomeric allosteric enzyme human glucokinase. , 2004, Structure.

[33]  Giuseppe d'Annunzio,et al.  Insights into the Structure and Regulation of Glucokinase from a Novel Mutation (V62M), Which Causes Maturity-onset Diabetes of the Young* , 2005, Journal of Biological Chemistry.

[34]  Alan Saghatelian,et al.  Erratum: Corrigendum: Assignment of protein function in the postgenomic era , 2005, Nature chemical biology.

[35]  Jan Balzarini,et al.  The molecular basis of resilience to the effect of the Lys103Asn mutation in non-nucleoside HIV-1 reverse transcriptase inhibitors studied by targeted molecular dynamics simulations. , 2005, Journal of the American Chemical Society.

[36]  Holger Gohlke,et al.  The Amber biomolecular simulation programs , 2005, J. Comput. Chem..

[37]  D. Nettleton,et al.  Glucose-induced conformational changes in glucokinase mediate allosteric regulation: transient kinetic analysis. , 2006, Biochemistry.

[38]  J. Lindsley,et al.  Whence cometh the allosterome? , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[39]  Jian Zhang,et al.  Conformational transition pathway in the allosteric process of human glucokinase , 2006, Proceedings of the National Academy of Sciences.

[40]  M. Karplus,et al.  The allosteric mechanism of yeast chorismate mutase: a dynamic analysis. , 2006, Journal of molecular biology.

[41]  J. Marcinkeviciene,et al.  A pre-steady state analysis of ligand binding to human glucokinase: evidence for a preexisting equilibrium. , 2007, Biochemistry.

[42]  N. Goodey,et al.  Allosteric regulation and catalysis emerge via a common route. , 2008, Nature chemical biology.

[43]  J. Vanderkooi,et al.  Sugar binding to recombinant wild-type and mutant glucokinase monitored by kinetic measurement and tryptophan fluorescence. , 2008, Biochemical Journal.

[44]  P. Njølstad,et al.  Catalytic activation of human glucokinase by substrate binding – residue contacts involved in the binding of D‐glucose to the super‐open form and conformational transitions , 2008, The FEBS journal.

[45]  Carlos Simmerling,et al.  A Partial Nudged Elastic Band Implementation for Use with Large or Explicitly Solvated Systems. , 2009, International journal of quantum chemistry.

[46]  G. Ferry,et al.  Binding kinetics of glucose and allosteric activators to human glucokinase reveal multiple conformational states. , 2009, Biochemistry.

[47]  Shaoxian Sun,et al.  Biochemical characterization of MODY2 glucokinase variants V62M and G72R reveals reduced enzymatic activities relative to wild type. , 2009, Biochemistry.

[48]  Antoine Taly,et al.  Nicotinic receptors: allosteric transitions and therapeutic targets in the nervous system , 2009, Nature Reviews Drug Discovery.

[49]  F. Matschinsky,et al.  Assessing the potential of glucokinase activators in diabetes therapy , 2009, Nature Reviews Drug Discovery.

[50]  N. Mirza,et al.  Developing analgesics by enhancing spinal inhibition after injury: GABAA receptor subtypes as novel targets. , 2009, Trends in pharmacological sciences.

[51]  M. Larion,et al.  Global fit analysis of glucose binding curves reveals a minimal model for kinetic cooperativity in human glucokinase. , 2010, Biochemistry.

[52]  Geng Wu,et al.  ASD: a comprehensive database of allosteric proteins and modulators , 2010, Nucleic Acids Res..

[53]  Shao-Yong Lu,et al.  Molecular Modeling and Molecular Dynamics Simulation Studies of the GSK3β/ATP/Substrate Complex: Understanding the Unique P+4 Primed Phosphorylation Specificity for GSK3β Substrates , 2011, J. Chem. Inf. Model..

[54]  Yong-Jun Jiang,et al.  Molecular modeling and molecular dynamics simulation studies on pyrrolopyrimidine-based α-helix mimetic as dual inhibitors of MDM2 and MDMX. , 2011, Journal of molecular graphics & modelling.

[55]  Shao-Yong Lu,et al.  Mechanism of kinase inactivation and nonbinding of FRATide to GSK3β due to K85M mutation: molecular dynamics simulation and normal mode analysis. , 2011, Biopolymers.

[56]  Shao-Yong Lu,et al.  Effect of double mutations K214/A–E215/Q of FRATide on GSK3β: insights from molecular dynamics simulation and normal mode analysis , 2011, Amino Acids.