項丹丹，楊曉芳，易干軍，陶海青，初元琦，李春雨

（廣東省農業(yè)科學院果樹研究所，農業(yè)農村部南亞熱帶果樹生物學與遺傳資源利用重點實驗室/廣東省熱帶亞熱帶果樹研究重點實驗室，廣州 510640）

香蕉作為世界貿易量最大的水果，也是世界十大主食之一。香蕉產業(yè)是我國熱帶亞熱帶地區(qū)農業(yè)支柱產業(yè)之一，該產業(yè)為我國國民經濟的增長作出了重要貢獻。然而，香蕉枯萎?。‵usarium wilt of banana）在我國廣東、廣西等香蕉主產區(qū)的蔓延，嚴重制約著我國香蕉產業(yè)的健康持續(xù)發(fā)展。香蕉枯萎病是一種由尖孢鐮刀菌古巴?；停‵usarium oxysporumf.sp.cubense,Foc）侵染引起維管束壞死的毀滅性土傳真菌病害。Foc遺傳多樣性豐富，根據致病香蕉品種的類型，可分為4個生理小種，其中危害性最大的是香蕉枯萎病菌熱帶4 號生理小種（Foctropical race 4,FocTR4），能侵染幾乎所有香蕉栽培品種[1]。近年來，由于FocTR4引起的枯萎病發(fā)展態(tài)勢嚴峻，我國香蕉產業(yè)收獲面積逐年下降[2]。2019 年，FocTR4 在哥倫比亞以及加勒比海地區(qū)的首次檢出預示著這種病害已成為全球香蕉產業(yè)的最大威脅[3]。

香蕉枯萎病的防治是世界性難題，目前主要采用培育抗性品種、生物防治、化學防治等手段進行綜合防控[3]?？梗停┎∑贩N的推廣正在香蕉枯萎病的防治中發(fā)揮關鍵作用，但近年來病原菌生理小種的多樣化加速了品種抗病性的喪失。生物防治發(fā)揮作用需要的時間較長、過程較為復雜，在病害防控中仍處于研發(fā)階段。化學殺菌劑具有高效、速效、使用方便、經濟效益高等優(yōu)點，在病蟲害綜合防治中占有重要地位。國內外學者開展了多種殺菌劑對Foc作用的試驗，如NEL等[4]表明1 μg/mL咪鮮胺可顯著抑制Foc菌絲生長；許文耀等[5]發(fā)現霉靈·溴菌腈對Foc孢子萌發(fā)有良好的抑制效果；除此之外，許多報道表明咪鮮胺及其混配劑對Foc有良好的抑菌活性[6-8]。但在田間環(huán)境下，室內試驗篩選的藥劑往往不能表現出相同的抑菌效果，難以在生產中得到應用與推廣[9]。因此，到目前為止還沒找到一種防治香蕉枯萎病的理想殺菌劑[10]。原因也在于香蕉枯萎病是一種維管束類型的病害，病菌以不同形式存在于土壤中，農藥施用于土壤后難以接觸到病菌。深入開展防治香蕉枯萎病化學藥劑的研究，篩選開發(fā)有效的藥劑，是香蕉枯萎病綜合防治體系中必須攻克的瓶頸。

三唑類殺菌劑具有優(yōu)良的保護、治療作用，其內吸性強，作用的病原菌產生抗藥性的風險低，持效期長，對環(huán)境與使用者安全，因此在殺菌劑領域占有重要地位。目前三唑類殺菌劑品種已有四十幾種，占市場上殺菌劑的1/3[11]。且丙環(huán)唑、戊唑醇、腈菌唑、丙硫菌唑等三唑類殺菌劑被陸續(xù)研發(fā)出來。三唑類殺菌劑屬于甾醇脫甲基抑制劑（sterol demethylation inhibitor,DMI）類殺菌劑，主要通過抑制甾醇生物合成中細胞色素P450 甾醇14α-脫甲基酶（cytochrome P450 sterol 14α-demethylase,CYP51）的活性來破壞病原菌細胞膜功能，最終導致細胞死亡[12]。該類殺菌劑屬于N-1-取代-1，2，4-三唑環(huán)的有機氮雜環(huán)類化合物，主要通過發(fā)揮內吸功能，隨蒸騰作用運輸至植物各部位。由于其獨特的作用方式及優(yōu)異的生物活性，三唑類殺菌劑已被廣泛應用于小麥、玉米、水稻、高粱等作物以及水果和蔬菜等的多種真菌病害防治中。但目前尚未見系統比較不同三唑類殺菌劑防治香蕉枯萎病效果的報道。因此，本研究基于常用的三唑類殺菌劑，開展殺菌劑對FocTR4 的抑菌效果研究，通過對三唑類殺菌劑與其靶標CYP51 的相互作用來探究它們對FocTR4 的抑制效果。本研究結果不但可為我國香蕉枯萎病的化學防治提供新策略，而且也為殺菌劑的科學合理應用提供了指導。

