Mathematical Modeling of the Canonical Aryl Hydrocarbon Receptor Pathway

Cette étude développe et calibre un modèle mécaniste d'équations différentielles ordinaires de la voie canonique du récepteur d'hydrocarbure aromatique en utilisant des données d'expression génique résolues dans le temps provenant de ligands divers, révélant que les réponses transcriptionnelles spécifiques aux ligands sont principalement encodées au niveau de la régulation transcriptionnelle plutôt qu'au niveau des événements de signalisation en amont.

L'essentiel

Imaginez que votre corps possède un système de sécurité spécialisé appelé le récepteur des hydrocarbures aromatiques (AhR). Considérez ce récepteur comme une serrure intelligente située sur une porte à l'intérieur de vos cellules. Sa fonction est de détecter lorsque des substances chimiques étrangères (comme des polluants ou des toxines) tentent de pénétrer. Lorsque la bonne « clé » (un ligand chimique) s'insère dans cette serrure, la porte s'ouvre et la cellule déclenche une chaîne d'événements spécifique pour faire face à l'intrus.

Les scientifiques savent depuis longtemps que différentes clés (substances chimiques) s'adaptent à la serrure avec des degrés de serrage variés et que certaines clés se dégradent plus rapidement que d'autres. Cependant, un grand mystère demeurait : une fois la porte ouverte, quelle étape précise du processus est à l'origine de la réaction différente de la cellule face à différentes clés ? La différence provient-elle du fait que la clé tourne la serrure différemment, ou bien parce que les instructions écrites sur le mur à l'intérieur de la pièce sont différentes ?

Pour résoudre ce mystère, les chercheurs ont construit une simulation mathématique — essentiellement un « simulateur de vol » numérique pour ce système de sécurité cellulaire. Ils ont utilisé des données réelles provenant de cellules immunitaires de souris exposées à trois types de clés chimiques très différents :

1. Le 3-méthylcholanthrène (un activateur puissant et connu).
2. L'indolo[3,2-b]carbazole (un composé naturel).
3. Le bisphénol A (un produit chimique courant des plastiques).

Ils ne se sont pas contentés de faire des hypothèses ; ils ont testé 528 versions différentes de leur simulation. Dans chaque version, ils ont ajusté un ou deux « réglages de vitesse » spécifiques (vitesses de réaction) pour déterminer quelle combinaison correspondait le mieux aux données réelles qu'ils avaient recueillies auprès des cellules.

La grande découverte :
Après avoir exécuté toutes ces simulations, les chercheurs ont constaté que le « secret » ne réside pas dans les premières étapes du processus (comme la façon dont la clé tourne la serrure ou la vitesse à laquelle la clé disparaît). Au contraire, la réponse unique à chaque substance chimique est déterminée à l'intérieur de la salle de contrôle, spécifiquement au moment où la cellule décide de la quantité d'un nouveau message (ARNm) à produire.

Pensez-y comme à une usine :

• L'ancienne théorie : Peut-être que la différence de production dépend de la vitesse à laquelle les camions de livraison (signaux en amont) arrivent.
• La nouvelle découverte : Les camions de livraison arrivent tous à peu près à la même vitesse. La véritable différence réside dans la décision du directeur concernant le nombre de copies du manuel d'instructions à imprimer une fois les camions arrivés. Le « directeur » est la partie de la cellule qui se fixe sur l'ADN (occupation du promoteur) et décide de l'intensité du travail en fonction de la clé chimique utilisée.

En résumé : L'article conclut que la manière unique dont vos cellules réagissent à différentes substances chimiques est principalement contrôlée par la façon dont l'interrupteur de la cellule active les instructions génétiques, et non par les étapes antérieures du signal traversant la cellule.

Les chercheurs ont partagé leur modèle numérique (comme du code logiciel open source) afin que d'autres scientifiques puissent l'utiliser pour vérifier si cette même règle de « décision du directeur » s'applique à d'autres types de cellules dans l'organisme.

