自抗基因導向的微生物天然產(chǎn)物挖掘

董 華,明燈明

(1.南京工業(yè)大學(xué) 藥學(xué)院,江蘇 南京 211800;2.南京工業(yè)大學(xué) 生物與制藥工程學(xué)院,江蘇 南京 211800)

天然產(chǎn)物(natural products,NPs)是源自生物體次級代謝的有機小分子,其高活性及成藥潛力一直備受人們關(guān)注,絕大多數抗癌、抗感染和抗菌藥物都來(lái)源于天然產(chǎn)物及其衍生物。天然產(chǎn)物主要存在于微生物和植物中,而微生物相較于植物因其較小的基因組和易培養的特性,所以人們對微生物來(lái)源的天然產(chǎn)物的研究更多且更深入。對天然產(chǎn)物生物合成基因的廣泛研究表明,編碼合成天然產(chǎn)物所需酶的基因其在基因組中通常緊密相連地排列在一起,形成所謂的生物合成基因簇(biosynthetic gene clusters,BGCs)。這些基因編碼天然產(chǎn)物合成過(guò)程所需要的不同的酶,共同負責特定天然產(chǎn)物的生物合成[1]。近年來(lái),基因組測序技術(shù)的快速發(fā)展徹底改變了天然產(chǎn)物的發(fā)現模式。借助生物信息學(xué)工具[2-3],人們已從微生物基因組中鑒定出了數以萬(wàn)計的BGCs,其數目遠遠超過(guò)了實(shí)驗表征的天然產(chǎn)物數量。在已鑒定的BGCs中,據估計只有不到10%與生物活性天然產(chǎn)物相關(guān),另外90%以上的基因簇產(chǎn)物未知[4]。這些BGCs的共同特征是轉錄沉默,從而形成所謂的生物合成暗物質(zhì)?，F代微生物基因組學(xué)的快速發(fā)展,催生了數據驅動(dòng)的天然產(chǎn)物挖掘方法,它使用模式識別技術(shù)在海量的基因組數據中尋找暗物質(zhì),挖掘新天然產(chǎn)物,從而加速基于天然產(chǎn)物的藥物發(fā)現。

大多數天然產(chǎn)物是生物體產(chǎn)生的,通過(guò)抑制必需的管家酶來(lái)殺死或限制競爭性生物體生長(cháng)。如果相同的酶靶標在生產(chǎn)者中是保守且必需的,那么天然產(chǎn)物對生產(chǎn)者也會(huì )有同樣的毒性。因此,生產(chǎn)者必須存在自我保護機制,以賦予自身對該合成天然產(chǎn)物的抗性[5]。編碼自抗性機制的基因是微生物生產(chǎn)者的特征性狀,而自抗性的表達通常不利于生產(chǎn)者自身的生長(cháng)與存活,因此自抗性的表達通常僅在天然產(chǎn)物生物合成時(shí)一同表達[6-7]。生產(chǎn)者的自抗性策略,一般包括其對天然產(chǎn)物本身的修飾和降解,對作用天然產(chǎn)物靶標的修飾、保護以及替換,以及廣泛利用轉運蛋白外排天然產(chǎn)物等[8]。在充分掌握天然產(chǎn)物與靶標相互作用的基礎上,人們可以通過(guò)計算識別抗性基因及其作用機制來(lái)甄別新的BGCs。

綜述了微生物生產(chǎn)者在天然產(chǎn)物生物合成過(guò)程中的自我抗性機制,整理了近年來(lái)文獻報道的抗性基因數據庫和相關(guān)計算工具。文章簡(jiǎn)述了通過(guò)識別自抗性基因從微生物中發(fā)現產(chǎn)生天然產(chǎn)物的潛在 BGCs的方法,并展望了開(kāi)發(fā)結合基因組注釋和蛋白質(zhì)結構進(jìn)行天然產(chǎn)物的挖掘方法。

