Investigating the Impact of the Parkinson’s-Associated GBA1 E326K Mutation on β-Glucocerebrosidase Dimerization and Interactome Dynamics Through an In Silico Approach

Heterozygous mutations or genetic variants in the GBA1 gene, which encodes for the β-glucocerebrosidase (GCase), a lysosomal hydrolase enzyme, may increase the risk of Parkinson’s disease (PD) onset. The heterozygous E326K form is one of the most common genetic risk factors for PD worldwide, but, to date, the underlying molecular mechanisms remain unclear. Here, we investigate the effect of the E326K on the structure, stability, dimerization process, and interaction mode with some proteins of the interactome of GCase using multiple molecular dynamics (MD) simulations at pH 5.5 and pH 7.0 to mimic the lysosomal and endoplasmic reticulum environments, respectively. The analysis of the MD trajectories highlights that the E326K mutation did not significantly alter the structural conformation of the catalytic dyad but significantly makes the structure of the dimeric complexes unstable, especially at lysosomal pH, potentially impacting the organization of the quaternary structure. Furthermore, the E326K mutation significantly impacts protein interactions by altering the binding mode with the activator Saposin C (SapC), reducing the binding affinity with the inhibitor α-Synuclein (α-Syn), and increasing the affinity for the Lysosomal integral membrane protein-2 (LIMP-2) transporter.