1 材料與方法

1.1 供試材料

供試菌株FocTR4（菌株Ⅱ5，VCG01213），保存于廣東省農業(yè)科學院果樹研究所香蕉遺傳改良中心，使用前將菌株于馬鈴薯葡萄糖瓊脂（potato dextrose agar,PDA）平板上進行活化培養(yǎng)（28 ℃，3 d）。丙硫菌唑、腈菌唑、丙環(huán)唑、戊唑醇標準品均購自德國默克公司。分別稱取對應量的不同三唑類殺菌劑原藥溶于二甲基亞砜（dimethyl sulfoxide,DMSO）中配制成10 mg/mL的母液，置于4 ℃冰箱中保存，備用。使用時將母液用DMSO稀釋到所需濃度，母液和稀釋后的藥劑每隔2 周重新配制。RNA提取試劑盒、cDNA 第一鏈合成試劑盒（HiScript Ⅱ1st Strand cDNA Synthesis Kit）和 以RNA 為 模 板的染料法實時熒光定量聚合酶鏈反應（real-time quantitative polymerase chain reaction, qRT-PCR）試 劑 盒（HiScript?ⅡOne Step qRT-PCR SYBR Green Kit），均購自南京諾唯贊生物科技股份有限公司；二喹啉甲酸蛋白濃度測定試劑盒（Bicinchoninic Acid Protein Assay Kit），購自廣州市紅爍林生物科技有限公司。

1.2 抑菌活性的測定

采用菌絲生長速率法測定不同三唑類殺菌劑對FocTR4 的抑菌活性。將PDA 加熱至完全融化，待其冷卻至40～50 ℃時加入相應濃度的丙硫菌唑、丙環(huán)唑、戊唑醇和腈菌唑，使得4 種三唑類殺菌劑在PDA 培養(yǎng)基中的終質量濃度分別為0.312 5、0.625、1.25、2.5、5、10、15、20 μg/mL，用已滅菌的打孔器在菌落邊緣打取直徑為5 mm的菌餅，菌餅倒置接種在冷卻至室溫的PDA平板上。用封口膜將平板密封，倒置放入28 ℃恒溫培養(yǎng)箱中培養(yǎng)。每個處理重復3次，同時以加入等體積的DMSO 為對照。培養(yǎng)7 d后利用十字交叉法測定各平板菌落的生長直徑并計算各濃度三唑類殺菌劑的抑制率，通過概率值分析法計算半最大效應濃度（concentration for 50%of maximal effect,EC50）值。菌絲生長抑制率計算公式如下：

1.3 菌絲形態(tài)觀察

菌株培養(yǎng)和處理同1.2節(jié)。于28 ℃的培養(yǎng)箱中倒置培養(yǎng)7 d 后，用牙簽挑取菌落邊緣少量菌絲于磷酸鹽緩沖液（phosphate buffer saline, PBS）（pH 7.4）中，用移液槍輕輕吹打菌絲至分散，并吸取少量菌液滴加于載玻片上，蓋上蓋玻片，在光學顯微鏡下觀察菌絲形態(tài)并拍照。

在平板邊緣打取直徑為8～10 mm的菌餅，將菌餅置于2.5%戊二醛-4.0%多聚甲醛固定液（pH 6.8）中，于4 ℃條件下進行前固定。3 h后用PBS（pH 6.8）沖洗樣品3 次，繼續(xù)用1%鋨酸固定并過夜，再次用PBS（pH 6.8）沖洗樣品3 次。用30%、50%、70%、90%、95%、100%、100%、100%乙醇梯度脫水，每個濃度梯度處理15～20 min。臨界點干燥、粘樣、噴金鍍膜后，于S-3400N-Ⅱ掃描電子顯微鏡（日本日立公司）下觀察、拍照。以上每個處理重復3次。

