Quantum chemical modelling of sequence encoded energetics and electrostatic landscapes governing early peptide association in α- synuclein and amyloid β

Sintesi Tecnica: Modellazione Chimico-Quantistica della Termodinamica Codificata dalla Sequenza e dei Paesaggi Elettrostatici nella Prima Associazione Peptidica

Problematica
Il ripiegamento errato e l'aggregazione delle proteine sono centrali nei disturbi neurodegenerativi come la malattia di Parkinson (PD) e la malattia di Alzheimer (AD), guidati dalla prima associazione molecolare di peptidi come l'$\alpha$-sinucleina e la beta-amiloide (A$\beta$). Mentre gli studi sperimentali hanno caratterizzato la cinetica e le strutture delle fibrille mature, i contributi energetici delle singole configurazioni molecolari durante le fasi iniziali, transitorie ed eterogenee di associazione rimangono scarsamente risolti. Esiste una lacuna critica nella comprensione di come l'ordine specifico dei residui (sequenza), distinto dalla composizione amminoacidica complessiva, modelli le proprietà molecolari che governano questi eventi di associazione iniziale. Questo studio mira a isolare il contributo dell'organizzazione della sequenza confrontando segmenti peptidici inclini all'aggregazione con controlli rimescolati (scrambled) con composizione corrispondente, determinando così come l'ordine dei residui influenzi l'energia, l'organizzazione strutturale e i paesaggi elettrostatici.

Metodologia
Lo studio impiega un flusso di lavoro computazionale coerente che integra la chimica quantistica semi-empirica e la teoria del funzionale della densità (DFT) per modellare gli stati monomerici e omodimerici.

• Sistemi: Segmenti inclini all'aggregazione ($\alpha$-sinucleina: VTGVTAVAQK; A$\beta$: LVFFAEDVGS) sono stati confrontati con controlli rimescolati a composizione corrispondente (AGATQTVVVK e EFLVGSDVAF, rispettivamente).
• Ottimizzazione della Geometria: Le strutture dei monomeri e degli omodimeri sono state ottimizzate utilizzando il metodo GFN2-xTB (Geometry, Frequency, Noncovalent extended Tight Binding, seconda generazione) implementato in xTB. Le strutture dei dimeri sono state inizialmente rilassate utilizzando il campo di forza GFN prima dell'ottimizzazione finale.
• Energia: Le energie di associazione ($\Delta E_{assoc}$) sono state calcolate come la differenza tra l'energia elettronica del dimero ottimizzato convergente tramite GFN2-xTB e il doppio dell'energia del monomero ottimizzato indipendentemente.
• Analisi della Struttura Elettronica: I calcoli DFT sono stati eseguiti su monomeri ottimizzati utilizzando GPAW con il funzionale di scambio-correlazione PBE e il formalismo Projector Augmented Wave (PAW) su una griglia in spazio reale (spaziatura di 0,20 Å).
• Descrittori:
  • Elettrostatica: Il potenziale elettrostatico molecolare (ESP) è stato mappato sulle superfici di densità elettronica (0,001 e Bohr$^{-3}$) ed analizzato tramite centratura della media per visualizzare i paesaggi relativi.
  • Distribuzione della Carica: L'analisi della popolazione Hirshfeld è stata utilizzata per quantificare le distribuzioni di carica atomica, inclusi i parametri statistici (media, mediana, deviazione standard) e le somme risolte per residuo.
  • Struttura/Sequenza: I descrittori includevano la decorazione della carica della sequenza (SCD), l'idropatia locale, la lunghezza della corsa idrofobica, il raggio di girazione ($R_g$), la distanza capo-coda e l'asfericità.
  • Interfaccia: I contatti intermolecolari sono stati quantificati tramite distanze tra atomi pesanti ($\le$ 4,5 Å e 5,0 Å) e contatti specifici N/O/S ($\le$ 3,5 Å).

Risultati Chiave
Lo studio rivela che il rimescolamento dei residui produce risposte distinte e specifiche per ogni peptide sia nei domini energetici che in quelli strutturali, nonostante la composizione amminoacidica conservata.

