
計算化学では、水モデルは水クラスター、液体の水、および明示的な溶媒を含む水溶液をシミュレートし、熱力学的に計算するために使用されます。モデルは、量子力学、分子力学、実験結果、およびこれらの組み合わせから決定されます。分子の特定の性質を模倣するために、多くの種類のモデルが開発されています。一般に、これらは次の3つの点で分類できます。(i)サイトと呼ばれる相互作用点の数、(ii)モデルが剛性か柔軟性か、(iii)モデルに分極効果が含まれているかどうか。
明示的な水モデルの代替として、連続体モデルとも呼ばれる暗黙的な溶媒和モデルを使用する方法があります。その例としては、COSMO溶媒和モデル、分極連続体モデル(PCM)、ハイブリッド溶媒和モデルなどがあります。[1]
シンプルな水モデル
剛体モデルは最も単純な水モデルと考えられており、非結合相互作用に依存しています。これらのモデルでは、結合相互作用は暗黙的にホロノミック制約によって扱われます。静電相互作用はクーロンの法則を使用してモデル化され、分散力と反発力はレナードジョーンズポテンシャルを使用してモデル化されます。[2] [3] TIP3P(3点の転送可能な分子間ポテンシャル)やTIP4Pなどのモデルのポテンシャルは次のように表されます。
ここで、静電定数k C は、分子モデリングで一般的に使用される単位では332.1 Å·kcal/(mol· e ²) の値を持ちます[要出典] ; [4] [5] [6] q iとq j は、電子の電荷に対する部分電荷です。 r ij は、2つの原子または荷電サイト間の距離です。AとB は、レナード・ジョーンズ・パラメータです。 荷電サイトは、原子上にある場合もあれば、ダミーサイト (孤立電子対など) 上にある場合もあります。 ほとんどの水モデルでは、レナード・ジョーンズ項は、酸素原子間の相互作用にのみ適用されます。
下の図は、3 ~ 6 サイトの水モデルの一般的な形状を示しています。正確な幾何学的パラメータ (OH 距離と HOH 角度) はモデルによって異なります。
2サイト
よく知られている3サイトSPCモデル(下記参照)に基づく水の2サイトモデルは、サイト再正規化分子流体理論を用いて水の誘電特性を予測できることが示されている。[7]
3サイト
3 サイト モデルには、水分子の 3 つの原子に対応する 3 つの相互作用点があります。各サイトには点電荷があり、酸素原子に対応するサイトには Lennard-Jones パラメータもあります。3 サイト モデルは計算効率が高いため、分子動力学シミュレーションの多くのアプリケーションで広く使用されています。ほとんどのモデルでは、実際の水分子と一致するリジッド ジオメトリを使用しています。例外は SPC モデルで、観測された角度 104.5° ではなく、理想的な四面体形状 (HOH 角度 109.47°) を想定しています。
以下の表に、いくつかの 3 サイト モデルのパラメータを示します。
SPC/E モデルは、ポテンシャル エネルギー関数に平均分極補正を追加します。
ここで、μは有効分極した水分子の電気双極子モーメント( SPC/Eモデルでは2.35 D )、μ 0は孤立した水分子の双極子モーメント(実験では1.85 D)、α iは等方性分極率定数で、その値は1.608 × 10 −40 F ·m 2。モデル内の電荷は一定なので、この補正により総エネルギーに 1.25 kcal/mol (5.22 kJ/mol) が追加されるだけです。SPC/E モデルでは、SPC モデルよりも密度と拡散定数が向上します。
CHARMM力場に実装されたTIP3Pモデルは、元のモデルをわずかに修正したバージョンです。違いはレナード・ジョーンズパラメータにあります。TIP3Pとは異なり、CHARMMバージョンのモデルでは、酸素に加えて水素原子にもレナード・ジョーンズパラメータを配置します。電荷は変更されません。[12] 3サイトモデル(TIP3P)は、比熱の計算において優れたパフォーマンスを発揮します。[13]
柔軟なSPC水モデル

