The current-voltage relation for electropores with conductivity gradients

In electroporation, an electric field transiently permeabilizes the cell membrane to gain access to the cytoplasm, and to deliver active agents such as DNA, proteins, and drug molecules. Past work suggests that the permeabilization is caused by the formation of aqueous, conducting pores on the lipid membrane, which are also known as electropores. The current-voltage relation across the membrane-bound pores is critical for understanding and predicting electroporation. In this work, we solve the Nernst-Planck equations in a geometry encompassing an isolated electropore to investigate this relation. In particular, we study cases where the intra- and extracellular electrical conductivities differ. We first derive an analytical solution, which is subsequently validated with a direct numerical simulation using a finite volume method. The main result of the current work is a formula for the effective pore resistance as a function of the pore radius, the membrane thickness, and the intra- and extracellular conductivities. This formula can be incorporated into whole-cell or planar-membrane electroporation models for system-level prediction and understanding.