Синхронизация (нейробиология)

В нейробиологии синхронизацией (от греч. συνχρόνος — одновременный) называют динамический режим, который характеризуется периодической одновременной активацией определенной популяции нейронов, или синхронизацию между локальными колебаниями двух или нескольких популяций нейронов.

Синхронизированная активность больших популяциях нейронов является главным механизмом образования макроколебаний, или так называемых ритмов головного мозга человека. Известно, что ритмы играют важную роль в процессах передачи и обработки информации в центральной нервной системе (сон, память, внимание, пространственная навигация и т. д.), однако общей теории относительно точного назначения и принципов работы ритмов пока нет.

С другой стороны, чрезмерная локальная синхронизация в некоторых сетях имеет патологические проявления, такие как тремор в болезни Паркинсона, эпилептические припадки, а также некоторые расстройства высшей мозговой функции (шизофрения, аутизм, и т. д.).

Важная роль синхронизации в процессах обработки информации в мозге делает её очень привлекательной темой исследований в медицинской и вычислительной нейронауке. Методы математического моделирования синхронизации (часто заимствованные из физики) в сочетании с доступностью вычислительных мощностей создают благоприятные условия для прогресса в направлении более глубокого и разностороннего понимания этого явления.

Вследствие того, что термин «синхронизация» является довольно распространённым и часто используется в различных отраслях науки, его точное определение может сильно варьировать в зависимости от конкретной сферы применения. В биологическом контексте удобно различать два вида синхронизации^[1]:

• локальная синхронизация, синхронизация активности отдельной популяции нейронов, которая приводит к образованию макроколебаний и может становиться заметной в сигнале ЭЭГ
• фазовая синхронизация, синхронизация фазы макроколебаний между удалёнными участками мозга или группами нейронов

Во всех секциях настоящей статьи, за исключением секции «Математическое моделирование», будем использовать именно вышеприведенное определение термина «фазовая синхронизация».

С другой стороны, в детальном математическом моделировании локальной синхронизации принято использовать термин «фазовая синхронизация» в более абстрактном смысле. А именно, два или несколько осцилляторов фазово синхронизированы, если соответствие их фаз не меняется со временем^[2], или, в наиболее общем определении, зависит по какой-либо функции. Причём амплитуды колебаний вышеупомянутых осцилляторов не обязательно должны быть одинаковы. Такими осцилляторами могут моделироваться как индивидуальные нейроны, так и сети или большие популяции нейронов. В секции «Математическое моделирование» описано несколько типов синхронизации.

Электроэнцефалография (ЭЭГ) — один из ключевых неинвазивных методов исследования головного мозга, который характеризуется высоким временным разрешением, но низкой пространственной разрешающей способностью, особенно в сравнении с МРТ или инвазивными электродными методами. Из-за низкого пространственного разрешения (обычно не более 100 каналов), полученный сигнал является результатом суперпозиции активности крупных популяций нейронов, а значит, является удобным методом для нахождения и исследования синхронизированных коллективных режимов в сетях головного мозга. Синхронизация достаточно большой популяции нейронов обычно производит глобальные колебания, которые и записываются электродами электроэнцефалографа.

В большинстве случаев сильная глобальная синхронизация больших популяций (синхронизация типа I) присуща состоянию бездействия или патологическому состоянию, поскольку динамика полностью синхронизированной сети недостаточно сложна для эффективной обработки информации. Во время нормальной работы (кроме состояния глубокого сна) небольшие локально-синхронизированные подсети генерируют колебания разных частот (синхронизация типа II), при этом глобальная синхронизация исчезает^[4]. Таким образом, спектр ЭЭГ сигнала содержит несколько важных составляющих, которые принято классифицировать по частоте, назначая каждому диапазону греческую букву. В таблице ниже приведено краткое описание каждого из ритмов согласно традиционной системе классификации, а также их соответствующие функции.

В следующих секциях более подробно рассмотрены свойства и пути образования некоторых ритмов.

