الانحراف التربيعي المتوسط ​​لمواقع الذرات

في المعلوماتية الحيوية ، يُعدّ الانحراف التربيعي المتوسط ​​للجذر لمواقع الذرات ، أو ببساطة الانحراف التربيعي المتوسط ​​للجذر (RMSD) ، مقياسًا لمتوسط ​​المسافة بين ذرات الجزيئات المتراكبة (عادةً ذرات الهيكل الأساسي) . [ 1 ] في دراسة التكوينات البروتينية الكروية، يُقاس التشابه في البنية ثلاثية الأبعاد عادةً باستخدام RMSD لإحداثيات ذرات الكربون ألفا (Cα) بعد تطبيق التراكب الأمثل للأجسام الصلبة.

عندما يتذبذب نظام ديناميكي حول موضع متوسط ​​محدد بدقة، يُشار إلى الانحراف التربيعي المتوسط ​​(RMSD) عن هذا المتوسط ​​بمرور الوقت باسم متوسط ​​التذبذب التربيعي (RMSF) . يمكن قياس حجم هذا التذبذب، على سبيل المثال باستخدام مطيافية موسباور أو الرنين المغناطيسي النووي ، ويمكن أن يوفر معلومات فيزيائية مهمة. يُعد مؤشر ليندمان طريقة لوضع متوسط ​​التذبذب التربيعي (RMSF) في سياق معلمات النظام.

إحدى الطرق الشائعة لمقارنة هياكل الجزيئات الحيوية أو الأجسام الصلبة هي نقل وتدوير أحد الهياكل بالنسبة للآخر لتقليل متوسط ​​الانحراف التربيعي (RMSD). قدم كوتسياس وآخرون اشتقاقًا بسيطًا، قائمًا على الكواترنيونات ، للتحويل الأمثل للجسم الصلب (الدوران-الانتقال) الذي يقلل من متوسط ​​الانحراف التربيعي بين مجموعتين من المتجهات. [ 2 ] أثبتوا أن طريقة الكواترنيونات مكافئة لخوارزمية كابش المعروفة . [ 3 ] يُعد الحل الذي قدمه كابش مثالًا على حل المسألة ذات الأبعاد d ، التي قدمها هيرلي وكاتيل. [ 4 ] نُشر حل الكواترنيونات لحساب الدوران الأمثل في ملحق ورقة بحثية لبيتيتجان. [ 5 ] تم توسيع حل الكواترنيونات هذا وحساب التساوي الأمثل في حالة الأبعاد d ليشمل المجموعات اللانهائية والحالة المستمرة في الملحق أ من ورقة بحثية أخرى لبيتيتجان. [ 6 ]

المعادلة

RمSد=1شمالأنا=1شمالدلتاأنا2{\displaystyle \mathrm {RMSD} ={\sqrt {{\frac {1}{N}}\sum _{i=1}^{N}\delta _{i}^{2}}}}

حيث تمثل δ i المسافة بين الذرة i وبنية مرجعية أو متوسط ​​موضع الذرات N المكافئة. ويتم حساب هذه المسافة عادةً لذرات العمود الفقري الثقيلة C و N و O و C α ، أو أحيانًا لذرات C α فقط .

عادةً ما يتم إجراء عملية تراكب صلبة تقلل من قيمة RMSD، ويتم إرجاع هذه القيمة الدنيا. بافتراض وجود مجموعتين منن{\displaystyle n}نقاطv{\displaystyle \mathbf {v} }وw{\displaystyle \mathbf {w} }، ويُعرَّف الانحراف المعياري للجذر التربيعي (RMSD) على النحو التالي:

RمSد(v،w)=1نأنا=1نvأنا-wأنا2=1نأنا=1ن((vأناx-wأناx)2+(vأناy-wأناy)2+(vأناz-wأناz)2){\displaystyle {\begin{align}\mathrm {RMSD} (\mathbf {v} ,\mathbf {w} )&={\sqrt {{\frac {1}{n}}\sum _{i=1}^{n}\|v_{i}-w_{i}\|^{2}}}\\&={\sqrt {{\frac {1}{n}}\sum _ {i=1}^{n}((v_{ix}-w_{ix})^{2}+(v_{iy}-w_{iy})^{2}+(v_{iz}-w_{iz})^{2}}})\end{محاذاة}}}

تُعبّر قيمة RMSD بوحدات الطول. الوحدة الأكثر استخدامًا في علم الأحياء البنيوي هي الأنجستروم (Å) والتي تساوي 10⁻¹⁰ متر  .

الاستخدامات

