Integrative In Silico Analysis of Synonymous Variants in ABCG2 and PPARA: Implications for Valproate Resistance in Epilepsy

Sintesi Tecnica: Analisi Integrativa In Silico di Varianti Sinonime in ABCG2 e PPARA

Problematica
Circa un terzo dei pazienti affetti da epilessia sviluppa resistenza ai farmaci antiepilettici, in particolare all'acido valproico (VPA). Il trasportatore ABCG2 e il suo regolatore trascrizionale PPARA sono implicati nella resistenza al VPA — specificamente attraverso una via in cui il VPA aumenta l'espressione di PPARA, il quale a sua volta incrementa l'espressione di ABCG2 per limitare la biodisponibilità del farmaco nel cervello — ma il contributo dei polimorfismi a singolo nucleotide sinonimi (sSNP) a questa resistenza rimane scarsamente compreso. Storicamente considerati "silenti" e funzionalmente neutri, gli sSNP sono sempre più riconosciuti come capaci di alterare l'espressione genica attraverso meccanismi quali l'alterazione dello splicing del pre-mRNA, i cambiamenti della struttura secondaria dell'mRNA e la modulazione della cinetica di traduzione. Questo studio affronta la lacuna nella comprensione di come gli sSNP in ABCG2 e PPARA possano contribuire meccanicisticamente alla variabilità interindividuale nella risposta al VPA.

Metodologia
Gli autori hanno condotto un'analisi integrativa in silico di 262 sSNP sinonimi in ABCG2 e 109 sSNP in PPARA, recuperati dal database dbSNP (riferimento GRCh38). Il flusso di lavoro ha valutato l'impatto funzionale attraverso quattro livelli regolatori:

1. Impatto sullo Splicing: Le varianti sono state analizzate utilizzando SpliceAI (una rete neurale profonda) per predire le alterazioni dei siti di splicing (guadagno/perdita di accettore/donatore) con una soglia di alta confidenza ($\Delta \ge 0,5$) e una soglia di sensibilità ($\Delta \ge 0,2$). I candidati promettenti sono stati ulteriormente valutati utilizzando Human Splicing Finder (HSF) per valutare gli effetti sugli enhancer dello splicing esonico (ESE), i silencer dello splicing esonico (ESS) e i siti di splicing criptici.
2. Struttura Secondaria dell'RNA: Lo strumento RNAsnp (Modalità 2 per trascritti lunghi, Modalità 1 per finestre localizzate) è stato utilizzato per predire i cambiamenti strutturali locali. La significatività è stata determinata tramite p-value ($\le 0,1$) e variazioni dell'energia libera minima ($\Delta\Delta G$).
3. Cinetica di Traduzione:
   • Bias dell'Uso dei Codoni: Sono stati calcolati i rapporti di Uso Relativo dei Codoni Sinonimi (RSCU) per identificare gli spostamenti nella preferenza dei codoni.
   • Efficienza di Traduzione: È stato applicato il framework dell'indice di adattamento del tRNA (tAI) utilizzando i numeri di copia dei geni tRNA (tGCN) da GtRNAdb e i parametri di vincolo derivati dal lievito ($s_{ij}$) per stimare le variazioni nell'efficienza di decodifica ($\Delta W_i$).
   • Classificazione della Decodifica: Le varianti sono state categorizzate in base alle transizioni tra i meccanismi di decodifica Watson-Crick, wobble e mista.
4. Ripiegamento Co-traduzionale: Le varianti sono state mappate su regioni sensibili note per influenzare il ripiegamento proteico, specificamente la regione N-terminale precoce (codoni 1–20) e le regioni 20–50 codoni a valle dei confini di dominio.
5. Analisi Statistica: È stata eseguita l'analisi di correlazione di Pearson per valutare la relazione tra i cambiamenti nel bias dell'uso dei codoni e l'efficienza di traduzione.

Risultati Chiave
L'analisi ha identificato un sottoinsieme di sSNP in entrambi i geni con effetti funzionali predetti:

