跳轉到

1994 年 CASP 結構預測競賽:盲測擂台的開幕之案#

副標題:當所有人都說自己能預言蛋白質的形狀,一對姊弟決定把考卷鎖進保險箱。 年表定位:第四幕之首,承接 BPTI 分子動力學與 QM/MM 之後,開啟結構預測的實證時代。

案發現場#

時間是 1994 年,地點是美國阿斯洛馬(Asilomar)海邊的會議廳。空氣裡瀰漫著一種奇特的焦慮。

案情背景有三重死結:

第一,序列爆炸而結構停滯。基因定序技術起飛,蛋白質胺基酸序列如雪片般湧入資料庫,能用 X 光晶體學與 NMR 解出的結構卻寥寥無幾。序列與結構之間的鴻溝,宛如滿屋子的指紋,卻沒有一枚對得上嫌犯。

第二,自吹自擂的推理者。八〇年代末,號稱能從序列預測結構的方法層出不窮:同源建模、穿線法(threading)、從頭預測(ab initio),每一派都宣稱自己的案子破了。但所有驗證都是「開卷考試」—— 拿已知結構做回測,參數調到過關為止。旁觀者根本分不清誰是神探、誰是江湖術士。

第三,Levinthal 悖論的陰影。1969 年,Cyrus Levinthal 提出一個令人心寒的計算:如果一條 100 個殘基的蛋白質每個殘基只有 3 種構形,全部搜尋需要 \(10^{47}\) 年,遠超宇宙年齡。但真實蛋白質卻在毫秒到秒級折疊成功。這意味著折疊絕非隨機搜尋,能量地形必有漏斗結構,而計算方法若不懂漏斗,就只能在指數迷宮裡打轉。

就在此時,兩位關鍵人物走進現場:John Moult 與他的姊妹?不,是學術意義上的「Moult 姊弟」—— John Moult 與 Esther Saul?更準確的說法是,主辦者 John Moult、Krzysztof Fidelis 與 Andriy Kryshtafovych。前兩人拍板了一個近乎無禮的提議:來一場盲測。

所謂盲測,就是預測者只能拿到胺基酸序列,而目標結構的實驗座標在評比結束前嚴格保密。實驗學家先把尚未發表的結構交給評審保管,計算學家限時繳卷,再由獨立裁判打分。這在當時是離經叛道的想法:等於要求所有名偵探蒙眼斷案。

偵查過程(含數學式/表格/理論)#

推理一:為何必須盲測?#

Moult 的推理像法醫重建現場:回測的污染無所不在。研究者無意間見過同源結構、調過參數、挑過好看的例子發表。這叫「數據洩漏」。唯一的解藥是時間之鎖:用尚未公開的結構當考題,讓作弊在物理上不可能。

於是 CASP(Critical Assessment of protein Structure Prediction)第一屆於 1994 年登場。流程堪稱嚴密的刑事鑑識程序:

階段 時間 動作 偵探比喻
徵集靶點 春季 實驗學家提交未發表結構序列 收集未曝光的命案卷宗
發布序列 夏季 公開靶點序列,隱藏結構 只給線索,不給答案
限時預測 數週至數月 各組提交三維座標模型 偵探限時繳交推理書
獨立評估 秋季 比對預測與真實結構 開棺驗屍
Asilomar 會議 12 月 公開放榜,檢討方法 偵探大會對質

推理二:如何打分?GDT_TS 的誕生背景#

最早的評分直覺是均方根偏差。給定預測結構與真實結構各有 \(N\) 個 Cα 原子座標,分別記為 \(r_i^{pred}\) 與 \(r_i^{exp}\) ,疊合後計算 \(RMSD\) :

\[ RMSD=\sqrt{\frac{1}{N}\sum_{i=1}^{N}|r_i^{pred}-r_i^{exp}|^2} \]

但辦案老手很快發現 \(RMSD\) 有致命弱點:它對局部大偏差極敏感。一條尾巴擺錯方向,整體分數就崩盤,無法反映「主體正確、局部走偏」的真實辦案品質。

