1
|
Sinclair M, Stein RA, Sheehan JH, Hawes EM, O’Brien RM, Tajkhorshid E, Claxton DP. Integrative analysis of pathogenic variants in glucose-6-phosphatase based on an AlphaFold2 model. PNAS NEXUS 2024; 3:pgae036. [PMID: 38328777 PMCID: PMC10849595 DOI: 10.1093/pnasnexus/pgae036] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/01/2023] [Accepted: 01/09/2024] [Indexed: 02/09/2024]
Abstract
Mediating the terminal reaction of gluconeogenesis and glycogenolysis, the integral membrane protein glucose-6-phosphate catalytic subunit 1 (G6PC1) regulates hepatic glucose production by catalyzing hydrolysis of glucose-6-phosphate (G6P) within the lumen of the endoplasmic reticulum. Consistent with its vital contribution to glucose homeostasis, inactivating mutations in G6PC1 causes glycogen storage disease (GSD) type 1a characterized by hepatomegaly and severe hypoglycemia. Despite its physiological importance, the structural basis of G6P binding to G6PC1 and the molecular disruptions induced by missense mutations within the active site that give rise to GSD type 1a are unknown. In this study, we determine the atomic interactions governing G6P binding as well as explore the perturbations imposed by disease-linked missense variants by subjecting an AlphaFold2 G6PC1 structural model to molecular dynamics simulations and in silico predictions of thermodynamic stability validated with robust in vitro and in situ biochemical assays. We identify a collection of side chains, including conserved residues from the signature phosphatidic acid phosphatase motif, that contribute to a hydrogen bonding and van der Waals network stabilizing G6P in the active site. The introduction of GSD type 1a mutations modified the thermodynamic landscape, altered side chain packing and substrate-binding interactions, and induced trapping of catalytic intermediates. Our results, which corroborate the high quality of the AF2 model as a guide for experimental design and to interpret outcomes, not only confirm the active-site structural organization but also identify previously unobserved mechanistic contributions of catalytic and noncatalytic side chains.
Collapse
Affiliation(s)
- Matt Sinclair
- Theoretical and Computational Biophysics Group, NIH Center for Macromolecular Modeling and Visualization, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
- Department of Biochemistry, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Richard A Stein
- Department of Molecular Physiology and Biophysics, Vanderbilt University, Nashville, TN 37232, USA
- Center for Applied Artificial Intelligence in Protein Dynamics, Vanderbilt University, Nashville, TN 37240, USA
| | - Jonathan H Sheehan
- Center for Structural Biology, Vanderbilt University, Nashville, TN 37240, USA
- Division of Infectious Diseases, Department of Internal Medicine, Washington University School of Medicine, St Louis, MO 63110, USA
| | - Emily M Hawes
- Department of Molecular Physiology and Biophysics, Vanderbilt University, Nashville, TN 37232, USA
| | - Richard M O’Brien
- Department of Molecular Physiology and Biophysics, Vanderbilt University, Nashville, TN 37232, USA
| | - Emad Tajkhorshid
- Theoretical and Computational Biophysics Group, NIH Center for Macromolecular Modeling and Visualization, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
- Department of Biochemistry, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
- Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Derek P Claxton
- Department of Molecular Physiology and Biophysics, Vanderbilt University, Nashville, TN 37232, USA
- Center for Applied Artificial Intelligence in Protein Dynamics, Vanderbilt University, Nashville, TN 37240, USA
- Center for Structural Biology, Vanderbilt University, Nashville, TN 37240, USA
| |
Collapse
|
2
|
Sinclair M, Stein RA, Sheehan JH, Hawes EM, O'Brien RM, Tajkhorshid E, Claxton DP. Molecular mechanisms of catalytic inhibition for active site mutations in glucose-6-phosphatase catalytic subunit 1 linked to glycogen storage disease. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.03.13.532485. [PMID: 36993754 PMCID: PMC10054992 DOI: 10.1101/2023.03.13.532485] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/19/2023]
Abstract
Mediating the terminal reaction of gluconeogenesis and glycogenolysis, the integral membrane protein G6PC1 regulates hepatic glucose production by catalyzing hydrolysis of glucose-6-phosphate (G6P) within the lumen of the endoplasmic reticulum. Consistent with its vital contribution to glucose homeostasis, inactivating mutations in G6PC1 cause glycogen storage disease (GSD) type 1a characterized by hepatomegaly and severe hypoglycemia. Despite its physiological importance, the structural basis of G6P binding to G6PC1 and the molecular disruptions induced by missense mutations within the active site that give rise to GSD type 1a are unknown. Exploiting a computational model of G6PC1 derived from the groundbreaking structure prediction algorithm AlphaFold2 (AF2), we combine molecular dynamics (MD) simulations and computational predictions of thermodynamic stability with a robust in vitro screening platform to define the atomic interactions governing G6P binding as well as explore the energetic perturbations imposed by disease-linked variants. We identify a collection of side chains, including conserved residues from the signature phosphatidic acid phosphatase motif, that contribute to a hydrogen bonding and van der Waals network stabilizing G6P in the active site. Introduction of GSD type 1a mutations into the G6PC1 sequence elicits changes in G6P binding energy, thermostability and structural properties, suggesting multiple pathways of catalytic impairment. Our results, which corroborate the high quality of the AF2 model as a guide for experimental design and to interpret outcomes, not only confirm active site structural organization but also suggest novel mechanistic contributions of catalytic and non-catalytic side chains.
Collapse
|
3
|
Tang Y, Cheng F, Feng Z, Jia G, Li C. Stereostructural Elucidation of Glucose Phosphorylation by Raman Optical Activity. J Phys Chem B 2019; 123:7794-7800. [PMID: 31335146 DOI: 10.1021/acs.jpcb.9b05968] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
Phosphorylation of glucose is the prime step in sugar metabolism and energy storage. Two key glucose phosphates are involved, that is, glucose 6-phosphate (G6P) and α-glucose 1-phosphate (αG1P). The chiral conformation of glucose, G6P, and αG1P plays an essential role in enzyme-mediated conversions. However, few techniques were able to give a direct view of the conformational changes from glucose to G6P and αG1P. Here, Raman optical activity (ROA) was used to elucidate the stereochemical evolution of glucose upon phosphorylation. ROA was found to be extremely sensitive to different phosphorylation sites. A characteristic ROA marker of (+)980 cm-1, originated from the phosphate group symmetric stretching vibration, is observed for αG1P with phosphorylation at chiral C1, while no corresponding ROA signal for G6P with phosphorylation at achiral C6 is observed. Phosphorylation-induced gauch-gauch (gg)/gauch-trans (gt) rotamer distribution changes can be sensitively probed by the sign of the ROA band around 1460 cm-1. A positive ROA band at 1465 cm-1 of glucose corresponds to a higher gt ratio, while a negative band at 1455 cm-1 of G6P suggests a dominant gg population, and the disappearance of this ROA band for αG1P indicates a nearly balanced gg/gt distribution. Meanwhile, the phosphorylation at C6 and C1 could cause dramatic reduction of the conformational flexibility of the adjacent C4-OH and C2-OH, respectively. These stereochemical changes revealed by ROA spectra offer a structural basis on the understanding of sugar phosphorylation from the perspective of chirality.
Collapse
Affiliation(s)
- Yuxuan Tang
- State Key Laboratory of Catalysis, Dalian Institute of Chemical Physics , Chinese Academy of Sciences , Zhongshan Road 457 , Dalian 116023 , China.,University of Chinese Academy of Sciences , No. 19A Yuquan Road , Beijing 100049 , China
| | - Feng Cheng
- State Key Laboratory of Catalysis, Dalian Institute of Chemical Physics , Chinese Academy of Sciences , Zhongshan Road 457 , Dalian 116023 , China
| | - Zhaochi Feng
- State Key Laboratory of Catalysis, Dalian Institute of Chemical Physics , Chinese Academy of Sciences , Zhongshan Road 457 , Dalian 116023 , China
| | - Guoqing Jia
- State Key Laboratory of Catalysis, Dalian Institute of Chemical Physics , Chinese Academy of Sciences , Zhongshan Road 457 , Dalian 116023 , China
| | - Can Li
- State Key Laboratory of Catalysis, Dalian Institute of Chemical Physics , Chinese Academy of Sciences , Zhongshan Road 457 , Dalian 116023 , China
| |
Collapse
|
4
|
Abstract
The anomeric specificity of D-glucose metabolism in intact hepatocytes remains a matter of debate. This issue was further investigated in the present study, which is based on the quantification of the alpha- and beta-anomers of the 13C-enriched isotopomers of D-glucose generated by rat liver cells exposed to either D-[1-13C] fructose or D-[2-13C] fructose in the presence of D2O. The D-[1-13C]glucose/D-[6-13C]glucose paired ratios found in the cells exposed to D-[1-13C] fructose and the D-[2-13C]glucose/D-[5-13C]glucose paired ratios found in the cells exposed to D-[2-13C] fructose yielded a paired beta/alpha ratio averaging (mean +/- S.E.M.) 79.3 +/- 6.1%. In the case of the isotopomers of D-glucose formed by gluconeogenesis, the D-[2-13C]glucose/D-[5-13C]glucose and D-[3-13C]glucose/D-[4-13C]glucose paired ratios found in cells exposed to D-[1-13C] fructose, as well as the D-[1-13C]glucose/D-[6-13C]glucose and D-[3-13C]glucose/D-[4-13C]glucose paired ratios found in cells exposed to D-[2-13C]fructose, yielded an alpha/beta paired ratio averaging 75.0 +/- 5.8%. Last, in the cells exposed to D-[2-13C]fructose, the beta/alpha ratio for the C2-deuterated isotopomers of D-[2-13C]glucose represented 78.9 +/- 3.7% of that for the C5-deuterated isotopomers of D-[5-13C]glucose. The three values representative of the anomeric specificity of D-glucose production by liver cells were not significantly different from one another, with an overall mean value of 76.9 +/- 3.6%. These findings unambiguously document that the anomeric specificity of phosphoglucoisomerase is operative in intact hepatocytes, resulting in a preferential output of the alpha-anomer of 13C-enriched D-glucose under the present experimental conditions.
Collapse
Affiliation(s)
- W J Malaisse
- Laboratory of Experimental Hormonology (ULB) and High Resolution NMR Centre (VUB), Brussels Free Universities, Brussels, Belgium.
| | | |
Collapse
|
5
|
Namba N, Watanabe F, Tokuda M, Mino M, Furuya E. A new method of quantitating serum and urinary levels of 1,5-anhydroglucitol in insulin-dependent diabetes mellitus. Diabetes Res Clin Pract 1994; 24:55-61. [PMID: 7924887 DOI: 10.1016/0168-8227(94)90086-8] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
A new method was developed for quantitating the serum and urinary levels of 1,5-anhydroglucitol (AG), a sensitive and informative marker of glycemic control. This method utilized a combination of ODS and pyranose oxidase-immobilized columns for HPLC, and monitored hydrogen peroxide production with an electrochemical detector. We applied this method to determine the serum and urinary AG levels in 15 patients with insulin-dependent diabetes mellitus (IDDM) as well as in control subjects. Baseline separation of AG from other sugars such as glucose and myoinositol was achieved. Quantitation of AG was achieved over the range from 0.2 ng to 0.3 micrograms based upon peak heights. The serum and urinary AG levels in the IDDM patients were 4.4 +/- 8.3 mg/l and 5.1 +/- 4.3 mg/day, respectively. We found that the urinary AG to serum AG ratio showed a linear correlation with the urinary glucose level in the IDDM patients (urinary glucose (y) vs. urinary AG to serum AG ratio (x): y = 9.071x-0.991; r = 0.968, P < 0.001). This method proved efficient and reliable for quantitating urinary AG. Since determination of both the AG and glucose levels in urine gives equivalent clinical information to the serum AG level, urinary monitoring could provide a valuable addition to the available methods for assessing the glycemic status of IDDM patients.
Collapse
Affiliation(s)
- N Namba
- Department of Pediatrics, Osaka Medical College, Japan
| | | | | | | | | |
Collapse
|
6
|
Okuyama H, Kamata S, Ishikawa S, Usui N, Okada A, Yamatodani A. Effects of synchronous pulsatile extracorporeal membrane oxygenation in an endotoxin-induced shock model: an experimental study. Artif Organs 1992; 16:477-84. [PMID: 10078296 DOI: 10.1111/j.1525-1594.1992.tb00327.x] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A synchronous pulsatile venoarterial device for extracorporeal membrane oxygenation (ECMO) was designed to accomplish more effective circulatory support for neonates. The effect of this device was studied using an endotoxin-induced shock model, compared with conventional nonpulsatile ECMO. Twenty puppies weighing 1.6 to 4.0 kg were given endotoxin (5 mg/kg) intravenously. Thirty minutes after the administration of endotoxin, 10 were placed on pulsatile ECMO, and the others were placed on nonpulsatile ECMO, and they were studied for an additional 180 min. Peak blood pressure, arterial pH, base excess, and renal blood flow were significantly higher in the pulsatile group than in the nonpulsatile group. Serum lactate and serum noradrenaline were significantly lower in the pulsatile group than in the nonpulsatile group. These results indicate that pulsatile ECMO may provide more effective cardiopulmonary support in the treatment of neonates with serious circulatory failure that has failed to be supported by nonpulsatile ECMO.
Collapse
Affiliation(s)
- H Okuyama
- Department of Pediatric Surgery, Osaka University Medical School, Japan
| | | | | | | | | | | |
Collapse
|
7
|
Takami H, Matsuda H, Kuki S, Nishimura M, Kawashima Y, Watari H, Furuya E, Tagawa K. Leakage of cytoplasmic enzymes from rat heart by the stress of cardiac beating after increase in cell membrane fragility by anoxia. Pflugers Arch 1990; 416:144-50. [PMID: 2352832 DOI: 10.1007/bf00370236] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
The effects of spontaneous beating after anoxia and the pumping stress induced by a left ventricular balloon on the leakage of myocardial enzymes from the isolated perfused rat heart were investigated. Beating of the heart was arrested by perfusion with high-K+ medium. When the beating was arrested during reoxygenation after anoxia, the leakage of lactate dehydrogenase (LDH) was significantly lower than during reoxygenation with spontaneous cardiac beating. After changing from K+ arrest to spontaneous beating by perfusion with low-K+ medium during reoxygenation, the leakage of LDH increased markedly. Imposition of left ventricular wall stress on the K(+)-arrested heart by repetitive passive distension during aerobic perfusion and after 20 min and 60 min of anoxia caused LDH leakages of 1.0, 4.6 and 21.0 units/g in 30 min, respectively. Under this mechanical stress, the release of LDH as a percentage of its total myocardial activity coincided well with that of cytoplasmic aspartate aminotransferase (AST), while the percentage release of mitochondrial AST was much less. These results appeared to indicate that the leakage of cytoplasmic enzymes during reoxygenation is accelerated by cardiac beating because of fragility of the cell membranes developing during the preceding anoxia.
Collapse
Affiliation(s)
- H Takami
- First Department of Surgery, Osaka University Medical School, Japan
| | | | | | | | | | | | | | | |
Collapse
|
8
|
Ramirez R, Zähner D, Marynissen G, Sener A, Malaisse WJ. Anomeric specificity of D-glucose phosphorylation by rat liver glucose-6-phosphatase. Biochem J 1989; 261:509-13. [PMID: 2549974 PMCID: PMC1138855 DOI: 10.1042/bj2610509] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
The anomeric specificity of D-glucose phosphorylation by hepatic glucose-6-phosphatase was examined in rat liver microsomes incubated in the presence of carbamoyl phosphate. At 10 degrees C, the Km for the equilibrated hexose and phosphate donor was close to 56 mM and 11 mM, respectively. The enzymic activity, which was increased in diabetic rats, was about 40% lower in untreated than in sonicated microsomes. No anomeric difference in affinity was found in sonicated microsomes. In untreated microsomes, however, the Km for beta-D-glucose was slightly lower than that for alpha-D-glucose. The maximal velocity was higher with beta- than alpha-D-glucose in both untreated and sonicated microsomes. These data indicate that the phosphotransferase activity of glucose-6-phosphatase cannot account for the higher rate of glycolysis and glycogen synthesis found in hepatocytes exposed to alpha- rather than beta-D-glucose.
Collapse
Affiliation(s)
- R Ramirez
- Laboratory of Experimental Medicine, Brussels Free University, Belgium
| | | | | | | | | |
Collapse
|
9
|
Achs MJ, Garfinkel D. Pancreatic islet discrimination of hexose anomers. I. Steady-state computer simulation. THE AMERICAN JOURNAL OF PHYSIOLOGY 1988; 255:E189-200. [PMID: 2970227 DOI: 10.1152/ajpendo.1988.255.2.e189] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
Pancreatic islets detect glucose level by phosphorylating it and converting the glycolytic rate to a signal to secrete insulin. Insulin secretion is greater from the alpha- than from the beta-anomer when the D-glucose level is below 22 mM. D-mannose behaves similarly but at nearly twofold higher concentrations. Two explanations have been proposed: 1) glucokinase, which has the same anomeric preference, is the principal hexose phosphorylating enzyme and limits glycolytic rate. 2) Phosphofructokinase limits glycolysis and hexokinase is the principal enzyme phosphorylating hexose; hexosediphosphate activators of phosphofructokinase are more readily synthesized from alpha-anomers of hexose phosphates. We have simulated both alternatives with a detailed anomerically specific model of the hexose-metabolizing glycolytic enzymes. The pathway preference for alpha-anomer of both hexoses was adequately reproduced with anomerically active limiting glucokinase. The other mechanism did not reproduce the observed pathway preference.
Collapse
Affiliation(s)
- M J Achs
- Department of Computer Science, University of Pennsylvania, Philadelphia 19104
| | | |
Collapse
|
10
|
Zähner D, Rasschaert J, Malaisse WJ. Anomeric specificity of liver glycogenolysis. BIOCHIMICA ET BIOPHYSICA ACTA 1987; 926:115-8. [PMID: 3663706 DOI: 10.1016/0304-4165(87)90226-1] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
In rat liver slices incubated in the absence of exogenous D-glucose, both the basal and glucagon-stimulated output of D-glucose resulted in the production of a greater relative amount of alpha-D-glucose than that found at anomeric equilibrium. Comparable results were obtained in isolated hepatocytes. In these experiments, the rate of glycogenolysis largely exceeded that of glycogen synthesis. These findings indicate that liver glycogenolysis represents an alpha-stereospecific process.
Collapse
Affiliation(s)
- D Zähner
- Laboratory of Experimental Medicine, Brussels Free University, Belgium
| | | | | |
Collapse
|