Guduchyadi Yoga对轻度认知障碍的多靶点计算验证：UHPLC-HRMS、网络药理学、分子对接、DFT、分子动力学与后动力学整合分析

《Results in Chemistry》：Chemical characterization and computational validation of Guduchyadi yoga in mild cognitive impairment: UHPLC-HRMS, network pharmacology, docking, DFT, dynamics, and post dynamics analysis.

　　

采用UHPLC-HRMS从Guduchyadi Yoga中鉴定并表征了368种植物化学物，其中55种生物活性化合物通过血脑屏障（BBB）和类药性筛选标准。分子对接（molecular docking）显示，异阿魏酸（isoferulic acid）、arctigenin和3,6,2-trimethoxyflavone与关键蛋白靶点AKT1、TNF和PPARG具有强而稳定的结合，其结合亲和力优于donepezil和lecanamab。密度泛函理论（DFT）计算揭示了先导化合物的电子性质和稳定性特征，表明其化学反应性与热力学稳定性的组合有利于治疗作用。分子动力学模拟（MD，500 ns）确认了在生理条件下通过氢键、疏水作用力和π-堆积（π-stacking）实现的动态稳定性和持久蛋白-配体相互作用。后动力学分析，包括主成分分析（PCA）和自由能景观（FEL），识别了复合物的稳定构象状态和低能盆地，尤其是PPARG-Arctigenin，表现出高稳定性。诸如CASP3（caspase-3）、TNF、IL6、AKT1、BCL2和TP53等中枢基因在早期神经变性中平衡凋亡和存活通路方面发挥关键作用，并被Guduchyadi Yoga植物成分有效调节。该研究突显了轻度认知障碍（MCI）中涉及凋亡调节、神经炎症、线粒体功能障碍和突触完整性的分子互作，定义了多靶点方法。计算见解支持Guduchyadi Yoga是一种科学上合理的多草药复方（polyherbal formulation），具有神经保护潜力，能够延缓或预防从MCI向痴呆的进展。该研究为进一步实验验证和神经退行性疾病转化临床研究提供了稳健的机制证据。

研究背景：轻度认知障碍（mild cognitive impairment, MCI）是正常认知老化与痴呆之间的关键过渡状态，其认知下降超过正常年龄相关变化但不满足痴呆标准。MCI影响65岁以上成人的15–20%，是阿尔茨海默病（Alzheimer's disease, AD）和其他痴呆的重要风险因素。其发病涉及内体失调、氧化应激、线粒体功能障碍、胆碱能缺陷和神经营养信号通路破坏等多层面改变，提示其可能为跨突触神经元-神经元断开综合征。当前胆碱酯酶抑制剂（donepezil、rivastigmine）或NMDA受体拮抗剂（memantine）仅提供症状缓解，未能阻止疾病进展，且尚无FDA批准的MCI药物。由于MCI病理复杂，多靶点策略被认为优于单靶点治疗。Guduchyadi Yoga是传统阿育吠陀多草药复方，作为Medhya Rasayana用于认知改善，并显示改善老年个体记忆和焦虑的临床效果；其主要成分Tinospora cordifolia（Guduchi）含生物碱、萜类、黄酮、酚类和皂苷等，具抗氧化、抗炎、胆碱酯酶抑制和神经递质调控等神经保护作用。然而，其与MCI相关分子靶点的相互作用机制仍不充分。尽管计算研究多聚焦AD，少关注MCI，因此研究人员开展整合计算研究以识别可调节MCI分子机制的多靶点植物成分。

关键技术方法：样本为GMP认证供应商IMPCL（Uttarakhand, India；批号19-ACH-LDA-539）提供的Guduchyadi Yoga，在AIIA, New Delhi提取。采用乙醇:水（70:30, v/v）提取后进行UHPLC-HRMS分析。以SwissADME和Molsoft筛选血脑屏障（BBB）通透性与Lipinski五规则；用SwissTargetPrediction预测化合物靶点，用NCBI Gene和GeneCards获取疾病靶点；用Cytoscape/STRING/CytoHubba构建蛋白-蛋白相互作用（PPI）网络，用DAVID进行基因本体（GO）和KEGG富集。用Schrodinger Glide进行分子对接，Gaussian 16以B3LYP/6-31G(d,p)进行DFT计算，Schrodinger Maestro进行500 ns MD模拟（TIP3P），并用MDAnalysis进行PCA和分子力学/广义玻恩表面积（MM-GBSA）分析。

