AlphaFold 蛋白質(zhì)結(jié)構(gòu)可視化5 行 Python 出 3D 結(jié)構(gòu)圖免費(fèi)完整指南【免費(fèi)下載鏈接】alphafoldOpen source code for AlphaFold 2.項(xiàng)目地址: https://gitcode.com/GitHub_Trending/al/alphafold用官方開(kāi)源的 AlphaFoldv2 的開(kāi)放實(shí)現(xiàn)把一段氨基酸序列喂進(jìn)去就能拿到一張可旋轉(zhuǎn)、可縮放、可著色的 3D 結(jié)構(gòu)圖。本文不繞彎子先讓你看到成品長(zhǎng)什么樣再給你一條最短跑通路徑最后按置信度著色、活性位點(diǎn)高亮、野生 vs 突變體對(duì)比三個(gè)場(chǎng)景把圖做有信息量。讀完你會(huì)獨(dú)立跑通渲染并知道每個(gè) pLDDT 閾值到底代表什么。先看圖這就是幾行代碼的產(chǎn)物左邊是實(shí)驗(yàn)測(cè)得的真實(shí)結(jié)構(gòu)右邊是模型預(yù)測(cè)出來(lái)的結(jié)構(gòu)兩幅幾乎重疊——這張對(duì)比圖背后的渲染就是后面用py3Dmol十幾行代碼能復(fù)現(xiàn)的東西。它不是渲染軟件的手工作品是預(yù)測(cè)結(jié)果的直接可視化。跑通最小示例從序列到 3D 結(jié)構(gòu)圖整條鏈路只有三步清洗序列 → 預(yù)測(cè)結(jié)構(gòu) → 渲染下面給你能跑通的最少代碼。先把序列洗干凈。AlphaFold 只認(rèn) 20 種標(biāo)準(zhǔn)氨基酸序列里不能有空格、不能太短from alphafold.notebooks import notebook_utils seq notebook_utils.clean_and_validate_single_sequence( input_sequenceMAAHKGAEHHHKAAEHHEQAAK..., # 換成你的序列 min_length16, max_length2500)跑完它會(huì)原樣返回干凈序列一旦混進(jìn)非氨基酸字符或長(zhǎng)度越界直接拋ValueError。接下來(lái)做結(jié)構(gòu)預(yù)測(cè)。真實(shí)環(huán)境需要下好模型參數(shù)和數(shù)據(jù)庫(kù)這里只寫(xiě)核心調(diào)用from alphafold.model import model, config, data cfg config.model_config(model_1) runner model.RunModel(cfg, data.get_model_haiku_params(model_1, params/)) pred runner.predict(runner.process_features(features, 0))這一步會(huì)把氨基酸坐標(biāo)算出來(lái)產(chǎn)物里最關(guān)鍵的是plddt——每個(gè)殘基的預(yù)測(cè)置信度。最后把坐標(biāo)交給py3Dmol變成可交互的 3D 圖import py3Dmol from alphafold.common import protein prot protein.from_prediction(features, pred) view py3Dmol.view(width800, height600) view.addModel(protein.to_pdb(prot), pdb) view.setStyle({cartoon: {color: spectrum}}) view.zoomTo(); view.show()上面這條彩虹絲帶圖就是view.show()后瀏覽器里能拖來(lái)拖去的那個(gè)窗口。把cartoon換成stick骨架棒狀或sphere空間填充球你看到的就是同一種結(jié)構(gòu)的三種畫(huà)法。讓圖更有信息量三個(gè)實(shí)用場(chǎng)景默認(rèn)渲染只能看個(gè)大概真正干活時(shí)要靠下面三招。① 按置信度著色pLDDT 染色。AlphaFold 把每個(gè)殘基的plddt寫(xiě)進(jìn)了 PDB 的 B 因子所以直接按b屬性映射四檔顏色即可bands {0: orange, 1: gold, 2: lightblue, 3: darkblue} view.setStyle({cartoon: {colorscheme: {prop: b, map: bands}}})深藍(lán)色區(qū)域pLDDT 90最可信橙色pLDDT 50基本是柔性無(wú)序區(qū)看結(jié)果時(shí)先盯深藍(lán)部分。② 高亮活性位點(diǎn)。想強(qiáng)調(diào)某幾個(gè)殘基用addStyle指定殘基號(hào)區(qū)間把它從卡通里挑出來(lái)view.addStyle({resi: (10, 19)}, {stick: {color: red, radius: 0.3}})③ 野生型 vs 突變體并排。用grid參數(shù)開(kāi)兩個(gè)子圖同屏對(duì)比構(gòu)象差異g py3Dmol.view(width800, height400, grid(1, 2)) g.addModel(wild_pdb, pdb, viewer(0, 0)) g.addModel(mut_pdb, pdb, viewer(0, 1)) g.setStyle({cartoon: {color: spectrum}}, viewer(0, 0)) g.zoomTo(); g.show()導(dǎo)出時(shí)記住一條view.png()存的是靜態(tài)死圖適合塞論文view._make_html()存的是保留交互的網(wǎng)頁(yè)適合發(fā)給同事自己拖。兩者別混用。踩坑速查表 現(xiàn)象 / 疑問(wèn)實(shí)際原因與對(duì)策ValueError: too short / too long單鏈16~2500殘基多鏈模型上限4000越界先截?cái)嗷驌Q模型圖上一段是橙色糊的那是pLDDT 50的低置信區(qū)別當(dāng)真是柔性無(wú)序區(qū)結(jié)構(gòu)有細(xì)小原子重疊AMBER relax 偶爾引入立體化學(xué)瑕疵不介意的話關(guān)掉run_relax預(yù)測(cè)卡很久MSA 數(shù)據(jù)庫(kù)搜索按分鐘到小時(shí)計(jì)序列越長(zhǎng)越久屬正常導(dǎo)出圖不能轉(zhuǎn)png()是靜態(tài)圖要可交互得導(dǎo)_make_html()配色不對(duì)pLDDT 藏在B 因子里染色務(wù)必用prop: b拿走就用序列先過(guò)clean_and_validate_single_sequence長(zhǎng)度守住 16~2500渲染前確認(rèn)plddt已寫(xiě)進(jìn) B 因子染色用prop: b論文用png()分享交互用_make_html()高置信區(qū)pLDDT 90優(yōu)先解讀低置信區(qū) 50只做參考官方完整示例見(jiàn) notebooks/AlphaFold.ipynb技術(shù)細(xì)節(jié)見(jiàn) docs/technical_note_v2.3.0.md【免費(fèi)下載鏈接】alphafoldOpen source code for AlphaFold 2.項(xiàng)目地址: https://gitcode.com/GitHub_Trending/al/alphafold創(chuàng)作聲明：本文部分內(nèi)容由AI輔助生成（AIGC），僅供參考