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

潘赐明，任 雪，韦姗姗，曹蓓苓，王 维，杨 艳，陈文慧，徐宏喜

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

潘赐明1, 2，任 雪2, 3，韦姗姗1*，曹蓓苓1，王 维1，杨 艳1，陈文慧1，徐宏喜2

1. 云南中医药大学，云南 昆明 650500 2. 上海中医药大学，上海 201203 3. 黑龙江中医药大学，黑龙江 哈尔滨 150040

基于中医药整合药理学和转录组学探究丹栀逍遥散防治焦虑肝郁化火证的作用机制。采用中医药整合药理学研究平台v2.0，检索并获取丹栀逍遥散的活性成分及靶标、焦虑肝郁化火证靶标信息，构建中药-成分-靶点-通路网络。将交集基因导入String数据库，构建交集靶点的蛋白互作网络，以确定最终需要验证的靶标。利用GEO数据库转录组学验证整合网络药理学结果，同时采用动物实验进行验证。共筛选出与焦虑肝郁化火证相关的靶标13个，丹栀逍遥散活性成分430种。整合药理学结果显示，丹栀逍遥散防治焦虑肝郁化火证主要通过调控α-氨基-3-羟基-5-甲基-4-异噁唑丙酸受体（α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid receptor，AMPAR）、蛋白激酶A（protein kinase A，PKA）、-甲基--天冬氨酸（-methyl--aspartic acid，NMDA）受体的激活及磷脂酰肌醇3-激酶（phosphatidylinositide 3-kinases，PI3K）/蛋白激酶B（protein kinase B，Akt）等信号通路。GEO微阵列芯片GSE4035数据集转录组学关键靶标雌激素受体α（estrogen receptor 1，）、腺苷酸环化酶 5（adenylate cyclase 5，）、谷氨酸离子型受体NMDA型亚基2A（glutamate ionotropic receptor NMDA type subunit 2A，）基因为差异表达基因（differential gene expression，DEGs）。实时荧光定量聚合酶链反应（quantitative real-time polymerase chain reaction，qRT-PCR）结果显示，mRNA表达上调和mRNA表达下调。丹栀逍遥散可能是通过上调mRNA表达和下调mRNA表达进而调节脑源性神经营养因子（brain-derived neurotrophic factor，BDNF）与其受体的结合和谷氨酸神经传递介导，进一步改善神经元生长、突触传导与突触的可塑性，从而发挥其生物学效应。

丹栀逍遥散；焦虑；肝郁化火证；整合药理学；转录组学

焦虑症（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）作为丹栀逍遥散活性成分的候选靶点，然后对关键靶标网络进行通路富集分析，选取值（＜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 丹栀逍遥散药物信息

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，）、谷氨酸离子型受体NMDA型亚基2A（GRIN2A-glutamate ionotropic receptor NMDA type subunit 2A Gene，）、脑源性神经营养因子（brain-derived neurotrophic factor，）的mRNA表达量，用2−ΔΔCt方法计算各基因mRNA相对表达量。以磷酸甘油醛脱氢酶（reduced glyceraldehyde-phosphate dehydrogenase，）为内参基因，所有实验共进行3次，引物合成信息见表2。

1.11 统计分析

2 结果

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

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

表2 引物信息

图1 丹栀逍遥散成分候选靶标的功能和通路分析气泡图

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

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

2.3 潜在靶点PPI分析

通过对13个靶点进行PPI分析，、、与其他靶标无关联性。其他10个靶标体现出较高的中心度，其中以为主，见图2。故在后续的验证中优先验证、、、、、、、、、。

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

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

图2 13种关键靶点基因的PPI分析

如图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和adj＜0.05阈值的ID有972个，在这阈值下，组2中高表达（log2FC为正）的数目有501个，组1中高表达（log2FC为负）的数目有471个；组间差异基因火山图见图5，热图见图6。

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

GSE4035数据集的差异表达基因与、、、、、、、、、得到3个交集基因，分别是、和，证明了上述分析的可靠性。

2.6 动物实验验证结果

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

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

图4 焦虑数据集主成分分析

上调基因用红色标记，下调基因用蓝色标记，无差异的基因用灰色标记

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

图8 qRT-PCR检测大脑海马组织中ADCY5和GRIN2A的mRNA表达量()

3 讨论

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

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

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

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

丹栀逍遥散化学成分分析，谷氨酸受体包括离子型受体和代谢性受体，离子型谷氨酸受体包括-甲基--天冬氨酸（-methyl--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表达增加磷酸化水平，并调节等基因的表达，从而发挥生物学效应[22]。

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

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

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

3.3 GEO数据库验证结果分析

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

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

3.4 动物实验验证结果分析

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

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

4 结论

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

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

