⚛️ biophysics

Kinase inhibitors can change protonation or tautomeric state upon binding

Mediante el uso de simulaciones de Electroestática de Continuo de Conformaciones Múltiples (MCCE), este estudio demuestra que los inhibidores de quinasas experimentan frecuentemente cambios dinámicos en sus estados de protonación y tautomería al unirse, revelando que las especies minoritarias pueden volverse dominantes en el complejo unido y destacando la necesidad crítica de considerar los conjuntos completos de protonación y tautomería para una predicción precisa de la energética de unión.

Autores originales: Ranepura, G. A., Chowdhury, S. I., Rosenzweig, E. A., Rustenburg, A. S., Lopez-Rios de Castro, R., Mao, J., Chodera, J. D., Singh, S., Gunner, M. R.

Publicado 2026-07-30

📖 5 min de lectura🧠 Análisis profundo

CC BY 4.0

Autores originales: Ranepura, G. A., Chowdhury, S. I., Rosenzweig, E. A., Rustenburg, A. S., Lopez-Rios de Castro, R., Mao, J., Chodera, J. D., Singh, S., Gunner, M. R.

Artículo original bajo licencia CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ⚕️ Esta es una explicación generada por IA de un preprint que no ha sido revisado por pares. No es consejo médico. No tome decisiones de salud basándose en este contenido. Leer descargo de responsabilidad completo

Imagina que estás intentando encajar una llave en una cerradura muy compleja y de alta tecnología. En el mundo de la medicina, esa cerradura es una proteína dentro de tu cuerpo, y la llave es una molécula de fármaco diseñada para detener una enfermedad. Durante años, los científicos han intentado predecir exactamente qué tan bien encajan estas llaves observando sus formas. Pero hay un giro oculto: tanto la llave como la cerradura están hechas de diminutos bloques de construcción que pueden cambiar su personalidad eléctrica. Algunas partes pueden capturar un protón (un pequeño núcleo de hidrógeno) y volverse con carga positiva, mientras que otras pueden perder uno y volverse negativas. También pueden reorganizar sus átomos internos en diferentes formas "tautoméricas", como una persona cambiando su peinado.

Estos cambios eléctricos son cruciales porque las cargas opuestas se atraen y las cargas iguales se repelen. Si una molécula de fármaco cambia su carga o forma en el momento en que toca la proteína, todo el encaje cambia. Sin embargo, la mayoría de los programas informáticos utilizados para diseñar fármacos son como maniquíes rígidos; asumen que el fármaco y la proteína mantienen un estado eléctrico fijo, generalmente aquel que tienen cuando flotan en el agua. Esto es un poco como intentar predecir cómo se pegará un imán a una nevera sin darse cuenta de que el imán podría invertir sus polos en el segundo en que se acerque. Comprender estos cambios dinámicos es vital para crear medicinas mejores y más efectivas, pero es increíblemente difícil de calcular porque hay tantas combinaciones posibles de cargas y formas para revisar.

Este artículo se sumerge en este mundo desordenado y dinámico para ver qué sucede realmente cuando los fármacos inhibidores de quinasas (una clase importante de combatientes contra el cáncer y la inflamación) se unen a sus proteínas diana. Los investigadores utilizaron una herramienta de simulación informática sofisticada llamada MCCE (Electroestática de Continuo de Multiconformadores) para actuar como un observador súper rápido y súper preciso. En lugar de adivinar una única "mejor" forma para el fármaco o la proteína, dejaron que la computadora explorara miles de diferentes posibilidades de carga y forma a la vez, calculando cuáles son las más probables de ocurrir. Observaron 18 fármacos diferentes uniéndose a 9 proteínas quinasas diferentes, simulando el entorno tanto en el agua de la célula como en lo profundo del bolsillo de unión de la proteína.

Los resultados revelaron algunos giros sorprendentes en la historia de la unión de fármacos. Primero, la idea de que un fármaco siempre se vuelve más "calmado" o menos cargado cuando sale del agua y entra en la proteína fue demostrada errónea en muchos casos. Aunque tiene sentido que un objeto cargado pierda energía al salir del agua (como una toalla mojada secándose), el interior de la proteína es a menudo una fiesta cargada que en realidad quiere retener esa carga eléctrica. En varios casos, el fármaco se volvió más cargado o mantuvo su alto nivel de carga una vez que estuvo unido, desafiando la expectativa de que se neutralizaría.

