《Brain Research Bulletin》:Integrating target prediction and transcriptomic analysis identifies the AKR1C2-KRT8 axis in rotenone-induced dopaminergic neuron injury
编辑推荐:
整合生物信息学与机器学习在鱼藤酮诱导的神经毒性中鉴定出一个六基因核心网络。研究人员发现AKR1C2在原代多巴胺能神经元中正向调控KRT8表达。AKR1C2敲低可挽救线粒体功能障碍、α-突触核蛋白(α-synuclein)聚集及电生理缺陷。KRT8重新表达逆转了
整合生物信息学与机器学习在鱼藤酮诱导的神经毒性中鉴定出一个六基因核心网络。研究人员发现AKR1C2在原代多巴胺能神经元中正向调控KRT8表达。AKR1C2敲低可挽救线粒体功能障碍、α-突触核蛋白(α-synuclein)聚集及电生理缺陷。KRT8重新表达逆转了AKR1C2敲低的保护效应,证实了二者的因果轴关系。
研究背景与立项依据
帕金森病(Parkinson's disease, PD)是第二常见的神经退行性疾病,核心病理为中脑黑质致密部多巴胺能神经元进行性不可逆丢失及以α-突触核蛋白(α-synuclein)为主要成分的路易小体沉积。现有多巴胺替代疗法仅缓解运动症状,无法延缓病程,瓶颈在于PD病因复杂,遗传易感、环境毒素暴露与衰老长期交互,下游分子网络未系统阐明。环境毒素如鱼藤酮(rotenone,线粒体复合体Ⅰ抑制剂)可模拟PD关键病理,但对鱼藤酮暴露后选择性损伤多巴胺能神经元的具体下游靶点与信号轴仍缺乏系统揭示。传统单组学或候选基因策略假阳性高、难以捕获致病网络,因此研究人员采用计算生物学与实验验证相整合的范式,从GEO(Gene Expression Omnibus)数据库PD患者转录组与鱼藤酮已知靶点交集出发,结合单细胞转录组(single-cell RNA sequencing, scRNA-seq)与体外原代神经元功能实验,系统鉴定并验证鱼藤酮-PD核心分子轴。该研究发表于《Brain Research Bulletin》。
主要关键技术方法
研究人员下载GEO数据库PD患者芯片数据集GSE7621、GSE20141、GSE49036及单细胞数据集GSE276684,并集ChEMBL、SwissTargetPrediction、SEA三库预测鱼藤酮靶点;经ComBat算法批次校正后,limma包筛选差异表达基因(DEGs,|logFC|>0.585,adj.P.Val<0.05);取鱼藤酮靶点与DEGs交集,用Random Forest、SVM(Support Vector Machine)、glmBoost等多算法机器学习十折交叉验证筛核心 Biomarker,构建Logistic回归列线图;CIBERSORT评估22类免疫细胞浸润;Seurat+SingleR注释单细胞群,scTenifoldKnk做虚拟敲除(virtual knockout)网络扰动;clusterProfiler做GO/KEGG/Reactome富集;CB-Dock 2做鱼藤酮与核心蛋白分子对接(ΔG5.0 kcal/mol为阈值);原代SD大鼠腹侧中脑多巴胺能神经元培养,10 nM鱼藤酮24 h造模,慢病毒介导AKR1C2/KRT8敲低或过表达及 Rescue 实验,CCK-8、Complex Ⅰ活性、ATP、qRT-PCR、Western blot、免疫荧光双标、全细胞膜片钳测动作电位(action potential, AP)验证。
研究结果
3.1 PD患者基因表达与免疫微环境变化
批次校正后三数据集PCA聚良好;DEGs火山图与热图显示PD显著紊乱。CIBERSORT表明不成熟B细胞、浆细胞样树突细胞、MDSC、中性粒细胞、巨噬细胞、效应记忆CD8 T细胞、NK等十余类免疫细胞在PD显著升高,提示免疫微环境激活。
3.2 靶点预测与机器学习筛选PD核心 Biomarker及诊断模型
三库并集得133鱼藤酮靶点,与DEGs交得10个交集基因;glmBoost+SVM模型AUC=0.732最优,单基因ROC筛出6个核心预测基因:AKR1C2、AKR1C3、CES1、CTSS、DRD2、HSPA1A;基于6基因建Logistic列线图评估PD风险。
3.3 核心 Biomarker 表达特征、免疫关联与功能通路
AKR1C2/3、CES1、CTSS、HSPA1A在PD上调,DRD2下调;AKR1C3、CTSS与多数免疫细胞正相关,DRD2多负相关。KEGG富集化学致癌-DNA加合物、花生四烯酸代谢、抗原加工提呈、凋亡、溶酶体、内质网蛋白加工等;Reactome富集胆汁酸代谢、角质化(Keratinization)、角化包膜形成;GO富集C21-甾体激素代谢、前列腺素代谢、醇代谢、先天免疫应答等。
3.4 核心靶点与鱼藤酮分子对接验证
6蛋白与鱼藤酮均形成稳定构象,Vina评分均7.7 kcal/mol,AKR1C3、AKR1C2、CTSS、DRD2结合尤紧密。
3.5 单细胞图谱揭示PD细胞状态改变并聚焦AKR1C2-KRT8调控轴
GSE276684分神经元与星形胶质细胞两大群;PD组神经元左区密度增、星形胶质细胞右下区富集。气泡图示AKR1C2在神经元最富集。scTenifoldKnk虚拟敲除AKR1C2仅扰约2%基因,KRT8为最显著下游之一;富集通路由角化包膜形成、IL-4/IL-13信号、GPCR(G protein-coupled receptor)配体结合、角质化主导。
3.6 鱼藤酮降低原代多巴胺能神经元活力、损线粒体功能并改变核心基因表达
0–400 nM梯度CCK-8选10 nM为造模浓度;10 nM 24 h显著降Complex Ⅰ活性与ATP;qRT-PCR示AKR1C2/3、CES1、CTSS、HSPA1A、KRT8 mRNA上调,DRD2下调;免疫荧光确认各核心蛋白与酪氨酸羟化酶(Tyrosine Hydroxylase, TH)共定位,AKR1C2与KRT8同细胞共表达。
3.7 AKR1C2-KRT8轴调控鱼藤酮处理神经元蛋白表达、线粒体功能与兴奋性
Western blot:鱼藤酮升AKR1C2/3、CES1、CTSS、HSPA1A、KRT8、α-synuclein,降DRD2、TH。OE-AKR1C2升KRT8/α-synuclein、降TH、损活力/Complex Ⅰ/ATP。shRNA-AKR1C2反向且降α-synuclein、升TH、改善功能;Forward Rescue(shRNA-AKR1C2+OE-KRT8)逆转保护效应。Reverse Rescue(OE-AKR1C2+shRNA-KRT8)回救过表达表型。膜片钳:鱼藤酮降AP幅度与频率;shRNA-AKR1C2恢幅度不恢频率,Rescue取消幅度恢复;OE-AKR1C2降频率不降幅度,Reverse Rescue恢频率。证实AKR1C2正向调KRT8,轴控线粒体、α-synuclein、TH及电生理。
讨论部分总结
研究人员指出,环境毒素暴露是散发性PD主因之一,鱼藤酮虽为经典模型,下游网络未明。本研究从多库靶点预测→机器学习6基因核心模块→分子对接直接结合证据→PD患者bulk与scRNA均示神经元富集→原代神经元证实AKR1C2-KRT8因果轴。AKR1C2属醛酮还原酶家族,参与神经甾体与神经活性醛代谢,其上调或为氧化应激适应反应,但驱动KRT8(中间丝蛋白)异常表达致细胞骨架重排,干扰ankyrin/spectrin对Na+/K+通道在轴丘锚定及微管依赖转运,从而降AP幅度/频率;同时KRT8异常促α-synuclein聚集。该轴连接环境暴露—代谢酶—细胞骨架—电生理—蛋白聚集,补入PD环境病因分子链条。局限:仅急性体外模型,缺体内动物与更多人脑验证,KRT8为保护还是致病需直接操纵明确。
结论部分翻译
从环境暴露病因学视角,本研究通过多维数据整合与体外实验验证,系统勾勒鱼藤酮暴露后多巴胺能神经元含多关键基因的分子网络,其中代谢酶AKR1C2对骨架蛋白KRT8的直接调控获实验确证。功能检测进一步表明该轴关键节点AKR1C2表达水平直接影响神经元存活与电生理活动,AKR1C2敲低部分挽救鱼藤酮诱导的功能损伤。工作揭示AKR1C2-KRT8新调控轴,为理解特定环境毒素如何影响神经元状态提供具体分子线索;所采用计算与实验整合范式亦为复杂病因下致病因子系统鉴定与验证提供方法学参考。