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

तकनीकी सारांश: प्रारंभिक पेप्टाइड जुड़ाव में अनुक्रम-कोडित ऊर्जा और इलेक्ट्रोस्टैटिक परिदृश्यों का क्वांटम रासायनिक मॉडलिंग

समस्या विवरण
प्रोटीन मिसफोल्डिंग और एकत्रीकरण (aggregation), पार्किंसंस रोग (PD) और अल्जाइमर रोग (AD) जैसे न्यूरोडीजेनेरेटिव विकारों के केंद्र में हैं, जो $\alpha$-सिन्यूक्लिन ( $\alpha$-synuclein) और अमाइलॉइड $\beta$ (A$\beta$) जैसे पेप्टाइड्स के प्रारंभिक आणविक जुड़ाव द्वारा संचालित होते हैं। जबकि प्रयोगात्मक अध्ययनों ने परिपक्व फाइब्रिल्स (fibrils) की गतिज (kinetics) और संरचनाओं को स्पष्ट किया है, जुड़ाव के क्षणिक, विषम प्रारंभिक चरणों के दौरान व्यक्तिगत आणविक विन्यासों के ऊर्जा योगदान अभी भी अपर्याूर्ण रूप से ज्ञात हैं। एक महत्वपूर्ण अंतराल यह समझने में मौजूद है कि कैसे अवशेषों (residues) का विशिष्ट रैखिक क्रम (अनुक्रम), समग्र अमीनो एसिड संरचना से अलग होकर, प्रारंभिक जुड़ाव घटनाओं को नियंत्रित करने वाले आणविक गुणों को आकार देता है। यह अध्ययन अनुक्रम व्यवस्था के योगदान को अलग करने के उद्देश्य से एकत्रीकरण-प्रवण (aggregation-prone) पेप्टाइड खंडों की तुलना संरचना-मिलान वाले स्कैम्बलड (scrambled) नियंत्रणों से करता है, जिससे यह निर्धारित किया जा सके कि अवशेष क्रम ऊर्जा, संरचनात्मक संगठन और इलेक्ट्रोस्टैटिक परिदृश्यों को कैसे प्रभावित करता है।

कार्यप्रणाली
यह अध्ययन मोनोमर और होमोडाइमर अवस्थाओं को मॉडल करने के लिए सेमी-एम्पिरिकल क्वांटम केमिस्ट्री और डेंसिटी फंक्शनल थ्योरी (DFT) को एकीकृत करते हुए एक सुसंगत कम्प्यूटेशनल वर्कफ़्लो का उपयोग करता है।

• प्रणालियाँ: एकत्रीकरण-प्रवण खंडों ($\alpha$-सिन्यूक्लिन: VTGVTAVAQK; A$\beta$: LVFFAEDVGS) की तुलना संरचना-मिलान वाले स्कैम्बलड नियंत्रणों (AGATQTVVVK और EFLVGSDVAF, क्रमशः) से की गई।
• ज्यामिति अनुकूलन (Geometry Optimization): xTB में कार्यान्वित GFN2-xTB विधि (Geometry, Frequency, Noncovalent extended Tight Binding, second generation) का उपयोग करके मोनोमर और होमोडाइमर संरचनाओं को अनुकूलित किया गया। डाइमर संरचनाओं को अंतिम अनुकूलन के लिए GFN फोर्स फील्ड का उपयोग करके प्रारंभ में शिथिल (relaxed) किया गया था।
• ऊर्जा (Energetics): जुड़ाव ऊर्जा ($\Delta E_{assoc}$) की गणना अनुकूलित डाइमर की अभिसरित (converged) GFN2-xTB इलेक्ट्रॉनिक ऊर्जा और दो स्वतंत्र रूप से अनुकूलित मोनोमर की ऊर्जा के अंतर के रूप में की गई।
• इलेक्ट्रॉनिक संरचना विश्लेषण: अनुकूलित मोनोमर्स पर वास्तविक-स्थान ग्रिड (0.20 Å स्पेसिंग) पर PBE एक्सचेंज-कोरिलेशन फंक्शनल और प्रोजेक्टर ऑगमेंटेड वेव (PAW) फॉर्मलिज्म का उपयोग करके GPAW के माध्यम से DFT गणनाएँ की गईं।
• डिस्क्रिप्टर (Descriptors):
  • इलेक्ट्रोस्टैटिक्स: मॉलिक्यूलर इलेक्ट्रोस्टैटिक पोटेंशियल (ESP) को इलेक्ट्रॉन घनत्व सतहों (0.001 e Bohr$^{-3}$) पर मैप किया गया और सापेक्ष परिदृश्यों को विज़ुअलाइज़ करने के लिए मीन-सेंटरिंग के माध्यम से विश्लेषित किया गया।
  • चार्ज वितरण: परमाणु-विशिष्ट चार्ज वितरण, जिसमें सांख्यिकीय मेट्रिक्स (माध्य, माध्यिका, मानक विचलन) और अवशेष-विशिष्ट योग शामिल हैं, को मापने के लिए हिरशफेल्ड (Hirshfeld) जनसंख्या विश्लेषण का उपयोग किया गया।
  • संरचनात्मक/अनुक्रम: डिस्क्रिप्टर्स में सीक्वेंस चार्ज डेकोरेशन (SCD), स्थानीय हाइड्रोपैथी, हाइड्रोफोबिक रन लेंथ, जिरोडियस ऑफ जॉयरेशन ($R_g$), एंड-टू-एंड दूरी और असफेरिसिटी (asphericity) शामिल थे।
  • इंटरफेस: इंटरमॉलिक्युलर संपर्कों को हेवी एटम दूरियों ($\le$ 4.5 Å और 5.0 Å) और विशिष्ट N/O/S संपर्कों ($\le$ 3.5 Å) द्वारा मात्रात्मक रूप से मापा गया।

