平衡常數的精確處理涉及活性而非濃度,活性係數理論、非理想溶液行為、現代統計力學定義以及生物學的多元結合平衡共同構成嚴謹的平衡熱力學。
嚴格熱力學平衡常數
基於活性:
K° = Π aᵢ^νᵢ
ν 為化學計量係數(產物正、反應物負)。活性 a = γc/c°(液體)、a = γP/P°(氣體)。
K° 為無因次量,僅依溫度。實驗測量 Kc 通常需要活性係數校正:
Kc = K° × Π γᵢ^(−νᵢ)
活性係數的計算
Debye-Hückel 極限律:log γ_± = −A|z⁺z⁻|√I
Davies 方程式(中等 I):
log γ_± = −A|z⁺z⁻| × (√I/(1+√I) − 0.3I)
Pitzer 方程式:高離子強度(>1 M),涉及多元相互作用參數,CSL(Pitzer)模型現代化學工程廣用。
統計力學的 K
從配分函數:
K° = exp(−ΔG°/RT) = Π (qᵢ^νᵢ / qᵢ°^νᵢ) × exp(−ΔE₀/RT)
q 為單個分子的配分函數(平動 × 轉動 × 振動 × 電子)。Ab initio 計算 K 達 ~1 kcal/mol 精度(CCSD(T) 級別)。
動力學與熱力學的橋樑
細緻平衡(detailed balance):每個正反應與逆反應在平衡時速率相等。
K = k_forward / k_reverse
從動力學測 k_f 與 k_r 求 K。應用:
- 催化劑改 k 但 k_f/k_r 不變 → K 不變
- 多步反應的 K_total = Π Kᵢ
多重平衡與耦合
生物系統含多重平衡:
E + S ⇌ ES → E + P(酵素)
HA ⇌ H⁺ + A⁻(緩衝)
Na⁺/K⁺ ATPase 多步循環
通量平衡分析(FBA):基因組規模代謝模型,假設穩態(非平衡,constant ΔG flux)。
生物大分子結合的精細分析
Scatchard 圖:r/[L] = nK − rK(r = bound/total receptor)
- 線性 → 單一結合位點
- 凹彎 → 多重位點或負協同
- 凸彎 → 正協同
Hill 圖:log(θ/(1−θ)) vs log[L]
- 斜率 = n(Hill 係數)
- 中點 = log K
等溫滴定量熱(ITC)
直接量測:
- N(化學計量)
- Kd(解離常數)
- ΔH(焓變)
- 計算 ΔG = −RT·lnKd、ΔS = (ΔH − ΔG)/T
金標準量測,但需 mg-mM 樣品。
離子液體與工業 K
離子液體(ionic liquids)作為「設計溶劑」改變反應 K:
- 提升選擇性(特定 K 增加)
- 降低能耗
- 環保(取代有機溶劑)
例:糖醇酸化在 [BMIM][BF₄] 中 K 增 100×。
遠平衡熱力學
Ilya Prigogine(1977 諾貝爾):生物學的 K 描述瞬時,但細胞遠離平衡。維持穩態靠:
- ATP 水解持續供能
- 開放系統:物質與能量持續流入流出
- 熵流(entropy flux):生物排熵到環境
Jarzynski 與 Crooks 定理
非平衡功 W 的指數平均:
⟨exp(−W/kT)⟩ = exp(−ΔG/kT)
這允許用「拉伸實驗」(single-molecule force spectroscopy)測 ΔG(與 K)。即使過程不可逆。
現代 SBDD 中的 K
基於結構藥物設計:
- FEP+ 預測 ΔΔG → 改 K 倍數
- MM/PBSA、MM/GBSA:快速但較不精確
- ABFE(absolute binding free energy):金標準
- Active learning 結合 FEP,加速優化
K 演化與基因組規模分析
演化平衡:突變-選擇平衡決定 allele 頻率(Hardy-Weinberg)。化學平衡概念延伸至演化(Kimura、Crow)。
極端條件下的 K
- 高溫高壓:地核、行星內部、超臨界水。K 可變 10⁵ 倍以上
- 低溫:深海、極區。K 隨 T 變化,影響嗜寒生物代謝
- 大氣壓力 vs 海洋深處:壓力影響 ΔV ≠ 0 反應(如 CaCO₃ 溶解度)
前沿:合成生物學的 K 設計
蛋白質工程改變受體 Kd(biased agonism)、抗體親和成熟(K 提升 1000×)、CRISPR sgRNA 親和力優化。AI 設計(AlphaFold、RoseTTAFold)正在改變我們設計 K 的能力。
