侯旭粲，成思凡，任兴德，李梦，王耘*

1.北京中医药大学 生命科学学院，北京 102488；2.北京中医药大学 中药学院，北京 102488

抑郁症是一类具有高患病率、高致残率、高自杀率和高复发率等特点的精神障碍疾病。抑郁患者经常陷入悲痛欲绝、自卑及悲观厌世情绪中，严重者可出现幻觉、妄情等精神异常状况，甚至出现自杀行为；而且，其发病情况和自杀事件已经开始呈现低龄化趋势，因此需要密切关注和研究。抑郁症不仅使患者承受精神上的痛苦，还有身体上的折磨，经常在抑郁人群中观察到各种躯体化症状，最常见的包括食欲下降、入睡困难、胸闷心慌、肢体无力和各种疼痛。有研究表明，抑郁症患者躯体化症状的发生率为65.0%～98.2%[1]，一些调查也显示了抑郁症与躯体化症状的紧密联系[2-4]。抑郁症对人们的精神世界和身体健康造成了双重威胁，因此，亟需发现能够有效治疗抑郁症的药物。

解郁丸在逍遥散和甘麦大枣汤的基础上研制而成，由柴胡、当归、白芍、郁金、茯苓、百合、合欢、甘草、小麦和大枣10味中药组成，是常用的抗抑郁方剂。目前，已有许多研究成果表明，解郁丸具有良好的抗抑郁功效。施桂兰等[5]通过实验观察解郁丸对慢性应激抑郁大鼠的行为，以及对下丘脑-垂体-肾上腺轴(hypothalamic-pituitary-adrenal axis，HPA轴)和免疫系统的影响，得到了解郁丸通过降低应激导致的过高浓度的皮质酮、促肾上腺皮质激素、肿瘤坏死因子-α和白细胞介素1β来逆转慢性应激大鼠HPA轴和免疫功能的紊乱，从而发挥抗抑郁功效。马荣等[6]通过小鼠实验，得出解郁丸具有良好的抗焦虑和催眠效果，及解郁丸可能通过调节不同脑区的5-羟色胺(5-hydroxytryptamine，5-HT)、去甲肾上腺素(noradrenaline，NE)等神经递质以发挥抗抑郁作用[7]。宋苗等[8]通过实验观察解郁丸对Wistar-Kyoto(WKY)大鼠抑郁样行为及脑源性神经营养因子(brain-derived neurotrophic factor，BDNF)表达的影响，得出了解郁丸能有效减少WKY大鼠的抑郁样行为，以及通过BDNF这一关键因子发挥抗抑郁作用的结论。除此之外，任素华等[9]通过观察48例轻、中度产后抑郁症患者经解郁丸和氟西汀治疗后的临床疗效和不良反应，得到了两者都具有明显治疗效果，但解郁丸不良反应明显较低的结论。杨静娟等[10]通过观察68例抑郁症患者经解郁丸和马普替林治疗后的临床疗效和不良反应，得到了两者疗效相当，但解郁丸不良反应较轻的结论。综上可知，解郁丸在实验和临床上都具有良好的抗抑郁功效，但由于中药多成分、多靶点的特性，到目前为止，解郁丸发挥抗抑郁功效的具体药理学机制尚不明确。

本研究基于实验室构建的“TCM语法系统”(traditional Chinese medicine，TCM)[11-12]平台和实体语法系统(entity grammar system，EGS)，结合蛋白质表达、细胞类型及其组织和器官的定位数据，获得了化学成分可以作用的抑郁症相关基因、细胞、组织和器官，得到了解郁丸抗抑郁的药理学机制及其化学成分的作用范围。与上述文献中的机制研究相比，本研究得到的结果更加明确，提供了解郁丸发挥抗抑郁效果的成分群和具体途径，这些分析为解郁丸的现代化应用提供了基础[13]。

1 材料和方法

1.1 数据准备

