Unraveling the molecular mechanisms of irisin in liver cancer through network toxicology and molecular docking strategies

Résumé Technique : Décryptage des mécanismes moléculaires de l'irisine dans le cancer du foie

Énoncé du problème
Le carcinome hépatocellulaire (CHC) demeure une cause majeure de mortalité liée au cancer à l'échelle mondiale, étant souvent diagnostiqué à des stades avancés où les options curatives sont limitées. Bien qu'il soit reconnu que l'exercice réduit les risques de maladie et améliore la qualité de vie, les mécanismes moléculaires spécifiques par lesquels les myokines induites par l'exercice, particulièrement l'irisine, modulent la pathogenèse du CHC restent mal caractérisés. La littérature existante présente des preuves contradictoires concernant le rôle de l'irisine dans le CHC : certaines études suggèrent qu'elle favorise la prolifération et la migration via la voie PI3K/Akt, tandis que d'autres ne rapportent aucun effet ou même une régulation à la baisse dans les tissus tumoraux. De plus, il n'existe pas de consensus sur les niveaux d'expression de l'irisine dans les tissus de CHC par rapport aux témoins normaux. Cette étude vise à résoudre ces incohérences en identifiant de manière exhaustive les cibles moléculaires et les voies de signalisation potentielles de l'irisine dans le CHC en utilisant une approche de biologie des systèmes.

Méthodologie
L'étude a employé un pipeline bioinformatique multi-étapes intégrant la toxicologie de réseau et le docking moléculaire :

1. Collecte des cibles :

   • Cibles de l'irisine : Les gènes cibles humains potentiels modulés par l'irisine ont été récupérés à partir des bases de données SwissTargetPrediction, ChEMBL et STITCH.
   • Cibles du CHC : Les gènes associés au cancer du foie ont été extraits de GeneCards et OMIM.
   • Expression différentielle : Les données d'ARN-seq de la cohorte TCGA-LIHC (374 échantillons tumoraux et 50 échantillons normaux) ont été analysées via DESeq2 pour identifier les gènes différentiellement exprimés (DEG) (seuil : adj. P < 0,05, |log₂FC| > 1).
   • Intersection : L'intersection des cibles de l'irisine, des gènes associés au CHC et des DEG dérivés de TCGA a produit 93 cibles potentielles.

2. Construction du réseau et identification des gènes pivots (hub genes) :

   • Un réseau d'interaction protéine-protéine (PPI) a été construit à l'aide de la base de données STRING (score de confiance > 0,4) et visualisé dans Cytoscape.
   • Analyse de module : Le plugin MCODE a identifié le sous-réseau le plus critique.
   • Sélection des hubs : Cinq algorithmes (Degree, Betweenness, Closeness, MNC, Radiality) au sein du plugin cytoHubba ont été utilisés pour classer les gènes. L'intersection des meilleurs gènes de ces algorithmes a défini les « gènes pivots » finaux.

3. Modélisation pronostique :

   • Un modèle de score de risque a été construit sur la base des niveaux d'expression des gènes pivots identifiés en utilisant la régression de Cox multivariée.
   • Les patients ont été stratifiés en groupes à haut risque et à faible risque. Des analyses de survie de Kaplan-Meier et des courbes ROC ont été utilisées pour évaluer la précision pronostique.
   • La régression de Cox multivariée a évalué le score de risque comme un facteur pronostique indépendant aux côtés des variables cliniques (âge, stade, sexe).

4. Enrichissement fonctionnel et de voies :

   • Des analyses d'ontologie génique (GO), KEGG, Reactome et d'ontologie des maladies (DO) ont été réalisées à l'aide des packages clusterProfiler et ReactomePA pour élucider les processus biologiques et les voies.

5. Validation par docking moléculaire :

   • Les structures protéiques pour les gènes pivots ont été prédites via AlphaFold 3.
   • Des simulations de docking moléculaire (AutoDock Vina) ont été menées pour évaluer l'affinité de liaison entre l'irisine et les protéines pivots. Les énergies de liaison ont été calculées pour évaluer la stabilité de l'interaction.

Résultats clés

