Identification of the novel inhibitors against M. tuberculosis ESX-1 secretion system EccA1 enzyme using virtual screening, docking and dynamics simulation techniques

ملخص تقني مفصل لورقة بحثية بعنوان: "تحديد المثبطات الجديدة ضد نظام إفراز ESX-1 لإنزيم EccA1 في المتفطرة السلية باستخدام تقنيات الفحص الافتراضي، والالتحام الجزيئي، ومحاكاة الديناميكا الجزيئية"

1. بيان المشكلة

• السياق الصحي العالمي: لا يزال مرض السل (TB)، وخاصة السلالات المقاومة للأدوية المتعددة (MDR) والسلالات شديدة المقاومة للأدوية (XDR)، يمثل أزمة صحية عالمية حرجة. وتواجه الأدوية الحالية (الأولية والثانوية) قيوداً بسبب تزايد المقاومة.
• تحديد الهدف: يعد نظام الإفراز ESX-1 في المتفطرة السلي (Mycobacterium tuberculosis) ضرورياً للفوعة (الضراوة) والبقاء على قيد الحياة داخل الفجوة الخلوية للمضيف. وبشكل محدد، يعد إنزيم EccA1 (المكون من 573 بقايا حمض أميني، وحوالي 62.4 كيلو دالتون) مكوناً حيوياً في هذا النظام، وهو مسؤول عن إفراز عوامل الفوعة (مثل EspAC وESAT-6 وCFP-10).
• الفجوة العلاجية: يتميز إنزيم EccA1 بكونه فريداً من نوعه في الممرض (غير موجود في البشر) ومحفوظاً تطورياً، مما يجعله هدفاً مثالياً لـ العلاج المضاد للفوعة. إن تثبيط EccA1 يمكن أن يمنع إفراز عوامل الفوعة دون ممارسة نفس الضغط الانتخابي الذي يؤدي للمقاومة كما تفعل الأدوية التقليدية القاتلة للبكتيريا. ومع ذلك، لا توجد حالياً مثبطات جزيئية صغيرة محددة لـ EccA1.

2. المنهجية

استخدمت الدراسة مساراً شاملاً لاكتشاف الأدوية عبر الحاسوب (in silico) باستخدام مجموعة برامج Schrödinger وأدوات حوسبة أخرى:

• النمذجة الهيكلية:

  • تم نمذجة هيكل EccA1 كامل الطول (P9WPH_9) باستخدام AlphaFold.
  • تم التحقق من صحة النموذج عبر PROCHECK (تحليل مخطط راماشاندران)، مما أكد سلامة الكيمياء الفراغية.
  • يتكون الإنزيم من نطاق N-terminal TPR (البقايا 1–272) ونطاق C-terminal ATPase (البقايا 273–573) الذي يحتوي على زخارف Walker A وWalker B.

• الفحص الافتراضي والالتحام الجزيئي:

  • المكتبة: تم فحص 7.5 مليون مركب شبيه بالدواء من قاعدة بيانات ZINC ضد جيب ATPase في الطرف الكربوكسيلي (C-terminal).
  • توليد الشبكة: تم إنشاء شبكة مستقبلات حول البقايا الرئيسية (Pro336، Thr338، Lys340، Arg429) باستخدام وحدة Glide.
  • الالتحام (Docking): تم إلتحام أفضل المرشحات باستخدام GLIDE (وضع XP). استُخدم الركيزة الطبيعية ADP ومثبطات p97 المعروفة (NMS873 و CB5083) كضوابط للتحليل المقارن.
  • تقليل الطاقة: تم تقليل طاقة المركبات باستخدام قوة المجال OPLS-3.

• حساب طاقة الارتباط الحرة:

  • تم إجراء حسابات MM-GBSA (ميكانيكا الجزيئات - نظرية السطح العام لـ Born المعممة) باستخدام وحدة Prime لتقدير طاقات الارتباط الحرة ($\Delta G_{bind}$).

