Growing, Buckling, and Swirling: motility from polymerization

Resumen Técnico: Crecimiento, Pandeo y Remolino: Motilidad a partir de la Polimerización

Planteamiento del Problema
El artículo aborda los mecanismos físicos subyacentes a una clase distinta de micro-nadadores impulsados por el crecimiento de filamentos del citoesqueleto anclados a la superficie. Mientras que los modelos tradicionales de micro-natación (por ejemplo, rotación flagelar o batido ciliar) suelen prescribir condiciones de contorno o velocidades de deslizamiento superficial, este trabajo investiga cómo las fuerzas activas emergen de manera autoconsistente a partir de la cinética microscópica de la polimerización, la reorientación de los filamentos y las interacciones hidrodinámicas. El estudio está motivado por sistemas biológicos como Acetobacter xylinum, que extruye fibras de celulosa que se pandean para generar motilidad, y coloides sintéticos propulsados por la polimerización de la actina. El desafío central es formular un marco teórico que resuelva simultáneamente la cinética de los filamentos, la dinámica de orientación y las interacciones fluido-estructura en un entorno de bajo número de Reynolds (estocástico/Stokesiano), donde la locomoción resultante está determinada por la dinámica autoorganizada de la alfombra de filamentos en lugar de condiciones de contorno predefinidas.

Metodología
Los autores desarrollan un marco continuo autoconsistente que acopla la cinética de polimerización microscópica con la dinámica de fluidos macroscópica:

1. Cinética de Campo Medio: El sistema se modela como una densa alfombra de filamentos que experimentan nucleación, crecimiento a velocidad constante ($V_g$) y depolimerización estocástica (catástrofe). Una función de distribución $\Psi(\ell, t)$ describe las longitudes de los filamentos, mientras que una función de distribución generalizada $\psi(y, \ell, p, t)$ da cuenta de la posición espacial, la longitud y la orientación $p$.
2. Separación de Escalas de Tiempo: El modelo asume que la polimerización y la depolimerización ocurren en escalas de tiempo más rápidas que la reorientación de los filamentos. Esto permite la eliminación adiabática de la dinámica de longitud, reduciendo el problema a la evolución de una distribución de orientación $\tilde{g}(y, p, t)$ y un campo de polarización $n(y, t)$.
3. Dinámica de los Filamentos: Los filamentos individuales se modelan como varillas rígidas ancladas a una superficie mediante resortes torsionales (rigidez $k_\alpha$). Su reorientación está gobernada por la ecuación de Jeffery (cizalladura del fluido) y torques de restauración, equilibrados contra las fuerzas activas generadas por la polimerización contra el fluido viscoso.
4. Acoplamiento Hidrodinámico: Las fuerzas activas se promedian para obtener un salto de tracción efectivo que actúa sobre una superficie desplazada $\partial S$ (ubicada en la longitud media de los filamentos $\bar{\ell}$ por encima del cuerpo rígido). Esta tracción impulsa el flujo de Stokes. El cuerpo rígido se trata como libre de fuerzas y torques, con sus velocidades traslacional ($U$) y angular ($\Omega$) determinadas de forma autoconsistente.
5. Implementación Numérica:
   • Cierre: Se utiliza un cierre de Bingham polar para aproximar los momentos de orientación de orden superior ($\langle pp \rangle, \langle ppp \rangle$) en términos del vector de polarización $n$, evitando el costo computacional de resolver la ecuación completa de Fokker-Planck.
   • Solucionador: Las ecuaciones de Stokes se resuelven utilizando un Método de Integral de Frontera (BIM) con un esquema de cuadratura de quinto orden y Cuadratura por Expansión (QBX) para interacciones casi singulares. El sistema se resuelve iterativamente mediante GMRES para manejar el acoplamiento implícito entre las tracciones superficiales y los gradientes de velocidad.

Resultados Clave

• Análisis de Estabilidad Lineal: En una geometría de semi-espacio, los autores demuestran que una alfombra uniforme de filamentos que crecen verticalmente es inestable ante perturbaciones de longitud de onda larga. Las fuerzas de compresión inducidas por la polimerización provocan una inestabilidad de pandeo que rompe la simetría, llevando a una inclinación espontánea de los filamentos y a la generación de flujos de gran escala. La inestabilidad está gobernada por la competencia entre los estreses activos desestabilizadores y los torques elásticos de restauración y las tasas de depolimerización.
• Dinámica de Esferas Fijas: Para una esfera rígida con nucleación uniforme, la inestabilidad conduce a un patrón de remolino caracterizado por dos defectos topológicos $+1$ antipodales. Esta configuración genera flujos rotacionales coherentes en el exterior de la esfera, análogos a los flujos citoplasmáticos en ovocitos de Drosophila.
• Dinámica de Esferas Libres (Rotación Espontánea): Cuando la esfera es libre de moverse, el flujo de remolino que rompe la simetría genera un torque hidrodinámico neto, causando que la esfera experimente una rotación de cuerpo rígido sostenida alrededor de un eje seleccionado espontáneamente.
• Nucleación Patroneada y Motilidad: Al romper la simetría espacial de los sitios de nucleación, los autores demuestran el control sobre el equilibrio de fuerzas y torques:
  • Banda Ecuatorial: Genera un torque neto pero fuerza neta cero, resultando en un giro puro.
  • Parche Polar: Genera tanto fuerza neta como torque, resultando en trayectorias circulares (traslación y rotación combinadas).
  • Simetría de Cuatro Pliegues: Un arreglo de cuatro parches cancela el torque neto mientras suma las fuerzas, produciendo una traslación dirigida y persistente.
  • Patrones Quirales: Un arreglo helicoidal de sitios de nucleación acopla la rotación con la traslación, produciendo trayectorias rectas con giro simultáneo, imitando mecanismos de propulsión quirales encontrados en la naturaleza (por ejemplo, el esperma de Ginkgo).
• Papel de la Geometría: Al extender el análisis a esferoides prolatos, los autores revelan que la forma del cuerpo altera significativamente las trayectorias. Debido a la anisotropía de arrastre, un esferoide con un patrón quiral traza una hélice circular en lugar de una línea recta, ya que el desalineamiento entre la tracción neta y el eje principal causa una precesión continua de la dirección de propulsión.

Significancia y Reivindicaciones
El artículo establece una base teórica general para la motilidad impulsada por la dinámica colectiva de filamentos en polimerización. Demuestra que la polimerización misma es una ruta genérica hacia la autopropulsión, capaz de generar comportamientos diversos (giro, motilidad dirigida, natación quiral) puramente a través de la autoorganización de la orientación y el crecimiento de los filamentos.

Los autores afirman que su marco proporciona un nuevo principio de diseño para micro-nadadores sintéticos, donde el "código de superficie" (el patronamiento espacial de los sitios de nucleación) dicta la estrategia de locomoción emergente. El trabajo cierra la brecha entre la cinética de polimerización microscópica y la hidrodinámica macroscópica, ofreciando una descripción continua unificada aplicable tanto a alfombras activas impulsadas por polimerización como por motores. Además, los resultados sugieren que velocidades de natación del orden de $0.1\text{--}1 \, \mu\text{m/s}$ son alcanzables mediante estos mecanismos. El artículo concluye señalando que, si bien el modelo actual asume una aproximación de "alfombra delgada", se requiere trabajo futuro para abordar regímenes donde los filamentos crecen hasta longitudes que exceden el tamaño del cuerpo, requiriendo enfoques híbridos discretos-continuos.