新最適化アルゴリズムによるタンパク質の折れたたみの研究

点击翻译按钮获取中文摘要

我々はこれまで、タンパク質の立体構造予測問題に、物理学の分野で開発された、徐冷モンテカルロ法やマルチカノニカル法適用することを提唱してきた。徐冷法ではモンテカルロ法とともに、分子動力学法でのシミュレーションが簡単に実現できるが、マルチカノニカル法の分子動力学法版は単純ではない。タンパク質の分野では、モンテカルロ法よりも分子動力学法が使われることが多いので、これは重要な問題である。よって、本年度は、まずマルチカノニカル分子動力学法を開発した。更には、ランジバン法及びハイブリッド法のマルチカノニカル法版も同時に開発した。これら3つの方法の有効性を調べた結果、ハイブリッド法が一番有効なことが示された。更に、アミノ酸数56のタンパク質ProteinGのシミュレーションを継続した。このタンパク質のX線回析実験によって得られた構造は1本のαヘリックスと4本のβストランドからなる1枚のβシート構造を持っている。昨年までは徐冷法によるシミュレーションを行ってきたが、本年度は1/k法によるシミュレーションを行った。初期構造は、徐冷法の時と同じく、完全にランダムな構造を採用した。この計算において、最小エネルギー構造がX線回析実験で得られた構造と同じ位置にαヘリックスを持つことで示した。また、このシミュレーションの結果から任意の温度のカノニカル分布が再現できるが、それから、%αヘリックス度と%βストランド度を残基番号の関数として計算した結果、X線実験で得られた構造と同じ位置にヘリックスやストランドができやすいという結果が得られた。