The model, developed by Chun-Hui Huang, Wen-Shuo Chen, and colleagues at National Taiwan University and National Taiwan Ocean University, is built on the classic Hodgkin-Huxley formalism, the biophysical framework that describes how voltage-gated ion channels generate action potentials. Rather than lumping brain regions into abstract nodes, the testbed resolves eleven distinct neuronal populations spanning cortex, thalamus, striatum, and the basal ganglia: cortical pyramidal cells, two classes of cortical interneurons, thalamocortical relay cells, reticular thalamic neurons, striatal fast-spiking interneurons, direct- and indirect-pathway medium spiny neurons, neurons of the external and internal globus pallidus, and subthalamic nucleus cells. Each population carries its own complement of sodium, potassium, calcium, and adaptation currents, allowing the simulation to reproduce complex discharge patterns such as the post-inhibitory rebound bursting that emerges from T-type calcium channels in the subthalamic nucleus and thalamus.

To emulate Parkinson’s disease, the researchers did not simply dial down a single parameter. Instead, they reconfigured the network to mimic the circuit-level consequences of dopamine depletion: cortical drive to the direct pathway was weakened, input to the indirect pathway was strengthened, and feedforward inhibition within the striatum was reduced. Background currents to the globus pallidus and subthalamic nucleus were also adjusted, with the subthalamic baseline dropping from 33 to 23 microamps. The calibrated parkinsonian state showed hallmark reorganization: beta-band power, oscillations in the 13 to 30 hertz range associated with rigidity and slowness of movement, rose nearly fourfold in the internal globus pallidus and fourfold in indirect-pathway striatal neurons, while firing rates of the two medium spiny neuron populations shifted dramatically upward.

Against this pathological backdrop, the team ran two interventions on identical networks. The first was conventional: 130-hertz deep brain stimulation of the subthalamic nucleus, delivered as symmetric, charge-balanced biphasic pulses with a pulse width of 0.06 milliseconds per phase. The second was an idealized periodic visual input, generated by a prescribed input unit that delivers rhythmic pulse trains directly to cortical pyramidal cells and thalamic relay cells. The authors are careful to note that this input unit is not a biophysical model of the retina or the lateral geniculate nucleus; it bypasses the sensory filtering that real visual signals undergo, and the archived identifier suggesting a 40-hertz condition is a legacy label—the implemented drive is a 10-hertz base-pulse train, with the 38 to 42 hertz band analyzed purely as an output measure.

The contrast between the two interventions was stark. Deep brain stimulation crushed the model’s internal measure of subthalamic beta-band synaptic current, reducing it from 0.0184 to 0.0018 arbitrary units—a suppression of roughly ninety percent, with a Cohen’s d effect size of 3.9 and a Holm-corrected p-value of 1.1 times ten to the minus fifth power. The periodic visual input did the opposite on this measure: subthalamic beta power actually rose by more than double. But the visual input produced something deep brain stimulation did not. Thalamic relay cell firing jumped from 8.6 to 35.3 hertz, previously silent reticular thalamic neurons woke up to 6.56 hertz, and cortical pyramidal firing nearly doubled from 4.86 to 10.33 hertz, accompanied by significant increases in broadband gamma power and narrowband power in the 38 to 42 hertz window. Deep brain stimulation, by contrast, left cortical firing statistically unchanged.

One of the most intriguing findings concerns what happened to the thalamus under deep brain stimulation. Spontaneous firing of thalamocortical relay cells fell to zero at every tested stimulation amplitude from 1 to 5 milliamps per square centimeter—a result that could, at first glance, suggest the stimulation had silenced the thalamus entirely, acting like a reversible lesion. The team dug deeper. In a full-resolution diagnostic run, they inspected 2.5 million membrane voltage samples from all ten modeled relay cells and found every sample finite and subthreshold, ranging from about -71 to -58 millivolts, with no crossings through zero millivolts. That pattern is inconsistent with sustained depolarization block, the pathological state in which neurons are stuck in a hyperpolarized inactivity. Instead, the data pointed to subthreshold quiescence: the relay cells remained alive and responsive, simply not firing spontaneously under the calibrated balance of excitatory and inhibitory inputs.

