Viscous vertex model for active epithelial tissues

إليك ملخص تقني مفصل للورقة البحثية بعنوان "نموذج الرأس المتقطع اللزج للأنسجة الظهارية النشطة" (Viscous vertex model for active epithelial tissues) لـ لين، وتيلي، وروبريخت.

1. بيان المشكلة

تتصرف الأنسجة الظهارية وتجمعات الخلايا الكثيفة كأدوات نشطة ذات لزوجة مرنة (viscoelastic). وبينما تصف نماذج الرؤوس المتقطعة (vertex models) القياسية ميكانيكا الأنسجة بفعالية باستخدام طاقات مرنة (المساحة والمحيط) والتوترات البينية، إلا أنها عادةً ما تنمذج التبديد من خلال الاحتكاك مع الركيزة فقط (الاحتكاك المحلي بين الرؤوس المتقطعة والمصفوفة خارج الخلوية الساكنة).

ومع ذلك، في العديد من السياقات البيولوجية — مثل الأنسجة الجنينية الطافية بحرية، أو العضويات (organoids)، أو الطبقات أحادية الطبقة المعلقة — يكون الاحتكاك مع الركيزة ضئيلاً للغاية. في هذه الأنظمة، تهيمن القوى اللزجة الداخلية (اللزوجة القشرية والكتلية) على عملية التبديد. وتفشل النماذج الحالية غالباً في:

1. معاملة اللزوجة عند الوصلات (البينية) واللزوجة الكتلية على قدم المساواة.
2. الحفاظ على الثبات الدوراني (الاستقلال عن الإطار المرجعي) عند تعريف القوى اللزجة.
3. البقاء متسقة رياضياً في حد انعدام احتكاك الركيزة (حيث يصبح النظام منفردًا/singular بسبب عدم وجود قيود على الحركة الانتقالية والدورانية العالمية).
4. توفير رابط واضح بين لزوجة الخلايا المجهرية والريولوجيا (علم التدفق) الكلية للنسيج.

2. المنهجية

قام المؤلفون بتطوير نموذج رأس متقطع لزج ثابت دورانياً يدمج مصدرين متميزين للتبديد:

• اللزوجة البينية بين الخلايا ($\eta^{(s)}$): التبديد الناتج عن الحركة النسبية للرؤوس على طول حواف الخلايا.
• اللزوجة الكتلية للخلية ($\eta^{(b)}$): التبديد الناتج عن الحركة النسبية بين رؤوس الخلية والمركز الهندسي للخلية.

خطوات الصياغة الرئيسية:

• معدلات التبديد: يحدد النموذج معدلات التبديد ($S_{interface}$ و $S_{bulk}$) بناءً على مربع معدلات الاستطالة المسقطة على طول الحواف والروابط بين الرأس والمركز، على التوالي. يؤدي هذا إلى توترات لزجة تتناسب مع فروق السرعة: $\Lambda^{(viscous)} \propto \mathbf{t} \cdot (\mathbf{v}_j - \mathbf{v}_i)$.
• الثبات الدوراني: على عكس الصياغات البديلة التي تقوم ببساطة بتخميد السرعات النسبية ($\dot{\mathbf{r}}_i - \dot{\mathbf{r}}_j$)، يقوم هذا النموذج بإسقاط السرعات على متجهات المماس الوحدة للحواف. يضمن هذا أن الدورات الصلبة العالمية لا تولد تبديداً زائفاً، مما يحقق موضوعية المادة (material objectivity).
• صيغة المصفوفة: يتم تجميع معادلات توازن القوة في نظام خطي $C \cdot \mathbf{v} = \mathbf{F}(t)$، حيث $C$ هي مصفوفة معامل الاحتكاك واللزوجة التي تجمع بين احتكاك الركيزة، واللزوجة البينية، واللزوجة الكتلية.
• التعامل مع الاحتكاك الصفري: لمعالجة حالة التفرد في المصفوفة $C$ عندما يكون احتكاك الركيزة $\gamma \to 0$، استخدم المؤلفون صيغة مضاعف لاغرانج. قاموا بتوسيع النظام إلى مصفوفة نقطة السرج $C_{ext}$ من خلال فرض قيود حركية عالمية (حفظ الزخم الخطي والزاوي الإجمالي). يعمل هذا على تنظيم النظام، مما يسمح بالمحاكاة عند انعدام الاحتكاك تماماً دون حدوث تباعد عددي.
• بروتوكول الريولوجيا: تم تقديم بروتوكول "قص الشريحة" (slab-shear) لاستخراج لزوجة النسيج الفعالة، متوسطة الحجم، عبر تطبيق معدلات انفعال القص وقياس الإجهاد اللزج الناتج.

