The subcellular distribution of T-type Ca2+ channels in interneurons of the lateral geniculate nucleus

Mechanism of an Intrinsic Oscillation in Rat Geniculate Interneurons

Depolarizing current injections produced a rhythmic bursting of action potentials - a bursting oscillation - in a set of local interneurons in the lateral geniculate nucleus (LGN) of rats. The current dynamics underlying this firing pattern have not been determined, though this cell type constitutes an important cellular component of thalamocortical circuitry, and contributes to both pathologic and non-pathologic brain states. We thus investigated the source of the bursting oscillation using pharmacological manipulations in LGN slices in vitro and in silico. 1. Selective blockade of calcium channel subtypes revealed that high-threshold calcium currentsI L andI P contributed strongly to the oscillation. 2. Increased extracellular K+ concentration (decreased K+currents) eliminated the oscillation. 3. Selective blockade of K+ channel subtypes demonstrated that the calcium-sensitive potassium current (I A H P ) was of primary importance. A morphologically simplified, multicompartment model of the thalamic interneuron characterized the oscillation as follows: 1. The low-threshold calcium currentI T provided the strong initial burst characteristic of the oscillation. 2. Alternating fluxes through high-threshold calcium channels andI A H P then provided the continuing oscillation's burst and interburst periods respectively. This interplay betweenI L andI A H P contrasts with the current dynamics underlying oscillations in thalamocortical and reticularis neurons, which primarily involveI T andI H , orI T andI A H P respectively. These findings thus point to a novel electrophysiological mechanism for generating intrinsic oscillations in a major thalamic cell type. Because local interneurons can sculpt the behavior of thalamocortical circuits, these results suggest new targets for the manipulation of ascending thalamocortical network activity.

Uncertainpy: A Python Toolbox for Uncertainty Quantification and Sensitivity Analysis in Computational Neuroscience

Computational models in neuroscience typically contain many parameters that are poorly constrained by experimental data. Uncertainty quantification and sensitivity analysis provide rigorous procedures to quantify how the model output depends on this parameter uncertainty. Unfortunately, the application of such methods is not yet standard within the field of neuroscience. Here we present Uncertainpy, an open-source Python toolbox, tailored to perform uncertainty quantification and sensitivity analysis of neuroscience models. Uncertainpy aims to make it quick and easy to get started with uncertainty analysis, without any need for detailed prior knowledge. The toolbox allows uncertainty quantification and sensitivity analysis to be performed on already existing models without needing to modify the model equations or model implementation. Uncertainpy bases its analysis on polynomial chaos expansions, which are more efficient than the more standard Monte-Carlo based approaches. Uncertainpy is tailored for neuroscience applications by its built-in capability for calculating characteristic features in the model output. The toolbox does not merely perform a point-to-point comparison of the "raw" model output (e.g., membrane voltage traces), but can also calculate the uncertainty and sensitivity of salient model response features such as spike timing, action potential width, average interspike interval, and other features relevant for various neural and neural network models. Uncertainpy comes with several common models and features built in, and including custom models and new features is easy. The aim of the current paper is to present Uncertainpy to the neuroscience community in a user-oriented manner. To demonstrate its broad applicability, we perform an uncertainty quantification and sensitivity analysis of three case studies relevant for neuroscience: the original Hodgkin-Huxley point-neuron model for action potential generation, a multi-compartmental model of a thalamic interneuron implemented in the NEURON simulator, and a sparsely connected recurrent network model implemented in the NEST simulator.

Biophysical network modeling of the dLGN circuit: Effects of cortical feedback on spatial response properties of relay cells

