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

Resumen Técnico: Desentrañando los Mecanismos Moleculares de la Irisina en el Cáncer de Hígado

Planteamiento del Problema
El carcinoma hepatocelular (CHC) sigue siendo una de las principales causas de mortalidad relacionada con el cáncer a nivel mundial, diagnosticándose a menudo en estadios avanzados donde las opciones curativas son limitadas. Aunque se sabe que el ejercicio reduce el riesgo de la enfermedad y mejora la calidad de vida, los mecanismos moleculares específicos mediante los cuales las miocinas inducidas por el ejercicio, particularmente la irisina, modulan la patogénesis del CHC, permanecen insuficientemente caracterizados. La literatura existente presenta evidencia contradictoria sobre el papel de la irisina en el CHC: algunos estudios sugieren que potencia la proliferación y la migración a través de la vía PI3K/Akt, mientras que otros reportan ningún efecto o incluso una regulación a la baja en los tejidos tumorales. Además, no existe un consenso sobre los niveles de expresión de la irisina en los tejidos de CHC en comparación con los controles normales. Este estudio tiene como objetivo resolver estas inconsistencias identificando de manera exhaustiva los potenciales objetivos moleculares y las vías de señalización de la irisina en el CHC utilizando un enfoque de biología de sistemas.

Metodología
El estudio empleó un flujo de trabajo bioinformático de múltiples pasos que integra la toxicología de redes y el acoplamiento molecular:

1. Recolección de Objetivos:

   • Objetivos de la Irisina: Los genes diana humanos potencialmente modulados por la irisina se recuperaron de las bases de datos SwissTargetPrediction, ChEMBL y STITCH.
   • Objetivos del CHC: Los genes asociados con el cáncer de hígado se curaron de GeneCards y OMIM.
   • Expresión Diferencial: Se analizaron los datos de RNA-seq de la cohorte TCGA-LIHC (374 muestras tumorales y 50 normales) utilizando DESeq2 para identificar genes diferencialmente expresados (DEGs) (corte: adj. P < 0.05, |log₂FC| > 1).
   • Intersección: La intersección de los objetivos de la irisina, los genes asociados al CHC y los DEGs derivados de TCGA arrojó 93 objetivos potenciales.

2. Construcción de la Red e Identificación de Genes Centrales (Hub Genes):

   • Se construyó una red de interacción proteína-proteína (PPI) utilizando la base de datos STRING (puntuación de confianza > 0.4) y se visualizó en Cytoscape.
   • Análisis de Módulos: El plugin MCODE identificó la subred más crítica.
   • Selección de Hubs: Se utilizaron cinco algoritmos (Grado, Interintermediación, Cercanía, MNC, Radialidad) dentro del plugin cytoHubba para clasificar los genes. La intersección de los genes principales de estos algoritmos definió los "genes centrales" finales.

3. Modelado Pronóstico:

   • Se construyó un modelo de puntuación de riesgo basado en los niveles de expresión de los genes centrales identificados mediante la regresión de Cox multivariante.
   • Los pacientes fueron estratificados en grupos de alto y bajo riesgo. Se utilizaron el análisis de supervivencia de Kaplan-Meier y las curvas ROC para evaluar la precisión pronóstica.
   • La regresión de Cox multivariante evaluó la puntuación de riesgo como un factor pronóstico independiente junto con variables clínicas (edad, estadio, género).

4. Enriquecimiento Funcional y de Vías:

   • Se realizaron análisis de Ontología Génica (GO), KEGG, Reactome y Ontología de Enfermedades (DO) utilizando los paquetes clusterProfiler y ReactomePA para elucidar procesos biológicos y vías.

5. Validación por Acoplamiento Molecular:

   • Las estructuras proteicas para los genes centrales fueron predichas utilizando AlphaFold 3.
   • Se realizaron simulaciones de acoplamiento molecular (AutoDock Vina) para evaluar la afinidad de unión entre la irisina y las proteínas centrales. Se calcularon las energías de unión para evaluar la estabilidad de la interacción.

Resultados Clave

