Deep Mechanistic Models reveal pathway-extrinsic drivers of mammary MAPK signalling heterogeneity
تقدم هذه الورقة النماذج الميكانيكية العميقة (DMMs)، وهي إطار عمل يدمج تعلم التمثيل شبه المُشرف مع نماذج المعادلات التفاضلية العادية للكشف عن أن عدم تجانس إشارات MAPK في الثدي مدفوع أساساً بعوامل خارجية المسار بدلاً من التباينات الجوهرية في المسار، بينما يستخدم في الوقت ذاته التباينات في النماذج لتحديد آليات بيولوجية جديدة.
تخيل جسدك كمدينة صاخبة حيث كل خلية هي مكتب صغير ومزدحم. تحتاج هذه المكاتب إلى التواصل مع بعضها البعض لتقرير ما يجب فعله: النمو، أو الانقسام، أو الراحة. وهي تفعل ذلك عن طريق إرسال رسائل عبر "مسارات الإشارات"، والتي تشبه شبكات معقدة من الأسلاك والمفاتيح داخل الخلية. أحد أهم هذه الشبكات هو مسار MAPK، وهو سلسلة قيادة تبدأ عندما تصطدم إشارة نمو (مثل رسالة من العمدة) بالباب الأمامي للخلية. عندما يفتح الباب، يركض سلسلة من الرسل عبر الخط، ويمررون الملحوظة حتى تصل إلى النواة، وهي "الرئيس" في الخلية، لتخبره ببدء العمل.
ولكن هنا يكمل اللغز: حتى لو كان لديك مجموعة من الخلايا التي هي توائم متطابقة جينياً، فإنها لا تتفاعل جميعها بنفس الطريقة مع نفس الرسالة. فبعضها يصبح متحمساً للغاية ويبدأ في النمو بسرعة؛ بينما يظل البعض الآخر هادئاً. لقد حاول العلماء فهم السبب باستخدام أداتين مختلفتين لحل هذا اللغز. إحدى الأداتين تشبه مخططاً تفصيلياً للأسلاك (النماذج الميكانيكية)، وهي رائعة لفهم كيفية عمل المفاتيح ولكنها سيئة في تفسير سبب كون مكتب ما أكثر فوضوية من غيره. أما الأداة الأخرى فهي تشبه كاميرا عالية التقنية تلتقط لقطة لكل شيء في المكتب وتجد الأنماط (تعلم الآلة)، لكنها لا تخبرك كيف تؤثر الفوضى على العمل. ولفترة طويلة، لم يستطع العلماء دمج هذين المنظورين بسهل للحصول على الصورة الكاملة.
هنا يأتي دور فريق من الباحثين بابتكار ذكي يسمى "النماذج الميكانيكية العميقة" (DMMs). فكر في هذا كـ "مترجم فائق الذكاء" يتحدث لغتي "المخطط" و"اللقطة". لقد بنوا نموذجاً حاسوبياً يأخذ كمية هائلة من البيانات من 63 خطاً خلوياً مختلفاً لسرطان الثدي — مثل التقاط صورة لكل عنصر في 63 مكتباً مختلفاً — ويستخدمها لتعلم القواعد الخفية لمسار MAPK. وبدلاً من مجرد التخمين، يتعلم النموذج من خلال محاولة التنبؤ بكيفية استجابة الخلايا لمختلف الأدوية وإشارات النمو، مع تعديل "مقابض" التحكم الداخلية لديه حتى تتطابق تنبؤاته مع الواقع.
كانت المفاجأة الكبيرة التي وجدوها هي أن الاختلافات بين هذه الخلايا لا تنتج غالسباً عن كون الأسلاك داخل مسار MAPK معطلة أو مختلفة. بدلاً من ذلك، يشير النموذج إلى أن المحركات الحقيقية هي عوامل "خارج المسار" (pathway-extrinsic) — وهي أشياء تحدث خارج سلسلة القيادة الرئيسية. الأمر يشبه إدراك أن مكتبين يتفاعلان بشكل مختلف ليس لأن أسلاكهما الداخلية مختلفة، بل لأن أحدهما لديه مدير نشط وصاخب جداً (تنشيط قاعدي عالٍ لـ ERBB2) والآخر لديه طاقم بناء صاخب بجواره يهز الأرضية (محور إشارات Ca²⁺/p38). لقد حدد النموذج أن هذه التأثيرات الخارجية هي ما يجعل سلوك الخلايا فريداً بهذا الشكل.
كما أظهر الباحثون أن عندما يخطئ نموذجهم، فإن ذلك يكون أمراً جيداً في الواقع. فهذه "الأخطاء" تعمل مثل كشاف ضوئي، حيث تشير مباشرة إلى طفرات نادرة أو برامج بيولوجية خفية لم يكن النموذج يعرفها. على سبيل المثال، واجه النموذج صعوبة في التنبؤ بسلوك الخلايا ذات طفرة معينة تسمى BRAFV600E، مما ساعد في تأكيد أن هذه الطفرة تغير كيفية استجابة الخلايا للأدوية. وبالمثل، فقد وضع علامة على مجموعة من الخلايا التي بدت وكأنها تمتلك برنامجاً "هيكلياً" فريداً (مرتبط بهيكل الخلية) وآلية مقاومة محتملة تتضمن بروتيناً يسمى AMPK.
باختصار، هذا البحث لا يقدم فقط أداة تنبؤ أفضل؛ بل يقدم طريقة جديدة للنظر إلى البيولوجيا. فمن خلال الجمع بين القواعد الصارمة للكيمياء وقوة تعلم الأنماط للذكاء الاصطناعي، يقترح المؤلفون أن التنوع الذي نراه في الخلايا لا ينبع من كون الآلات الأساسية مختلفة، بل من البيئة الفريدة والإشارات الإضافية التي تتلقاها كل خلية. كما وجدوا أنه إذا قمت بتغذية النموذج بمعلومات محددة عن "الباب الأمامي" للخلية (مستويات EGFR)، فإن الخريطة الداخلية للنموذج تتغير تماماً، مما يظهر أن نفس البيانات يمكن فهمها إما كقصة جزيئية أو كقصة على مستوى النظام، اعتماداً على كيفية طرح السؤال. وبينما لا يزال لدى النموذج بعض الفجوات — مثل عدم القدرة على التنبؤ بدقة بكيفية تفاعل الخلايا مع تركيبات معينة من الأدوية — إلا أنه نجح في تحويل تلك الفجوات إلى خارطة طريق لاكتشاف بيولوجيا جديدة، مما يثبت أن معرفة ما لا تفهمه قد تكون بنفس قيمة معرفة ما تفهمه.
ملخص تقني: النماذج الميكانيكية العميقة تكشف عن الدوافع الخارجية للمسارات المسببة لتباين إشارات MAPK في خلايا الثدي
بيان المشكلة تُحكم الاستجابات الخلوية للمؤثرات البيئية بواسطة مسارات الإشارات، مثل شلال EGFR/MAPK، حيث تشكل الديناميكيات مصير الخلية حتى ضمن المجموعات المتطابقة جينياً. ومع ذلك، فإن النهج الحسابية الحالية لنمذجة هذا التباين منقسمة إلى مسارين: النماذج الميكانيكية (عادةً المعادلات التفاضلية العادية، ODEs) التي ترمز طوبولوجيا المسار وحركيته، لكنها تواجه صعوبة في تمثيل التباين الناجم عن عوامل خارج المسار المُمَذْج صراحةً. ومن ناحية أخرى، طرق تعلم التمثيل (مثل PCA، والمشفرات التلقائية/autoencoders) التي ترسم خرائط للبيانات الجينومية عالية الأبعاد إلى فضاءات منخفضة الأبعاد، لكنها تفتقر إلى التفسير الميكانيكي، وتفشل في شرح كيفية تنفيذ حالات خلوية معينة لوظائف بيولوجية. إن التوفيق بين هذين المنظورين لتفسير تباين الإشارات بطريقة تعتمد على البيانات وتكون مفسرة ميكانيكياً في آن واحد ظل تحدياً كبيراً.
المنهجية: النماذج الميكانيكية العميقة (DMMs) قدم المؤلفون النماذج الميكانيكية العميقة (DMMs)، وهو إطار عمل يربط بين تعلم التمثيل شبه المُوجه ونموذج معادلة تفاضلية عادية (ODE) ثابت لمسار إشارات EGFR/MAPK. تعمل البنية وفقاً لما يلي:
المُشفر (Encoder): تقوم شبكة عصبية (بنية مشفر تلقائي) بترميز القياسات الجزيئية عالية الأبعاد (النسخ، البروتينات، أو الفسفرة البروتينية) إلى تمثيل كامن (latent embedding) منخفض الأبعاد.
الرأس الميكانيكي (Mechanistic Head): يقوم هذا التمثيل الكامن بتحديد الانحرافات في المعاملات الحركية لنموذج ODE ثابت. يحاكي نموذج ODE ديناميكيات مستقبلات تيروزين كيناز (EGFR، ERBB2، ERBB3) وأحداث الفسفرة في المصب (MEK، ERK، p90rsk) تحت تأثير تنبيهات مختلفة (تحفيز EGF، مثبطات MEK، مثبطات RTK).
المفكك (Decoder): يقوم مفكك موازٍ بإعادة بناء الميزات الجزيئية الأساسية من التمثيل الكامن.
التدريب: يتم تدريب النظام بالكامل باستخدام التحسين القائم على التدرج (المُنفذ عبر JAX وAMICI). تجمع دالة الخسارة بين الملاءمة الميكانيكية لبيانات الإشارات الديناميكية (خاضعة للإشراف) وإعادة بناء الميزات (غير خاضعة للإشراف)، مما يضمن أن التمثيل الكامن يلتقط المعلومات ذات الصلة بكل من ديناميكيات المسار والحالة الجزيئية.
على عكس المعادلات التفاضلية العادية العصبية (Neural ODEs) أو المعادلات التفاضلية الشاملة (Universal Differential Equations)، فإن الـ ODEs في نماذج DMMs ثابتة ومفسرة ميكانيكياً، وتعمل في الفضاء المرصود بدلاً من الفضاء الكامن.
النتائج الرئيسية طُبق إطار العمل على بيانات CyTOF متعددة الإرسال من 63 خطاً خلوياً لسرطان الثدي (مجموعة التدريب: 58 خطاً؛ مجموعة الاختبار: 5 خطوط مستبعدة).
التعميم والأداء: أظهرت نماذج DMMs قدرة قوية على التعميم على الخطوط الخلوية المستبعدة، محققةً خطأً جذرياً متوسط المربع (RMSE) منخفض في بيانات الاختبار مقارنة ببيانات التدريب. كما تفوقت النماذج المدربة على مدخلات متعددة الأنماط (التي تجمع بين النسخ، والبروتينات، وCyTOF) على تلك المدربة على أنماط مفردة، وتفوقت على نماذج الانحدار الخطي الأساسية.
تحديد دوافع التباين: كشف تحليل الحساسية لانحرافات المعاملات المستنتجة أن معظم تباين الإشارات مدفوع بـ عوامل خارجية للمسار بدلاً من التباين في مكونات MAPK الأساسية. تم تحديد محورين رئيسيين:
تنشيط ERBB2 الأساسي: معامل يتحكم في الفسفرة الذاتية لـ ERBB2، ويرتبط بقوة بمستويات p-ERBB2 الأساسية وتنشيط ERK المستدام.
محور إشارات Ca²⁺/p38: معامل يتحكم في انتقال الإشارة من ERK إلى p90rsk. ارتبط هذا المحور بمؤشرات إشارات الكالسيوم، وPLCG، وp38، مما يشير إلى أن التداخل بين مسار MAPK وإشارات الكالسيوم ينظم تضخيم الإشارة إلى p90rsk.
التكامل الميكانيكي مقابل تعزيز الميزات: أدى دمج مستويات نسخ EGFR مباشرة في نموذج الـ ODE (عبر إعادة تحجيم معدلات التصنيع) إلى إعادة تشكيل التمثيل الكامن. وبينما أدى ذلك إلى إلغاء الفصل بين الأنواع الفرعية اللمعية (luminal) والأساسية (basal) في الفضاء الكامن (لأن مستويات EGFR أصبحت محسوبة صراحةً بواسطة الآلية)، إلا أنه حافظ على القدرة التنبؤية. أثبت هذا أن DMMs يمكنها التبديل بين "عدسة جزيئية" (حيث يقود EGFR الأنماط الظاهرية عبر الآلية) و"عدسة نظمية" (حيث يعمل EGFR كعلامة للنوع الفرعي) بقدر متساوٍ من الأداء التنبئي.
فشل النموذج كأداة للاكتشاف: تم تحليل التناقضات بين تنبؤات النموذج والبيانات التجريبية (فشل النموذج) بشكل منهجي لتحديد الرؤى البيولوجية:
الطفرات النادرة: رُبط الخطأ العالي في خطوط خلوية معينة (مثل HCC2218، DU4475) بطفرات نادرة (حذف MET exon-14، وBRAF V600E) تغير استجابات الإشارات، والتي فشلت نماذج الانحدار القياسية في التقاطها دون وجود أولويات ميكانيكية.
البرامج المتكررة: أشارت الأخطاء المنهجية في مجموعات من الخطوط الخلوية إلى برامج متكررة، بما في ذلك محور مقاومة مثبطات BRAF MEK المرتبط بـ AMPK (المرتبط بانخفاض PRKAA2) وبرنامج هيكلي/إفرازي (مرتبط بطفرات Rho GTPase وجينات مثل TAGLN2 وGSN) الذي يقود التباين في إشارات ERK-to-p90rsk.
القيود والعقبات: أشار تحليل حجم مجموعة التدريب إلى عوائد متناقصة بعد 30 خطاً خلوياً، مما يوحي بأن السقف التنبئي الحالي محدود بـ بنية النموذج الميكانيكي وليس بكمية البيانات. أدى تعزيز الميزات إلى تحسين الأداء قليلاً، لكن التكامل الميكانيكي لفرضيات بيولوجية محددة (مثل فسفرة ERBB2 المدفوعة بالليجاند) لم يؤدِ دائماً إلى مكاسب، حيث عمل كمرشح للفرضيات.
الأهمية والادعاءات يزعم البحث أن نماذج DMMs تقدم إطار عمل عام لدمج الآلية مع التعلم، مما يتيح:
القابلية للتفسير: توفير تفسير ميكانيكي لكيفية تنفيذ الحالات الخلوية المستمدة من البيانات لوظائفها، مما يسد الفجوة بين الارتباط الإحصائي والسببية البيولوجية.
اكتشاف التنظيم الخارجي: النجاح في إرجاع التباين إلى عوامل خارجية للمسار (مثل التداخل بين Ca²⁺/p38) التي لم تستطع النماذج الميكانيكية الثابتة وحدها التقاطها.
التحليل المنهجي للفشل: تحويل فشل النموذج إلى مسار منهجي لاكتشاف البيولوجيا المرشحة، مثل الطفرات النادرة والبرامج الإشارية المتكررة، من خلال رصد التناقضات التي تغفل عنها النماذج القياسية.
الاستقلالية عن النمط (Modality Agnosticism): إطار العمل قابل للتوسع طبيعياً لأنماط أخرى (مثل التصوير) ويمكنه التوفيق بين الأوصاف الجزيئية والنظمية لنفس البيانات دون فرض خيار بينهما.
يخلص المؤلفون إلى أنه بينما تواجه نماذج DMMs حالياً قيوداً في بنية النموذج (تحديداً الطبيعة الخشنة لـ ODEs)، إلا أنها تمثل خطوة مهمة نحو توحيد تعلم التمثيل القائم على البيانات مع النمذجة البيولوجية الميكانيكية.