[1] 世界中医药学会联合会, 中华中医药学会. 国际中医临床实践指南焦虑症(2020-10-11) [J]. 世界中医药, 2021, 16(8): 1188-1191.

[2] Yu W, Singh S S, Calhoun S,. Generalized anxiety disorder in urban China: Prevalence, awareness, and disease burden [J]., 2018, 234: 89-96.

[3] Saragih I D, Tonapa S I, Saragih I S,. Global prevalence of mental health problems among healthcare workers during the Covid-19 pandemic: A systematic review and meta-analysis [J]., 2021, 121: 104002.

[4] 刘璇. 后疫情时代的情绪稳定技术[J]. 东方养生, 2022(6): 138-139.

[5] 王剑峰. 琥珀香附宽心汤联合盐酸帕罗西汀治疗肝郁气滞型广泛性焦虑症的疗效研究 [J]. 中国实用医药, 2020, 15(16): 148-149.

[6] 耿东. 基于现代文献及临床调研的广泛性焦虑中医证候及证候要素研究 [D]. 北京: 北京中医药大学, 2012.

[7] 王佳丽, 吴清明. 针刺联合加味丹栀逍遥散治疗气郁化火型广泛性焦虑症的疗效观察 [J]. 中医药临床杂志, 2021, 33(7): 1334-1337.

[8] 施丹. 丹栀逍遥散治疗广泛性焦虑症(障碍)肝郁化火型临床研究 [J]. 中国中医药现代远程教育, 2018, 16(17): 106-108.

[9] Xu H Y, Zhang Y Q, Liu Z M,. ETCM: An encyclopaedia of traditional Chinese medicine [J]., 2019, 47(D1): D976-D982.

[10] Xue R C, Fang Z, Zhang M X,. TCMID: Traditional Chinese Medicine integrative database for herb molecular mechanism analysis [J]., 2013, 41(Database issue): D1089-D1095.

[11] 刘志超, 陈金红, 高晟玮, 等. 基于整合药理学策略的黄芪-丹参抗心肌梗死作用机制研究 [J]. 中国中医基础医学杂志, 2021, 27(10): 1644-1648.

[12] 中国中西医结合学会精神疾病专业委员会. 广泛性焦虑障碍中医证候辨证分型及量化分级标准专家共识 [J]. 中医杂志, 2022, 63(11): 1096-1100.

[13] 刘志超, 高晟玮, 陈金红, 等. 基于整合药理学探究当归四逆汤抗寒凝心脉型心肌梗死的作用机制 [J]. 中成药, 2022, 44(4): 1121-1125.

[14] Davis S, Meltzer P S. GEOquery: A bridge between the gene expression omnibus (GEO) and BioConductor [J]., 2007, 23(14): 1846-1847.

[15] Gentleman R, Carey V J, Zhang J. Meta-data resources and tools in bioconductor [A] //[M]. New York: Springer-Verlag, 2005: 113-133.

[16] Gu Z G, Eils R, Schlesner M. Complex heatmaps reveal patterns and correlations in multidimensional genomic data [J]., 2016, 32(18): 2847-2849.

[17] 吕佳. 焦虑症采用丹栀逍遥散治疗的效果 [J]. 中国医药指南, 2021, 19(32): 30-32.

[18] 王玉, 杨雪, 靳晓杰, 等. 基于中医药整合药理学平台、GEO数据库芯片及分子对接探讨大黄抗肝癌的作用机制 [J]. 中草药, 2020, 51(20): 5207-5219.

[19] 宋文悦, 高明奇, 梁泽怀, 等. 作用于-甲基--天门冬氨酸受体的抗抑郁药物的研究进展 [J]. 中国新药杂志, 2021, 30(15): 1396-1400.

[20] 周珊珊, 刘晓朦, Elona Khasanova, 等. 大麻二酚通过调节脑内AMPA受体表达改善痴呆模型大鼠认知损害 [J]. 中风与神经疾病杂志, 2022, 39(7): 596-600.

[21] 胡勇箭. 代谢型谷氨酸2型受体和γ-氨基丁酸B受体信号通路研究 [D]. 武汉: 华中科技大学, 2018.

[22] Wang W F, Wang Y C, Guo Q T,. Valerian essential oil for treating insomniathe serotonergic synapse pathway [J]., 2022, 9: 927434.

[23] Xu K, Wang M Y, Zhou W,. Chronic-ribose and-mannose overload induce depressive/anxiety-like behavior and spatial memory impairment in mice [J]., 2021, 11: 90.

[24] Pritchett D, Hasan S, Tam S K E,. D-amino acid oxidase knockout (Dao−/−) mice show enhanced short-term memory performance and heightened anxiety, but no sleep or circadian rhythm disruption [J]., 2015, 41(9): 1167-1179.

