Truncated mutants of beta-glucosidase 2 (GBA2) are localized in the mitochondrial matrix and cause mitochondrial fragmentation

The proximity interactome of PML isoforms I and II under fatty acid stress

Promyelocytic leukemia (PML) protein forms the scaffold for PML nuclear bodies (PML NB) that reorganize into Lipid-Associated PML Structures (LAPS) under fatty acid stress. We determined how the fatty acid oleate alters the interactome of PMLI or PMLII by expressing fusions with the ascorbate peroxidase APEX2 in U2OS cells. The resultant interactome included ESCRT and COPII transport protein nodes. Proximity ligation assay (PLA) revealed that COPII proteins SEC23B, SEC24A and USO1 preferentially associated with PML NBs. Nuclear localization of USO1, but not SEC23B and SEC24A, was reduced in PML knockout cells and restored by PMLII expression. Thus, proximity-labelling methods identified COPII transport protein interactions with PML NBs that are disrupted by fatty acid stress.

Hereditary spastic paraparesis type 46 (SPG46): new GBA2 variants in a large Italian case series and review of the literature

Hereditary spastic paraparesis (HSP) is a group of central nervous system diseases primarily affecting the spinal upper motor neurons, with different inheritance patterns and phenotypes. SPG46 is a rare, early-onset and autosomal recessive HSP, linked to biallelic GBA2 mutations. About thirty families have been described worldwide, with different phenotypes like complicated HSP, recessive cerebellar ataxia or Marinesco-Sjögren Syndrome. Herein, we report five SPG46 patients harbouring five novel GBA2 mutations, the largest series described in Italy so far. Probands were enrolled in five different centres and underwent neurological examination, clinical cognitive assessment, column imaging for scoliosis assessment, ophthalmologic examination, brain imaging, GBA2 activity in peripheral blood cells and genetic testing. Their phenotype was consistent with HSP, with notable features like upper gaze palsy and movement disorders. We review demographic, genetic, biochemical and clinical information from all documented cases in the existing literature, focusing on the global distribution of cases, the features of the syndrome, its variable presentation, new potential identifying features and the significance of measuring GBA2 enzyme activity.

Structural basis for inhibition of a GH116 β-glucosidase and its missense mutants by GBA2 inhibitors: Crystallographic and quantum chemical study

The crystal structure of the Thermoanaerobacterium xylanolyticum in glycoside hydrolase family 116 (TxGH116) β-glucosidase provides a structural model for human GBA2 glucosylceramidase, an enzyme defective in hereditary spastic paraplegia and a potential therapeutic target for treating Gaucher disease. To assess the therapeutic potential of known inhibitors, the X-ray structure of TxGH116 in complex with isofagomine (IFG) was determined at 2.0 Å resolution and showed the IFG bound in a relaxed chair conformation. The binding of IFG and 7 other iminosugar inhibitors to wild-type and mutant enzymes (Asp508His and Arg786His) mimicking GBA2 pathogenic variants was then evaluated computationally by two-layered ONIOM calculations (at the B3LYP:PM7 level). Calculations showed that six charged residues, Glu441, Asp452, His507, Asp593, Glu777, and Arg786 influence inhibitor binding most. His507, Glu777 and Arg786, form strong hydrogen bonds with the inhibitors (∼1.4-1.6 Å). Thus, the missense mutation of one of these residues in Arg786His has a greater effect on the interaction energies for all inhibitors compared to Asp508His. In line with the experimental data for the inhibitors that have been tested, the favorable interaction energy between the inhibitors and the TxGH116 protein followed the trend: isofagomine > 1-deoxynojirimycin > glucoimidazole > N-butyl-deoxynojirimycin ≈ N-nonyl-deoxynojirimycin > conduritol B epoxide ≈ azepane 1 > azepane 2. The obtained structural and energetic properties and comparison to the GBA2 model can lead to understanding of structural requirement for inhibitor binding in GH116 to aid the design of high potency GBA2 inhibitors.

Systematic Functional and Computational Analysis of Glucose-Binding Residues in Glycoside Hydrolase Family GH116

Glycoside hydrolases (GH) bind tightly to the sugar moiety at the glycosidic bond being hydrolyzed to stabilize its transition state conformation. We endeavored to assess the importance of glucose-binding residues in GH family 116 (GH116) β-glucosidases, which include human β-glucosylceramidase 2 (GBA2), by mutagenesis followed by kinetic characterization, X-ray crystallography, and ONIOM calculations on Thermoanaerobacterium xylanolyticum TxGH116, the structural model for GH116 enzymes. Mutations of residues that bind at the glucose C3OH and C4OH caused 27–196-fold increases in KM for p-nitrophenyl-β-D-glucoside, and significant decreases in the kcat, up to 5000-fold. At the C6OH binding residues, mutations of E777 decreased the kcat/KM by over 60,000-fold, while R786 mutants increased both the KM (40-fold) and kcat (2–4-fold). The crystal structures of R786A and R786K suggested a larger entrance to the active site could facilitate their faster rates. ONIOM binding energy calculations identified D452, H507, E777, and R786, along with the catalytic residues E441 and D593, as strong electrostatic contributors to glucose binding with predicted interaction energies > 15 kcal mol−1, consistent with the effects of the D452, H507, E777 and R786 mutations on enzyme kinetics. The relative importance of GH116 active site residues in substrate binding and catalysis identified in this work improves the prospects for the design of inhibitors for GBA2 and the engineering of GH116 enzymes for hydrolytic and synthetic applications.

