Computationally guided modifications of CviUPO to improve catalytic activity

Sintesi Tecnica: Modifiche della CviUPO guidate dal calcolo computazionale per migliorare l'attività catalitica

Definizione del Problema
Le perossigenasi aspecifiche (UPO) sono biocatalizzatori promettenti capaci di ossifunzionalizzare selettivamente gli idrocarburi saturi utilizzando solo il perossido di idrogeno ($H_2O_2$) come co-substrato. Tuttavia, la loro applicazione industriale è ostacolata dall'inattivazione enzimatica indotta dal perossido, un processo probabilmente innescato dai radicali ossidrilici che ossidano gli amminoacidi nella tasca attiva. Sebbene l'ingegneria mirata degli enzimi sia essenziale per mitigare questo fenomeno e migliorarne le prestazioni, l'influenza specifica degli amminoacidi del sito attivo sul ciclo catalitico rimane scarsamente compresa. Questo studio affronta la necessità di elucidare come le modifiche ai residui chiave del centro attivo influenzino le barriere energetiche e l'efficienza catalitica di Collariella virescens UPO (CviUPO).

Metodologia
Gli autori hanno impiegato un workflow computazionale sinergico combinando simulazioni di Dinamica Molecolare (MD) e calcoli di Meccanica Quantistica/Meccanica Molecolare (QM/MM) per investigare due specifiche mutazioni in CviUPO:

1. Glu162 $\to$ Asp162: Sostituendo l'acido glutammico (coinvolto nell'astrazione del protone) con l'acido aspartico, più corto, per testare se un aumento del volume della tasca migliori l'attività.
2. Cys19 $\to$ His19: Sostituendo la cisteina che ancora l'eme con l'istidina per investigare una potenziale resistenza al bleaching dell'eme o uno spostamento verso un comportamento simile a quello della perossidasi.

Lo studio ha utilizzato le seguenti tecniche computazionali:

• Simulazioni MD Classiche: Sono state eseguite simulazioni di 120 ns su CviUPO nativa, modificata con Asp e modificata con His per valutare la stabilità strutturale, il volume della tasca attiva e la dinamica di idratazione. La stabilità del substrato (etilbenzene/ETBE) è stata testata tramite 100 ns di MD dopo il docking.
• Simulazioni QM/MM: Utilizzando il software ORCA con il funzionale B3LYP (base set def2-TZVP) e il campo di forza CHARMM36, lo studio ha modellato il ciclo catalitico fino alla formazione del Composto I (il complesso ferro-oxo reattivo). Ciò ha incluso ottimizzazioni della geometria e analisi delle frequenze per gli stati Low Spin (LS), Intermediate Spin (IS) e High Spin (HS).
• Simulazioni Nudged Elastic Band (NEB): Sono state condotte calcoli CI-NEB per mappare i percorsi di minima energia e identificare le barriere energetiche tra i key intermedi (Stato Iniziale $\to$ Composto 0 $\to$ Composto I).
• Scansioni di Superficie Energetica Rilassata: Eseguite per l'enzima nativo e l'enzima modificato con His insieme al substrato nella tasca attiva per investigare la seconda metà del ciclo catalitico (Composto I $\to$ Composto II $\to$ Prodotto).

Risultati Chiave

• Comportamento Dinamico e Idratazione: Tutte e tre le varianti enzimatiche sono rimaste stabili per 120 ns. Sebbene le modifiche con Asp e His abbiano comportato volumi della tasca attiva leggermente inferiori rispetto all'enzima nativo, la densità media delle molecole d'acqua all'interno delle tasche è rimasta comparabile in tutte le varianti. Notevolmente, l'enzima modificato con His ha mostrato una maggiore variazione temporale nel volume della tasca, suggerendo una potenziale instabilità strutturale.
• Formazione del Composto I (Nativo vs. Modificato):
  • Modifica Asp: Le simulazioni NEB hanno rivelato una barriera energetica significativamente aumentata (circa 3 eV) per la transizione dallo stato iniziale al Composto 0. Ciò è attribuito al residuo di acido aspartico più corto, che aumenta la distanza che il protone deve percorrere da $H_2O_2$, ostacolando la deprotonazione. Di conseguenza, la modifica Asp è prevista come deleteria per l'attività catalitica.
  • Modifica His: I risultati NEB hanno indicato l'assenza di una barriera energetica tra lo stato iniziale e il Composto 0, suggerendo una facile formazione del complesso ferro-idroperossido. Tuttavia, la barriera tra il Composto 0 e il Composto I era più alta (1,25–1,63 eV) rispetto all'enzima nativo. Inoltre, lo stato finale (Composto I) non era energeticamente inferiore al Composto 0, sollevando la possibilità che il ciclo catalitico possa rimanere intrappolato in un minimo locale.
• Stati di Spin: Per gli enzimi modificati, le energie calcolate per gli stati LS, IS e HS erano quasi identiche, indicando che piccoli cambiamenti nell'ambiente elettronico locale influenzano significativamente la configurazione elettronica preferita. Al contrario, l'enzima nativo mostrava una separazione più distinta tra gli stati di spin, con lo stato IS che generalmente presentava i valori di energia più bassi.
• Interazione con il Substrato: Le simulazioni MD hanno confermato che l'ETBE è stabile nelle tasche attive degli enzimi nativo e modificato con His, ma instabile nella variante modificata con Asp.
• Percorso di Reazione con il Substrato: Le scansioni di superficie rilassata per l'enzima modificato con His hanno mostato una barriera inferiore per la formazione del Composto II (0,67 eV) rispetto all'enzima nativo (1,37 eV). Tuttavia, la scansione ha suggerito una deviazione meccanicistica in cui una molecola d'acqua, piuttosto che un diretto rebound del gruppo ossidrile, potrebbe facilitare la formazione del prodotto. Ciò solleva la possibilità che la modifica con His alteri il meccanismo dell'enzima, potenzialmente convertendolo da perossigenasi a perossidasi, sebbene ciò non sia stato confermato definitivamente.

Significato e Rivendicazioni
L'articolo conclude che, sebbene le simulazioni computazionali abbiano avuto successo nel chiarire i processi fondamentali del ciclo catalitico di CviUPO e abbiano evidenziato i ruoli critici dell'idratazione della tasca attiva e degli stati di spin, esse sono insufficienti per determinare definitivamente se le modifiche testate aumentino l'attività catalitica complessiva.

• La modifica Asp è computazionalmente prevista come sfavorevole a causa delle elevate barriere energetiche associate al trasferimento protonico.
• La modifica His presenta un compromesso complesso: facilita la formazione iniziale del complesso ferro-idroperossido, ma introduce barriere più elevate nelle fasi successive del ciclo e può alterare interamente il meccanismo di reazione.
  L'autore sottolinea che le simulazioni da sole non possono determinare in modo conclusivo la specificità del substrato o confermare se la modifica con His converta l'enzima in una perossidasi. Essi affermano che un approccio sinergico, che connetta le previsioni teoriche con la verifica sperimentale, è necessario per un'efficace ingegneria enzimatica. I tentativi sperimentali preliminari di esprimere l'enzima modificato con His sono stati difficili, producendo scarse quantità di enzima e caricamento dell'eme, il che sottolinea ulteriormente la necessità di una validazione sperimentale.

Lo studio funge da dimostrazione di come i calcoli chimici quantistici di alto livello, se integrati con la dinamica molecolare, possano fornire approfondimenti profondi sulle interazioni elettroniche e sui paesaggi energetici, guidando futuri sforzi di progettazione razionale pur riconoscendo i limiti delle sole previsioni teoriche.