Identification of the novel inhibitors against M. tuberculosis ESX-1 secretion system EccA1 enzyme using virtual screening, docking and dynamics simulation techniques

Este estudio identifica cinco compuestos de la base de datos ZINC como potenciales inhibidores antivirulencia de la enzima EccA1 del sistema de secreción ESX-1 de *M. tuberculosis*, mediante un enfoque computacional que incluyó modelado estructural, cribado virtual, análisis de propiedades farmacocinéticas y simulaciones de dinámica molecular que confirmaron su estabilidad y afinidad de unión.

Lo esencial

Imagina que la tuberculosis (Mycobacterium tuberculosis) es como un ladrón muy astuto que se mete en tu casa (tu cuerpo) y se esconde en una habitación especial (los macrófagos, que son las células de defensa). Para sobrevivir y robar recursos, este ladrón necesita abrir una caja fuerte secreta llamada sistema de secreción ESX-1. Dentro de esa caja fuerte hay una pieza clave, una especie de "motor" o "bomba" llamada enzima EccA1.

Si logramos apagar ese motor, el ladrón no puede abrir la caja, no puede esconderse bien y, finalmente, muere o es expulsado. El problema es que los ladrones actuales (las bacterias resistentes a los medicamentos) son muy fuertes y las llaves antiguas (los antibióticos tradicionales) ya no les funcionan.

Aquí es donde entra este estudio, que es como una búsqueda de nuevas llaves maestras usando una computadora muy potente.

1. El Mapa del Tesoro (Modelado de la Enzima)

Los científicos primero crearon un mapa digital 3D de la enzima EccA1. Imagina que es como construir una réplica exacta de la cerradura de la caja fuerte del ladrón en un videojuego. Saben que esta cerradura tiene dos partes: una que sirve de soporte y otra que es el motor (donde entra la energía, el ATP).

2. La Búsqueda de Llaves (Virtual Screening)

En lugar de fabricar millones de llaves físicas y probarlas una por una (lo cual tardaría años), usaron una base de datos gigante llamada ZINC, que es como un almacén digital con millones de llaves potenciales (compuestos químicos).

• La prueba: Lanzaron estas llaves virtuales contra la cerradura digital (la enzima) para ver cuáles encajaban mejor.
• El hallazgo: De entre millones, encontraron 5 llaves especiales (llamadas Z1, Z2, Z3, Z4 y Z5) que encajaban de forma increíblemente perfecta, mucho mejor que la llave original que usa la bacteria (el ADP) o incluso mejor que otras llaves de prueba que ya conocemos (como la CB5083 y la NMS873, que son medicamentos para el cáncer que se están probando).

3. La Prueba de Fuego (Simulación de Dinámica)

Encontrar una llave que encaja en un dibujo no es suficiente; hay que ver si funciona cuando la cerradura se mueve.

• La analogía: Imagina que pones la llave en la cerradura y luego agitas la puerta violentamente durante 100 segundos (en el mundo de la computadora, esto es una simulación de 100 nanosegundos).
• El resultado: Las 5 llaves nuevas (Z1-Z5) se mantuvieron firmes. No se salieron, no se rompieron. La cerradura se estabilizó. En cambio, algunas de las llaves de prueba antiguas (CB5083 y NMS873) se tambalearon y no encajaron tan bien en esta cerradura específica.
• Conclusión: Las nuevas llaves son muy estables y probablemente bloquearán el motor de la bacteria de manera efectiva.

4. ¿Son seguras para el humano? (Propiedades ADMET)

Antes de dar una llave a alguien, hay que asegurarse de que no sea tóxica.

• Los científicos revisaron si estas 5 llaves nuevas podrían ser absorbidas por el cuerpo humano, si el hígado las procesaría bien y si cumplirían con las reglas básicas para ser un medicamento oral (como el famoso "Regla de los 5" de Lipinski, que es como una lista de verificación de seguridad).
• Veredicto: ¡Sí! Cumplen con todos los requisitos. Parecen ser "buenos candidatos" para convertirse en pastillas reales.

En resumen

Este estudio es como si un grupo de ingenieros hubiera diseñado 5 nuevos dispositivos de bloqueo digitales para detener a un ladrón microscópico.

