Quantum chemical modelling of sequence encoded energetics and electrostatic landscapes governing early peptide association in α- synuclein and amyloid β

技術要約：初期ペプチド会合における配列エンコードされたエネルギーおよび静電ランドスケープの量子化学モデリング

問題提起
タンパク質のミスフォールディングと凝集は、パーキンソン病（PD）やアルツハイマー病（AD）などの神経変性疾患の中心的な要因であり、$\alpha$-シヌクレインやアミロイド$\beta$（A$\beta$）といったペプチドの初期の分子会合によって駆動される。実験的研究により成熟した線維の動力学や構造は特徴付けられているが、一過性で不均一な初期会合段階における個々の分子構成のエネルギー的寄与については、依然として解明が進んでいない。特定の残基の線形順序（配列）が、全体のアミノ酸組成とは別に、初期会合イベントを支配する分子特性をどのように形成するかを理解することには、重要なギャップが存在する。本研究では、凝集しやすいペプチドセグメントを組成が一致するスクランブルド・コントロールと比較することで、配列による寄与を単離し、残基の順序がエネルギー、構造組織、および静電ランドスケープにどのように影響するかを決定することを目的とする。

方法論
本研究は、セミエンピリカル量子化学と密度汎関数理論（DFT）を統合した一貫した計算ワークフローを用いて、モノマーおよびホモダイマー状態をモデリングする。

• 系: 凝集しやすい配列セグメント（$\alpha$-シヌクレイン: VTGVTAVAQK; A$\beta$: LVFFAEDVGS）を、組成が一致するスクランブルド・コントロール（AGATQTVVVK および EFLVGSDVAF）と比較した。
• 構造最適化: モノマーおよびホモダイマー構造は、xTBに実装されたGFN2-xTB法（Geometry, Frequency, Noncovalent extended Tight Binding, second generation）を用いて最適化した。ダイマー構造は、最終的な最適化の前に、最初にGFN力場を用いて緩和させた。
• エネルギー論: 会合エネルギー（$\Delta E_{assoc}$）は、最適化されたダイマーの収束したGFN2-xTB電子エネルギーと、独立して最適化されたモノマーのエネルギーの2倍との差として計算された。
• 電子構造解析: 最適化されたモノマーに対し、実空間グリッド（0.20 Å 間隔）上でPBE交換相関汎関数と投影増加平面波（PAW）形式を用いたGPAWによるDFT計算を行った。
• 記述子:
  • 静電特性: 分子静電ポテンシャル（ESP）を電子密度表面（0.001 e Bohr$^{-3}$）上にマッピングし、相対的なランドスケープを可視化するために平均中心化を用いて解析した。
  • 電荷分布: ヒルシュフェルト（Hirshfeld）原子分極解析を用い、統計的指標（平均、中央値、標準偏差）および残基ごとの総和を含む、原子分解された電荷分布を定量化した。
  • 構造/配列: 記述子には、配列電荷デコレーション（SCD）、局所的な疎水性、疎水性ラン長、回転半径（$R_g$）、端間距離、および非球形性が含まれる。
  • 界面: 重原子距離（$\le$ 4.5 Å および 5.0 Å）および特定のN/O/S接触（$\le$ 3.5 Å）による分子間接触を定量化した。

主な結果
本研究は、残基の再配置が、アミノ酸組成が保存されているにもかかわらず、エネルギーおよび構造の両領域において、ペプチド特異的な異なる応答を生じさせることを明らかにしている。

