Next generation Human Neocortical Neurosolver model with biophysically realistic somatic, dendritic and synaptic dynamics

ملخص تقني: نموذج المحلل العصبي للقشرة المخية البشرية من الجيل التالي (HNN)

بيان المشكلة
يُعد المحلل العصبي للقشرة المخية البشرية (HNN) أداة مفتوحة المصدر مصممة لتفسير إشارات تخطيط كهربائية الدماغ (EEG) وتخطيط مغناطيسية الدماغ (MEG) البشرية من خلال سد الفجوة بين التسجيلات واسعة النطاق وديناميكيات الدوات دقيقة النطاق. ومع ذلك، اعتمد نموذج HNN الأصلي على معلمات فيزيائية حيوية مستمدة أساساً من دراسات القوارض والحيوانات غير البشرية، نظراً لأن الديناميكيات العصبية البشرية لم تكن مشخصة جيداً وقت تطويره. ونتيجة لذلك، افتقر النموذج الأصلي إلى قيود محددة مستمدة من الأنسجة البشرية، لا سيما فيما يتعلق بديناميكيات الكالسيوم (Ca2+) في الشجيرات العصبية وحركية المشابك العصبية النوعية، مما حد من قدرته على الفصل بدقة بين المساهمات الجسدية، والمشبكية، والشجيرية لإشارات الـ E-/MEG البشرية.

المنهجية
يقدم المؤلفون نموذجاً محدثاً بشكل كبير للقشرة المخية، يُسمى duecker_ET_model مدمجاً في مجموعة برمجيات HNN. تضمنت المنهجية ما يلي:

1. ضبط المعلمات الفيزيائية الحيوية: تم ضبط الخلايا العصبية الهرمية (الطبقة 2/3 والطبقة 5) والخلايا العصبية البينية باستخدام بيانات فيزيولوجية كهربائية في المختبر (in vitro) من أنسجة بشرية.
   • الخلايا العصبية الهرمية من النوع (L5 ET) - الناقلة خارج التلفنسيفال: تم نمذجتها باستخدام مورفولوجيا مخفضة (9 حجرات) ولكن تم ضبطها لإعادة إنتاج 12 مقياساً فيزيولوجياً كهربائياً جسدياً من خلايا L5ET البشرية. ومن الأهمية بمكان أن النموذج يدمج قنوات شجيرية نشطة لمحاكاة نبضات الكالسيوم (Ca2+) الشجيرية، وهي ميزة كانت غير مقيدة سابقاً في HNN.
   • الخلايا العصبية الهرمية (L2/3): تم ضبطها لإعادة إنتاج الديناميكيات الجسدية والنبضات الشجيرية القصيرة الملاحظة في خلايا L2/3 البشرية.
   • الخلايا العصبية البينية: تم ضبط نموذج أحادي الحجرة ليتوافق مع ديناميكيات الخلايا عديمة الأشواك البشرية، مع دمج تيار Ih لإعادة إنتاج نسب الـ (sag ratios).
2. تحديثات المشابك العصبية: أدخل النموذج آليات مشبكية واقعية بيولوجياً، بما في ذلك:
   • نموذج جديد لمشبك NMDA يتضمن تيارات الكالسيوم المشبكية وحصار المغنيسيوم (Mg2+) المعتمد على الجهد.
   • حركية مُحدثة لمستقبلات GABA-A لتثبيط الشجيرات الشبيه بـ SOM+ (ارتفاع/اضمحلال أبطأ).
   • نموذج مشبك GABA-B واقعي بيولوجياً يتضمن اقتران البروتين (G-protein coupling).