To test whether this quiescence depended on one particular parameter, the researchers ran a matched-seed sensitivity analysis on the inhibitory coupling from the internal globus pallidus to the thalamus, varying it by twenty percent in either direction. Across all thirty runs, thalamic spontaneous firing remained at zero, and all 750,000 saved voltage samples stayed subthreshold, between -72 and -55 millivolts. A further three-parameter sensitivity analysis—varying a GABA decay constant, the subthalamic-pallidal coupling conductance, and the subthalamic injected current, each by twenty percent—left all four principal outcomes intact, with beta suppression under deep brain stimulation remaining between 88.9 and 95.8 percent. These checks do not prove the model is robust to every conceivable perturbation, but they show the headline results are not artifacts of a single finely tuned value.

The dose-response story added another layer of nuance. Rather than a clean monotonic relationship between stimulation intensity and beta suppression, the amplitude series revealed a non-monotonic response: at 1 milliamp per square centimeter, subthalamic beta actually increased above the unstimulated parkinsonian baseline, while at 2 through 5 milliamps the mean beta values settled into a narrow band corresponding to over ninety percent suppression. Because no amplitudes between zero and one milliamp were sampled, the exact transition point where thalamic firing collapses remains unresolved. The authors also stress that the model’s current density is not directly translatable to clinical device current, which depends on electrode geometry, volume conduction, and axonal recruitment in real tissue.

The study is equally candid about its boundaries. The subthalamic beta measure is derived from outgoing synaptic current rather than a simulated local field potential, and the model omits the hyper-direct cortico-subthalamic pathway now considered a key driver of beta oscillations in human Parkinson’s disease—perhaps explaining why the canonical parkinsonian rise in subthalamic beta appeared in the pallidum and striatum instead. The simplified non-REM sleep representation failed to reproduce the pathological beta that clinical recordings show can persist during sleep, and the cortical module, with its simplified horizontal inhibition, produced stronger high-beta synchrony than would be expected clinically. The authors frame these as defined targets for future iterations rather than fatal flaws.

What emerges is not a treatment recommendation but a set of sharply defined, testable hypotheses. Deep brain stimulation and rhythmic visual input, in this calibrated network, produce fundamentally different response profiles: one suppresses pathological subthalamic beta while leaving cortical output largely untouched, the other recruits thalamic and cortical activity vigorously without touching the beta measure. Whether visual recruitment preserves information processing, whether thalamic quiescence under stimulation is compatible with faithful relay, and whether either profile translates to motor benefit in patients are questions the model cannot answer. But by placing both interventions on the same digital substrate—complete with wakefulness, sleep, and anesthesia states—the testbed gives experimentalists a concrete map of where to look next, and a reminder that in the brain, suppressing one circuit’s noise and amplifying another’s signal are not the same thing.

Subject of Research: Computational modeling of deep brain stimulation versus periodic visual input in the Parkinsonian cortico-basal ganglia-thalamic circuit

Article Title: An in silico cortico-basal ganglia-thalamic testbed for comparing visual and deep brain stimulation in Parkinson’s disease

Article References: Huang, C.-H., Chen, W.-S., Lai, J.-R., Cho, C.-H., Sung, Y.-T., Chang, C.-T., & Lin, C.-W. (2026). An in silico cortico-basal ganglia-thalamic testbed for comparing visual and deep brain stimulation in Parkinson’s disease. Frontiers in Computational Neuroscience, 20, Article 1950516. https://doi.org/10.3389/fncom.2026.1950516

Image Credits: AI Generated

DOI: 10.3389/fncom.2026.1950516

Keywords: Parkinson's disease, deep brain stimulation, computational neuroscience, Hodgkin-Huxley model, basal ganglia, thalamus, beta oscillations, visual stimulation, neuromodulation, cortico-basal ganglia-thalamic loop, in silico modeling, gamma oscillations