Vagus nerve stimulation (VNS). Image from: http://heatherdane.com/
Computational MIDA head model , specifically realized for neurostimulation investigations.
Vagus nerve stimulation (VNS) was approved in 1997 by the US Food and Drug Administration (FDA)  as an invasive neuromodulator approach for the treatment of epilepsy in anti-epileptic drug (AED) resistant subjects. As the vagus nerve (VN) innervates many organs, it is candidate for many new potentially therapeutically relevant applications of selective neurostimulation.
The VN (like many other large nerves in the human body) is composed by multiple functional units assembling many myelinated A- and B- axons of different diameters combined with small unmyelinated C-fibers. Hereafter, VNS approaches for therapeutically relevant applications require approaches providing high fiber selectivity. Electrode arrays with optimized stimulation waveforms can be used to provide, in theory, such selectivity. Computational models that feature simplified or realistic VN models embedded in realistic human anatomical models along realistic trajectories, together with electrophysiological models of axonal fibers capturing the electric-neuronal interaction, are fundamental for the computationally assisted formulation of new VNS protocols, the design of electrode arrays, the optimization of stimulation waveforms, and the prediction of setup-specific axonal fiber recruitment.
Functionalized VN models, featuring head models, nerve trajectories with an arbitrary numbers of electrophysiological axonal models, can now be created in Sim4Life with the recently added T-NEURO feature in combination with electromagnetic (EM) simulations in the low-frequency range.
Sim4Life, with its expanded simulative abilities of T-NEURO, permits the investigation of the mechanisms of interaction between EM fields and the electrical activity of neuronal membranes. A wide range of applications ranging from neurostimulation investigations (e.g. transcranial electric (TES) and magnetic (TMS) stimulation), safety (e.g. against peripheral nerve stimulation (PNS) in magnetic resonance imaging (MRI)), as well as the computationally assisted development of electroceuticals or neuroprostetic devices, is now possible. Anatomically realistic compartmentalized electrophysiological representations of axons and neurons modeled using NEURON  libraries according to Sim4Life pre-defined biophysical models or available in the web-depository databases – can be positioned within the computational human head/body models to predict the physiological response of nerve and neurons to applied electric (E-) or magnetic (M-) fields. Highly customizable solutions for modeling, execution of simulations, and post-processing analysis for numerical optimization are also provided in Sim4Life.
2D section of Helmers’s VN model  (top); simplified 3D model (middle); anatomically realistic VN model that follows the trajectory of the VN in the MIDA model [2,5] (bottom).
Cross-sectional 2D VN models that feature, e.g., epineurium, perineurium, and fascicles, can be created from medical images or scratched in Sim4Life (see Figure top) with the graphical user interface (GUI). 3D models of the VN can be created, for example, by extruding realistic 2D nerve cross sections along user-defined trajectories (see Figure middle), or along anatomical nerve trajectories within computational human models (see Figure bottom) Electrode geometries can be imported in Sim4life as CAD files or created as parameterized model objects even using available templates entities (helical, spiral, etc.).
Sim4Life’s electro-quasi-static current dominated (EQSCD) solver can be used to set up low frequency EM simulations for either homogeneous or anisotropic dielectric tissue parameters. Electrode voltages may be assigned as Dirichlet boundary conditions.
Sim4Life’s post-processing tools permit the analysis and visualization of exposure-related quantities such as E-field distribution, input currents at electrodes as well as neurostimulation-related quantities as Eigenvalues and Eigenvectors of the E-field Jacobian to identify regions of neurostimulation on the basis of the ‘activation function concept’ .
Fascicles populated with axonal trajectories modeled as lines.
Axonal trajectories, created in terms of splines either using the GUI CAD features or numerically using the Python interface can be positioned in the relevant structures (i.e., fascicles) according to user-defined rules (number of axons, length, displacement, etc.) or following realistic nerve trajectories in computational human models. They can be created as well as random spline trajectories within specified entities (for example within fascicles) using the Sim4Life IMSafe tool for implant safety evaluations.
Complete electrophysiological axonal models are created assigning fiber diameters and pre-defined biophysical models (SENN , Sweeney , MOTOR ) or custom (as hoc files) to the axonal trajectories. Nerves can be functionalized with arbitrary number of axons and arbitrary biophysical descriptions to mime realistic mixture of A-, B-, and C- fibers.
T-Neuro simulations can be executed serially or in parallel. Each simulation can include multiple independent electric sources of neurostimulation (as in the case of electrode arrays), each with its own stimulation waveform.
Point or line sensors can be defined to record transmembrane potential or currents information, to be used for post-pro analyses. The implemented ‘titration procedure’ is fundamental to identify threshold values of E-field or stimulation related quantities (such as input currents) for initiation of action potentials (APs) within each individual axon in the model.
Visualization of the activation function along fiber trajectories (top); examples of fiber recruitment curves for three types of axons populations (bottom).
Sim4Life T-Neuro is conceived to assist in the computational assisted investigation of neurostimulation, the design and optimization of electroceuticals and neurostimulators, and for studies about the mechanisms of electromagnetic-neuronal interaction. In this example the benefits of using T-Neuro to investigate fiber recruitment on realistic or simplified models of Vagus Nerve have been illustrated. Similar procedures can be applied to model arbitrary complex nerve stimulation in arbitrary computational human or animal bodies. Sim4Life provides unique and unprecedented features to investigate the complex domain of neurostimulation, affirming itself as the leader tool in the market for neurostimulation investigations
|Sim4Life → Computable Human Phantoms → ViP3.0|
|Sim4Life → Tissue Models → T-NEURO|
|Sim4Life → Physics Models → P-EM-QS|
|Sim4Life → Modules → IMSAFE|
|Sim4Life → Framework → PYTHON|