Real-time adaptive microstimulation increases reliability of electrically evoked cortical potentials

Activity-dependent recruitment of inhibition and excitation in the awake mammalian cortex during electrical stimulation

Electrical stimulation is an effective tool for mapping and altering brain connectivity, with applications ranging from treating pharmacology-resistant neurological disorders to providing sensory feedback for neural prostheses. Paramount to the success of these applications is the ability to manipulate electrical currents to precisely control evoked neural activity patterns. However, little is known about stimulation-evoked responses in inhibitory neurons nor how stimulation-evoked activity patterns depend on ongoing neural activity. In this study, we used 2-photon imaging and cell-type specific labeling to measure single-cell responses of excitatory and inhibitory neurons to electrical stimuli in the visual cortex of awake mice. Our data revealed strong interactions between electrical stimulation and pre-stimulus activity of single neurons in awake animals and distinct recruitment and response patterns for excitatory and inhibitory neurons. This work demonstrates the importance of cell-type-specific labeling of neurons in future studies.

Local Neuronal Responses to Intracortical Microstimulation in Rats' Barrel Cortex Are Dependent on Behavioral Context

The goal of cortical neuroprosthetics is to imprint sensory information as precisely as possible directly into cortical networks. Sensory processing, however, is dependent on the behavioral context. Therefore, a specific behavioral context may alter stimulation effects and, thus, perception. In this study, we reported how passive vs. active touch, i.e., the presence or absence of whisker movements, affects local field potential (LFP) responses to microstimulation in the barrel cortex in head-fixed behaving rats trained to move their whiskers voluntarily. The LFP responses to single-current pulses consisted of a short negative deflection corresponding to a volley of spike activity followed by a positive deflection lasting ~100 ms, corresponding to long-lasting suppression of spikes. Active touch had a characteristic effect on this response pattern. While the first phase including the negative peak remained stable, the later parts consisting of the positive peak were considerably suppressed. The stable phase varied systematically with the distance of the electrode from the stimulation site, pointing to saturation of neuronal responses to electrical stimulation in an intensity-dependent way. Our results suggest that modulatory effects known from normal sensory processing affect the response to cortical microstimulation as well. The network response to microstimulation is highly amenable to the behavioral state and must be considered for future approaches to imprint sensory signals into cortical circuits with neuroprostheses.

State-aware detection of sensory stimuli in the cortex of the awake mouse

Cortical responses to sensory inputs vary across repeated presentations of identical stimuli, but how this trial-to-trial variability impacts detection of sensory inputs is not fully understood. Using multi-channel local field potential (LFP) recordings in primary somatosensory cortex (S1) of the awake mouse, we optimized a data-driven cortical state classifier to predict single-trial sensory-evoked responses, based on features of the spontaneous, ongoing LFP recorded across cortical layers. Our findings show that, by utilizing an ongoing prediction of the sensory response generated by this state classifier, an ideal observer improves overall detection accuracy and generates robust detection of sensory inputs across various states of ongoing cortical activity in the awake brain, which could have implications for variability in the performance of detection tasks across brain states.

Establishing the link between neural activity and behavior is a central goal of neuroscience. One context in which to examine this link is in a sensory detection task, in which an animal is trained to report the presence of a barely perceptible sensory stimulus. In such tasks, both sensory responses in the brain and behavioral responses are highly variable. A simple hypothesis, originating in signal detection theory, is that perceived inputs generate neural activity that cross some threshold for detection. According to this hypothesis, sensory response variability would predict behavioral variability, but previous studies have not born out this prediction. Further complicating the picture, sensory response variability is partially dependent on the ongoing state of cortical activity, and we wondered whether this could resolve the mismatch between response variability and behavioral variability. Here, we use a computational approach to study an adaptive observer that utilizes an ongoing prediction of sensory responsiveness to detect sensory inputs. This observer has higher overall accuracy than the standard ideal observer. Moreover, because of the adaptation, the observer breaks the direct link between neural and behavioral variability, which could resolve discrepancies arising in past studies. We suggest new experiments to test our theory.

Sensory percepts induced by microwire array and DBS microstimulation in human sensory thalamus

BACKGROUND

Microstimulation in human sensory thalamus (ventrocaudal, VC) results in focal sensory percepts in the hand and arm which may provide an alternative target site (to somatosensory cortex) for the input of prosthetic sensory information. Sensory feedback to facilitate motor function may require simultaneous or timed responses across separate digits to recreate perceptions of slip as well as encoding of intensity variations in pressure or touch.