Кроме вышеупомянутых канонических ритмов, сигнал ЭЭГ может также содержать высокочастотные компоненты с частотой свыше 400 Гц. В большинстве случаев такие компоненты имеют очень низкую энергию и часто игнорируются или воспринимаются как шум. Однако недавние исследования показали, что высокочастотные колебания в сигнале ЭЭГ, которые можно наблюдать в соматосенсорной коре головного мозга в ответ на стимуляцию периферийного нерва, могут быть результатом очень точной синхронизации нейронных всплесков в соответствующих сетях^[9]. В экспериментах с одновременным получением ЭЭГ и значения внешнеклеточного сигнала было показано, что форма высокочастотных компонентов сигнала ЭЭГ содержит информацию о временной структуре всплеска спайков индивидуальных нейронов. Таким образом, высокочастотные ритмы ЭЭГ могут пролить свет на динамику клеточного уровня без вмешательства в биологические ткани^[9].

Синхронизированная активность нейронов является базой для генерации периодических моторных команд для ритмических движений. Ритмичные команды, в свою очередь, производятся группой связанных нейронов, которые формируют сеть под названием центральный генератор упорядоченной активности (ЦГУА). При активации такие нейронные сети могут генерировать сложные ритмические моторные команды, даже при отсутствии сигналов обратной связи, которые несут информацию о конкретных временных интервалах. Ходьба, плавание или дыхание являются примерами работы ЦГУА^[10]. ЦГУА наиболее исследованы у низших животных, таких как минога, однако также существуют свидетельства о наличии ЦГУА у людей^[11].

Фазовая синхронизация (ФС) колебаний между различными зонами мозга имеет ряд важных функций, которые являются критическими для эффективной работы памяти. Важной функцией является улучшение передачи информации между удалёнными участками мозга во время ФС^[1]. Оптимальная передача информации может происходить благодаря ФС следующим образом: высокий уровень локального потенциала обычно свидетельствует о высоком уровне активности в локальной популяции в данный момент времени. В то же время чувствительность нейронов этой популяции также повышена, потому что мембраны нейронов в среднем более деполяризованы, чем когда локальный потенциал имеет низкие значения. Таким образом, если две удаленные популяции нейронов имеют синхронизированные коллективные колебания, то момент повышенной активности одной популяции совпадет с высокой возбудимостью другой, обеспечивая высокую вероятность передачи информации^[12].

ФС также может способствовать усилению синаптической пластичности. Синаптическая пластичность, которая зависит от времени импульса (англ. Spike-Timing Dependent Plasticity, STDP) влечёт долговременную потенциацию синапсов, если принимающий нейрон срабатывает не более чем через 10-20 мс после входного импульса. ФС может обеспечить достаточно узкие интервалы срабатывания нейронов между двумя группами и таким образом увеличить или уменьшить силу связи между ними^[1].

По распространённой системе классификации, существуют два основных вида памяти:

• Кратковременная память
• Долговременная память

В работе кратковременной памяти ФС наблюдается между префронтальной корой и височной долей мозга (на частоте θ-ритма) во время чтения, записи и удержания краткосрочных воспоминаний, а также между фронтальными и париетальными зонами коры (на γ и β-частотах) во время удержания воспоминаний. Запись и чтение визуальных объектов из декларативной памяти (тип долговременной памяти) сопровождаются фазовой синхронизацией колебаний между передними и задними зонами мозга на θ, γ, δ и β-частотах. Другие виды синхронизации также принимают участие в презентации объектов в памяти^[1].

Результаты экспериментов с одновременной записью электрической активности в гиппокампе и миндалине свидетельствуют об усилении синхронизации на θ-частотах между этими участками после формации и во время действия условного рефлекса страха^[13].

Во время глубокой фазы сна, а также при некоторых видах анестезии на энцефалограммах обычно наблюдаются колебания низкой частоты (0,3 — 1 Гц) и дельта-колебания. Такие ритмы часто ассоциируются с отсутствием сознания, а также часто используются как индикатор наличия сознания у пациента под наркозом. На клеточном уровне осцилляции выглядят как периодическая смена состояний активного срабатывания нейронов и тишины, каждый из которых длится до одной секунды^[14]. Одним из возможных механизмов генерации медленных колебаний (до 1 Гц) является совместная работа соответствующих сетей таламуса и неокортекса следующим образом:

• активность нейронов в коре инициирует внутренние осцилляции в таламусе
• колебания в таламусе усиливают синхронизацию тормозных интернейронов
• активация тормозных интернейронов подавляет активность нейронов кортекса, и наступает период тишины
• цикл повторяется

Другие механизмы образования низкочастотных колебаний хорошо освещены в работе, например, Г. Тимофеева^[14].

