Predicted and experimental protein-ligand coordinates with distance labels and evaluation splits
تقدم هذه الورقة PLI-Parallax، وهو عبارة عن مجموعة بيانات شاملة تدمج بين هياكل البروتين واللجين المتوقعة والتجريبية مع تسميات دقة مُعايرة وتقسيمات تقييم صارمة، مما يتيح تقدير موثوقية التنبؤ من خلال توافق الطرق حتى في غياب الحقيقة التجريبية الأرضية.
في العالم المجهري لعلم الأحياء، تعتمد الحياة على سلسلة مستمرة من المصافحات. فالبروتينات، وهي الآلات الجزيئية المعقدة التي تبني خلايانا وتديرها، يجب أن تتعرف على شركاء كيميائيين محددين وترتبط بهم، والذين يُطلق عليهم غالبًا اسم "الربيطات" (ligands)، للقيام بمهامها. لعقود من الزمن، اعتمد العلماء على تقنية بلورات الأشعة السينية لالتقاط صور عالية الدقة لهذه المصافحات، عبر تجميد البروتين وشريكه في مكانه لرؤية كيفية تطابقهما تمامًا. هذه الصور هي المعيار الذهبي، لكن التقاطها صعب ومكلف. ونتيجة لذلك، بالنسبة للغالبية العظمى من أزواج البروتين والربيطة التي يعلم العلماء بوجودها، لا توجد صورة فوتوغرافية. ولملء هذه الفجوة، لجأ الباحثون إلى برامج حاسوبية تحاول التنبؤ بشكل هذه المعقدات من الصفر. هذه الأدوات قوية، لكنها تعاني من نقطة ضعف كبيرة: يمكنها إنتاج شكل ما، لكنها غالبًا لا تستطيع إخبارك بمدى الثقة في هذا الشكل. قد ينتج الحاسوب نموذجًا يبدو مثاليًا ولكنه خاطئ تمامًا، أو نموذجًا يبدو فوضويًا ولكنه صحيح في الواقع، مما يترك العلماء دون وسيلة موثوقة للحكم على النتائج.
يعالج مورد جديد يسمى "PLI-Parallax" حالة عدم اليقين هذه عن طريق تحويل مشكلة التنبؤ إلى اختبار للتوافق. فبدلاً من الاعتماد على برنامج حاسوبي واحد لإعلان صحة النتيجة، قام هذا المشروع بتشغيل أربعة محركات تنبؤ مختلفة على نفس المجموعة من أزواج البروتين والربيطة. تضمنت هذه المحركات أداتين حديثتين للتنبؤ بالبنية تعتمدان على التعلم من قواعد بيانات ضخمة للهياكل البيولوجية المعروفة، ومحرك إرساء (docking) تقليدي قائم على الفيزياء يحسب كيفية تطابق الذرات بناءً على القوى الفيزيائية. طبق الباحثون هذه الأدوات على مجموعة هائلة تضم أكثر من 51,000 نظام. ومن الأهمية بمكان أنهم قسموا هذه المجموعة إلى مجموعتين: المجموعة الأولى، التي تحتوي على ما يقرب من 20,000 نظام، تألفت من أزواج توجد لها بالفعل صور تجريبية حقيقية؛ مما سمح للفريق بالتحقق من التنبؤات الحاسوبية مقابل الحقيقة المعروفة. أما المجموعة الثانية، التي تضم أكثر من 31,000 نظام، فقد تألفت من أزواج لا توجد لها صور تجريبية. وفي هذه الحالات المجهولة، لم يستطع الباحثون التحقق من الإجابة مباشرة، لذا استخدموا الدروس المستفادة من المجموعة الأولى لتقدير مدى موثوقية التنبؤات.
الاكتشاف الجوهري لهذا العمل هو أنه عندما تتفق طرق التنبؤ المختلفة مع بعضها البعض، فإن ذلك يعد إشارة قوية على أن النتيجة من المرجح أن تكون صحيحة. ومن خلال مقارنة مواقع الذرات عبر النماذج الحاسوبية المختلفة، وجد الباحثون أن الأنظمة التي أنتجت فيها النماذج أشكالاً متشابهة كانت أكثر عرضة للمطابقة مع الواقع التجريبي من الأنظمة التي اختلفت فيها النماذج. لقد استخدموا البيانات التجريبية المعروفة لمعايرة هذه الإشارة، مما أدى إلى إنشاء نظام تسجيل يمنح تقديرًا للموثوقية لكل تنبؤ على حدة. وهذا يعني أنه بالنسبة لآلاف أزواج البروتين والربيطة التي لم تُجرَ عليها أي تجربة من قبل، أصبح لدى العلماء الآن وسيلة لتقييم جودة التنبؤ. لا يوفر هذا المورد الإحداثيات المتوقعة فحسب، بل يوفر أيضًا خريطة للمسافات بين الذرات ودرجة ثقة تخبر المستخدم بمدى إمكانية الوثوق في الهندسة البنائية.
أودع المشروع أكثر من 300 مليون سجل مسافة، حيث رصد المسافات الدقيقة بين كل ذرة في الربيطة وكل جزء من البروتين. وتسمح هذه السجلات للعلماء الآخرين بإعادة فحص البيانات وفقًا لمعاييرهم الخاصة. كما نظم الباحثون البيانات في مجموعات اختبار محددة لضمان أن النتائج لم تكن مجرد حفظ لأمثلة قديمة؛ حيث فصلوا بعناية بين بيانات التدريب وبيانات الاختبار لضمان اختبار النماذج على بروتينات ومواد كيميائية جديدة وغير مرئية سابقًا. وكشف هذا الإعداد الصارم أنه بينما تعد أدوات التنبؤ مثيرة للإعجاب، إلا أنها ليست مثالية. فقد كان أداء الأدوات أفضل في البروتينات الصغيرة والمواد الكيميائية البسيطة، بينما واجهت صعوبة أكبر مع الهياكل الضخمة والمعقدة أو تلك التي تتضمن أيونات معدنية. وقد ثبت أن التوافق بين الطرق المختلفة هو مؤشر أكثر موثوقية للدقة من درجات الثقة الداخلية التي توفرها أي أداة منفردة.
لا يدعي هذا المشروع حل مشكلة التنبؤ بالبروتين والربيطة، ولا يشير إلى أن النماذج الحاسوبية معصومة من الخطأ. بدلاً من ذلك، فإنه يقدم إطارًا عمليًا للتعامل مع عدم اليقين. ومن خلال توفير مجموعة بيانات ضخمة حيث يتم قياس وتوصيف موثوقية كل تنبؤ بشكل صريح، يسمح "PLI-Parallax" للباحثين بتصفية التخمينات منخفضة الجودة والتركيز على التخمينات عالية الثقة. إنه يحول مجموعة من الأشكال الخام غير المتحقق منها إلى مجموعة بيانات مهيكلة وقابلة للاستخدام حيث يتم تحديد مستوى الثقة بوضوح. وبالنسبة للعلماء الذين يدرسون كيفية تفاعل الأدوية مع الجسم، أو كيفية عمل الإنزيمات، فإن هذا يعني أنه يمكنهم الآن استخدام هذه البنى المتوقعة مع فهم أوضح لحدودها. إن هذا العمل يمثل جسرًا بين العالم المعروف للهياكل التجريبية والمجال الشاسع وغير المستكشف للبيولوجيا المتوقعة، موفرًا الأدوات اللازمة للإبحار في المجهول بدقة أكبر.
ملخص تقني لـ PLI-Parallax
بيان المشكلة تنتج أدوات التنبؤ بالبنية (مثل نماذج التكوين المشترك ومحركات الإرساء) الآن هندسات البروتين-الربيطة للأزواج التي تفتقر إلى الدقة التجريبية. ومع ذلك، فإن هذه الأدوات عادةً ما توفر درجات ثقة غير معايرة مقابل الدقة الملحوظة؛ فعلى سبيل المثال، غالبًا ما تتركز قيم الثقة في نموذج Boltz-2 فوق 0.9 دون تمييز بين التوضعات الصحيحة والخاطئة. وبينما يُعد الاتفاق بين الطرق المستقلة إشارة على الموثوقية، لم تكن هناك آلية لاستخلاص القيمة الكمية لمستوى محدد من الاتفاق. الموارد الحالية إما تفتقر إلى المركبات المتنبأ بها، أو لا تربط التنبؤات بالمركبات المقاسة للتقييم، أو تفشل في توفير تقديرات موثوقية لكل نظام تمت معايرتها مقابل الحقيقة التجريبية المقاسة.
المنهجية يقدم المؤلف PLI-Parallax، وهو مجموعة بيانات تودع إحداثيات البروتين-الربيطة المتنبأ بها والتجريبية جنباً إلى جنب مع تسميات المسافة وتقسيمات التقييم. يتضمن البناء ثلاث مراحل:
اختيار المصدر وتحديد الطبقات:
طبقة البلورات (Crystal Tier): 19,350 مركبًا تجريبيًا تم اختيارها من BioLiP2 وProtein Data Bank (PDB). تعمل هذه كحقيقة مطلقة للمعايرة.
طبقة المتن (Corpus Tier): 31,746 زوجًا تم إجراء عملية إرساء متقاطع لها (cross-docked) بناءً على مستقبلات AlphaFold حيث لا يوجد مركب تجريبي. تعمل هذه كهدف للتعلم الخاضع للإشراف حيث تغيب الحقيقة المطلقة.
كون البروتين (Protein Universe): موروث من معيار ذهبي للتفاعلات بين بروتين الإنسان وبروتين آخر (Bernett et al.)، تم تقسيمه للتحكم في تسرب التسلسل ولكن لم يتم تحسينه لتنوع موقع الارتباط.
تكوينات المعلم (Teacher Configurations): تم تشغيل أربع تكوينات "معلم" متميزة على مدخلات مشتركة:
Chai-1: خطوة إعادة تدوير جذع واحدة، مع تفعيل تضمينات ESM.
Boltz-2 (Single-Sequence): بدون شرط المحاذاة، خطوة إعادة تدوير واحدة.
Boltz-2 (Alignment-Conditioned): محاذات ColFold MMseqs2، خطوة إعادة تدوير واحدة في طبقة البلورات وثلاث خطوات في طبقة المتن التجريبية.
smina: محرك إرساء قائم على الفيزياء (نسخة من AutoDock Vina) يستخدم صندوقًا مكعبًا بمقدار 22 أنغستروم.
ملاحظة: تتنبأ أذرع التكوين المشترك ببنيات المستقبل الخاصة بها، بينما يستخدم ذراع الإرساء المستقبل المودع.
اشتقاق التسميات والتعليق:
جداول المسافة: مشتقة من الإحداثيات عند حد قطع 15.0 أنغستروم، وتسجل الحد الأدنى لمسافة الذرة الثقيلة (dmin) ومسافة Cα لأزواج ذرات البقايا والربيطة. تُشتق أعلام الاتصال المنطقية (Boolean) عند 4 و5 و8 أنغستروم.
التعليق على الموثوقية: يتم حساب الاتفاق بين المعلمين كمتوسط مؤشر جاكارد الزوجي لمجموعات بقايا اتصال 4 أنغستروم بين Chai-1 وBoltz-2 (تسلسل واحد) وذراع smina المعاد توليده. يقوم الانحدار الأيزوتوني (Isotonic regression) بربط هذا الاتفاق بالدقة الملحوظة في طبقة البلورات. كما يتم حساب فاصل زائف (split-conformal interval) لتوفير دقة تسمية متنبأ بها لكل نظام مع حد ثقة.
تقسيمات التقييم: تم توفير 646 تكوينًا (تدريب/تحقق/اختبار) عبر سبعة عائلات (C1, C2p, C2l, C3, C3comp, temporal, lead-optimisation). تستخدم هذه التقسيمات التجميع (MMseqs2 للبروتينات، وButina للربيطات) للتحكم في تسرب التسلسل والهيكل، مع الإبلاغ عن درجات مصنف ثنائي العينة للتحقق من الفصل.
المساهمات الرئيسية وسجلات البيانات
النطاق: تغطي المجموعة 51,096 نظامًا، مما ينتج 307,314,646 سجل مسافة من بقايا إلى ذرة الربيطة.
التنسيق: تُخزن البيانات في مخزن إحداثيات HDF5 (5.01 جيجابايت) وجداول Parquet (1.14 جيجابايت). مخزن الإحداثيات ممركز المتوسط ومكمم بدقة 0.01 أنغستروم.
المعايرة: تحمل 30,567 نظامًا (تقاطع المعلمين الثلاثة) حقل دقة تسمية متنبأ به مع فاصل زائف.
النزاهة: تتضمن مجموعة البيانات بيانًا (manifest)، وبيانات وصفية بتنسيق Croissant، وقاموس حقول (235 حقلًا) لضمان قابلية التكرار والقراءة الآلية.
النتائج والتحقق التقني
الاتفاق مقابل الدقة: في طبقة البلورات، يترابط الاتفاق بين المعلمين مع الدقة الملحوظة. على سبيل المثال، ترتفع دقة مجموعات الاتصال مقابل الحقيقة المطلقة من 0.2175 (مع معلم واحد) إلى 0.9356 (مع ثلاثة معلمين).
أداء المعايرة: يحقق الانحدار الأيزوتوني متوسط خطأ مطلق (MAE) قدره 0.0805 على تقسيم المعايرة المحجوز، وهو ما يتفوق بشكل كبير على الخط المرجعي الثابت (MAE 0.1752).
الصلاحية الفيزيائية: باستخدام PoseBusters، حقق ذراع إرساء smina المعاد توليده كسر صلاحية قدره 0.835 في طبقة البلورات، متفوقًا على التحكم التجريبي (0.633) وأذرع التكوين المشترك (0.248–0.532). ويُعزى انخفاض صلاحية أذرع التكوين المشترك إلى أخطاء الإدراك في استعادة الروابط/التكافؤ من الإحداثيات المتنبأ بها بدلاً من الفشل الهندسي.
مقاييس الدقة: في طبقة البلورات، وضع Chai-1 نسبة 32.1% من الربيطات ضمن 2 أنغستروم RMSD، بينما وضع smina نسبة 28.9%. ومع ذلك، استعاد smina المزيد من بقايا الاتصال (تداخل جاكارد 0.5592) مقارنة بـ Chai-1 (0.4442)، مما يوضح أن المقاييس المختلفة ترتب المعلمين بترتيب مختلف.
الفصل في التقسيم: نجحت التقسيمات المقدمة في فصل الكيانات التدريبية والاختبارية. في طبقة المتن، تتبع تقسيمات محور البروتين خط الأساس المعتمد على الطول فقط (AUC ~0.50)، مما يشير إلى أن تقسيم البروتينات لا يجعلها قابلة للتمييز بطبيعتها بما يتجاوز الحجم، بينما تظهر تقسيمات محور الربيطة فصلًا قويًا (AUC 0.82–0.90).
الأهمية والادعاءات يضع البحث PLI-Parallax ليس كمعيار لإعلان "أفضل" طريقة، بل كمورد لتمكين معايرة موثوقية التنبؤ.
سد الفجوة: يعالج نقص مقاييس الدقة المعايرة لكل نظام للمركبات المتنبأ بها من خلال استخدام طبقة البلورات لـ "تحديد قيمة مستوى معين من الاتفاق".
الإشراف على البيانات غير المقاسة: من خلال نقل المعايرة من طبقة البلورات إلى طبقة المتن، تتيح المورد استخدام طبقة المتن كإشراف مرجح بالاتفاق، بدلاً من كونها مجرد هندسة غير موسومة.
نطاق متواضع: يذكر المؤلف صراحة أن تكوينات المعلم موجهة نحو الإنتاجية وليست محسنة لدقة الهدف الواحد. تمثل الإحداثيات ما تنتجه هذه الإعدادات المحددة، وليس الحد الأقصى النظري للطرق.
القيود: يقر المورد بأن تقاطع المعلمين الثلاثة منحاز نحو الربيطات الأصغر وغير المعدنية ويقل فيه تمثيل روابط الواجهة. علاوة على ذلك، فإن فاصل الموثوقية هو حد هامشي تمت معايرته في طبقة البلورات؛ وبسبب تغير المتغيرات (covariate shift)، فإنه يرتب أنظمة المتن حسب الجودة المتوقعة ولكنه لا يوفر حدود دقة معتمدة لأنظمة المتونة الفردية.
قابلية التكرار: يشير المؤلف إلى أن إعادة توليد التنبؤات بنفس البذور لا يؤدي إلى إحداثيات متطابقة بسبب عشوائية النموذج؛ التمثيل المودع هو سحبة واحدة من فضاء التكوين، حيث يعمل عدم الاتفاق عبر الطرق كوسيط للتباين.