Influence of Autapses on Synchronization in Neural Networks With Chemical Synapses

A great deal of research has been devoted on the investigation of neural dynamics in various network topologies. However, only a few studies have focused on the influence of autapses, synapses from a neuron onto itself via closed loops, on neural synchronization. Here, we build a random network with adaptive exponential integrate-and-fire neurons coupled with chemical synapses, equipped with autapses, to study the effect of the latter on synchronous behavior. We consider time delay in the conductance of the pre-synaptic neuron for excitatory and inhibitory connections. Interestingly, in neural networks consisting of both excitatory and inhibitory neurons, we uncover that synchronous behavior depends on their synapse type. Our results provide evidence on the synchronous and desynchronous activities that emerge in random neural networks with chemical, inhibitory and excitatory synapses where neurons are equipped with autapses.

Adaptive integrate-and-fire model reproduces the dynamics of olfactory receptor neuron responses in a moth

In order to understand how olfactory stimuli are encoded and processed in the brain, it is important to build a computational model for olfactory receptor neurons (ORNs). Here, we present a simple and reliable mathematical model of a moth ORN generating spikes. The model incorporates a simplified description of the chemical kinetics leading to olfactory receptor activation and action potential generation. We show that an adaptive spike threshold regulated by prior spike history is an effective mechanism for reproducing the typical phasic-tonic time course of ORN responses. Our model reproduces the response dynamics of individual neurons to a fluctuating stimulus that approximates odorant fluctuations in nature. The parameters of the spike threshold are essential for reproducing the response heterogeneity in ORNs. The model provides a valuable tool for efficient simulations of olfactory circuits.

Models for Spiking Neurons: Integrate-and-Fire Units and Relaxation Oscillators

Relaxation oscillators are nonlinear electronic circuits that produce a repetitive non-sinusoidal waveform when sufficient voltage is applied. In this fashion, they are reminiscent of integrate-and-fire neuron models, except that they also include components with hysteresis, and thus require no threshold rule to determine when an impulse has occurred or to return the voltage to its reset value. Here, I discuss the pros and cons of teaching elementary neurophysiology using first-order linear integrate-and-fire neurons versus relaxation oscillator circuits. I suggest that the shortcomings of both types of models are useful in order to foster a critical understanding of the neurophysiology underlying the firing dynamics of biological neurons.

Decoding of the spike timing of primary afferents during voluntary arm movements in monkeys

Understanding the mechanisms of encoding forelimb kinematics in the activity of peripheral afferents is essential for developing a somatosensory neuroprosthesis. To investigate whether the spike timing of dorsal root ganglion (DRG) neurons could be estimated from the forelimb kinematics of behaving monkeys, we implanted two multi-electrode arrays chronically in the DRGs at the level of the cervical segments in two monkeys. Neuronal activity during voluntary reach-to-grasp movements were recorded simultaneously with the trajectories of hand/arm movements, which were tracked in three-dimensional space using a motion capture system. Sixteen and 13 neurons, including muscle spindles, skin receptors, and tendon organ afferents, were recorded in the two monkeys, respectively. We were able to reconstruct forelimb joint kinematics from the temporal firing pattern of a subset of DRG neurons using sparse linear regression (SLiR) analysis, suggesting that DRG neuronal ensembles encoded information about joint kinematics. Furthermore, we estimated the spike timing of the DRG neuronal ensembles from joint kinematics using an integrate-and-fire model (IF) incorporating the SLiR algorithm. The temporal change of firing frequency of a subpopulation of neurons was reconstructed precisely from forelimb kinematics using the SLiR. The estimated firing pattern of the DRG neuronal ensembles encoded forelimb joint angles and velocities as precisely as the originally recorded neuronal activity. These results suggest that a simple model can be used to generate an accurate estimate of the spike timing of DRG neuronal ensembles from forelimb joint kinematics, and is useful for designing a proprioceptive decoder in a brain machine interface.

Cancer risk at low doses of ionizing radiation: artificial neural networks inference from atomic bomb survivors