1.4 菌絲電導率的測定

菌株培養(yǎng)和處理同1.2 節(jié)。參考張俊[13]的方法測定菌絲細胞膜的通透性以確定菌絲電導率。用滅菌的打孔器在菌落邊緣打取菌餅并放入配制好的馬鈴薯葡萄糖液體（potato dextrose broth,PDB）培養(yǎng)基中，于28 ℃、180 r/min 振蕩培養(yǎng)箱中培養(yǎng)7 d后，用5 層擦鏡紙過濾。將稱取的0.15 g 菌絲放入已提前裝好等量雙蒸水的小燒杯中，加入不同質量濃度（0.25、1.25 μg/mL）的丙硫菌唑、丙環(huán)唑、戊唑醇、腈菌唑，以加入等體積的DMSO 為對照。分別在5、15、30、60、90、120、150、180、240、300、360 min時用電導儀測定電導率，最后將菌絲煮沸15 min，冷卻至室溫后測定最終電導率，每個處理設置3 個重復。菌絲相對電導率計算公式如下：

1.5 丙二醛（malondialdehyde,MDA）含量測定

在超凈工作臺中，將相應濃度的三唑類殺菌劑加入冷卻至40～50 ℃的PDA 培養(yǎng)基中，使不同殺菌劑的終質量濃度分別為0.25、1.25 μg/mL，同時分別以加入等體積的DMSO 為對照。待平板冷卻凝固，將已提前滅菌的玻璃紙平鋪于PDA 表面，待玻璃紙表面水分吹干后接種菌餅，隨后密封置于28 ℃恒溫培養(yǎng)箱中培養(yǎng)7 d 后收集菌絲。采用MDA 含量檢測試劑盒（微量法）（南京建成生物工程研究所）進行MDA 含量的測定：分別稱取對照組和處理組菌絲樣品0.10 g，加入1 mL提取液進行冰浴勻漿，于8 000 r/min 下離心10 min，收集上清液，根據試劑盒說明書進行操作，反應結束后，在450 nm 波長下測定各孔的吸光度值，計算樣品中MDA含量。

1.6 CYP51 酶活性的測定

菌株培養(yǎng)和處理同1.2節(jié)。根據細胞色素P450芳香化酶檢測試劑盒（上海江萊生物科技有限公司）說明書要求進行CYP51酶活性的測定。準確稱取0.10 g 菌絲樣品，加入1 mL 提取液，利用超聲波破碎菌絲細胞（冰浴，功率200 W，超聲3 s 停10 s，重復30 次），于4 ℃、12 000 r/min 下離心10 min，收集上清液，置于冰上待測。加入50 μL 終止液終止反應，在450 nm 波長下依次測定各孔的吸光度值。用標準品濃度與吸光度值計算出標準曲線的直線回歸方程。在標準曲線上計算出樣品經過稀釋后的蛋白質濃度，再乘以稀釋倍數，即為樣品的實際濃度。

1.7 基因表達量的測定

參照RNA 提取試劑盒說明書提取菌絲的總RNA，以RNase-free 水為對照，檢測RNA 的濃度與純度；采用1%瓊脂糖凝膠電泳檢測RNA的完整性；采用cDNA第一鏈合成試劑盒對RNA進行反轉錄，具體試劑添加和反應條件按照說明書要求進行，反應終止后得到的cDNA 于－80 ℃冰箱中保存，備用。

所用引物（表1）由生工生物工程（上海）股份有限公司合成，20 μL qRT-PCR體系：0.4 μL 10 μmol/L正、反向引物，10 μL 2×通用型高靈敏度染料法定量PCR 預混液（Cham Q Universal SYBR qPCR Master Mix），2 μL cDNA（提前用ddH2O 稀釋4 倍），7.2 μL ddH2O。反應擴增程序：95 ℃預變性30 s；95 ℃變性10 s，60 ℃延伸10 s，40 個循環(huán)。溶解曲線制作：95 ℃變性15 s，60 ℃復性60 s，升溫至95 ℃保持15 s，結束整個過程。目的基因的相對表達量采用2－ΔΔCT法進行計算。

