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

ملخص تقني: النمذجة الكيميائية الكمومية للطاقات والآفاق الكهروستاتيكية المشفرة بالتسلسل في الارتباط المبكر للببتيدات

بيان المشكلة
يعد سوء طي البروتين وتجمعه من الأمور المركزية في الاضطرابات العصبية التنكسية مثل مرض باركنسون (PD) ومرض ألزهايمر (AD)، والتي تقودها عمليات الارتباط الجزيئي المبكرة للببتيدات مثل ألفا-سينوكلين (α-synuclein) وأميلويد بيتا (Aβ). وبينما شخصت الدراسات التجريبية الحركية والهياكل للألياف الناضجة، لا تزال المساهمات الطاقية للتكوينات الجزيئية الفردية خلال المراحل الانتقالية غير المتجانسة والمبكرة للارتباط غير واضحة بدقة. وتوجد فجوة حرجة في فهم كيفية تأثير الترتيب الخطي المحدد للبقايا (التسلسل)، بمعزل عن التركيب الكلي للأحماض الأمينية، في تشكيل الخصائص الجزيئية التي تحكم أحداث الارتباط الأولية هذه. تهدف هذه الدراسة إلى عزل مساهمة ترتيب التسلسل من خلال مقارنة قطع الببتيد القابلة للتجمع مع ضوابط مبعثرة متطابقة في التركيب، وذلك لتحديد كيفية تأثير ترتيب البقايا على الطاقات، والتنظيم الهيكلي، والآفاق الكهروستاتيكية.

المنهجية
تستخدم الدراسة سير عمل حوسبي متسق يدمج الكيمياء الكمومية شبه التجريبية ونظرية الدالة الوظيفية للكثافة (DFT) لنمذجة حالات المونومر (الجزيء المنفرد) والهوموديمر (الجزيء المزدوج المتماثل).

• الأنظمة: تمت مقارنة مقاطع معرضة للتجمع (ألفا-سينوكلين: VTGVTAVAQK؛ Aβ: LVFFAEDVGS) مقابل ضوابط مبعثرة متطابقة في التركيب (AGATQTVVVK و EFLVGSDVAF، على التوالي).
• تحسين الهندسة: تم تحسين هياكل المونومر والهوموديمر باستخدام طريقة GFN2-xTB (الهندسة، التردد، الارتباط غير التساهمي الممتد، الجيل الثاني) المطبقة في xTB. تم استرخاء هياكل الديمر أولياً باستخدام مجال القوة GFN قبل التحسين النهائي.
• الطاقات: تم حساب طاقات الارتباط ($\Delta E_{assoc}$) كفرق بين الطاقة الإلكترونية المتجمعة لـ GFN2-xTB للديمر المحسن وضِعف طاقة المونومر المحسن بشكل مستقل.
• تحليل البنية الإلكترونية: أجريت حسابات DFT على المونومرات المحسنة باستخدام GPAW مع دالة التبادل والارتباط PBE وصيغة الموجة المستوية (PAW) على شبكة في الفضاء الحقيقي (بفاصل 0.20 أنغستروم).
• الموصِّفات:
  • الكهروستاتيكا: تم رسم جهد الكهروستاتيكا الجزيئي (ESP) على أسطح كثافة الإلكترونات (0.001 e Bohr$^{-3}$) وتحليله عبر تمركُز المتوسط لتصور الآفاق النسبية.
  • توزيع الشحنة: استُخدم تحليل تعداد هيرشفيلد (Hirshfeld) لتحديد كمية توزيعات الشحنة على مستوى الذرة، بما في ذلك المقاييس الإحصائية (المتوسط، الوسيط، الانحراف المعياري) والمجاميع المحددة لكل بقايا.
  • الهيكل/التسلسل: شملت الموصِّفات زخرفة شحنة التسلسل (SCD)، والرهاب المائي المحلي، وطول سلسلة الرهاب المائي، ونصف قطر الدوران ($R_g$)، والمسافة من الطرف إلى الطرف، وعدم التماثل.
  • الواجهة: تم قياس الاتصالات بين الجزيئات عبر مسافات الذرات الثقيلة ($\le$ 4.5 Å و 5.0 Å) واتصالات N/O/S المحددة ($\le$ 3.5 Å).

النتائج الرئيسية
تكشف الدراسة أن إعادة ترتيب البقايا تنتج استجابات متميزة ومحددة لكل ببتيد في كل من المجالات الطاقية والهيكلية، رغم ثبات تركيب الأحماض الأمينية.