Cancer risk at low doses of ionizing radiation remains poorly defined because of ambiguity in the quantitative link to doses below 0.2 Sv in atomic bomb survivors in Hiroshima and Nagasaki arising from limitations in the statistical power and information available on overall radiation dose. To deal with these difficulties, a novel nonparametric statistics based on the 'integrate-and-fire' algorithm of artificial neural networks was developed and tested in cancer databases established by the Radiation Effects Research Foundation. The analysis revealed unique features at low doses that could not be accounted for by nominal exposure dose, including (i) the presence of a threshold that varied with organ, gender and age at exposure, and (ii) a small but significant bumping increase in cancer risk at low doses in Nagasaki that probably reflects internal exposure to (239)Pu. The threshold was distinct from the canonical definition of zero effect in that it was manifested as negative excess relative risk, or suppression of background cancer rates. Such a unique tissue response at low doses of radiation exposure has been implicated in the context of the molecular basis of radiation-environment interplay in favor of recently emerging experimental evidence on DNA double-strand break repair pathway choice and its epigenetic memory by histone marking.

Patterns of interval correlations in neural oscillators with adaptation

Neural firing is often subject to negative feedback by adaptation currents. These currents can induce strong correlations among the time intervals between spikes. Here we study analytically the interval correlations of a broad class of noisy neural oscillators with spike-triggered adaptation of arbitrary strength and time scale. Our weak-noise theory provides a general relation between the correlations and the phase-response curve (PRC) of the oscillator, proves anti-correlations between neighboring intervals for adapting neurons with type I PRC and identifies a single order parameter that determines the qualitative pattern of correlations. Monotonically decaying or oscillating correlation structures can be related to qualitatively different voltage traces after spiking, which can be explained by the phase plane geometry. At high firing rates, the long-term variability of the spike train associated with the cumulative interval correlations becomes small, independent of model details. Our results are verified by comparison with stochastic simulations of the exponential, leaky, and generalized integrate-and-fire models with adaptation.

How adaptation currents change threshold, gain, and variability of neuronal spiking

Many types of neurons exhibit spike rate adaptation, mediated by intrinsic slow K(+) currents, which effectively inhibit neuronal responses. How these adaptation currents change the relationship between in vivo like fluctuating synaptic input, spike rate output, and the spike train statistics, however, is not well understood. In this computational study we show that an adaptation current that primarily depends on the subthreshold membrane voltage changes the neuronal input-output relationship (I-O curve) subtractively, thereby increasing the response threshold, and decreases its slope (response gain) for low spike rates. A spike-dependent adaptation current alters the I-O curve divisively, thus reducing the response gain. Both types of an adaptation current naturally increase the mean interspike interval (ISI), but they can affect ISI variability in opposite ways. A subthreshold current always causes an increase of variability while a spike-triggered current decreases high variability caused by fluctuation-dominated inputs and increases low variability when the average input is large. The effects on I-O curves match those caused by synaptic inhibition in networks with asynchronous irregular activity, for which we find subtractive and divisive changes caused by external and recurrent inhibition, respectively. Synaptic inhibition, however, always increases the ISI variability. We analytically derive expressions for the I-O curve and ISI variability, which demonstrate the robustness of our results. Furthermore, we show how the biophysical parameters of slow K(+) conductances contribute to the two different types of an adaptation current and find that Ca(2+)-activated K(+) currents are effectively captured by a simple spike-dependent description, while muscarine-sensitive or Na(+)-activated K(+) currents show a dominant subthreshold component.

Discharge regularity in the turtle posterior crista: comparisons between experiment and theory

Intra-axonal recordings were made from bouton fibers near their termination in the turtle posterior crista. Spike discharge, miniature excitatory postsynaptic potentials (mEPSPs), and afterhyperpolarizations (AHPs) were monitored during resting activity in both regularly and irregularly discharging units. Quantal size (qsize) and quantal rate (qrate) were estimated by shot-noise theory. Theoretically, the ratio, σV/(dμV/dt), between synaptic noise (σV) and the slope of the mean voltage trajectory (dμV/dt) near threshold crossing should determine discharge regularity. AHPs are deeper and more prolonged in regular units; as a result, dμV/dt is larger, the more regular the discharge. The qsize is larger and qrate smaller in irregular units; these oppositely directed trends lead to little variation in σV with discharge regularity. Of the two variables, dμV/dt is much more influential than the nearly constant σV in determining regularity. Sinusoidal canal-duct indentations at 0.3 Hz led to modulations in spike discharge and synaptic voltage. Gain, the ratio between the amplitudes of the two modulations, and phase leads re indentation of both modulations are larger in irregular units. Gain variations parallel the sensitivity of the postsynaptic spike encoder, the set of conductances that converts synaptic input into spike discharge. Phase variations reflect both synaptic inputs to the encoder and postsynaptic processes. Experimental data were interpreted using a stochastic integrate-and-fire model. Advantages of an irregular discharge include an enhanced encoder gain and the prevention of nonlinear phase locking. Regular and irregular units are more efficient, respectively, in the encoding of low- and high-frequency head rotations, respectively.

