A thermodynamic framework for mapping elastic recoil mechanism across the human proteome

技術要約：ヒトプロテオームにおける弾性回帰メカニズムのマッピングのための熱力学的枠組み

問題提起
弾性タンパク質は生物学的力学において不可欠であり、組織や分子集合体が機械的エネルギーを貯蔵・回収することを可能にしている。エラスチン（脊椎動物の結合組織）やレジリン（節足動物のクチクラ）のような典型的な弾性タンパク質のメカニズムはよく解明されているが、ヒトプロテオームの大部分は弾性機能に関して未探索のままである。配列相同性に依存する現在の同定手法は不十分である。例えば、ヒトのエラスチンには明確なパラログが存在せず、ハエのレジリンには予測されるヒトのオーソログが存在しない。さらに、弾性タンパク質はしばしば反復的、あるいは固有の無秩序構造や架橋構造を持つため、従来の構造解析や遺伝学的解析を困難にしており、これは弾性機能が著しく過小評価されていることを示唆している。この分野には、候補となる弾性タンパク質を特定し、さらに決定的な点として、それらの根底にある熱力学的回帰メカニズム（例：溶媒駆動型 vs バックボーン駆動型）を区別するための一般的な戦略が欠けている。

手法
著者らは、バックボーンのコンフォメーション・エントロピー（$\Delta S_{bb}$）と溶媒エントロピー（$\Delta S_{sol}$）のどちらが弾性回帰に相対的に寄与しているかに基づいて、タンパク質をマッピングする、配列由来の熱力学的枠組みを開発した。

1. 熱力学的パラメータ:

   • バックボーン・エントロピー ($\Delta S_{bb}$): D'Aquinoらによるスケールを用いて推定。ラマチャンドラン・プロットの分布に基づき、伸長に伴うコンフォメーション自由度の損失を計算する。柔軟な残基（グリシン、アラニン）は高い $\Delta S_{bb}$ に寄与する。
   • 溶媒エントロピー ($\Delta S_{sol}$): Schauperlらによるスケールを用い、分子動力学トラジェクトリに適用されたグリッド非均一溶媒理論（GIST）から導出。これは、アミノ酸側鎖の周囲における水分子の秩序化を定量化する。疎水性残基は、露出に伴う溶媒エントロピーの変化に大きく寄与する。
   • マッピング: タンパク質は、x軸を $\Delta S_{bb}$、y軸を $\Delta S_{sol}$ とする2D「エントロピー・マップ」上に投影される。

2. データソースと処理:

   • 本研究では、AlphaFold DBからのヒト参照プロテオーム（23,587エントリー）を分析した。
   • タンパク質は、pLDDTスコアおよびDSSP二次構造注釈に基づき、「折り畳み構造あり（folded）」または「無秩序（disordered）」に分類された。
   • 解析は、構造化されたタンパク質内の機械的特徴が隠蔽されるのを防ぐため、全タンパク質配列および特定された固有無秩序領域（最小50残基、pLDDT < 0.7）に対して行われた。

3. タンパク質言語モデル（PLM）との統合:

   • 単純なアミノ酸組成を超えた配列文法を捉えるため、著者らはESM2（650Mパラメータ）タンパク質言語モデルを用いて、無秩序領域を埋め込み（embedding）化した。
   • レイヤー24、30、33からの埋め込みを連結し、複合特徴ベクトルを作成した。
   • 近傍（neighborhood）は、配列相同性が低い場合でも類似の潜在空間表現を持つタンパク質を特定するために、k近傍グラフ（kNN）とライデン・コミュニティ検出を用いて定義された。

主な結果