يُستخدم عادةً متوسط ​​الانحراف التربيعي الجذري (RMSD) كمقياس كمي للتشابه بين بنيتين بروتينيتين أو أكثر. على سبيل المثال، تستخدم مسابقة CASP للتنبؤ ببنية البروتين متوسط ​​الانحراف التربيعي الجذري (RMSD) كأحد معايير تقييم مدى تطابق البنية المُقدمة مع البنية المستهدفة المعروفة. وبالتالي، كلما انخفضت قيمة متوسط ​​الانحراف التربيعي الجذري (RMSD)، كان النموذج أفضل مقارنةً بالبنية المستهدفة. [ 7 ]

كما يستخدم بعض العلماء الذين يدرسون طي البروتين عن طريق المحاكاة الحاسوبية RMSD كإحداثية تفاعل لتحديد مكان البروتين بين الحالة المطوية والحالة غير المطوية.

تُعدّ دراسة متوسط ​​الانحراف التربيعي الجذري (RMSD) للجزيئات العضوية الصغيرة (والتي تُسمى عادةً بالروابط عند ارتباطها بالجزيئات الكبيرة، مثل البروتينات) شائعة في سياق الالتحام الجزيئي ، [ 1 ] وكذلك في طرق أخرى لدراسة بنية الروابط عند ارتباطها بالجزيئات الكبيرة. تجدر الإشارة إلى أنه في حالة الروابط (على عكس البروتينات، كما ذُكر أعلاه)، لا يتم عادةً تطبيق بنيتها قبل حساب متوسط ​​الانحراف التربيعي الجذري.

يُعدّ RMSD أيضًا أحد المقاييس العديدة التي تم اقتراحها لتحديد التشابه التطوري بين البروتينات، بالإضافة إلى جودة محاذاة التسلسلات. [ 8 ] [ 9 ]

انظر أيضاً

مراجع

  1. 1 2 "الالتحام الجزيئي، وتقدير طاقات الربط الحرة، ومجال القوة شبه التجريبي لبرنامج AutoDock" . موقع سيباستيان راشكا الإلكتروني . 26-06-2014 . تاريخ الاسترجاع: 07-06-2016 .
  2. كوتسياس إي إيه، سيوك سي، ديل كيه إيه (2004). "استخدام الكواترنيونات لحساب متوسط ​​الانحراف التربيعي الجذري". مجلة الكيمياء الحاسوبية . 25 (15): 1849-1857 . doi : 10.1002/jcc.20110 . PMID 15376254. S2CID 18224579 .  
  3. 1 2 كابش دبليو (1976). "حل لأفضل دوران لربط مجموعتين من المتجهات". أكتا كريستالوغرافيكا . 32 (5): 922-923 . Bibcode : 1976AcCrA..32..922K . doi : 10.1107/S0567739476001873 .
  4. هيرلي جيه آر، كاتيل آر بي (1962). "برنامج بروكروستس: إنتاج دوران مباشر لاختبار بنية عاملية مفترضة". العلوم السلوكية . 7 (2): 258-262 . doi : 10.1002/bs.3830070216 .
  5. بيتيجان م (1999). "حول قياسات التناظر الكمي للجذر التربيعي المتوسط ​​والتناظر الكمي" (ملف PDF) . مجلة الفيزياء الرياضية . 40 (9): 4587-4595 . Bibcode : 1999JMP....40.4587P . doi : 10.1063/1.532988 .
  6. بيتيتجان م (2002). "الخلائط الكيرالية" (ملف PDF) . مجلة الفيزياء الرياضية . 43 (8): 185-192 . Bibcode : 2002JMP....43.4147P . doi : 10.1063/1.1484559 .
  7. جامبر جيه، إيفانز آر، بريتزل إيه، وآخرون (2021). "تنبؤ دقيق للغاية ببنية البروتين باستخدام ألفا فولد" . نيتشر . 596 (7873): 583-589 . doi : 10.1038/s41586-021-03819-2 . PMC 8371605 .  
  8. جويت إيه آي، هوانغ سي سي، فيرين تي إي (2003). "MINRMS: خوارزمية فعالة لتحديد تشابه بنية البروتين باستخدام متوسط ​​الجذر التربيعي للمسافة" (ملف PDF) . المعلوماتية الحيوية . 19 (5): 625-634 . doi : 10.1093/bioinformatics/btg035 . PMID 12651721 . {{cite journal}}: صيانة CS1: أسماء متعددة: قائمة المؤلفين ( رابط )
  9. أرموغوم ف، موريتي س، كيدواس ف، نوتردام س (2006). "مؤشر الانحراف المعياري الجذري التربيعي (iRMSD): مقياس محلي لدقة محاذاة التسلسل باستخدام المعلومات الهيكلية" (ملف PDF) . المعلوماتية الحيوية . 22 (14): e35–39. doi : 10.1093/bioinformatics/btl218 . PMID 16873492 . 

للمزيد من القراءة