Reverse-engineered models reveal differential membrane properties of autonomic and cutaneous unmyelinated fibers

Unmyelinated C-fibers constitute the vast majority of axons in peripheral nerves and play key roles in homeostasis and signaling pain. However, little is known about their ion channel expression, which controls their firing properties. Also, because of their small diameters (~ 1 μm), it has not been possible to characterize their membrane properties using voltage clamp. We developed a novel library of isoform-specific ion channel models to serve as the basis functions of our C-fiber models. We then developed a particle swarm optimization (PSO) framework that used the isoform-specific ion channel models to reverse engineer C-fiber membrane properties from measured autonomic and cutaneous C-fiber conduction responses. Our C-fiber models reproduced experimental conduction velocity, chronaxie, action potential duration, intracellular threshold, and paired pulse recovery cycle. The models also matched experimental activity-dependent slowing, a property not included in model optimization. We found that simple conduction responses, characterizing the action potential, were controlled by similar membrane properties in both the autonomic and cutaneous C-fiber models, but complicated conduction response, characterizing the afterpotenials, were controlled by differential membrane properties. The unmyelinated C-fiber models constitute important tools to study autonomic signaling, assess the mechanisms of pain, and design bioelectronic devices. Additionally, the novel reverse engineering approach can be applied to generate models of other neurons where voltage clamp data are not available.