संरचना और संरचना में झिल्ली लिपिड विविधता सेलुलर प्रक्रियाओं में एक महत्वपूर्ण योगदानकर्ता है और बीमारी का एक मार्कर हो सकता है। आणविक गतिशीलता सिमुलेशन हमें परमाणु संकल्प पर बायोमोलेक्यूल्स के साथ झिल्ली और उनकी बातचीत का अध्ययन करने की अनुमति देते हैं। यहां, हम जटिल झिल्ली प्रणालियों के निर्माण, चलाने और विश्लेषण करने के लिए एक प्रोटोकॉल प्रदान करते हैं।
लिपिड कोशिका झिल्ली के संरचनात्मक निर्माण खंड हैं; लिपिड प्रजातियां सेल ऑर्गेनेल और जीवों में भिन्न होती हैं। इस किस्म के परिणामस्वरूप झिल्ली में विभिन्न यांत्रिक और संरचनात्मक गुण होते हैं जो सीधे इस इंटरफ़ेस पर होने वाले अणुओं और प्रक्रियाओं को प्रभावित करते हैं। लिपिड संरचना गतिशील है और सेल सिग्नलिंग प्रक्रियाओं को संशोधित करने के लिए काम कर सकती है। बायोमोलेक्यूल्स के बीच बातचीत की भविष्यवाणी करने और प्रयोगात्मक अवलोकनों को आणविक अंतर्दृष्टि प्रदान करने के लिए कम्प्यूटेशनल दृष्टिकोण का तेजी से उपयोग किया जाता है। आणविक गतिशीलता (एमडी) सांख्यिकीय यांत्रिकी पर आधारित एक तकनीक है जो उन पर कार्य करने वाले बलों के आधार पर परमाणुओं के आंदोलन की भविष्यवाणी करती है। एमडी सिमुलेशन का उपयोग बायोमोलेक्यूल्स की बातचीत को चिह्नित करने के लिए किया जा सकता है। यहां, हम संक्षेप में तकनीक का परिचय देते हैं, शुरुआती लोगों के लिए व्यावहारिक कदमों की रूपरेखा तैयार करते हैं जो लिपिड बाइलेयर का अनुकरण करने में रुचि रखते हैं, शुरुआती अनुकूल सॉफ्टवेयर के साथ प्रोटोकॉल का प्रदर्शन करते हैं, और विकल्पों, चुनौतियों और प्रक्रिया के महत्वपूर्ण विचारों पर चर्चा करते हैं। विशेष रूप से, हम सिमुलेशन में उपयुक्त हाइड्रोफोबिक और यांत्रिक वातावरण को पकड़ने के लिए रुचि के सेल झिल्ली को मॉडल करने के लिए जटिल लिपिड मिश्रण का उपयोग करने की प्रासंगिकता पर जोर देते हैं। हम कुछ उदाहरणों पर भी चर्चा करते हैं जहां झिल्ली संरचना और गुण अन्य बायोमोलेक्यूल्स के साथ बाइलेयर की बातचीत को संशोधित करते हैं।
लिपिड झिल्ली के प्रमुख घटक हैं, जो कोशिकाओं के लिए सीमाएं प्रदान करते हैं और इंट्रासेल्युलर कंपार्टमेंटलाइजेशन 1,2,3 को सक्षम करते हैं। लिपिड एम्फीफिलिक होते हैं, जिसमें एक ध्रुवीय सिर समूह और दो हाइड्रोफोबिक फैटी एसिड पूंछ होते हैं; ये पानी के साथ हाइड्रोफोबिक श्रृंखलाओं के संपर्क को कम करने के लिए एक बाईलेयर में स्वयं इकट्ठा होते हैं 3,4. हाइड्रोफिलिक हेड ग्रुप और हाइड्रोफोबिक पूंछ के विभिन्न संयोजनों के परिणामस्वरूप जैविक झिल्ली में लिपिड के विभिन्न वर्ग होते हैं, जैसे ग्लिसरोफॉस्फोलिपिड्स, स्फिंगोलिपिड्स और स्टेरोल्स (चित्रा 1)1,5,6। ग्लिसरोफॉस्फोलिपिड्स यूकेरियोटिक कोशिका झिल्ली के प्राथमिक निर्माण खंड हैं जो ग्लिसरोफॉस्फेट, लंबी श्रृंखला फैटी एसिड और कम आणविक भारके सिर समूहों से बने होते हैं। लिपिड नामकरण सिर समूहों में अंतर पर आधारित है; उदाहरणों में फॉस्फेटिडिल-कोलीन (पीसी), फॉस्फेटिडिल-इथेनॉलमाइन (पीई), फॉस्फेटिडिल-सेरीन (पीएस), फॉस्फेटिडिल-ग्लिसरॉल (पीजी), फॉस्फेटिडिल-इनोसिटोल (पीआई), या असंशोधित फॉस्फेटिडिक एसिड (पीए) 5,6 शामिल हैं। हाइड्रोफोबिक पूंछ के लिए, रीढ़ की संरचना के साथ संतृप्ति की लंबाई और डिग्री भिन्न होती है। संभावित संयोजन कई हैं, जिसके परिणामस्वरूप स्तनधारी कोशिकाओं में हजारों लिपिड प्रजातियांहैं। झिल्ली लिपिड संरचना में परिवर्तन विभिन्न यांत्रिक और संरचनात्मक झिल्ली गुणों को जन्म देता है जो अभिन्न झिल्ली प्रोटीन और परिधीय प्रोटीन 2,6 दोनों की गतिविधि को प्रभावित करते हैं।
चित्र 1. प्रतिनिधि लिपिड संरचनाएं। फैटी एसिड पूंछ नीले बक्से में दिखाए जाते हैं, नारंगी में सामान्य लिपिड हेड समूह, और बैंगनी में नमूना बैकबोन। कृपया इस आंकड़े का एक बड़ा संस्करण देखने के लिए यहाँ क्लिक करें.
लिपिड सेलुलर प्रक्रियाओं में सक्रिय खिलाड़ी हैं, सिग्नलिंग कैस्केड में प्रोटीन सक्रियण, और स्वस्थ सेल होमियोस्टैसिस 8,9। परिवर्तित लिपिड गतिशीलता संक्रमण का परिणाम है या रोग10,11,12,13,14,15 के रोगजनन के मार्कर हो सकते हैं। कोशिका के लिए बाधाओं के रूप में, झिल्ली लिपिड का अध्ययन और छोटे अणुओं के पारगम्यता में उनकी भूमिका दवा वितरण प्रणाली और झिल्ली विघटन तंत्र16,17 के लिए प्रासंगिक है। रासायनिक विविधता और ऑर्गेनेल, ऊतकों और जीवों में लिपिड प्रजातियों के विभिन्न अनुपात जटिलझिल्ली गतिशीलता को जन्म देते हैं। इसलिए लिपिड बाइलेयर के मॉडलिंग अध्ययन में इन विशेषताओं को बनाए रखना महत्वपूर्ण है, खासकर जब एक अध्ययन का लक्ष्य झिल्ली के साथ अन्य बायोमोलेक्यूल्स की बातचीत की जांच करना है। एक मॉडल में विचार करने के लिए लिपिड प्रजातियां जीव और रुचि के सेलुलर डिब्बे पर निर्भर करती हैं। उदाहरण के लिए, पीजी लिपिड प्रकाश संश्लेषक बटेरिया18 में इलेक्ट्रॉन हस्तांतरण के लिए महत्वपूर्ण हैं, जबकि फॉस्फोराइलेटेड इनोसिटोल लिपिड (पीआईपी) प्लाज्मा झिल्ली (पीएम) गतिशीलता और स्तनधारी कोशिकाओं में सिग्नलिंग कैस्केड 19,20 में प्रमुख खिलाड़ी हैं। कोशिका के अंदर, पीएम, एंडोप्लाज्मिक रेटिकुलम (ईआर), गोल्गी और माइटोकॉन्ड्रियल झिल्ली में अद्वितीय लिपिड बहुतायत होते हैं जो उनके कार्य को प्रभावित करते हैं। उदाहरण के लिए, ईआर लिपिड बायोजेनेसिस के लिए केंद्र है और कोलेस्ट्रॉल को पीएम और गोल्गी तक पहुंचाता है; इसमें पीसी और पीई की प्रचुरता के साथ एक उच्च लिपिड विविधता होती है, लेकिन कम स्टेरोल सामग्री, जो झिल्ली तरलता21,22,23,24 को बढ़ावा देती है। इसके विपरीत, पीएम जीव25 के आधार पर सैकड़ों और यहां तक कि हजारों लिपिड प्रजातियों को शामिल करता है, इसमें स्फिंगोलिपिड्स और कोलेस्ट्रॉल के उच्च स्तर होते हैं जो इसे सेल24 में अन्य झिल्ली की तुलना में एक विशिष्ट कठोरता देते हैं। पत्रक विषमता को पीएम जैसी झिल्लियों के लिए माना जाना चाहिए, जिसमें स्फिंगोमाइलिन, पीसी और कोलेस्ट्रॉल से भरपूर एक बाहरी पत्रक होता है, और पीई, पीआई और पीएस में समृद्ध एक आंतरिक पत्रक होता है जो सिग्नलिंग कैस्केड24 के लिए महत्वपूर्ण होता है। अंत में, लिपिड विविधता भी माइक्रो-डोमेन के गठन को प्रेरित करती है जो पैकिंग और आंतरिक क्रम में भिन्न होती है, जिसे लिपिड राफ्ट24,26 के रूप में जाना जाता है; ये पार्श्व विषमता प्रदर्शित करते हैं, सेलुलर सिग्नलिंग26 में महत्वपूर्ण भूमिका निभाने के लिए परिकल्पित हैं, और उनकी क्षणिक प्रकृति के कारण अध्ययन करना कठिन है।
फ्लोरोस्कोपी, स्पेक्ट्रोस्कोपी, और मॉडल झिल्ली प्रणालियों जैसे विशाल यूनिलामेलर पुटिकाओं (जीयूवी) जैसी प्रयोगात्मक तकनीकों का उपयोग झिल्ली के साथ बायोमोलेक्यूल्स की बातचीत की जांच करने के लिए किया गया है। हालांकि, इसमें शामिल घटकों की जटिल और गतिशील प्रकृति को अकेले प्रयोगात्मक तरीकों से पकड़ना मुश्किल है। उदाहरण के लिए, प्रोटीन के ट्रांसमेम्ब्रेन डोमेन की इमेजिंग, ऐसे अध्ययनों में उपयोग की जाने वाली झिल्ली की जटिलता और ब्याज की प्रक्रिया के दौरान मध्यवर्ती या क्षणिक अवस्थाओं की पहचान27,28,29 पर सीमाएं हैं। 1980 के दशक में लिपिड मोनोलेयर औरबाइलेयर के आणविक सिमुलेशन के आगमन के बाद से, लिपिड-प्रोटीन सिस्टम और उनकी बातचीत को अब आणविक स्तर पर निर्धारित किया जा सकता है। आणविक गतिशीलता (एमडी) सिमुलेशन एक सामान्य कम्प्यूटेशनल तकनीक है जो कणों के आंदोलन की भविष्यवाणी उनके इंटरमॉलिक्युलर बलों के आधार पर करती है। एक योजक इंटरैक्शन क्षमता सिस्टम30 के कणों के बीच बंधुआ और गैर-बंधुआ बातचीत का वर्णन करती है। इन इंटरैक्शन को मॉडल करने के लिए उपयोग किए जाने वाले मापदंडों के सेट को सिमुलेशन फोर्सफील्ड (एफएफ) कहा जाता है। इन मापदंडों को प्रारंभिक गणनाओं, अर्ध-अनुभवजन्य और क्वांटम यांत्रिक गणनाओं से प्राप्त किया जाता है, और एक्स-रे और इलेक्ट्रॉन विवर्तन प्रयोगों, एनएमआर, इन्फ्रारेड, रमन और न्यूट्रॉन स्पेक्ट्रोस्कोपी से डेटा को पुन: पेश करने के लिए अनुकूलित किया जाताहै।
एमडी सिमुलेशन का उपयोग संकल्प32,33,34 के विभिन्न स्तरों पर सिस्टम का अध्ययन करने के लिए किया जा सकता है। सिस्टम जो विशिष्ट बायोमोलेक्यूलर इंटरैक्शन, हाइड्रोजन बॉन्ड और अन्य उच्च-रिज़ॉल्यूशन विवरणों को चिह्नित करने का लक्ष्य रखते हैं, उनका अध्ययन ऑल-एटम (एए) सिमुलेशन के साथ किया जाता है। इसके विपरीत, मोटे-दाने वाले (सीजी) सिमुलेशन कम्प्यूटेशनल लागत को कम करने और बड़े पैमाने पर गतिशीलता33 की जांच करने के लिए बड़े कार्यात्मक समूहों में गांठ परमाणुओं को जोड़ते हैं। इन दोनों के बीच में स्थित हैं संयुक्त-परमाणु (यूए) सिमुलेशन, जहां हाइड्रोजन परमाणुओं को उनके संबंधित भारी परमाणुओं के साथ जोड़ा जाता है ताकि गणना33,35 में तेजी लाई जा सके। एमडी सिमुलेशन लिपिड झिल्ली की गतिशीलता और अन्य अणुओं के साथ उनकी बातचीत की खोज के लिए एक शक्तिशाली उपकरण हैं और झिल्ली इंटरफ़ेस पर रुचि की प्रक्रियाओं के लिए आणविक स्तर के तंत्र प्रदान करने के लिए काम कर सकते हैं। इसके अतिरिक्त, एमडी सिमुलेशन प्रयोगात्मक लक्ष्यों को कम करने और सूक्ष्म इंटरैक्शन के आधार पर किसी दिए गए सिस्टम के मैक्रोमोलेक्यूलर गुणों की भविष्यवाणी करने के लिए काम कर सकते हैं।
संक्षेप में, प्रारंभिक निर्देशांक, वेग और निरंतर तापमान और दबाव जैसी स्थितियों के एक सेट को देखते हुए, प्रत्येक कण की स्थिति और वेग की गणना इंटरैक्शन पोटेंशियल और न्यूटन के गति के नियम के संख्यात्मक एकीकरण के माध्यम से की जाती है। यह पुनरावृत्ति रूप से दोहराया जाता है, जिससे सिमुलेशन प्रक्षेपवक्र30 उत्पन्न होता है। ये गणना एमडी इंजन के साथ की जाती है; कई ओपन-सोर्स पैकेजों के बीच, GROMACS36 सबसे अधिक इस्तेमाल किए जाने वाले इंजनों में से एक है और जिसे हम यहां वर्णित करते हैं। इसमें विश्लेषण और सिस्टम के प्रारंभिक निर्देशांक के निर्माण के लिए उपकरण भी शामिल हैं जिन्हें सिम्युलेटेडकिया जा सकता है। अन्य एमडी इंजनों में एनएएमडी38 शामिल हैं; चार्मम39, और एम्बर40, जिसे उपयोगकर्ता किसी दिए गए सिस्टम के कम्प्यूटेशनल प्रदर्शन के आधार पर अपने विवेक पर चुन सकता है। सिमुलेशन के दौरान प्रक्षेपपथ की कल्पना करना और साथ ही परिणामों के विश्लेषण और व्याख्या के लिए महत्वपूर्ण है। विभिन्न प्रकार के उपकरण उपलब्ध हैं; यहां हम दृश्य आणविक गतिशीलता (वीएमडी) पर चर्चा करते हैं जो सुविधाओं की एक विस्तृत श्रृंखला प्रदान करता है, जिसमें विस्तृत ड्राइंग और रंग विधियों के साथ तीन आयामी (3-डी) विज़ुअलाइज़ेशन, वॉल्यूमेट्रिक डेटा विज़ुअलाइज़ेशन, एमडी सिमुलेशन सिस्टम के निर्माण, तैयारी और विश्लेषण और सिस्टम आकार पर कोई सीमा नहीं होने के साथ प्रक्षेपवक्र-फिल्म बनाना शामिल है, अगर मेमोरी41,42,43 उपलब्ध है।
सिस्टम घटकों के बीच अनुमानित गतिशीलता की सटीकता सीधे प्रक्षेपवक्र के प्रसार के लिए चुने गए एफएफ से प्रभावित होती है। अनुभवजन्य एफएफ परमेट्राइजेशन प्रयासों को कुछ शोध समूहों द्वारा आगे बढ़ाया जाता है। एमडी के लिए सबसे स्थापित और आम एफएफ में चार्मम39, एम्बर 40, मार्टिनी 44, ओपीएलएस45 और एसआईआरएएच 46 शामिल हैं। ऑल-एटम एडिटिव CHARMM36 (सी 36) बल क्षेत्र47 का व्यापक रूप से झिल्ली प्रणालियों के एए एमडी के लिए उपयोग किया जाता है क्योंकि यह प्रयोगात्मक संरचनात्मक डेटा को सटीक रूप से पुन: पेश करता है। यह मूल रूप से चार्मम समुदाय द्वारा विकसित किया गया था, और यह ग्रोमाक्स और एनएएमडी जैसे कई एमडी इंजनों के साथ संगत है। सामान्य एफएफ में सुधार के बावजूद, पैरामीटर सेटों में सुधार करने के लिए निरंतर प्रयास किया जाता है ताकि उन भविष्यवाणियों की अनुमति मिल सके जो प्रयोगात्मक अवलोकनों को बारीकी से पुन: पेश करते हैं, जो अध्ययनकी विशेष प्रणालियों में रुचियों से प्रेरित हैं।
लिपिड झिल्ली का अनुकरण करते समय एक चुनौती सिमुलेशन प्रक्षेपवक्र की लंबाई निर्धारित कर रही है। यह काफी हद तक विश्लेषण किए जाने वाले मैट्रिक्स और उस प्रक्रिया पर निर्भर करता है जिसे कोई चिह्नित करना चाहता है। आमतौर पर, जटिल लिपिड मिश्रण को संतुलन तक पहुंचने के लिए अधिक समय की आवश्यकता होती है, क्योंकि अधिक प्रजातियों को झिल्ली तल पर फैलने और एक स्थिर पार्श्व संगठन तक पहुंचने के लिए पर्याप्त समय होना चाहिए। एक सिमुलेशन को संतुलन में कहा जाता है जब ब्याज की संपत्ति एक पठार पर पहुंच जाती है और एक स्थिर मूल्य के बारे में उतार-चढ़ाव करती है। रुचि के गुणों और अंतःक्रियाओं पर उचित सांख्यिकीय विश्लेषण करने के लिए कम से कम 100-200 एनएस समतुल्य प्रक्षेपवक्र प्राप्त करना आम बात है। लिपिड मिश्रण और शोध प्रश्न की जटिलता के आधार पर 200-500 एनएस के बीच झिल्ली-केवल सिमुलेशन चलाना आम है। प्रोटीन-लिपिड इंटरैक्शन को आमतौर पर 500-2000 एनएस के बीच लंबे सिमुलेशन समय की आवश्यकता होती है। झिल्ली प्रणालियों के साथ नमूनाकरण और अवलोकन योग्य गतिशीलता में तेजी लाने के लिए कुछ दृष्टिकोण हैं: (i) अत्यधिक मोबाइल झिल्ली मिमेटिक (एचएमएमएम) मॉडल, जो नमूना करण में तेजी लाने के लिए कार्बनिक विलायक के साथ झिल्लीमें लिपिड के अंत कार्बन को प्रतिस्थापित करता है; और (ii) हाइड्रोजन द्रव्यमान पुनर्विभाजन (एचएमआर), जो हाइड्रोजन परमाणुओं के साथ एक प्रणाली के भीतर भारी परमाणुओं के द्रव्यमान के एक अंश को जोड़ता है ताकि एक बड़े सिमुलेशन टाइमस्टेप51 के उपयोग की अनुमति मिल सके।
निम्नलिखित प्रोटोकॉल एए एमडी का उपयोग करके यथार्थवादी झिल्ली मॉडल बनाने, चलाने और विश्लेषण करने के लिए एक शुरुआती-अनुकूल दृष्टिकोण पर चर्चा करता है। एमडी सिमुलेशन की प्रकृति को देखते हुए, परिणामों की प्रजनन क्षमता और उचित सांख्यिकीय विश्लेषण के लिए कई प्रक्षेपपथ चलाए जाने चाहिए। यह वर्तमान प्रथा है कि रुचि की प्रति प्रणाली में कम से कम तीन प्रतिकृतियां चलाई जाती हैं। एक बार जब जीव और रुचि की प्रक्रिया के लिए लिपिड प्रजातियों का चयन किया जाता है, तो झिल्ली-केवल प्रणाली के सिमुलेशन प्रक्षेपवक्र के निर्माण, चलाने और विश्लेषण करने के लिए बुनियादी चरणों को चित्र 2 में रेखांकित और सारांशित किया जाता है।
चित्र 2. एमडी सिमुलेशन चलाने के लिए योजनाबद्ध। नारंगी बक्से प्रोटोकॉल में वर्णित तीन मुख्य चरणों के अनुरूप हैं। नीचे सिमुलेशन प्रक्रिया का वर्कफ़्लो है। सिस्टम सेटअप के दौरान, सॉल्वेटेड मेम्ब्रेन सिस्टम के प्रारंभिक निर्देशांक वाले सिस्टम को चार्म-जीयूआई मेम्ब्रेन बिल्डर जैसे सिस्टम इनपुट जनरेटर के साथ बनाया जाता है। इनपुट फ़ाइलों को उच्च-प्रदर्शन कंप्यूटिंग क्लस्टर में स्थानांतरित करने के बाद, सिमुलेशन प्रक्षेपवक्र को एमडी इंजन का उपयोग करके प्रचारित किया जाता है, जैसे कि GROMACS। प्रक्षेपवक्र विश्लेषण कंप्यूटर क्लस्टर या विज़ुअलाइज़ेशन के साथ एक स्थानीय कार्य स्टेशन पर किया जा सकता है। विश्लेषण तब किया जाता है, या तो अंतर्निहित विश्लेषण कोड जैसे GROMACS और VMD के साथ पैकेज का उपयोग करके, या बैश स्क्रिप्ट या विभिन्न पायथन पुस्तकालयों का उपयोग करके। कृपया इस आंकड़े का एक बड़ा संस्करण देखने के लिए यहाँ क्लिक करें.
प्रायोगिक तकनीक क्रायो-इलेक्ट्रॉन माइक्रोस्कोपी (क्रायो-ईएम) 58, फ्लोरेसेंस तकनीक और परमाणु बल माइक्रोस्कोपी (एएफएम) 59 का उपयोग करके उच्च रिज़ॉल्यूशन पर बायोमोलेक्यूल्स की कल्पना कर…
The authors have nothing to disclose.
लेखक ों ने इस पांडुलिपि के लेखन के दौरान अपने सिमुलेशन प्रक्षेपवक्र और चर्चाओं के लिए जिनहुई ली और रिकार्डो एक्स रामिरेज़ को धन्यवाद दिया। ओ.सी. को बफ़ेलो प्रेसिडेंशियल फैलोशिप विश्वविद्यालय और नेशनल इंस्टीट्यूट ऑफ हेल्थ की पहल फॉर मैक्सिममाइज़िंग स्टूडेंट डेवलपमेंट ट्रेनिंग ग्रांट 1T32GM144920-01 द्वारा मार्गरीटा एल डुबोकोविच (पीआई) को सम्मानित किया गया था।
Anaconda3 | Anaconda Inc (Python & related libraries) | N/A | |
CHARMM-GUI.org | Im lab, Lehigh University | N/A | |
GROMACS | GROMACS development team | N/A | |
Linux HPC Cluster | UB CCR | N/A | |
MATLAB | MathWorks | N/A | |
VMD | Theoretical and Computational Biophysics Group | N/A |