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

Technische Samenvatting: Kwantumchemische Modellering van Sequentie-gecodeerde Energetica en Elektrostatische Landschappen in Vroege Peptide-associatie

Probleemstelling
Eiwitmisvalling en aggregatie staan centraal bij neurodegeneratieve aandoeningen zoals Parkinson (PD) en Alzheimer (AD), gedreven door de vroege moleculaire associatie van peptiden zoals $\alpha$-synucleïne en amyloïd $\beta$ (A$\beta$). Hoewel experimentele studies de kinetiek en structuren van volwassen fibrillen hebben gekarakteriseerd, blijven de energetische bijdragen van individuele moleculaire configuraties tijdens de transiënte, heterogene vroege stadia van associatie slecht opgelost. Er bestaat een kritische kloof in het begrip van hoe de specifieke lineaire volgorde van residuen (sequentie), onderscheiden van de algehele aminozuurcompositie, de moleculaire eigenschappen vormgeeft die deze initiële associatiegebeurtenissen beheersen. Deze studie beoogt de bijdrage van de sequentie-arrangement te isoleren door aggregatiegevoelige peptidefragmenten te vergelijken met compositie-gematchte gescrambelde controles, om zo te bepalen hoe de volgorde van residuen de energetica, structurele organisatie en elektrostatische landschappen beïnvloedt.

Methodologie
De studie maakt gebruik van een consistente computationele workflow die semi-empirische kwantumchemie en dichtheidsfunctionaaltheorie (DFT) integreert om monomeer- en homodimeerstaten te modelleren.

• Systemen: Aggregatiegevoelige segmenten ($\alpha$-synucleïne: VTGVTAVAQK; A$\beta$: LVFFAEDVGS) werden vergeleken met compositie-gematchte gescrambelde controles (AGATQTVVVK en EFLVGSDVAF, respectievelijk).
• Geometrie Optimalisatie: Monomeer- en homodimeerstructuren werden geoptimaliseerd met de GFN2-xTB methode (Geometry, Frequency, Noncovalent extended Tight Binding, tweede generatie) geïmplementeerd in xTB. Dimerstructuren werden aanvankelijk ontspannen met het GFN-krachtveld voor de uiteindelijke optimalisatie.
• Energetica: Associatie-energieën ($\Delta E_{assoc}$) werden berekend als het verschil tussen de geconvergeerde GFN2-xTB elektronische energie van het geoptimaliseerde dimeer en twee keer de energie van de onafhankelijk geoptimaliseerde monomeer.
• Elektronische Structuur Analyse: DFT-berekeningen werden uitgevoerd op geoptimaliseerde monomeren met behulp van GPAW met de PBE exchange-correlation functionaal en Projector Augmented Wave (PAW) formalisme op een real-space grid (0,20 Å afstand).
• Descriptoren:
  • Elektrostatica: Het Moleculaire Elektrostatische Potentiaal (ESP) werd gemapt op elektronendichtheidsoppervlakken (0,001 e Bohr$^{-3}$) en geanalyseerd via middelcentrering om relatieve landschappen te visualiseren.
  • Ladingverdeling: Hirshfeld populatieanalyse werd gebruikt om atoom-opgeloste ladingverdelingen te kwantificeren, inclusief statistische metrieken (gemiddelde, mediaan, standaarddeviatie) en residue-opgeloste sommen.
  • Structureel/Sequentie: Descriptoren omvatten sequence charge decoration (SCD), lokale hydropathie, hydrofobe run-lengte, radius van gyratie ($R_g$), end-to-end afstand en asfericiteit.
  • Interface: Intermoleculaire contacten werden gekwantificeerd door zware atoomafstanden ($\le$ 4,5 Å en 5,0 Å) en specifieke N/O/S contacten ($\le$ 3,5 Å).

Belangrijkste Resultaten
De studie laat zien dat de herschikking van residuen leidt tot onderscheidende, peptide-specifieke responsen in zowel energetische als structurele domeinen, ondanks een geconserveerde aminozuurcompositie.

