Modeling a Ca2+ Channel/BKCa Channel Complex at the Single-Complex Level

BK_Ca-channel activity often affects the firing properties of neurons, the shapes of neuronal action potentials (APs), and in some cases the extent of neurotransmitter release. It has become clear that BK_Ca channels often form complexes with voltage-gated Ca²⁺ channels (CaV channels) such that when a CaV channel is activated, the ensuing influx of Ca²⁺ activates its closely associated BK_Ca channel. Thus, in modeling the electrical properties of neurons, it would be useful to have quantitative models of CaV/BK_Ca complexes. Furthermore, in a population of CaV/BK_Ca complexes, all BK_Ca channels are not exposed to the same Ca²⁺ concentration at the same time. Thus, stochastic rather than deterministic models are required. To date, however, no such models have been described. Here, however, I present a stochastic model of a CaV2.1/BK_Ca(α-only) complex, as might be found in a central nerve terminal. The CaV2.1/BK_Ca model is based on kinetic modeling of its two component channels at physiological temperature. Surprisingly, The CaV2.1/BK_Ca model predicts that although the CaV channel will open nearly every time during a typical cortical AP, its associated BK_Ca channel is expected to open in only 30% of trials, and this percentage is very sensitive to the duration of the AP, the distance between the two channels in the complex, and the presence of fast internal Ca²⁺ buffers. Also, the model predicts that the kinetics of the BK_Ca currents of a population of CaV2.1/BK_Ca complexes will not be limited by the kinetics of the CaV2.1 channel, and during a train of APs, the current response of the complex is expected to faithfully follow even very rapid trains. Aside from providing insight into how these complexes are likely to behave in vivo, the models presented here could also be of use more generally as components of higher-level models of neural function.