Electrical stimulation in the human cochlea: a computational study based on high-resolution micro-CT scans

Background: Many detailed features of the cochlear anatomy have not been included in existing 3D cochlear models, including the microstructures inside the modiolar bone, which in turn determines the path of auditory nerve fibers (ANFs).

Method: We captured the intricate modiolar microstructures in a 3D human cochlea model reconstructed from μCT scans. A new algorithm was developed to reconstruct ANFs running through the microstructures within the model. Using the finite element method, we calculated the electrical potential as well as its first and second spatial derivatives along each ANF elicited by the cochlear implant electrodes. Simulation results of electrical potential was validated against intracochlear potential measurements. Comparison was then made with a simplified model without the microstructures within the cochlea.

Results: When the stimulus was delivered from an electrode located deeper in the apex, the extent of the auditory nerve influenced by a higher electric potential grew larger; at the same time, the maximal potential value at the auditory nerve also became larger. The electric potential decayed at a faster rate toward the base of the cochlea than toward the apex. Compared to the cochlear model incorporating the modiolar microstructures, the simplified version resulted in relatively small differences in electric potential. However, in terms of the first and second derivatives of electric potential along the fibers, which are relevant for the initiation of action potentials, the two models exhibited large differences: maxima in both derivatives with the detailed model were larger by a factor of 1.5 (first derivative) and 2 (second derivative) in the exemplary fibers. More importantly, these maxima occurred at different locations, and opposite signs were found for the values of second derivatives between the two models at parts along the fibers. Hence, while one model predicts depolarization and spike initiation at a given location, the other may instead predict a hyperpolarization.

Conclusions: Although a cochlear model with fewer details seems sufficient for analysing the current spread in the cochlear ducts, a detailed-segmented cochlear model is required for the reconstruction of ANF trajectories through the modiolus, as well as the prediction of firing thresholds and spike initiation sites.