ROOTS: An Algorithm to Generate Biologically Realistic Cortical Axons and an Application to Electroceutical Modeling

Evolving characterization of the human hyperdirect pathway

The hyperdirect pathway (HDP) represents the main glutamatergic input to the subthalamic nucleus (STN), through which the motor and prefrontal cerebral cortex can modulate basal ganglia activity. Further, direct activation of the motor HDP is thought to be an important component of therapeutic deep brain stimulation (DBS), mediating the disruption of pathological oscillations. Alternatively, unintended recruitment of the prefrontal HDP may partly explain some cognitive side effects of DBS therapy. Previous work describing the HDP has focused on non-human primate (NHP) histological pathway tracings, diffusion-weighted MRI analysis of human white matter, and electrophysiology studies involving paired cortical recordings with DBS. However, none of these approaches alone yields a complete understanding of the complexities of the HDP. As such, we propose that generative modeling methods hold promise to bridge anatomy and physiology results, from both NHPs and humans, into a more detailed representation of the human HDP. Nonetheless, numerous features of the HDP remain to be experimentally described before model-based methods can simulate corticosubthalamic activity with a high degree of scientific detail. Therefore, the goals of this review are to examine the experimental evidence for HDP projections from across the primate neocortex and discuss new data which are required to improve the utility of anatomical and biophysical models of the human corticosubthalamic system.

Subthalamic deep brain stimulation of an anatomically detailed model of the human hyperdirect pathway

The motor hyperdirect pathway (HDP) is considered a key target in the treatment of Parkinson's disease with subthalamic deep brain stimulation (DBS). This hypothesis is partially derived from the association of HDP activation with evoked potentials (EPs) generated in the motor cortex and subthalamic nucleus (STN) after a DBS pulse. However, the biophysical details of how and when DBS-induced action potentials (APs) in HDP neurons reach their terminations in the cortex or STN remain unclear. Therefore, we used an anatomically detailed representation of the motor HDP, as well as the internal capsule (IC), in a model of human subthalamic DBS to explore AP activation and transmission in the HDP and IC. Our results show that small diameter HDP axons exhibited AP initiation in their subthalamic terminal arbor, which resulted in relatively long transmission latencies to cortex (∼3.5-8 ms). Alternatively, large diameter HDP axons were most likely to be directly activated in the capsular region, which resulted in short transmission times to the cortex (∼1-3 ms). However, those large diameter HDP antidromic APs would be indistinguishable from any other IC axons that were also activated by the stimulus. Conversely, DBS-induced APs in both small and large diameter HDP axons reached their synaptic boutons in the STN with similar timings, but both spanned a wide temporal range (∼0.5-5 ms). We also found that using anodic or bipolar stimulation helped to bias activation of the HDP over the IC. These computational results provide useful information for linking HDP activation with EP recordings in clinical experiments.NEW & NOTEWORTHY We used biophysical models to study pathway recruitment and conduction latencies of the hyperdirect pathway (HDP) in response to subthalamic deep brain stimulation (DBS). The model system allowed us to assess the influence of increased anatomical realism on pathway activity and the possibility of identifying HDP activity in evoked potentials (EPs) recorded in either the subthalamic nucleus (STN) or cortex. The model predicts that HDP activation is accentuated by complex axonal branching in the STN.

Histology-driven model of the macaque motor hyperdirect pathway

Emerging appreciation for the hyperdirect pathway (HDP) as an important cortical glutamatergic input to the subthalamic nucleus (STN) has motivated a wide range of recent investigations on its role in motor control, as well as the mechanisms of subthalamic deep brain stimulation (DBS). However, the pathway anatomy and terminal arbor morphometry by which the HDP links cortical and subthalamic activity are incompletely understood. One critical hindrance to advancing understanding is the lack of anatomically detailed population models which can help explain how HDP pathway anatomy and neuronal biophysics give rise to spatiotemporal patterns of stimulus-response activity observed in vivo. Therefore, the goal of this study was to establish a population model of motor HDP axons through application of generative algorithms constrained by recent histology and imaging data. The products of this effort include a de novo macaque brain atlas, detailed statistical analysis of histological reconstructions of macaque motor HDP axons, and the generation of 10,000 morphometrically constrained synthetic motor HDP axons. The synthetic HDP axons exhibited a 3.8% mean error with respect to parametric distributions of the fiber target volume, total length, number of bifurcations, bifurcation angles, meander angles, and segment lengths measured in BDA-labeled HDP axon reconstructions. As such, this large population of synthetic motor HDP axons represents an anatomically based foundation for biophysical simulations that can be coupled to electrophysiological and/or behavioral measurements, with the goal of better understanding the role of the HDP in motor system activity.

Deep brain stimulation of terminating axons

BACKGROUND

Deep brain stimulation (DBS) of the subthalamic region is an established treatment for the motor symptoms of Parkinson's disease. Several types of neural elements reside in the subthalamic region, including subthalamic nucleus (STN) neurons, fibers of passage, and terminating afferents. Recent studies suggest that direct activation of a specific population of subthalamic afferents, known as the hyperdirect pathway, may be responsible for some of the therapeutic effects of subthalamic DBS.

OBJECTIVE

The goal of this study was to quantify how axon termination affects neural excitability from DBS. We evaluated how adjusting different stimulation parameters influenced the relative excitability of terminating axons (TAs) compared to fibers of passage (FOPs).

METHODS

We used finite element electric field models of DBS, coupled to multi-compartment cable models of axons, to calculate activation thresholds for populations of TAs and FOPs. These generalized models were used to evaluate the response to anodic vs. cathodic stimulation, with short vs. long stimulus pulses.

RESULTS

Terminating axons generally exhibited lower thresholds than fibers of passage across all tested parameters. Short pulse widths accentuated the relative excitability of TAs over FOPs.

CONCLUSION(S)

Our computational results demonstrate a hyperexcitability of terminating axons to DBS that is robust to variation in the stimulation parameters, as well as the axon model parameters.