基于網(wǎng)絡藥理學和分子對接探討苦參堿在奶牛生產(chǎn)中的抗炎作用機制

趙梓軒 楊嘉睿 蔣林樹 童津津

(北京農(nóng)學院動物科學技術學院，奶牛營養(yǎng)學北京市重點實驗室，北京102206)

隨著居民生活水平的不斷提高，畜產(chǎn)品的安全質(zhì)量問題日益受到重視，這對我國畜牧業(yè)實現(xiàn)綠色、健康和高質(zhì)量發(fā)展提出了更高要求。在過去的很長一段時間里，獸用抗生素對動物疾病的防治、養(yǎng)殖效益的提高以及畜禽產(chǎn)品的高效供應起到了舉足輕重的作用；然而，抗生素使用不當?shù)膯栴}也日益凸顯，對畜產(chǎn)品的質(zhì)量、公共衛(wèi)生、生態(tài)環(huán)境以及動物健康造成了危害[1]，阻礙了我國畜牧業(yè)大力推進綠色健康養(yǎng)殖發(fā)展的步伐。所以，尋找安全有效并且不產(chǎn)生耐藥性的抗生素替代物成為當務之急。

苦參(SophoraflavescensAit.)在我國已經(jīng)具有兩千多年的使用歷史，作為傳統(tǒng)的中藥材其作用主要有清熱利尿和祛濕殺蟲等。據(jù)報道，苦參的主要生物活性成分為苦參堿，其同時也是發(fā)揮抗菌、抗病毒、抗腫瘤以及抗炎等功能的最主要的活性物質(zhì)[2]。研究發(fā)現(xiàn)，不僅是苦參，還包括苦豆、廣豆根等其他豆科植物的中藥材，它們的有效活性成分均為苦參堿，都能發(fā)揮殺菌消炎、抵抗病毒和腫瘤、殺蟲以及調(diào)節(jié)機體免疫力等多種生理功能[3]。此外，研究表明，苦參堿對斷奶仔豬腹瀉[4]、牛乳頭狀瘤病毒[5]以及小鼠膠原誘導性關節(jié)炎[6]等疾病的治療均有較好的效果?？鄥A在奶牛臨床應用中也取得了較好的療效。據(jù)報道，苦參堿對預防和治療奶牛子宮內(nèi)膜炎的效果顯著[3]，對奶牛乳腺上皮細胞(bMECs)增殖、凋亡及抗氧化能力有一定的影響[7]。盧金霞[8]研究發(fā)現(xiàn)，苦參堿對引起奶牛乳房炎的常見致病菌金黃色葡萄球菌具有顯著的抑制作用，可通過干預其對bMECs的黏附作用，從而發(fā)揮苦參堿對乳腺的保護作用。

雖然中藥在疾病預防和治療方面發(fā)揮著舉足輕重的作用，但是由于人們在其發(fā)揮療效的作用靶點及機制方面尚未闡明，嚴重限制了中藥的重新定位及作為新獸藥的開發(fā)利用[9]。網(wǎng)絡藥理學借助計算機技術、系統(tǒng)生物學和多向藥理學等多學科的交叉融合，形成了以網(wǎng)絡靶標為核心的技術體系。據(jù)報道，通過網(wǎng)絡藥理學方法對藥物的組成成分以及作用機制進行研究，可從生物網(wǎng)絡穩(wěn)態(tài)的角度來揭示藥物與疾病之間的關系[10]。網(wǎng)絡藥理學的應用可以提高疾病的治療效率，減少機體對相關藥物的不良反應，提升中藥治療疾病的成功率，節(jié)約藥品研發(fā)成本[11]。因此，本研究通過網(wǎng)絡藥理學的研究方法，探討苦參堿在奶牛生產(chǎn)中發(fā)揮抗炎作用的機制，以期為苦參堿在奶牛生產(chǎn)中的臨床應用提供參考。

1 材料與方法

1.1 苦參堿相關靶點的獲取

