Efficient Proton Relay Orchestrated by Covalent Bond Switching of Active Amino Acids in Protein Channels

技術要約：タンパク質チャネルにおける活性アミノ酸の共有結合切替による効率的なプロトンリレーのオーケストレーション

問題提起
タンパク質は、質量輸送やエネルギー変換（しばしばプロトン移動を伴う）を含む多様な機能を実行するために、20種類の標準的なアミノ酸に依存している。ペプチド結合や非共有結合的相互作用を通じた構造的安定性はよく理解されているが、アミノ酸の側鎖における末端のX–H共有結合（ここでXはO、S、N、またはC）の役割は、従来、静的なものと見なされてきた。共有結合の切替（例：hHv1のGlu119やM2チャネルのHis）によって媒介されるプロトン移動の特定の事例は特定されているものの、20種類すべてのアミノ酸残基が量子力学的な動的結合切替を通じてどのようにプロトンと相互作用するかについての系統的な理解は欠けている。本研究が取り組む根本的な問いは、アミノ酸を、共有結合切替を通じてプロトン移動を媒介する固有の能力に基づいて分類するための普遍的なメカニズムが存在するかどうかである。

手法
著者らは、ヒト電圧依存性プロトンチャネル（hHv1）内での系統的な変異シミュレーションを用い、鍵となる残基であるGlu119を他の19種類のアミノ酸へと置換した。

• シミュレーション技術: プロトンリレーのダイナミクスをモデル化するために、量子力学/分子力学（QM/MM）分子動力学（MD）シミュレーションを用いた。
• 自由エネルギー計算: 異なる残基におけるX–H結合切替の自由エネルギー障壁（$\Delta G^\ddagger$）を算出するために、拡張適応バイアス力（eABF）法を適用した。自発的な切替が観察されなかった残基については、これらの障壁を推定するために仮説的な経路を構築した。
• 記述子解析: 観察された活性の物理的起源を特定するため、局所的な求核性指数（$N_k$, $f^-_k$）、静電ポテンシャル極値（ESPE）、原子電荷、および結合双極子モーメントを含む11の物理化学的記述子について、絶対スペアマン順位相関を算出した。
• 一般化: 知見は呼吸複合体Iにおいて検証され、3つの異なるタンパク質チャネル（hHv1、Otop1 Cドメイン、およびインフルエンザAウイルスM2）において、活性残基と不活性残基のラジアル分布を解析することで検討された。

主な結果

1. アミノ酸の分類: 本研究は、共有結合切替能に基づく20種類のアミノ酸の二値分類を確立した。
   • 活性残基（13種類）: Glu, Asp, Ser, Tyr, Thr, Cys, Asn, His, Lys, Trp, Arg, Gln, および Phe。これらは側鎖がO–H、S–H、N–H、または$sp^2$ C–Hで終端している。これらは、低エネルギー障壁（$\Delta G^\ddagger < 20$ kJ mol$^{-1}$）のX–H結合切替を介してプロトンリレーを媒介する。
   • 不活性残基（7種類）: Leu, Ala, Val, Gly, Ile, Met, および Pro。これらは$sp^3$ C–H結合で終端している。これらは結合切替を行わず、代わりに隣接する水分子のワイヤーを介してプロトン移動をサポートする。
2. メカニズムのバリエーション:
   • O–H/S–H 切替: GluやAspのような残基は、ほぼ障壁のない切替を示す。Cys（S–H）は、ヒドロキシル含有残基のメカニズムを模倣する。
   • N–H 切替: Lys, Trp, Asn, Gln, および Argは、単一の窒素部位で切替を行う。イミダゾール環を含むHisは、2つの窒素部位での協同的なメカニズムを利用しており、これが効率的な長距離リレーを可能にしている可能性がある。
   • C–H 切替: フェニルアラニン（Phe、$sp^2$ C–H）はC–H結合切替を介して転送を媒介するが、メチオニン（Met、$sp^3$ C–H）は水ワイヤーを形成するための回転再配列に依存する。
3. 活性の決定要因: 結合切替のエネルギー障壁は、主にプロトン受容原子の局所的な求核性（$N_k$, $f^-_k$）および**静電ポテンシャル極値（ESPE）**によって支配される。
   • 活性残基は負のESPE（$< -0.02$ a.u.）と高い求核性を持ち、ヒドロニウムイオン（$H_3O^+$）を惹きつけ、エネルギー障壁を低下させる。
   • 不活性残基は正のESPE（$> 0.02$ a.u.）と低い求核性を持ち、結合切替を抑制するエネルギー障壁を生じさせる。
4. 空間分布: hHv1、Otop1、およびM2チャネルにおいて、活性残基はプロトンリレーを媒介するために中心孔内に優先的に濃縮されている一方、不活性残基は構造的安定化を提供するために周辺領域に集中している。この分布関数関係の詳細は、タンパク質チャネルの具体的なアーキテクチャによって異なる。
5. 一般化可能性: X–H結合切替メカニズムは、呼吸複合体I（His292およびLys282において）でも確認され、これらの部位を他の活性残基（Gln, Asn）へ変異させても持続することから、生物学的システムにおける効率的なプロトン移動の一般的なメカニズムであることが示唆された。

意義および主張
本論文は、電荷、極性、または疎水性に基づく従来の分類を超え、共有結合切替活性に基づくアミノ酸の新しい分類基準を確立したと主張している。著者らは、本研究が生命がいかにして20種類のアミノ酸を利用しているかという基礎的な物理学的洞察を提供すると断言している。具体的には、13種類の活性残基は、低障壁の量子力学的動的結合切替を通じて効率的なプロトン移動を媒介するように構造的・化学的に準備されている一方で、7種類の不活性残基は安定化の役割を果たす。静電ポテンシャル極値を主要な記述子として特定したことは、バイオエナジェティクス（生体エネルギー論）における酵素やチャネルのプロトン移動メカニズムを理解するための予測ツールを提供するものである。