प्रमुख परिणाम
अध्ययन से पता चलता है कि अवशेष पुनर्व्यवस्था, संरक्षित अमीनो एसिड संरचना के बावजूद, ऊर्जावान और संरचनात्मक डोमेन में विशिष्ट, पेप्टाइड-विशिष्ट प्रतिक्रियाएं उत्पन्न करती है।

1. मोनोमर ऊर्जा: अवशेष क्रम ने अलग किए गए मोनोमर्स की इलेक्ट्रॉनिक ऊर्जा को महत्वपूर्ण रूप से बदल दिया। A$\beta$ के लिए, स्कैम्बलड नियंत्रण, एकत्रीकण-प्रवण अनुक्रम की तुलना में ~45.33 kJ/mol कम ऊर्जा वाला था। इसके विपरीत, $\alpha$-सिन्यूक्लिन के लिए, एकत्रीकरण-प्रवण अनुक्रम स्कैम्बलड नियंत्रण की तुलना में ~151.12 kJ/mol कम ऊर्जा वाला था।
2. जुड़ाव ऊर्जा बनाम इंटरफेस पैकिंग:
   • A$\beta$: दोनों (एकत्रीकरण-प्रवण अनुक्रम [115.617 kJ/mol] और स्कैम्बलड नियंत्रण [232.552 kJ/mol]) ने सकारात्मक जुड़ाव ऊर्जा प्रदर्शित की, जो इंगित करता कि अनुकूलित डाइमर स्वतंत्र रूप से अनुकूलित मोनोमर्स की तुलना में उच्च इलेक्ट्रॉनिक ऊर्जा वाले थे। हालांकि, एकत्रीकरण-प्रवण अनुक्रम ने स्कैम्बलड नियंत्रण की तुलना में कम सकारात्मक जुड़ाव ऊर्जा (116.935 kJ/mol कम) दिखाई। विरोधाभासी रूप से, इस अनुक्रम ने नियंत्रण की तुलना में कम हेवी एटम संपर्क (40 बनाम 81) और कम संपर्क घनत्व बनाया।
   • $\alpha$-सिन्यूक्लिन: स्कैम्बलड नियंत्रण ने एकत्रीकरण-प्रवण अनुक्रम (-172.819 kJ/mol) की तुलना में अधिक अनुकूल (अधिक ऋणात्मक) जुड़ाव ऊर्जा (-494.371 kJ/mol) प्रदर्शित की, जिसका अंतर 321.552 kJ/mol था। हालांकि, एकत्रीकरण-प्रवण अनुक्रम ने अधिक विस्तृत इंटरफेस बनाया, जिसमें 166 हेवी एटम संपर्क (105 बनाम 166) और उच्च N/O/S संपर्क घनत्व (15 बनाम 5) था।
   • निष्कर्ष: जुड़ाव ऊर्जा और इंटरफेस पैकिंग (संपर्क संख्या/घनत्व) रैखिक रूप से भिन्न नहीं होते हैं; वे जुड़ाव के पूरक, न कि समतुल्य माप हैं।