• 明確な熱力学的シグネチャ: エントロピー・マップは、典型的な弾性タンパク質を正常に分離した。レジリンは、高い $\Delta S_{bb}$ と、エラスチンよりも高い $\Delta S_{sol}$ の領域にマッピングされた。これは、バックボーンのコンフォメーション・エントロピーによって駆動されるゴムのようなメカニズムであるが、チロシンや電荷を持つ残基の高含有量によって溶媒エントロピーが複雑化していることと一致している。エラスチンは、より低い $\Delta S_{bb}$ と、レジリンよりも低い $\Delta S_{sol}$ の領域にマッピングされた。注目すべきは、エラスチンのこの低い溶媒エントロピーは、その疎水性回帰メカニズムがより高い溶媒エントロピー変化を反映するという当初の仮説に反していた点である。この不一致は、アラニンに富む足場領域が全体的な疎水性を減少させていることや、レジリンにおける電荷を持つ残基の存在といった、特定の配列特性に起因するとされた。
• 候補弾性タンパク質の発見:
  • 皮膚プロテオーム: 皮膚組織に富むタンパク質は、エントロピー・マップの低い領域（$\Delta S_{sol} < 12$ cal/mol·K）に過剰に集積している。これには、表皮における架橋マトリックスや毛髪を形成するケラチン関連タンパク質（KRTAP）や、後期角化包層タンパク質が含まれる。
  • 無秩序領域: 無秩序領域の解析により、エラスチンの近傍には、タイプIIケラチンの末端ドメイン、フィブリリン-2領域、およびムチン-19が含まれることが明らかになった。これらのドメインは、中間径フィラメントや細胞外マトリックスなどの大きな構造体の中に組み込まれており、コンプライアンスの高いスペーサーやエントロピー的なバネとして機能している可能性がある。
  • レジリン近傍: レジリンの近傍には、高チロシン含有のKRTAP（おそらくジチロシンによる架橋）や、相分離凝縮体に関与するRNA結合タンパク質（FUS、EWSR1など）が含まれており、柔軟で一時的に架橋されたポリマーネットワークに対する共通の要求を示唆している。
• 組成のみのモデルの限界: Gly/Pro組成フレームワーク（Rauscherら）との比較では、重複は不完全であった。多くの皮膚関連ドメインはグリシンに富んでいるが、自律的な弾性体としての自己組み立てに必要なプロリン含有量を欠いており、配列文法とパターニングが材料の挙動を決定する重要な要因であることを示している。
• PLMによる洞察: ESM2の埋め込みは、同様の回帰メカニズムを共有するタンパク質が潜在空間内で明確なコミュニティを形成することを示した。
  • エラスチン、レジリン、およびチチンのPEVK領域は、それぞれ異なる熱力学的メカニズムを反映して、別々のコミュニティを占めている。
  • エラスチンの近傍には、細胞外マトリックスタンパク質（フィブリリン、コラーゲン）だけでなく、調節用足場や長い無秩序リンカーを持つRNA結合タンパク質も含まれており、弾性的な配列文法が多様な細胞コンテキストで再利用されていることを示唆している。
  • シャッブリング（配列の入れ替え）実験により、レジリンの近傍は配列のランダム化に対して頑健である（低複雑性組成に起因する）のに対し、エラスチンの近傍は残基の順序（配列文法）に対して非常に敏感であることが示された。

意義と主張
本論文は、弾性タンパク質を探索し、メカニズム的に分類するための一般的な戦略を提供すると主張している。ヒトプロテオームを、バックボーン・エントロピーと溶媒エントロピーによって定義される熱力学的ランドスケープ上にマッピングすることで、配列相同性のみでは検出できなかった多数の候補弾性ドメインを特定している。

主な貢献は以下の通りである：

1. 弾性の範囲の拡大: 本知見は、弾性機能がエラスチンやレジリンのような特殊な材料を超えて、細胞骨格、細胞外マトリックス、および動的な分子機械の構成要素を含む、現在認識されているよりもはるかに広く普及していることを示唆している。
2. メカニズムによる分類: 本フレームワークにより、候補タンパク質を、予測される回帰メカニズム（疎水性/溶媒駆動型 vs バックボーン/ゴム状）に基づいて分類することが可能になり、それらの潜在的な機械的挙動に関する洞察が得られる。
3. 文法と熱力学の統合: タンパク質言語モデルと熱力学的推定を組み合わせることで、弾性挙動が単なるアミノ酸組成だけでなく、配列文法や進化的なコンテキストによってエンコードされていることを実証した。これは、「弾性的」な配列モチーフが、構造的な揺らぎを緩衝したり、高分子集合体における間隔を調節したりするなど、多様な機能のために転用され得ることを明らかにし、専用の弾性タンパク質と、エントロピー的な復元力を伴う固有無秩序領域との境界を曖昧にしている。

著者らは、現在のモデルには限界（例：明示的なコンフォメーション・アンサンブル・サンプリングを行わずに、伸長によるアクセシビリティの変化を仮定している点）があるとしつつも、ヒトプロテオームにおける弾性の分子起源を解明するための基礎的な枠組みを提供するものであると結論付けている。