溶酶体鞘脂降解通路中GCase1抑制诱导的双稳态与滞后：从Gaucher disease到Parkinson’s disease及α-synuclein纤维化

《Mathematical Biosciences》：Parkinson’s disease, Gaucher disease, and bistability in the lysosomal degradation pathway

　　

鞘脂（sphingolipids）是细胞膜的内在组分，其被转运至溶酶体（lysosomes）降解为神经酰胺（ceramide），后者在高尔基体（Golgi）中用于合成新的鞘脂。尽管鞘脂代谢非常复杂，但降解通路相对简单，可将鞘脂转化为神经酰胺和鞘氨醇（sphingosine）。该通路中酶基因的多态性与多种溶酶体贮积症（lysosomal storage disorders, LSDs）相关，其中一种酶，葡萄糖脑苷脂酶（glucosylceramidase, GCase1），与帕金森病（Parkinson’s disease, PD）、路易体痴呆（Lewy-body dementia, LBD）以及戈谢病（Gaucher disease）相关。研究人员构建了鞘脂降解通路的数学模型，并表明在GCase1的V_max值较大范围内，该通路表现出双稳态（bi-stability），即两个稳态且二者之间存在滞后（hysteresis）。此外，已知由单体（monomers）和寡聚体（oligomers）形成突触核蛋白纤维（synuclein fibrils）的过程受降解通路影响，并反过来抑制GCase1。研究人员表明存在两个稳态，一个具有低纤维生成，一个具有高纤维生成，并伴有滞后。存在两个稳态且具有滞后这一事实，为药物发现和给药既带来挑战也带来机会。研究人员的工作对衰老（aging）具有意义，因为已知GCase1基因的表达水平随年龄下降。

论文解读：

研究背景：鞘脂（sphingolipids）是细胞膜内在组分，参与膜通透性、细胞信号和膜流动性调控。鞘脂在内质网和高尔基体（Golgi）由神经酰胺（ceramide）合成，转运至细胞膜，再进入溶酶体（lysosomes）和内质网降解回神经酰胺。溶酶体降解通路中，exo-α-sialidase（SIAL）、beta-galactosidase（GLB2）、glucosylceramidase（GCase1）和acid ceramidase（ACR）依次将鞘脂转化为神经酰胺和鞘氨醇（sphingosine），再循环至内质网和高尔基体。该通路酶缺陷或基因多态性与多种溶酶体贮积症相关，如SIAL与神经退行性疾病、GLB2与Tay-Sachs、GCase1与Gaucher disease、ACR与Farber disease。GCase1基因GBA1涉及帕金森病（Parkinson’s disease, PD）和路易体痴呆（dementia with Lewy bodies, DLB），患者血液GBA1 mRNA降低；GCase1活性随年龄下降，可能解释PD多为晚发。小鼠证据提示衰老导致的鞘脂代谢扰动，包括GCase1降低，可能加速PD病理过程。PD特征是黑质致密部（substantia nigra pars compacta, SNc）多巴胺能神经元退化，线粒体氧化应激和溶酶体功能障碍与神经元退化和Parkinsonism相关。GBA1突变是Parkinsonism主要遗传风险因素，也导致Gaucher disease，使GCase1底物GluCer积累并升高血浆GluSph。临床已描述Gaucher disease与PD联系，但机制不清。α-synuclein积聚形成Lewy bodies。Mazzulli et al.提出高GluCer激活α-synuclein组装，组装体反过来抑制GCase1，构成“双向致病环路”。因此，研究人员构建数学模型，以探究该降解通路中GCase1抑制是否导致双稳态和滞后，并评估其对药物设计、给药和衰老的意义。

关键技术方法：研究人员建立溶酶体鞘脂降解通路的常微分方程质量平衡模型；变量为底物浓度（μM），多数酶采用Michaelis-Menten动力学，K_m取自文献，V_max按底物浓度与测量范围一致来选取；GCase1速度含GluSph抑制项和synuclein寡聚体抑制项，ACR2速度对GluCer呈S型，synuclein速度含GluCer刺激项；线性速率常数作为复杂转运与代谢过程的占位。模拟通过改变GCase1 V_max、GM3输入和Cer输入，并随机扰动除GCase1外所有酶V_max 1000次检验稳健性。实验数据采用Dekker et al.和Rolfs et al.中Gaucher患者与健康对照的GluCer、GluSph数据重绘分析。

研究结果：

