教材：讀懂 Nico 的 Cas9 變體篩選漏斗（v5 / v6）

0 一頁總覽

Nico 這台漏斗，本質是候選削減機，不是「物理真理機」。SpCas9（1368 殘基）的多點突變庫組合爆炸（10⁴–10⁵⁺），濕實驗一輪只做得起 ~10–50 個，所以用便宜、會看核酸、可以有雜訊的篩子在頂層把庫富集，把貴的物理留給少數存活者。

層	處理量	在算什麼	你的類比
Phase 1 序列	10⁴–10⁵	序列合不合理（幾何/親和/演化/摺疊/溶解）	便宜 descriptor 初篩
Phase 2 結構	100–500 摺疊	疊出 all-atom 複合體 + 清座標	建結構 + minimize
Phase 3 打分	10–100	結合能 ΔG_bind（三條路）	隱溶劑單點 / AIMD 系綜
Phase 4 FEP	1–5	精確相對 ΔΔG_bind	RPA / CCSD(T) 級把關

貫穿全文的重點：漏斗每一個定量分數都是「平衡結合自由能 ΔG_bind」（外加摺疊穩定 ΔΔG_fold）。但 Cas9 的專一性由「動力學」決定（HNH 構形檢查點 / kinetic proofreading），而動力學是態函數 ΔG 數學上碰不到的一維。→ 漏斗是很棒的「結合可行性富集器」，但不是專一性引擎。

1 名詞速查

PAM：Cas9 動手前先找的短 DNA 標記（NGG）。找到才開始解 DNA。

R-loop：sgRNA 與 target DNA 從 PAM 端像拉鏈配對、頂開另一股，形成的 RNA:DNA 雜合。

RNP：蛋白+核酸複合體。Cas9 = 蛋白+sgRNA+DNA(+Mg²⁺)，缺核酸就沒生物學。

Inverse folding：AF 問「序列→結構」；反向摺疊問「給定固定骨架，哪些序列會摺成它」。

ΔΔG_fold vs ΔΔG_bind（最重要）：前者=蛋白自己還折不折得起來；後者=複合體抓不抓得牢。是熱力學循環不同邊、可號相反。

pLDDT / PAE：摺疊器對自己座標的信心（learned error bars），不是實驗準度、不是能量。

隱溶劑 / PB：把顯性水換成介電連續體+離子屏蔽近似溶劑靜電（=VASPsol/COSMO 那套）。

介電 override（pdie 4–10）：MM-PBSA 裡溶質內部介電（非溶劑 80），蛋白預設 ~1–2，對帶電骨架調高。是旋鈕不是常數。

Kinetic proofreading：把不可逆的切割藏在慢的構形檢查點後，讓錯基質先脫離。專一性 ∝ k_cleave/k_off，不是 ∝ 1/K_d。

I → D state：HNH 核酸酶域要從中間態 docking 進催化活性態才會切。

2 核心問題與設計哲學

為什麼需要漏斗：20 個位點就有 ~10⁵ 雙突變、~10⁷ 三突變，濕實驗做不起百萬個。漏斗唯一的任務：從天文數字的庫裡便宜地挑出「值得真的合成」的 ~10–50 個。

「Information Economy」每往下一層，每個變體更貴、但變體更少，總成本 =（便宜×很多）+（昂貴×很少）有界。一個把 90% 明顯垃圾便宜刷掉的雜訊篩子，就算不準也省下巨量下游成本。

階層篩選漏斗，從很多候選收斂到少數 — 階層漏斗：頂層便宜工具處理 10⁵ 個變體做富集，底層昂貴物理只跑少數存活者——與你熟的多層 DFT 篩選（PBE 掃全部、HSE/GW 只留給存活者）同構。

三個設計直覺（讀懂 Nico 每個選擇的鑰匙）：

便宜雜訊在前、昂貴精準在後：早期誤殺好變體便宜可容忍；到濕實驗的 false positive 才貴。
正交多票否決：每個閘門攔不同失敗模式，要全過才前進——不是相加成一個總能量。
富集不是真理：頂層工具是「提高存活池好貨比例」，不是最終排名。用「有沒有 MD 準」評它就搞錯重點。

3 Phase 1：序列級篩子（10⁴–10⁵ → 100–1000）

