Structure-Informed Stress-Testing of the Standard Numbering Scheme for Class B β-Lactamases Using AlphaFold2
تؤكد هذه الدراسة متانة مخطط الترقيم القياسي "غاليني-غاراو" (BBL) لإنزيمات بيتا-لاكتاماز من الفئة B عبر 87 عائلة تم نمذجتها حديثاً باستخدام AlphaFold2، مما يؤكد قابليته للتطبيق العام مع تحديد مجموعة فرعية صغيرة من المتغيرات ذات التباينات الهيكلية التي تستوجب المتابعة التجريبية.
طورت البكتيريا دفاعاً هائلاً ضد المضادات الحيوية التي نعتمد عليها لعلاج العدوى. ومن بين أكثر أسلحتها فعالية هي إنزيمات تُسمى "بيتا-لاكتاماز"، والتي تعمل مثل مقصات جزيئية، حيث تقطع الحلقات الكيميائية للبنسلين والأدوية ذات الصلة قبل أن تتمكن من قتل البكتيريا. وبينما تستخدم بعض هذه الإنزيمات أداة كيميائية محددة للقطع، فإن مجموعة خطيرة بشكل خاص تُعرف باسم "ميتالو-بيتا-لاكتاماز" تستخدم ذرات الزنك للقيام بالمهمة نفسها. وتعد هذه الإنزيمات المعتمدة على الزنك مصدر قلق كبير في الطب الحديث لأنها تستطيع تدمير حتى أقوى مضاداتنا الحيوية التي نعتبرها الملاذ الأخير، مما يترك الأطباء بخيارات قليلة لعلاج حالات العدوى الشديدة.
لدراسة هذه الإنزيمات وتصميم أدوية جديدة لإيقافها، يحتاج العلماء إلى لغة مشتركة. ولأن هذه الإنزيمات تأتي بأشكال وأحجام مختلفة، فقد ابتكر الباحثون نظام ترقيم قياسياً. يعمل هذا النظام مثل خريطة عالمية، مما يسمح للعلماء بالإشارة إلى نفس البقعة تماماً في إنزيمات مختلفة، بغض النظر عن طول أو قصر السلسلة البروتينية. وقد رُسمت هذه الخريطة منذ عقود، بناءً على عدد صغير جداً من الهياكل الفيزيائية التي تمكن العلماء من التقاطها في المختبر. واليوم، ومع ذلك، نحن نعرف مئات الإنزيمات هذه من البكتيريا الموجودة في التربة والمياه والمستشفيات، لكننا لم نرَ أشكالها الفيزيائية قط. والسؤال الكبير هو ما إذا كانت الخريطة القديمة، التي رُسمت من بضعة أمثلة فقط، لا تزال صالحة للعمل مع هذا العالم الواسع وغير المرئي من الإنزيمات الجديدة.
لقلبت فريق من الباحثين لاختبار هذه الخريطة باستخدام أداة جديدة قوية تسمى "ألفافولد" (AlphaFold)، والتي تستخدم الذكاء الاصطناعي للتنبؤ بالشكل ثلاثي الأبعاد للبروتينات بناءً على شفرتها الجينية. وبدلاً من الانتظار لسنوات لالتقاط الهياكل الفيزيائية في المختبر، استخدم الفريق هذا البرنامج الحاسوبي لبناء نماذج لـ 87 فئة مختلفة من هذه الإنزيمات التي تستخدم الزنك والتي لم يسبق رؤيتها من قبل. بدأوا بإثبات نجاح طريقتهم، حيث قارنوا نماذجهم الحاسوبية لثمانية إنزيمات معروفة مقابل هياكلها الفعلية المقاسة تجريبياً. وكانت تنبؤات الكمبيوتر دقيقة للغاية، حيث طابقت الأشكال الحقيقية بدقة عالية لدرجة أن الفرق كان أقل من نصف عرض ذرة واحدة.
بعد التحقق من صحة طريقتهم، طبق الباحثون منهجهم على 87 عائلة غير مرئية. وتحققوا مما إذا كان نظام الترقيم القياسي يحدد بشكل صحيح المواقع الحرجة حيث تلتصق ذرات الزنك بالإنزيم. وفي الغالبية العظمى من الحالات، صمدت الخريطة القديمة بشكل مثالي. أظهرت النماذج الحاسوبية أن مواقع التصاق الزنك تتطابق تماماً كما يتوقع النظام القياسي، حتى بالنسبة للإنزيمات التي لم تُشاهد قط في المختبر. وهذا يشير إلى أن نظام الترقيم، الذي أُنشئ قبل أكثر من عشرين عاماً، قوي بما يكفي لتوجيه العلماء عبر العائلة الأكبر والأكثر تنوعاً من الإنزيمات المكتشفة منذ ذلك الحين.
ومع ذلك، وجدت الدراسة أيضاً بعض الاستثناءات المثيرة للاهتمام. ففي إحدى عشرة حالة بدا فيها أن الخريطة قد فشلت، اكتشف الباحثون سبباً بيولوجياً لهذا الاختلاف، مثل انتماء الإنزيم إلى مجموعة فرعية محددة تمتلك طبيعياً ترتيباً مختلفاً قليلاً. ولكن في عائلتين محددتين، كشفت النماذج الحاسوبية عن أمر أكثر غموضاً؛ فبينما كانت هذه الإنزيمات لا تزال تمتلك تسلسل الحروف الصحيح لبناء مواقع التصاق الزنك، أظهرت أشكالها ثلاثية الأبعاد المتوقعة أن المنطقة النشطة كانت ملتوية أو مرتبة بشكل مختلف عما اقترحته الخريطة القياسية. وتمثل هاتان العائلتان، بالإضافة إلى عائلة أخرى أظهرت نمطاً مشابهاً، مجموعة صغيرة من الإنزيمات حيث قد لا يتطابق الهيكل الفيزيائي مع التوقعات القديمة.
ويخلص الباحثون إلى أنه بينما يعد نظام الترقيم القياسي دليلاً موثوقاً لمعظم هذه الإنزيمات، فإن هذه الحالات الاستثنائية القليلة هي حالات مهمة. فهي تمثل مجموعة صغيرة من المتغيرات حيث قد يبدو الموقع النشط، وهو قلب وظيفة الإنزيم، مختلفاً في الفضاء ثلاثي الأبعاد عما كنا نعتقد. ولأن النماذج الحاسوبية قد حددت هذه الاختلافات، فإن هذه الدراسة تحدد هذه الإنزيمات كأهداف ذات أولوية عالية للأعمال المختبرية المستقبلية. وسيعرف العلماء الآن كيف يركزون جهودهم التجريبية على هذه العائلات القليلة لتأكيد أشكالها الفريدة، مما يضمن بقاء فهمنا لهذه الدفاعات البكتيرية دقيقاً بينما نواصل اكتشاف أنواع جديدة في البرية.
ملخص تقني: اختبار الإجهاد القائم على البنية لنظام الترقيم القياسي لإنزيمات بيتا-لاكتاماز من الفئة B باستخدام AlphaFold2
بيان المشكلة وُضع نظام الترقيم القياسي لإنزيمات بيتا-لاكتاماز من الفئة B (الميتالو-بيتا-لاكتاماز، MBLs)، المعروف باسم مخطط BBL، في عام 2001 وتم تحديثه في عام 2004 بناءً على مجموعة محدودة مكونة من ثمانية هياكل بلورية تم حلها تجريبياً. ومنذ ذلك الحين، تم تصنيف مئات التسلسلات الإضافية من الفئة B في قواعد بيانات مثل قاعدة بيانات البيتا-لاكتاماز (BLDB)، إلا أن الغالبية العظمى منها تفتقر إلى هياكل تجريبية. وبناءً على ذلك، يظل من غير المتحقق ما إذا كانت مراسلات البقايا (residue correspondences) المحددة بواسطة مخطط BBL — لا سيما تلك التي تحكم أنماط تنسيق الزنك الكلاسيكية وعراوي الموقع النشط التشخيصية للفئات الفرعية — تظل صحيحة عبر هذا الفضاء المتسع والمتنوع من التسلسلات. تعالج هذه الدراسة الفجوة بين الأساس الهيكلي الأصلي لنظام الترقيم والواقع الحديث للتسلسل الميتاجينومي العميق، الذي كشف عن تنوع في إنزيمات MBL يتجاوز بكثير العينة الأصلية.
المنهجية استخدم المؤلفون سير عمل يجمع بين التنبؤ بالبنية عبر التعلم العميق والتحقق الصارم وتحليل نقل البقايا:
بناء مجموعة البيانات: من قاعدة بيانات BLDB، حدد المؤلفون 87 عائلة من الفئة B (تشمل الفئات الفرعية B1 وB2 وB3) التي لا تمتلك حالياً هيكلاً تم حله تجريبياً. تضمن ذلك إجراء مطابقة مع مدخلات PDB والأدبيات العلمية لاستبعاد العائلات ذات الهياكل المحلولة حديثاً أو أخطاء التوصيف.
التحقق من صحة خط الأنابيب (Pipeline Validation): قبل التنبؤ بـ 87 هدفاً، تم التحقق من صحة خط أنابيب ColabFold (AlphaFold2/MMseqs2) مقابل مجموعة مرجعية تاريخية مكونة من ثمانية إنزيمات (BcII, CcrA, IMP-1, VIM-2, BlaB, CphA, L1, FEZ-1) وهي الهياكل المستخدمة لتعريف مخطط BBL. تم دمج النماذج المتوقعة مع نظيراتها التجريبية، مما أسفر عن انحرافات جذر التربيع لمتوسط المربعات (RMSD) تتراوح بين 0.35 و0.46 أنجستروم، مما أكد دقة خط الأنابيب.
التعامل مع الببتيد الإشاري: تمت معالجة التسلسلات المستهدفة لإزالة الببتيدات الإشارية عند الطرف النتروجيني (N-terminal) باستخدام DeepSig وSignalP 6.0. تم الاحتفاظ بالتسلسلات ذات النداءات الغامضة للببتيد الإشاري بطولها الكامل مع وضع علامة عليها كقصور في الدراسة.
التنبؤ بالبنية: تم نمذجة 8 l 87 عائلة مستهدفة باستخدام ColoreFold v1.6.3 بالإعدادات القياسية (بدون قوالب، ثلاث دورات تكرارية، واسترخاء Amber).
نقل البقايا والتحليل:
على أساس التسلسل: تم مسح أنماط ربط الزنك الكلاسيكية (مثل HxHxD للفئة B1) عبر الـ 87 تسلسلاً ناضجاً.
على أساس البنية: بالنسبة للعائلات التي كان فيها التوافق مع مراجع الفئة الفرعية غامضاً أو أنتج فجوات غير رتيبة، تم دمج الهياكل المتوقعة فوق الهياكل المرجعية التاريخية في برنامج PyMOL لرسم خرائط مواقع ربط الزنك مباشرة.
تحليل العروة (Loop Analysis): تم تقييم وجود أو غياب عروة الموقع النشط التشخيصية للفئة الفرعية (BBL 61–65) عبر الفحص المباشر للمحاذات المرجعية والتحليل الوكيل الحسابي.
النتائج الرئيسية
التحقق من الصحة: نجح خط أنابيب AlphaFold2 في إعادة إنتاج الطيات التجريبية لمجموعة الثمانية إنزيمات المرجعية بدقة عالية (RMSD < 0.5 Å)، مما يؤكد صلاحية استخدامه لاختبار مخططات الترقيم في العائلات غير المحلولة.
حفظ الأنماط (Motif Conservation): من بين 87 عائلة، أظهرت 71 عائلة (82%) تطابقاً تاماً مع نمط ربط الزنك الكلاسيكي المتوقع لفئتها الفرعية.
حل الاستثناءات: من بين 16 عائلة افتقرت في البداية إلى تطابق النمط الدقيق، تم حل 11 حالة بناءً على أسباب بيولوجية محددة:
تأكدت ثلاث منها كأعضاء في الفئة الفرعية B2 مع استبدال Asn-for-His المتوقع في موضع Zn1 الأول.
تأكدت سبع منها كأعضاء في الفئة الفرعية B3 مع استبدال Gln-for-His المسموح به في موضع Zn1 الأول.
واحدة (SPR) طابقت مجموعة "B4" المقترحة والمؤقتة سابقاً.
التباعد الهيكلي: لم تتمكن عائلتان، BB13 وBB32، من التمركز بثقة عبر محاذاة التسلسل، وأظهرتا تباعداً هيكلياً كبيراً عند الدمج (RMSDs بلغت 3.0 Å و3.7 Å على التوالي). وعلى الرغم من احتفاظهما بأنماط ربط الزنك في التسلسل، إلا أن عراوي الموقع النشط لديهما أظهرت تباعداً هيكلياً موضعياً حال دون رسم خرائط ترقيم ربط الزنك بدقة.
تأكيد العروة: أكدت الدراسة أن العروة (BBL 61–65) موجودة في B1 وغائبة في B3، وهو ما يتوافق مع المخطط الأصلي، رغم أن التحليل الوكيل الحسابي عبر المجموعة المستهدفة قد تعثر بسبب عمليات الإدراج والحذف (indels) غير المرتبطة.
مقاييس الثقة: أظهرت غالبية النماذج المتوقعة ثقة عالية (متوسط pLDDT > 90). كانت BB32 هي النموذج الأقل ثقة (pLDDT 81.3)، مع تركز الثقة المنخفضة في أجزاء معينة من السلسلة التي تقابل التباعد في الموقع النشط.
الأهمية والادعاءات يزعم البحث أن مخطط ترقيم BBL، الذي صيغ في الأصل قبل عصر التسلسل الميتاجينومي أو التنبؤ بالبنية عبر التعلم العميق، يعمم بشكل قوي على فضاء التسلسلات الأكبر والأكثر تنوعاً المتاح الآن. وتثبت الدراسة أن تنبؤات AlphaFold2 التي تم التحقق من صحتها يمكنها فعلياً اختبار ما إذا كان نظام الترقيم المبني على عينة هيكلية صغيرة يظل صالحاً للإنزيمات غير المحلولة.
تتمثل المساهمة الأساسية في تحديد مجموعة فرعية صغيرة من المتغيرات (تحديداً BB13 وBB32) التي تبدو مواقعها النشطة متباعدة هيكلياً عن الترتيب الكلاسيكي رغم احتفاظها بنمط التسلسل. ويجادل المؤلفون بأن هذه النتائج تسلط الضال على قصور محدد في أدوات إعادة الترقيم القائمة على التسلسل فقط (مثل تلك القائمة على نماذج ماركوف المخفية للملفات التعريفية)، والتي لا يمكنها اكتشاف التباعد الهيكلي الثلاثي الأبعاد الموضعي. ومن خلال تحديد هذه الاستثناءات، توفر الدراسة قائمة ذات أولوية للإنزيمات التي تستحق التوصيف الهيكلي التجريبي المباشر.
يضع المؤلفون هذا العمل كنموذج عام لتوسيع نطاق أنظمة التصنيف والترقيم المعتمدة على البنية عبر عائلات البروتينات حيث يتأخر التوصيف الهيكلي عن اكتشاف التسلسلات. كما صرحوا بوضوح أن نتائجهم تدعم الاستمرارية في استخدام مخطط BBL مع التنبيه إلى متغيرات محددة تتطلب متابعة تجريبية لحل بنية مواقعها النشطة.