Mutation-induced conformational remodeling explains gain- and loss-of-function phenotypes in Euglena gracilis phytoene synthase

Resumen Técnico: Remodelación Conformacional Inducida por Mutaciones en la Fitoeno Sintasa de Euglena gracilis

Planteamiento del Problema
La fitoeno sintasa (CrtB) cataliza el primer paso limitante de la velocidad en la biosíntesis de carotenoides, lo que la convierte en un objetivo crítico para la ingeniería metabólica. Si bien el cribado experimental en Euglena gracilis ha identificado mutaciones específicas (Y197A, E121H y L349F) que dan lugar a fenotipos contrastantes de ganancia de función (producción mejorada de carotenoides) y pérdida de función (reducción de la producción), los mecanismos estructurales que impulsan estas diferencias siguen sin estar claros. Los modelos estructurales estáticos y las observaciones experimentales por sí solos no logran explicar cómo las sustituciones de un solo residuo alteran la actividad enzimática. El estudio aborda la brecha en la comprensión de la base estructural dinámica que vincula mutaciones específicas con la alteración de la productividad de carotenoides.

Metodología
Los autores emplearon simulaciones de Dinámica Molecular (MD) de la forma apo (libre de sustrato) de la CrtB de E. gracilis para investigar la dinámica conformacional de la enzima silvestre y tres variantes: Y197A, E121H y L349F. El enfoque computacional incluyó:

• Análisis Energético: Cálculos de MM-GBSA y MM-PBSA para evaluar las perturbaciones energéticas dependientes de la mutación y la estabilidad ($\Delta\Delta G$).
• Caracterización Electrostática: Análisis de PROPKA para determinar los cambios en los puntos isoeléctricos (pI) y la carga neta a pH 7.0.
• Métricas Estructurales Globales: Análisis del Área de Superficie Accesible al Solvente (SASA) y el Radio de Giro ($R_g$) para evaluar la compacidad global y el muestreo conformacional.
• Dinámica Local: Análisis de la Fluctuación de la Raíz Cuadrada Media (RMSF) y $\Delta$RMSF para mapear los cambios en la movilidad específica de los residuos.
• Arquitectura del Sitio Activo: Examen detallado de la accesibilidad al solvente, volúmenes de los bolsillos de unión y redes de enlaces de hidrógeno dentro de los motivos ricos en Asp (esenciales para la coordinación de Mg²⁺) y la solapa hidrofóbica.

Resultados Clave
Las simulaciones revelaron que las mutaciones inducen una remodelación distinta del paisaje de energía conformacional, afectando las redes de interacción local y la arquitectura del sitio activo sin causar cambios estructurales globales sustanciales:

1. Desplazamientos Energéticos y Electrostáticos:

   • E121H exhibió la mayor redistribución de la dinámica conformacional, con valores positivos de $\Delta\Delta G$ que indican una desestabilización energética en relación con la proteína nativa. También mostró el pI predicho más alto y una ligera carga positiva a pH neutro.
   • Y197A mostró cambios energéticos casi neutros o favorables.
   • L349F indujo perturbaciones energéticas modestas, pero mostró una desestabilización pronunciada en las predicciones de la energía libre de plegamiento de PROPKA.

2. Muestreo Conformacional y Flexibilidad:

   • Variantes de ganancia de función (Y197A, E121H): Ambas retuvieron solapas hidrofóbicas relativamente accesibles. E121H mostró la distribución de $R_g$ más amplia, mientras que Y197A y E121H exhibieron distribuciones de SASA más amplias que la proteína nativa, lo que sugiere un muestreo conformacional mejorado manteniendo la compacidad general.
   • Variante de pérdida de función (L349F): Se asemejó estrechamente a la proteína nativa en las distribuciones de SASA y $R_g$ globales, pero mostró una constricción localizada y dramática de la solapa del bolsillo hidrofóbico (reducción de volumen de aproximadamente el 35%).

3. Redes de Interacción Local y Remodelación del Sitio Activo:

   • Solapa Hidrofóbica: La mutación L349F causó una reducción significativa en el volumen del bolsillo de la solapa y en la accesibilidad al solvente, consistente con una conformación más cerrada. En contraste, las variantes de ganancia de función mantuvieron una mayor exposición de la solapa.
   • Motivos Ricos en Asp: Todos los mutantes mostraron un aumento en la ocupación de enlaces de hidrógeno en la región rica en Asp 1. Notablemente, L349F mostró la ocupación más alta a pesar de su rendimiento catalítico reducido, lo que sugiere que la eficiencia catalítica depende de un equilibrio entre estabilidad y flexibilidad, más que de una estabilización local máxima.
   • Rigidez en L349F: El análisis estructural del mutante L349F reveló que el residuo Tyr127 adoptó una conformación más rígida, formando una red de interacción persistente que involucra interacciones $\pi$-catión y $\pi$-anión con Arg154 y Glu158. Esto contrasta con la proteína nativa y las variantes de ganancia de función, donde Tyr127 permaneció conformacionalmente flexible.

Significancia y Afirmaciones
El artículo afirma proporcionar un marco estructural que vincula la remodelación conformacional inducida por mutaciones alrededor del centro catalítico con las diferencias observadas experimentalmente en la producción de carotenoides. Los autores sostienen que los fenotipos contrastantes surgen no meramente de las perturbaciones estructurales locales en los sitios de mutación, sino de la redistribución del muestreo conformacional alrededor de motivos funcionalmente importantes, específicamente la solapa hidrofóbica y el motivo adyacente rico en Asp 1.

El estudio destaca que el fenotipo de pérdida de función en L349F se caracteriza por una "constricción pronunciada" de la solapa hidrofóbica y una "red de interacción rígida", mientras que las variantes de ganancia de función retienen una "solapa hidrofóbica relativamente accesible" y una mayor flexibilidad. Estos hallazgos sugieren que la organización estructural dinámica es un determinante crítico de la función de CrtB. Los autores concluyen que estos conocimientos ofrecen una base para la ingeniería racional de las fitoeno sintasas, aunque señalan explícitamente que se requieren simulaciones futuras que incorporen sustratos (GGPP) e iones Mg²⁺ para elucidar completamente cómo estos cambios conformacionales intrínsecos influyen en el reconocimiento del sustrato y la catálisis.