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

Résumé Technique : Modélisation de la Chimie Quantique des Énergétiques et des Paysages Électrostatiques Codés par la Séquence dans l'Association Précoce des Peptides

Énoncé du Problème
Le repliement incorrect et l'agrégation des protéines sont au cœur des troubles neurodégénératifs tels que la maladie de Parkinson (MP) et la maladie d'Alzheimer (MA), pilotés par l'association moléculaire précoce de peptides comme l'$\alpha$-synucléine et la protéine amyloïde $\beta$ (A$\beta$). Bien que les études expérimentales aient caractérisé la cinétique et les structures des fibrilles matures, les contributions énergétiques des configurations moléculaires individuelles lors des étapes transitoires et hétérogènes de l'association précoce restent mal résolues. Une lacune critique existe dans la compréhension de la manière dont l'ordre spécifique des résidus (la séquence), distinct de la composition globale en acides aminés, façonne les propriétés moléculaires régissant ces événements d'association initiaux. Cette étude vise à isoler la contribution de l'arrangement de la séquence en comparant des segments peptidiques sujets à l'agrégation avec des contrôles mélangés (scrambled) de composition identique, afin de déterminer comment l'ordre des résidus influence l'énergie, l'organisation structurelle et les paysages électrostatiques.

Méthodologie
L'étude emploie un flux de travail computationnel cohérent intégrant la chimie quantique semi-empirique et la théorie de la fonctionnelle de la densité (DFT) pour modéliser les états monomère et homodimère.

• Systèmes : Des segments sujets à l'agrégation ($\alpha$-synucléine : VTGVTAVAQK ; A$\beta$ : LVFFAEDVGS) ont été comparés à des contrôles mélangés de composition identique (AGATQTVVVK et EFLVGSDVAF, respectivement).
• Optimisation de la Géométrie : Les structures monomères et homodimères ont été optimisées via la méthode GFN2-xTB (Geometry, Frequency, Noncovalent extended Tight Binding, deuxième génération) implémentée dans xTB. Les structures de dimères ont été initialement relaxées en utilisant le champ de force GFN avant l'optimisation finale.
• Énergétique : Les énergies d'association ($\Delta E_{assoc}$) ont été calculées comme la différence entre l'énergie électronique du dimère optimisé par GFN2-xTB et deux fois l'énergie du monomère optimisé indépendamment.
• Analyse de la Structure Électronique : Les calculs DFT ont été effectués sur les monomères optimisés avec GPAW, utilisant la fonctionnelle d'échange-corrélation PBE et le formalisme de l'onde projetée (PAW) sur une grille d'espace réel (espacement de 0,20 Å).
• Descripteurs :
  • Électrostatique : Le potentiel électrostatique moléculaire (ESP) a été cartographié sur les surfaces de densité électronique (0,001 e Bohr$^{-3}$) et analysé via une centrage de la moyenne pour visualiser les paysages relatifs.
  • Distribution de Charge : L'analyse de population de Hirshfeld a été utilisée pour quantifier les distributions de charges atomiques, incluant les métriques statistiques (moyenne, médiane, écart-type) et les sommes par résidu.
  • Structure/Séquence : Les descripteurs comprenaient la décoration de charge de séquence (SCD), l'hydropathie locale, la longueur de la course hydrophobe, le rayon de giration ($R_g$), la distance bout-à-bout et l'asphéricité.
  • Interface : Les contacts intermoléculaires ont été quantifiés par les distances entre atomes lourds ($\le$ 4,5 Å et 5,0 Å) et les contacts spécifiques N/O/S ($\le$ 3,5 Å).

Résultats Clés
L'étude révèle que le réarrangement des résidus produit des réponses distinctes et spécifiques aux peptides, tant dans les domaines énergétiques que structurels, malgré une composition en acides aminés conservée.

