thumb||alt=Ilustração da propagação do potencial de ação em um neurônio ao longo de seu axônio.|Ilustração da propagação do potencial de ação em um neurônio. Modelos de disparos neuronais são modelos matemáticos que descrevem os padrões com os quais potenciais de ação são iniciados e propagados nos neurônios. Tais modelos compreendem descrições matemáticas das propriedades das células no sistema nervoso que possuem potenciais elétricos. Grosso modo, trata-se de descrições, em linguagem matemática, da comunicação entre os neurônios. Esses modelos são utilizados nos estudos da dinâmica de redes neurais e fundamentam a neurociência computacional e a neuromatemática. Variados modelos são utilizados em pesquisas de neurociência, com diversas implicações no modo como as redes neurais são descritas e simuladas.
A quantidade de modelos de disparos neuronais é alta. Essa situação está em parte associada a características de diferentes experimentos neurocientíficos. De modo geral, pode-se distinguir os modelos em algumas famílias, conforme sua abordagem. Abaixo, um panorama é feito de alguns desses modelos.
Modelos elétricos de entrada-saída
Hodgkin-Huxley
thumb|right|260x260px|Potencial de membrana (em milivolts) do modelo Hodgkin-Huxley com uma corrente injetada variando de −5 a 12 nanoampères.|alt=O gráfico do potencial de membrana v(t) passando por três estágios: uma fase de equilíbrio, uma fase de um disparo único e uma fase de ciclo limite, conforme a variação da corrente injetada de -5nA a 12nA. O modelo de Hodgkin-Huxley, ou modelo baseado em condutância, é um modelo matemático que descreve como potenciais de ação nos neurônios são iniciados e propagados. É um conjunto de equações diferenciais ordinárias não-lineares, e, portanto, é um modelo de tempo contínuo, que se aproxima das características elétricas das células excitáveis, tais como neurônios e miócitos cardíacos. Alan Lloyd Hodgkin e Andrew Huxley descreveram o modelo em 1952 para explicar os mecanismos iônicos subjacentes à iniciação e propagação dos potenciais de ação no axônio gigante de uma lula. É um dos principais modelos e um dos grandes avanços da biofísica no século XX. Rendeu aos seus autores o prêmio Nobel de Fisiologia ou Medicina de 1963.
É composto por uma série de quatro equações diferenciais não lineares, com suas variáveis de estado descrevendo o potencial de membrana, a ativação de correntes de íon sódio e de íon potássio e a inativação da corrente de íon sódio. Denota-se a relação tensão-corrente, com várias correntes dependentes de voltagem, que carregam a capacitância da membrana com
.
Cada corrente é dada por
,
em que é o inverso da resistência elétrica, ou condutância, que pode ser expandida em termos da sua média constante e as frações dos canais de ativação e inativação e , respectivamente, que determinam quantos íons podem fluir através de canais de membrana disponíveis. Tal expansão é
e as frações seguem a cinética de primeira ordem
com uma dinâmica similar para , em que se pode usar tanto e ou e para definir a fração do canal.
O modelo pode ser expandido para incluir outras correntes iônicas. Geralmente, essas consideram a adição de correntes de íons Ca2+ e Na+ para o interior do neurônio e diversas variedades de correntes de K+ para o fora do neurônio, incluindo uma corrente de "vazamento".
Integra-e-dispara O modelo integra-e-dispara (em inglês, integrate-and-fire ou I&F) é um dos mais utilizados pela neurociência computacional. É um dos primeiros modelos de disparos neuronais criados, tendo sido proposto por Louis Lapicque em 1907. Não é um modelo explícito, isto é, não descreve o mecanismo que gera o potencial de ação de um neurônio. Pode ser descrito como
,
sendo a corrente de entrada e o potencial da membrana. A fórmula é a derivada no tempo da lei da capacitância, . Quando uma corrente é aplicada, a voltagem da membrana aumenta até atingir uma constante , quanto então o neurônio dispara, como uma função delta de Dirac, e a voltagem retorna ao valor de potencial de repouso. Assim, a frequência de disparo cresce linearmente sem restrições com o aumento da corrente.
Integra-e-dispara adaptativo Para lidar com a uni-dimensionalidade do IF, é possível adicionar uma segunda equação linear descrevendo uma dinâmica de ativação para uma corrente de potássio de limite alto para permitir uma adaptação de frequência de disparos no modelo. Assim, com
a cada disparo aumenta-se o valor de ativação de K, via função delta de Dirac . e denotam, respectivamente, as primeiras derivadas temporais das variáveis e .
Integra-e-dispara-ou-rajada Outra melhoria feita a partir do modelo inicial de integra-e-dispara é o modelo integra-e-dispara-ou-rajada ou I&FB:
,
considerando-se que caso , então .
e denotam, respectivamente, as primeiras derivadas temporais das variáveis e . , , , e são parâmetros para descrever a dinâmica da corrente canal de cálcio T, a intativação dessa corrente, e uma função degrau de Heaviside. A segunda variável permite modelar rajadas.
Ressoa-e-dispara Enquanto o modelo integra-e-dispara é uni-dimensional, o modelo ressoa-e-dispara é uma forma análoga em versão bi-dimensional, em que representa o potencial de membrana, e é descrito por
,
considerando-se que caso , então .
denota a primeira derivada temporal de . , , e são parâmetros, e é uma função arbitrária que descreve o reinicio do sistema após o disparo.
Integra-e-dispara quadrático O modelo integra-e-dispara quadrático (ou QIF) é uma alternativa aos modelos de integra-e-dispara com vazamento. Em Latham et al. o potencial de membrana é modelado por e sua dinâmica é descrita como
,
considerando-se que caso , então .
denota a primeira deriava temporal de . e são os valores de repouso e de limiar máximo do potencial de membrana. De acordo com Izhikevich, “este deve ser o modelo de escolha quando alguém simula redes de integradores em grande escala”.
Modelos estocásticos Modelos nesta categoria são generalizações dos casos integra-e-dispara com um certo nível de estocasticidade. A estocasticidade pode aparecer na modelagem de neurônios por dois motivos. O primeiro ocorre devido à abertura e ao fechamento de canais iônicos, que resultam em uma pequena variabilidade no valor do potencial de membrana e no instante de tempo dos disparos. O outro motivo decorre do fato que, para um neurônio inserido em uma rede cortical, é difícil controlar de qual outro neurônio um sinal de estímulo chega, dado que a maioria dos impulsos são oriundos de neurônios em alguma outra região do cérebro não observada.
Galves-Löcherbach
thumb|Visualização 3D do modelo de Galves-Löcherbach com 4.000 neurônios (4 camadas com uma população de neurônios inibitórios e uma população de neurônios excitatórios cada) em 180 intervalos de tempo. O modelo de Galves-Löcherbach consiste de rede de neurônios idealizados que interagem por eventos discretos esporádicos quase instantâneos (picos ou disparos). A cada instante, cada neurônio dispara independentemente, com uma probabilidade dependente do histórico de disparos dos demais neurônios.
Dessa forma, a probabilidade de um dado neurônio disparar em um certo tempo , de acordo com o modelo probabilístico, é dada por
,
em que é um peso sináptico que representa o aumento do potencial de ação do neurônio devido ao disparo do neurônio , é uma função que modela o vazamento de potencial e é o momento de disparo mais recente do neurônio antes do tempo em questão , descrito por
.
Modelos didáticos Os modelos dessa categoria são simplificações que descrevem qualitativamente o potencial de membrana como função de uma entrada. Enquanto são mais simples de implementar quando comparado com modelos mais realistas como o de Hodgkin-Huxley, apresentam uma robustez maior do que os modelos simples de integra-e-dispara. Geralmente são baseados em sistemas de equações diferenciais que ajustam medidas reais, suas principais aplicações são didáticas e em simulações com uma grande quantidade de neurônios.
FitzHugh-Nagumo
alt=Trata-se do gráfico da evolução temporal da variável x. Estão presentes 4 ciclos completos e parte de um outro, em que a variável altera rapidamente por valores entre 2 e -2.|miniaturadaimagem|259x259px|Gráfico da evolução temporal de no modelo de FitzHugh-Nagumo. O modelo de FitzHugh-Nagumo, originalmente chamado de oscilador Bonhoeffer–van der Pol, é um exemplo de sistema dinâmico excitatório-oscilatório com duas variáveis de estado. , é a variável relacionada ao potencial de membrana, enquanto é responsável pela acomodação e resistência do sistema. Sua dinâmica é descrita por
,
em que e representam, respectivamente, as primeiras derivadas temporais das variáveis e .
Hindmarsh-Rose
alt=Gráfico da evolução temporal da variável x. Constam 2 sequências de rajadas de disparos rápidos e de elevada amplitude, apresentam um considerável tempo de acomodação entre as mesmas.|miniaturadaimagem|260x260px|Gráfico da evolução temporal de no modelo de Hindmarsh-Rose. O modelo de Hindmarsh-Rose é um sistema dinâmico baseado em três variáveis de estado capaz de modelar disparos únicos e em rajadas. refere-se ao potencial de membrana, é a taxa de transporte dos íons sódio e potássio feito por canais rápidos, enquanto o transporte de outros íons, como o cálcio, que é feito através de canais lentos, é representado por . O sistema que descreve esse modelo é
,
em que , e representam, respectivamente, as primeiras derivadas temporais das variáveis , e .
Morris-Lecar
alt=Gráfico da evolução temporal da variável V(t) com correntes excitatórias crescentes. Conforme a corrente injetada aumenta, surgem os disparos.|miniaturadaimagem|260x260px|Gráfico da evolução temporal de no modelo de Morris-Lecar com diferentes correntes de excitação. O modelo de Morris-Lecar é um modelo bidimensional desenvolvido para reproduzir oscilações na fibra muscular de uma cirripedia. As variáveis são o potencial de membrana, , e a variável de recuperação, , que indica a probabilidade de o canal de íon potássio estar conduzindo corrente. O modelo pode considerado uma simplificação do modelo de Hodgkin-Huxley e é baseado no sistema
,
em que
e e denotam as primeiras derivadas temporais das variáveis e .
Rulkov
miniaturadaimagem|261x261px|Gráfico da evolução temporal de no modelo de Rulkov em regime de disparo único (a) e de rajadas (b).|alt=Dois gráficos demonstrando os regimes de disparo único e de rajada descritos pelo modelo. O mapa de Rulkov é um modelo de disparo neuronal baseado em um sistema dinâmico discreto. Foi proposto por Nikolai F. Rulkov em 2001. representa o potencial de membrana do neurônio e , por mais que não tenha significado biológico explícito, pode ser feita alguma analogia com uma variável de controle simulando um canal iônico. Seu comportamento é descrito pelo sistema
,
em que varia conforme as condições
.
Izhikevich thumb|Visualização 3D do modelo de Izhikevich com 4.000 neurônios (4 camadas com uma população de neurônios inibitórios e uma população de neurônios excitatórios cada) em 180 intervalos de tempo. O matemático Eugene Izhikevich propôs um modelo a partir do qual podem ser reproduzidos diversos padrões de disparo observados no córtex motor de ratos. O sistema que descreve o modelo é baseado na redução, por metodologias de bifurcações de outros modelos mais biofisicamente apropriados, como o de Hodgkin-Huxley, no sistema bidimensional de equações diferenciais ordinárias
com a inclusão de um mecanismo de reinício após um disparo conforme se mV, então ,
em que e denotam as primeiras derivadas temporais das variáveis e , que são respectivamente, o potencial de membrana e a variável de recuperação da membrana.
Wilson