[論文レビュー] Flexibility of $\alpha$-helices: Results of a statistical analysis of database protein structures

Eldon Emberly, Ranjan Mukhopadhyay+2 more|arXiv (Cornell University)|Sep 25, 2002

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

ひとこと要約

本研究は、タンパク質データバンク（PDB）の399個のタンパク質構造を用いた主成分分析（PCA）により、αヘリックスの固有の柔軟性を定量的に評価した。3つの主要モード—2つのほぼデジェネレートした曲げモードと1つのねじりモード—が同定され、それらの振幅は独立したガウス分布に従い、準熱的平衡と整合的である。水素結合と排他的体積を組み込んだ単純なばねモデルが、これらのモードをうまく再現でき、低エネルギーの集団運動を組み込むことで、タンパク質構造の精緻化と設計の向上が可能になった。

ABSTRACT

$\alpha$-helices stand out as common and relatively invariant secondary structural elements of proteins. However, $\alpha$-helices are not rigid bodies and their deformations can be significant in protein function ({\it e.g.} coiled coils). To quantify the flexibility of $\alpha$-helices we have performed a structural principal-component analysis of helices of different lengths from a representative set of protein folds in the Protein Data Bank. We find three dominant modes of flexibility: two degenerate bend modes and one twist mode. The data are consistent with independent Gaussian distributions for each mode. The mode eigenvalues, which measure flexibility, follow simple scaling forms as a function of helix length. The dominant bend and twist modes and their harmonics are reproduced by a simple spring model, which incorporates hydrogen-bonding and excluded volume. As an application, we examine the amount of bend and twist in helices making up several coiled-coil proteins. Incorporation of $\alpha$-helix flexibility into structure refinement and design is discussed.

研究の動機と目的

• 剛体仮定を越えてαヘリックスの固有の柔軟性を定量すること。
• αヘリックス内の集団的変形が、熱的平衡と整合的な統計的パターンに従うかどうかを特定すること。
• タンパク質構造モデリングに使用可能な、支配的で単純な物理的モデルの開発。
• コイルドコイルのような生物学的に重要な系におけるヘリックスの柔軟性の関連性を評価すること。
• αヘリックスのソフトで集団的なモードを組み込むことで、より正確なタンパク質構造の精緻化と設計を可能にすること。

提案手法

• SCOPに登録された399の代表的タンパク質構造を用い、L = 10–25アミノ酸残基のαヘリックスを抽出して、構造的主成分分析（PCA）を実施した。
• Cα原子座標を理想ヘリックス（1回転あたり3.6アミノ酸、上昇量1.5 Å/アミノ酸）に逐次整列させ、平均ヘリックス構造を反復計算で算出し、10−4 Å/アミノ酸以内に収束させた。
• 整列済みヘリックス間のCα座標のずれから、3L × 3Lの構造的共分散行列を構築し、固有値と固有ベクトルを求めるために対角化した。
• エネルギー関数の数値的2階微分を用いて、水素結合強度と排他的体積の2つのパラメータを持つばねモデルを構築し、主成分分析で得られた主要モードを再現した。
• PCAで得られた正規モードとモデルの正規モードを比較することでモデルの妥当性を検証した。低エネルギー領域のモードでは、良好な一致が得られた。
• コイルドコイルタンパク質（ロイシン・ジッパー、Mntトランスクリプション調節ドメイン、線溶たんぱく質）にモデルを適用し、観察された変形に寄与するモードの寄与度を評価した。

実験結果

リサーチクエスチョン

• RQ1 天然タンパク質構造における、長さが異なるαヘリックスの支配的集団モードの柔軟性は何か？
• RQ2 これらの柔軟モードの振幅は、熱的平衡と整合的な統計的分布に従うか？
• RQ3 水素結合と排他的体積を含む単純なばねモデルは、実際のタンパク質構造のPCAで観察された主要モードを再現できるか？
• RQ4 特に長めのヘリックスを持つタンパク質において、曲げモードとねじりモードの寄与度はどのように異なるか？
• RQ5 これらのソフトモードを組み込むことで、タンパク質構造の精緻化と設計がどの程度向上できるか？

主な発見

• 3つの支配的モード—2つのほぼデジェネレートした曲げモードと1つのねじりモード—が同定され、αヘリックスの構造的分散の大部分を占める。
• これらのモードの振幅は、独立したガウス分布に従い、特徴的な温度T∗ ≈ 室温と整合的な準熱的平衡状態を示している。
• 支配的モードの固有値は、ヘリックス長さに従って単純なべき乗則に従って変化し、サイズ依存性が明確に示された。
• 水素結合と排他的体積の2つのパラメータのみを有する最小限のばねモデルが、低エネルギー正規モード（高調波を含む）をうまく再現できた。
• コイルドコイルタンパク質では、短いヘリックス（例：ロイシン・ジッパー、Mnt）は第1曲げモードと第1ねじりモードでよく記述できるが、長いヘリックス（例：線溶たんぱく質）では高次の曲げモードの寄与が顕著に現れた。
• これらの結果から、タンパク質モデリングにこれらのソフトモードを組み込むことで、全アトムの自由度を必要とせずとも、エネルギー最適化と構造の現実性の両方を向上させられると示唆された。