| 研究課題/領域番号 |
21K12130
|
|---|
| 研究機関 | 立命館大学 |
|---|
研究代表者 |
菊地 武司 立命館大学, 生命科学部, 教授 (90195206)
|
|---|
| 研究期間 (年度) |
2021-04-01 – 2024-03-31
|
| キーワード | 結合自由エネルギー / ファインマン自由エネルギー不等式 / SARS-CoV-2 |
|---|
| 研究実績の概要 |
本研究の目的は、ファインマンの自由エネルギー不等式に基づく相対的結合自由エネルギーの予測法(Method based on free energy inequality of Feynman (Method based on FEIF) に新たに溶媒効果を取り入れることである。これまでの方法（T. Ashida and T. Kikuchi, J. Comut-Aided Mol. Des. 27, 479 (2013),T. Ashida and T. Kikuchi, J. Comut-Aided Mol. Des.34, 647 (2020)）では、分子動力学シミュレーションの部分は陽溶媒を用いていたが、ファインマンの自由エネルギー不等式による相対的結合自由エネルギー計算の部分は陰溶媒を用いた計算、すなわちMM-GB/SA法を用いている。そして、自由エネルギー摂動法（FEP法）、MM-PB(GB)/SA法、線形相互作用エネルギー（LIE）法と同程度の結果が得られている。本研究では、この部分の計算をreference interaction site model (RISM)理論を導入を目指す。それが成功すれば溶媒和自由エネルギーの計算精度が高まり、比較的計算コストは低いにも関わらず、タンパク質-リガンド結合自由エネルギーの予測精度の高い方法の開発が可能となる。現在FEIF法にRISM理論のFEIE法への実装を進め、Pim-1キナーゼの阻害剤の結合自由エネルギー計算への応用を始めている。新しいFEIF法をSARS-CoV-2に関連するACE阻害剤ドメイン－リガンド系へ応用し、Covid-19治療薬候補の提案を試みる予定である。またさらなる精度の向上を目指しZwanzigの相対的自由エネルギーの式の展開式の次数を上げた計算も進めている。
|
| 現在までの達成度 (区分) |
現在までの達成度 (区分)
2: おおむね順調に進展している
理由
現在、FEIF法に溶媒効果を有効に取り込むため、RISM理論の実装を進めている。FEIF法では、まずタンパク質とリガンドの複合体のMDシミュレーションを陽溶媒を用いて行う。本研究では、300K、1atm、の環境下、20psecの昇温過程、20psecの平衡過程を経て10.2nsecのシミュレーションを行う。この際、リガンドとして、非摂動系リガンドと摂動系リガンドを対象とする。MDシミュレーションの後、FEIF法の手続きを行う。試験的に対象としてPim-1キナーゼの阻害剤を取り上げる。Pim-1キナーゼは、セリン／スレオニンキナーゼの一種であり、様々な組織で発現していることが知られ、その阻害剤はがん治療薬として期待されている。現在のところ、Pim-1キナーゼ（PDBコード：３BGQ）に対し、FEIF法にRISM計算を取り入れて、タンパク質ーリガンド相対的結合自由エネルギー計算を試みている。RISMモデルにはAMBERモジュールを採用した。これまでに、相対誤差が1.86kcal/molから0.12kcal/molへと大幅な改善がなされた。リガンドはtriazolo[4, 3-b]pyridazineを基本骨格とする。この傾向が様々なリガンドにおいて見られるか、パラメータ調整を含め検討を進めている。また別の角度からも精度を上げるため、Zwanzigの相対的自由エネルギーの式の展開式の第２項まで取り込むことも進めているところである。これが成功すれば、FEIF法においてより高い摂動の次数が求められ、精度は大幅に上がる。
|
| 今後の研究の推進方策 |
すでにPim-1キナーゼ阻害剤については、いくつかのリガンドにおいて、予測した結合自由エネルギーの実験値との決定係数（相関係数の2乗）が0.7-0.8程度であることを確認している。この数値は本方法が成功していると見なせることを示している。現在、一般的に良好な結果が得られるか、パラメータ調整および確認を行っている。またRISM理論を導入することにより計算コストが増大することが懸念されており、その場合は水素原子などの効果を重原子に繰り込むunited atom modelを導入することを検討している。またZwanzigの相対的自由エネルギーの式の展開式の高次項を取り込むことも検討中である。これが成功すれば、FEIF法においてより高い摂動の次数が求められ、精度は大幅に上がる。
平行してSARS-CoV-2阻害剤のリード化合物の探索を試みる予定である。まずSARS-CoV-2タンパク質とangiotensin converting enzyme 2(ACE2)との相互作用を阻害する化合物の探索を試みる。SARS-CoV-2に感染した場合、コロナウィルスタンパク質とACE2が結合することにより、コロナウィルスが細胞にとりこまれ、症状が重篤化することが知られている。SARS-CoV-2タンパク質-ACE2複合体構造として、Protein Data Bank(PDB)に登録されている、PDBコードが6VW1であるものを用いる。その界面部分に結合可能な化合物を探索し、探索した化合物の結合自由エネルギーを本方法を用いて計算する。また別の候補はSARS-CoV-2ウィルスの3CLプロテアーゼの阻害剤である。このタンパク質とリガンドの複合体もPDBに登録されている（PDBコード；7K0G）。このタンパク質も本研究の対象とする予定である。
|
| 次年度使用額が生じた理由 |
当該年度の当初の予想金額が、実際の金額より若干多かったため、その分を次年度にまわすこととする。
|
