Computational modeling of neurons: intensity-duration relationship of extracellular electrical stimulation for changes in intracellular calcium.

In many instances of extensive nerve damage, the injured nerve never adequately heals, leaving lack of nerve function. Electrical stimulation (ES) has been shown to increase the rate and orient the direction of neurite growth, and is a promising therapy. However, the mechanism in which ES affects neuronal growth is not understood, making it difficult to compare existing ES protocols or to design and optimize new protocols. We hypothesize that ES acts by elevating intracellular calcium concentration ([Ca(2+)]i) via opening voltage-dependent Ca(2+) channels (VDCCs). In this work, we have created a computer model to estimate the ES Ca(2+) relationship. Using COMSOL Multiphysics, we modeled a small dorsal root ganglion (DRG) neuron that includes one Na(+) channel, two K(+) channels, and three VDCCs to estimate [Ca(2+)]i in the soma and growth cone. As expected, the results show that an ES that generates action potentials (APs) can efficiently raise the [Ca(2+)]i of neurons. More interestingly, our simulation results show that sub-AP ES can efficiently raise neuronal [Ca(2+)]i and that specific high-voltage ES can preferentially raise [Ca(2+)]i in the growth cone. The intensities and durations of ES on modeled growth cone calcium rise are consistent with directionality and orientation of growth cones experimentally shown by others. Finally, this model provides a basis to design experimental ES pulse parameters, including duration, intensity, pulse-train frequency, and pulse-train duration to efficiently raise [Ca(2+)]i in neuronal somas or growth cones.