自抗性通常是基因編碼的,可以通過(guò)多種機制表現,包括現有基因的過(guò)度表達、復制和由點(diǎn)突變獲得全新基因。新的測序技術(shù)和計算方法,使得人們可以在海量的細菌和真菌基因組中快速地鑒定和表征天然產(chǎn)物生物合成基因簇,同時(shí)挖掘出越來(lái)越多的自抗性基因,從而對各種抗性作用機制獲得比較深入的認識[9]。細菌(真菌)已經(jīng)進(jìn)化出多種復雜的機制,可以針對自身的活性小分子進(jìn)行自我防御;很多時(shí)候它們同時(shí)擁有多種抗性機制,以確保完全抵御它們產(chǎn)生的生物活性分子傷害自身[10]。這些機制一般包括:外排天然產(chǎn)物、產(chǎn)物修飾、產(chǎn)物降解、天然產(chǎn)物靶標修飾等[11](圖1)。本文重點(diǎn)介紹微生物在合成天然產(chǎn)物過(guò)程中激發(fā)的典型的自抗機制類(lèi)型,并為每個(gè)類(lèi)型提供了具體示例。

圖1 自抗性機制[11]

1.1 外排泵

外排泵與所有原核或真核細胞對有毒化合物的抗性相關(guān),它們存在于所有生物中,從人類(lèi)(對抗癌藥物的抗性)到細菌(對抗生素的抗性)[12]。對于產(chǎn)生天然產(chǎn)物的微生物生產(chǎn)者來(lái)說(shuō),天然產(chǎn)物的外排是一種普遍的自抗機制,通常會(huì )與其他機制一起發(fā)生,例如,天然產(chǎn)物或靶標的修飾[10]。細菌通常擁有多藥外排泵,從細胞中主動(dòng)輸出活性小分子,已經(jīng)發(fā)現了一些外排泵,它們基因的突變和過(guò)度表達會(huì )增加有毒物質(zhì)從細胞內流出,使細胞內毒物濃度保持在較低水平,防止達到抑制濃度[13]。

研究發(fā)現,外排泵可分為5個(gè)結構家族(表1)。腺嘌呤核苷三磷酸結合盒(adenosine triphosphate binding cassette,ABC)轉運蛋白,利用腺嘌呤核苷三磷酸(adenosine triphosphate,ATP)的結合和水解來(lái)轉運底物,它們是一個(gè)廣泛的蛋白質(zhì)家族,承擔原核生物和真核生物中物質(zhì)的輸入和輸出[14]。主要協(xié)助轉運蛋白超家族(the major facilitator superfamily,MFS)包括一組轉運蛋白,使用質(zhì)子動(dòng)力作為運輸底物的能源[15]。多藥與有毒物質(zhì)外排轉運蛋白(the multidrug and toxic compound extrusion,MATE)分為3個(gè)亞家族,基于它們的氨基酸序列相似性,分為多藥耐藥蛋白NorM(multidrug resistance protein NorM)、DNA損傷誘導蛋白DinF(DNA damage-inducible protein F)和真核亞家族[16-17]。MATE轉運蛋白在功能上更通用,可以使用 H+或 Na+電化學(xué)梯度輸出藥物[18]。葡萄球菌多藥耐藥轉運蛋白(the staphylococcal multi-resistance,SMR)使用質(zhì)子動(dòng)力來(lái)排出有毒化合物[19]?？剐越Y瘤和細胞分裂超家族(the resistance nodulation and cell division superfamily,RND)外排泵廣泛分布在革蘭氏陰性微生物中,該復合物由3種不同的蛋白質(zhì)形成,分別是可利用質(zhì)子動(dòng)力排出底物且嵌入細胞內膜中的二級主動(dòng)外排泵蛋白、外膜蛋白和膜融合蛋白,將位于兩個(gè)膜上的蛋白質(zhì)跨周質(zhì)空間連接起來(lái)[20]。

表1 5種外排泵

