[論文レビュー] CSAW: a dynamical model of protein folding

Kerson Huang|arXiv (Cornell University)|Jan 12, 2006

Protein Structure and Dynamics参考文献 4被引用数 4

ひとこと要約

CSAW（条件付き自己回避ウォーク）は、自己回避ウォークとエネルギーに基づく受容ルールを組み合わせることで疎水性相互作用と水素結合をモデル化する、確率的でモンテカルロベースの動的モデルである。このモデルは、疎水性クラスターの急速な収縮による溶融グロブュール状態への遷移と、その後のゆっくりとした安定したネイティブ構造への移行といった、主要な折りたたみダイナミクスを再現し、実験的な回転半径データと定性的に一致する。

ABSTRACT

CSAW (conditioned self-avoiding walk) is a model of protein folding that combines SAW (self-avoiding walk) with Monte-Carlo. It simulates the Brownian motion of a chain molecule in the presence of interactions, both among chain residues, and with the environment. In a first model that includes the hydrophobic effect and hydrogen bonding, a chain of 30 residues folds into a native state with stable secondary and tertiary structures. The process starts with a rapid collapse into an intermediate "molten globule", which slowly decays into the native state afer a relatively long quiescent period. The behavior of the radius of gyration mimics experimental data.

研究の動機と目的

• 疎水性力や水素結合といった主要な物理的力の本質を捉えつつ、計算効率が高く物理的に解釈可能なタンパク質折りたたみモデルの開発。
• 疎水性相互作用と水素結合が三次元構造および二次構造の形成にどのように共同で寄与するかの解明。
• 溶融グロブュール中間体を含む折りたたみの動的経路をシミュレートし、実験的観察と比較すること。
• 骨格、側鎖、水素結合といった構造的・エネルギー的特徴を段階的に組み込み、それらの役割を体系的に調査可能にするための手法の構築。
• フル分子動力学とは比較して計算コストが低く、折りたたみメカニズムを研究するための柔軟なフレームワークの提供。

提案手法

• タンパク質鎖を表すために、二面角で定義される接続された「クランク」からなる骨格モデルを用いる。
• 硬い球体の排除による空間的重なり回避を実現するため、条件付き自己回避ウォーク（SAW）を導入する。
• 熱的揺らぎとエネルギー最小化を模倣するために、メトロポリス受容基準を用いたモンテカルロ更新を実装する。
• 構成エネルギー関数を導入：E = -g₁K₁ - g₂K₂、ここでK₁は疎水性残基間の接触数、K₂は水素結合数を表す。
• 疎水性と水素結合を調整可能なパラメータ（g₁, g₂）として扱い、正規化のためkBT = 1を用いる。
• 幾何的および角度的制約を用いて水素結合をモデル化：O–H距離 ≤ 2 Å かつ反平行な結合配置。

実験結果

リサーチクエスチョン

• RQ1 疎水性相互作用のみがタンパク質鎖において三次元構造を形成するプロセスをどのように駆動するか。
• RQ2 折りたたみ過程において、水素結合はαヘリックスのような二次構造を安定化させる役割を果たすか。
• RQ3 このモデルは、実験的に観察された溶融グロブュール中間体状態を再現できるか。
• RQ4 回転半径は時間とともにどのように変化するか。実験的折りたたみキネティクスと一致するか。
• RQ5 溶融グロブュール状態からネイティブ状態への遷移は、核化制御の一次相転移の特徴を示すか。

主な発見

• 疎水性相互作用のみで、明確な三次元構造が形成されるが、顕著な構造的フラクチュエーションを示す。
• 水素結合の追加により、三次元構造が安定化され、αヘリックスの形成が誘導される。
• 鎖は急速に溶融グロブュール状態へ収縮し、その後ゆっくりとネイティブ状態へ移行する。これは実験的折りたたみキネティクスを模倣している。
• 回転半径の時間発展には、収縮過程で一時的な拡張が観察され、これはレイノルズ膨張に類似した挙動を示唆する可能性がある。
• 溶融グロブュール状態は、最終的なネイティブ状態と類似したコンパクトさと構造的特徴（例：ヘリックス、ループ）を示しており、妥当な折りたたみ中間体である可能性が示唆される。
• 溶融グロブュール状態からネイティブ状態への遷移は、ランダムなプロセスであるため、複数回のシミュレーションで一貫性がない。これは、核化駆動の確率的プロセスであることを示している。