⚛️ biophysics

Extending the MARTINI 3 Coarse-Grained Forcefield to Polypeptoids

यह अध्ययन 19 सामान्य पेप्टॉइड अवशेषों के लिए पहला MARTINI 3-अनुकूल कोर्स-ग्रेन्ड फोर्सफील्ड प्रस्तुत करता है, जिसे ऑल-एटम संदर्भ सिमुलेशन से प्राप्त किया गया है और मार्टिनाइज़2 (martinize2) टूल में एकीकृत किया गया है, ताकि पेप्टॉइड संरचनाओं, स्व-संयोजन और अंतःक्रियाओं के कुशल और सटीक मेसोस्केल मॉडलिंग को सक्षम बनाया जा सके।

मूल लेखक: Wang, J., Yu, Z., Zhao, M.

प्रकाशित 2026-04-11

📖 5 मिनट में पढ़ें🧠 गहराई से पढ़ें

CC BY 4.0

मूल लेखक: Wang, J., Yu, Z., Zhao, M.

मूल पेपर CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/) के तहत लाइसेंस किया गया है। ⚕️ यह एक ऐसे प्रीप्रिंट की AI से तैयार की गई व्याख्या है जिसकी अभी सहकर्मी समीक्षा नहीं हुई है। यह चिकित्सकीय सलाह नहीं है। इस सामग्री के आधार पर स्वास्थ्य संबंधी फैसले न लें। पूरा डिस्क्लेमर पढ़ें

कल्पना कीजिए कि आप यह समझने की कोशिश कर रहे हैं कि एक जटिल मशीन कैसे काम करती है। आपके पास दो विकल्प हैं: या तो आप हर एक पेंच, गियर और तार को देख सकते हैं (एक "ऑल-एटम" दृश्य), या आप मशीन को बड़े, कार्यात्मक ब्लॉकों के एक संग्रह के रूप में देख सकते हैं (एक "कोर्स-ग्रेन्ड" दृश्य)।

यह शोध पत्र एक विशिष्ट प्रकार के सिंथेटिक पदार्थ जिसे पॉलीपेप्टॉइड्स (Polypeptoids) कहा जाता है, के लिए एक नया, बेहतर "कार्यात्मक ब्लॉक" बनाने के बारे में है।

यहाँ शोधकर्ताओं ने जो किया, उसे सरल भाषा में समझाया गया है:

1. समस्या: "स्लो-मोशन" सामग्री

पॉलीपेप्टॉइड्स प्राकृतिक प्रोटीनों की तरह होते हैं। ये स्मार्ट सामग्री, ड्रग डिलीवरी सिस्टम और सूक्ष्म नैनोवायर बनाने के लिए उपयोग किए जाने वाले सिंथेटिक चेन हैं।

चुनौती: प्राकृतिक प्रोटीनों के विपरीत, पॉलीपेप्टॉइड्स बहुत लचीले और घुमावदार होते हैं। उनके बैकबोन (रीढ़) में एक विशेष "हिंज" (कब्जा) होता है जो आगे-पीछे झूल सकता है (जैसे एक दरवाजा खुलता और बंद होता है)।
सिमुलेशन का दुस्वप्न: इन सामग्रियों का अध्ययन करने के लिए कंप्यूटर पर, वैज्ञानिक आमतौर पर "ऑल-एटम" मॉडल का उपयोग करते हैं। यह एक ऐसी फिल्म देखने जैसा है जहाँ हर एक फ्रेम को हाथ से बनाया जाता है। यह अविश्वसनीय रूप से सटीक है, लेकिन इसमें बहुत समय लगता है। क्योंकि पॉलीपेप्टॉइड के हिंज बहुत धीरे-धीरे घूमते हैं, एक कंप्यूटर सिमुलेशन को अणु के एक बार मुड़ने तक देखने के लिए हफ्तों तक चलना पड़ सकता है। यह घोंखे की दौड़ को गिनने के लिए उसके हर कदम को गिनने जैसा है।

