基于中医药整合药理学和转录组学探究丹栀逍遥散防治焦虑肝郁化火证的作用机制

焦虑症（ anxiety ）是以广泛和持续性焦虑或反复发作的惊恐不安为主要特征，常伴有自主神经紊乱及运动性不安的疾病。临床常分为广泛性焦虑与惊恐障碍2 种主要形式。广泛性焦虑包括精神性焦虑、躯体性焦虑、心理性警觉、运动性不安和睡眠障碍；惊恐障碍患者主要在日常活动中，突然出现强烈恐惧的典型表现[1] 。中国国家卫生与健康调查统计数据显示，在2012 —2013 年我国城市广泛性焦虑障碍（36 086 例焦虑患者）的发病率为5.3%[2] 。在新型冠状病毒肺炎流行期间，医护人员的焦虑患病率高达40%[3] 。在后疫情时代，新型冠状病毒疫情对人们的生活、出行、学习以及自由活动产生了较大影响，从而产生负面情绪[4] 。因此，对于焦虑症的研究，特别是寻找中医药防治焦虑症的有效方药并阐明其机制是当前的关键和热点，甚至在未来几年将会是研究热点。

焦虑症实证层面的发病机制在于久郁化火，诚如朱丹溪在《丹溪心法 · 六郁》中所言：“气血冲和，万病不生，一有怫郁，诸病生焉”“凡郁病必先气病”；刘河间曰：“六气皆从火化……五志所伤皆化为热”。可见“气- 郁- 火”是焦虑症发生发展的病理过程，焦虑症实证患者临床症状表现为一派阳亢行为，在精神上表现为提心吊胆、主观过虑，但不能控制；行为上表现为患者躯体性主观和客观不适，如心悸、眩晕、多汗、呼吸困难等，以及运动性不安；除此之外，还有睡眠障碍和心理警觉等[5] 。焦虑症患者从临床表现来看更符合“阴不制阳”的机制，“郁”则是发病的一个不可或缺的枢纽环节。焦虑症虚证为主层面的发病机制在于心神耗损，久病必虚，《灵枢邪客》云：“心者，五脏六腑之大主也，精神之所舍也。”焦虑患者长期的急躁易怒暗损人体的津液，其中心阴首当其冲，《灵枢经· 口问》亦指出：“悲哀愁忧则心动，心动则五脏六腑皆摇”。心血不足、心阴亏耗则心失所养，阴不制阳，阴阳皆损，表现出心情郁闷、心慌多梦等。

焦虑证初期主要以肝郁气滞为主，渐则痰凝、血郁、化热，后期以气血不足、阴虚火旺为主，病位主要在心、肝、肾，尤以肝为甚，气机失调是发病的关键。有研究显示焦虑证候要素频率高于 5% 的有火热、气滞、阴虚[6] 。丹栀逍遥散由牡丹皮、栀子、茯苓、白术、白芍、当归、甘草和柴胡组成，出自明代薛己《内科摘要· 卷下十一· 各症方药》，具有清热泻火、建中补脾、养血柔肝和疏肝调达的功效，临床上丹栀逍遥散已经广泛运用在肝郁化火型广泛性焦虑症（障碍）的治疗中，研究结果显示可有效改善患者的焦虑症状[7-8] 。

1 材料与方法

1.1 丹栀逍遥散化学成分库的建立

利用中医药整合药理学平台 [9] （http://www. tcmip.cn/ ，TCMIP v2.0 ）的中药材数据库和中药成分数据库，检索“牡丹皮”“栀子”“柴胡”“白芍”“白术”“茯苓”“当归”“甘草”所包含的化学成分信息，并建立“丹栀逍遥散化学成分库”。

1.2 药物靶标预测

利用 TCMIP v2.0 ，设置Tanimoto Score 为0.8 （高准确度），提取“丹栀逍遥散化学成分库”中药物成分的化学分子指纹，预测化学成分的候选靶点谱，并将其与Drug Bank （http://www.drugbank.ca/ ）数据库中获美国食品药品监督管理局批准的抗焦虑药物指纹进行比对，将选定评分＞0.8 相似性药物靶点作为丹栀逍遥散药物组合的候选靶点[10] 。

1.3 疾病 、 证候靶标确定