3. المساهمات الرئيسية

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

4. النتائج

أ. لزوجة النسيج متوسطة الحجم (Coarse-Grained)

• اللزوجة قصيرة المدى: في غياب إعادة ترتيب الخلايا (انتقالات T1)، تتناسب لزوجة النسيج قصيرة المدى $\eta_{tissue}^{(ST)}$ خطياً مع مجموع اللزوجات البينية ($\eta_s$) والكتلية ($\eta_b$).
  • النتيجة التحليلية: بالنسبة لشبكة سداسية، $\eta_{tissue}^{(ST)} = \frac{\sqrt{3}}{4}\eta$ (بافتراض أن $\eta_s = \eta_b = \eta$).
• اللزوجة طويلة المدى: تحت تأثير القص المستمر الذي يسمح بإعادة ترتيب الخلايا، تتناسب اللزوجة طويلة المدى أيضاً خطياً مع لزوجات الخلايا، ولكنها تكون عموماً أقل من قيمة المدى القصير بسبب استرخاء الإجهاد عبر انتقالات T1.
• نطاقات القص (Shear Banding): في وجود احتكاك عالٍ مع الركيزة، يتمركز تدفق القص بالقرب من الحدود (نطاقات القص). في حد انعدام الاحتكاك، يختفي هذا النطاق ويصبح ملف التدفق منتظماً، مما يشبه السائل النيوتوني.

ب. ديناميكيات الأنسجة النشطة

تم تطبيق النموذج على سيناريوهين نشطين:

1. الجر القطبي النشط (هجرة الحافة القائدة):
   • تؤدي اللزوجة المنخفضة إلى تدفقات تشبه الاضطراب مع وجود العديد من العيوب المتحركة.
   • تعمل زيادة اللزوجة على تعزيز الارتباط المكاني لتوجهات الخلايا، وتقليل عدد العيوب، وتعزيز نمط دوران متماسك عالمي (دوامي) مع عيوب $+1/2$ و $-1/2$ محددة جيداً.
2. الإجهادات النيماتية النشطة (المعتمدة على الشكل):
   • تؤدي اللزوجة المنخفضة إلى خلايا مستطيلة، ووريدات متعددة الخلايا، وكثافة عالية للعيوب.
   • تعمل اللزوجة العالية على تثبيت العيوب الطوبولوجية المحددة جيداً وتعزيز الارتباطات المكانية طويلة المدى في توجه الخلايا، مما يطابق الملاحظات التجريبية في الأنسجة النيماتية النشطة.
   • تعمل اللزوجة كمعامل تحكم، مما يدفع بالانتقال من "الاضطراب النشط" إلى الحركة المتماسكة والمنظمة.

ج. الاستقرار العددي

• يتناسب شرط المصفوفة (condition number) للمصفوفة الموسعة $C_{ext}$ مع $N_c^{2.5}$ (حيث $N_c$ هو عدد الخلايا)، مما يجعلها قابلة للمعالجة عددياً لأنظمة تصل إلى $\sim 6 \times 10^5$ خلية.
• نجح الإطار في محاكة الأنسجة الطافية بحرية دون احتكاك اصطناعي، مما يلتقط الأنماكن "الرخوة" (التشوهات ذات التبديد الصفري) التي لا يتم كبحها إلا عندما تكون اللزوجة الكتلية غير صفرية.

5. الأهمية

يوفر هذا العمل جسراً حاسماً بين النماذج الخلوية المنفصلة و النظريات المستمرة للمادة النيماتية النشطة.

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

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