Phenylalanine modification in plasma-driven biocatalysis revealed by solvent accessibility and reactive dynamics in combination with protein mass spectrometry

Resumen Técnico: Modificación de la Fenilalanina en la Biocatálisis Impulsada por Plasma

Planteamiento del Problema
La biocatálisis ofrece una alternativa respetuosa con el medio ambiente a la catálisis química convencional, pero enfrenta desafíos significativos relacionados con el suministro de co-sustratos reactivos, particularmente el peróxido de hidrógeno ($H_2O_2$). Aunque es esencial para la rotación catalítica, las altas concentraciones de $H_2O_2$ y otras especies generadas por plasma (PGS) pueden inactivar las enzimas. Las fuentes de plasma no térmico (por ejemplo, descarga de barrera dieléctrica) pueden suministrar estas especies en concentraciones controladas, pero las interacciones específicas entre las PGS (tales como $H_2O_2$, $OH$, $O$, $O_2$, $NO$y$H$) y las estructuras enzimáticas aún se comprenden poco a nivel atomístico. Comprender estas interacciones es crítico para prevenir la desactivación enzimática manteniendo la eficiencia catalítica.

Metodología
El estudio emplea un enfoque computacional y experimental combinado para investigar las interacciones entre las PGS y tres enzimas: AaeUPO (una peroxigenasa inespecífica modelo), CviUPO (la principal enzima modelo para este estudio) y GapA (una enzima glucolítica inactivada por $H_2O_2$).

1. Tamizaje Computacional (SASA y MD):

   • Análisis del Área de Superficie Accesible al Solvente (SASA): Se utilizó un paquete personalizado construido sobre VMD y LAMMPS para mapear las energías de interacción de las PGS con la superficie accesible al solvente de las enzimas. Esto implicó colocar moléculas de prueba en los puntos de SASA (a 1.4 Å de la superficie) y realizar minimizaciones de energía restringidas para identificar regiones de fuerte interacción.
   • Dinámica Molecular Reactiva (rMD): Utilizando el campo de fuerza ReaxFF (parametrización de biomoléculas por Monti et al.) dentro de LAMMPS, el estudio simuló la ruptura y formación dinámica de enlaces. Las simulaciones se realizaron a 300 K, 400 K y 500 K, tanto en vacío como en solvente explícito, para evaluar la estabilidad de los sitios de interacción predichos y el efecto del desplegamiento proteico.
   • Simulaciones de Alta Concentración: Para imitar condiciones de difusión realistas, se realizaron simulaciones con concentraciones elevadas de PGS (por ejemplo, 100 radicales $OH$) distribuidos aleatoriamente alrededor de la enzima.

2. Validación Experimental:

   • Tratamiento con Plasma: La proteína recombinante CviUPO (rCviUPO) fue tratada con una fuente de plasma de descarga de barrera dieléctrica (DBD) durante 5 minutos.
   • Espectrometría de Masas (MS): Las proteínas tratadas fueron reducidas, alquiladas y digeridas con tripsina. Se realizó un análisis LC-MS/MS para detectar modificaciones específicas de aminoácidos, centrándose en residuos con energías de interacción predichas altas (por ejemplo, Fenilalanina, Metionina, Cisteína, Arginina, Lisina).

Resultados Clave

• Perfiles de Interacción: El análisis de SASA reveló que las PGS interactúan más fuertemente con residuos cargados y polares. Específicamente, los radicales $OH$ mostraron fuertes interacciones con la Lisina (Lys), Fenilalanina (Phe) y Arginina (Arg). El $H_2O_2$ interactuó fuertemente con la Serina (Ser) y el Glutamato (Glu).
• Efectos de la Temperatura: A medida que la temperatura aumentó (simulando el desplegamiento), los perfiles de interacción cambiaron. A 400 K y 500 K, los residuos previamente ocultos se volvieron accesibles, aumentando la probabilidad de modificación en el sitio activo (por ejemplo, cerca del centro de hemo en CviUPO).
• Rol del Solvente: En las simulaciones con solvente, la mayoría de las especies reactivas reaccionaron con las moléculas de agua en lugar de con la superficie de la proteína, lo que sugiere que el solvente actúa como una barrera protectora. Sin embargo, las simulaciones de alta concentración confirmaron que, con el tiempo, las especies aún se unen a los residuos predichos por el tamizaje SASA.
• Confirmación Experimental: La espectrometría de masas de CviUPO tratada con plasma confirmó la oxidación significativa de los residuos de Fenilalanina (Phe51, Phe67, Phe208) y de Metionina (Met42, Met220).
  • La modificación de la Phe concuerda con la predicción de SASA de una fuerte interacción de $OH$ con los anillos aromáticos expuestos.
  • La modificación de la Met ocurrió a pesar de las bajas puntuaciones de interacción de SASA, lo que los autores atribuyen a que los residuos de Met están ocultos pero siguen siendo accesibles mediante difusión o entornos locales específicos no capturados totalmente por el modelo estático de SASA.
  • Las modificaciones a Arg y Lys se detectaron solo a intensidades bajas, mientras que las modificaciones de Cys no fueron observadas, probablemente debido al contexto estructural específico (por ejemplo, Cys19 está covalentemente unida al hemo).

Significancia y Reivindicaciones
El artículo sostiene que la combinación del tamizaje de interacción basado en SASA y la dinámica molecular reactiva es una herramienta viable para predecir posibles sitios de modificación en enzimas expuestas al plasma. El estudio demuestra que:

1. Utilidad Predictiva: El enfoque computacional identificó con éxito la Fenilalanina como un objetivo primario para la oxidación inducida por plasma, un hallazgo que fue validado experimentalmente. Esto sugiere que el método puede dirigir los esfuerzos experimentales que consumen mucho tiempo hacia residuos específicos.
2. Información Mecanística: Los resultados resaltan que la accesibilidad al solvente y la ruta de difusión de las especies reactivas son factores críticos; la difusión prolongada a través de la fase líquida aumenta la probabilidad de que las especies reaccionen con el solvente en lugar de con la enzima.
3. Integración Metodológica: El trabajo subraya el valor de integrar el tamizaje computacional con la espectrometría de masas. Si bien el método SASA proporcionó un "escaneo grueso" que priorizó con éxito la Phe para el análisis, los autores señalan que una alta energía de interacción no garantiza la unión, y la validación experimental sigue siendo necesaria.

Los autores concluyen que, si bien el enfoque SASA/MD es efectivo para el tamizaje inicial, el trabajo futuro debe centrarse en expandir el rango de modificaciones probadas y aplicar el modelo a otros sistemas para verificar su validez general. El estudio no pretende haber resuelto el problema de la inactivación enzimática, sino que proporciona un marco para comprender y predecir los objetivos moleculares específicos del daño inducido por plasma.