在抗藥性外排泵系統中,研究得最多的系統之一是吖啶抗性蛋白A/B-外膜因子TolC(AcrAB-TolC)系統(圖2),屬于RND轉運蛋白家族,該泵組件包括外膜通道TolC、位于內膜的二級轉運蛋白AcrB和連接二者的周質(zhì)膜融合蛋白AcrA[21]。RND 泵起到質(zhì)子反向轉運蛋白的作用,AcrAB-TolC能夠以矢量方式運輸各種化學(xué)相似性很小的化合物,因此對廣譜抗生素具有抗性[22],能夠運輸多種底物,對四環(huán)素、氯霉素、某些β-內酰胺類(lèi)、新生霉素、夫西地酸和氟喹諾酮類(lèi)藥物產(chǎn)生抗性[23]。它是少數已經(jīng)獲得復合物完整結構的外排泵之一,在許多的革蘭氏陰性菌中都發(fā)現了RND的同源復合物,包括動(dòng)物和植物的病原體中[20],是生物代謝物和抗菌化合物的輸出者,在細菌抗性中發(fā)揮重要的作用。

圖2 AcrAB-TolC系統的結構[24](PDB:5ng5)

1.2 天然產(chǎn)物的修飾

給天然產(chǎn)物小分子結構上添加修飾(例如,磷酸鹽、乙酸鹽)是一種較為普遍且靈活的自抗性機制?；瘜W(xué)修飾包括O-?；?、N-?；?、O-磷酸化、O-核苷酸化、O-核糖基化、O-糖基化和硫醇轉移等。這些共價(jià)修飾策略都需要活性的共底物,包括ATP、乙酰輔酶A、煙酰胺腺嘌呤二核苷酸(NAD+)、二磷酸尿苷葡萄糖(UDPG)或谷胱甘肽。因此,這些酶通常只有在細胞質(zhì)中具有活性[25]。

博萊霉素BLM(bleomycin)是鏈霉菌(Streptomycesverticillus)產(chǎn)生的糖肽類(lèi)衍生抗生素?；騜lmB是在博萊霉素生物合成基因簇內發(fā)現的抗性元件,blmB編碼一種N-乙酰轉移酶[26],命名為BlmB,它通過(guò)將乙?；鶊F從乙酰輔酶A轉移到β-氨基丙氨酸部分的伯胺,即BLM金屬離子結合域的軸向配體,賦予宿主自抗性(圖3)。雖然乙?；煌耆柚菇饘倥cBLM結合,但可以抑制金屬螯合配位物與分子氧配位和還原分子氧,從而防止形成活化的BLM。

圖3 博萊霉素自抗性機制[26]

1.3 天然產(chǎn)物的降解

天然產(chǎn)物中經(jīng)常包含酰胺鍵和酯鍵等易被分解的化學(xué)鍵,其完整性對生物活性至關(guān)重要。生物體中有些酶可以靶向切割這樣的鍵,為生物體提供了一種新的自我抗性手段[25]。

β-內酰胺酶是其中最為典型的一種降解酶,它的作用是水解β-內酰胺類(lèi)抗生素中的β-內酰胺鍵。多種微生物體中都有編碼β-內酰胺酶的基因[27]。β-內酰胺酶本身并不是必需的細菌蛋白質(zhì),據推測它可能是從一些會(huì )產(chǎn)生β-內酰胺類(lèi)代謝物的細菌,如鏈霉菌中的青霉素結合蛋白(penicillin-binding proteins,PBPs)進(jìn)化而來(lái)的,因為這些細菌必須要有一些針對β-內酰胺類(lèi)抗生素的自我保護策略[28-29]?；诎被嵝蛄械南嗨菩?β-內酰胺酶分為A、B、C、D四類(lèi)[30],其中A 類(lèi)(青霉素水解)、C 類(lèi)(頭孢菌素水解)和 D 類(lèi)(苯唑西林水解)β-內酰胺酶的活性位點(diǎn)都是絲氨酸,稱(chēng)為絲氨酸-β-內酰胺酶,它們從β-內酰胺的細胞靶標PBPs進(jìn)化而來(lái),并在鍵裂解過(guò)程中使用催化絲氨酸[31](圖4(a))。B 類(lèi)β-內酰胺酶需要兩個(gè)鋅離子的活性,稱(chēng)為金屬-β-內酰胺酶,它們利用活性位點(diǎn)鋅離子來(lái)協(xié)調親核氫氧化物以促進(jìn)鍵斷裂[32](圖4(b))。因此,β-內酰胺酶采用兩種主要的分子策略來(lái)水解裂解抗生素的β-內酰胺環(huán):通過(guò)活性位點(diǎn)絲氨酸親核試劑的作用以及通過(guò)Zn2+中心激活水[30, 33-34]。