1. Energia del Monomero: L'ordine dei residui ha alterato significativamente l'energia elettronica dei monomeri isolati. Per l'A$\beta$, il controllo rimescolato era inferiore in energia di circa 45,33 kJ/mol rispetto alla sequenza incline all'aggregazione. Al contrario, per l'$\alpha$-sinucleina, la sequenza incline all'aggregazione era inferiore in energia di circa 151,12 kJ/mol rispetto al suo controllo rimescolato.
2. Energia di Associazione vs. Impacchettamento dell'Interfaccia:
   • A$\beta$: Sia la sequenza incline all'aggregazione (115,617 kJ/mol) che il controllo rimescolato (232,552 kJ/mol) hanno mostrato energie di associazione positive, indicando che i dimeri ottimizzati avevano un'energia elettronica superiore rispetto ai monomeri ottimizzati indipendentemente. Tuttavia, la sequenza incline all'aggregazione mostrava un'energia di associazione positiva inferiore (di 116,935 kJ/mol) rispetto al controllo rimescolato. Paradossalmente, questa sequenza formava meno contatti tra atomi pesanti (40 contro 81) e aveva una densità di contatto inferiore rispetto al controllo.
   • $\alpha$-Sinucleina: Il controllo rimescolato ha presentato un'energia di associazione più favorevole (più negativa) (-494,371 kJ/mol) rispetto alla sequenza incline all'aggregazione (-172,819 kJ/mol), con una differenza di 321,552 kJ/mol. Tuttavia, la sequenza incline all'aggregazione ha formato un'interfaccia più estesa, con 166 contatti tra atomi pesanti (rispetto a 105) e una maggiore densità di contatti stretti N/O/S (15 contro 5).
   • Conclusione: L'energia di associazione e l'impacchettamento dell'interfaccia (numero di contatti/densità) non variano linearmente; rappresentano misure complementari, non equivalenti, dell'associazione.
3. Paesaggi Elettrostatici (ESP): Il rimescolamento dei residui ha riorganizzato la distribuzione spaziale dell'ESP. Sebbene le coppie di A$\beta$ e $\alpha$-sinucleina condividessero stati di carica globali simili, l'organizzazione spaziale specifica delle regioni di potenziale positivo e negativo è cambiata significativamente tra le sequenze inclini all'aggregazione e quelle rimescolate. Ciò indica che l'ordine della sequenza detta la presentazione 3D del carattere elettrostatico.
4. Distribuzioni di Carica Hirshfeld: L'analisi della carica atomica ha mostrato che il rimescolamento dei residui ha ridistribuito principalmente le popolazioni di carica atomica piuttosto che spostare uniformemente le magnitudo di carica globale. Nell'$\alpha$-sinucleina, la proporzione di atomi positivi/negativi è rimasta quasi costante, ma le distribuzioni di probabilità dei valori di carica sono cambiate. Nell'A$\beta$, è stato osservato uno spostamento modesto nell'equilibrio tra atomi positivi e negativi insieme alla ridistribuzione.
5. Descrittori Sequenziali e Strutturali:
   • A$\beta$: Il rimescolamento ha alterato la decorazione della carica della sequenza (SCD), la lunghezza della corsa idrofobica (5 residui contro 3) e la compattezza molecolare, con il peptide incline all'aggregazione che presentava un $R_g$ di circa 10,1 Å rispetto ai circa 11,0 Å del controllo rimescolato.
   • $\alpha$-Sinucleina: Il rimescolamento ha alterato principalmente l'idropatia locale (la media massima della idropatia su 3 residui è aumentata da ~2,6 a ~4,2 nel controllo), mentre SCD, $R_g$ e asfericità sono rimasti comparabili.

Significato e Rivendicazioni
L'articolo stabilisce che l'ordine dei residui è un determinante fondamentale della prima associazione peptidica, indipendentemente dalla composizione amminoacidica. Lo studio dimostra che composizioni identiche possono generare caratteristiche energetiche, strutturali ed elettrostatiche distinte solo attraverso l'organizzazione molecolare dipendente dalla sequenza.

Le rivendicazioni chiave includono:

• La Sequenza come Determinante: La disposizione dei residui, non solo la composizione amminoacidica, governa l'equilibrio delle interazioni intermolecolari, influenzando sia la termodinamica dell'associazione che la geometria dell'interfaccia risultante.
• Risposte Specifiche per Peptide: Le conseguenze molecolari del rimescolamento dei residui non sono universali; l'$\alpha$-sinucleina e l'A$\beta$ mostrano risposte contrastanti in cui i cambiamenti nell'ordine della sequenza portano a diverse combinazioni di favorevolezza energetica e impacchettamento dell'interfaccia.
• Descrittori Complementari: Lo studio evidenzia che l'energia di associazione e le metriche di contatto dell'interfaccia catturano aspetti diversi dell'associazione molecolare. Un'interfaccia più popolata non corrisponde necessariamente a un'energia di associazione più favorevole, come osservato nelle tendenze divergenti tra i due sistemi peptidici.
• Riorganizzazione Elettrostatica: Identiche composizioni amminoacidiche possono produrre ambienti elettrostatici distinti a causa della riorganizzazione della schematizzazione della carica e della struttura 3D, fornendo una base molecolare per i comportamenti di associazione specifici della sequenza.

Il lavoro fornisce un quadro quantitativo per identificare i determinanti molecolari codificati dalla sequenza, suggerendo che il "codice" per la prima aggregazione risiede nella specifica disposizione lineare dei residui, che detta i risultanti paesaggi energetici ed elettrostatici.