<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE root>
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:ali="http://www.niso.org/schemas/ali/1.0/" article-type="oration" dtd-version="1.2" xml:lang="en"><front><journal-meta><journal-id journal-id-type="publisher-id">Genes &amp; Cells</journal-id><journal-title-group><journal-title xml:lang="en">Genes &amp; Cells</journal-title><trans-title-group xml:lang="ru"><trans-title>Гены и Клетки</trans-title></trans-title-group><trans-title-group xml:lang="zh"><trans-title>Genes and Cells</trans-title></trans-title-group></journal-title-group><issn publication-format="print">2313-1829</issn><issn publication-format="electronic">2500-2562</issn><publisher><publisher-name xml:lang="en">Human Stem Cells Institute</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="publisher-id">623481</article-id><article-id pub-id-type="doi">10.17816/gc623481</article-id><article-categories><subj-group subj-group-type="toc-heading" xml:lang="en"><subject>Conference proceedings</subject></subj-group><subj-group subj-group-type="toc-heading" xml:lang="ru"><subject>Материалы конференции</subject></subj-group><subj-group subj-group-type="article-type"><subject>Conference Report, Theses of Report</subject></subj-group></article-categories><title-group><article-title xml:lang="en">Bursting activity in the reduced mean-field model of neuron-glial interaction</article-title><trans-title-group xml:lang="ru"><trans-title>Пачечная активность в редуцированной среднеполевой модели нейрон-глиального взаимодействия</trans-title></trans-title-group></title-group><contrib-group><contrib contrib-type="author"><name-alternatives><name xml:lang="en"><surname>Olenin</surname><given-names>S. M.</given-names></name><name xml:lang="ru"><surname>Оленин</surname><given-names>С. М.</given-names></name></name-alternatives><address><country country="RU">Russian Federation</country></address><email>tatiana.levanova@itmm.unn.ru</email><xref ref-type="aff" rid="aff1"/></contrib><contrib contrib-type="author"><name-alternatives><name xml:lang="en"><surname>Levanova</surname><given-names>T. A.</given-names></name><name xml:lang="ru"><surname>Леванова</surname><given-names>Т. А.</given-names></name></name-alternatives><address><country country="RU">Russian Federation</country></address><email>tatiana.levanova@itmm.unn.ru</email><xref ref-type="aff" rid="aff1"/></contrib><contrib contrib-type="author"><name-alternatives><name xml:lang="en"><surname>Stasenko</surname><given-names>S. V.</given-names></name><name xml:lang="ru"><surname>Стасенко</surname><given-names>С. В.</given-names></name></name-alternatives><address><country country="RU">Russian Federation</country></address><email>tatiana.levanova@itmm.unn.ru</email><xref ref-type="aff" rid="aff1"/></contrib></contrib-group><aff-alternatives id="aff1"><aff><institution xml:lang="en">National Research Lobachevsky State University of Nizhny Novgorod</institution></aff><aff><institution xml:lang="ru">Национальный исследовательский Нижегородский государственный университет имени Н.И. Лобачевского</institution></aff></aff-alternatives><pub-date date-type="pub" iso-8601-date="2023-12-15" publication-format="electronic"><day>15</day><month>12</month><year>2023</year></pub-date><volume>18</volume><issue>4</issue><issue-title xml:lang="en"/><issue-title xml:lang="ru"/><fpage>870</fpage><lpage>873</lpage><history><date date-type="received" iso-8601-date="2023-11-16"><day>16</day><month>11</month><year>2023</year></date><date date-type="accepted" iso-8601-date="2023-11-18"><day>18</day><month>11</month><year>2023</year></date></history><permissions><copyright-statement xml:lang="en">Copyright ©; 2023, Eco-Vector</copyright-statement><copyright-statement xml:lang="ru">Copyright ©; 2023, Эко-Вектор</copyright-statement><copyright-year>2023</copyright-year><copyright-holder xml:lang="en">Eco-Vector</copyright-holder><copyright-holder xml:lang="ru">Эко-Вектор</copyright-holder><ali:free_to_read xmlns:ali="http://www.niso.org/schemas/ali/1.0/" start_date="2027-02-20"/><license><ali:license_ref xmlns:ali="http://www.niso.org/schemas/ali/1.0/">https://creativecommons.org/licenses/by-nc-nd/4.0/</ali:license_ref></license></permissions><self-uri xlink:href="https://genescells.ru/2313-1829/article/view/623481">https://genescells.ru/2313-1829/article/view/623481</self-uri><abstract xml:lang="en"><p>The study of neural activity synchronization in the brain is a central focus of modern neurobiology and neurodynamics. In a healthy brain, cognitive functions necessitate the accurate integration of neural activity at specific spatiotemporal scales. Some neurological and psychiatric disorders exhibit observable variations in patterns of synchronous population activity. One of the most significant and compelling synchronous population activity patterns is the bursting activity that participates in various informational and physiological processes, including epilepsy. Several mathematical models have been proposed to explain the mechanisms of formation of the bursting activity. The most intriguing and scientifically sound models enable incorporation of astrocytic modulation of neural activity [1].</p> <p>In this paper, we present a novel phenomenological model for accurately replicating the bursting neural population activity. Our model builds upon the Tsodyks-Markram [2] model and incorporates vital aspects of neuron-glial interaction through a tripartite synapse. The presented model is a simplified version of the previously proposed model [3]. During astrocyte activity (which lasts seconds), we can set the probability of neurotransmitter release, <italic>u</italic>, which is determined in fractions of a second. This enables the proposed model to be written in the following manner:</p> <p><italic>r</italic>=–<italic>E</italic>+<italic>α</italic> ln(1+exp((<italic>Ju</italic>(<italic>Y</italic>)<italic>xE</italic>+<italic>I</italic><sub>0</sub>)/<italic>α</italic>)),</p> <p>=(1–<italic>x</italic>)/<italic>τ</italic><italic><sub>D</sub></italic>–<italic>u</italic>(<italic>Y</italic>)<italic>xE</italic>,</p> <p>=–<italic>y</italic>/<italic>τ</italic><italic><sub>Y</sub></italic>+<italic>βσ</italic><italic><sub>y</sub></italic>(<italic>X</italic>).</p> <p>Here, <italic>E</italic>(<italic>t</italic>) represents the average neuronal activity of the excitatory population, while <italic>x</italic>(<italic>t</italic>) denotes the quantity of available neurotransmitter. Additionally, <italic>y</italic>(<italic>t</italic>) describes the concentration of gliotransmitters, which are released as a result of biochemical reactions during neuron-astrocyte interactions. The function describes the alteration in the likelihood of neurotransmitter release when gliotransmitters are present</p> <p><italic>u</italic>(<italic>y</italic>)=<italic>u</italic><sub>0</sub>+(∆<italic>u</italic><sub>0</sub>)/(1+exp(–50(<italic>y</italic>–<italic>y</italic><italic><sub>thr</sub></italic>)).</p> <p>Here <italic>u</italic><sub>0</sub> represents the probability of neurotransmitter release without astrocytic influence; ∆<italic>u</italic><sub>0</sub> indicates the change in the probability of neurotransmitter release caused by the interaction of the gliotransmitter with the presynaptic terminal, and <italic>y</italic><italic><sub>thr</sub></italic> represents the threshold value that determines the change in the probability of neurotransmitter release due to the effect of a gliotransmitter. The function describes the influence of neurotransmitter on the concentration of the gliotransmitter</p> <p><italic>σ</italic><italic><sub>y</sub></italic>(<italic>x</italic>)=1/(1+exp(–20(<italic>x</italic>–<italic>x</italic><italic><sub>thr</sub></italic>)),</p> <p>where <italic>x</italic><italic><sub>thr</sub></italic> is the astrocyte activation threshold.</p> <p>Note that the proposed model does not account for synaptic depression’s mechanism. Therefore, various dynamical regimes’ formation in the model is solely controlled by astrocytic dynamics.</p> <p>In the study, control parameters <italic>I</italic><sub>0</sub> and <italic>u</italic><sub>0</sub> were selected while the remaining parameters were held constant. Specifically, <italic>τ</italic>=0.013, <italic>τ</italic><italic><sub>D</sub></italic>=0.08, <italic>α</italic>=1.58, <italic>J</italic>=3.07. Changes in the concentration of neurotransmitters and gliotransmitters were characterized by parameters ∆<italic>u</italic><sub>0</sub>=0.305, <italic>τ</italic><italic><sub>ν</sub></italic>=3.3, <italic>β</italic>=0.3, <italic>x</italic><italic><sub>xthr</sub></italic>=0.75, <italic>y</italic><italic><sub>thr</sub></italic>=0.4. The research was conducted using computer modeling and numerical methods in nonlinear dynamics.</p> <p>The proposed model exhibits a diverse range of population activity patterns, such as spiking and bursting regular and chaotic activity. Chaotic dynamics zones are identified in the control parameter plane. An account is given of the mechanisms behind the emergence of these activity patterns through bifurcation analysis. It is shown that the occurrence of chaotic activity in the system can be associated both with the cascade of period doubling bifurcations and with the further development of chaos, as a result of which a homoclinic attractor appears in the system according to the Shilnikov scenario. It was also shown that multistability is observed in the system in a certain range of parameters.</p> <p>Note that the emergence of multistability and bursting activity in the model is independent of the intricacy of neuron and glial cell dynamics. These dynamics are instead influenced by the feedback loop between the presynaptic neuron and the glial cell. The effects of neuron-like dynamics and neural-glial interaction demonstrated are generally applicable, without specifying particular characteristics of neural-glial interaction, network architecture, or single-cell neuron dynamics. It should be noted that our proposed model solely concerns the potentiation of synaptic transmission by astrocytes.</p> <p>In summary, the proposed phenomenological model for population activity can replicate various patterns of neuron activity found in a wide range of dynamic memory and information processing studies. This research can potentially lead to the development of new, effective treatment methods for neurological diseases that involve neuron-glial interaction. Another potential application for these results is in the development of an optimized live chip that has preset functions. This requires a deeper understanding of rhythmogenesis in neural networks and brain function.</p></abstract><trans-abstract xml:lang="ru"><p>Исследование синхронизации нейронной активности в мозге является одним из ключевых направлений современной нейробиологии и нейродинамики. В здоровом мозге когнитивные процессы требуют точной интеграции нейронной активности в определённых пространственно-временных масштабах. Изменение синхронизации можно наблюдать при некоторых неврологических и психических расстройствах. Одним из наиболее интересных и важных паттернов синхронной популяционной активности является пачечная активность, которая может лежать в основе как нормальных физиологических, так и патологических процессов (к примеру, эпилепсии). К настоящему времени предложен ряд математических моделей, описывающих механизмы формирования пачечной активности. Наиболее интересные и биологически правдоподобные модели позволяют учитывать астроцитарную модуляцию нейронной активности [1].</p> <p>В данной работе мы предлагаем новую феноменологическую модель, которая позволяет воспроизводить пачечную популяционную активность нейронов. Предлагаемая модель основана на модели Цодыкса–Маркрама [2] и учитывает особенности взаимодействия нейронов и астроцитов через трёхчастный синапс. Модель является упрощением ранее предложенной модели [3]. На временах активности астроцита (порядка секунд) мы можем зафиксировать значение вероятности высвобождения нейротрансмиттера, <italic>u</italic>, динамика которой определяется на временах порядка миллисекунд. Это позволяет нам записать предлагаемую модель в следующем виде:</p> <p><italic>r</italic>=–<italic>E</italic>+<italic>α</italic> ln(1+exp((<italic>Ju</italic>(<italic>Y</italic>)<italic>xE</italic>+<italic>I</italic><sub>0</sub>)/<italic>α</italic>)),</p> <p>=(1–<italic>x</italic>)/<italic>τ</italic><italic><sub>D</sub></italic>–<italic>u</italic>(<italic>Y</italic>)<italic>xE</italic>,</p> <p>=–<italic>y</italic>/<italic>τ</italic><italic><sub>Y</sub></italic>+<italic>βσ</italic><italic><sub>y</sub></italic>(<italic>X</italic>).</p> <p>Здесь <italic>E</italic>(<italic>t</italic>) — средняя нейронная активность возбуждающей популяции. Переменная <italic>x</italic>(<italic>t</italic>) моделирует количество доступного нейротрансмиттера, а <italic>y</italic>(<italic>t</italic>) — описывает концентрацию глиотрансмиттера, высвобождаемого в результате биохимических реакций при нейрон-астроцитарном взаимодействии. Изменение вероятности высвобождения нейротрансмиттера в присутствии глиотрансмиттера описывается функцией</p> <p><italic>u</italic>(<italic>y</italic>)=<italic>u</italic><sub>0</sub>+(∆<italic>u</italic><sub>0</sub>)/(1+exp(–50(<italic>y</italic>–<italic>y</italic><italic><sub>thr</sub></italic>)).</p> <p>Здесь <italic>u</italic><sub>0</sub> — вероятность высвобождения нейротрансмиттера в отсутствие астроцитарного влияния; ∆<italic>u</italic><sub>0</sub> — изменение вероятности высвобождения нейротрансмиттера за счёт действия глиотрансмиттера на пресинаптическую терминаль, <italic>y</italic><italic><sub>thr</sub></italic> — пороговое значение, определяющее изменение вероятности высвобождения нейротрансмиттера за счёт воздействия глиотрансмиттера. Влияние нейротрансмиттера на концетрацию глиотрансмиттера описывается функцией</p> <p><italic>σ</italic><italic><sub>y</sub></italic>(<italic>x</italic>)=1/(1+exp(–20(<italic>x</italic>–<italic>x</italic><italic><sub>thr</sub></italic>)),</p> <p>где <italic>x</italic><italic><sub>thr</sub></italic> — порог активации астроцитов.</p> <p>Отметим, что предлагаемая модель не включает механизм синаптической депрессии, и формирование различных динамических режимов в модели управляется исключительно астроцитарной динамикой.</p> <p>В работе в качестве управляющих параметров были выбраны параметры <italic>I</italic><sub>0</sub> и <italic>u</italic><sub>0</sub>. Значения остальных параметров были зафиксированы следующим образом: <italic>τ</italic>=0,013, <italic>τ</italic><italic><sub>D</sub></italic>=0,08, <italic>α</italic>=1,58, <italic>J</italic>=3,07. Параметры изменения концентрации нейротрансмиттера и глиотрансмиттера: ∆<italic>u</italic><sub>0</sub>=0,305, <italic>τ</italic><italic><sub>ν</sub></italic>=3,3, <italic>β</italic>=0,3, <italic>x</italic><italic><sub>xthr</sub></italic>=0,75, <italic>y</italic><italic><sub>thr</sub></italic>=0,4. Исследование проводится с использованием компьютерного моделирования и численных методов нелинейной динамики.</p> <p>Показано, что предложенная модель воспроизводит богатый набор паттернов популяционной активности, в том числе спайковые и пачечные регулярные и хаотические режимы. На плоскости управляющих параметров найдены области, в которых система демонстрирует хаотическую динамику. С помощью бифуркационного анализа объяснены механизмы возникновения указанных паттернов активности. Показано, что возникновение хаотической активности в системе может быть связано как со сценарием удвоения периода, так и с дальнейшим развитием хаоса, в результате чего в системе появляется гомоклинический аттрактор по сценарию Шильникова. Показано также, что в системе в определённом диапазоне параметров наблюдается мультистабильность.</p> <p>Отметим, что появление мультистабильности и пачечной активности в модели не зависит от сложности локальной динамики нейронов и глиальных клеток. Эти типы динамики определяются наличием петли обратной связи между пресинаптическим нейроном и глиальной клеткой. Продемонстрированные эффекты нейроноподобной динамики и нейроно-глиального взаимодействия являются достаточно общими, поскольку они не подразумевают специфических характеристик нейрон-глиального взаимодействия, конкретной архитектуры нейронной сети или динамики отдельных нейронов. Следует отметить, что в предложенной модели мы рассматриваем только усиление (потенциацию) синаптической передачи астроцитом.</p> <p>Подводя итог, можно сказать, что предложенная феноменологическая модель популяционной активности может использоваться для воспроизведения различных паттернов популяционной активности нейронов в широком диапазоне исследований динамической памяти и обработки информации. Одним из возможных приложений таких исследований является разработка новых эффективных методов лечения неврологических заболеваний, связанных с нейрон-глиальным взаимодействием. Другая область, где эти результаты могут быть полезны, касается создания эффективного живого чипа с заранее заданными функциями, что требует лучшего понимания ритмогенеза в нейронных сетях и функционирования мозга.</p></trans-abstract><kwd-group xml:lang="en"><kwd>actrocyte</kwd><kwd>tripartite synapse</kwd><kwd>population activity</kwd><kwd>dynamical chaos</kwd><kwd>short-term plasticity</kwd></kwd-group><kwd-group xml:lang="ru"><kwd>астроциты</kwd><kwd>трёхчастный синапс</kwd><kwd>популяционная активность</kwd><kwd>динамический хаос</kwd><kwd>кратковременная пластичность</kwd></kwd-group><funding-group><funding-statement xml:lang="en">The study was supported by RSF grant No. 19-72-10128</funding-statement><funding-statement xml:lang="ru">Работа поддержана грантом Российского научного фонда 19-72-10128</funding-statement></funding-group></article-meta></front><body></body><back><ref-list><ref id="B1"><label>1.</label><citation-alternatives><mixed-citation xml:lang="en">Lazarevich IA, Stasenko SV, Kazantsev VB. Synaptic multistability and network synchronization induced by the neuron–glial interaction in the brain. Jetp Lett. 2017;105:210–213. doi: 10.1134/S0021364017030092</mixed-citation><mixed-citation xml:lang="ru">Lazarevich I.A., Stasenko S.V., Kazantsev V.B. Synaptic multistability and network synchronization induced by the neuron–glial interaction in the brain // Jetp Lett. 2017. Vol. 105. P. 210–213. doi: 10.1134/S0021364017030092</mixed-citation></citation-alternatives></ref><ref id="B2"><label>2.</label><citation-alternatives><mixed-citation xml:lang="en">Tsodyks M, Pawelzik K, Markram H. Neural networks with dynamic synapses. Neural Comput. 1998;10(4):821–835. doi: 10.1162/089976698300017502</mixed-citation><mixed-citation xml:lang="ru">Tsodyks M., Pawelzik K., Markram H. Neural networks with dynamic synapses // Neural Comput. 1998. Vol. 10, N 4. P. 821. doi: 10.1162/089976698300017502</mixed-citation></citation-alternatives></ref><ref id="B3"><label>3.</label><citation-alternatives><mixed-citation xml:lang="en">Barabash N, Levanova T, Stasenko S. Rhythmogenesis in the mean field model of the neuron-glial network. Eur Phys J Spec Top. 2023;232:529–534. doi: 10.1140/epjs/s11734-023-00778-9</mixed-citation><mixed-citation xml:lang="ru">Barabash N., Levanova T., Stasenko S. Rhythmogenesis in the mean field model of the neuron-glial network // Eur Phys J Spec Top. 2023. Vol. 232. P. 529–534. doi: 10.1140/epjs/s11734-023-00778-9</mixed-citation></citation-alternatives></ref></ref-list></back></article>