柔軟な単純点電荷水モデル(または柔軟なSPC水モデル)は、 3サイトSPC水モデルの再パラメータ化です。[14] [15] SPCモデルは剛性ですが、柔軟なSPCモデルは柔軟です。ToukanとRahmanのモデルでは、O-H伸縮が非調和化されているため、動的挙動がよく説明されています。これは、分極を考慮しない最も正確な3中心水モデルの1つです。分子動力学シミュレーションでは、水の正確な密度と誘電率を与えます。[16]
柔軟な SPC は、 MDynaMixおよびAbaloneプログラムに実装されています。
その他のモデル
- ファーガソン(フレキシブルSPC)[17]
- CVFF(フレキシブル)
- MG(柔軟かつ解離性)[18]
- KKYポテンシャル(柔軟なモデル)[19]
- BLXL(スミア帯電電位)[20]
4サイト
4 サイトモデルは、3 サイトモデルの HOH 角の二等分線に沿って酸素の近くにダミー原子を 1 つ追加することにより、4 つの相互作用点を持ちます (図ではMで示されています)。ダミー原子は負の電荷のみを持ちます。このモデルは、水分子の周りの静電分布を改善します。このアプローチを使用した最初のモデルは、1933 年に発表された Bernal–Fowler モデルであり、[21]これはおそらく最も古い水モデルでもあります。ただし、BF モデルは、密度や蒸発熱などの水のバルク特性をうまく再現しないため、歴史的な関心事にすぎません。これはパラメーター化方法の結果です。現代のコンピューターが利用可能になった後に開発された新しいモデルは、メトロポリス モンテ カルロまたは分子動力学シミュレーションを実行し、バルク特性が十分に再現されるまでパラメーターを調整することによってパラメーター化されました。
1983年に初めて発表されたTIP4Pモデルは、計算化学ソフトウェアパッケージに広く実装されており、生体分子システムのシミュレーションによく使用されています。その後、TIP4Pモデルは特定の用途向けに再パラメータ化されてきました。TIP4P-Ewモデルはエワルド和法で使用し、TIP4P/Iceは固体水氷のシミュレーションに、TIP4P/2005は凝縮水の状態図全体をシミュレーションするための一般的なパラメータ化であり、TIP4PQ/2005は同様のモデルですが、シミュレーションに量子効果が含まれている場合に固体水と液体水の特性を正確に記述するように設計されているモデルです。[22]
4 サイト水モデルのほとんどは、自由水分子のものと一致する OH 距離と HOH 角度を使用します。例外は OPC モデルで、水分子の基本的な C 2v 分子対称性以外の幾何学的制約は課されません。代わりに、点電荷とその位置は、水分子の静電気を最もよく表すように最適化されます。OPC は、一般的に使用されているいくつかの剛性nサイト水モデルよりも正確に、包括的なバルク特性を再現します。OPC モデルは、AMBER力場内に実装されています。
その他:
- q-TIP4P/F(フレキシブル)[29]
- TIP4P/2005f(フレキシブル)[30]
5サイト
5サイトモデルは、負電荷を、酸素原子の孤立電子対を表すダミー原子( Lと表記)に配置し、四面体のような形状にする。このタイプの初期のモデルは、1971年に提案されたベン・ナイムとスティリンガーのBNSモデルであり[要出典] 、すぐに1974年にスティリンガーとラーマンのST2モデルに引き継がれた。[31]主に計算コストが高いため、5サイトモデルは、2000年にマホニーとジョルゲンセンのTIP5Pモデルが発表されるまではあまり開発されなかった。[32]以前のモデルと比較すると、TIP5Pモデルでは、水二量体の形状が改善され、中性子回折による実験的な動径分布関数と水の最大密度の温度をよりよく再現する、より「四面体」の水構造となっている。 TIP5P-E モデルは、エワルド和で使用するために TIP5P を再パラメータ化したものです。
ただし、BNS モデルと ST2 モデルでは、静電項にクーロンの法則を直接使用するのではなく、スイッチング関数S ( r )を乗じて短距離で縮小された修正バージョンを使用していることに注意してください。
したがって、R LおよびR Uパラメータは BNS と ST2 にのみ適用されます。
6サイト
もともと水/氷系の研究のために設計された、4サイトモデルと5サイトモデルのすべてのサイトを組み合わせた6サイトモデルが、Nadaとvan der Eerdenによって開発されました。[34]このモデルは、周期静電条件(エワルド和)で使用すると融点が非常に高かったため[35]、クーロン相互作用を推定するためにエワルド法を使用して最適化された修正版が後に公開されました[36]。
他の
- 生体分子シミュレーションにおける溶質挙動に対する明示的な溶質モデルの効果も広く研究されている。明示的な水モデルは折り畳まれていないペプチドの特定の溶媒和とダイナミクスに影響を与えるが、折り畳まれたペプチドの立体配座挙動と柔軟性は損なわれないことが示された。[37]
- MBモデル。メルセデスベンツのロゴに似た、2次元システムにおける水の特徴を再現する、より抽象的なモデル。これは「現実の」(つまり3次元の)システムのシミュレーションには使用されないが、定性的研究や教育目的には役立つ。[38]
- 粗視化モデル。水の1サイトモデルと2サイトモデルも開発されている。[39]粗視化モデルでは、各サイトは複数の水分子を表すことができる。
- 多体モデル。トレーニングセット構成を使用して構築された水モデルは量子力学的に解かれ、機械学習プロトコルを使用してポテンシャルエネルギー面を抽出します。これらのポテンシャルエネルギー面はMDシミュレーションに入力され、凝縮相システムの物理的特性を前例のない精度で計算します。[40]
- 多体モデルのもう一つの分類法[41]は、基礎となる静電気の展開に基づくものであり、例えばSCME(Single Center Multipole Expansion)モデル[42]などがある。
計算コスト
水シミュレーションの計算コストは、水モデル内の相互作用サイトの数に応じて増加します。CPU 時間は、計算する必要のある原子間距離の数にほぼ比例します。3 サイト モデルの場合、水分子の各ペアに 9 つの距離が必要です (1 つの分子のすべての原子と他の分子のすべての原子、つまり 3 × 3)。4 サイト モデルの場合、10 の距離が必要です (すべての荷電サイトとすべての荷電サイト、および O–O 相互作用、つまり 3 × 3 + 1)。5 サイト モデルの場合、17 の距離が必要です (4 × 4 + 1)。最後に、6 サイト モデルの場合、26 の距離が必要です (5 × 5 + 1)。
分子動力学で剛体水モデルを使用する場合、制約アルゴリズムを使用して構造を制約したままにしておくことに関連する追加コストが発生します(ただし、結合長が制約されている場合、多くの場合、時間ステップを増やすことができます)。
参照
参考文献
- ^ Skyner RE, McDonagh JL, Groom CR, van Mourik T, Mitchell JB (2015 年 3 月). 「溶液自由エネルギーの計算方法と溶液中のシステムのモデリングに関するレビュー」(PDF) .物理化学化学物理学. 17 (9): 6174–91. Bibcode :2015PCCP...17.6174S. doi : 10.1039/C5CP00288E . PMID 25660403.
- ^ Allen MP、Tildesley DJ (1989)。液体のコンピュータシミュレーション。Clarendon Press。ISBN 978-0-19-855645-9。
- ^ Kirby BJ. マイクロおよびナノスケールの流体力学:マイクロ流体デバイスにおける輸送。
- ^ Swails JM、Roitberg AE (2013)。「prmtop ファイル of {A}mber」(PDF)。
- ^ Swails JM (2013). 一定pHでの複雑な生物システムの自由エネルギーシミュレーション(PDF)。フロリダ大学。
- ^ Case DA、Walker RC、Cheatham III TE、Simmerling CL、Roitberg A、Merz KM、他。 (2019年4月)。 「Amber 2019 リファレンス マニュアル (Amber18 および AmberTools19 をカバー)」(PDF)。
- ^ Dyer KM、Perkyns JS、Stell G、 Pettitt BM (2009)。「サイト再正規化分子流体理論:水の2サイトモデルの有用性について」。分子物理学。107 (4–6): 423–431。Bibcode :2009MolPh.107..423D。doi : 10.1080/ 00268970902845313。PMC 2777734。PMID 19920881。
- ^ Jorgensen, William L. (1981). 「液体の量子力学と統計力学の 研究。10. 水、アルコール、エーテルの転送可能な分子間ポテンシャル関数。液体の水への応用」。アメリカ化学会誌。103 (2)。アメリカ化学会 (ACS): 335–340。doi :10.1021/ja00392a016。ISSN 0002-7863 。
- ^ HJC Berendsen、JPM Postma、WF van Gunsteren、および J. Hermans、「分子間力」、B. Pullman 編集 (ライデル、ドルドレヒト、1981 年)、p. 331.
- ^ ab Jorgensen WL、Chandrasekhar J、Madura JD、Impey RW、Klein ML (1983)。「液体の水をシミュレートするための単純なポテンシャル関数の比較」。The Journal of Chemical Physics。79 ( 2): 926–935。Bibcode :1983JChPh..79..926J。doi :10.1063/1.445869。
- ^ ベレンセン HJ、グリゲラ JR、ストラツマ TP (1987)。 「有効ペアポテンシャルの欠落項」。物理化学ジャーナル。91 (24): 6269–6271。土井:10.1021/j100308a038。
- ^ MacKerell AD、Bashford D、Bellott M、Dunbrack RL、Evanseck JD、Field MJ、他 (1998 年 4 月)。「タンパク質の分子モデリングとダイナミクス研究のための全原子経験的ポテンシャル」。The Journal of Physical Chemistry B。102 (18): 3586–616。doi :10.1021/jp973084f。PMID 24889800 。
- ^ Mao Y, Zhang Y (2012). 「剛性水モデルの熱伝導率、せん断粘度、比熱」. Chemical Physics Letters . 542 :37–41. Bibcode :2012CPL...542...37M. doi :10.1016/j.cplett.2012.05.044.
- ^ Toukan K, Rahman A (1985年3月). 「水中の原子運動の分子動力学研究」. Physical Review B. 31 ( 5): 2643–2648. Bibcode :1985PhRvB..31.2643T. doi :10.1103/PhysRevB.31.2643. PMID 9936106.
- ^ ベレンセン HJ、グリゲラ JR、ストラツマ TP (1987)。 「有効ペアポテンシャルの欠落項」。物理化学ジャーナル。91 (24): 6269–6271。土井:10.1021/j100308a038。
- ^ Praprotnik M, Janezic D, Mavri J (2004). 「水の振動スペクトルの温度依存性: 分子動力学シミュレーション研究」Journal of Physical Chemistry A . 108 (50): 11056–11062. Bibcode :2004JPCA..10811056P. doi :10.1021/jp046158d.
- ^ Ferguson, David M. (1995年4月). 「柔軟な水モデルのパラメータ化と評価」. Journal of Computational Chemistry . 16 (4): 501–511. doi :10.1002/jcc.540160413. S2CID 206038409. 2021年7月28日閲覧。
- ^ MGモデル Archived 2016-03-04 at the Wayback Machine .
- ^ 熊谷 暢、河村 健、横川 哲也 (1994). 「H2O の原子間ポテンシャルモデル: 水と氷の多形への応用」.分子シミュレーション. 12 (3–6). Informa UK Limited: 177–186. doi :10.1080/08927029408023028. ISSN 0892-7022.
- ^ Burnham CJ、Li J、Xantheas SS、Leslie M (1999)。「第一原理による Thole 型全原子分極水モデルのパラメータ化と水クラスター (n=2–21) と氷 Ih のフォノンスペクトルの研究への応用」。The Journal of Chemical Physics。110 ( 9 ): 4566–4581。Bibcode :1999JChPh.110.4566B。doi :10.1063/1.478797。
- ^ ab Bernal JD、Fowler RH (1933)。「水とイオン溶液の理論、特に水素イオンとヒドロキシルイオンについて」。The Journal of Chemical Physics。1 ( 8): 515。Bibcode : 1933JChPh...1..515B。doi :10.1063/1.1749327。
- ^ McBride, C.; Vega, C.; Noya, EG; Ramirez, R.; Sese', LM (2009). 「氷相における量子寄与: 水の新しい経験的モデルへの道—TIP4PQ/2005」. J. Chem. Phys . 131 (2): 024506. arXiv : 0906.3967 . Bibcode :2009JChPh.131b4506M. doi :10.1063/1.3175694. PMID 19604003. S2CID 15505037.
- ^ Jorgensen (1982). 「液体の水と水溶液のシミュレーションのための改訂版TIPS」. The Journal of Chemical Physics . 77 (8): 4156–4163. Bibcode :1982JChPh..77.4156J. doi : 10.1063/1.444325 .
- ^ Horn HW、Swope WC、Pitera JW、Madura JD、Dick TJ、Hura GL、Head-Gordon T (2004 年 5 月)。「生体分子シミュレーションのための改良型 4 サイト水モデルの開発: TIP4P-Ew」。The Journal of Chemical Physics。120 ( 20): 9665–78。Bibcode :2004JChPh.120.9665H。doi : 10.1063 /1.1683075。PMID 15267980。S2CID 39545298 。
- ^ Abascal JL、Sanz E、García Fernández R、 Vega C (2005 年 6 月)。「氷と非晶質水の研究のための潜在的モデル: TIP4P/Ice」。The Journal of Chemical Physics。122 ( 23): 234511。Bibcode : 2005JChPh.122w4511A。doi : 10.1063/1.1931662。PMID 16008466。S2CID 8382245。
- ^ Abascal JL、Vega C (2005 年12 月)。「水の凝縮相の汎用モデル: TIP4P/2005」。The Journal of Chemical Physics。123 ( 23): 234505。Bibcode :2005JChPh.123w4505A。doi : 10.1063/ 1.2121687。PMID 16392929。S2CID 9757894 。
- ^ Izadi S、Anandakrishnan R 、 Onufriev AV (2014 年 11 月)。「水モデルの構築: 異なるアプローチ」。The Journal of Physical Chemistry Letters。5 ( 21 ) : 3863–3871。arXiv : 1408.1679。Bibcode :2014arXiv1408.1679I。doi : 10.1021 /jz501780a。PMC 4226301。PMID 25400877。
- ^ Piana S、Donchev AG、 Robustelli P、Shaw DE (2015 年 4 月)。「水分散相互作用は、無秩序なタンパク質状態のシミュレーションによる構造特性に強く影響する」。The Journal of Physical Chemistry B。119 ( 16 ) : 5113–23。doi : 10.1021 /jp508971m。PMID 25764013。
- ^ Habershon, S.; Markland, TE; Manolopoulos, DE (2009). 「柔軟な水モデルのダイナミクスにおける競合する量子効果」J. Chem. Phys . 131 (2): 024501. arXiv : 1011.1047 . Bibcode :2009JChPh.131b4501H. doi :10.1063/1.3167790. PMID 19603998. S2CID 9095938.
- ^ Gonzalez, MA; Abascal, JJF (2011). 「TIP4P/2005 に基づく水の柔軟なモデル」J. Chem. Phys . 135 (22): 224516. Bibcode :2011JChPh.135v4516G. doi :10.1063/1.3663219. PMID 22168712.
- ^ abc Stillinger FH、Rahman A (1974)。 「分子動力学による液体の水の改良シミュレーション」。The Journal of Chemical Physics。60 ( 4): 1545–1557。Bibcode :1974JChPh..60.1545S。doi :10.1063/1.1681229。S2CID 96035805 。
- ^ ab Mahoney MW、 Jorgensen WL (2000)。「液体の水の5サイトモデルと、剛性の非分極ポテンシャル関数による密度異常の再現」。The Journal of Chemical Physics。112 ( 20): 8910–8922。Bibcode : 2000JChPh.112.8910M。doi :10.1063/1.481505。S2CID 16367148。
- ^ Rick SW (2004 年 4 月). 「Ewald 和で使用するための 5 サイト水ポテンシャル (TIP5P) の再最適化」. The Journal of Chemical Physics . 120 (13): 6085–93. Bibcode :2004JChPh.120.6085R. doi :10.1063/1.1652434. PMID 15267492.
- ^ Nada, H. (2003). 「融点付近の氷と水のシミュレーションのための分子間ポテンシャルモデル:H 2 Oの6サイトモデル」. The Journal of Chemical Physics . 118 (16): 7401. Bibcode :2003JChPh.118.7401N. doi :10.1063/1.1562610.
- ^ Abascal JL、Fernández RG、 Vega C、Carignano MA (2006 年 10 月)。「水の 6 サイトポテンシャルモデルの融点」。The Journal of Chemical Physics。125 ( 16): 166101。Bibcode : 2006JChPh.125p6101A。doi : 10.1063/ 1.2360276。PMID 17092145。S2CID 33883071 。
- ^ Nada H (2016年12月). 「2Oと分子動力学シミュレーション」. The Journal of Chemical Physics . 145 (24): 244706. Bibcode :2016JChPh.145x4706N. doi :10.1063/1.4973000. PMID 28049310.
- ^ Florová P、Sklenovský P、Banáš P、Otyepka M (2010 年 11 月)。「明示的な水モデルは、折り畳まれたペプチドの立体配座挙動と柔軟性はそのままで、折り畳まれていないペプチドの特定の溶媒和とダイナミクスに影響を与えます」。Journal of Chemical Theory and Computation。6 ( 11): 3569–79。doi : 10.1021/ ct1003687。PMID 26617103 。
- ^ Silverstein KA, Haymet AD, Dill KA (1998). 「水と疎水性効果の簡単なモデル」アメリカ化学会誌. 120 (13): 3166–3175. doi :10.1021/ja973029k.
- ^ Izvekov S, Voth GA (2005 年 10 月). 「液体状態のシステムのマルチスケール粗視化」. The Journal of Chemical Physics . 123 (13). AIP Publishing: 134105. Bibcode :2005JChPh.123m4105I. doi :10.1063/1.2038787. PMID 16223273.
- ^ Medders GR 、Paesani F (2015 年 3 月)。「第一原理多体分子動力学による液体水の赤外線およびラマン分光法」。Journal of Chemical Theory and Computation。11 ( 3): 1145–54。doi :10.1021/ct501131j。PMID 26579763。
- ^ Cisneros GA、Wikfeldt KT、Ojamäe L、Lu J、Xu Y、Torabifard H、他。 (2016年7月)。 「水中の分子相互作用のモデリング: ペアワイズから多体のポテンシャル エネルギー関数まで」。化学レビュー。116 (13): 7501–28。土井:10.1021/acs.chemrev.5b00644。PMC 5450669。PMID 27186804。
- ^ Wikfeldt KT、Batista ER、Vila FD、Jónsson H (2013 年10 月)。「単一中心多極展開に基づく転送可能な H2O 相互作用ポテンシャル: SCME」。 物理化学化学物理学。15 (39): 16542–56。arXiv : 1306.0327。Bibcode : 2013PCCP ... 1516542W。doi : 10.1039/c3cp52097h。PMID 23949215。S2CID 15215071 。