[25] Comings D E, Muhleman D, Johnson P,. Potential role of the estrogen receptor gene () in anxiety [J]., 1999, 4(4): 374-377.

[26] Sundermann E E, Maki P M, Bishop J R. A review of estrogen receptor alpha gene () polymorphisms, mood, and cognition [J]., 2010, 17(4): 874-886.

[27] Ryan J, Ancelin M L. Polymorphisms of estrogen receptors and risk of depression [J]., 2012, 72(13): 1725-1738.

[28] 陶晶, 仇剑崟, 顾红亮, 等. ESR1基因与5-HTT基因多态性的交互作用对女性围绝经期抑郁障碍的影响 [J]. 上海交通大学学报: 医学版, 2012, 32(6): 766-770.

Discussion on mechanism of Danzhi Xiaoyao Powder in prevention and treatment of anxiety with liver depression transforming into fire syndrome based on integrative pharmacology of traditional Chinese medicine and transcriptomics

PAN Ci-ming1, 2, REN Xue2, 3, WEI Shan-shan1, CAO Bei-ling1, WANG Wei1, YANG Yan1, CHEN Wen-hui1, XU Hong-xi2

1. Yunnan University of Traditional Chinese Medicine, Kunming 650500, China 2. Shanghai University of Traditional Chinese Medicine, Shanghai 201203, China 3. Heilongjiang University of Traditional Chinese Medicine, Harbin 150040, China

To explore the mechanism of Danzhi Xiaoyao Powder (丹栀逍遥散) in the prevention and treatment of anxiety with liver depression transforming into fire syndrome based on integrated pharmacology of traditional Chinese medicine and transcriptomics.The integrated pharmacology research platform of traditional Chinese medicine v2.0 was used to retrieve and obtain the active ingredients and target information of Danzhi Xiaoyao Powder, and the target information of anxiety with liver depression transforming into fire syndrome, and construct the network of traditional Chinese medicine-component-target-pathway. The intersection genes were imported into the String database, and the protein interaction network of intersection targets was constructed to determine the target that needed to be verified. The transcriptomic of GEO database were used to verify the integrated network pharmacology results, and animal experiments were used to verify the results.A total of 13 targets related to anxiety with liver depression transforming into fire syndrome and 430 kinds of active ingredients of Danzhi Xiaoyao Powder were screened out, and the integrated pharmacological results showed that Danzhi Xiaoyao Powder prevented and treated anxiety with liver depression transforming into fire syndrome by mainly regulating the activation of α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid receptor (AMPAR), protein kinase A (PKA),-methyl--aspartic acid (NMDA) receptor, phosphatidylinositide 3-kinases (PI3K) of protein kinase B (Akt) and other signal pathways of anxiety liver depression and fire syndrome. The GEO microarray GSE4035 data set, the key targets of transcriptomics estrogen receptor 1 (), adenylate cyclase 5 () and glutamate ionotropic receptor NMDA type subunit 2A () genes were differential gene expression (DEGs). Quantitative real-time polymerase chain reaction (qRT-PCR) results showed thatmRNA expression was up-regulated andmRNA expression was down-regulated.Danzhi Xiaoxiao Powder may regulate brain-derived neurotrophic factor (BDNF) by up-regulating the expression ofmRNA and down-regulating the expression ofmRNA, thus regulating the binding of BDNF and its receptors and the neurotransmission of glutamate, further improve neuron growth, synaptic conduction and synaptic plasticity, so as to exert its biological effects.

Danzhi Xiaoyao Powder; anxiety; liver depression transforming into fire syndrome; integrative pharmacology; transcriptomics

R285

A

0253 - 2670(2023)12 - 3911 - 10

10.7501/j.issn.0253-2670.2023.12.018

2023-02-23

国家自然科学基金项目（81703964）；云南省科技厅科技计划项目（202101AZ070001-021）；云南省教育厅基金项目：基于GC-NLRP3炎症小体通路研究丹栀逍遥散干预肝郁化火型焦虑行为大鼠的作用机制（2022Y332）；云南省科技厅科技人才与平台计划—云南省中医药学分子生物学重点实验室[2019DG016BK-基于整合网络药理学探究丹栀逍遥散防治郁化火型焦虑行为的作用机制（30172110174），BK-丹栀逍遥散通过GC-NLRP3信号通路干预焦虑行为大鼠的机制研究（30172110179）]

潘赐明（1993—），男，在读硕士研究生，研究方向为中西医结合防治肝胆及肿瘤疾病的应用基础研究。E-mail: panciming@ynutcm.edu.cn

通信作者：韦姗姗（1987—），女，云南昆明人，副教授，博士学位，主要从事证候的生物学基础研究工作。E-mail: weiss814@163.com

[责任编辑 潘明佳]