При нормальных условиях, глубокий сон играет важную роль в процессах консолидации памяти. Согласно популярной теории, медленные коллективные колебания регулируют ре-активацию и транспорт воспоминаний от гиппокампа к кортексу. Во время фазы быстрого сна (повышенная тета-активность и локальная синаптическая пластичность), ре-активированные воспоминания в кортексе закрепляются на синаптическом уровне, то есть откладываются в долговременную память^[15][16].

Внимание — направленность психической деятельности человека или животного и её сосредоточенность в определённый момент на объектах или явлениях, которые имеют определённое значение при одновременном абстрагировании от других, в результате чего они отображаются полнее, яснее и глубже, чем другие^[17]. Во время концентрации внимания на конкретном объекте или его свойстве, повышается скорость обработки информации, уменьшается время реакции, повышается уровень точности, чувствительности к небольшим изменениям в стимуле и восприятие контраста.

Во-первых, внимание усиливает активность нейронов, которые представляют объект внимания или определенную характеристику (цвет, ориентация). Так, например, нейроны зрительной коры, рецептивное поле которых содержит объект внимания, сильнее синхронизированы с локальным потенциалом (на γ-частоте), чем нейроны, которые реагируют на другие объекты вне внимания. Такой способ селективного усиления синхронизации предположительно является альтернативой повышению частоты срабатывания нейронов (она обычно не меняется в зависимости от направления внимания)^[18].

Во-вторых, внимание также регулирует информационное сообщение между отдалёнными частями мозга, усиливая связь между группами нейронов, которые несут важную для животного информацию. Многие исследования приходят к выводу что такое усиление связи происходит за счет селективной синхронизации между зонами коры^[19].

Во время нормальной работы мозга, локальная и межзонная синхронизации нейронной активности играют важную роль в таких ключевых процессах, как внимание, память, моторика и т. д. Однако важным является также и определенный баланс синхронизации и десинхронизации^[20]. Чрезмерная или недостаточная синхронизация может быть причиной многих патологических процессов, в числе которых эпилепсия, тремор, шизофрения, деменция и некоторые другие. Таким образом, важным направлением развития нейронауки является поиск эффективных путей обезвреживания патологической синхронизации и понимание принципов её создания.

Эпилепсия характеризуется внезапным началом синхронизации в относительно больших нейронных сетях, в результате чего становится невозможной эффективная обработка информации, что в свою очередь препятствует нормальной работе мозга. Эпилептические припадки могут сопровождаться кратковременной потерей сознания или более сложными сенсорными и моторными проявлениями^[20]. Наиболее проблематичным следствием эпилепсии является непредсказуемость припадков, вследствие чего людям, страдающим эпилепсией, может быть на законодательном уровне запрещено управлять транспортными средствами и выполнять некоторые другие занятия, требующие постоянного внимания.

Точные нейронные механизмы эпилептической синхронизации пока в совершенстве не известны. Однако применение теории динамических систем анализа и моделирования сигнала, который генерируется такими сетями, даёт определённые результаты в направлении понимания феномена. Любая нейронная сеть может быть описана как (нелинейная) динамическая система с одним или несколькими атракторами. (Аттрактор — это множество точек в фазовом пространстве системы, в котором сходятся траектории системы. Другими словами, это тип поведения, когда система возвращается в исходное состояние после небольшого отклонения в результате внешнего вмешательства. Простейшим примером аттрактора является предельный цикл. На таком атракторе система начинает описывать периодические колебания). Согласно одной из теорий, переход от нормальной к эпилептической активности сети может быть двух видов^[21]:

• В результате постепенной деформации аттрактора нормальной активности в аттрактор эпилептической активности (через соответствующую бифуркацию), вследствие изменения параметров системы (например, внеклеточной концентрации ионов кальция или магния), или внешнего стимула. Такая модель соответствует типам припадков, когда изменение состояния протекает постепенно, в несколько фаз.
• Вследствие резкого «перескакивания» с одного нормального аттрактора на аттрактор синхронизированной активности, что было вызвано случайной флуктуацией или внешним стимулом (например, вспышками света). В отличие от первого случая, здесь система одновременно имеет два или несколько аттракторов, включая нормальные и эпилептические (такую систему называют мультистабильной). В этом случае, в мозге здорового человека такие аттракторы могут быть расположены дальше друг от друга, что делает невозможным случайный переход в состояние эпилептического припадка).

В зависимости от типа эпилепсии, могут быть применены различные подходы к моделированию динамических процессов, приводящих к массовой синхронизации нейронов^[14][21][22].

