Predicting Binding Affinities for the Binding Domain of Hyperpolarization-Activated Cyclic Nucleotide-Gated Channel Isoforms Using Free-Energy Perturbation

Mediante simulaciones de dinámica molecular y el método de perturbación de energía libre, este estudio calculó las energías de unión absolutas del AMPc a los dominios de unión de nucleótidos cíclicos de las isoformas 1-4 de los canales HCN, proporcionando así una comprensión más profunda de las diferencias en la sensibilidad y activación de estos canales.

Lo esencial

Imagina que tu corazón y tu cerebro son como una orquesta gigante que necesita mantener un ritmo constante para funcionar bien. En esta orquesta, hay unos "directores de orquesta" moleculares llamados canales HCN. Su trabajo es encender y apagar la electricidad en las células para que el corazón lata y el cerebro piense al ritmo correcto.

Pero estos directores tienen un problema: a veces se vuelven demasiado estrictos y no dejan que la música fluya, lo que puede causar arritmias o problemas neurológicos. Para ayudarlos, existe un "mensajero" químico llamado cAMP (una pequeña molécula) que llega y les dice: "¡Relájense, abran las puertas y dejen pasar la música!".

El problema es que no todos los directores de orquesta (hay cuatro versiones diferentes llamadas HCN1, HCN2, HCN3 y HCN4) reaccionan igual a este mensajero. Algunos son muy sensibles y se abren con un susurro, mientras que otros son duros de pelar y necesitan un grito para moverse. Los científicos quieren saber por qué pasa esto, para poder diseñar medicamentos que solo arreglen al director que está fallando, sin molestar a los otros.

¿Qué hicieron los científicos en este estudio?

En lugar de solo mirar fotos estáticas de estos canales (como si fueran estatuas), los investigadores usaron una supercomputadora para crear una película de alta velocidad de cómo se mueven y actúan en tiempo real.

1. La Simulación (El Laboratorio Virtual):
   Imagina que tienes cuatro modelos de plastilina de estos canales (uno para cada versión). Los científicos metieron a cada uno en una "piscina" virtual llena de agua y sal (simulando el interior de una célula) y les lanzaron el mensajero cAMP. Luego, observaron durante un tiempo muy largo (en la vida real, esto equivale a microsegundos, pero en la computadora es mucho tiempo) para ver cómo el mensajero se sentaba en su silla y cómo el canal reaccionaba.

2. La Prueba de Fuerza (Perturbación de Energía Libre):
   Para medir exactamente qué tan fuerte es la unión entre el mensajero y el canal, usaron una técnica matemática muy sofisticada llamada Perturbación de Energía Libre (FEP).

   • La analogía: Imagina que intentas sacar un imán pegado a una nevera. A veces es fácil, a veces cuesta mucho. Los científicos hicieron un cálculo matemático para "despegar" virtualmente al mensajero del canal, paso a paso, y midieron cuánta "fuerza" (energía) se necesitaba para separarlos. Cuanta más fuerza se necesita para separarlos, más fuerte es la unión.

¿Qué descubrieron?

Sus resultados son como encontrar las claves de la cerradura para cada tipo de puerta:

• Los "Amigables" (HCN1 y HCN3): Estos dos canales se llevan muy bien con el mensajero cAMP. La unión es fuerte y estable. Es como si el mensajero y el canal fueran viejos amigos que se abrazan con fuerza. Esto significa que estos canales son muy sensibles y se activan fácilmente.
• El "Tímido" (HCN2): Este canal es el que tiene la unión más débil. Es como si el mensajero intentara saludarlo, pero el canal se aleja un poco. Necesita mucho más esfuerzo para activarse.
• El "Intermedio" (HCN4): Está en medio, ni tan fuerte como los primeros ni tan débil como el segundo.

El Secreto de la Química (Los Aminoácidos)

