1977 年 BPTI 蛋白質 MD：蛋白會動的目擊證詞#

報案人：靜態晶體結構。案情：X 光照片裡的蛋白像標本，沒人見過它動。偵探 McCammon、Gelin、Karplus 讓 BPTI 跑了 8.8 皮秒，第一次看見蛋白呼吸。

案發現場#

1977 年的結構生物學是一間停屍間。

X 光給出漂亮的蛋白質結構，BPTI 僅 58 個殘基，緊湊明確，堪稱模範被害人。但照片全是凍結平均， B 因子只敢暗示晃動，無人知原子級動態真相。

謎團有三：酶如何找底物，抗體環為何柔軟，別構如何遠距傳遞，靜態結構一概沉默。另有 Levinthal 佯謬，構象空間天文數字，蛋白卻毫秒折疊，地形必有玄機。

McCammon 等人鋌而走險：把 BPTI 放進真空近似，用 Verlet 硬跑數千步，看蛋白會不會散架。 8.8 皮秒今天連眨眼都不夠，當時卻是創世。

偵查從力場驗屍開始，勢能寫成可加形式：

\[ U = U_{bond} + U_{angle} + U_{dihedral} + U_{vdW} + U_{elec} \]

其中 \(U_{bond}\) 為鍵伸縮， \(U_{angle}\) 為角彎曲， \(U_{dihedral}\) 為扭角， \(U_{vdW}\) 為凡得瓦， \(U_{elec}\) 為庫侖靜電。

牛頓方程由 Verlet 推進：

\[ m_i \ddot r_i = -\nabla_i U \]

其中 \(m_i\) 為原子質量， \(r_i\) 為原子座標， \(U\) 為上述力場。

現場筆錄戲劇性：BPTI 沒有散架，偏差維持 1 至 2 埃，穩定卻處處顫動。側鏈翻轉、氫鍵開合、表面大幅擺動、內部剛硬。

均方位移與晶體 B 因子趨勢一致，抖動不是雜訊而是真相。 8.8 皮秒只見熱振動，不見折疊，偵探坦承此限，卻預言更長軌道將揭示功能運動。此程式後來長成 CHARMM。

1977 年《Nature》短文如驚雷，蛋白會動從猜測變證詞， MD 殺入生物學。遺產有三：動態結構觀確立，功能藏在漲落裡；CHARMM 與 AMBER 時代開啟；時間尺度軍備賽起跑，從皮秒邁向毫秒，起點正是這 8.8 皮秒。

BPTI 用顫抖證明：靜態即是假象。

關鍵證物一：BPTI 檔案。

關鍵證物二：MD 啟動虛擬碼。

# r 為座標，v 為速度，U 為力場函數
v = init_maxwell(T=300)
for step in range(9000):
    f = -grad(U, r)
    r, v = verlet_step(r, v, f, h=1e-15)

延伸閱讀伏筆：Nosé-Hoover 將為 MD 補上溫度控制拼圖。