在 TCMIP v2.0 的“疾病相关分子库”中，以“anxiety ”作为疾病关键词进行检索，检索出的所有靶点作为焦虑症的候选靶点集。

2020 年《国际中医临床实践指南· 焦虑症》将肝郁化火证临床表现定义为情绪不宁、郁闷烦躁、胸胁胀痛、脘闷嗳气、不思饮食、大便不调，或见急躁易怒、口苦而干，或头痛、目赤、耳鸣，或嘈杂吞酸、大便秘结，舌质红、苔黄，脉弦或弦数[11] 。《广泛性焦虑障碍中医证候辨证分型及量化分级标准专家共识》将肝郁化火证临床表现定义为主症烦躁易怒，善太息，咽中不适，如物梗阻，敏感多疑；次症为头痛，面红目赤，口苦咽干，两肋胀满，脘腹不适，痞满纳差，失眠多梦，女子月经不调，舌质红或淡红、苔黄或苔白，脉弦数[12] 。总结以上指南和专家共识，最后确定肝郁化火证临床表型为烦躁、易怒、恐慌、紧张、担忧、紧张性头痛、头晕、呼吸急促、心动过速、坐卧不宁、出汗、上腹不适、注意力困难 。在 TCMIP v2.0 平台数据库，以 irritation 、 anger 、 panic 、 nervousness 、 worry 、 tension headache 、 dizziness 、 shortness of breath 、 tachycardia 、 restlessness 、 sweating 、 upper abdominal discomfort 、 difficulty paying attention 为关键词，获取肝郁化火证证候相关靶标。

1.4 蛋白质相互作用 （ protein-protein interaction ，PPI ） 网络构建 、 分析与可视化

通过 TCMIP v2.0 的“网络分析模块”对PPI 网络进行拓扑分析，选取中药靶点- 疾病基因的关键节点（hubs ）作为丹栀逍遥散活性成分的候选靶点，然后对关键靶标网络进行通路富集分析，选取P 值（P ＜0.05 ） 排名前20 的通路[13] ，构建完整的“中药- 核心靶标- 关键通路”多层次关联网络，并进行可视化数据分析[13] 。

1.5 核心靶点的微阵列芯片验证

为检验中医药整合药理学得到的核心靶点的可靠性，以“ anxiety ”为关键词，检索美国国立生物技术信息中心的GEO 数据库（https://www.ncbi.nlm. nih.gov/geo/ ），筛选anxiety 相关的微阵列芯片数据集。采用R 软件（3.6.3 版本）进行统计分析与可视化：利用GEOquery 包（2.54.1 版本）[14] 下载数据，limma 包（3.42.2 版本）[15] 进行差异表达基因分析，umap 包（0.2.7.0 版本）进行UMAP （uniform manifold approximation and projection ）分析，ggplot2 包（3.3.3 版本）进行可视化，ComplexHeatmap 包（2.2.0 版本）[16] 绘制热图 。

1.6 实验动物

9 只雄性SFP 级SD 大鼠，体质量（200 ±20 ）g ，11 周龄。动物来源于北京斯贝福生物技术科技有限公司，生产许可证编号：SCXK （京）2019-0010 。在云南中医药大学动物房适应性饲养3 d ，自由饮水和摄取食物（普通饲料来源于云南楚商科技有限公司，生产批号20210711 ）。行为学实验前1 d 搬至行为学实验室，适应24 h ，保持与动物房相同的温度（25 ℃ ）和湿度（60% ～70% ）。 动物实验经云南中医药大学动物实验伦理委员会批准，批准文号 R-062021G106 。

1.7 药 物

丹栀逍遥散组成药物的饮片购买于云南中医药大学汇贤中医医院，由主管中药师闫永启主任鉴定，均符合《中国药典》 2020 年版规定。8 味中药，每付共计78 g ，见表1 。 将全部饮片浸泡在 10 倍体积的纯水中，浸泡 1 h ，先武火煮开，后改用文火，煮 1 h 后滤过药液。然后加入 10 倍体积的纯水继续二煎， 1 h 后混合 2 次滤过后的药液。最后再将药物浓缩成含生药 1.56 g/mL 的水煎液，将煮好的水煎液放置冰箱保鲜层，每天给药前 2 h 取出当天剂量，放置常温环境待用。

