論文タイトル
De Novo Design of Protein Switches with Diffusion-Based Ensemble Sampling
出典
要旨
刺激に応じて複数の立体構造を可逆的に行き来する「タンパク質スイッチ」のde novo設計は、単一の安定構造を狙う既存パイプラインでは扱いにくいテーマです。本論文は、事前学習済みタンパク質拡散モデルに推論時制御を加えることで、参照構造ひとつだけからスイッチに必要な代替構造のアンサンブルを自動生成する手法Diff-Switchを提案しています。
用語
- タンパク質スイッチ: 単一のアミノ酸配列が2つ(以上)の安定な立体構造をとり、リガンド結合や翻訳後修飾などの刺激に応じて可逆的に変換するタンパク質。シグナル伝達や制御の中核を担います
- collective variable(CV): 分子動力学の強化サンプリング分野で、系の複雑な状態を要約する少数の変数。本論文ではドメイン間の重心距離などが該当し、スイッチの開閉状態を区別する軸として使われます
- メタダイナミクス: CV空間上に、既に訪れた領域を避けるよう履歴依存のバイアスを蓄積し、自由エネルギー障壁を越えた探索を促す強化サンプリング手法
解説など
背景
スイッチのde novo設計には、
- 相互変換可能な複数のコンフォメーションを定義すること
- その平衡を偏らせるトリガーを指定すること
- これら全てをコードする単一の配列を特定すること
という3つの要件を同時に満たす必要があります。ProteinMPNN-MSDやDynamicMPNNのようなmulti-state inverse folding法は3.には対応できますが、1.の「複数バックボーンのアンサンブル」自体が既に手元にあることを前提としており、このアンサンブルをどう用意するかが設計全体のボトルネックになっていました。
本論文の特徴
このボトルネックに対し、直近の2つの研究がどう向き合ってきたかを見ると、Diff-Switchの立ち位置がはっきりします。
- Praetorius et al.(2023年, Science)はヒンジ型スイッチの設計で、密に詰まったヘリカルリピートタンパク質を閉状態の出発点とし、開状態は「シフト・回転・接合・Rosetta緩和」という一連の手作業構造操作で人力構築していました。
- Buckley et al.(2025年, bioRxiv)はリン酸化誘導型スイッチの設計で、リン酸化部位が埋没した閉状態のフォールドから、リン酸化部位を含むヘリックス末端を「マスク」してProteinMPNN-MSDに配列予測を任せる、という半自動的な工夫をとっています。
両者に共通するのは、代替コンフォメーションを得るために手作業の構造操作が不可欠だった点です。Diff-Switchは、この「代替バックボーンをどう作るか」という工程そのものを拡散モデルへの推論時制御によって自動化しようとしています。
その仕組みは、「各ドメインの局所構造は参照構造に忠実に保ちつつ、ドメイン同士の大域的な配置だけを多様化させる」という2段構えのサンプリングです。
拡散モデルの事前分布を、twisted diffusion samplerでサンプリングします。これだけだと各ドメインの形は保たれますが、拡散モデルが最も「自然」とみなす配置に毎回収束してしまい、スイッチに必要な別の配置にはほとんど到達しません。
そこで導入されるのがメタダイナミクスです。ドメイン間の重心距離のような1次元のCVを定義し、その上にwell-temperedメタダイナミクスの履歴依存バイアスを構築します。手順はイテレーションを繰り返す形で進みます。既にたくさん探索済みのCV領域には小さい山が、未探索領域には元の高さの山が置かれるため、サンプラーは徐々に「まだ訪れていない配置」へ押し出されていきます。
十分な回数のイテレーションを終えたら、蓄積されたアンサンブルの中から参照構造から最も乖離したバックボーンを選び、そのペアをmulti-state inverse foldingに渡して候補配列を生成し、最後にAlphaFold3でフィルタリングします。パイプラインの後段は既存手法そのままで、Diff-Switchが介入するのはバックボーン生成の段階だけという設計です。
結果
筆者らは、閉状態のみが既知でリン酸化に応じて開くスイッチのde novo設計に取り組んでいます。de novo設計されたαβフォールドを出発点に、リン酸化部位を含むヘリックスとスキャフォールドの重心距離をCVとしてDiff-Switchを適用し、Calibyで16配列を設計してAlphaFold3で評価しました。結果、Diff-Switchによる設計は16配列中4配列(25%)が意図通りのスイッチ挙動(リン酸化なしで閉じ、ありで開く)を示したのに対し、閉状態のみを入力としたbaseline設計は0配列でした。

