electrostatic catalysis
結論から言うと、酵素電場による触媒促進（electrostatic catalysis）は、 反応座標に沿った電場ベクトル E⃗ が、基質の双極子変化 Δμ⃗ や電荷再配置 Δq⃗ に仕事をすることで、活性化自由エネルギー ΔG‡ を下げる、という数式モデルで記述できます。
以下に、現在もっとも広く使われている「Warshel モデル」を中心に、数式で整理します。
1. 基本式：電場が活性化自由エネルギーを下げる
酵素内電場 E⃗ が反応座標に沿って整列すると、 基質の双極子モーメント変化 Δμ⃗ に対して仕事をし、 活性化自由エネルギーを下げます。
中心式（Warshel 1998–2014）
意味：
E⃗：酵素が作る内部電場（vector）
Δμ⃗：反応前後で変化する双極子モーメント
ΔΔG‡：電場による活性化自由エネルギーの低下量
電場と双極子が同じ方向なら、活性化障壁が下がる。
⚡ 2. 電荷再配置を含むより一般的な式
反応では双極子だけでなく、部分電荷の再配置も起こるため、 より一般的には以下の式が使われます。
qi：基質の原子 i の電荷
Δϕi：酵素内電場によるポテンシャル変化
これは、電場が基質の電荷分布に対して行う仕事の総和です。
3. 反応座標に沿った電場成分だけを取り出すモデル
反応座標を n^ とすると、 電場の有効成分は以下で与えられます。
そして、
ここで Δμ∥ は反応座標方向の双極子変化。
4. 酵素の「電場整列（electrostatic preorganization）」の数式
Warshel が強調したのは、 酵素は反応に最適化された電場を“事前に整列させている”という点。
そのため、溶液反応との差は：
酵素内電場は溶媒よりも方向性が強く、反応座標に沿っているため、 触媒促進が起こる。
5. QM/MM で使われる「電場寄与の分解式」
QM/MM では、活性化自由エネルギーの電場寄与を以下のように分解します。
ここで
ΔGelec‡：電場による静電安定化
ΔGpol‡：偏極による安定化
ΔGsolv‡：水和構造の再編成
電場項は
で計算されます。
6. 実際の酵素での数値例（代表的な値）
DHFR（ジヒドロ葉酸還元酵素） → 電場による障壁低下：4–6 kcal/mol
セリンプロテアーゼ → 求核攻撃方向の電場整列：3–5 kcal/mol
キナーゼ（ATPリン酸転移） → Mg²⁺＋Asp/Glu による電場：5–10 kcal/mol
これらは QM/MM と実験速度論の両方で支持されています。
まとめ（数式での本質）
酵素触媒の本質は、 反応座標に沿って電場を整列させ、双極子変化に仕事をすること。
その数式モデルは以下の3つに集約されます：
(1) 双極子モデル
(2) 電荷再配置モデル
(3) 反応座標方向の有効電場
これらが、2020年代の酵素電場モデルの標準的な数式です。
次に深掘りするなら
electrostatic_preorganizationの実例
QM_MMで電場寄与をどう分解するか
酵素内電場のベクトル場をどう計算するか
どれを進めたい？
2026年7月22日 | カテゴリー：その他 |
