Recently Haas et al. (J Neurophysiol 96: 3305-3313, 2006), observed a novel form of spike timing dependent plasticity (iSTDP) in GABAergic synaptic couplings in layer II of the entorhinal cortex. Depending on the relative timings of the presynaptic input at time t (pre) and the postsynaptic excitation at time t (post), the synapse is strengthened (Deltat = t(post) - t(pre) > 0) or weakened (Deltat < 0). The temporal dynamic range of the observed STDP rule was found to lie in the higher gamma frequency band (> or =40 Hz), a frequency range important for several vital neuronal tasks. In this paper we study the function of this novel form of iSTDP in the synchronization of the inhibitory neuronal network. In particular we consider a network of two unidirectionally coupled interneurons (UCI) and two mutually coupled interneurons (MCI), in the presence of heterogeneity in the intrinsic firing rates of each coupled neuron. Using the method of spike time response curve (STRC), we show how iSTDP influences the dynamics of the coupled neurons, such that the pair synchronizes under moderately large heterogeneity in the firing rates. Using the general properties of the STRC for a Type-1 neuron model (Ermentrout, Neural Comput 8:979-1001, 1996) and the observed iSTDP we determine conditions on the initial configuration of the UCI network that would result in 1:1 in-phase synchrony between the two coupled neurons. We then demonstrate a similar enhancement of synchrony in the MCI with dynamic synaptic modulation. For the MCI we also consider heterogeneity introduced in the network through the synaptic parameters: the synaptic decay time of mutual inhibition and the self inhibition synaptic strength. We show that the MCI exhibits enhanced synchrony in the presence of all the above mentioned sources of heterogeneity and the mechanism for this enhanced synchrony is similar to the case of the UCI.

No abstract

We conclude with a discussion on the utility of computational models for the development of new epilepsy treatment protocols.

Channelrhodopsins-2 (ChR2) are a class of light sensitive proteins that offer the ability to use light stimulation to regulate neural activity with millisecond precision. In order to address the limitations in the efficacy of the wild-type ChR2 (ChRwt) to achieve this objective, new variants of ChR2 that exhibit fast mon-exponential photocurrent decay characteristics have been recently developed and validated. In this paper, we investigate whether the framework of transition rate model with 4 states, primarily developed to mimic the biexponential photocurrent decay kinetics of ChRwt, as opposed to the low complexity 3 state model, is warranted to mimic the mono-exponential photocurrent decay kinetics of the newly developed fast ChR2 variants: ChETA (Gunaydin et al., Nature Neurosci. 13:387-392, 2010) and ChRET/TC (Berndt et al., Proc. Natl. Acad. Sci. 108:7595-7600, 2011). We begin by estimating the parameters of the 3-state and 4-state models from experimental data on the photocurrent kinetics of ChRwt, ChETA, and ChRET/TC. We then incorporate these models into a fast-spiking interneuron model (Wang and Buzsaki, J. Neurosci. 16:6402-6413, 1996) and a hippocampal pyramidal cell model (Golomb et al., J. Neurophysiol. 96:1912-1926, 2006) and investigate the extent to which the experimentally observed neural response to various optostimulation protocols can be captured by these models. We demonstrate that for all ChR2 variants investigated, the 4 state model implementation is better able to capture neural response consistent with experiments across wide range of optostimulation protocol. We conclude by analytically investigating the conditions under which the characteristic specific to the 3-state model, namely the monoexponential photocurrent decay of the newly developed variants of ChR2, can occur in the framework of the 4-state model.

An understanding of the in vivo spatial emergence of abnormal brain activity during spontaneous seizure onset is critical to future early seizure detection and closed-loop seizure prevention therapies. In this study, we use Granger causality (GC) to determine the strength and direction of relationships between local field potentials (LFPs) recorded from bilateral microelectrode arrays in an intermittent spontaneous seizure model of chronic temporal lobe epilepsy before, during, and after Racine grade partial onset generalized seizures. Our results indicate distinct patterns of directional GC relationships within the hippocampus, specifically from the CA1 subfield to the dentate gryus, prior to and during seizure onset. Our results suggest sequential and hierarchical temporal relationships between the CA1 and dentate gyrus within and across hippocampal hemispheres during seizure. Additionally, our analysis suggests a reversal in the direction of GC relationships during seizure, from an abnormal pattern to more anatomically expected pattern. This reversal correlates well with the observed behavioral transition from tonic to clonic seizure in time-locked video. These findings highlight the utility of GC to reveal dynamic directional temporal relationships between multichannel LFP recordings from multiple brain regions during unprovoked spontaneous seizures.

The motor pathway responsible for the complex vocalizations of songbirds has been extensively characterized, both in terms of intrinsic and synaptic physiology in vitro and in terms of the spatiotemporal patterns of neural activity in vivo. However, the relationship between the neural architecture of the song motor pathway and the acoustic features of birdsong is not well understood. Using a computational model of the song motor pathway and the songbird vocal organ, we investigate the relationship between song production and the neural connectivity of nucleus HVc (used as a proper name) and the robust nucleus of the archistriatum (RA). Drawing on recent experimental observations, our neural model contains a population of sequentially bursting HVc neurons driving the activity of a population of RA neurons. An important focus of our investigations is the contribution of intrinsic circuitry within RA to the acoustic output of the model. We find that the inclusion of inhibitory interneurons in the model can substantially influence the features of song syllables, and we illustrate the potential for subharmonic behavior in RA in response to forcing by HVc neurons. Our results demonstrate the association of specific acoustic features with specific neural connectivities and support the view that intrinsic circuitry within RA may play a critical role in generating the features of birdsong.

We automated a traditionally labor-intensive, yet widely-used capsule production system.

We compare the performance of three support vector machine (SVM) types: weighted SVM, one-class SVM and support vector data description (SVDD) for the application of seizure detection in an animal model of chronic epilepsy. Large EEG datasets (273 h and 91 h respectively, with a sampling rate of 1 kHz) from two groups of rats with chronic epilepsy were used in this study. For each of these EEG datasets, we extracted three energy-based seizure detection features: mean energy, mean curve length and wavelet energy. Using these features we performed twofold cross-validation to obtain the performance statistics: sensitivity (S), specificity (K) and detection latency (tau) as a function of control parameters for the given SVM. Optimal control parameters for each SVM type that produced the best seizure detection statistics were then identified using two independent strategies. Performance of each SVM type is ranked based on the overall seizure detection performance through an optimality index metric (O). We found that SVDD not only performed better than the other SVM types in terms of highest value of the mean optimality index metric (O⁻) but also gave a more reliable performance across the two EEG datasets.