Solvation Free Energy Descriptors from Density Functional Theory for Unified Retention Prediction of Flavonoids and Isoflavonoids

कल्पना कीजिए कि आप पहेली के टुकड़ों के एक विशाल ढेर को छाँटने की कोशिश कर रहे हैं। कुछ टुकड़े लगभग एक जैसे दिखते हैं, लेकिन वे अलग-अलग चित्रों के हैं। रसायन विज्ञान की दुनिया में, वैज्ञानिक अक्सर इस समस्या का सामना करते हैं जिन्हें फ्लेवोनोइड्स (flavonoids) कहा जाता है। ये पौधों में पाए जाने वाले प्राकृतिक यौगिक हैं जो कैंसर से लड़ सकते हैं या एंटीऑक्सीडेंट के रूप में कार्य कर सकते हैं। इनका अध्ययन करने के लिए, वैज्ञानिक हाई-परफॉर्मेंस लिक्विड क्रोमैटोग्राफ (HPLC) नामक मशीन का उपयोग करते हैं। HPLC को एक लंबी, घुमावदार फिसलने वाली स्लाइड (slide) के रूप में सोचें। आप मिश्रण के अणुओं को ऊपर की ओर डालते हैं, और वे नीचे फिसलते हैं। कुछ तेजी से फिसलते हैं, और कुछ फंस जाते हैं और धीरे-धीरे फिसलते हैं। किसी अणु को नीचे पहुँचने में लगने वाला समय उसका "रिटेंशन टाइम" (retention time) है।

आमतौर पर, यह जानने के लिए कि कौन सा अणु कौन सा है, वैज्ञानिकों को एक संदर्भ मानक (reference standard) की आवश्यकता होती है—एक आदर्श, ज्ञात नमूना जिसकी तुलना की जा सके। लेकिन कई दुर्लभ पौधों के यौगिकों के लिए, ऐसे संदर्भ नमूने अभी तक मौजूद नहीं हैं। इसलिए, वैज्ञानिक यह अनुमान लगाने की कोशिश करते हैं कि एक अणु HPLC स्लाइड पर कितनी तेजी से फिसलेगा, केवल उसके आकार और विद्युत गुणों को देखकर। इसे क्वांटिटेटिव स्ट्रक्चर-रिटेंशन रिलेशनशिप (QSRR) कहा जाता है। बड़ी चुनौती यह रही है कि अणुओं के थोड़े अलग आकार (जैसे कि एक अलग स्थान पर रिंग जुड़ी होना) के लिए आमतौर पर पूरी तरह से अलग भविष्यवाणी नियमों की आवश्यकता होती है। यह ऐसा है जैसे कारों के लिए एक नियम पुस्तिका हो और मोटरसाइकिलों के लिए बिल्कुल अलग, भले ही दोनों एक ही सड़क पर चलती हों। यह शोध पत्र पूछता है: क्या हम एक एकल, सार्वभौमिक नियम पा सकते हैं जो दोनों के लिए काम करे?

---

महान आणविक स्लाइड: दो आकारों के लिए एक नियम पुस्तिका

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

टीम ने पाया कि उन्हें अणु के हर छोटे मोड़ और घुमाव को ट्रैक करने की आवश्यकता नहीं है। इसके बजाय, उन्होंने पाया कि कुल विलायक मुक्त ऊर्जा (total solvation free energy) जादुई कुंजी है। इसे समझने के लिए, कल्पना करें कि अणु एक यात्री है जो दो अलग-अलग भीड़ों के माध्यम से गुजरने की कोशिश कर रहा है: पानी के अणुओं की भीड़ और मेथनॉल के अणुओं की भीड़। "विलायक मुक्त ऊर्जा" मूल रूप से इस बात का स्कोर है कि यात्री उस भीड़ में कितना फिट बैठता है या उसे पसंद करता है या नहीं। यदि यात्री अच्छी तरह से फिट बैठता है, तो स्कोर कम होता है; यदि वह अजीब और दबाव डालने वाला महसूस करता है, तो स्कोर अधिक होता है। शोधकर्ताओं ने गणना की कि इन अणुओं को पानी बनाम मेथनॉल में रहने के लिए कितनी ऊर्जा की आवश्यकता होती है।

दो-चरणीय नृत्य: सही मुद्रा खोजना

इन ऊर्जा स्कोरों की गणना करना कठिन है क्योंकि अणु लचीले होते हैं; वे हिलते और मुड़ते हैं। यदि आप कंप्यूटर को अणु का आकार खोजने के लिए कहते हैं, तो वह एक "लोकल मिनिमम" (local minimum) में फंस सकता है—ऊर्जा परिदृश्य में एक आरामदायक छोटा गड्ढा जो वास्तव में सबसे निचला बिंदु नहीं है। यह एक हाइकर (पर्वतारोही) की तरह है जो एक आरामदायक गुफा पाकर रुक जाता है, यह महसूस किए बिना कि अगली पहाड़ी के ठीक नीचे एक गहरी घाटी है।

इसे हल करने के लिए, लेखकों ने एक चतुर दो-स्तरीय रणनीति का उपयोग किया:

1. एक कच्चा स्केच (PM6): सबसे पहले, उन्होंने अणु की "बी-रिंग" को एक त्वरित स्केच कलाकार की तरह स्कैन करने के लिए एक तेज़, अर्ध-अनुभवजन्य विधि (PM6) का उपयोग किया। उन्होंने रिंग को घुमाया ताकि सामान्य ऊर्जा परिदृश्य और सबसे गहरे गड्ढे (ग्लोबल मिनिमम) को देखा जा सके। यह सस्ता और तेज़ था, जिससे यह सुनिश्चित हुआ कि वे गलत गुफा में न फंस जाएं।
2. एक हाई-डेफिनिशन फोटो (DFT): एक बार जब वे सही घाटी पा लेते, तो वे डेंसिटी फंक्शनल थ्योरी (DFT) नामक एक अत्यंत सटीक विधि पर स्विच कर देते थे, जो B3LYP-D3/6-31G(d) स्तर पर थी। इसने आकार को परिष्कृत किया, जिससे उन्हें अणु की सबसे स्थिर मुद्रा की एक हाई-डेफिनिशन फोटो मिली।

उन्होंने यह भी देखा कि अणु कहाँ हाथ मिला सकता है (हाइड्रोजन बॉन्ड बना सकता है) पानी और मेथनॉल के साथ। मॉलिक्यूलर इलेक्ट्रोस्टैटिक पोटेंशियल (ESP) विश्लेषण नामक एक तकनीक का उपयोग करके, उन्होंने उन सटीक स्थानों की पहचान की जहाँ अणु के किसी विलायक अणु को पकड़ने की सबसे अधिक संभावना है। फिर उन्होंने इन विशिष्ट "हाथ मिलाने" (handshakes) की ऊर्जा की गणना की और उन्हें कुल स्कोर में जोड़ दिया।

जादुई सूत्र

इन दो ऊर्जा स्कोरों के साथ—पानी के लिए (EW) और मेथनॉल के लिए (EM)—टीम ने यह अनुमान लगाने के लिए कि अणु स्लाइड पर कितनी देर तक रहेगा, एक सरल गणितीय समीकरण बनाया:

tR = −0.8486 EW + 0.9396 EM + 40.55

यहाँ, tR रिटेंशन टाइम (स्लाइड पर अणु के रुकने का समय) है।

• EW के लिए नकारात्मक चिह्न का अर्थ है कि यदि कोई अणु पानी को बहुत पसंद करता है (कम ऊर्जा), तो वह स्लाइड से तेज़ी से फिसल जाता है (कम रिटेंशन टाइम)।
• EM के लिए सकारात्मक चिह्न का अर्थ है कि यदि कोई अणु मेथनॉल को पसंद करता है, तो वह अधिक समय तक टिका रहता है।

परिणाम प्रभावशाली थे। मॉडल ने डेटा के उतार-चढ़ाव की 92.8% व्याख्या की (R² = 0.928)। जब उन्होंने एक समय में एक अणु को छोड़कर परीक्षण किया कि क्या मॉडल सही ढंग से अनुमान लगा सकता है, तो स्कोर अभी भी एक मजबूत 0.872 था। मॉडल बनाने के लिए उपयोग किए गए 11 अणुओं के लिए, औसत त्रुटि केवल 2.80% थी। इसका मतलब है कि मॉडल एक ही समीकरण का उपयोग करके फ्लेवोन और आइसोफ्लेवोन दोनों के स्लाइड समय की सफलतापूर्वक भविष्यवाणी कर सकता है, जिससे यह सिद्ध होता है कि यदि आप सही ऊर्जा स्कोर देखते हैं, तो बी-रिंग की स्थिति नियमों को तोड़ती नहीं है।

वह एक अणु जिसने नियमों को तोड़ा

हालाँकि, कहानी पूरी तरह से परफेक्ट नहीं है। शोधकर्ताओं ने नियोबैविसोफ्लेवोन (neobavaisoflavone) नामक एक अणु पर मॉडल का परीक्षण किया, जिसमें एक भारी, गैर-ध्रुवीय समूह (एक आइसोपेंटेनिल समूह) जुड़ा हुआ है। यह अणु एक ऐसे यात्री की तरह है जिसने वसा (grease) से बना एक विशाल, भारी बैकपैक पहना हुआ है। मानक ऊर्जा गणनाएँ (जो विद्युत आकर्षण और हाइड्रोजन बॉन्डिंग पर ध्यान केंद्रित करती हैं) इस "वसा" के इंटरैक्शन को पूरी तरह से नहीं समझा सकतीं। इसके परिणामस्वरूप, नियोबैविसोफ्लेवोन के लिए मॉडल की भविष्यवाणी 26.4% गलत थी।

इस विफलता ने वास्तव में टीम को कुछ महत्वपूर्ण सिखाया: मॉडल तब पूरी तरह से काम करता है जब अणु का व्यवहार विद्युत बलों और हाइड्रोजन बॉन्डिंग द्वारा संचालित होता है। लेकिन जब कोई अणु "चिकनाई" वाले हाइड्रोफोबिक प्रभावों और डिस्पर्शन बलों (जैसे कि वह भारी बैकपैक) से हावी होता है, तो वर्तमान विधि विफल हो जाती है। ऐसा नहीं है कि गणित गलत है; बल्कि यह है कि विशिष्ट प्रकार के इंटरैक्शन को ऊर्जा स्कोर में शामिल नहीं किया गया था।

मुख्य निष्कर्ष

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