2. समाधान: "MARTINI 3" टूलकिट

वैज्ञानिकों के पास एक लोकप्रिय टूलकिट है जिसे MARTINI कहा जाता है, जो अणुओं को सरल बनाता है। अणुओं के हर परमाणु को बनाने के बजाय, यह उन्हें "बीड्स" (मोतियों की तरह) में समूहबद्ध करता है (जैसे लेगो ब्लॉक्स)। यह सिमुलेशन को हजारों गुना तेज़ बना देता है।

अंतराल: अब तक, पॉलीपेप्टॉइड्स के लिए कोई अच्छा "लेगो सेट" नहीं था जो टूलकिट के नवीनतम संस्करण (MARTINI 3) के साथ काम कर सके। मौजूदा सेट या तो बहुत भारी थे या अन्य सामग्रियों (जैसे सेल मेम्ब्रेन) के साथ फिट नहीं बैठते थे जिनका अध्ययन वैज्ञानिक इन पॉलीपेप्टॉइड्स के साथ करना चाहते थे।

3. सफलता: नया लेगो सेट बनाना

अलाबामा विश्वविद्यालय की टीम ने पॉलीपेप्टॉइड्स के लिए पहला MARTINI 3-संगत मॉडल बनाया। उन्होंने इसे इस प्रकार किया:

ब्लूप्रिंट (मैपिंग): उन्हें यह तय करना था कि एक जटिल पॉलीपेप्टॉइड अणु को सरल बीड्स में कैसे बदला जाए। उन्होंने अणुओं को समूहबद्ध करने के विभिन्न तरीकों का परीक्षण किया। कल्पना कीजिए कि आप एक सूटकेस पैक कर रहे हैं: आप कपड़ों को कसकर मोड़ सकते हैं (मैपिंग 1) या उन्हें ढीला भर सकते हैं (मैपिंग 2)। उन्होंने वह "कसकर फोल्ड" वाला तरीका खोजा जो अणु के आकार को सटीक रखते हुए उसकी लचीलेपन को खोए बिना काम करता है।
प्रशिक्षण (स्लो-मोशन से सीखना): यह सुनिश्चित करने के लिए कि उनका नया "फास्ट" मॉडल सटीक है, उन्होंने पहले सुपर-स्लो, सुपर-एक्यूरेट "ऑल-एटम" सिमुलेशन चलाए। उन्होंने पैरेलल बायस मेटाडायनामिक्स नामक एक विशेष ट्रिक का उपयोग किया (इसे एक सहायक मार्गदर्शक के रूप में सोचें जो अणु को उसके धीमे घुमावों के माध्यम से धकेलता है ताकि वे देख सकें कि इसके सभी संभावित आकार क्या हो सकते हैं)।
अनुवाद (बोल्ट्ज़मैन इन्वर्जन): उन्होंने उन धीमे सिमुलेशन से डेटा लिया और गणितीय रूप से उसे अपने नए तेज़ मॉडल के नियमों में "अनुवादित" किया। इसने यह सुनिश्चित किया कि जब तेज़ मॉडल चले, तो वह बिल्कुल उसी तरह व्यवहार करे जैसे धीमा, सटीक मॉडल करता है।

4. परिणाम: गति बनाम सटीकता

उन्होंने अपने नए मॉडल का परीक्षण चार तरीकों से "गोल्ड स्टैंडर्ड" (धीमे सिमुलेशन) के विरुद्ध किया:

आकार: क्या श्रृंखला सही ढंग से मुड़ती है या फैलती है? हाँ।
घुमाव: क्या बैकबोन सही ढंग से घूमता है? हाँ।
चिपकना: क्या दो श्रृंखलाएं सही तरीके से आपस में जुड़ती हैं? हाँ।
निर्माण: क्या वे आपस में मिलने पर सही आकार बनाती हैं? हाँ।

जादुई संख्या: नया मॉडल पुराने तरीके की तुलना में 57 गुना तेज़ है।