1.8 动物 分组 、 给药及模型制备

9 只大鼠适应性喂养 3 d ，按照随机数字表法分组，分为丹栀逍遥散组、模型组和对照组，每组 3 只。 3 组均给予正常普通饲料和自由饮水，每天 10:00 时，丹栀逍遥散组大鼠给予丹栀逍遥散水提物（生药 15.6 g/kg ），对照组和模型组大鼠给予等量的生理盐水（ 10 mL/kg ），共 2 周。除了对照组外，每天下午随机时间开始对丹栀逍遥散组和模型组大鼠连续 2 周给予随机慢性束缚应激（每天随机使用特制束缚架将模型组和丹栀逍遥散组大鼠固定在内不能活动，束缚架如圆筒状，大鼠头部和四肢从圆筒孔中暴露出来，每天 1 次，每次 2 h ）进行造模。每隔 3 d 测 1 次体质量。第 14 天 19:00 时开始进行高架十字迷宫实验，每只大鼠 5 min ，按时间顺序进行。

1.9 高架十字迷宫实验

采用高架十字迷宫实验测试动物焦虑程度。高架十字迷宫实验器材由 2 个开放臂（ 0.35 m × 0.60 m ）、 2 个封闭臂（ 0.35 m × 0.60 m ）以及交叉连接部（ 0.05 m × 0.05 m ）的开阔平台组成。封闭臂有 0.15 m 高的黑色挡板，十字迷宫平台离地面 0.5 m 。保持昼夜节律的光暗调整，实验室在实验过程中保持无灯光状态，开始将大鼠面向开臂放在高架十字迷宫的交叉连接部中央，同时开始软件记录大鼠在高架十字迷宫实验平台的运动轨迹以及在不同部位的时间和次数。 ① 进入开放臂次数（ open arm entry ， OE ）：即进入到任一开放臂的次数，以大鼠 4 只爪子均进入到臂内为准，中途 1 只爪子从该臂中完全退出则为该次进入活动完成。 ② 进入开放臂时间（ open arm time ， OT ）：即进入开放臂的时间（ s ）。 ③ 进入封闭臂次数（ close arm entry ， CE ）：即进入任一封闭臂的次数，以大鼠 4 只爪子均进入到臂内为准。 ④ 进入封闭臂时间（ close arm time ， CT ）：即进入封闭臂的时间（ s ）。由以上指标计算进入开放臂次数比例（ OE% ）、开放臂停留时间比例（ OT% ）。上述记录由 SuperMaze 动物行为学视频分析系统完成，以排除测试过程的人为干扰因素。每次测试结束后清除粪便，并用 75% 乙醇清理迷宫，以掩盖气味。待乙醇蒸发后方可开始下 1 只动物实验。

OE% ＝ OE/(OE ＋ CE)

OT% ＝ OT/(OT ＋ CT)

1.10 qRT-PCR 验证

1.10.1 主要仪器和材料 CFX Connect Real-Time System 荧光定量 PCR 仪（美国 Bio-rad ）， D3024R 型离心机（美国 SCILOGEX ）， Analytikjena-Easycycler 普通 PCR 仪（德国 Analytik Jena ）， Scandrop100 型超微量核酸蛋白测定仪（德国 Analytik Jena ），移液器（美国 Bio-rad ）， M2 全自动加样仪（四川科劲生物科技有限公司） 。 TUREscript 1st Stand cDNA SYNTHESIS Kit 反转录试剂盒（北京艾德莱生物科技有限公司，货号 pc70 ）； 2 × SYBR® Green 预混液（成都丹凤科技有限公司，货号 DF001 ）。

1.10.2 qRT-PCR 实验 各组大鼠麻醉后取脑组织，在冰块上分离海马组织和前额叶皮质，然后将海马组织、前额叶皮质和剩余脑组织分别装于 1.5 mL 离心管，迅速在液氮中冻存，放置于 −80 ℃冰箱保存，委托北京诺禾致源科技股份有限公司完成测序验证。