中藥系統(tǒng)藥理學數(shù)據(jù)庫與分析平臺(TCMSP，https://tcmsp-e.com/tcmsp.php)是一個查找藥物、靶點和疾病之間的關系的有效中藥系統(tǒng)藥理平臺[12]。本文以“matrine”為關鍵詞進行搜索，將獲取的苦參堿有效化合物作用靶點信息輸入UniProt數(shù)據(jù)庫(https：//www.uniprot.org)和PubChem(https://pubchem.ncbi.nlm.nih.gov)數(shù)據(jù)庫(即有機小分子生物活性數(shù)據(jù)庫，是一種化學模組的數(shù)據(jù)庫[13])。然后通過PharmMapper(http://www.lilab-ecust.cn/pharmmapper/submitfile.html)預測并得到苦參堿相關的靶點。

1.2 炎癥相關靶點的獲取

GeneCards是一個可檢索和集成的數(shù)據(jù)庫，可以在大量的基因與疾病、變異與蛋白質(zhì)和細胞與通路之間進行有效鏈接[14]。DisGeNET是一個綜合和規(guī)范化各種與疾病有關的基因和有關變異的數(shù)據(jù)庫，其中包含了科學文獻，同時涵蓋了所有的人類疾病以及正常與不正常的特點[15]。以“inflammation”為關鍵詞，分別在上述2個數(shù)據(jù)庫進行搜索，獲取的交集為炎癥相關靶點。

1.3 苦參堿抗炎靶點可視化分析

將上述獲得的成分靶點、炎癥相關靶點上傳于Venny(https://bioinfogp.cnb.csic.es/tools/venny/)網(wǎng)絡在線服務器中繪制韋恩圖，獲得苦參堿-炎癥靶點的交集，并進行可視化處理。

1.4 蛋白質(zhì)-蛋白質(zhì)相互作用(protein-protein interaction，PPI)網(wǎng)絡構建

登錄STRING數(shù)據(jù)庫上傳苦參堿抗炎靶點，分析苦參堿的相關靶點并構建PPI網(wǎng)絡。STRING數(shù)據(jù)庫是現(xiàn)階段信息量最豐富，研究蛋白相互作用的數(shù)據(jù)庫之一，該數(shù)據(jù)庫對已知和預測的蛋白質(zhì)互作信息進行評估與整合，形成了覆蓋大于1 100生物體的綜合蛋白質(zhì)網(wǎng)絡。該數(shù)據(jù)庫綜合多個數(shù)據(jù)庫的相關試驗數(shù)據(jù)，同時也含有生物信息學等研究方法預測的相關研究結果。STRING數(shù)據(jù)庫依據(jù)不同證據(jù)通道的概率計算，結合隨機交互概率校正，最終綜合得分為0.4時，為中等可信度；為0.7時，為高可信度；為0.9時，為最高可信度[16]。本研究選用了評分為0.4的PPI網(wǎng)絡，即中等可信度網(wǎng)絡。登錄STRING數(shù)據(jù)庫導入苦參堿發(fā)揮抗炎作用的潛在靶點并選擇“Bostaurus”物種，獲取靶點蛋白相互作用關系；required interaction score選擇“medium confidence(0.400)”并導出tsv文件；上傳tsv文件至Cytoscape 3.9.1并制作PPI網(wǎng)絡圖，選擇Cytoscape 3.9.1插件CytoNCA，對PPI網(wǎng)絡進行拓撲分析并得到關鍵核心靶點，通過CytoNCA插件計算節(jié)點的介數(shù)，接近中心性和子網(wǎng)絡中心性，疊加3個指標，篩選出前10%作為關鍵核心靶點[17]。

1.5 基因本體(gene ontology，GO)功能和京都基因與基因組百科全書(Kyoto encyclopedia of genes and genomes，KEGG)通路富集分析

登錄David數(shù)據(jù)庫(https://david.ncifcrf.gov/)輸入苦參堿抗炎靶點，選擇“Bostaurus”物種，select identifier設置為official gene symbol，list type設置為gene list，進行GO功能和KEGG通路富集分析，導出txt文件。將靶點數(shù)目進行降序排列，選擇生物學過程(biological process，BP)、細胞組分(cellular component，CC)和分子功能(molecular function，MF)進行富集分析。根據(jù)微生信在線作圖平臺(http://www.bioinformatics.com.cn/)對相關數(shù)據(jù)進行分析并制作直方圖和KEGG氣泡圖。

1.6 苦參堿-炎癥靶點-信號通路網(wǎng)絡構建

Cytoscape源于系統(tǒng)生物學，其將生物分子交互網(wǎng)絡、高通量基因表達數(shù)據(jù)和其他的分子狀態(tài)信息綜合分析整合，最終將數(shù)據(jù)結合成可視化網(wǎng)絡[18]。本研究通過Cytoscape 3.9.1軟件(https://cytoscape.org/)導入苦參堿抗炎靶點以及相關通路信息，建立苦參堿-炎癥靶點-信號通路網(wǎng)絡圖。

1.7 分子對接

利用TCMSP數(shù)據(jù)庫(https://tcmsp-e.com/tcmsp.php)搜索苦參堿，復制苦參堿的InChIKey，登錄PubChem數(shù)據(jù)庫(https://pubchem.ncbi.nlm.nih.gov/)，搜索苦參堿InChIKey，下載配體苦參堿的2D結構，保存成SDF格式；利用Chem 3D軟件轉(zhuǎn)換為mol2格式文件，導入Autodock Tools 1.5.6軟件進行加全氫、設置為配體(自動分布電荷)、檢測配體的Torsion tree中的扭轉(zhuǎn)鍵、選擇扭轉(zhuǎn)鍵等處理，保存文件為pdbqt。登錄UniProt數(shù)據(jù)庫(https://www.uniprot.org)復制核心靶點的，選擇條件應滿足“已驗證，物種人”。再登錄RCSB PDB數(shù)據(jù)庫(https://www.rcsb.org/)，粘貼對應靶點Entry Name，選擇蛋白質(zhì)晶體分辨率介于2.0～3.0，下載核心靶點，即受體白蛋白(ALB)、腫瘤壞死因子(TNF)、MYC原癌基因堿性螺旋-環(huán)-螺旋轉(zhuǎn)錄因子(MYC)、半胱氨酸-天冬氨酸蛋白酶3(CASP3)、周期蛋白D1(CCND1)、CD44和白細胞介素-6(IL-6)的PDB格式文件。同時在AutodockTools 1.5.6軟件中對核心靶點蛋白采用除去水、加入氫等處理方法。首先在AutodockTools 1.5.6軟件中設置Grid，選擇對應的配體文件和受體文件，選擇半柔性對接，設置GridBox，包裹住全部受體，保存為GPF格式文件，運行RUN AutoGrid，設置對接參數(shù)及運算方法，保存為DPF格式文件；然后設置Dock，依次選擇配體文件、受體文件等，運行Run Autodock，保存為DPF格式文件；最終對運算結果進行分析并選用自由結合能最低的結合構象作為最佳構象。利用PyMOL 2.5.2軟件進行可視化處理。

1.8 細胞分組

采用DMEM-F12、10%胎牛血清(FBS)和1%青-鏈霉素培養(yǎng)基，在溫度為37 ℃的5% CO2細胞培養(yǎng)箱中培養(yǎng)bMECs。未經(jīng)處理正常培養(yǎng)的bMECs為對照組(CON組)；將無乳鏈球菌與bMECs按感染復數(shù)(MOI)50∶1、共孵育6 h為奶牛乳腺炎癥模型組(TRT組)[19-20]；參照已發(fā)表的文獻[21]，利用50 μg/mL的苦參堿添加無乳鏈球菌與bMECs炎癥模型共孵育6 h，為苦參堿組(MAT組)。每組3個重復。

1.9 Western blot

參照本團隊前期報道Western blot的研究方法[22-23]，檢測CASP3和CCND1蛋白表達量的變化?？贵wCASP3(#9662)和CCND1(#55506)均購自于Cell Signaling公司。使用Image J 軟件進行灰度值分析。

1.10 數(shù)據(jù)統(tǒng)計分析

試驗數(shù)據(jù)最初使用Excel 2019進行整理，應用軟件SPSS 20.0進行單因素方差分析(one-way ANOVA)，試驗結果用“平均值±標準差”表示，P<0.05表示差異顯著，P<0.01表示差異極顯著。

2 結果與分析

2.1 苦參堿炎癥靶點獲取結果

本研究發(fā)現(xiàn)，苦參堿在TCMSP數(shù)據(jù)庫中共有10個靶點，通過PharmMapper數(shù)據(jù)庫獲得276個靶點，兩者取并集共獲得282個苦參堿靶點。通過GeneCards數(shù)據(jù)庫獲得10 898個炎癥靶點，通過DisGeNET數(shù)據(jù)庫獲得467個靶點，兩者取并集共獲得10 900個炎癥靶點。

2.2 苦參堿炎癥靶點可視化分析

通過生物信息在線作圖平臺，分別輸入苦參堿靶點和炎癥靶點通過韋恩圖對苦參堿和炎癥靶點進行可視化(圖1)。

2.3 苦參堿炎癥靶點網(wǎng)絡構建及關鍵靶點的篩選

將獲得的苦參堿炎癥靶點共71個輸入STRING數(shù)據(jù)庫，限制物種為“Bostaurus”，獲取蛋白互作PPI網(wǎng)絡(圖2)。結果表明，節(jié)點數(shù)為71，邊數(shù)為180，預期的邊數(shù)為91，平均節(jié)點度值為5.07，平均局部聚類系數(shù)為0.47，PPI富集P值為2.22×10-16，保存為tsv文件。將文件導入Cytoscape 3.9.1軟件，根據(jù)其度值(degree)選擇排名前30繪制柱狀圖(圖3)，由圖可見degree排名前5的靶點主要包括ALB、TNF、MYC、CASP3及CCND1。CytoNCA插件計算節(jié)點中心性，篩選核心靶點。CytoNCA包含8種中心性，分別為介數(shù)、接近中心性、自由度、特征向量中心性、平均最短路徑長度、網(wǎng)絡中心性、子網(wǎng)絡中心性和信息中心性。本研究選取介數(shù)、接近中心性和子網(wǎng)絡中心性的前10%的節(jié)點，并將3個中心性重疊，從而獲得具有關鍵靶點的子網(wǎng)絡。由表1和圖4可知，苦參堿抗炎核心靶點為ALB、TNF、MYC、CASP3、CCND1、CD44和IL-6。

圖1 苦參堿炎癥靶點韋恩圖

ACAT2：乙酰輔酶A乙酰轉(zhuǎn)移酶2 acetyl-CoA acetyltransferase 2；ACO1：順烏頭酸酶1 aconitase 1；IDH：異檸檬酸脫氫酶 isocitrate dehydrogenase；HSD17B4：17β-羥基類固醇4型脫氫酶 hydroxysteroid 17-beta dehydrogenase 4；FABP2：脂肪酸結合蛋白2 fatty acid-binding protein 2；MTAP：S-甲基-5′-硫代腺苷磷酸化酶 S-methyl-5′-thioadenosine phosphorylase；GAD：谷氨酸脫羧酶 glutamate decarboxylase；ALDH18A1：乙醛脫氫酶18家族成員A1 aldehyde dehydrogenase 18 family member A1；ODC1：鳥氨酸脫羧酶1 ornithine decarboxylase 1；HPD：4-羥基苯丙酮酸雙加氧酶 4-hydroxyphenylpyruvate dioxygenase；PYGL：肝糖原磷酸化酶 glycogen phosphorylase, liver form；DGKA：二?；视图っ甫?diacylglycerol kinase alpha；CANT1：鈣激活核苷酸酶1 calcium activated nucleotidase 1；PRM1：魚精蛋白1 protamine 1；SYNCRIP：突觸結合蛋白結合細胞質(zhì)RNA相互作用蛋白 synaptotagmin binding cytoplasmic RNA interacting protein；SPR54：信號識別顆粒54 signal recognition particle 54；PDK3：丙酮酸脫氫酶激酶同工酶3 pyruvate dehydrogenase kinase 3；HPSE：乙酰肝素酶 heparanase；PSAP：鞘脂激活蛋白原 prosaposin；BMP7：骨形成蛋白7 bone morphogenetic protein 7；CCNA2 周期蛋白A2 cyclin A2；CAMK2A：鈣/鈣調(diào)蛋白依賴性蛋白激酶激酶2α calcium/calmodulin dependent protein kinase 2 alpha；BCR：斷裂點簇集群區(qū) breakpoint cluster region；ALK：間變性淋巴瘤激酶 anaplastic lymphoma kinase；ALB：白蛋白 albumin；MYC：MYC原癌基因堿性螺旋-環(huán)-螺旋轉(zhuǎn)錄因子 MYC proto-oncogene, basic helix-loop-helix transcription factor；CUL1：滯蛋白1 Cullin 1；SERPINA5：絲氨酸蛋白酶抑制劑家族A成員5 Serpin family A member 5；TPR：核蛋白易位啟動子區(qū)域 nucleoprotein translocation promoter region；UTRN：肌營養(yǎng)相關蛋白 Utrophin；TAB1：轉(zhuǎn)化生長因子β活化激酶結合蛋白1 transforming growth factor-β-activated kinase-binding protein 1；MMP2：基質(zhì)金屬蛋白酶2 matrix metalloproteinase 2；TTR：轉(zhuǎn)甲狀腺素蛋白 transthyretin；RAC1：Ras相關的C3肉毒桿菌毒素底物1 Ras-related C3 botulinum toxin substrate 1；MAPK6：絲裂原活化蛋白激酶6 mitogen-activated protein kinase 6；ESR1 雌激素受體1 estrogen receptor 1；CCND1：周期蛋白D1 cyclin D1；HBB：血紅蛋白β亞基 hemoglobin subunit beta；C8A：補體C8α鏈 complement C8 alpha chain；MAPKAPK2：絲裂原激活蛋白激酶激活蛋白激酶2 mitogen-activated protein kinase activated protein kinase 2；PRKCQ：蛋白激酶Cθ protein kinase C theta；CD3E：T細胞受體復合物CD3ε亞基 CD3 epsilon subunit of T-cell receptor complex；RELA：轉(zhuǎn)錄因子p65 transcription factor p65；TNF：腫瘤壞死因子 tumor necrosis factor；CASP3：半胱氨酸-天冬氨酸蛋白酶3 Caspase-3；TJP1：緊密連接蛋白1 tight junction protein 1；DPP4：二肽基肽酶-4 dipeptidyl peptidase-4；TLL1：類Tolloid蛋白1 Tolloid like 1；INSR：胰島素受體 insulin receptor；MUC1：黏蛋白1 mucin 1；ICAM1：細胞間黏附分子1 intercellular adhesion molecule 1；IL-6：白細胞介素-6 interleukin-6；ACBD7：乙酰輔酶A結合域包含7 acyl-CoA binding domain containing 7；BCL9：B細胞淋巴瘤9 B-cell lymphoma 9；PRKD2：絲氨酸/蘇氨酸蛋白激酶D2 Serine/threonine protein kinase D2；HMGB2：高遷移率族蛋白B2 high mobility group box 2；GTF2E1：通用轉(zhuǎn)錄因子2E亞基1 general transcription factor 2E subunit 1；HDAC8：組蛋白脫乙?；? histone deacetylase 8；ANXA1：膜聯(lián)蛋白A1 annexin A1；NMT1：N-十四烷酰轉(zhuǎn)移酶1 N-myristoyltransferase 1；PTGS1：前列腺素內(nèi)過氧化物合酶1 prostaglandin-endoperoxide synthase 1；PTPRJ：蛋白酪氨酸磷酸酶受體J protein tyrosine phosphatase receptor J；VPS4B：液泡蛋白分選蛋白4B vacuolar protein sorting 4 homolog B；KLRD1：殺傷細胞凝集素樣受體D1 killer cell lectin like receptor D1；CSF2RB：集落刺激因子2受體β亞基 colony stimulating factor 2 receptor subunit beta；ALOX12：花生四烯酸12-脂氧合酶 arachidonate 12-lipoxygenase；SMAP1：基質(zhì)膜相關蛋白1 stromal membrane-associated protein 1。圖3、圖4和圖7同 the same as Fig.3, Fig.4 and Fig.7。

圖3 苦參堿炎癥潛在靶點度值排序

表1 苦參堿炎癥核心靶點

通過拓撲分析，確定苦參堿對奶牛炎癥作用相關靶點的PPI網(wǎng)絡，顏色越深，其度值越高。

2.4 GO功能和KEGG通路富集分析

登錄David數(shù)據(jù)庫上傳71個苦參堿抗炎靶點，該物種限定為“Bostaurus”，并進行GO功能和KEGG通路富集分析。GO功能富集分析結果獲得62種BP、16種CC和12種MF。選取閾值P<0.01的條目，在微生信(http://www.bioinformatics.com.cn/)中導入GO term、Count、P-value以及分類相關數(shù)據(jù)生成BP、CC和MF三合一圖(圖5)。其中，BP主要包括細胞對過氧化氫的反應(cellular response to hydrogen peroxide)、一氧化氮(NO)生物合成過程的正調(diào)控(positive regulation of nitric oxide biosynthetic process)、脂多糖細胞反應(cellular response to lipopolysaccharide)、異型細胞黏附的正調(diào)控(positive regulation of heterotypic cell-cell adhesion)以及白細胞介素-8產(chǎn)生的正調(diào)控(positive regulation of interleukin-8 production)等；CC主要包括細胞質(zhì)(cytoplasm)和細胞外間隙(extracellular space)等；MF主要包括蛋白激酶結合(protein kinase binding)、ATP結合(ATP binding)和質(zhì)膜(membrane)等。

圖5 GO功能富集分析

KEGG通路富集分析(圖6)共獲得22條信號通路，主要包括腫瘤壞死因子信號通路(TNF signaling pathway)、絲裂原活化蛋白激酶信號通路(MAPK signaling pathway)、白細胞介素-17信號通路(IL-17 signaling pathway)、核因子-κB信號通路(NF-κB signaling pathway)以及Toll樣受體信號通路(Toll-like receptor signaling pathway)等通路，與苦參堿發(fā)揮抗炎作用密切相關。

圖6 KEGG通路富集分析

2.5 苦參堿-炎癥靶點-信號通路網(wǎng)絡構建

將David通路注釋得到的數(shù)據(jù)，登錄Cytoscape 3.9.1軟件上傳相關數(shù)據(jù)，構建苦參堿-炎癥靶點-信號通路網(wǎng)絡圖(圖7)。Network Analyzer結果顯示，網(wǎng)絡節(jié)點數(shù)為94，網(wǎng)絡邊數(shù)為182，網(wǎng)絡密度為0.04，網(wǎng)絡直徑為4，網(wǎng)絡半徑為2，網(wǎng)絡中心度為0.74，網(wǎng)絡異質(zhì)性為1.92，特征路徑長度為2.35。

2.6 分子對接結果

將苦參堿相關活性成分與潛在作用靶點拓撲分析網(wǎng)絡圖中連接度最高的成分苦參堿與7個核心靶點進行對接，即受體為ALB、TNF、MYC、CASP3、CCND1、CD44和IL-6，配體為苦參堿，分子對接結果見表2和圖8。結合自由能用來評估成分和靶點之間的結合；如果結合能小于0，則表明活性成分和靶點可以自由結合。研究表明，若兩者之間結合自由能較低，則兩者之間親和力較強[21]。通常認為，結合自由能<-5 kcal/mol表示活性成分與作用靶點之間親和力良好[20]。本研究結果表明，ALB、TNF、MYC、CASP3、CCND1、CD44和IL-6皆可與苦參堿穩(wěn)定結合，提示苦參堿與上述靶點之間的相互作用可能為苦參堿發(fā)揮抗炎的重要作用機制之一，需要進一步試驗驗證。

2.7 苦參堿對無乳鏈球菌誘導的bMECs相關蛋白表達量的影響

如圖9所示，與CON組相比，TRT組CASP3蛋白表達量極顯著升高(P<0.01)，CCND1蛋白表達量極顯著降低(P<0.01)；與TRT組相比，MAT組CASP3蛋白表達量極顯著降低(P<0.01)，CCND1蛋白表達量極顯著提高(P<0.01)。

Matrine：苦參堿；Metabolic pathways：代謝通路；MAPK signaling pathway：絲裂原活化蛋白激酶信號通路；IL-17 signaling pathway：白細胞介素-17信號通路；HIF-1 signaling pathway：缺氧誘導因子-1信號通路；Hematopoietic cell lineage：造血細胞譜系；Glutathione metabolism：谷胱甘肽代謝；Citrate cycle (TCA cycle)：檸檬酸循環(huán)(三羧酸循環(huán))；Cellular senescence：細胞衰老；Carbon metabolism：碳代謝；C-type lectin receptor signaling pathway C型凝集素受體信號通路；Butanoate metabolism：丁酸鹽代謝；Biosynthesis of amino acids：氨基酸生物合成；Adherens junction：黏著連接；2-oxocarboxylic acid metabolism；2-氧羧酸代謝；Wnt signaling pathway Wnt信號通路；Toll-like receptor signaling pathway Toll樣受體信號通路；TNF signaling pathway：腫瘤壞死因子信號通路；Th17 cell differentiation：Th17細胞分化；TGF-beta signaling pathway：轉(zhuǎn)化生長因子-β信號通路；T cell receptor signaling pathway：T細胞受體信號通路；NF-kappa B signaling pathway：核因子-κB信號通路；Neurotrophin signaling pathway：神經(jīng)營養(yǎng)因子信號通路。

表2 苦參堿與核心靶點的分子對接結果

3 討 論

3.1 苦參堿抗炎關鍵靶點分析

本研究發(fā)現(xiàn)，苦參堿在奶牛生產(chǎn)中發(fā)揮抗炎的核心靶點為ALB、TNF、MYC、CASP3、CCND1、CD44和IL-6。據(jù)報道，由巨噬細胞和單核細胞等特定細胞產(chǎn)生的促炎因子——腫瘤壞死因子-α(TNF-α)，用于刺激炎癥發(fā)生，并調(diào)節(jié)細胞中引起壞死和凋亡的一系列信號調(diào)控。這種炎性因子在抵抗感染、癌癥等過程中發(fā)揮一定的作用[24]。有研究表明，苦參堿對金黃色葡萄球菌脂磷壁酸(LTA)誘導的小鼠子宮內(nèi)膜炎具有一定的保護作用，且呈劑量依賴性降低促炎細胞因子[TNF-α和白細胞介素-1(IL-1)]的表達。此外，研究發(fā)現(xiàn)，在LTA刺激的牛子宮內(nèi)膜上皮細胞(bEECs)的炎癥模型中，苦參堿處理可抑制Toll樣受體(TLR)2的表達及其下游核因子NF-κB的活化，從而降低炎癥損傷[25]。于曉紅等[26]建立了脂多糖(LPS)誘導的牛子宮內(nèi)膜上皮細胞炎癥模型，結果表明，苦參堿能顯著降低TNF-α、IL-1、IL-6和白細胞介素-8(IL-8)mRNA的表達，具有較強的抗炎作用。據(jù)報道，IL-6是一種多功能細胞因子，通過其獨特的受體系統(tǒng)具有多種生物活性，能夠調(diào)節(jié)免疫應答、急性期應答和炎癥反應[27]。此外，與本研究結果一致，丁園園等[28]研究發(fā)現(xiàn)，苦參堿通過調(diào)節(jié)病毒復制、凋亡和炎癥反應作用于TNF-α、IL-6和CASP3靶點，從而達到治療COVID-19的目的。本研究通過Western blot結果證明，苦參堿可通過對CASP3和CCND1蛋白的調(diào)節(jié)發(fā)揮抗炎作用。據(jù)報道，清肝化瘀顆粒的有效成分苦參堿以劑量依賴方式抑制了肝癌細胞的增殖[29]，對TP53和CCND1等核心靶點具有調(diào)節(jié)作用，對肝癌細胞氧化應激脅迫標志物產(chǎn)生影響，引起細胞凋亡。綜上可知，本研究進一步闡明了苦參堿抗炎作用的核心靶點及其作用機制。

a～g分別代表苦參堿與ALB、TNF、MYC、CASP3、CCND1、CD44和IL-6的對接示意圖。

**表示P<0.01，****表示P<0.000 1。

3.2 苦參堿炎癥GO功能和KEGG關鍵通路分析

GO功能富集分析結果表明，苦參堿發(fā)揮抗炎作用的主要靶點是通過過氧化氫反應、NO生物合成過程和IL-8等途徑。這與先前的研究結果一致，苦參堿可作用于白細胞介素-1β(IL-1β)、IL-17淋巴細胞或單核巨噬細胞釋放的炎性細胞因子，通過Nrf2信號通路發(fā)揮作用，與提高機體抗氧化能力及降低機體氧化損傷密切有關[30]。NO作為機體重要的信號傳遞因子，對中樞神經(jīng)系統(tǒng)、免疫系統(tǒng)及生殖系統(tǒng)等都具有重要作用，可廣泛參與哺乳動物的生理調(diào)節(jié)。張想軍等[31]發(fā)現(xiàn)，25～100 μg/mL的苦參堿可顯著促進bMECs的增殖，并可抑制金黃色葡萄球菌入侵感染bMECs，此外苦參堿可通過干預抑制NO的分泌和NF-κB的活化發(fā)揮其抗炎作用。而且，研究表明細胞因子干擾素-γ(INF-γ)和IL-1可誘導巨噬細胞表達大量的誘導型一氧化氮合酶(iNOS)蛋白并產(chǎn)生NO，NO作為一個免疫功能分子從而發(fā)揮抗炎抗病毒作用[32]。IL-8作為一種細胞因子，在炎癥反應發(fā)生時，具有重要的調(diào)節(jié)作用[33]，研究證明，苦參堿可顯著降低血清中IL-8含量，有效抑制中晚期舌癌術后患者炎癥反應和腫瘤轉(zhuǎn)移[21]。由此可見，苦參堿對炎癥反應具有顯著的抑制作用。

KEGG通路富集分析顯示，苦參堿在奶牛生產(chǎn)中發(fā)揮抗炎作用的機制與多種信號通路等密切相關。據(jù)報道，MAPK信號通路作為哺乳動物細胞中重要的信號傳遞方式之一，可將細胞表面的信號通過刺激傳導給細胞及其細胞核，參與細胞的增殖、分化以及細胞凋亡等生物學過程[34]。研究發(fā)現(xiàn)，LPS與細胞膜表面的TLR結合激活轉(zhuǎn)錄因子NF-κB，活化的NF-κB轉(zhuǎn)移到細胞核內(nèi)，從而啟動炎癥因子如TNF-α、白細胞介素-19(IL-19)和IL-6等基因的轉(zhuǎn)錄[35]。IL-17信號依賴的基因有一種獨特的、組織特異性的模式，是奶牛體內(nèi)重要的生理功能的基礎[36]。有研究表明，苦參堿通過抑制絲裂原活化蛋白激酶(MKK)/p38 MAPK信號通路抑制炎癥反應，對氧化低密度脂蛋白(ox-LDL)暴露巨噬細胞具有良好的抗炎作用[37]。據(jù)報道，苦參堿可通過抑制TNF-α、IL-6和IL-1β的生成，減少炎癥浸潤和脫髓鞘，緩解自身免疫性腦脊髓炎(EAE)的發(fā)生與發(fā)展。此外，苦參堿還通過高遷移率族蛋白B1(HMGB1)/TLR4/NF-κB信號通路激活發(fā)揮一定的抗炎作用[38]。此外，有研究發(fā)現(xiàn)，苦參堿可通過調(diào)節(jié)Th17所誘導的角質(zhì)形成細胞，抑制白細胞介素-17A(IL-17A)、白細胞介素-21受體(IL-21R)以及白細胞介素-22受體1(IL-22R1)的表達，對由Th17介導的免疫性皮膚病中發(fā)揮較好的抗炎作用[39]。然而，目前關于苦參堿對奶牛抗炎作用相關的通路研究報道較少，本研究結果為進一步闡明苦參堿在奶牛抗炎作用中的機制研究提供了依據(jù)。

4 結 論

本研究通過網(wǎng)絡藥理學、分子對接技術及多種生物信息學工具，對苦參堿多靶點、多途徑防治奶牛炎癥的互作網(wǎng)絡進行分析?？鄥A發(fā)揮抗炎作用的核心靶點為ALB、TNF、MYC、CASP3、CCND1、CD44和IL-6。KEGG通路富集分析結果表明，苦參堿可通過MAPK信號通路、TNF信號通路、IL-17信號通路、Toll樣受體信號通路和NF-κB信號通路等發(fā)揮抗炎作用。本研究結果為苦參堿進一步在治療奶牛炎癥相關疾病中的應用提供了理論依據(jù)。