Solvation Free Energy Descriptors from Density Functional Theory for Unified Retention Prediction of Flavonoids and Isoflavonoids
تُنشئ هذه الدراسة نموذجًا موحدًا لعلاقة البنية والخاصية الكروماتوغرافية (QSRR) قادرًا على التنبؤ بالاحتجاز الكروماتوغرافي لكل من الفلافونات والأيزوفلافونات باستخدام مجموعة دنيا مكونة من وصفي طاقة حرة للذوبان، وهما وصفيان مفسران فيزيائيًا ومستمدان من نظرية الدالة الوظيفية للكثافة، مما يتغلب بذلك على القيد التقليدي المتمثل في النماذج الخاصة بالهياكل البنائية.
البحث الأصلي مرخَّص بموجب CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). هذا شرح مولَّده بالذكاء الاصطناعي للبحث أدناه. لم يكتبه المؤلفون ولم يصادقوا عليه. وللتحقق من الدقة التقنية، يرجى الرجوع إلى البحث الأصلي. اقرأ إخلاء المسؤولية الكامل
تخيل أنك تحاول فرز كومة ضخمة من قطع الأحجية المختلطة. بعض القطع تبدو متطابقة تقريبًا، لكنها تنتمي إلى صور مختلفة. في عالم الكيمياء، غالبًا ما يواجه العلماء هذه المشكلة مع جزيئات تسمى الفلافونيدات (flavonoids). هذه مركبات طبيعية توجد في النباتات ويمكنها محاربة السرطان أو العمل كمضادات للأكسدة. لدراستها، يستخدم العلماء جهازًا يسمى الكروماتوغرافيا السائلة عالية الأداء (HPLC). فكر في جهاز الـ HPLC كأنه منزلق طويل ومتعرج. تصب مزيجًا من الجزيئات في الأعلى، وهي تنزلق للأسفل. بعضها ينزلق بسرعة، وبعضها يعلق وينزلق ببطء. الوقت الذي تستغرقه الجزيئة للوصول إلى الأسفل هو "زمن الاستبقاء" (retention time) الخاص بها.
عادةً، لكي يعرف العلماء أي جزيء هو المقصود، يحتاجون إلى معيار مرجعي — عينة مثالية معروفة للمقارنة بها. ولكن بالنسبة للعديد من المركبات النباتية النادرة، فإن هذه العينات المرجعية لا توجد بعد. لذا، يحاول العلماء التنبؤ بمدى سرعة انزلاق الجزيء على منزلق الـ HPLC بمجرد النظر إلى شكله وخصائصه الكهربائية. وهذا ما يسمى بالعلاقة الكمية بين البنية وزمن الاستبقاء (QSRR). كان التحدي الكبير هو أن الجزيئات ذات الأشكال المختلفة قليلاً (مثل وجود حلقة متصلة في مكان مختلف) تتطلب عادةً قواعد تنبؤ مختلفة تمامًا. الأمر يشبه امتلاك كتاب قواعد للسيارات وكتاب قواعد مختلف تمامًا للدراجات النارية، رغم أن كليهما يسير على نفس الطريق. هذا البحث يسأل: هل يمكننا إيجًااد قاعدة واحدة عالمية تعمل لكليهما؟
المنزلق الجزيئي العظيم: كتاب قواعد واحد لشكلين مختلفين
في هذه الدراسة، تناول الباحثون ييهان شينغ، وهايتيان ليو، وجيان يان لي، وتشانغهاي سون من جامعة جياموسي، لغز التنبؤ بكيفية سلوك الفلافونات (flavones) والأيزوفلافونات (isoflavones) على منزلق الكروماتوغرافيا. هاتان المجموعتان من الجزيئات تشبهان التوائم الذين انفصلوا عند الولادة: فهما يتشاركان في نفس الهيكل الأساسي ونفس الأجزاء، لكن "حلقة B" (وهي حلقة محددة من الذرات) متصلة في مكان مختلف قليًّا. وبسبب هذا الاختلاف الطفيف، تتغير مساحاتهما الكهربائية، مما يجعل من الصعب استخدام نموذج تنبؤ واحد لكليهما.
اكتشف الفريق أنهم ليسوا بحاجة لتتبع كل التواء وانعطاف صغير في الجزيء. بدلاً من ذلك، وجدوا أن طاقة التذوب الحرة الإجمالية (total solvation free energy) هي المفتاح السحري. لفهم ذلك، تخيل أن الجزيء هو مسافر يحاول التحرك عبر حشدين مختلفين: حشد من جزيئات الماء وحشد من جزيئات الميثانول. "طاقة التذوب الحرة" هي ببساطة درجة توضح مدى حب المسافر لهذا الحشد أو عدم رغبته فيه. إذا كان المسافر يتناسب جيدًا مع الحشد، تكون الدرجة منخفضة؛ وإذا كان يشعر بالارتباك والضغط، تكون الدرجة عالية. قام الباحثون بحساب مقدار الطاقة التي يتطلبها بقاء هذه الجزيئات في الماء مقابل الميثانول.
رقصة الخطوتين: إيجاد الوضعية المثالية
حساب درجات الطاقة هذه أمر صعب لأن الجزيئات مرنة؛ فهي تتمايل وتلتوي. إذا طلبت من حاسوب إيجاد شكل الجزيء، فقد يعلق في "نهاية صغرى محلية" (local minimum) — وهي فجوة مريحة صغيرة في مشهد الطاقة ليست هي النقطة الأدنى فعليًا. الأمر يشبه متسلقًا يجد كهفًا مريحًا ويتوقف عنده، دون أن يدرك أن هناك واديًا أعمق خلف التل التالي مباشرة.
لحل هذه المشكلة، استخدم المؤلفون استراتيجية من مستويين:
- المخطط الأولي (PM6): أولاً، استخدموا طريقة شبه تجريبية سريعة (PM6) لمسح "حلقة B" للجزيء مثل رسام سريع. قاموا بتدوير الحلقة لرؤية مشهد الطاقة العام وإيجاد أعمق وادي (النهاية الصغرى العالمية). كانت هذه الطريقة رخيصة وسريعة، مما ضمن عدم تعثرهم في الكهف الخطأ.
- الصورة عالية الدقة (DFT): بمجرد العثور على الوادي الصحيح، انتقلوا إلى طريقة دقيقة للغاية تسمى نظرية الكثافة الوظيفية (DFT) عند مستوى B3LYP-D3/6-31G(d). قامت هذه الطوة بتنقيح الشكل، مما منحهم صورة عالية الدقة للوضعية الأكثر استقرارًا للجزيء.
نظروا أيضًا عن كثب إلى الأماكن التي يمكن للجزيء فيها "الإمساك بالأيدي" (تكوين روابط هيدروجينية) مع الماء والميثانول. باستخدام تقنية تسمى تحليل الجهد الكهروستاتيكي الجزيئي (ESP)، حددوا النقاط الدقيقة التي من المرجح أن يمسك فيها الجزيء بجزيء المذيب. ثم قاموا بحساب طاقة هذه "المصافحات" المحددة وأضافوها إلى الدرجة الإجمالية.
الصيغة السحرية
باستخدام درجتي الطاقة هاتين — واحدة للماء (EW) وواحدة للميثانول (EM) — بنى الفريق معادلة رياضية بسيطة للتنبؤ بالوقت الذي سيقضيه الجزيء على المنزلق:
tR = −0.8486 EW + 0.9396 EM + 40.55
هنا، tR هو زمن الاستبقاء (الوقت الذي يقضيه الجزيء على المنزلق).
- الإشارة السالبة لـ EW تعني أنه إذا كان الجزيء يحب الماء حقًا (طاقة منخفضة)، فإنه سينزلق من المنزلق بشكل أسرع (زمن استبقاء أقصر).
- الإشارة الموجبة لـ EM تعني أنه إذا كان الجزيء يحب الميثانول، فإنه سيبقى لفترة أطول.
كانت النتائج مبهرة. لقد فسر النموذج 92.8% من التباين في البيانات (R² = 0.928). وعندما اختبروه عن طريق ترك جزيء واحد في كل مرة لمعرفة ما إذا كان النموذج يمكنه تخمين الجزيء بشكل صحيح، كانت النتيجة لا تزال قوية عند 0.872. وبالنسبة لـ 11 جزيئًا استُخدمت لبناء النموذج، كان متوسط الخطأ 2.80% فقط. وهذا يعني أن النموذج استطاع بنجاح التنبؤ بزمن المنزلق لكل من الفلافونات والأيزوفلافونات باستخدام نفس المعادلة، مما يثبت أن موقع حلقة B لا يكسر القواعد إذا نظرت إلى درجات الطاقة الصحيحة.
الجزيء الوحيد الذي كسر القواعد
ومع ذلك، فإن القصة ليست مثالية تمامًا. اختبر الباحثون النموذج على جزيء يسمى نيوبافافلافون (neobavaisoflavone)، والذي يحتوي على مجموعة ضخمة وغير قطبية مرتبطة به (مجموعة إيزوبنتينيل). هذا الجزيء يشبه مسافرًا يرتدي حقيبة ظهر ضخمة وثقيلة مصنوعة من الشحم. لم تستطع حسابات الطاقة القياسية (التي تركز على التجاذبات الكهربائية والروابط الهيدروجينية) تفسير تفاعلات "الشحم" بشكل كامل. ونتيجة لذلك، كان تنبؤ النموذج لـ نيوبافلافون بعيدًا بنسبة 26.4%.
هذا الفشل علم الفريق شيئًا مهمًا: النموذج يعمل بشكل مثالي عندما يكون سلوك الجزيء مدفوعًا بالقوى الكهربائية والروابط الهيدروجينية. ولكن عندما يهيمن على الجزيء التأثيرات "الدهنية" الكارهة للماء وقوى التشتت (مثل تلك الحقيبة الثقيلة)، يصطدم المنهج الحالي بحائط مسدود. ليس الخطأ في الرياضيات، بل في أن نوع التفاعل المحدد لم يكن مدرجًا في درجة الطاقة.
الخلاصة
يظهر هذا البحث أننا لسنا بحاجة إلى كتاب قواعد مختلف لكل شكل مختلف قليلاً من الفلافونيدات. فمن خلال التركيز على تكلفة الطاقة الإجمالية للانتقال بين الماء والميثانول، يمكننا إنشاء نموذج موحد يتنبأ بكيفية سلوك هذه المركبات الطبيعية في المختبر. وبينما يعد النموذج دقيقًا للغاية لمعظم الجزيئات القياسية، إلا أنه يذكرنا بأن الكيمياء معقدة: إذا كان للجزيء "حقيبة ظهر" ضخمة وغير قطبية، فيجب أن نكون حذرين بشأن كيفية حساب رحلته. يوفر هذا النهج أداة قوية لتحديد المركبات النباتية دون الحاجة إلى عينة فيزيائية لكل واحدة منها، مما يساعد العلماء على كشف أسرار كيمياء الطبيعة بشكل أسرع من أي وقت مضى.
غارق في أبحاث مجالك؟
تصلك نشرة يومية بأحدث الأبحاث المطابقة لكلماتك البحثية المفتاحية — مع ملخصات تقنية، بلغتك.