A nonlinear theory for chemotactic fronts of mixed populations

Résumé technique : Une théorie non linéaire pour les fronts chimiotactiques de populations mixtes

Énoncé du problème
La migration collective de populations cellulaires hétérogènes est fondamentale dans les processus biologiques tels que le développement, la réponse immunitaire et les métastases cancéreuses. Bien que des travaux expérimentaux et théoriques récents aient établi que les interactions asymétriques avec des gradients chimiques auto-générés façonnent la distribution spatiale de types cellulaires distincts, les principes régissant l'organisation spatiale robuste de ces populations hétérogènes restent mal compris. Les études théoriques antérieures se sont largement concentrées sur des populations homogènes ou ont reposé sur des modèles linéarisés qui capturent la vitesse du front, mais échouent à caractériser pleinement la structure spatiale du front migrant. Plus précisément, il existe un manque de compréhension globale de la manière dont les interactions chimiques non réciproques (où un type de cellule consomme un chimioattractant et un autre le détecte) façonnent la co-migration de populations de cellules eucaryotes.

Méthodologie
Les auteurs s'appuient sur le modèle capteur/consommateur (SC) proposé par Ucar et al. [10], qui décrit la co-migration de cellules dendritiques (consommateurs, $\rho$) et de cellules T (capteurs, $\eta$) guidées par un chimioattractant auto-généré ($a$). Les équations directrices constituent un système de Keller–Segel modifié impliquant des termes de réaction, d'advection et de diffusion.

L'étude emploie une approche combinée de calibration quantitative et d'analyse asymptotique :

1. Calibration du modèle : Les auteurs calibrent les équations d'ondes progressives non dimensionnalisées par rapport aux données microfluidiques expérimentales concernant les densités spatiales des cellules dendritiques et des cellules T. En utilisant l'inférence bayésienne, ils déduisent les distributions postérieures des paramètres du modèle ($\chi_\rho, \kappa, D_\eta, D_\rho$) sans mesure directe du profil du chimioattractant.
2. Analyse asymptotique : Pour dériver une théorie non linéaire générale, les auteurs effectuent une analyse asymptotique systématique des équations différentielles ordinaires couplées décrivant l'onde progressive. Ils identifient des « limites asymptotiques distinguées » basées sur les magnitudes relatives de la sensibilité chimiotactique du consommateur ($\chi_\rho$) et de la diffusivité de la cellule consommatrice ($D_\rho$).
3. Classification des régimes : L'analyse catégorise le système en quatre régimes dynamiques distincts (A, B, C, D) basés sur ces paramètres. Pour les régimes où un équilibre physique unique ne domine pas (C et D), les auteurs utilisent des développements asymptotiques appariés pour résoudre les couches limites et connecter les solutions internes et externes.
4. Quantification du mélange : L'organisation spatiale est quantifiée par une métrique de mélange ($\mu$), définie comme la probabilité qu'une cellule capteuse soit située en avant d'une cellule consommatrice, dérivée de la distribution cumulative des cellules capteuses par rapport à la densité des consommateurs.

Contributions clés et résultats

• Théorie analytique non linéaire : L'article dérive une théorie analytique non linéaire complète pour les collectifs cellulaires hétérogènes. Contrairement aux linéarisations précédentes, ce cadre capture la structure spatiale complète du front progressif, liant directement les paramètres du modèle aux profils de densité des deux populations cellulaires.
• Quatre régimes dynamiques distincts : L'analyse révèle quatre régations distinctes régissant la morphologie du front du consommateur :
  • Régime A ($\chi_\rho \ll 1, D_\rho \gg \chi_\rho$) : Caractérisé par des interfaces diffuses et une signalisation à longue portée. Ici, la densité du consommateur est proportionnelle au carré du signal ($\rho \sim S^2$).
  • Régime B ($\chi_\rho, D_\rho \gg 1$) : Caractérisé par des interfaces diffuses mais une signalisation à courte portée. La densité et le signal sont fortement couplés ($\rho \sim S$).
  • Régime C ($D_\rho \ll \chi_\rho \ll 1$) : Caractérisé par des interfaces nettes et une signalisation à longue portée (décroissance algébrique du signal en amont du front).
  • Régime D ($\chi_\rho \gg 1, D_\rho \ll 1$) : Caractérisé par des interfaces nettes et une signalisation à courte portée (décroissance exponentielle).
• Calibration aux données Dendritiques/Cellules T : Lorsqu'il est calibré aux données expérimentales, le système se trouve opérer dans le Régime A. L'analyse prédit que ce régime est défini par une signalisation intermédiaire à longue portée du chimioattractant générée par une forte consommation du chimioattractant par les cellules dendritiques, équilibrée par une faible force chimiotactique.
• Mécanismes de robustesse : La théorie démontre que le système dendritique/cellule T opère dans un régime de paramètres qui équilibre le mélange intercellulaire avec la localisation des cellules T au front de tête. Les auteurs montrent que ce régime spécifique (Régime A) permet une population migrante mixte cohérente tout en permettant aux cellules capteuses de se localiser à l'avant.
• Impact de l'hétérogénéité : L'étude quantifie comment les variations de la diffusivité cellulaire ($D_\eta$) et de la sensibilité chimiotactique ($\kappa$) affectent la colocalisation. Dans les régimes diffusifs (A et B), le mélange est robuste aux perturbations. Dans les régimes à interface nette (C et D), le mélange est hautement sensible aux changements de paramètres, menant souvent à la séparation des populations de capteurs et de consommateurs.

Signification
L'article affirme fournir un cadre complet pour comprendre comment les interactions chimiques non réciproques façonnent une migration collective robuste dans les populations cellulaires hétérogènes. En allant au-delà des approximations linéaires, les auteurs démontrent que l'information spatiale encodée dans les profils de densité des populations mixtes est suffisante pour récupérer tous les paramètres du modèle, y compris le profil du chimioattractant non mesuré.

La signification réside dans l'identification de principes physiques spécifiques — à savoir l'équilibre entre la consommation du chimioattractant, la diffusion et la sensibilité chimiotactique — qui permettent une co-migration stable. La théorie souligne qu'une signalisation intermédiaire à longue portée (Régime A) est critique pour maintenir un collectif cohérent tout en permettant la localisation fonctionnelle de types cellulaires spécifiques (par exemple, les cellules T à l'avant). Ce cadre offre des principes généraux applicables à un large éventail de systèmes biologiques et synthétiques présentant des interactions chimiques asymétriques, fournissant une base mathématique pour prédire comment l'organisation spatiale propre à chaque trait émerge dans les collectifs migrants.