1. モノマーのエネルギー論: 残基の順序は、単離されたモノマーの電子エネルギーを著しく変化させた。A$\beta$の場合、スクランブルド・コントロールは凝集しやすい配列よりもエネルギーが約45.33 kJ/mol低かった。逆に、$\alpha$-シヌクレインの場合、凝集しやすい配列の方がスクランブルド・コントロールよりもエネルギーが約151.12 kJ/mol低かった。
2. 会合エネルギー対界面パッキング:
   • A$\beta$: 凝集しやすい配列（115.617 kJ/mol）とスクランブルド・コントロール（232.552 kJ/mol）の両方が正の会合エネルギーを示しており、これは最適化されたダイマーが、独立して最適化されたモノマーよりも電子エネルギーが高いことを示している。しかし、凝集しやすい配列は、スクランブルド・コントロールよりも低い正の会合エネルギー（116.935 kJ/mol低い）を示した。逆説的ではあるが、この配列は重原子接触数において、コントロール（81）に対し、より少ない接触数（40）を示し、接触密度も低かった。
   • $\alpha$-シヌクレイン: スクランブルド・コントロールは、凝集しやすい配列（-172.819 kJ/mol）よりも有利な（より負の）会導エネルギー（-494.371 kJ/mol）を示した。しかし、凝集しやすい配列は、より広範な界面を形成し、166個の重原子接触（105個に対して）と、より高い密度の近接N/O/S接触（5個に対して15個）を有していた。
   • 結論: 会合エネルギーと界面パッキング（接触数/密度）は線形には変化せず、これらは相互に補完的なものであり、等価ではないことを示している。
3. 静電ランドスケープ (ESP): 残基の再配置は、ESPの空間分布を再編成した。A$\beta$と$\alpha$-シヌクレインのペアは同様のグローバルな電荷状態を共有していたが、正および負のポテンシャル領域の特定の空間配置は、凝集しやすい配列とスクランブルド配列の間で大きく変化した。これは、配列の順序が静電特性の3次元的な提示を決定していることを示している。
4. ヒルシュフェルト電荷分布: 原子電荷解析により、残基の再配置は、グローバルな電荷の大きさを一様にシフトさせるのではなく、主に原子電荷の分布を再分配することが示された。$\alpha$-シヌクレインでは、正/負の原子の割合はほぼ一定であったが、電荷値の確率分布が変化した。A$\beta$では、正/負の原子のバランスにおける緩やかなシフトとともに、再分配が観察された。
5. 配列および構造記述子:
   • A$\beta$: 再配置は、配列電荷デコレーション（SCD）、疎水性ラン長（3残基に対し5残基）、および分子のコンパクトネスを変化させた。凝集しやすいペプチドの$R_g$は約10.1 Åであったのに対し、スクランブルド・コントロールは約11.0 Åであった。
   • $\alpha$-シヌクレイン: 再配置は、主に局所的な疎水性（最大3残基平均疎水性がコントロールの約2.6から約4.2へ増加）を変化させたが、SCD、$R_g$、および非球形性は同程度であった。

意義および主張
本論文は、アミノ酸組成とは独立して、残基の順序が初期ペプチド会合の根本的な決定因子であることを確立している。本研究は、同一の組成であっても、配列依存的な分子組織化を通じて、異なるエネルギー的、構造的、および静電的な特性を生み出し得ることを示している。

主な主張は以下の通りである：

• 決定因子としての配列: 残基の配置が、アミノ酸組成だけでなく、分子間相互作用のバランスを支配し、その結果として生じる熱力学および界面の幾何学に影響を与える。
• ペプチド特異的な応答: 残基の再配置による分子レベルの帰結は普遍的なものではない。$\alpha$-シヌクレインとA$\beta$は対照的な応答を示し、配列の変更がエネルギー的な有利さと界面パッキングの異なる組み合わせをもたらす。
• 補完的な記述子: 本研究は、会合エネルギーと界面接触メトリクスが、分子会合の異なる側面を捉えていることを強調している。より多くの界面接触が存在しても、必ずしも会合エネルギーが有利になるとは限らず、これは二つのペプチド系で見られた相反する傾向によって示されている。
• 静電の再編成: 同一のアミノミノ酸組成であっても、電荷パターニングと3次元構造の再編成により、異なる静電環境が生じ得る。これは、配列特異的な会合挙動の分子的な根拠を提供している。

本研究は、配列エンコードされた分子決定因子を特定するための定量的枠組みを提供しており、初期凝集の「コード」は、結果として生じるエネルギーおよび静電ランドスケープを規定する、残基の特定の線形順序の中に存在することを示唆している。