Determination of active forces in actomyosin systems as inverse source problems for the Stokes equation

ملخص تقني: تحديد القوى النشطة في أنظمة الأكتوميوسين كمسائل مصدر عكسية لمعادلة ستوكس

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

تتناول الدراسة إطارين فيزيائيين متميزين:

1. الأنظمة المحتواة: شبكات الأكتوميوسين المحتواة داخل قطرات زيت في الماء، والتي تُنمذج باستخدام شروط روبين (الانزلاق) الحدية.
2. الأنظمة الحرة (الكتلية): الشبكات المعلقة بحرية دون احتواء مرئي، والتي تُنمذج باستخدام شروط ديريكليه غير المتجانسة.

تتمثل إحدى التحديات الحرجة المحددة في عدم اكتمال البيانات التجريبية؛ وتحديداً، حقل الضغط $p$ غير مقاس، وتُفترض اللزوجة $\eta$ ثابتة بناءً على الحجم الضئيل لشبكة الخيوط مقارنة بالسوائل المحيطة.

المنهجية
يستخدم المؤلفون معادلة ستوكس كنموذج أمامي، مع دمج شرط عدم الانضغاط ($\nabla \cdot v = 0$) ولزوجة متغيرة مكانياً $\eta$. تم بناء الإطار الرياضي بناءً على الخطوات التالية:

1. تحليل المسألة الأمامية: يحلل المؤلفون بدقة جودة المسألة (well-posedness) لمسألة ستوكس لكلا نوعي الشروط الحدية (روبين وديريكليه). كما يحددون فضاءات الدوال المناسبة (مثل $H^1_{\parallel}(\Omega)^d$ لشروط روبين) ويثبتون الوجود الفريد والاستقرار لزوج السرعة والضغط $(v, p)$ لقوة معطاة $f$ باستخدام شرط inf-sup ومتراجحة كورن.
2. صياغة المسألة العكسية: تُعرف المسألة العكسية بأنها إيجاد $f$ بحيث يطابق المؤثر الأمامي $A$ (الذي يربط القوة بالسرعة) البيانات المقاسة $v_{meas}$. ونظراً لغياب بيانات الضغط، تُعامل المسألة كعملية تحديد للقوة بناءً على السرعة فقط.
3. التفكيك النظري: يحلل البحث تفكيك هلمهولتز للقوة إلى مكونات خالية من التباعد (solenoidal) وخالية من الدوران (irrotational). وقد تبين أنه بالنسبة للزوجة الثابتة، فإن المكون الخالي من التباعد للقوة يتوافق مع الحد المعتمد على السرعة في معادلة ستوكس، بينما يتوافق المكون الخالي من الدوران مع تدرج الضغط. وبناءً على ذلك، ومن دون بيانات الضغط، يمكن إعادة بناء المكون الخالي من التباعد فقط من القوة النشطة بشكل فريد.
4. اشتقاق الملحق (Adjoint): لتمكين الحل العددي، اشتق المؤلفون مشتقات فريشيه للمؤثرات الأمامية وملحقاتها المقابلة (ملحقات فضاء باناخ وفضاء هيلبرت). ويعد هذا الاشتقاق ضرورياً للتنظيم القائم على التدرج (gradient-based regularization).
5. التنظيم العددي: يتم حل المسألة العكسية باستخدام طريقة تكرار لاندويبر (Landweber iteration). تقوم الخوارزمية بتقليل الفرق بين السرعة المقاسة والسرعة المحاكية. ولضمان الاستقرار ضد الضجيج، يتم إنهاء التكرار باستخدام مبدأ التباين (discrepancy principle). وتكون التخمين الأولي عبارة عن إعادة بناء "ساذجة" مشتقة مباشرة من حد السرعة في معادلة ستوكس ($f_0 = -\nabla \cdot \eta D(v_{meas})$).

النتائج الرئيسية
تم التحقق من صحة المنهجية باستخدام كل من البيانات الاصطناعية والقياسات التجريبية:

• البيانات الاصطناعية: أُجريت عمليات إعادة البناء على نطاقات ثنائية الأبعاد تمثل كل من القطرات (روبين) والعينات الحرة (ديريكليه).
  • للبيانات الخالية من الضجيج: نجحت عملية تكرار لاندويبر في استعادة الحقيقة الأرضية (ground truth) للقوة، بما في ذلك أنماط الدوامات المعقدة، مع خطأ نسبي منخفض ($\tilde{\epsilon} \approx 0.3$ لروبين، و $\approx 0.73$ لديريكليه).
  • للبيانات المشوبة بالضجيج: أظهرت الطريقة متانة ملحوية. فحتى مع مستويات ضجيج نسبية تصل إلى $\tilde{\delta} = 1.37$، ظل خطأ إعادة البناء النسبي محدوداً ($\tilde{\epsilon} \approx 0.29$)، مما حافظ على البنية العامة والقدر لشدة القوة، رغم فقدان التفاصيل الدقيقة عند مستويات الضجيج العالية.
• البيانات التجريبية: طُبقت الخوارزمية على بيانات PIV لشبكات الأكتوميوسين في القطرات وعينات حرة.
  • القطرات: أظهر حقل القوة المعاد بناؤه مكوناً خالياً من التباعد يتوافق مع اتجاه التدفق الملاحظ داخل القطرة.
  • الأنظمة الحرة: أُجريت عمليات إعادة البناء على شبكات غير محتواة. ويشير المؤلفون إلى أن إعادة البناء في الأنظمة الحرة تعد أكثر صعوبة بسبب احتمال أهمية المكون المفقود المرتبط بالدوران (الضغط) والغموض حول ما إذا كانت التدفقات الملاحظة تنبع فقط من نشاط الأكتوميوسين أم من التفاعلات بين السائل والمحيط.

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

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

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