本年度はまず、Linear Response Path Following (LRPF)法とThree-Dimensional Reference Interaction Site Model (3D-RISM)法を組み合わせた新たなシミュレーション手法を開発するために、タンパク質が水やイオン等の溶媒種から受ける摂動を3D-RISM法により得られる溶媒種の分布から計算し、LRPF法に組み込むためのソースコードを作成した。続いて、本年度で開発したシミュレーション手法(以下、LRPF/3D-RISM法と呼ぶ)の有用性を検証するために、ユビキチンの尿素変性に適用した。先行研究から、尿素変性のメカニズムとしてαヘリックスよりもβシート構造の方から優先的に変性が起こることが報告されている。LRPF/3D-RISM法をユビキチンの尿素変性に適用した結果、5本のうち4本のトラジェクトリで300 ns以内に顕著な構造変化を引き起こすことに成功した。一方で比較のために実行した、尿素水溶液環境を模した通常の分子動力学シミュレーションでは400 nsのシミュレーションでも目立った構造変化は見られなかった。顕著な構造変化が見られたLRPF/3D-RISMシミュレーションの4本のトラジェクトリにおける二次構造の変化の過程の詳細を調べたところ、βシート構造の部分はシミュレーションの進行と共にランダムコイル構造へと変化した一方で、αヘリックス構造の部分にはほとんど変化が見られず、シミュレーション時間全体にわたって構造がほぼ維持されていた。以上の結果は、尿素水溶液中ではβシート構造の方から優先的に変性が起こるという尿素変性のメカニズムをLRPF/3D-RISMシミュレーションが再現していることを示唆している。得られた研究成果は現在論文として国際誌に投稿準備中である。
|