表1 實時熒光定量PCR引物序列Table 1 Primer sequences for qRT-PCR

1.8 同源建模和分子對接

選擇與FocTR4 的CYP51（TR4-CYP51）蛋白質序列一致性大于50% 的蛋白質序列pdb ID 6CR2[14]、5V5Z[15]、5ESI[16]和5JLC[16]進 行 多 模 板 建模[17]。使用Modeller 9.20 軟件[18]多模板建模工具預測CYP51 的三維結構，獲得1 000 個結構，通過DOPE 軟件打分選擇最佳模型；用AutoDock Tools制備受體pdbqt 文件，加氫和分配電荷。采用AutoDock Vina工具進行分子對接，選擇血紅素鐵原子為對接中心，設置對接盒子為15 ?×15 ?×15 ?。根據模板晶體結構中配體的結合構象，主要按照如下2 個標準選擇最優(yōu)的結合構象：第一是配體結合在4 個芳香氨基酸殘基Tyr137、Phe131、Phe230 和Phe511 組成的結合口袋中；第二是三氮唑環(huán)上的N原子與血紅素鐵接近，N-FE 接近垂直于血紅素平面。使用PyMOL軟件和Ligplot軟件對結果進行分析和對接構象的可視化。

1.9 親和動力學分析

前期的研究中，我們已通過原核表達系統獲得了純度較高的CYP51 蛋白（純度≥90%）[17]。采用表面等離子共振（surface plasmon resonance,SPR）技術，通過Biacore T200軟件（瑞士GE公司）進行多循環(huán)動力學測定，分析不同三唑類殺菌劑與CYP51之間的相互作用。采用氨基偶聯反應將CYP51蛋白共價連接于CM5傳感芯片表面。將三唑類殺菌劑稀釋成濃度梯度125.000、62.500、31.250、15.625 nmol/L，緩沖液為含5%DMSO的HBS-EP［0.01 mol/LN-（2-羥乙基）哌嗪-N′-（2-乙磺酸）（HEPES），0.15 mol/L NaCl，0.003 mol/L 乙二胺四乙酸（EDTA），0.05%（體積比）表面活性劑P20，pH 7.4］。依次進行親合力與動力學測試，每個濃度設置1 個循環(huán)，結合120 s，解離300 s，設置溫度25 ℃，流速30 μL/min。Biacore T200 Evaluation 軟件用于數據分析，同時使用1∶1穩(wěn)態(tài)親和力模型進行動力學/親和力擬合分析。

1.10 數據處理與分析

采用GraphPad Prism 7.00和Statistics 19軟件對數據進行處理，對照組和不同三唑類殺菌劑處理組采用單因素方差分析方法，并選擇Tukey’s 多重比較法進行檢驗分析。

2 結果與分析

2.1 不同三唑類殺菌劑對Foc TR4生物活性的影響

采用菌絲生長速率法測定了4種三唑類殺菌劑對香蕉枯萎病菌FocTR4 的抑制作用。結果表明：對照組菌絲生長迅速，各藥劑處理組菌絲生長緩慢；對FocTR4的抑制效果存在濃度依賴效應，隨著三唑類殺菌劑濃度的升高，菌落直徑逐漸減小（圖1A）。在相同濃度下，不同殺菌劑對FocTR4 抑制效果存在顯著差異，其中丙硫菌唑的抑菌活性最高，EC50為0.79 μg/mL（R2=0.95）；戊唑醇和丙環(huán)唑對FocTR4 的抑制作用低于丙硫菌唑，EC50分別為1.33 μg/mL（R2=0.95）和5.32 μg/mL（R2=0.95）；腈菌唑的抑菌活性遠低于其他3 種殺菌劑，其EC50為10.13 μg/mL（R2=0.95）（圖1A～B）。

利用光學顯微鏡觀察殺菌劑處理對FocTR4菌絲形態(tài)的影響發(fā)現，對照組菌絲表面光滑，生長粗細均勻，分枝正常；而處理組菌絲呈現不規(guī)則分枝增多、生長粗細不均等畸形現象（圖1C）。進一步通過掃描電子顯微鏡觀察發(fā)現，三唑類殺菌劑處理可引起FocTR4 菌絲出現表面干癟、凹陷和扁平化現象，其中丙硫菌唑對FocTR4 菌絲的形態(tài)影響最為明顯。光學顯微鏡下觀察到腈菌唑引起FocTR4菌絲分枝增多，但掃描電子顯微鏡觀察到其對菌絲的顯微結構無明顯影響（圖1D）。以上結果表明，這4種殺菌劑均可影響FocTR4 菌絲的正常生長發(fā)育，且不同殺菌劑引起的菌絲畸形程度存在差異。

