Beyond HCD: complementary fragmentation chemistries expand HLA class I immunopeptidome discovery

技術要約：HCDを超えて – 相補的な断片化化学がHLAクラスI免疫ペプチドームの探索を拡大する

問題提起
質量分析ベースの免疫ペプチドミクスは、ペプチド同定における決定的なボトルネックに直面している。タンパク質を複数のペプチドから推論するプロテオミクスとは異なり、免疫ペプチドミクスでは個々のペプチドを分析物かつ治療標的として扱うため、同定されなかったリガンドは、見逃された標的を意味する。HLAクラスIのリガンドは通常、短く（8〜12アミノ酸）、非トリプチックであり、電荷状態が低いことが多く、断片化が不完全な場合もしばしばある。こうした課題にもかかわらず、現在のワークフローのほとんどは、単一の活性化法であるHCD（Higher-Energy Collisional Dissociation）に依存している。比較プロテオミクスにおいては、単一の断片化法ではすべてのペプチドクラスを捉えることはできないことが長らく確立されているが、免疫ペプチドミクスにおいては、代替手法の価値が体系的に評価されてこなかった。さらに、ペプチド・スペクトル一致（PSM）の再スコアリングに使用されるディープラーニングモデルは、ほぼ排他的にHCDスペクトルに基づいて学習されており、CID（Collision-Induced Dissociation）やIRMPD（Infrared Multiphotion Dissociation）といった代替的な化学的手法の価値を過小評価するバイアスを生む可能性がある。

方法論
本研究では、Orbitrap Tribuid Apex MultiOmicsプラットフォームを用い、赤外線レーザーによるIRMPDへのアクセスを可能にする標準的な手法とともに、3つの断片化法（HCD、イオントラップCID、IRMPD）を評価した。

• サンプルセット： 4種類のHLAタイピングされたヒト細胞株（EBV変換された2つのBリンパ芽球様細胞株：BLCL064およびBLCL004、ならびに2つの癌細胞株：HCT116およびMDA-MB-231）をテクニカル・トリプリケートで分析した。
• 取得： 各サンプルはHLAクラスI免疫親和性精製を経てLC-MS/MSにより分析された。各細胞株に対し、特定のエネルギー設定（例：HCDはNCE 30%、CIDは36%、IRMPDはレーザー出力50%）を用いて、3つの断片化法（HCD、CID、IRMPD）が適用された。
• データベース検索および再スコアリング： 生データは、標準的なヒトSwiss-Protプロテオミクスと、Ribo-seq由来のクリプティック（非標準）データベースの2つに対して検索された。極めて重要な点として、公平な比較を保証するために、本研究ではAlphaPeptDeepディープラーニングモデルを断片化法および細胞株の組み合わせごとに個別にファインチューニングした。これにより、HCD用に訓練されたモデルに内在するバイアスに対処した。再スコアリング後、Percolatorによる1% PSMレベルの偽発見率（FDR）フィルタリングが行われた。
• 分析： 研究では、HCD、CID、IRMPDの3つの断法について、同定収量、フラグメントイオンのラダーの完全性、アレル表現型、および塩基性残基（リシン、アルギニン、ヒスチジン）を持つペプチドと持たないペプチドのスペクトル品質（検索スコアおよびイオン電流）を評価した。

主要な貢献および結果

1. 手法特異的な再スコアリングの検証： 各化学的手法に対してAlphaPeptDeepをファインチューニングすることは、すべての手法においてユニークなペプチド同定数を大幅に増加させた。その増加率は、標準データベースで中央値+17.1%、クリプティックデータベースで+54.5%であった。これは、代替的な化学的手法を「HCD型の定規」でスコアリングすることを防ぐために、手法特異的なモデルが必須であることを裏付けた。
2. 相補的、かつ冗長ではないレパートリー： 3つの手法は、全体的なサイズや組成において同様の免疫ペプチドミクスをもたらしたが、それらは冗長ではなかった。
   • HCDは、最も完全なフラグメントイオンのラダーを生成し、塩基性C末端アンカーを持つアレル（特にA3スーパータイプのメンバーであるA11:01およびA33:03）に対して優れていた。
   • CIDおよびIRMPDは、塩基性C末端アンカーを欠くアレルにおいて、より大きなシェアを回収した。これらは、KRH（リシン、アルギニン、ヒスチジン）を含まないペプチドにおいて、より高い検索スコアとb/yイオン電流を実現した。
3. スペクトル品質と塩基性： 本研究は、根本的なトレードオフを特定した。すなわち、HCDは信号をより完全なバックボーンラダー全体に分散させるのに対し、CIDおよびIRMPDは信号をより少数の、より強度の高いフラグメントへと集中させる。この集中効果は、可動プロトンモデルの下で断片化するKRHフリーペプチドにおいて、CIDおよびIRMPDに特有のスペクトル品質上の優位性（高いハイパースコアおよび総イオン電流）をもたらした。
4. 装置時間の最適化： 3回のLC-MS/MS実行という固定された予算において、3つの手法（各1回ずつ）に実行を分散させることで、単一の手法を3回繰り返すよりも多くのユニークなペプチドが得られた（+0.4%から+10.9%の利得）。
5. 新規リガンドの発見： 癌細胞株において、免疫エピトープデータベース（IEDB）には存在しない239個の予測HLAバインダーのうち、55%が単一の手法によってのみ検出された。このうち9個は癌関連抗原に由来しており、そのうち6個はHCDのみのワークフローでは見逃されていた。

意義および主張
本論文は、CIDおよびIRMPDがHCD単独では同定されない可能性のある配列を提供するか、あるいは化学的に定義されたサブセットのリガンド（特に塩基性残基を欠くもの）に対して優れたスペクトル品質を提供する、と主張している。これらの手法の価値は、HCDに取って代わるものではなく、HCDと並んで存在するものである。

著者らは、これらの化学的手法の潜在能力を最大限に引き出すには、現在欠けている以下のインフラストラクチャが必要であると強調している。

• HCDからファインチューニングされたものではなく、CIDおよびIRMPDの免疫ペプチドミクスデータに対してネイティブに訓練された、フラグメント強度および保持時間モデル。
• 固定の設定ではなく、前駆体の特性（例：塩基性）に基づいて活性化法を割り当てる取得スキーム。
• データ非依存取得（DIA）ワークフローにおける評価。

こうしたインフラが構築されるまでの間、本論文の実践的な推奨事項は控えめながらも実行可能なものである。すなわち、装置の時間が許すのであれば、複数の断片化手法に実行を分散させるべきである。HLA型が既知である場合は、A2スーパータイプのアレル（疎水性アンカー）が支配的なサンプルにはCIDまたはIRMPDを、A3スーパーテールのアルレル（塩基性アンカー）が支配的なサンプルにはHCDを重視して割り当てるべきである。このアプローチにより、ユニークなHLAリガンドの回収を最大化し、潜在的な治療標的のランドスケープを拡大することができる。