|
|
| Line 1: |
Line 1: |
| A '''biological neuron model''' (also known as '''spiking neuron model''') is a mathematical description of the properties of nerve cells, or [[neurons]], that is designed to accurately describe and predict biological processes. This is in contrast to the [[artificial neuron]], which aims for computational effectiveness, although these goals sometimes overlap.
| | Whatever your type could be, Iunderwearstore.com has everything you need! We offer in a large selection of mens underwear to steam up the bedroom. Shop at our large catalog, we assured you that clients attention is our utmost priority. Mens underwear area is growing day by day due to its growing demand. Now men are too getting style savvy and they want the variety due to their innerwear selection. This is the reason why we have got a huge variety for sexy underwear. If you desire to identify further on [https://www.youtube.com/watch?v=__j-NLVT9gw gay porn star], there are many libraries you should consider investigating. Their your ultimate on the web destination for everything attractive. <br><br>Please spend some time and pay particular focus on our site Iunderwearstore.com to truly have a good examination about our products. These Mens underwear are specially made to provide mobility and extra support. It allows you with regular breathing than some other mens underwear obtainable in industry. We"ve got the largest collection for guys underwear and if you are thinking that from where to start your collection then all the trendy and trendy Mens underwear are now online. There are numerous reasons to pride ourselves but we really care for our customers and just concern about their advantages. We"re often determined to supply our client the very best services possible. Browse here at [https://www.youtube.com/watch?v=__j-NLVT9gw gay sex] to discover the meaning behind this viewpoint. <br><br>At 1underwearstore.com, we also deal in Boy Underwear <br><br>, Fighter Briefs, Hot Underwear, Male Underwear, and Calvin Klein Underwear. We all know with such a massive collection you might get overcome. Click here [https://www.youtube.com/watch?v=ozVCMbdo0T0 date online] to learn the meaning behind this view. As mentioned previously we know exactly what a customer want and are extremely cautious about customer satisfaction. We set up our site to remove the confusion and make shopping for guys underwear as easy as possible. As a man if you know exactly what males underwear you wants then you are at the right location. Choose according to your need and budget. <br><br>More information please visit: http://www.1underwearstore.com/.<br><br>If you liked this post and you would like to get additional information about [http://fallaciousaccor38.xtgem.com government health insurance] kindly stop by our web page. |
| | |
| ==Artificial neuron abstraction==
| |
| The most basic model of a neuron consists of an input with some [[synaptic weight]] vector and an [[activation function]] or [[transfer function]] inside the neuron determining output. This is the basic structure used in artificial neurons, which in a [[neural network]] often looks like
| |
| | |
| <math>y_i = \phi\left( \sum_j w_{ij} x_j \right)</math>
| |
| | |
| where {{math|''y''<sub>''i''</sub>}} is the output of the {{math|''i''}} th neuron, {{math|''x''<sub>''j''</sub>}} is the {{math|''j''}}th input neuron signal, {{math|''w''<sub>''ij''</sub>}} is the synaptic weight (or strength of connection) between the neurons {{math|''i''}} and {{math|''j''}}, and {{math|''φ''}} is the [[activation function]]. Some of the earliest biological models took this form until kinetic models such as the [[Hodgkin-Huxley model]] became dominant.
| |
| | |
| ==Biological abstraction==
| |
| In the case of modelling a biological neuron, physical analogues are used in place of abstractions such as "weight" and "transfer function". The input to a neuron is often described by an ion [[Electric current|current]] through the [[cell membrane]] that occurs when [[neurotransmitters]] cause an activation of [[ion channels]] in the cell. We describe this by a physical time-dependent current {{math|''I''(''t'')}}. The cell itself is bound by an insulating cell membrane with a concentration of charged ions on either side that determines a [[capacitance]] {{math|''C''<sub>m</sub>}}. Finally, a neuron responds to such a signal with a change in [[voltage]], or an electrical potential energy difference between the cell and its surroundings, which is observed to sometimes result in a [[voltage spike]] called an [[action potential]]. This quantity, then, is the quantity of interest and is given by {{math|''V''<sub>m</sub>}}.
| |
| | |
| ===Integrate-and-fire===
| |
| One of the earliest models of a neuron was first investigated in 1907 by [[Louis Lapicque]].<ref>{{cite journal |last=Abbott |first=L.F. |authorlink= |coauthors= |year=1999 |month= |title=Lapique's introduction of the integrate-and-fire model neuron (1907) |journal=Brain Research Bulletin |volume=50 |issue=5/6 |pages=303–304 |id= |url=http://neurotheory.columbia.edu/~larry/AbbottBrResBul99.pdf |accessdate= 2007-11-24 |quote= |doi=10.1016/S0361-9230(99)00161-6 |pmid=10643408}}</ref> A neuron is represented in time by
| |
| | |
| :<math>I(t)=C_\mathrm{m} \frac{d V_\mathrm{m}}{d t}</math>
| |
| | |
| which is just the time [[derivative]] of the law of [[Capacitor#Capacitance|capacitance]], {{math|''Q'' {{=}} ''CV''}}. When an input current is applied, the membrane voltage increases with time until it reaches a constant threshold {{math|''V''<sub>th</sub>}}, at which point a [[Dirac delta function|delta function]] spike occurs and the voltage is reset to its resting potential, after which the model continues to run. The ''firing frequency'' of the model thus increases linearly without bound as input current increases.
| |
| | |
| The model can be made more accurate by introducing a [[refractory period (physiology)|refractory period]] {{math|''t''<sub>ref</sub>}} that limits the firing frequency of a neuron by preventing it from firing during that period. Through some calculus involving a [[Fourier transform]], the firing frequency as a function of a constant input current thus looks like
| |
| | |
| :<math>\,\! f(I)= \frac{I} {C_\mathrm{m} V_\mathrm{th} + t_\mathrm{ref} I}</math>.
| |
| | |
| A remaining shortcoming of this model is that it implements no time-dependent memory. If the model receives a below-threshold signal at some time, it will retain that voltage boost forever until it fires again. This characteristic is clearly not in line with observed neuronal behavior.
| |
| | |
| ===Leaky integrate-and-fire===
| |
| In the leaky integrate-and-fire model, the memory problem is solved by adding a "leak" term to the membrane potential, reflecting the diffusion of ions that occurs through the membrane when some equilibrium is not reached in the cell. The model looks like
| |
| | |
| :<math>I(t)-\frac{V_\mathrm{m} (t)}{R_\mathrm{m}} = C_\mathrm{m} \frac{d V_\mathrm{m} (t)}{d t}</math>
| |
| | |
| where {{math|''R''<sub>m</sub>}} is the membrane resistance, as we find it is not a perfect insulator as assumed previously. This forces the input current to exceed some threshold {{math|''I''<sub>th</sub> {{=}} ''V''<sub>th</sub> / ''R''<sub>m</sub>}} in order to cause the cell to fire, else it will simply leak out any change in potential. The firing frequency thus looks like
| |
| | |
| :<math>f(I) =
| |
| \begin{cases}
| |
| 0, & I \le I_\mathrm{th} \\
| |
| {[} t_\mathrm{ref}-R_\mathrm{m} C_\mathrm{m} \log(1-\tfrac{V_\mathrm{th}}{I R_\mathrm{m}}) {]}^{-1}, & I > I_\mathrm{th}
| |
| \end{cases} </math>
| |
| | |
| which converges for large input currents to the previous leak-free model with refractory period.<ref name="Koch98">{{cite book|last1=Koch|first1=Christof|title=Methods in neuronal modeling : from ions to networks|year=1999|pages=687|publisher=MIT Press|location=Cambridge, Mass.|isbn=0-262-11231-0|url=http://www.klab.caltech.edu/MNM/|edition=2nd|last2=Segev|first2=Idan}}</ref>
| |
| | |
| ===Exponential integrate-and-fire===
| |
| {{main|Exponential integrate-and-fire}}
| |
| In the [[Exponential Integrate-and-Fire]], spike generation is exponential, following the equation:
| |
| | |
| :<math> \frac{dX}{dt} = \Delta_T \exp \left( \frac{X - X_T} {\Delta_T} \right) </math>.
| |
| | |
| where <math>X</math> is the membrane potential, <math>X_T</math> is the membrane potential threshold, and <math>\Delta_T</math> is the sharpness of action potential initiation, usually around 1 mV for cortical pyramidal neurons. Once the membrane potential crosses <math>X_T</math>, it diverges to infinity in finite time.<ref name="Ostojic09">{{cite pmid| 19692598 }}</ref>
| |
| | |
| ===Hodgkin-Huxley===
| |
| {{main|Hodgkin-Huxley model}}
| |
| The most successful and widely used models of neurons have been based on the Markov kinetic model developed from Hodgkin and Huxley's 1952 work based on data from the [[squid giant axon]]. We note as before our voltage-current relationship, this time generalized to include multiple voltage-dependent currents:
| |
| | |
| :<math>C_\mathrm{m} \frac{d V(t)}{d t} = -\sum_i I_i (t, V)</math>.
| |
| | |
| Each current is given by [[Ohm's Law]] as
| |
| | |
| :<math>I(t,V) = g(t,V)\cdot(V-V_\mathrm{eq})</math>
| |
| | |
| where {{math|''g''(''t'',''V'')}} is the [[Electrical conductance|conductance]], or inverse resistance, which can be expanded in terms of its constant average {{math|''ḡ''}} <!-- \bar{g} --> and the activation and inactivation fractions {{math|''m''}} and {{math|''h''}}, respectively, that determine how many ions can flow through available membrane channels. This expansion is given by
| |
| | |
| :<math>g(t,V)=\bar{g}\cdot m(t,V)^p \cdot h(t,V)^q</math>
| |
| | |
| and our fractions follow the first-order kinetics
| |
| | |
| :<math>\frac{d m(t,V)}{d t} = \frac{m_\infty(V)-m(t,V)}{\tau_\mathrm{m} (V)} = \alpha_\mathrm{m} (V)\cdot(1-m) - \beta_\mathrm{m} (V)\cdot m</math>
| |
| | |
| with similar dynamics for {{math|''h''}}, where we can use either {{math|''τ''}} and {{math|''m''<sub>∞</sub>}} or {{math|''α''}} and {{math|''β''}} to define our gate fractions.
| |
| | |
| With such a form, all that remains is to individually investigate each current one wants to include. Typically, these include inward Ca<sup>2+</sup> and Na<sup>+</sup> input currents and several varieties of K<sup>+</sup> outward currents, including a "leak" current. The end result can be at the small end 20 parameters which one must estimate or measure for an accurate model, and for complex systems of neurons not easily tractable by computer. Careful simplifications of the Hodgkin-Huxley model are therefore needed.
| |
| | |
| ===FitzHugh-Nagumo===
| |
| {{main|FitzHugh-Nagumo model}}
| |
| Sweeping simplifications to Hodgkin-Huxley were introduced by FitzHugh and Nagumo in 1961 and 1962. Seeking to describe "regenerative self-excitation" by a nonlinear positive-feedback membrane voltage and recovery by a linear negative-feedback gate voltage, they developed the model described by
| |
| | |
| :<math>\begin{array}{rcl}
| |
| \dfrac{d V}{d t} &=& V-V^3 - w + I_\mathrm{ext} \\ \\
| |
| \tau \dfrac{d w}{d t} &=& V-a-b w
| |
| \end{array}</math>
| |
| | |
| where we again have a membrane-like voltage and input current with a slower general gate voltage {{math|''w''}} and experimentally-determined parameters {{math|''a'' {{=}} -0.7, ''b'' {{=}} 0.8, ''τ'' {{=}} 1/0.08}}. Although not clearly derivable from biology, the model allows for a simplified, immediately available dynamic, without being a trivial simplification.<ref>{{cite doi|10.4249/scholarpedia.1349}}</ref>
| |
| | |
| ===Morris-Lecar===
| |
| {{main|Morris-Lecar model}}
| |
| In 1981 Morris and Lecar combined Hodgkin-Huxley and FitzHugh-Nagumo into a voltage-gated calcium channel model with a delayed-rectifier potassium channel, represented by
| |
| | |
| :<math>\begin{array}{rcl} | |
| C\dfrac{d V}{d t} &=& -I_\mathrm{ion}(V,w) + I \\ \\
| |
| \dfrac{d w}{d t} &=& \phi \cdot \dfrac{w_{\infty} - w}{\tau_{w}}
| |
| \end{array}</math>
| |
| | |
| where <math>I_\mathrm{ion}(V,w) = \bar{g}_\mathrm{Ca} m_{\infty}\cdot(V-V_\mathrm{Ca}) + \bar{g}_\mathrm{K} w\cdot(V-V_\mathrm{K}) + \bar{g}_\mathrm{L}\cdot(V-V_\mathrm{L})</math>.<ref name="Koch98"/>
| |
| | |
| ===Hindmarsh-Rose===
| |
| {{main|Hindmarsh-Rose model}}
| |
| Building upon the FitzHugh-Nagumo model, Hindmarsh and Rose proposed in 1984 a model of neuronal activity described by three coupled first order differential equations:
| |
| | |
| :<math>\begin{array}{rcl}
| |
| \dfrac{d x}{d t} &=& y+3x^2-x^3-z+I \\ \\
| |
| \dfrac{d y}{d t} &=& 1-5x^2-y \\ \\
| |
| \dfrac{d z}{d t} &=& r\cdot (4(x + \tfrac{8}{5})-z)
| |
| \end{array}</math>
| |
| | |
| with {{math|''r''<sup>2</sup> {{=}} ''x''<sup>2</sup> + ''y''<sup>2</sup> + ''z''<sup>2</sup>}}, and {{math|''r'' ≈ 10<sup>-2</sup> }} so that the {{math|''z''}} variable only changes very slowly. This extra mathematical complexity allows a great variety of dynamic behaviors for the membrane potential, described by the {{math|''x''}} variable of the model, which include chaotic dynamics. This makes the Hindmarsh-Rose neuron model very useful, because being still simple, allows a good qualitative description of the many different patterns of the action potential observed in experiments.
| |
| | |
| ==Expanded neuron models==
| |
| While the success of integrating and kinetic models is undisputed, much has to be determined experimentally before accurate predictions can be made. The theory of neuron integration and firing (response to inputs) is therefore expanded by accounting for the nonideal conditions of cell structure.
| |
| | |
| ===Cable theory===
| |
| {{see also|Cable theory}}
| |
| Cable theory describes the [[dendrites|dendritic arbor]] as a cylindrical structure undergoing a regular pattern of [[Bifurcation theory|bifurcation]], like branches in a tree. For a single cylinder or an entire tree, the input conductance at the base (where the tree meets the cell body, or any such boundary) is defined as
| |
| | |
| :<math>G_{in} = \frac{G_\infty \tanh(L) + G_L}{1+(G_L / G_\infty )\tanh(L)}</math>,
| |
| | |
| where {{math|''L''}} is the electrotonic length of the cylinder which depends on its length, diameter, and resistance. A simple recursive algorithm scales linearly with the number of branches and can be used to calculate the effective conductance of the tree. This is given by
| |
| | |
| :<math>\,\! G_D = G_m A_D \tanh(L_D) / L_D</math> | |
| | |
| where {{math|''A''<sub>''D''</sub> {{=}} ''πld''}} is the total surface area of the tree of total length {{math|''l''}}, and {{math|''L''<sub>''D''</sub>}} is its total electrotonic length. For an entire neuron in which the cell body conductance is {{math|''G''<sub>''S''</sub>}} and the membrane conductance per unit area is {{math|''G''<sub>''md''</sub> {{=}} ''G''<sub>''m''</sub> / ''A''}}, we find the total neuron conductance {{math|''G''<sub>''N''</sub>}} for {{math|''n''}} dendrite trees by adding up all tree and soma conductances, given by
| |
| | |
| :<math>G_N = G_S + \sum_{j=1}^n A_{D_j} F_{dga_j}</math>,
| |
| | |
| where we can find the general correction factor {{math|''F''<sub>''dga''</sub>}} experimentally by noting {{math|''G''<sub>''D''</sub> {{=}} ''G''<sub>''md''</sub>''A''<sub>''D''</sub>''F''<sub>''dga''</sub>}}.
| |
| | |
| ===Compartmental models===
| |
| {{see also|Multi-compartment model}}
| |
| The cable model makes a number of simplifications to give closed analytic results, namely that the dendritic arbor must branch in diminishing pairs in a fixed pattern. A compartmental model allows for any desired tree topology with arbitrary branches and lengths, but makes simplifications in the interactions between branches to compensate. Thus, the two models give complementary results, neither of which is necessarily more accurate.
| |
| | |
| Each individual piece, or compartment, of a dendrite is modeled by a straight cylinder of arbitrary length {{math|''l''}} and diameter {{math|''d''}} which connects with fixed resistance to any number of branching cylinders. We define the conductance ratio of the {{math|''i''}}th cylinder as {{math|''B''<sub>''i''</sub> {{=}} ''G''<sub>''i''</sub> / ''G''<sub>∞</sub>}}, where <math>G_\infty=\tfrac{\pi d^{3/2}}{2\sqrt{R_i R_m}}</math> and {{math|''R''<sub>''i''</sub>}} is the resistance between the current compartment and the next. We obtain a series of equations for conductance ratios in and out of a compartment by making corrections to the normal dynamic {{math|''B''<sub>out,''i''</sub> {{=}} ''B''<sub>in,''i+1''</sub>}}, as
| |
| | |
| *<math>B_{\mathrm{out},i} = \frac{B_{\mathrm{in},i+1}(d_{i+1}/d_i)^{3/2} }{ \sqrt{R_{\mathrm{m},i+1}/R_{\mathrm{m},i}} }</math>
| |
| *<math>B_{\mathrm{in},i} = \frac{ B_{\mathrm{out},i} + \tanh X_i }{ 1+B_{\mathrm{out},i}\tanh X_i }</math>
| |
| *<math>B_\mathrm{out,par} = \frac{B_\mathrm{in,dau1} (d_\mathrm{dau1}/d_\mathrm{par})^{3/2}} {\sqrt{R_\mathrm{m,dau1}/R_\mathrm{m,par}}} + \frac{B_\mathrm{in,dau2} (d_\mathrm{dau2}/d_\mathrm{par})^{3/2}} {\sqrt{R_\mathrm{m,dau2}/R_\mathrm{m,par}}} + \ldots</math>
| |
| | |
| where the last equation deals with ''parents'' and ''daughters'' at branches, and <math>X_i = \tfrac{l_i \sqrt{4R_i}}{\sqrt{d_i R_m}}</math>. We can iterate these equations through the tree until we get the point where the dendrites connect to the cell body (soma), where the conductance ratio is {{math|''B''<sub>in,stem</sub>}}. Then our total neuron conductance is given by
| |
| | |
| :<math>G_N = \frac{A_\mathrm{soma}}{R_\mathrm{m,soma}} + \sum_j B_{\mathrm{in,stem},j} G_{\infty,j}</math>.
| |
| | |
| ===Synaptic transmission===
| |
| {{see also|Neurotransmission}}
| |
| The response of a neuron to individual neurotransmitters can be modeled as an extension of the classical Hodgkin-Huxley model with both standard and nonstandard kinetic currents. Four neurotransmitters primarily have influence in the CNS. [[AMPA receptor|AMPA/kainate receptors]] are fast [[Excitatory synapse|excitatory]] mediators while [[NMDA receptor]]s mediate considerably slower currents. Fast [[Inhibitory synapse|inhibitory]] currents go through [[GABAA receptor|GABA<sub>A</sub> receptors]], while [[GABAB receptor|GABA<sub>B</sub> receptors]] mediate by secondary ''G''-protein-activated potassium channels. This range of mediation produces the following current dynamics:
| |
| | |
| *<math>I_\mathrm{AMPA}(t,V) = \bar{g}_\mathrm{AMPA} \cdot [O] \cdot (V(t)-E_\mathrm{AMPA})</math>
| |
| *<math>I_\mathrm{NMDA}(t,V) = \bar{g}_\mathrm{NMDA} \cdot B(V) \cdot [O] \cdot (V(t)-E_\mathrm{NMDA})</math>
| |
| *<math>I_\mathrm{GABA_A}(t,V) = \bar{g}_\mathrm{GABA_A} \cdot ([O_1]+[O_2]) \cdot (V(t)-E_\mathrm{Cl})</math>
| |
| *<math>I_\mathrm{GABA_B}(t,V) = \bar{g}_\mathrm{GABA_B} \cdot \tfrac{[G]^n}{[G]^n+K_\mathrm{d}} \cdot (V(t)-E_\mathrm{K})</math>
| |
| | |
| where {{math|''ḡ''}} <!-- \bar{g} --> is the maximal<ref name="Hod52">{{cite pmid| 12991237 }}</ref><ref name="Koch98"/> conductance (around 1[[Siemens (unit)|S]]) and {{math|''E''}} is the equilibrium potential of the given ion or transmitter (AMDA, NMDA, [[Chlorine|Cl]], or [[Potassium|K]]), while {{math|[''O'']}} describes the fraction of receptors that are open. For NMDA, there is a significant effect of ''magnesium block'' that depends [[sigmoid]]ally on the concentration of intracellular magnesium by {{math|''B''(''V'')}}. For GABA<sub>B</sub>, {{math|[''G'']}} is the concentration of the ''G''-protein, and {{math|''K''<sub>d</sub>}} describes the dissociation of ''G'' in binding to the potassium gates.
| |
| | |
| The dynamics of this more complicated model have been well-studied experimentally and produce important results in terms of very quick [[synaptic plasticity|synaptic potentiation and depression]], that is, fast, short-term learning.
| |
| | |
| ===Other conditions===
| |
| The models above are still idealizations. Corrections must be made for the increased membrane surface area given by numerous dendritic spines, temperatures significantly hotter than room-temperature experimental data, and nonuniformity in the cell's internal structure.<ref name="Koch98"/> Certain observed effects do not fit into some of these models. For instance, the temperature cycling (with minimal net temperature increase) of the cell membrane during action potential propagation not compatible with models which rely on modeling the membrane as a resistance which must dissipate energy when current flows through it. The transient thickening of the cell membrane during action potential propagation is also not predicted by these models, nor is the changing capacitance and voltage spike that results from this thickening incorporated into these models. The action of some anesthetics such as inert gases is problematic for these models as well. New models, such as the [[soliton model]] attempt to explain these phenomena, but are less developed than older models and have yet to be widely applied.
| |
| | |
| ==References==
| |
| {{reflist}}
| |
| | |
| {{DEFAULTSORT:Biological Neuron Model}}
| |
| [[Category:Neuroscience]]
| |
| [[Category:Biophysics]]
| |
| [[Category:Computational neuroscience]]
| |