看到“schrodinger 薛定諤”這個(gè)關(guān)鍵詞很多人的第一反應(yīng)是那只既死又活的貓。但我作為一個(gè)常年跑計(jì)算模擬的人腦子里閃過(guò)的其實(shí)是另一個(gè)東西Schr?dinger Suite業(yè)內(nèi)直接叫“薛定諤”一套做分子模擬、藥物設(shè)計(jì)、材料科學(xué)的高通量計(jì)算平臺(tái)。薛定諤方程是量子力學(xué)的基石而 Schr?dinger 這家公司把這套物理模型做成了工程化軟件讓人可以直接在圖形界面里做分子對(duì)接、虛擬篩選、自由能微擾計(jì)算算得上是藥物發(fā)現(xiàn)領(lǐng)域繞不開(kāi)的“重型裝備”。這篇文章不聊量子力學(xué)科普就圍繞我在實(shí)際項(xiàng)目里用 Schr?dinger 套件做分子對(duì)接和虛擬篩選的經(jīng)驗(yàn)展開(kāi)。我會(huì)從最初的軟件選型邏輯講起一步一步把結(jié)構(gòu)準(zhǔn)備、Glide 對(duì)接參數(shù)、MM-GBSA 重打分、FEP 邊界這些核心環(huán)節(jié)拆開(kāi)再把實(shí)操中踩過(guò)的坑、總結(jié)的排查方法一并放出來(lái)。適合正在學(xué)計(jì)算化學(xué)、生物信息學(xué)或者剛拿到薛定諤 license 準(zhǔn)備跑第一個(gè)真實(shí)課題的人參考。1. 先說(shuō)清楚此薛定諤非彼薛定諤1.1 從量子力學(xué)的薛定諤到計(jì)算平臺(tái)Schr?dinger 這個(gè)名字源于物理學(xué)家 Erwin Schr?dinger他提出的薛定諤方程描述了微觀粒子的波函數(shù)演化?，F(xiàn)代分子模擬的底層邏輯本質(zhì)上還是解這個(gè)方程——不過(guò)不是解析解而是通過(guò)各種近似方法分子力學(xué)、半經(jīng)驗(yàn)方法、從頭算、DFT在計(jì)算機(jī)上逼近真實(shí)體系的能量和運(yùn)動(dòng)狀態(tài)。Schr?dinger 這家公司成立于1990年核心就是用物理模型去預(yù)測(cè)分子的性質(zhì)和行為把學(xué)術(shù)圈里散落的算法工具包裝成一套工業(yè)級(jí)工作流。很多人第一次接觸“薛定諤”是在藥企實(shí)習(xí)或者計(jì)算化學(xué)課程里打開(kāi) Maestro 界面看到一堆模塊菜單。整套套件里最常見(jiàn)的幾個(gè)模塊是Glide分子對(duì)接、Prime蛋白結(jié)構(gòu)預(yù)測(cè)與 MM-GBSA 重打分、LigPrep配體準(zhǔn)備、Desmond分子動(dòng)力學(xué)模擬、FEP自由能微擾、Canvas化學(xué)信息學(xué)與構(gòu)效關(guān)系分析。它們不是獨(dú)立的小工具而是共享同一套力場(chǎng)參數(shù)、同一套文件格式、同一個(gè)圖形界面這意味著從受體準(zhǔn)備到結(jié)果分析整個(gè)流程都是串起來(lái)的不太需要在不同軟件之間玩“文件搬家”。1.2 誰(shuí)在什么場(chǎng)景下會(huì)用到它我見(jiàn)過(guò)三類人會(huì)真正用到這套軟件。第一類是藥物研發(fā)團(tuán)隊(duì)的計(jì)算化學(xué)家他們拿它做靶點(diǎn)的虛擬篩選、先導(dǎo)化合物優(yōu)化配合濕實(shí)驗(yàn)驗(yàn)證第二類是結(jié)構(gòu)生物學(xué)出身的研究者手里有晶體結(jié)構(gòu)或冷凍電鏡結(jié)構(gòu)想快速看這個(gè)靶點(diǎn)能裝下什么形狀的分子第三類是跟計(jì)算沾邊的學(xué)生和跨界研究者比如做農(nóng)藥設(shè)計(jì)、材料篩選、甚至酶催化機(jī)理的人也會(huì)借用對(duì)接和動(dòng)力學(xué)模擬來(lái)補(bǔ)充實(shí)驗(yàn)解釋。這些場(chǎng)景有一個(gè)共同點(diǎn)核心需求都是“預(yù)測(cè)”。結(jié)晶一個(gè)蛋白可能花掉幾個(gè)月合成一個(gè)化合物又需要時(shí)間和經(jīng)費(fèi)而計(jì)算可以通過(guò)打分函數(shù)先篩一遍把候選范圍從百萬(wàn)級(jí)壓到幾十個(gè)。所以 Schr?dinger 解決的不是“算出真相”而是“用較少資源找到最可能正確的方向”。這點(diǎn)非常重要后面所有參數(shù)選擇、流程設(shè)計(jì)都是圍繞它來(lái)的。2. 用 Schr?dinger 搭一條分子對(duì)接流水線2.1 結(jié)構(gòu)準(zhǔn)備Protein Prep Wizard 和 LigPrep別跳過(guò)這一步我第一次用薛定諤跑對(duì)接的時(shí)候以為拿一個(gè) PDB 結(jié)構(gòu)直接 Glide 就能出結(jié)果結(jié)果對(duì)出來(lái)的打分爛得離譜。后來(lái)才明白PDB 文件里存的是實(shí)驗(yàn)解析的原子坐標(biāo)但這個(gè)坐標(biāo)系里其實(shí)缺了很多化學(xué)上必須的信息氫原子基本都沒(méi)有組氨酸的質(zhì)子化狀態(tài)不明確Asp、Glu 等殘基該帶幾個(gè)質(zhì)子也沒(méi)定活性口袋里可能還混著結(jié)晶緩沖液分子、金屬離子、二聚體界面上的水。Protein Prep Wizard蛋白準(zhǔn)備向?qū)ё龅木褪沁@些“補(bǔ)齊”工作。實(shí)際操作里我會(huì)依次做三件事。第一給蛋白加氫并分配正確的質(zhì)子化狀態(tài)第二用 Prime 補(bǔ)缺失的殘基和側(cè)鏈同時(shí)處理掉不合理的原子位置沖突第三做一輪能量最小化讓整個(gè)結(jié)構(gòu)在力場(chǎng)下達(dá)到一個(gè)穩(wěn)定構(gòu)象。這里要特別提醒最小化的原子位移限制不要設(shè)太大默認(rèn) 0.3 ? 就夠了否則會(huì)把活性位點(diǎn)壓變形后面的對(duì)接全廢。配體端的準(zhǔn)備用 LigPrep。它會(huì)生成配體的 3D 構(gòu)象計(jì)算特定 pH 下的質(zhì)子化狀態(tài)還能產(chǎn)生合理的立體異構(gòu)體和環(huán)構(gòu)象。這一步很多人不重視直接拿 2D 結(jié)構(gòu)的 smiles 去對(duì)接結(jié)果發(fā)現(xiàn) Glide 報(bào)錯(cuò)或者產(chǎn)出一堆高能構(gòu)象。常見(jiàn)做法是設(shè)定 pH 7.0 ± 2.0讓 Epik 模塊按生理?xiàng)l件計(jì)算所有可電離基團(tuán)的質(zhì)子化狀態(tài)再對(duì)每個(gè)輸入結(jié)構(gòu)生成最多 32 個(gè)低能 3D 構(gòu)象。聽(tīng)起來(lái)簡(jiǎn)單但這一步直接決定后續(xù)對(duì)接能不能找到“對(duì)的姿勢(shì)”。提示結(jié)構(gòu)準(zhǔn)備是整個(gè)流程的地基。模板結(jié)構(gòu)里如果帶著共晶配體Prep Wizard 默認(rèn)會(huì)保留它不要盲目刪掉。這個(gè)配體往往是網(wǎng)格生成的定位點(diǎn)也是判斷對(duì)接結(jié)果是否合理的錨定參照。2.2 Glide 分子對(duì)接參數(shù)怎么選Glide 是 Schr?dinger 的招牌對(duì)接工具它把配體放進(jìn)受體結(jié)合口袋搜索所有可能的結(jié)合構(gòu)象然后用 GlideScore 打分。具體用哪種精度取決于你要解決的問(wèn)題。Glide 提供了三檔模式HTVS高通量虛擬篩選、SP標(biāo)準(zhǔn)精度、XP超高精度。HTVS 最粗糙也最快適合動(dòng)輒百萬(wàn)級(jí)的超大庫(kù)初篩SP 是日常最常用的精度兼顧速度和可靠性XP 會(huì)對(duì)配體-受體之間的形狀互補(bǔ)、疏水作用做更精細(xì)的懲罰和獎(jiǎng)勵(lì)但極其耗時(shí)適合對(duì)幾十個(gè)命中做深度分析。網(wǎng)格生成這一步最容易被忽略。Glide 需要先定義一個(gè)受體網(wǎng)格Receptor Grid網(wǎng)格的中心通常放在共晶配體的質(zhì)心上這樣能保證結(jié)合口袋被完整覆蓋。網(wǎng)格盒子的大小不需要貪大默認(rèn)的 15 ? × 15 ? × 15 ? 一般夠用盒子里還可以指定某些殘基為柔性殘基讓側(cè)鏈在對(duì)接時(shí)微調(diào)但每加一個(gè)柔性殘基計(jì)算量會(huì)指數(shù)增加所以一般只對(duì)活性位點(diǎn)里最關(guān)鍵的幾個(gè)殘基開(kāi)啟。如果做一個(gè)小規(guī)模的篩選我的習(xí)慣是先用 SP 跑一遍把打分前 20% 的配體拿出來(lái)再用 XP 重新對(duì)接這兩輪的結(jié)果綜合起來(lái)決定下一步實(shí)驗(yàn)。很多人上來(lái)直接跑 XP一個(gè)配體跑十幾分鐘換來(lái)的是被噪聲淹沒(méi)的分?jǐn)?shù)純屬浪費(fèi)時(shí)間。2.3 打分函數(shù)怎么選、怎么看對(duì)接完成后你會(huì)在 Maestro 的結(jié)果表里看到一列 GlideScore。它表示的是配體結(jié)合到受體上時(shí)的一種“經(jīng)驗(yàn)性評(píng)分”不是真實(shí)的結(jié)合自由能但數(shù)值越負(fù)說(shuō)明預(yù)測(cè)的結(jié)合越強(qiáng)。它由好幾項(xiàng)加在一起靜電相互作用、范德華力、氫鍵、疏水接觸、溶劑化效應(yīng)、以及配體內(nèi)部的應(yīng)變能。理解“打分函數(shù)不是能量”這一點(diǎn)很重要。真實(shí)結(jié)合自由能要考慮熵、去溶劑化、蛋白構(gòu)象變化這些在對(duì)接打分里要么是近似項(xiàng)要么根本沒(méi)有。所以分?jǐn)?shù)差 0.3 分根本不算差距我自己一般以 1 分以上作為篩選閾值。另一點(diǎn)是分?jǐn)?shù)高更負(fù)不代表活性好因?yàn)榇蚍趾瘮?shù)可能被某些“表面友好”的分子誤導(dǎo)比如分子太大、埋在口袋里強(qiáng)行填滿空隙導(dǎo)致范德華項(xiàng)虛高。我見(jiàn)過(guò)不少人直接按 GlideScore 排序取前 50 個(gè)去做活性測(cè)試結(jié)果命中率慘淡。我的做法是先按分?jǐn)?shù)看排名再逐個(gè)目視檢查結(jié)合姿勢(shì)——看配體有沒(méi)有完全離開(kāi)口袋、有沒(méi)有嚴(yán)重的原子碰撞、極性基團(tuán)有沒(méi)有和水或缺電子區(qū)域形成合理相互作用。視覺(jué)檢查能過(guò)濾掉 30% 到 50% 的“分?jǐn)?shù)好看但姿勢(shì)離譜”的候選。3. 進(jìn)階從虛擬篩選到結(jié)合自由能計(jì)算3.1 虛擬篩選流程設(shè)計(jì)與富集率如果你面對(duì)的化合物庫(kù)有幾萬(wàn)甚至上百萬(wàn)個(gè)分子那就不是一個(gè)個(gè)跑 Glide 的事了而是一個(gè)系統(tǒng)流程設(shè)計(jì)問(wèn)題。我習(xí)慣的虛擬篩選管線分四步。第一步是預(yù)過(guò)濾。所有配體先過(guò)一遍物理化學(xué)性質(zhì)過(guò)濾器比如分子量、AlogP、氫鍵供體/受體數(shù)量、可旋轉(zhuǎn)鍵數(shù)參考類藥五規(guī)則。再跑一遍 ADMET 預(yù)測(cè)把有明顯毒性風(fēng)險(xiǎn)或透膜性差的分子提前踢掉。這一步能把庫(kù)縮到原來(lái)的三分之一甚至更少。第二步是構(gòu)象和異構(gòu)體準(zhǔn)備。用 LigPrep 批量生成 3D 結(jié)構(gòu)。注意控制立體異構(gòu)體數(shù)量沒(méi)有手性中心的分子不要生成無(wú)謂的異構(gòu)體否則后面對(duì)接時(shí)間直接翻倍。第三步是分級(jí)對(duì)接。所有化合物先跑 HTVS按分?jǐn)?shù)保留前 20%30%剩下的跑 SPSP 結(jié)果里取前 1000 到 2000 個(gè)跑 XP。每一級(jí)都沒(méi)有必要把上一級(jí)的全部結(jié)果都拿過(guò)來(lái)篩的就是“低分基本不值得細(xì)看”的思路。第四步是重打分和聚類。用 Prime MM-GBSA 對(duì) XP 命中的前 200 個(gè)分子重新打分再對(duì)化學(xué)結(jié)構(gòu)聚類選每一類里分?jǐn)?shù)最高的一兩個(gè)。這樣最終留下的幾十個(gè)化合物在結(jié)構(gòu)和理化性質(zhì)上都有代表性和差異避免 50 個(gè)命中全是同一個(gè)母核的不同尾巴。富集率是評(píng)價(jià)篩選流程優(yōu)劣的關(guān)鍵指標(biāo)。操作上可以準(zhǔn)備一組已知活性化合物和一組誘餌分子混入候選庫(kù)中跑一遍完整流程看看活性化合物是否比誘餌排在更前面。如果活性化合物大部分排在庫(kù)的頭部說(shuō)明流程的富集能力靠譜如果活性化合物排名跟隨機(jī)差不多那即使表面分?jǐn)?shù)好看也說(shuō)明你的受體結(jié)構(gòu)或準(zhǔn)備步驟有問(wèn)題。3.2 Prime MM-GBSA 的批量重打分GlideScore 是快速篩選的機(jī)槍但到了決賽圈你會(huì)需要精度更高一點(diǎn)的評(píng)價(jià)。Prime MM-GBSA 計(jì)算的是配體結(jié)合前后的能量差綜合了分子力學(xué)能量MM、連續(xù)溶劑化模型GBSA里的極性項(xiàng)和非極性項(xiàng)。它的核心優(yōu)勢(shì)在于部分考慮了溶劑化效應(yīng)這是對(duì)接打分里處理得比較粗的部分。實(shí)際操作不復(fù)雜在 Glide 的結(jié)果列表里選中一批配體右鍵選擇 Prime MM-GBSA設(shè)置力場(chǎng)用 OPLS4、溶劑模型用 VSGB其他保持默認(rèn)然后提交任務(wù)。它會(huì)自動(dòng)對(duì)每個(gè)配體做一個(gè)小規(guī)模的構(gòu)象采樣并給出一個(gè) delta G bind 的估計(jì)值。通常可以把 Glide 的排名和 MM-GBSA 的排名做一個(gè)交叉驗(yàn)證兩者都靠前的分子優(yōu)先做實(shí)驗(yàn)。但這個(gè)手段也有明顯的邊界。它仍是一個(gè)“端態(tài)”計(jì)算忽略了配體結(jié)合過(guò)程中的熵效應(yīng)與通路上的中間態(tài)因此它輸出的數(shù)值適合做排序不適合當(dāng)真實(shí)自由能解讀。比如對(duì)比 A、B 兩個(gè)類似物說(shuō)“B 比 A 預(yù)估強(qiáng) 5 kcal/mol”是合理的但說(shuō)“A 的結(jié)合自由能就是 -48.6 kcal/mol”就沒(méi)有物理意義因?yàn)榻^對(duì)數(shù)值受力場(chǎng)參數(shù)和模型簡(jiǎn)化影響太大。在我自己的流程里MM-GBSA 永遠(yuǎn)是“相對(duì)比較”工具絕不單獨(dú)用它拍板一個(gè)化合物能不能進(jìn)合成清單。3.3 FEP 的思路和適用邊界如果項(xiàng)目推進(jìn)到先導(dǎo)化合物優(yōu)化階段常見(jiàn)的問(wèn)題是在某個(gè)核心骨架上我想把這個(gè)苯環(huán)換成吡啶把甲基換成乙基或者在這里加一個(gè)氟原子活性會(huì)變好還是變壞這時(shí) MM-GBSA 的精度已經(jīng)不夠了就需要引入自由能微擾FEP。FEP 的核心思路是用一系列中間態(tài)連接兩個(gè)結(jié)構(gòu)高度相似的配體計(jì)算從一個(gè)配體“演化”到另一個(gè)配體過(guò)程中的自由能差值。因?yàn)閮蓚€(gè)分子結(jié)構(gòu)接近很多誤差會(huì)在求差的過(guò)程中相互抵消所以它預(yù)測(cè)的相對(duì)結(jié)合自由能精度可以做得非常高誤差通常在 1 kcal/mol 以內(nèi)遠(yuǎn)好于 MM-GBSA也更適合指導(dǎo)化學(xué)家往哪個(gè)方向修飾。但 FEP 不是拿來(lái)隨便跑的。它的第一個(gè)硬性要求是配體之間結(jié)構(gòu)足夠相似最好是只在局部官能團(tuán)上有差異如果你拿一個(gè)完全不同的骨架去做 FEP中間態(tài)很難收斂結(jié)果就失去了意義。第二個(gè)要求是計(jì)算資源充足每個(gè)配體需要搭一張合理的圖哪些配體之間建微擾邊每個(gè)邊都要跑分子動(dòng)力學(xué)模擬動(dòng)輒需要 GPU 集群跑幾天。第三個(gè)要求是你得有一個(gè)可靠的共晶結(jié)構(gòu)或?qū)咏Y(jié)構(gòu)作為起點(diǎn)起點(diǎn)蛋白質(zhì)構(gòu)象如果不對(duì)后面一切都白搭。我自己項(xiàng)目里的分工是這樣的初篩用 Glide中篩用 MM-GBSA最后的十來(lái)個(gè)類似物決定具體合成順序時(shí)才上 FEP。這個(gè)三級(jí)遞進(jìn)既保證了效率也把最貴的計(jì)算留給了最重要的決策點(diǎn)不會(huì)出現(xiàn)“GPU 集群跑了一禮拜結(jié)果發(fā)現(xiàn)模擬的體系根本不是活性構(gòu)象”這種慘劇。4. 實(shí)操中我踩過(guò)的坑4.1 別拿原始 PDB 直接對(duì)接我剛開(kāi)始做激酶項(xiàng)目時(shí)直接從 PDB 下載了一個(gè)分辨率不錯(cuò)的晶體結(jié)構(gòu)沒(méi)有用 Prep Wizard 處理只是刪掉了水分子和其他鏈就生成了網(wǎng)格開(kāi)始對(duì)接。結(jié)果篩出來(lái)一個(gè)打分極高的化合物合成出來(lái)卻完全沒(méi)活性。事后回頭檢查才發(fā)現(xiàn)原始結(jié)構(gòu)里 DFG 基序附近有一個(gè)關(guān)鍵殘基的側(cè)鏈密度不完整PDB 里存的坐標(biāo)是扭曲的Prep Wizard 能自動(dòng)檢測(cè)并修正而我跳過(guò)了這一步。正確操作是每次下載 PDB 后先看序列信息和配體信息再用 Protein Prep Wizard 按默認(rèn)流程走一遍重點(diǎn)觀察 Log 窗口里有沒(méi)有“missing residues”或“alternate conformations”提示。如果有缺失長(zhǎng)度超過(guò)五六個(gè)殘基的 loop看它離活性位點(diǎn)遠(yuǎn)不遠(yuǎn)遠(yuǎn)的話直接保持缺口不要強(qiáng)行建模近的話要考慮換一個(gè)更完整的結(jié)構(gòu)因?yàn)槟Ｐ脱a(bǔ)出來(lái)的 loop 不確定性很大對(duì)接結(jié)果可信度低。4.2 質(zhì)子化狀態(tài)和異構(gòu)體的問(wèn)題配體端最容易翻車的是異構(gòu)體數(shù)量失控。有一次我想篩一個(gè)含有兩個(gè)手性中心的母核LigPrep 跑完提醒我生成了 4 個(gè)立體異構(gòu)體各 32 個(gè)構(gòu)象總共 128 個(gè)文件對(duì)接時(shí)間直接翻了 16 倍。后來(lái)我學(xué)會(huì)了在有明確實(shí)驗(yàn)數(shù)據(jù)的前提下用 chiral 標(biāo)記指定只有目標(biāo)手性構(gòu)型參與計(jì)算或者用 LigPrep 的“retain specified stereoisomers”功能把不需要的異構(gòu)體刪掉。另一件常被忽略的事是水分子的處理。有些活性位點(diǎn)里會(huì)有一個(gè)水分子它同時(shí)在蛋白和配體之間搭橋就好像兩個(gè)陌生人之間的傳話人處于關(guān)鍵的位置。準(zhǔn)備受體時(shí)如果無(wú)條件刪除所有水分子某些配體可能就少了一個(gè)重要的極性相互作用打分被低估。我的做法是保留所有太陽(yáng)下看起來(lái)協(xié)調(diào)的水分子與蛋白極性殘基距離在 3.0 ? 左右且不與非極性殘基沖突的生成網(wǎng)格時(shí)再把它們當(dāng)作受體的一部分參與打分。4.3 資源消耗與并行調(diào)優(yōu)Schr?dinger 的任務(wù)調(diào)度邏輯跟普通程序不太一樣很多新手會(huì)誤以為把所有 CPU 核心都填滿就能加速。實(shí)際上單個(gè) Glide 對(duì)接任務(wù)本身對(duì)多核利用有限提速的關(guān)鍵是“同時(shí)跑多個(gè)任務(wù)”。我習(xí)慣在 Linux 服務(wù)器上把配體庫(kù)拆成幾個(gè)子文件每個(gè)子文件單獨(dú)提交一個(gè) Glide 任務(wù)利用批處理腳本把它們并行調(diào)度起來(lái)。Desmond 和 FEP 則剛好相反它們強(qiáng)烈依賴 GPU 加速?zèng)]有 GPU 的話跑一個(gè) 100 ns 的分子動(dòng)力學(xué)模擬可能要等一周而用單張中端 GPU 能縮短到一天以內(nèi)。所以如果你準(zhǔn)備長(zhǎng)期跑薛定諤的自由能計(jì)算別省 GPU如果只是跑對(duì)接和虛擬篩選CPU 就夠用優(yōu)先保證內(nèi)存不要爆掉就行。許可證連接也是一個(gè)隱蔽的坑。Schr?dinger 的 license 是網(wǎng)絡(luò)浮動(dòng)式的如果服務(wù)器網(wǎng)絡(luò)不穩(wěn)定或者 license 服務(wù)沒(méi)啟動(dòng)啟動(dòng) Maestro 或命令行工具時(shí)會(huì)一直卡在“could not checkout license”的報(bào)錯(cuò)。排查時(shí)先看 license 服務(wù)器日志再確認(rèn)本機(jī)環(huán)境變量 SCHROD_LICENSE_FILE 是否指向正確的端口和地址通常都能解決。5. 常見(jiàn)問(wèn)題速查現(xiàn)象可能原因處理方法Prep Wizard 報(bào)錯(cuò)缺少大量殘基PDB 文件本身分辨率差或截?cái)鄵Q用分辨率更高的結(jié)構(gòu)或刪除無(wú)序區(qū)域盡力補(bǔ)全關(guān)鍵殘基Glide 網(wǎng)格生成失敗共晶配體坐標(biāo)不完整或原子類型沖突檢查配體是否含金屬離子或非標(biāo)準(zhǔn)殘基編號(hào)必要時(shí)重新用 LigPrep 處理LigPrep 產(chǎn)出的構(gòu)象太多立體異構(gòu)體數(shù)量未限制指定手性構(gòu)型控制每個(gè)結(jié)構(gòu)的構(gòu)象生成上限對(duì)接分?jǐn)?shù)很好但結(jié)合姿勢(shì)奇怪配體沒(méi)進(jìn)入口袋或原子碰撞嚴(yán)重用 Maestro 檢查受力分布關(guān)掉不合理的柔性殘基設(shè)置MM-GBSA 結(jié)果不收斂配體電荷設(shè)置不合理或結(jié)構(gòu)構(gòu)象過(guò)少重新檢查配體質(zhì)子化狀態(tài)生成更多的初始構(gòu)象FEP 自由能結(jié)果方差很大中間態(tài)不夠密或者模擬時(shí)間不夠增加兩個(gè)相鄰配體間的中間態(tài)數(shù)量延長(zhǎng)平衡時(shí)長(zhǎng)檢查是否選了結(jié)構(gòu)差異過(guò)大的配體對(duì)還有一個(gè)值得單獨(dú)說(shuō)的技巧批量任務(wù)報(bào)錯(cuò)時(shí)優(yōu)先看 job 目錄下的 write.log 文件。Schr?dinger 的圖形界面經(jīng)常把錯(cuò)誤吞掉只顯示一個(gè)紅色的 “Job failed”但具體的失敗原因一定寫在日志里。學(xué)會(huì)從日志的第一條 error 開(kāi)始處理比反復(fù)重跑要有效得多。實(shí)際項(xiàng)目里我第一次排查一個(gè)批量對(duì)接任務(wù)失敗的問(wèn)題就是因?yàn)槟硞€(gè)配體包含未被力場(chǎng)定義的鹵素原子類型日志里一行 “cannot type atom” 瞬間讓我定位到問(wèn)題換了參數(shù)版本后整個(gè)隊(duì)列就順了。另外建議定期備份 Maestro 中構(gòu)建好的自定義網(wǎng)格和結(jié)果視圖。薛定諤的項(xiàng)目文件結(jié)構(gòu)里除了 .mae 主文件還有配套的 .maegz、.log、.grd 等文件。只拷貝主文件回家換個(gè)機(jī)器繼續(xù)看結(jié)果時(shí)經(jīng)常發(fā)現(xiàn)結(jié)構(gòu)缺失就是因?yàn)闆](méi)有把同一目錄下的附屬文件一起打包。把這些文件統(tǒng)一放進(jìn)一個(gè)項(xiàng)目文件夾里算是一個(gè)做了就會(huì)感謝自己的好習(xí)慣。這套軟件不便宜學(xué)習(xí)曲線也不平緩但只要理順了“準(zhǔn)備結(jié)構(gòu)—生成網(wǎng)格—分級(jí)對(duì)接—重打分—自由能計(jì)算”這條主線它就能成為決策效率的強(qiáng)杠桿。我個(gè)人這幾年最深的體會(huì)是薛定諤的價(jià)值不在于幫你跑出更高的分?jǐn)?shù)而在于逼著你想清楚每一步到底在算什么、誤差在哪里、結(jié)論能支撐多大范圍的判斷。想明白這層軟件選哪家、參數(shù)怎么調(diào)反而都是次要問(wèn)題了。