1. طاقة المونومر: غير ترتيب البقايا بشكل كبير الطاقة الإلكترونية للمونومرات المعزولة. بالنسبة لـ Aβ، كان التحكم المبعثر أقل في الطاقة بمقدار ~45.33 كيلوجول/مول مقارنة بالتسلسل المعرض للتجمع. وعلى العكس من ذلك، بالنسبة لـ ألفا-سينوكلين، كان التسلسل المعرض للتجمع أقل في الطاقة بمقدار ~151.12 كيلوجول/مول مقارنة بالتحكم المبعثر الخاص به.
2. طاقة الارتباط مقابل تعبئة الواجهة:
   • Aβ: أظهر كل من تسلسل Aβ المعرض للتجمع (115.617 كيلوجول/مول) والتحكم المبعثر (232.552 كيلوجول/مول) طاقات ارتباط موجبة، مما يشير إلى أن الديمرات المحسنة كانت أعلى في الطاقة الإلكترونية من المونومرات المحسنة بشكل مستقل. ومع ذلك، أظهر التسلسل المعرض للتجمع طاقة ارتباط موجبة أقل (بمقدار 116.935 كيلوجول/مول) من التحكم المبعثر. ومن المفارقات أن هذا التسلسل شكل اتصالات ذرية ثقيلة أقل (40 مقابل 81) وكثافة اتصالات أقل من التحكم.
   • ألفا-سينوكلين: أظهر التحكم المبعثر طاقة ارتباط أكثر تفضيلاً (أكثر سالبية) بلغت (-494.371 كيلوجول/مول) مقارنة بالتسلسل المعرض للتجمع (-172.819 كيلوجول/مول)، بفارق قدره 322.552 كيلوجول/مول. ومع ذلك، شكل التسلسل المعرض للتجمع واجهة أكثر اتساعاً، مع 166 اتصالاً من الذرات الثقيلة (مقابل 105) وكثافة أعلى من اتصالات N/O/S الوثيقة (15 مقابل 5).
   • الاستنتاج: طاقة الارتباط ومقاييس تعبئة الواجهة لا تتغير خطياً؛ فهما يمثلان مقاييس متكاملة وليست متكافئة للارتباط.
3. الآفاق الكهروستاتيكية (ESP): أعادت إعادة ترتيب البقية تنظيم التوزيع المكاني لـ ESP. وبينما تشاركت أزواج Aβ وألفا-سينوكلين في حالات شحنة عالمية متشابهة، فإن التنظيم المكاني لمناطق الجهد الموجب والسالب قد تغير بشكل كبير بين التسلسلات المعرضة للتجمع والمبعثرة. وهذا يشير إلى أن ترتيب التسلسل يحدد العرض ثلاثي الأبعاد للطابع الكهروستاتيكي.
4. توزيعات شحنة هيرشفيلد: أظهر تحليل الشحنة الذرية أن إعادة ترتيب البقايا أدت أساساً إلى إعادة توزيع تعدادات الشحنة الذرية بدلاً من إزاحة مقادير الشحنة العالمية بشكل موحد. في ألفا-سينوكلين، ظلت نسبة الذرات الموجبة/السالبة ثابتة تقريباً، لكن توزيعات قيم الشحنة قد تغيرت. وفي Aβ، لوحظت إزاحة طفيفة في توازن الذرات الموجبة/السالبة جنباً إلى جنب مع إعادة التوزيع.
5. الموصِّفات الهيكلية والتسلسلية:
   • Aβ: غيرت إعادة الترتيب من زخرفة شحنة التسلسل (SCD)، وطول سلسلة الرهاب المائي (5 بقايا مقابل 3)، وتراص الجزيء، حيث أظهر الببتيد المعرض للتجمع نصف قطر دوران ($R_g$) يبلغ حوالي 10.1 أنغستروم مقارنة بحوالي 11.0 أنغستروم للتحكم المبعثر.
   • ألفا-سينوكلين: غيرت إعادة الترتيب بشكل أساسي الرهاب المائي المحلي (ارتفع متوسط الرهاب المائي لـ 3 بقايا من ~2.6 إلى ~4.2 في التحكم)، بينما ظلت SCD و $R_g$ وعدم التماثل متقاربة.

الأهمية والادعاءات
تثبت الورقة أن ترتيب البقايا هو محدد أساسي للارتباط المبكر للببتيد، بشكل مستقل عن تركيب الأحماض الأمينية. وتوضح الدراسة أن الأحماض الأمينية ذات التركيب المتطابق يمكن أن تولد خصائص طاقية وهيكلية وكهروستاتيكية متميزة فقط من خلال التنظيم الجزيئي المعتمد على التسلسل.

تشمل الادعاءات الرئيسية ما يلي:

• التسلسل كمحدد: إن ترتيب البقايا، وليس مجرد التركيب، هو ما يحكم توازن التفاعلات بين الجزيئات، مما يؤثر على كل من ديناميكا الحرارية للارتباط وهندسة الواجهة الناتجة.
• الاستجابات النوعية للببتيد: إن النتائج الجزيئية لإعادة ترتيب البقايا ليست عالمية؛ حيث تظهر ألفا-سينوكلين و Aβ استجابات متضاربة حيث تؤدي التغييرات في ترتيب التسلسل إلى مزيج مختلف من التفضيل الطاقي وتعبئة الواجهة.
• الموصِّفات المتكاملة: تسلط الدراسة الضوء على أن طاقة الارتباط ومقاييس اتصال الواجهة يلتقطان جوانب مختلفة من الارتباط الجزيئي. فواجهة أكثر ازدحاماً لا تعني بالضرورة طاقة ارتباط أكثر تفضيلاً، كما لوحظ في الاتجاهات المتباعدة لكلا النظامين الببتيديين.
• إعادة التنظيم الكهروستاتيكي: يمكن للتركيبات المتطابقة من الأحماض الأمينية أن تنتج بيئات كهروستاتيكية متميزة بسبب إعادة تنظيم أنماط الشحنة والهيكل ثلاثي الأبعاد، مما يوفر أساساً جزيئياً لسلوكيات الارتباط المعتمدة على التسلسل.

يوفر هذا العمل إطاراً كمياً لتحديد المحددات الجزيئية المشفرة بالتسلسل، مما يشير إلى أن "الشفرة" للتجمع المبكر تكمن في الترتيب الخطي المحدد للبقايا الذي يحدد الآفاق الطاقية والكهروستاتيكية الناتجة.