Résumé technique

Résumé technique : Modélisation mathématique de la voie canonique du récepteur des hydrocarbures aromatiques

Énoncé du problème
Le récepteur des hydrocarbures aromatiques (AhR) fonctionne comme un facteur de transcription activé par un ligand, essentiel pour la détection des xénobiotiques, le développement, l'immunité et l'homéostasie tissulaire. Bien qu'il soit établi que différents ligands présentent des affinités de liaison et des cinétiques de dégradation variables, et que ces différences subtiles influencent la signalisation en aval, les étapes de réaction biochimique spécifiques au sein de la voie canonique de l'AhR responsables des réponses transcriptionnelles distinctes propres à chaque ligand restent obscures. Le défi central abordé par cette étude consiste à identifier quels paramètres mécanistiques au sein de la voie codent la spécificité des ligands.

Méthodologie
Pour répondre à cette question, les auteurs ont développé un modèle mécanistique d'équations différentielles ordinaires (EDO) représentant la voie canonique de l'AhR. L'étude a utilisé des données de PCR quantitative (qPCR) résolues dans le temps mesurant les niveaux d'ARNm de Cyp1a1 et d'Ahrr dans des macrophages dérivés de moelle osseuse de souris. Ces cellules ont été exposées à trois ligands structurellement diversifiés, pertinents sur le plan environnemental et possédant une activité immunomodulatrice connue : le 3-méthylcholanthrène, l'indolo[3,2-b]carbazole et le bisphénol A.

Le processus d'étalonnage du modèle a fait appel à une optimisation globale sous une vraisemblance hétéroscédastique pour ajuster les données expérimentales. Afin de disséquer les sources de la spécificité des ligands, les auteurs ont évalué 528 modèles candidats. Ces modèles ont été construits en permettant à une ou deux constantes de vitesse de réaction impliquant le ligand de varier, testant ainsi différentes hypothèses concernant l'origine des effets spécifiques aux ligands au sein de la voie. La sélection de modèles a été réalisée à l'aide de critères basés sur Akaike.

Contributions et résultats clés
La contribution principale de ce travail réside dans le développement et l'étalonnage d'un cadre quantitatif reliant les paramètres biochimiques spécifiques aux ligands aux résultats transcriptionnels. L'analyse de sélection de modèles a révélé un régime dynamique dominant où les réponses transcriptionnelles spécifiques aux ligands sont principalement gouvernées par deux facteurs :

1. Occupation du promoteur : L'ampleur selon laquelle le complexe récepteur lié au ligand occupe le promoteur du gène cible.
2. Synthèse d'ARNm du gène cible : Le taux d'initiation de la transcription suite à la liaison au promoteur.

Les résultats indiquent que les réponses transcriptionnelles distinctes observées pour différents ligands sont codées au niveau de la régulation transcriptionnelle plutôt que d'être pilotées par des événements de signalisation en amont. Par conséquent, les variations des affinités de liaison et des cinétiques de dégradation seules sont insuffisantes pour expliquer l'ensemble du spectre des réponses spécifiques aux ligands sans prendre en compte ces étapes de régulation en aval.

Importance et affirmations
L'article affirme que son modèle mécanistique fournit une explication robuste de la manière dont la spécificité des ligands est atteinte au sein de la voie canonique de l'AhR, déplaçant l'accent des variations de signalisation en amont vers les mécanismes de régulation transcriptionnelle. Pour assurer la reproductibilité et faciliter la recherche ultérieure, le modèle résultant est disponible aux formats SBML et PEtab. Les auteurs déclarent que cette ressource permet de futures investigations sur la question de savoir si ces effets spécifiques aux ligands sont conservés ou réorganisés entre différents types cellulaires, sans formuler d'affirmations plus larges concernant des applications thérapeutiques spécifiques ou des propositions expérimentales non testées.