Computing in carbon Basic elements of neuroelectronics -- membranes - - PowerPoint PPT Presentation

computing in carbon
SMART_READER_LITE
LIVE PREVIEW

Computing in carbon Basic elements of neuroelectronics -- membranes - - PowerPoint PPT Presentation

Computing in carbon Basic elements of neuroelectronics -- membranes -- ion channels -- wiring Elementary neuron models -- conductance based -- modelers alternatives Wires -- signal propagation -- processing in dendrites Wiring neurons


slide-1
SLIDE 1

Basic elements of neuroelectronics Elementary neuron models

  • - conductance based
  • - modelers’ alternatives
  • - membranes
  • - ion channels
  • - wiring

Wires

  • - signal propagation
  • - processing in dendrites

Wiring neurons together

  • - synapses
  • - long term plasticity
  • - short term plasticity

Computing in carbon

slide-2
SLIDE 2

Equivalent circuit model

slide-3
SLIDE 3

Membrane patch

slide-4
SLIDE 4

Ohm’s law: Capacitor: C = Q/V Kirchhoff:

The passive membrane

slide-5
SLIDE 5

Energetics: qV ~ kBT V ~ 25mV

Movement of ions through ion channels

slide-6
SLIDE 6

Na+, Ca2+ K+

Ions move down their concentration gradient until opposed by electrostatic forces

Nernst:

The equilibrium potential

slide-7
SLIDE 7

Different ion channels have associated conductances. A given conductance tends to move the membrane potential toward the equilibrium potential for that ion V > E  positive current will flow outward V < E  positive current will flow inward ENa ~ 50mV ECa ~ 150mV EK ~

  • 80mV

ECl ~

  • 60mV

depolarizing depolarizing hyperpolarizing shunting

V

Vrest ENa EK more polarized

Each ion type travels through independently

slide-8
SLIDE 8

Several I-V curves in parallel: New equivalent circuit:

Parallel paths for ions to cross membrane

slide-9
SLIDE 9

Neurons are excitable

slide-10
SLIDE 10
  • Voltage dependent
  • transmitter dependent (synaptic)
  • Ca dependent

Excitability arises from ion channel nonlinearity

slide-11
SLIDE 11

Persistent conductance

K channel: open probability increases when depolarized

PK ~ n4 n is open probability 1 – n is closed probability Transitions between states

  • ccur at voltage dependent

rates

C  O O  C

n describes a subunit

The ion channel is a cool molecular machine

slide-12
SLIDE 12

Gate acts as in previous case

PNa ~ m3h

Additional gate can block channel when open

m and h have opposite voltage dependences: depolarization increases m, activation hyperpolarization increases h, deinactivation m is activation variable h is inactivation variable

Transient conductances

slide-13
SLIDE 13

We can rewrite: where

Dynamics of activation and inactivation

slide-14
SLIDE 14

Dynamics of activation and inactivation

slide-15
SLIDE 15
  • Ohm’s law:

and Kirchhoff’s law

Capacitative current Ionic currents Externally applied current

Putting it together

slide-16
SLIDE 16

The Hodgkin-Huxley equation

slide-17
SLIDE 17

Na ~ m3h K ~ m3h

EK ENa

Anatomy of a spike

slide-18
SLIDE 18

EK ENa

Runaway +ve feedback Double whammy

Anatomy of a spike

slide-19
SLIDE 19

Where to from here?

Hodgkin-Huxley Biophysical realism Molecular considerations Geometry Simplified models Analytical tractability

slide-20
SLIDE 20

Ion channel stochasticity

slide-21
SLIDE 21

approach to macroscopic description

Microscopic models for ion channel fluctuations

slide-22
SLIDE 22

Different from the continuous model: interdependence between inactivation and activation transitions to inactivation state 5 can occur only from 2,3 and 4 k1, k2, k3 are constant, not voltage dependent

Transient conductances

slide-23
SLIDE 23

Like a passive membrane: but with the additional rule that when V  VT, a spike is fired and V  Vreset. EL is the resting potential of the “cell”.

The integrate-and-fire neuron

slide-24
SLIDE 24

V

Vmax Vreset Vth

f(V)

Vrest

f(V) = -V + exp([V-Vth]/D)