因此 CASP 後來發展出 GDT(Global Distance Test)及其總分 \(GDT\_TS\) 。概念是:把預測結構往真實結構上疊合,統計落在不同距離截斷內的殘基比例,再取平均。若記 \(P_d\) 為 Cα 距離小於 \(d\) 埃的殘基百分比,則有:

\[ GDT\_TS=\frac{P_1+P_2+P_4+P_8}{4} \]

其中下標數字代表 1、2、4、8 埃四檔截斷。 \(GDT\_TS\) 介於 0 到 100 之間,越高越好。它像寬容的陪審團:主體正確即給高分,不因一處瑕疵判死刑。

推理三:TBM 與 FM 的分類學#

隨著競賽累積卷宗,評審發現不能把所有案子混為一談,必須按「線索多寡」分案:

類別 全稱 線索狀況 典型方法 難度
TBM Template-Based Modeling 資料庫中有同源模板 同源建模、穿線法 較易,如認親案
FM Free Modeling 無可用模板 從頭折疊、物理勢能 極難,如無頭案
TBM/FM 邊界 遠源同源 模板微弱 綜合策略 懸案區

這個分類至關重要。它承認一個殘酷事實:有模板的案子靠「認人」,沒模板的案子靠「推理」。CASP1 的結果顯示,TBM 已初具人形,FM 則慘不忍睹—— 最好的 FM 模型 \(GDT\_TS\) 往往不到 30 分,幾乎等同亂猜。這份成績單像一盆冷水,澆醒了整個領域。

推理四:Levinthal 悖論的理論回響#

CASP 的慘淡 FM 成績,反過來印證了 Levinthal 悖論的理論重量。若折疊真是隨機搜尋,計算機理應束手無策;但蛋白質秒級折疊證明能量漏斗存在。偵探的結論是:預測方法必須內建漏斗物理—— 疏水塌縮、局部片段偏好、協同成核,缺一不可。這為五年後的 Rosetta 片段組裝埋下了伏筆。

結案報告#

CASP1 於 1994 年底放榜。沒有全壘打,沒有神探橫空出世。但它完成了三件比破案更重要的事:

其一,建立了盲測的司法制度。從此結構預測有了獨立的法院,任何方法都必須接受同一批密卷的考驗。弄虛作假的空間被壓縮到最小。

其二,摸清了領域的底線。TBM 可用、FM 不行,這個診斷讓後續十年資源投向正確的戰場:片段庫、知識勢能、採樣算法。

其三,啟動了兩年一屆的連載探案。CASP2(1996)、CASP3(1998)一路延續,直到 2018 年 AlphaFold 登場、2020 年 AlphaFold 2 封神。這部連載的起點,正是 1994 年阿斯洛馬海邊的那間會議廳。

用一句話結案:Moult 姊弟沒有抓住兇手,但他們蓋了一座法院,讓真正的兇手—— 蛋白質折疊之謎—— 無所遁形。

證據與工具#

  • 關鍵公式一: \(RMSD\) 疊合偏差,用於早期評估,公式見偵查過程。
  • 關鍵公式二: \(GDT\_TS\) 四截斷平均分,CASP 的標準量刑尺,公式見偵查過程。
  • 關鍵分類表:TBM 對 FM,見上表;辦案前先問有無模板。
  • 關鍵理論:Levinthal 悖論,搜尋空間隨鏈長指數爆炸, \(3^{100}\) 量級的迷宮。
  • 工具與數據:PDB(蛋白質資料庫)提供靶點來源;DALI、TM-score、 \(GDT\_TS\) 為後續評分工具箱。
  • 辦案心法:盲測即正義。任何未經盲測的準確率宣稱,都應視為嫌疑證詞。
  • 延伸卷宗:下一案是 1999-Rosetta蛋白質預測.md,看 Baker 如何用片段組裝回應 FM 慘案。