3. इलेक्ट्रोस्टैटिक परिदृश्य (ESP): अवशेष पुनर्व्यवस्था ने ESP के स्थानिक वितरण को पुनर्गठित किया। जबकि A$\beta$ और $\alpha$-सिन्यूक्लिन युग्मों ने समान वैश्विक चार्ज अवस्था साझा की, धनात्मक और ऋणात्मक क्षमता क्षेत्रों का विशिष्ट स्थानिक संगठन स्कैम्बलड अनुक्रमों के बीच काफी बदल गया। यह इंगित करता है कि अनुक्रम क्रम इलेक्ट्रोस्टैटिक चरित्र के त्रि-आयामी प्रस्तुतीकरण को निर्धारित करता है।
4. हिरशफेल्ड चार्ज वितरण: परमाणु चार्ज विश्लेषण ने दिखाया कि अवशेष पुनर्व्यवस्था ने वैश्विक चार्ज परिमाणों को समान रूप से स्थानांतरित करने के बजाय मुख्य रूप से परमाणु चार्ज आबादी को पुनर्वितरित किया। $\alpha$-सिन्यूक्लिन में, धनात्मक/ऋणात्मक परमाणुओं का अनुपात लगभग स्थिर रहा, लेकिन चार्ज मानों के संभाव्यता वितरण बदल गए। A$\beta$ में, पुनर्वितरण के साथ धनात्मक/ऋणात्मक परमाणुओं के संतुलन में एक मामूली बदलाव देखा गया।
5. अनुक्रम और संरचनात्मक डिस्क्रिप्टर:
   • A$\beta$: पुनर्व्यवस्था ने सीक्वेंस चार्ज डेकोरेशन (SCD), हाइड्रोफोबिक रन लेंथ (5 अवशेष बनाम 3), और आणविक सघनता को परिवर्तित किया, जिसमें एकत्रीकरण-प्रवण पेप्टाइड का $R_g$ लगभग 10.1 Å था, जबकि स्कैम्बलड नियंत्रण के लिए यह लगभग 11.0 Å था।
   • $\alpha$-सिन्यूक्लिन: पुनर्व्यवस्था ने मुख्य रूप से स्थानीय हाइड्रोपैथी को परिवर्तित किया (अधिकतम 3-अवशेष माध्य हाइड्रोपैथी ~2.6 से बढ़कर नियंत्रण में ~4.2 हो गई), जबकि SCD, $R_g$, और असफेरिसिटी तुलनीय रहे।

महत्व और दावे
यह शोध स्थापित करता है कि अवशेष क्रम, अमीनो एसिड संरचना से स्वतंत्र, प्रारंभिक पेप्टाइड जुड़ाव का एक मौलिक निर्धारक है। अध्ययन दर्शाता है कि समान संरचनाएं केवल अनुक्रम-निर्भर आणविक संगठन के माध्यम से विशिष्ट ऊर्जावान, संरचनात्मक और इलेक्ट्रोस्टैटिक विशेषताओं को उत्पन्न कर सकती हैं।

मुख्य दावे निम्नलिखित हैं:

• अनुक्रम एक निर्धारक के रूप में: अवशेषों का क्रम, न कि केवल संरचना, अंतः आणविक अंतःक्रियाओं के संतुलन को नियंत्रित करता है, जो परिणामी इंटरफेस की ज्यामिति और थर्मोडायनामिक्स दोनों को प्रभावित करता है।
• पेप्टाइड-विशिष्ट प्रतिक्रियाएं: अवशेष पुनर्व्यवस्था के आणविक परिणाम सार्वभौमिक नहीं हैं; $\alpha$-सिन्यूक्लिन और A$\beta$ विपरीत प्रतिक्रियाएं प्रदर्शित करते हैं जहाँ अनुक्रम क्रम में परिवर्तन ऊर्जावान अनुकूलता और इंटरफेस पैकिंग के विभिन्न संयोजनों की ओर ले जाता है।
• पूरक डिस्क्रिप्टर: अध्ययन रेखांकित करता है कि जुड़ाव ऊर्जा और इंटरफेस संपर्क मेट्रिक्स आणविक जुड़ाव के विभिन्न पहलुओं को कैप्चर करते हैं। एक अधिक आबादी वाला इंटरफेस आवश्यक रूप से अधिक अनुकूल जुड़ाव ऊर्जा के साथ सहसंबंधित नहीं होता है, जैसा कि दोनों पेप्टाइड प्रणालियों में देखे गए विचलित रुझानों से पता चलता है।
• इलेक्ट्रोस्टैटिक पुनर्गठन: समान अमीनो एसिड संरचनाएं चार्ज पैटर्न और 3D संरचना के पुनर्गठन के कारण विशिष्ट इलेक्ट्रोस्टैटिक वातावरण प्रदान कर सकती हैं, जो अनुक्रम-विशिष्ट जुड़ाव व्यवहारों के लिए एक आणविक आधार प्रदान करती है।

यह कार्य उन मात्रात्मक ढांचे को प्रदान करता है जो अनुक्रम-एन्कोडेड आणविक निर्धारकों की पहचान करने में मदद करता है, यह सुझाव देते हुए कि प्रारंभिक एकत्रीकरण का "कोड" अवशेषों के विशिष्ट रैखिक क्रम में निहित है जो परिणामी ऊर्जावान और इलेक्ट्रोस्टैटिक परिदृश्यों को निर्धारित करता है।