Morphogenesis driven by nematic defects in active biological networks

Résumé Technique : Morphogenèse pilotée par les défauts nématiques dans les réseaux biologiques actifs

Énoncé du Problème
L'article traite du défi fondamental de la modélisation de la morphogenèse cellulaire, plus précisément de la manière dont les tissus biologiques passent d'états désordonnés à des structures ordonnées aux formes spécifiques. Bien qu'il soit établi que les filaments du cytosquelette peuvent s'organiser en phases nématiques présentant un ordre orientationnel partiel, le rôle mécanique précis des défauts topologiques au sein de ces phases pour piloter la déformation tissulaire et la formation de motifs reste peu clair. Les modèles existants reposent souvent sur des approximations linéaires ou des gradients de morphogènes prédéfinis, qui peuvent échouer à capturer l'interaction complexe et non linéaire entre l'ordre nématique, l'élasticité et la croissance induite par la contrainte observée chez les organismes régénératifs comme Hydra. Les auteurs visent à développer un cadre continu traitant les tissus vivants comme des réseaux biologiques actifs (RBA) où les changements de forme émergent du couplage entre l'ordre nématique, l'élasticité et la relaxation adaptative, sans prescrire le motif de croissance a priori.

Méthodologie
Les auteurs proposent un cadre continu quasi-statique modélisant les tissus vivants comme des réseaux biologiques actifs composés de chaînes de filaments réticulés capables de former une phase nématique. Les hypothèses fondamentales et la formulation mathématique comprennent :

• Cinématique et Croissance : Le tissu est modélisé comme localement incompressible. La morphogenèse est décrite comme une séquence d'états d'équilibre. La croissance n'est pas prescrite mais est pilotée par la contrainte mécanique locale. Le tenseur de gradient de déformation $\mathbf{F}$ est décomposé de manière multiplicative en un tenseur de pas de longueur nématique $\mathbf{L}$, une distorsion élastique $\mathbf{F}_e$, un tenseur de croissance $\mathbf{F}_g$ et un tenseur de relaxation $\mathbf{F}_{rel}$.
• Formulation Énergétique : L'énergie libre totale $\mathcal{F}$ comprend deux contributions principales :
  1. Énergie Nématique ($\mathcal{F}_n$) : Basée sur la théorie continue d'Ericksen, ce terme rend compte des variations spatiales du champ de directeur $\mathbf{n}$ et du paramètre d'ordre scalaire $s$. Crucialement, $s$ est autorisé à varier spatialement, s'annulant aux défauts topologiques pour régulariser les singularités et permettre des transitions locales vers une phase isotrope.
  2. Énergie Élastique ($\mathcal{F}_e$) : Une fonction d'énergie de déformation néo-classique étendue pour inclure la croissance et la relaxation. Elle dépend de la partie élastique de la déformation, capturant la réponse du réseau de fibres aux distorsions.
