動力學與自由能：當單一姿勢不夠用時

蛋白不是雕像

晶體結構凍結了一種構象，但真實的蛋白在不停抖動。分子動力學（MD）模擬這種運動：它把蛋白、配體和顯式水放進一個盒子，施加描述所有原子間作用力的力場，並以極小的時間步長向前推進，觀察體系的演化。MD 揭示出單張靜態圖所掩蓋的構象柔性——擺開的側鏈、把守口袋的環、以及剛性對接看不到的誘導契合。

MD 也是詢問水分子的誠實方式。通過觀察水分子停留在哪裡、被結合得有多緊，你可以辨認出值得作為目標的可被取代的水，以及必須保留的結構水。而且由於 MD 報告漲落，它直接關乎結合熵——結合中常被忽視的那一半，它來自配體鎖定時柔性與無序的變化。

把結合視為自由能平衡

親和力由結合自由能 ΔG 決定，它綜合了焓（你畫出的交互作用）和熵（有序度、柔性和水）。評分函數只粗略地估計 ΔG，這正是它排序差的原因。要把活性預測到足以指導合成的程度，你需要一種紮根於真實統計熱力學、並計入蛋白運動和顯式溶劑的方法。

自由能微擾（FEP）就是這種方法。FEP 不去計算單個分子的絕對 ΔG（那非常困難），而是計算兩個相近類似物之間結合自由能的差值：在口袋中和在水中，全程用 MD，把一個逐漸「變形」為另一個。由於兩條路徑共享大部分結構，誤差相互抵消，現代 FEP 能把相對活性預測到約 1 kcal/mol——足以為下一步該做哪個類似物排序。

Why a difference is easier than an absolute:

  ligand_A  --(alchemical morph in WATER)-->  ligand_B    : dG_water
  ligand_A  --(alchemical morph in POCKET)--> ligand_B    : dG_bound

  ddG_binding(A->B) = dG_bound - dG_water

A simple edit -- swap H for F, add a methyl -- changes little,
so the shared parts cancel and only the change is 'paid for'.

FEP 通過一個熱力學循環計算差值（ΔΔG），於是大部分困難的物理量相互抵消。

用好 FEP，不浪費它

FEP 準確但昂貴，而且它的優劣完全取決於你起始的姿勢。它在正確的地方才有回報——對一個已知系列進行後期、細緻的優化——而如果你讓它去比較差異極大的分子，或依賴一個你並不確定的姿勢，就會浪費算力。把它當作一道精密過濾，並通過用你已經測量過的化合物來核對它，保持它的誠實。

在先導優化中使用 FEP，此時你有一個可信的結合姿勢，並在固定骨架上探索小的修改。
讓微擾小且化學上合理——H→F、加一個甲基、單個環替換——以保持循環可靠。
把計算錨定到已測量的類似物上，以便用真實資料核對預測。
對虛擬想法排序，合成排在前面的幾個，把結果反饋回來，再重複——FEP 是設計循環內的一道過濾，而不是實驗室的替代品。