提取对照组、丹栀逍遥散组和模型组大鼠海马组织的总 RNA ，采用 SYBR GREEN I 检测对照组、模型组和丹栀逍遥散组大鼠海马组织 RNA 样品 3 个目的基因腺苷酸环化酶 5 （ adenylate cyclase 5 ， ADCY5 ）、谷氨酸离子型受体 NMDA 型亚基 2A （ GRIN2A-glutamate ionotropic receptor NMDA type subunit 2A Gene ， GRIN2A ） 、脑源性神经营养因子（ brain-derived neurotrophic factor ， BDNF ）的 mRNA 表达量，用 2−ΔΔCt 方法计算各基因 mRNA 相对表达量。以磷酸甘油醛脱氢酶（ reduced glyceraldehyde-phosphate dehydrogenase ， GAPDH ）为内参基因，所有实验共进行 3 次，引物合成信息见表 2 。

1.11 统计分析

采用 SPSS 26.0 软件进行分析，计量资料以 表示。两两比较若方差齐则采用 LSD 检验，若方差不齐，使用 Dunnett T3 检验方法。组间比较若方差齐采用单因素方差 One way ANOVA 分析；若方差不齐采用 Dunnett-t 检验。 P ＜ 0.05 表示差异有统计学差异。数据图采用生信工具 - 仙桃学术（ https://www.xiantao.love/products ）绘制。

2 结果

2.1 丹栀逍遥散成分靶标预测及分析

通过 TCMIPV2.0 平台对丹栀逍遥散化学成分的分析结果可知，共包含 430 种成分，候选靶标 520 个。对候选靶标进行基因功能分析和通路富集分析，发现其功能大多与谷氨酸受体信号通路、谷氨酸代谢过程、离子型谷氨酸受体信号通路、视黄醇代谢过程、 L - 谷氨酸转运、兴奋性突触后电位、 γ- 氨基丁 酸信号通路正向调节密切相关，排名前 20 的功能和通路如图 1 所示。

2.2 焦虑疾病靶标及相关基因

对全疾病相关潜在靶标合并去重后，共获得 209 个，其中与丹栀逍遥散密切相关的共同靶基因共 13 个，分别是长链酯酰辅酶 A 合成酶 4 （ Acyl-CoA synthetase long chain family member 4 ， ACSL4 ）、 ADCY5 、 醛脱氢酶 5 家族成员 A1 （ aldehyde dehydrogenase 5 family member A1 ， ALDH5A1 ）、 ATP 酶 Na+/K+ 转运亚基 α3 （ ATPase Na+/K+transporting subunit alpha 3 ， ATP1A3 ）、钙敏感受体（ calcium sensing receptor ， CaSR ）、儿茶酚 -O- 甲基转移酶（ catecholO-methyltransferase ， COMT ）、 D - 氨基酸氧化酶（ D -amino acid oxidase ， DAO ）、 GRIN2A 、 GNAS （ GNAS complex locus ）、 PI3K 催化亚基 p110alpha 的编码基因（ PIK3CA ）、转化生长因子受体 2 （ transforming growth factor beta receptor 2 ， TGFBR2 ）、雌激素受体 1 （ estrogen receptor 1 ， ESR1 ）、核受体亚家族 3 C 组成员 1 （ nuclear receptor subfamily 3 group C member 1 ， NR3C1 ） 。

2.3 潜在靶点 PPI 分析

通过对 13 个靶点进行 PPI 分析， PIK3CA 、 ACSL4 、 ALDH5A1 与其他靶标无关联性。其他 10 个靶标体现出较高的中心度，其中以 ESR1 为主，见图 2 。故在后续的验证中优先验证 ADCY5 、 ATP1A3 、 CASR 、 COMT 、 DAO 、 GRIN2A 、 GNAS 、 TGFBR2 、 ESR1 、 NR3C1 。

2.4 丹栀逍遥散防治焦虑的 “ 中药 - 靶点 - 通路 ” 多维网络分析

TCMIP v2.0 平台采用人工智能、数据挖掘、网络计算及可视化等方法和技术，形成 7 大整合药理学分析模块，包括疾病相关分子集及其功能挖掘、证候相关分子挖掘及功能分析、中药（含方剂）靶标预测及功能分析、中药药性相关分子挖掘及功能分析、 组方用药规律分析、中医药关联网络挖掘和反向查找中药（含方剂）。丹栀逍遥散防治焦虑的“中药 - 靶点 - 通路”多维网络图是基于“中医药关联网络挖掘”构建，根据中药方剂、中药材、中药成分、核心靶标、 通路、疾病和证候之间的关联性，可任意选择不同类型的节点，构建中药多维关联网络 [17] 。

