neurons to networks - university of washingtonfaculty.washington.edu/etsb/amath342/materials... ·...

22
Neurons to networks

Upload: others

Post on 15-Aug-2020

2 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: Neurons to networks - University of Washingtonfaculty.washington.edu/etsb/AMATH342/materials... · max τ τ t t P st Pee − − ... If synapse receives input spikes at rate r, then

Neurons to networks

Page 2: Neurons to networks - University of Washingtonfaculty.washington.edu/etsb/AMATH342/materials... · max τ τ t t P st Pee − − ... If synapse receives input spikes at rate r, then

How do synapses transform inputs?

Page 3: Neurons to networks - University of Washingtonfaculty.washington.edu/etsb/AMATH342/materials... · max τ τ t t P st Pee − − ... If synapse receives input spikes at rate r, then

Excitatory synapse

Input spike

! Neurotransmitter release

binds to/opens Na channels

Change in synaptic conductance

! Na+ influx

! Depolarization due to EPSP (excitatory postsynaptic potential)

E.g. AMPA synapse

Vocab: Depolarization means make V less neg = more positive

Page 4: Neurons to networks - University of Washingtonfaculty.washington.edu/etsb/AMATH342/materials... · max τ τ t t P st Pee − − ... If synapse receives input spikes at rate r, then

Input spike

! Neurotransmitter release

binds to/opens K channels

Change in synaptic conductance

K+ leaves cell

! Hyperpolarization due to IPSP (inhibitory postsynaptic potential)

Inhibitory synapse

Vocab: hyperpolarization means make V more negative

Page 5: Neurons to networks - University of Washingtonfaculty.washington.edu/etsb/AMATH342/materials... · max τ τ t t P st Pee − − ... If synapse receives input spikes at rate r, then

Modeling a synaptic input to a (RC) neuron

V

Synaptic conductance

R

Page 6: Neurons to networks - University of Washingtonfaculty.washington.edu/etsb/AMATH342/materials... · max τ τ t t P st Pee − − ... If synapse receives input spikes at rate r, then

Modeling a synaptic input to a (RC) neuron

V

Synaptic conductance

R

Probability of transmitter release given an input spike

Probability of postsynaptic channel opening(= fraction of channels opened)

Page 7: Neurons to networks - University of Washingtonfaculty.washington.edu/etsb/AMATH342/materials... · max τ τ t t P st Pee − − ... If synapse receives input spikes at rate r, then

Basic synapse model

Assume Prel = 1 (for now)

What does a single spike input do to Ps?

Kinetic model:

sssss PPdtdP

