golgi: an open-source graphical platform for image-to-recruitment modeling of peripheral nerve stimulation

Computational models of peripheral nerve stimulation--coupling finite-element bioelectric fields to biophysical axon models--have become essential for designing electrodes and waveforms for neuromodulation therapies. Yet the established open tools are code-first and assume substantial modeling expertise, and several depend on commercial finite-element solvers, placing realistic nerve modeling out of reach for many experimentalists and clinicians. We present golgi, an open-source platform that takes a peripheral nerve from image to stimulated fiber population through a single graphical interface, with an equivalent scriptable Python API and command-line interface for batch studies. golgi integrates the full pipeline: image segmentation (or import of surfaces or masks), automated multi-region tetrahedral meshing, anisotropic finite-element solution of the extracellular field with explicit perineurium contact impedance, generation of realistic fiber populations and their three-dimensional trajectories (straight, or curved via a quasi-static streamline solver), and biophysical activation thresholds through interchangeable NEURON and a GPU-accelerated surrogate backend. We demonstrate golgi on extruded multifascicular swine and human cervical vagus nerves and on real three-dimensional, branching human and rabbit vagus nerves reconstructed from micro-computed tomography. It reproduces the physiological fiber-diameter recruitment order, quantifies fascicular selectivity and current steering with a multi-contact cuff electrode, and resolves anatomically defined nerve branches. Using this branch resolution, we find that selectively engaging a vagal cardiac branch from a proximal cuff depends on anatomy. In the rabbit, whose cardiac fibers are predominantly small and whose superior cardiac branch forms a discrete, spatially segregated tract, current steering isolates even its small cardiac (B-type) fibers; in the human cervical vagus only the large myelinated fibers are separable, because the high thresholds of the small cardiac fibers force stimulus currents that also recruit off-target fibers. Every study can be exported as an integrity-hashed, self-contained bundle whose finite-element-to-recruitment provenance is verifiable byte-for-byte with a single command--a reproducibility guarantee absent from existing tools. By combining non-specialist usability, anatomical realism, and verifiable reproducibility in one open package, golgi lowers the barrier to in-silico peripheral nerve stimulation modeling. golgi is freely available as open-source software. Author summaryElectrical stimulation of peripheral nerves treats a growing range of conditions, from epilepsy to inflammatory and cardiovascular disease. Deciding where to place an electrode and how to shape the stimulus increasingly relies on computer models that combine the electric field around the electrode with detailed models of how individual nerve fibers respond. We found that existing software for this is powerful but primarily designed for expert modelers: it generally requires programming, substantial modeling expertise, and sometimes expensive commercial software, which can limit its adoption by experimentalists and clinicians. We built golgi to remove that barrier. With golgi, a user can go from a nerve image all the way to predicted fiber recruitment through a single point-and-click interface, while advanced users keep full scripting control. golgi builds anatomically realistic nerve models, simulates how different fiber types and fascicles are recruited, and lets users compare electrode designs. Using golgi, we also found that whether a small but clinically important nerve branch--such as the cardiac branch of the vagus nerve--can be selectively stimulated depends on its anatomy. In a rabbit nerve, where this branch forms a discrete, spatially separated bundle and its fibers are mostly small, even its small fibers can be targeted from a cuff on the main trunk; in the human, only the large fibers can be reached selectively, because the small cardiac fibers are harder to excite and the stronger currents needed to reach them also activate off-target fibers. Critically, every result can be packaged so that anyone else can verify it reproduces exactly--making peripheral nerve stimulation models easier to build, share, and trust.