Despite half-a-century of research since the seminal work of Hubel and Wiesel, the role of the dorsal lateral geniculate nucleus (dLGN) in shaping the visual signals is not properly understood. Placed on route from retina to primary visual cortex in the early visual pathway, a striking feature of the dLGN circuit is that both the relay cells (RCs) and interneurons (INs) not only receive feedforward input from retinal ganglion cells, but also a prominent feedback from cells in layer 6 of visual cortex. This feedback has been proposed to affect synchronicity and other temporal properties of the RC firing. It has also been seen to affect spatial properties such as the center-surround antagonism of thalamic receptive fields, i.e., the suppression of the response to very large stimuli compared to smaller, more optimal stimuli. Here we explore the spatial effects of cortical feedback on the RC response by means of a a comprehensive network model with biophysically detailed, single-compartment and multicompartment neuron models of RCs, INs and a population of orientation-selective layer 6 simple cells, consisting of pyramidal cells (PY). We have considered two different arrangements of synaptic feedback from the ON and OFF zones in the visual cortex to the dLGN: phase-reversed (‘push-pull’) and phase-matched (‘push-push’), as well as different spatial extents of the corticothalamic projection pattern. Our simulation results support that a phase-reversed arrangement provides a more effective way for cortical feedback to provide the increased center-surround antagonism seen in experiments both for flashing spots and, even more prominently, for patch gratings. This implies that ON-center RCs receive direct excitation from OFF-dominated cortical cells and indirect inhibitory feedback from ON-dominated cortical cells. The increased center-surround antagonism in the model is accompanied by spatial focusing, i.e., the maximum RC response occurs for smaller stimuli when feedback is present.

The functional role of the dorsal lateral geniculate nucleus (dLGN), placed on route from retina to primary visual cortex in the early visual pathway, is still poorly understood. A striking feature of the dLGN circuit is that dLGN cells not only receive feedforward input from the retina, but also a prominent feedback from cells in the visual cortex. It has been seen in experiments that cortical feedback modifies the spatial properties of dLGN cells in response to visual stimuli. In particular, it has been shown to increase the center-surround antagonism for flashing-spot and patch-grating visual stimuli, i.e., the suppression of responses to very large stimuli compared to smaller stimuli. Here we investigate the putative mechanisms behind this feature by means of a comprehensive network model of biophysically detailed neuron models for RCs and INs in the dLGN and orientation-selective cortical cells providing the feedback. Our results support that the experimentally observed feedback effects may be due to a phase-reversed (‘push-pull’) arrangement of the cortical feedback where ON-symmetry RCs receive (indirect) inhibitory feedback from ON-dominated cortical cell and excitation from OFF-dominated cortical cells, and vice versa for OFF-symmetry RCs.

Effect of Ionic Diffusion on Extracellular Potentials in Neural Tissue

Recorded potentials in the extracellular space (ECS) of the brain is a standard measure of population activity in neural tissue. Computational models that simulate the relationship between the ECS potential and its underlying neurophysiological processes are commonly used in the interpretation of such measurements. Standard methods, such as volume-conductor theory and current-source density theory, assume that diffusion has a negligible effect on the ECS potential, at least in the range of frequencies picked up by most recording systems. This assumption remains to be verified. We here present a hybrid simulation framework that accounts for diffusive effects on the ECS potential. The framework uses (1) the NEURON simulator to compute the activity and ionic output currents from multicompartmental neuron models, and (2) the electrodiffusive Kirchhoff-Nernst-Planck framework to simulate the resulting dynamics of the potential and ion concentrations in the ECS, accounting for the effect of electrical migration as well as diffusion. Using this framework, we explore the effect that ECS diffusion has on the electrical potential surrounding a small population of 10 pyramidal neurons. The neural model was tuned so that simulations over ∼100 seconds of biological time led to shifts in ECS concentrations by a few millimolars, similar to what has been seen in experiments. By comparing simulations where ECS diffusion was absent with simulations where ECS diffusion was included, we made the following key findings: (i) ECS diffusion shifted the local potential by up to ∼0.2 mV. (ii) The power spectral density (PSD) of the diffusion-evoked potential shifts followed a 1/f2 power law. (iii) Diffusion effects dominated the PSD of the ECS potential for frequencies up to several hertz. In scenarios with large, but physiologically realistic ECS concentration gradients, diffusion was thus found to affect the ECS potential well within the frequency range picked up in experimental recordings.

