Multiscale order, flocking and phenotypic hysteresis in the cellular Potts model of epithelia

ملخص تقني: النظام متعدد المقاييس، التجمع (Flocking)، والتباطؤ المظهري في نموذج "بوتس" الخلوي للأنسجة الطلائية

بيان المشكلة
تتناول الورقة البحثية السؤال الجوهري حول كيفية تأثير كل من الهجرة الخلوية الجماعية والتنظيم المكاني للأنسجة على بعضهما البعض في الأنسجة الطلائية. وبينما قدمت الدراسات الحديثة مفهوم "النظام متعدد المقاييس" (MsO) — وهو التعايش بين أطوار غير متوازنة مختلفة عبر مقاييس طولية مختلفة — إلا أن الآليات التي تقود تطور النظام النيماتي (nematic order) واسع النطاق والطبيعة المحددة للتنظيم "السداسي النيماتي" (hexanematic) لا تزال غير واضحة. علاوة على ذلك، فإن العلاقة بين الهجرة الموجهة (التجمع/Flocking) وآليات التراصف الذاتي الكامنة في الأنسجة، مقابل الهجرة غير الموجهة القائمة على التداخل (intercalation)، تتطلب مزيدًا من التوضيح. يهدف المؤلفون إلى رسم مخطط الطور الكامل لهجرة الأنسجة الطلائية وربط أنماط الهجرة بتمثلات متميزة للنظام متعدد المقاييس.

المنهجية
يستخدم المؤلفون عمليات محاكاة عددية واسعة النطاق باستخدام نموذج "بوتس الخلوي" (Cellular Potts Model - CPM)، وهو نموذج قائم على الشبكة حيث تُمثل الخلايا كمجموعات متصلة من البكسلات. وتشمل السمات المنهجية الرئيسية ما يلي:

• تنفيذ النموذج: تم تنفيذ نموذج CPM على وحدات معالجة الرسومات (GPUs) باستخدام لغة CUDA، مما سمح بإجراء عمليات محاكاة ذات نطاق غير مسبوق (يصل إلى شبكات $L=2^{13}$، مما يمثل ما يقرب من 700,000 خلية).
• دمج الهيكل الخلوي للأكتين: يتضمن النموذج "نموذج أكتين" (Act-model) صريحًا حيث يمثل المجال $\phi_i$ كثافة خيوط الأكتين في لاميليوبوديوم (lamellipodium) الخلية. هذا المجال يدفع سلوك "التدحرج" (tread-milling)، حيث تترجم البكسلات نحو المناطق ذات كثافة الأكتين العالية، مما يولد الحركة.
• ديناميكيات هاميلتون (Hamiltonian Dynamics): يتطور النظام عبر خطوات "مونت كارلو" لآلية "متروبوليس-هاستينغز" بناءً على دالة هاميلتون $H$ التي تتضمن تفاعلات السطح الجانبية ($J$)، وقيود المساحة ($\lambda_{cell}$)، وطاقة تجميع الهيكل الخلوي ($\lambda_{cyto}$).
• استكشاف فضاء المعلمات: يقوم المؤلفون بتغيير معلمات نشاط الهيكل الخلوي ($\lambda_{cyto}$ و $\Phi$) ودرجة الحرارة الفعالة/الضجيج ($T$) بشكل منهجي عبر 901 تكوينًا مختلفًا.
• تحليل معلمات النظام: لتوصيف الأطوار، يحسب المؤلفون ثلاثة مجالات توجيهية:
  • $\psi_1$: اتجاه المحور القطبي (مرتبط بحركة الخلية).
  • $\gamma_2$: نظام الشكل النيماتي (Nematic shape order) ($p=2$).
  • $\gamma_6$: نظام الشكل السداسي (Hexatic shape order) ($p=6$).
    يتم تجريد هذه المجالات خشنًا عبر مقاييس طولية $\ell$ لتحديد معلمات النظام المعتمدة على المقياس $\Gamma_p(\ell)$ وأسس اضمحلالها $\eta_p$.
• التصنيف: يتم استخدام خوارزمية تصنيف الغابة العشوائية لاتخاذ القرار (RDFCA) بالاقتران مع خوارزمية "المكعبات المارة" (Marching Cubes) لرسم حدود الطور في فضاء المعلمات ثلاثي الأبعاد.