圖4 β-內酰胺酶作用機制[25]

1.4 分子靶標的修飾

自抗性還可以通過(guò)突變或修飾天然產(chǎn)物靶標蛋白來(lái)實(shí)現。一般情況下,這種自抗機制有很大的可能性是通過(guò)特定基因(管家酶)中的點(diǎn)突變產(chǎn)生的,從而使生產(chǎn)者獲得相對快速和簡(jiǎn)單的抗性,而且這種抗性對生產(chǎn)者適應性的影響最小。生產(chǎn)者會(huì )重新編碼出與管家酶功能等效的自抗酶,它屬于管家酶的變體;自抗酶在序列上與管家酶高度相似、仍保持自身原有的活性,但突變使自抗酶對天然產(chǎn)物不敏感。

地霉素(dityromycin)是一種肽抗生素,從土壤微生物鏈霉菌Streptomycessp.AM-2504的培養液中分離出來(lái)。結構研究表明,地霉素以30S 核糖體亞基中的核糖體蛋白 S12為靶標,通過(guò)阻斷tRNA來(lái)抑制蛋白質(zhì)合成[35]。宿主菌株對次級代謝產(chǎn)物地霉素的自抗性源于核糖體的特定修飾:位于抗生素結合位點(diǎn) S12 的蛋白高度保守區域有2個(gè)氨基酸替換,分別是Val36Thr和Arg59Lys[36],這些突變使地霉素對 30S 核糖體亞基的親和力顯著(zhù)降低,從而保護宿主菌株免受抗生素的毒性作用。

由于抗生素的廣泛應用,所以抗生素耐藥性機制研究最為普遍。目前,現有的抗生素抗性基因數據庫通常是網(wǎng)絡(luò )在線(xiàn)平臺,提供耐藥性相關(guān)的參考以及用于抗性基因序列分析和注釋的網(wǎng)絡(luò )工具。這些平臺接收用戶(hù)的核苷酸或蛋白質(zhì)序列查詢(xún)并返回抗性基因預測結果,包含帶有置信度值的注釋以及指向外部資源的鏈接[37]。下文描述了當前部分可用的抗生素抗性基因數據庫(表2)。

表2 抗性基因數據庫

在抗性基因相關(guān)的數據庫中,抗生素抗性基因數據庫(antibiotic resistance genes database,ARDB)是第一個(gè)建立的手動(dòng)管理的抗性基因序列資源,其中每個(gè)基因都注釋了抗性類(lèi)型、機制和本體[38]。目前不再更新,但所有的ARDB序列都集成到了抗生素耐藥性綜合數據庫(the comprehensive antibiotic resistance database,CARD)中。

CARD提供了與抗菌素耐藥性分子基礎相關(guān)的數據和模型,是抗性基因及其產(chǎn)物和相關(guān)表型的生物信息學(xué)數據庫[39]。目前擁有4 834個(gè)抗性基因本體、3 339個(gè)參考序列、1 788個(gè)單核苷酸多態(tài)性、2 774個(gè)已報道的抗性基因、3 385個(gè)抗菌檢測模型。CARD集成了許多預測和分析工具,包括用于分析和注釋序列的局部序列比對基本檢索工具(basic local alignment search tool,BLAST),用于預測抗性基因的抗性基因鑒定器(resistance gene identifier,RGI)。CARD 2020中最顯著(zhù)的更新是一個(gè)新的抗性突變體模塊,它提供了對來(lái)自88個(gè)病原體、9 560條染色體、21 362個(gè)質(zhì)粒、102 181個(gè)全基因組測序組裝和222 011個(gè)計算機預測分析的全基因組耐藥變體,RGI分析工具的存在也增強了CARD在預測基因組或宏基因組數據中抗生素耐藥性方面的效用[40]。

