How Spatially Modulated Activity Reshapes Active Polymer Conformations

技術要約：空間的に変調された活性が能動的ポリマーのコンフォメーションをどのように再形成するか

問題提起
内部的な力によって平衡状態から外れた能動的ポリマー（アクティブポリマー）は、受動的な鎖とは異なるコンフォメーション特性を示す。これまでの研究では、一定の力強度を持つ接線方向の能動的な力がポリマーにグロビュール（球状）のような相を誘起することが確立されているが、「空間的に変調された」活性の具体的な役割については、未だ十分に解明されていない。既存の解析モデルの多くは、均一な力分布を仮定しているか、あるいは系統的な解析的枠組みを欠いた数値シミュレーションに依存している。本研究は、接線方向の推進力の空間パターン（波長およびモード）が、半柔軟性ポリマーのコンフォメーション統計（特に回転半径 $R_G$ および末端間距離 $R_E$）にどのように影響するかという空白を埋めるものである。

手法
著者らは、解析的摂動論と数値シミュレーションを組み合わせた二角的なアプローチを採用している。

1. 解析モデル:

   • フレームワーク: 有限の曲げ剛性を持つ連続体ラウズモデル（Rouse model）を用いて、半柔軟性ポリマーを記述する。運動方程式には、輪郭 $s$ に沿って変化する接線方向の能動的力 $f(s)$ が含まれる。
   • 展開: 弱活性（小さなペクレ数 $Pe_m$）および小さな持続長（persistence length）の極限における系統的な展開を行う。モデルでは、ポリマーの座標と能動的な力を展開するために、受動的ラウズモデルの固有関数（余弦基底）を利用する。
   • 力プロファイル: 能動的な力は、単一モードの正弦波変調 $f(s) \propto \cos(\pi m s / N)$ としてモデル化される。ここで、$m$ は力のモードを表す整数である。
   • 導出: 著者らは、能動的な力の振幅の2次まで、モード相関 $\langle r_i \cdot r_j \rangle_\infty$ の解析的な表現を導出する。これらの相関は、$R_G^2$ および $R_E^2$ を計算するために用いられる。導出における重要な側面は、モード和の収束を保証するために小さな曲げ剛性を導入することであり、これにより、能動的な力を持つ柔軟なラウズモデルで見られる発散を防いでいる。

2. 数値シミュレーション:

   • モデル: 3次元における半柔軟なビーズ・スプリング鎖のランジュバン動力学シミュレーション。
   • 相互作用: 隣接するモノマーは、鎖の連結性と固定された結合長を維持するために、FENEおよびWCAポテンシャルを介して相互作用する。曲げ剛性は、余弦ベースのポテンシャル ($U_b(\theta) = \kappa(1 + \cos \theta)$) を通じて導入される。
   • 活性: 各ビーズの力強度が、解析モデルで定義された正弦波プロファイルに従うように、接線方向の自己推進が実装されている。
   • パラメータ: 理論的予測を様々なレジーム（有限のペクレ数を含む、摂動展開が形式的に制限される領域）において検証するため、ポリマー長 ($N$)、力のモード ($m$)、および能動的な力の強さ ($f_m$) を変化させてシミュレーションを行う。

主要な貢献と結果

• モード依存の転移: 本研究は、活性の空間モード $m$ によって駆動される、ポリマーのコンフォメーションにおける根本的な転移を特定している。

  • 低モード ($m \le 1$): 一様または低モードの強制力（長波長）は、ポリマーの収縮を誘起し、コンパクトなグロビュール様のコンフォメーションへと導く。これは、一定の能動的な力に関する先行研究と一致する。
  • 高モード ($m > 1$): 高モードの強制力は、鎖に沿って交互に伸びたセグメントと圧縮されたセグメントを生成する。これらの領域が十分に多くなると、全体的な効果としてポリマーの膨潤が生じ、伸長したコンフォメーションとなる。
  • 閾値: 臨界モード $m^*$ が、収縮レジームと膨張レジームを分かつ。解析モデルは $m^* = 4$ と予測するが、シミュレーションでは $m^* = 2$ が観察される。

• サイズ指標のデカップリング（分離）: 能動的システムにおける回転半径 ($R_G$) と末端間距離 ($R_E$) の標準的な関係の崩壊が、重要な発見として示されている。

  • $R_G$ の挙動: $R_G^2$ への主要な補正項は、能動的な力の振幅の2次式である。したがって、弱力極限において $R_G$ は力プロファイルの符号（方向）に依存せず、一般に低モードでは減少し、高モードでは増加する。
  • $R_E$ の挙動: $R_E^2$ への主要な補正項は、力のモードのパリティ（奇偶）に依存する。奇数モードの場合、補正項は力の振幅に対して線形であり、$R_E$ は力パターンの向き（例：力が端に向かって内向きか外向きか）に対して敏感である。これにより、$R_G$ ではコンパクトであるが $R_E$ では伸長しているという、受動的な平衡ポリマーでは不可能な現象が可能となる。

• 制御パラメータ ($Pe_m$): 著者らは、モード依存のペクレ数 $Pe_m = f_m b \sqrt{N} / k_B T$ を定義している。彼らは、このパラメータがポリマーサイズ $N$ に依存することを強調しており、これは、同じ局所的な力強度において、より長いポリマーの方が活性に対して敏感であることを意味する。このサイズ依存性は、摂動展開の妥当性が、固定された力ではなく、固定された小さな $Pe_m$ に対して制限されることを示唆している。

• 崩壊の性質: 本論文は、活性による収縮と、引力的相互作用による平衡崩壊を区別している。解析モデルにおいて、接線方向の活性による $R_G$ の減少は、力の強さに対して2次的であるが、保存的な引力的相互作用は線形依存性を示す。このことは、活性による収縮が、単純な有効引力とは異なる非平衡現象であることを示唆している。

意義と主張
本論文は、空間的に変化する活性を持つ能動的ポリマーを理解するための、体系的な解析的枠組みを提供すると主張している。この領域は、これまで数値的な観察に支配されてきた。主な意義は、活性の空間的なパターニングが、平衡から遠く離れたポリマーのコンフォメーションを制御する強力なメカニズムであることを示した点にある。

著者らは以下の点を強調している：

1. パターンが重要である: 能動的な力の空間波長は、ポリマーを収縮させるか膨張させるかを決定する上で、その大きさと同じくらい重要である。
2. 非普遍性: 能動的ポリマーは、単一の有効なサイズ指標では記述できない。$R_G$ と $R_E$ は、同じ駆動力に対して異なる応答を示す。特に、力の方向やパリティに対する感度の面で顕著である。
3. 堅牢性: 弱力での導出であるにもかかわらず、定性的な傾向（収縮 vs 膨張）および転移閾値の存在に関する解析的予測は、より強い活性で行われた数値シミュレーションに対しても成立している。

結論として、現在の理論は弱力レジームおよびガウス鎖に限定されているものの、より複雑なパターン、より強い駆動、およびトポロジー的に制約されたフィラメント（環状構造など）の研究のための基礎を築くものであるとしている。