• Identificación de Objetivos: El estudio identificó 93 objetivos potenciales que vinculan la irisina con el CHC.
• Genes Centrales (Hub Genes): Se identificaron ocho genes centrales: RACGAP1, BIRC5, NEK2, UBE2C, TYMS, EZH2, AURKB y CDK1. Los ocho genes fueron significativamente sobreexpresados en los tejidos de CHC y asociados con un mal pronóstico (Hazard Ratio > 1). Aunque AURKA fue resaltado inicialmente en el análisis de intersección, fue excluido del conjunto final de genes centrales, mientras que los otros ocho fueron retenidos.
• Valor Pronóstico: La firma de riesgo basada en la irisina sirvió como un factor pronóstico independiente para pacientes con CHC, distinto de la edad, el estadio y el género. Los pacientes de alto riesgo exhibieron una supervivencia global significativamente peor (p < 0.001). El modelo mostró una capacidad predictiva moderada con valores de AUC de 0.76 (1 año), 0.66 (3 años) y 0.582 (5 años).
• Mecanismos Funcionales: El análisis de enriquecimiento reveló que los objetivos relacionados con la irisina están potencialmente involucrados en cuatro mecanismos de la inducción de CHC por la irisina:
  • Estabilidad Mitótica: Regulación de la unión del huso del microtúculo al cinetocoro y la segregación de las cromátidas hermanas.
  • Desregulación del Ciclo Celular: Participación en la transición de la fase G2/M y la degradación de proteínas del ciclo celular mediada por APC/C.
  • Actividad Oncogénica: Mejora de la actividad de la proteína serina/treonina quinasa y la unión a chaperonas.
  • Viabilidad Celular: Promoción del crecimiento tumoral y resistencia a agentes quimioterapéuticos (por ejemplo, resistencia al antifolato).
• Acoplamiento Molecular: Las ocho proteínas centrales demostraron fuertes afinidades de unión con la irisina. Siete proteínas exhibieron energías de unión inferiores a −7.2 kcal/mol, mientras que RACGAP1 mostró una energía de unión de −6.8 kcal/mol, lo que indica una robusta estabilidad de unión y una alta afinidad por todos los objetivos.

Significancia y Reivindicaciones
Los autores afirman que este estudio proporciona una perspectiva a nivel de sistemas sobre el papel de la irisina en el CHC, yendo más allá de los informes inconsistentes de genes individuales hacia una comprensión basada en redes. Las contribuciones primarias son:

1. Identificación de Objetivos Centrales: El estudio descubre ocho genes centrales específicos (RACGAP1, BIRC5, NEK2, UBE2C, TYMS, EZH2, AURKB, CDK1) que son críticos en la progresión del CHC mediada por la irisina.
2. Información Mecanística: Identifica cuatro mecanismos potenciales: inestabilidad cromosómica/aneuploidía durante la mitosis, desregulación del ciclo celular, mejora de la actividad de proteínas oncogénicas y aumento de la viabilidad de las células cancerosas.
3. Herramienta Pronóstica: La firma de 8 genes basada en la irisina fue validada como un indicador pronóstico independiente para la supervivencia de los pacientes.
4. Potencial Terapéutico: Las fuertes afinidades de unión confirmadas por el acoplamiento molecular sugieren que estas proteínas centrales son objetivos viables para futuras estrategias terapéuticas.

Limitaciones y Modestia
Los autores reconocen explícitamente que este es un estudio computacional. Señalan que:

• Los conjuntos de datos de objetivos de diversas bases de datos pueden sufrir de heterogeneidad debido a los diferentes algoritmos y estándares de predicción.
• Las interacciones de unión son predichas in silico y requieren validación experimental para confirmar la afinidad real.
• No todos los mecanismos identificados han sido verificados experimentalmente en el contexto del CHC inducido por la irisina.
• El estudio hace un llamado a una rigurosa validación experimental para elucidar completamente el papel preciso de la irisina en la patogénesis del CHC.

En conclusión, el artículo postula que la irisina potencialmente regula la progresión del CHC a través de vías mitóticas y del ciclo celular específicas mediadas por ocho genes centrales, ofreciendo un nuevo marco para comprender la hepatocarcinogénesis e identificar posibles objetivos para la prevención y el tratamiento, a la espera de la confirmación experimental.