उपमा: यदि पुराने तरीके को आणविक गतिविधि के एक दिन के सिमुलेशन के लिए 57 घंटे लगते, तो नए तरीके में केवल 1 घंटा लगेगा। यह वैज्ञानिकों को इन सामग्रियों को असेंबल होते और आपस में क्रिया करते हुए वास्तविक समय में देखने की अनुमति देता है, न कि परिणामों के लिए महीनों इंतजार करने की।

5. यह क्यों मायने रखता है

यह केवल एक स्पीड अपग्रेड नहीं है; यह भविष्य की कुंजी है।

तर्कसंगत डिजाइन (Rational Design): वैज्ञानिक अब कंप्यूटर पर नए पॉलीपेप्टॉइड सामग्रियों को डिजाइन कर सकते हैं, तुरंत उनका परीक्षण कर सकते हैं, और लैब में रसायन मिलाने से पहले ही भविष्यवाणी कर सकते हैं कि वे कैसे व्यवहार करेंगे।
सार्वभौमिक भाषा: क्योंकि यह मॉडल मानक MARTINI 3 इकोसिस्टम में फिट बैठता है, आप अब एक ही सिमुलेशन में सेल मेम्ब्रेन, दवाओं या अन्य प्रोटीनों के साथ आसानी से पॉलीपेप्टॉइड्स को मिला सकते हैं। यह एक सार्वभौमिक भाषा बोलने जैसा है जो विभिन्न प्रकार के अणुओं को कंप्यूटर में एक-दूसरे से "बात" करने की अनुमति देती है।

सारांश

शोधकर्ताओं ने एक धीमी, कठिन अध्ययन वाली सामग्री (पॉलीपेप्टॉइड्स) ली, उसका एक सरल लेकिन अत्यधिक सटीक "लेगो" संस्करण बनाया, और इसे दुनिया के सबसे लोकप्रिय सिमुलेशन टूलकिट के अनुकूल बनाया। परिणाम एक ऐसा टूल है जो 57 गुना तेज़ है, जिससे वैज्ञानिक बहुत तेज़ी से अगली पीढ़ी की जीवन रक्षक दवाओं और स्मार्ट सामग्रियों को डिजाइन कर सकते हैं।

यहाँ शोध पत्र "Extending the MARTINI 3 Coarse-Grained Forcefield to Polypeptoids" का विस्तृत तकनीकी सारांश दिया गया है।

1. समस्या विवरण (Problem Statement)

पॉलीपेप्टॉइड्स (polypeptoids - poly-N-substituted glycines) सिंथेटिक पेप्टिडोमिमेटिक्स हैं जहाँ पार्श्व श्रृंखलाएं (side chains) $\alpha$ -कार्बन के बजाय बैकबोन एमाइड नाइट्रोजन से जुड़ी होती हैं। यह संरचनात्मक अंतर बैकबोन हाइड्रोजन बॉन्डिंग को समाप्त कर देता है, जिसके परिणामस्वरूप अद्वितीय विन्यास संबंधी लचीलापन (conformational flexibility) और धीमी cis/trans $\omega$ -डाइहेड्रल आइसोमेराइज़ेशन (isomerization) प्राप्त होता है। हालाँकि ऑल-एटम (AA) फोर्स फील्ड मौजूद हैं (जैसे, CHARMM-संगत MoSiC-CGenFF-NTOID), वे कम्प्यूटेशनल रूप से महंगे हैं और छोटे सिस्टम तथा कम समय सीमा (सब-माइक्रोसेकंड) तक ही सीमित हैं।

