KM_MA2020_GoC

KM_MA2020_GoC#

class braincell.channel.KM_MA2020_GoC(size, g_max=Quantity(0.25, 'mS / cm^2'), temp=Quantity(303.15, 'K'), name=None)#

M-type potassium current of the Golgi cell model.

Slow, non-inactivating M-type potassium current imported from the cerebellar Golgi cell model of Masoli et al. (2020) [2]. A single first-order n gate 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 to 0.25 mS/cm2, which is exactly the source mechanism’s gkbar = 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 the celsius = 30 (degC) written in the source mechanism’s PARAMETER block.

  • name (str | None) – Optional channel name.

See also

KM_MA2020_GrC

Granule-cell port of the same mechanism.

KM_RI2021_SC

Stellate-cell port of the same mechanism.

Kv1p1_MA2020_GoC

Low-threshold Kv1.1 current of the same Golgi cell model.

Notes

Ported from GoC/channel/KM_MA20_GoC.mod. That file, the granule port GrC/channel/KM_MA20_GrC.mod and the stellate port SC/channel/KM_RI21_SC.mod are byte-identical apart from their SUFFIX line, and so are the three BrainCell classes: the rate constants above are shared verbatim with KM_MA2020_GrC and KM_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 .mod source 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.0 at a reference of 22 degrees Celsius, so HH.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 .mod file instead multiplies Q10 into alp_n and bet_n, which divides its tau_n by the same factor. The two forms are algebraically identical, but it means f_n_tau() returns the q10-free time constant rather than the mechanism’s tau_n.

Provenance, and what the header does not establish. The .mod header carries Author: A. Fontana and CoAuthor: T.Nieus. That credit line is copy-pasted verbatim across every cell-type port of this mechanism and names people unrelated to the Golgi-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 Golgi-cell paper [2] names the model this parameterisation was imported from, not the origin of the equations. Reference [2] is the Golgi-cell paper specifically – the companion granule-cell paper of the same year belongs to KM_MA2020_GrC, not to this class.

Conductance default. 0.25 mS/cm2 is the deposit’s tuned value, carried across from the .mod file. It is not a value printed in either cited paper.

Import deviations. The original mechanism’s NMODL TABLE over [-100, 30] mV, covering n_inf and tau_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 from derivimplicit to cnexp; with one independent gate ODE that substitution is exact.

References

root_type#

alias of Potassium