Resfam數據庫包含一系列精選的蛋白質(zhì)家族,而定義每個(gè)蛋白質(zhì)家族均使用了與某種抗性功能相關(guān)的蛋白質(zhì)的序列。這些序列經(jīng)一種叫隱馬爾可夫模型(hidden markov model,HMM)的統計算法即生成了具有某種抗性功能的蛋白質(zhì)代表模型。其使用抗性蛋白數據來(lái)自CARD、內酰胺酶工程數據庫(lactamase engineering database,LacED)和其他精選的β-內酰胺酶Lahey等。相對于其他一些主要關(guān)注病原體相關(guān)抗性基因的平臺,該數據庫還提供了抗性基因所處的環(huán)境概況的描述。與基于BLAST的搜索結果相比,其檢索結果對于土壤和人類(lèi)腸道微生物群中的抗性基因鑒定率提高了64%[41-43]。另一方面,由于用戶(hù)需要考慮到HMM方法可能會(huì )以較差的特異性為代價(jià),因此可能會(huì )產(chǎn)生大量的假陽(yáng)性預測,有時(shí)可能無(wú)法區分具有密切相關(guān)功能的抗性基因[37]。

ResFinder數據庫是在A(yíng)RDB數據庫的基礎上,附加收錄通過(guò)水平轉移獲得的耐藥基因數據而成的。此庫分為15種耐藥小庫、3 097條耐藥基因序列,分類(lèi)全面,更新及時(shí),目前已經(jīng)更新到4.1版本,并且支持本地化安裝使用。目前,除了可以查找分析已知的耐藥基因序列外,該數據庫可使用一種基于全基因組測序數據預測耐藥性染色體突變的工具(PointFinder)來(lái)識別染色體上的靶基因的點(diǎn)突變。這些突變信息的識別,對于研究耐藥突變的發(fā)生機制和新型抗生素的開(kāi)發(fā)等方面具有重要的意義。

由于β-內酰胺酶的范圍和臨床重要性,絕大多數專(zhuān)業(yè)數據庫都關(guān)注了該家族。β-內酰胺酶數據庫(beta-lactamase database,BLDB)包含所有類(lèi)別的β-內酰胺酶的信息,包括從文獻、國際生物技術(shù)中心信息檢索數據庫(national center for biotechnology information,NCBI)和蛋白質(zhì)數據庫PDB中的蛋白質(zhì)三維結構收集的抗性模式。該數據庫結合了動(dòng)力學(xué)和突變信息,為理解β-內酰胺酶的結構和功能之間的關(guān)系奠定了基礎。除了能夠訪(fǎng)問(wèn)結合配體的蛋白質(zhì)結構和理化特性的信息外,該數據庫還提供與外部蛋白質(zhì)和核苷酸資源的鏈接[44]。

產(chǎn)生天然產(chǎn)物的酶促途徑會(huì )在連續的生物合成基因簇中編碼。每個(gè)基因簇都編碼至少一個(gè)與其同源天然產(chǎn)物相關(guān)的抗性基因[51]?？剐曰虺霈F的位置通常與合成基因簇的位置非?？拷?存在于基因簇中或者出現在連續的基因簇附近,確保當天然產(chǎn)物開(kāi)始在細胞內積累時(shí),抗性基因能及時(shí)發(fā)揮作用進(jìn)行防御[52]。出現在基因簇中的抗性基因,其發(fā)生的位置呈無(wú)規則分布,目前還沒(méi)有研究揭示抗性基因出現的相對位置對天然產(chǎn)物有什么具體的影響。自 2017 年以來(lái),抗生素抗性目標搜尋器(antibiotic resistant target seeker,ARTS)促進(jìn)了靶向(自抗)基因組挖掘的方法,以便通過(guò)將管家酶和已知抗性基因與 BGC 鄰近、復制和水平基因轉移事件快速聯(lián)系起來(lái),優(yōu)先考慮生產(chǎn)具有假定新作用模式的抗生素的菌株[53-54],體現了自抗基因對天然產(chǎn)物挖掘的重要作用。