Болезнь Паркинсона — это медленно прогрессирующее двигательное возбуждение, которое сопровождается нарушением речи, скованностью, повышением тонуса мышц, тремором и ригидностью мышц, замедлением движений. Тремор, который характеризуется ритмичными неконтролируемыми движениями конечностей, является очень проблематичным симптомом, так как, при определенных условиях, делает невозможным выполнение больным многих повседневных действий, таких как управление авто, удержание открытой ёмкости с жидкостью, рукописание и т. д.

Патологический процесс при этом заболевании характеризуется деградацией дофаминергических нейронов чёрного вещества, что влечёт за собой усиление активности и синхронизации нейронов в субталамических ядрах (часть базальных ганглий) и приводит к чрезмерному торможению таламо-кортикальной связи.

На настоящий момент точные механизмы генерации тремора у пациентов с болезнью Паркинсона неизвестны. Однако всё больше данных корреляционного характера свидетельствуют о связи между абнормально синхронизированной активностью в определённых частях мозга и мышечным тремором^[23][24]. Так, в соответствии с ведущей гипотезой, тремор генерируется группой нейронов в базальных ядрах, которые работают как центральные осцилляторы и продуцируют периодические импульсы определённым группам мышц. Причём важным является как общий уровень активности центральных осцилляторов, так и уровень их фазовой синхронизации^[23].

Шизофрения — это когнитивное расстройство, главными симптомами которого являются галлюцинации, сильно искажённое восприятие реальности и общая дезорганизация мышления и действий. Как и в случае с многими сложными ментальными заболеваниями, патофизиология шизофрении до сих пор остаётся до конца не известной^[25].

Как известно, осцилляторная активность сетей головного мозга связана с такими важными процессами, как селективность к внешним стимулам с помощью внимания, памяти и сознания. Так как симптомы шизофрении связаны именно с расстройствами вышеупомянутых процессов, есть основания считать, что анормально низкая β и γ-синхронизация, которая наблюдается у больных, является важным фактором в патофизиологии заболевания. Современные эмпирические данные свидетельствуют о нелокальном и системном характере дисфункции головного мозга, которая возникает при шизофрении и может включать в себя расстройство в динамике нейронных сетей в префронтальных и теменных зонах, а также ухудшение связей между отдалёнными зонами^[25][26].

Аутизм, заболевание, которое подавляет социальную и коммуникативную способности пациентов, также связано с уменьшенной осиляторной активностью на γ-частотах. Больные аутизмом страдают суженным фокусом внимания и не способны достаточно хорошо интегрировать фрагменты информации в единое целое, что может быть следствием недостаточно эффективной передачи информации между зонами мозга^[25][26].

Одним из распространённых типов фармацевтического подавления чрезмерной синхронизации является использование разнообразных блокаторов ионных каналов, вследствие чего уменьшается возбудимость целевых нейронов или же угнетаются синаптические связи. Например, увеличение внеклеточной концентрации Mg²⁺ может прекращать синхронизированную колебательную активность в моделях эпилепсии. С одной стороны, это происходит благодаря способности ионов магния блокировать НМДА-рецепторы в синаптических связях между нейронами. С другой стороны, магний и некоторые другие 2- и 3-валентные ионы (например, Ca²⁺) могут влиять на локальный заряд вокруг натриевых и калиевых потенциалчувствительных каналов, снижая активационный порог и некоторые динамические характеристики нервной клетки^[27][28].

Фармацевтических средства также могут действовать через связывание и деактивацию определённых нейротрансмиттеров.

Почти треть пациентов с эпилепсией не показывают положительной динамики в результате лечения антиэпилептическими медикаментами^[29]. Однако в некоторых случаях (около 7-8 %), удаётся достаточно точно локализовать источник анормальной эпилептической синхронизации в мозге, который удаляется с помощью оперативного вмешательства. Современные хирургические методы также позволяют проведение операций по удалению чрезмерно синхронизированных популяций в глубоком сером веществе (бледный шар, субталамические ядра) для устранения симптомов многих моторных заболеваний, включая болезнь Паркинсона^[30].

