新型冠状病毒肺炎（corona virus disease 2019,COVID-19），简称“新冠肺炎”，已发展为全球性流行病，目前还未有对COVID-19的特效药，因此研究预防或治疗新冠肺炎的有效药物，迫在眉睫。芳香中药自古以来在防治瘟疫的抗争中发挥着重要作用，芳香疗法是中医药疗法的重要组成部分，《景岳全书》卷十三《瘟疫》论述避疫法时，就提到“用降真香烧焚，大解邪秽，小儿带之，能解诸邪，最验”[1]。芳香中药在现代防治病毒性疾病中也发挥着重要作用，如鱼腥草、野菊花等，多具有清热解毒、芳香化湿功效，这些药材大多有抗病毒、抗菌、消炎解热及调节免疫功能等作用，尤以抗病毒作用最为显著。

经过整合芳香疗法行业专家组综合文献研究、专家共识等研究结果提出一套用于防疫新冠肺炎的芳香中药挥发油复方，通过香薰、喷洒、吸嗅等方式运用于江西省中医院防治新冠肺炎全过程，病人经鼻腔吸入后感觉良好，从疫情初期江西省中医院到该院确诊人数清零，医护人员零感染，该挥发油作用功不可没。该复方作为天然药物产物，来源丰富，毒副作用小，适宜性佳，在变异性强、传染率高的病毒性疾病的防治中显示出独特优势。本文结合现代网络药理学和分子对接技术对芳香中药挥发油在预防及治疗方面的作用机制进行预测和探析，为研究探索新冠肺炎的预防治疗提供现代医学理论基础。

芳香中药挥发油由10味药协同发挥扶正驱秽、芳香化湿、疏肝解郁之功。其中，4个核心药味挥发油，肉桂挥发油、丁香挥发油、荆芥挥发油、香薷挥发油均有抗菌抗病毒作用，对肺部感染有较好疗效；此外，还有镇静、降温、提高免疫力等辅助功效，对治疗新冠肺炎发挥主要治疗作用。月桂挥发油、雪松挥发油、橙花挥发油、迷迭香挥发油、真薰衣草挥发油、八角茴香挥发油除大都能抗菌外，部分具有抗病毒作用，且具有抗焦虑、抗炎等作用，可以帮助核心药味进行治疗，预防其它感染，并改善患者情绪，有助于提高患者依从性。

1、材料与方法

1.1 基于文献获取复方中各药挥发油成分

