AlphaFold3がタンパク質-DNA複合体について教えてくれる最も有用なことは、各原子がどこにあるかではなく、各原子対がどれだけ接触しやすいかだとしたらどうだろう。
その問いを出発点に、華東師範大学と北京大学の研究チームはContactSeekを開発した。本日Natureに掲載されたこの計算フレームワークは、AlphaFold3の接触確率出力を活用し、DNA塩基エディターの特異性を劇的に向上させる。その結果は、タンパク質工学における広く信じられた前提―複合体の精密な三次元構造を知ることが、機能の理解と改善への最善の道である―を覆すものだ。研究チームは、AlphaFold3の二次的な予測（複合体内のすべての残基とすべてのヌクレオチド間のペアワイズ接触確率の行列）が、予測三次元構造では完全に見逃されていた相互作用の違いを捉えることを発見した。
十分ではなかった構造
接触確率（CP）はAlphaFold3における情報のサイドチャネルである。おなじみの原子座標を含む.cifファイルと並んで、すべてのAF3実行は`_confidences.json`ファイルを生成する。このファイルには、入力内のすべてのトークンペアに対する0から1の値の行列が含まれる。Cas9残基K1020と特定のDNAヌクレオチド間の値0.95は、それらの2つの部位が接触しているとAF3が95%の確信を持っていることを意味する。これはアンサンブル量であり、単一のフォールディング状態ではなく、AlphaFold3の拡散サンプリング実行全体の分布を反映している。
第一著者のHaowei Meng氏らはまず、dIプロファイリングと呼ばれる手法を用いてアデニン塩基エディター（ABE）のゲノム全体のオフターゲットプロファイルをマッピングし、複数のガイドRNAにおけるCas9-TadA8eエディターのオフターゲットDNA部位の包括的なカタログを作成した。その後、それらのオフターゲットDNA配列を対応するオンターゲット配列とともにAlphaFold3に入力し、Cas9、シングルガイドRNA、および標的DNAからなる数千の予測三元複合体を生成した。
研究者らが、三次元構造の偏差（RMSDで測定）とCP変化のどちらがオンターゲット複合体とオフターゲット複合体をよりよく区別できるかを比較したところ、結果は明確だった。CPの方が系統的により高感度だった。オンターゲット複合体とオフターゲット複合体の間のRMSD変化は、確実に解析するには小さすぎることが多く、予測のノイズ範囲内に収まるサブオングストロームレベルのずれだった。対照的に、CP変化は明確で局所化された差異を示し、Cas9の既知の機能ドメインに直接マッピングされた。特定のスペーサー位置における一塩基のミスマッチは、全体的な三次元構造が本質的に変化していないように見える場合でも、DNA-RNAハイブリッドのその領域と接触するCas9残基に特徴的なCP変化シグネチャを生み出した。
これは物理的に理にかなっている。タンパク質-DNA-RNA複合体は動的な集合体であり、複数のコンフォメーション状態をサンプリングする。静的構造は一つのスナップショットを捉えるが、接触確率はアンサンブルを符号化する：特定の界面がどの程度の頻度で形成されるか、その安定性、そして配列の摂動にどのように応答するか。オフターゲットDNA配列は側鎖の位置を大きく変えることはないかもしれないが、その側鎖が安定化接触を形成する確率を大幅に低下させる可能性がある――これは単一の構造では明らかにできないが、CP行列なら可能な差異である。
ContactSeekパイプライン
この洞察をもとに、研究チームはContactSeekを構築した。これはCP行列を実用的な工学的ターゲットに変換するモジュラー計算フレームワークである。パイプラインは3段階で動作する。
第一に、各オフターゲットDNA配列について、ContactSeekはデルタCP行列（オフターゲットCP行列とオンターゲットCP行列の差）を計算する。これらのデルタCP値は、DNA配列が意図した標的から逸脱したときに、どの残基-ヌクレオチド接触が弱まるか強まるかを捉える。
第二に、コンセンサス接触領域（CCR）を特定する。これらは、Cas9タンパク質内の隣接する残基の連続領域であり、その接触確率がオフターゲット配列のセット全体にわたって相関した形で変化する。CCRとは本質的に、核酸との相互作用がオフターゲット結合によって集合的に摂動を受けるタンパク質ドメインであり、工学的介入の有望な部位となる。
第三に、ContactSeekは各CCR内の残基を、接触強調相関（contact enhancement correlation）と呼ばれる指標でランク付けする。これは本質的に、残基のCP変化がシーケンシングで測定された実験的オフターゲット編集シグナルとどの程度密接に連動するかを示す。接触確率の変化がオフターゲット編集率と強く相関する残基ほど、特異性決定因子である可能性が高い。
計算からより優れたエディターへ
研究チームはContactSeekを広く使用されているABE8eアデニン塩基エディターに適用した。これは触媒能が向上したTadA8eデアミナーゼをCas9ニッカーゼに融合させたものである。ContactSeekは2つの重要な特異性決定残基を指名した：Cas9 REC2ドメインのK1020と、TadA8eデアミナーゼ自体のH29である。
K1020DはCas9におけるトップヒットだった。Cas9残基K1020は、非標的DNA鎖近くのREC2ドメインのループに位置し、そのDNAとの接触確率は、オフターゲット配列を解析した際に他のほとんどの残基よりも大きく変化した。分子動力学シミュレーションにより、K1020D置換がループを微妙に再配置し、非標的鎖との接触を減少させることが確認された。その結果、オンターゲット相互作用はほぼそのままに、オフターゲット結合を不均衡に penalize する効果が生まれる。
デアミナーゼ側では、ContactSeekはTadA8e H29Dをトップの特異性向上変異として特定した。H29はTadA8eの活性部位ポケット近くに位置し、RNAガイドおよびDNA基質との接触確率がオフターゲット配列全体で系統的に変化していた。H29D置換は、不完全に位置する基質とのデアミナーゼの無差別な係合を低減した。
複合変異体（ABE8e-DD: Cas9 K1020D + TadA8e H29D）は顕著な結果を示した。ゲノム全体のオフターゲットプロファイリングによれば、ABE8e-DDは複数のガイドRNAにおいて、親株ABE8eと比較して総オフターゲット編集を98-99%削減し、一方でオンターゲット編集効率は維持した。既知のいくつかの高忠実度ABE変異体（ABE8e-V28C、ABE8e-V106W、ABE8e-N108Q）との直接比較において、ABE8e-DDは、ターゲットアンプリコンシーケンシング、ゲノム全体のdIプロファイリング、および直交Rループアッセイによって測定されたすべての指標で優れた特異性を達成した。
アデニンエディターを超えて
ContactSeekはABEに限定されない。研究チームはこのフレームワークをまったく異なるシステム、LbCas12aベースのシトシン塩基エディター（CBE）に適用することで、そのモジュール性を実証した。Cas12aは異なるガイドRNA構造とDNA切断機構を持つ別のCRISPRエフェクターである。Cas12a-APOBEC3A CBEのオフターゲット配列を同じパイプラインに入力することで、Cas12a R284EとAPOBEC3A H29Dが特異性決定残基として特定された。得られた変異体はオフターゲット編集を大幅に減少させ、ContactSeekがCRISPRエフェクターファミリー全体に一般化可能であることが確認された。
このフレームワークはRNAオフターゲット編集の低減にも有効であることが示された。特にTadA8eおよびAPOBEC3Aデアミナーゼを使用する塩基エディターは、RNA転写産物を無差別に編集し、トランスクリプトーム内で意図しないA-to-IまたはC-to-U変化を引き起こす可能性がある。ContactSeekはRNAオフターゲットに対するアプローチを適応させ、トランスクリプトーム全体のRNA編集部位を高活性・中活性・低活性のグループに分類し、グループ間のデルタCPを計算し、RNAオフターゲット活性と相関するデアミナーゼ内のCCRを特定した。得られたTadA8e H29D変異は、DNAオンターゲット活性を犠牲にすることなく、RNAオフターゲット編集を同時に低減した。
特異性工学の新しいパラダイム
ContactSeekのより広範な意義は塩基編集を超える。長年にわたり、タンパク質工学者は構造情報（結晶構造、クライオ電子顕微鏡マップ、そして最近ではAlphaFold2による予測構造）に依存して合理的設計を導いてきた。暗黙の前提は、より高い構造解像度が常により優れた工学的洞察をもたらすというものだった。ContactSeekは、統計的アンサンブル特性である接触確率が、動的分子認識に依存するタスクにおいて静的構造を凌駕できることを示すことで、その前提に挑戦する。
これはゲノム編集において特に重要である。特異的活性と無差別活性の差は、タンパク質-DNA複合体がミスマッチ、バルジ、その他の完全標的からの逸脱にどのように応答するかにかかっていることが多い。結晶構造は理想的な状態を示す。接触確率は、基質が不完全である場合にシステムがどのように振る舞うかを示す。そして特異性工学にとって、まさにそれが重要なのである。
このアプローチは実験的にも効率的である。塩基エディターの特異性を改善する従来の戦略は、大規模な変異スクリーニング、指向性進化キャンペーン、または構造に基づく合理的設計（一度に一つの変異）に依存してきた。ContactSeekに必要なのは、オフターゲットシーケンシングデータ（ますます日常的になりつつある）とAF3予測（急速に手頃な価格になりつつある）のみである。数千の可能な変異から、一回の計算パスで少数の高信頼度候補へと探索空間を絞り込む。
コードとデータがGitHubとZenodoでオープンに公開されているContactSeekは、AI駆動の接触分析がどのようにしてより安全なゲノム編集ツールの開発を加速できるかのテンプレートを提供する。分野が塩基編集の治療応用に向かうにつれ、オンターゲット効力を維持しながらオフターゲット活性を系統的に排除する能力は極めて重要になる。ContactSeekが示唆するのは、前進の道はより鮮明な構造スナップショットを必要とするのではなく、代わりに確率に目を向けることかもしれないということだ。
雅子 訳