基因組測序技術(shù)和生物信息分析方法的最新進(jìn)展表明,人們使用傳統方法僅僅發(fā)現了微生物的天然產(chǎn)物的一小部分,嚴重低估了微生物次級代謝的化學(xué)廣度。將天然產(chǎn)物與其合成基因簇連接起來(lái)的研究方法,以及不斷積累的關(guān)于生物合成過(guò)程的邏輯知識,催生了基于基因組的天然產(chǎn)物挖掘這一新領(lǐng)域[11]。Tang[55]等研究表明,事先預判生物合成基因簇中的抗性基因,可以在研究發(fā)現天然產(chǎn)物的結構和作用機制之前確定它們的分子靶標[55]。這種挖掘策略不僅可以識別已知的抗性基因,而且可以發(fā)現潛在地靶向尚未表征的天然產(chǎn)物的新抗性基因。

根據已知活性的天然產(chǎn)物及其靶標定位到生物合成基因簇。Tang等[55]最早使用了這種挖掘策略,他們發(fā)現了一組脂肪酸合酶Ⅱ(fatty acid synthase Ⅱ,FASⅡ)抑制劑。通過(guò)篩查 86 個(gè)鹽孢菌屬菌株的基因組,確定了管家基因的重復拷貝,其中一些與FASⅡ基因有關(guān)。特別是tlmE基因,編碼飽和脂肪酸合成酶B/F(FabB/F)的同源物,它與合成FASⅡ抑制劑platensimycin和platencin的BGC的特征自抗蛋白PtmP3和PtnP3 有高度的序列相似性,所以tlmE被推定為自抗性基因。將該自抗基因相關(guān)聯(lián)的BGC異源表達,證明該BGC與已知的FASⅡ抑制劑硫代乳霉素TLM(thiolactomycin)[56]以及3種類(lèi)似化合物的產(chǎn)生有關(guān)。他們遵循類(lèi)似的方法從Streptomycesafghaniensis中識別出相關(guān)的 BGC,編碼兩個(gè)FabB/F 同源物,TtmE和TtmJ,推定為自抗性基因。異源表達表明該基因簇可以產(chǎn)生一系列新的TLM類(lèi)似物。Tang等[55]證實(shí)了推定的抗性基因TtmE和TtmJ發(fā)揮自我抗性。

針對選定的管家酶靶標,尋找特定功能的天然產(chǎn)物。Yan等[57]在最近的一項研究中發(fā)展了自抗基因導向的基因組挖掘技術(shù),以新的作用模式在真菌基因組數據庫中篩查新的除草劑先導化合物[57]?？紤]到該化合物的作用靶標是管家酶乙酰羥酸脫水酶(dihydroxy-acid dehydratase,DHAD),通過(guò)尋找與已知植物的DHAD相匹配的真菌密切同源物,獲得了含有編碼DHAD重復拷貝的BGC,推定該重復拷貝是自抗性基因。研究確定了一個(gè)保守的基因簇,其中編碼了一段DHAD的重復拷貝基因AstD,與所有真菌中存在的保守的管家酶DHAD具有 60% 的相似度。然后將來(lái)自土曲霉的生物合成核心基因在異源宿主中重組表達產(chǎn)生化合物,產(chǎn)生一種已知的倍半萜烯天然產(chǎn)物,生物學(xué)靶標未知。Yan等[57]證明了該化合物是DHAD的競爭性抑制劑,而重復拷貝的AstD確實(shí)對化合物具有抗性,是該化合物的自抗性基因。

