何鈺，許海燕，單進軍，趙霞，沈存思，紀建建，林麗麗，楊瑞，謝彤，汪受傳

(1.南京中醫(yī)藥大學，江蘇省兒童呼吸疾病(中醫(yī)藥)重點實驗室，江蘇 南京 210023；2.南通市疾病預防控制中心微生物檢驗科，江蘇 南通 226007)

新型冠狀病毒(SARS-CoV-2)感染誘發(fā)的肺炎(COVID-19)引發(fā)嚴重的全球性公共衛(wèi)生事件。截至北京時間2020年4月2日13:00，全球NCP累積確診病例937 567例，死亡47 256例，涉及180個國家和地區(qū)[1]。

中國成功遏制了疫情蔓延，取得階段性成果，為國際社會提供重要借鑒，其中中醫(yī)藥發(fā)揮的力量不可小覷，特別是阻止病情繼續(xù)擴大以及縮短疾病治療時程等方面方面發(fā)揮了極其獨特的積極作用。金欣口服液是在麻杏石甘湯的基礎上加減所得，大量研究結果顯示其治療病毒性肺炎效果顯著[2-3]。中藥方劑成分復雜，具有多靶多效特點[4]，注重從整體觀上對病證辨證施治。網絡藥理學技術選取信號節(jié)點進行多靶點藥物分子設計，與中醫(yī)藥多靶點、整體效應相符合[5-7]。本研究旨在利用網絡藥理學探究金欣口服液治療COVID-19的可能性和潛在靶點，為金欣口服液更廣泛地應用于臨床提供理論依據。

1 方法

1.1 篩選類藥性成分及其作用靶點

