Predictions of energy efficient Berger-Levy model neurons with constraints

More specifically, in our contribution, in line with [1] and [2], a neuron is modeled as a memory-less constant communication channel with a Gamma conditional probability distribution function (PDF) [1]. In this setting, the channel input and output are the excitatory postsynaptic potential intensity, $\lambda$, and the inter spike interval (ISI), $t$, with PDFs $f_{\Lambda }\left(\lambda \right)$ and $f_T\left(t\right)$, respectively. We then formulate two new constraints: First, we impose a lower bound $t_{\mathsf{\text{min}}}$ on the duration $t$of ISIs. The rational for this is to account for a maximal firing rate. Second, we consider a peak energy expenditure constraint per ISI as compared to only bounding the expected energy expenditure. This translates into an upper bound $t_{\mathsf{\text{max}}}$ on the ISI duration. We then derive the $f_T\left(t\right)$ (corresponding to valid $f_{\Lambda }\left(\lambda \right)$) of a bpj-maximizing neuron for the original unconstrained setting from [1] and in the presence of the above two constraints for different expected ISIs. (Details omitted here for brevity.) Figure 1 shows three $f_T\left(t\right)$s obtained in the unconstrained (dashed curves) and constrained settings (solid curves) for $t_{\mathsf{\text{min}}}=1$and $t_{\mathsf{\text{max}}}=5$. While the constrained and unconstrained solutions have the same mean, the shape of their $f_T\left(t\right)$ differ. For comparison with experimental data, we computed the coefficient of variation (CV) as a function of the mean ISI as an "observable" (Figure 2), which is easier to measure experimentally than the full distribution $f_T\left(t\right)$. Interestingly, the CV is predicted i) to be lower in the constrained setting, and ii) to increase and then decrease with the mean ISI while it only decreases in the unconstrained setting. Thus, we demonstrated that constraints can affect predictions based on bpj-maximization, and should be explicitly taken into account. Ongoing work makes these predictions more quantitative via simulating biophysically realistic model neurons.