¿Por qué son diferentes si se ven casi iguales? Los científicos miraron los "ladrillos" que componen estos canales (llamados aminoácidos) y descubrieron dos personajes clave:

1. El Arginina (El Abrazo Fuerte): En los canales HCN1 y HCN3, hay un aminoácido llamado arginina que actúa como un imán muy fuerte. Se une al mensajero con un abrazo de oso.
2. El Glutamato (El Plan B): En los canales HCN2 y HCN4, ese "abrazo fuerte" no funciona igual. En su lugar, usan otro aminoácido llamado glutamato para intentar mantener la unión, pero no es tan efectivo como el primero.

Es como si en una fiesta, algunos invitados (HCN1/3) tuvieran una mano de hierro para agarrar a sus amigos, mientras que otros (HCN2/4) tuvieran una mano de goma que se estira y se suelta más fácil.

¿Por qué es importante esto?

Hasta ahora, solo tenemos un medicamento aprobado que afecta a estos canales (Ivabradine), pero es como un martillo: sirve para todo, pero no es preciso.

Este estudio es como tener un manual de instrucciones detallado para cada tipo de canal. Ahora, los científicos saben exactamente qué "ladrillos" (aminoácidos) hacen que un canal sea sensible y otro no. Esto les permite diseñar medicamentos de precisión (como llaves maestras) que puedan:

• Apagar un canal que está causando epilepsia sin afectar al corazón.
• O arreglar un canal del corazón sin afectar al cerebro.

En resumen, los científicos usaron supercomputadoras para entender la "química del abrazo" entre un mensajero y sus receptores, descubriendo que pequeños cambios en la forma de abrazar explican por qué nuestros órganos funcionan de manera diferente. ¡Y eso es un gran paso para curar enfermedades con menos efectos secundarios!

Resumen técnico

Resumen Técnico: Predicción de Afinidades de Unión para Dominios de Unión de Canales HCN Mediante Perturbación de Energía Libre

1. El Problema

Los canales activados por hiperpolarización y nucleótidos cíclicos (HCN) son canales iónicos de Na+/K+ regulados por voltaje y modulados por nucleótidos cíclicos, esenciales para la actividad eléctrica rítmica espontánea en el corazón y el cerebro. Su disfunción se vincula a arritmias, epilepsia y otros trastornos neurológicos.
A pesar de su importancia terapéutica, el desarrollo de fármacos selectivos para los diferentes isoformas de HCN (HCN1 a HCN4) se ha visto obstaculizado por la falta de comprensión detallada de:

• La dinámica conformacional de los canales.
• Los mecanismos precisos de activación dual (por voltaje y por el ligando cAMP).
• Las diferencias específicas en la sensibilidad y afinidad de unión al cAMP entre las distintas isoformas, lo cual es crucial para diseñar fármacos con selectividad isoform-específica.

2. Metodología

Los autores emplearon una combinación de simulaciones de dinámica molecular (MD) de todo átomo y cálculos de perturbación de energía libre (FEP) para determinar las energías libres de unión absolutas del cAMP a los dominios de unión de nucleótidos cíclicos (CNBD) de las isoformas HCN1-4.

• Sistemas de Simulación: Se construyeron cuatro sistemas basados en estructuras experimentales (Cryo-EM para HCN1, cristalografía de rayos X para HCN2 y HCN4) y un modelo de AlphaFold para HCN3. Se aislaron los CNBDs (monómeros) y se solvataron en condiciones fisiológicas (0.15 M NaCl, agua TIP3P).
• Dinámica Molecular (MD): Se realizaron simulaciones de 100 ns por réplica (10 réplicas por isoforma, total 1 µs por isoforma) utilizando el campo de fuerzas CHARMM36m y NAMD 2.14. Se aplicaron restricciones iniciales para estabilizar el ligando antes de su liberación.
• Selección de Réplicas: Se filtraron las réplicas para el análisis FEP basándose en criterios estrictos de estabilidad: bajo desplazamiento del centro de masa del ligando (< 2 Å), baja RMSD del ligando (< 3 Å), distancia de "tapa" (lid distance) estable y formación de enlaces de hidrógeno estables en el fondo de la bolsa de unión.
• Perturbación de Energía Libre (FEP): Se calcularon las energías libres de unión relativas utilizando el método de Perturbación de Energía Libre (FEP) con el método de Bennett Acceptance Ratio (BAR). El proceso implicó la aniquilación alquímica del ligando en 20 ventanas, utilizando potenciales de núcleo blando para evitar problemas de extremos en las interacciones de van der Waals.
• Análisis por Residuo: Se calcularon las energías libres de interacción por residuo para identificar los contribuyentes clave a la unión.

