The choice of structural template can create apparent links between predicted protein stability and longevity in tetrapods

技術的要約：構造テンプレート、タンパク質の安定性、および四肢動物の長寿

問題提起
比較ゲノミクスは、種間のアミノ酸変異をスコアリングするために、構造ベースの予測手法にますます依存しており、多くの場合、単一の構造テンプレート（例：ヒトやモデル生物のバックボーン）を適用してタンパク質の安定性を推論している。このような手法は分子レベルの安定性と長寿のような表現型特性を結びつけることを目的としているが、構造テンプレート自体の選択が、種の特徴との間に偽の相関を生み出している可能性については、これまでほとんど注意が払われてこなかった。本研究は、再構築された祖先と現生種との間のアミノ酸の違いによる予測される折り畳み安定性が、四肢動物の長寿と真に関連しているのか、あるいは観察された相関がテンプレート選択によるアーティファクト（人工的な産物）であるのかを調査するものである。具体的には、「長寿な種はより高いタンパク質安定性を維持している」という仮説と、「長寿な種（有効集団サイズが小さい）は、わずかに不安定化させる置換を蓄積している可能性がある」とする「ほぼ中立理論」という、対立する仮説を検証する。

方法論
著者らは、性質の異なる2つの独立した四肢動物の系統進化について比較研究を行った：

1. 胎盤哺乳類： 57種（1,755遺伝子）。
2. 爬虫類および両生類： 40種（1,910遺伝子）。

実験設計およびワークフロー：

• 祖先再構築： 現代のリファレンス種と比較する代わりに、IQ-TREEを用いて祖先タンパク質配列（哺乳類については胎盤哺乳類の祖先、爬虫類・両生類については四肢動物の祖先）を再構築した。
• 安定性予測： 膨大な安定性データセットで学習されたディープラーニング予測器である ThermoMPNN を用いて、祖先と各現生種との間のアミノ酸の違いを評価した。すべての違いについて、折り畳み自由エネルギーの変化（$\Delta\Delta G$）を算出した。
• 構造テンプレート： 予測は単一の構造バックボーン上にマッピングされた。
  • 研究1（哺乳類）： 主要テンプレートは ヒト のAlphaFoldモデルとした。確認テストとして、オポッサム（事後検証用）および カモノハシ （事前指定された確認用）のバックボーンを使用した。
  • 研究2（爬虫類・両生類）： 主要テンプレートは、時間的な等距離性を考慮して選ばれた ゼブラフィッシュ のAlphaFoldモデルとした。確認テストとして、シーラカンス および ヒト のバックボーンを使用した。
• 指標： 主要指標（$M_s$）は、祖先からのアミノ酸の違いあたりの平均予測 $\Delta\Delta G$ である。二次指標は、部位あたりの純 $\Delta\Delta G$ および部位あたりの分岐度（$R_s$）である。
• 統計解析： 系統学的一般化最小二乗法（PGLS）回帰を用い、成体体重を調整した上で、対数変換された最大寿命と安定性指標を関連付けた。解析は、固定された修正案に従った事前登録計画に基づき、有意性テストのためにブートストラップ・リサンプリングを行った。

主な結果

1. 部位レベルでの予測器の妥当性： ThermoMPNN予測器は、期待される部位レベルのシグナルを正確に捉えた：観察されたアミノ酸の違いは、ランダムな代替配列よりも不安定化しにくかった（浄化選択と一致）。また、埋没した部位は、露出した部位よりも不安定化に対して高い感受性を示した。
2. 長寿とのロバストな関連性の欠如：
   • 哺乳類： ヒトのテンプレートを用いた場合、主に霊長類が「最も不安定化していない」ことに起因する弱い負の傾きが観察された。しかし、この相関は外群テンプレート（オポッサム：$P=0.595$; カモノハシ：$P=0.701$）を使用すると消失した。「霊長類のオフセット」（他の胎盤哺乳類からの偏差）は、ヒトテンプレートでの $-1.96$SD から、カモノハシテンプレートでの$-0.62$ SD へと減少した。
   • 爬虫類・両生類： ゼブラフィッシュのテンプレートでは関連は見られなかった（$P=0.799$）。ヒトのテンプレートを用いた場合は、名目上の負の傾きが得られたが（$P=0.012$）、これは進化速度の遅い系統（カメとワニ）が、進化の速い有鱗類から分離していることに起因していた。分岐度や系統を制御すると、この相関は消失した。
3. テンプレート・バイアス： テンプレートの選択が、見かけ上の相関を生じさせていた。サンプル内のテンプレート（例：哺乳類におけるヒト）を使用すると、そのテンプレート自身の系統が有利になり、それらがより安定しているように見える。たとえ「等距離」のテンプレート（爬虫類におけるゼブラフィッシュ）であっても、それが配列空間において等距離でない場合、バックボーンとの適合性が高い、進化の遅い系統を有利に扱うというバイアスが生じ得る。
4. 二次的なコントロール： 長寿な哺乳類は進化が遅い（部位あたりの分岐度が低い）ことが確認されたが、これはアミノ酸の違いの平均的な予測安定性の明確な差にはつながらなかった。

主要な貢献

• テンプレート・バイアスの特定： 本論文は、単一の構造テンプレートの選択が、予測されるタンパク質安定性と長寿との間に、偽の種間相関を生み出す可能性があることを示した。
• 方法論的フレームワーク： 著者らは、比較構造解析におけるこのバイアスを検出し、軽減するための4つの具体的なチェック項目を提案している：
  1. 現代のリファレンスではなく、再構築された祖先からの差異を測定すること。
  2. 外群テンプレートを使用し、別の無関係な外群テンプレートで結果を確認すること。
  3. 見かけ上の相関が、テンプレートの系統または進化の遅い系統に従っていないかをテストすること。
  4. 特定の系統への依存性を、ブートストラップテスト、系統共変量、または系統内解析を用いてテストすること。
• 安定性と長寿の関連に関する「負の結果」： 解析された特定のデータセットと指標において、長寿な四肢動物が、タンパク質の折り畳み安定性（または祖先からの不安定化の低さ）が高いという強固な証拠は見出されなかった。

意義と主張
本論文は、タンパク質の安定性と長寿の関連を示唆するこれまでの知見が、方法論的な選択、具体的には、しばしばヒト中心となる単一の構造テンプレートの使用によるアーティファクトである可能性があると主張している。著者らは、中程度の効果を排除することはできないものの、複数のテンプレート制御を用いた厳格な事前登録解析の結果、予測されるタンパク質安定性と長寿との間の強固な関連を支持する証拠は得られなかったと断言している。

本研究は、構造ベースの予測器は、特定のバックボーンの文脈において変異を評価するものであり、そのバックボーンには、それが設計された配列の「刻印」が残っていることを強調している。したがって、このような比較は、テンプレートに対する感度の検証が必要である。著者らは、これらの防護策を適用した場合、長寿な種が（祖先からの平均 $\Delta\Delta G$ の観点から）系統的に安定したタンパク質を保持している、あるいは長寿な種が系統的に不安定化させる置換を蓄積しているという仮説を、データは支持しないと結論付けている。この結果は、長寿に関連する適応が、孤立した鎖の平均的な折り畳み安定性ではなく、タンパク質界面、複合体形成、または調節機構において作用している可能性を示唆している。