1. Énergétique des Monomères : L'ordre des résidus a modifié de manière significative l'énergie électronique des monomères isolés. Pour l'A$\beta$, le contrôle mélangé était plus bas en énergie d'environ 45,33 kJ/mol par rapport à la séquence sujette à l'agrégation. Inversement, pour l'$\alpha$-synucléine, la séquence sujette à l'agrégation était plus basse en énergie d'environ 151,12 kJ/mol par rapport à son contrôle mélangé.
2. Énergie d'Association vs Empilement de l'Interface :
   • A$\beta$ : La séquence sujette à l'agrégation (115,617 kJ/mol) ainsi que le contrôle mélangé (232,552 kJ/mol) présentaient tous deux des énergies d'association positives, indiquant que les dimères optimisés étaient plus hauts en énergie électronique que les monomères optimisés indépendamment. Cependant, la séquence sujète à l'agrégation présentait une énergie d'association positive plus faible (de 116,935 kJ/mol) que le contrôle mélangé. Paradoxalement, cette séquence formait moins de contacts entre atomes lourds (40 contre 81) et présentait une densité de contact plus faible que le contrôle.
   • $\alpha$-Synucléine : Le contrôle mélangé présentait une énergie d'association plus favorable (plus négative) (-494,371 kJ/mol) que la séquence sujette à l'agrégation (-172,819 kJ/mol), soit une différence de 321,552 kJ/mol. Cependant, la séquence sujette à l'agrégation formait une interface plus étendue, avec 166 contacts d'atomes lourds (contre 105) et une densité plus élevée de contacts proches N/O/S (15 contre 5).
   • Conclusion : L'énergie d'association et l'empilement de l'interface (nombre de contacts) ne varient pas linéairement ; ils représentent des mesures complémentaires, et non équivalentes, de l'association.
3. Paysages Électrostatiques (ESP) : Le réarrangement des résidus a réorganisé la distribution spatiale de l'ESP. Bien que les paires d'A$\beta$ et d'$\alpha$-synucléine partagent des états de charge globaux similaires, l'organisation spatiale spécifique des régions de potentiel positif et négatif a considérablement changé entre les séquences sujettes à l'agrégation et les séquences mélangées. Cela indique que l'ordre de la séquence dicte la présentation 3D du caractère électrostatique.
4. Distributions de Charge de Hirshfeld : L'analyse de la charge atomique a montré que le réarrangement des résidus redistribuait principalement les populations de charges atomiques plutôt que de déplacer uniformément les magnitudes de charge globale. Pour l'$\alpha$-synucléine, la proportion d'atomes positifs/négatifs est restée presque constante, mais les distributions de probabilité des valeurs de charge ont décalé. Pour l'A$\beta$, un léger décalage dans l'équilibre des atomes positifs/négatifs a été observé parallèlement à la redistribution.
5. Descripteurs de Séquence et de Structure :
   • A$\beta$ : Le réarrangement a modifié la décoration de charge de séquence (SCD), la longueur de la course hydrophobe (5 résidus contre 3) et la compacité moléculaire, le peptide sujet à l'agrégation présentant un $R_g$ d'environ 10,1 Å contre environ 11,0 Å pour le contrôle mélangé.
   • $\alpha$-Synucléine : Le réarrangement a principalement modifié l'hydropathie locale (l'hydropathie moyenne sur 3 résidus est passée d'environ 2,6 à environ 4,2 pour le contrôle), tandis que la SCD, le $R_g$ et l'asphéricité sont restés comparables.

Signification et Revendications
L'article établit que l'ordre des résidus est un déterminant fondamental de l'association précoce des peptides, indépendamment de la composition en acides aminés. L'étude démontre que des compositions identiques peuvent générer des caractéristiques énergétiques, structurelles et électrostatiques distinctes uniquement par l'organisation moléculaire dépendante de la séquence.

Les principales revendications sont :

• La Séquence comme Déterminant : L'arrangement des résidus, et non seulement la composition, régit l'équilibre des interactions intermoléculaires, influençant à la fois la thermodynamique de l'association et la géométrie de l'interface résultante.
• Réponses Spécifiques aux Peptides : Les conséquences moléculaires du réarrangement des résidus ne sont pas universelles ; l'$\alpha$-synucléine et l'A$\beta$ présentent des réponses contrastées où les changements d'ordre de séquence mènent à différentes combinaisons de favorabilité énergétique et d'empilement d'interface.
• Descripteurs Complémentaires : L'étude souligne que l'énergie d'association et les métriques de contact d'interface capturent différents aspects de l'association moléculaire. Une interface plus peuplée ne corrèle pas nécessairement avec une énergie d'association plus favorable, comme on le voit dans les tendances divergentes des deux systèmes peptidiques.
• Réorganisation Électrostatique : Des compositions d'acides aminés identiques peuvent produire des environnements électrostatiques distincts en raison de la réorganisation de la structuration des charges et de la structure 3D, fournissant une base moléculaire aux comportements d'association spécifiques à la séquence.

Ce travail fournit un cadre quantitatif pour identifier les déterminants moléculaires codés par la séquence, suggérant que le "code" de l'agrégation précoce réside dans l'ordre linéaire spécifique des résidus qui dicte les paysages énergétiques et électrostatiques résultants.