Morphogenesis driven by nematic defects in active biological networks

Resumen Técnico: Morfogénesis Impulsada por Defectos Nemáticos en Redes Biológicas Activas

Planteamiento del Problema
El artículo aborda el desafío fundamental de modelar la morfogénesis celular, específicamente cómo los tejidos biológicos transicionan de estados desordenados a estructuras ordenadas con formas específicas. Si bien se ha establecido que los filamentos del citoesqueleto pueden organizarse en fases nemáticas con un orden orientacional parcial, el papel mecánico preciso de los defectos topológicos dentro de estas fases para impulsar la deformación tisular y la formación de patrones sigue sin estar claro. Los modelos existentes suelen depender de aproximaciones lineales o gradientes de morfógenos predefinidos, los cuales pueden fallar al capturar la compleja interacción no lineal entre el orden nemático, la elasticidad y el crecimiento impulsado por el estrés observada en organismos regenerativos como Hydra. Los autores pretenden desarrollar un marco continuo que trate a los tejidos vivos como redes biológicas activas (ABN, por sus siglas en inglés) donde los cambios de forma emergen del acoplamiento entre el orden nemático, la elasticidad y la relajación adaptativa, sin prescribir el patrón de crecimiento a priori.

Metodología
Los autores proponen un marco continuo cuasi-estático que modela los tejidos vivos como redes biológicas activas compuestas de cadenas de filamentos reticulados capaces de formar una fase nemática. Los supuestos centrales y la formulación matemática incluyen:

• Cinemática y Crecimiento: El tejido se modela como localmente incompresible. La morfogénesis se describe como una secuencia de estados de equilibrio. El crecimiento no es prescrito, sino que es impulsado por el estrés mecánico local. El gradiente de deformación $\mathbf{F}$ se descompone multiplicativamente en un tensor de paso de longitud nemática $\mathbf{L}$, una distorsión elástica $\mathbf{F}_e$, un tensor de crecimiento $\mathbf{F}_g$ y un tensor de relajación $\mathbf{F}_{rel}$.
• Formulación Energética: La energía libre total $\mathcal{F}$ comprende dos contribuciones principales:
  1. Energía Nemática ($\mathcal{F}_n$): Basada en la teoría continua de Ericksen, este término da cuenta de las variaciones espaciales en el campo del director $\mathbf{n}$ y el parámetro de orden escalar $s$. Crucialmente, $s$ puede variar espacialmente, desapareciendo en los defectos topológicos para regularizar las singularidades y permitir transiciones locales a una fase isotrópica.
  2. Energía Elástica ($\mathcal{F}_e$): Una función de energía de deformación neo-clásica extendida para incluir crecimiento y relajación. Depende de la parte elástica de la deformación, capturando la respuesta de la red de fibras ante las distorsiones.
• Mecanismos de Acoplamiento: El modelo integra cuatro mecanismos: (a) estrés mecánico cerca de los defectos topológicos, (b) orden inducido por actividad (aumento de $s$), (c) crecimiento impulsado por el estrés (donde se añade nuevo material proporcional a la densidad de energía elástica) y (d) relajación estructural (actualización de la configuración libre de estrés para aliviar los estreses internos).
• Implementación Numérica: El problema no lineal acoplado se resuelve utilizando una formulación de elementos finitos mixtos implementada en FEniCSx. Para manejar la simetría cabeza-cola de los directores nemáticos y los asociados cortes de rama para defectos de medio entero, el campo del director se reformula utilizando un campo de tensor de rango 2 simétrico $\mathbf{M} = \mathbf{n} \otimes \mathbf{n}$. Las simulaciones emplean un método de Newton con búsqueda de línea y carga incremental para rastrear la evolución desde un esferoide isotrópico inicial hacia morfologías maduras.

Contribuciones Clave

1. Marco Continuo para ABNs: El artículo introduce un modelo matemático que acopla el orden nemático, la elasticidad y el crecimiento/relajación impulsados por el estrés, yendo más allá de las teorías de cáscara lineales para capturar resultados morfogenéticos no lineales.
2. Morfogénesis Impulsada por Defectos: El estudio demuestra que la topología de los defectos nemáticos (específicamente su carga y disposición espacial) actúa como el motor primario de la evolución de la forma. El modelo muestra que los defectos no son meramente singularidades, sino centros organizadores que localizan el estrés y dirigen el crecimiento.
3. Reproducción de Morfologías Complejas: Al variar la configuración inicial de defectos, el modelo reproduce con éxito tres resultados morfogenéticos distintos observados en la regeneración de Hydra:
   • Regeneración Uniaxial (Cabeza-Pie): Dos defectos $+1$ en polos opuestos impulsan la elongación y la formación de una boca (ruptura apical) y un pie.
   • Desarrollo Bicéfalo (Dos Cabezas): Dos defectos $+1$ en el hemisferio superior, equilibrados por defectos $-1/2$, conducen a la formación de dos cabezas distintas.
   • Formación de Tentáculos: Una configuración con defectos $+1$ en las puntas de los tentáculos y defectos $-1/2$ en sus bases impulsa el crecimiento de los tentáculos.
4. Mecanismo de Crecimiento Impulsado por Estrés: Las simulaciones muestran explícitamente que el crecimiento se localiza en regiones de alto estrés elástico (cerca de los defectos), mientras que la relajación permite la ampliación de aperturas (bocas) donde el estrés puede aliviarse geométricamente, contrastando con el estrés persistente en el pie o las puntas de los tentáculos.

Resultados
Las simulaciones de elementos finitos de la regeneración de Hydra revelan que:

• Localización del Estrés: Las distorsiones en el campo del director nemático cerca de los defectos topológicos generan altos estreses mecánicos. Estos estreses impulsan la acumulación de nuevo material (crecimiento) en su vecindad.
• Resultados Asimétricos: En la configuración uniaxial, la relajación del estrés permite que el polo apical (boca) se rompa y se amplíe, mientras que el polo basal (pie) permanece cerrado y estructuralmente estresado.
• Roles de los Defectos: Los defectos $+1$ impulsan la formación de protuberancias y la ruptura del tejido. Los defectos $-1/2$ actúan principalmente como estabilizadores estructurales o "equilibradores topológicos" en etapas tempranas, aunque en la configuración de tentáculos, la acumulación de estrés en estos defectos eventualmente desencadena el crecimiento localizado en las bases de los tentáculos.
• Evolución Morfológica: El modelo captura la transición de un fragmento de tejido esférico e isotrópico hacia estructuras elongadas complejas con características específicas (poros orales, tentáculos) únicamente a través de la evolución del parámetro de orden nemático y el resultante crecimiento impulsado por el estrés.

Significancia y Reivindicaciones
Los autores afirman que su trabajo proporciona una base teórica robusta para comprender cómo los defectos topológicos en redes nemáticas activas pueden impulsar procesos morfogenéticos complejos sin la necesidad de patrones de crecimiento preespecificados o gradientes de morfógenos externos. Al tratar el tejido como una red biológica activa donde el crecimiento es una respuesta directa al estrés mecánico generado por las distorsiones nemáticas, el modelo cierra la brecha entre la organización microscópica del citoesqueleto y la forma macroscópica del tejido. Los resultados se alinean con las observaciones experimentales en Hydra, sugiriendo que la "memoria" de la polaridad del eje corporal y la emergencia de estructuras específicas están codificadas en la carga topológica y la disposición de los defectos nemáticos. El artículo posiciona este marco como un paso hacia la comprensión de los impulsores mecánicos intrínsecos de la formación de patrones en tejidos vivos, distintos de los modelos puramente bioquímicos de señalización.