3. Contribuciones Clave

• Cálculo de Energías de Unión Absolutas: Se determinaron por primera vez, mediante FEP, las energías libres de unión absolutas del cAMP a los CNBD de las cuatro isoformas de HCN en un marco comparativo directo.
• Análisis Conformacional Detallado: Se caracterizó la estabilidad de la "tapa" (lid) del dominio de unión y la dinámica de los hélices B y C, revelando variaciones en la estabilidad del dominio entre isoformas.
• Identificación de Residuos Críticos: Se cuantificó la contribución energética de residuos específicos (conservados y no conservados), destacando el papel de la arginina y el glutamato conservados en la modulación de la afinidad.
• Validación de Modelos: El estudio validó el uso de modelos de AlphaFold (para HCN3) combinados con MD y FEP para estudiar sistemas donde la estructura experimental completa no está disponible.

4. Resultados

• Jerarquía de Afinidad de Unión: Los resultados de FEP revelaron un orden claro de afinidad de unión al cAMP:
  1. HCN1 y HCN3: Presentan la mayor afinidad de unión comparable (ΔG ≈ -11.1 y -11.7 kcal/mol, respectivamente).
  2. HCN4: Afinidad intermedia (ΔG ≈ -8.9 kcal/mol).
  3. HCN2: La menor afinidad de unión (ΔG ≈ -6.9 kcal/mol).
• Mecanismos de Interacción:
  • HCN1 y HCN3: Muestran una fuerte interacción con un residuo de arginina conservado (offset +47), que actúa como un anclaje principal.
  • HCN2 y HCN4: Pierden esta interacción fuerte con la arginina. Como compensación, interactúan más fuertemente con un residuo de glutamato conservado (offset +38), aunque la energía total de unión es menor que en HCN1/3.
• Estabilidad del Dominio: Se observó que HCN2 y HCN4 mantienen una distancia de "tapa" y un ángulo entre hélices más estables que HCN1 y HCN3, lo que sugiere diferencias en la flexibilidad del dominio de unión que podrían influir en la cinética de activación.

5. Significado e Impacto

Este estudio proporciona una base molecular para entender por qué las isoformas de HCN difieren en su sensibilidad al cAMP, un factor determinante en su función fisiológica y patológica.

• Diseño de Fármacos: Los hallazgos ofrecen una guía racional para el diseño de fármacos selectivos. Al identificar que la arginina (offset +47) es crítica para la alta afinidad en HCN1/3, los investigadores pueden diseñar ligandos que exploten o eviten esta interacción para modular selectivamente isoformas específicas (por ejemplo, para tratar arritmias sin afectar la función neuronal).
• Avance Metodológico: Demuestra la viabilidad de combinar modelos predictivos (AlphaFold) con simulaciones avanzadas de energía libre para estudiar sistemas de canales iónicos complejos, superando la limitación de la falta de estructuras experimentales para todas las isoformas.
• Comprensión Mecanística: Refuta la idea de que las diferencias de sensibilidad se deben únicamente a la secuencia de aminoácidos, sugiriendo que la dinámica conformacional y las compensaciones energéticas entre residuos específicos juegan un papel igual de importante.