L8 · 35 دقيقة
التحقّق: إعادة الالتحام وROC
RMSD وAUC وEF وBEDROC ومجموعات الخداع.
أهداف التعلّم
- تعريف redocking وcross-docking وعتبة RMSD كاختبار هندسي لا كاختبار نشاط.
- استخدام ROC وAUC وEF وBEDROC وdecoys لاختبار الإثراء دون ادّعاء «مثبط جيد».
- تفسير لماذا الدرجة الجيدة لا تستلزم تثبيطاً.
طبقات التحقّق ليست واحدة
التحقق من docking أسئلة منفصلة. (1) هل نُعيد هيئة معقد معروف في بنيته (redocking / self-docking)؟ (2) هل نُعيد ليغند في بنية بلورية أخرى لنفس البروتين (cross-docking) — اختبار مرونة وأمانة الشبكة؟ (3) هل نُثري actives بين decoys (ROC، AUC، EF، BEDROC)؟ (4) هل الترتيب يوافق SAR داخل سلسلة؟ النجاح في (1) شائع نسبياً ولا يستلزم (3) ولا (4). الفشل في (1) يضعف الثقة في فرز لاحق، لكنه قد يعني بروتوناً خطأ لا عجزاً أبدياً للدالة. «جيد في الالتحام» بلا تحديد الطبقة جملة فارغة.
RMSD لإعادة الالتحام
جذر متوسط مربعات انحراف ذرات الليجند (عادة الثقيلة) عن المرجع بعد مواءمة مناسبة.
- Nعدد الذرات الداخلة في الحساب
- \mathbf{r}_{i}موضع الذرة في الهيئة المحسوبة · Å
- \mathbf{r}_{i}^{\mathrm{ref}}موضع الذرة في المرجع (البلورة عادة) · Å
عتبة 2 Å شائعة لنجاح redocking، وهي عرف لا قانون طبيعة. RMSD صغير مع تفاعلات خاطئة (قلب مانح/مستقبل) ليس نجاحاً كيميائياً. RMSD كبير بسبب تناظر رباعي السطوح في tert-butyl قد يكون نجاحاً كيميائياً.
تطبيق في التصميم الدوائي: احسب RMSD بعد مواءمة البروتين لا بعد تحريك الليجند حرّاً في الفضاء، وحدد الذرات (بدون H عادة، وبحذر مع تناظر).
مقاييس الإثراء
- ROC: منحنى معدل الإيجاب الحقيقي مقابل الإيجاب الزائف وأنت تُنزّل عتبة الدرجة.
- AUC: مساحة تحت ROC. 0.5 = عشواء. قريب من 1 إثراء شامل، وقد يخفي فشلاً في أعلى 1% إن كان التحسن في الوسط.
- EF (enrichment factor) عند x%: كم مرة أغنى الجزء العلوي بالـ actives مقارنة بالعشواء. يناسب سؤال الفرز العملي.
- BEDROC: يوزن أعلى القائمة أكثر (Boltzmann-enhanced). مفيد حين تهمّك القمة فقط.
- Decoys: جزيئات يُفترض أنها غير نشطة ومتشابهة فيزيوكيميائياً مع الـ actives حتى لا يُميّز النموذج بالحجم أو logP وحدهما. DUD-E مثال معروف لمجموعات من هذا النوع — استخدمها بوعي لانحيازها هي أيضاً.
| Redocking | Cross-docking | |
|---|---|---|
| المستقبل | بنية المعقد نفسه بعد إزالة الليجند. | بنية أخرى (apo أو ليغند مختلف). |
| ما يُختبر | بحث + دالة في جيب «متعاون». | نقل الهيئة عبر مرونة وحالة بلورية. |
| التفاؤل | مُفرط كمؤشر لفرز حقيقي. | أقسى وأقرب لواقع مكتبة على بنية واحدة. |
لماذا؟ — لماذا decoys المتطابقة فيزيوكيميائياً ألزم من «جزيئات عشوائية من ZINC»؟
إن كان الـ actives أثقل وأشد كراهة للماء من الخلفية، فأي دالة تحب الدفن hydrophobic ستُثري بلا فهم للجيب. المطابقة في MW وlogP وعدد الروابط القابلة للدوران تُجبر التمييز على هندسة التفاعلات. حتى ذلك الحين يبقى انحياز الهيكل: decoys ليست ضمان عدم نشاط تجريبي.
ماذا لو؟ — ماذا لو كان AUC = 0.9 لأن كل الـ actives من سلسلة واحدة وكل decoys بعيدة كيمياء؟
النموذج قد يميّز هيكلاً لا جيباً. هذا تسرب تقييم (analogue bias). اختبر بسلسلة مختلفة، أو بتقسيم زمني/هيكلي. أعلن AUC مع وصف تنوع الـ actives وإلا فالرقم مسرح.
تمرين
redocking RMSD = 1.5 Å، ولا cross-docking، ولا decoys، والمؤلف يستنتج أن الفرز سيكتشف مثبطات نانومولارية. حدّد الطبقة التي نجحت والدعوى الزائدة.
مواضيع مرتبطة