1. Identificaron la cerradura exacta que usa el ladrón.
2. Encontraron 5 llaves virtuales que encajan mejor que cualquier otra conocida.
3. Probaron que estas llaves aguantan el movimiento y no se sueltan.
4. Verificaron que son seguras para el cuerpo humano.

¿Qué sigue?
Ahora, los científicos necesitan salir del mundo digital y entrar al laboratorio real. Tienen que probar estas 5 llaves en un plato de Petri (con bacterias reales) y luego en animales para confirmar que realmente matan a la tuberculosis. Si todo sale bien, podríamos tener una nueva clase de medicamentos para curar a las personas que ya no responden a los tratamientos actuales. Es una esperanza muy prometedora en la lucha contra una enfermedad antigua y resistente.

Resumen técnico

Título: Identificación de nuevos inhibidores contra la enzima EccA1 del sistema de secreción ESX-1 de M. tuberculosis mediante cribado virtual, acoplamiento y simulaciones de dinámica molecular.

1. El Problema

La tuberculosis (TB), causada por Mycobacterium tuberculosis, sigue siendo una amenaza global de salud, con un aumento alarmante de casos de TB multirresistente (MDR) y extremadamente resistente (XDR). Las terapias actuales enfrentan desafíos significativos debido a la resistencia a los fármacos. El sistema de secreción ESX-1 es esencial para la virulencia de la bacteria, permitiendo la secreción de factores de virulencia (como EspAC) y la supervivencia dentro del fagosoma del huésped.
La enzima EccA1 es un componente clave de este sistema. Es una ATPasa única de la bacteria (ausente en humanos), lo que la convierte en un objetivo terapéutico ideal para desarrollar fármacos antivirulencia. Sin embargo, carecía de inhibidores específicos diseñados para su sitio activo ATPasa. El objetivo de este estudio fue identificar y validar compuestos pequeños que puedan inhibir la función de EccA1, bloqueando así la virulencia bacteriana.

2. Metodología

Los autores emplearon un enfoque integral de descubrimiento de fármacos basado en la estructura (SBDD):

• Modelado Estructural: Se reconstruyó la estructura de la enzima EccA1 de longitud completa (573 residuos, ~62.4 kDa) utilizando el servidor AlphaFold. El modelo se validó mediante análisis de Ramachandran (PROCHECK), confirmando la presencia de un dominio N-terminal TPR y un dominio C-terminal ATPasa tipo CbxX/CfqX.
• Cribado Virtual y Acoplamiento Molecular:
  • Se realizó un cribado virtual de la base de datos ZINC (aprox. 7.5 millones de compuestos) contra la bolsa ATPasa de EccA1.
  • Se utilizaron los módulos Glide (modo XP) y Prime de la suite Schrödinger para el acoplamiento molecular.
  • Se seleccionaron los 5 mejores compuestos candidatos (denominados Z1 a Z5) basándose en sus puntuaciones de acoplamiento y energías de unión.
  • Para comparación, se incluyeron el sustrato natural ADP y dos inhibidores clínicos conocidos de la ATPasa p97: NMS873 y CB5083.
• Cálculo de Energía de Unión: Se utilizó el método MM-GBSA (Área de Superficie Generalizada Born basada en Mecánica Molecular) para calcular la energía libre de unión ($\Delta G_{bin}$) de los complejos proteína-ligando.
• Simulación de Dinámica Molecular (MD): Se realizaron simulaciones de 100 ns utilizando el módulo Desmond (Schrödinger) con el campo de fuerzas OPLS3. Se analizaron los complejos de EccA1 con:
  1. Los compuestos Z1-Z5.
  2. El sustrato ADP.
  3. Los fármacos NMS873 y CB5083.
  4. La proteína apó (sin ligando).
• Análisis de Propiedades Farmacocinéticas (ADMET): Se evaluaron las propiedades de absorción, distribución, metabolismo y excreción, así como el cumplimiento de la Regla de los Cinco de Lipinski, utilizando el módulo QikProp.

3. Contribuciones Clave