圖1 不同三唑類殺菌劑對Foc TR4的抑制效果及菌絲形態(tài)Fig.1 Inhibitory effects of different triazole fungicides against Foc TR4 and hyphae morphology

2.2 不同三唑類殺菌劑對Foc TR4 細胞膜通透性的影響

電導率可作為病原菌菌絲細胞膜損傷程度的評價指標，電導率越大，表明細胞內電解質滲漏越嚴重，細胞膜的損傷越厲害。因此，本研究通過測定不同時間點FocTR4菌絲在雙蒸水中的相對電導率來表征三唑類殺菌劑對菌絲細胞膜通透性的影響，相對電導率越大，表示菌絲細胞膜通透性越高。如圖2所示，隨著時間的延長，對照組和處理組菌絲的相對電導率均逐漸升高，且不同時間點各處理組菌絲電導率均顯著高于對照組，在處理時間點為360 min 時，對照組的電導率為20.69%，經0.25、1.25 μg/mL丙硫菌唑處理的菌絲相對電導率分別達到59.16%和61.42%（圖2A），相同條件下戊唑醇處理的菌絲相對電導率分別為32.03%和51.37%（圖2B），丙環(huán)唑處理組相對電導率分別為30.55%和47.65%（圖2C），腈菌唑處理組的相對電導率分別為30.70%和41.55%（圖2D）。以上結果表明，4種三唑類殺菌劑均可引起FocTR4菌絲細胞內的電解質外滲，提高相對電導率，且濃度越高，對菌絲細胞膜通透性的影響越大。除此之外，不同殺菌劑對菌絲細胞膜通透性的影響也存在差異，表現為丙硫菌唑＞戊唑醇＞丙環(huán)唑＞腈菌唑。

圖2 不同三唑類殺菌劑對Foc TR4菌絲細胞相對電導率的影響Fig.2 Effects of different triazole fungicides on the relative conductivity of Foc TR4 hypha cells

MDA 是細胞內部生物大分子過氧化的產物之一，也是衡量細胞膜損傷程度的一個重要指標，MDA 含量越高，表示細胞損傷越嚴重。與相對電導率結果一致，不同三唑類殺菌劑處理可導致FocTR4 菌絲細胞內MDA 含量顯著上升（圖3），0.25、1.25 μg/mL 丙硫菌唑處理組菌絲的MDA 含量分別為對照組（5.43 nmol/mL）的1.72 倍（9.33 nmol/mL，P＜0.001）和2.50 倍（13.57 nmol/mL，P＜0.001）；0.25、1.25 μg/mL戊唑醇處理組菌絲的MDA含量分別為對照組的1.53 倍（8.33 nmol/mL，P＜0.001）和2.24倍（12.18 nmol/mL，P＜0.001）；0.25、1.25 μg/mL丙環(huán)唑處理組菌絲的MDA 含量分別為對照組MDA 的1.40 倍（7.62 nmol/mL，P＜0.001）和1.96 倍（10.65 nmol/mL，P＜0.001）；0.25、1.25 μg/mL腈菌唑處理組菌絲的MDA 含量分別為對照組的1.40 倍（7.58 nmol/mL，P＜0.001）和1.93 倍（10.48 nmol/mL，P＜0.001）。以上結果表明，三唑類殺菌劑可促進MDA在FocTR4菌絲細胞中的積累，從而造成菌絲不同程度的膜脂過氧化損傷。

圖3 不同三唑類殺菌劑對Foc TR4 菌絲細胞內丙二醛含量的影響Fig.3 Effects of different triazole fungicides on concentrations of MDA in Foc TR4 hypha cells

2.3 三唑類殺菌劑與靶標CYP51 的相互作用