研究结果：3.1 植物化学分析与计算机药理学评价：研究人员制备Guduchyadi Yoga水醇提取物，经UHPLC-HRMS正、负离子模式鉴定，分别暂定109和259种植物化学物，共368种，并用于后续网络药理学、分子对接和MD分析。3.2 生物活性化合物和疾病靶点筛选：368种化合物经SwissADME和Molsoft筛选，仅55种满足BBB通透性和Lipinski五规则；SwissTargetPrediction去重后得到1204个化合物靶点，GeneCards和NCBI Gene以“Mild-Cognitive Impairment”检索并去重后得到16035个疾病靶点，Venn交集得到942个中枢基因（hub genes）。3.3 基因表达和通路富集分析：通过STRING构建PPI并经Cytoscape CytoHubba分析，得到前15位核心靶点AKT1、TP53、TNF、SRC、IL6、EGFR、CTNNB1、STAT3、ESR1、HSP90AA1、CASP3、JUN、BCL2、HIF1A和PPARG。GO富集显示其与蛋白结合、酶结合、蛋白激酶结合、DNA结合转录激活活性等相关，定位涉及核、胞质、细胞质和核质；生物过程涉及凋亡、炎症反应、氧化应激信号和突触传递。KEGG富集涉及AD、PI3K-AKT、MAPK、TNF、凋亡、HIF-1、雌激素和Wnt/β-catenin信号通路。功能解释提示AKT1、STAT3、BCL2、EGFR和ESR1相关通路正调节可能促进神经元存活、减少凋亡、维持突触完整性并减轻淀粉样β和tau病理；TP53、TNF、IL6过度激活及SRC、CTNNB1失调则加剧神经炎症和凋亡。3.4 分子对接：AKT1–isoferulic acid复合物对接分数为?8.010 kcal/mol，与GLN61、TYR151形成氢键并与TYR59发生π–π堆积，疏水作用涉及TYR119、SER60和LEU57。TNF–3,6,2′-trimethoxyflavone为?7.891 kcal/mol，与GLN61、TYR151形成氢键并与TYR59发生π–π堆积。PPARG–arctigenin为?7.647 kcal/mol，与HIE449、TYR473和SER289形成氢键。donepezil和lecanamab在相同条件下分别为?5.635和?6.121，植物化合物的对接分数更有利，但对接分数不能单独预测生物效能或临床效果。3.5 DFT：isoferulic acid电离势为5.68 eV、能隙4.13 eV、亲电性ω=3.16 eV，提示较高化学反应性和亲电性；arctigenin电离势8.45 eV、能隙12.01 eV、亲电性ω=0.50 eV，提示高化学硬度和稳定性；3,6,2′-trimethoxyflavone能隙4.27 eV、亲电性ω=3.27 eV，提示强电子接受能力。结果表明3,6,2′-trimethoxyflavone和isoferulic acid反应性较高，arctigenin最稳定且化学惰性较强。3.6 分子动力学：500 ns MD显示TNF–3,6,2′-trimethoxyflavone中蛋白骨架RMSD约1.2–3.2 ?，配体稳定停留于预测活性位点，GLN61、TYR151和TYR59运动受限，氢键、疏水相互作用和水桥持续存在。PPARG–arctigenin的RMSD约1.0–3.0 ?，HIE449、TYR473和SER289持续参与氢键和疏水作用，结合口袋刚性增强，配体扭转应变小。AKT1–isoferulic acid蛋白骨架RMSD约4.0 ?，提示构象柔性但未失稳，配体稳定，GLN61和TYR151形成氢键，TYR59提供疏水和π–π锚定。所有复合物的回旋半径、溶剂可及表面积和配体内部氢键稳定，支持对接预测。3.7 后动力学分析：PCA显示AKT1复合物有两个主要簇，构象景观受限；PPARG轨迹紧凑并有两个主要构象状态，稳定性高；TNF轨迹较拉长但端点仍聚类。FEL显示AKT1有多个浅能峰，构象采样活跃但中等稳定；PPARG具有更宽更深的低能盆地和较少高能波动，最稳定；TNF为窄浅盆地，灵活性有限但中等稳定。MM-GBSA显示AKT1–isoferulic acid的ΔG_bind为?74.87±4.94 kcal/mol，PPARG–arctigenin为?74.81±5.24 kcal/mol，TNF–3,6,2′-trimethoxyflavone为?74.87±5.01 kcal/mol，表明三者具有强结合能和构象稳定性，支持进一步实验研究。