基于中藥系統(tǒng)藥理數據庫和分析平臺(TCMSP，http://lsp.nwu.edu.cn/tcmsp.php)查詢配伍藥物麻黃(Ephedra herba)、苦杏仁(Amygdalus communis vas)、黃芩(Scutellariae radix)、葶藶子(Lepidii semen descurainiae semen)、桑白皮(Mori cortex)、前胡(Peucedani radix)和虎杖(Polygoni cuspidati rhizoma et radix)的物質組成分，以口服生物利用度(Oral bioavailability,OB)≥30%、成藥相似性(Drug likeindex,DL)≥0.18為標準篩選合適的類藥成分。通過PubChem獲取化學結構標準的SMILES信息以及InChIKeys信息。以ClassyFire化學結構分類器(https://cfb.fiehnlab.ucdavis.edu)對得到的化學結構進行結構分類。數據可視化通過sunburst R包實現。

1.2 疾病對應靶點的收集

疾病對應的靶點通過DisGeNET數據庫(https://www.disgenet.org/)收集。以疾病對應的MeSH標準主題詞為關鍵詞，在數據庫中搜索疾病對應的靶點。所選擇MeSH標準主題詞及其對應的疾病ID號分別為：嚴重急性呼吸綜合征C1175175(Severe acute respiratory syndrome,SARS)、中東呼吸癥候群冠狀病毒感染癥C3694279(Middle east respiratory syndrome,MERS)、呼吸道合胞病毒C0035235(Respiratory syncytial virus infection,RSV)和肺炎C0032285(Pneumonia)。

1.3 潛在治療靶點互作網絡的構建

所有的靶點基因和靶點蛋白通過Uniport(https://www.uniprot.org/)數據庫進行注釋標準化。將金欣口服液中類藥性成分相關的靶點和疾病相關的靶點進行匹配，篩選出具有交互性的治療靶點基因。將篩選得到的靶點基因輸入STRING數據庫(https://string-db.org/)進行蛋白質互作分析，構建蛋白質相互作用網絡圖。數據可視化通過venn R包實現。

1.4 基因本體分析和代謝通路分析

篩選得到的基因靶點輸入DAVID數據庫(https://david.ncifcrf.gov/)進行差異基因的基因本體功能注釋和代謝通路分析，特征地識別和分析潛在靶點。從結果中挖掘FDR值小于0.05，且基因本體注釋語句與肺炎直接相關的基因本體和代謝通路。構建化合物→基因→基因本體通路的網絡構建，通過CytoScape軟件(https://cytoscape.org)實現。

2 結果

2.1 金欣口服液的效應成分分析

根據TCMSP數據庫查詢到金欣口服液成方中七種中藥的物質組成分共962種，以OB≥30%且DL≥0.18為篩選標準得到129個類藥性成分。韋恩圖的結果表明，金欣口服液中類藥性的成分主要集中在前胡、桑白皮、麻黃以及苦杏仁中(圖1A～B)?；瘜W結構分類結果表明，類藥性成分主要分布在黃酮大類和脂類成分，黃酮大類包括黃烷酮類、黃酮類以及甲基化黃酮類，脂質類成分主要包括固醇及甾體類、異戊烯醇脂、脂肪酰類(圖1C～D)。篩選的具體結果見表1。

表1 金欣口服液物質成分化學分類結果

2.2 金欣口服液治療新型冠狀病毒的潛在靶點

通過TCMSP數據庫篩選出類藥成分對應的靶點共268個。新型冠狀病毒是一種單鏈RNA正鏈包膜β冠狀病毒，與SARS、MERS相似，基因組序列分別具有79%和50%的同源性[8]，因此選擇SARS和MERS所對應的疾病為關鍵詞在DisGeNET搜索疾病相關靶點。由于金欣口服液在治療RSV肺炎方面療效顯著，因此我們將RSV肺炎也作為關鍵詞納入篩選標準。共挖掘到疾病相關靶點基因805個。整合類藥成分相關靶點基因和疾病相關基因，挖掘出金欣口服液可能用于治療新型冠狀病毒肺炎的潛在靶點共75個(圖2)。

2.3 潛在治療靶點的蛋白質互作分析

75個靶點基因輸入數據庫，以互作得分>7.0作為篩選標準進行蛋白質互作分析。共得到75個節(jié)點和527條相互作用線，蛋白質相互作用富集P<0.01。采用K-means聚類算法對蛋白質相互作用網絡進行聚類分析，聚類結果以Uniport關鍵詞和Reactome信號通路注釋條目表示(圖3，表2)。

2.4 基因本體分析和KEGG通路分析

將靶基因輸入DAVIDA數據庫進行基因本體分析和代謝通路分析。以富集的FDR值<0.05篩選基因本體(Gene ontology,GO)數據和KEGG代謝通路數據。共挖掘到生物過程(Biological process,BP)注釋條目45個、細胞組成(Cellular component,CC)條目12個、分子功能條目(Molecular function,MF)19個、KEGG代謝通路分析(KEGG pathway)條目66個，從中挖掘與抗病毒免疫以及炎癥應答相關的注釋條目繪制氣泡圖(圖4)，并進一步構建化合物→基因→富集通路圖(圖5)。

3 討論

新型冠狀病毒肺炎在中醫(yī)理論上屬于疫病、瘟疫的范疇[9]，臨床上以清熱、解毒、化濕為主。金欣口服液在麻杏石甘湯的基礎上加減所得，全方宣肺開閉、清熱解毒、化痰止咳、除濕,治療病毒性肺炎有顯著療效[10]。本研究利用網絡藥理學共篩選出金欣口服液中129種類藥性成分，包括苯丙素類83種、脂類成分27種、苯環(huán)型化合物6種、有機氧化合物7種等。報道指出天然植物中的黃酮類成分高效低毒，具有良好的體內外抗多種病毒活性，能有效預防和治療病毒性感染[11]。

篩選到的類藥成分-疾病交互的靶點75個，進行進一步的蛋白質互作的關聯網絡分析。根據相似度算法可將關聯網絡中各個節(jié)點聚類成不同社區(qū)(Community)。實驗結果顯示：炎癥因子相關的靶點聚類在同一個社區(qū)內，如IL-1A、IL-1B、IL-10等(圖3中紅色靶點)，信號轉導通路相關靶點聚集在同一個社區(qū)內，如MAPK和JUN等(圖3中藍色靶點)。對蛋白質互作靶點進行富集分析，Uniport關鍵詞富集分析發(fā)現，互作主要富集在細胞因子、蛋白質-病毒相互作用、炎癥應答等。Reactome信號通路數據庫收集了關于蛋白質相互作用的信息。Reactome信號通路富集結果顯示，靶點基因主要富集到白細胞介素和免疫炎癥因子通路上，表明金欣口服液中的類藥成分有可能通過抗病毒免疫發(fā)揮抗COVID-19的作用。

利用75個靶點基因進行GO分析和KEGG代謝通路分析，發(fā)現基因靶點富集在抗病毒免疫以及抗炎的生物過程和代謝通路上。如KEGG代謝通路富集在T細胞受體細胞通路、B細胞信號受體通路、Toll樣受體信號通路、NOD樣信號受體通路等。BP條目主要富集在炎癥響應通路、氧化還原過程、NF-κB信號通路、MAPK信號通路等。此外，富集到G-蛋白偶聯受體生物與趨化因子表達相關。肺部在受到病毒侵襲時，往往會激發(fā)機體的固有免疫或獲得性免疫[12-15]。本研究結果顯示，金欣口服液中的類藥性成分可能通過發(fā)揮抗病毒免疫和抗炎作用，從而具有潛在的治療COVID-19的療效。

通過化合物→基因→通路注釋關聯網絡分析，節(jié)點度最高的兩個靶點基因P35354和P23219分別對應于為前列腺素G/H合酶2和前列腺素G/H合酶1，節(jié)點度最高的兩個類藥化合物MOL000098、MOL000422分別對應于槲皮素和山奈酚，由此說明金欣口服液通過抗炎發(fā)揮抗病毒功效。本次抗疫情過程中，提到“炎癥風暴”是新冠狀病毒患者從輕癥向重癥轉化的重要節(jié)點。炎癥風暴引起機體免疫損傷，引起體液中多種細胞因子產生，是引起多器官衰竭的重要原因[16]。由網絡藥理學結果顯示，金欣口服液具有抗炎抗病毒免疫的作用，提示金欣口服液在間接抗病毒方面，用于新型冠狀病毒肺炎患者，恢復機體整體功能具有一定的可行性。

網絡藥理學通過挖掘數據庫中的信息，從整體、多靶角度闡述中藥復方可能作用的環(huán)節(jié)，為中醫(yī)藥科學研究提供指導。本研究基于網絡藥理學挖掘出金欣口服液中可能發(fā)揮用于治療COVID-19的基因靶點，為后期驗證工作提供相應的理論依據。需要指出的是，網絡藥理學的研究多基于數據庫進行數據挖掘，而數據的信息往往來源于文獻，因此一些熱點化合物和靶點基因在分析過程會被富集，而一些強效的特異性化合物信息可能會丟失。在后期的研究中，應基于復方特色和化合物結構特點，選擇性開展驗證工作。