Ion-Modulated Polyelectrolyte Complexation of DNA and Polyacrylic Acid from Molecular Dynamics Simulations

技術要約：DNAとポリアクリル酸のイオン変調による高分子電解質複合体形成

問題提起
同電荷を持つ高分子電解質間の複合体形成は、同様の電荷を持つ高分子間に反発が生じると予測する従来の静電的な直感に挑戦するものである。多価陽イオンによるDNAの凝縮（例：多価陽イオンによるDNAの凝縮）のように、カウンターイオンの相関や電荷のパッチ状分布が引力を媒介するメカニズムは既知であるが、化学的に性質の異なる同電荷の高分子間の会合については、依然として理解が進んでいない。特に、半剛直なバイオポリマー（DNA）と柔軟なアニオン性合成ポリマー（ポリアクリル酸、PAA）の界面におけるイオンの原子レベルでの組織化については、体系的な検討が行われてこなかった。本研究は、水溶液中において、カウンターイオンの原子価とイオンサイズが、二つの負電荷を持つポリマー間の実効的な引力および架橋をどのように変調するかという理解の空白を埋めるものである。

手法
著者らは、3つの異なるイオン環境（単価のNaCl、および二価のMgCl₂およびCaCl₂）におけるDNA–PAA相互作用を調査するために、原子レベルの分子動力学（MD）シミュレーションを用いた。

• システム構築： シミュレーションでは、24塩基対のB型DNA二重らせん（純電荷 −46 e⁻）と、50%の脱プロトン化を受けた40モノマーの単一PAA鎖（純電荷 −−20 e⁻）を用いた。これらのシステムは、生理的イオン強度0.15 MのTIP3P水モデル中に溶媒和させた。
• 力場： DNAにはAMBER parmBSC1力場を用い、PAAにはGeneralized Amber Force Field (GAFF) を用いた。イオンパラメータは、TIP3P水に対して最適化されたJoung and Cheathamのセットから取得した。
• シミュレーションプロトコル： 生成ランは、300 K、1 barのNPTアンサンブルで行われた。DNA単独のシステムは500 ns、DNA–PAA複合体は1 μsのシミュレーションを行い、統計的な信頼性を確保するために各塩条件下で3つの独立したレプリカを実行した。
• 解析フレームワーク：
  • 状態分類： 重原子の距離カットオフ（≤ 0.6 nm）および少なくとも2つの水素結合の存在に基づき、「結合（bound）」状態を定義した。これは、持続的な複合体形成と一時的な遭遇を区別するために、100 ns以上の持続性を必要とする。
  • 構造指標： 結合状態および非結合状態に応じて、回転半径（Rg）、慣性主モーメント、および非球形度（asphericity）を含む、平方平均二乗偏差（RMSD）を計算した。
  • 相互作用解析： 短距離相互作用エネルギー（Lennard-JonesおよびCoulombic）を分解した。イオン媒介接触は、[DNA–ion–PAA]のトリプレットを0.35 nmのカットオフ内で特定することによって定量化した。また、DNAのリン酸基およびPAAのカルボキシレート周囲のイオンの動径分布関数（RDF）を算出した。

主要な結果
シミュレーションにより、イオンの原子価と水和特性に依存した、結合強度と構造的組織化の明確な階層性が明らかになった。

1. 複合体の安定性と結合率：

   • CaCl₂： 最も安定した複合体を生成し、平均結合率（$f_{bound}$）は0.89 ± 0.06であった。DNA–PAA複合体は、シミュレーション期間のほぼ全域にわたって結合状態を維持した。
   • MgCl₂： 中間的な安定性（$f_{bound}$ = 0.51 ± 0.21）を示し、より動的で断続的な解離が見られた。
   • NaCl： 最も高い変動性（$f_{bound}$ = 0.53 ± 0.28）を示し、結合状態と非結合状態の間を頻繁に行き来し、一部のレプリカでは長い解離期間が見られた。

2. イオンの配位と架橋メカニズム：

   • Ca²⁺： 持続的な分子架橋剤として機能した。RDFは、DNAのリン酸骨格およびPAAのカルボキシレートの両方に対する強い内圏配位（0.25–0.30 nm）を示した。Ca²⁺は平均2.17個の架橋を形成し、結合時間の約72%において少なくとも1つの架橋を維持した。
   • Mg²⁺： 弱く一時的な架橋（平均1.08個の架橋、結合時間の約44%）を示した。その強固に結合した水和シェルは外圏配位を優先させ、安定したバックボーンへのロックよりも、より動的なグルーブ関連の相互作用をもたらした。
   • Na⁺： 主に静電遮蔽を提供し、架橋形成は無視できる程度であった（結合時間の約10%）。NaClにおける会合は、直接的なイオン架橋ではなく、遮蔽と弱いファンデルワールス接触によって駆動されていた。

3. 構造応答：

   • DNA： CaCl₂において、結合したDNAは非結合セグメント（Rg ≈ 2.5 nm）やDNA単独の場合（Rg ≈ 2.5 nm）と比較して、より大きな回転半径（Rg ≈ 3.0 nm）を示し、強力な架橋と整列によるDNAフレームワークの拡張を示した。対照的に、MgCl₂およびNaClでは、状態依存的なDNAの構造変化は最小限であった。
   • PAA： CaCl₂において、PAAは複合体内で収縮し、より凝縮したコンフォメーションを採用した。
   • 配向： Ca²⁺は明確な方向的整列を誘導し、DNAの長軸とPAAの短軸が安定した低角の構成を採用した（$\cos\theta$ 分布のピーク）。Mg²⁺およびNa⁺は、より柔軟で秩序の低い会合を示した。

4. 結合モード：

   • Ca²⁺： 持続的なバックボーンへの結合が支配的であり、架橋はほぼ排他的にリン酸骨格に局在していた。
   • Mg²⁺： 主要溝（major groove）および副溝（minor groove）内での一時的な相互作用を促進し、PAAがCa²⁺系よりもDNA表面へ深く浸透することを可能にした。
   • Na⁺： 接触はまばらであり、主にバックボーンに限定され、安定した架橋は見られなかった。

意義と主張
本論文は、イオン特有の性質（原子価、水和シェルの不安定性、および移動度）がいかにして、同電荷のポリマー電解質の自己組織化を制御するかを理解するための、分子レベルの枠組みを提供すると主張している。著者らは、多価イオンは互換性があるものではないことを強調している。すなわち、不安定な水和シェルと柔軟な配位を持つCa²⁺は、真の架橋剤として機能し、複合体を安定化させ、ポリマーのトポロジーを再編成するが、Mg²⁺は一時的な相互作用を通じて中程度の安定化を提供し、Na⁺は遮蔽に依存するというものである。

本研究は、このような動的な系を正確に特徴付けるためには、結合状態と非結合状態を分離することが不可ло欠であることを強調している。単純な時間平均の指標では、強力な塩依存性の挙動を見落とす可能性があるためである。さらに、結果は、強力な局所的架橋（Ca²⁺に見られるような）が、パラドキシカルにDNA鎖のグローバルなサイズ（拡張）を増加させつつ、関連するポリマーを収縮させることを示しており、これは複合体を単なる凝縮としてのみ解釈する単純な見方に疑問を投げかけるものである。これらの知見は、多価イオン環境における遺伝子デリバリー、バイオセンシング、および応答性バイオマテリアルの設計に関連する、イオン応答性の核酸–ポリマー集合体の理解のための基礎を提供する。