讨论部分总结：MCI是正常衰老与痴呆尤其AD之间的过渡阶段。通路富集和中枢基因分析显示CASP3、TNF、IL6、AKT1、BCL2、PPARG、β-catenin和TP53存在于AD、亨廷顿病、神经退行和多疾病通路及肌萎缩侧索硬化中，提示与痴呆进展相关的分子通路可能在MCI阶段已参与。CASP3为半胱氨酸依赖、天冬氨酸特异性蛋白酶，参与凋亡及神经可塑性和长时程抑制；AD中Aβ激活CASP3，裂解APP和tau，促进突触毒性，并与NLRP3炎症小体信号相关。TNF-α由小胶质细胞产生，经TNFR1和TNFR2信号，增加APP和BACE1表达，激活NF-κB，CSF中TNF-α升高与MCI向痴呆更快进展相关。IL-6经IL-6R/gp130激活JAK/STAT3和MAPK，AD脑和血浆中升高与认知损害及海马体积缩小负相关，并通过cdk5/p35导致tau磷酸化。AKT在AD中PI3K/Akt通路下调，导致tau过度磷酸化、GSK-3β过度活化和葡萄糖转运体表达下降；AKT可促进CREB依赖的神经保护基因转录并抑制促凋亡蛋白。BCL2家族抗凋亡蛋白控制钙稳态和线粒体外膜通透化，防止神经元死亡。β-catenin是Wnt信号主要效应分子，AD中受DKK1和GSK-3β抑制，导致突触丢失，并影响帕金森病、亨廷顿病和ALS。TP53和PPARG等也参与神经退行机制。研究强调MCI与多种神经退行疾病共享分子特征，涉及凋亡介质、促炎细胞因子和代偿性生存通路的互作。该工作流程与联合网络药理学、分子对接和MD优先植物成分-靶点对的研究一致，计算结果提示Guduchyadi Yoga化合物值得作为神经退行疾病相关通路的多靶点调节剂进一步研究。

结论翻译：本计算研究突出MCI相关通路的预测关系，包括凋亡调节因子（CASP3、TP53、BAX）、炎性细胞因子（TNF-α、IL-6）和代偿性促生存通路（PI3K/Akt、BCL2、PPARG、Wnt/β-catenin）。这些中枢基因在多个神经退行通路中的识别提示，与痴呆进展相关的病理级联可能在MCI阶段已参与。重要的是，Guduchyadi Yoga植物成分在计算条件下对这些靶点显示有利预测相互作用，提示其作为凋亡、炎症和神经元存活的多靶点调节剂值得进一步研究。这些计算见解为传统多草药复方作为MCI相关多靶点调节剂的进一步研究提供了机制依据。进一步实验和临床验证对于将这些发现转化为MCI和早期AD的有效治疗策略至关重要。未进行基于重对接的回顾性验证和配体姿态RMSD计算，因此预测对接姿态需要独立基于结构或实验验证。更广泛地，这些计算发现需要生化测定、动物研究和临床研究验证，然后才能确立治疗应用。