Computing in carbon Basic elements of neuroelectronics -- membranes - - PDF document

computing in carbon
SMART_READER_LITE
LIVE PREVIEW

Computing in carbon Basic elements of neuroelectronics -- membranes - - PDF document

2/1/2017 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


slide-1
SLIDE 1

2/1/2017 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

Equivalent circuit model

slide-2
SLIDE 2

2/1/2017 2

Membrane patch

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

The passive membrane

slide-3
SLIDE 3

2/1/2017 3 Energetics: qV ~ kBT V ~ 25mV

Movement of ions through ion channels

Na+, Ca2+ K+

Ions move down their concentration gradient until opposed by electrostatic forces

Nernst:

The equilibrium potential

slide-4
SLIDE 4

2/1/2017 4

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

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

Parallel paths for ions to cross membrane

slide-5
SLIDE 5

2/1/2017 5

Neurons are excitable

  • Voltage dependent
  • transmitter dependent (synaptic)
  • Ca dependent

Excitability arises from ion channel nonlinearity

slide-6
SLIDE 6

2/1/2017 6

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

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-7
SLIDE 7

2/1/2017 7

We can rewrite: where

Dynamics of activation and inactivation Dynamics of activation and inactivation

slide-8
SLIDE 8

2/1/2017 8

  • Ohm’s law:

and Kirchhoff’s law

Capacitative current Ionic currents Externally applied current

Putting it together The Hodgkin-Huxley equation

slide-9
SLIDE 9

2/1/2017 9 Na ~ m3h K ~ m3h

EK ENa

Anatomy of a spike

EK ENa

Runaway +ve feedback Double whammy

Anatomy of a spike

slide-10
SLIDE 10

2/1/2017 10

Where to from here?

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

Ion channel stochasticity

slide-11
SLIDE 11

2/1/2017 11

approach to macroscopic description

Microscopic models for ion channel fluctuations

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-12
SLIDE 12

2/1/2017 12

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 V

Vmax Vreset Vth

f(V)

Vrest

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

Exponential integrate-and-fire neuron

slide-13
SLIDE 13

2/1/2017 13 Vth Vrest

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

The theta neuron

Ermentrout and Kopell Vspike

  • 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-14
SLIDE 14

2/1/2017 14

Two-dimensional models

V w

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

Truccolo and Brown, Paninski, Pillow, Simoncelli

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

The generalized linear model

slide-15
SLIDE 15

2/1/2017 15

Dendritic computation

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

Dendritic computation

slide-16
SLIDE 16

2/1/2017 16 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

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

Linear cables

slide-17
SLIDE 17

2/1/2017 17 ri rm cm For a length L of membrane cable: ri  ri L rm  rm / L cm  cm L

Axial and membrane resistance

(1) (2)

The cable equation

x+dx x

slide-18
SLIDE 18

2/1/2017 18 (1) (2) (1) 

  • r

where Time constant Space constant

The cable equation General solution: filter and impulse response

Exponential decay Diffusive spread

slide-19
SLIDE 19

2/1/2017 19

Current injection at x=0, T  ∞

Voltage decays exponentially away from source

 Electrotonic length

Properties of passive cables

slide-20
SLIDE 20

2/1/2017 20 Johnson and Wu

Electrotonic length

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

Properties of passive cables

slide-21
SLIDE 21

2/1/2017 21 Johnson and Wu  Current can escape through additional pathways: speeds up decay

Voltage rise time

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

Properties of passive cables

slide-22
SLIDE 22

2/1/2017 22  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 Step response: pulse travels

Conduction velocity

slide-23
SLIDE 23

2/1/2017 23

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

Conduction velocity

Finite cables Active channels

Other factors

slide-24
SLIDE 24

2/1/2017 24

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

New cable equation for each dendritic compartment

Active cables

slide-25
SLIDE 25

2/1/2017 25

Genesis, NEURON

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

Passive computations

London and Hausser, 2005

slide-26
SLIDE 26

2/1/2017 26

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