通過(guò)自抗基因定位到生物合成基因簇進(jìn)而確定天然產(chǎn)物的生物活性。Panter等[58]使用相似的方法,用五肽重復蛋白(pentapeptide repeat proteins,PRP)作為靶標從粘細菌中發(fā)現了新型拓撲異構酶抑制劑。Baumann等[59]先是使用人們以往表征的模板序列,在其粘細菌基因組數據庫中搜索獲得了相似蛋白PRP;有研究表明該蛋白賦予針對拓撲異構酶的化合物的自抗性。隨后,他們在Pyxidicoccusfallax基因組中的聚酮生物合成基因簇中鑒定出了編碼 PRP的基因,該基因簇編碼合成兩個(gè)新的天然產(chǎn)物,pyxidicycline A和B。Panter等[58]研究證實(shí)了pyxidicycline A和B都是拓撲異構酶的選擇性抑制劑。盡管這里的抗性基因并不是靶向拓撲異構酶的直接重復拷貝,但該實(shí)驗驗證了使用自抗性基因發(fā)現具有所需生物活性的新天然產(chǎn)物的可能性[58]。

這些研究使用的共同策略是,使用所需的生物合成核心酶在基因組數據庫中搜索攜帶編碼推定的自抗性基因的BGC,然后注釋每個(gè)基因的功能,從而識別抗性基因。研究表明,這種通過(guò)自抗基因引導挖掘發(fā)現天然產(chǎn)物的方法確實(shí)可以發(fā)現具有新作用模式的天然產(chǎn)物[60]。

本文綜述的自抗基因導向的天然產(chǎn)物挖掘策略,其原理是抗菌化合物的生產(chǎn)者必須具有某種形式的固有抗性機制,以確保自我保護免受活性化合物的傷害。這些自我抗性基因通常出現在BGC中與代謝產(chǎn)物組裝的核心基因一起,這就使得這些抗性基因能夠與代謝產(chǎn)物一起表達,為天然產(chǎn)物生產(chǎn)者提供立即免疫。正如本文第3節中所描述的例子,這些抗性基因(ttmE、astD等)的功能、在基因組中出現的位置,都可以作為一種潛在的BGC合成的天然產(chǎn)物的生物活性和作用機制的挖掘策略。

然而,目前自抗基因導向的天然產(chǎn)物挖掘還面臨不少的挑戰?，F有的生物合成知識還比較有限,常常缺乏對于謝產(chǎn)物組裝的核心基因功能的完整表征,這使得研究人員們對于BGCs的預測有時(shí)可能不準確,對于BGC中自抗基因的預判也不準確[5]。例如,實(shí)際的天然產(chǎn)物可能并不是按照預測的方式參與自我抗性作用,也可能預測的管家酶同源物不參與自抗機制,而是作為生物合成所必需的合成酶而存在。這就需要大規模的拓展現有的BGCs知識庫,在比較基因組學(xué)的基礎上系統地研究BGCs的進(jìn)化規律,克服可能存在的基因組數據中BGCs元件碎片化和不完整性的認識問(wèn)題。隨著(zhù)微生物基因組序列信息數量的指數增長(cháng)和BGCs工作元件數量的不斷積累,相信計算機識別BGCs 的能力和精確度也會(huì )不斷增強。

準確地預測自抗基因,特別是對于那些可能在生物合成途徑中發(fā)揮關(guān)鍵作用的酶,正確地解析其結構與功能,這無(wú)疑是自抗基因導向的天然產(chǎn)物挖掘所面臨的一項挑戰任務(wù)。隨著(zhù)生物合成酶學(xué)知識的積累、計算模擬工具的改進(jìn),對自抗基因的預測準確度也有望得到快速提升。特別是,近兩年計算結構生物的重大突破,包括蛋白質(zhì)結構預測工具AlphaFold2和RoseTTAFold的開(kāi)源等,使得人們可以在更為準確的結構基礎上解析自抗基因的功能,分析其與小分子和其他蛋白質(zhì)相互作用。這為系統地研究自抗基因的功能提供了前所未有的機遇。有望基于自抗基因導向挖掘天然產(chǎn)物的方法,在不久的未來(lái)可能開(kāi)發(fā)出集序列、結構和其他特性為一體的工具,應用于天然產(chǎn)物的生物合成挖掘,并為新藥的研發(fā)儲備資源。