### Slides - University of Reddit

```Neural Basis of Cognition
Lecture 7
Learning and Memory: Part III
LTP/LTD
• In the last lecture, we briefly explored a
candidate biochemical pathway for the induction
of LTP/LTD
• Questions:
– Can we model the induction of LTP/LTD?
– Is this pathway plausible? i.e. if we model this
pathway, will our model neuron exhibit LTP/LTD as
observed in real neurons?
– If a neural network is equipped with such rules, can it
“learn”?
LTP/LTD
• Can we model the induction of LTP/LTD? Yes.
• Doing so will require some background.
A Mathematical Detour
• Next, we are going to consider the HodgkinHuxley model of the neuron. In order to do
that, we will first review several topics:
– Ordinary differential equations
– Equilibrium potential
– Electrical circuits
• These tools will give us the ability to
understand the basics of the Hodgkin-Huxley
model
Differential Equations
• Given a function f(x), the average change in
the function over a time h starting at a point x
is “rise over run,” that is, the change in f(x)
divided by the change in x
– [f(x+h)-f(x)]/h
Differential Equations
• As we decrease h more and more until it
approaches 0, the secant line becomes a line
tangent to f(x) at the point x.
Differential Equations
• If we do not choose the point x at which we
want to find the tangent line of f(x) and leave
x as a variable, we can find the general
equation for the tangent line of f(x) at any
point x; this equation is the derivative of f(x),
denoted f’(x) or df(x)/dx
Differential Equations
• Example: f(x) = x^2
[f(x+h) – f(x)] / h
= [(x+h)^2 – x^2] / h
= [(x^2 + 2xh + h^2) – x^2] / h
= (2xh + h^2) / h
= 2x + h
lim (h->0) 2x+h = 2x
So, the derivative of x^2 is 2x: f’(x) = 2x
Differential Equations
• What does the derivative mean? The derivative of a
function describes how quickly the function changes
with time.
• If f(x) = x^2, f’(x) = 2x; this means that at any point x,
f(x) is increasing by an amount equal to twice that x
• Therefore, if we know that some quantity x increases
twice as fast as its current value x, we know that the
function that can give the value of that quantity at any
point in time has derivative 2x, and therefore the
function describing the change in that quantity over
time is x^2.
Differential Equations
• The process of finding a function from its
derivative is called antidifferentiation, or, more
commonly, integration.
• When describing physical systems, it is often
easier to observe how a system changes based on
its current state and to write that down an
equation, called a differential equation;
integrating this equation will lead to an equation
that can describe the state of the system at any
point in time.
Equilibrium potential
• The concentrations of an ion inside and outside of a cell are
maintained by a combination of electrical and chemical
• The Nernst equation can be used to calculate the potential
of an ion of charge z across a membrane. This potential is
determined using the concentration of the ion both inside
and outside the cell when it is at rest:
– E = (RT/zF) ln([outside]/[inside])
•
•
•
•
R = universal gas constant, 8.314 Joules / (K mol)
T = temperature in Kelvin
z = charge of ion
F = Faraday’s constant, 96485 Coulombs/mol
• When the membrane potential is equal to E, there is no net
flux of this ion across the membrane
Equilibrium Potential
• Using the Nernst equation, we can calculate
the equilibrium potentials of Na+ and K+ in a
neuron:
– If, at resting potential, K+ has concentrations 5mM
outside and 140 mM inside the cell,
• EK = RT/F*ln(5/140) = -0.089 V = -89 mV
– If, at resting potential, Na+ has concentrations 12
mM inside and 140mM outside the cell,
• ENa = RT/F*ln(140/12) = 0.066 V = 66 mV
Electrical Circuits
• An electrical circuit is a closed loop through
which current (electricity) runs.
– “The voltage between two points is a short name for
the electrical force that would drive an electric
current between those points.”
– “Electric current means, depending on the context, a
flow of electric charge (a phenomenon) or the rate of
flow of electric charge (a quantity). This flowing
electric charge is typically carried by moving electrons,
in a conductor such as wire….”
Electrical Circuits
• “A resistor is a two-terminal electronic component that
produces a voltage across its terminals that
is proportional to the electric current passing through
it.”
• “A capacitor is a passive electronic
component consisting of a pair
of conductors separated by a dielectric (insulator).
When a voltage exists across the conductors,
an electric field is present in the dielectric. This field
stores energy.”
– “An ideal capacitor is wholly characterized by a constant
capacitance C, defined as the ratio of charge ±Q on each
conductor to the voltage V between them: C = Q/V”
Hodgkin-Huxley Model
In 1952, Hodgkin, Huxley, and Katz, based on their
experimental work with the squid giant axon, published
their model describing how action potentials in neurons
are initiated and propagated. They received the 1963
Nobel prize in Physiology or Medicine for their work.
Their model has since been significantly expanded.
The Hodgkin-Huxley model is based on an
“equivalent circuit” for the neuron. The
membrane is modeled as a capacitor, the
voltage-gated channels as resistors, and the
membrane for each ion as voltage sources.
Hodgkin-Huxley Model
• The current through each of the
resistors (i.e. the current through
each of the types of channels) is
proportional to the difference
between the voltage drop across the
membrane and the relevant ion’s
reversal potential, (Vm – Ek), where k
is the channel type: ik = -gk(Vm-Ek), so
Im = -Σk gk(Vm-Ek)
• The proportionality factor gk is the
conductance of the channel (the
conductance of the resistor in the
circuit). However, the conductance is
generally not a constant and instead
changes nonlinearly with the
potential of the neuron.
Hodgkin-Huxley Model
•
•
•
Recall that the capacitance of a capacitor
(which is a constant) is described by the
equation C = Q/V, rearrange, and
differentiate with respect to time:
CmVm = Q
Cm dVm/dt = dQ/dt
and the change in charge over time is current, so
Cm dVm/dt = Im
and Im was defined above, so
Cm dVm/dt = -Σk gk(Vm-Ek)
This is the fundamental equation describing
the Hodgkin-Huxley equivalent circuit.
Further complexity arises from the fact that
gk is also dependent on voltage, which can
make this differential equation more difficult
to solve.
LTP/LTD
• Is this pathway plausible? Yes.
LTP/LTD
•
“Calcium Time Course as a Signal for Spike-Timing–Dependent Plasticity” by Rubin
et al. 2005
– Two-compartment model, for the dendrite and the soma
– “The mechanism for plasticity in our model involves a biophysically plausible calcium
detection system that responds to calcium and then changes the strength of the synapse
accordingly. In the model, three detector agents (P, A, V) respond to the instantaneous
calcium level in the spine compartment. The interactions of these three agents, together
with two others (D, B) that they influence, act to track the calcium time course in the
spine (see Fig. 3A). More specifically, different calcium time courses lead to different
time courses of P and D, which compete to influence a plasticity variable W. This variable
W is used as a measure of the sign and magnitude of synaptic strength changes from
baseline. Note that this scheme is significantly different from a detection of peak
calcium levels, in which the change in W would be determined by how large spine
calcium becomes during an appropriate set of spikes. The interactions between agents
within our detector system qualitatively resemble the pathways influencing the
regulation of Ca2/calmodulin-dependent protein kinase II (CaMKII).”
Rubin, J. E., R. C. Gerkin, et al. (2005). "Calcium time course as a signal for spike-timing-dependent plasticity." J Neurophysiol 93(5): 2600-13.
LTP/LTD
“Calcium has been proposed as a
postsynaptic signal underlying synaptic spiketiming-dependent plasticity (STDP). We
examine this hypothesis with computational
modeling based on experimental results from
hippocampal cultures… in which pairs and
triplets of pre- and postsynaptic spikes induce
potentiation and depression in a temporally
asymmetric way. Specifically, we present a set
of model biochemical detectors, based on
plausible molecular pathways, which make
direct use of the time course of the calcium
signal to reproduce these experimental STDP
results…. Simulations of our model are also
shown to be consistent with classical LTP and
LTD induced by several presynaptic
Rubin, J. E., R. C. Gerkin, et al. (2005). "Calcium time course as a signal for spike-timing-dependent plasticity." J Neurophysiol 93(5): 2600-13.
LTP/LTD
• If a neural network is equipped with such rules, can it
“learn”? Yes.
LTP/LTD
•
“Learning Real-World Stimuli in a Neural Network with Spike-Driven Synaptic
Dynamics” by Brader, Senn, Fusi 2007
– “We present a model of spike-driven synaptic plasticity inspired by experimental observations
and motivated by the desire to build an electronic hardware device that can learn to classify
complex stimuli in a semisupervised fashion. During training, patterns of activity are
sequentially imposed on the input neurons, and an additional instructor signal drives the
output neurons toward the desired activity. The network is made of integrate-and-fire neurons
with constant leak and a floor. The synapses are bistable, and they are modified by the arrival
of presynaptic spikes. The sign of the change is determined by both the depolarization and the
state of a variable that integrates the postsynaptic action potentials. Following the training
phase, the instructor signal is removed, and the output neurons are driven purely by the
activity of the input neurons weighted by the plastic synapses. In the absence of stimulation,
the synapses preserve their internal state indefinitely.Memories are also very robust to the
disruptive action of spontaneous activity. A network of 2000 input neurons is shown to be able
to classify correctly a large number (thousands) of highly overlapping patterns (300 classes of
preprocessed Latex characters, 30 patterns per class, and a subset of the NIST characters data
set) and to generalize with performances that are better than or comparable to those of
artificial neural networks. Finally we show that the synaptic dynamics is compatible with many
of the experimental observations on the induction of long-term modifications (spike-timingdependent plasticity and its dependence on both the postsynaptic depolarization and the
frequency of pre- and postsynaptic neurons).”
Brader, J. M., W. Senn, et al. (2007). "Learning real-world stimuli in a neural network with spike-driven synaptic dynamics." Neural Comput 19(11): 2881-912.
LTP/LTD
“A network of 2000 input neurons is
shown to be able to classify
correctly a large number
(thousands) of highly overlapping
patterns (300 classes of
preprocessed Latex characters, 30
patterns per class, and a subset of
the NIST characters data set) and to
generalize with performances that
are better than or comparable to
those of artificial neural networks.”
Brader, J. M., W. Senn, et al. (2007). "Learning real-world stimuli in a neural network with spike-driven synaptic dynamics." Neural Comput 19(11): 2881-912.
LTP/LTD
“Figure 6: The full Latex data set containing 293
classes. (a) Percentage of nonclassified patterns in the
training set as a function of the number of classes for
different numbers of output units per class. Results
are shown for 1(∗), 2, 5, 10, and 40(+) outputs per
class using the abstract rule (points) and for 1, 2, and
5 outputs per class using the spike-driven network
(squares, triangles, and circles, respectively). In all
cases, the percentage of misclassified patterns is less
than 0.1%. (b) Percentage of nonclassified patterns as
a function of the number of output units per class for
different numbers of classes (abstract rule). Note the
logarithmic scale. (c) Percentage of nonclassified
patterns as a function of number of classes for
generalization on a test set (abstract rule). (d) Same
as for c but showing percentage of misclassified
patterns.”
Brader, J. M., W. Senn, et al. (2007). "Learning real-world stimuli in a neural network with spike-driven synaptic dynamics." Neural Comput 19(11): 2881-912.
```