三唑類殺菌劑通過與靶標CYP51 作用抑制真菌生長，其與靶標酶的結合能力能夠在一定程度上反映他們的抑菌能力，因此，本研究通過明確不同三唑類殺菌劑處理對細胞色素P450酶活性、CYP51基因表達量的影響以及不同三唑類殺菌劑與CYP51的結合作用模式，探究不同三唑類殺菌劑對CYP51的影響。

如圖4A 所示，不同濃度的4 種殺菌劑對FocTR4 的細胞色素P450 酶活性均具有一定的誘導作用，且存在濃度依賴效應。其中，與對照組相比，0.25 μg/mL丙硫菌唑和0.25 μg/mL戊唑醇可誘導細胞色素P450 酶活性顯著上升（P＜0.05），細胞色素P450酶活性分別為對照組的1.70倍和1.45倍；當處理組的質量濃度為1.25 μg/mL 時，4 種殺菌劑均可顯著誘導細胞色素P450 酶活性上升（P＜0.05）。FocTR4 含 有3 個CYP51基 因，分 別 為CYP51-1（FOIG_04112）、CYP51-2（FOIG_10483）和CYP51-3（FOIG_00890）。采用實時熒光定量PCR 測定不同三唑類殺菌劑對基因表達量的影響，結果（圖4）表明：4 種殺菌劑均可顯著上調CYP51-1（圖4B）和CYP51-3（圖4D）基因的表達量（P＜0.05），而對CYP51-2基因的表達量均無顯著影響（圖4C）。CYP51-1基因表達量的上調范圍在0.25 μg/mL處理下為1.29～5.66 倍，在1.25 μg/mL 處理下為3.04～16.08 倍。CYP51-3基因表達量上調范圍在0.25 μg/mL 處理下為1.26～6.30 倍，在1.25 μg/mL 處理下為3.50～13.93 倍。與對細胞色素P450 酶活性的誘導作用不同，在相同處理濃度下，丙硫菌唑對CYP51-1和CYP51-3基因表達水平的影響顯著小于戊唑醇和丙環(huán)唑。

圖4 不同三唑類殺菌劑對細胞色素P450酶活性和靶標基因CYP51相對表達量的影響Fig.4 Effects of different triazole fungicides on cytochrome P450 enzyme activity and relative expression levels of target gene CYP51

為進一步研究不同三唑類殺菌劑與CYP51的作用模式，本研究采用Modeller 9.20 軟件對CYP51-3蛋白質序列進行同源建模（圖5A），利用AutoDock Vina 工具分析不同三唑類殺菌劑與CYP51 蛋白質的作用模式（圖5B）。采用盲對接將每個配體進行5次對接，每次產生20個構象，窮盡性為100。通過分析它們與CYP51的結合構象發(fā)現，這4種殺菌劑均位于同樣的結合口袋，共同作用的疏水殘基有Ala308、Phe511、Tyr123、Phe131、Ile374 和Tyr137（圖6）。結合自由能反映了配體和受體間的結合程度，結合自由能越小，分子間結合就越緊密。通過分子對接得到丙硫菌唑、戊唑醇、腈菌唑和丙環(huán)唑同源建模結構的結合自由能分別為－33.79、－34.92、－34.08、－36.59 kJ/mol，上述4 種殺菌劑三氮唑環(huán)上的N 原子與CYP51 的血紅素鐵之間的距離分別為2.66、2.81、2.67、2.96 ?（表2）。

圖5 Foc TR4中CYP51蛋白質的同源建模Fig.5 Homology modeling of CYP51 protein from Foc TR4

圖6 分子對接預測4種三唑類殺菌劑與CYP51的結合模式Fig.6 Binding modes of four triazole fungicides with CYP51 predicted by molecular docking

為驗證不同三唑類殺菌劑與CYP51 的互作模式，分別對4 種殺菌劑與CYP51 蛋白質進行SPR 親和力測試。如圖7 所示，CYP51 蛋白質的響應值（response unit,RU）隨殺菌劑濃度的升高而上升，表現出很好的濃度依賴效應。4種殺菌劑與CYP51的結合呈現出2 種模式：丙硫菌唑和腈菌唑在進樣幾秒鐘內即可以取得穩(wěn)態(tài)反應，進樣完畢藥物解離迅速，響應信號值立即回歸基線，呈現出“快上快下”的結合模式；而戊唑醇和丙環(huán)唑與CYP51的結合呈現出“慢上快下”結合模式。這些結果證實了CYP51 蛋白偶聯芯片具有良好的活性和特異性。利用1∶1穩(wěn)態(tài)親和力模型對4種殺菌劑與CYP51的結合情況進行擬合，得到了4 種殺菌劑與CYP51 蛋白質的親和動力常數（KD）（表2），分別為2.87×10－7mol/L（丙硫菌唑）、2.42×10－9mol/L（戊唑醇）、7.86×10－8mol/L（丙環(huán)唑）和2.15×10－6mol/L（腈菌唑）。

