吉布斯自由能

吉布斯自由能在現代生物物理、計算化學、藥物開發中扮演定量化的核心角色。其精確計算涉及統計力學、自由能微擾、增強抽樣等先進方法。

統計熱力學基礎

G 與配分函數的關係（正則系綜）：
G = −kBT × ln Q + PV

其中 Q 為系統的配分函數（Σ exp(−Ei/kT)）。

反應的 ΔG° 來自反應物與產物的配分函數差。Sackur-Tetrode 方程式（單原子氣體）：
G_m = −RT·ln(qN/N!) + PV

複雜分子需考慮平動 + 轉動 + 振動 + 電子配分函數。

自由能計算方法

1. 熱力學積分（TI）：
   ΔG = ∫₀¹ ⟨∂U/∂λ⟩_λ dλ
   λ 為耦合參數，從 0（A 態）到 1（B 態）

2. 自由能微擾（FEP）：
   ΔG = −kT·ln⟨exp(−ΔU/kT)⟩_A

3. Bennett Acceptance Ratio（BAR）：
   組合 A→B 與 B→A 計算，提升統計效率

4. Umbrella Sampling：
   施加偏壓力沿反應座標抽樣，WHAM 重建 PMF

5. Metadynamics：
   動態填充已訪問狀態的位能，加速稀有事件抽樣

6. Replica Exchange MD（REMD）：
   多溫度副本交換，跨越能量障壁

藥物結合自由能

FEP+（Schrödinger）、Amber TI 等商業/學術工具達 ~1 kcal/mol 預測精度（~3× 親和力）。藥物優化常用 FEP 篩選千萬量級的衍生物。

相對 FEP：
ΔΔG = ΔG_bind(A) − ΔG_bind(B)

通常更精確（系統誤差抵消）。

生物大分子摺疊熱力學

蛋白質 ΔG_folding：
ΔG = ΔH − TΔS

典型值（小蛋白質）：

• ΔG = −5 to −15 kcal/mol（不大）
• ΔH = −50 to −200 kcal/mol
• TΔS = −50 to −180 kcal/mol

H/S 補償（H/S compensation）：相關蛋白質中 ΔH 與 TΔS 線性相關，使 ΔG 變化小。原因爭議（部分為實驗誤差伴隨偏差）。

等溫滴定量熱（ITC）

直接測量 ΔH，搭配 Kd 求 ΔG，計算 −TΔS：
ΔG = −RT·lnKd
ΔS = (ΔH − ΔG)/T

ITC 是熱力學「金標準」量測，但需 mg-mM 量級樣品。

焓-熵補償現象

相關反應系列常顯示 ΔH 與 TΔS 線性相關（slope ≈ 1）。原因：

• 部分為數據處理 artifact（窄 ΔG 範圍）
• 部分為物理：強相互作用同時降焓增有序
• 部分為水熵：疏水效應的 H/S 補償

影響藥物優化策略：「焓驅動」結合通常結合特異性更高。

遠離平衡的熱力學

Ilya Prigogine（1977 諾貝爾化學獎）擴展熱力學至非平衡開放系統：

• 耗散結構（dissipative structures）：穩定遠平衡狀態靠能量輸入
• 生命本質：低熵狀態維持靠 +ΔG 食物 + 排熱
• Jarzynski 等式：⟨exp(−W/kT)⟩ = exp(−ΔG/kT)，從非平衡功量測 ΔG

膜能量學

Mitchell（1961、1978 諾貝爾化學獎）化學滲透理論：
Δμ_H⁺ = ΔG = RT·ln([H⁺]_out/[H⁺]_in) + FΔψ
= −2.303·RT·ΔpH + FΔψ

線粒體 ΔpH ≈ 0.5、Δψ ≈ 150 mV → Δμ_H⁺ ≈ −22 kJ/mol

足以推動 ATP 合成（−30 kJ/mol，需 3-4 個 H⁺ 通過 ATP 合酶）。

酵素催化熱力學

酵素降低 ΔG‡（活化自由能），不改變 ΔG（總反應自由能）。

Kramers 理論（1940）：
k = (ω₀·ω‡)/(2πγ) × exp(−ΔG‡/kT)

含摩擦項 γ，描述溶劑黏度依賴。生物酵素中內部摩擦影響反應速率。

ΔG° 修正

生物化學 ΔG°' 採用：

• pH 7（H⁺ 活性固定，不計於 K）
• T = 25°C 或 37°C
• 標準濃度 1 M（H₂O 除外，活性 = 1）
• 離子強度修正（debye-Hückel）

前沿主題

1. 單分子拉伸：AFM 拉伸蛋白質，量測 ΔG_unfolding
2. Jarzynski 與 Crooks 定理：非平衡量測 ΔG
3. 量子熱力學：奈米尺度的工作與熱定義
4. 資訊熱力學：信息與熵的等價（Maxwell's demon）
5. 生命起源：自由能流如何驅動有機分子自組裝
6. 代謝網路熱力學：基因組規模 ΔG 通量平衡分析
7. CRISPR 結合：sgRNA-DNA ΔG 預測脫靶活性