βα −−= )1(

closed open

open closed

s

s

⎯→⎯

⎯→⎯

β

α

Opening rate Closing rate

Fraction of channels closed

Fraction of channels open

Where:

Page 8: Neurons to networks - University of Washingtonfaculty.washington.edu/etsb/AMATH342/materials... · max τ τ t t P st Pee − − ... If synapse receives input spikes at rate r, then

Synaptic filters

A difference of exponentials model better fits biological data for GABA,NMDA synapse types

Exponential

Difference of exponentials

Page 9: Neurons to networks - University of Washingtonfaculty.washington.edu/etsb/AMATH342/materials... · max τ τ t t P st Pee − − ... If synapse receives input spikes at rate r, then

Simplified synaptic models

Difference of exponentials:

Alpha function:

peak

t

peaks ettP τ

τ

⋅= const)(

)(const)( 21max

ττtt

s eePtP−−

−⋅=t

Ps

Page 10: Neurons to networks - University of Washingtonfaculty.washington.edu/etsb/AMATH342/materials... · max τ τ t t P st Pee − − ... If synapse receives input spikes at rate r, then

What happens with a sequence of input spikes?

Input

(Markram & Tsodyks, 1998)

Short-term facilitation

Short-term depression

• Biological synapses are dynamic!

Page 11: Neurons to networks - University of Washingtonfaculty.washington.edu/etsb/AMATH342/materials... · max τ τ t t P st Pee − − ... If synapse receives input spikes at rate r, then

Short-term synaptic plasticity: describe this via P_relRecall definition of synaptic conductance:

Idea: Specify how Prel changes as a

function of consecutive input spikes

srelss PPgg max,=

relrel

P PPdtdP

−= 0τBetween input spikes, Prel decays exponentially back to P0

relDrel PfP → depression: decrement Prel

If input spike:

Abbott et al 1997

Input

[sketch on board]

Page 12: Neurons to networks - University of Washingtonfaculty.washington.edu/etsb/AMATH342/materials... · max τ τ t t P st Pee − − ... If synapse receives input spikes at rate r, then
Page 13: Neurons to networks - University of Washingtonfaculty.washington.edu/etsb/AMATH342/materials... · max τ τ t t P st Pee − − ... If synapse receives input spikes at rate r, then

Lab exercise:

Write a code that implements the Abbott et al mechanism for synaptic depression.

Drive the synapse with spikes occurring regularly at 20 Hz, as in Fig. 1A of Abbott et al ’97. Can you reproduce that figure?

Hint: this should be a few lines of code.

Page 14: Neurons to networks - University of Washingtonfaculty.washington.edu/etsb/AMATH342/materials... · max τ τ t t P st Pee − − ... If synapse receives input spikes at rate r, then

Key result (a few lines of calculation, see (7) in paper [ESB notes]):

If synapse receives input spikes at rate r, then the steady state value of

P_rel (r) ~ 1/r

Impact of synaptic depression

Page 15: Neurons to networks - University of Washingtonfaculty.washington.edu/etsb/AMATH342/materials... · max τ τ t t P st Pee − − ... If synapse receives input spikes at rate r, then

Consequences of synaptic depression: steady state

depression

Input Rate

average release probability

average transmission rate

At steady state,

Constant synaptic input for Wide range of inputs!

Steady-state gain control

Page 16: Neurons to networks - University of Washingtonfaculty.washington.edu/etsb/AMATH342/materials... · max τ τ t t P st Pee − − ... If synapse receives input spikes at rate r, then

Consequences of synaptic depression: dynamic response

Steady-state transmission rates are similar for different rates

Transient inputs are amplified relative to steady-state inputs

Input rates r

In fact: an equal-percent change from baseline gives an equal transient response.Who cares? Abbott et al 97: Neuron gets inputs from 1000’s of upstream cells, each of which fires at 1-200 Hz. How can we be responsive to all?

Synaptic depression yields “gain control” to satisfy this. vs. Adaptation or inhibition, does so in an INPUT-SPECIFIC way!

Page 17: Neurons to networks - University of Washingtonfaculty.washington.edu/etsb/AMATH342/materials... · max τ τ t t P st Pee − − ... If synapse receives input spikes at rate r, then

Extending the model to include facilitationRecall definition of synaptic conductance:

srelss PPgg max,=

Between input spikes, Prel still decays exponentially back to P0

relDrel PfP →

)1( relFrelrel PfPP −+→

depression: decrement Prel

facilitation: increment Prel

If input spike:

Abbott et al 1997

Page 18: Neurons to networks - University of Washingtonfaculty.washington.edu/etsb/AMATH342/materials... · max τ τ t t P st Pee − − ... If synapse receives input spikes at rate r, then

Effects of synaptic facilitation & depression

facilitation depression

Input Rate Input Rateaverage release probability

average transmission rate

Page 19: Neurons to networks - University of Washingtonfaculty.washington.edu/etsb/AMATH342/materials... · max τ τ t t P st Pee − − ... If synapse receives input spikes at rate r, then

More detailed plasticity model

ER

I

• Incoming spike activates a fraction of recovered resources R to become effective (E)

• Amount of effective resources E govern size of gs(t)• Effective resources rapidly inactivate (msec)• Inactivated recover (100s of msec)• E(t) determines postsynaptic current

Markram and Tsodyks, PNAS 1997

Page 20: Neurons to networks - University of Washingtonfaculty.washington.edu/etsb/AMATH342/materials... · max τ τ t t P st Pee − − ... If synapse receives input spikes at rate r, then

Match between model and experiment

Page 21: Neurons to networks - University of Washingtonfaculty.washington.edu/etsb/AMATH342/materials... · max τ τ t t P st Pee − − ... If synapse receives input spikes at rate r, then

Synapses provide an additional layer of dynamics and computation, beyond that occurring in single neurons!

Summary:

Page 22: Neurons to networks - University of Washingtonfaculty.washington.edu/etsb/AMATH342/materials... · max τ τ t t P st Pee − − ... If synapse receives input spikes at rate r, then

Network computation via simplified “firing rate” models

xt ythtStimulus Function

W, connection weight matrix

W_ij = weight from j to i neuron (or neural population) ifires with rate

inputi = stimi(t) +X

j

Wijrj(t)<latexit sha1_base64="vWXcx+HQS5jBNPKg+Iuk0eE0UAI=">AAACEnicbVBNSwMxFMzWr1q/qh69BItQFcpWBPUgFL14rODaQrss2TTbpk12l+StUJb+By/+FS8eVLx68ua/Md32oK0DgWHmPV5m/FhwDbb9beUWFpeWV/KrhbX1jc2t4vbOvY4SRZlDIxGppk80EzxkDnAQrBkrRqQvWMMfXI/9xgNTmkfhHQxj5krSDXnAKQEjecUjHsYJeBxfYg1cerwMh/gYt3UivT5ueCnvj7Dy+kb2iiW7YmfA86Q6JSU0Rd0rfrU7EU0kC4EKonWrasfgpkQBp4KNCu1Es5jQAemylqEhkUy7aZZphA+M0sFBpMwLAWfq742USK2H0jeTkkBPz3pj8T+vlUBw7qZZbBbSyaEgERgiPC4Id7hiFMTQEEIVN3/FtEcUoWBqLJgSqrOR54lzUrmo2LenpdrVtI082kP7qIyq6AzV0A2qIwdR9Iie0St6s56sF+vd+piM5qzpzi76A+vzB/l3nIc=</latexit><latexit sha1_base64="vWXcx+HQS5jBNPKg+Iuk0eE0UAI=">AAACEnicbVBNSwMxFMzWr1q/qh69BItQFcpWBPUgFL14rODaQrss2TTbpk12l+StUJb+By/+FS8eVLx68ua/Md32oK0DgWHmPV5m/FhwDbb9beUWFpeWV/KrhbX1jc2t4vbOvY4SRZlDIxGppk80EzxkDnAQrBkrRqQvWMMfXI/9xgNTmkfhHQxj5krSDXnAKQEjecUjHsYJeBxfYg1cerwMh/gYt3UivT5ueCnvj7Dy+kb2iiW7YmfA86Q6JSU0Rd0rfrU7EU0kC4EKonWrasfgpkQBp4KNCu1Es5jQAemylqEhkUy7aZZphA+M0sFBpMwLAWfq742USK2H0jeTkkBPz3pj8T+vlUBw7qZZbBbSyaEgERgiPC4Id7hiFMTQEEIVN3/FtEcUoWBqLJgSqrOR54lzUrmo2LenpdrVtI082kP7qIyq6AzV0A2qIwdR9Iie0St6s56sF+vd+piM5qzpzi76A+vzB/l3nIc=</latexit><latexit sha1_base64="vWXcx+HQS5jBNPKg+Iuk0eE0UAI=">AAACEnicbVBNSwMxFMzWr1q/qh69BItQFcpWBPUgFL14rODaQrss2TTbpk12l+StUJb+By/+FS8eVLx68ua/Md32oK0DgWHmPV5m/FhwDbb9beUWFpeWV/KrhbX1jc2t4vbOvY4SRZlDIxGppk80EzxkDnAQrBkrRqQvWMMfXI/9xgNTmkfhHQxj5krSDXnAKQEjecUjHsYJeBxfYg1cerwMh/gYt3UivT5ueCnvj7Dy+kb2iiW7YmfA86Q6JSU0Rd0rfrU7EU0kC4EKonWrasfgpkQBp4KNCu1Es5jQAemylqEhkUy7aZZphA+M0sFBpMwLAWfq742USK2H0jeTkkBPz3pj8T+vlUBw7qZZbBbSyaEgERgiPC4Id7hiFMTQEEIVN3/FtEcUoWBqLJgSqrOR54lzUrmo2LenpdrVtI082kP7qIyq6AzV0A2qIwdR9Iie0St6s56sF+vd+piM5qzpzi76A+vzB/l3nIc=</latexit>

neuron (or neural population) ireceives input

ri(t)<latexit sha1_base64="TDxq+riMmEuvG1mwEN4UMd0gCh8=">AAAB7HicbVBNS8NAEJ34WetX1aOXxSLUS0lFUG9FLx4rmLbQhrLZbtq1m03YnQgl9D948aDi1R/kzX/jts1BWx8MPN6bYWZekEhh0HW/nZXVtfWNzcJWcXtnd2+/dHDYNHGqGfdYLGPdDqjhUijuoUDJ24nmNAokbwWj26nfeuLaiFg94DjhfkQHSoSCUbRSU/dEBc96pbJbdWcgy6SWkzLkaPRKX91+zNKIK2SSGtOpuQn6GdUomOSTYjc1PKFsRAe8Y6miETd+Nrt2Qk6t0idhrG0pJDP190RGI2PGUWA7I4pDs+hNxf+8TorhlZ8JlaTIFZsvClNJMCbT10lfaM5Qji2hTAt7K2FDqilDG1DRhlBbfHmZeOfV66p7f1Gu3+RpFOAYTqACNbiEOtxBAzxg8AjP8ApvTuy8OO/Ox7x1xclnjuAPnM8fYceOiA==</latexit><latexit sha1_base64="TDxq+riMmEuvG1mwEN4UMd0gCh8=">AAAB7HicbVBNS8NAEJ34WetX1aOXxSLUS0lFUG9FLx4rmLbQhrLZbtq1m03YnQgl9D948aDi1R/kzX/jts1BWx8MPN6bYWZekEhh0HW/nZXVtfWNzcJWcXtnd2+/dHDYNHGqGfdYLGPdDqjhUijuoUDJ24nmNAokbwWj26nfeuLaiFg94DjhfkQHSoSCUbRSU/dEBc96pbJbdWcgy6SWkzLkaPRKX91+zNKIK2SSGtOpuQn6GdUomOSTYjc1PKFsRAe8Y6miETd+Nrt2Qk6t0idhrG0pJDP190RGI2PGUWA7I4pDs+hNxf+8TorhlZ8JlaTIFZsvClNJMCbT10lfaM5Qji2hTAt7K2FDqilDG1DRhlBbfHmZeOfV66p7f1Gu3+RpFOAYTqACNbiEOtxBAzxg8AjP8ApvTuy8OO/Ox7x1xclnjuAPnM8fYceOiA==</latexit><latexit sha1_base64="TDxq+riMmEuvG1mwEN4UMd0gCh8=">AAAB7HicbVBNS8NAEJ34WetX1aOXxSLUS0lFUG9FLx4rmLbQhrLZbtq1m03YnQgl9D948aDi1R/kzX/jts1BWx8MPN6bYWZekEhh0HW/nZXVtfWNzcJWcXtnd2+/dHDYNHGqGfdYLGPdDqjhUijuoUDJ24nmNAokbwWj26nfeuLaiFg94DjhfkQHSoSCUbRSU/dEBc96pbJbdWcgy6SWkzLkaPRKX91+zNKIK2SSGtOpuQn6GdUomOSTYjc1PKFsRAe8Y6miETd+Nrt2Qk6t0idhrG0pJDP190RGI2PGUWA7I4pDs+hNxf+8TorhlZ8JlaTIFZsvClNJMCbT10lfaM5Qji2hTAt7K2FDqilDG1DRhlBbfHmZeOfV66p7f1Gu3+RpFOAYTqACNbiEOtxBAzxg8AjP8ApvTuy8OO/Ox7x1xclnjuAPnM8fYceOiA==</latexit>

⌧dridt

= f(inputi)� ri<latexit sha1_base64="DRH0qSCSPSgX/UqOZX8tqSzhi5w=">AAACEXicbVDLSsNAFJ3UV62vqEs3g0VoF5ZEBHUhFN24rGBtoS1hMpm0QycPZm6EEvINbvwVNy5U3Lpz5984TbPQ1gMDh3Pu5c45biy4Asv6NkpLyyura+X1ysbm1vaOubt3r6JEUtamkYhk1yWKCR6yNnAQrBtLRgJXsI47vp76nQcmFY/CO5jEbBCQYch9TgloyTHrfSAJ7vuS0NTD0uFZ6kGGL7FfwzyME3B4HR9PDcesWg0rB14kdkGqqEDLMb/6XkSTgIVABVGqZ1sxDFIigVPBsko/USwmdEyGrKdpSAKmBmkeKcNHWvGwH0n9QsC5+nsjJYFSk8DVkwGBkZr3puJ/Xi8B/3yQ5tFYSGeH/ERgiPC0H+xxySiIiSaESq7/iumI6HpAt1jRJdjzkRdJ+6Rx0bBuT6vNq6KNMjpAh6iGbHSGmugGtVAbUfSIntErejOejBfj3fiYjZaMYmcf/YHx+QOLjZxP</latexit><latexit sha1_base64="DRH0qSCSPSgX/UqOZX8tqSzhi5w=">AAACEXicbVDLSsNAFJ3UV62vqEs3g0VoF5ZEBHUhFN24rGBtoS1hMpm0QycPZm6EEvINbvwVNy5U3Lpz5984TbPQ1gMDh3Pu5c45biy4Asv6NkpLyyura+X1ysbm1vaOubt3r6JEUtamkYhk1yWKCR6yNnAQrBtLRgJXsI47vp76nQcmFY/CO5jEbBCQYch9TgloyTHrfSAJ7vuS0NTD0uFZ6kGGL7FfwzyME3B4HR9PDcesWg0rB14kdkGqqEDLMb/6XkSTgIVABVGqZ1sxDFIigVPBsko/USwmdEyGrKdpSAKmBmkeKcNHWvGwH0n9QsC5+nsjJYFSk8DVkwGBkZr3puJ/Xi8B/3yQ5tFYSGeH/ERgiPC0H+xxySiIiSaESq7/iumI6HpAt1jRJdjzkRdJ+6Rx0bBuT6vNq6KNMjpAh6iGbHSGmugGtVAbUfSIntErejOejBfj3fiYjZaMYmcf/YHx+QOLjZxP</latexit><latexit sha1_base64="DRH0qSCSPSgX/UqOZX8tqSzhi5w=">AAACEXicbVDLSsNAFJ3UV62vqEs3g0VoF5ZEBHUhFN24rGBtoS1hMpm0QycPZm6EEvINbvwVNy5U3Lpz5984TbPQ1gMDh3Pu5c45biy4Asv6NkpLyyura+X1ysbm1vaOubt3r6JEUtamkYhk1yWKCR6yNnAQrBtLRgJXsI47vp76nQcmFY/CO5jEbBCQYch9TgloyTHrfSAJ7vuS0NTD0uFZ6kGGL7FfwzyME3B4HR9PDcesWg0rB14kdkGqqEDLMb/6XkSTgIVABVGqZ1sxDFIigVPBsko/USwmdEyGrKdpSAKmBmkeKcNHWvGwH0n9QsC5+nsjJYFSk8DVkwGBkZr3puJ/Xi8B/3yQ5tFYSGeH/ERgiPC0H+xxySiIiSaESq7/iumI6HpAt1jRJdjzkRdJ+6Rx0bBuT6vNq6KNMjpAh6iGbHSGmugGtVAbUfSIntErejOejBfj3fiYjZaMYmcf/YHx+QOLjZxP</latexit>

DYNAMICS: rates approach steady states f(input)f(input)

input