基于网络的转录组学识别阿尔茨海默病中海马关键靶点及其被芹菜素、木犀草素和小檗碱的调节

《Scientific African》：Network-Based Transcriptomics Identifies Key Hippocampal Targets in Alzheimer’s Disease and Their Modulation by Apigenin, Luteolin, and Berberine

　　

阿尔茨海默病（AD）是一种进行性神经退行性疾病，以认知功能下降和突触功能障碍为特征，海马是最易受影响的脑区之一。在本研究中，研究人员对公开可用的海马微阵列数据集进行了基于网络的转录组学分析，以识别与AD病理相关的差异表达基因（DEGs）和失调的生物通路。功能富集分析揭示了多巴胺能突触信号、钙信号、MAPK通路、胶质细胞分化和神经炎症过程的改变。蛋白-蛋白相互作用网络分析突出显示了几个高度连接的基因，包括DRD4、PTAFR、CAMK1G、RIMS2和CACNG3，表明它们在海马分子网络中的核心作用。为了探索这些优先靶点的潜在结构-功能关系，研究人员使用三种天然黄酮类化合物——芹菜素、木犀草素和小檗碱进行了分子对接和分子动力学模拟。比较分析表明，木犀草素在多个激酶靶点上表现出相对更稳定的结合，表现为在整个模拟轨迹中结构波动减小、构象紧凑和持续的分子间相互作用。这些计算发现被定性解释，同时承认了方法学局限性和缺乏实验验证。总体而言，这种整合方法将转录组网络扰动与结构相互作用模式联系起来，并为与AD相关信号通路相关的多靶点分子相互作用提供了产生假设的见解。研究结果强调了结合基于网络的转录组学与计算机结构分析的价值，同时强调了未来需要实验研究来确认所提出的机制。

阿尔茨海默病（Alzheimer's disease, AD）是一种进行性神经退行性疾病，其病理特征包括细胞外淀粉样蛋白-β（Aβ）沉积、细胞内过度磷酸化Tau蛋白、突触丢失和持续性神经炎症激活。海马是AD早期受影响的关键脑区，但其分子网络的系统层面组织仍不完全清楚，限制了多靶点治疗策略的开发。天然黄酮类化合物如芹菜素（Apigenin）、木犀草素（Luteolin）和小檗碱（Berberine）具有神经保护、抗炎和抗氧化活性，但它们在AD海马复杂分子网络中的系统相互作用尚未得到充分阐明。本研究旨在通过网络转录组学与结构分析相结合，探索这些化合物与AD海马关键靶点的相互作用，为多靶点干预提供理论依据。论文发表在《Scientific African》。

研究人员主要使用了以下关键方法：从GEO数据库（GSE1297）获取海马微阵列数据，包含对照、早期、中度和重度AD四个阶段的样本，其中对照组9例、早期7例、中度8例、重度7例。采用limma包进行差异表达分析，筛选标准为调整后P<0.05且|log?FC|≥1。通过Enrichr和Metascape进行功能富集分析，使用蛋白-蛋白相互作用（PPI）网络确定枢纽基因，基于网络拓扑指标（度和介数中心性）优先选择了19个靶点。随后对这三种化合物与靶点进行分子对接（AutoDock 4），对结合能最优的复合物进行300 ns分子动力学模拟（CHARMM36力场，CHARMM-GUI构建系统），并计算MM/PBSA结合自由能。

**Identification and Prioritization of Dysregulated Hippocampal Genes Across Alzheimer’s Disease Stages**：通过差异表达分析，在重度AD与对照组比较中识别出显著失调的基因，包括下调的DRD4、PTAFR、MEOX2，以及上调的CAMK1G、CACNG3、RIMS2。在重度与早期AD比较中，仅DRD4显著上调，提示其可能参与晚期多巴胺能功能障碍。

**Pathway and Functional Enrichment Analysis**：KEGG通路富集显示催产素信号通路、钙信号通路、胰岛素分泌、多巴胺能突触和谷氨酸能突触等通路显著富集。GO分析显示定位、生物过程调节、刺激应答等过程受扰，神经递质转运调节、细胞突起组织、微管过程等形成密集互作簇。PPI网络分析发现PKP4、PRKCB和APC构成核心功能模块，参与细胞骨架组织和信号转导。

