Imperial College London


Faculty of EngineeringDepartment of Bioengineering

Professor of Neurotechnology



+44 (0)20 7594 1533s.schultz Website




4.11Royal School of MinesSouth Kensington Campus






BibTex format

author = {Lubba, CH and Guen, YL and Jarvis, S and Jones, NS and Cork, SC and Eftekhar, A and Schultz, SR},
doi = {10.1101/196451},
publisher = {Cold Spring Harbor Laboratory},
title = {Multiscale simulation of peripheral neural signaling},
url = {},
year = {2017}

RIS format (EndNote, RefMan)

AB - <jats:title>Abstract</jats:title><jats:p>Bioelectronic Medicines that modulate the activity patterns on peripheral nerves have promise as a new way of treating diverse medical conditions from epilepsy to rheumatism. Progress in the field builds upon time consuming and expensive experiments in living organisms to evaluate spontaneous activity patterns, stimulation efficiency, and organ responses. To reduce experimentation load and allow for a faster, more detailed analysis of both recording from and stimulation of peripheral nerves, adaptable computational models incorporating insights won in experiments will be of great help. We present a peripheral nerve simulator that combines biophysical axon models and numerically solved and idealized extracellular space models in one environment. Two different scales of abstraction were merged. On the one hand we modeled the extracellular space in a finite element solver as a three dimensional resistive continuum governed by the electro-quasistatic approximation of the Maxwell equations. Potential distributions were precomputed for different media (homogeneous, nerve in saline, nerve in cuff). Axons, on the other hand, were modeled at a higher level of abstraction as one dimensional chains of compartments; each consisting of lumped linear elements and, for some, channels with non-linear dynamics. Unmyelinated fibres were based on the Hodgkin-Huxley model; for myelinated fibers, we instead adapted the model proposed by McIntyre et al. in 2002 to smaller diameters. To obtain realistic axon shapes, an iterative algorithm positioned fibers along the nerve with variable tortuosity, with tortuosity values fit to imaged trajectories. We validated our model with data from the stimulated rat vagus nerve. Simulation results predicted that tortuosity leads to differentiation in recorded signal shapes, with unmyelinated axons being the most affected. Tortuosity was further shown to increase the stimulation threshold. The
AU - Lubba,CH
AU - Guen,YL
AU - Jarvis,S
AU - Jones,NS
AU - Cork,SC
AU - Eftekhar,A
AU - Schultz,SR
DO - 10.1101/196451
PB - Cold Spring Harbor Laboratory
PY - 2017///
TI - Multiscale simulation of peripheral neural signaling
UR -
ER -