全部「只吃序列、不折疊、一次前向推論」。六個正交閘門：

Gate 1a — LigandMPNN 主篩

做什麼：把整個 Cas9+核酸骨架固定，逐殘基問「這位置放這胺基酸合不合理」。關鍵：它的圖包含 DNA/RNA/離子，會懲罰撞到 guide 或 target 的序列。

為什麼 Nico 用：Cas9 是 RNP，最要命的殘基都朝著核酸。v3/v4 的通用 ESM-IF1 對核酸是盲的。換成 LigandMPNN 是 v5 最重要的工具決策。

物理類比：對固定 host lattice（含 adsorbate 原子）打 site-occupancy 分數。ProteinMPNN 像把核酸從 supercell 刪掉再算。

⚠ 它給的是 pseudo-energy / log-likelihood，沒校準到 kcal/mol，只能排序、不能當 ΔG 加總。

Gate 1b — DeePNAP

做什麼：ProNAB 訓練的深度迴歸器，不建結構直接從序列預測結合的變化量 ΔΔG_bind / ΔK_d（是變化，不是絕對 K_d）。CPU 秒級。

為什麼：唯一直接指名「結合」目標量的、又便宜，當獨立一票。

物理類比：ML surrogate / MLP 取代昂貴能量評估——像用迴歸從組成猜吸附能。

⚠ domain-of-applicability：多點 Cas9 突變很可能遠離 ProNAB 分布，外推像把 MLP 用在沒 fit 過的組態。

Gate 1c — ESM-2 + thermoMPNN v6 新增 + CamSol + EvoEF2

ESM-2：蛋白語言模型，打「演化語法合理性」（只看序列，是 prior 不是 measurement）。

thermoMPNN（v6 新增）：Megascale 27 萬筆訓練，預測 ΔΔG_fold（摺疊穩定度）。ThermoMPNN-D 做雙突變。

CamSol：溶解度/聚集；EvoEF2：經驗力場式 scaffold 穩定度（thermoMPNN 的非 ML、可解釋兄弟）。

為什麼 v6 加 thermoMPNN：補漏洞——一個突變可以「紙上結合漂亮但根本折不起來」。把「會不會摺」和「結不結合」兩種失敗模式分開。

⚠ 鐵則：ΔΔG_fold ≠ ΔΔG_bind。thermoMPNN 給的是「蛋白自己穩不穩」，不是「抓不抓得牢」。

4 Phase 2：結構級（100–1000 → 100–500 摺疊）

為什麼需要 3D：Phase 1 全是「序列進、分數出」，看不到「這顆突變側鏈實際擺哪、跟哪個 DNA 磷酸撞到」。類比：光有化學式 Fe₃Ni 排不出能量，得先有 POSCAR。

Gate 2a — AlphaFold3 / Chai-1 / Boltz-2（會看核酸的摺疊器）

做什麼：輸入蛋白+sgRNA+DNA+Mg²⁺，吐整個 RNP 的 all-atom 座標 + 信心圖（pLDDT/PAE）。~1–3 min/複合體。

為什麼：AF3 介面最準當 primary；但 AF3 權重非商用、叢集部署有法律摩擦，故備開源退路 Chai-1 / Boltz-2。「primary + 可互換」不被單一授權綁死。

物理類比：學出來的結構生成器（對應你「從頭建初始結構」）。pLDDT/PAE 像收斂/誤差診斷。

Gate 2b — 鬆弛作為「選擇」：FoldX vs Mini-MD v6 改動

做什麼：ML 座標是「生的」，有 clash，直接算連續介質靜電會爆炸。v6 二選一：FoldX（經驗、快轉側鏈、CPU 秒級）或 Mini-MD（OpenMM/GROMACS 最小化+短 NVT，顯式 ff14SB/OL15/OL3 真物理鬆弛）。

物理類比（你的主場）：Mini-MD 就是你跑生產 MD 前的 minimize+equilibrate；跳過就像對重疊原子做單點，SCF 不收斂。用 Mini-MD 還能與後面打分共用力場、能量面自洽。

5 Phase 3：打分（10–100）——v6 把物理升為主線

Path A — 靜態 MM-PBSA v6 primary

做什麼：用一張鬆弛快照估結合能：ΔG_bind ≈ E_MM(Coulomb+vdW) + G_solv(PB)。能抓 ML 抹平的斷裂鹽橋。