**Selection of Key AD-Associated Genes for Docking**：基于通路参与和功能相关性，选择了18个基因（如DRD4、RIMS2、CACNG3、CAMK1G、PTAFR、PRKCB、MAP3K9等）进行分子对接。

**Molecular Docking Results and Interaction Analysis**：木犀草素、芹菜素和小檗碱与19个靶点均形成结合，其中CAMK1G、PRKCB和MAP3K9的结合能最有利（ΔG_bind约-7.9至-8.2 kcal/mol）。通过重结晶配体对接验证了结合位点的准确性。

**Molecular Dynamics Simulation Results**：300 ns模拟显示所有复合物在50-100 ns内达到平衡。CAMK1G-木犀草素复合物的平均RMSD为0.23±0.03 nm，优于其他配体。Rg分析表明蛋白质紧密性保持稳定，氢键分析显示持续相互作用。

**Structural Dynamics and Interaction Analysis of Ligand–Protein Complexes**：溶剂可及表面积（SASA）分析显示木犀草素复合物具有更低的SASA值，表明更紧密的疏水堆积。距离分析显示木犀草素复合物具有更小的最小接触距离（0.24-0.32 nm）。残基柔性分析（RMSF）和残基ASA表明木犀草素结合使催化环和调控环区域刚性增强。氢键占有率分析表明木犀草素与关键残基形成持续氢键（如CAMK1G的Thr193、Glu204、Gln277）。

**Principal Component Analysis and Free Energy Landscape of Protein–Ligand Complexes**：主成分分析和自由能景观（FEL）显示木犀草素复合物处于单一且深的能量盆地，表明构象稳定且受限，而芹菜素和小檗碱复合物探索更宽的能量盆地。

**MM/PBSA Binding Free Energy Analysis**：木犀草素与CAMK1G、PRKCB和MAP3K9的总结合自由能分别为-9.24±3.43、-4.83±4.24和-6.28±1.69 kcal/mol，主要由范德华力和静电相互作用驱动，且极性溶剂化惩罚适中，显示出能量协同效应。小檗碱虽有时表现出高瞬时结合能，但波动大；芹菜素结合能较弱。

讨论部分指出，本研究强调了AD海马中多巴胺能、钙依赖和神经炎症信号通路的失调，DRD4、CAMK1G、CACNG3等基因的上调以及PTAFR、MEOX2的下调，与以往大规模转录组学研究一致。木犀草素在分子动力学模拟中表现出最优的稳定性和结合能，可能与其中枢神经系统生物利用度和脑渗透性有关。已有实验证据支持木犀草素改善空间记忆、降低氧化应激和神经炎症、抑制Aβ聚集和Tau过度磷酸化。本研究为网络转录组学与结构分析的整合提供了框架，但需进一步实验验证。

**研究结论**：在本研究中，研究人员通过全面的转录组学和网络分析，识别了与海马功能障碍相关的关键AD相关基因和通路，包括多巴胺能突触、钙信号、MAPK级联、神经递质受体活性、胶质细胞分化和神经炎症。PPI网络分析进一步突出了CAMK1G、PRKCB和MAP3K9作为具有潜在调控作用的核心枢纽蛋白。这些19个优先靶点通过分子对接和分子动力学模拟与天然黄酮类化合物（木犀草素、芹菜素和小檗碱）进行了系统评估，以探究其结合亲和力、稳定性和功能相关性。其中，木犀草素始终表现出最强的相互作用（ΔG_bind：-6.28至-9.24 kcal/mol）、降低的RMSD和RMSF、紧凑的SASA、持久的氢键以及能量有利的构象景观（PCA/FEL），表明催化环和调控环的稳健稳定化。这些结构和动态特征表明，木犀草素能够调节关键的AD相关激酶，减轻神经炎症信号，保护突触完整性，并可能抑制淀粉样蛋白和Tau病理。总体而言，这些发现建立了转录组来源的AD靶点与黄酮类化合物介导的多靶点调节之间的明确机制联系，为未来阿尔茨海默病的实验验证和治疗探索提供了强有力的依据。本研究代表了第一个整合通路富集、基因优先排序和AD中黄酮-激酶相互作用的深入计算机表征的框架。