如图3 所示，图中五角星表示方剂丹栀逍遥散，正方形表示中药，圆形表示靶标，椭圆形表示肝郁化火证临床表型，三角形表示丹栀逍遥散防治焦虑肝郁化火证的相关通路。通过网络图可以一目了然地发现丹栀逍遥散防治焦虑肝郁化火证的关键网络靶标及其作用机制，有助于揭示其科学内涵。

2.5 核心靶点微阵列芯片验证结果

最终选定微阵列芯片数据集 GSE4035 。该数据集基于 GPL1261 [Mouse430_2] Affymetrix Mouse Genome 430 2.0 Array 阵列平台，包含 12 例焦虑 RNA 样本（组 2 ）和 12 例健康 RNA 样本（组 1 ）。

由图 4 可见，基因表达谱区分度良好，保证了后续差异表达分析的准确性。去掉空值后的 ID 总数为 20 814 个，满足 |log2(FC)| ＞ 1 和 P adj ＜ 0.05 阈值的 ID 有 972 个，在这阈值下，组 2 中高表达（ log2FC 为正）的数目有 501 个，组 1 中高表达（ log2FC 为负）的数目有 471 个；组间差异基因火山图见图 5 ，热图见图 6 。

GSE4035 数据集的差异表达基因与ADCY5 、ATP1A3 、CASR 、COMT 、DAO 、GRIN2A 、GNAS 、TGFBR2 、ESR1 、NR3C1 得到3 个交集基因，分别是ESR1 、GRIN2A 和ADCY5 ，证明了上述分析的可靠性。

2.6 动物实验验证结果

2.6.1 一般情况和体质量 造模期间，丹栀逍遥散组与对照组的大鼠皮毛光泽顺滑，二便正常，日常活动灵敏。模型组大鼠变得较为暴躁，皮毛光泽暗淡。对照组大鼠体质量增长较快，丹栀逍遥散组和模型组大鼠的体质量增长较缓，与对照组相比差异均具有统计学意义（ P ＜ 0.05 ）。

2.6.2 高架十字迷宫实验 如图7 所示，在OT% 方面，模型组低于对照组，差异具有统计学意义（P ＜0.01 ）；丹栀逍遥散组高于模型组，差异不具有统计学意义（P ＝0.069 ）。在OE% 方面，模型组高于对 照组，差异具有统计学意义（ P ＜ 0.01 ）；丹栀逍遥散组低于模型组，差异具有统计学意义（ P ＜ 0.01 ）。

2.6.3 qRT-PCR 实验结果 样品中 ESR1 mRNA 表达丰度低，有些样本没有 C t 值，无法进行后续计算。与对照组比较，模型组 ADCY5 mRNA 表达量升高，但差异无统计学意义（P ＝0.056 7 ），GRIN2A mRNA 表达量显著降低（ P ＜ 0.05 ）；与模型组比较，丹栀逍遥散组 ADCY5 mRNA 表达量显著降低（ P ＜ 0.05 ）， GRIN2A mRNA 表达量显著升高（ P ＜ 0.001 ）（图 8 ）。表明这 2 个核心靶点基因的相对定量值与整合网络药理学的数据趋势一致，证明整合网络药理学和 GEO 数据挖掘结果的可靠且重复性好。

3 讨论

3.1 焦虑肝郁化火证中医病机与治疗策略

焦虑属于中医“郁”的范畴。南宋的陈言提出“脏气不行，郁而生涎……皆七气所生所成”，郁是疾病发生的罪魁祸首，郁既有病的概念，也是证的范畴；明清时期，“郁”多被单列一门，《医学正传》首用“郁证”作为病名进行论述[5]。但对于“郁”的治疗不仅是疏散之法，有研究显示古代医家治疗郁多以“培土”为基础，辅以疏肝，酌用清热、行气、安神之品[6]。

