Data availability
The cryo-EM maps and atomic coordinates for the eight UGCG structures reported in this study have been deposited in the Electron Microscopy Data Bank and PDB, respectively, under the following accession codes: apo UGCG, EMD-80823 and PDB 26QF; UGCG bound to UDP-glucose, EMD-80832 and PDB 26QS; UGCG bound to UDP, EMD-80833 and PDB 26QT; UGCG bound to UDP-glucose and PS, EMD-80829 and PDB 26QM; UGCG bound to UDP and C6-ceramide, EMD-80824 and PDB 26QG; UGCG bound to miglustat, EMD-80827 and PDB 26QJ; UGCG bound to ibiglustat, EMD-80826 and PDB 26QI; and UGCG bound to eliglustat, EMD-80825 and PDB 26QH. The reference structure of Chlorella virus hyaluronan synthase used for structural comparison is available from the PDB under accession code 7SP7. Structural homology searches were performed against the PDB and AlphaFold Database. For molecular dynamics simulations, simulation input files, including initial system coordinate, topology and force-field files, equilibration configuration files and run scripts, as well as final output coordinate files, are available at Zenodo (https://doi.org/10.5281/zenodo.20809831)45. For sequence data sources: the amino acid sequences of UGCG analysed in this study were obtained from the NCBI protein database. The specific accession numbers for the species used in the phylogenetic analysis are as follows: Prototheria: Ornithorhynchus anatinus (platypus, XP_028909707.1), Tachyglossus aculeatus (echidna, XP_038626167.1); Metatheria: Monodelphis domestica (opossum, XP_001365801.1), Notamacropus eugenii (tammar wallaby, XP_072454167.1), Phascolarctos cinereus (koala, XP_020828114.1), Vombatus ursinus (wombat, XP_027695337.1), Sarcophilus harrisii (Tasmanian devil, XP_003761468.2); Afrotheria: Loxodonta africana (African elephant, XP_003407853.1), Trichechus manatus (manatee, XP_004372220.1), Orycteropus afer (aardvark, XP_007934906.1), Echinops telfairi (tenrec, XP_012862260.1); Xenarthra: Dasypus novemcinctus (armadillo, XP_058158524.1); Laurasiatheria: Erinaceus europaeus (hedgehog, XP_007529571.1), Myotis lucifugus (bat, XP_023609045.1), Equus caballus (Horse, XP_001490953.3), Tursiops truncatus (dolphin/whale, XP_073662506.1), Canis lupus familiaris (dog, NP_001165706.1), Felis catus (cat, XP_011286799.1), Manis javanica (pangolin, XP_017515964.1); Euarchontoglires: Mus musculus (mouse, NP_035803.1), Oryctolagus cuniculus (rabbit, XP_069911685.1), Tupaia belangeri (tree shrew, ELW70583.1), Carlito syrichta (tarsier, XP_008064016.1), Macaca mulatta (monkey, NP_001244652.1) and Homo sapiens (human, NP_003349.1). Source data are provided with this paper.
References
Simons, K. & Ikonen, E. Functional rafts in cell membranes. Nature 387, 569–572 (1997).
Article ADS CAS PubMed Google Scholar
Simons, K. & Toomre, D. Lipid rafts and signal transduction. Nat. Rev. Mol. Cell Biol. 1, 31–39 (2000).
Article CAS PubMed Google Scholar
Hakomori, S. Glycosphingolipids in cellular interaction, differentiation, and oncogenesis. Annu. Rev. Biochem. 50, 733–764 (1981).
Article CAS PubMed Google Scholar
Kasahara, K. & Sanai, Y. Functional roles of glycosphingolipids in signal transduction via lipid rafts. Glycoconj. J. 17, 153–162 (2000).
Article CAS PubMed Google Scholar
Merrill, A. H. Jr. Sphingolipid and glycosphingolipid metabolic pathways in the era of sphingolipidomics. Chem. Rev. 111, 6387–6422 (2011).
Article CAS PubMed PubMed Central Google Scholar
Hannun, Y. A. & Obeid, L. M. Sphingolipids and their metabolism in physiology and disease. Nat. Rev. Mol. Cell Biol. 19, 175–191 (2018).
Article CAS PubMed Google Scholar
Platt, F. M. Emptying the stores: lysosomal diseases and therapeutic strategies. Nat. Rev. Drug Discov. 17, 133–150 (2018).
Article CAS PubMed Google Scholar
Platt, F. M. Sphingolipid lysosomal storage disorders. Nature 510, 68–75 (2014).
Article ADS CAS PubMed Google Scholar
Ichikawa, S., Sakiyama, H., Suzuki, G., Hidari, K. I. & Hirabayashi, Y. Expression cloning of a cDNA for human ceramide glucosyltransferase that catalyzes the first glycosylation step of glycosphingolipid synthesis. Proc. Natl Acad. Sci. USA 93, 12654 (1996).
Article ADS CAS PubMed PubMed Central Google Scholar
Ichikawa, S. & Hirabayashi, Y. Glucosylceramide synthase and glycosphingolipid synthesis. Trends Cell Biol. 8, 198–202 (1998).
Article CAS PubMed Google Scholar
Grabowski, G. A. Phenotype, diagnosis, and treatment of Gaucher’s disease. Lancet 372, 1263–1271 (2008).
Article CAS PubMed Google Scholar
Germain, D. P. Fabry disease. Orphanet. J. Rare Dis. 5, 30 (2010).
Article PubMed PubMed Central Google Scholar
Sidransky, E. et al. Multicenter analysis of glucocerebrosidase mutations in Parkinson’s disease. N. Engl. J. Med. 361, 1651–1661 (2009).
Article CAS PubMed PubMed Central Google Scholar
Mazzulli, J. R. et al. Gaucher disease glucocerebrosidase and α-synuclein form a bidirectional pathogenic loop in synucleinopathies. Cell 146, 37–52 (2011).
Article CAS PubMed PubMed Central Google Scholar
Zunke, F. et al. Reversible conformational conversion of α-synuclein into toxic assemblies by glucosylceramide. Neuron 97, 92–107 (2018).
Article CAS PubMed Google Scholar
Cox, T. M. et al. Eliglustat compared with imiglucerase in patients with Gaucher’s disease type 1 stabilised on enzyme replacement therapy: a phase 3, randomised, open-label, non-inferiority trial. Lancet 385, 2355–2362 (2015).
Article CAS PubMed Google Scholar
Mistry, P. K. et al. Effect of oral eliglustat on splenomegaly in patients with Gaucher disease type 1: the ENGAGE randomized clinical trial. JAMA 313, 695–706 (2015).
Article CAS PubMed PubMed Central Google Scholar
Butters, T. D., Dwek, R. A. & Platt, F. M. Imino sugar inhibitors for treating the lysosomal glycosphingolipidoses. Glycobiology 15, 43r–52r (2005).
Article CAS PubMed Google Scholar
Marks, D. L. et al. Oligomerization and topology of the Golgi membrane protein glucosylceramide synthase. J. Biol. Chem. 274, 451–456 (1999).
Article CAS PubMed Google Scholar
Lairson, L. L., Henrissat, B., Davies, G. J. & Withers, S. G. Glycosyltransferases: structures, functions, and mechanisms. Annu. Rev. Biochem. 77, 521–555 (2008).
Article CAS PubMed Google Scholar
Breton, C., Snajdrová, L., Jeanneau, C., Koca, J. & Imberty, A. Structures and mechanisms of glycosyltransferases. Glycobiology 16, 29r–37r (2006).
Article CAS PubMed Google Scholar
Basu, S., Kaufman, B. & Roseman, S. Enzymatic synthesis of glucocerebroside by a glucosyltransferase from embryonic chicken brain. J. Biol. Chem. 248, 1388–1394 (1973).
Article CAS PubMed Google Scholar
Peterschmitt, M. J. et al. Pharmacokinetics, pharmacodynamics, safety, and tolerability of oral venglustat in healthy volunteers. Clin. Pharmacol. Drug. Dev. 10, 86–98 (2021).
Article CAS PubMed Google Scholar
Platt, F. M., Neises, G. R., Dwek, R. A. & Butters, T. D. N-butyldeoxynojirimycin is a novel inhibitor of glycolipid biosynthesis. J. Biol. Chem. 269, 8362–8365 (1994).
Article CAS PubMed Google Scholar
Wild, R. et al. Structure of the yeast oligosaccharyltransferase complex gives insight into eukaryotic N-glycosylation. Science 359, 545–550 (2018).
Article ADS CAS PubMed Google Scholar
van Kempen, M. et al. Fast and accurate protein structure search with Foldseek. Nat. Biotechnol. 42, 243–246 (2024).
Article ADS PubMed Google Scholar
Maloney, F. P. et al. Structure, substrate recognition and initiation of hyaluronan synthase. Nature 604, 195–201 (2022).
Article ADS CAS PubMed PubMed Central Google Scholar
Coutinho, P. M., Deleury, E., Davies, G. J. & Henrissat, B. An evolving hierarchical family classification for glycosyltransferases. J. Mol. Biol. 328, 307–317 (2003).
Article CAS PubMed Google Scholar
Finstermeier, K. et al. A mitogenomic phylogeny of living primates. PLoS ONE 8, e69504 (2013).
Article ADS CAS PubMed PubMed Central Google Scholar
Steiper, M. E. & Young, N. M. Primate molecular divergence dates. Mol. Phylogenet. Evol. 41, 384–394 (2006).
Article CAS PubMed Google Scholar
Andersson, U. et al. Improved outcome of N-butyldeoxygalactonojirimycin-mediated substrate reduction therapy in a mouse model of Sandhoff disease. Neurobiol. Dis. 16, 506–515 (2004).
Article CAS PubMed Google Scholar
Cai, H. et al. An improved fluorescent tag and its nanobodies for membrane protein expression, stability assay, and purification. Commun. Biol. 3, 753 (2020).
Article CAS PubMed PubMed Central Google Scholar
Punjani, A., Rubinstein, J. L., Fleet, D. J. & Brubaker, M. A. cryoSPARC: algorithms for rapid unsupervised cryo-EM structure determination. Nat. Methods 14, 290–296 (2017).
Article CAS PubMed Google Scholar
Bepler, T. et al. Positive-unlabeled convolutional neural networks for particle picking in cryo-electron micrographs. Nat. Methods 16, 1153–1160 (2019).
Article CAS PubMed PubMed Central Google Scholar
Abramson, J. et al. Accurate structure prediction of biomolecular interactions with AlphaFold 3. Nature 630, 493–500 (2024).
Article ADS CAS PubMed PubMed Central Google Scholar
Pettersen, E. F. et al. UCSF ChimeraX: structure visualization for researchers, educators, and developers. Protein Sci. 30, 70–82 (2021).
Article CAS PubMed Google Scholar
Emsley, P., Lohkamp, B., Scott, W. G. & Cowtan, K. Features and development of Coot. Acta Crystallogr. D Biol. Crystallogr. 66, 486–501 (2010).
Article ADS CAS PubMed PubMed Central Google Scholar
Adams, P. D. et al. PHENIX: a comprehensive Python-based system for macromolecular structure solution. Acta Crystallogr. D Biol. Crystallogr. 66, 213–221 (2010).
Article ADS CAS PubMed PubMed Central Google Scholar
O’Leary, N. A. et al. Reference sequence (RefSeq) database at NCBI: current status, taxonomic expansion, and functional annotation. Nucleic Acids Res. 44, D733–D745 (2016).
Article PubMed Google Scholar
The UniProt Consortium. UniProt: the Universal Protein Knowledgebase in 2023. Nucleic Acids Res. 51, D523–D531 (2023).
Article Google Scholar
Edgar, R. C. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 32, 1792–1797 (2004).
Article CAS PubMed PubMed Central Google Scholar
Jo, S., Kim, T., Iyer, V. G. & Im, W. CHARMM-GUI: a web-based graphical user interface for CHARMM. J. Comput. Chem. 29, 1859–1865 (2008).
Article ADS CAS PubMed Google Scholar
Abraham, M. J. et al. GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX 1–2, 19–25 (2015).
Article ADS Google Scholar
Valdés-Tresanco, M. S., Valdés-Tresanco, M. E., Valiente, P. A. & Moreno, E. gmx_MMPBSA: a new tool to perform end-state free energy calculations with GROMACS. J. Chem. Theory Comput. 17, 6281–6291 (2021).
Article PubMed Google Scholar
Wu, C. & Jin, S. Simulation files for UGCG. Zenodo https://doi.org/10.5281/zenodo.20809831 (2026).
Download references
Acknowledgements
The cryo-EM data were collected at the Shanghai Advanced Center for Electron Microscopy, Shanghai Institute of Materia Medica, Chinese Academy of Sciences.
Funding
This work was supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (grant nos. XDB0830000 and XDB37030103 to H.E.X.); National Natural Science Foundation of China (grant nos. 825B2010 to J.J.W.; 32301016 to C.W.; 32130022 and 82495184 to H.E.X.; 82121005 to H.E.X. and Y.J.; and 82404881 to Q.Y.); The National Key R&D Program of China (grant no. 2022YFC2703105 to H.E.X.); National Key R&D Program ‘Strategic Scientific and Technological Innovation Cooperation’ Key Project (grant no. 2022YFE0203600) released by the Ministry of Science and Technology; Shanghai Municipal Science and Technology Major Project (grant no. 2019SHZDZX02 to H.E.X.); Program of Shanghai Academic/Technology Research Leader (grant no. 22XD1425200 to Y.J.) and the Shanghai Oriental Talents Program (Y.J.).
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature thanks Yu Cao, who co-reviewed with Kexin Hu; Binks W. Wattenberg, who co-reviewed with Usha Mahawar; and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data figures and tables
Extended Data Fig. 1 Chain-length-dependent apparent competition and lateral-access model for ceramide engagement.
a, Apparent competition assay using unlabeled C6–C22 ceramides. Lower C6-NBD-glucosylceramide formation indicates stronger competition with the fluorescent acceptor substrate. C8-ceramide showed the strongest apparent competition, followed by C6-ceramide, whereas longer-chain ceramides showed weaker effects under the same delivery conditions. b, Docked/relaxed ceramide poses in the UGCG acceptor tunnel, comparing C6-, C8- and C16-ceramide. C16-ceramide was used as a representative long-chain ceramide. c, Membrane-coupled lateral-access model. The ceramide headgroup and proximal acyl segment enter the catalytic tunnel, whereas the distal acyl chain remains associated with the bilayer. Insets show representative docked/relaxed poses of C8- and C16-ceramide. Data in a are mean ± s.e.m.; statistical significance was assessed relative to DMSO. For statistical comparisons in a, n = 3 independent experiments. Significance thresholds are defined as follows: n.s., P ≥ 0.05; *P < 0.05; **P < 0.01; ***P < 0.001; ****P < 0.0001. Exact P values are provided in the Source Data.
Source data
Extended Data Fig. 2 Functional and structural characterization of a non-productive phospholipid-occupied state.
a, Lipid inhibition screen of UGCG activity using 18:0–18:1 PC, PE, PG, PS and DAG under matched assay conditions. PS showed the strongest inhibition among the tested lipids, whereas PC and DAG had little effect. Ibiglustat was included as a positive inhibition control. Data are mean ± s.e.m.; statistical significance was assessed relative to DMSO. For statistical comparisons in a, n = 3 independent experiments. Significance thresholds are defined as follows: n.s., P ≥ 0.05; *P < 0.05; **P < 0.01; ***P < 0.001; ****P < 0.0001. Exact P values are provided in the Source Data. b, Dose–response analysis of PS inhibition. The apparent IC50 for 18:0–18:1 PS was 112.2 ± 6.8 μM. c, Cryo-EM density and model of the UDP-glucose-bound, phospholipid-occupied UGCG state. The phospholipid-like density is located in the membrane-facing ceramide-access tunnel. d, Isolated density for the bound phospholipid, shown in two orientations, with the fitted PS model. e, Ligand RMSD of the fitted PS molecule from three independent 500-ns MD simulations, supporting the stability of the PS-bound configuration in the ceramide-access tunnel.
Source data
Extended Data Fig. 3 Ligand-responsive remodeling of the C-terminal β-hairpin.
a, Electrostatic surface representations of apo and UDP-glucose-bound UGCG, showing opening and occlusion of the cytosolic donor-binding pocket. b, Cartoon representations highlighting the C-terminal β-hairpin in the apo and UDP-glucose-bound states. The β-hairpin adopts a more compact conformation over the donor-binding pocket upon UDP-glucose binding. c, Structural superposition of apo and UDP-glucose-bound UGCG showing an approximately 5 Å displacement of the β-hairpin tip toward the donor-binding pocket. The hinge-proximal residue T372 is distal to UDP-glucose, indicating that the effect of the T372P mutation is not caused by loss of a direct ligand contact. d, UDP-glucose-dependent activity of wild-type UGCG and T372P. The T372P substitution increases the apparent Km for UDP-glucose and reduces apparent kcat, supporting a role for local β-hairpin flexibility in productive donor-pocket remodeling. Data are mean ± s.e.m. from independent experiments.
Source data
Extended Data Fig. 4 Molecular dynamics analysis of C-terminal β-hairpin remodeling.
a, Backbone RMSD of the C-terminal β-hairpin in apo and UDP-glucose-bound UGCG from three independent 500-ns MD simulations. Apo UGCG shows greater β-hairpin conformational variability, whereas the UDP-glucose-bound state maintains a more restrained β-hairpin conformation. b, Schematic showing the distance measurement used to compare β-hairpin positioning relative to the catalytic center. The distance was defined between the backbone centroid of the selected β-hairpin segment, residues Ile373–Tyr380, and the backbone atoms of Ser145. c, Violin plot showing the distribution of the β-hairpin centroid–S145 distance from the MD trajectories. UDP-glucose binding shifts the β-hairpin toward S145, consistent with a more compact, inward conformation over the donor-binding pocket. Dotted lines indicate the median and quartiles. Data are mean ± s.e.m. from n = 3 independent experiments.
Extended Data Fig. 5 Cellular and dynamic analysis of the Tyr196 modulatory checkpoint.
a, Cell-based C6-NBD-ceramide glucosylation assay in cells expressing human UGCG WT, human UGCG Y196S, mouse UGCG WT or mouse UGCG S196Y. Formation of C6-NBD-glucosylceramide was quantified and normalized to human UGCG WT. Human Y196S increased cellular glucosylation activity, whereas the reciprocal mouse S196Y substitution reduced activity relative to mouse WT, supporting a conserved activity-dampening effect of tyrosine at position 196. Data are mean ± s.e.m.; statistical significance was assessed relative to the corresponding WT enzyme. For statistical comparisons in a, n = 3 independent experiments. Significance thresholds are defined as follows: n.s., P ≥ 0.05; *P < 0.05; **P < 0.01; ***P < 0.001; ****P < 0.0001. Exact P values are provided in the Source Data. b, Definition of the Tyr196 orientation angle relative to the membrane-normal Z axis. Representative apo and UDP/C6-ceramide-bound UGCG structures illustrate the upward/occluded and downward/open Tyr196 conformations, respectively. The inset shows local contacts associated with the two states, including the Tyr196–R280 contact in the upward/occluded conformation and the Tyr196–E295 contact in the downward/open conformation. c, Tyr196 angle distributions from three independent 500-ns MD simulations across eight UGCG states. Apo, UDP-bound, UDP-glucose-bound and miglustat-bound states predominantly sample the upward/occluded conformation, whereas UDP/C6-ceramide-bound, UDP-glucose/PS-occupied, eliglustat-bound and ibiglustat-bound states preferentially sample the downward/open conformation. Angles >90° correspond to upward/occluded orientations, and angles <90° correspond to downward/open orientations.
Source data
Extended Data Fig. 6 Functional validation of inhibitor-binding modes.
a, Dose–response curves for eliglustat inhibition of wild-type UGCG and selected eliglustat-pocket mutants, including Y196A and Y196S, retaining sufficient basal activity for IC50 determination. b, ΔpIC50 values for eliglustat-pocket mutants relative to wild type. ΔpIC50 was calculated as pIC50 (mutant) − pIC50 (WT); negative values indicate reduced inhibitor potency. Dashed lines indicate 10-fold and 100-fold decreases in apparent potency. c, Dose–response curves for ibiglustat inhibition of wild-type UGCG and selected ibiglustat-pocket mutants, including Y196A and Y196S. d, ΔpIC50 values for ibiglustat-pocket mutants relative to wild type, calculated as in b. e, Dose-response curves for miglustat inhibition of wild-type UGCG and M211A, an activity-retaining probe adjacent to the miglustat-binding environment. f, Eliglustat dose–response curves measured at low and high C6-NBD-ceramide concentrations. The rightward shift at higher C6-NBD-ceramide concentration is consistent with competition at or near the acceptor site. g, Miglustat dose–response curves measured at low and high UDP-glucose concentrations. The rightward shift at higher UDP-glucose concentration is consistent with competition at the donor-binding pocket. Data are mean ± s.e.m. For b and d, statistical significance was assessed relative to WT pIC50 values. For statistical comparisons in b and d, n = 3 independent experiments. Significance thresholds are defined as follows: n.s., P ≥ 0.05; *P < 0.05; **P < 0.01; ***P < 0.001; ****P < 0.0001. Exact P values are provided in the Source Data. Y196A and Y196S had little effect on eliglustat potency but modestly reduced ibiglustat potency.
Source data
Extended Data Fig. 7 Working model for coordinated donor- and acceptor-site engagement by UGCG.
Schematic summary of UGCG substrate engagement, catalysis and inhibition at the Golgi membrane. The apo, UDP-glucose-bound and UDP/C6-ceramide-bound states correspond to experimentally observed structural states. UDP-glucose accesses the cytosolic donor pocket, whereas ceramide engages the membrane-facing acceptor tunnel through lateral access from the bilayer. Miglustat primarily blocks donor-pocket engagement, whereas eliglustat and ibiglustat block acceptor-tunnel engagement. The catalytic assembly and product-like state are shown as inferred models based on the observed structural snapshots. Tyr196 adopts an upward/occluded conformation in the apo and UDP-glucose-bound states and a downward/open conformation in the acceptor-engaged and inferred catalytic states. UDP and glucosylceramide release are depicted without assigning a release order. Enhanced dashed arrows indicate possible or modelled transitions; the order of substrate binding and product release remains unresolved.
Full size table
Full size table
Supplementary information
Source data
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
Reprints and permissions
About this article
Cite this article
Wu, C., Jin, S., Xu, J. et al. Primate-specific regulation of the human glycosphingolipid gatekeeper UGCG. Nature (2026). https://doi.org/10.1038/s41586-026-10927-4
Download citation
Received:
Accepted:
Published:
Version of record:
DOI: https://doi.org/10.1038/s41586-026-10927-4