We review different aspects of the simulation of spiking neural networks. We start by reviewing the different types of simulation strategies and algorithms that are currently implemented. We next review the precision of those simulation strategies, in particular in cases where plasticity depends on the exact timing of the spikes. We overview different simulators and simulation environments presently available (restricted to those freely available, open source and documented). For each simulation tool, its advantages and pitfalls are reviewed, with an aim to allow the reader to identify which simulator is appropriate for a given task. Finally, we provide a series of benchmark simulations of different types of networks of spiking neurons, including Hodgkin-Huxley type, integrate-andfire models, interacting with current-based or conductance-based synapses, using clock-driven or NIH-PA Author ManuscriptNIH-PA Author Manuscript NIH-PA Author Manuscript event-driven integration strategies. The same set of models are implemented on the different simulators, and the codes are made available. The ultimate goal of this review is to provide a resource to facilitate identifying the appropriate integration strategy and simulation tool to use for a given modeling problem related to spiking neural networks.

We designed a model-based analysis to predict the occurrence of population patterns in distributed spiking activity. Using a maximum entropy principle with a Markovian assumption, we obtain a model that accounts for both spatial and temporal pairwise correlations among neurons. This model is tested on data generated with a Glauber spin-glass system and is shown to correctly predict the occurrence probabilities of spatiotemporal patterns significantly better than Ising models only based on spatial correlations. This increase of predictability was also observed on experimental data recorded in parietal cortex during slow-wave sleep. This approach can also be used to generate surrogates that reproduce the spatial and temporal correlations of a given data set.

Plasticity of cortical responses in vivo involves activity-dependent changes at synapses, but the manner in which different forms of synaptic plasticity act together to create functional changes in neurons remains unknown. We found that spike timing-induced receptive field plasticity of visual cortex neurons in mice is anchored by increases in the synaptic strength of identified spines. This is accompanied by a decrease in the strength of adjacent spines on a slower time scale. The locally coordinated potentiation and depression of spines involves prominent AMPA receptor redistribution via targeted expression of the immediate early gene product Arc. Hebbian strengthening of activated synapses and heterosynaptic weakening of adjacent synapses thus cooperatively orchestrate cell-wide plasticity of functional neuronal responses.

As in other sensory modalities, one function of the somatosensory system is to detect coherence and contrast in the environment. To investigate the neural bases of these computations, we applied different spatiotemporal patterns of stimuli to rat whiskers while recording multiple neurons in the barrel cortex. Model-based analysis of the responses revealed different coding schemes according to the level of input correlation. With uncorrelated stimuli on 24 whiskers, we identified two distinct functional categories of neurons, analogous in the temporal domain to simple and complex cells of the primary visual cortex. With correlated stimuli, however, a complementary coding scheme emerged: two distinct cell populations, similar to reinforcing and antagonist neurons described in the higher visual area MT, responded specifically to correlations. We suggest that similar context-dependent coexisting coding strategies may be present in other sensory systems to adapt sensory integration to specific stimulus statistics.