मौजूदा कोर्स-ग्रेन्ड (CG) मॉडल मुख्य रूप से MARTINI 2 फ्रेमवर्क या विशेष डिजाइनों तक सीमित हैं जिनमें आधुनिक बायोमॉलिकुलर इकोसिस्टम के साथ अनुकूलता का अभाव है। MARTINI 3 फोर्स फील्ड बेहतर इंटरेक्शन स्कीम्स और बीड डेफिनिशन प्रदान करता है, लेकिन इसे अभी तक पॉलीपेप्टॉइड्स के लिए विस्तारित नहीं किया गया है। एक MARTINI 3-संगत मॉडल विकसित करना पॉलीपेप्टॉइड स्व-संयोजन (self-assembly), झिल्ली अंतःक्रिया (membrane interactions), और बड़े पैमाने पर नैनोस्ट्रक्चर निर्माण के मेसोस्केल सिमुलेशन को सक्षम करने के लिए महत्वपूर्ण है, जो वर्तमान में AA सिमुलेशन के लिए अप्राप्य हैं।

2. कार्यप्रणाली (Methodology)

लेखकों ने एक कठोर पैरामीट्राइजेशन वर्कफ़्लो का उपयोग करके एक बॉटम-अप CG फोर्स फील्ड विकसित किया:

संदर्भ सिमुलेशन (All-Atom):
- सिस्टम: एसिटाइल (Acetyl) और एमीनोमिथाइल (Aminomethyl) समूहों के साथ कैप किए गए होमोपैप्टॉइड चेन (4 से 20 अवशेष)।
- फोर्स फील्ड: MoSiC-CGenFF-NTOID (CHARMM-संगत)।
- एन्हांस्ड सैंपलिंग: बैकबोन $\omega$ -डाइहेड्रल cis/trans आइसोमेराइज़ेशन की उच्च सक्रियण बाधा (~63 kJ/mol) को दूर करने के लिए पैरेलल बायस मेटाडायनामिक्स (PBMetaD) का उपयोग किया गया। इसने कन्वर्ज्ड सैंपलिंग सुनिश्चित की, जो मानक MD के साथ कठिन है।
मैपिंग रणनीति (Mapping Strategy):
- आणविक आयतन को संरक्षित करने के लिए सेंटर-ऑफ-जियोमेट्री (COG) मैपिंग योजना को अपनाया गया (जो MARTINI 3 के लिए मानक है), न कि सेंटर-ऑफ-मास को।
- बैकबोन: चार मैपिंग स्कीमा का मूल्यांकन किया गया। मैपिंग 1 को चुना गया, जिसने केवल भारी परमाणुओं (N, C $\alpha$ , C, O) का उपयोग करके C–C $\alpha$ बॉन्ड के साथ ज्यामितीय केंद्र को रखा। इसने हार्मोनिक पोटेंशियल के लिए उपयुक्त एक यूनिमोडल, सममित बॉन्ड लंबाई वितरण प्रदान किया।
- साइड चेन्स: इन्हें बिना नए बीड प्रकार पेश किए, मानक MARTINI 3 बीड प्रकारों (बैकबोन के लिए P2, साइड चेन्स के लिए विभिन्न वर्ग) में मैप किया गया, जिससे मौजूदा रासायनिक ऑन्टोलॉजी के साथ अनुकूलता बनी रही।
पैरामीट्राइजेशन (Parameterization):
- बॉन्डेड इंटरेक्शन: PBMetaD प्रक्षेपवक्र (trajectories) के संभाव्यता वितरण के डायरेक्ट बोल्ट्ज़मैन इनवर्जन (DBI) के माध्यम से प्राप्त किया गया।
  - एकीकृत दृष्टिकोण (Unified Approach): अधिकांश बैकबोन टर्म्स (बॉन्ड्स, एंगल्स, डाइहेड्रल्स) को MARTINI 3 के यूनिफाइड रिप्रेजेंटेशन के अनुकूल बनाने के लिए अवशेषों (residues) के बीच औसत किया गया।
  - क्लस्टर्ड दृष्टिकोण (Clustered Approach): BB–BB $^+$ –SC1 $^+$ एंगल ने cis/trans बाधाओं के कारण महत्वपूर्ण अवशेष निर्भरता दिखाई। लेखकों ने 18 अवशेषों को 5 विशिष्ट क्लस्टर्स में समूहबद्ध करने के लिए k-means clustering का उपयोग किया, और प्रत्येक के लिए विशिष्ट पैरामीटर निर्धारित किए।
  - विशेष मामले: बैकबोन से सीधे जुड़े आवेशित साइड चेन्स वाले अवशेषों (जैसे, NAsp, Nae) को संरचनात्मक पतन (conformational collapse) को रोकने के लिए विशिष्ट इम्प्रापर डाइहेड्रल रिफाइनमेंट की आवश्यकता थी।