数据来源于课题组建立的多功能平台“TCM语法系统”，它是以实体语法系统为工具，综合各大数据库和文献构建的多功能平台，具体数据库包括：包含中药化学成分信息的中药系统药理学数据库与分析平台(Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform，TCMSP)[14]、中药信息数据[15](Traditional Chinese Medicine Information Database，TCM-ID)等数据库；包含化学成分靶点信息的DrugBank[16]、Stitch[17]等数据库；包含蛋白质相互作用信息的String[18]、蛋白质相互作用数据库[19](Database of Interacting Proteins，DIP)等数据库；包含疾病相关基因信息的在线人类孟德尔遗传数据[20](Online Mendelian Inheritance in Man，OMIM)、治疗靶点数据[21](Therapeutic Target Database，TTD)等数据库。

蛋白质的表达及其表达量的数据来源于文献[22]和人类蛋白质图谱(The Human Protein Atlas，HPA)数据库[23]。文献中的数据是哺乳动物基因组的功能注释5(functional annotation of the mammalian genome，FANTOM5)项目的一部分。此项目专注于哺乳动物基因组的功能注释和转录调控网络的表征，对人体主要器官、多种原代细胞类型、超过200种癌细胞系、30种细胞分化过程和小鼠发育时间过程的RNA样品使用基因表达的帽分析(cap analysis of gene expression，CAGE)方法测序，为进一步分析提供了方便。HPA数据库旨在使用各种组学技术的整合来绘制细胞、组织和器官中的所有人类蛋白质，包括基于抗体的成像，基于质谱的蛋白质组学、转录组学和系统生物学等技术。2个数据库的数据完备，并且基本包含本研究所涉及蛋白质的表达数据，所以使用这2种数据源互相补充完善，以获得可以表达抑郁症相关基因的细胞。

细胞所在组织和器官的数据来源于HPA数据库，用来表征上一部分中所涉及的细胞在组织和器官层面的定位。

1.2 方法

1.2.1获得解郁丸化学成分作用的抑郁症相关基因 解郁丸包括10味中药，每味中药都含有众多化学成分，但有些成分进入体内后不能被有效地吸收利用，也有些成分不适合被开发为药物。因此，本研究根据Lipinski规则和口服生物利用度(oral bioavailability，OB)对解郁丸的化学成分进行初步筛选。具体筛选规则为：分子量<500；氢键供体数目<5；氢键受体数目<10；脂水分配系数<5；OB值≥30%。

所用方法为实体语法系统，实体语法系统是适用于绝大多数复杂系统建模的一种形式语法系统，已经得到了广泛应用和验证[24-26]。实体语法系统用五元组G=(VN，VT，F，P，S) 表示，其中VN和VT指系统中的基本元素，分别代表非末端字符集和末端字符集，V=VN∪VT；F代表基本元素之间的结构关系；P代表一种推理规则，即从已知实体推理出新实体的规则；S代表推理的起点，即系统的初始状态。如果不区分VN和VT，则G=(V，F，P，S)。

V代表网络中各种类型的节点，即中药名称、中药的化学成分和抑郁症相关基因等；F指网络中的各种节点之间的关系，即include或interaction；P为基于上述实体推断解郁丸成分与抑郁症相关基因的作用关系的推理规则，类似于基本的推理公式；S为系统的初始状态，即解郁丸包含的中药、每味中药包含的化学成分、化学成分的靶点、蛋白相互作用信息以及抑郁症相关基因。具体可用下列公式表示：

V=V1∪V2∪V3∪V4∪V5

(1)

V1是解郁丸中所有中药的集合，V2是每味中药包含的所有化学成分的集合，V3是所有化学成分的靶点的集合，V4是蛋白质的集合，V5是抑郁症相关基因的集合。其中V3⊂V4，V5⊂V4，V3和V5之间有交集。这是因为化学成分的靶点和抑郁症的相关基因都是蛋白质，且存在某些化学成分直接作用于抑郁症相关基因。

F={ include1 (A，B，C，D，E，F，G，H，I，J)，
include2 (K，L，T)，link1 (L，M，T)，
link2(M，N，T，O)，dis-target (P，Q，T，R)，
min-distance (L，Q，T，R) }