Biophysical Network Modelling of the dLGN Circuit: Different Effects of Triadic and Axonal Inhibition on Visual Responses of Relay Cells

Despite its prominent placement between the retina and primary visual cortex in the early visual pathway, the role of the dorsal lateral geniculate nucleus (dLGN) in molding and regulating the visual signals entering the brain is still poorly understood. A striking feature of the dLGN circuit is that relay cells (RCs) and interneurons (INs) form so-called triadic synapses, where an IN dendritic terminal can be simultaneously postsynaptic to a retinal ganglion cell (GC) input and presynaptic to an RC dendrite, allowing for so-called triadic inhibition. Taking advantage of a recently developed biophysically detailed multicompartmental model for an IN, we here investigate putative effects of these different inhibitory actions of INs, i.e., triadic inhibition and standard axonal inhibition, on the response properties of RCs. We compute and investigate so-called area-response curves, that is, trial-averaged visual spike responses vs. spot size, for circular flashing spots in a network of RCs and INs. The model parameters are grossly tuned to give results in qualitative accordance with previous in vivo data of responses to such stimuli for cat GCs and RCs. We particularly investigate how the model ingredients affect salient response properties such as the receptive-field center size of RCs and INs, maximal responses and center-surround antagonisms. For example, while triadic inhibition not involving firing of IN action potentials was found to provide only a non-linear gain control of the conversion of input spikes to output spikes by RCs, axonal inhibition was in contrast found to substantially affect the receptive-field center size: the larger the inhibition, the more the RC center size shrinks compared to the GC providing the feedforward excitation. Thus, a possible role of the different inhibitory actions from INs to RCs in the dLGN circuit is to provide separate mechanisms for overall gain control (direct triadic inhibition) and regulation of spatial resolution (axonal inhibition) of visual signals sent to cortex.

Anatomical localization of Cav3.1 calcium channels and electrophysiological effects of T-type calcium channel blockade in the motor thalamus of MPTP-treated monkeys

Conventional anti-Parkinsonian dopamine replacement therapy is often complicated by side effects that limit the use of these medications. There is a continuing need to develop nondopaminergic approaches to treat Parkinsonism. One such approach is to use medications that normalize dopamine depletion-related firing abnormalities in the basal ganglia-thalamocortical circuitry. In this study, we assessed the potential of a specific T-type calcium channel blocker (ML218) to eliminate pathologic burst patterns of firing in the basal ganglia-receiving territory of the motor thalamus in Parkinsonian monkeys. We also carried out an anatomical study, demonstrating that the immunoreactivity for T-type calcium channels is strongly expressed in the motor thalamus in normal and 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP)-treated monkeys. At the electron microscopic level, dendrites accounted for >90% of all tissue elements that were immunoreactive for voltage-gated calcium channel, type 3.2-containing T-type calcium channels in normal and Parkinsonian monkeys. Subsequent in vivo electrophysiologic studies in awake MPTP-treated Parkinsonian monkeys demonstrated that intrathalamic microinjections of ML218 (0.5 μl of a 2.5-mM solution, injected at 0.1-0.2 μl/min) partially normalized the thalamic activity by reducing the proportion of rebound bursts and increasing the proportion of spikes in non-rebound bursts. The drug also attenuated oscillatory activity in the 3-13-Hz frequency range and increased gamma frequency oscillations. However, ML218 did not normalize Parkinsonism-related changes in firing rates and oscillatory activity in the beta frequency range. Whereas the described changes are promising, a more complete assessment of the cellular and behavioral effects of ML218 (or similar drugs) is needed for a full appraisal of their anti-Parkinsonian potential.