Electric discharge of electrocytes: Modelling, analysis and simulation

•A PNP (or PDE) model is formulated for the mechanism of electrocyte discharge.•Time variation of the electric fields is essential for the conservation of total current.•PNP model reduces to an ordinary differential equations model at leading order by asymptotic analysis.•Delay in contact leads to interesting dynamics and new features, compared to the non-contact case investigated previously. In this paper, we investigate the electric discharge of electrocytes by extending our previous work on the generation of electric potential. We first give a complete formulation of a single cell unit consisting of an electrocyte and a resistor, based on a Poisson-Nernst-Planck (PNP) system with various membrane currents as interfacial conditions for the electrocyte and a Maxwell’s model for the resistor. Our previous work can be treated as a special case with an infinite resistor (or open circuit). Using asymptotic analysis, we simplify our PNP system and reduce it to an ordinary differential equation (ODE) based model. Unlike the case of an infinite resistor, our numerical simulations of the new model reveal several distinct features. A finite current is generated, which leads to non-constant electric potentials in the bulk of intracellular and extracellular regions. Furthermore, the current induces an additional action potential (AP) at the non-innervated membrane, contrary to the case of an open circuit where an AP is generated only at the innervated membrane. The voltage drop inside the electrocyte is caused by an internal resistance due to mobile ions. We show that our single cell model can be used as the basis for a system with stacked electrocytes and the total current during the discharge of an electric eel can be estimated by using our model.