Chemically Informed and Leakage-Aware QSPR Modeling: The 4L-QSPR Framework Applied to Aqueous Solubility Prediction
تقدم هذه الورقة إطار عمل 4L-QSPR، وهو سير عمل مدرك للتسرب يفصل بين التجميع الكيميائي المستقل عن النهاية وبين اختيار الواصفات الخاضع للإشراف وتحسين النموذج المتداخل لتحقيق تنبؤات بالذوبانية المائية قابلة للتعميم كيميائياً وتفسيرية بمتوسط خطأ مطلق قدره 0.459 وحدة logS.
في عالم اكتشاف الأدوية والهندسة الكيميائية، يواجه العلماء باستمرار لغزًا عمليًا: ما مدى قدرة جزيء جديد على الذوبان في الماء؟ هذه الخاصية، المعروفة باسم الذوبانية المائية، هي التي تحدد ما إذا كان الدواء سيُمتص في الجسم، أو كيف سيتصرف كيميائي في البيئة، أو ما إذا كان يمكن معالجة مادة ما بفعالية. وللإجابة على ذلك، يستخدم الباحثون طريقة تسمى نمذجة العلاقة الكمية بين البنية والخواص. وببساخات، تحاول هذه الطريقة التنبؤ بسلوك المادة الكيميائية من خلال تحليل شكلها وبنيتها، تمامًا كما قد يتنبأ ميكانيكي بكيفية تعامل السيارة مع الطريق بناءً على تصميم محركها وتوزيع وزنها. والهدف هو بناء نموذج حاسوبي يمكنه النظر إلى جزيء جديد لم يسبق رؤيته، ليخبرنا ما إذا كان سيذوب بسهولة أم سيظل صلبًا بعناد. ومع ذلك، عانت هذه النماذج لعقود من عيب خفي؛ فعندما يختبر العلماء توقعاتهم، غالبًا ما يقسمون بياناتهم عشوائيًا، مما يضع الجزيئات المتشابهة كيميائيًا في كل من مجموعة التدريب ومجموعة الاختبار. وهذا يخلق حالة لا يتعلم فيها الحاسوب قواعد الكيمياء حقًا، بل يقوم ببساطة بحفظ الإجابات على أسئلة رآها بالفعل بأشكال مختلفة قليلاً. والنتيجة هي نموذج يبدو رائعًا على الورق ولكنه يفشل عندما يواجه نوعًا جديدًا حقًا من الجزيئات.
لقد اقترح باحث يدعى بيوتر سيسيوسكي وسيلة جديدة لإصلاح هذه المشكلة، حيث قدم إطار عمل مصمم لمنع الحواسيب من الاعتماد على التشابهات الكيميائية. وفي دراسة ركزت على التنبؤ بالذوبانية المائية، طور سيسيوسكي سير عمل مكونًا من أربع خطوات يجبر النموذج على إثبات فهمه للكيمياء الأساسية بدلاً من مجرد التعرف على الأنماط المألوفة. ويتمثل جوهر هذه الطريقة الجديدة في قاعدة صارمة: قبل أن يرى الحاسوب بيانات الذوبانية، يجب عليه أولاً تنظيم جميع الجزيئات في عائلات متميزة بناءً على شكلها البنيوي فقط. ويحدث هذا التنظيم بشكل مستقل عن مفتاح الإجابة. وبمجرد تثبيت هذه العائلات، يُسمح للحاسوب فقط بالتعلم من عائلة واحدة في كل مرة، بينما يتم اختباره على عائلات مختلفة تمامًا لم يسبق له مواجهتها. وهذا يضمن أنه عندما يصدر النموذج تنبؤًا، فإنه يطبق حقًا معرفته على منطقة كيميائية جديدة، وليس مجرد الاستكمال بين جيران يعرفهم بالفعل.
طبقت الدراسة هذا النهج الصارم على مجموعة معروفة تضم أكثر من 1,100 جزيء فريد، وهي مجموعة بيانات تُستخدم غالبًا كمعيار للتنبؤ بالذوبانية. قام الباحث أولاً ببناء مكتبة مرجعية ضخمة تضم مئات الآلاف من الهياكل الكيميائية لإنشاء خريطة مستقرة للمجال الكيميائي. وباستخدام هذه الخريطة، قام بتجميع جزيئات الاختبار في 81 عنقودًا متميزًا بناءً على مدى تشابه أشكالها، متجاهلاً أي معلومات حول مدى ذوبانها. ثم استخدم خوارزمية تعلم آلي قوية لبناء نموذج تنبؤ، ولكن مع قيد حاسم: كان على الحاسوب أن يتعلم من بعض هذه العناقيد ويُختبر على عناقيد أخرى لم يرها أثناء التدريب. ولجعل النموذج أكثر قوة، دمج الباحث نوعين مختلفين من المعلومات الكيميائية؛ حيث وصفت مجموعة من البيانات الشكل ثنائي الأبعاد للارتباط والاتصال، بينما وصفت المجموعة الأخرى كيفية تفاعل الجزيء مع الماء على المستوى الفيزيائي، بما في ذلك الطاقة وتأثيرات الإذابة. ومن خلال تغذية كلا النوعين من المعلومات في النموذج، سمح له برؤية الجزيء من منظورين مختلفين في آن واحد.
أظهرت النت النتائج أن هذا النهج الأكثر صرامة والواعي بـ "تسرب البيانات" لم يجعل النموذج أسوأ؛ بل على العكس، أنتج تنبؤات دقيقة وموثوقة للغاية. فعند اختباره عبر ترتيبات عشوائية مختلفة للبيانات، حقق النموذج باستمرار متوسط خطأ يبلغ حوالي 0.46 وحدة logS (وهي مقياسة قياسية للذوبانية)، مع درجة عالية من الثقة الإحصائية. والأهم من ذلك، ظل النموذج مستقرًا بغض النظر عن كيفية خلط البيانات، مما أثبت أن نجاحه لم يكن ضربة حظ ناتجة عن تقسيم عشوائي معين. واختار الحاسوب مجموعة صغيرة يمكن إدارتها تتكون من حوالي 22 ميزة رئيسية لاتخاذ قراراته. وكانت هذه الميزات مزيجًا من الواصفات الهيكلية ومصطلحات التفاعل الفيزيائي، مما يشير إلى أن النموذج قد نجح في الجمع بين شكل الجزيء وسلوكه الفيزيائي في الماء. وقد استبعدت الدراسة صراحةً فكرة إمكانية تحقيق أداء عالٍ بمجرد الاعتماد على تقسيمات البيانات العشوائية، مظهرةً بدلاً من ذلك أن الفهم الكيميائي الحقيقي يتطلب بروتوكولات تحقق تحترم العائلات الطبيعية للجزيئات.
يبرهن هذا العمل على أنه من الممكن بناء نماذج تنبؤية تكون دقيقة وصادقة بشأن حدودها. فمن خلال فصل مهمة تجميع الجزيئات عن مهمة التنبؤ بخصائصها، أنشأ الباحث نظامًا أقل عرضة للخداع بالتشابهات السطحية. لم يكتفِ النموذج النهائي بحفظ قائمة من الإجابات، بل حدد نمطًا ثابتًا من السمات الهيكلية والفيزيائية التي تحكم الذوبانية. ويقدم هذا النهج مخططًا للدراسات المستقبلية في الكيمياء وعلوم المواد، مشيرًا إلى أن الطريقة التي نختبر بها نماذجنا لا تقل أهمية عن النماذج نفسها. فإذا أردنا الوثوق بالحاسوب للتنبؤ بسلوك دواء جديد أو مادة جديدة، يجب أن نتأكد من أنه قد تم اختباره على أرض كيميائية لم يمشِ عليها من قبل. وتخلص الدراسة إلى أنه بينما تكون التكلفة الحسابية لإنشاء هذه الأوصاف الفيزيائية التفصيلية أعلى، فإن النماذج الناتجة هي أكثر موثوقية للتطبيقات الواقعية حيث يمكن أن يكون التخمين الخاطئ مكلفًا.
ملخص تقني: نمذجة علاقة البنية بالخاصية (QSPR) القائمة على المعرفة الكيميائية والوعي بالتسرب (4L-QSPR)
بيان المشكلة غالبًا ما تعاني نماذج العلاقة بين البنية والخاصية الكمية (QSPR) من تقديرات أداء متفائلة للغاية بسبب "التسرب الهيكلي" (structural leakage). ففي بروتوكولات التقسيم العشوائي القياسية، تظهر المركبات المتشابهة هيكليًا (مثل السلاسل المتجانسة أو نظائر الهياكل الأساسية) في كل من مجموعات التدريب والاختبار. وهذا يسمح للنماذج بالاستكمال الداخلي (interpolation) ضمن جيران كيميائيين مألوفين بدلاً من التعميم على مركبات متميزة هيكليًا. وبينما توجد استراتيجيات مثل تقسيم الهياكل الأساسية (scaffold splitting) أو التقسيم القائم على التكتل (cluster-based splitting)، إلا أنها غالبًا ما تعتمد على تعريفات تعسفية أو تعتمد على الخاصية المستهدفة، مما قد يؤدي إلى إدخال التحيز. هناك فجوة منهجية لإطار عمل للتحقق يحدد مجموعات ذات معنى كيميائي بشكل مستقل عن النهاية الطرفية (endpoint)، أو ترتيب الواصفات، أو بقايا النموذج، مما يضمن أن مقاييس التحقق تقيس التعميم الكيميائي الحقيقي.
المنهجية: إطار عمل 4L-QSPR يقدم البحث إطار عمل 4L-QSPR، وهو سير عمل مكون من أربع طبقات مصمم للفصل بين التنظيم الهيكلي المستقل عن النهاية الطرفية وتحسين النموذج المعتمد على النهاية الطرفية. تم تطبيق إطار العمل على مجموعة بيانات ذوبانية مائية منقحة (ESOL) تتكون من 1,111 هوية جزيئية فريدة.
الطبقة الأولى (L1): التجميع المستقل عن النهاية الطرفية. تؤسس L1 تسلسلًا هرميًا هيكليًا قابلًا لإعادة الإنتاج دون استخدام قيم الذوبانية. وهي تستخدم مسافة هجينة ثابتة للتشابه الجزيئي تجمع بين:
التشابه الهيكلي: مسافة تانيموتو (Tanimoto distance) المشتقة من بصمات مورغان (Morgan fingerprints) بنصف قطر 2 وبطول 2048 بت.
التشابه الفيزيائي الكيميائي: مسافة جيب التمام (Cosine distance) في فضاء معياري من ثمانية واصفات RDKit (مثل TPSA، وMolLogP، وNumHBD).
الأوزان: نسبة ثابتة قدرها 0.80 (بصمة) إلى 0.20 (واصفات).
التقسيم: إجراء أخذ عينات MaxMin قائم على المعالم (landmark-based) على كون مرجعي يضم حوالي 342 ألف جزيء لتوليد تكتل هرمي. بالنسبة لدراسة ESOL، تم اختيار التقسيم hier_level_100 الذي نتج عنه 81 مجموعة هيكلية بمتوسط حجم يبلغ حوالي 13.7 جزيء.
الطبقة الثانية (L2): الاختيار المسبق للواصفات الخاضع للإشراف. تقوم L2 بعملية تنظيف واختيار مسبق محدد للمنظم (regressor) والخاص بالهدف. فهي تقوم بتصفية الواصفات غير الصالحة، أو الثابتة، أو المتكررة، وتطبق إجراء اختيار الاستقرار. بالنسبة لدراسة ESOL، تضمن ذلك توليد ثلاثة فضاءات من الواصفات:
PaDEL: واصفات هيكلية ثنائية الأبعاد.
المشتقة من COSMO: مصطلحات فيزيائية كيميائية من حسابات COSMO-RS (معاملات النشاط، الجهود الكيميائية، طاقات التفاعل في الماء/الهكسان/الأوكتانول الرطب).
المدمجة: دمج واصفات PaDEL وCOSMO. تقوم L2 بتقليص هذه الفضاءات إلى مجمعات مرشحة مدمجة (أفضل 100) للطبقة التالية، باستخدام القيم المستهدفة فقط للتصفية، وليس لتعريف المجموعات.
الطبقة الثالثة (L3): التحسين المتداخل المدرك للمجموعات (DOOIT2). تستخدم L3 خوارزمية DOOIT2 (التحسين متعدد الأهداف مع تقليم الميزات المتكرر) ضمن مخطط تحقق StratifiedGroupKFold (SGKF) متداخل.
هيكل التحقق: تستخدم الحلقة الخارجية المجموعات الهيكلية من L1 كمتغير التجميع لضمان عدم وجود تداخل هيكلي بين طيات (folds) التدريب والاختبار. تتعامل الحلقة الداخلية مع ضبط المعلمات الفائقة وتقليم الواصفات.
التحسين: باستخدام Optuna مع خوارزمية NSGA-II، يقوم الإجراء بتحسين هدفين: تقليل متوسط الخطأ المطلق (MAE) وتقليل تعقيد النموذج.
التقليم المتكرر: بدءًا من مجمع L2 المختار، يتم إزالة الواصفات الأقل أهمية بشكل تكراري بناءً على الأداء في الحلقة الداخلية حتى الوصول إلى الحد الأدنى لعدد الواصفات (المحدد بـ 4).
التقييم النهائي: الطية الخارجية محجوزة بصرامة؛ حيث تُسجل مقاييس الأداء فقط بعد انتهاء الحلقة الداخلية من تحديد مجموعة الواصفات والمعلمات الفائقة.
الطبقة الرابعة (L4): محجوزة للنشر والمراقبة في المستقبل (لم يتم استكشافها في هذه الدراسة).
النتائج الرئيسية تم تقييم إطار العمل باستخدام منظم LightGBM على فضاء واصفات COSMO-PaDEL المدمج. أكدت الدراسة على المتانة من خلال تشغيل سير عمل L3 الكامل باستخدام ثلاث بذور عشوائية مستقلة (7، 42، 123).
مقاييس الأداء:
متوسط البذور: MAE = 0.459 ± 0.004 وحدة logS؛ R² = 0.897 ± 0.007.
تشغيل نموذجي (البذرة 42): MAE الطية الخارجية = 0.460 ± 0.046؛ R² = 0.895 ± 0.031.
عدد الواصفات: اختار النموذج مجموعة مدمجة بلغت في المتوسط 22.2 ± 5.8 واصفًا.
قابلية تفسير الواصفات: شملت الواصفات الأكثر تكرارًا في الاختيار عبر جميع الطيات والبذور مزيجًا من:
المشتقة من COSMO: مصطلحات الإذابة (مثل ln_gamma_water وln_gamma_wet_octanol)، والجهود الكيميائية، وتباينات طاقة التفاعل (مثل ΔH_int_wet_octanol_minus_water).
PaDEL: واصفات كره الماء (مثل XLogP)، والواصفات الطوبولوجية/الرسومية (مثل SpMax3_Bhm وpiPC4)، والواصفات الإلكترونية (مثل MLFER_E). يدعم هذا التكرار الفرضية القائلة بأن الجمع بين الطوبولوجيا الهيكلية وفيزياء الإذابة الصريحة ينتج نماذج قوية وقابلة للتفسير كيميائيًا.
التحقق من L1: تم التحقق من تجميع L1 نوعيًا مقابل تصنيف دروند (Durand) للمذيبات، مما أظهر أن التسلسل الهرمي المستقل عن النهاية الطرفية قد صنف بنجاح العائلات المرتبطة كيميائيًا (مثل الأمينات الحلقية، والنتريلات، والكحولات) دون الوصول إلى بيانات الذوبانية.
الأهمية والادعاءات يزعم البحث أن المساهمة الرئيسية لهذا العمل هي منهجية وليست مجرد إنشاء متنبئ عالمي للذوبانية أو دخولًا في قائمة المتصدرين.
التحقق الواعي بالتسرب: تظهر الدراسة أنه يمكن تحقيق نماذج QSPR تنافسية وقابلة للتفسير تحت بروتوكولات تحقق أكثر صرامة تعتمد على المجموعات الكيميائية التي تمنع التسرب الهيكلي. وتشير النتائج إلى أن الأداء العالي تحت التقسيم العشوائي غالبًا ما يخفي عدم قدرة النموذج على التعميم إلى هياكل كيميائية متميزة.
الفصل الصريح بين المهام: من خلال فصل التجميع الهيكلي (L1) عن التعلم الخاضع للإشراف (L2/L3)، يوفر إطار العمل حدودًا شفافة بين التقسيم الكيميائي وملاءمة النموذج. وهذا يضمن أن مقاييس التحقق تعكس التعميم إلى مناطق هيكلية متميزة بدلاً من مجرد الاستكمال الداخلي.
تكامل الواصفات: تؤكد النتائج أن نمذجة الذوبانية المائية تستفيد من فضاء واصفات هجين. ويشير تكرار كل من مصطلحات الإذابة المشتقة من COSM والواصفات الهيكلية من PaDEL إلى أن أي مصدر منهما بمفرده ليس كافيًا للتعميم الأمثل تحت التحقق الصارم.
قابلية إعادة الإنتاج: تم تصميم سير العمل ليكون حتميًا فيما يتعلق بتقسيم L1 ومتينًا تجاه التباينات العشوائية في تدريب النموذج، كما يتضح من التباين المنخفض عبر البذور العشوائية.
يخلص المؤلف إلى أن التحقق من QSPR يجب أن يُعامل كمشكلة تقسيم كيميائي، وليس مجرد مهمة إعادة أخذ عينات إحصائية. يوفر إطار عمل 4L-QSPR مسارًا عمليًا لتطوير نماذج ليست دقيقة فحسب، بل وشفافة أيضًا فيما يتعلق بقدرتها على التعميم عبر الفضاء الكيميائي.