Structure-guided prioritization of Andrographis paniculata diterpenoids at the PCSK9-LDLR interface

Resumen Técnico: Priorización de diterpenoides de Andrographis paniculata guiada por la estructura en la interfaz PCSK9-LDLR

Planteamiento del problema
La interacción entre la proproteína convertasa subtilisina/kexina tipo 9 (PCSK9) y el receptor de lipoproteínas de baja densidad (LDLR) es una diana terapéutica validada para reducir el colesterol LDL. Sin embargo, la interrupción de esta interacción mediante moléculas pequeñas es estructuralmente desafiante. La interfaz de unión es amplia, expuesta al solvente y carece de un bolsillo profundo canónico, lo que puede hacer que las puntuaciones de acoplamiento (docking) estándar sean potencialmente engañosas, ya que pueden premiar la acomodación superficial inespecífica en lugar de una competencia productiva. Los diterpenoides de Andrographis paniculata han mostrado efectos moduladores de lípidos en etapas preclínicas, pero sus dianas moleculares directas aún no se han definido. Este estudio aborda la necesidad de una estrategia de priorización que distinga entre la afinidad nominal de acoplamiento y el compromiso estructural relevante de la interfaz.

Metodología
El estudio empleó un flujo de trabajo computacional por etapas y limitado por evidencia para evaluar seis diterpenoides de A. paniculata (L01–L06):

1. Preparación de ligando y receptor: Los ligandos se prepararon desde PubChem preservando la estereoquímica. El complejo PCSK9-LDLR (PDB 3BPS) sirvió como referencia estructural. Se definió prospectivamente un conjunto específico de 19 residuos de PCSK9 (153–446) involucrados en el reconocimiento de LDLR como la diana de la interfaz.
2. Acoplamiento con semillas triplicadas: Se utilizó AutoDock Vina para realizar el acoplamiento por triplicado para cada ligando utilizando diferentes semillas aleatorias para evaluar la convergencia. La clasificación de la afinidad se basó en las medianas de las puntuaciones, mientras que el mapeo de interacción se centró en el solapamiento con los residuos de la interfaz predefinidos.
3. Perfilado ADMET: Se realizaron predicciones multiplataforma (SwissADME, pkCSM, ProTox) para evaluar la semejanza con fármacos, la solubilidad, la absorción y la toxicidad.
4. Calificación del modelo del receptor: Se rechazó un flujo de trabajo inicial de cadena experimental debido a brechas en las coordenadas. Se generó un modelo comparativo de la secuencia de los residuos 153–446 de PCSK9 mediante SWISS-MODEL, validado mediante dinámica molecular (MD) apo de corta escala temporal.
5. Dinámica molecular emparejada: Los compuestos clasificados por afinidad (L05, isoandrographolide) y el compuesto clasificado por interfaz (L06, andrograpanin) fueron sometidos a simulaciones de MD no restringidas emparejadas de 1 ns. Las métricas incluyeron RMSD de la columna vertebral de la proteína, RMSF a nivel de residuo, RMSD del ligando, persistencia de puentes de hidrógeno y ocupación de contacto a nivel de residuo.

Resultos clave

• Afinidad de acoplamiento vs. Compromiso de la interfaz: Los seis compuestos exhibieron un rango estrecho de afinidad (−5.6 a −6.3 kcal mol⁻¹). L05 (isoandrographolide) tuvo la puntuación mediana más favorable (−6.3 kcal mol⁻¹). Sin embargo, el mapeo de contactos reveló una divergencia crítica: L05 se solapaba con un solo residuo de la interfaz predefinida (ILE154), mientras que L06 (andrograpanin) contactaba de forma reproducible con ocho residuos de la interfaz (SER153, PRO155, ASP238, ALA239, ILE369, PHE379, VAL380, SER381) que abarcan segmentos proximales, centrales y distales de la superficie de unión.
• Estabilidad dinámica y ocupación: En las simulaciones de 1 ns emparejadas, ambos complejos fueron técnicamente estables sin disrupción de la estructura global. No obstante, L06 demostró un compromiso dinámico superior:
  • Estabilidad del ligando: L06 mostró un RMSD de ligando alineado con la proteína menor (0.0961 nm) en comparación con L05 (0.1355 nm).
  • Puentes de hidrógeno: L06 mantuvo puentes de hidrógeno en el 58.4% de los fotogramas frente al 47.5% de L05.
  • Ocupación de residuos: L06 mantuvo una ocupación del 100% en PRO155, ALA239 e ILE369, y una alta ocupación (>94%) en SER153, PHE379 e ILE154. En contraste, la ocupación de L05 se localizó principalmente en ILE154 (85.1%), con un contacto insignificante en los residuos distales.
• Perfiles ADMET: L02 (neoandrographolide) fue despriorizado debido a la hepatotoxicidad y la responsabilidad hERG predichas, a pesar de tener una puntuación de acoplamiento favorable. L06 mostró una alta absorción intestinal predicha, pero una mayor lipofilicidad y menor solubilidad. L05 presentó un perfil moderado, pero fue señalado como un sustrato potencial de la P-glicoproteína.

Significancia y Reivindicaciones
El artículo sostiene que la amplitud de la interfaz y la persistencia dinámica proporcionan una base más informativa desde el punto de vista mecánico para la selección de candidatos que las pequeñas diferencias en la afinidad de acoplamiento por sí solas.

• Diferenciación de candidatos: El estudio distingue con éxito entre un candidato "clasificado por afinidad" (L05) y un candidato "clasificado por interfaz" (L06). Mientras que L05 se une fuertemente de manera localizada, L06 compromete la superficie funcional de forma más exhaustiva, imitando la topología requerida para potencialmente interferir con el reconocimiento PCSK9-LDLR.
• Priorización para pruebas experimentales: Los autores identifican a la andrograpanina (L06) como el andamiaje (scaffold) más fuerte dirigido por la interfaz para las pruebas experimentales. Esta priorización se basa en la convergencia de las poses de acoplamiento, la amplitud del contacto con los residuos y la retención dinámica de dichos contactos, en lugar de en los cálculos de energía libre de unión.
• Rigor metodológico: El estudio enfatiza la importancia de la definición prospectiva de la interfaz, el uso de semillas triplicadas para asegurar la reproducibilidad de la pose y el rechazo de modelos de receptores inválidos (debido a brechas en las coordenadas) para evitar artefactos de conectividad artificial.
• Alcance modesto: Los autores declaran explícitamente que estos hallazgos representan una priorización basada en la estructura, no una evidencia de inhibición de la interacción PCSK9-LDLR. Debido a los tiempos de simulación cortos (1 ns), la falta de cálculos de energía libre y la ausencia de validación bioquímica, los resultados se interpretan como generadores de hipótesis. El estudio no afirma que L06 sea un fármaco, sino que es un andamiaje líder digno de validación biofísica y funcional (por ejemplo, ensayos de competición, SPR, estudios de reciclaje en hepatocitos).