A Simple Model
We start with the most basic model possible so we can focus on introducing you to the core components of the NEURON application and navigating the graphical user interface (GUI).
Part A: “Ball and Stick” model
To introduce the sequence of events involved in setting up a model, we will create an extremely simple neuron consisting solely of a soma and a dendrite called a “Ball and Stick” model.
STEP A1: Launch CellBuilder
CellBuilder is a graphical interface used to generate the basic morphology and biophysics of a neuron.
From the Build menu on the Main Menu window, select CellBuilder and read the descriptions provided.
STEP A2: Create morphology with CellBuilder
First we will specify the two components of our “ball and stick” model: a soma and a dendrite.
- First, let’s click on Topology
- The components will be displayed graphically in the white canvas region - a
soma
is created automatically - Create the dendrite component by clicking anywhere on the canvas once
- This generates a
dend
section as specified by Basename which is set by default todend
(Hint: click the button to change)
STEP A3: Specify morphological dimensions
The dimensions of the components we created can now be defined in the Geometry window.
NEURON refers to each component as a section. Each section is divided into one or more segments. A segment is the core fundamental unit consisting of a cylinder of uniform electrotonic characteristics based on cable theory.
The soma
and dendrite
are each cylindrical sections of uniform dimensions. We can now specify exactly what dimensions we would like them to have.
Section | Ref | Length (µm) | Diameter (µm) |
---|---|---|---|
Soma | soma |
20 | 20 |
Dendrite | dend |
1000 | 5 |
- Click on Geometry
all
should be selected by default - if not, selectall
in the side panel- Ensure Specify Strategy is ticked.
- Select length
L
, diameterdiam
from Distinct values over subset - Select
d_lambda
from Spatial grid (This parameter ensures segments are dynamically distributed - see ref 1).
What is Spatial Grid?
The Spatial Grid here refers to your model and is a computational issue rather than a biological one. A section is divided into one or more segments of uniform conductance properties. A finer Spatial Grid will therefore produce a more accurate result but this will also increase computational time.
The Spatial Grid can be defined three ways:
- The
nseg
button sets the number of segments. More info... - The
d_X
parameter specifies the spacing of segments by physical length, in µm - The
d_lambda
parameter specifies the spacing of segments by electrotonic length, expressed as a fraction of the AC length constant at 100 Hz for a cylindrical cable with the same diameter, Ra, and cm. This is often the best choice to use.
It is best practice to make the number of segments equal an odd number with nseg
. Why? This will ensure there is always a segment at the centre of the section (0.5).
nseg | section centre |
---|---|
1 | 0.5 |
2 | 0.33, 0.66 |
3 | 0.25, 0.5, 0.75 |
If your code calls for a segment that doesn't exist, then NEURON will round to the nearest segment - this may impact the accuracy of your results.
Now, we need to specify our lengths and diameters:
- Unselect Specify Strategy
d-lambda
is by default set to 0.1 which is suitable for most models- In the number box next to
soma.L (um)
, enter20
as per our table above - In the number box next to
dend.L (um)
, enter1000
as per our table above - In the number box next to
soma.diam (um)
, enter20
as per our table above - In the number box next to
dend.diam (um)
, enter5
as per our table above - so that the final result looks like:
Ball and Stick
We have now created a neuron with a soma that is a cylinder 20µm long, 20µm wide, with one dendrite protruding out 1000µm long, 5µm wide. Certainly lives up to the name "Ball and Stick"!
STEP A4: Specify biophysical characteristics
The next part of defining the neuron is to specify the biophysical characteristics such as the electrical resistance within the cytoplasm (Ra
), the capacitance of the plasma membrane (cm
), membrane ion channels, buffers and pumps. NEURON specifies ion channel density by setting the maximum combined conductance (in Siemens) of those channels in that particular segment.
- Click on Biophysics
- Ensure Specify Strategy is selected.
Ra
andcm
are uniform in this particular model, so we select theall
subset and then click on theRa
andcm
checkboxes.- We will add a Hodgkin-Huxley (
hh
) conductance tosoma
anddend
by selecting each then click onhh
Now we will add our values:
- Now deselect Specify Strategy
- Select
Ra
underall
and enter160
- The default for
cm
is 1 µF/cm^2 which is fine
Inserting ionic conductances
The generic hh
conductance is actually composed of three conductances:
- voltage-gated Na+ channel component
- voltage-gated K+ channel component
- the combined passive leak conductance
Note: If you have no hh
you will need to add the passive leak conductance by selecting pas
in specify strategy.
We will now insert HH into the soma and dendrite.
For the dendrite, we will enter a reduced hh
which means altering the standard HH conductances to be 10% of their initial values.
- Select
hh
underdend
and enter 10% of the Na+ (gnabar_hh
) and K+ (gnakbar_hh
) conductance values - No change to leak current of HH (
gl_hh
) - For the equilibrium potential (
el_hh
), change this to-64mV
Save Me
You can save this to a file called bs_cell.ses.
- Set your working directory using the File -> working dir command from the Main Menu window
- Enter the file path in the window or double-click on folders as they appear. The
../
will return you to the parent directory. - Select Move To which sets you up in the correct directory.
- Under the File -> Save session command, enter the name of the file in the top window (the bottom window is used as a filter only)
STEP A5: Running a simple simulation
We are now ready to load the specifications of our model into the NEURON simulator. We will use the Continuous Create method which is very useful when initially testing a model as it updates parameters in real time.
- Click Continuous Create in the CellBuilder window
- We will need to create a stimulus which is known as a Point Process so from the Main Menu window, select Tools -> Point Processes -> Managers -> Point Manager
- In the PointProcessManager window, click on SelectPointProcess then select IClamp
- We will accept the default of
soma(0.5)
which means the stimulus has been placed in the middle of the soma - We will insert a
0.6nA
current of1ms
pulse width starting at del=5ms
(a delay allows for initialization of any conductances with dependencies) so enter the values as shown.
- Now from the Main Menu window, select Tools->RunControl. This is our stimulus parameter window and is where the simulation is launched.
Run control window parameters
- The following options can be set:
Init(mV)
: determines the voltage we start at (generally keep this at the resting membrane potential you will expect from your ion channels you have placed in – the further away it is from that, the longer your cell will take to reach an equilibrium at the start-
Tstop
: controls the duration of our experiment, and -
dt
andPoints plotted/ms
: controls our temporal resolution of the experiments.
- We will accept the default values, except for the time, so enter
20
forTstop(ms)
. - To see an output of the simulation, from the Main Menu, select Graph -> Voltage axis
- Now in the RunControl window, click Init & Run
Have a look in the graph and you should see your neuron respond to the current.
Now we will rerun the simulator to see the effect of increasing the current injection amplitude.
- In the PointProcess window, change
0.6
to1.0
in theamp(nA)
- Now rerun in the RunControl window, by clicking Init & Run and we get a lovely AP.
STEP A6 Save Me
You can save this to a file called bs_iclamprig.ses. (Saving to a new file separates our simulation from our cell.)
Part B: Simple neuron model
We have covered the basics of creating a simple model and simulation. We can now turn our “Ball and Stick” model into a “Simple neuron” by replacing our single dendrite with a tree of branching apical dendrites with the following characteristics:
Section | Ref | Length (µm) | Diameter (µm) | Biophysics |
---|---|---|---|---|
Apical dendrite | ap |
400 | 2 | reduced hh |
Apical dendrite branch | ap[1] |
300 | 1 | reduced hh |
Apical dendrite branch | ap[2] |
300 | 1 | reduced hh |
STEP B1: Reloading sessions
A saved session can now be reloaded with all the parameters and windows already set.
- Close down the NEURON application via File -> Quit
- Restart with
nrngui
- You may need to return to your working directory with File -> recent dir
- Under File -> Load session, select the bs_cell.ses session created above
STEP B2: Modifying neurons
We will now modify the topology of the “ball and stick” neuron.
- From the CellBuilder, select Topology
- Select Delete Section then click on
dend
to remove this
- Click on Basename and type in
ap
and Accept - Click on Make Section
- Create two apical dendrite components by clicking anywhere on the canvas (twice!)
- Add a third dendrite as a branch by click and hold from the joint of the first two dendrites, then drag and release. Errors can be rectified by selecting the commands as required.
- This generates
ap
(root),ap[1]
andap[2]
STEP B3: Creating subsets
- As they share common properties, we will group the dendrite components so click on Subsets
- Click on Select Subtree
- Click on the
ap
segment (root of the subtree) - Click on New Section List
- In the popup window, enter
apicals
and Accept - This new list should appear in the sets.
STEP B4: Modifying parameters - Geometry
We will need to specify the morphological characteristics of the new dendrites
- Click on Geometry
- Ensure Specify Strategy is unticked
- Enter the lengths and diameters of the
ap
sections as per the table
STEP B5: Modifying parameters - Biophysics
For all the dendrites, we need to enter a reduced HH which means altering the standard HH conductances to be 10% of their initial values.
- Click on Biophysics
- Ensure Specify Strategy is ticked
- Select
apicals
then tickhh
- Now untick Specify Strategy
- Select
hh
underapicals
and enter 10% of the Na+ (gnabar_hh
) and K+ (gnakbar_hh
) conductance values - No change to leak current of HH (
gl_hh
) - For the equilibrium potential (
el_hh
), change this to-64mV
STEP B6: Save Me
You can save this to a file called simple_cell.ses.
B. EXPERIMENTS
Now it’s time for you to complete the simple neuron. To check you’ve got the hang of it, go back to the CellBuilder and add a basal dendrite and an axon with the following characteristics:
Section | Ref | Length (µm) | Diameter (µm) | Biophysics |
---|---|---|---|---|
Basal dendrite | bas |
200 | 3 | pas (e_pas =-65mV) |
Axon | axon |
800 | 1 | hh |
Don’t forget to save your changes to simple_cell.ses