Transcriptomic characterization of tissues from patients and subsequent pathway analyses reveal biological pathways that are implicated in spastic ataxia

BACKGROUND

Spastic ataxias (SAs) encompass a group of rare and severe neurodegenerative diseases, characterized by an overlap between ataxia and spastic paraplegia clinical features. They have been associated with pathogenic variants in a number of genes, including GBA2. This gene codes for the non-lysososomal β-glucosylceramidase, which is involved in sphingolipid metabolism through its catalytic role in the degradation of glucosylceramide. However, the mechanism by which GBA2 variants lead to the development of SA is still unclear.

METHODS

In this work, we perform next-generation RNA-sequencing (RNA-seq), in an attempt to discover differentially expressed genes (DEGs) in lymphoblastoid, fibroblast cell lines and induced pluripotent stem cell-derived neurons derived from patients with SA, homozygous for the GBA2 c.1780G > C missense variant. We further exploit DEGs in pathway analyses in order to elucidate candidate molecular mechanisms that are implicated in the development of the GBA2 gene-associated SA.

RESULTS

Our data reveal a total of 5217 genes with significantly altered expression between patient and control tested tissues. Furthermore, the most significant extracted pathways are presented and discussed for their possible role in the pathogenesis of the disease. Among them are the oxidative stress, neuroinflammation, sphingolipid signaling and metabolism, PI3K-Akt and MAPK signaling pathways.

CONCLUSIONS

Overall, our work examines for the first time the transcriptome profiles of GBA2-associated SA patients and suggests pathways and pathway synergies that could possibly have a role in SA pathogenesis. Lastly, it provides a list of DEGs and pathways that could be further validated towards the discovery of disease biomarkers.

Effect of Expression of Human Glucosylceramidase 2 Isoforms on Lipid Profiles in COS-7 Cells

Glucosylceramide (GlcCer) is a major membrane lipid and the precursor of gangliosides. GlcCer is mainly degraded by two enzymes, lysosomal acid β-glucosidase (GBA) and nonlysosomal β-glucosidase (GBA2), which may have different isoforms because of alternative splicing. To understand which GBA2 isoforms are active and how they affect glycosphingolipid levels in cells, we expressed nine human GBA2 isoforms in COS-7 cells, confirmed their expression by qRT-PCR and Western blotting, and assayed their activity to hydrolyze 4-methylumbelliferyl-β-D-glucopyranoside (4MUG) in cell extracts. Human GBA2 isoform 1 showed high activity, while the other isoforms had activity similar to the background. Comparison of sphingolipid levels by ultra-high resolution/accurate mass spectrometry (UHRAMS) analysis showed that isoform 1 overexpression increased ceramide and decreased hexosylceramide levels. Comparison of ratios of glucosylceramides to the corresponding ceramides in the extracts indicated that GBA2 isoform 1 has broad specificity for the lipid component of glucosylceramide, suggesting that only one GBA2 isoform 1 is active and affects sphingolipid levels in the cell. Our study provides new insights into how increased breakdown of GlcCer affects cellular lipid metabolic networks.

Analyzing Gene Expression Profiles from Ataxia and Spasticity Phenotypes to Reveal Spastic Ataxia Related Pathways

Spastic ataxia (SA) is a group of rare neurodegenerative diseases, characterized by mixed features of generalized ataxia and spasticity. The pathogenetic mechanisms that drive the development of the majority of these diseases remain unclear, although a number of studies have highlighted the involvement of mitochondrial and lipid metabolism, as well as calcium signaling. Our group has previously published the GBA2 c.1780G > C (p.Asp594His) missense variant as the cause of spastic ataxia in a Cypriot consanguineous family, and more recently the biochemical characterization of this variant in patients’ lymphoblastoid cell lines. GBA2 is a crucial enzyme of sphingolipid metabolism. However, it is unknown if GBA2 has additional functions and therefore additional pathways may be involved in the disease development. The current study introduces bioinformatics approaches to better understand the pathogenetic mechanisms of the disease. We analyzed publicly available human gene expression datasets of diseases presented with ‘ataxia’ or ‘spasticity’ in their clinical phenotype and we performed pathway analysis in order to: (a) search for candidate perturbed pathways of SA; and (b) evaluate the role of sphingolipid signaling pathway and sphingolipid metabolism in the disease development, through the identification of differentially expressed genes in patients compared to controls. Our results demonstrate consistent differential expression of genes that participate in the sphingolipid pathways and highlight alterations in the pathway level that might be associated with the disease phenotype. Through enrichment analysis, we discuss additional pathways that are connected to sphingolipid pathways, such as PI3K-Akt signaling, MAPK signaling, calcium signaling, and lipid and carbohydrate metabolism as the most enriched for ataxia and spasticity phenotypes.