1. Monomeer Energetica: De volgorde van residuen veranderde de elektronische energie van de geïsoleerde monomeren significant. Voor A$\beta$ was de gescrambelde controle ongeveer 45,33 kJ/mol lager in energie dan de aggregatiegevoelige sequentie. Daarentegen was de $\alpha$-synucleïne aggregatiegevoelige sequentie ongeveer 151,12 kJ/mol lager in energie dan de gescrambelde controle.
2. Associatie Energetica versus Interface Pakking:
   • A$\beta$: Zowel de aggregatiegevoelige sequentie (115,617 kJ/mol) als de gescrambelde controle (232,552 kJ/mol) vertoonden positieve associatie-energieën, wat aangeeft dat de geoptimaliseerde dimeren een hogere elektronische energie hadden dan de onafhankelijk geoptimaliseerde monomeren. De aggregatiegevoelige sequentie vertoonde echter een lagere positieve associatie-energie (met 116,935 kJ/mol) dan de gescrambelde controle. Paradoxaal genoeg vormde deze sequentie minder zware atoomcontacten (40 vs. 81) en had een lagere contactdichtheid dan de controle.
   • $\alpha$-Synucleïne: De gescrambelde controle vertoonde een gunstigere (negatievere) associatie-energie (-494,371 kJ/mol) dan de aggregatiegevoelige sequentie (-172,819 kJ/mol), een verschil van 321,552 kJ/mol. Echter, de aggregatiegevoelige sequentie vormde een uitgebreider interface, met 166 zware atoomcontacten (vs. 105) en een hogere dichtheid van nauwe N/O/S contacten (15 vs. 5).
   • Conclusie: Associatie-energie en interface-pakking (contactaantal/dichtheid) variëren niet lineair; ze zijn complementaire, en niet equivalente, maten van associatie.
3. Elektrostatische Landschappen (ESP): De herschikking van residuen reorganiseerde de ruimtelijke distributie van de ESP. Hoewel de A$\beta$- en $\alpha$-synucleïne-paren vergelijkbare globale ladingstoestanden deelden, verschoof de specifieke ruimtelijke organisatie van positieve en negatieve potentiaalgebieden aanzienlijk tussen de aggregatiegevoelige en gescrambelde sequenties. Dit geeft aan dat de sequentievolgorde de 3D-presentatie van het elektrostatische karakter dicteert.
4. Hirshfeld Ladingverdelingen: Atoom-ladinganalyse toonde aan dat de herschikking van residuen primair de atomaire ladingpopulaties herverdeelde in plaats van de globale ladingmaten uniform te verschuiven. In $\alpha$-synucleïne bleef de proportie positieve/negatieve atomen bijna constant, maar de waarschijnlijkheidsverdelingen van de ladingswaarden verschoven. In A$\beta$ werd naast de herverdeling ook een bescheiden verschuiving in de balans van positieve/negatieve atomen waargenomen.
5. Sequentie en Structurele Descriptoren:
   • A$\beta$: Herschikking veranderde de sequence charge decoration (SCD), de hydrofobe run-lengte (5 residuen vs. 3) en de moleculaire compactheid, waarbij de aggregatiegevoelige peptide een $R_g$ van ongeveer 10,1 Å vertoonde vergeleken met ongeveer 11,0 Å voor de gescrambelde controle.
   • $\alpha$-Synucleïne: Herschikking veranderde primair de lokale hydropathie (de maximale 3-residue gemiddelde hydropathie steeg van ~2,6 naar ~4,2 in de controle), terwijl SCD, $R_g$ en asfericiteit vergelijkbaar bleven.

Betekenis en Claims
Het artikel stelt vast dat de volgorde van residuen een fundamentele determinant is van vroege peptide-associatie, onafhankelijk van de aminozuurcompositie. De studie demonstreert dat identieke composities unieke energetische, structurele en elektrostatische kenmerken kunnen genereren, enkel door de sequentie-afhankelijke moleculaire organisatie.

Kernclaims zijn onder meer:

• Sequentie als Determinant: De arrangement van residuen, en niet alleen de compositie, beheerst de balans van intermoleculaire interacties, wat zowel de thermodynamica van associatie als de geometrie van de resulterende interface beïnvloedt.
• Peptide-Specifieke Responsen: De moleculaire consequenties van residu-herschikking zijn niet universeel; $\alpha$-synucleïne en A$\beta$ vertonen contrasterende responsen waarbij veranderingen in sequentievolgorde leiden tot verschillende combinaties van energetische gunstigheid en interface-pakking.
• Complementaire Descriptoren: De studie benadrukt dat associatie-energie en interface-contactmetrieken verschillende aspecten van moleculaire associatie vastleggen. Een meer gepopuleerde interface correleert niet noodzakelijkerwijs met een gunstigere associatie-energie, zoals te zien is aan de divergente trends tussen de twee peptide-systemen.
• Elektrostatische Reorganisatie: Identieke aminozuurcomposities kunnen unieke elektrostatische omgevingen opleveren vanwege de reorganisatie van ladingpatronen en 3D-structuur, wat een moleculair fundament biedt voor sequentie-specifieke associatiegedragingen.

Het werk biedt een kwantitatief kader voor het identificeren van sequentie-gecodeerde moleculaire determinanten, en suggereert dat de "code" voor vroege aggregatie ligt in de specifieke lineaire ordening van residuen, die de resulterende energetische en elektrostatische landschappen dicteert.