Модели нейронных вычислений - это попытки абстрактно и математически объяснить основные принципы, лежащие в основе обработки информации в биологических нервных системах или их функциональных компонентах. Эта статья призвана предоставить обзор наиболее точных моделей нейробиологических вычислений, а также инструментов, обычно используемых для их построения и анализа.
Вступление
Из-за сложности поведения нервной системы соответствующие границы экспериментальной ошибки не определены, но относительные достоинства различных моделей конкретной подсистемы можно сравнивать в зависимости от того, насколько точно они воспроизводят поведение в реальном мире или реагируют на определенные входные сигналы. . В тесно связанной области вычислительной нейроэтологии практика заключается в том, чтобы включить среду в модель таким образом, чтобы цикл был замкнутым . В тех случаях, когда конкурирующие модели недоступны или когда были измерены или количественно оценены только общие ответы, четко сформулированная модель может помочь ученому в разработке экспериментов для исследования биохимических механизмов или сетевых подключений.
Во всех случаях, кроме простейших, математические уравнения, составляющие основу модели, не могут быть решены точно. Тем не менее компьютерные технологии, иногда в виде специализированных программных или аппаратных архитектур, позволяют ученым выполнять итерационные вычисления и искать вероятные решения. Компьютерный чип или робот, который может взаимодействовать с природной средой способами, подобными исходному организму, является одним из воплощений полезной модели. Однако высшим мерилом успеха является способность делать проверяемые прогнозы.
Общие критерии оценки моделей
Скорость обработки информации
Скорость обработки информации в биологических нейронных системах ограничена скоростью, с которой потенциал действия может распространяться по нервному волокну. Эта скорость проводимости колеблется от 1 м / с до более 100 м / с и обычно увеличивается с увеличением диаметра нейронального отростка. Медленные по шкале времени биологически значимых событий, продиктованные скоростью звука или силой тяжести, нервная система в подавляющем большинстве предпочитает параллельные вычисления последовательным в критичных по времени приложениях.
Надежность
Модель считается устойчивой, если она продолжает давать те же результаты вычислений при изменениях входных или рабочих параметров, вызванных шумом. Например, направление движения, вычисленное надежным детектором движения , не изменится при небольших изменениях яркости , контраста или дрожания скорости. Для простых математических моделей нейрона, например, зависимость паттернов спайков от задержки сигнала намного слабее, чем зависимость от изменений «весов» межнейронных связей. [1]
Получить управление
Это относится к принципу, согласно которому реакция нервной системы должна оставаться в определенных пределах, даже если входные данные из окружающей среды резко меняются. Например, при переходе между солнечным днем и безлунной ночью сетчатка изменяет соотношение между уровнем света и выходом нейронов более чем в 1 раз.так что сигналы, посылаемые на более поздние стадии зрительной системы, всегда остаются в гораздо более узком диапазоне амплитуд. [2] [3] [4]
Линейность против нелинейности
Линейная система является одним ответ которого в заданной единице измерения, к набору входов , рассматриваемых одновременно, является суммой его ответов в связи с входами , рассматриваемых по отдельности.
Линейные системы легче анализировать математически, и они являются убедительным допущением во многих моделях, включая нейрон Мак-Каллока и Питтса, модели кодирования населения и простые нейроны, часто используемые в искусственных нейронных сетях . Линейность может возникать в основных элементах нейронной цепи, таких как реакция постсинаптического нейрона, или как возникающее свойство комбинации нелинейных подсхем. [5] Хотя линейность часто рассматривается как неправильная, в недавних работах предполагалось, что в некоторых случаях она может быть биофизически правдоподобной. [6] [7]
Примеры
Вычислительная нейронная модель может быть ограничена уровнем биохимической передачи сигналов в отдельных нейронах или может описывать весь организм в его окружающей среде. Примеры здесь сгруппированы в соответствии с их объемом.
Модели передачи информации в нейронах
Наиболее широко используемые модели передачи информации в биологических нейронах основаны на аналогиях с электрическими цепями. Решаемые уравнения представляют собой зависящие от времени дифференциальные уравнения с электродинамическими переменными, такими как ток, проводимость или сопротивление, емкость и напряжение.
Модель Ходжкина – Хаксли и ее производные.
Модель Ходжкина-Хаксли, широко известная как одно из величайших достижений биофизики 20-го века, описывает, как потенциалы действия в нейронах инициируются и распространяются в аксонах через потенциалзависимые ионные каналы . Это набор нелинейных обыкновенных дифференциальных уравнений, которые были введены Аланом Ллойдом Ходжкином и Эндрю Хаксли в 1952 году для объяснения результатов экспериментов по фиксации напряжения на гигантском аксоне кальмара . Аналитических решений не существует, но алгоритм Левенберга – Марквардта , модифицированный алгоритм Гаусса – Ньютона , часто используется для согласования этих уравнений с данными фиксации напряжения.
Модель ФитцХью-Нагумо является упрощение модели Ходжкина-Хаксли. Модель Хиндмарша – Роуза - это расширение, которое описывает всплески нейрональных спайков. Модель Морриса-Лекара является модификацией, которая не генерирует спайков, но описывает распространение медленных волн, которое участвует в тормозных синаптических механизмах центральных генераторов паттернов .
Солитоны
Модель солитон является альтернативой модели Ходжкина-Хаксли , что требования объяснить , как потенциалы действия инициируются и проводятся в форме определенных видов одиночного звука (или плотности ) импульсов , которые могут быть смоделированы как солитоны вдоль аксонов , основанный на термодинамической теория распространения нервного импульса.
Передаточные функции и линейные фильтры
Этот подход, основанный на теории управления и обработке сигналов , рассматривает нейроны и синапсы как инвариантные во времени объекты, которые производят выходные данные, которые представляют собой линейные комбинации входных сигналов, часто изображаемых как синусоидальные волны с четко определенными временными или пространственными частотами.
Все поведение нейрона или синапса закодировано в передаточной функции , несмотря на отсутствие знаний о точном механизме, лежащем в основе. Это заставляет высокоразвитую математику решать проблему передачи информации.
Сопутствующая таксономия линейных фильтров оказывается полезной для характеристики нейронных схем. Оба с низкими и высокочастотными фильтрами постулируются существовать в некоторой форме в сенсорных системах, так как они действуют , чтобы предотвратить потерю информации в условиях высокой и низкой контрастность, соответственно.
Действительно, измерения передаточных функций нейронов сетчатки подковообразного краба в соответствии с анализом линейных систем показывают, что они устраняют краткосрочные колебания входных сигналов, оставляя только долгосрочные тенденции, подобно фильтрам нижних частот. Эти животные не могут видеть малоконтрастные объекты без помощи оптических искажений, вызванных подводными течениями. [8] [9]
Модели вычислений в сенсорных системах
Боковое торможение в сетчатке: уравнения Хартлина – Рэтлиффа
В сетчатке возбужденный нервный рецептор может подавлять активность окружающих нейронов в области, называемой тормозным полем. Этот эффект, известный как латеральное торможение , увеличивает контраст и резкость зрительного отклика, но приводит к эпифеномену полос Маха . Это часто иллюстрируется оптической иллюзией светлых или темных полос рядом с резкой границей между двумя областями на изображении разной яркости.
Модель Hartline-Ratliff описывает взаимодействия внутри группы p- фоторецепторных клеток . [10] Предполагая, что эти взаимодействия являются линейными , они предложили следующую зависимость для скорости реакции в установившемся режиме. данного p -го фоторецептора с точки зрения скорости реакции в установившемся режимеиз j окружающих рецепторов:
.
Здесь,
возбуждение целевого p-го рецептора от сенсорной трансдукции
- связанный порог срабатывания ячейки, а
- коэффициент ингибирующего взаимодействия между p -м и j- м рецепторами. Ингибирующее взаимодействие уменьшается по мере удаления от целевого p-го рецептора.
Взаимная корреляция в локализации звука: модель Джеффресса
Согласно Джеффрессу , [11] для того , чтобы вычислить местоположение источника звука в пространстве от интерауральных разницы во время , слуховая система опирается на линии задержки : индуцированный сигнал от ипсилатерального слухового рецептора к конкретному нейрону задерживаются на то же самое время так как исходный звук распространяется в пространстве от одного уха до другого. Каждая постсинаптическая клетка задерживается по-разному и, таким образом, специфична для определенной разницы во времени между ушами. Эта теория эквивалентна математической процедуре взаимной корреляции .
Следуя Фишеру и Андерсону [12], ответ постсинаптического нейрона на сигналы из левого и правого уха определяется выражением
где
а также
представляет функцию задержки. Это не совсем правильно, и для того, чтобы расставить символы по порядку, нужен ясный взгляд.
В «сипухе» были обнаружены конструкции, соответствующие механизмам типа Джеффресса. [13]
Взаимная корреляция для обнаружения движения: модель Хассенштейна – Райхардта
Детектор движения должен удовлетворять трем основным требованиям: пара входов, асимметрия и нелинейность. [14] Операция взаимной корреляции, реализованная асимметрично на ответах пары фоторецепторов, удовлетворяет этим минимальным критериям и, кроме того, предсказывает особенности, которые наблюдались в ответе нейронов дольчатой пластинки у двукрылых насекомых. [15]
Основное уравнение для ответа:
Модель HR предсказывает пик отклика на определенной входной временной частоте. Концептуально подобная модель Барлоу – Левика несовершенна в том смысле, что стимула, предъявляемого только к одному рецептору из пары, достаточно для генерации ответа. Это отличается от модели HR, которая требует, чтобы два коррелированных сигнала доставлялись в упорядоченном по времени образом. Однако модель HR не показывает насыщения отклика при высоких контрастах, что наблюдается в эксперименте. Расширения модели Барлоу-Левика могут компенсировать это несоответствие. [16]
Модель Уотсона – Аумады для оценки движения человека
Это использует взаимную корреляцию как в пространственном, так и во временном направлениях и связано с концепцией оптического потока . [17]
Антихеббовская адаптация: пластичность, зависящая от времени всплеска
- Цунопулос, Т; Kim, Y; Oertel, D; Трассел Л.О. (2004). «Клеточно-специфическая пластичность, зависящая от времени спайков в дорсальном ядре улитки». Nat Neurosci . 7 (7): 719–725. DOI : 10.1038 / nn1272 . PMID 15208632 . S2CID 17774457 .
- Робертс, Патрик Д.; Портфорс, Кристина В. (2008). «Принципы проектирования сенсорной обработки в мозжечковых структурах». Биологическая кибернетика . 98 (6): 491–507. DOI : 10.1007 / s00422-008-0217-1 . PMID 18491162 . S2CID 14393814 .
Модели сенсорно-моторной связи
Нейрофизиологические метрономы: нейронные цепи для генерации паттернов
Взаимно тормозящие процессы являются объединяющим мотивом всех центральных генераторов паттернов . Это было продемонстрировано на стоматогастральной нервной системе раков и омаров. [18] На основе STG были построены осциллирующие сети с двумя и тремя ячейками, которые поддаются математическому анализу и которые простым образом зависят от силы синапсов и общей активности, предположительно от ручек этих вещей. [19] Используемая математика - это теория динамических систем .
Обратная связь и управление: модели управления полетом на лету
Считается, что управление полетом у мухи обеспечивается за счет сигналов от зрительной системы, а также от жужжальцев , пары органов, похожих на ручки, которые измеряют угловую скорость. Интегрированные компьютерные модели дрозофилы , лишенные нейронных схем, но основанные на общих рекомендациях, данных теорией управления и данными о привязанных полетах мух, были построены для исследования деталей управления полетом. [20] [21]
Сенсорно-двигательный контроль мозжечка
Теория тензорных сетей - это теория функции мозжечка, которая обеспечивает математическую модель преобразования сенсорных пространственно-временных координат в моторные координаты и наоборот нейронными сетями мозжечка . Теория была разработана Андрасом Пеллионисом и Родольфо Ллинасом в 1980-х годах как геометризация функций мозга (особенно центральной нервной системы ) с использованием тензоров . [22] [23]
Подходы и инструменты программного моделирования
Нейронные сети
В этом подходе сила и тип, возбуждающий или тормозящий, синаптических связей представлены величиной и знаком весов, то есть числовыми коэффициентами перед входами к конкретному нейрону. Ответ-й нейрон задается суммой нелинейных, обычно « сигмоидальных » функций входов как:
.
Затем этот ответ вводится в другие нейроны и так далее. Цель состоит в том, чтобы оптимизировать веса нейронов для вывода желаемого ответа на выходном слое, соответствующего набору заданных входов на входном уровне. Эта оптимизация весов нейронов часто выполняется с использованием алгоритма обратного распространения ошибки и метода оптимизации, такого как градиентный спуск или метод оптимизации Ньютона . Обратное распространение сравнивает выходные данные сети с ожидаемыми выходными данными обучающих данных, а затем обновляет веса каждого нейрона, чтобы минимизировать вклад этого отдельного нейрона в общую ошибку сети.
Генетические алгоритмы
Генетические алгоритмы используются для развития нейронных (а иногда и телесных) свойств в модельной системе мозг-тело-окружающая среда, чтобы продемонстрировать желаемые поведенческие характеристики. Затем разработанные агенты могут быть подвергнуты детальному анализу, чтобы раскрыть их принципы действия. Эволюционные подходы особенно полезны для исследования пространства возможных решений данной поведенческой задачи, потому что эти подходы минимизируют априорные предположения о том, как должно быть реализовано данное поведение. Они также могут быть полезны для изучения различных способов завершения вычислительной модели нейроэтологии, когда для интересующей биологической системы доступны только частичные нейронные схемы. [24]
НЕЙРОН
Программное обеспечение NEURON , разработанное в Университете Дьюка, представляет собой симуляционную среду для моделирования отдельных нейронов и сетей нейронов. [25] Среда NEURON - это автономная среда, позволяющая взаимодействовать через графический интерфейс или сценарии с помощью hoc или python . Механизм моделирования NEURON основан на модели типа Ходжкина – Хаксли с использованием формулировки Борга – Грэма. Несколько примеров моделей, написанных на NEURON, доступны в онлайн-базе данных ModelDB. [26]
Воплощение в электронном оборудовании
Кремниевые нейроны на основе проводимости
Нервные системы отличаются от большинства вычислительных устройств на основе кремния тем, что они напоминают аналоговые компьютеры (не процессоры цифровых данных ) и массивно- параллельные процессоры, а не последовательные процессоры. Для точного моделирования нервной системы в режиме реального времени требуется альтернативное оборудование.
Наиболее реалистичные на сегодняшний день схемы используют аналоговые свойства существующей цифровой электроники (работающей в нестандартных условиях) для реализации in silico моделей типа Ходжкина – Хаксли . [27] [28]
Ретиноморфные чипы
[29]
Смотрите также
- Когнитивная архитектура
- Познавательная карта
- Вычислительная нейробиология
- Восприятие движения
- Нейронное кодирование
- Нейронный коррелят
- Нейронное декодирование
- Нейроэтология
- Нейроинформатика
- Количественные модели потенциала действия
- Пиковая нейронная сеть
- Системная неврология
Рекомендации
- ^ Цейнар, Павел; Вышата, Олдржих; Валиш, Мартин; Прохазка, Алеш (2019). «Сложное поведение простой модели нейронного осциллятора в коре головного мозга человека». IEEE Transactions по нейронным системам и реабилитационной инженерии . 27 (3): 337–347. DOI : 10.1109 / TNSRE.2018.2883618 . PMID 30507514 . S2CID 54527064 .
- ^ Николас Дж. Прибе и Дэвид Ферстер (2002). «Новый механизм управления усилением нейронов». 35 (4). Нейрон. С. 602–604.
- Перейти ↑ Klein, SA, Carney, T., Barghout-Stein, L., & Tyler, CW (1997, июнь). Семь моделей маскировки. В Electronic Imaging'97 (стр. 13–24). Международное общество оптики и фотоники.
- Перейти ↑ Barghout-Stein, Lauren. О различиях между маскировкой периферического и фовеального рисунка. Дисс. Калифорнийский университет, Беркли, 1999.
- ^ Мольнар, Алеша; Сюэ, Хайн-Анн; Роска, Ботонд; Верблин, Франк С. (2009). «Торможение кроссовера в сетчатке: схема, которая компенсирует нелинейную выпрямляющую синаптическую передачу» . Журнал вычислительной нейробиологии . 27 (3): 569–590. DOI : 10.1007 / s10827-009-0170-6 . ISSN 0929-5313 . PMC 2766457 . PMID 19636690 .
- ^ Сингх, Чандан; Леви, Уильям Б. (13 июля 2017 г.). «Пирамидный нейрон консенсусного слоя V может выдерживать межимпульсное интервальное кодирование» . PLOS One . 12 (7): e0180839. arXiv : 1609.08213 . Bibcode : 2017PLoSO..1280839S . DOI : 10.1371 / journal.pone.0180839 . ISSN 1932-6203 . PMC 5509228 . PMID 28704450 .
- ^ Кэш, Сидней; Юсте, Рафаэль (1 января 1998 г.). «Суммирование входных данных культивированными пирамидальными нейронами является линейным и независимым от положения» . Журнал неврологии . 18 (1): 10–15. DOI : 10.1523 / JNEUROSCI.18-01-00010.1998 . ISSN 0270-6474 . PMC 6793421 . PMID 9412481 .
- ^ Роберт Б. Барлоу младший; Рамкришна Пракаш; Эдуардо Солессио (1993). «Нейронная сеть сетчатки Limulus: от компьютера к поведению» . Американский зоолог . 33 : 66–78. DOI : 10.1093 / ICB / 33.1.66 .
- ^ Роберт Б. Барлоу, Джеймс М. Хитт и Фредерик А. Додж (2001). «Видение Limulus в морской среде». Биологический бюллетень . 200 (2): 169–176. CiteSeerX 10.1.1.116.5190 . DOI : 10.2307 / 1543311 . JSTOR 1543311 . PMID 11341579 .
- ^ К.П. Хаделер и Д. Кун (1987). «Стационарные состояния модели Хартлина – Рэтлиффа». Биологическая кибернетика . 56 (5–6): 411–417. DOI : 10.1007 / BF00319520 . S2CID 8710876 .
- ^ Джеффресс, Лос-Анджелес (1948). «Теория места локализации звука». Журнал сравнительной и физиологической психологии . 41 (1): 35–39. DOI : 10.1037 / h0061495 . PMID 18904764 .
- ^ Фишер, Брайан Дж .; Андерсон, Чарльз Х. (2004). «Расчетная модель локализации звука у сипухи». Нейрокомпьютеры . 58–60: 1007–1012. DOI : 10.1016 / j.neucom.2004.01.159 .
- ^ Кэтрин Э. Карр, 1993. Модели линий задержки локализации звука в сипухе "Американский зоолог" Vol. 33, № 1 79–85
- ^ Borst A, Egelhaaf M., 1989. Принципы визуального обнаружения движения. «Тенденции в нейронауках» 12 (8): 297–306
- ^ Joesch, M .; и другие. (2008). «Свойства реакции чувствительных к движению зрительных интернейронов в лобулярной пластинке Drosophila melanogaster». Curr. Биол . 18 (5): 368–374. DOI : 10.1016 / j.cub.2008.02.022 . PMID 18328703 . S2CID 18873331 .
- ^ Гонсало Г. де Polavieja, 2006. Нейронные алгоритмыкоторые обнаруживают временной порядок событий «Neural Исчисление» 18 (2006) 2102-2121
- ^ Эндрю Б. Уотсон и Альберт Дж. Ахумада младший, 1985. Модель человеческого зрительного восприятия движения "J. Opt. Soc. Am. A" 2 (2) 322–341
- ^ Майкл П. Нусбаум и Марк П. Бенхаккер, Подход малых систем к генерации двигательных паттернов, Nature 417, 343–350 (16 мая 2002 г.)
- ↑ Cristina Soto-Treviño, Kurt A. Thoroughman и Eve Marder, LF Abbott, 2006. Зависимая от активности модификация тормозных синапсов в моделях ритмических нейронных сетей Nature Vol 4 No 3 2102–2121
- ^ «Модель Grand Unified Fly (GUF)» .
- ^ http://www.mendeley.com/download/public/2464051/3652638122/d3bd7957efd2c8a011afb0687dfb6943731cb6d0/dl.pdf [ постоянная мертвая ссылка ]
- ^ Пеллионис А., Ллинас Р. (1980). «Тензорный подход к геометрии функции мозга: координация мозжечка с помощью метрического тензора» (PDF) . Неврология . 5 (7): 1125–1136. DOI : 10.1016 / 0306-4522 (80) 90191-8 . PMID 6967569 . S2CID 17303132 .CS1 maint: несколько имен: список авторов ( ссылка )
- ^ Пеллионис А., Ллинас Р. (1985). «Тензорная сетевая теория метаорганизации функциональной геометрии в центральной нервной системе». Неврология . 16 (2): 245–273. DOI : 10.1016 / 0306-4522 (85) 90001-6 . PMID 4080158 . S2CID 10747593 .CS1 maint: несколько имен: список авторов ( ссылка )
- ^ Пиво, Рэндалл; Чиль, Гилель (4 марта 2008 г.). «Вычислительная нейроэтология» . Scholarpedia . 3 (3): 5307. Bibcode : 2008SchpJ ... 3.5307B . DOI : 10,4249 / scholarpedia.5307 .
- ^ «NEURON - для эмпирического моделирования нейронов и нейронных сетей» .
- ^ Макдугал Р.А., Морс TM, Карневале Т., Маренко Л., Ван Р., Мильоре М., Миллер П.Л., Шеперд GM, Хайнс М.Л. Двадцать лет ModelDB и далее: создание необходимых инструментов моделирования для будущего нейробиологии. J Comput Neurosci. 2017; 42 (1): 1–10.
- ^ L. Alvadoa, J. Tomasa, S. Saghia, S. Renauda, T. Balb, A. Destexheb, G. Le Masson, 2004. Аппаратные вычисления моделей нейронов на основе проводимости. Нейрокомпьютеры 58–60 (2004) 109–115
- ^ Индивери, Джакомо; Дуглас, Родни; Смит, Лесли (29 марта 2008 г.). «Кремниевые нейроны» . Scholarpedia . 3 (3): 1887. Bibcode : 2008SchpJ ... 3.1887I . DOI : 10,4249 / scholarpedia.1887 .
- ^ Квабена Боахен, "Ретиноморфный чип с параллельными путями: кодирование УВЕЛИЧЕНИЕ, ВКЛЮЧЕНИЕ, УМЕНЬШЕНИЕ и ВЫКЛЮЧЕНИЕ визуальных сигналов", Аналоговые интегральные схемы и обработка сигналов, 30, 121–135, 2002
Внешние ссылки
- Нейродинамика в NSI - Интернет-страница Патрика Д. Робертса в Институте неврологических наук
- Лаборатория Дикинсона - веб-страница группы Дикинсона в Калифорнийском технологическом институте, изучающей управление полетом у дрозофилы.