Computational Design and In Silico Evaluation of Indole-2-Carboxamide Derivatives as Mycobacterium tuberculosis Membrane Protein (MmpL3) Inhibitors Using 2D-QSAR, Molecular Docking, and Molecular Dynamics Simulation
استخدمت هذه الدراسة نهجاً حوسبياً متكاملاً، شمل نماذج العلاقة بين البنية والنشاط الكمي ثنائية الأبعاد (2D-QSAR)، والرسو الجزيئي، ومحاكاة الديناميكا الجزيئية، لتصميم وتقييم مشتقات إندول-2-كربوكساميد كمثبطات محتملة لبروتين MmpL3، مما أدى إلى تحديد المركب 17c كمرشح واعد يتميز بألفة ارتباط فائقة وخصائص حركية دوائية ملائمة رغم الحاجة إلى التحقق التجريبي مستقبلاً.
المؤلفون الأصليون:Thomas Aondofa Nyijime, Gideon Adamu Shallangwa, Adamu Uzairu, Abdullahi Bello Umar, Muhammad Tukur Ibrahim
لا يزال السل أحد أكثر التحديات الصحية استمراراً في العالم، وهو عدوى بكتيرية لا تزال تحصد الأرواح رغم وجود علاجات قائمة. البكتيريا المسؤولة عن ذلك، وهي "المتفطرة السلية" (Mycobacterium tuberculosis)، قد طورت آلية دفاع هائلة: جدار خلوي شمعي سميك يحميها من البيئة ومن العديد من الأدوية. ولبناء هذا الجدار، تعتمد البكتيريا على آلة جزيئية محددة تسمى MmpL3، والتي تعمل كـناقل متخصص، حيث تقوم بنقل اللبنات الأساسية عبر غشاء الخلية. وإذا توقف هذا الناقل، فلن تتمكن البكتيريا من بناء غلافها الواقي، وتموت الخلية. لعقود من الزمن، سعى العلماء لإيجاد طرق جديدة لتعطيل هذه الآلة، آملين في العثور على سلاح فعال ضد السلالات التي أصبحت مقاومة للأدوية الحالية. وغالباً ما يتجه البحث نحو هياكل كيميائية تُعرف باسم "إندول-2-كربوكساميد" (indole-2-carboxamides)، وهي عائلة من المركبات التي أظهرت وعداً في تعطيل هذه العملية الحيوية.
في دراسة حديثة، شرع الباحثون في تحسين هذه الهياكل الكيميائية باستخدام محاكاة حاسوبية قوية بدلاً من أنابيب الاختبار. بدأوا بمجموعة مكونة من ستة وأربعين جزيئاً موجوداً من مركبات "إندول-2-كربوكساميد"، وحللوا كيف ترتبط أشكالها وخصائصها الإلكترونية بقدرتها على إيقاف ناقل MmpL3. ومن خلال رسم هذه العلاقات، بنوا نموذجاً تنبؤياً يمكنه تخمين مدى فعالية جزيء جديد غير مختبر بناءً على تصميمه. عمل هذا النموذج كمخطط، حيث وجه الفريق لتعديل مركب قوي بشكل خاص، يُعرف بالمركب 17، لإنشاء أربع نسخ جديدة محسنة. كان الهدف هو تعديل أطراف الجزيء بحيث يتناسب بشكل أكثر إحكاماً مع الموقع النشط للناقل، تماماً مثل ضبط مفتاح ليدور في القفل بسلاسة أكبر.
كشفت عمليات المحاكاة الحاسوبية أن أحد التصاميم الجديدة، والذي أُطلق عليه اسم المركب 17c، برز بشكل ملحوظ. فعندما وضع الباحثون هذا الجزيء في نموذج افتراضي لبروتين MmpL3، استقر في جيب الارتباط بقوة تفوق كلاً من المركب الأصلي وأدوية السل القياسية مثل "آيزونيازيد" و"إيثامبوتول". أظهرت المحاكاة أن المركب 17c شكل شبكة معقدة من الروابط مع البروتين، بما في ذلك التجاذبات الكهربائية والاتصالات الكارهة للماء التي ثبتته في مكانه بإحكام. وللتأكد من أن هذا لم يكن مجرد لقطة ثابتة، أجرى الفريق محاكاة ديناميكية استمرت مائتي نانو ثانية، لمراقبة كيفية سلوك الجزيء والبروتين بمرور الوقت. وخلال هذه الفترة الممتدة، ظل المركب الجديد مستقراً، يتحرك ضمن نطاق ضيق ويحافظ على قبضته على الهدف، بينما أظهرت الأدوية المرجعية حركة أكبر وارتباطاً أضعف.
علاقت تحليل إضافي الطاقة المطلوبة لإبقاء الدواء مرتبطاً بالبروتين، وهو مقياس لمدى قوة تمسكهما. أشارت النتائج إلى أن المركب 17c يتطلب طاقة أكبر بكثير للانفصال عن البروتين مقارنة بالأدوية القديمة، مما يشير إلى تفاعل أقوى وأكثر ديمومة. كما تحقق الباحثون من ملف السلامة المتوقع لهذا الجزيء الجديد. أشارت النماذج الحاسوبية إلى أنه سيتم امتصاصه جيداً من قبل الأمعاء البشرية، وهي سمة حاسمة للدواء الذي يؤخذ عن طريق الفم، وأنه لن يسبب على الأرجح ضرراً جينياً فورياً أو تحسساً جلدياً. وبينما أكدت الدراسة أن هذا الجزيء هو مرشح واعد للغاية في المجال الرقمي، فإن المؤلفين يلاحظون بحذر أن هذه النتائج تستند بالكامل إلى عمليات المحاكاة. الاختبار الحقيقي لن يأتي إلا عندما يتم تصنيع هذه المركبات في المختبر واختبارها ضد البكتيريا الحية للتأكد من أن تنبؤات الكمبيوتر تصمد في العالم الحقيقي. وحتى ذلك الحين، يظل المركب 17c نموذجاً أولياً مهذباً ومحسناً نظرياً، ويقدم اتجاهاً واضحاً للعمل التجريبي المستقبلي في مكافحة السل المقاوم للأدوية.
ملخص تقني: التصميم الحسابي والتقييم في السيليكوو (In Silico) لمشتقات إندول-2-كاربوكساميد كمثبطات لـ MmpL3
بيان المشكلة لا يزال السل (TB) يمثل تحدياً صحياً عالمياً حرجاً، يتفاقم بسبب الظهور السريع للسلالات المقاومة للأدوية المتعددة (MDR) والسلالات المقاومة للأدوية الشاملة (XDR) من المتفطرة السلية (Mycobacterium tuberculosis). أصبحت العلاجات الحالية من الخط الأول غير فعالة بشكل متزايد ضد هذه السلالات المقاومة، مما يستدعي تحديد أهداف استراتيجية وعلاجات جديدة. وقد برز بروتين الغشاء الميكوبكتيري الكبير 3 (MmpL3)، وهو ناقل ضروري مسؤول عن تصدير الترهالوز مونومايكولات (TMM) لتخليق حمض الميكوليك، كهدف واعد. وبينما أظهرت مشتقات إندول-2-كاربوكساميد نشاطاً تثبيطياً قوياً لـ MmpL3، إلا أن العديد من المرشحات الحالية تعاني من قيود مثل الخصائص الفارماكوكينيتيكية الضعيفة أو السمية الخلوية. علاوة على ذلك، تعتمد العديد من الدراسات الحالية على نهج حسابي أحادي الطريقة، مما قد يحد من دقة التنبؤ. وهناك حاجة إلى دراسات حسابية متكاملة وشاملة لتقييم مشتقات إندول-2-كاربوكساميد بشكل منهجي من حيث ألفة الارتباط، والخصائص الإلكترونية، والسلوك الفارماكوكينيتيكي لتوجيه تصميم عوامل مضادة للسل محسنة.
المنهجية استخدمت الدراسة استراتيجية حسابية موحدة تدمج تقنيات متعددة لتصميم وتقييم مشتقات إندول-2-كاربوكساميد:
مجموعة البيانات والتحسين: تم استخدام مجموعة بيانات مكونة من 46 مشتقاً من إندول-2-كاربوكساميد ذات قيم IC₅₀ تجريبية مسجلة. تم تحسين البنى الجزيئية باستخدام نظرية الكثافة الوظيفية (DFT) عند مستوى B3LYP/6-31G*.
نمذجة 2D-QSAR: تم تقسيم مجموعة البيانات إلى مجموعتين: تدريب (70%) واختبار (30%). استُخدم نهج الخوارزمية الجينية - الانحدار الخطي المتعدد (GA-MLR) لتطوير نموذج 2D-QSAR. تم التحقق من قوة النموذج باستخدام التحقق المتبادل "ترك واحد خارج" (LOO-CV)، والتحقق الخارجي، واختبارات العشوائية Y.
الإرساء الجزيئي (Molecular Docking): أُجريت دراسات الإرساء ضد بروتين MmpL3 (PDB ID: 7NVH) باستخدام Molegro Virtual Docker. تم التحقق من البروتوكول عن طريق إعادة إرساء اللجين المشترك (Av0). حللت الدراسة تفاعلات الارتباط لمجموعة البيانات وأربعة نظائر مصممة حديثاً (17a–17d) بناءً على المركب الأفضل أداءً 17.
توصيف ADMET: تم تقييم خصائص الشبه بالدواء والفارماكوكينيتيكا باستخدام SwissADME وpkCSM، لتقييم معايير مثل الامتصاص المعوي، ونفاذية الحاجز الدموي الدماغي (BBB)، والاستقرار الأيضي، والسمية.
محاكاة الديناميكا الجزيئية (MD Simulations): أُجريت عمليات محاكاة ديناميكا جزيئية لمدة 200 نانو ثانية باستخدام وحدة Desmond مع حقل القوة OPLS4 لتقييم الاستقرار الديناميكي للمعقدات المتكونة من المركب الرائد (17)، وأفضل النظائر المصممة (17c)، والأدوية المرجعية (آيزونيازيد وإيثامبوتول).
تحليل MM-GBSA: تم حساب طاقات الارتباط الحرة باستخدام طريقة ميكانيكا الجزيئات - مساحة سطح غيبس الموجهة بالنموذج المستمر (MM-GBSA) لتقدير استقرار الارتباط.
المساهمات والنتائج الرئيسية
أداء نموذج QSAR: أظهر نموذج 2D-QSAR المطور قدرة تنبؤية مرضية مع معامل تحديد داخلي (R2) قدره 0.646، ومعامل تحقق متبادل (QLOO2) قدره 0.527، ومعامل تنبؤ خارجي R2 قدره 0.636. حدد النموذج واصفات رئيسية (GATS6v, GATS5i, SpMax7_Bhe) تحكم النشاط، مما يشير إلى أن الضخامة الفراغية، والجهد التأيني، والتفرع الجزيئي تؤثر بشكل كبير على قوة التثبيط.
الإرساء الهيكلي والتحسين: حدد الإرساء الجزيئي المركب 17 كركيزة ملائمة بدرجة MolDock بلغت -146.758. وبناءً على ذلك، تم تصميم أربعة نظائر. ومن بينها، أظهر المركب 17c أعلى ألفة ارتباط بدرجة MolDock بلغت -159.019، متفوقاً على المركب الأصلي والأدوية المرجعية (آيزونيازيد: -121.305؛ إيثامبوتول: -90.576). كشف تحليل التفاعل أن 17c يشكل تفاعلات كهروستاتيكية، وروابط هيدروجينية، وتفاعلات كارهة للماء واسعة النطاق داخل جيب ارتباط MmpL3.
ADMET والتشابه مع الدواء: استوفت جميع المركبات التي تم فحصها، بما في ذلك النظائر المصممة، بشكل عام قاعدة "ليبنشكي" للقواعد الخمس، مما يشير إلى توافق جيد مع الخصائص الدوائية الفموية. أشارت تنبؤات ADMET إلى امتصاص معوي بشري عالٍ (>90%) لجميع المركبات. ومن الجدير بالذكر أن المركب 17c أظهر امتصاصاً عالياً (93.82%) وكان يُتوقع أن يكون ركيزة لإنزيم CYP3A4، ولكنه غير مطفر وغير مسبب لحساسية الجلد.
الاستقرار الديناميكي (محاكاة MD): كشفت محاكاة MD لمدة 200 نانو ثانية أن معقد 17c-MmpL3 حافظ على استقرار هيكلي عالٍ، حيث تراوحت قيم RMSD للعمود الفقري للبروتين بين 2.4 و2.8 آنجستروم، وقيم RMSD لللجين بين 2.4 و4.8 آنجستروم. كان هذا أكثر استقراراً من المعقد الناتج عن المركب الأصلي 17، الذي أظهر تقلبات أعلى في البداية. حافظ المركب 17c على روابط هيدروجينية مستمرة وتفاعلات كارهة للماء مع بقايا الأحماض الأمينية الرئيسية (مثل GLU-424، ARG-448، ILE-422).
طاقة الارتباط الحرة (MM-GBSA): أكد تحليل MM-GBSA استقرار الارتباط الفائق للمركب 17c، حيث حقق طاقة ارتباط حرة (ΔGbind) قدرها -76.27 كيلو كالوري/مول. وقد تفوق هذا بشكل كبير على المركب 17 (-68.59 كيلو كالوري/مول) والأدوية المرجعية آيزونيازيد (-22.04 كيلو كالوري/مول) وإيثامبوتول (-36.88 كيلو كالوري/مول). وأشار تحليل تفكك الطاقة إلى أن تفاعلات فان دير فالس والتفاعلات الكارهة للماء كانت القوى الدافعة الأساسية لهذا الاستقرار.
الأهمية والادعاءات يدعي البحث أن النهج الحسابي المتكامل نجح في تحديد المركب 17c كمرشح رائد واعد لتثبيط MmpL3. وتفترض الدراسة أن التعديلات الهيكلية التي أُدخلت على 17c عززت من ألفة الارتباط والاستقرار الديناميكي مقارنة بالركيزة الأصلية والأدوية المضادة للسل القياسية. ويؤكد المؤلفون أن المركب 17c يمتلك خصائص ارتباط جيدة، وخصائص فارماكوكينيتيكية متوقعة (امتصاص معوي عالٍ)، وملف سمية آمن، مما يجعله مرشحاً قوياً لمزيد من التطوير.
القصور والتواضع يقر المؤلفون صراحةً بأن الدراسة محدودة بغياب التحقق التجريبي. فجميع النتائج تعتمد على التنبؤات الحاسوبية (in silico). وبناءً على ذلك، يخلص البحث بتواضع إلى أنه بينما يظهر المركب 17c إمكانات تثبيط واعدة حسابياً، فإن الدراسات المخبرية (in vitro) والدراسات الحيوية (in vivo) ضرورية تماماً لتأكيد النشاط البيولوجي المتوقع، والإمكانات العلاجية، والسلوك الفارماكوكينيتيكي قبل التطبيق السريري. كما وُصف الأداء التنبئي لنموذج QSAR بأنه "مرضٍ" ولكنه "متوسط"، ويُعزى ذلك جزئياً إلى حجم مجموعة البيانات البالغ 46 مركباً.