• توصيف ADMET:

  • تم تقييم الخصائص الحركية الدوائية (الامتصاص، التوزيع، الاستقلاب، الإخراج) والسمية باستخدام QikProp.
  • تم تقييم الامتثال لـ قاعدة ليبينسكي للخمسة (Lipinski's Rule of Five).

• محاكاة الديناميكا الجزيئية (MD):

  • إعداد النظام: تم إجراء عمليات محاكاة لمدة 100 نانو ثانية على apo-EccA1 والمعقدات مع أفضل 5 مركبات من ZINC (Z1–Z5)، بالإضافة إلى ADP وNMS873 وCB5083.
  • الظروف: مجموعة NPT (300 كلفن، 1 بار) باستخدام وحدة Desmond مع قوة المجال OPLS3 وماء TIP3P.
  • مقاييس التحليل: تم تحليل متوسط الانحراف للجذر التربيعي (RMSD)، ونصف قطر الدوران (Rg)، ومتوسط الانحراف للجذر التربيعي للتقلبات (RMSF)، ومساحة السطح المتاحة للمذيب (SASA)، وأنماط الروابط الهيدروجينية لتقييم الاستقرار والتغيرات الهيكلية.

3. المساهمات الرئيسية

• التحقق من صحة الهدف الجديد: تأكيد أن جيب ATPase لإنزيم EccA1 هو هدف قابل للتطبيق لتثبيط الجزيئات الصغيرة، وهو متميز عن نظائره البشرية.
• اكتشاف المركبات الرائدة: تحديد خمسة مركبات جديدة من ZINC (Z1–Z5) تتميز بخصائص ارتباط متفوقة مقارنة بالركيزة الطبيعية (ADP) ومثبطات p97 الموجودة.
• التحقق الشامل: تقديم نهج تحقق متعدد الطبقات يجمع بين الالتحام الساكن، وحسابات طاقة الارتباط الديناميكية الحرارية، وتحليل الاستقرار الديناميكي (100 نانو ثانية MD).
• التحليل المقارن: أظهرت النتائج أنه بينما تمتلك NMS873 وCB5083 (مثبطات p97) بعض الألفة، فإن مركبات ZINC المكتشفة حديثاً تظهر استقراراً وتفاعلات أكثر تحديداً داخل جيب EccA1.

4. النتائج الرئيسية

• الألفة الارتباطية ودرجات الالتحام:

  • أظهرت أفضل 5 مركبات طاقات ارتباط أعلى بكثير من ADP ($\Delta G \approx -35.0$ كيلو كالوري/مول).
  • سجل المركب Z3 أقوى طاقة ارتباط (-55.83 كيلو كالوري/مول)، يليه Z4 (-52.33)، ثم Z2 (-49.56)، ثم Z5 (-44.44)، وأخيراً Z1 (-43.45).
  • تراوحت درجات الالتحام لـ Z1–Z5 بين -10.88 و -11.33 كيلو كالوري/مول، متفوقة على ADP (-8.3)، وNMS873 (-5.0)، وCB5083 (-5.2).

• آليات التفاعل:

  • تتشابه مركبات Z1–Z5 هيكلياً مع النوكليوزيدات البيورينية (حلقات تشبه الأدينين، وسكريات الريبوز) وتُكون روابط هيدروجينية واسعة مع بقايا موتييف Walker A (Gly339، Lys340، Thr341، Thr342) وبقايا رئيسية أخرى (Ile294، Tyr466، Arg519، Arg522).
  • شكل المركب Z3 روابط هيدروجينية ثلاثية التفرع مع Thr341 وThr342 وArg522.
  • في المقابل، أظهرت NMS873 وCB5083 عدداً أقل من الروابط الهيدروجينية وتلاؤماً أقل مثالية في جيب EccA1 مقارنة بمركبات ZINC.

• استقرار الديناميكا الجزيئية (100 نانو ثانية):

  • RMSD: أظهر معقد Z5-EccA1 أكثر ملفات RMSD استقراراً (تذبذب حوالي 5–6 Å)، مما يشير إلى هيئة صلبة ومنخفضة الطاقة. أظهر apo-EccA1 مرونة عالية.
  • Rg (التماسك): حافظ معقد ADP وZ5 على أكثر الهياكل تماسكاً (Rg مستقر)، بينما أظهرت NMS873 وCB5083 تقلبات أولية كبيرة قبل الاستقرار، مما يشير إلى ارتباط أولي أضعف.
  • SASA: أظهر معقد Z5 أقل مساحة سطح متاحة للمذيب وأكثرها استقراراً (~28,000 Ų)، مما يشير إلى واجهة ارتباط متماسكة ومستقرة للغاية.
  • الروابط الهيدروجينية: حافظ معقد ADP على أعلى عدد من الروابط الهيدروجينية (8–16). ومن بين المثبطات، أظهر Z5 أكثر ملف تعريف ثابت للروابط الهيدروجينية (3–5 روابط) مع تقلب منخفض، بينما أظهر Z2 وZ3 استقراراً أقل في الروابط الهيدروجينية بمرور الوقت.

• ADMET والتشابه الدوائي:

  • استوفت جميع مركبات ZINC (Z1–Z5) قاعدة ليبينسكي للخمسة.
  • أظهرت توقعات امتصاص فموي جيدة (30–55% استهلاك بشري فموي) وملفات تعريف ذوبان مقبولة.
  • أظهرت اختراق منخفض للحاجز الدموي الدماغي (QPlogBB < -1.921)، وهو أمر مرغوب فيه عموماً لأدوية السل التي تستهدف البكتيريا داخل الخلايا دون التسبب في سمية عصبية.

5. الأهمية والاستنتاج

• إمكانات مضادات الفوعة: نجحت الدراسة في تحديد Z1–Z5 كمثبطات قوية وجديدة لـ ATPase الخاص بـ EccA1، مما يقدم استراتيجية جديدة لمكافحة السل المقاوم للأدوية المتعددة (MDR) والشديد المقاومة (XDR) عبر استهداف الفوعة بدلاً من الحيوية، مما قد يقلل من تطوير المقاومة.
• تحسين المركبات الرائدة: تمثل المركبات، وخاصة Z3 (أعلى طاقة ارتباط) وZ5 (أعلى استقرار ديناميكي)، مرشحات رائدة قوية لمزيد من التطوير.
• الاتجاهات المستقبلية: بينما تعتبر البيانات المستمدة من الحاسوب (in silico) قوية، يشدد المؤلفون على ضرورة إجراء اختبارات إنزيمية مخبرية (in vitro) ونماذج حيوانية (in vivo) للتحقق من القدرة التثبيطية والفعالية العلاجية لهذه المركبات قبل تطبيقها السريري.

يملأ هذا العمل الفجوة بين البيولوجيا الهيكلية واكتشاف الأدوية، مما يوفر إطار عمل للتصميم العقلاني لاستهداف نظام إفراز ESX-1 في المتفطرة السلية.