OBJECTIVES

To determine the feasibility of evoking sensory percepts on separate digits with variable intensity through either a microwire array or deep brain stimulation (DBS) electrode, recreating "natural" and scalable percepts relating to the arm and hand.

METHODS

We compared microstimulation within ventrocaudal sensory thalamus through either a 16-channel microwire array (∼400 kΩ per channel) or a 4-channel DBS electrode (∼1.2 kΩ per contact) for percept location, size, intensity, and quality sensation, during thalamic DBS electrode placement in patients with essential tremor.

RESULTS

Percepts in small hand or finger regions were evoked by microstimulation through individual microwires and in 5/6 patients sensation on different digits could be perceived from stimulation through separate microwires. Microstimulation through DBS electrode contacts evoked sensations over larger areas in 5/5 patients, and the apparent intensity of the perceived response could be modulated with stimulation amplitude. The perceived naturalness of the sensation depended both on the pattern of stimulation as well as intensity of the stimulation.

CONCLUSIONS

Producing consistent evoked perceptions across separate digits within sensory thalamus is a feasible concept and a compact alternative to somatosensory cortex microstimulation for prosthetic sensory feedback. This approach will require a multi-element low impedance electrode with a sufficient stimulation range to evoke variable intensities of perception and a predictable spread of contacts to engage separate digits.

State-Dependent Decoding Algorithms Improve the Performance of a Bidirectional BMI in Anesthetized Rats

Brain-machine interfaces (BMIs) promise to improve the quality of life of patients suffering from sensory and motor disabilities by creating a direct communication channel between the brain and the external world. Yet, their performance is currently limited by the relatively small amount of information that can be decoded from neural activity recorded form the brain. We have recently proposed that such decoding performance may be improved when using state-dependent decoding algorithms that predict and discount the large component of the trial-to-trial variability of neural activity which is due to the dependence of neural responses on the network's current internal state. Here we tested this idea by using a bidirectional BMI to investigate the gain in performance arising from using a state-dependent decoding algorithm. This BMI, implemented in anesthetized rats, controlled the movement of a dynamical system using neural activity decoded from motor cortex and fed back to the brain the dynamical system's position by electrically microstimulating somatosensory cortex. We found that using state-dependent algorithms that tracked the dynamics of ongoing activity led to an increase in the amount of information extracted form neural activity by 22%, with a consequently increase in all of the indices measuring the BMI's performance in controlling the dynamical system. This suggests that state-dependent decoding algorithms may be used to enhance BMIs at moderate computational cost.

Implications of the Dependence of Neuronal Activity on Neural Network States for the Design of Brain-Machine Interfaces

Brain-machine interfaces (BMIs) can improve the quality of life of patients with sensory and motor disabilities by both decoding motor intentions expressed by neural activity, and by encoding artificially sensed information into patterns of neural activity elicited by causal interventions on the neural tissue. Yet, current BMIs can exchange relatively small amounts of information with the brain. This problem has proved difficult to overcome by simply increasing the number of recording or stimulating electrodes, because trial-to-trial variability of neural activity partly arises from intrinsic factors (collectively known as the network state) that include ongoing spontaneous activity and neuromodulation, and so is shared among neurons. Here we review recent progress in characterizing the state dependence of neural responses, and in particular of how neural responses depend on endogenous slow fluctuations of network excitability. We then elaborate on how this knowledge may be used to increase the amount of information that BMIs exchange with brain. Knowledge of network state can be used to fine-tune the stimulation pattern that should reliably elicit a target neural response used to encode information in the brain, and to discount part of the trial-by-trial variability of neural responses, so that they can be decoded more accurately.

Implantable neurotechnologies: electrical stimulation and applications

Neural stimulation using injected electrical charge is widely used both in functional therapies and as an experimental tool for neuroscience applications. Electrical pulses can induce excitation of targeted neural pathways that aid in the treatment of neural disorders or dysfunction of the central and peripheral nervous system. In this review, we summarize the recent trends in the field of electrical stimulation for therapeutic interventions of nervous system disorders, such as for the restoration of brain, eye, ear, spinal cord, nerve and muscle function. Neural prosthetic applications are discussed, and functional electrical stimulation parameters for treating such disorders are reviewed. Important considerations for implantable packaging and enhancing device reliability are also discussed. Neural stimulators are expected to play a profound role in implantable neural devices that treat disorders and help restore functions in injured or disabled nervous system.