• Identificación de Nuevos Inhibidores: Se descubrieron cinco compuestos de la base de datos ZINC (Z1-Z5) que muestran una afinidad de unión superior a la del sustrato natural ADP y a los fármacos de referencia (NMS873 y CB5083) contra EccA1.
• Validación Computacional Rigurosa: No solo se realizó el acoplamiento estático, sino que se validó la estabilidad de los complejos mediante simulaciones de dinámica molecular de larga duración (100 ns), analizando parámetros críticos como RMSD, RMSF, Radio de Giro (Rg) y Área Superficial Accesible al Solvente (SASA).
• Análisis Comparativo de Mecanismos: Se demostró que, a diferencia de los inhibidores de p97 (NMS873/CB5083) que no se adaptan bien a la bolsa de EccA1, los compuestos Z1-Z5 forman interacciones de hidrógeno estables y específicas con residuos clave del sitio activo (como la Motivo Walker A y B).
• Perfil Farmacocinético: Se estableció que los compuestos candidatos cumplen con los criterios de "droga-like" y tienen propiedades ADMET favorables para su desarrollo futuro.

4. Resultados Principales

• Energías de Unión: Los compuestos Z1-Z5 mostraron energías de unión muy favorables, oscilando entre -43.45 y -55.83 kcal/mol. En comparación, el ADP tuvo -35.00 kcal/mol, NMS873 -48.68 kcal/mol y CB5083 -50.88 kcal/mol. Esto indica que los nuevos compuestos tienen una mayor afinidad teórica que el sustrato natural.
• Interacciones Moleculares: Los compuestos Z1-Z5, que poseen estructuras similares a nucleósidos de purina, forman múltiples enlaces de hidrógeno con residuos críticos de la ATPasa (ej. Tyr466, Ile294, Thr341, Arg519, Arg522).
• Estabilidad en Dinámica Molecular (100 ns):
  • RMSD (Desviación Cuadrática Media): Los complejos EccA1-Z5 y EccA1-ADP mostraron los perfiles de RMSD más estables, indicando que el ligando se mantiene firmemente en la bolsa de unión. Los compuestos Z1-Z4 mostraron fluctuaciones moderadas pero estables.
  • RMSF (Fluctuación Cuadrática Media): La unión de los inhibidores redujo la flexibilidad de los residuos del sitio activo en comparación con la proteína apó, sugiriendo una estabilización conformacional.
  • SASA (Área Superficial Accesible): El complejo EccA1-Z5 mostró la mayor reducción en SASA, indicando la conformación más compacta y estable.
  • Enlaces de Hidrógeno: El complejo EccA1-ADP mantuvo el mayor número de enlaces de hidrógeno durante la simulación, seguido de cerca por Z5. Los fármacos NMS873 y CB5083 mostraron un perfil de unión más débil e inestable en EccA1, confirmando que no son inhibidores óptimos para este objetivo específico.
• Propiedades ADMET: Todos los compuestos Z1-Z5 cumplieron con la Regla de los Cinco de Lipinski, mostraron buena solubilidad y permeabilidad celular predicha, y no violaron las reglas de toxicidad estándar.

5. Significancia

Este estudio proporciona una base sólida para el desarrollo de nuevos fármacos antivirulencia contra la tuberculosis, específicamente dirigidos a la resistencia MDR/XDR.

• Especificidad: Al dirigirse a EccA1, una enzima exclusiva de M. tuberculosis, se minimiza el riesgo de toxicidad en humanos.
• Mecanismo de Acción: La inhibición de EccA1 bloquea la secreción de factores de virulencia, desactivando la capacidad de la bacteria para sobrevivir en el huésped sin necesariamente matarla, lo que podría reducir la presión selectiva para el desarrollo de resistencia.
• Avance Futuro: Los compuestos Z1-Z5, especialmente Z5 (que mostró la mayor estabilidad), se presentan como candidatos prometedores para ensayos in vitro e in vivo. El estudio demuestra la viabilidad de utilizar la dinámica molecular y el cribado virtual para rediseñar inhibidores de ATPasas bacterianas, ofreciendo una vía alternativa para combatir la TB resistente a los medicamentos.

En conclusión, la investigación valida computacionalmente a los compuestos Z1-Z5 como potentes inhibidores de EccA1, estableciendo un punto de partida crucial para la generación de la próxima generación de terapias antituberculosas.