為什麼 v6 升 primary：最響亮的訊號——把可解釋、能歸因到殘基的物理能量分解推到預設，把黑盒 DL 降為備選。

物理類比：力場版隱溶劑結合能 = DFT + implicit solvation（VASPsol）算吸附能。單快照的問題跟你算吸附能一樣。

⚠ 對 dense protein-RNA ranking 相關性天花板約 r≈0.56；是「top-20% 粗篩」不是精排。

Path B — PDIScore + AlphaBridge（DL 幾何打分）v6 alternate

做什麼：完全不算物理。PDIScore=graph 讀「原子–核苷酸距離」；AlphaBridge=直接從 AF3 的 PAE/pLDDT 抽介面剛性。極快、可吃 unrelaxed。

物理類比：材料界的 ML descriptor（像用 d-band center 迴歸先排 HEA）。

⚠ PDIScore 論文自報 ranking PCC≈0.50——能刷爛貨、抓大趨勢，不能分辨兩個接近的候選誰較好。

Path C — 系綜 MD + MM-PBSA premium

做什麼：跑 20–50 ns 顯性溶劑 MD，沿途多快照做 MM-PBSA 再平均。對帶電骨架把內部介電 override 到 4–10。

物理類比：單快照 vs 系綜 = 0 K 靜態能量 vs 有限溫度自由能。Path C 像跑 AIMD/TI 取系綜平均，把熵與構型漲落納入。

⚠ 介電 override 一定要做 pdie 2/4/8 敏感度掃描報成誤差棒（等同收斂測試）。

6 Phase 4：Alchemical FEP（1–5）——唯一的精確物理

概念（你會很有共鳴）：沿一個非物理的 λ 座標把 WT 側鏈煉金術般 morph 成 mutant，在複合體與 apo 各做一次，靠熱力學循環讓算不出的邊抵消，得到精確 ΔΔG_bind = ΔG_mut(複合體) − ΔG_mut(apo)。pmx 建 hybrid 拓撲，MBAR 把所有 λ 窗最優合成一個數+誤差（=WHAM 的 binless 繼承者）。

物理類比：就是你認得的熱力學循環（Hess/Born-Haber）+ TI / adiabatic switching。

⚠ 一個要更正的常見誤解（很重要） v6 明確聲稱它能處理會改變淨電荷的胺基酸突變（例如 Ala→Asp），做法是同時做 co-alchemical 反離子變換維持電中性——也就是說 Nico 已經想到 FEP 最容易出錯的那個 case，不是忽略它。
所以對他 push 的不是「你沒做 charge-changing」（他有），而是：「counter-ion co-transformation 對 Cas9 這種中和骨架接觸鹼性殘基的突變，收斂/可靠度到底如何？值得 stress-test。」——更精準、也更站得住。

⚠ 範圍是蛋白胺基酸突變（不做 DNA/RNA 鹼基 alchemy）。而且它是 end-state 平衡量：再精確，對速率、障礙、構形檢查點一無所知。

7 v3/v4 → v5 → v6：讀出 Nico 的意圖

版本	動作	透露的意圖
v3/v4	通用工具（ESM-IF1）+ 有 off-target specificity MD tier	泛用 + 有嘗試碰專一性
v5	換核酸感知工具、模組化打分、AF3 primary；拿掉 specificity tier；DL 推論當 primary	「現代化、快、模組化」——但唯一碰專一性的層消失了
v6	加 thermoMPNN；Gate 2b 做成 FoldX/Mini-MD 選擇；把靜態 MM-PBSA 升為 primary	預設越來越物理——把物理邊界往漏斗上方推

讀出的核心Nico 是個熱力學腦、會自我修正的設計者——ML 在頂層負責通量，物理在底層負責信任。唯一一直沒回來的，是「專一性 / 判別」這一維：這不是疏忽，是這台「熱力學優先」漏斗依它自己一致的邏輯，留下未量的那一軸。（v6 的「geometric specificity」只是目標宣示、沒有對應模組。）

8 最關鍵的一維：kinetics（構形檢查點）

Cas9 不是「結合=切」的開關，是多步機器。辨識 on/off-target 主要發生在第 3 步（構形檢查點），不是結合那步：

