| 研究概要 |
我々は、最近、実空間グリッドを基底とするハイブリッド型の第一原理分子法(QM/MM法)と新規な溶液論(エネルギー表示の理論)を組み合わせることによって、効率良く、正確に溶液内の化学過程の自由エネルギー変化を計算する方法を開発した(QM/MM-ER法)。本研究課題の主たる目的は、酵素内の様々な化学過程の自由エネルギー変化を計算するためにQM/MM-ER法を拡張することである。本年度は、タンパク質1b4v中に埋め込まれた補酵素FADの還元反応に伴う自由エネルギー変化をQM/MM-ER法を拡張することによつて計算した。この方法の大きな特徴は、FADに付加する余剰電子1個を溶質と見なすことにある。水、アポタンパク、FADが混合溶媒を形成するとして、これらの溶媒について、エネルギー分布関数を構成し、各溶媒からの還元自由エネルギーに対する寄与を計算した。水分子、アポタンパク、FAD、ADPによる還元自由エネルギーへの寄与は、それぞれ、-57.5kcal/mol,-149.7 kcal/mol,-31.9 kcal/mol,55.1 kcal/molと計算された。また、バルクの水の寄与は、21.0 kcal/molである。FADの還元自由エネルギーは、これらの寄与の単純な和として、-163.1 kcal/molと計算される。このような大きな還元の自由エネルギー変化は、タンパク質全体が+3価に帯電している為である。この計算値は、報告されている実験値、-90.6kcal/molから大きくずれている。実験では、タンパク質の電荷を中和するようなカウンターイオンが含まれていることが、その差の原因と考えられる。
また、本年度の後半では、酵素活性サイトのプロトン化状態を決定するための方法論の定式化を行い、プログラム開発を行った。10個の残基からなる最小のタンパク質、シニョリン(iuao)をテスト系として、QM/MM-ERシミュレーションを行った。シニョリンのような小さなタンパク質でも気相中と溶液中の平均構造に大きな違いがあり、それ故に参照系と溶液系のエネルギー分布関数が十分にオーバーラップしないという結果がでた。
|