• Identification des cibles : L'étude a identifié 93 cibles potentielles reliant l'irisine au CHC.
• Gènes pivots : Huit gènes pivots centraux ont été identifiés : RACGAP1, BIRC5, NEK2, UBE2C, TYMS, EZH2, AURKB et CDK1. Ces huit gènes ont tous été significativement surexprimés dans les tissus de CHC et associés à un mauvais pronostic (Hazard Ratio > 1). Bien qu'AURKA ait été initialement mise en évidence dans l'analyse d'intersection, elle a été exclue de l'ensemble final des gènes pivots, tandis que les huit autres ont été conservées.
• Valeur pronostique : La signature de risque basée sur l'irisine a servi d'indicateur pronostique indépendant pour les patients atteints de CHC, distinct de l'âge, du stade et du sexe. Les patients à haut risque présentaient une survie globale significativement plus faible (p < 0,001). Le modèle a montré une capacité prédictive modérée avec des valeurs AUC de 0,76 (1 an), 0,66 (3 ans) et 0,582 (5 ans).
• Mécanismes fonctionnels : L'analyse d'enrichissement a révélé que les cibles liées à l'irisine sont potentiellement impliquées dans quatre mécanismes de l'effet de l'irisine sur le CHC :
  • Stabilité mitotique : Régulation de l'attachement des microtubules du fuseau aux kinétochores et de la ségrégation des chromatides sœurs.
  • Dysrégulation du cycle cellulaire : Implication dans la transition G2/M et la dégradation médiée par l'APC/C des protéines du cycle cellulaire.
  • Activité oncogénique : Amélioration de l'activité des protéines sérine/thréonine kinase et de la liaison des chaperons.
  • Viabilité cellulaire : Promotion de la croissance tumorale et résistance aux agents chimiothérapeutiques (ex: résistance aux antifolates).
• Docking moléculaire : Les huit protéines pivots ont démontré de fortes affinités de liaison avec l'irisine. Sept protéines présentaient des énergies de liaison inférieures à −7,2 kcal/mol, tandis que RACGAP1 présentait une énergie de liaison de −6,8 kcal/mol, indiquant une stabilité de liaison robuste et une haute affinité pour tous les cibles.

Signification et Revendications
Les auteurs affirment que cette étude fournit une perspective systémique sur le rôle de l'irisine dans le CHC, dépassant les rapports incohérents sur des gènes uniques pour proposer une compréhension basée sur les réseaux. Les principales contributions sont :

1. Identification des cibles centrales : L'étude révèle huit gènes pivots spécifiques (RACGAP1, BIRC5, NEK2, UBE2C, TYMS, EZH2, AURKB, CDK1) qui sont critiques dans la progression du CHC médiée par l'irisine.
2. Aperçu mécanistique : Elle identifie quatre mécanismes potentiels : l'instabilité chromosomique/aneuploïdie pendant la mitose, la dysrégulation du cycle cellulaire, l'amélioration de l'activité des protéines oncogéniques et l'augmentation de la viabilité des cellules cancéreuses.
3. Outil pronostique : La signature à 8 gènes liée à l'irisine est validée comme un indicateur pronostique indépendant de la survie des patients.
4. Potentiel thérapeutique : Les fortes affinités de liaison confirmées par le docking moléculaire suggèrent que ces protéines pivots sont des cibles viables pour de futures stratégies thérapeutiques.

Limites et Modestie
Les auteurs reconnaissent explicitement qu'il s'agit d'une étude computationnelle. Ils notent que :

• Les ensembles de données de cibles provenant de diverses bases de données peuvent souffrir d'hétérogénéité en raison de différences d'algorithmes et de standards de prédiction.
• Les interactions de liaison sont prédites in silico et nécessitent une validation expérimentale pour confirmer l'affinité réelle.
• Tous les mécanismes identifiés n'ont pas été vérifiés expérimentalement dans le contexte de l'induction du CHC par l'irisine.
• L'étude appelle à une validation expérimentale rigoureuse pour élucider pleinement le rôle précis de l'irisine dans la pathogenèse du CHC.

En conclusion, l'article postule que l'irisine régule potentiellement la progression du CHC via des voies mitotiques et du cycle cellulaire spécifiques médiées par huit gènes pivots, offrant un nouveau cadre pour comprendre l'hépatocarcinogenèse et identifier des cibles potentielles pour la prévention et le traitement, sous réserve de confirmation expérimentale.