辨識 on/off-target 主要發生在第 3 步（構形檢查點），不是第 1–2 步（結合）。

HNH 域 on-target docking 切割 vs off-target 卡住脫離 — 左：完全配對時 HNH 域 docking 進活性態、切割（綠）。右：有 mismatch 時同一個域卡在抬起的不活性態、切不下去而脫離（紅）。高保真變體對 on/off-target 的**結合都跟野生型差不多**，靠的正是這一步的差別（Chen 2017；Dagdas 2017）。

一句話講死：ΔG_bind 是態函數，只看 end-state「已結合 vs 未結合」的差，數學上不含任何「路徑上的障礙高度」。所以無論 FEP 收斂到多化學精度，都變不出 HNH docking 的活化能障 ΔG^‡——這是結構性的不可能，不是精度不夠。就像你無法從 ΔG_reaction 反推速率。

深井（結合）與高障礙（動力學）的能量地景 — funnel 算的是井深（ΔG_bind）；Cas9 成敗由**障礙高度**（ΔG^‡，HNH docking）決定。off-target 是被困在深井裡的 kinetic trap——穩定結合＝深井，但跨不過障礙去切。MM-PBSA 算井深，不算障礙。

為什麼這麼容易被漏（對任何人都一樣） ① 可觀測性偏誤：結合看得見摸得著，直覺「抓得牢=作用強」。② 工具成熟度落差：平衡結合工具鏈幾十年成熟、可高通量；kinetic barrier 沒有現成高通量工具。③ 歷史包袱：v3/v4 那個 specificity tier 被 v5 拿掉，而它剛好是唯一碰專一性的層。
就像早期只用 ΔG_ads 當 HER descriptor（好算、可掃），但真正限速常是 barrier——NEB 又貴又難自動化。Cas9 funnel 走的是一模一樣的路。

★ 我們實際跑了什麼（把漏斗上游接上真硬體）

前面都在講「該怎麼讀這台漏斗」。這一節是進度：漏斗上游的幾個閘門，我們已經在國網 Nano5 GPU 上實際跑起來了。更關鍵的是——Gate 2a 摺出來的結構，剛好給了上一節 kinetics 論點一個我們自己算出來的證據。

你在國網跑的真軌跡 — WT SpCas9 的 2 ns MD（Cα trace，40 幀，H100）。顏色＝該殘基的 RMSF（藍=剛硬 → 紅=最會動）。紅色那一團就是 HNH 檢查點域，明顯比其餘骨架更會擺動。可拖曳旋轉；apo 蛋白，動作幅度是上限。

Gate 1c — ESM-2（演化適應度）已跑 ✓

算了什麼：14 個變體面板（催化殘基 + eSpCas9 + HF1 + HypaCas9）逐一打演化分，外加 REC3 區（660–700）整段 site-saturation 掃描（779 變體）+ 熱圖。

結果：催化殘基最保守（−10 ~ −13，動不得，符合預期）；eSpCas9 的 K1003A/R1060A、HypaCas9 的 H698A 最可容忍（−4 ~ −6）→ 演化上不介意，是好候選。本機 CPU 與 H100 叢集逐位元完全重現。

物理類比：先用便宜 descriptor 把明顯不合理的組成刷掉——是 prior，不是 measurement。

Gate 1a — LigandMPNN（核酸感知的序列打分）已跑 ✓

算了什麼：把整個 RNP（蛋白＋sgRNA＋DNA）餵進去，對蛋白逐位打分，graph 裡含核酸當 context——這正是 ESM-2 看不到的那一塊。跟 ESM-2 同一批變體對照。

結果：兩者都抓到催化殘基該罰（sanity check 過）；工程專一性變體被判「結構上可容忍」（同 thermoMPNN）。誠實修正：我原本猜「它看得到 DNA → 會更罰 nt-groove 突變」——沒成立，因為它量的是「序列對骨架的契合」，不是「DNA 結合強度」；拿掉一個接觸 DNA 的 Lys 不會讓蛋白 fold 變差。

⚠ 印證教材前面說的：LigandMPNN 的加值是「不會設計出撞到 guide/target 的序列」，不是直接量結合。

Gate 2a — Chai-1（摺 WT RNP）已跑 ✓