Развитие компьютерных технологий открыло новую эру в лечении пациентов, страдающих болезнью Паркинсона. Всё чаще применяется новый метод глубокой стимуляции мозга (ГСМ). Метод ГСМ заключается в использовании импульсов электрического стимулятора, встроенного поверхностно в организм пациента. Через введённый в мозг специальный электрод импульсы подаются в базальные ганглии, в частности, в субталамические ядра, и выполняют высокочастотную стимуляцию групп нейронов, что приводит к заметному улучшению перенесения болезни. Вследствие существенного ослабления симптомов болезни, пациенты могут вернуться к нормальной повседневной деятельности и значительно уменьшить объём употребления лекарств. Кроме болезни Паркинсона, показаниями к операции являются врожденные дистония и эссенциальный тремор^[31]. Также продолжаются клинические испытания для применения ГСМ в лечении депрессии.

Несмотря на большой успех ГСМ в лечении многих болезней, технология пока находится на раннем этапе развития, и полный отказ от хирургических вмешательств на пользу ГСМ является темой дебатов^[30].

Методы математического моделирования физических и химических процессов нашли широкое применение в биофизике и нейронауке. Начиная от моделирования стохастической динамики отдельных ионных каналов и заканчивая сетями со сложной неоднородной топологией и миллионами нейронов, модели помогают учёным понять механизмы функционирования нервной системы на различных уровнях абстракции, а значит и найти новые пути лечения соответствующих заболеваний, и разрабатывать более умные алгоритмы и вычислительные машины. Главным назначением математической модели является устранение незначительных факторов какого-то процесса или наблюдения и кристаллизация именно того минимального набора механизмов, который является достаточным для объяснения данных в конкретном контексте. Правильно построенная модель позволяет тестировать интересные гипотезы и делать важные предсказания относительно работы той или иной системы. В контексте синхронизации, моделирование позволяет применять аналитические и вычислительные методы для оценки динамических режимов сети. Например, бифуркационный анализ детальной модели нейронной сети может помочь определить максимально допустимые значения физиологических параметров, которые не приведут к переходу сети в патологически синхронизированный режим.

Биологический нейрон можно моделировать на различных уровнях абстракции, но суть многих моделей заключается в выделении квазипериодических свойств его динамики, то есть изображении нейрона как осциллятора определённого уровня сложности.

Модель соединённых осцилляторов Курамото^[32] — одна из наиболее абстрактных моделей, предназначенных для исследования синхронизации колебаний в нейронных сетях (осцилляторы имеют только одну степень свободы — фазу колебаний, игнорируя амплитуду)^[33]. Взаимодействия между такими осцилляторами описываются простой алгебраической формой (такой, как синус) и коллективно генерируют определённую картину глобальной активности. Модель Курамото и её расширения (такие как, например, добавление возможности формировать связи различной силы) широко используются для изучения колебательных процессов в мозге^[34]. В частности, её удобно использовать для изучения механизмов синхронизации в группе связанных нейронов.

В простейшем виде, фаза (${\displaystyle \theta _{i}}$) каждого из N нейронов задаётся следующим образом:

${\displaystyle {\dot {\theta _{i}}}=\omega _{i}+\sum _{j=1}^{N}\sigma _{ij}a_{ij}\sin(\theta _{j}-\theta _{i}),\qquad i=1\ldots N}$

где ${\displaystyle \omega _{i}}$ это частота колебания i-го осциллятора, ${\displaystyle \sigma _{ij}}$ и ${\displaystyle a_{ij}}$ — сила соединения и элемент матрицы соединений соответственно.

Благодаря простоте определения, можно относительно быстро симулировать динамику крупных популяций таких нейронов, а также находить аналитические решения в определённых случаях.

Спаечные модели нейронов используются в случаях, когда нельзя ограничиться только фазовой или частотной абстракцией нейрона, то есть когда необходимо учитывать амплитуду колебаний, и/или другие, более детальные аспекты работы биологического нейрона. Характерным свойством спаечных моделей является способность генерировать активность, приближённую к потенциалу действия (спайку). Такие сети бывают очень полезными, когда необходимо сделать предсказание относительно значений определённых биологических параметров, потому что в их математическое определение часто включены аналоги биологических параметров (например, динамика ионных каналов или концентрация определенных веществ во внеклеточной среде). Однако главным недостатком таких моделей является сложность и нелинейность их математического определения, что ограничивает возможности аналитических решений и требует существенных затрат компьютерного времени для вычислительных решений.

Обычно один такой нейрон определяется системой нелинейных дифференциальных уравнений. В вычислительной нейронауке наиболее популярными являются следующие модели:

• экспоненциальная импульсная модель (модель англ. integrate-and-fire: «интегрировать и возбуждаться», «интегрировать и сработать») (AdEx)^[35]. Является простейшей моделью, которая способна продуцировать все основные динамические режимы биологического нейрона, такие как однократное срабатывание, регулярное срабатывание в ответ на стимуляцию постоянным током, пачечная активность и т. д. Широко используется в исследованиях синхронизации у крупных гетерогенных популяций нейронов. Описывается системой двух дифференциальных уравнений.
• Модель Хиндмарша — Роуза (HR) — одна из самых простых непрерывных моделей, которая способна генерировать пачечную, а также хаотичную активности. Описывается системой трёх дифференциальных уравнений^[36].
• Модель Ходжкина — Хаксли (HH) и производные — семейство детальных биологически-реалистичных моделей нейронов, которые включают в себя кинетику натриевых, калиевых, кальциевых и других ионных каналов. Различные модификации способны очень точно имитировать форму потенциала действия настоящих нейронов. Система обычно описывается как минимум четырьмя нелинейными дифференциальными уравнениями. Несмотря на сложность математического определения, HH-подобные модели успешно используются для исследований эффектов синхронизации в больших популяциях нейронов^[37].

Существует большое количество различных определений явления синхронизации, которые могут быть применены в том или ином контексте. Отметим некоторые основные виды синхронизации

• Полная синхронизация — самый простой с математической точки зрения тип синхронизации, который характеризуется полной сходимостью соответствующих переменных всех нейронов в сети. Например, если i-й нейрон описывается вектором переменных ${\displaystyle \mathbf {x} _{i}(t)}$, то говорят что сеть из N нейронов полностью синхронизирована, если её состояние соответствует уравнению ${\displaystyle [\mathbf {x} _{1}(t)=\mathbf {x} _{2}(t)=\dots =\mathbf {x} _{N}(t)]}$. Состояние полной синхронизации может быть достигнуто только в сети с идентичными нейронами.
• Фазовая синхронизация накладывает определённые ограничения на фазовые отношения между осцилляторами, игнорируя отношение амплитуд. Если в системе можно определить понятие фазы, то фазовая синхронизация сети таких осцилляторов означает, что разность фаз между элементами сети должна или вообще не меняться со временем, или находиться в определённых конечных границах^[38]. Такой подход позволяет исследовать синхронизацию сетей с осцилляторами, параметры которых могут отличаться.
• Обобщенная синхронизация подобна фазовой, только отношение между фазами должны описываться определённой функцией.
• Пачечная синхронизация — синхронизация нейронов во время пачечной активности. Такая синхронизация не требует точного совпадения моментов спайков нейронов, но требует синхронизированного срабатывания пачек спайков (ритмических, ограниченных во времени всплесков активности с генерацией нескольких потенциалов действия)^[39].

Полная синхронизация

Фазовая синхронизация (с временной задержкой)

Пачечная синхронизация

Простота математического определения состояния полной синхронизации делает возможным применение интересных аналитических методов для исследования этого вида синхронизации. Глобальная функция стабильности синхронизированного состояния (англ. Master Stability Function, MSF), например, позволяет определять стабильность состояния полной синхронизации для сети идентичных нейронов^[40]. Другими словами, рассчитав MSF для сети нужных моделей нейронов, можно с точностью сказать, существуют ли такие топологии соединения нейронов, при которых данная сеть будет способна синхронизироваться (то есть будет иметь стабильное состояние полной синхронизации). Красота метода заключается в том, что MSF надо рассчитывать лишь один раз для каждого конкретного типа модели и набора параметров, а из результата можно делать выводы о стабильности синхронизированного состояния для любого количества нейронов и любой топологии (при соблюдении условия одинаковой суммарной силы входного сигнала для нейронов).

Основной процесс анализа синхронизации системы с помощью MSF можно разбить на несколько шагов. Рассмотрим сеть идентичных связанных нейронов, каждый из которых описывается определённой системой диф. уравнений (например, модель Ходжкина — Хаксли). Тогда сеть с N нейронами можно записать в таком виде:

${\displaystyle {\dot {\mathbf {x} }}_{i}=\mathbf {f} (\mathbf {x} _{i})+\sum _{j}G_{ij}\mathbf {h} (\mathbf {x} _{j}),\qquad (1)}$

