1.チトクロム酸化酵素における第一の電子受容体であるCuA部位について、電子構造と力場を非経験的分子軌道法により解析した。主目的は、分子動力学シミュレーションを行うためのポテンシャル関数を構築することにある。関連して、この部位が特異な複核銅錯体を形成していることの機能的な意味について考察を試みた。CuA部位のアミノ酸配位子については、側鎖部分を切出し水素で終端したモデル分子で近似した。座標は月原らのX線構造解析データに水素原子を付加した。この際に、非経験的分子軌道法で最適化したアミノ酸側鎖モデルの構造を参照した。銅原子にはStevensらの有効内殻ポテンシャルと基底関数、他原子には4-31G基底を用いた。イオン化ポテンシャルの計算にはRHFとROHF法を、構造最適化にはRHF、UHF、CASSCF法を用いた。比較のために、類似のアミノ酸配位子を持ち、かつ銅中心が単核であるような一連のタンパク質(アズリンなど)についても同様の計算を行った。 2.水中の酸電離プロトン移動の分子的機構について、塩酸とフッ化水素酸を比較しながら議論した。これらは良く知られたハロゲン酸だが、酸の強さの差異の原因、プロトン移動の量子的断熱性、多重プロトン移動機構など未解決の問題が多く残っている。ここで得られた知見と方法論は、膜タンパクにおけるプロトンポンプの経路や機構の議論に密接に関連する。 3.タンパク質シミュレーションにおける従来の標準的相互作用ポテンシャルモデルは、アミノ酸側鎖の電子分極効果を無視しているという点で不十分である。この効果を取り入れるために、「電子陰性度平衡化法」が有望視されているが、動力学計算に応用した場合に電荷およびエネルギーの発散がしばしば起こるという問題がある。これについて根本的な解決となり得る新手法を開発し、水溶媒中の電子移動反応系へ応用した。
|