Pose: a pure-NumPy, chirality-aware, structural engine for building and manipulating the molecular structures of proteins, nucleic acids, and small organic molecules
تُعد Pose مكتبة بايثون خفيفة الوزن تعتمد كلياً على NumPy، وتوفر إطار عمل موحداً ومدركاً للتماثل المرآتي لبناء ومعالجة هياكل البروتينات، والأحماض النووية، والجزيئات الصغيرة، مع تنفيذ دقيق لعدة حقول قوة ودوال تسجيل بأداء يضاهي الأدوات الراسخة مثل OpenMM وPyRosetta وAutoDock Vina.
غالباً ما يواجه العلماء الذين يدرسون شكل وسلوك اللبنات الأساسية للحياة واقعاً محبطاً: فمن أجل فهم بروتين واحد، يتعين عليهم التعامل مع عشرات البرامج المختلفة. كل برنامج يتحدث لغته الخاصة، ويستخدم تنسيقات ملفاته الخاصة، ويتبع قواعده الخاصة. قد يقوم الباحث ببناء جزيء في أداة ما، ثم يحاكي حركته في أداة أخرى، ثم يقيم استقراره في أداة ثالثة. وفي كل مرة ينتقل فيها بين هذه الأدوات، يجب ترجمة البنية، وفي هذه الترجمة، يمكن أن تتسلل أخطاء دقيقة. والأسوأ من ذلك، أن معظم هذه الأدوات تجد صعوبة في التعامل مع الجزيئات التي تعد صوراً مرآتية للنسخ الطبيعية. ففي الطبيعة، تُبنى البروتينات من أحماض أمينية يسارية اليد، لكن العلماء يحتاجون أحياناً إلى إنشاء نسخ يمينية لتصنيع أدوية تدوم لفترة أطول في الجسم. وغالباً ما تتعامل البرمجيات الحالية مع هذه الصور المرآتية كأجسام غريبة أو تعطلها تماماً، مما يضطر الباحثين إلى استخدام حلول بديلة تثير الشك وعدم اليقين.
هنا يأتي دور "بوز" (Pose)، وهي أداة جديدة صُممت لدمج هذه العملية المجزأة في مساحة واحدة شفافة. "بوز"، التي ابتكرتها ساري سابان من جامعة الملك عبد العزيز، هي مكتبة برمجية يمكنها بناء وتحريك وتحليل البروتينات والحمض النووي (DNA) والجزيئات الصغيرة الشبيهة بالأدوية دون الحاجة إلى أي محركات خارجية معقدة. لقد بُنيت بالكامل على مجموعة أدوات رياضية قياسية تسمى "نمباي" (NumPy)، مما يعني أنها تعمل على أي جهاز كمبيوتر تقريباً دون متطلبات التثبيت الثقيلة لبرمجيات النمذجة الجزيئية التقليدية. الميزة الأبرز لـ "بوز" هي قدرتها على التعامل مع الجزيئات اليسارية واليمينية بتماثل تام. ففي العالم الطبيعي، يكون الحمض الأميني يساري اليد وصورته المرآتية يمينية اليد متطابقين كيميائياً في كل شيء باستثناء اتجاههما في الفراغ. ويحترم "بوز" هذه القاعدة من خلال تطبيق مجموعة واحدة ومتسقة من الجداول الرياضية لكليهما، مما يضمن أن تحصل الجزيئة المرآتية على نفس درجة الطاقة تماماً مثل الأصلية. وتسمح هذه الاتساقية للباحثين بخلط أجزاء يسارية ويمينية في سلسلة واحدة دون أن ترتبك البرمجية أو تعطي نتائج متناقضة.
اختبر الباحثون "بوز" من خلال مقارنة حساباتها مع خمسة برامج قياسية ومعتمدة في الصناعة. فقد أدخلوا نفس الهياكل الجزيئية في "بوز" والبرامج المرجعية ليروا مدى تطابق الأرقام. وعند استخدام إحداثيات متطابقة، أعاد "بوز" إنتاج حسابات الطاقة للبرامج المرجعية بدقة متناهية. وبالنسبة لمجالات القوة المستخدمة لمحاكاة الحركة الفيزيائية، كان الفرق غالباً أقل من جزء من مليون من الكيلوجول لكل مول، وهو هامش ضئيل جداً يُرجح أنه ناتج عن طريقة معالجة أجهزة الكمبيوتر للأرقام العشرية وليس خللاً في المنطق. وعند اختبار وظائف التقييم المستخدمة للتنبؤ بمدى قدرة الدواء على الارتباط بالبروتين، طابق "بوز" نتائج الأدوات المرجعية في حدود بضعة أجزاء من الألف من السعرة الحرارية لكل مول. والأهم من ذلك، عندما قام الباحثون بعكس هيكل ما لإنشاء نسخة يمينية اليد، ظلت درجات الطاقة في "بوز" دون تغيير، مما أكد أن البرنامج يعامل الصورة المرآتية تماماً كما تفعل الطبيعة. وفي المقابل، أظهر أحد البرامج المرجعية، وهو "بيروسيتتا" (PyRosetta)، اختلافات كبيرة في درجاته للهياكل المعكوسة، مما كشف عن عدم تماثل في كيفية تعامله مع البيانات.
وبعيداً عن مجرد حساب الطاقة، تم اختبار "بوز" في قدرته على إعادة بناء هياكل معقدة من الصفر. أخذ الباحثون الزوايا بين الذرات في بروتينات حقيقية وطلبوا من "بوز" إعادة بناء الشكل الكامل من تلك الزوايا وحدها. وكانت النتيجة هيكلاً يختلف عن الأصلي بأقل من جزء من مئة من الأنجستروم، وهي مسافة أصغر من عرض ذرة واحدة. وقد ظل هذا المستوى من الدقة ثابتاً حتى بالنسبة لسلسلة قصيرة من الحمض النووي تحتوي على أحماض أمينية يمينية اليد، وهي مهمة لا تستطيع العديد من الأدوات الأخرى القيام بها بموثوقية. كما أظهرت البرمجية قدرتها على أداء مهام ديناميكية، مثل محاكاة حركة الذرات بمرور الوقت أو إيجاد الترتيب الأكثر استقراراً للسلاسل الجانبية. وبينما استغرقت هذه المحاكاة وقتاً أطول للتشغيل مقارنة بالبرامج المرجعية المتخصصة وعالية السرعة، إلا أنها أنتجت نتائج لا يمكن تمييزها إحصائياً في منتصف النطاق. فعلى سبيل المثال، عند ربط الجزيئات الدوائية في جيوب البروتين، نجح "بوز" في إيجاد الموضع الصحيح في ست حالات من أصل عشر تجارب، محققاً أداءً يماثل الأدوات الراسخة المستخدمة على نطاق واسع في اكتشاف الأدوية.
لا يدعي هذا العمل أنه الطريقة الأسرع لمحاكاة الجزيئات؛ بل يقر المؤلفون بأن نهج البرمجيات الصرفة الخاص بهم أبطأ بكثير من المحركات المجمعة (compiled engines) عالية التحسين التي تستخدمها شركات الأدوية الكبرى. بدلاً من ذلك، تكمن قيمة "بوز" في وضوحه ومرونته. ولأن كل عملية حسابية مكتوبة كسطر برمجي مقروء بدلاً من كونها مخفية داخل برنامج "صندوق أسود"، يمكن للباحث رؤية كيفية اشتقاق درجة الطاقة بدقة. ويمكنه تعديل حد واحد، أو استبدال قاعدة، أو إضافة نوع جديد من الجزيئات دون الحاجة إلى إعادة تجميع مجموعة برمجية ضخمة. وهذا الوضوح يجعل منه أداة مثالية لتعليم الطلاب كيفية عمل القوى الجزيئية، أو لتطوير وظائف طاقة جديدة، أو لتصميم جزيئات مبتكرة تمزج بين الأجزاء اليسارية واليمينية. ومن خلال توحيد التعامل مع مختلف مجالات القوة وضمان التعامل مع الصور المرآتية بأمانة رياضية، يوفر "بوز" أساساً نظيفاً وموثوقاً للجيل القادم من التصميم الجزيئي.
ملخص تقني لـ "Pose: محرك بنيوي يعتمد على NumPy الصرف ويدعم مراعاة الكيرالية (التماثل المرآتي)"
بيان المشكلة تتطلب بروتوكولات النمذجة الجزيئية بشكل متكرر عبور الحدود بين أنظمة برمجية متباينة (على سبيل المثال، بناء الهياكل في أداة ما، والمحاكاة باستخدام OpenMM، والتقييم باستخدام Rosetta، والالتحام باستخدام AutoDock Vina). يعمل كل برنامج بآليات، وتنسيقات ملفات، ومتطلبات تثبيت، واصطلاحات داخلية متميزة. توجد فجوة حرجة في المعالجة المتسقة للكيرالية (Chirality)؛ فبينما تعد الأحماض الأمينية من النوع D والأحماض غير القياسية ضرورية لتصميم الأدوية ومقاومة التحلل البروتيني، تعاني العديد من الأدوات في التعامل معها بشكل أصلي. وتحديداً، فإن دالات الطاقة التي تمت معايرتها لبقايا النوع L (مثل شبكات CMAP أو مكتبات الروتامير) غالباً ما تعطي طاقات غير صحيحة للصور المرآتية ما لم يتم عكس زوايا الالتواء صراحةً، وهي عملية لا يتم تنفيذها بشكل موحد عبر المحركات المختلفة. علاوة على ذلك، فإن إعادة تنفيذ هذه الدوال أمر صعب التحقق منه دون مقارنة شفافة، بنداً ببند، مقابل البرامج المرجعية.
المنهجية يقدم المؤلفون Pose، وهو مكتبة بايثون مصممة كمحرك بنيوي يعتمد كلياً على NumPy دون أي امتدادات مجمعة (compiled extensions). تعتمد المكتبة في تشغيلها حصرياً على NumPy، وتمنح الأولوية للشفافية والقابلية للتوسع على حساب السرعة الحسابية الخام.
البنية والتصميم: يمثل Pose الحالات الجزيئية (الببتيدات المتعددة، الأحماض النووية، والجزيئات الصغيرة) كقواميس تحتوي على الإحداثيات الكارتيزية، وسجلات الذرات، ورسوم الروابط البيانية، ونماذج البقايا (templates). تم تنفيذ جميع حدود الطاقة كعمليات مصفوفات NumPy قابلة للقراءة، مما يسمح للمستخدمين بفحص أو تعديل أو إعادة وزن الحدود الفردية مباشرة في بايثون.
اصطلاح الكيرالية: يطبق Pose اصطلاح كيرالية موحداً حيث يؤدي استخدام رموز البقايا ذات الأحرف الصغيرة في التسلسل إلى بناء بقايا من النوع D. ويتم تحقيق ذلك عن طريق عكس نموذج النوع L عبر المستوى xy. ومن الأهمية بمكان أنه بالنسبة لبقايا النوع D، يقرأ Pose جميع الجداول المعتمدة على اليدوية (مكتبات الروتامير، خرائط راماشاندران، شبكات CMAP) عند زوايا التواء معكوسة. وهذا يضمن أن تحصل البنية وصورتها المرآتية على طاقات متطابقة، وهي خاصية تحقق منها المؤلفون رياضياً وتجريبياً.
حقول القوة ودالات التقييم: أعاد Pose تنفيذ ثلاثة حقول قوة (OpenFF Sage 2.3.0، وAMBER ff19SB، وCHARMM36) ودالتين للتقييم (Rosetta REF15، وAutoDock Vina).
حقول القوة: تتضمن الروابط التوافقية والزوايا، حدود Urey-Bradley، التواءات فورير، لندارد-جونز، الكهروسكونية كولوم، وتصحيحات CMAP ثنائية التكعيب (bicubic). يتم اشتقاق الشحنات من Sage، أو مجموعات AM1-BCC مسبقة الحساب، أو نسخة NumPy معاد تنفيذها لنموذج OpenFF-GNN.
دالات التقييم: تم تنفيذ REF15 باستخدام 19 حداً موزوناً، بالاعتماد على مكتبة روتامير تعتمد على الهيكل العظمي لعام 2010. كما تم تنفيذ AutoDock Vina باستخدام حدودها الخمسة المتبادلة.
بروتوكولات أخذ العينات (Sampling): تتضمن المكتبة طرقاً للتقليل (minimization) (مثل المنحدر المتدرج، L-BFGS، FIRE)، والديناميكا الجزيئية (Verlet، BAOAB Langevin)، والتلدين المحاكى (simulated annealing)، وتعبئة السلسلة الجانبية، والالتحام (docking).
البقايا غير القياسية: يدعم Pose 26 حمضاً أمينياً، بما في ذلك ستة أحماض أمينية غير قياسية (مثل سيلينوميثيونين، وفوسفوتيروسين). يتم اشتقاق النماذج ومكتبات الروتامير لهذه الأحماض من أصولها القياسية أو باستخدام بروتوكولات Rosetta.
المساهمات الرئيسية
معالجة موحدة للكيرالية: يوفر Pose إطاراً متسقاً لبناء وتقييم السلاسل مختلطة الكيرالية (L/D) والبقايا غير القياسية، مما يضمن التماثل المرآتي في حسابات الطاقة حيث تفشل البرامج المرجعية غالباً.
إعادة تنفيذ شفافة: من خلال استخدام NumPy الصرف، يتيح Pose للباحثين تتبع كل عملية حسابية للطاقة من معادلتها الرياضية وصولاً إلى الكود، مما يسهل تطوير المنهجيات والتعليم.
التكامل عبر البروتوكولات: تمكن المكتبة البروتوكولات التي تمتد عبر حقول قوة ودالات تقييم متعددة دون الحاجة لتحويلات ملفات وسيطة، مما يقلل من الأخطاء المتعلقة بالتنسيقات.
مجموعة التحقق: يقدم المؤلفون مجموعة اختبار شاملة تقارن Pose بكل من OpenMM وPyRosetta وAutoDock Vina.
النتائج تحقق المؤلفون من صحة Pose من خلال أربعة اختبارات مرجعية:
الاختبار المرجعي 1 (الاتفاق في الطاقة): على إحداثيات متطابقة، طابق Pose البرامج المرجعية بدقة عالية:
طابق OpenFF Sage 2.3.0 برنامج OpenMM/OpenFF Interchange بدقة 8.8×10−5 كيلوجول/مول.
طابق ff19SB وCHARMM36 برنامج OpenMM بدقة تصل إلى 6.5×10−4 كيلوجول/مول.
طابق AutoDock Vina البرنامج بدقة 0.013 كيلو كالوري/مول.
طابق Rosetta REF15 برنامج PyRosetta بدقة 2×10−4 REU في 25 من أصل 30 بنية. وعُزيت التناقضات إلى بقع برولين معينة عند الطرف النتروجيني (N-terminal) في Rosetta والتي لا يحاكيها Pose بشكل كامل.
الاختبار المرجحي 2 (تماثل الكيرالية): أدى الانعكاس المرآتي للهياكل المكونة من النوع L بالكامل والنوع D بالكامل في Pose إلى تغييرات في الطاقة لا تتجاوز 4×10−12 كيلوجول/مول (ضجيج النقطة العائمة)، وظلت الدرجات دون تغيير. وفي المقابل، تغيرت درجات REF15 في PyRosetta للمجمعات المرآتية بمقدار يصل إلى 31.4 REU بسبب عدم التماثل في كيفية تقييم الجلايسين وغيره من البقايا في السلاسل المرآتية.
الاختبار المرجعي 3 (إعادة بناء الهيكل): أعاد Pose بناء البروتينات من زوايا الالتواء بجذر متوسط مربع (Cα RMSD) قدره 0.0097 أنغستروم مقارنة بـ PyRosetta. ونجح في إعادة بناء جزء من غراميسيدين A يحتوي على ست بقايا من النوع D بجذر متوسط مربع قدره 0.455 أنغستروم.
الاختبار المرجعي 4 (بروتوكولات أخذ العينات): أنتجت طرق النمذجة في Pose (التعبئة، التلدين، التقليل، الالتحام) نتائج مماثلة للمراجع في الأداء الوسيط، وإن كانت بمعدلات نجاح أقل بكثير في الالتحام (7/10 مقابل 10/10 للبروتين-بروتين؛ 6/10 مقابل 8/10 للمستقبل-الربيطة) وأوقات تنفيذ أبطأ.
الأهمية والادعاءات يزعم البحث أن Pose يثبت بنجاح أن مكتبة تعتمد على NumPy الصرف يمكنها إعادة تنفيذ ميكانيكا جزيئية معقدة ودالات تقييم بدقة مع الحفاظ على تماثل مرآتي صارم للأنظمة الكيرالية. ويؤكد المؤلفون أن Pose مصمم للشفافية والقابلية للتوسع وليس للسرعة عالية الإنتاجية. وهو يعمل كأداة للباحثين الذين يطورون حدود طاقة جديدة، وللمعلمين الذين يشرحون آليات الحساب، وللمستخدمين الذين يعملون مع أنظمة مختلطة الكيرالية حيث قد تفتقر الحزم المجمعة القياسية إلى الدعم أو الاتساق.
يشير المؤلفون بتواضع إلى أنه بينما يطابق Pose البرامج المرجعية عند وجود إحداثيات متطابقة، فإن الاختبارات المرجعية هي اختبارات اتفاق وليست دقة تنبؤية مقابل البيانات التجريبية. كما يقرون بأن Pose أبطأ من المحركات المجمعة (بمعامل يتراوح بين 2 إلى 2,000 مرة)، وأن بعض بروتوكولات أخذ العينات (مثل الطي من الصفر) لم تصل إلى الأشكال الأصلية، وهو ما يتفق مع قيود الطرق المستخدمة عند تنفيذها في بيئة غير محسنة. المساهمة الأساسية هي توفير بيئة واحدة قابلة للفحص للبروتوكولات التي تعبر حقول القوة، ودالات التقييم، واصطلاحات الكيرالية.