نمذجة إمكانيات العمل في الألياف العصبية الميالينية The modelling of the action potentials in myelinated nerve fibres

شارك:
المجلة: Biomechanics and Modeling in Mechanobiology، المجلد: 25، العدد: 1
DOI: https://doi.org/10.1007/s10237-025-02030-w
PMID: https://pubmed.ncbi.nlm.nih.gov/41504943
تاريخ النشر: 2026-01-08
المؤلف: Kert Tamm وآخرون
الموضوع الرئيسي: تنظيم ووظيفة قنوات الأيونات

نظرة عامة

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

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

مقدمة

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

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

مناقشة

تتناول قسم المناقشة في ورقة البحث نموذج ليبرشتاين لانتشار جهد الفعل (AP) في المحاور غير الميالينية، مع التأكيد على أهمية تضمين الحث جنبًا إلى جنب مع آلية أيونات هودجكين-هاكسلي (HH). يعدل نموذج ليبرشتاين الإطار الكلاسيكي لـ HH من خلال دمج الحث، مما يؤدي إلى معادلة حاكمة هيبرولية تضمن انتشار الإشارة السببية، مما يتناقض مع الطبيعة البارابولية لمعادلة الكابل التقليدية التي تسمح بسرعات غير سببية غير محدودة. تعتبر هيبرولية النموذج حاسمة لوصف ديناميات جهود الفعل بدقة، بما في ذلك سلوكها عند العتبة، والانهيار عند الاصطدام، وفترة الاستجابة.

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

Journal: Biomechanics and Modeling in Mechanobiology, Volume: 25, Issue: 1
DOI: https://doi.org/10.1007/s10237-025-02030-w
PMID: https://pubmed.ncbi.nlm.nih.gov/41504943
Publication Date: 2026-01-08
Author(s): Kert Tamm et al.
Primary Topic: Ion channel regulation and function

Overview

In 2025, the authors undertook a significant revision of their initial manuscript, resulting in a bifurcation into two distinct parts. The first part concentrates on modeling action potential (AP) propagation in unmyelinated axons, while the second part aims to incorporate the effects of myelination. The revision was prompted by feedback indicating that the original presentation was confusing for readers with experimental backgrounds in nerve processes. Consequently, the authors re-evaluated their numerical simulations, aligning parameters more closely with those used by Hodgkin and Huxley in their seminal 1952 paper, which introduced the classical Hodgkin-Huxley (HH) model.

The revised model posits that AP propagation is influenced not only by the HH ion mechanisms but also by the axon’s capacitance and inductance. The authors specifically modified the selection logic for the specific inductance value in their numerical examples, now choosing it to match the experimentally observed AP propagation velocity reported by Hodgkin and Huxley, rather than aiming for a predetermined velocity based on a fixed axon radius. This foundational work sets the stage for further modifications to the Lieberstein model to account for the effects of myelination on AP dynamics, acknowledging the significant impact of myelin sheaths on propagation velocity as revealed by experimental studies.

Introduction

The introduction of the paper discusses the Hodgkin-Huxley (HH) model, which effectively describes the dynamics of action potentials (AP) in unmyelinated axons by accounting for sodium and potassium ion currents. While the HH model provides a detailed representation of ion current physics, it simplifies the propagation of AP by neglecting inductance, relying instead on a diffusion-type cable equation. The authors highlight that advancements in computational resources since the model’s inception in the 1950s allow for more complex modeling approaches. Lieberstein’s alternative model incorporates inductance, enabling a more accurate representation of propagating AP signals without the additional assumptions required by the HH model.

The paper also addresses the differences in AP propagation between myelinated and unmyelinated axons, noting that myelinated axons exhibit higher AP velocities due to the insulating properties of the myelin sheath. The authors propose a phenomenological model that integrates the effects of the myelin sheath on AP propagation, emphasizing the importance of structural parameters such as the myelin length ratio (µ-ratio) and the g-ratio. By incorporating inductance into the governing equations, the model aims to explain the nonlinear changes in AP velocity, supported by experimental findings. The subsequent sections of the paper will detail the model for unmyelinated axons, present experimental results on myelinated axons, and ultimately propose a comprehensive model for AP propagation influenced by myelin sheath structure.

Discussion

The discussion section of the research paper elaborates on the Lieberstein model for action potential (AP) propagation in unmyelinated axons, emphasizing the significance of including inductance alongside the Hodgkin-Huxley (HH) ion mechanism. The Lieberstein model modifies the classical HH framework by incorporating inductance, resulting in a hyperbolic governing equation that ensures causal signal propagation, contrasting with the parabolic nature of the traditional cable equation which allows for non-causal infinite velocity. The model’s hyperbolicity is crucial for accurately describing the dynamics of APs, including their threshold behavior, annihilation upon collision, and the refractory period.

Numerical simulations based on the Lieberstein model demonstrate its ability to replicate key physiological characteristics of APs, such as the existence of a threshold for signal generation and the rapid dissipation of sub-threshold disturbances. The results align with previous findings and validate the model’s applicability for future modifications to include myelinated axons. The paper also outlines the structure of the model, the numerical methods employed, and the parameters used for simulations, ultimately reinforcing the model’s relevance in understanding nerve pulse propagation and its potential for further exploration in neurophysiology.