النتائج الرئيسية
تكشف عمليات المحاكاة عن مخطط طور غني يتميز بأنواع متعددة من النظام التوجيهي التي تتعايش أو تهيمن في مقاييس طولية مختلفة:

1. النظام السداسي النيماتي (Hexanematic Order): عند مستويات متوسطة من نشاط الهيكل الخلوي، يظهر النظام نظامًا "سداسيًا نيماتيًا". عند المقاييس الطولية المماثلة لحجم الخلية ($\ell \approx \ell_{cell}$)، تتبنى الخلايا مورفولوجيا سداسية تقريبًا، مما يؤدي إلى هيمنة النظام السداسي ($\Gamma_6 > \Gamma_2$). ومع ذلك، عند المقاييس الطولية الأكبر، يصبح النظام النيماتي هو المهيمن ($\Gamma_2 > \Gamma_6$). يحدث هذا الانتقال عند مقياس طولي محدد $\ell^\times$.
2. انتقال التجمع (Flocking Transition): مع زيادة معدل بلمرة الأكتين ($\lambda_{cyto}$)، ينتقل النظام إلى حالة تجمع طويلة المدى. ويتميز هذا بوجود نظام طويل المدى في المحور القطبي ($\eta_1 = 0$).
3. ظهور النظام النيماتي: بالتزامن مع بدء عملية التجمع، يظهر نظام نيماتي شبه طويل المدى عند مقاييس طولية أكبر بكثير من حجم الخلية نتيجة للتفاعل بين الحركة الموجهة والتفاعلات الجانبية.
4. سلوك الأسس: تحدد الدراسة مسارًا محددًا حيث يؤدي زيادة $\lambda_{cyto}$ إلى انتقال يصبح فيه النظام النيماتي طويل المدى بالكامل ($\eta_2 = 0$) بينما يظل النظام السداسي شبه طويل المدى ($\eta_6 \approx 1$) ومحصورًا في مقاييس طولية قصيرة، بشكل مستقل عن أي زيادة إضافية في نشاط الهيكل الخلوي.
5. التباطؤ المظهري (Phenotypic Hysteresis): عندما يكون الضجيج ($T$) منخفضًا بما يكفي للسماح بالتبلور، فإن دورة معامل نشاط الهيكل الخلوي $\lambda_{cyto}$ عبر مرحلة الانصهار تؤدي إلى حلقة تباطؤ (hysteresis loop). يظهر النظام حالات مختلفة (خلايا متحركة مقابل خلايا قرص عسل غير متحركة) اعتمادًا على اتجاه مسح المعامل، وهو ما يشبه الانتقال الظهاري-المتوسطي (EMT).

الأهمية والادعاءات
تدعي الورقة توفير إطار كمي لفهم التغذية الراجعة بين الهجرة الجماعية والتنظيم المكاني في الأنسجة الطلائية. ومن خلال الجمع بين مورفولوجيا الخلية ونشاط الهيكل الخلوي، يثبت المؤلفون أن:

• التنظيم متعدد المقاييس: يمكن تفسير التنظيم المعقد للأنسجة الطلائية من خلال لغة هيدروديناميكا البلورات السائلة، وإن كان ذلك بأبعاد شاذة تتجاوز حدود الفيزياء الإحصائية للتوازن.
• آلية التجمع: تشير الدراسة إلى وجود رابط سببي بين ظهور النظام السداسي النيماتي وظهور التجمع، حيث يتطور النظام النيماتي عند المقاييس الكبيرة بينما يستمر النظام السداسي عند المقاييس الصغيرة.
• اللدونة المظهرية: يوفر رصد التباطؤ في النموذج آلية محتملة لفهم اللدونة المظهرية، وتحديدًا كيف يمكن للخلايا الاحتفاظ بالنمط الميزنشيمي حتى بعد سحب المحفزات، وهي ظاهرة ذات صلة بسرطان الانتشار (metastasis).

يخلص المؤلفون إلى أنه بينما يلتقط نموذجهم الفيزياء الأساسية لهذه الانتقالات، فإن الفهم الكامل لللدونة المظهرية سيتطلب مستقبلاً دمج مسارات كيميائية حيوية محددة تتجاوز دورة بلمرة الأكتين. ويعمل هذا العمل كجسر بين نظرية المادة النشطة والملاحظات البيولوجية في الطبقات الخلوية المعاد تكوينها.