Estimation of diffusion coefficients from voltammetric signals by support vector and gaussian process regression

BACKGROUND

Support vector regression (SVR) and Gaussian process regression (GPR) were used for the analysis of electroanalytical experimental data to estimate diffusion coefficients.

RESULTS

For simulated cyclic voltammograms based on the EC, Eqr, and EqrC mechanisms these regression algorithms in combination with nonlinear kernel/covariance functions yielded diffusion coefficients with higher accuracy as compared to the standard approach of calculating diffusion coefficients relying on the Nicholson-Shain equation. The level of accuracy achieved by SVR and GPR is virtually independent of the rate constants governing the respective reaction steps. Further, the reduction of high-dimensional voltammetric signals by manual selection of typical voltammetric peak features decreased the performance of both regression algorithms compared to a reduction by downsampling or principal component analysis. After training on simulated data sets, diffusion coefficients were estimated by the regression algorithms for experimental data comprising voltammetric signals for three organometallic complexes.

CONCLUSIONS

Estimated diffusion coefficients closely matched the values determined by the parameter fitting method, but reduced the required computational time considerably for one of the reaction mechanisms. The automated processing of voltammograms according to the regression algorithms yields better results than the conventional analysis of peak-related data.

A tensor-product-kernel framework for multiscale neural activity decoding and control

Brain machine interfaces (BMIs) have attracted intense attention as a promising technology for directly interfacing computers or prostheses with the brain's motor and sensory areas, thereby bypassing the body. The availability of multiscale neural recordings including spike trains and local field potentials (LFPs) brings potential opportunities to enhance computational modeling by enriching the characterization of the neural system state. However, heterogeneity on data type (spike timing versus continuous amplitude signals) and spatiotemporal scale complicates the model integration of multiscale neural activity. In this paper, we propose a tensor-product-kernel-based framework to integrate the multiscale activity and exploit the complementary information available in multiscale neural activity. This provides a common mathematical framework for incorporating signals from different domains. The approach is applied to the problem of neural decoding and control. For neural decoding, the framework is able to identify the nonlinear functional relationship between the multiscale neural responses and the stimuli using general purpose kernel adaptive filtering. In a sensory stimulation experiment, the tensor-product-kernel decoder outperforms decoders that use only a single neural data type. In addition, an adaptive inverse controller for delivering electrical microstimulation patterns that utilizes the tensor-product kernel achieves promising results in emulating the responses to natural stimulation.

System identification of the nonlinear dynamics in the thalamocortical circuit in response to patterned thalamic microstimulation in vivo

OBJECTIVE

Nonlinear system identification approaches were used to develop a dynamical model of the network level response to patterns of microstimulation in vivo.

APPROACH

The thalamocortical circuit of the rodent vibrissa pathway was the model system, with voltage sensitive dye imaging capturing the cortical response to patterns of stimulation delivered from a single electrode in the ventral posteromedial thalamus. The results of simple paired stimulus experiments formed the basis for the development of a phenomenological model explicitly containing nonlinear elements observed experimentally. The phenomenological model was fit using datasets obtained with impulse train inputs, Poisson-distributed in time and uniformly varying in amplitude.

MAIN RESULTS

The phenomenological model explained 58% of the variance in the cortical response to out of sample patterns of thalamic microstimulation. Furthermore, while fit on trial-averaged data, the phenomenological model reproduced single trial response properties when simulated with noise added into the system during stimulus presentation. The simulations indicate that the single trial response properties were dependent on the relative sensitivity of the static nonlinearities in the two stages of the model, and ultimately suggest that electrical stimulation activates local circuitry through linear recruitment, but that this activity propagates in a highly nonlinear fashion to downstream targets.

SIGNIFICANCE

The development of nonlinear dynamical models of neural circuitry will guide information delivery for sensory prosthesis applications, and more generally reveal properties of population coding within neural circuits.

Flexible charge balanced stimulator with 5.6 fC accuracy for 140 nC injections