где ${\displaystyle i,j\in 1,2,\dots ,N}$; ${\displaystyle \mathbf {x} _{i}\in \mathbb {R} ^{D}}$ это D-мерный вектор переменных i-го нейрона (вектор состояния); ${\displaystyle \mathbf {f} :\mathbb {R} ^{D}\to \mathbb {R} ^{D}}$ и ${\displaystyle \mathbf {h} :\mathbb {R} ^{D}\to \mathbb {R} ^{D}}$ векторные поля, которые описывают соответственно (нелинейную) внутреннюю динамику (обычное мат. определение нейрона) и функцию связи между элементами; ${\displaystyle G_{ij}}$ элемент матрицы ${\displaystyle \mathbf {G} \in \mathbb {R} ^{N\times N}}$ , который содержит силы синаптических связей между нейронами.

Итак, общая размерность фазового пространства всей сети будет равна ${\displaystyle N\cdot D}$. Часто нужно анализировать сеть с большим количеством нейронов (${\displaystyle N>1000}$), и тогда прямой анализ стабильности синхронизированного состояния сети (с помощью, например, экспонент Ляпунова) становится слишком громоздкой вычислительной задачей. С другой стороны, согласно формализму MSF, для определения стабильности синхронизированного состояния вышеприведённой системы, достаточно рассчитать экспоненты Ляпунова для гораздо меньшей, линеаризованной D-мерной системы, которая задаётся уравнением (предположим, что матрица связей имеет исключительно реальные собственные числа):

${\displaystyle {\dot {\mathbf {\zeta } }}=\left[D\mathbf {f} (\mathbf {x} _{s})+\alpha D\mathbf {h} (\mathbf {x} _{s})\right]\mathbf {\zeta } ,\qquad (2)}$

где ${\displaystyle D\mathbf {\cdot } (\mathbf {x} _{s})}$ якобиан векторного поля, рассчитанный вдоль синхронизированного решения ${\displaystyle \mathbf {x} _{s}}$ (под синхронизированным решением здесь имеется в виду вычислительное решение уравнения одного нейрона, соединённого с самим собой), а ${\displaystyle \alpha }$ — специальный скалярный параметр.

Опустив детали, имеем следующий алгоритм:

1. Рассчитать максимальную экспоненту Ляпунова (${\displaystyle \lambda _{max}}$) небольшой линейной системы (2) для нужного диапазона параметра ${\displaystyle \alpha }$;
2. Построить график зависимости ${\displaystyle \lambda _{max}}$ от параметра ${\displaystyle \alpha }$ и определить диапазоны параметра, где ${\displaystyle \lambda _{max}}$ имеет отрицательный знак;
3. Найти собственные числа матрицы связей ${\displaystyle \mathbf {G} }$
4. Если все собственные числа матрицы связей (за исключением того, которое равняется строковой сумме матрицы, см.^[40]) попадают в диапазон ${\displaystyle \alpha }$, что соответствует отрицательной ${\displaystyle \lambda _{max}}$ — сеть имеет стабильное состояние полной синхронизации

Пример графика MSF для HR нейронов, а также видео процесса полной синхронизации небольшой сети (согласно предсказанию MSF) приведены справа.

• Buzsáki, György. Rhythms of the Brain (англ.). — Oxford University Press, 2006. — ISBN 978-0-19-530106-9.
• Скляров Олег Павлович. Нарушения ритма речи, возникающие в результате синхронизации, индуцированной в связанных осцилляторах // Техническая акустика : журнал. — 2005. — № 5. — С. 420—423.
• Бушов Юрий Валентинович, Светлик Михаил Васильевич. Фазовые взаимодействия между высокои низкочастотными ритмами ЭЭГ и их роль в когнитивных процессах // Вестник Томского государственного университета. Биология : журнал. — 2012. — № 3. — С. 137—145.
• Цукерман В.Д. Математическая модель фазового кодирования событий в мозге // Математическая биология и биоинформатика. — 2006. — Т. 1, № 1. — С. 97—107.
• Бушов Ю.В., Светлик М.В., Крутенкова Е.П. Гамма- активность коры головного мозга: связь с интеллектом и точностью восприятия времени // Физиология человека. — Академический научно-издательский, производственно-полиграфический и книгораспространительский центр РАН «Издательство „Наука“» (Москва), 2010. — Т. 36, № 4. — С. 15—21. — ISSN 0131-1646.