Background: Deep brain stimulation is an established therapy for several neurological disorders; however, its effects on neuronal activity vary across brain regions and depend on stimulation settings. Understanding these variable responses can aid in the development of physiologically-informed stimulation paradigms in existing or prospective indications. Objective: Provide experimental and computational insights into the brain-region-specific and frequency-dependent effects of extracellular stimulation on neuronal activity. Methods: In patients with movement disorders, single-neuron recordings were acquired from the subthalamic nucleus, substantia nigra pars reticulata, ventral intermediate nucleus, or reticular thalamus during microstimulation across various frequencies (1e100 Hz) to assess single-pulse and frequencyresponse functions. Moreover, a biophysically-realistic computational framework was developed which generated postsynaptic responses under the assumption that electrical stimuli simultaneously activated all convergent presynaptic inputs to stimulation target neurons. The framework took into consideration the relative distributions of excitatory/inhibitory afferent inputs to model site-specific responses, which were in turn embedded within a model of short-term synaptic plasticity to account for stimulation frequency-dependence. Results: We demonstrated microstimulation-evoked excitatory neuronal responses in thalamic structures (which have predominantly excitatory inputs) and inhibitory responses in basal ganglia structures (predominantly inhibitory inputs); however, higher stimulation frequencies led to a loss of site-specificity and convergence towards neuronal suppression. The model confirmed that site-specific responses could be simulated by accounting for local neuroanatomical/microcircuit properties, while suppression of neuronal activity during high-frequency stimulation was mediated by short-term synaptic depression. Conclusions: Brain-region-specific and frequency-dependant neuronal responses could be simulated by considering neuroanatomical (local microcircuitry) and neurophysiological (short-term plasticity) properties.

Multiplexing refers to the simultaneous encoding of two or more signals. Neurons have been shown to multiplex, but different stimuli require different multiplexing strategies. Whereas the frequency and amplitude of periodic stimuli can be encoded by the timing and rate of the same spikes, natural scenes, which comprise areas over which intensity varies gradually and sparse edges where intensity changes abruptly, require a different multiplexing strategy. Recording in vivo from neurons in primary somatosensory cortex during tactile stimulation, we found that stimulus onset and offset (edges) evoked highly synchronized spiking, whereas other spikes in the same neurons occurred asynchronously. Stimulus intensity modulated the rate of asynchronous spiking, but did not affect the timing of synchronous spikes. From this, we hypothesized that spikes driven by high- and low-contrast stimulus features can be distinguished on the basis of their synchronization, and that differentially synchronized spiking can thus be used to form multiplexed representations. Applying a Bayesian decoding method, we verified that information about high- and low-contrast features can be recovered from an ensemble of model neurons receiving common input. Equally good decoding was achieved by distinguishing synchronous from asynchronous spikes and applying reverse correlation methods separately to each spike type. This result, which we verified with patch clamp recordings in vitro, demonstrates that neurons receiving common input can use the rate of asynchronous spiking to encode the intensity of low-contrast features while using the timing of synchronous spikes to encode the occurrence of high-contrast features. We refer to this strategy as synchrony-division multiplexing.