L9 · 28 دقيقة
RMSD — المعادلة وسوء الاستخدام
مسطّح لا يعني متقارباً.
أهداف التعلّم
- كتابة معادلة RMSD وشرح المواءمة والذرات الداخلة.
- رفض تفسير RMSD المسطح كتقارب أو كثبات ارتباط.
ماذا يقيس RMSD؟
RMSD مسافة بين هيئتين بعد إزالة انزياح ودوران جسميّين عادة (مواءمة). في MD يُقاس غالباً مقابل الإطار الأول أو مقابل البلورة، لذرات العمود الفقري أو لليجند. هو مقياس هندسي واحد: لا طاقة، لا إنتروبي، لا احتمال وادٍ. مسار يقفز بين هيئتين متباعدتين قد يُظهر RMSD «مزدوجاً»؛ مسار يزحف ببطء في اتجاه واحد قد يبدو مسطحاً في نافذة قصيرة. التسطيح يعني أن المسافة إلى المرجع لم تتغير كثيراً في تلك النافذة — لا أن المعاينة اكتملت، ولا أن الليجند «مرتبط بقوة».
RMSD مقابل مرجع
جذر متوسط مربعات انحراف الذرات المختارة بعد المواءمة إن طُبّقت.
- Nعدد الذرات الداخلة (مثلاً Cα أو ليغند ثقيل)
- \mathbf{r}_{i}(t)موضع الذرة i في الإطار t · nm أو Å
- \mathbf{r}_{i}^{\mathrm{ref}}الموضع في الهيئة المرجعية · nm أو Å
اختيار N يغيّر القصة: Cα يُخفي حلقة ليغند؛ ليغند وحده يُخفي تفكك نطاق بروتيني. المواءمة على البروتين ثم RMSD لليجند تقيس انزياح الليجند في الجيب — وهذا أقرب لسؤال الالتحام من RMSD لكل النظام بما فيه الماء.
تطبيق في التصميم الدوائي: حدد المرجع والذرات والمواءمة في كل رسم. بلا ذلك الرسم ليس معلومة.
gmx rms -s md.tpr -f md.xtc -o rmsd.xvg
الغرض: حساب RMSD من مسار مقابل بنية في tpr أو مرجع تُحدده تفاعلياً.
المُدخل: md.tpr (طوبولوجيا/مرجع) وmd.xtc (مسار).
المُخرج: rmsd.xvg زمن مقابل RMSD. اختيار المجموعات تفاعلي (Backbone، ligand، إلخ).
المعنى: رقم كل إطار بعد المواءمة الافتراضية للمجموعة. تفسيره العلمي ليس داخل الأمر.
أخطاء شائعة: مواءمة على ماء، أو مجموعة ليغند غير معرّفة في الفهرس، أو مسار غير مصحح لقفزات PBC قبل التحليل البصري (يُعالج بأدوات تحويل المسار حسب الحاجة).
لماذا؟ — لماذا المرجع = الإطار الأول قرار غير بريء؟
الإطار الأول قد يكون ما بعد الاتزان مباشرة — لا البلورة ولا الهيئة النشطة. RMSD المنخفض يعني «لم نبتعد عن انطلاقنا» وهذا نجاح بروتوكول أحياناً وعمى أحياناً. المرجع البلوري يقيس الانحراف عن التجربة. المرجعان مشروعان إن سُمّيا.
ماذا لو؟ — ماذا لو نشرت RMSD لليجند دون مواءمة البروتين، فظهر قفز كل بضعة نانوثوانٍ؟
قد يكون قفز PBC أو دوران جسم صلب للمعقد كله، لا تفككاً. صحّح الحدود الدورية وواءم على الجيب ثم أعد الرسم قبل أن تكتب «unstable binding».
تمرين
RMSD Cα مسطح عند 0.25 nm لـ 200 ns، وRMSD الليجند بعد مواءمة الجيب يرتفع فجأة إلى 0.8 nm في النسخة 2 فقط. ماذا تستنتج وماذا لا تستنتج؟