• Disruzione dello Splicing: SpliceAI ha identificato 16 sSNP significativi in ABCG2 e 3 in PPARA. Le varianti degne di nota includono ABCG2 c.1134A>G (guadagno/perdita duale di donatore) e c.786A>T (creazione di un sito donatore criptico e grave interruzione di ESE/ESS). L'analisi HSF ha confermato che queste varianti spesso interrompono i motivi regolatori, portando potenzialmente all'exon skipping o a frameshift. In PPARA, c.333C>T e c.126C>T sono stati predetti come creatori di siti donatori criptici. L'analisi HSF ha confermato che queste varianti spesso interrompono i motivi regolatori, portando potenzialmente all'exon skipping o a frameshift.
• Struttura Secondaria dell'mRNA: In ABCG2, 20 varianti (7,6%) hanno mostrato una significativa disruzione strutturale ($p \le 0,1$). La più destabilizzante è stata c.480G>A ($\Delta\Delta G = +6,6$ kcal/mol), mentre c.1044A>G ha mostrato una significativa stabilizzazione ($\Delta\Delta G = -4,0$ kcal/mol). In PPARA, 4 varianti hanno soddisfatto la soglia di significatività, con c.1323G>A che causa un riarrangiamento globale della struttura dell'RNA, interrompendo un'interazione elicoidale a lungo raggio che si estende per 240 nucleotidi.
• Efficienza di Traduzione e Bias dei Codoni:
  • RSCU: In ABCG2, 29 varianti hanno mostrato una preferenza aumentata (rapporto >2) e 56 una preferenza diminuita (rapporto <0,5). Tendenze simili sono state osservate in PPARA.
  • Efficienza ($\Delta W_i$): In ABCG2, il 45,8% delle varianti ha ridotto l'efficienza di traduzione, mentre il 40,8% l'ha aumentata. In PPRA, il 51,3% ha ridotto l'efficienza.
  • Correlazione: È stata trovata una moderata correlazione positiva tra i cambiamenti nel bias dell'uso dei codoni ($\Delta$RSCU) e l'efficienza di traduzione ($\Delta W_i$) per entrambi, ABCG2 ($r=0,478$) e PPARA ($r=0,567$).
• Ripiegamento Co-traduzionale: Sono state identificate varianti in regioni di ripiegamento sensibili: 11 nella regione N-terminale precoce e 25 a valle dei confini di dominio per ABCG2; 11 e 8 rispettivamente per PPARA.

Significatività e Rivendicazioni
L'articolo sostiene che le varianti sinonime in ABCG2 e PPARA non sono funzionalmente neutre, ma possono influenzare la resistenza al valproato attraverso molteplici meccanismi molecolari:

1. Fedeltà dello Splicing: La disruzione degli elementi regolatori dello splicing può portare a isoforme di trascritto aberranti, risultando potenzialmente in trasportatori o recettori non funzionali.
2. Stabilità dell'mRNA: La destabilizzazione strutturale può accelerare la degradazione del trascritto, riducendo l'abbondanza proteica, mentre la stabilizzazione potrebbe alterare i livelli di espressione.
3. Dinamiche di Traduzione: I cambiamenti nell'uso dei codoni e nell'efficienza di decodifica possono alterare i tassi di allungamento ribosomiale, causando potenzialmente il pausing del ribosoma che interrompe il ripiegamento proteico co-traduzionale.

Gli autori propongono un quadro meccanicistico in cui gli sSNP in PPARA potrebbero compromettere la sua capacità di regolare positivamente ABCG2 in risposta al VPA, o gli sSNP in ABCG2 potrebbero alterare direttamente la funzione o la stabilità del trasportatore. Queste perturbazioni convergono nell'alterare la farmacocinetica del valproato alla barriera emato-encefalica.

Conclusione
Lo studio conclude che un sottoinsieme biologicamente significativo di varianti sinonime in ABCG2 e PPARA presenta perturbazioni predette sufficienti a contribuire alla resistenza al valproato. Gli autori raccomandano l'integrazione sistematica del profilo degli sSNP nello screening farmacogenomico per migliorare la medicina di precisione nell'epilessia. Essi sottolineano che questi risultati supportano un quadro di valutazione variante per variante che considera le dimensioni strutturali, di splicing e di traduzione, piuttosto che assumere la neutralità basata esclusivamente sulla conservazione amminoacidica. L'articolo nota che questi risultati richiedono una futura validazione sperimentale e studi di associazione clinica per confermare le relazioni genotipo-fenotipo.