Swimming-limited aggregation of bacteria in liquid crystals

Technische Samenvatting: Zwem-beperkte Aggregatie van Bacteriën in Vloeibare Kristallen

Probleemstelling
Aggregatie en fragmentatie zijn fundamentele processen in de techniek en de natuur, die doorgaans worden beheerst door botsingssuccesratio's in diffusiebeperkte of reactiebeperkte regimes. Hoewel bacteriële aggregatie in isotrope vloeistoffen goed gedocumenteerd is (bijv. biofilmvorming gedreven door depletie-aantrekking of door stroming geïnduceerde streamers), blijft de dynamica van mobiele bacteriën in sterk anisotrope omgevingen minder goed begrepen. Specifiek presenteren de collectieve dynamica van mobiele Escherichia coli in nematische vloeibare kristallen (LC's) een uniek systeem waarbij de bacteriële beweging door de nematische director wordt beperkt tot één dimensie (1D). Eerdere studies hebben aangetoond dat LC's passieve colloïdale aggregatie kunnen induceren via elastische krachten en dat mobiele bacteriën ketens kunnen vormen in LC's, maar het fysieke mechanisme dat de snelheid en stabiliteit van deze actieve ketens bepaalt, met name de relatie tussen de ketenlengte en de zwemsnelheid, miste een kwantitatieve verklaring.

Methodologie
De studie combineert experimentele observatie met theoretische modellering en numerieke simulatie:

1. Experimentele Opstelling: Mobiele E. coli werden ondergedompeld in een biocompatibele nematische LC (disodiumcromoglycaat, DSCG) ingesloten in dunne glazen kamers ($h \approx 40 \pm 10$ $\mu$m). De LC-uitlijning werd voorgeschreven door nano-krassen, waardoor bacteriën langs 1D-banen parallel aan de director zwommen. Fluorescentiemicroscopie werd gebruikt om individuele bacteriën en ketens (tot 16 cellen, hoewel de analyse zich concentreerde op $n \le 8$) te volgen over intervallen van 5–7 minuten, 2–30 minuten na het mengen.
2. Theoretische Modellering:
   • Minimaal Model: Een probabilistisch model werd ontwikkeld voor een centrale bacterie en zijn twee dichtstbijzijnde buren. Het berekent de waarschijnlijkheid en de mediane ontmoetings-tijd op basis van relatieve snelheden, uitgaande van een lognormale verdeling van zwemsnelheden en waarbij hydrodynamische en elastische krachten worden genegeerd voor de berekening van de ontmoetingskans.
   • Agent-Based Simulaties: Een 1D "droge" mechanische koppelingsmodellen werd geïmplementeerd voor 10.000 agenten. In dit model zwemmen bacteriën deterministisch tot ze botsen. Bij het samensmelten wordt de snelheid van de nieuwe keten bepaald door de conservering van viskeuze weerstand (het balanceren van totale stuwkracht en totale weerstand), wat resulteert in een ketensnelheid die de gemiddelde snelheid van de constituerende componenten is. Het model incorporeert een lognormale verdeling van initiële snelheden en staat zowel frontale als achterwaartse botsingen toe.
3. Data-analyse: Experimentele gegevens (ketenlengtes, snelheden) werden vergeleken met de output van simulaties bij een specifieke "dynamische leeftijd" ($t^*$), gedefinieerd als de tijd waarop het totale aantal ketens in de simulatie daalt tot 45% van het initiële aantal bacteriën, wat overeenkomt met de gemiddelde ketenlengte ($\langle n \rangle \approx 2,2$) die in experimenten werd waargenomen.

Belangrijkste Resultaten

1. Actieve versus Passieve Aggregatietijdschalen: De studie toont aan dat actieve aggregatie plaatsvindt op de orde van grootte van minuten, wat aanzienlijk sneller is dan de uren die nodig zijn voor passieve aggregatie die uitsluitend wordt gedreven door LC-gemedieerde elastische aantrekking. Dit bevestigt dat intrinsieke bacteriële mobiliteit en variaties in relatieve snelheid de primaire drijfveren zijn van het aggregatieproces.
2. Lengte-Snelheid Correlatie: In tegenstelling tot eenvoudige krachtbalansmodellen die voorspellen dat de ketensnelheid de gemiddelde waarde van de componenten zou moeten zijn (en dus onafhankelijk van de lengte als componenten willekeurig gesorteerd zijn), onthullen experimenten een positieve correlatie: langere ketens zwemmen gemiddeld sneller.
3. Mechanisme van Zelfsortering: Het "zwem-beperkte aggregatiemechanisme" wordt geïdentificeerd als de oorzaak. Snellere bacteriën komen vaker tegen met buren (kortere mediane ontmoetingstijden) en worden bij voorkeur downstream opgenomen in langere ketens. Deze dynamische zelfsortering leidt tot een verrijking van snellere zwemmers in langere ketens.
4. Model Validatie: De agent-gebaseerde simulaties, gebruikmakend van een lognormale snelheidsverdeling met een gemiddelde $V \approx 9,18$ $\mu$m s$^{-1}$ en een standaarddeviatie $\sigma_V \approx 7,26$ $\mu$m s$^{-1}$, reproduceren succesvol de experimentele snelheid-lengte correlatie. Het model vangt ook de volledige vorm van de snelheidsverdelingen voor ketens van verschillende lengtes op.
5. Parametersensitiviteit: De snelheid-lengte correlatie wordt primair gedreven door de variantie in zwemsnelheden ($\sigma_V$) en directionele bias. Variaties in weerstandscoëfficiënten, cel-lengtes of initiële ruimtelijke arrangementen hebben een verwaarloosbaar effect op de correlatie vergeleken met de snelheidvariabiliteit. De correlatie verzwakt naarmate het systeem ouder wordt.

Betekenis en Claims
Het artikel vestigt een onderscheidend regime van "zwem-beperkte aggregatie", waarbij de evolutie van actieve materie wordt beperkt door de intrinsieke snelheidsverdeling van de agenten. De auteurs beweren dat:

• De waargenomen toename van de ketensnelheid met de lengte geen schending is van de krachtbalans, maar een gevolg van een dynamisch zelfsorteringsproces gedreven door relatieve snelheden.
• Variaties tussen cellen in zwemsnelheid essentieel zijn voor deze zelforganisatie; zonder heterogeniteit zou het sorteringsmechanisme falen.
• Dit mechanisme een kwantitatief kader biedt voor het begrijpen van bacteriële zelforganisatie in anisotrope media, wat suggereert dat motiliteit-geïnduceerde fase-scheiding in 1D kan worden versterkt en gestabiliseerd door nematische elasticiteit.
• Deze bevindingen bieden een generieke route voor motiliteit-gebaseerde zelforganisatie in beperkte omgevingen, waar fysieke interacties individuele snelheidvariaties omzetten in ruimtelijke en collectieve heterogeniteit.

De studie erkent beperkingen en merkt op dat het "droge" model hydrodynamische interacties, temporele snelheidsschommelingen en nematische elasticiteit verwaarloost, wat de kwantitatieve details (zoals de exacte helling van de snelheid-lengte curve) kan beïnvloeden. Echter, de kernbevinding — dat actieve aggregatie wordt gedreven door snelheid-afhankelijke ontmoetingsfrequenties die leiden tot zelfsortering — is robuust over de gepresenteerde experimentele en theoretische kaders heen.