圖7 不同三唑類殺菌劑與CYP51的結合動力學曲線Fig.7 Binding kinetic curves of different triazole fungicides with CYP51

表2 不同三唑類殺菌劑與CYP51的分子對接和親和動力學參數Table 2 Molecular docking and affinity-kinetic parameters of different triazole fungicides with CYP51

3 討論與結論

三唑類殺菌劑因其優(yōu)良的保護、治療作用及內吸性功能已被廣泛用于多種植物病原菌的防治。但目前尚無關于不同三唑類殺菌劑防治香蕉枯萎病效果差異比較的系統性報道。本研究測定了4種常用三唑類殺菌劑對香蕉枯萎病菌FocTR4 的抑菌活性，結果表明：不同三唑類殺菌劑對FocTR4菌絲生長抑制效果差異明顯，其中丙硫菌唑的抑菌活性最高，EC50為0.79 μg/mL，而相同條件下腈菌唑EC50為10.13 μg/mL，二者抑菌活性相差近12倍。由此可見三唑類殺菌劑對FocTR4 的作用具有一定的選擇性。細胞膜具有多種生物功能，可為細胞的生命活動提供相對穩(wěn)定的內部環(huán)境，是真菌自我保護的屏障。病原菌細胞膜受損可導致細胞膜屏障作用被破壞，細胞內物質外泄，外源物質也可以隨意進入細胞，最終導致病原菌形態(tài)發(fā)生改變。三唑類殺菌劑通過抑制真菌甾醇生物合成中CYP51活性，引起麥角甾醇合成受阻以及甾醇前體物質的積累，從而破壞細胞膜的正常功能，在不同程度上對細胞膜的通透性、完整性以及細胞的形態(tài)等各方面產生影響[19]，最終導致細胞死亡。本研究發(fā)現，使用三唑類殺菌劑尤其是丙硫菌唑處理可導致FocTR4菌絲出現表面干癟、凹陷和扁平化現象。電導率和MDA含量是衡量細胞膜通透性和細胞膜損傷程度的重要指標。與丙環(huán)唑對核盤菌[20]，戊唑醇、三唑醇對小麥紋枯病菌（Rhizoctonia cerealis）菌絲細胞膜的影響結果[21]一致，丙硫菌唑、戊唑醇、丙環(huán)唑和腈菌唑均可以引起FocTR4 菌絲細胞內電解質外滲，使細胞膜通透性提高，表明三唑類殺菌劑對FocTR4 的生長影響不僅表現在對菌絲微觀形態(tài)的改變上，也體現在對菌絲細胞膜通透性的改變上，從而引起細胞內一系列生理生化反應，對菌絲生長產生不利影響。同時，相較于其他3種三唑類殺菌劑，丙硫菌唑對FocTR4的作用效果是最優(yōu)的。