3. استراتيجية التحسين: لضبط النموذج على البيانات التجريبية، استخدم المؤلفون استراتيجية تطور مصفوفة التغاير (CMA-ES). تم تطوير دالة خسارة مبتكرة لتحسين كل من شكل الموجة (عبر معامل ارتباط بيرسون) والسعة (عبر جذر متوسط مربع الخطأ)، مع فرض عقوبات على الحلول التي تطابق السعة ولكن تفشل في التقاط شكل الموجة. كما تم تطبيق عقوبات إضافية لضمان أنماط إطلاق نبضات بيولوجية معقولة (على سبيل المثال، منع الإطلاق المفرط لخلايا L5ET مقارنة بـ L2/3 في نوافذ الاستجابة المبكرة).
4. أهداف المحاكاة: تم اختبار النموذج المحدث مقابل:
   • الجهود المرتبطة بالحدث (ERPs): الاستجابات اللمسية في القشرة الحسية الجسدية الأولية (SI).
   • التذبذبات: مركب ألفا-بيتا لإيقاع "مو" (mu rhythm) الحسي الجسدي وتذبذبات غاما (آلية PING).

النتائج الرئيسية

• الدقة الخلوية: نجح نموذج duecker_ET_model في إعادة إنتاج مجموعة واسعة من البيانات البشرية في المختبر، بما في ذلك جهود الغشاء أثناء الراحة، وأشكال جهد الفعل، وترددات الرنين، ونبضات الكالسيوم الشجيرية في خلايا L5ET.
• توليد ثنائي القطب (Dipole Generation): أدى إدراج نبضات الكالسيوم الشجيرية وتحديث المقاومة المحورية (التي انخفضت من 200 إلى 100 أوم.سم في L5ET) إلى تغيير عملية توليد ثنائي القطب بشكل كبير. يحقق النموذج ملاءمة وثيقة لأشكال موجات MEG التجريبية (RMSE = 10.82، معامل بيرسون r = 0.99) مع عامل قياس قدره 1,000، مما يعني وجود حوالي 200,000 خلية عصبية هرمية نشطة، مقارنة بالنموذج السابق الذي كان يتطلب حوالي 600,000 خلية عصبية (عامل قياس 3,000).
• الفصل الآلي: يسمح النموذج المحدث بفصل المساهمات في موجة ERP:
  • الذروة الموجبة الأولية (M25) تنشأ من التحفيز القريب والتثبيط الشجيري القمي.
  • الانحراف السالب (M35) مدفوع بالتثبيط حول الجسدي (perisomatic).
  • المنخفض البارز (M70) يقابل التحفيز البعيد الذي يطلق نبضات الكالسيوم في خلايا L5ET.
  • المكونات الموجبة المتأخرة (M100, M135) ناتجة عن تسلسلات إيقاعية من الاستثارة والتثبيط.
• الديناميكيات التذبذبية: يعيد النموذج إنتاج الطبيعة "الانفجارية" العابرة لإيقاعات ألفا وبيتا "مو"، وتماثل قمم وقيعان الموجات الملاحظة في بيانات MEG. كما نجح في محاكاة تذبذبات غاما المستمرة والعابرة عبر آليات PING، مع تعديلات في المعلمات تعكس ديناميكيات الانفجار المتميزة لنموذج L5ET الجديد مقارنة بالإصدار السابق.

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

1. القيود الخاصة بالبشر: تجاوز المعلمات المستمدة من الحيوانات لاستخدام البيانات البشرية في المختبر مباشرة للضبط الجسدي والشجيري.
2. الديناميكيات الشجيرية: كونه، حسب علم المؤلفين، أول إطار عمل فيزيائي حيوي لنمذجة E-/MEG يدمج صراحةً ديناميكيات الكالسي_م الشجيرية في كل من الخلايا الهرمية L2/3 وL5، مما يسمح للباحثين باستقصاء كيفية مساهمة هذه الآليات المحددة في الإشارات واسعة النطاق.
3. التفسير متعدد المقاييس: توفير أداة تم التحقق من صحتها للفصل بين المساهمات المحددة للجهود الجسدية، والنبضات الشجيرية، والتيارات المشبكية في توقيعات E-/MEG، مما يسهل اختبار الفرضيات بدقة أكبر فيما يتعلق بالمولدات العصبية للنشاط الدماغي البشري.

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