Modeling Dendrites and Synapses

In the last lecture we learned how to make a compartmental model of a neuron, including models of voltage activated channels. This time, we will talk about what goes on in the dendrites and how this relates to the behavior of more complex cells and networks. The first step is to model synaptically (chemically) activated channels.

The variable conductance in the generic compartment diagram could also represent a synaptically activated cconductance, usually in a dendritic compartment.

This diagram shows a synapse with an action potential propagating to the pre-synaptic terminal, causing neurotransmitter release, and the resulting post-synaptic conductance change, G(t). Current flow through this conductance can lead to a post-synaptic potential (PSP).

There is a lot of biochemistry and molecular biology involved in this behavior. Fortunately, we can often use an empirical fit to the observed behavior, rather than modeling it in detail.

Typically, the conductance change from a quantum of neurotransmitter follows a linear rise and exponential decay, so it is often modeled with a so-called "alpha" function with a single time constant, tau.

	G(t) = const * t * exp(-t/tau)

The current due to this conductance is I = G·(Ek - Vm), which may be into or out of the cell, depending on the size of the ionic equilbrium potential Ek, relative to the membrane potential. Here we have adopted the convention that a positive current flows into the cell. Thus, if the Ek is large (Na or Ca), the current will be into the cell, and it will be a depolarizing (excitatory) synapse. If it is large and negative (Ek < Vm, as for potassium), it will be an inhibitory synapse. So, we can use the same model for both types of synapses.

Other models use a difference of two exponentials, with two time constants. We'll use the alpha function form in a GENESIS tutorial simulation called "Neuron" in order to understand both temporal and spatial summation of synaptic inputs.

Although we won't do it here, this simulation will also let you experiment with "silent inhibition" (shunting inhibition), in which there can be an inhibitory effect when there is no PSP, or even a slightly depolarizing one. Do you know how this can happen?

Another thing that we can do in this simulation is to scale the conductance by a "synaptic weight" factor. This would let us experiment with the effect of multiple inputs in synchrony, or with learning (synaptic plasticity) by varying the weight.

Another time constant that is relevant for determining the properties of the cell is the membrane time constant, Rm·Cm. This represents the time that it takes to charge up and to discharge the membrane capacitance, and will also affect the duration of a post-synaptic potential. We could also use the Neuron simulation (but won't here) to see what happens if you have a small time constant for synaptic activation and a large membrane time constant.

Single Neuron Behavior - Temporal Summation

This simulation has a help menu that can show a schematic diagram of the model neuron and the types of inputs that we can give to it.

We have a simple model neuron with a soma and two dendrite compartments that are connected by axial resistances, as well as the option of putting any number of passive cable compartments between these two dendrite compartments. These passive compartments are like the ones shown, but have no variable conductances. This will let us see what happens if we have spatially separated inputs to the neuron. The soma has Hodgkin-Huxley voltage activated channels like the ones we used in the last simulation, and the dendrite compartments have both excitatory and inhibitory synaptically activated channels which respond to spikes applied at the synapses. Not only can we change the various parameters of the cell by popping up menus, but we can apply different types of stimulation. For example, we can inject pulses of current into any of the compartments, or connect spike trains to any of the synapses with a specified weighting. We'll start by applying a burst of spikes 10 milliseconds apart to the Excitatory synapse in the first dendrite section.

The spike train starts after 10 milliseconds. {UPPER WINDOW}

We can see that each spike causes a rise and decay in the conductance of this channel. {MIDDLE WINDOW}

However, the conductance decays before the postsynaptic potential can build up to a very high level. {LOWER WINDOW}

The persistence of a PSP is determined not only by the decay constant of the conductance, but also by Rm·Cm, the rate at which the voltage across the membrane capacitance is discharged through the membrane resistance. In this simulation, tau for the decay of the excitatory synaptic conductance is 3 msec, and Rm·Cm is 5 msec. Both of these are less than the spike interval. So, let's see what happens if we make the input spikes come four times as fast.

This is done in overlay mode, with the new results plotted in blue, so we can compare them with the previous results. (However, the plots of activation parameters and soma conductances were not overlayed, in order to keep the plots simple.)

Now, we are getting a good buildup of conductance as each increase adds to the previous one, and we see action potentials both in the dendrite compartment {LOWER LEFT} and in the soma {LOWER RIGHT}. Why do we see action potentials in the dendrite compartment? You might notice that the action potentials are a little higher in the soma. This is because the PSP in the dendrite is propagated to the soma, where the voltage activated channels cause the action potentials, and these are propagated back to the dendrite through an axial resistance. In the dendrite compartment, we're seeing a superposition of the small PSP produced here and the large action potential that is produced in the soma. If we were to look in a more distant dendrite section, the peaks would be much more attenuated. We'll do this in the next simulation run.

Single Neuron Behavior - Spatial Summation

Now we'll explore the effects of summing spatially separated inputs to the neuron. Although it isn't shown on this diagram, I've added 5 passive cable sections between the two dendrite compartments, so dendrite compartment #2 is much further from the soma. With the dimensions and resistances that are used here, cable theory predicts that a voltage in this compartment will be attenuated by about 50% by the time it gets to the soma. This distance is 0.7 of the "length constant", so the attenuation is exp(-0.7). Do you know how to calculate this? If you are interested, take a look at the optional material on Cable Theory of Passive Propagation in Dendrites.

You may want to consult the neuron inputs diagram to remind yourself of the location of the three compartments.

I've set the interval for both sources of spike trains to 2 msecs, but I'll start applying Source A to the excitatory input of dendrite #2 right away, and wait 60 msec before applying Source B to the inhibitory input of dendrite #1.

Since the excitatory input is further away, we'll give it a fighting chance by giving it a higher synaptic weight. The weight of 40 means that we are scaling the conductance as if 40 synapses were being made at this point in the dendrite, and they are all getting the same input. The inhibitory input to dendrite #1 gets half this weight. We use this synaptic weighting because the input from a single synapse is rarely enough to fire an action potential, unless the cell is a motor neuron.

The excitatory spike train (in red) starts right away (Inputs to cell - UPPER LEFT), and causes a large buildup of conductance in the excitatory channel of the distal dendrite compartment {UPPER RIGHT}. This leads to a large increase in the membrane potential in this compartment - way above threshold.

This passively propagates to the soma where it produces action potentials. These are shown in slightly attenuated form in the nearby dendrite #1 compartment (LOWER LEFT), and in much more attenuated form in the distant dendrite #2 compartment (LOWER RIGHT). In this compartment, we can see the postsynaptic potential clearly, with the attenuated somatic action potentials superimposed.

After about 60 msec, the action potentials are disappearing. Let's look at our inputs and see what's happening. Now we are getting spikes to the inhibitory channel dendrite #1 {UPPER LEFT}. It produces a fairly small buildup of conductance {MIDDLE LEFT}, which produces a modest hyperpolaring influence in the compartment {LOWER LEFT}. But, since this compartment is close to the soma, it has a stronger effect than the depolarizing PSP in the distal dendrite compartment, and it stops the production of action potentials in the soma.

If we like, we can change the display to plot the soma membrane potential and see this directly.

We've seen two properties of neurons that are relevant to understanding the ways which neurons compute:

Next, we will look at the behavior of more complex neurons, and conclude with a model of a cortical network.


Back to Lectures on Computational Neuroscience

GENESIS Home Page


Dave Beeman, University of Colorado, Boulder
Mon Nov 19 15:50:52 MST 2001