Además, el estudio encontró que los fármacos son camaleones. Un fármaco puede flotar en el agua con una forma o carga específica, pero en el momento en que se bloquea en la proteína, a menudo cambia a una forma o estado de carga completamente diferente, y a veces más rara. De hecho, para algunos fármacos, la versión "minoritaria" que apenas existe en el agua se convierte en la versión "mayoritaria" una vez dentro de la proteína. Esto significa que si solo observas el fármaco en un tubo de ensayo, podrías estar diseñando para un fantasma que no existe en el objetivo real.

Curiosamente, mientras los fármacos estaban ocupados cambiando sus personalidades eléctricas, las proteínas fueron sorprendentemente tranquilas. Incluso cuando un fármaco altamente cargado se une a ellas, la carga eléctrica general de la proteína apenas se movió. Es como si la proteína fuera una multitud masiva y bien amortiguada; incluso si una persona (el fármaco) entra gritando y agitando una bandera, la multitud no entra en pánico ni cambia su estado de ánimo colectivo. La proteína gestiona estos cambios realizando pequeños ajustes en solo unos pocos puntos específicos cerca del sitio de unión, en lugar de cambiar toda su identidad eléctrica. Los investigadores también verificaron si la "pose" de la proteína (específicamente un bucle llamado DFG, que puede estar "dentro" o "fuera") dictaba estos cambios, pero descubrieron que los cambios de carga ocurrían independientemente de la pose de la proteína.

En resumen, este artículo sugiere que para predecir con precisión qué tan bien funcionará un fármaco, no podemos limitarnos a mirar el fármaco de forma aislada o asumir que permanece igual una vez que entra en el cuerpo. El fármaco y la proteína están en un baile dinámico donde se ajustan constantemente sus estados eléctricos para encajar el uno con el otro. Ignorar estos cambios es como intentar predecir los movimientos de un compañero de baile sin verlo moverse. Al tener en cuenta estos cambios fluidos de carga y forma, los científicos pueden construir mejores modelos para diseñar la próxima generación de medicinas que salvan vidas.

Resumen Técnico: Los inhibidores de quinasas pueden cambiar el estado de protonación o de tautómeros tras la unión

Planteamiento del problema
La predicción precisa de la afinidad de unión proteína-ligando ( $\Delta G$ ) es un pilar fundamental en el diseño de fármacos asistido por computadora (CADD). Sin embargo, los flujos de trabajo estándar suelen fallar al no considerar la distribución dinámica de los estados de protonación y tautomería tanto del ligando como de la proteína diana. Al pH fisiológico, residuos como la histidina, el aspartato, el glutamato, la arginina y la lisina existen en conjuntos de estados de ionización, y los ligandos poseen múltiples formas de protonación y tautomería accesibles. Los métodos actuales dependen frecuentemente de la asignación heurística de estados de protonación fijos o requieren un muestreo computacionalmente prohibitivo (por ejemplo, MD de pH constante) que es incompatible con los procesos de alto rendimiento. En consecuencia, la asignación errónea de estos estados puede conducir a errores significativos en las predicciones de la energía de unión, particularmente para los inhibidores de quinasas donde la dependencia del pH y los estados de ionización son críticos para la función y el diseño de fármacos.

Metodología
Los autores emplearon el método de Electroestática de Continuo de Múltiples Conformaciones (MCCE), un enfoque de Monte Carlo de Gran Canónica (GCMC) basado en la Electroestática de Continuo (CE), para investigar las distribuciones de protonación y tautomería en complejos quinasa-inhibidor.

