ModelDB is moving. Check out our new site at The corresponding page is

Schiz.-linked gene effects on intrinsic single-neuron excitability (Maki-Marttunen et al. 2016)

 Download zip file 
Help downloading and running models
Python scripts for running NEURON simulations that model a layer V pyramidal cell with certain genetic variants implemented. The genes included are obtained from genome-wide association studies of schizophrenia.
1 . Mäki-Marttunen T, Halnes G, Devor A, Witoelar A, Bettella F, Djurovic S, Wang Y, Einevoll GT, Andreassen OA, Dale AM (2016) Functional Effects of Schizophrenia-Linked Genetic Variants on Intrinsic Single-Neuron Excitability: A Modeling Study. Biol Psychiatry Cogn Neurosci Neuroimaging 1:49-59 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Neuron or other electrically excitable cell;
Brain Region(s)/Organism:
Cell Type(s): Neocortex L5/6 pyramidal GLU cell;
Channel(s): I Na,p; I Na,t; I L high threshold; I T low threshold; I K; I K,leak; I M; I h; I K,Ca; I A, slow; Ca pump;
Gap Junctions:
Gene(s): Nav1.1 SCN1A; Nav1.7 SCN9A; Cav3.3 CACNA1I; Cav1.3 CACNA1D; Cav1.2 CACNA1C; Kv2.1 KCNB1; HCN1;
Simulation Environment: NEURON; Python;
Model Concept(s): Coincidence Detection; Active Dendrites; Detailed Neuronal Models; Schizophrenia;
Implementer(s): Maki-Marttunen, Tuomo [tuomomm at];
Search NeuronDB for information about:  Neocortex L5/6 pyramidal GLU cell; I Na,p; I Na,t; I L high threshold; I T low threshold; I K; I K,leak; I M; I h; I K,Ca; I A, slow; Ca pump;
Ca_HVA.mod *
Ca_LVAst.mod *
CaDynamics_E2.mod *
epsp.mod *
Ih.mod *
Im.mod *
K_Pst.mod *
K_Tst.mod *
SK_E2.mod *
SKv3_1.mod *
# Python-scripts for reproducing Figures 2-6 in the article 
# "Functional effects of schizophrenia-linked genetic variants on intrinsic single-neuron excitability: A modeling study."
# The scripts can be run using Python 2.7.8 with an interface to NEURON 7.3.
# If LaTeX installed, change the variable 'useLatex' to True for cleaner figure texts.
# All scripts are open for distribution and reuse through Creative Commons Attribution 2.0 Generic license
# (CC BY 2.0). The NEURON interface commands use the cell model and simulation of BAC firing by Etay Hay (Hay et al. 2011)
# The included .asc and .hoc files are copied as such from ModelDB entry 139653 with the author's permission, and the .mod
# files are copied with minor changes (hard-coded parameters are replaced with parameters whose values may be changed).
# Tuomo Maki-Marttunen, Jan 2015
# To run the python scripts, save them to the same directory with the .mod files, and make sure directories
# 'morphologies' and 'models' exist (in more detail, make sure that files 'morphologies/cell1.asc',
# 'morphologies/cell2.asc', 'models/L5PC_template.hoc', and 'models/L5PCbiophys3.hoc' are accessible). Then,
# follow the procedures below:
# 1) Run the .mod file complier:


# Before it is possible to draw any of figures, the mutations have to be scaled. To do this,
# do one of the following:
# 2a)
# python
# python
# Running calculates the f-I curve and membrane potential limit cycle for the control neuron.
# Running finds the scaling factors for each variant separately.
# This is a computationally heavy task, so if you can parallelize it to different computing nodes,
# you can alternatively run on each machine the following:
# 2b)
# python i
# where i=0,...,57.
# After 2a) or 2b), run
# python
# to collect the scalings into a single file, scalings.sav
# 2c) Alternatively, download the file scalings.sav, where the scalings are calculated in advance.
# 3) You are now ready to draw Figures 2-5 from the article:


# On a tested standard computer, Figure 2 took around one minute to calculate, Figure 3
# was finished in 10 hours, Figure 4 in 15 minutes, and Figure 5 in five hours.
# 4) To draw Figure 6, first calculate the spiking thresholds for uniformly distributed
# synaptic inputs on apical dendrite:


# (The second script takes some time to calculate, it can be distributed to different
# computers as was done in 2b with
# 5) Now, Figure 6 can be drawn:


# This is the figure with heaviest computational load, due to the fine resolution of ISI.
# It took around 6 days to finish the Figure 6 on the tested machine.

Loading data, please wait...