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

Simulations of modulation of HCN channels in L5PCs (Mäki-Marttunen and Mäki-Marttunen, 2022)

 Download zip file 
Help downloading and running models
"... In this work, we build upon existing biophysically detailed models of thick-tufted layer V pyramidal cells and model the effects of over- and under-expression of Ih channels as well as their neuromodulation by dopamine (gain of Ih function) and acetylcholine (loss of Ih function). We show that Ih channels facilitate the action potentials of layer V pyramidal cells in response to proximal dendritic stimulus while they hinder the action potentials in response to distal dendritic stimulus at the apical dendrite. We also show that the inhibitory action of the Ih channels in layer V pyramidal cells is due to the interactions between Ih channels and a hot zone of low voltage-activated Ca2+ channels at the apical dendrite. Our simulations suggest that a combination of Ih-enhancing neuromodulation at the proximal apical dendrite and Ih-inhibiting modulation at the distal apical dendrite can increase the layer V pyramidal excitability more than any of the two neuromodulators alone..."
1 . Mäki-Marttunen T, Mäki-Marttunen V (2022) Excitatory and inhibitory effects of HCN channel modulation on excitability of layer V pyramidal cells Plos Comp Biol [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 layer 5 pyramidal cell;
Gap Junctions:
Receptor(s): AMPA; NMDA; Gaba;
Transmitter(s): Acetylcholine; Dopamine; Glutamate; Gaba;
Simulation Environment: NEURON;
Model Concept(s): Neuromodulation;
Implementer(s): Maki-Marttunen, Tuomo [tuomomm at];
Search NeuronDB for information about:  AMPA; NMDA; Gaba; Acetylcholine; Dopamine; Gaba; Glutamate;
TITLE GABA receptor

Conductance-based GABA synaptic current. Copied from AMPANMDA.mod, retained only the AMPA-form current


        RANGE gGABAmax
	RANGE E_Cl, tau_sGABA
        RANGE i, i_GABA, g_GABA, sGABA



	gGABAmax = 0.01  (uS)

	E_Cl = -80       (mV)
	tau_sGABA = 2   (ms)



        v (mV)
        i (nA)
	i_GABA (nA)
        g_GABA (uS)


        sGABA       : GABA state variable to construct the single-exponential profile - decays with conductance tau_sGABA


	sGABA = 0


        SOLVE state METHOD cnexp
        g_GABA = gGABAmax*sGABA          :compute time varying conductance
        i_GABA = g_GABA*(v-E_Cl) :compute the GABA driving force based on the time varying conductance, membrane potential, and GABA reversal
	i = i_GABA


        sGABA' = -sGABA/tau_sGABA

:NET_RECEIVE (weight, Pv, Pr, u, tsyn (ms)){
NET_RECEIVE (weight){
        sGABA = sGABA + 1
        if (sGABA > 1) { :Do not allow larger values than 1
          sGABA = 1

Loading data, please wait...