Integrated Machine Learning, Molecular Docking, and Molecular Dynamics Simulations Design Novel Peptide Inhibitors of Rabies Virus Glycoprotein
تستخدم هذه الدراسة نهجاً حوسبياً متكاملاً يجمع بين التعلم الآلي، والرسو الجزيئي، ومحاكاة الديناميات الجزيئية لتصميم والتحقق من صحة مثبطات ببتيدية جديدة غير سامة مشتقة من سم "ألفا-بونغاروتوكسين"، حيث حددت المرشح P40 كعامل علاجي واعد يستهدف البروتين السكري لفيروس داء الكلب.
تخيل جسم الإنسان كمدينة صاخبة، والفيروسات كلصوص متسللين يحاولون الاقتحام. ولكي يتمكن هؤلاء اللصوص من الدخول، يحتاجون إلى مفتاح رئيسي لفتح الباب الأمامي. في حالة فيروس داء الكلب، هذا "الباب الأمامي" هو مستقبل محدد على خلايانا العصبية، و"المفتاح الرئيسي" هو بروتين على سطح الفيروس يسمى البروتين السكري (أو RABV-G). هذا المفتاح فعال للغاية لدرجة أنه يحاكي شكل سم الأفعى، مما يخدع دفاعات الجسم ويسمح للفيروس بالتسلل مباشرة إلى الجهاز العصبي. وبمجرد دخوله، يسبب الفيروس مرضاً يكاد يكون مميتاً في معظم الحالات، مما يجعل السباق لإيجاد طريقة لتعطيل ذلك القفل مسألة حياة أو موت. لقد عرف العلماء منذ زمن طويل أنه إذا استطاعوا تصميم "مفتاح مزيف" أو "سدادة" تناسب قفل الفيروس تماماً ولكنها لا تفتح الباب، فيمكنهم إيقاف العدوى قبل أن تبدأ. وهنا يلتقي عالم علوم الحاسوب بالبيولوجيا: استخدام خوارزميات قوية لتصميم قطع جزيئية صغيرة، تسمى الببتيدات، تعمل كـ "سدادات" مثالية.
في هذه الدراسة، قرر فريق من الباحثين من معهد هافكين في مومباي لعب لعبة "تيتريس" الجزيئية القصوى. بدأوا بمخطط من مصدر مختلف تماماً: سم موجود في سم أفعى "الكرايت"، المعروف باسم "ألفا-بونغاروتوكسين" (α-bungarotoxin). يشتهر سم الأفعى هذا بقدرته على الارتباط بنفس المستقبلات العصبية التي يستخدمها فيروس داء الكلب. أدرك الباحثون أن الفيروس وسم الأفعى يتنافسان أساساً على نفس الموقع. لذا، طرحوا سؤالاً ذكياً: "ماذا لو استطعنا أخذ الجزء من سم الأفعى الذي يمسك بالقفل، ولكن مع تجريده من السم القاتل، وتحويله إلى درع غير ضار؟"
باستخدام مزيج من الأدوات الحاسوبية المتقدمة، صمم الفريق مكتبة من هذه "الدروع غير الضارة". أولاً، استخدموا إطار عمل هجين يجمع بين التعلم الآلي والتعلم العميق — فكر في الأمر كمرشح رقمي فائق الذكاء — لفحص آلاف التصاميم المحتملة. تم تدريب هذا المرشح على رصد واستبعاد أي تصميم قد يظل ساماً أو ضاراً للبشر، لضمان انتقال المرشحين الآمنين فقط إلى المرحلة التالية. ثم استخدموا "الرسو الجزيئي" (molecular docking)، وهو ما يشبه محاكاة افتراضية عالية السرعة حيث حاولوا ملاءمة كل ببتيد في قفل الفيروس لمعرفة أي منها يتلاحم بقوة أكبر.
كانت نتائج هذا الفحص الرقمي مثيرة للاهتمام. فقد حدد الباحثون مرشحاً بارزاً، وهو ببتيد أطلقوا عليه اسم P40 (بالتسلسل MKMKCHTTATMKIV). في محاكاتهم الحاسوبية، ناسب P40 البروتين السكري لفيروس داء الكلب بشكل أفضل من أي تصميم آخر اختبروه. لقد شكل رابطة قوية ومستقرة، مما أدى فعلياً إلى "تعطيل" القفل. وللتأكد تماماً من أن هذا لم يكن مجرد مصادفة في المحاكاة، أجروا محاكاة "الديناميكا الجزيئية" لمدة 100 نانو ثانية. تخيل هذا كفيلم عالي السرعة يظهر تفاعل الببتيد مع الفيروس بمرور الوقت؛ أظهر الفيلم أن P40 حافظ على موقعه، وظل ملتصقاً بقوة بالفيروس دون اهتزاز أو سقوط، حتى مع تحرك واهتزاز الجزيئات.
تشير الدراسة إلى أن هذا الببتيد المهندس، والمستمد من سم أفعى ولكن بعد تجريده من خطورته، يمكن أن يكون سلاحاً جديداً قوياً ضد داء الكلب. وبينما تؤكد الورقة البحثية أن هذه النتائج تعتمد حالياً على النماذج والمحاكاة الحاسوبية، فإن النتائج تشير إلى أن P40 يعد نموذجاً واعداً للغاية. إنها خطوة مفعمة بالأمل نحو مستقبل قد نمتلك فيه دواءً قائماً على الببتيدات يمكنه إيقاف فيروس داء الكلب في مساره، محولاً القفل المميت إلى لغز غير ضار لم يعد الفيروس قادراً على حله.
ملخص تقني: تصميم مثبطات ببتيدية مبتكرة لبروتين الجليكوبروتين الخاص بفيروس داء الكلب باستخدام التعلم الآلي المتكامل، والرسو الجزيئي، ومحاكاة الديناميكا الجزيئية
بيان المشكلة لا يزال داء الكلب أحد أكثر الأمراض الفيروسية حيوانية المنشأ فتكاً في العالم، حيث يتسبب في حوالي 55,000 حالة وفاة سنوياً، وتحدث معظمها في آسيا وأفريقيا. يسببه فيروس داء الكلب (RABV)، وهو فيروس RNA عصبي، حيث يعمل بروتين الجليكوبروتين (RABV-G) كمحدد رئيسي للتعرف على الخلايا المضيفة والدخول إليها. يتفاعل RABV-G مع مستقبلات المضيف، بما في ذلك مستقبل الأسيتيل كولين النيكوتيني (nAChR)، وهي آلية مشتركة مع السم العصبي "ألفا-بانغاراتوكسين" (α-bungarotoxin). وبينما يعد RABV-G هدفاً تم التحقق منه لتطوير مضادات الفيروسات وتصميم اللقاحات، إلا أن هناك حاجة ماسة لوكلاء علاجيين مبتكرين يمكنهم استهداف التفاعل بين الفيروس والمستقبل بشكل محدد دون التسبب في السمية العصبية المرتبطة بالسموم الطبيعية. تعالج هذه الدراسة تحدي تصميم مثبطات ببتيدية آمنة وغير سامة تحاكي أنماط ربط المستقبلات لـ α-bungarotoxin ولكنها تفتقر إلى خصائصه السامة.
المنهجية استخدمت الدراسة مسار عمل حاسوبي متعدد المراحل يدمج هندسة الببتيدات، والتعلم الآلي (ML)، والتعلم العميق (DL)، ومحاكاة البيولوجيا الهيكلية:
نمذجة الهدف: تم إنشاء نموذج هيكلي ثلاثي الأبعاد لبروتين RABV-G باستخدام SWISS-MODEL (القالب 6lgx.1.A) وتم التحقق من صحته عبر تحليل مخطط راماشاندران، مما أكد وجود 92.5% من البقايا في المناطق المفضلة.
تصميم مكتبة الببتيدات: بدءاً من التسلسل الأحماض الأمينية لـ α-bungarotoxin، تم إنشاء مكتبة من نظائر الببتيد. تضمنت استراتيجية التصميم تجزئة تسلسل السم وتطبيق طفرات نقطية منهجية ضمن مجموعات فيزيائية كيميائية محددة (حمضية، قاعدية، كارهة للماء، إلخ) للحفاظ على أنماط الربط مع القضاء على السمية العصبية.
فحص السمية (إطار عمل ML-DL): تم استخدام إطار عمل هجين يجمع بين التعلم الخاضع للإشراف (الغابة العشوائية، وآلات ناقلات الدعم، وK-أقرب جار) والتعلم العميق (الشبكات العصبية العميقة مع المشفرات التلقائية) لفحص الببتيدات بناءً على الواصفات الفيزيائية الكيميائية (الكراهية للماء، الشحنة، مؤشر بومان). قام هذا النهج المزدوج بتصفية 163 مرشحاً محتملاً من السمية من أصل مجموعة أولية، مما أبقى 40 ببتيداً غير سام.
التوصيف الفيزيائي الكيميائي ومضادات الفيروسات: تم تحليل المرشحين الأربعين الـ 40 من حيث الاستقرار، والذوبان، والنشاط المضاد للفيروسات المتوقع (IC50) باستخدام أدوات مثل ProtParam وAVP-IC50Pred.
الرسو الجزيئي: تم إجراء رسو بروتين-ببتيد عالي الغموض باستخدام خوادم HADDOCK وHDOCK لتحديد أوضاع الربط وتقارب الارتباط مع جيب ربط RABV-G. تم تحديد البقايا الرئيسية في جيب الربط (مثل Glu272، Glu279، Asp283) عبر DoGSiteScorer.
محاكاة الديناميكا الجزيئية (MDS): خضع أفضل المرشحين لمحاكاة MDS لمدة 100 نانو ثانية باستخدام وحدة Desmond. تم تقييم الاستقرار من خلال تحليل الانحراف في متوسط مربع الجذر (RMSD)، وتقلب متوسط مربع الجذر (RMSF)، وتحليل الروابط الهيدروجينية، وتحليل المكونات الرئيسية (PCA)، ومشهد طاقة الحرية (FEL)، ومصفوفة الارتباط المتقاطع الديناميكي (DCCM).
النتائج الرئيسية
تحديد المرشح الرائد: من بين المرشحين الذين تم تقييمهم، برز P40 (MKMKCHTTATMKIV) كأكثر المثبطات واعدية. أظهر أعلى درجة HADDOCK بلغت -90.0 ± 3.8، وطاقة ارتباط حرة (ΔG) قدرها -10.3 كيلو كالوري/مول، وثابت تفكك متوقع (Kd) قدره 2.7×10−8 مول. يحدد الملخص صراحةً P40 كأفضل مرشح، يليه P3 (CLDLGYTIVCHTTATSPI) بدرجة -86.3، ثم P25 (MKIVCHTTMKSPMK) بدرجة -86.0، وP21 (MKIVCHTTATSPMKV) بدرجة -75.0.
المرشحون الثانويون: أظهرت الببتيدات P3 وP25 وP21 أيضاً درجات رسو وتقاربات ربط تنافسية. ومن الجدير بالذكر أن P25 أظهر درجة HADDOCK متفوقة (-86.0) مقارنة بـ P21 (-75.0)، مما يجعله في المرتبة الثانية أو الثالثة كأقوى رابط اعتماداً على المعيار المحدد (HADDOCK مقابل طاقة الارتباط)، ولكنه يتفوق باستمرار في تقارب الرسو على P21.
آلية التفاعل: كشف تحليل الرسو والتفاعل أن الببتيدات الرائدة ترتبط بجيب محدد في RABV-G يتضمن بقايا رئيسية هي: Glu272، Glu279، Asp283، Lys276، وArg297. يتم تثبيت الارتباط من خلال مزيج من التفاعلات الكهروستاتيكية (بين بقايا الببتيد القاعدية مثل Lys/His والبقايا البروتينية الحمضية) والاتصالات الكارهة للماء.
الاستقرار الديناميكي: أشارت نتائج MDS إلى أن معقد P40-RABV-G حقق استقراراً هيكلياً سريعاً، محافظاً على قيم RMSD منخفضة (2–4 Å) وقيم RMSF طوال محاكاة الـ 100 نانو ثانية. أظهر تحليل مشهد طاقة الحرية أن P40 شكل حوض طاقة عميقاً ومحدداً جيداً، مما يشير إلى حالة تشكيلية مستقرة للغاية مقارنة بالبروتين غير المرتبط أو المعقدات مع الببتيدات الأخرى.
علاقة البنية بالنشاط (SAR): حددت الدراسة أن الببتيدات ذات نقاط التعادل الكهربائي المتوسطة إلى العالية، والتي تحتوي على زخارف CHTT محددة، وطول يترا-وح 14–16 حمضاً أمينياً، أظهرت أقوى إمكانات مضادة للفيروسات متوقعة.
الأهمية والادعاءات يدعي البحث نجاح استراتيجية حاسوبية متكاملة للتصميم de novo لمثبطات ببتيدية ضد بروتين الجليكوبروتين لفيروس داء الكلب. تكمن أهمية العمل في:
السلامة بالتصميم: النجاح في القضاء على السمية العصبية من تسلسلات مشتقة من α-bungarotoxin مع الاحتفاظ بالقدرة على استهداف مستقبل الفيروس.
التكامل المنهجي: الاستخدام الفعال لإطار عمل هجين (ML-DL) لفحص السمية، مقترناً بالتحقق الهيكلي الصارم ومحاكاة الديناميكا، مما يوفر مساراً قوياً لاكتشاف الببتيدات المضادة للفيروسات.
الإمكانات العلاجية: يشير تحديد P40 (MKMKCHTTATMKIV) كمرشح ارتباط مستقر وعالي الألفة إلى أنه مرشح رائد قابل للتطبيق لتطوير علاجات ببتيدية لداء الكلب. يشير الملخص إلى أن الببتيدات المشتقة من α-bungarotoxin يمكن أن تستهدف RABV-G وأن P40 (المصحح من الخطأ المطبعي "P440" في الملخص) يمكن أن يكون مرشحاً رائداً واعداً.
التوجه المستقبلي: يذكر المؤلفون أن هذه النتائج توفر "أساساً مناسباً" للتحقق التجريبي المستقبلي، وتحديداً من خلال الإشارة إلى الحاجة إلى اختبارات الارتباط in vitro، ودراسات النشاط المضاد للفيروسات، والتحقيقات in vivo لتأكيد التنبؤات الحاسوبية.
تخلص الدراسة إلى أن الببتيدات المهندسة المشتقة من α-bungarotoxin، وتحديداً P40، تمتلك إمكانات كبيرة لعرقلة دخول الخلايا المضيفة بوساطة RABV-G، مما يوفر مساراً جديداً لعلاج داء الكلب.