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

Technisches Resümee: Quantenchemische Modellierung der sequenzkodierten Energetik und elektrostatischen Landschaften in der frühen Peptidassoziation

Problemstellung
Proteinfehlfaltung und -aggregation stehen im Zentrum neurodegenerativer Erkrankungen wie der Parkinson-Krankheit (PD) und der Alzheimer-Krankheit (AD), die durch die frühe molekulare Assoziation von Peptiden wie $\alpha$-Synuclein und Amyloid-$\beta$ (A$\beta$) vorangetrieben werden. Während experimentelle Studien die Kinetik und Strukturen reifer Fibrillen charakterisiert haben, bleiben die energetischen Beiträge einzelner molekularer Konfigurationen während der transienten, heterogenen frühen Stadien der Assoziation unzureichend aufgelöst. Es besteht eine kritische Lücke im Verständnis darüber, wie die spezifische lineare Anordnung der Reste (Sequenz) – im Gegensatz zur gesamten Aminosäurezusammensetzung – die molekularen Eigenschaften beeinflusst, welche die anfänglichen Assoziationsereignisse steuern. Diese Studie zielt darauf ab, den Beitrag der Sequenzanordnung zu isolieren, indem aggregationsanfällige Peptidsegmente mit kompositionsgleichen, geskrambelten Kontrollen verglichen werden, um zu bestimmen, wie die Reihenfolge der Reste die Energetik, die strukturelle Organisation und die elektrostatischen Landschaften beeinflusst.

Methodik
Die Studie verwendet einen konsistenten computergestützten Workflow, der semiempirische Quantenchemie und Dichtefunktionaltheorie (DFT) integriert, um Monomer- und Homodimerzustände zu modellieren.

• Systeme: Aggregationsanfällige Segmente ($\alpha$-Synuclein: VTGVTAVAQK; A$\beta$: LVFFAEDVGS) wurden gegenüber kompositionsgleichen geskrambelten Kontrollen (AGATQTVVVK und EFLVGSDVAF, jeweils) verglichen.
• Geometrieoptimierung: Monomer- und Homodimerstrukturen wurden mittels der GFN2-xTB-Methode (Geometry, Frequency, Noncovalent extended Tight Binding, zweite Generation), implementiert in xTB, optimiert. Dimerstrukturen wurden initial mittels des GFN-Kraftfelds relaxiert, bevor die abschließende Optimierung erfolgte.
• Energetik: Assoziationsenergien ($\Delta E_{assoc}$) wurden als Differenz zwischen der konvergierten GFN2-xTB-Elektronenenergie des optimierten Dimers und der zweifachen Energie des unabhängig optimierten Monomers berechnet.
• Elektronische Strukturanalyse: DFT-Berechnungen wurden an optimierten Monomeren mittels GPAW mit dem PBE-Austausch-Korrelationsfunktional und dem Projector Augmented Wave (PAW)-Formalismus auf einem Realraum-Gitter (0,20 Å Abstand) durchgeführt.
• Deskriptoren:
  • Elektrostatik: Das molekulare elektrostatische Potenzial (ESP) wurde auf Elektronendichteflächen (0,001 e Bohr$^{-3}$) abgebildet und mittels Mittelzentrierung analysiert, um relative Landschaften zu visualisieren.
  • Ladungsverteilung: Die Hirshfeld-Populationsanalyse wurde verwendet, um atomaufgelöste Ladungsverteilungen zu quantifizieren, einschließlich statistischer Metriken (Mittelwert, Median, Standardabweichung) und residuenaufgelöster Summen.
  • Struktur/Sequenz: Zu den Deskriptoren gehörten die Sequenz-Ladungsdekoration (SCD), lokale Hydropathie, hydrophobe Run-Länge, der Trägheitsradius ($R_g$), der End-zu-Ende-Abstand und die Asphärizität.
  • Interface: Intermolekulare Kontakte wurden über Schweratome-Distanzen ($\le$ 4,5 Å und 5,0 Å) sowie spezifische N/O/S-Kontakte ($\le$ 3,5 Å) quantifiziert.

Kernergebnisse
Die Studie zeigt, dass die Umordnung der Reste distinkte, peptidspezifische Antworten in sowohl energetischen als auch strukturellen Domänen hervorruft, trotz konservierter Aminosäurezusammensetzung.

