Morphogenesis driven by nematic defects in active biological networks

Resumo Técnico: Morfogênese Impulsionada por Defeitos Nemáticos em Redes Biológicas Ativas

Declaração do Problema
O artigo aborda o desafio fundamental de modelar a morfogênese celular, especificamente como tecidos biológicos transitam de estados desordenados para estruturas ordenadas com formas específicas. Embora esteja estabelecido que os filamentos do citoesqueleto podem se organizar em fases nemáticas com ordem orientacional parcial, o papel mecânico preciso dos defeitos topológicos dentro dessas fases na condução da deformação tecidual e da formação de padrões permanece incerto. Os modelos existentes frequentemente dependem de aproximações lineares ou gradientes de morfógenos predefinidos, que podem falhar em capturar a complexa interação não linear entre ordem nemática, elasticidade e crescimento impulsionado por tensão observada em organismos em regeneração, como a Hydra. Os autores visam desenvolver um arcabouço contínuo que trate tecidos vivos como redes biológicas ativas (ABNs), onde mudanças de forma emergem do acoplamento entre ordem nemática, elasticidade e relaxação adaptativa, sem prescrever o padrão de crescimento a priori.

Metodologia
Os autores propõem um arcabouço contínuo quase-estático, modelando tecidos vivos como redes biológicas ativas compostas por cadeias de filamentos reticulados capazes de formar uma fase nemática. As principais suposições e a formulação matemática incluem:

• Cinemática e Crescimento: O tecido é modelado como localmente incompressível. A morfogênese é descrita como uma sequência de estados de equilíbrio. O crescimento não é prescrito, mas é impulsionado pelo estresse mecânico local. O gradiente de deformação $\mathbf{F}$ é decomposto multiplicativamente em um tensor de comprimento de passo nemático $\mathbf{L}$, uma distorção elástica $\mathbf{F}_e$, um tensor de crescimento $\mathbf{F}_g$ e um tensor de relaxação $\mathbf{F}_{rel}$.
• Formulação Energética: A energia livre total $\mathcal{F}$ compreende duas contribuições principais:
  1. Energia Nemática ($\mathcal{F}_n$): Baseada na teoria contínua de Ericksen, este termo contabiliza as variações espaciais no campo do diretor $\mathbf{n}$ e no parâmetro de ordem escalar $s$. Crucialmente, $s$ pode variar espacialmente, desaparecendo em defeitos topológicos para regularizar singularidades e permitir transições locais para uma fase isotrópica.
  2. Energia Elástica ($\mathcal{F}_e$): Uma função de energia de deformação neoclássica estendida para incluir crescimento e relaxação. Ela depende da parte elástica da deformação, capturando a resposta da rede de fibras às distorções.
• Mecanismos de Acoplamento: O modelo integra quatro mecanismos: (a) estresse mecânico próximo a defeitos topológicos, (b) ordenamento induzido pela atividade (aumento de $s$), (c) crescimento impulsionado pelo estresse (onde novo material é adicionado proporcionalmente à densidade de energia elástica) e (d) relaxação estrutural (atualizando a configuração livre de estresse para aliviar tensões internas).
• Implementação Numérica: O problema não linear acoplado é resolvido usando uma formulação de elementos finitos mistos implementada no FEniCSx. Para lidar com a simetria cabeça-cauda dos diretores nemáticos e as associações de cortes de ramo para defeitos de meio inteiro, o campo do diretor é reformulado usando um campo de tensor simétrico de posto 2 $\mathbf{M} = \mathbf{n} \otimes \mathbf{n}$. As simulações empregam um método de Newton com busca de linha e carregamento incremental para rastrear a evolução de um esferoide isotrópico inicial para morfologias maduras.

Principais Contribuições

1. Arcabouço Contínuo para ABNs: O artigo introduz um modelo matemático que acopla ordem nemática, elasticidade e crescimento/relaxação impulsionados por estresse, indo além das teorias de casca lineares para capturar resultados morfogenéticos não lineares.
2. Morfogênese Impulsionada por Defeitos: O estudo demonstra que a topologia dos defeitos nemáticos (especificamente sua carga e arranjo espacial) atua como o principal motor da evolução da forma. O modelo mostra que os defeitos não são meramente singularidades, mas centros organizadores que localizam o estresse e direcionam o crescimento.
3. Reprodução de Morfologias Complexas: Ao variar a configuração inicial de defeitos, o modelo reproduz com sucesso três resultados morfogenéticos distintos observados na regeneração de Hydra:
   • Regeneração Uniaxial (Cabeça-Pé): Dois defeitos $+1$ nos polos opostos impulsionam o alongamento e a formação de uma boca (ruptura apical) e um pé.
   • Desenvolvimento Bicefálico (Dois Cabeças): Dois defeitos $+1$ no hemisfério superior, equilibrados por defeitos $-1/2$, levam à formação de duas cabeças distintas.
   • Formação de Tentáculos: Uma configuração com defeitos $+1$ nas pontas dos tentáculos e defeitos $-1/2$ em suas bases impulsiona o crescimento dos tentáculos.
4. Mecanismo de Crescimento Impulsionado por Estresse: As simulações mostram explicitamente que o crescimento é localizado em regiões de alto estresse elástico (próximo aos defeitos), enquanto a relaxação permite o alargamento de aberturas (bocas) onde o estresse pode ser geometricamente aliviado, contrastando com o estresse persistente no pé ou nas pontas dos tentáculos.

Resultados
Simulações de elementos finitos de regeneração de Hydra revelam que:

• Localização de Estresse: Distorções no diretor nemático perto de defeitos topológicos geram altos estresses mecânicos. Esses estresses impulsionam o acúmulo de novo material (crescimento) em suas proximidades.
• Resultados Assimétricos: Na configuração uniaxial, a relaxação de estresse permite que o polo apical (boca) se rompa e se alargue, enquanto o polo basal (pé) permanece fechado e estruturalmente estressado.
• Papéis dos Defeitos: Defeitos $+1$ impulsionam a formação de protuberâncias e a ruptura do tecido. Defeitos $-1/2$ atuam principalmente como estabilizadores estruturais ou "balanceadores topológicos" em estágios iniciais, embora na configuração de tentáculos, o acúmulo de estresse nesses defeitos acabe por desencadear o crescimento localizado nas bases dos tentáculos.
• Evolução Morfológica: O modelo captura a transição de um fragmento de tecido esférico e isotrópico para estruturas complexas e alongadas com características específicas (poros orais, tentáculos) exclusivamente através da evolução do parâmetro de ordem nemática e do resultante crescimento impulsionado por estresse.

Significância e Alegações
Os autores afirmam que seu trabalho fornece uma base teórica robusta para entender como defeitos topológicos em redes nemáticas ativas podem impulsionar processos morfogenéticos complexos sem a necessidade de padrões de crescimento pré-especificados ou gradientes de morfógenos externos. Ao tratar o tecido como uma rede biológica ativa onde o crescimento é uma resposta direta ao estresse mecânico gerado por distorções nemáticas, o modelo preenche a lacuna entre a organização microscópica do citoesqueleto e a forma tecidual macroscópica. Os resultados alinham-se com observações experimentais em Hydra, sugerindo que a "memória" da polaridade do eixo corporal e a emergência de estruturas específicas estão codificadas na carga topológica e no arranjo dos defeitos nemáticos. O artigo posiciona este arcabouço como um passo em direção à compreensão dos drivers mecânicos intrínsecos da formação de padrões em tecidos vivos, distintos de modelos puramente bioquímicos de sinalização.