Selección del sistema: El estudio analizó 37 estructuras cristalinas experimentales que comprenden 9 dominios de quinasas únicos (ABL1, ALK, CDK6, DDR1, EGFR, JAK3, MEK, MET, VEGFR) unidos a 18 inhibidores aprobados por la FDA.
Flujo de trabajo: Se realizaron tres cálculos independientes de MCCE para cada sistema: (1) la proteína apo (sin el ligando), (2) el complejo holo-proteína, y (3) el ligando aislado en solución.
Muestreo: El método muestrea rotámeros de las cadenas laterales, confórmeros del ligando y estados de protonación/tautomería. Las autoenergías y las interacciones de pares se calcularon utilizando la ecuación de Poisson-Boltzmann (resuelta mediante Delphi) con una constante dieléctrica de la proteína de 4 y un solvente de 80.
Análisis: Se computaron las probabilidades ponderadas por Boltzmann de los microestados para determinar la carga neta promedio y los estados dominantes de protonación/tautomería tanto de la proteína como del ligando en solución frente al estado unido. El estudio examinó específicamente los cambios de carga ( $\Delta q$ ) y las poblaciones tautoméricas tras la unión.

Contribuciones clave y resultados

Carga dinámica del ligando: Contrario a la suposición de que los ligandos pierden carga al entrar en el bolsillo hidrofóbico (debido a las penalizaciones por desolvatación), el estudio encontró que las cargas de los inhibidores individuales son altamente dinámicas. Mientras que algunos inhibidores se vuelven menos cargados al unirse, otros mantienen o incluso aumentan su carga neta. Por ejemplo, la Crizotinib pasa de un estado dominante de +2 en solución a un estado de +1 cuando se une a ALK, mientras que la Lapatinib pasa de un estado predominantemente +1 en solución a un estado neutro cuando se une a EGFR.
Estabilidad de la carga proteica: A pesar de la unión de ligandos cargados, la carga neta de los dominios de quinasa se mantuvo mayoritariamente estable. La carga neta de la proteína está efectivamente "amortiguada" por un gran conjunto de microestados con energías similares, lo que permite la redistribución de los estados de protonación sin cambios significativos en la carga neta total.
Cambios de protonación localizados: Los cambios de protonación significativos ( $\Delta q > 0.2$ ) se restringieron a una pequeña fracción de los residuos, predominantemente Lisina y Cisteína, a menudo ubicados cerca del sitio de unión pero ocasionalmente hasta 10 Å de distancia (por ejemplo, LYS1205 en ALK). La gran mayoría de los residuos no mostraron cambios en su estado de protonación entre las formas apo y holo.
Cambios tautoméricos: El estudio identificó casos en los que el tautómero de menor energía en solución no es la especie dominante en el estado unido. En algunos casos, las especies minoritarias en solución se convierten en el componente dominante en el complejo unido, con diferenciales de energía libre entre tautómeros que alcanzan los 3–6 kcal/mol.
Independencia conformacional: La magnitud de los cambios de carga tras la unión no mostró una fuerte dependencia de la conformación del bucle DFG (DFG-in vs. DFG-out). Ni la proteína ni el ligando exhibieron un patrón de cambio de carga predecible basado únicamente en el estado DFG.

Significado y afirmaciones
El artículo afirma que la predicción precisa de la energética de unión proteína-ligando requiere considerar el conjunto completo de estados de protonación y tautomería para ambos socios, en lugar de depender de las propiedades del estado de solución o de asignaciones fijas.

Necesidad del conjunto (Ensemble): Los hallazgos enfatizan que las "especies minoritarias" en solución pueden convertirse en dominantes en el estado unido, lo que requiere una consideración de todo el conjunto para evitar grandes errores en la predicación de la afinidad.
Asignación específica del objetivo: Los estados de protonación y tautomería deben tratarse como propiedades específicas del sitio de unión. No pueden extrapolarse del estado de solución del ligando, de la conformación global de la quinasa o de un único modelo fijo.
Eficiencia metodológica: El estudio demuestra que el MCCE proporciona una alternativa escalable y eficiente a la MD de pH constante para capturar estos efectos en sistemas farmacológicamente relevantes, permitiendo la identificación de diferenciales de energía libre significativos (3–6 kcal/mol) asociados con estados tautoméricos específicos.

Los autores concluyen que, si bien la carga neta de la proteína es robusta a la unión del ligando, el estado de protonación específico del ligando es altamente dependiente del contexto y debe determinarse dentro del complejo para asegurar un modelado físico fiable.

Más como este