Multiscale autoregressive modeling generates ordered protein conformational trajectories
تقدم الورقة البحثية نموذج TEMPO، وهو نموذج تلقائي التكرار هرمي يولد بكفاءة مسارات تشكيل بروتينية طويلة ومرتبة مع تقليل تراكم الأخطاء بشكل كبير ودقة عالية تجاه بيانات الديناميكا الجزيئية، متفوقاً بذلك على الأساليب الحالية مثل AlphaFlow.
تُعد البروتينات الآلات العاملة في الحياة، وهي هياكل جزيئية دقيقة تنطوي لتتخذ أشكالًا محددة لأداء مهام مثل حمل الأكسجين، أو محاربة العدوى، أو إرسال الإشارات بين الخلايا. ومع ذلك، فإن هذه الآلات ليست تماثيل ساكنة؛ فهي تتحرك وتتنفس وتغير أشكالها باستمرار لإنجاز المهمة. قد يبدأ البروتين في وضعية مغلقة، ثم ينفتح ليمسك بجزيء، ثم ينغلق بسرعة ليغيره. لقد عرف العلماء منذ زمن طويل أن فهم هذه التحركات لا يقل أهمية عن معرفة الشكل النهائي، لأن المسار الذي يسلكه البروتين للانتقال من حالة إلى أخرى غالبًا ما يحمل المفتاح لكيفية عمله أو كيفية إيقافه بواسطة دواء ما. ولعقود من الزمن، حاول الباحثون محاكاة هذه التحركات باستخدام أجهزة كمبيوتر قوية تحسب فيزياء كل ذرة، لكن هذه المحاكاة بطيئة ومكلفة للغاية، حيث تستغرق غالبًا أسابيع لتوليد بضع ثوانٍ فقط من حركة بروتين واحد. وقد جعل هذا القصور من الصعب دراسة كيفية تحرك البروتينات المختلفة أو مقارنة مسارات العديد من الجزيئات جنبًا إلى جنب.
ولحل هذه المشكلة، طور فريق من الباحثين أداة ذكاء اصطناعي جديدة تسمى TEMPO، يمكنها توليد تسلسلات طويلة ومنظمة من تحركات البروتين في جزء بسيط من الوقت الذي كان يستغرقه الأمر سابقًا. فبدلاً من محاولة حساب كل قوة فيزيائية بين الذرات، يتعلم TEMPO من آلاف المحاكاة الحاسوبية الموجودة للتعرف على الأنماط التي تتطور بها البروتينات طبيعيًا بمرور الوقت. يعمل النظام مثل كاميرا ذات سرعتين: فهو يتنبأ أولاً بالتحولات الكبرى والبطيئة في شكل البروتين، والتي تعمل كنقاط ارتكاز، ثم يملأ الفراغات بالاهتزازات السريعة والأصغر التي تحدث بين تلك النقاط. يسمح هذا النهج للنموذج بإنشاء فيلم مستمر وسلس لدورة حياة البروتين دون الأخطاء التي تتراكم عادةً عندما تحاول أجهزة الكمبيوتر التنبؤ خطوة بخطوة عبر فترات طويلة. وقد اختبر الباحثون هذه الطريقة على مجموعة متنوعة من البروتينات، بما في ذلك الإنزيمات التي تساعد الخلايا على إنتاج الطاقة وعائلة كبيرة من البروتينات تسمى "كينازات" (kinases) التي تنظم العديد من العمليات الخلوية.
تظهر النتائج أن TEMPO يمكنه إعادة إنشاء الخصائص الإحصائية لتحركات البروتين بدقة ملحوية. وعندما قارن الفريق التسلسلات المولدة بالمعيار الذهبي للمحاكاة القائمة على الفيزياء، طابق النموذج الجديد أنماط الحركة المتوقعة بشكل شبه مثالي. فعلى سبيل المثال، تنبأ النموذج بشكل صحيح بمقدار اهتزاز الأجزاء المختلفة من البروتين وكيف سيتغير الشكل العام، محققًا مستوى من الاتفاق يتجاوز بكثير الطرق السابقة. وفي اختبار محدد يتعلق ببروتين يفتح ويغلق مثل صدفة المحار، نجح النموذج في توليد مسار مستمر من الحالة المغلقة إلى الحالة المفتوحة، مظهرًا تقدمًا منطقيًا في تغيرات التلامس يحاكي ما يراه العلماء في التجارب الحقيقية. كما أن سرعة هذه الطريقة الجديدة تمثل طفرة كبرى؛ حيث استغرق توليد تسلسل من 400 إطار أقل من دقيقة واحدة، وهي مهمة كانت ستستغرق ساعات أو أيامًا باستخدام التقنيات القديمة.
وإلى جانب السرعة والدقة، أظهرت الدراسة أن TEMPO يمكن توجيهه لإنتاج أنواع محددة من التحركات. فقد أظهر الباحثون أنه يمكنهم توجيه النموذج لإنتاج أشكال محلية محددة للبروتينات، مثل تكوين معين لنمط كيميائي يُعرف باسم DFG، وهو هدف حاسم للعديد من الأدوية. ومن خلال تطبيق قواعد بسيطة أثناء عملية التوليد، تمكنوا من توجيه النموذج لإنشاء هياكل تبدو مثل حالات "التشغيل" أو "الإيقاف" لهذه البروتينات، ونجحوا في ذلك لعشرات من الكينازات المختلفة. يشير هذا إلى أن النموذج قد تعلم القواعد الأساسية لمرونة البروتين جيدًا بما يكفي لتوجيهه نحو نتائج محددة دون الحاجة إلى إعادة تدريبه لكل بروتين جديد. كما طبق الفريق هذه الأداة على مجموعة من البروتينات تسمى "مستقبلات البروتين G المقترنة" (G protein-coupled receptors)، وهي منغرسة في أغشية الخلايا وتشارك في استشعار البيئة، مما أظهر أن النموذج يمكنه نقل معرفته إلى بروتينات لم يرها من قبل.
تؤكد الدراسة أن البروتينات تتحرك بطريقة هرمية، حيث توفر التغيرات الكبيرة والبطيئة الإطار العام للاهتزازات المحلية الأسرع. ومن خلال الفصل بين هذين المستويين من الحركة، تمكن الباحثون من منع الكمبيوتر من الانحراف عن المسار، وهي مشكلة شائعة حيث تتراكم الأخطاء الصغيرة بمرور الوقت لتخلق شكلًا نهائيًا غير منطقي. وبينما يركز الإصدار الحالي من الأداة على الهيكل الأساسي للبروتين ولا يشمل بعد الماء أو الأيونات أو الروابط الكيميائية المحيطة، إلا أنه يوفر إطارًا قويًا جديدًا لمقارنة كيفية تحرك البروتينات المختلفة. وتفتح هذه القدرة الباب أمام الفحص السريع لكيفية تأثير الطفرات أو الأدوية على السلوك الديناميكي للبروتينات، مما يوفر عدسة جديدة يمكن من خلالها رؤية الآلات الجزيئية للحياة. إن هذا العمل لا يحل محل الحاجة إلى عمليات محاكاة فيزيائية مفصلة، بل يقدم مكملاً سريعًا وموثوقًا يسمح للعلماء باستكشاف المشهد الشاسع لحركة البروتين بوضوح وسرعة كانا بعيدي المنال في السابق.
تعتمد وظيفة البروتين على الانتقالات بين الحالات التشكيلية، مما يستلزم فهماً للتسلسلات الهيكلية المرتبة بدلاً من مجرد النقاط النهائية الساكنة أو المجموعات غير المنظمة. وبينما توفر محاكاة الديناميكا الجزيئية (MD) مسارات دقيقة ذرية، إلا أنها مكلفة حاسوبياً، مما يحد من المقارنات المنهجية عبر مجموعات كبيرة من البروتينات، أو الطفرات، أو حالات الربيطة (ligand). وتتفوق النماذج التوليدية العميقة الحالية (مثل AlphaFlow وBioEmu) في أخذ عينات من المجموعات التشكيلية المتنوعة، لكنها غالباً ما تفشل في الحفاظ على التبعيات الزمنية المطلوبة للمسارات التسلسلية ذات المعنى. ومن ناحية أخرى، تواجه المولدات التسلسلية تحدي تراكم الأخطاء أثناء عمليات التوليد الطويلة (rollouts)، حيث تتراكم الأخطاء الهندسية مع كل تنبؤ يصبح هو السياق للتنبؤ التالي. علاوة على ذلك، فإن حركات البروتين هي بطبيعتها متعددة المقاييس، حيث تتضمن تقلبات محلية سريعة مدمجة ضمن إعادة ترتيب جماعية أبطأ؛ ويواجه آلية التحديث أحادية المقياس صعوبة في الموازنة بين التباين قصير المدى والتطور مستقر طويل المدى.
المنهجية: TEMPO
قدم المؤلفون TEMPO، وهو نموذج تكراري هرمي تم تدريبه مباشرة على مسارات الديناميكا الجزيئية (MD) لتوليد تشكيلات هيكل العمود الفقري للبروتين المرتبة عند مستويين متميزين من الدقة الزمنية.
البنية والتمثيل
تمثيل العمود الفقري: يتم تمثيل كل ثمالة (residue) بإطار عمود فقري SE(3) (دوران عبر كواتيرنيون وحدة، وإزاحة) وزوايا الالتواء (ϕ,ψ,ω) المشفرة كأزواج جيب-جيب تمام (sine-cosine) لتجنب الانقطاعات.
البنية العصبية: يستخدم النموذج بنية مشتركة لكلا المقياسين، حيث يوظف انتباه النقاط الثابتة (Invariant Point Attention - IPA) للترميز المكاني لأطر SE(3) والانتباه الذاتي متعدد الرؤوس للتفاعلات بين الثمالات. ويتم نمذجة التبعيات الزمنية باستخدام طبقات الوحدة المتكررة ذات البوابة (GRU).
صياغة الديناميكيات: يتم تحديد التطور التشكيلي باستخدام تحديث منفصل مستوحى من المعادلات التفاضلية العشوائية (SDE): Xt+Δt=Xt+fθ(ctxt)Δt+σΔtϵt. التنبؤ مشروط بسياق مكون من إطارين، مما يجعل العملية غير ماركوفية من الدرجة الأولى.
الإطار الهرمي
يقوم TEMPO بتفكيك عملية التوليد إلى مقياسين للحد من الانجراف طويل الأمد:
نموذج المقياس البطيء: يقوم بتقديم تشكيلات "المرساة" (anchor) عند خطوة أخذ عينات قدرها 20 نانو ثانية (مطابقة لخطوة بيانات التدريب). تنظم هذه المراسي التطور طويل المدى للمسار.
نموذج المقياس السريع: يولد الإطارات المتوسطة عند خطوة قدرها 1 نانو ثانية. هذا النموذج مشروط بسياق المقياس البطيء (تشكيلات المرساة) لحل التباين الهيكلي بين المراسي.
عملية الاستدلال: أثناء الاستدلال، يولد نموذج المقياس البطيء تسلسل المراسي أولاً. بعد ذلك، يولد نموذج المقياس السريع الـ 19 إطاراً متوسطاً لكل مقطع قدره 20 نانو ثانية، مما يشكل مساراً مستمراً ومرتباً.
التدريب والاشتراط
بيانات التدريب: تم التدريب على 320,000 مسار من 4,000 بروتين عشوائي من مجموعة mdCATH.
دالة الخسارة: تجمع بين إعادة بناء متوسط مربع الخطأ (MSE) وعقوبة التصادم الفراغي (عقوبة ReLU للذرات غير المتجاورة التي تبعد أقل من 1.2 أنجستروم).
آليات الاشتراط:
المشروط بالنقاط النهائية: يقوم مشفر قابل للتعلم برسم خريطة لهيكل مستهدف (X∗) إلى تضمين لكل ثمالة، والذي يُضاف إلى الحالة الكامنة لتوجيه عملية التوليد نحو نقطة نهاية محددة.
التوجيه أثناء الاستدلال: للتحكم في الزخارف المحلية (مثل هندسة DFG)، تقوم عملية توجيه بتنقية ثمالات محددة بعد كل خطوة تكرارية باستخدام هدف متعدد الحدود دون تحديث معاملات النموذج.
النتائج الرئيسية
دقة وكفاءة المجموعة
على بروتينات mdCATH وATLAS المحجوبة، أعاد TEMPO إنتاج ملاحظات المجموعة المستمدة من الديناميكا الجزيئية:
الارتباطات: حقق ارتباطات بيرسون بلغت 0.91 لمتوسط جذر متوسط المربعات (RMSD) المزدوج للمجموعة، و0.89 لمتوسط جذر متوسط المربعات (RMSF) على مستوى الثمالة في ATLAS.
تحليل المكونات الرئيسية (PCA): اتفقت اتجاهات المكون الرئيسي الرائدة مع الديناميكا الجزيئية في 76% من حالات ATLAS، متفوقاً بشكل كبير على AlphaFlow (44%).
الكفاءة: تطلب توليد 400 إطار أقل من دقيقة واحدة، أي أسرع بحوالي 750 مرة من AlphaFlow تحت الإعدادات المعلنة.
تقليل الانجراف: مقارنة بنموذج تكراري أحادي المقياس، قلل التسلسل الهرمي متعدد المقاييس من خطأ RMSD طويل الأمد من 8.62 أنجستروم إلى 1.78 أنجستروم لعملية توليد مكونة من 400 إطار.
التطبيق 1: انفتاح أدينيلات كيناز (ADK)
في 10 عمليات توليد مشروطة بالنقاط النهائية من الحالة المغلقة (1AKE) إلى الحالة المفتوحة (4AKE):
شكلت التشكيلات المولدة ممرًا مستمرًا في مساحة RMSD.
حدد التنسيق الزمني (tICA) أن الإحداثي المشترك من المغلق إلى المفتوح (tIC1) أظهر تقدماً متسقاً مع ترتيب الإطارات.
كشف التحليل عن تبادل اتصالات تدريجي: تناقصت الاتصالات الخاصة بالحالة المغلقة بينما تراكمت الاتصالات الخاصة بالحالة المفتوحة، مع تمركز الحركة في نطاقي NMP وLID.
التطبيق 2: توجيه DFG في الكيناز
تم تطبيق التوجيه أثناء الاستدلال على 100 كيناز لتوليد هندسات بديلة لعمود فقري DFG المحلي:
بالنسبة لـ 68 كيناز، نجحت الطريقة في إنتاج إطارات تشبه حالة DFG-in وحالة DG-out.
شغلت الإطارات المقبولة مناطق متميزة في مخطط راماشاندران لـ DFG-Phe، مما أظهر القدرة على التحكم في الهندسة المحلية للعمود الفقري دون إعادة تدريب.
التطبيق التطبيق 3: مرونة الكيناز وانتقال GPCR
مرونة الكيناز: عبر 252 كيناز، أظهرت ملفات التقلب المولدة توافقاً متغيراً مع عوامل B البلورية (متوسط سبيرمان ρ=0.358؛ و22% مع ρ≥0.5).
انتقال GPCR: تم تطبيق النموذج، المدرب فقط على البروتينات القابلة للذوبان (mdCATH)، على ستة مستقبلات GPCR (الفئات A–C) دون تحديثات محددة. تداخلت التشكيلات المولدة مع المناطق البارزة لإسقاطات الديناميكا الجزيئية الخاصة بالمستقبلات فقط، مع تركز التباين في المناطق المرنة (النهاية N، والحلقات) بينما ظل قلب عبر الغشاء مستقراً.
الأهمية والادعاءات
يدعي البحث أن TEMPO يوفر إطاراً لتوليد ومقارنة تسلسلات تشكيلية بروتينية مرتبة، معالِجاً التحدي المحدد المتمثل في تراكم الأخطاء في التوليد طويل المدى. تكمن أهميته الأساسية في:
الاستقرار الهرمي: إثبات أن فصل التوليد إلى مقاييس بطيئة (جماعية) وسريعة (محلية) يقلل بشكل كبير من الانجراف مع الحفاظ على الإحصائيات التشكيلية الشبيهة بالديناميكا الجزيئية.
القابلية للتحكم: إظهار أن النموذج يدعم التوليد المشروط بالنقاط النهائية (لمسارات الانتقال) والتوجيه أثناء الاستدلال (للتلاعب بالزخارف المحلية).
الكفاءة: تقديم بديل سريع للمحاكاة القائمة على الفيزياء لمقارنة التغيرات التشكيلية المرتبة عبر مجموعات كبيرة من البروتينات.
يشير المؤلفون صراحةً إلى القيود: النموذج يركز على العمود الفقري (مما يستبعد السلاسل الجانبية، والربيطات، والمذيب)، وتسميات الوقت تعكس خطوات أخذ العينات وليس وقتاً حركياً معايراً، كما أن عمليات التوليد المشروطة أو الموجهة ليست عينات غير منحازة. يضع هذا العمل TEMPO كمكمل، وليس بديلاً، للمحاكاة القائمة على الفيزياء والتحقق التجريبي.