Hydration Free Energies of Linear Alkanes: Systematic Deviations in Common Water Models and Their Correction

技術的サマリー：直鎖アルカンの水和自由エネルギー：一般的な水モデルにおける体系的な偏差とその補正

問題提起
一般的な分子力場は、疎水性溶質の水和自由エネルギー（$\Delta G_{hyd}$）を体系的に過大評価しており、シミュレーションにおいて疎水効果が誇張される結果を招いている。この不一致は、界面活性剤の臨界ミセル濃度の過小評価や、油 - 水界面吸着自由エネルギーの過大評価として現れる。アルカンは非極性溶媒和を理解するための基本的なモデル系として機能するが、標準的なパラメータ化（具体的には、アルカン用の TraPPE-UA 力場と、標準的なローレンツ - ベルトロ混合則を適用した各種水モデルの組み合わせ）は、実験的な溶解度データを正確に再現できない。先行研究ではこれらの偏差が指摘されていたが、アルカン鎖長（メタンからエイコサンまで）の広い範囲および複数の水モデルに適用可能な体系的な補正は欠如していた。

手法
著者らは、等温等圧アンサンブル内における自由エネルギー摂動（FEP）を用いて、メタン（$C_1$）からエイコサン（$C_{20}H_{42}$）までの直鎖アルカンの水和自由エネルギーを計算した。本研究では、アルカン用に TraPPE-UA 団結原子力場を採用し、以下の 4 つの異なる水モデルと組み合わせて使用した。

• 3 サイトモデル: SPC/E および OPC3。
• 4 サイトモデル: TIP4P/2005 および OPC。

シミュレーションは、化学的変換中の特異性を回避するためにソフトコア・レナード・ジョーンズポテンシャルを用いた LAMMPS により実施された。結合パラメータ $\lambda$ は 40 ウィンドウで変化した。本研究では、初期比較のために標準的なローレンツ - ベルトロ混合則を採用し、シフトされていないポテンシャルに対して長距離尾部補正を適用した。

体系的な過大評価を補正するため、著者らは水和自由エネルギーの分解に基づく再パラメータ化戦略を採用した。
$$\Delta G_{hyd} = \Delta G_{cavity} + \Delta H_{att}$$
ここで、$\Delta G_{cavity}$ は溶質サイズの空洞を形成する自由エネルギーであり、$\Delta H_{att}$ は引力性溶質 - 溶媒相互作用（レナード・ジョーンズ相互作用エネルギー $\Delta U_{att}$ として近似される）を表す。確立された $\Delta G_{cavity}$ の値を利用することで、著者らは実験的な目標 $\Delta H_{att}$ を導き出し、300 K においてこれらの目標値に一致するようにアルカン - 水間のレナード・ジョーンズ井戸深さパラメータ（$\epsilon$）を調整した。空洞自由エネルギーは $\sigma$ に非常に敏感であるため、距離パラメータ（$\sigma$）は一定に保たれた。

さらに、本研究では以下の評価を行った。

1. TIP4P/2005 および SPC/E 水モデルと組み合わせた一般 Amber 力場（GAFF）。
2. 以前に提案された HH-アルカンモデル（Ashbaugh ら）を $C_{20}$ まで拡張したもの。
3. シフトされたレナード・ジョーンズポテンシャル（カットオフでポテンシャルをゼロにする）を適用した場合と、尾部補正を伴うシフトされていないポテンシャルとの比較による影響。

主要な結果

• 体系的な過大評価: 試験されたすべての水モデル（SPC/E、OPC3、TIP4P/2005、OPC）は、TraPPE-UA と組み合わせて、実験データおよびグループ寄与推定値に対して $\Delta G_{hyd}$ を過大評価した。偏差は 3 サイトモデルでより顕著であり、アルカン鎖長が増加するにつれて増大した。
• 再パラメータ化の成功: アルカン - 水間のレナード・ジョーンズ井戸深さ（$\epsilon$）をローレンツ - ベルトロ値に対して約 5% 増加させることで調整したところ、300 K におけるシミュレーション結果が実験値と見事に一致した。
  • 再パラメータ化されたモデルは、290–350 K の温度範囲にわたって実験データと優れた一致を維持した。
  • 再パラメータ化された TIP4P/2005、SPC/E、および OPC3 モデルは、それぞれ 0.07、0.10、および 0.06 $kT$ の二乗平均平方根偏差（RMSD）を示した。
• GAFF の性能: TIP4P/2005 または SPC/E 水モデルと組み合わせた GAFF 力場は、元の TraPPE-UA/TIP4P/2005 組み合わせと比較して実験値からの偏差が小さかったが、それでも体系的な傾向（TIP4P/2005 との組み合わせにおける小分子アルカンの過大評価、および大分子アルカンの過小評価）を示した。
• HH-アルカンモデル: 当初ブタンまでしか検証されていなかった HH-アルカンモデルをエイコサンまで拡張した結果、$C_{20}$ までの実験データおよびグループ寄与データを再現する能力が確認された。
• ポテンシャルのシフト: 本研究は、球状カットオフにおける値によってレナード・ジョーンズポテンシャルをシフトさせることが、$\Delta G_{hyd}$ に体系的な偏差を導入することを示した。シフトされたポテンシャルに尾部補正を加えることで数学的にシフトされていない結果を回復させることは可能であるが、厳密に補正されない限り、シフト行為そのものが誤差を増幅させることを示しており、化学的自由エネルギーがポテンシャルの切断方式に対していかに敏感であるかを浮き彫りにした。

意義と主張
本論文は、一般的な力場における疎水和の体系的な過大評価が、アルカン - 水間のレナード・ジョーンズ井戸深さパラメータの一歩の調整によって補正可能な頑健な現象であると主張している。著者らは以下の点を強調している。

1. 補正の転送性: $\epsilon$ の約 5% の増加は、異なる水モデル（3 サイトおよび 4 サイトの両方）で一貫しており、非極性溶質に対する一般化可能な補正メカニズムを示唆している。
2. 混合則の限界: 純成分熱力学に基づいてパラメータ化された力場は、特定の相互作用の較正なしに経験的な混合則（ローレンツ - ベルトロなど）を適用された場合、混合物において性能が低下することを結果は浮き彫りにしている。
3. 方法論的洞察: 本研究は、標準的な混合則のみに依存するのではなく、空洞形成自由エネルギーおよび実験的溶解度データに対して相互作用パラメータを直接較正することにより、水和自由エネルギーの予測を改善するための実践的な道筋を提供している。
4. ポテンシャルのシフト: この研究は、自由エネルギー計算におけるシフトされたポテンシャルの使用に関する警告として機能し、そのようなシフトが水和熱力学に重大な体系的誤差を導入し得ることを実証している。

著者らは、非極性種（アルカンなど）の正確な水和自由エネルギーを達成するためには、空洞形成自由エネルギーを用いて相互作用パラメータを直接較正することが可能であり、それによって疎水性の自己集合および溶媒和現象に関する分子シミュレーションの信頼性が向上すると結論づけている。