算了什麼：把完整 WT SpCas9 + sgRNA + target/nontarget DNA（1526 tokens）一起摺成 all-atom 結構，H200 上跑，5 個 model、零 inter-chain clash，逐域做 pLDDT 品管。

結果（兩極）：scaffold 很硬 —— REC lobe / RuvC 核心 pLDDT 78–85、sgRNA 75；但 HNH 域只有 63、target DNA 只有 64，最低信心片段是 HNH 的 L1 linker（765–809，pLDDT 46）。

這就是給 kinetics 論點的證據：摺疊器最沒把握的地方（HNH 域 + L1/L2 linker + target DNA），正好就是上一節那個構形檢查點的機械。換句話說——一張靜態摺出來的結構，本身就解不出 HNH 檢查點該長怎樣。這不是模型爛，是「檢查點是動力學、不是單一構型」的直接體現：你要補的那一維，連摺疊器自己都告訴你它不確定。

Gate 1c — thermoMPNN（ΔΔG_fold）已跑 ✓

結果：全 27,360 條突變掃完。所有文獻高保真變體（eSpCas9 / HF1 / HypaCas9）的 ΔΔG_fold 都 < ~1 kcal/mol——全部 fold-neutral、都通過穩定度閘（本來就該如此,它們是真實可用的變體）。最去穩定的高信心突變 D499P/V/K（+4.6）也合理。

一個呼應主軸的點：thermoMPNN 無法區分這幾個變體（全擠在 0~+1）——穩定度跟結合 ΔG 一樣,不是判別專一性的那一維。它正確地讓它們全過閘,真正的差別在動力學。

⚠ 鐵則不變：這給的是 ΔΔG_fold（會不會摺），不是 ΔΔG_bind、更不是專一性。

Gate 2b + 短 MD — OpenMM（2 ns, H100）已跑 ✓

結果：WT 蛋白 scaffold 顯式溶劑 MD（ff14SB/TIP3P,235k atoms,2 ns）。HNH 是最會動的域（Cα RMSF 1.86 Å,其餘域 0.8–1.3）；逐殘基 RMSF 與 pLDDT 反相關（Spearman −0.36）——越低信心的地方越軟。

第三條獨立證據：ESM-2 說 HNH 演化上可容忍、Chai 說它結構信心低、MD 說它最會動——三個獨立方法指向同一件事：HNH 檢查點區域本質上是動態的,正是 kinetics 那一維的機械。

⚠ 誠實 caveat：這是 apo 蛋白（沒放核酸）,完整 RNP 會束縛部分運動,絕對 RMSF 是上限;但相對模式（HNH 最柔）與文獻一致。完整 RNP 建置要 tleap 處理核酸端點,是下一步。

Phase 1 共識篩 — 三個序列閘門合起來（全長 25,992 突變）已跑 ✓

問的問題：教材前面（§2）說漏斗靠「正交多票否決」。那這些閘門真的正交嗎？把 ESM-2 + LigandMPNN + thermoMPNN 在全 1368 位點 × 19 突變（25,992 個）上合起來，實測它們的一致度。

結果：三者中度正相關（Spearman 0.32–0.64）——共享一個「這突變破不破壞」的訊號，但沒到冗餘。兩個結構型的（LigandMPNN／thermoMPNN）最像（0.64）；演化型的 ESM-2 最獨立（0.32–0.40）。→ 「正交」要修正成「互補」：不是完全獨立，但每個都帶不同資訊，多票確實有意義。

方法教訓：一開始只用 14 個變體看，得到 ρ≈0.13（看似完全不相關）——那是小樣本騙人。放大到 25,992 個突變才看到真相 0.4。這正是「統計 finding 在 robustness check 之前別下結論」的活教材。

漏斗吐出的實際名單：三票都放行的可改位點（F1037、H420、S297… 表面／柔性）；三票都否決的不變核心（W883 在 HNH、埋在核心的 Leu→帶電… 結構錨點）。域層次上，REC1／REC3（DNA 辨識葉）最保守、最動不得，PI/CTD 最可改——跟生物學完全一致。這就是漏斗 Phase 1 真正該吐出的東西：一張可改性地圖。

真實數據（不是示意圖）

每殘基 RMSF — **MD 每殘基柔性** — 2 ns 軌跡的 Cα RMSF。HNH 域（橘底）明顯高出其餘域；越柔的地方，正是 Chai 摺疊信心最低處。

