Phenylalanine modification in plasma-driven biocatalysis revealed by solvent accessibility and reactive dynamics in combination with protein mass spectrometry

技術要約：プラズマ駆動型バイオカタリシスにおけるフェニルアラニンの修飾

問題提起
バイオカタリシスは、従来の化学触媒に代わる環境に優しい選択肢を提供するが、反応性の高い共基質、特に過酸化水素（$H_2O_2$）の供給に関して重大な課題に直面している。$H_2O_2$は触媒回転に不可欠であるが、高濃度の$H_2O_2$やその他のプラズマ生成種（PGS）は酵素を不活性化させる可能性がある。非熱プラズマ源（誘電体バリア放電など）は、これらの種を制御された濃度で供給できるが、PGS（$H_2O_2$、$OH$、$O$、$O_2$、$NO$、$H$など）と酵素構造との間の具体的な相互作用は、原子レベルでは依然として解明が進んでいない。酵素の失活を防ぎつつ触媒効率を維持するためには、これらの相互作用を理解することが極めて重要である。

手法
本研究では、計算科学的手法と実験的手法を組み合わせ、3つの酵素（モデル非特異的ペルオキシゲナーゼであるAaeUPO、本研究の主要なモデル酵素であるCviUPO、および$H_2O_2$によって不活性化される解糖系酵素であるGapA）とPGSとの相互作用を調査した。

1. 計算機スクリーニング（SASAおよびMD）：

   • 溶媒接近可能表面積（SASA）解析： VMDおよびLAMMPS上に構築されたカスタムパッケージを用い、酵素の溶媒接近可能表面に対するPGSとの相互作用エネルギーをマッピングした。これには、表面から1.4 Åの位置にプローブ分子を配置し、強い相互作用領域を特定するための拘束エネルギー最小化を行うプロセスが含まれる。
   • 反応性分子動力学（rMD）： LAMMPS内でReaxFF力場（Montiらによる生体分子パラメータ化）を用い、動的な結合の切断と形成をシミュレートした。シミュレーションは、タンパク質のアンフォールディングの影響を評価するため、真空および明示的溶媒中で300 K、400 K、500 Kの温度にて実施された。
   • 高濃度シミュレーション： 現実的な拡散条件を模倣するため、酵素の周囲に100個の$OH$ラジカルをランダムに分布させた高濃度のPGSを用いたシミュレーションを行った。

2. 実験的検証：

   • プラズマ処理： 組換えCviUPO（rCviUPO）を誘電体バリア放電（DBD）プラズマ源で5分間処理した。
   • 質量分析（MS）： 処理されたタンパク質を還元、アルキル化し、トリプシンで消化した。LC-MS/MS解析を行い、予測された相互作用エネルギーが高いアミノ酸残基（例：フェニルアラニン、メチオニン、システイン、アルギニン、リジン）に焦点を当てて、特定のアミノ酸修飾を検出した。

主な結果

• 相互作用プロファイル： SASA解析により、PGSは電荷を持つ残基や極性残基と最も強く相互作用することが明らかになった。具体的には、$OH$ラジカルはリシン（Lys）、フェニルアラニン（Phe）、アルギニン（Arg）と強く相互作用した。$H_2O_2$はセリン（Ser）およびグルタミン酸（Glu）と強く相互作用した。
• 温度効果： 温度の上昇（アンフォールディングを模倣）に伴い、相互作用プロファイルが変化した。400 Kおよび500 Kでは、以前は埋もれていた残基が露出したことで、活性部位（CviUPOのヘム中心付近など）での修飾の可能性が高まった。
• 溶媒の役割： 溶媒を含むシミュレーションでは、大部分の反応性種がタンパク質表面ではなく水分子と反応しており、溶媒が保護障壁として機能していることが示唆された。しかし、高濃度シミュレーションでは、時間の経過とともに、種が依然としてSASAスクリーニングで予測された残基に結合することが確認された。
• 実験的確認： プラズマ処理されたCviUPOの質量分析により、フェニルアラニン残基（Phe51, Phe67, Phe208）およびメチオニン残基（Met42, Met220）の顕著な酸化が確認された。
  • フェニルアラニンの修飾は、露出した芳香環に対する$OH$の強い相互作用というSASAの予測と一致している。
  • メチオニンの修飾は、SASAの相互作用スコアが低かったにもかかわらず発生したが、これは著者らによれば、メチオニン残基が埋もれていても、拡散や静的なSASAモデルでは捉えきれない局所的な環境を通じてアクセス可能であったためである。
  • アルギニンおよびリシンの修飾は低強度でしか検出されず、システインの修飾は見られなかった。これは、特定の構造的コンテキスト（例：Cys19はヘムと共有結合している）によるものと考えられる。

意義および主張
本論文は、SASAに基づく相互作用スクリーニングと反応性分子動力学の組み合わせが、プラズマに曝露された酵素における潜在的な修飾部位を予測するための有効なツールであることを主張している。本研究は以下のことを示している：

1. 予測的有用性： 計算手法は、プラズマ誘起酸化の主要な標的としてフェニルアラニンを特定することに成功し、この知見は実験的に検証された。これは、本手法が、時間のかかる実験的努力を特定の残基へと直接導けることを示唆している。
2. メカニズム的洞察： 結果は、溶媒の接近可能性と反応種の拡散経路が極めて重要であることを強調している。液相を通じた拡散が長時間続くと、種が酵素ではなく溶媒と反応する確率が高まる。
3. 手法の統合： 本研究は、計算機スクリーニングと質量分析を統合することの価値を強調している。SASA法は「粗いスキャン」を提供し、フェニルアラニスの解析を優先させることに成功したが、著者らは、高い相互作用エネルギーが必ずしも結合を保証するわけではなく、実験的検証が依然として必要であると述べている。

著者らは、SASA/MDアプローチは初期スクリーニングには効果的であるが、今後の課題は、テストされる修飾の範囲を拡大し、モデルを他の系に適用してその一般的な妥当性を検証することであると結論付けている。本研究は、酵素の不活性化問題を解決したと主張するものではなく、プラズマ誘起損傷の特定の分子標的を理解し予測するための枠組みを提供するものである。