КАРТОЧКА ПРОЕКТА ФУНДАМЕНТАЛЬНЫХ И ПОИСКОВЫХ НАУЧНЫХ ИССЛЕДОВАНИЙ,
ПОДДЕРЖАННОГО РОССИЙСКИМ НАУЧНЫМ ФОНДОМ
Информация подготовлена на основании данных из Информационно-аналитической системы РНФ, содержательная часть представлена в авторской редакции. Все права принадлежат авторам, использование или перепечатка материалов допустима только с предварительного согласия авторов.
ОБЩИЕ СВЕДЕНИЯ
Номер проекта 25-25-00219
НазваниеПостроение математической модели физиологических состояний таламокортикальной системы головного мозга
Руководитель Сысоева Марина Вячеславовна, Доктор физико-математических наук
Организация финансирования, регион федеральное государственное автономное образовательное учреждение высшего образования "Санкт-Петербургский политехнический университет Петра Великого" , г Санкт-Петербург
Конкурс №102 - Конкурс 2025 года «Проведение фундаментальных научных исследований и поисковых научных исследований малыми отдельными научными группами»
Область знания, основной код классификатора 05 - Фундаментальные исследования для медицины; 05-106 - Нейробиология
Ключевые слова Математическое моделирование, анализ сигналов, анализ связанности, реконструкция моделей по временным рядам, состояния мозга, сон, бодрствование, эпилепсия, таламокортикальная система
Код ГРНТИ29.03.77
ИНФОРМАЦИЯ ИЗ ЗАЯВКИ
Ожидаемые результаты
По результатам проекта будут построены математические модели в виде систем обыкновенных дифференциальных уравнений без шума и с шумом (динамические и стохастические) сигналов локальных потенциалов мозга, соответствующих трём физиологическим состояниям мозга: пассивному бодрствованию, поверхностному и глубокому медленному сну. Будут использованы два подхода к моделированию: 1) по экспериментальным сигналам (реконструкция моделей по временным рядам), то есть модели главным образом феноменологические, в первую очередь отражающие особенности наблюдаемой активности, а не биологические принципы её генерации, и 2) из первых принципов, то есть на основе представлений биофизики и физиологии о строении таламокортикальной системы, в том числе о типах клеток, об архитектуре их связей, о строении синапсов и т.п.
Хотя феноменологические модели не претендуют на описание процессов в таламокортикальной системе с точки зрения биологии, построение даже такого оператора эволюции имеет огромные преимущества перед подходами, основанными на спектральном, вейвлетном анализе или моделировании формы сигнала, а именно эти подходы используются в настоящее время для решения задач разметки типов состояний по ЭЭГ, МЭГ, сигналам локальных потенциалов. Во-первых, феноменологические модели позволяют определять момент изменения оператора эволюции, когда сама динамика ещё не успела перестроиться (в терминах нелинейной динамики этому соответствует затягивание петли устойчивости), что важно для понимания и предсказания переходов между состояниями. Во-вторых, модели позволяют количественно охарактеризовать активность по значениям параметров и форме нелинейных функций, то есть дают измеримый критерий для диагностирования типов активности. Наконец, они могут быть использованы для диагностики связанности, как это было показано не только с использованием эмпирических моделей в виде отображений последования [Sysoeva M.V., Sitnikova E., Sysoev I.V., Bezruchko B.P., van Luijtelaar G. Application of adaptive nonlinear Granger causality: Disclosing network changes before and after absence seizure onset in a genetic rat model. Journal of Neuroscience Methods, 2014; 226: 33–41. DOI: 10.1016/j.jneumeth.2014.01.028], но и более сложных моделей в виде дифференциальных уравнений [Sysoev I.V. & Bezruchko B. P. Noise robust approach to reconstruction of van der Pol-like oscillators and its application to Granger causality. Chaos, 2021; 31: 083118. DOI: 10.1063/5.0056901].
Во второй части проекта будут построены модели из первых принципов, представляющие собою классические мезомасштабные модели, где большие группы нейронов одного типа и близкие по параметрам моделируются одной клеткой, но при этом каждая структура мозга представлена не одним или несколькими уравнениями для ионных токов, как в моделях нейронных масс, а ансамблем осцилляторов, что позволяет естественным образом учесть связи внутри структуры мозга и смоделировать популяционную вариабельность (изменчивость структуры связей и итоговой динамики в популяции различных особей одного происхождения, например, крыс одной и той же генетической линии). Основная идея — построить модели таким образом, чтобы все три рассматриваемых физиологических состояния были режимами одной и той же модели при различных параметрах или коэффициентах связи между узлами. Таким образом будет сделан важный шаг от моделирования отдельных функций (или патологических состояний) подсистем мозга к моделированию подсистемы в целом. При этом мы также намерены воспроизвести популяционную вариабельность — в разных моделях, отличающихся деталями матрицы связей при общих основных принципах организации будут варьировать распределения длительностей фаз сна и бодрствования, значения параметров, при которых они достигаются, вероятности переходов.
Итого, проект направлен на то, чтобы на новом уровне достоверности и новыми средствами решить важные проблемы обработки и воспроизведения сигналов электромагнитной активности мозга для целей диагностики различных состояний и прогнозирования переходов между ними.
ОТЧЁТНЫЕ МАТЕРИАЛЫ
Аннотация результатов, полученных в 2025 году
В результате ручной и автоматической разметки (в том числе специально адаптированными методами) большого объёма многоканальных экспериментальных записей локальных потенциалов из нескольких структур коры и таламуса получен большой набор референсных сигналов для моделирования тета- и дельта-ритмов в таламокортикальной системе. Используя ранее разработанные подходы к реконструкции маломерных динамических систем и их ансамблей, по наблюдаемым сигналам реконструированы модели в виде обобщённых осцилляторов ван дер Поля, в том числе связанных. Для этого были использованы, во-первых, временные ряды радиотехнических генераторов (модельных нейронов), во-вторых ряды тета-ритма, выделенные из записей крыс линии Wistar.
Используя математические модели нейронов в виде уравнений Ходжкина-Хаксли были построены математические модели генерации дельта- и тета-ритма. В зависимости от типа клетки для четырёх типов нейронов: пирамид и интернейронов коры, таламокортикальных клеток и клеток ретикулярного ядра таламуса уравнения были дополнены соответствующими токами, связи между уравнениями были организованы с помощью дополнительных дифференциальных уравнений. Все параметры моделей отдельных нейронов и синапсов были взяты из литературы, где они измерены в экспериментах. Для генерации ритмов в модели из небольшого числа нейронов была привлечена ключевая идея о преобразовании сетью внешнего высокочастотного сигнала гамма-диапазона (в нашем случае для простоты — сигнала нейрона-пейсмекера из коры больших полушарий). Генерация моделей на основе случайных матриц связей, построенных в соответствии с анатомическими особенностями таламокортикальной системы, а затем отбор моделей по принципу соответствия спектра и режимов колебаний позволили получить модели, качественно воспроизводящие высокочастотный дельта-ритм. Построение моделей в виде многосвязных цепочек нейронов — последовательных преобразователей частоты и/или передатчиков сигнала позволили получить как основную частоту дельта-ритма, так и тета-ритма. При этом важно, что для соответствующего преобразования частоты необходимы оказались нейроны всех четырёх рассмотренных типов. Это означает, что в рамках использованного модельного приближения, основанного на экспериментально измеренных параметрах и специально адаптированных уравнениях, тета- и дельта-ритмы могут возникнуть только как результат работы таламокортикальной системы в целом, а не в пределах какой-то её отдельной области.
Публикации
1.
Капустников А.А., Сысоев И.В., Сысоева М.В.
Network motifs in talamocortical system model responsible for delta-rhythm generation
Chaos, Solitons & Fractals, Kapustnikov A.A., Sysoev I.V., Sysoeva M.V. Network motifs in talamocortical system model responsible for delta-rhythm generation. Chaos, Solitons & Fractals. 2025. V.200. 117081. (год публикации - 2025)
https://doi.org/10.1016/j.chaos.2025.117081
2. Капустников А.А., Сысоева М.В., Сысоев И.В. Модель частотного преобразователя в таламокортикальной системе Молодёжная петербургская школа-конференция инженеров-педагогов : Сборник материалов научной конференции, Капустников А.А., Сысоева М.В., Сысоев И.В. Модель частотного преобразователя в таламокортикальной системе // Молодёжная петербургская школа-конференция инженеров-педагогов : Сборник материалов научной конференции, Санкт-Петербург, 22-23 апреля 2025 года. С. 34-35. (год публикации - 2025)
3. Долинина А.Ю. Исследование структуры локальных потенциалов мозга у крысмоделей в состоянии бодрствования Молодёжная петербургская школа-конференция инженеров-педагогов : Сборник материалов научной конференции, Долинина А.Ю. Исследование структуры локальных потенциалов мозга у крысмоделей в состоянии бодрствованиясистеме // Молодёжная петербургская школа-конференция инженеров-педагогов : Сборник материалов научной конференции, Санкт-Петербург, 22-23 апреля 2025 года. С. 25-27. (год публикации - 2025)
4. Капустников А.А., Сысоев И.В., Сысоева М.В. Преобразование частот с помощью простых сетевых мотивов таламо-кортикальной системы Нелинейные дни в Саратове для молодых - 2025 : Материалы XXXI Всероссийской научной конференции, Капустников А.А., Сысоев И.В., Сысоева М.В. Преобразование частот с помощью простых сетевых мотивов таламо-кортикальной системы // Нелинейные дни в Саратове для молодых - 2025 : Материалы XXXI Всероссийской научной конференции, Саратов, 26-30 мая 2025 года. С. 61-62. (год публикации - 2025)