Exponential integrate-and-fire neuron

slide-25
SLIDE 25

Vth Vrest

dq/dt = 1 – cos q + (1+ cos q) I(t)

The theta neuron

Ermentrout and Kopell Vspike

slide-26
SLIDE 26
  • determine f from the linearized HH equations
  • fit a threshold
  • paste in the spike shape and AHP

Kernel f for subthreshold response  replaces leaky integrator Kernel for spikes  replaces “line”

Gerstner and Kistler

The spike response model

slide-27
SLIDE 27

Two-dimensional models

V w

Simple™ model: V’ = -aV + bV2 - cW W’ = -dW + eV

slide-28
SLIDE 28

Truccolo and Brown, Paninski, Pillow, Simoncelli

  • general definitions for k and h
  • robust maximum likelihood fitting procedure

The generalized linear model

slide-29
SLIDE 29

Dendritic computation

slide-30
SLIDE 30

Passive contributions to computation Active contributions to computation Dendrites as computational elements: Examples

Dendritic computation

slide-31
SLIDE 31

r Vm = Im Rm Current flows uniformly out through the cell: Im = I0/4pr2 Input resistance is defined as RN = Vm(t∞)/I0 = Rm/4pr2 Injecting current I0

Geometry matters

slide-32
SLIDE 32

rm and ri are the membrane and axial resistances, i.e. the resistances of a thin slice of the cylinder

Linear cables

slide-33
SLIDE 33

ri rm cm For a length L of membrane cable: ri  ri L rm  rm / L cm  cm L

Axial and membrane resistance

slide-34
SLIDE 34

(1) (2)

The cable equation

x+dx x

slide-35
SLIDE 35

(1) (2) (1) 

  • r

where Time constant Space constant

The cable equation

slide-36
SLIDE 36

General solution: filter and impulse response

Exponential decay Diffusive spread

slide-37
SLIDE 37

Current injection at x=0, T  ∞

Voltage decays exponentially away from source

slide-38
SLIDE 38

 Electrotonic length

Properties of passive cables

slide-39
SLIDE 39

Johnson and Wu

Electrotonic length

slide-40
SLIDE 40

 Electrotonic length Current can escape through additional pathways: speeds up decay

Properties of passive cables

slide-41
SLIDE 41

Johnson and Wu Current can escape through additional pathways: speeds up decay

Voltage rise time

slide-42
SLIDE 42

 Electrotonic length Current can escape through additional pathways: speeds up decay  Cable diameter affects input resistance

Properties of passive cables

slide-43
SLIDE 43

 Electrotonic length Current can escape through additional pathways: speeds up decay  Cable diameter affects input resistance  Cable diameter affects transmission velocity

Properties of passive cables

slide-44
SLIDE 44

Step response: pulse travels

Conduction velocity

slide-45
SLIDE 45

www.physiol.usyd.edu/au/~daved/teaching/cv.html

Conduction velocity

slide-46
SLIDE 46

Finite cables Active channels

Other factors

slide-47
SLIDE 47

Impedance matching: If a3/2 = d1

3/2 + d2 3/2

can collapse to an equivalent cylinder with length given by electrotonic length

Rall model

slide-48
SLIDE 48

New cable equation for each dendritic compartment

Active cables

slide-49
SLIDE 49

Genesis, NEURON

Who’ll be my Rall model, now that my Rall model is gone

slide-50
SLIDE 50

Passive computations

London and Hausser, 2005

slide-51
SLIDE 51

Enthusiastically recommended references

  • Johnson and Wu, Foundations of Cellular Physiology, Chap 4

The classic textbook of biophysics and neurophysiology: lots of problems to work through. Good for HH, ion channels, cable theory.

  • Koch, Biophysics of Computation

Insightful compendium of ion channel contributions to neuronal computation

  • Izhikevich, Dynamical Systems in Neuroscience

An excellent primer on dynamical systems theory, applied to neuronal models

  • Magee, Dendritic integration of excitatory synaptic input,

Nature Reviews Neuroscience, 2000 Review of interesting issues in dendritic integration

  • London and Hausser, Dendritic Computation,

Annual Reviews in Neuroscience, 2005 Review of the possible computational space of dendritic processing