• Mécanismes de Couplage : Le modèle intègre quatre mécanismes : (a) la contrainte mécanique à proximité des défauts topologiques, (b) l'ordonnancement induit par l'activité (augmentation de $s$), (c) la croissance induite par la contrainte (où du nouveau matériau est ajouté proportionnellement à la densité d'énergie élastique), et (d) la relaxation structurelle (mise à jour de la configuration sans contrainte pour relâcher les contraintes internes).
• Implémentation Numérique : Le problème non linéaire couplé est résolu à l'aide d'une formulation d'éléments finis mixtes implémentée dans FEniCSx. Pour gérer la symétrie tête-queue des directeurs nématiques et les coupures de branche associées aux défauts de demi-entier, le champ de directeur est reformulé à l'aide d'un champ de tenseur symétrique de rang 2 $\mathbf{M} = \mathbf{n} \otimes \mathbf{n}$. Les simulations utilisent une méthode de Newton avec recherche linéaire et chargement incrémentiel pour suivre l'évolution d'un sphéroïde isotrope initial vers des morphologies matures.

Contributions Clés

1. Cadre Continu pour les RBA : L'article introduit un modèle mathématique qui couple l'ordre nématique, l'élasticité et la croissance/relaxation induites par la contrainte, allant au-delà des théories de coquilles linéaires pour capturer des résultats morphogénétiques non linéaires.
2. Morphogenèse Pilotée par les Défauts : L'étude démontre que la topologie des défauts nématiques (spécifiquement leur charge et leur arrangement spatial) agit comme le principal moteur de l'évolution de la forme. Le modèle montre que les défauts ne sont pas de simples singularités mais des centres organisateurs qui localisent la contrainte et dirigent la croissance.
3. Reproduction de Morphologies Complexes : En faisant varier la configuration initiale des défauts, le modèle reproduit avec succès trois résultats morphogénétiques distincts observés lors de la régénération d'Hydra :
   • Régénération Uniaxiale (Tête-Pied) : Deux défauts $+1$ aux pôles opposés pilotent l'élongation et la formation d'une bouche (rupture apicale) et d'un pied.
   • Développement Bicephale (Deux Têtes) : Deux défauts $+1$ dans l'hémisphère supérieur, équilibrés par des défauts $-1/2$, conduisent à la formation de deux têtes distinctes.
   • Formation de Tentacules : Une configuration avec des défauts $+1$ aux extrémités des tentacules et des défauts $-1/2$ à leurs bases pilote la croissance des tentacules.
4. Mécanisme de Croissance Induit par la Contrainte : Les simulations montrent explicitement que la croissance est localisée dans les régions de forte contrainte élastique (près des défauts), tandis que la relaxation permet l'élargissement des ouvertures (bouches) là où la contrainte peut être géométriquement relâchée, contrastant avec la contrainte persistante au pied ou aux extrémités des tentacules.

Résultats
Les simulations par éléments finis de la régénération d'Hydra révèlent que :

• Localisation de la Contrainte : Les distorsions du directeur nématique près des défauts topologiques génèrent des contraintes mécaniques élevées. Ces contraintes pilotent l'accumulation de nouveau matériau (croissance) à proximité de ces derniers.
• Résultats Asymétriques : Dans la configuration uniaxiale, la relaxation de la contrainte permet au pôle apical (bouche) de se rompre et de s'élargir, tandis que le pôle basal (pied) reste fermé et structurellement contraint.
• Rôles des Défauts : Les défauts $+1$ pilotent la formation de protrusions et la rupture tissulaire. Les défauts $-1/2$ agissent principalement comme des stabilisateurs structurels ou des « équilibreurs topologiques » lors des phases précoces, bien que dans la configuration tentaculaire, l'accumulation de contrainte à ces défauts finisse par déclencher une croissance localisée à la base des tentacules.
• Évolution Morphologique : Le modèle capture la transition d'un fragment de tissu sphérique et isotrope vers des structures complexes et allongées dotées de caractéristiques spécifiques (pores oraux, tentacules) uniquement par l'évolution du paramètre d'ordre nématique et la contrainte de croissance résultante.

Signification et Revendications
Les auteurs affirment que leur travail fournit une base théorique robuste pour comprendre comment les défauts topologiques dans les réseaux nématiques actifs peuvent piloter des processus morphogénétiques complexes sans nécessiter de modèles de croissance pré-spécifiés ou de gradients de morphogènes externes. En traitant le tissu comme un réseau biologique actif où la croissance est une réponse directe à la contrainte mécanique générée par les distorsions nématiques, le modèle comble le fossé entre l'organisation microscopique du cytosquelette et la forme macroscopique du tissu. Les résultats concordent avec les observations expérimentales chez Hydra, suggérant que la « mémoire » de la polarité de l'axe corporel et l'émergence de structures spécifiques sont encodées dans la charge topologique et l'arrangement des défauts nématiques. L'article positionne ce cadre comme une étape vers la compréhension des moteurs mécaniques intrinsèques de la formation de motifs dans les tissus vivants, distincts des modèles purement biochimiques de signalisation.