RMSF vs pLDDT — **兩個獨立方法收斂** — 每殘基 MD 柔性 vs Chai 結構信心，反相關 **Spearman ρ=−0.36**。物理（會動）與 ML（沒把握）指向同一批殘基＝HNH 檢查點區。

**thermoMPNN 專一性面板** — 所有文獻高保真變體的 ΔΔG_fold 都擠在 0～+1 kcal/mol、全部 fold-neutral。穩定度分不出它們＝穩定度不是判別軸。

ESM-2 REC3 掃描熱圖 — **ESM-2 REC3（660–700）逐位掃描** — 779 個單點變體的演化適應度熱圖。這段正是 HF1/Hypa 專一性突變聚集處。

**全長每域可改性** — 三閘門共識容忍度（全 25,992 突變）。差異其實不大（0.47–0.55），但 **REC1／REC3（DNA 辨識葉）最保守**、PI/CTD 與 bridge helix 最可改——設計時該避開 vs 可動的域一目了然。

相關上線頁：OpenMM GPU MD demo · 為什麼結合自由能漏斗量錯了東西（kinetics 論點）。所有圖為國網 Nano5 實算輸出，非示意。

MD Nico 的 MD benchmark：在信任 MD 之前，先讓它重現「已知答案」

§★ 是「我們已經跑的上游閘門」（apo 蛋白）。這一節是 Nico 剛傳來的兩份文件——一份 07-07 會議整理，一份給你的 MD benchmark protocol（對齊 Skeens 2024）。一句話：在拿 MD 去排「未知的」設計變體之前，先用它重現「已知答案」的系統，確認參數是對的。這正是你熟的——選 functional 前先算已知 lattice constant / formation energy 對標。

A. 核心邏輯：先驗證（benchmark），再設計

不能一開始就用 MD 去排你不知道答案的新變體——你怎麼知道它算得準？protocol 的作法是：先跑一批文獻上答案很清楚的 Cas9 專一性變體（eSpCas9 / HF1 / HypaCas9 / evoCas9），它們都有已發表的 MD ＋實驗結果。如果我們的 MD 能重現這些已知趨勢 → 參數可信 → 才拿去算新的設計變體。

MD 方法驗證的校準迴圈 — protocol 的骨架＝一個**校準迴圈**：先對「已知答案」的專一性變體重現趨勢，過了才相信參數、拿去設計新變體。跟你 DFT 先對已知量校 functional 一模一樣。

B. 兩條軸（最重要的觀念，直接接上 §8 kinetics）

protocol 一開頭就把要量的東西切成兩條軸，而且明說要分清楚：

Cas9 MD 的兩條驗證軸 — 左 **Axis A（HNH 動力學 / fidelity）**：HNH 域大幅、慢速的構形擺動——就是 §8 那個**構形檢查點**。對齊 Skeens 2024 + smFRET，但慢（要 multi-µs）。右 **Axis B（介面 engagement）**：R-loop 穩定度、non-target 股外露、介面接觸——**局部、快**（數百 ns 收斂），是 base / prime editing 設計真正要的軸。

軸	量什麼	快慢	對這專案
Axis A HNH 動力學	HNH 大幅構形擺動、allosteric network（對齊 Skeens 2024 + smFRET）	慢：multi-µs（Skeens ~16 µs/系統）	只當參數校準；你的預算通常到不了
Axis B 介面 engagement	R-loop 穩定度、non-target 股外露、介面接觸	快：~數百 ns 收斂	設計真正要的軸（BE / PE engagement）

接回 §8 的關鍵：Axis A 就是你之前講的 kinetics / HNH 構形檢查點那一維——慢、難、要 µs。Axis B 是介面 / 結合，快、可做。所以 Nico 這份 protocol 其實是務實地承認「kinetics（Axis A）短 MD 抓不到」，先在介面軸（Axis B）把 MD 做扎實。短 MD 重現不了 Axis A 是意料中、可接受的——因為設計本來就靠 Axis B。

C. 要算的系統（benchmark mutation set）——這就是「在算啥」

protocol 的核心是一張對標突變清單：每個系統都有明確功能、已發表的 MD／實驗、和要看的觀測量。FAST＝介面 / 局部，數百 ns 可得；SLOW＝HNH allostery，要 multi-µs。

