KM_MA2020_GrC#
- class braincell.channel.KM_MA2020_GrC(size, g_max=Quantity(0.25, 'mS / cm^2'), temp=Quantity(303.15, 'K'), name=None)#
M-type potassium current of the granule cell model.
Slow, non-inactivating M-type potassium current imported from the cerebellar granule cell model of Masoli et al. (2020) [2]. A single first-order
ngate of power 1 drives an ohmic current:\[\begin{split}\begin{aligned} n_\infty &= \frac{1}{1 + \exp(-(V + 35) / 6)} \\ \alpha_n &= 0.0033 \, \exp((V + 30) / 40) \\ \beta_n &= 0.0033 \, \exp(-(V + 30) / 20) \\ \tau_n &= \frac{1}{\alpha_n + \beta_n} \end{aligned}\end{split}\]where \(V\) is in millivolts, the rates are per millisecond and \(\tau_n\) is in milliseconds. This class applies no voltage shift, and the reversal potential comes from the potassium ion object rather than from the class.
- Parameters:
size (
int|Sequence[int] |integer|Sequence[integer]) – Channel state shape.g_max (
Array|ndarray|bool|number|bool|int|float|complex|Quantity|Callable) – Maximal conductance density. Defaults to0.25 mS/cm2, which is exactly the source mechanism’sgkbar = 0.00025 mho/cm2.temp (
Array|ndarray|bool|number|bool|int|float|complex|Quantity) – Absolute temperature driving the gate’s q10 factor, default 30 degrees Celsius. This matches thecelsius = 30 (degC)written in the source mechanism’sPARAMETERblock.
See also
KM_MA2020_GoCGolgi-cell port of the same mechanism.
KM_RI2021_SCStellate-cell port of the same mechanism.
Kir2p3_MA2020_GrCInward rectifier of the same granule cell model, sharing the same origin paper.
Notes
Ported from
GrC/channel/KM_MA20_GrC.mod. That file, the Golgi portGoC/channel/KM_MA20_GoC.modand the stellate portSC/channel/KM_RI21_SC.modare byte-identical apart from theirSUFFIXline, and so are the three BrainCell classes: the rate constants above are shared verbatim withKM_MA2020_GoCandKM_RI2021_SC. What differs is only the deposit each was imported from, and therefore the model paper cited below.The mechanism does not use the steady state implied by its own rates. Its
n_inf = a_n/(a_n + b_n)line is commented out in the.modsource and replaced by the explicit Boltzmann shown above, so \(n_\infty\) and \(\tau_n\) are independent expressions here.Where the q10 factor is applied. The gate declares
q10 = 3.0at a reference of 22 degrees Celsius, soHH.compute_derivative()scales the whole \((n_\infty - n)/\tau_n\) term by \(\phi = 3^{(T - 22)/10}\) (about 2.41 at the default 30 degrees Celsius). The.modfile instead multipliesQ10intoalp_nandbet_n, which divides itstau_nby the same factor. The two forms are algebraically identical, but it meansf_n_tau()returns the q10-free time constant rather than the mechanism’stau_n.Provenance, and what the header does not establish. The
.modheader carriesAuthor: A. FontanaandCoAuthor: T.Nieus. That credit line is copy-pasted verbatim across every cell-type port of this mechanism and names people unrelated to the granule-cell key, so it is not treated as a citation here. The kinetics originate in the cerebellar granule cell model of D’Angelo et al. (2001) [1]; the granule-cell paper [2] names the model this parameterisation was imported from, not the origin of the equations. Reference [2] is the granule-cell paper specifically – the companion Golgi-cell paper of the same year belongs toKM_MA2020_GoC, not to this class.Conductance default.
0.25 mS/cm2is the deposit’s tuned value, carried across from the.modfile. It is not a value printed in either cited paper.Import deviations. The original mechanism’s NMODL
TABLEover[-100, 30] mV, coveringn_infandtau_n, is not reproduced: both expressions are evaluated per call. NEURON clamped tabulated values to the boundary outside that window, so any BrainCell-versus-NEURON divergence below -100 mV or above 30 mV is expected rather than a port error. The integration method was also changed fromderivimplicittocnexp; with one independent gate ODE that substitution is exact.References