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

Sintesi Tecnica: Rimodellamento Conformazionale Indotto da Mutazione in la Fitoene Sintasi di Euglena gracilis

Problema
La fitoene sintasi (CrtB) catalizza il primo passaggio, limitante per velocità, nella biosintesi dei carotenoidi, rendendola un bersaglio crito per l'ingegneria metabolica. Sebbene lo screening sperimentale in Euglena gracilis abbia identificato mutazioni specifiche (Y197A, E121H e L349F) che danno origine a fenotipi contrastanti di guadagno di funzione (produzione di carotenoidi aumentata) e perdita di funzione (produzione ridotta), i meccanismi strutturali che guidano queste differenze rimangono poco chiari. I modelli strutturali statici e le osservazioni sperimentali da soli non riescono a spiegare come le sostituzioni di singoli residui alterino l'attività enzimatica. Lo studio affronta la lacuna nella comprensione della base strutturale dinamica che lega le specifiche mutazioni alla produttività dei carotenoidi.

Metodologia
Gli autori hanno impiegato simulazioni di Dinamica Molecolare (MD) della forma apo (priva di substrato) di CrtB di E. gracilis per investigare la dinamica conformazionale dell'enzima wild-type e di tre varianti: Y197A, E121H e L349F. L'approccio computazionale ha incluso:

• Analisi Energetica: Calcoli MM-GBSA e MM-PBSA per valutare le perturbazioni energetiche dipendenti dalla mutazione e la stabilità ($\Delta\Delta G$).
• Caratterizzazione Elettrostatica: Analisi PROPKA per determinare gli spostamenti nei punti isoelettrici (pI) e la carica netta a pH 7,0.
• Metriche Strutturali Globali: Analisi della Superficie Accessibile al Solvente (SASA) e del Raggio di Giro ($R_g$) per valutare la compattezza globale e il campionamento conformazionale.
• Dinamica Locale: Analisi della Deviazione Standard della Fluttuazione (RMSF) e $\Delta$RMSF per mappare i cambiamenti di mobilità specifici per residuo.
• Architettura del Sito Attivo: Esame dettagliato dell'accessibilità al solvente, dei volumi delle tasche di legame e delle reti di legami a idrogeno all'interno dei motivi ricchi di Asp (essenziali per la coordinazione di Mg²⁺) e del flap idrofobico.

Risultati Chiave
Le simulazioni hanno rivelato che le mutazioni inducono un distinto rimodellamento del panorama energetico conformazionale, influenzando le reti di interazione locale e l'architettura del sito attivo senza causare cambiamenti strutturali globali sostanziali:

1. Spostamenti Energetici ed Elettrostatici:

   • E121H ha esibito la più grande ridistribuzione della dinamica conformazionale, con valori di $\Delta\Delta G$ positivi che indicano una destabilizzazione energetica rispetto alla proteina nativa. Ha inoltre mostrato il pI predetto più alto e una leggera carica positiva a pH neutro.
   • Y197A ha mostrato cambiamenti energetici quasi neutri o favorevoli.
   • L349F ha indotto modeste perturbazioni energetiche ma ha mostrato una pronunciata destabilizzazione nelle previsioni dell'energia libera di ripiegamento PROPKA.

2. Campionamento Conformazionale e Flessibilità:

   • Varianti di guadagno di funzione (Y197A, E121H): Entrambe hanno mantenuto flap idrofobici relativamente accessibili. E121H ha mostrato la distribuzione di $R_g$ più ampia, mentre Y197A ed E121H hanno esibito distribuzioni di SASA più ampie rispetto alla proteina nativa, suggerendo un potenziamento del campionamento conformazionale pur mantenendo la compattezza complessiva.
   • Variante di perdita di funzione (L349F): Ha somigliato strettamente alla proteina nativa nelle distribuzioni globali di SASA e $R_g$, ma ha mostrato una drastica e localizzata costrizione del flap idrofobico (riduzione del volume di circa il 35%).

3. Reti di Interazione Locale e Rimodellamento del Sito Attivo:

   • Flap Idrofobico: La mutazione L349F ha causato una significativa riduzione del volume della tasca del flap e dell'accessibilità al solvente, coerentemente con una conformazione più chiusa. Al contrario, le varianti di guadagno di funzione hanno mantenuto una maggiore esposizione del flap.
   • Motivi Ricchi di Asp: Tutti i mutanti hanno mostrato un aumento dell'occupanza dei legami a idrogeno nel motivo Asp-ricco 1. Notevolmente, L349F ha esibito l'occupanza più elevata nonostante la ridotta performance catalitica, suggerendo che l'efficienza catalitica dipenda da un equilibrio tra stabilità e flessibilità piuttosto che dalla massimale stabilizzazione locale.
   • Rigidificazione in L349F: L'analisi strutturale del mutante L349F ha rivelato che il residuo Tyr127 ha adottato una conformazione più rigida, formando una rete di interazione persistente che coinvolge interazioni $\pi$-catione e $\pi$-anione con Arg154 e Glu158. Ciò contrasta con la proteina nativa e le varianti di guadagno di funzione, dove Tyr127 rimaneva conformazionalmente flessibile.

Significatività e Rivendicazioni
L'articolo sostiene di fornire un quadro strutturale che collega il rimodellamento conformazionale indotto dalla mutazione attorno al centro catalitico alle differenze osservate sperimentalmente nella produzione di carotenoidi. Gli autori affermano che i fenotipi contrastanti derivano non solo da perturbazioni strutturali locali nei siti di mutazione, ma dalla ridistribuzione del campionamento conformazionale attorno a motivi funzionalmente importanti, specificamente il flap idrofobico e il motivo adiacente Asp-ricco 1.

Lo studio evidenzia che il fenotipo di perdita di funzione in L349F è caratterizzato da una "pronunciata costrizione" del flap idrofobico e da una "rete di interazione rigida", mentre le varianti di guadagno di funzione mantengono un "flap idrofobico relativamente accessibile" e una maggiore flessibilità. Questi risultati suggeriscono che l'organizzazione strutturale dinamica sia un determinante critico della funzione di CrtB. Gli autori concludono che queste intuizioni offrono una base per l'ingegneria razionale delle fitoene sintasi, sebbene specifichino esplicitamente che sono necessarie simulazioni future che includano i substrati (GGPP) e gli ioni Mg²⁺ per elucidare pienamente come questi cambiamenti conformazionali intrinseci influenzino il riconoscimento del substrato e la catalisi.