通过cnki(https：//www-cnki-net-443.webvpn.jxutcm.edu.cn/）和pubmed(https：//www.ncbi.nlm.nih.gov/pubmed/）分别检索“肉桂挥发油”[2,3,4]“丁香挥发油”[5,6,7]“荆芥挥发油”[8,9,10,11]“香薷挥发油”[12,13]“月桂挥发油”[14,15]“雪松挥发油”[16]“橙花挥发油”[17]“迷迭香挥发油”[18,19]“真薰衣草挥发油”[20,21,22,23]“八角茴香挥发油”[24,25]文献，整理文献中所指各个挥发油中含量较高化学成分。

1.2获取中药挥发油成分靶点和新型冠状病毒肺炎靶点及潜在作用靶点预测

通过Swiss Target Prediction(http：//www.swisstargetprediction.ch/）在线平台收集挥发油成分的靶点并加以整理。通过NCBI(https：//www.ncbi.nlm.nih.gov/）、GenCards(https：//www.genecards.org/）等数据库分别收集新型冠状病毒肺炎相关靶点。将挥发油处方中10味芳香药物的挥发油成分靶点与新冠肺炎潜在靶点共同导入ImageGP(http：//www.ehbio.com/）平台，运用韦恩图功能映射取交集，得到的交集靶点，即为普济避瘟扶正熏香挥发油防疫新冠肺炎潜在作用靶点。

1.3交集靶点的基因本体（GO）及京都基因与基因组百科全书（KEGG）分析

在R软件（Version3.5.2)安装Bioconductor软件包“org.Hs.eg.db”并运行，将关键交集靶点转换成entrez ID。然后在R软件安装“cluster Profiler”包，根据已转换的entrez ID，以P<0.05，进行关键靶基因GO与KEGG功能富集分析，并将结果输出。

1.4 构建中药-挥发油成分-靶点-关键通路关系网络

采用Cytoscape3.7.2软件对上述中药，交集靶点，作用与的交集靶点挥发油成分，以及重点通路构建网络PPI图，在网络中，Degree值代表某一节点与其它节点的相关数目，Degree值越大，则越多节点与该节点相关联。

1.5 基于人体血管紧张素转换酶II(ACE2）的分子对接

文献报道，SARS-CoV-2与严重急性呼吸综合征冠状病毒（SARS-CoV，非典）两种病毒感染细胞的途径均通过其表达的S-蛋白与人体内血管紧张素转换酶II(ACE2）的结合，导致病毒侵入机体而致病[26]。因此从PDB(http：//www.rcsb.org/）数据库下载ACE2(PDB ID:1R42）蛋白质3D结构作为分子对接靶标。通过Discovery Studio软件，采用基于CHARMm的半柔性对接程序CDOCKER对挥发油成分和ACE2进行分子对接，且将结果与已注册用于新冠肺炎治疗的临床试验化学药[27]瑞德西韦（remdesivir）、阿比朵尔（arbidol）、氯喹（chloroquine）、恩曲他滨（emtricitabine）、达芦那韦（darunavir）和ACE2的对接结果进行比较。

2、结果

2.1 芳香中药挥发油成分

基于1.1中所述方法，得到芳香中药挥发油成分共131个。见表1。

表1芳香中药复方中各挥发油含量较高成分

2.2 芳香中药复方挥发油和新冠肺炎靶点及复方潜在作用靶点预测

通过“1.2”所述方法，得到丁香挥发油成分522个靶点，荆芥428个，肉桂213个，香薷662个，八角茴香614，橙花212个，迷迭香287个，雪松274个，月桂391个，真薰衣草450个。得到新型冠状病毒肺炎相关靶点348个。得到关键交集靶点，即普济避瘟扶正熏香挥发油防疫新冠肺炎潜在作用靶点65个，如图1。

图1芳香中药复方挥发油和新型冠状病毒肺炎靶点韦恩图

2.3 交集靶点的GO分析及KEGG通路富集分析

为进一步探讨上述65个交集靶点在防疫新冠肺炎时的作用机制，通过Rstudio对65个的靶点进行GO富集分析，见图2。GO-BP分析共得到1327条生物学过程，选前20条生物学过程。主要涉及对脂多糖的反应（Response to lipopolysaccharide）、对细菌起源分子的反应（Response to molecule of bacterial origin）、细胞对生物刺激的反应（Cellular response to biotic stimulus）、细胞对脂多糖的反应（Cellular response to lipopolysaccharide）、白细胞迁移（Leukocyte migration）、细胞对细菌起源分子的反应（Cellular response to molecule of bacterial origin）、血小板活化（Platelet activation）、凋亡信号通路的调节（Regulation of apoptotic signaling pathway）、凝血（Blood coagulation）、止血（Hemostasis）等。GO-CC分析共得到27个细胞组成，根据P<0.01筛选出20个结果。主要包括膜筏（Membrane raft）、膜微区（Membrane microdomain）、膜区（Membrane region）、囊泡腔（Vesicle lumen）、细胞质囊泡腔（Cytoplasmic vesicle lumen）、粘着斑（Focal adhesion）、细胞-基质粘附连接（Cell-substrate adherens junction）、细胞-底物连接（Cell-substrate junction）等。GO-MF分析共得到94条有关分子功能的结果，取结果前20。主要涉及蛋白质丝氨酸/苏氨酸激酶活性（Protein serine/threonine kinase activity）、磷酸酶结合（Phosphatase binding）、蛋白磷酸酶结合（Protein phosphatase binding）、MAP激酶活性（MAP kinase activity）、磷脂酰肌醇双磷酸激酶活性（Phosphatidylinositol bisphosphate kinase activity）、蛋白激酶C活性（Protein kinase C activity）等。

图2交集靶点的GO富集分析

KEGG分析结果显示共有161条通路途径，取结果前20条通路作为重点通路，可以发现糖尿病并发症中的AGE-RAGE信号通路（AGE-RAGE signaling pathway in diabetic complications）、卡波西氏肉瘤相关疱疹病毒感染（Kaposi sarcomaassociated herpesvirus infection）、细胞凋亡（Apoptosis）、甲型流感（Influenza A）、百日咳（Pertussis）等通路与疾病发展有关。为了反应KEGG术语之间的关系，使用ClusterProfiler包对前生物学功能进行重构，图中更大的红色节点代表更明显的途径。我们可以看到在这些通路途径中CXCL8、NFKB1、PIK3CA、PIK3CB、PIK3CD、MAPK1、MAPK3、RELA等靶点比较重要，它们可以同时在多种通路中发挥作用，如图3。靶点与主要通路关系图中靶点代码见表2。

表2靶点代码对照表

图3 KEGG富集分析

2.4 中药-挥发油成分-靶点-关键通路关系网络

通过1.4所述方法，得到中药-挥发油成分-靶点-关键通路关系网络图，如图4。分析显示Degree>10的挥发油成分有42个，见表3。Degree>30的靶点共有13个，分别是PTGS1、JAK1、MAPK14、NOS2、PARP1和MAPK3等；Degree>18的通路有4条，分别为卡波西氏肉瘤相关疱疹病毒感染（Kaposi sarcoma-associated herpesvirus infection）、人巨细胞病毒感染（Human cytomegalovirus infection）、糖尿病并发症中的AGE-RAGE信号通路（AGE-RAGE signaling pathway in diabetic complications）、细胞凋亡（Apoptosis）。网络关系图表明以上结果是芳香中药复方挥发油防疫新冠肺炎的主要挥发油成分，靶点和通路。

图4中药-挥发油成分-靶点-关键通路关系网络图（红色-中药，绿色-挥发油成分，蓝色-靶点，黄色-通路注：挥发油成分以CAS号替代）

表3 Degree>10的挥发油成分

2.5 挥发油成分-靶标分子对接

选取Degree>10的42个挥发油成分作为对接成分。通过1.5所述方法，将挥发油成分、临床试验化学药与ACE2进行分子对接，结果见表4，分子对接模式图见图5。

图5分子对接模式图

表4挥发油成分和部分化学药与ACE2水解酶的结合能

一般认为，当化合物以配位与疾病受体结合时，配体与受体结合的构象稳定时能量越低，发生的作用可能性越大[28]。分子对接结果显示较多挥发油成分和ACE2的结合能低于-20 kJ•mol-1，表明这些挥发油成分与ACE2有较好的结合活性；同时该方中部分挥发油成分和ACE2的结合能可与化学药相媲美甚至超过化学药，且还可以发挥出单一化学药成分所没有的有助于防治新冠肺炎的作用。

3、讨论

自疫情发生以来，中医药在预防疫病方面积累了大量宝贵经验，除了在消灭传染源、切断传播途径方面的措施外，以芳香药用植物挥发油为原料，可用于预防疾病，杀菌除秽，还可以预防呼吸道疾病，减轻焦虑，促进中枢神经兴奋[29]。所以本研究采用多味中药挥发油组成的方剂应对新冠肺炎，探索其抗病毒、抗炎、祛痰等功效，为芳香疗法在防疫新冠肺炎疾病全过程发挥着积极作用理论提供依据，也为共同防疫提供新思路。

预防期，芳香中药复方可提高免疫，抵抗病毒外侵。从KEGG富集分析结果可知，HIF-1信号通路（HIF-1 signaling pathway）可能是芳香中药复方干涉新型冠状病毒肺炎涉及的主要通路之一。王明亮等[30]研究调查发现致肺癌途径包含HIF-1α等，可通过HIF-1通路促进非小细胞肺癌细胞增殖和侵袭，抑制细胞凋亡。故推测本方可能通过抑制HIF-1通路，帮助受损细胞正常凋亡来控制疾病的蔓延和发展。根据成分-靶点-信号通路网络图，NOS2可能是重要靶点之一，Scott B.Drutman等[31]研究结果显示NOS2缺乏症是患者感染巨细胞病毒（CMV）致死原因，因此可以推测NOS2在人体抗病毒机制中发挥重要作用。

一般认为，当化合物以配位与疾病受体结合时，配体与受体结合的构象稳定时能量越低，发生的作用可能性越大[28]。分子对接结果显示较多挥发油成分和ACE2的结合能低于-20 kJ•mol-1，表明这些挥发油成分与ACE2有较好的结合活性；同时该方中部分挥发油成分和ACE2的结合能可与化学药相媲美甚至超过化学药，且还可以发挥出单一化学药成分所没有的有助于防治新冠肺炎的作用。

表4挥发油成分和部分化学药与ACE2水解酶的结合能

治疗期，芳香中药复方可清肺祛痰，消除炎症外侵。富集分析显示多条通路主要富集在抗病毒感染和修复肺部损伤等功能上，起到抗病毒抗炎作用。其中AGEs-RAGE信号通路能激活多种“损伤反应基因”的多效转录调节因子NF-κB（核转录因子-κB），进入细胞核，可引起大量黏附分子、生长因子、促炎细胞因子的表达与释放[32]，故本方可能可以通过作用该通路上的靶点抑制AGEsRAGE信号通路表达，从而抑制炎症的发生。人巨细胞病毒属于β亚科的DNA疱疹病毒，肺脏为人巨细胞病毒最易感染的器官，容易引起肺炎。在病理生理学中，引起肺炎的主要原因是细支气管炎尤其是毛细支气管周围炎，易使小支气管粘膜壁的充血和水肿，常常可导致细支气管部分性或者完全性梗阻[33]。

从成分-靶点-信号通路网络图分析，重要靶标分别为：PTGS1、JAK1、NOS2、MAPK14、PARP1、EGFR等。PTGS也称为环氧合酶（COX），是前列腺素生物合成中的关键酶，COX-1通常存在于身体的各个区域，不仅包括胃，而且包括任何炎症部位，COX-1将花生四烯酸转化为产物LXs,LX结构它们在体内均显示出强大的抗炎和促分解作用[34]。MAPK14可在炎症状态下调控促炎因子的产生，可协调细胞对压力和炎症的反应，这在调节神经元细胞死亡中起着重要作用[35,36]。由此可见芳香中药复方挥发油可能通过调控促炎、抗炎因子的产生，两种因子相互协调达到保护机体的作用，靶点EGFR存在于通路Human cytomegalovirus infection、HIF-1 signaling pathway、PD-L1 expression and PD-1 checkpoint pathway in cancer中，与靶点MUC1有关，李洪晓等[37]通过研究发现患者痰液MUC1、Muc5ac表达与气道炎症相关，可以通过消除炎症的发生达到祛痰的效果，因此，EGFR可能通过MUC1对新冠病毒感染后肺部产生的痰瘀阻肺发挥积极的疗效。

KEGG分析发现靶点ACE主要通过恰加斯病（美国锥虫病）[Chagas disease(american trypanosomiasis）]，肾素-血管紧张素系统（Reninangiotensin system），肾素分泌（Renin secretion）等通路途径发挥作用；ACE参与催化血管紧张素I转化为生理活性肽血管紧张素II的酶。而新型冠状病毒（SARS-Co V-2）结合血管紧张素转化酶2(ACE2）受体激活经典的肾素-血管紧张素系统（RAS）调节途径，作用于肺及肺外靶器官，引发多脏器损伤[38]。因此，将挥发油成分与ACE2进行分子对接，结果显示棕榈酸（Palmitic acid）、亚油酸（Linoleic acid）、茴香酸（Anisic acid）与ACE2的结合能比目前推荐使用的化学药瑞德西韦（remdesivir）等的结合能更有优势，挥发油成分是否可以通过作用于ACE2受体调节多条信号通路起到治疗COVID-19的作用有待临床一线医疗工作者的进一步验证。

康复期，保护神经，修复损伤。芳香中药挥发油的生物学过程不仅富集在免疫系统、内分泌系统、信号转导等的一系列生物过程信号通路上，用于对抗SARS-CoV-2感染导致强烈的免疫反应和炎症风暴，因此大量细胞因子被激活，同时还富集在神经系统上具有保护神经，调节情志等功能，能够缓解疫情带来的抑郁焦躁情绪。JAK1、PARP1可能是本方干预新冠肺炎涉及的重要靶点，有文献表明JAK1的mRNA水平与疾病的预后和免疫浸润水平相关[37]，而PARP1充当检测DNA损伤的第一响应者，与几种DNA修复过程的调控有关，可通过调节染色质结构以及与多种DNA修复因子相互作用和修饰而有助于DNA修复[38]。综上，芳香中药挥发油可通过多靶点共同协助在新冠肺炎康复期发挥积极作用。

目前新冠肺炎感染人数虽然已基本得到控制，但是我们尚未研发出特效药来全面治疗新冠肺炎，所以现在我们最关键的任务就是能够提高治愈率，降低致死率，加强治愈患者的愈后调理作用。建议是否可以综合现阶段治愈COVID-19患者的中药处方，筛选得到与冠状病毒感染中的关键蛋白结合较好的的芳香中药复方挥发油成分，根据辨证论治考虑其所在芳香中药，建立一些普适的芳香中药复方或精油处方，而芳香中药精油处方不仅在预防和治疗期对新型冠状病毒引起的肺部炎症有着干预改善作用。而且用药方便简单易行，经济实惠，同时能提高患者后期生活质量，再结合推荐的其他药物在新冠肺炎全过程进行综合治疗，为预防治疗新型冠状病毒肺炎及恢复治疗提供新的治疗方法和研究方向。

参考文献：

[1]王文远.古代中国防疫思想与方法及其现代应用研究[D].南京:南京中医药大学,2011.

[2]张小飞,冯玲玲,伍振峰,等.四川产肉桂挥发油化学成分分析及药效学研究[J].中国医药工业杂志,2016,47(9):1183-1187.

[3]吴修富.肉桂提取物的主要化学成分及药理活性研究进展[J].中国药房,2015,26(24):3454-3456.

[4]侯小涛,郝二伟,秦健峰,等.肉桂的化学成分、药理作用及质量标志物(Q-marker)的预测分析[J].中草药,2018,49(1):20-34.

[5]美丽,朱懿敏,罗晶,等.丁香化学成分、药效及临床应用研究进展[J].中国实验方剂学杂志,2019,25(15):222-227.

[6]牛彪,金川,梁剑平,等.牛至、香茅、丁香精油化学成分及体外抑菌活性研究[J].食品研究与开发,2020,41(3):46-52.

[7]李京晶,籍保平,周峰,等.丁香和肉桂挥发油的提取、主要成分测定及其抗菌活性研究[J].食品科学,2006,27(8):64-68.

[8]王凤,温桃群,桑文涛,等.荆芥挥发油化学成分及药理作用研究现状[J].中南药学,2017,15(3):312-318.

[9]袁岸,刘淇,饶志粒,等.基于网络药理学的荆芥挥发油主要成分抗炎机制研究[J].中国药理学通报,2020,36(1):97-103.

[10]邓寒霜,李筱玲.土荆芥挥发油成分及其活性研究进展[J].陕西农业科学,2018,64(8):93-96.

[11]黄晓巍,刘玥欣,刘轶蔷,等.荆芥化学成分及药理作用研究进展[J].吉林中医药,2017,37(8):817-819.

[12]金哲,马林,崔成哲,等.东北野生香薷挥发油香味成分分析及在卷烟中的应用[J].中国农学通报,2008,24(7):93-96.

[13]马国娟,李晖,陈勇,等.芳香解表类中药香薷、薄荷、荆芥穗、防风中挥发性成分的气质联用分析[J].江苏中医药,2009,41(2):57-59.

[14]苗红林,翁新楚.月桂叶香气成分的初步分析[J].上海大学学报(自然科学版),2009,15(3):326-330.

[15]林正奎,华映芳,龚国萍,等.月桂叶精油成分及逐月动态变化[J].JINTEGRPLANTBIOL,1990(11):878-882.

[16]扶巧梅,彭映辉,熊国红,等.两种松科植物精油对蚊虫的熏杀活性及其化学成分分析[J].中国生物防治学报,2013,29(3):370-375.

[17]饶先立,郭宏霞,孙胜南.巴西苦橙精油和苦橙花精油的挥发性成分分析及卷烟应用效果[J].化工技术与开发,2019,48(4):34-39.

[18]徐勇,姚雷,张艳玲,等.三种迷迭香植物学性状和精油成分研究[J].上海交通大学学报(农业科学版),2006,24(5):429-434.

[19]周永红,纪良霞.广西产迷迭香挥发油化学成分的分析[J].林产化工通讯,2004,38(6):34-36.

[20]万传星,朱丽莉,刘文杰.薰衣草精油化学成分及抗菌活性研究(英文)[J].塔里木大学学报,2008,20(2):40-43.

[21]胡喜兰,韩照祥,刘玉芬,等.薰衣草挥发油的抗氧化活性测定和成分比较[J].食品科技,2006(9):115-117.

[22]张健,蔡宝国,章苏宁,等.薰衣草精油化学成分的GC-MS分析比较[J].食品工业,2007(5):52-54.

[23]王新玲,娜·卡斯木,胡君萍,等.薰衣草不同部位中挥发油化学成分的比较[J].华西药学杂志,2010,25(3):361-362.

[24]杨靖,李瑞丽,陈芝飞,等.八角茴香油的超临界CO_2萃取及分析[J].中国调味品,2010,35(12):96-98.

[25]权美平.八角茴香精油的成分分析及生物活性研究进展[J].中国调味品,2017,42(1):164-166.

[26]宗阳,丁美林,马世堂,等.以血管紧张素转换酶Ⅱ(ACE2)为受体挖掘治疗新型冠状病毒肺炎(COVID-19)潜在中药及单体成分[J].中草药,2020,51(5):1123-1129.

[27]张超,陈姝冰,张洁,等.浅析注册用于新冠肺炎治疗的临床试验药物[J].药学学报,2020,55(3):355-365.

[28]徐森楠,庄莉,翟园园,等.基于网络药理学研究二至丸防治骨质疏松症的物质基础与作用机制[J].中国药学杂志,2018,53(22):1913-1920.

[29]芳香疗法有讲究[J].国内外香化信息,2003(3):10.

[32]杨超茅,杨志新,马晓玲.AGEs-RAGE信号通路在糖尿病肾病中的作用机制及中医药研究进展[J].中医学报,2019,34(9):1864-1868.

[33]赵爱玲,王雪琴,王江涛.人巨细胞病毒感染与小婴儿肺炎的临床分析[J].医药论坛杂志,2019,40(5):37-38.

[36]刘杰,刘广学,尚明英,等.华细辛和北细辛HPLC特征图谱识别与抗炎靶点及抗炎成分分析[J].中国中药杂志,2020,45(6):1374-1383.

罗晶,殷宁,黄小英,张华,李慧婷,任桂林,柳小莉,杨明.芳香中药防治新型冠状病毒肺炎作用机制探究[J].江西中医药大学学报,2020(03):66-73.