《Genes》：Integrative Identification of Candidate Protein Targets and Compounds for Dystonia Using Mendelian Randomization, Single-Cell RNA Sequencing, and Network Pharmacology
编辑推荐：
背景：肌张力障碍是一种发病机制不明的严重神经系统疾病。当前治疗选择有限，难以阻止疾病进展，亟需开发新的靶向治疗药物。方法：研究人员开展了蛋白质组范围孟德尔随机化（MR）研究及敏感性分析，评估肌张力障碍与蛋白质之间的因果关系；对肌张力障碍相关蛋白进行基因本体（G
背景：肌张力障碍是一种发病机制不明的严重神经系统疾病。当前治疗选择有限，难以阻止疾病进展，亟需开发新的靶向治疗药物。方法：研究人员开展了蛋白质组范围孟德尔随机化（MR）研究及敏感性分析，评估肌张力障碍与蛋白质之间的因果关系；对肌张力障碍相关蛋白进行基因本体（GO）与京都基因与基因组百科全书（KEGG）富集分析；构建蛋白互作（PPI）网络，并利用单细胞测序数据识别核心基因在特定脑神经元中的表达；对核心基因进行药物富集分析，并采用网络药理学与分子对接方法识别肌张力障碍的潜在药物。结果：研究识别出51个蛋白质与肌张力障碍风险存在遗传预测关联，符合潜在因果效应。GO与KEGG富集分析显示这些蛋白参与转化生长因子-β刺激的细胞应答及细胞因子-细胞因子受体相互作用。PPI网络在51个相关蛋白中表现出21个群落关系。单细胞RNA脑簇特异性注释显示肿瘤坏死因子（TNF）在小胶质细胞中高表达。药物富集分析识别出5个中药单体（芍药苷、青蒿琥酯、人参皂苷Rh1、补骨脂素、二水槲皮素）作为分子对接候选。其中芍药苷-TNF、二水槲皮素-TNF、青蒿琥酯-TNF结合能最高（?9.1 kcal/mol）。结论：分子对接分析提示芍药苷、二水槲皮素、青蒿琥酯等中药单体可作为未来药物开发的候选。
研究背景与立项依据
肌张力障碍（dystonia）是以不自主肌肉收缩导致异常运动与姿势为特征的神经系统疾病，在运动障碍中患病率仅次于帕金森病与特发性震颤。其病理生理机制复杂，现有药物（抗胆碱药、苯二氮?类、肌松剂、左旋多巴）疗效有限且不良反应明显，新药发现成功率持续下降。循环蛋白质作为可成药基因产物，是潜在治疗靶点；蛋白质定量性状位点（pQTL）驱动的孟德尔随机化（MR）可规避混杂与反向因果，适合靶点发现。同时，网络药理学（NP）与分子对接为中药单体（TCM monomers）的靶向筛选提供新范式。鉴于肌张力障碍的MR与中药单体研究稀缺，研究人员整合蛋白质组范围MR、单细胞RNA测序（scRNA-seq）、PPI网络、NP与分子对接，系统挖掘候选蛋白靶点及中药单体，论文发表于《Genes》。
主要技术方法概括
研究人员采用欧洲裔人群汇总数据：血浆pQTL来自UK Biobank-PPP（54,219人，2,923蛋白，筛选cis-pQTL，F>10，LD r2<0.001，p<5×10??）；肌张力障碍结局来自Finngen R10（ICD-10 G24，1,647例/399,663对照）。以TwoSampleMR包行两样本MR（IVW、MR-Egger、加权中位数等），辅以Cochran Q、MR-PRESSO、Radial MR、Steiger过滤与BH-FDR校正；clusterProfiler行GO/KEGG富集；STRING与Cytoscape（cytoHubba）建PPI并提核心基因；PanglaoDB取人脑scRNA-seq两数据集验表达；DSigDB做药物富集，Autodock Vina行分子对接（结合能
研究结果
3.1 MR估计与敏感性分析
对941个合格血浆蛋白行MR，IVW初筛51个与肌张力障碍名义相关（p<0.05，FDR后q>0.05）。每SD蛋白水平升高，OR介于MINK1的0.337至SPRY2的3.943，TNF为2.051。25个蛋白（≥3 IVs）多法一致；敏感性分析（异质性、水平多效性、方向性）支持可靠性；26个蛋白因IVs<3仅用IVW。
3.2 功能与通路富集
51蛋白GO得117项（FDR<0.05），生物过程（BP）富集于急性期反应、转化生长因子-β应答、急性炎症应答；分子功能（MF）涉细胞外基质结构、糖胺聚糖结合；细胞组分（CC）归胶原外基质、质膜外侧。KEGG最显著为细胞因子-细胞因子受体相互作用、氮代谢、脂肪细胞因子信号。
3.3 PPI网络与核心基因
51蛋白曼哈顿图标出CD101、ITGB6、IL6R、TNF、F2、ACAN等；PPI网络含21个群落；cytoHubba提10核心基因，TNF、F2、ACAN为关键。
3.4 单细胞RNA注释
两脑scRNA-seq数据集显示TNF于小胶质细胞高表达；IFNGR1脑特异且在多免疫细胞表达；F2、ACAN、MBL12、IL6R、CD40、POSTN、COL9A1、SERPINA1低或无声。
3.5 药物富集与分子对接
10核心基因药物富集得617药（FDR<0.05），含5中药单体（芍药苷、青蒿琥酯、人参皂苷Rh1、补骨脂素、二水槲皮素）。对接TNF与SERPINA1，结合能?6.7至?9.1 kcal/mol；芍药苷-TNF、二水槲皮素-TNF、青蒿琥酯-TNF最优（?9.1 kcal/mol），人参皂苷Rh1-TNF（?6.8）、补骨脂素-TNF（?6.7）、二水槲皮素-SERPINA1（?6.7）较弱。
讨论部分总结
研究人员指出本研究首次将蛋白质组MR、scRNA-seq、NP整合用于肌张力障碍靶点发现；Steiger与cis-pQTL控方向偏倚。富集指向TGF-β应答与细胞因子通路，PPI核心TNF于小胶质高表，契合中枢神经炎症机制（TNF经TNFR1促M1极化、NF-κB自放大环路）。芍药苷、青蒿琥酯、二水槲皮素均有抑制微胶质TNF报道，与MR“高TNF增风险”方向一致。局限：仅欧洲裔、未用Priority index、FDR后无显著蛋白（故51关联为初步）、对接仅为计算机筛查非亲和力定量、缺动物验证、生物制剂跨血脑屏障（BBB）难。意义在于提供“因果蛋白→小胶质TNF→中药单体”可药靶框架。
结论部分翻译
研究识别出51个蛋白质与肌张力障碍风险存在名义显著遗传预测关联。分子对接分析提示芍药苷、二水槲皮素、青蒿琥酯等中药来源单体可与TNF蛋白相互作用，这些化合物是需进一步实验验证的候选剂，可为肌张力障碍治疗靶点开发提供新线索；结果凸显中药单体作为肌张力障碍候选靶向干预的潜力。
