FAST++ Paper Atlas · 1976

最佳剛體旋轉

Wolfgang Kabsch, Acta Cryst. A 32, 922–923 (1976)

從已知點對應出發，逐步理解去中心、最佳正交旋轉、平移與 RMSD，並用可操作的 3D 點雲觀察雜訊和離群點。

teaching-model

透明星圖校準台：配對已知，姿勢未知

這是教學比喻，不是論文情節。兩張透明星圖上的星星已經一一編號；任務不是猜誰對誰，而是移動其中一張，讓同號星星盡量重合。

  1. 第一章：同號星星早已配好

    校準員收到兩張有相同編號的星圖，因此不用猜 1 號應該對哪一顆；真正未知的是整張圖的姿勢。

    拆回論文語言： Kabsch 問題假設點對應已知，求的是使加權平方誤差最小的正交旋轉。

    Kabsch 1976, p. 922, eqs. (1)–(3).

  2. 第二章：先把兩張圖的中心對齊

    校準員先找出每張星圖的平均中心，把兩張圖都移到同一個工作台原點，暫時消掉整體位移。

    拆回論文語言： 兩組點各自減去質心後，平移與旋轉可以分開處理。

    Kabsch 1976, p. 922.

  3. 第三章：只能轉，不能把透明片拉長

    透明片可以整張轉動，卻不能拉長、縮短或翻成鏡像；校準員尋找總誤差最小的角度。

    拆回論文語言： 旋轉矩陣受正交限制，保持距離與角度；論文推導最佳旋轉並處理秩不足情況。

    Kabsch 1976, pp. 922–923, eqs. (4)–(12).

  4. 第四章：把中心放回正確位置

    轉好之後，校準員再把透明片平移回目標星圖的位置，最後檢查每對星星還差多遠。

    拆回論文語言： 最佳旋轉後以兩個質心恢復平移；RMSD 濃縮已配對點的殘差，但不會判斷配對是否正確。

    Kabsch 1976, pp. 922–923.

Kabsch 的四站白話地圖

只追蹤一件事：在點對應已知時，如何用一次剛體運動把兩組點疊得最好。

  1. 1. 先確認配對是題目給的

    Kabsch 不負責從兩個蛋白質中找出殘基配對；它接手的是已經配好的座標。

    這一步只記住： 先配對，再疊合；這篇只解後半段。

    Kabsch 1976, p. 922.

  2. 2. 各自減去質心

    把每一點減去該組的平均位置，兩組點便都以原點為中心。

    這一步只記住： 去中心是在暫時拿掉平移。

    Kabsch 1976, p. 922.

  3. 3. 找最佳正交旋轉

    在不拉伸幾何的條件下，找出讓所有點對平方距離總和最小的旋轉。

    這一步只記住： 最佳是全體誤差的折衷，不保證每一點都完美。

    Kabsch 1976, pp. 922–923.

  4. 4. 恢復平移，再讀 RMSD

    旋轉完成後用兩個質心算回平移，最後以 RMSD 描述配對點的平均殘差尺度。

    這一步只記住： 低 RMSD 只回答已選點對疊得近不近。

    Kabsch 1976, pp. 922–923.

理解檢查

  1. Kabsch 1976 會自行猜殘基配對嗎？

    • 不會；它假設點對應已知
    • 會，而且不需要座標
    • 只在 RMSD 為 0 時會

    它解的是已知對應下的最佳剛體旋轉，不是完整的結構搜尋。

  2. 為什麼先減去質心？

    • 先把平移與旋轉拆開
    • 把所有點縮成同樣長度
    • 刪除離群點

    去中心讓兩組點共用原點，之後可專心求旋轉。

  3. 低 RMSD 能證明殘基配對一定正確嗎？

    • 不能；它只摘要目前點對的殘差
    • 能，任何低值都證明同源
    • 能，但只限三個點

    短配對或錯誤挑選都可能得到好看的 RMSD，仍要看 aligned length、coverage 與配對合理性。

完成任務： 不用看上文，用四句話說完：輸入是什麼、為何去中心、旋轉不能做什麼、RMSD 沒回答什麼。

paper-fact

問題設定：對應已知，姿勢未知

給定兩組一一對應的向量 xₙ 與 yₙ，以及每一對的權重 wₙ，目標是找一個不拉伸、不縮放的正交矩陣 U，使旋轉後的平方偏差總和最小。平移可先用兩組點的質心消去。

這一步不負責猜測哪個殘基對哪個殘基。若對應本身錯了，最佳旋轉仍只能對錯誤的配對做最小平方折衷。

E(U) = ½ Σₙ wₙ ‖U xₙ − yₙ‖², subject to UᵀU = I

誤差 E 是所有已知點對距離平方的加權總和；U 必須保持長度與角度。

Kabsch 1976, p. 922, eqs. (1)–(3).

paper-fact

四步讀懂方法

第一步，把兩組點各自減去質心，讓平移問題消失。第二步，從配對點形成關係矩陣 R 與對稱矩陣 S。

第三步，以正交限制和拉格朗日乘子求出能使誤差最小的旋轉 U。第四步，把目標質心減去旋轉後的來源質心，取回平移量。

論文也處理秩不足與平面情況；這提醒我們：點若幾乎共線或共面，有些旋轉方向會缺少足夠資訊。

Kabsch 1976, pp. 922–923, eqs. (4)–(12).

project-reading

如何讀 RMSD

RMSD 是最佳疊合後殘差的平方平均根。它把整體誤差濃縮成一個長度單位，但不會告訴你誤差集中在哪一點，也不會自動判斷配對是否合理。

少數離群點可能把整個最小平方解拉向自己。互動實驗的離群點預設正是用來看這個限制。

project-reading

它在後續論文中的位置

CE、FAST 與 TM-align 都必須在某個階段把已選出的殘基對疊合；Kabsch 問題是其中的幾何子問題。DALI 則先用不受剛體旋轉影響的內部距離矩陣搜尋結構關係。

小例子：被搬動的三角形

A 有三個點 (0,0,0)、(2,0,0)、(0,1,0)。B 是同一三角形先旋轉 90°，再平移到 (4,−2,1) 附近。三個對應已知。

  1. 分別計算 A、B 質心，並從每個點扣除質心。
  2. 由中心化點建立關係矩陣，求保持長度與角度的最佳旋轉。
  3. 將旋轉套回 B，再用兩個質心恢復平移。

沒有雜訊時，三對點可重合，疊合後 RMSD 接近 0；加入雜訊或離群點後，結果變成所有點的最小平方折衷。

名詞表

Centroid／質心
一組點各座標的平均位置。
Orthogonal rotation／正交旋轉
保持距離與角度、且不產生鏡射的剛體旋轉。
RMSD
已配對點距離平方的平均值再開根號。

互動實驗

啟用 JavaScript 後可操作互動實驗。

來源與限制

  • 不會自行尋找殘基對應，也不能作蛋白質搜尋。
  • 最小平方解會受離群點影響；單一剛體旋轉也不能描述柔性或多域運動。
  • 點雲是合成資料，不是 PDB 生物結論。