1. Monomer-Energetik: Die Reihenfolge der Reste veränderte die Elektronenenergie der isolierten Monomere signifikant. Für A$\beta$ war die geskrambelte Kontrolle um ~45,33 kJ/mol niedriger in der Energie als die aggregationsanfällige Sequenz. Im Gegensatz dazu war für $\alpha$-Synuclein die aggregationsanfällige Sequenz um ~151,12 kJ/mol niedriger in der Energie als seine geskrambelte Kontrolle.
2. Assoziationsenergetik vs. Interface-Packing:
   • A$\beta$: Sowohl die aggregationsanfällige Sequenz (115,617 kJ/mol) als auch die geskrambelte Kontrolle (232,552 kJ/mol) wiesen positive Assoziationsenergien auf, was bedeutet, dass die optimierten Dimere eine höhere Elektronenenergie aufwiesen als die unabhängig optimierten Monomere. Die aggregationsanfällige Sequenz zeigte jedoch eine geringere positive Assoziationsenergie (um 116,935 kJ/mol) als die geskrambelte Kontrolle. Paradoxerweise bildete diese Sequenz weniger Schweratom-Kontakte (40 vs. 81) und hatte eine geringere Kontaktdichte als die Kontrolle.
   • $\alpha$-Synuclein: Die geskrambelte Kontrolle wies eine günstigere (negativere) Assoziationsenergie auf (-494,371 kJ/mol) als die aggregationsanfällige Sequenz (-172,819 kJ/mol), mit einer Differenz von 321,552 kJ/mol. Dennoch bildete die aggregationsanfällige Sequenz ein ausgedehnteres Interface mit 166 Schweratom-Kontakten (vs. 105) und einer höheren Dichte an engen N/O/S-Kontakten (15 vs. 5).
   • Fazit: Assoziationsenergie und Interface-Packing (Kontaktanzahl/Dichte) variieren nicht linear; sie stellen komplementäre, nicht äquivalente Maße der Assoziation dar.
3. Elektrostatische Landschaften (ESP): Die Umordnung der Reste reorganisierte die räumliche Verteilung des ESP. Während A$\beta$- und $\alpha$-Synuclein-Paare ähnliche globale Ladungszustände teilten, verschob sich die spezifische räumliche Organisation der positiven und negativen Potenzialregionen signifikant zwischen aggregationsanfälligen und geskrambelten Sequenzen. Dies deutet darauf hin, dass die Sequenzordnung die 3D-Präsentation des elektrostatischen Charakters diktiert.
4. Hirshfeld-Ladungsverteilungen: Die Analyse der Atomladungen zeigte, dass die Umordnung der Reste primär die Verteilung der atomaren Ladungspopulationen veränderte, statt die globalen Ladungsmagnituden einheitlich zu verschieben. Bei $\alpha$-Synuclein blieb der Anteil der positiven/negativen Atome nahezu konstant, aber die Wahrscheinlichkeitsverteilungen der Ladungswerte verschoben sich. Bei A$\beta$ wurde neben der Umverteilung auch eine moderate Verschiebung im Gleichgewicht der positiven/negativen Atome beobachtet.
5. Sequenz- und Strukturdeskriptoren:
   • A$\beta$: Die Umordnung veränderte die Sequenz-Ladungsdekoration (SCD), die hydrophobe Run-Länge (5 Reste vs. 3) und die molekulare Kompaktheit, wobei das aggregationsanfällige Peptid einen $R_g$ von etwa 10,1 Å aufwies im Vergleich zu etwa 11,0 Å für die geskrambelte Kontrolle.
   • $\alpha$-Synuclein: Die Umordnung veränderte primär die lokale Hydropathie (die maximale 3-Rest-Mittelhydropathie stieg von ~2,6 auf ~4,2 in der Kontrolle), während SCD, $R_g$ und Asphärizität vergleichbar blieben.

Bedeutung und Behauptungen
Die Arbeit etabliert, dass die Reihenfolge der Reste ein fundamentaler Determinant der frühen Peptidassoziation ist, unabhängig von der Aminosäurezusammensetzung. Die Studie demonstriert, dass identische Zusammensetzungen allein durch die sequenzabhängige molekulare Organisation distinkte energetische, strukturelle und elektrostatische Charakteristika erzeugen können.

Zentrale Behauptungen sind:

• Sequenz als Determinant: Die Anordnung der Reste, nicht nur die Zusammensetzung, steuert das Gleichgewicht der intermolekularen Wechselwirkungen und beeinflusst sowohl die Thermodynamik der Assoziation als auch die Geometrie des resultierenden Interfaces.
• Peptidspezifische Antworten: Die molekularen Konsequenzen der Residuen-Umordnung sind nicht universell; $\alpha$-Synuclein und A$\beta$ zeigen kontrastierende Reaktionen, bei denen Änderungen der Sequenzordnung zu unterschiedlichen Kombinationen aus energetischer Begünstigung und Interface-Packing führen.
• Komplementäre Deskriptoren: Die Studie hebt hervor, dass Assoziationsenergie und Interface-Kontaktmetriken unterschiedliche Aspekte der molekularen Assoziation erfassen. Ein stärker besetztes Interface korreliert nicht zwangsläufig mit einer günstigeren Assoziationsenergie, wie die divergierenden Trends in den beiden Peptidsystemen zeigen.
• Elektrostatische Reorganisation: Identische Aminosäurezusammensetzungen können aufgrund der Reorganisation der Ladungsmuster und der 3D-Struktur distinkte elektrostatische Umgebungen hervorbringen, was eine molekulare Basis für sequenzspezifische Assoziationsverhalten liefert.

Die Arbeit liefert einen quantitativen Rahmen zur Identifizierung sequenzkodierter molekularer Determinanten und legt nahe, dass der „Code“ für die frühe Aggregation in der spezifischen linearen Ordnung der Reste liegt, welche die resultierenden energetischen und elektrostatischen Landschaften diktiert.