How spike generation mechanisms determine the neuronal response to fluctuating inputs

This study examines the ability of neurons to track temporally varying inputs, namely by investigating how the instantaneous firing rate of a neuron is modulated by a noisy input with a small sinusoidal component with frequency (f). Using numerical simulations of conductance-based neurons and analytical calculations of one-variable nonlinear integrate-and-fire neurons, we characterized the dependence of this modulation on f. For sufficiently high noise, the neuron acts as a low-pass filter. The modulation amplitude is approximately constant for frequencies up to a cutoff frequency, fc, after which it decays. The cutoff frequency increases almost linearly with the firing rate. For higher frequencies, the modulation amplitude decays as C/falpha, where the power alpha depends on the spike initiation mechanism. For conductance-based models, alpha = 1, and the prefactor C depends solely on the average firing rate and a spike "slope factor," which determines the sharpness of the spike initiation. These results are attributable to the fact that near threshold, the sodium activation variable can be approximated by an exponential function. Using this feature, we propose a simplified one-variable model, the "exponential integrate-and-fire neuron," as an approximation of a conductance-based model. We show that this model reproduces the dynamics of a simple conductance-based model extremely well. Our study shows how an intrinsic neuronal property (the characteristics of fast sodium channels) determines the speed with which neurons can track changes in input.

Mathematical models of cochlear nucleus onset neurons: I. Point neuron with many weak synaptic inputs

The cochlear nucleus (CN) presents a unique opportunity for quantitatively studying input-output transformations by neurons because it gives rise to a variety of different response types from a relatively homogeneous input source, the auditory nerve (AN). Particularly interesting among CN neurons are Onset (On) neurons, which have a prominent response to the onset of sustained sounds followed by little or no response in the steady-state. On neurons contrast sharply with their AN inputs, which respond vigorously throughout stimuli. On neurons can entrain to stimuli (firing once per cycle of a periodic stimulus) at up to 1000 Hz, unlike their AN inputs. To understand the mechanisms underlying these response patterns, we tested whether an integrate-to-threshold point-neuron model with a fixed refractory period can account for On discharge patterns for tones, systematically examining the effect of membrane time constant and the number and strength of the exclusively excitatory AN synaptic inputs. To produce both onset responses to high-frequency tone bursts and entrainment to a broad range of low-frequency tones, the model must have a short time constant ( approximately 0.125 ms) and a large number (>100) of weak synaptic inputs, properties that are consistent with the electrical properties and anatomy of On-responding cells. With these parameters, the model acts like a coincidence detector with a threshold-like relationship between the instantaneous discharge rates of the output and the inputs. Onset responses to high-frequency tone bursts result because the threshold effect enhances the initial response of the AN inputs and suppresses their relatively lower sustained response. However, when the model entrains across a broad range of frequencies, it also produces short interspike intervals at the onset of high-frequency tone bursts, a response pattern not found in all types of On neurons. These results show a tradeoff, that may be a general property of many neurons, between following rapid stimulus fluctuations and responding without short interspike intervals at the onset of sustained stimuli.

How noise contributes to contrast invariance of orientation tuning in cat visual cortex

The width of the orientation tuning curves of the spike response of neurons in V1 is invariant to contrast. This property constrains the possible mechanisms underlying orientation selectivity. It has been suggested that noise circumvents the iceberg effect that would prevent contrast invariance in the purely feedforward mechanism. Here we investigate systematically how noise contributes to the contrast invariance of orientation tuning curves in V1. We study three models of increasing complexity: a simple threshold-linear firing rate model, a leaky integrate-and-fire model, and a conductance-based model. We show that the noise transmutes the threshold nonlinearity of the input-output relationships into an approximate power law without a threshold within some firing rate range. This implies that, under certain conditions which are derived here, the tuning of the neuron output is approximately contrast invariant. In particular we show that this mechanism for contrast invariance requires that the neuron firing rate must not be too large and that increasing or lowering the contrast too much destroys this invariance. We also show that if this mechanism operates in V1, the spike response, R, and average voltage response V of the neurons in V1 should vary with the contrast, C, according to R(C)gamma proportional to V(C)gamma. The exponent gamma can be estimated from the amount by which the spike tuning curve is sharpened with respect to the voltage tuning curves of the neurons. This prediction does not depend on the specifics of the model and can be tested experimentally.