《BioMed Research International》：Integrative Transcriptomic and Network Analyses Reveal Common Molecular Signatures and Turmeric-Derived Common Therapeutic Candidates for Neurological Diseases
编辑推荐：
摘要
背景：神经系统疾病（Neurological diseases, NDs），包括阿尔茨海默病（Alzheimer′s disease, AD）、帕金森病（Parkinson′s disease, PD）、卒中（stroke）和脑肿瘤，给全球健康造成了沉重
摘要
背景：神经系统疾病（Neurological diseases, NDs），包括阿尔茨海默病（Alzheimer′s disease, AD）、帕金森病（Parkinson′s disease, PD）、卒中（stroke）和脑肿瘤，给全球健康造成了沉重负担。尽管这些疾病具有不同的临床表现和病因学，但越来越多的证据表明它们共享某些分子机制。然而，其共同的分子特征以及姜黄来源生物活性化合物的治疗潜力尚未被系统研究。
目的：本研究旨在利用整合生物信息学框架，识别AD、PD、卒中和胶质母细胞瘤（glioblastoma, GBM）之间的共享分子机制和关键枢纽基因（hub genes），并探究姜黄来源化合物作为潜在的多靶点共同治疗候选物。
材料与方法：从基因表达综合数据库（Gene Expression Omnibus, GEO）获取基因表达数据集（GSE22491、GSE5281、GSE122709和GSE4290）。通过整合生物信息学分析识别共同宿主差异表达基因（common host differentially expressed genes, cHDEGs）。随后构建蛋白质-蛋白质相互作用（protein–protein interaction, PPI）网络以识别共同枢纽基因（common hub genes, cHubGs），并进行功能富集及转录因子（transcription factor, TF）和miRNA介导的调控网络分析。进一步通过分子对接和500纳秒分子动力学模拟（molecular dynamics simulations, MDS）研究所识别cHubGs与姜黄来源化合物之间的相互作用。
结果：共识别出81个cHDEGs和4个cHubGs（CCNB1、CDK1、CYCS和STAT3）。功能富集和调控网络分析揭示了共享的生物学过程和信号通路，以及与ND发病机制潜在相关的关键TFs和miRNAs。在77个通过毒性筛选的姜黄来源化合物中，去甲氧基姜黄素（demethoxycurcumin）、环姜黄素（cyclocurcumin）和(1E,4E)-1-(4-羟基-3-甲氧基苯基)-5-(4-羟基苯基)-1,4-戊二烯-3-酮表现出最强的预测结合亲和力。所选化合物-靶点复合物在MDS过程中也表现出稳定的行为。
结论：本研究为AD、PD、卒中和GBM之间的共享分子特征提供了见解，并将姜黄来源化合物优先列为潜在的多靶点治疗候选物。这些计算发现需要通过独立数据集和实验研究进行进一步验证。
神经系统疾病共享分子机制与姜黄素多靶点治疗潜力的整合研究解读
一、研究背景与科学问题
神经系统疾病（Neurological diseases, NDs）涵盖所有影响神经系统的病理状态，包括阿尔茨海默病（Alzheimer′s disease, AD）、帕金森病（Parkinson′s disease, PD）、卒中和脑肿瘤等。2021年全球数据显示，神经系统疾病已成为第二大死亡原因和伤残调整生命年（disability-adjusted life years, DALYs）的首要来源。AD患者全球约5500万，预计2050年将增至1.52亿；PD患者2019年约851万；GBM患者5年生存率仅略高于5%。尽管这些疾病具有截然不同的临床表现和病因学特征，但近年研究日益表明它们可能共享部分趋同的分子和细胞机制，如异常蛋白聚集、神经血管耦合失调、免疫相关基因表达异常等。然而，此前研究多聚焦于单一疾病或有限疾病组合，尚未系统性地跨多种神经系统疾病识别共同的转录组特征并关联潜在治疗候选物。此外，天然产物因其结构多样性和多效生物学活性，成为多靶点治疗药物的重要来源，其中姜黄（Curcuma longa）具有抗氧化、抗炎和神经保护等广泛药理特性。因此，研究人员开展了本项整合性研究。
二、研究内容与主要结论
研究人员采用整合转录组学和网络药理学框架，系统研究了AD、PD、缺血性卒中和GBM四种神经系统疾病的共享分子特征，并评估了姜黄来源化合物作为潜在多靶点治疗候选物的可行性。研究识别出81个共同宿主差异表达基因（cHDEGs）和4个共同枢纽基因（cHubGs）：CCNB1、CDK1、CYCS和STAT3。功能富集分析揭示了低氧反应、蛋白磷酸化、脂质代谢紊乱等共享生物学过程；调控网络分析鉴定了E2F1、SP1等关键转录因子及hsa-miR-92a-3p、hsa-miR-21-5p等关键miRNAs。分子对接和500纳秒分子动力学模拟筛选出去甲氧基姜黄素、环姜黄素和(1E,4E)-1-(4-羟基-3-甲氧基苯基)-5-(4-羟基苯基)-1,4-戊二烯-3-酮三种具有最强结合亲和力和稳定复合物构象的候选化合物。这些发现为跨疾病共享分子机制的深入理解及多靶点治疗策略的开发提供了理论基础。
三、关键技术与方法
研究人员从基因表达综合数据库（Gene Expression Omnibus, GEO）获取了四个独立数据集：GSE22491（PD，10例病例/8例对照）、GSE5281（AD，87例病例/74例对照）、GSE122709（缺血性卒中，10例病例/5例对照）和GSE4290（脑肿瘤，157例病例/23例对照）。主要技术方法包括：利用GEO2R和limma工具进行差异表达基因（DEGs）鉴定；通过STRING数据库和Cytoscape的cytoHubba插件（8种算法）构建PPI网络并识别cHubGs；利用TRRUST和miRTarBase数据库经ShinyGO工具构建TF和miRNA调控网络；通过GO和KEGG进行功能富集分析；采用SwissADME和pkCSM评估药物相似性和ADMET特性；使用ProTox 3.0进行毒性预测；利用PyRx/AutoDock Vina进行分子对接；使用Gromacs 2023（CHARMM36力场）进行500纳秒分子动力学模拟。
四、研究结果
4.1 神经系统疾病DEGs的识别：研究人员从四个数据集中分别鉴定出9751、5342、11048和10373个DEGs，通过Venn图分析获得22143个不同DEGs，其中0.366%为共同基因，最终识别出81个cHDEGs。
4.2 cHDEGs的PPI分析：通过STRING数据库构建PPI网络，应用8种cytoHubba算法（closeness、radiality、betweenness、stress、degree、EPC、MNC、MCC），经UpSet图分析，4个基因（CCNB1、CDK1、CYCS、STAT3）被所有算法一致识别为cHubGs。
4.3 cHDEGs的功能和通路富集分析：GO富集分析显示，生物学过程（BP）主要富集于蛋白磷酸化、低氧反应、肽基丝氨酸磷酸化等；分子功能（MF）主要涉及蛋白激酶活性、ATP结合、核苷酸结合等；细胞组分（CC）主要涉及细胞质、细胞溶质、线粒体、cyclin B1–CDK1复合物、凋亡体等。KEGG通路分析显示富集于癌症通路、乙型肝炎、FoxO信号通路、p53信号通路、病毒致癌、凋亡等多物种等通路。
4.4 cHDEGs与其调控因子的相互作用分析：TF调控网络识别出SP1、STAT3、MYC、TP53、BRCA1、E2F1、YBX1、PTTG1、NR4A1和E2F3等关键转录因子，其中E2F1最为显著并调控包括CDK1和CCNB1在内的9个cHDEGs。miRNA调控网络识别出hsa-miR-92a-3p（关联17个基因，含3个hub基因）、hsa-miR-21-5p（最显著，关联12个基因，含STAT3）等关键miRNAs。
4.5 药物相似性、ADMET和毒性分析：从IMPPAT数据库收集164种姜黄来源化合物，经Lipinski五规则、ADMET预测和毒性筛选后，77种化合物保留用于后续分析。
4.6 分子对接分析：针对CDK1（PDB: 6GU2）、STAT3（PDB: 6QHD）、CCNB1（PDB: 4YC3）和CYCS（PDB: 5EXQ）四种靶蛋白进行分子对接。初始选定的四个复合物为CDK1+CID10469828（?8.4 kcal/mol）、STAT3+CID69879809（?8.1 kcal/mol）、CCNB1+CID5469424（?8.9 kcal/mol）和CYCS+CID69879809（?9.1 kcal/mol）。因CID69879809同时与STAT3和CYCS结合，选择亲和力更优的CYCS复合物，最终三个复合物进入MDS分析。
4.7 分子动力学模拟分析：对三个复合物进行500纳秒MDS分析。RG分析显示CYCS+CID69879809保持较紧凑结构（<1.4 nm）；RMSD分析显示CCNB1+CID5469424稳定性较差，CYCS+CID69879809高度稳定；RMSF分析显示CYCS复合物残基波动较低，CCNB1复合物波动较大；SASA分析显示CYCS复合物溶剂可及表面积显著较低（约60 nm2），表明更紧凑的构象和强蛋白-配体相互作用。PCA分析显示CYCS复合物前三个主成分占78%方差，表现出优异的聚类分离；DCCM分析揭示了各复合物的协同和反协同结构域运动；FEL分析显示CDK1复合物为单峰能量景观（单一收敛结合模式），CCNB1和CYCS复合物为双峰能量景观（两个动力学隔离的结合模式）；SSA分析显示各复合物以弯曲和卷曲/环区域为主。
五、讨论与结论
研究人员在讨论中指出，尽管姜黄素类化合物具有令人鼓舞的分子相互作用结果，但其临床转化应用需谨慎考虑生物利用度低、吸收和溶解度问题、快速代谢及药代动力学不确定性等挑战。计算预测结果不应被解读为治疗功效或临床适用性的直接证据，而应视为产生假设的工具，需要通过生化、细胞、动物、药代动力学及最终临床研究进行验证。研究识别的cHubGs——CCNB1、CDK1、CYCS和STAT3——在细胞周期调控、线粒体凋亡信号和炎症信号转导中发挥核心作用，可能代表连接神经退行性变、脑肿瘤进展和缺血性损伤的关键分子调控节点。功能富集分析表明低氧应激、异常蛋白磷酸化和脂质代谢紊乱是四种疾病共享的病理机制；TF和miRNA调控网络揭示了转录和转录后水平的多层调控。研究结论认为，本研究为AD、PD、卒中和GBM的共享分子特征提供了系统性见解，并将去甲氧基姜黄素、环姜黄素和(1E,4E)-1-(4-羟基-3-甲氧基苯基)-5-(4-羟基苯基)-1,4-戊二烯-3-酮优先列为潜在多靶点治疗候选物，但这些计算发现仍需通过独立数据集和实验研究进一步验证。