Electrical stimulations of neuronal structures must ensure net injected charges to be zero for biological safety and voltage compliance reasons. We present a novel architecture of general purpose biphasic constant current stimulator that exhibits less than 5.6 fC error while injecting 140 nC charges using 1.4 mA currents. The floating current sources and conveyor switch based system can operate in monopolar or bipolar modes. Anodic-first or cathodic-first pulses with optional inter-phase delays have been demonstrated with zero quiescent current requirements at the analog front-end. The architecture eliminates blocking capacitors, electrode shorting and complex feedbacks. Bench-top and in-vivo measurement results have been presented with emulated electrode impedances (resistor-capacitor network), Ag-AgCl electrodes in saline and in-vivo (acute) peripheral nerve stimulations in anesthetized rats.

Adaptive inverse control of neural spatiotemporal spike patterns with a reproducing kernel Hilbert space (RKHS) framework

The precise control of spiking in a population of neurons via applied electrical stimulation is a challenge due to the sparseness of spiking responses and neural system plasticity. We pose neural stimulation as a system control problem where the system input is a multidimensional time-varying signal representing the stimulation, and the output is a set of spike trains; the goal is to drive the output such that the elicited population spiking activity is as close as possible to some desired activity, where closeness is defined by a cost function. If the neural system can be described by a time-invariant (homogeneous) model, then offline procedures can be used to derive the control procedure; however, for arbitrary neural systems this is not tractable. Furthermore, standard control methodologies are not suited to directly operate on spike trains that represent both the target and elicited system response. In this paper, we propose a multiple-input multiple-output (MIMO) adaptive inverse control scheme that operates on spike trains in a reproducing kernel Hilbert space (RKHS). The control scheme uses an inverse controller to approximate the inverse of the neural circuit. The proposed control system takes advantage of the precise timing of the neural events by using a Schoenberg kernel defined directly in the space of spike trains. The Schoenberg kernel maps the spike train to an RKHS and allows linear algorithm to control the nonlinear neural system without the danger of converging to local minima. During operation, the adaptation of the controller minimizes a difference defined in the spike train RKHS between the system and the target response and keeps the inverse controller close to the inverse of the current neural circuit, which enables adapting to neural perturbations. The results on a realistic synthetic neural circuit show that the inverse controller based on the Schoenberg kernel outperforms the decoding accuracy of other models based on the conventional rate representation of neural signal (i.e., spikernel and generalized linear model). Moreover, after a significant perturbation of the neuron circuit, the control scheme can successfully drive the elicited responses close to the original target responses.

Optimizing microstimulation using a reinforcement learning framework

The ability to provide sensory feedback is desired to enhance the functionality of neuroprosthetics. Somatosensory feedback provides closed-loop control to the motor system, which is lacking in feedforward neuroprosthetics. In the case of existing somatosensory function, a template of the natural response can be used as a template of desired response elicited by electrical microstimulation. In the case of no initial training data, microstimulation parameters that produce responses close to the template must be selected in an online manner. We propose using reinforcement learning as a framework to balance the exploration of the parameter space and the continued selection of promising parameters for further stimulation. This approach avoids an explicit model of the neural response from stimulation. We explore a preliminary architecture--treating the task as a k-armed bandit--using offline data recorded for natural touch and thalamic microstimulation, and we examine the methods efficiency in exploring the parameter space while concentrating on promising parameter forms. The best matching stimulation parameters, from k = 68 different forms, are selected by the reinforcement learning algorithm consistently after 334 realizations.

Subspace matching thalamic microstimulation to tactile evoked potentials in rat somatosensory cortex

We show experimental results that the evoked local field potentials of the rat somatosensory cortex from natural tactile touch of forepaw digits and matched thalamic microstimulation can be qualitatively and quantitively similar. In ongoing efforts to optimize the microstimulation settings (e.g., location, amplitude, etc.) to match the natural response, we investigate whether subspace projection methods, specifically the eigenface approach proposed in the computer vision community (Turk and Pentland 1991 [1]), can be used to choose the parameters of microstimulation such that the response matches a single tactile touch realization. Since the evoked potentials from multiple electrodes are high dimensional spatio-temporal data, the subspace projections improve computational efficiency and can reduce the effect of noisy realizations. In computing the PCA projections we use the peristimulus averages instead of the realizations. The dataset is pruned of unreliable stimulation types. A new subspace is computed for the pruned stimulation type, and is used to estimate a sequence of microstimulations to best match the natural responses. This microstimulation sequence is applied in vivo and quantitative analysis shows that per realization matching does statistically better than choosing randomly from the pruned subset.