《Computational Biology and Chemistry》：Elucidating the Mechanisms of Aristolochic Acid-Induced Upper Tract Urothelial Carcinoma: A Multi-Omics Approach Combining Bioinformatics and Computational Modeling
編輯推薦：
為系統(tǒng)闡明馬兜鈴酸(AA)誘導(dǎo)上尿路尿路上皮癌(UTUC)的致癌機(jī)制，研究人員整合網(wǎng)絡(luò)毒理學(xué)、機(jī)器學(xué)習(xí)、分子對接與分子動力學(xué)模擬等計算生物學(xué)方法，深入探索了AA致UTUC的潛在分子作用靶點(diǎn)與通路網(wǎng)絡(luò)。研究發(fā)現(xiàn)AA可能通過直接結(jié)合EGFR、PTGS2等核心靶點(diǎn)，協(xié)同經(jīng)典DNA加合物介導(dǎo)的基因毒性通路，共同驅(qū)動UTUC的發(fā)生發(fā)展。該研究為理解AA相關(guān)的致癌機(jī)制提供了新視角，并為防治策略的開發(fā)奠定了理論基礎(chǔ)。
在全球健康意識提升的浪潮下，天然藥物的使用日益廣泛，但其中的“隱形殺手”不容忽視。馬兜鈴酸，這種存在于某些傳統(tǒng)草藥中的天然化合物，已被世界衛(wèi)生組織國際癌癥研究機(jī)構(gòu)明確列為I類人類致癌物，與一種惡性程度較高的癌癥——上尿路尿路上皮癌密切相關(guān)。長久以來，科學(xué)家們普遍認(rèn)為，馬兜鈴酸的“罪證”在于它能與細(xì)胞的“生命藍(lán)圖”DNA結(jié)合，形成難以清除的“烙印”（DNA加合物），從而引發(fā)基因突變，最終導(dǎo)致細(xì)胞癌變。然而，隨著研究的深入，一個謎團(tuán)逐漸浮現(xiàn)：除了這套廣為人知的、直接攻擊DNA的“基因毒性”劇本，馬兜鈴酸是否還扮演著其他不為人知的角色？它能否像一把“萬能鑰匙”，直接“插”入細(xì)胞信號通路中的關(guān)鍵蛋白質(zhì)“鎖孔”，開啟或關(guān)閉致癌開關(guān)，與基因毒性機(jī)制“里應(yīng)外合”，共同導(dǎo)演癌癥這出悲劇？這個超越單一機(jī)制的、多靶點(diǎn)的致癌網(wǎng)絡(luò)，正是當(dāng)前研究的空白與挑戰(zhàn)。
為了回答這個復(fù)雜的問題，Tongpeng Liu、Yu Yao、Yuan Xu、Haojun Zhang、Lijiang Sun、Guiming Zhang這組研究人員展開了一場深入的“虛擬偵查”。他們繞開了傳統(tǒng)耗時費(fèi)力的實(shí)驗(yàn)摸索，而是巧妙地運(yùn)用了一系列強(qiáng)大的計算生物學(xué)工具，對馬兜鈴酸的致癌“暗網(wǎng)”進(jìn)行了系統(tǒng)性排查與模擬推演。他們的研究成果最終發(fā)表在《Computational Biology and Chemistry》期刊上，為我們理解馬兜鈴酸的“多重身份”提供了全新的線索。
為開展研究，作者團(tuán)隊(duì)主要運(yùn)用了以下幾項(xiàng)關(guān)鍵技術(shù)：首先，通過整合網(wǎng)絡(luò)毒理學(xué)方法，從ChEMBL、STITCH、SwissTargetPrediction、OMIM、GeneCards等公共數(shù)據(jù)庫中篩選并交叉比對，獲得與馬兜鈴酸及上尿路尿路上皮癌共同相關(guān)的潛在靶點(diǎn)。其次，利用Cytoscape軟件構(gòu)建和分析蛋白質(zhì)-蛋白質(zhì)相互作用網(wǎng)絡(luò)，并結(jié)合多種拓?fù)鋵W(xué)算法（如MCC、MNC、Degree、Closeness）鑒定核心靶點(diǎn)。接著，采用機(jī)器學(xué)習(xí)模型（基于LASSO回歸和隨機(jī)森林的嵌套交叉驗(yàn)證策略）對核心靶點(diǎn)的預(yù)測效能進(jìn)行驗(yàn)證和解釋，所使用的訓(xùn)練與驗(yàn)證數(shù)據(jù)來源于公共基因表達(dá)數(shù)據(jù)集GEO GSE166716。然后，運(yùn)用分子對接技術(shù)（AutoDock Vina）模擬馬兜鈴酸與核心靶點(diǎn)蛋白的結(jié)合模式與親和力。最后，通過分子動力學(xué)模擬（GROMACS軟件，100納秒）和分子力學(xué)/泊松-玻爾茲曼表面積結(jié)合自由能計算，深入評估和驗(yàn)證關(guān)鍵復(fù)合物的動態(tài)穩(wěn)定性和熱力學(xué)結(jié)合強(qiáng)度。
研究結(jié)果：
1. 與AA誘導(dǎo)UTUC相關(guān)靶點(diǎn)的篩選
通過數(shù)據(jù)庫交叉比對，研究人員從570個馬兜鈴酸相關(guān)靶點(diǎn)和1410個UTUC相關(guān)靶點(diǎn)中，精確定位了97個兩者共享的潛在關(guān)鍵靶點(diǎn)，提示這些靶點(diǎn)可能參與了AA誘導(dǎo)UTUC的病理過程。
2. GO和KEGG富集分析
對97個潛在靶點(diǎn)的功能富集分析揭示，它們顯著富集于脂質(zhì)代謝、外源性物質(zhì)解毒、氧化還原酶活性等生物學(xué)過程，以及癌癥相關(guān)通路（如PI3K-Akt信號通路）、細(xì)胞色素P450介導(dǎo)的藥物代謝和細(xì)胞凋亡等關(guān)鍵通路。這提示AA可能通過擾亂細(xì)胞代謝、壓力應(yīng)答和生存死亡平衡來驅(qū)動腫瘤發(fā)生。
3. PPI網(wǎng)絡(luò)構(gòu)建與核心靶點(diǎn)鑒定
基于97個靶點(diǎn)構(gòu)建的蛋白質(zhì)相互作用網(wǎng)絡(luò)具有高連接密度。通過多種拓?fù)鋵W(xué)算法和交叉驗(yàn)證，最終鑒定出六個共識核心基因：CASP3、EGFR、PARP1、CYCS、PTGS2和HSP90AA1。它們在網(wǎng)絡(luò)中處于樞紐地位，可能是AA發(fā)揮作用的關(guān)鍵節(jié)點(diǎn)。
4. 基于機(jī)器學(xué)習(xí)的核心基因篩選與模型可解釋性
利用公開數(shù)據(jù)集GSE166716，研究人員通過嵌套交叉驗(yàn)證的隨機(jī)森林模型對六個候選基因進(jìn)行驗(yàn)證。模型顯示出卓越的區(qū)分性能。穩(wěn)定性分析進(jìn)一步將核心基因鎖定為五個：CASP3、EGFR、PARP1、PTGS2和HSP90AA1。SHAP分析表明，CASP3的高表達(dá)與“正常”類別相關(guān)，而PARP1和HSP90AA1的高表達(dá)則傾向于推動模型預(yù)測為“癌癥”類別，這為這些基因在UTUC中的可能作用提供了計算層面的支持。
5. AA與核心靶點(diǎn)蛋白的分子對接
分子對接預(yù)測顯示，馬兜鈴酸與五個核心靶點(diǎn)蛋白（CASP3、EGFR、PARP1、PTGS2、HSP90AA1）均具有較高的結(jié)合親和力，其中與PTGS2（結(jié)合能為-9.3 kcal/mol）和EGFR（結(jié)合能為-8.2 kcal/mol）的結(jié)合最為穩(wěn)定。詳細(xì)的相互作用模式圖展示了氫鍵、疏水相互作用等多種分子間作用力，從原子層面揭示了可能的結(jié)合機(jī)制。
6. AA與核心靶點(diǎn)蛋白的分子動力學(xué)模擬
以結(jié)合親和力預(yù)測較高的AA-EGFR復(fù)合物為代表，進(jìn)行了長達(dá)100納秒的分子動力學(xué)模擬。結(jié)果顯示，復(fù)合物的結(jié)構(gòu)整體穩(wěn)定，構(gòu)象緊湊，溶劑可及表面積呈下降趨勢，回轉(zhuǎn)半徑保持穩(wěn)定，均方根偏差在平衡后波動較小。這些動力學(xué)指標(biāo)共同證實(shí)了AA與EGFR在模擬時間尺度內(nèi)能夠形成穩(wěn)定的復(fù)合物。自由能景觀圖揭示了復(fù)合物存在多個穩(wěn)定的低能量構(gòu)象狀態(tài)。
7. AA-EGFR復(fù)合物的結(jié)合自由能計算
基于分子動力學(xué)模擬軌跡，采用MM/PBSA方法對AA-EGFR復(fù)合物的結(jié)合自由能進(jìn)行了定量計算。結(jié)果顯示，其結(jié)合自由能ΔGbind為-55.68 ± 1.00 kcal/mol，表明兩者結(jié)合是一個高度自發(fā)的熱力學(xué)過程。能量分解顯示，強(qiáng)大的靜電相互作用是結(jié)合的主要驅(qū)動力，盡管被極性溶劑化能抵消了大部分，但范德華力和非極性溶劑化貢獻(xiàn)提供了重要的穩(wěn)定作用。這從熱力學(xué)角度為兩者的強(qiáng)結(jié)合提供了定量證據(jù)。
研究結(jié)論與重要意義：
本研究通過整合網(wǎng)絡(luò)毒理學(xué)、機(jī)器學(xué)習(xí)、分子對接與分子動力學(xué)模擬等多層次計算策略，首次系統(tǒng)描繪了馬兜鈴酸誘導(dǎo)上尿路尿路上皮癌的潛在多靶點(diǎn)致癌網(wǎng)絡(luò)。研究不僅重申了馬兜鈴酸通過形成DNA加合物（基因毒性）這一經(jīng)典致癌途徑的核心地位，更重要的是，提出并提供了計算證據(jù)支持一個全新的協(xié)同機(jī)制假說：馬兜鈴酸或其代謝物可能像一把“特洛伊木馬”，直接、高親和力地結(jié)合并潛在調(diào)控EGFR、PTGS2、PARP1等關(guān)鍵信號蛋白，從而干擾細(xì)胞增殖、凋亡、DNA修復(fù)和炎癥反應(yīng)等核心通路。
這種直接的“非基因毒性”作用，可能與它造成的DNA損傷和突變（基因毒性）產(chǎn)生“協(xié)同放大”效應(yīng)。例如，DNA損傷激活的生存信號可能因馬兜鈴酸對EGFR的直接作用而被進(jìn)一步加強(qiáng)；凋亡執(zhí)行蛋白CASP3被結(jié)合可能讓受損細(xì)胞“逃過一死”，積累更多突變；與促炎酶PTGS2結(jié)合可能持續(xù)“煽風(fēng)點(diǎn)火”，營造利于腫瘤生長的慢性炎癥微環(huán)境。這種基因毒性與非基因毒性雙管齊下的模式，可能更全面地解釋了為何馬兜鈴酸相關(guān)的UTUC往往表現(xiàn)出更高的侵襲性。
該研究的結(jié)論具有重要的理論意義和潛在的轉(zhuǎn)化價值。在理論上，它突破了將馬兜鈴酸致癌性僅歸于DNA加合物的單一視角，為理解復(fù)雜致癌物的多模式作用機(jī)制提供了新的分析框架和假設(shè)。在應(yīng)用上，所鑒定的核心靶點(diǎn)（如EGFR、PARP1、PTGS2）本身已是重要的抗癌藥物靶點(diǎn)，這為開發(fā)針對馬兜鈴酸暴露人群的化學(xué)預(yù)防策略或聯(lián)合治療方案提供了新的線索和思路。當(dāng)然，這些基于計算的發(fā)現(xiàn)仍需后續(xù)濕實(shí)驗(yàn)的嚴(yán)格驗(yàn)證，以最終確定馬兜鈴酸與這些蛋白的直接相互作用及其在致癌過程中的確切功能角色。這項(xiàng)研究猶如一幅精心繪制的“嫌犯關(guān)系網(wǎng)絡(luò)圖”，為后續(xù)的實(shí)驗(yàn)偵探們指明了深入調(diào)查的關(guān)鍵方向。