系統	突變（域）	機制／要驗的（文獻）	快慢
WT	—	基準：複合體穩定度、R-loop 完整、HNH 構形（Palermo 2017 / Skeens 2024 / smFRET）	both
eSpCas9	K848A / K1003A / R1060A（non-target 股溝槽）	中和穩定被頂開 non-target 股的正電溝槽 → 去穩 off-target R-loop（Slaymaker 2016）	FAST · BE 相關
SpCas9-HF1	N497A / R661A / Q695A / Q926A（Rec3＋HNH 旁）	移除對 target 股骨架的氫鍵 → 去掉「多餘結合能」→ 較不容忍 mismatch（Kleinstiver 2016）	FAST
HypaCas9	N692A / M694A / Q695A / H698A（Rec3 α37）	Rec3 / α37 突變 → allosteric 降低 HNH 柔性 → 增強 proofreading（Chen 2017）	SLOW
evoCas9	M495V / Y515N / K526E / R661Q（Rec3）	高專一性但 HNH 動力學 ≈ WT（不同機制）→ 當鑑別對照（Casini 2018）	SLOW
WT + PAM-distal mismatch	DNA 端 protospacer 18–20 錯配（非蛋白突變）	off-target R-loop 去穩 → 你的 Δ(on−off) engagement 軸（Singh 2016）	FAST-ish
dCas9	D10A + H840A（催化死）	結合 / 形 R-loop 如 WT，只是不切 → 穩定度對照	FAST

建議跑序（protocol §6，先做便宜 FAST 的）：① WT（必須先過）→ ② eSpCas9 & HF1（介面 benchmark）→ ③ WT+mismatch（Δ on−off）→ ④ dCas9 / nickase（穩定度對照）→（預算到 µs 才做）⑤ HypaCas9 + evoCas9（重現 Skeens 的 HNH 柔性模式）。

D. Protocol 要點（怎麼算）

起始結構：4UN3（Anders 2014 晶體，HNH 在 checkpoint 態，Skeens 就用這個）。不要用 4OO8（HNH inactive、R-loop 不完整）。→ 要跟廖庭尉拿這個結構＋補 loop。

力場：要跟 Skeens apples-to-apples 用 ff99SBnmr2（NMR-tuned backbone）；一般用 ff19SB 也行。核酸：RNA OL3、DNA bsc1 或 OL15（先從 Skeens SI 確認哪個）；水 TIP3P；離子 Joung–Cheatham；active site 放 Mg²⁺。

平衡：restrained 分段（固定溶質、放水鬆弛 → 慢慢升溫、約束慢慢放）→ production 前 ~200 ns 當額外平衡丟掉（Skeens 對這種大系統的明確 cutoff）。

取樣（你負擔不起 µs 時）：多條短 replica（4–8 × 250–500 ns）而非一條長跑；報 Δ(變異−WT) ＋ replica 誤差棒——差異要大於 replica 散布才可信。

觀測量：穩定度（RMSD、R-loop、RNA:DNA hybrid 配對）；Axis B（介面接觸／氫鍵佔據、non-target 股 RMSF/SASA、per-residue MM-GBSA）；Axis A（HNH ΔRMSF、Jensen–Shannon distance、PCA）。

通過標準＝重現「相對趨勢的方向」，不是絕對數字WT 穩、eSpCas9 non-target 溝槽接觸變弱、HF1 target 骨架接觸變弱、WT+mismatch R-loop 較不穩。方向對就算過。

E. 我們現在正在算什麼

§★ 那條 MD 是 apo 蛋白（沒核酸）。現在照這份 protocol，我已經把 WT 三元複合體（system #1，正對照）建好、在國網 Nano5 上跑起來了：

WT Cas9–sgRNA–DNA 三元 MD 已上機 ✓

建了什麼：從完整 WT RNP（蛋白＋sgRNA＋target/non-target DNA，完整 R-loop）出發，泡水（TIP3P＋0.15 M NaCl）、amber14（ff14SB＋OL15＋OL3）→ 272,791 個原子。解掉了之前卡住的 RNA 5′ 端力場模板問題（把懸空的 5′ 磷酸改成 5′-OH）。本機先建好＋最小化＋跑 300 步驗證穩定，才上機。