通过焦虑肝郁化火证的临床证候表型可以看出，焦虑肝郁化火证患者实际临床中往往呈现虚实夹杂，本虚标实，阴虚阳亢的病机，譬如运动性紧张、坐卧不宁、紧张性头痛、颤抖、心动过速或呼吸急促、无法放松属于虚阳外亢的临床表现；担忧、注意困难、头重脚轻、出汗、头晕、口干属于正虚的临床表现。丹栀逍遥散组方和功效划分恰好符合焦虑肝郁化火证的整体病机，其中丹皮、栀子清热泻火治疗阳亢于外，茯苓、炒白术、炙甘草建中补脾治疗脾气亏虚，当归、白芍、柴胡养血柔肝补肝体而疏肝用。总体而言，焦虑肝郁化火证涉及主要病位在肝、脾，主要病性虚、实，丹栀逍遥散的组方符合焦虑肝郁化火证的发病机制。本研究借助TCMIP v2.0平台，从中药成分、证候要素、疾病靶点的多层次研究以及从药物自身出发，较为全面地展现了中药-证候-疾病-靶标-通路的衔接性研究，形成中医药现代研发体系，从分子角度解释了中医药疗效内涵[18]。

3.2 整合网络药理学结果分析

丹栀逍遥散化学成分分析，谷氨酸受体包括离子型受体和代谢性受体，离子型谷氨酸受体包括 N - 甲基-D- 天冬氨酸（N-methyl-D-aspartic acid ，NMDA ）受体、α- 氨基-3- 羟基-5- 甲基-4- 异噁唑丙酸（α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid ，AMPA ）受体和Kainate 受体；应激状态时，NMDARs 会被过度激活，Ca2+ 通道反复开放，过多的Ca2+ 内流使细胞内Ca2+ 超载，引起神经细胞死亡[19] 。 应激状态下，上调的AMPAR 增加了脑内谷氨酸兴奋性毒性，出现细胞内钙累积，导致细胞稳态的破坏和脑组织的结构损伤[20] 。代谢型谷氨酸2 型受体和γ- 氨基丁酸B 受体属于C 家族G 蛋白偶联受体，分别是兴奋性神经递质谷氨酸和抑制性神经递质γ- 氨基丁酸的代谢型受体，γ- 氨基丁酸B 受体在中枢神经系统中介导缓慢的抑制性神经传导和突触可塑性调控，其活性的失调会导致细胞功能紊乱、损伤甚至细胞凋亡[21] 。

“中药 - 靶点- 通路”多维网络分析显示，丹栀逍遥散通过调控AMPA 受体、蛋白激酶A （protein kinase A ，PKA ）、AMPA 受体的转运、胰高血糖素信号中的PKA 激活、NMDA 受体激活、磷脂酰肌醇3- 激酶（phosphatidylinositide 3-kinases ，PI3K ）/ 蛋白激酶B （protein kinase B ，Akt ）等信号通路，对蛋白质结合、胞质溶胶、质膜、细胞质、线粒体等通路进行干预。

ADCY5 是一种腺苷酸环化酶，通过合成第2 信使环磷酸腺苷（cyclic adenosine monophosphate ，cAMP ）来介导G 蛋白偶联受体的信号传导。此外，cAMP 也是一种典型的第2 信使，其通过浓度的变化响应细胞外信号与细胞表面受体的结合，调节细胞内酶和非酶蛋白的活性，激活依赖于cAMP 的PKA 表达增加磷酸化水平，并调节BDNF 等基因的表达，从而发挥生物学效应[22] 。

GRIN2A 是重要的谷氨酸离子向受体NMDA 型亚基，焦虑样行为与海马中GRIN2A 的上调有关，可通过调节与Ca2+/ 钙调蛋白依赖性蛋白激酶II （calcium-calmodulin-dependent protein kinase II ，CaMKII ）的结合来影响突触的可塑性[23] 。DAO 是一种降解D- 丝氨酸的酶，是NMDAR 的主要内源性共激动剂，灭活或缺失DAO 将会加剧焦虑[24] 。