3. Inhibition of GCase1 by glucosylsphingosine：通过关闭红色相互作用、开启GluSph对GCase1的绿色抑制，模拟标准GM3和Cer输入均为50 μM/hr。结果在GCase1 V_max约100至450 μM/hr范围内，GluCer和GluSph均有两个稳态和滞后；V_max低于100 μM/hr仅高稳态，高于450 μM/hr仅低稳态。降低输入至40 μM/hr，双稳态区域左移并缩小至约68至255 μM/hr；降低至20 μM/hr，高稳态消失，仅低稳态。随机扰动其他酶V_max 1000次仍见双稳态和滞后。Dekker et al.和Rolfs et al.数据中Gaucher患者与对照的GluSph/GluCer呈分离簇，支持两个分离稳态。

4. Synuclein：通过关闭绿色抑制、开启红色相互作用，模拟标准输入。结果在GCase1 V_max约177至280 μM/hr范围内，GluCer和synuclein fibrils均有两个稳态和滞后；高V_max仅低稳态，低V_max仅高稳态。输入升高时转变点右移且双稳态区域增大；输入降低时区域减小；在90、100、110 μM/hr输入下均存在高、低稳态及二者共存的V_max范围。极低输入时仅一个正常稳态。双稳态和滞后在独立同时变化所有V_max参数时稳健。

讨论部分总结：鞘脂代谢难研究，涉及复杂生化及高尔基体、细胞膜、内体、内质网和溶酶体间转运。研究人员聚焦降解通路，展示两种GCase1抑制均可导致双稳态和滞后。酶K_m来自文献，V_max代表基因表达水平且随时间与个体变异大，其他参数为估计，因此双稳态发生的GCase1 V_max范围系统依赖。同时包含红、绿箭头也产生双稳态；ACR2 V_max变化亦可产生双稳态。Gaucher disease自1882年已知，患者表型广泛；1965、1966年发现其均存在GCase1有害变异。非神经元型Gaucher disease常表现为器官肿大、血细胞减少和骨病，GBA1突变使巨噬细胞积累GluCer等糖脂并形成“Gaucher cells”，浸润器官引起器官肿大、炎症等。神经元型Gaucher disease可出现脑干功能障碍、肝脾肿大和双侧第6神经麻痹等，GluCer在溶酶体异常积累导致内质网钙释放过多。PD与Gaucher disease联系较新：Sidranski于2005年提示GBA1突变与synucleinopathy可能存在临床联系，Mazzulli et al.提出机制；现认为GBA1多态性是PD主要风险因素。高GluCer刺激synuclein单体向寡聚体转化，寡聚体抑制GCase1使GluCer进一步升高，最终形成synuclein纤维；纤维具神经毒性，与SNc多巴胺神经元死亡相关。因此低GluCer稳态称“正常”，高GluCer稳态称“疾病”。双稳态在代谢模型中较少见，Ha、Satin和Sherman的糖尿病葡萄糖-胰岛素动力学模型存在稳定“健康”和“病理”稳态，糖酵解通量开关也是双稳态与滞后例子。正反馈可产生不稳定、双稳态和滞后。药物挑战与机会：若健康人处于低稳态，V_max降至200仍可稳定，低于177才切换至病理稳态；一旦进入病理稳态，药物需将V_max升至280以上方可使病理稳态消失并恢复。低Cer输入时仅正常稳态，提示抑制ceramide合成可能有用；联合限制ceramide输入和增加GCase1 V_max可能有益。GCase1活性随年龄下降，PD和DLB患者血液GBA1 mRNA降低，老年患者可能需持续药物保持GCase1 V_max高。若GCase1 V_max下降与SNc神经元死亡相关，提高GBA1表达和/或GCase1活性的药物可能改善或延缓PD，也可能对痴呆如此。

研究结论翻译：最近发现GCase1活性随年龄下降，使这一情形更加复杂和有趣，因为GCase1 V_max低值可能并非突然的特发性下降，而是衰老导致的下降。血液GBA1 mRNA在PD患者和Lewy-body dementia患者中降低，这使老年患者很可能需要持续药物治疗以维持GCase1 V_max处于高水平。PD的近端原因是SNc多巴胺神经元死亡。如果GCase1 V_max下降与SNc神经元死亡相关，那么通过药物提高GBA1表达水平和/或GCase1活性水平，可能改善或至少延缓PD，且对痴呆也可能如此。