- नॉनबॉन्डेड इंटरेक्शन: मुख्य रूप से मानक MARTINI 3 लाइब्रेरी से लिए गए, जो अत्यधिक एकत्रीकरण (aggregation) से बचने के लिए पुन: संतुलित इंटरेक्शन मैट्रिक्स का लाभ उठाते हैं।
एकीकरण (Integration): वर्कफ़्लो को martinize2 टूल में एकीकृत किया गया, जिससे स्वचालित टोपोलॉजी जनरेशन संभव हुआ।

3. प्रमुख योगदान (Key Contributions)

प्रथम MARTINI 3 पेप्टॉइड फोर्स फील्ड: विकसित किया गया एक ट्रांसफरेबल CG मॉडल जो 19 सामान्य रूप से उपयोग किए जाने वाले अवशेष प्रकारों को कवर करता है और आधुनिक MARTINI 3 फ्रेमवर्क के साथ पूरी तरह से संगत है।
दुर्लभ घटनाओं की मजबूत सैंपलिंग: धीमी cis/trans $\omega$ -डाइहेड्रल ट्रांजिशन, जो पॉलीपेप्टॉइड लचीलेपन की एक परिभाषित विशेषता है, के लिए कन्वर्ज्ड AA संदर्भ डेटा उत्पन्न करने हेतु PBMetaD का उपयोग किया गया।
एकीकृत फिर भी लचीला पैरामीट्राइजेशन: अधिकांश इंटरेक्शन के लिए एक एकीकृत बैकबोन पैरामीटर सेट बनाया, जबकि अवशेष-विशिष्ट एंगल विविधताओं के लिए क्लस्टरिंग रणनीति का उपयोग किया, जिससे रासायनिक सटीकता के साथ ट्रांसफरेबिलिटी का संतुलन बना रहा।
स्वचालित वर्कफ़्लो: मॉडल को martinize2 टूल में एकीकृत किया गया और GitHub के माध्यम से रिलीज़ किया गया, जिससे शोध समुदाय के लिए आसानी से MARTINI 3 पेप्टॉइड संरचनाएं और टोपोलॉजी बनाना संभव हो गया।
प्रदर्शन लाभ: ऑल-एटम सिमुलेशन की तुलना में ~57-गुना कम्प्यूटेशनल दक्षता प्राप्त की गई, जबकि संरचनात्मक और ऊष्मप्रवैगिकी (thermodynamic) निष्ठा बनाए रखी गई।

4. परिणाम और सत्यापन (Results and Validation)

मॉडल को कई अवलोकनों (observables) पर AA संदर्भ डेटा के विरुद्ध कड़ाई से सत्यापित किया गया:

पॉलिमर स्केलिंग (रेडियस ऑफ गाइरेशन, $R_g$ ):
- CG मॉडल ने लघु श्रृंखलाओं ( $N \le 10$ ) के लिए AA $R_g$ रुझानों को अधिकांश अवशेषों के लिए <30% विचलन के साथ पुनरुत्पादित किया।
- अपवाद: NAsp (एसिडिक) ने एकल बीड पर चार्ज लोकलाइजेशन के कारण CG में महत्वपूर्ण ओवर-कॉम्पैक्शन दिखाया। इसे एक विशिष्ट अवशेष-विशिष्ट इम्प्रापर डाइहेड्रल (NAsp*) पेश करके सुधारा गया, जिससे लघु श्रृंखलाओं के लिए बेहतर सहमति प्राप्त हुई।
बैकबोन कन्फर्मेशनल लैंडस्केप:
- बैकबोन डाइहेड्रल्स ( $\psi_1, \psi_2$ ) के फ्री एनर्जी सर्फेस (FES), CG सिमुलेशन से प्राप्त, PBMetaD-एन्हांस्ड AA सिमुलेशन से प्राप्त व्यापक, कन्वर्ज्ड परिदृश्यों से मेल खाते हैं।
- मॉडल ने सुलभ कन्फर्मेशनल स्पेस को सफलतापूर्वक कैप्चर किया, जिसमें अनबायस्ड AA सिमुलेशन में प्रतिबंधित क्षेत्र भी शामिल हैं।
डाइमेराइजेशन थर्मोडायनामिक्स:
- हाइड्रोफोबिक (NVal), पोलर (NCys), केशनिक (NArg), और एरोमैटिक (NPhe) टेट्रामर्स के डाइमेराइजेशन के लिए पोटेंशियल ऑफ मीन फोर्स (PMF) ने AA डेटा के साथ उत्कृष्ट सहमति दिखाई।
- डाइमेराइजेशन फ्री एनर्जी विभिन्न प्रणालियों में 1.69 $k_BT$ से कम भिन्न थी, जो अवशेष-विशिष्ट जुड़ाव व्यवहार के सटीक कैप्चर को दर्शाती है।
अनुक्रम-निर्भर एकत्रीकरण (Sequence-Dependent Aggregation):
- 50-चेन हेटरोट्राइमर सिस्टम (FXF, FKF, KFF, XFF) के सिमुलेशन ने प्रायोगिक रुझानों को सफलतापूर्वक पुनरुत्पादित किया: FXF (फेनिल-लाइसिन-फेनिल एनालॉग) ने रैखिक नैनोवायर बनाए, जबकि अन्य गोलाकार संरचनाएं बनाईं।
- हालांकि CG मॉडल ने क्लस्टर आकार को थोड़ा अधिक आंका (एक ज्ञात MARTINI प्रवृत्ति), इसने अनुक्रम-निर्भर एकत्रीकरण पदानुक्रम और ओरिएंटेशनल ऑर्डरिंग को सही ढंग से कैप्चर किया।

5. महत्व (Significance)

यह कार्य पॉलीपेप्टॉइड्स के मेसोस्केल सिमुलेशन के लिए एक ट्रांसफरेबल और कम्प्यूटेशनल रूप से कुशल फ्रेमवर्क स्थापित करता है। परमाणु विवरण और बड़े पैमाने की गतिशीलता के बीच की खाई को पाटकर, नया फोर्स फील्ड निम्नलिखित को सक्षम बनाता है:

तर्कसंगत डिजाइन (Rational Design): अगली पीढ़ी के अनुक्रम-विशिष्ट कार्यात्मक सामग्री, दवा वितरण प्रणाली और बायोमिमेटिक नैनोस्ट्रक्चर के डिजाइन को सुगम बनाना।
स्केलेबिलिटी (Scalability): शोधकर्ताओं को बड़े पैमाने के पेप्टॉइड असेंबली, झिल्ली अंतःक्रिया और स्व-संयोजन प्रक्रियाओं को सिम्युलेट करने की अनुमति देना, जो पहले कम्प्यूटेशनल रूप से निषेधात्मक थे।
इकोसिस्टम एकीकरण (Ecosystem Integration): सिंथेटिक फोल्डामर्स और अनुक्रम-परिभाषित पॉलिमर को शामिल करने के लिए MARTINI 3 इकोसिस्टम की बहुमुखी प्रतिभा का विस्तार करना, जो गैर-प्राकृतिक बायो-पॉलिमर के मल्टीस्केल मॉडलिंग के भविष्य के विकास का मार्ग प्रशस्त करता है।