目前，有關不同三唑類殺菌劑抑菌活性差異作用機制，尤其是不同分子結構與其靶標的互作模式仍缺乏進一步研究。本研究中，三唑類殺菌劑處理可引起FocTR4菌絲中細胞色素P450酶活性上升，表明4 種三唑類殺菌劑對細胞色素P450 酶具有誘導作用。動植物或真菌中都不只含有1個CYP51基因[22]。如煙曲霉（Aspergillus fumigatus）[23]以及稻瘟病菌（Magnaporthe oryzae）[24]中均含有2 個同源CYP51基因，尖孢鐮刀菌（Fusarium oxysporum）[19]、禾谷鐮刀菌（Fusarium graminearum）[25]以及黃曲霉（Aspergillus flavus）[26]中則均含有3 個CYP51基因。對細胞色素P450 基因的異源表達和調控的研究結果證明，P450 有廣泛的底物誘導特性[27]，如戊唑醇可誘導禾谷鐮刀菌[25]和亞洲鐮刀菌（Fusarium asiaticum）[28]的CYP51A和CYP51B基因的過量表達；三唑類殺菌劑可誘導稻瘟病菌[24]和指狀青霉菌（Penicillium digitatum）[29]中CYP51A基因表達量上調。同樣地，本研究中，與對照組相比，4種三唑類殺菌 劑 可誘 導FocTR4 中CYP51-1和CYP51-3的 表達量顯著性上調，表明靶標基因CYP51的較高表達是4 種三唑類殺菌劑誘導的結果。許多研究指出，CYP51 在病原菌對三唑類殺菌劑的敏感性中發(fā)揮著重要作用。如與野生型菌株相比，ΔCYP51A的稻瘟病菌和禾谷鐮刀菌菌株對三唑類殺菌劑的敏感性顯著增加；稻瘟病菌的CYP51A基因被敲除后，其對多種三唑類殺菌劑敏感性提高[24-25]。以病原菌靶標CYP51的三維結構為基礎，解析其與三唑類殺菌劑互作的分子機制，有助于進一步闡明FocTR4 對三唑類殺菌劑敏感性差異的機制。FocTR4 的3 個CYP51 中，CYP51-3 和 鐮 刀 菌Fusarium musae的CYP51B 氨基酸序列一致性高達98.67%[30]，在已報道的絲狀真菌中，CYP51B被認為主要負責編碼甾醇14α-脫甲基酶[31]?；诖?，本研究選取了CYP51-3氨基酸序列進行同源建模，并進一步對4 種殺菌劑與靶標CYP51-3的分子對接模式進行分析。結果表明，4 種殺菌劑與CYP51 的結合自由能分別為－33.79 kJ/mol（丙硫菌唑）、－34.92 kJ/mol（戊唑醇）、－34.08 kJ/mol（腈菌唑）和－36.59 kJ/mol（丙環(huán)唑），其中CYP51與丙環(huán)唑的親和力最強，與丙硫菌唑的親和力最弱，但不同殺菌劑之間的結合自由能無顯著性差異。因此本研究進一步采用SPR手段對4種殺菌劑與CYP51蛋白進行親和力檢測。結果發(fā)現，4種殺菌劑均顯示出較好的結合力。與分子對接結果相一致，與戊唑醇和丙環(huán)唑相比，丙硫菌唑與CYP51的結合模式不同，表現為“快上快下”模式，且結合能力較弱，其親和動力常數為2.87×10－7mol/L。該結果與對CYP51基因表達量的研究結論一致，丙硫菌唑對CYP51-1和CYP51-3基因表達的誘導水平顯著小于戊唑醇和丙環(huán)唑。4種殺菌劑中，丙硫菌唑的作用效果優(yōu)于其他3種殺菌劑。推測原因可能是與其他三唑類殺菌劑相比，丙硫菌唑與CYP51的作用方式存在特異性，戊唑醇和丙環(huán)唑是非競爭性抑制劑，而丙硫菌唑是競爭性抑制劑。PARKER等[32]的研究結果表明，當戊唑醇、三唑醇、丙環(huán)唑與禾生球腔菌（Mycosphaerella graminicola）的CYP51（MgCYP51）緊密結合時，產生強的Ⅱ型差異光譜（峰值在423～429 nm 波長處，波谷在406～409 nm波長處），表明形成了典型的低自旋六配位體復合物。而丙硫菌唑和MgCYP51蛋白質的結合具有不同的吸收峰，且結合能力比戊唑醇等其他殺菌劑低800多倍[32]，該研究結果也很好地證實了上述推測。

與其他三唑類殺菌劑相比，丙硫菌唑具有更廣譜的殺菌活性。硫酮結構的引入賦予了丙硫菌唑獨特的性能以及良好的內吸活性。因此，丙硫菌唑對于香蕉枯萎病的防治效果值得進一步的探究，這對于擴大丙硫菌唑的應用范圍和高效防治香蕉枯萎病具有重要意義。在此基礎上，本課題組下一步的研究工作將對FocTR4 的CYP51基因進行突變，探究其與三唑類殺菌劑互作中發(fā)揮作用的關鍵氨基酸位點，以期為香蕉枯萎病的田間防治提供理論支持。