تذبذبات متعددة الاستقرار في الشبكات القشرية مع فئتين من التثبيط Multi-stable oscillations in cortical networks with two classes of inhibition

شارك:
المجلة: PLoS Computational Biology، المجلد: 22، العدد: 6
DOI: https://doi.org/10.1371/journal.pcbi.1014391
PMID: https://pubmed.ncbi.nlm.nih.gov/42263113
تاريخ النشر: 2026-06-09
المؤلف: Arnab Dey Sarkar وآخرون
الموضوع الرئيسي: ديناميات الأعصاب ووظيفة الدماغ

نظرة عامة

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

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

مقدمة

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

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

الطرق

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

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

النتائج

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

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

المناقشة

في هذه الدراسة، طور المؤلفون نموذج شبكة عصبية نابضة يتكون من مجموعات مثيرة (E) وثلاث مجموعات مثبطة (I، S، VIP) للتحقيق في ظهور تذبذبات إيقاعية متميزة. تتكون كل مجموعة من 400 خلية عصبية من نوع التكامل والنار، مع اتصالات مشبكية تتميز بديناميات تتلاشى بشكل أسي. يكشف النموذج عن دورتين ثابتتين، تُسمى “كبيرة” (15 هرتز) و”صغيرة” (16 هرتز)، تختلفان بشكل كبير في التماسك والسعة. تظهر الدورة الكبيرة تزامنًا أعلى بين المجموعات، كما تم قياسه بواسطة معامل ترتيب كوراموتو، بينما تظهر الدورة الصغيرة نشاطًا أضعف، خاصة في مجموعة S. تسلط الدراسة الضوء على تعددية الاستقرار في النظام، حيث يمكن أن تتواجد كلا الإيقاعات تحت نفس ظروف المعلمات، مما يبرز دور تهيئة الشبكة في تحديد الديناميات المرصودة.

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

Journal: PLoS Computational Biology, Volume: 22, Issue: 6
DOI: https://doi.org/10.1371/journal.pcbi.1014391
PMID: https://pubmed.ncbi.nlm.nih.gov/42263113
Publication Date: 2026-06-09
Author(s): Arnab Dey Sarkar et al.
Primary Topic: Neural dynamics and brain function

Overview

In this research, the authors challenge the classical perspective on cortical rhythms, which primarily attributes the generation of gamma- and beta-band oscillations to the interactions between excitatory pyramidal neurons (E) and inhibitory parvalbumin-expressing interneurons (I). They emphasize the significance of various inhibitory interneuron subtypes, particularly somatostatin-expressing neurons (S), which not only receive excitatory inputs from E cells but also inhibit both E and I populations. Additionally, the S cells are influenced by vasoactive intestinal peptide (VIP) neurons, which integrate inputs from other cortical regions.

To investigate these dynamics, the authors develop a spiking network model that incorporates populations of E, I, and S neurons, alongside VIP modulation. This model is simplified into a system of nine differential equations that capture the mean membrane potential, firing rate, and synaptic conductance for each neuronal population. Through this reduced model, the study identifies a broad spectrum of parameters that support multiple coexisting rhythms. Utilizing nonlinear dynamics tools, the authors analyze the contributions of the two classes of inhibition and VIP modulation, shedding light on their roles in shaping the characteristics of cortical oscillations.

Introduction

The introduction of the research paper discusses the role of neuronal oscillations in the cortex, emphasizing their distinct frequency bands and functions in cognitive processes. Gamma rhythms (30-80 Hz) are linked to local circuit synchronization and sensory processing, primarily mediated by parvalbumin (PV) interneurons. In contrast, beta rhythms (13-30 Hz) are divided into beta1 (13-20 Hz) and beta2 (20-30 Hz), with beta1 associated with long-range communication and motor control, while beta2 relates to local circuitry and inhibitory feedback. The interplay between different inhibitory interneuron classes, particularly PV and somatostatin (SOM) interneurons, is highlighted as crucial for generating and modulating these oscillations.

Recent studies have shown that the dynamics of excitatory (E) and inhibitory (I) neuron networks, particularly involving PV and SOM, are essential for understanding cortical rhythms. Computational models have demonstrated that beta1 and beta2 rhythms arise from distinct network dynamics, with PV enhancing beta2 through strong inhibition and SOM modulating phase coupling between the two bands. The paper aims to investigate the emergence of multi-rhythmicity in cortical networks using a mathematical model that incorporates these inhibitory populations. By employing a mean-field reduction and bifurcation analysis, the study seeks to elucidate how the interactions among E, PV, and SOM populations influence network oscillations, stability, and the coexistence of distinct rhythmic states, thereby contributing to a deeper understanding of cortical dynamics.

Methods

The “Materials and Methods” section of the research paper outlines the experimental design and procedures employed to investigate the research questions. It details the specific materials used, including any reagents, equipment, and biological samples, ensuring reproducibility of the study. The methodology encompasses the techniques applied for data collection and analysis, such as statistical tests, experimental protocols, and any computational models utilized.

Additionally, the section may describe the sampling methods, controls, and any ethical considerations taken into account during the research. By providing a comprehensive overview of the methods, this section aims to enable other researchers to replicate the study and validate the findings presented in subsequent sections of the paper.

Results

The “Results” section of the research paper presents key findings derived from the conducted experiments or analyses. It details the outcomes of the study, including any significant statistical results, observed trends, and comparisons with existing literature. The section may include visual representations such as graphs or tables to illustrate the data effectively.

The results highlight the main hypotheses tested and provide evidence supporting or refuting them. Specific metrics, such as p-values or confidence intervals, may be reported to quantify the significance of the findings. Overall, this section serves to communicate the empirical contributions of the research, laying the groundwork for subsequent discussion and interpretation.

Discussion

In this study, the authors developed a spiking neural network model comprising excitatory (E) and three inhibitory populations (I, S, VIP) to investigate the emergence of distinct rhythmic oscillations. Each population consists of 400 quadratic integrate-and-fire neurons, with synaptic connections characterized by exponentially decaying dynamics. The model reveals two stable limit cycles, termed “big” (15 Hz) and “small” (16 Hz), which differ significantly in coherence and amplitude. The big limit cycle exhibits higher synchronization among the populations, as quantified by the Kuramoto order parameter, while the small limit cycle shows weaker activity, particularly in the S population. The study highlights the multistability of the system, where both rhythms can coexist under the same parameter conditions, emphasizing the role of network initialization in determining the observed dynamics.

To further understand the underlying mechanisms of these oscillations, the authors constructed a mean-field model that simplifies the dynamics of the spiking network. This model captures the essential behavior of the network as the number of neurons approaches infinity, revealing bifurcations that dictate the stability of oscillatory states. The analysis indicates that varying the excitatory drive ($\mu_e$) and the inhibitory influences ($\gamma$ and $\lambda$) can lead to transitions between different rhythmic states, including bi-rhythmicity. The findings suggest that the interplay between distinct inhibitory populations is crucial for generating multiple oscillatory frequencies, which may have implications for understanding neural dynamics in biological systems.