怎麼跑：restrained 分段平衡 → NPT 300 K production，4 條 replica 平行（perturbed seeds），可 checkpoint 續跑。

會得到：complex／介面 RMSD、逐殘基 RMSF、介面接觸佔據率 → 先確認 WT 穩（RMSD 進平台、R-loop 沒散），過了才上 mutant 做 Δ(變異−WT)。

⚠ 誠實 caveat：起始用 Chai 摺的結構（不是 4UN3 晶體）、沒放 active-site Mg²⁺、先跑 100 ns 當穩定度先導。正式 Skeens benchmark 要換 4UN3＋補 loop、擴到 250–500 ns × replica。

所以我可以怎麼回 Nico（他問 "any support helpful?"） 你現在看得懂這兩份文件了，可以這樣回：
① 確認理解：「protocol 讀懂了——先用 Skeens 的已知專一性變體對標 MD、分 Axis A(HNH, µs)／Axis B(介面, ns) 兩軸，我先在 Axis B 做。」
② 報進度：「WT 三元已建好、在 Nano5 上跑 4 條 replica（順手解掉 RNA 5′ 端力場模板問題）。」
③ 亮出你的優勢：「Nano5 的 H100 很空、幾乎不排隊——replica／掃參數這種需要大量平行的計算，我這邊跑很快。」（這正是你在合作裡的實質貢獻。）
④ 問幾個關鍵澄清：(a) Skeens SI 裡 DNA 到底用 bsc1 還是 OL15？(b) 第一輪設計要不要先只靠 Axis B（介面）、Axis A 之後再補？(c) 4UN3 需不需要廖庭尉那邊先補好 loop 的版本？

9 怎麼跟 Nico 講（先肯定 → 說好終點 → 共構側枝）

真心肯定 funnel：它把 10⁵→10–50 的濾漏做得又快又物理，v6 往 static MM-PBSA 主線漂移是往更嚴謹的正確方向。這層要留、做得好。
提案 a：pre-register 驗證終點：濕實驗前先講好「拿什麼當成功判準」。若讀的是編輯專一性（on/off cleavage ratio、GUIDE-seq 這類 kinetic/functional 終點），就明確承認 ΔΔG_bind 排名是必要非充分條件。
提案 b：加一條可選 kinetic 側枝：不打掉重來，只跑漏斗底部 ~10–50 個候選：算 HNH I→D docking 的能障 ΔG^‡（NEB/string/metadynamics），或用 MSM 估構形轉換速率。等於用現代工具把 v3/v4 被拿掉的 specificity tier 低成本接回來。

為什麼對方難拒絕：pre-register 是攤開假設的好習慣；側枝是加法、可選、只跑少數候選（成本可控、不動主流程）。而且把對話錨定在你最強的主場——free energy、barrier、MD、MSM。你是來「貢獻缺的那塊物理」的合作者，不是來挑錯的。

附帶：這個「binding vs kinetics 的 scoring-axis」在 AI-Protein-Interface dashboard 裡已被登記成一個待裁決的 open decision（D-006/H-011，待 T-011）——你開這話題，是接一個已經在桌上的問題。

10 可直接照唸的講稿

「你的 funnel 我很欣賞——v5 換成會看核酸的工具、v6 又把靜態 MM-PBSA 升為主打分，方向完全對，它是一台很好的結合可行性富集器，能把 10⁵ 個突變體有效壓到能做的幾十個。

我想補的是它結構上看不到、但恰好是 Cas9 成敗關鍵的一維：專一性其實由動力學決定——off-target 會穩穩結合卻切不下去（HNH 構形檢查點），高保真變體 on/off 的結合都跟 WT 一樣，靠的是結合後的速率步驟。所以 ΔΔG_bind 排名相近的兩個變體，真實專一性可以差很多。態函數 ΔG 數學上不含 barrier，再精確的 FEP 也補不出來。

兩個小提案，都不動你的主流程：(a) 濕實驗前先講好我們比的是編輯專一性（kinetic/functional 終點），把 ΔG_bind 定位成必要非充分；(b) 只在漏斗底部 10–50 個候選掛一條可選側枝算 HNH docking 的能障 / 用 MSM 估速率——等於把 v3/v4 那個被拿掉的 specificity tier 用現代工具低成本接回來。這塊正好是我的物理本行，我來做。」