Swimming-limited aggregation of bacteria in liquid crystals

技術要約：液晶中における細菌の遊泳制限的凝集

問題提起
凝集と断片化は、エンジニアリングや自然界における基本的なプロセスであり、通常、拡散制限または反応制限のレジームにおける衝突成功率によって支配される。等方性流体中での細菌の凝集（枯渇引力や流れに誘起されたストリーマーによるバイオフィルム形成など）はよく知られているが、高度に異方的な環境における運動性細菌のダイナミクスについては、未だ十分に理解されていない。具体的には、ネマチック液晶（LC）中における運動性大腸菌（Escherichia coli）の集団ダイナミクスは、ネマチックディレクターによって細菌の運動が1次元（1D）に制約されるという、ユニークな系を提示している。先行研究では、液晶が弾性力を通じて受動的なコロイド凝集を誘起することや、運動性細菌が液晶中で鎖状構造を形成することが示されているが、これらの活性鎖の速度と安定性、特に鎖の長さと遊泳速度の関係を規定する物理的メカニズムについては、定量的な説明が欠けていた。

手法
本研究は、実験的観察と理論的モデリングおよび数値シミュレーションを組み合わせている：

1. 実験セットアップ： 運動性の大腸菌を、薄いガラスチャンバー（$h \approx 40 \pm 10$ $\mu$m）内に閉じ込められた生体適合性ネマチック液晶（ジソジウム クロモグリケート、DSCG）に浸した。液晶の配向はナノスクラッチによって規定され、細菌はディレクターに平行な1Dレーンに沿って泳ぐよう強制された。蛍光顕微鏡を用いて、混合から2〜30分後、5〜7分間にわたり、個々の細菌および鎖（最大16細胞まで観察したが、解析は $n \le 8$ に焦点を当てた）を追跡した。
2. 理論的モデリング：
   • 最小モデル： 中心となる細菌とその両隣の近接細胞に関する確率モデルを開発した。これは、相対速度に基づき、対面（エンカウンター）の確率と中央値となる時間を計算するものであり、対面確率の計算においては、対面時の速度分布を対数正規分布と仮定し、流体力学的および弾性的力を無視している。
   • エージェントベース・シミュレーション： 10,000のエージェントを用いた1D「ドライ」機械的結合モデルを実装した。このモデルでは、細菌は衝突するまで決定論的に泳行する。合体時、新しい鎖の速度は粘性抵抗の保存（総推力と総抵抗のバランス）によって決定され、その結果、鎖の速度は構成要素の速度の平均となる。モデルは初期速度の対数正規分布を取り入れ、正面衝突と追突の両方を許容している。
3. データ解析： 実験データ（鎖の長さ、速度）を、シミュレーションにおける特定の「動的年齢」（$t^*$、シミュレーション内の全鎖数が初期細菌数の45%に減少した時点と定義され、実験で観察された平均鎖長 $\langle n \rangle \approx 2.2$ と一致する）と比較した。

主要な結果

1. 活性凝集と受動的凝集のタイムスケール： 本研究は、活性凝集が分単位で行われることを示しており、これは液晶媒介の弾性引力のみによって駆動される受動的凝集に要する数時間よりも大幅に速い。これにより、細菌固有の運動性と相対速度の変動が、凝集プロセスを駆動する主要な要因であることが確認された。
2. 長さと速度の相関： 鎖の速度は構成要素の平均となる（成分がランダムにソートされている場合、長さには依存しない）という単純な力バランスモデルの予測に反して、実験では正の相関が見られた。すなわち、長い鎖ほど平均的に速く泳ぐ。
3. 自己ソーティングのメカニズム： 「遊泳制限的凝集（swimming-limited aggregation）」メカニズムがその原因として特定された。速い細菌は、近接者とより頻繁に遭遇し（中央値としてのエンカウンター時間が短い）、下流においてより長い鎖へと優先的に組み込まれる。この動的な自己ソーティング（self-sorting）により、長い鎖の中に速い泳ぎ手が濃縮される。
4. モデルの検証： 平均 $V \approx 9.18$ $\mu$m s$^{-1}$、標準偏差 $\sigma_V \approx 7.26$ $\mu$m s$^{-1}$ の対数正規速度分布を用いたエージェントベース・シミュレーションは、実験で見られた速度と長さの相関を再現することに成功した。また、モデルは異なる長さの鎖に対する速度分布の全形状も捉えている。
5. パラメータ感度： 速度と長さの相関は、主に遊泳速度の分散（$\sigma_V$）と方向のバイアスによって駆動される。ドラッグ係数、細胞長、または初期の空間配置の変化は、速度の変動と比較して、相関に対して無視できる程度の影響しか与えない。この相関は、システムが経年変化するにつれて弱まる。

意義と主張
本論文は、エージェントの固有の速度分布によって制約される「遊泳制限的凝集」という明確なレジームを確立している。著者らは以下のことを主張している：

• 観察された鎖の長さと速度の増加は、力のバランスへの違反ではなく、相対速度によって駆動される動的な自己ソーティングの結果である。
• 細胞間の遊泳速度の差異は、この活性的な自己組織化に不可欠である。速度の不均一性がなければ、ソーティング機構は機能しない。
• このメカニズムは、異方性媒体中における細菌の自己組織化を理解するための定量的な枠組みを提供し、1Dにおける運動誘起相分離がネマチック弾性によって強化・安定化され得ることを示唆している。
• これらの知見は、物理的な相互作用が個々の速度の差異を空間的および集団的な不均一性へと変換する、閉じ込められた環境における運動性に基づく自己組織化への一般的な経路を提示している。

本研究は、限界についても言及しており、「ドライ」モデルは流体力学的相互作用、速度の時間的変動、およびネマチック弾性を無視しており、これらが（速度と長さの曲線の正確な傾きなどの）定量的な詳細に影響を与える可能性があるとしている。しかし、活性凝集が速度依存的な遭遇率によって駆動され、それが自己ソーティングにつながるという核心的な発見は、提示された実験的および理論的枠組みの両方において強固である。