(2)

在F的集合中，A～J代表解郁丸中的10味中药，A～J∈V1；K代表10味中药中的任意一味，L代表中药K包含的所有成分，L∈V2；T∈{include，interaction }；M代表中药成分的靶点，N代表与M有相互作用关系的蛋白质，M∈V3，N∈V4；O代表M经过O步数作用于N，R的含义与之相同；Q代表抑郁症相关基因，P代表与Q有相互作用关系的蛋白质，P∈V4，Q∈V5；min-distance (L，Q，T，R) 代表从中药成分到疾病基因所需的最小步数。

P=P1∪P2∪P3∪P4，

P1={ include2 (K，L，T)，link1 (L，M，T) ⟹ draw(L，M，T，2) }，

P2={ draw(L，M，T，S)，link2 (M，N，T，O)，O=S+1 ⟹ draw(L，N，T，O) }，

P3={ draw(L，N，T，O)，dis-target (P，Q，T，R)，R=O+1 ⟹ draw(L，Q，T，R) }，

P4={ #min{ F：distance(L，Q，T，X)}=Y，distance(L，Q，T，Y) ⟹ min-distance(L，Q，T，Y) }

(3)

规则P1指从成分靶点数据库中匹配到解郁丸中的所有化学成分的靶点，并将推理步数记为2；规则P2指获得与P1中得到的化学成分靶点具有相互作用关系的蛋白质，并将推理步数加1；规则P3指通过多步推理，最终得到作用于抑郁症相关基因的化学成分，并将推理步数加1；规则P4指获得从成分到抑郁症相关基因的最小推理步数。

S=S1∪S2，

S1={ include1 (A，B，C，D，E，F，G，H，I，J)，include2 (K，L)，link1 (L，M)，link2 (M，N，O) }，

S2={ dis-target (P，Q，T) }

(4)

S1是用于推理的初始信息，包括解郁丸的10味中药、每味中药的化学成分、化学成分的靶点和蛋白相互作用信息。S2是推理的终点，即抑郁症相关基因。通过系统构建及推理，可以得到解郁丸化学成分作用的抑郁症相关基因，最终使用Cytoscape_v 3.7.0软件将解析结果以网络的形式可视化。

1.2.2分析抑郁症相关基因在细胞中的表达情况 在生物体的发育过程中，细胞通过分化形成了形态结构和生理功能不同的细胞，而细胞分化的实质是基因的选择性表达。同一种基因可能在不同类型的细胞中表达，因此不同的细胞会具有相同类型的蛋白质，这种情况在生物体中普遍存在。获得了解郁丸的化学成分作用的抑郁症相关基因后，本研究从FANTOM5项目和HPA数据库中收集了这些基因在不同类型细胞中的表达情况，随后对数据进行了筛选，筛选标准为：基因表达程度为中等表达水平及以上。

1.2.3分析筛选出的细胞的组织和器官定位 细胞是生物体基本的结构和功能单位，不同类型的细胞行使不同的功能，形态相似、功能相同的细胞聚集在一起，并与细胞间质共同构成组织，组织是处于细胞和器官之间的结构，它不一定具备某种特定的功能。几种不同类型的组织经过发育分化并相互结合后，可以构成具有一定形态和特定功能的结构，即器官。虽然组织和器官处于不同层面，但本研究对细胞归类时既涉及组织又涉及器官。这是因为虽然器官具有某种特定的功能，更加适合研究和分析，但对于某些细胞来说，分类到组织层次更合适。例如脂肪细胞和肌肉细胞，它们在机体中普遍存在，分别被划分到脂肪组织和肌肉组织会比划分到某个具体的器官更能体现其功能特性。本研究的分类标准主要参考HPA数据库对人体组织和器官的分类。

2 结果与讨论

2.1 解郁丸化学成分作用的抑郁症相关基因

在“TCM语法系统”中检索到的解郁丸成分共有1293个，其中白芍77个、百合84个、柴胡288个、大枣126个、当归125个、茯苓34个、甘草20个、郁金222个、合欢20个、小麦37个。检索到的抑郁症相关基因共有34个。

通过使用实体语法系统解析解郁丸化学成分所作用的抑郁症相关基因，得到共有181个化学成分可以直接作用于与抑郁症相关的9个基因，分别为5-羟色胺受体2A(HTR2A)、钠依赖性去甲肾上腺素转运蛋白(SLC6A2)、毒蕈碱乙酰胆碱受体M2(CHRM2)、β-2肾上腺素能受体(ADRB2)、胺氧化酶[含黄素]A(MAOA)、钠依赖性5-羟色胺转运蛋白(SLC6A4)、甘氨酸受体亚基α-1(GLRA1)、D(2)多巴胺受体(DRD2)和瞬时受体电位阳离子通道亚家族V成员1(TRPV1)。另外，还有116个化学成分通过多个步骤间接作用于11个抑郁症相关基因，分别为腺苷受体A2a(ADORA2A)、离子型的谷氨酸受体、红藻氨酸4(GRIK4)、伤害感受肽受体(OPRL1)、促肾上腺皮质激素释放因子受体1(CRHR1)、腺苷受体A3(ADORA3)、Ⅰ型神经降压素受体(NTSR1)、神经介素-K受体(TACR3)、物质K受体(TACR2)、环AMP反应元件结合蛋白1(CREB1)、溶质载体家族52-核黄素转运蛋白-成员2(SLC52A2)、神经元PAS结构域蛋白2(NPAS2)。由于中药多成分的特性，某些成分可能存在于多味中药，并在进入体内后发挥协同作用。例如α-cubebol同时存在于郁金和甘草中，进入体内后以效果叠加的形式共同发挥作用。所识别到的297个解郁丸成分共同发挥抗抑郁作用。本研究用图1直观地体现了10味中药、中药的化学成分以及抑郁症相关基因之间的关系。

该网络图直观地显示了解郁丸的化学成分可以作用的抑郁症相关基因及作用关系，解释了解郁丸抗抑郁的药理学作用机制。

关于解郁丸的化学成分靶向的20个抑郁症相关基因，UniProt数据库[27]和GeneCards数据库[28]对它们的功能做出了详尽的解释。HTR2A是5-HT(血清素)的G蛋白偶联受体，能够影响神经活动、感知、认知和情绪，在行为调节中发挥作用。SLC6A2能够通过其高亲和力的钠依赖性再摄取进入突触前末梢终止去甲肾上腺素的作用。CHRM2能够介导各种细胞反应，包括腺苷酸环化酶的抑制、磷酸肌醇的分解和通过G蛋白的作用调节钾离子通道。ADRB2能以较高的亲和力与肾上腺素结合。MAOA在中枢神经系统和外周组织中的神经活性和血管活性胺的代谢中具有重要功能，可以优先氧化生物胺，例如5-HT、NE和肾上腺素。SLC6A4可以将5-HT分子从突触间隙转运回突触前终端以重新利用，从而调节5-HT能信号传导，能够终止5-HT的作用并以钠依赖性方式回收它。GLRA1的异戊二烯通道被甘氨酸激活后，在神经元兴奋性的下调中起重要作用，

注：方形节点代表解郁丸中的10味中药；椭圆节点代表中药的化学成分；V形节点代表抑郁症相关基因；黄色的边代表2个节点之间的相互作用关系；红色的边表示2个节点之间的包含关系。图1 中药-成分-抑郁症相关基因的关系

有助于产生抑制性突触后电流。DRD2为多巴胺受体D2，多巴胺的不足与抑郁症有着密切的联系。TRPV1作为离子型内源性大麻素受体时具有中枢神经调节作用，通过影响AMPA受体的内吞作用，触发海马和伏隔核中内源性大麻素介导的一种长期抑郁症。ADORA2A在许多生物功能中起重要作用，如心律和循环、免疫功能、疼痛调节和睡眠。GRIK4参与L-谷氨酸介导的中枢神经系统兴奋。OPRL1在调节伤害感受和疼痛感觉中发挥作用，对外界应激做出反应。CRHR1在对刺激引起的焦虑的反应中发挥作用。ADORA3参与多种细胞内信号传导途径和生理功能，与神经保护功能相关。NTSR1可以介导神经降压素的多种功能，如低血压、高血糖、体温过低、抗伤害作用和分泌的调节。TACR3可以参与多巴胺代谢过程的调节。TACR2参与乙酰胆碱分泌和神经传递的正调节。CREB1与昼夜节律的同步和脂肪细胞的分化有关。SLC52A2可作为4-羟基丁酸的受体，而4-羟基丁酸是一种神经药物，可用于治疗抑郁症。NPAS2是一种转录激活因子，是形成生物钟的核心组成部分，也是各种生理功能的重要调节因子，包括新陈代谢、睡眠、体温、血压、内分泌和免疫等。

上述20个基因，都与抑郁症有着密切联系，它们大多与神经递质相关，如5-羟色胺、去甲肾上腺素、乙酰胆碱、多巴胺、甘氨酸和谷氨酸。有的与肾上腺素、生物节律和对外界应激的反应相关。这些都在抑郁症的发生发展过程中起着重要作用，解郁丸的化学成分可以通过作用于这些基因综合调节它们介导的生物过程，体现了中药方剂治疗疾病的独特优势和中医遣方用药的整体性思维。

2.2 抑郁症相关基因在细胞中的表达情况

本研究一共分析了18个抑郁症相关基因，缺少靶基因甘氨酸受体亚基α-1(GLRA1)和瞬时受体电位阳离子通道亚家族V成员1(TRPV1)的数据。这是因为GLRA1在HPA数据库中没有检测到表达量，而在FANTOM5项目的数据源中检测到在小脑、胰腺和视网膜中有相对较大的表达量，但是可靠性还未验证，因此本研究不予考虑；TRPV1在mRNA水平上检测到有表达，但2个数据源中均没有其在蛋白水平表达的数据，这种情况是比较常见的，因为从基因到蛋白要经历转录、翻译和翻译后修饰等过程，在这些过程中会有很多因素导致蛋白不能表达出来，而且基因检测和蛋白检测所用方法不同，灵敏度有差异，也会导致上述情况的出现。因此本研究最后分析了18个抑郁症相关基因。

根据方法中设置的筛选条件，如表1所示，得到了可以在中等及中等以上水平表达抑郁症相关基因的细胞。这意味着解郁丸成分进入体内后可以通过靶向抑郁症相关基因而作用于表达这些基因的所有细胞。此部分的研究结果表明，解郁丸中的化学成分可以通过作用于多种类型的细胞，系统地调节抑郁症患者的机体状态。这种治疗方式更加体现了中医药在治疗疾病时的整体性思维。

表1 部分抑郁症相关基因在细胞中的表达情况

2.3 细胞的组织和器官定位

如表2所示，根据方法中设置的分类标准，细胞中抑郁相关基因的表达情况，以及细胞的组织和器官的定位，可以获得在每个组织或器官中表达的抑郁相关基因。

解郁丸成分通过作用于抑郁症相关基因而发挥抗抑郁作用，而靶基因作为蛋白质是可以在多种细胞中表达的，因此，解郁丸成分进入体内后会作用于表达这些基因的细胞，作用效果因靶基因涉及的生化反应及细胞所在的组织和器官的功能的不同而不同。从表2中可以看出，解郁丸成分更多地作用于免疫组织、脂肪组织、心脏和神经系统(包括大脑、小脑、神经细胞和神经干细胞)。

表2 组织和器官中表达的抑郁症相关基因

以抑郁症患者的免疫情况为例,抑郁症常被描述为应激相关疾病，而应激反应主要是以HPA轴的持续激活为核心的心理和生理反应。HPA轴的持续激活会导致体内循环糖皮质激素水平升高，而其浓度过高时会造成脑和其他器官的损伤，使突触数量减少，而且糖皮质激素通常被认为是免疫抑制剂，它能够抑制细胞免疫。另外，压力应激实验表明，肾上腺素的血浆浓度与淋巴细胞和单核细胞的特异性免疫功能呈负相关[29]。据研究，遭受慢性心理应激的人群和抑郁症患者体内的前炎性细胞因子表达水平增加[30-32]，而抗炎性细胞因子表达水平降低[30-33]，过度表达的前炎性细胞因子与5-HT损耗、HPA轴的功能失调有关，并且可以通过调节P38基因和核转录因子-κB(Nuclear factor-kappaB，NFκ-B)信号通路引起神经元的损伤，其浓度变化还会被脑部认为是压力源，从而强化对HPA轴的激活，形成恶性循环。

综上可知，抑郁症患者常伴随着免疫功能低下和前炎性细胞因子表达增加，而这又反过来使抑郁情况加重。鉴于抑郁症的复杂性，单纯进行抗抑郁治疗很难达到理想的效果，因此需要在抗抑郁的同时提高机体的免疫水平。解郁丸的化学成分不仅可以作用于神经系统中的抑郁症相关基因，还有很大一部分作用于机体的免疫组织，通过同时调整抑郁状态和机体免疫功能协同改善抑郁症状，体现了中药方剂治疗复杂性疾病时的独特功效和强大功能。解郁丸成分对其他组织的作用也提示其具有治疗抑郁症与其他疾病的合并症的潜能，例如，抑郁症与肥胖的共病[34-36]，抑郁症与心脏病的共病[37-39]等。

由于从抑郁症相关基因到细胞定位再到组织和器官定位的过程比较复杂，为了方便查看，该研究通过图的形式可视化三者之间的关系(见图2)，并提取子网(见图3)作为简要描述。

注：V形节点代表与抑郁症相关的基因；椭圆形节点代表可以表达这些基因并满足筛选标准的细胞；三角形节点代表这些细胞所在的组织或器官；橙色的边表示2个节点之间的表达关系；紫色的边表示2个节点之间的定位关系；下同。图2 靶基因-细胞-组织/器官的关系

图3 靶基因-细胞-组织/器官关系图的子网络

图3直观地说明了基因可以同时在多个细胞中表达的情况，并且存在多个功能相似的细胞属于同一组织或器官的情况。子网清楚地体现了解郁丸的化学成分可以作用的基因、细胞、组织和器官，明确了它们在人体中的作用范围。

3 总结与展望

本研究以实体语法系统为框架，综合利用信息学、中药学和生物学数据及工具，通过推理得到了可以靶向抑郁症相关基因的297个解郁丸成分，并得出了解郁丸抗抑郁的药理学机制，即这些化学成分可以直接或间接地作用于抑郁症相关基因。随后根据蛋白表达数据筛选了可以表达上述抑郁症相关基因的细胞，最后得到这些细胞的组织和器官定位，从而明确了解郁丸化学成分在人体中的作用范围。根据解郁丸成分影响的组织和器官，提出了解郁丸在抗抑郁的同时也能够调整机体的免疫功能，并且具有治疗一些抑郁症与其他疾病的合并症的潜能，为解郁丸的进一步研究提供了方向。

今后的研究可以从本研究中提供的作用于抑郁症相关基因的化学成分入手，通过对成分各方面性质的评价，筛选出可以作为抗抑郁新药的先导化合物的成分。另外，根据解郁丸成分作用的组织和器官，可以为一些抑郁症与其他疾病共病的特殊人群开发个性化治疗药物，例如，可以筛选出同时作用于脂肪组织和抑郁症相关基因的解郁丸成分，从而获得治疗抑郁合并肥胖疾病的成分群，然后可以针对这些成分进行后续的筛选和药物设计。总之，本研究结果为解郁丸的应用和之后药物开发提供了依据。