雌激素受体 1 （estrogen receptor 1 ，ESR1 ）二核苷酸重复等位基因与90 项症状清单（symptom check list 90 ，SCL-90 ）9 个亚分数和总分之间的关联结果显示，其与SCL-90 评分中焦虑、恐惧焦虑和总症状有显著关联[25] 。已有的研究表明，ESR1 变异与认知结果之间具有密切关系，初步证据表明ESR1 基因在某些情绪中的作用，发现ESR1 多态性与焦虑或抑郁症状关联显著[26] 。

本研究旨在探讨丹栀逍遥散治疗焦虑的本质，完善和补充现有的研究。目前有诸多学者对丹栀逍遥散治疗焦虑展开了深入研究，在生物学机制方面丹栀逍遥散可以显著降低焦虑大鼠炎症反应、机制细胞凋亡和和减轻护神经元损伤以及提高线粒体功能、提高脑源性神经营养因子和蛋白激酶蛋白表达，促进脑源性神经营养因子分泌，从而改善焦虑样行为。初步阐明了丹栀逍遥散治疗焦虑症的作用机制和主要成分协同作用的机制，同时支撑了丹栀逍遥散在防治疗焦虑症方面的应用。

3.3 GEO 数据库验证结果分析

GEO 数据库中转录组学验证结果显示，ESR1 、ADCY5 、GRIN2A 与焦虑肝郁化火证密切相关。雌激素在很大程度上通过其主要的ESR1 和ESR2 的细胞内激活来发挥其生物学作用。因此，雌激素受体的遗传变异可能会改变雌激素信号，从而影响女性患者抑郁症的易感性[27] 。也有研究显示，5- 羟色胺基因与ESR1 基因多态交互作用可能增加女性围绝经期抑郁障碍患病的风险[28] 。进一步验证了ESR1 基因是核心靶点基因之一。

ADCY5 可以通过调节BDNF 基因的表达进而对谷氨酸损伤的神经元具有保护作用，BDNF 可增加谷氨酸损伤皮层神经元NMDAR1 mRNA 及蛋白的表达，增强其NMDA 受体活性。同时GRIN2A 是谷氨酸离子向受体NMDA 型重要亚基，亦影响着NMDA 受体的活性。因此，推测丹栀逍遥散可能是通过下调ADCY5 和上调GRIN2A 基因的表达水平进一步上调BDNF 表达和增加NMDA 活性，发挥其生物学效应。

3.4 动物 实验验证 结果分析

高架十字迷宫实验结果显示，在 OT% 方面，虽然丹栀逍遥散组高于模型组，差异不具有统计学意义（P ＝0.069 ），但通过统计结果可以看出具有明显的趋势，其原因可能与样本量少有关。在OE% 方面，丹栀逍遥散组与模型组差异均具有统计学意义，说明丹栀逍遥散能改善焦虑模型大鼠进入开放臂的次数。

核心靶点基因验证结果显示， ADCY5 和GRIN2A mRNA 表达量模型组与对照组、丹栀逍遥散组与模型组差异均具有统计学意义；核心靶点经qPCR 实验验证，丹栀逍遥散能降低ADCY5 mRNA 表达（P ＜0.05 ），升高GRIN2A mRNA 的表达（P ＜0.001 ）。

4 结论

本研究对丹栀逍遥散组方候选靶标与焦虑症进行关联并以网络可视化形式展现，并预测其可能靶标和分子信号通路。丹栀逍遥散调控的主要通路有 AMPA 受体的激活、PKA 激活、胰高血糖素信号中的PKA 激活、NMDA 受体激活、PI3K/Akt 等。中医药整合药理学和GEO 数据集获取交集基因，包括ESR1 、ADCY5 、GRIN2A 基因；丹栀逍遥散可能通过上调ADCY5 mRNA 表达和下调GRIN2A mRNA 表达进一步调节BDNF 表达和NMDA 活性，进而调节BDNF 与其受体的结合和NMDA 活性，从而调节BDNF 与其受体结合，进一步改善神经元生长、突触传导与突触的可塑性，从而发挥其生物学效应。

利益冲突 所有作者均声明不存在利益冲突

参考文献（略）

来 源：潘赐明，任 雪，韦姗姗，曹蓓苓，王 维，杨 艳，陈文慧，徐宏喜．基于中医药整合药理学和转录组学探究丹栀逍遥散防治焦虑肝郁化火证的作用机制 [J]. 中草药, 2023, 54(12):3911-3920.