Integrated Machine Learning and Structure-Based Computational Prioritization of KIF11 Inhibitor Candidates for CNS-Oriented Drug Discovery
تقدم هذه الدراسة سير عمل حوسبي متكامل يعتمد على التعلم الآلي والأساس الهيكلي، نجح في تحديد خمسة مرشحات متنوعة هيكلياً لتثبيط بروتين KIF11 ذات قدرة تنبؤية عالية ونافذية عالية للحاجز الدموي الدماغي لاستخدامها في اكتشاف أدوية الأورام الموجهة للجهاز العصبي المركزي.
السرطان هو مرض يتميز بالنمو غير المنضبط، حيث تنقسم الخلايا دون توقف. ولإيقاف هذه العملية الجامحة، يبحث العلماء غالبًا عن طرق لتعطيل الآلية التي تسمح للخلايا بالانقسام إلى اثنتين. أحد الأجزاء الحيوية في هذه الآلية هو بروتين يسمى KIF11، والذي يعمل كمحرك جزيئي لتنظيم السقالات الداخلية اللازمة لانقسام الخلايا. وعندما يتم حظر هذا البروتين، لا تستطيع الخلية الانقسام بشكل صحيح وتموت في النهاية. وهذا يجعل KIF11 هدفًا واعدًا لأدوية السرطان الجديدة. ومع ذلك، فإن علاج أنواع السرطان التي تنشأ في الدماغ، مثل الورم الأرومي الدبقي، يمثل تحديًا فريدًا. فالدماغ محمي بحاجز بيولوجي يمنع معظم المواد من الدخول، مما يعني أن الدواء يجب أن يُصمم خصيًا ليعبر هذا الجدار ليصل إلى الورم. إن العثور على جزيء يكون قويًا بما يكفي لإيقاف السرطان وصغيرًا بما يكفي لدخول الدماغ هو عملية توازن صعبة للباحثين.
في دراسة حديثة، طور فريق من العلماء من الهند استراتيجية جديدة تعتمد على الحاسوب لإيجاد عقاقير محتملة يمكنها حل هذه المشكلة. فبدلاً من اختبار آلاف المواد الكيميائية في المختبر واحدة تلو الأخرى، وهو أمر بطيء ومكلف، استخدموا برامج حاسوبية متقدمة لمحاكاة كيفية تفاعل الجزيئات المختلفة مع بروتين KIF11. بدأ الباحثون بجمع مجموعة ضخمة من البيانات الموجودة حول كيفية تأثير المواد الكيميائية المختلفة على KIF11، وقاموا بتنقية المعلومات لضمان دقتها واتساقها. ثم علموا نظامًا حاسوبيًا كيفية التعرف على الأنماط في أشكال هذه الجزيئات وخصائصها الكيميائية، مما سمح للنظام بالتنبؤ بأي المواد الكيميائية الجديدة وغير المختبرة التي قد تكون مثبطات فعالة. ولضمان موثوقية تنبؤاتهم، استخدموا طريقة صارمة حيث تعلم الحاسوب من مجموعة من المواد الكيميائية وتم اختباره على مجموعة منفصلة تمامًا لم يسبق له رؤيتها من قبل، مما منع النظام من مجرد حفظ الإجابات.
نجح النموذج الحاسوبي في تحديد العديد من المرشحين الواعدين، لكن الباحثين لم يتوقفوا عند هذا الحد. فقد أرادوا معرفة ما إذا كان بإمكان هذه الجزيئات أن تتناسب بالفعل داخل الجيب المحدد في بروتين KIF11 حيث يتعين عليها الارتباط للقيام بعملها. وباستخدام تقنية تسمى "الإرساء الجزيئي" (molecular docking)، قاموا بمحاكاة التفاعل الفيزيائي بين أفضل المرشحين وبنية البروتين. هذه العملية تشبه محاولة إدخال مفتاح في قفل لمعرفة ما إذا كان الشكل متطابقًا، ولكن يتم ذلك بالكامل داخل الحاسوب. أظهرت عمليات المحاكاة أن الجزيئات المختارة يمكنها بالفعل الاستقرار في المكان الصحيح على البروتين. كما أجرى الفريق فحوصات لمعرفة ما إذا كانت هذه الجزيئات ستتمكن من عبور الحاجز الدموي الدماغي والوصول إلى الدماغ، وكذلك ما إذا كانت ستكون آمنة للتعامل معها من قبل الجسم. ووجدوا أن أفضل المرشحين لديهم احتمالية عالية جدًا لدخول الدماغ، وهو مطلب حاسم لعلاج أورام الدماغ.
من بين الجزيئات التي تم تحديدها، برز واحد كأكثرها واعدًا. هذا المركب، المعروف بالرمز CHEMBL392148، تم التنبؤ بأنه سيكون فعالاً للغاية في إيقاف بروتين KIF11 ولديه فرصة ممتازة لعبور الدماغ. نظر الباحثون أيضًا عن كثب في التركيب الكيميائي لهذا الجزيء وغيره للتأكد من أنها لا تحتوي على أجزاء معروفة بسميتها أو عدم استقرارها. وبينما أظهر أحد المرشحين علامة تحذير صغيرة تتعلق بسمة كيميائية محددة، بدا الخيار الأول نظيفًا وآمنًا بناءً على النماذج الحاسوبية. كما سلطت الدراسة الضوء على مصدر قلق محتمل يتعلق بالسلامة: حيث أظهرت جميع المرشحات الخمسة الأفضل إشارة إلى أنها قد تؤثر على النظام الكهربائي للقلب، وهي مخاطرة تتطلب الاختبار في مختبر حقيقي قبل أي استخدام بشري.
يؤكد الباحثون أن هذه النتائج هي نقطة انطلاق وليست منتجًا نهائيًا. فالجزيئات التي حددوها لا تزال مجرد تنبؤات حاسوبية. الخطوة التالية ستكون تصنيع هذه المواد الكيميائية في المختبر واختبارها على بروتينات وخلايا حقيقية لمعرفة ما إذا كانت ستتصرف كما توقع الحاسوب. توفر الدراسة مسارًا واضحًا ومنطقيًا للمضي قدمًا، حيث تقدم قائمة مختصرة من الهياكل الكيميائية المحددة التي تستحق المزيد من البحث. ومن خلال الجمع بين التعلم الآلي والتحليل الهيكلي المفصل، نجح الفريق في تضييق نطاق البحث عن دواء لسرطان الدماغ إلى حفنة من المرشحين عالي الجودة، مما وفر الوقت والموارد مع التركيز الجهود على الحلول الأكثر احتمالية.
ملخص تقني: التكامل بين التعلم الآلي والترتيب الحسابي القائم على البنية لمرشحات مثبطات KIF11 من أجل اكتشاف أدوية موجهة للجهاز العصبي المركزي
بيان المشكلة لا يزال الورم الأرومي الدبقي (Glioblastoma - GBM) يمثل تحدياً كبيراً في مجال طب الأورام نظراً لطبيعته العدوانية، ومقاومته العلاجية، ومحدودية فعالية العلاج. وبينما يعد بروتين Kinesin Family Member 11 (KIF11/Eg5) هدفاً معتمداً للعلاج المضاد للانقسام الفتيلي عبر تعطيل تجميع المغزل ثنائي القطب، فإن تحديد مثبطات جديدة قادرة على عبور الحاجز الدموي الدماغي (BBB) من أجل اكتشاف أدوية موجهة للجهاز العصبي المركزي يتطلب استراتيجيات حسابية صارمة. تعالج هذه الدراسة الحاجة إلى تقليص الحيز الكيميائي للاختبار التجريبي من خلال دمج التعلم الآلي (ML) مع الطرق القائمة على البنية لتحديد أولويات مثبطات KIF11 ذات النفاذية المتوقعة للجهاز العصبي المركزي.
المنهجية استخدمت الدراسة سير عمل حسابي منهجي متعدد المراحل:
تنسيق البيانات وتقسيمها:
تم استرجاع مجموعة بيانات مكونة من 1,817 سجلاً لنشاط KIF11 من ChEMBL (الهدف CHEMBL4581).
طُبقت معايير تنسيق صارمة: الاحتفاظ بقياسات IC50 الدقيقة، وإزالة التحذيرات المتعلقة بـ "النطاق غير النموذجي"، واستبعاد المقايسات الخاصة بالطفرات، وإزالة المقايسات القائمة على الخلايا لضمان الحصول على نهاية بيوكيميائية على مستوى الهدف.
تكونت مجموعة البيانات المنسقة النهائية من 831 مركباً فريداً (نطاق pIC50 يتراوح بين 4.046 و9.854).
تم تقسيم مجموعة البيانات إلى مجموعة تدريب (654 مركباً) ومجموعة اختبار خارجية مغلقة (177 مركباً) قبل تطوير النموذج. ومن الأهمية بمكان أن التقسيم كان منفصلاً من حيث الهيكل الكيميائي (scaffold-disjoint) باستخدام هياكل Bemis–Murcko، مما يضمن عدم وجود هياكل مشتركة بين مجموعتي التدريب والاختبار لمنع تسرب السلاسل الكيميائية.
نمذجة العلاقة الكمية بين البنية والنشاط (QSAR) والتحقق من صحتها:
الموصِّفات (Descriptors): تم توليد أحد عشر موصِّفاً فيزيائياً كيميائياً وبصمات جزيئية بنظام ECFP4 (Extended Connectivity Fingerprint 4) بسعة 2048 بت.
الخوارزميات: تم تقييم الغابة العشوائية (RF)، والأشجار الإضافية (ET)، وتعزيز التدرج المستند إلى المدرج التكراري (HGB) باستخدام تمثيلات هجينة، وتمثيلات تعتمد فقط على الموصِّفات، وتمثيلات تعتمد فقط على البصمات.
استراتيجية التحقق:
التحقق المتقاطع الواعي بالهيكل (Scaffold-aware cross-validation): لضمان بقاء المركبات التي تشترك في نفس الهيكل ضمن نفس الطية.
التحقق المتقاطع المتداخل (Nested scaffold cross-validation): لفصل اختيار النموذج عن تقدير الأداء لتوفير تقدير داخلي محافظ.
العشوائية لمتغير Y (Y-randomization): أُجريت 10 دورات لاستبعاد الارتباطات العرضية.
نطاق التطبيق (AD): تم تحديده باستخدام تشابه Morgan ECFP4 Tanimato (العتبة: 0.5701) لتقييم تغطية الحيز الكيميائي.
التحقق الخارجي المغلق: تم تطبيق النموذج النهائي المجمد على الـ 177 مركباً المستقلة دون مزيد من التحسين.
الفحص الافتراضي والتحليل القائم على البنية:
تنبأ نموذج QSAR بقيم pIC50 للمجموعة الخارجية.
حددت عملية اختيار قائمة على التنوع (خوارزمية جشعة، تشابه Max Tanimoto < 0.70) 20 مرشحاً عالي الإمكانات.
الإرساء الجزيئي (Molecular Docking): أُجري باستخدام AutoDock Vina 1.2.7 ضد بروتين KIF11 البشري (PDB ID: 2WOG). تم التحقق من البروتوكول عن طريق إعادة إرساء اللجين المشترك في البلورة STLC. ورغم عدم تحقيق إعادة إنتاج وضعية البلورة بدقة (< 2 Å RMSD) بسبب تناظر اللجين وحدود أخذ العينات، إلا أن الإرساء نجح في تحديد موقع اللجينات داخل جيب الارتباط المحدد تجريبياً.
توصيف التفاعل: تم تحليل البقايا (residues) ضمن نطاق 5 أنجستروم من وضعيات الإرساء لتحديد أنماط التفاعل الشائعة.
توصيف ADMET والتشابه الدوائي:
خضع خمسة من أفضل المرشحين لتنبؤ شامل لـ ADMET باستخدام ADMET-AI، مع التركيز على اختراق الحاجز الدموي الدماغي (BBB)، والتوافر الحيوي، ونشاط البروتين السكري P-glycoprotein، وتفاعلات CYP450، ونهايات السمية (hERG، AMES، DILI، إلخ).
تم تطبيق فحص التنبيهات الهيكلية (PAINS, BRENK, NIH) ومقاييس التشابه الدوائي (QED, Lipinski rules).
أُجريت محاكاة استرخاء المذيب الأولية لمدة 1 نانو ثانية (1 ns MD simulation) لمعقد ممثل (KIF11–CHEMBL195232) لتقييم الاستقرار الهيكلي، والتي أُبلغ عنها بشكل منفصل في المعلومات التكميلية.
النتائج الرئيسية
أداء النموذج:
تفوقت النماذج الهجينة (التي تجمع بين الموصِّفات وECFP4) على النماذج التي تعتمد فقط على الموصِّفات. حقق نموذج الأشجار الإضافية الهجين أفضل أداء أساسي (R2 = 0.5985).
التحقق المتقاطع المتداخل للهيكل: حقق التحقق الداخلي الصارم R2 قدره 0.5471، وRMSE قدره 0.9021، وMAE قدره 0.7040.
التحقق الخارجي المغلق: حققت مجموعة الاختبار المستقلة R2 قدره 0.2671، وRMSE قدره 0.8754، وMAE قدره 0.7384. يعكس انخفاض R2 الخارجي التحدي المتمثل في التنبؤ بالنشاط للهياكل الكيميائية المتميزة وحقيقة أن 48% فقط من المركبات الخارجية وقعت ضمن نطاق التطبيق المحدد مسبقاً.
العشوائية لمتغير Y: أنتجت النماذج العشوائية قيم R2 سالبة (المتوسط -0.3593)، مما يؤكد أن النموذج يلتقط علاقات حقيقية بين البنية والنشاط وليس مجرد ارتباطات عرضية.
المرشحات ذات الأولوية: تم تحديد خمسة مركبات متنوعة هيكلياً بناءً على درجة متكاملة تجمع بين تنبؤات QSAR، ودرجات الإرساء، وتشابه AD، وملفات ADMET التعريفية:
CHEMBL392143: أعلى قدرة متوقعة (pIC50 = 7.75؛ IC50 ≈ 17.77 نانومولار)، وأعلى احتمالية لاختراق الحاجز الدموي الدماغي (0.98)، ولا توجد تنبيهات هيكلية.
CHEMBL1086410: أفضل درجة إرساء (-9.882 كيلو كالوري/مول) واحتمالية عالية لاختراق الحاجز الدموي الدماغي (0.92).
CHEMBL383402: أعلى تشابه AD (0.80) ولكنه احتوى على تنبيه واحد من نوع BRENK (سلسلة أليفاتية طويلة).
CHEMBL400042: أعلى قيمة QED (0.87) وملف فيزيائي كيميائي متوازن.
تفاعلات الإرساء: احتلت جميع المركبات ذات الأولوية جيب ارتباط STLC. لوحظت تفاعلات متكررة مع البقايا Arg119، وGlu116، وTyr211، وLeu214، وAla218، وArg221.
نتائج ADMET: أظهرت جميع المرشحات الخمسة احتمالات عالية لاختراق الحاجز الدموي الدماوي (0.91–0.98) وتوافراً حيوياً عن طريق الفم. ومع ذلك، لوحظت إشارة حسابية متسقة لمسؤولية hERG عبر جميع المرشحين الخمسة (الاحتمالات 0.71–0.96)، مما يشير إلى خطر محتمل على سلامة القلب يتطلب تحققاً تجريبياً.
الأهمية والادعاءات يدعي البحث تقديم نقطة انطلاق حسابية عقلانية لاكتشاف مثبطات KIF11 المصممة خصيصاً لأورام الجهاز العصبي المركزي. وتؤكد الدراسة أن المركبات المحددة هي مرشحات ذات أولوية حسابية، وليست مثبطات مؤكدة تجريبياً.
تشمل المساهمات الرئيسية ما يلي:
تطوير سير عمل QSA robust (قوي) وواعٍ بالهيكل يقلل من تسرب السلاكل الكيميائية ويوفر إطار تحقق خارجي صارم.
دمج نهج يعتمد على اللجين (QSAR) ونهج يعتمد على البنية (الإرساء) مع توصيف ADMET للترشيح من أجل نفاذية الجهاز العصبي المركزي.
تحديد خمسة كيانات كيميائية (بقيادة CHEMBL392148) تستوفي معايير متعددة لتثبيط KIF11 واختراق الحاجز الدموي الدماغي.
يصرح المؤلفون صراحةً أن الدراسة لا تثبت النشاط التثبيطي أو السلامة. كما تم تسليط الضوء على مسؤولية hERG المكتشفة ومحدودية بروتوكول الإرساء (عدم القدرة على إعادة إنتاج الوضعيات البلورية الدقيقة) كجوانب حرجة تتطلب تحقراً تجريبياً. يعمل هذا العمل كإطار لتوليد الفرضيات لعمليات التقييم البيوكيميائي اللاحقة، والتقييم الخلوي في نماذج GBM، والدراسات الحركية الدوائية.