《Letters in Drug Design & Discovery》：Network pharmacology-based identification and molecular docking with dynamics simulation of phytochemicals from
Centella asiatica,
Panax ginseng, and
Rosmarinus officinalis as potential multi-target therapeutics against Ischemic Stroke
编辑推荐：
背景：缺血性卒中（ischemic stroke）是突发脑供血丧失导致神经元死亡和功能损害的全球主要神经病学挑战。尽管再灌注治疗已有进展，缺血性卒中的复杂病理生理、神经炎症、氧化应激和血管功能障碍使靶向药物治疗效果有限。传统草药已显示神经保护、抗炎和抗氧化作用
背景：缺血性卒中（ischemic stroke）是突发脑供血丧失导致神经元死亡和功能损害的全球主要神经病学挑战。尽管再灌注治疗已有进展，缺血性卒中的复杂病理生理、神经炎症、氧化应激和血管功能障碍使靶向药物治疗效果有限。传统草药已显示神经保护、抗炎和抗氧化作用，但其在缺血性卒中中的分子机制尚未充分明确。方法：本研究采用计算方法，包括网络药理学（network pharmacology）、蛋白-蛋白相互作用（protein-protein interaction, PPI）、基因本体（Gene Ontology, GO）富集、KEGG（Kyoto Encyclopedia of Genes and Genomes）通路分析、化合物-靶点网络构建、药代动力学表征、分子对接（molecular docking）和分子动力学模拟（molecular dynamics, MD）（n = 3），以鉴定这些草药中的重要生物活性化合物和分子靶点。结果：研究发现498个草药相关靶点和99个缺血性卒中相关靶点，重叠靶点9个。PPI网络分析形成含28个节点的网络，提示较高相互作用潜力。KEGG关键通路为肾素分泌、血管平滑肌收缩和神经配体-受体相互作用。药代动力学分析显示Carnosic acid、Asiatic acid和Ginsenosides具有良好血脑屏障（blood-brain barrier, BBB）通透性。结论：在所研究的植物化学物中，Carnosic acid对关键靶点显示最佳结合亲和力。进一步，MD证实这些相互作用的稳定性；然而MM/PBSA（molecular mechanics Poisson–Boltzmann surface area）预测Carnosic acid与DAPK1和ESR1结合有利，与PRKCH结合边际，与PDE4D结合较弱。最终，Carnosic acid似乎是具有最佳神经保护潜力的缺血性卒中候选物，尚需进一步实验验证。
缺血性卒中是全球重要致死和长期致残原因，主要由血栓或栓塞导致脑血流阻断，约占全部卒中的85%。再灌注治疗如重组组织型纤溶酶原激活剂（recombinant tissue plasminogen activator, rt-PA）静脉溶栓和血管内机械取栓受时间窗窄和禁忌证多限制，多数患者仍缺乏有效神经保护治疗。缺血性卒中病理为多机制级联，包括兴奋性毒性、离子失衡、线粒体功能障碍、氧化/亚硝化应激、神经炎症、血脑屏障（blood-brain barrier, BBB）破坏和程序性细胞死亡；单靶点策略不足，需多靶点治疗。积雪草、人参和迷迭香具神经保护、抗炎、抗氧化和血管保护证据，但三者共享靶点、关键富集通路及先导植物化学物缺乏系统计算分析，故研究人员开展本研究。研究整合网络药理学、分子对接和分子动力学模拟，鉴定9个共享靶点及关键通路，并确定Carnosic acid为最佳多靶点先导物；意义在于为缺血性卒中多靶点治疗提供计算依据和候选分子，仍需实验验证。
关键技术方法概括（不超过250字）：研究人员检索DisGeNET、OMIM、GeneCards和NCBI Gene获得99个缺血性卒中靶点，检索IMPPAT获得498个草药靶点，以InteractiveVenn取交集9个；用STRING构建PPI网络，用GeneMANIA扩展功能网络，并行GO、Reactome和KEGG富集；用Cytoscape v3.9.1构建化合物-靶点网络（compound-target network, CTN）；用Molinspiration等进行ADMET（吸收、分布、代谢、排泄和毒性）与BBB预测；对ESR1（PDB 1GWQ）、PRKCH（8FP3）、DAPK1（5AV2）和PDE4D（6NJJ）行Glide SP/XP对接及再对接验证；用GROMACS行300 ns×3重复MD（n = 3），分析RMSD（均方根偏差）、RMSF（均方根涨落）、Rg（回旋半径）、SASA（溶剂可及表面积）和氢键；用gmx_MMPBSA计算MM/PBSA。未涉及实验样本队列来源。
结果：3.1 共同分子靶点鉴定：通过多数据库Venn交集研究得出9个共同靶点，包括HDAC1、PRKCH、EDNRA、ALOX5AP、HCAR2、DAPK1、ACE、ESR1和PDE4D，提示表观遗传、血管张力、炎症、凋亡、肾素-血管紧张素系统、激素信号和第二信使调控等靶点覆盖。3.2 PPI网络分析：STRING置信度≥0.4构建28个节点、94条边、平均度6.71、富集p = 0.001的网络；ACE、ESR1和DAPK1具高中心性，提示关键治疗节点。3.3 GeneMANIA相互作用功能网络：共表达43.19%、物理相互作用29.98%、共定位23.74%、共享蛋白域1.84%、遗传相互作用1.26%，提示共调控、物理接触和亚细胞共定位。3.4 GO和通路分析：BP富集血管收缩、激素调节系统动脉压和血管相关平滑肌收缩；MF富集转录辅阻遏物活性和组蛋白脱乙酰酶结合；CC富集组蛋白脱乙酰酶复合物；Reactome富集肽配体结合受体信号；KEGG富集肾素分泌、血管平滑肌收缩和神经活性配体-受体相互作用。3.5 化合物-靶点网络分析：Cytoscape构建7个化合物与9个靶点网络，Carnosic acid和Ginsenosides靶覆盖最广，支持多药理作用。3.6 药理学性质：ADMET显示Carnosic acid分子量332.44 Da、miLogP 4.60、TPSA 77.75、无Lipinski违规、HIA 0.96；除Ginsenosides 0.91外其余BBB评分0.50–0.74；Ginseng tetrapeptide最不类药（miLogP -5.25、TPSA 252.32、HIA 0.66、2项违规）。3.7 参考抑制剂再对接验证：ZTW、Y3I、KMP和KR7再对接得分分别为-8.01、-9.652、-8.262和-7.323 kcal/mol，RMSD在2 ?内并复现关键残基接触，验证对接流程稳健。3.8 分子对接结果：7个植物化学物对4个靶点行Glide XP对接，Carnosic acid对ESR1、PRKCH、DAPK1和PDE4D得分分别为-8.8、-10.3、-9.5和-9.1 kcal/mol，均为最佳，结合MM/GBSA为-49.05至-57.16 kcal/mol，被鉴定为先导物。3.9 Carnosic acid相互作用分析：3.9.1 Carnosic acid-PRKCH：最高对接得分-10.3 kcal/mol，Asn484氢键1.58/2.04 ?，疏水残基接触，未形成铰链氢键；3.9.2 Carnosic acid-DAPK1：得分-9.5 kcal/mol，Lys42氢键和盐桥2.27/4.75 ?，疏水接触；3.9.3 Carnosic acid-ESR1：得分-8.8 kcal/mol，Leu346氢键1.75/1.81 ?，Phe404 π-π堆积5.06 ?，疏水接触；3.9.4 Carnosic acid-PDE4D：得分-9.1 kcal/mol，Asn375氢键1.73/2.03 ?，Phe599 π-π堆积5.17 ?，疏水接触。3.10 分子动力学模拟分析：4系统×3重复×300 ns共3.6 μs；各复合物RMSD稳定约0.10–0.35 nm，无解离或展开；结合口袋RMSF低，Rg稳定，氢键持续；5AV2最稳定（RMSD 0.156±0.002 nm，至少一个氢键占98.4±2.4%帧），6NJJ氢键较少（79.0±10.6%帧）。3.11 结合自由能分析—MMPBSA：DAPK1 ΔTOTAL -27.28±4.19 kcal/mol、ESR1 -20.99±3.56 kcal/mol、PRKCH -4.85±4.06 kcal/mol、PDE4D +12.31±8.00 kcal/mol；范德华和非极性溶剂化为主要驱动，MD平均MM/PBSA较单构象MM/GBSA更可靠，最可能靶点为DAPK1和ESR1。
讨论总结：作者认为ACE为共享靶点并与肾素分泌通路相关，脑局部肾素-血管紧张素系统过度激活参与神经炎症、BBB破坏、凋亡和氧化应激；PRKCH静态对接最佳但MD平均MM/PBSA仅边际且未复现Y3I铰链氢键，故不能确认其为最高亲和靶点；DAPK1在对接、MM/GBSA、MM/PBSA、RMSD和氢键持续性上最一致；Carnosic acid药代特征利于中枢神经系统，但需关注CYP450（细胞色素P450酶）抑制风险；GeneMANIA提示共表达和共定位的网络交叉；9个靶点涉及兴奋毒性/凋亡、血管张力/血栓炎症及转录/激素调节。局限性包括结果均为计算预测，需体外和体内验证；数据库注释偏倚、MD采样有限、MM/PBSA忽略构象熵、ADMET和BBB为模型值等。研究结论翻译：当前网络药理学和分子模拟分析鉴定出缺血性卒中病理与积雪草、人参和迷迭香药理空间重叠的9个特异分子靶点，并为这些传统草药的多靶点治疗机制提供系统和计算证据。富集分析表明肾素分泌通路、血管平滑肌收缩和神经配体-受体相互作用通路是这些靶点介导缺血性卒中病理的主要机制轴。对7种植物化学物的药代筛选确认Carnosic acid具有最有利于中枢神经系统药物开发的类药性质。对四个关键蛋白靶点（ESR1、PRKCH、DAPK1和PDE4D）的分子对接在所有四个靶点中均鉴定Carnosic acid为先导化合物，结合得分分别为-8.8、-10.3、-9.5和-9.1 kcal/mol，在每个靶点类别中记录到最佳亲和力。300 ns分子动力学模拟支持所有Carnosic acid-蛋白复合物的结构稳定性和动态相互作用完整性，在分子水平显示紧密稳定结合模式。然而MM/PBSA表明对DAPK1和ESR1结合有利，对PRKCH边际，对PDE4D较不利。总体，这些结果将迷迭香来源Carnosic acid定位为有前景的多靶点神经保护植物化学物，能够同时调节PRKCH介导的凋亡、DAPK1介导的兴奋性毒性、ESR1介导的神经保护信号和PDE4D介导的神经炎症，代表一个经计算验证、值得在体外和体内实验中研究作为缺血性卒中新型治疗分子的先导分子。
