Circuits that contain the Current : I Potassium

Re-display model names without descriptions
    Models   Description
1. 3D model of the olfactory bulb (Migliore et al. 2014)
This entry contains a link to a full HD version of movie 1 and the NEURON code of the paper: "Distributed organization of a brain microcircuit analysed by three-dimensional modeling: the olfactory bulb" by M Migliore, F Cavarretta, ML Hines, and GM Shepherd.
2. 3D olfactory bulb: operators (Migliore et al, 2015)
"... Using a 3D model of mitral and granule cell interactions supported by experimental findings, combined with a matrix-based representation of glomerular operations, we identify the mechanisms for forming one or more glomerular units in response to a given odor, how and to what extent the glomerular units interfere or interact with each other during learning, their computational role within the olfactory bulb microcircuit, and how their actions can be formalized into a theoretical framework in which the olfactory bulb can be considered to contain "odor operators" unique to each individual. ..."
3. A 1000 cell network model for Lateral Amygdala (Kim et al. 2013)
1000 Cell Lateral Amygdala model for investigation of plasticity and memory storage during Pavlovian Conditioning.
4. A Model Circuit of Thalamocortical Convergence (Behuret et al. 2013)
“… Using dynamic-clamp techniques in thalamic slices in vitro, we combined theoretical and experimental approaches to implement a realistic hybrid retino-thalamo-cortical pathway mixing biological cells and simulated circuits. … The study of the impact of the simulated cortical input on the global retinocortical signal transfer efficiency revealed a novel control mechanism resulting from the collective resonance of all thalamic relay neurons. We show here that the transfer efficiency of sensory input transmission depends on three key features: i) the number of thalamocortical cells involved in the many-to-one convergence from thalamus to cortex, ii) the statistics of the corticothalamic synaptic bombardment and iii) the level of correlation imposed between converging thalamic relay cells. In particular, our results demonstrate counterintuitively that the retinocortical signal transfer efficiency increases when the level of correlation across thalamic cells decreases. …”
5. A Moth MGC Model-A HH network with quantitative rate reduction (Buckley & Nowotny 2011)
We provide the model used in Buckley & Nowotny (2011). It consists of a network of Hodgkin Huxley neurons coupled by slow GABA_B synapses which is run alongside a quantitative reduction described in the associated paper.
6. A multilayer cortical model to study seizure propagation across microdomains (Basu et al. 2015)
A realistic neural network was used to simulate a region of neocortex to obtain extracellular LFPs from ‘virtual micro-electrodes’ and produce test data for comparison with multisite microelectrode recordings. A model was implemented in the GENESIS neurosimulator. A simulated region of cortex was represented by layers 2/3, 5/6 (interneurons and pyramidal cells) and layer 4 stelate cells, spaced at 25 µm in each horizontal direction. Pyramidal cells received AMPA and NMDA inputs from neighboring cells at the basal and apical dendrites. The LFP data was generated by simulating 16-site electrode array with the help of ‘efield’ objects arranged at the predetermined positions with respect to the surface of the simulated network. The LFP for the model is derived from a weighted average of the current sources summed over all cellular compartments. Cell models were taken from from Traub et al. (2005) J Neurophysiol 93(4):2194-232.
7. A network model of tail withdrawal in Aplysia (White et al 1993)
The contributions of monosynaptic and polysynaptic circuitry to the tail-withdrawal reflex in the marine mollusk Aplysia californica were assessed by the use of physiologically based neural network models. Effects of monosynaptic circuitry were examined by the use of a two-layer network model with four sensory neurons in the input layer and one motor neuron in the output layer. Results of these simulations indicated that the monosynaptic circuit could not account fully for long-duration responses of tail motor neurons elicited by tail stimulation. A three-layer network model was constructed by interposing a layer of two excitatory interneurons between the input and output layers of the two-layer network model. The three-layer model could account for long-duration responses in motor neurons. Sensory neurons are a known site of plasticity in Aplysia. Synaptic plasticity at more than one locus modified dramatically the input-output relationship of the three-layer network model. This feature gave the model redundancy in its plastic properties and points to the possibility of distributed memory in the circuitry mediating withdrawal reflexes in Aplysia. Please see paper for more results and details.
8. A simulation method for the firing sequences of motor units (Jiang et al 2006)
" ... a novel model based on the Hodgkin–Huxley (HH) system is proposed, which has the ability to simulate the complex neurodynamics of the firing sequences of motor neurons. The model is presented at the cellular level and network level, and some simulation results from a simple 3-neuron network are presented to demonstrate its applications." See paper for more and details.
9. A single column thalamocortical network model (Traub et al 2005)
To better understand population phenomena in thalamocortical neuronal ensembles, we have constructed a preliminary network model with 3,560 multicompartment neurons (containing soma, branching dendrites, and a portion of axon). Types of neurons included superficial pyramids (with regular spiking [RS] and fast rhythmic bursting [FRB] firing behaviors); RS spiny stellates; fast spiking (FS) interneurons, with basket-type and axoaxonic types of connectivity, and located in superficial and deep cortical layers; low threshold spiking (LTS) interneurons, that contacted principal cell dendrites; deep pyramids, that could have RS or intrinsic bursting (IB) firing behaviors, and endowed either with non-tufted apical dendrites or with long tufted apical dendrites; thalamocortical relay (TCR) cells; and nucleus reticularis (nRT) cells. To the extent possible, both electrophysiology and synaptic connectivity were based on published data, although many arbitrary choices were necessary.
10. A two-layer biophysical olfactory bulb model of cholinergic neuromodulation (Li and Cleland 2013)
This is a two-layer biophysical olfactory bulb (OB) network model to study cholinergic neuromodulation. Simulations show that nicotinic receptor activation sharpens mitral cell receptive field, while muscarinic receptor activation enhances network synchrony and gamma oscillations. This general model suggests that the roles of nicotinic and muscarinic receptors in OB are both distinct and complementary to one another, together regulating the effects of ascending cholinergic inputs on olfactory bulb transformations.
11. A unified thalamic model of multiple distinct oscillations (Li, Henriquez and Fröhlich 2017)
We present a unified model of the thalamus that is capable of independently generating multiple distinct oscillations (delta, spindle, alpha and gamma oscillations) under different levels of acetylcholine (ACh) and norepinephrine (NE) modulation corresponding to different physiological conditions (deep sleep, light sleep, relaxed wakefulness and attention). The model also shows that entrainment of thalamic oscillations is state-dependent.
12. Activity patterns in a subthalamopallidal network of the basal ganglia model (Terman et al 2002)
"Based on recent experimental data, we have developed a conductance-based computational network model of the subthalamic nucleus and the external segment of the globus pallidus in the indirect pathway of the basal ganglia. Computer simulations and analysis of this model illuminate the roles of the coupling architecture of the network, and associated synaptic conductances, in modulating the activity patterns displayed by this network. Depending on the relationships of these coupling parameters, the network can support three general classes of sustained firing patterns: clustering, propagating waves, and repetitive spiking that may show little regularity or correlation. ...". Terman's XPP code and a partial implementation by Taylor Malone in NEURON and python are included.
13. Axonal gap junctions produce fast oscillations in cerebellar Purkinje cells (Traub et al. 2008)
Examines how electrical coupling between proximal axons produces fast oscillations in cerebellar Purkinje cells. Traub RD, Middleton SJ, Knopfel T, Whittington MA (2008) Model of very fast (>75 Hz) network oscillations generated by electrical coupling between the proximal axons of cerebellar Purkinje cells. European Journal of Neuroscience.
14. Basal ganglia-thalamic network model for deep brain stimulation (So et al. 2011)
This is a model of the basal ganglia-thalamic network, modified from the Rubin and Terman model (High frequency stimulation of the Subthalamic Nucleus, Rubin and Terman 2004). We subsequently used this model to investigate the effectiveness of STN and GPi DBS as well as lesion when various proportions of local cells and fibers of passage were activated or silenced. The BG network exhibited characteristics consistent with published experimental data, both on the level of single cells and on the network level. Perhaps most notably, and in contrast to the original RT model, the changes in the thalamic error index with changes in the DBS frequency matched well the changes in clinical symptoms with changes in DBS frequency.
15. Biophysically realistic neural modeling of the MEG mu rhythm (Jones et al. 2009)
"Variations in cortical oscillations in the alpha (7–14 Hz) and beta (15–29 Hz) range have been correlated with attention, working memory, and stimulus detection. The mu rhythm recorded with magnetoencephalography (MEG) is a prominent oscillation generated by Rolandic cortex containing alpha and beta bands. Despite its prominence, the neural mechanisms regulating mu are unknown. We characterized the ongoing MEG mu rhythm from a localized source in the finger representation of primary somatosensory (SI) cortex. Subjects showed variation in the relative expression of mu-alpha or mu-beta, which were nonoverlapping for roughly 50% of their respective durations on single trials. To delineate the origins of this rhythm, a biophysically principled computational neural model of SI was developed, with distinct laminae, inhibitory and excitatory neurons, and feedforward (FF, representative of lemniscal thalamic drive) and feedback (FB, representative of higher-order cortical drive or input from nonlemniscal thalamic nuclei) inputs defined by the laminar location of their postsynaptic effects. ..."
16. Burst induced synaptic plasticity in Apysia sensorimotor neurons (Phares et al 2003)
The Aplysia sensorimotor synapse is a key site of plasticity for several simple forms of learning. Intracellular stimulation of sensory neurons to fire a burst of action potentials at 10 Hz for 1 sec led to significant homosynaptic depression of postsynaptic responses. During the burst, the steady-state depressed phase of the postsynaptic response, which was only 20% of the initial EPSP of the burst, still contributed to firing the motor neuron. To explore the functional contribution of transient homosynaptic depression to the response of the motor neuron, computer simulations of the sensorimotor synapse with and without depression were compared. Depression allowed the motor neuron to produce graded responses over a wide range of presynaptic input strength. Thus, synaptic depression increased the dynamic range of the sensorimotor synapse and can, in principle, have a profound effect on information processing. Please see paper for results and details.
17. Bursting and oscillations in RD1 Retina driven by AII Amacrine Neuron (Choi et al. 2014)
"In many forms of retinal degeneration, photoreceptors die but inner retinal circuits remain intact. In the rd1 mouse, an established model for blinding retinal diseases, spontaneous activity in the coupled network of AII amacrine and ON cone bipolar cells leads to rhythmic bursting of ganglion cells. Since such activity could impair retinal and/or cortical responses to restored photoreceptor function, understanding its nature is important for developing treatments of retinal pathologies. Here we analyzed a compartmental model of the wild-type mouse AII amacrine cell to predict that the cell's intrinsic membrane properties, specifically, interacting fast Na and slow, M-type K conductances, would allow its membrane potential to oscillate when light-evoked excitatory synaptic inputs were withdrawn following photoreceptor degeneration. ..."
18. Ca+/HCN channel-dependent persistent activity in multiscale model of neocortex (Neymotin et al 2016)
"Neuronal persistent activity has been primarily assessed in terms of electrical mechanisms, without attention to the complex array of molecular events that also control cell excitability. We developed a multiscale neocortical model proceeding from the molecular to the network level to assess the contributions of calcium regulation of hyperpolarization-activated cyclic nucleotide-gated (HCN) channels in providing additional and complementary support of continuing activation in the network. ..."
19. CA1 network model for place cell dynamics (Turi et al 2019)
Biophysical model of CA1 hippocampal region. The model simulates place cells/fields and explores the place cell dynamics as function of VIP+ interneurons.
20. CA1 pyramidal cell: reconstructed axonal arbor and failures at weak gap junctions (Vladimirov 2011)
Model of pyramidal CA1 cells connected by gap junctions in their axons. Cell geometry is based on anatomical reconstruction of rat CA1 cell (NeuroMorpho.Org ID: NMO_00927) with long axonal arbor. Model init_2cells.hoc shows failures of second spike propagation in a spike doublet, depending on conductance of an axonal gap junction. Model init_ring.hoc shows that spike failure result in reentrant oscillations of a spike in a loop of axons connected by gap junctions, where one gap junction is weak. The paper shows that in random networks of axons connected by gap junctions, oscillations are driven by single pacemaker loop of axons. The shortest loop, around which a spike can travel, is the most likely pacemaker. This principle allows us to predict the frequency of oscillations from network connectivity and visa versa. We propose that this type of oscillations corresponds to so-called fast ripples in epileptic hippocampus.
21. Cerebellar granular layer (Maex and De Schutter 1998)
Circuit model of the granular layer representing a one-dimensional array of single-compartmental granule cells (grcs) and Golgi cells (Gocs). This paper examines the effects of feedback inhibition (grc -> Goc -> grc) versus feedforward inhibition (mossy fibre -> Goc -> grc) on synchronization and oscillatory behaviour.
22. Changes of ionic concentrations during seizure transitions (Gentiletti et al. 2016)
"... In order to investigate the respective roles of synaptic interactions and nonsynaptic mechanisms in seizure transitions, we developed a computational model of hippocampal cells, involving the extracellular space, realistic dynamics of Na+, K+, Ca2+ and Cl - ions, glial uptake and extracellular diffusion mechanisms. We show that the network behavior with fixed ionic concentrations may be quite different from the neurons’ behavior when more detailed modeling of ionic dynamics is included. In particular, we show that in the extended model strong discharge of inhibitory interneurons may result in long lasting accumulation of extracellular K+, which sustains the depolarization of the principal cells and causes their pathological discharges. ..."
23. Classic model of the Tritonia Swim CPG (Getting, 1989)
Classic model developed by Petter Getting of the 3-cell core CPG (DSI, C2, and VSI-B) mediating escape swimming in Tritonia diomedea. Cells use a hybrid integrate-and-fire scheme pioneered by Peter Getting. Each model cell is reconstructed from extensive physiological measurements to precisely mimic I-F curves, synaptic waveforms, and functional connectivity. **However, continued physiological measurements show that Getting may have inadvertently incorporated modulatory and or polysynaptic effects -- the properties of this model do *not* match physiological measurements in rested preparations.** This simulation reconstructs the Getting model as reported in: Getting (1989) 'Reconstruction of small neural networks' In Methods in Neural Modeling, 1st ed, p. 171-196. See also, an earlier version of this model reported in Getting (1983). Every attempt has been made to replicate the 1989 model as precisely as possible.
24. Collection of simulated data from a thalamocortical network model (Glabska, Chintaluri, Wojcik 2017)
"A major challenge in experimental data analysis is the validation of analytical methods in a fully controlled scenario where the justification of the interpretation can be made directly and not just by plausibility. ... One solution is to use simulations of realistic models to generate ground truth data. In neuroscience, creating such data requires plausible models of neural activity, access to high performance computers, expertise and time to prepare and run the simulations, and to process the output. To facilitate such validation tests of analytical methods we provide rich data sets including intracellular voltage traces, transmembrane currents, morphologies, and spike times. ... The data were generated using the largest publicly available multicompartmental model of thalamocortical network (Traub et al. 2005), with activity evoked by different thalamic stimuli."
25. Competing oscillator 5-cell circuit and Parameterscape plotting (Gutierrez et al. 2013)
Our 5-cell model consists of competing fast and slow oscillators connected to a hub neuron with electrical and inhibitory synapses. Motivated by the Stomatogastric Ganglion (STG) circuit in the crab, we explored the patterns of coordination in the network as a function of the electrical coupling and inhibitory synapse strengths with the help of a novel visualization method that we call the "Parameterscape." The code submitted here will allow you to run circuit simulations and to produce a Parameterscape with the results.
26. Competition for AP initiation sites in a circuit controlling simple learning (Cruz et al. 2007)
"The spatial and temporal patterns of action potential initiations were studied in a behaving leech preparation to determine the basis of increased firing that accompanies sensitization, a form of non-associative learning requiring the S-interneurons. ... The S-interneurons, one in each ganglion and linked by electrical synapses with both neighbors to form a chain, are interposed between sensory and motor neurons. ... the single site with the largest initiation rate, the S-cell in the stimulated segment, suppressed initiations in adjacent ganglia. Experiments showed this was both because (1) it received the earliest, greatest input and (2) the delayed synaptic input to the adjacent S-cells coincided with the action potential refractory period. A compartmental model of the S-cell and its inputs showed that a simple, intrinsic mechanism of inexcitability after each action potential may account for suppression of impulse initiations. Thus, a non-synaptic competition between neurons alters synaptic integration in the chain. In one mode, inputs to different sites sum independently, whereas in another, synaptic input to a single site precisely specifies the overall pattern of activity."
27. Composite spiking network/neural field model of Parkinsons (Kerr et al 2013)
This code implements a composite model of Parkinson's disease (PD). The composite model consists of a leaky integrate-and-fire spiking neuronal network model being driven by output from a neural field model (instead of the more usual white noise drive). Three different sets of parameters were used for the field model: one with basal ganglia parameters based on data from healthy individuals, one based on data from individuals with PD, and one purely thalamocortical model. The aim of this model is to explore how the different dynamical patterns in each each of these field models affects the activity in the network model.
28. Computational analysis of NN activity and spatial reach of sharp wave-ripples (Canakci et al 2017)
Network oscillations of different frequencies, durations and amplitudes are hypothesized to coordinate information processing and transfer across brain areas. Among these oscillations, hippocampal sharp wave-ripple complexes (SPW-Rs) are one of the most prominent. SPW-Rs occurring in the hippocampus are suggested to play essential roles in memory consolidation as well as information transfer to the neocortex. To-date, most of the knowledge about SPW-Rs comes from experimental studies averaging responses from neuronal populations monitored by conventional microelectrodes. In this work, we investigate spatiotemporal characteristics of SPW-Rs and how microelectrode size and distance influence SPW-R recordings using a biophysical model of hippocampus. We also explore contributions from neuronal spikes and synaptic potentials to SPW-Rs based on two different types of network activity. Our study suggests that neuronal spikes from pyramidal cells contribute significantly to ripples while high amplitude sharp waves mainly arise from synaptic activity. Our simulations on spatial reach of SPW-Rs show that the amplitudes of sharp waves and ripples exhibit a steep decrease with distance from the network and this effect is more prominent for smaller area electrodes. Furthermore, the amplitude of the signal decreases strongly with increasing electrode surface area as a result of averaging. The relative decrease is more pronounced when the recording electrode is closer to the source of the activity. Through simulations of field potentials across a high-density microelectrode array, we demonstrate the importance of finding the ideal spatial resolution for capturing SPW-Rs with great sensitivity. Our work provides insights on contributions from spikes and synaptic potentials to SPW-Rs and describes the effect of measurement configuration on LFPs to guide experimental studies towards improved SPW-R recordings.
29. Computational aspects of feedback in neural circuits (Maass et al 2006)
It had previously been shown that generic cortical microcircuit models can perform complex real-time computations on continuous input streams, provided that these computations can be carried out with a rapidly fading memory. We investigate ... the computational capability of such circuits in the more realistic case where not only readout neurons, but in addition a few neurons within the circuit have been trained for specific tasks. This is essentially equivalent to the case where the output of trained readout neurons is fed back into the circuit. We show that this new model overcomes the limitation of a rapidly fading memory. In fact, we prove that in the idealized case without noise it can carry out any conceiv- able digital or analog computation on time-varying inputs. But even with noise the resulting computational model can perform a large class of biologically relevant real-time computations that require a non-fading memory. ... In particular we show that ... generic cortical microcircuits with feedback provide a new model for working memory that is consistent with a large set of biological constraints. See paper for more and details.
30. Computational Model of a Central Pattern Generator (Cataldo et al 2006)
The buccal ganglia of Aplysia contain a central pattern generator (CPG) that mediates rhythmic movements of the foregut during feeding. This CPG is a multifunctional circuit and generates at least two types of buccal motor patterns (BMPs), one that mediates ingestion (iBMP) and another that mediates rejection (rBMP). The present study used a computational approach to examine the ways in which an ensemble of identified cells and synaptic connections function as a CPG. Hodgkin-Huxley-type models were developed that mimicked the biophysical properties of these cells and synaptic connections. The results suggest that the currently identified ensemble of cells is inadequate to produce rhythmic neural activity and that several key elements of the CPG remain to be identified.
31. Computer model of clonazepam`s effect in thalamic slice (Lytton 1997)
Demonstration of the effect of a minor pharmacological synaptic change at the network level. Clonazepam, a benzodiazepine, enhances inhibition but is paradoxically useful for certain types of seizures. This simulation shows how inhibition of inhibitory cells (the RE cells) produces this counter-intuitive effect.
32. Conductance-based model of Layer-4 in the barrel cortex (Argaman et Golomb 2017)
Layer 4 in the mouse barrel cortex includes hundreds of inhibitory PV neurons and thousands of excitatory neurons. Despite this fact, its dynamical state is similar to a balanced state of large neuronal circuits.
33. Convergence regulates synchronization-dependent AP transfer in feedforward NNs (Sailamul et al 2017)
We study how synchronization-dependent spike transfer can be affected by the structure of convergent feedforward wiring. We implemented computer simulations of model neural networks: a source and a target layer connected with different types of convergent wiring rules. In the Gaussian-Gaussian (GG) model, both the connection probability and the strength are given as Gaussian distribution as a function of spatial distance. In the Uniform-Constant (UC) and Uniform-Exponential (UE) models, the connection probability density is a uniform constant within a certain range, but the connection strength is set as a constant value or an exponentially decaying function, respectively. Then we examined how the spike transfer function is modulated under these conditions, while static or synchronized input patterns were introduced to simulate different levels of feedforward spike synchronization. We observed that the synchronization-dependent modulation of the transfer function appeared noticeably different for each convergence condition. The modulation of the spike transfer function was largest in the UC model, and smallest in the UE model. Our analysis showed that this difference was induced by the different spike weight distributions that was generated from convergent synapses in each model. Our results suggest that the structure of the feedforward convergence is a crucial factor for correlation-dependent spike control, thus must be considered important to understand the mechanism of information transfer in the brain.
34. Cortex-Basal Ganglia-Thalamus network model (Kumaravelu et al. 2016)
" ... We developed a biophysical network model comprising of the closed loop cortical-basal ganglia-thalamus circuit representing the healthy and parkinsonian rat brain. The network properties of the model were validated by comparing responses evoked in basal ganglia (BG) nuclei by cortical (CTX) stimulation to published experimental results. A key emergent property of the model was generation of low-frequency network oscillations. Consistent with their putative pathological role, low-frequency oscillations in model BG neurons were exaggerated in the parkinsonian state compared to the healthy condition. ..."
35. Current Dipole in Laminar Neocortex (Lee et al. 2013)
Laminar neocortical model in NEURON/Python, adapted from Jones et al 2009.
36. Deconstruction of cortical evoked potentials generated by subthalamic DBS (Kumaravelu et al 2018)
"... High frequency deep brain stimulation (DBS) of the subthalamic nucleus (STN) suppresses parkinsonian motor symptoms and modulates cortical activity. ... Cortical evoked potentials (cEP) generated by STN DBS reflect the response of cortex to subcortical stimulation, and the goal was to determine the neural origin of cEP using a two-step approach. First, we recorded cEP over ipsilateral primary motor cortex during different frequencies of STN DBS in awake healthy and unilateral 6-OHDA lesioned parkinsonian rats. Second, we used a biophysically-based model of the thalamocortical network to deconstruct the neural origin of the cEP. The in vivo cEP included short (R1), intermediate (R2) and long-latency (R3) responses. Model-based cortical responses to simulated STN DBS matched remarkably well the in vivo responses. R1 was generated by antidromic activation of layer 5 pyramidal neurons, while recurrent activation of layer 5 pyramidal neurons via excitatory axon collaterals reproduced R2. R3 was generated by polysynaptic activation of layer 2/3 pyramidal neurons via the cortico-thalamic-cortical pathway. Antidromic activation of the hyperdirect pathway and subsequent intracortical and cortico-thalamo-cortical synaptic interactions were sufficient to generate cEP by STN DBS, and orthodromic activation through basal ganglia-thalamus-cortex pathways was not required. These results demonstrate the utility of cEP to determine the neural elements activated by STN DBS that might modulate cortical activity and contribute to the suppression of parkinsonian symptoms."
37. Dentate Gyrus Feed-forward inhibition (Ferrante et al. 2009)
In this paper, the model was used to show how that FFI can change a steeply sigmoidal input-output (I/O) curve into a double-sigmoid typical of buffer systems.
38. Dentate gyrus network model (Santhakumar et al 2005)
Mossy cell loss and mossy fiber sprouting are two characteristic consequences of repeated seizures and head trauma. However, their precise contributions to the hyperexcitable state are not well understood. Because it is difficult, and frequently impossible, to independently examine using experimental techniques whether it is the loss of mossy cells or the sprouting of mossy fibers that leads to dentate hyperexcitability, we built a biophysically realistic and anatomically representative computational model of the dentate gyrus to examine this question. The 527-cell model, containing granule, mossy, basket, and hilar cells with axonal projections to the perforant-path termination zone, showed that even weak mossy fiber sprouting (10-15% of the strong sprouting observed in the pilocarpine model of epilepsy) resulted in the spread of seizure-like activity to the adjacent model hippocampal laminae after focal stimulation of the perforant path. See reference for more and details.
39. Dentate gyrus network model (Tejada et al 2014)
" ... Here we adapted an existing computational model of the dentate gyrus (J Neurophysiol 93: 437-453, 2005) by replacing the reduced granule cell models with morphologically detailed models coming from (3D) reconstructions of mature cells. ... Different fractions of the mature granule cell models were replaced by morphologically reconstructed models of newborn dentate granule cells from animals with PILO-induced Status Epilepticus, which have apical dendritic alterations and spine loss, and control animals, which do not have these alterations. This complex arrangement of cells and processes allowed us to study the combined effect of mossy fiber sprouting, altered apical dendritic tree and dendritic spine loss in newborn granule cells on the excitability of the dentate gyrus model. Our simulations suggest that alterations in the apical dendritic tree and dendritic spine loss in newborn granule cells have opposing effects on the excitability of the dentate gyrus after Status Epilepticus. Apical dendritic alterations potentiate the increase of excitability provoked by mossy fiber sprouting while spine loss curtails this increase. "
40. Dentate gyrus network model pattern separation and granule cell scaling in epilepsy (Yim et al 2015)
The dentate gyrus (DG) is thought to enable efficient hippocampal memory acquisition via pattern separation. With patterns defined as spatiotemporally distributed action potential sequences, the principal DG output neurons (granule cells, GCs), presumably sparsen and separate similar input patterns from the perforant path (PP). In electrophysiological experiments, we have demonstrated that during temporal lobe epilepsy (TLE), GCs downscale their excitability by transcriptional upregulation of ‘leak’ channels. Here we studied whether this cell type-specific intrinsic plasticity is in a position to homeostatically adjust DG network function. We modified an established conductance-based computer model of the DG network such that it realizes a spatiotemporal pattern separation task, and quantified its performance with and without the experimentally constrained leaky GC phenotype. ...
41. Duration-tuned neurons from the inferior colliculus of vertebrates (Aubie et al. 2012)
These models reproduce the responses of duration-tuned neurons in the auditory midbrain of the big brown bat, the rat, the mouse and the frog (Aubie et al. 2012). They are written in the Python interface to NEURON and a subset of the figures from Aubie et al. (2012) are pre-set in (raw data is generated and a separate graphing program must be used to visualize the results).
42. Dynamic cortical interlaminar interactions (Carracedo et al. 2013)
"... Here we demonstrate the mechanism underlying a purely neocortical delta rhythm generator and show a remarkable laminar, cell subtype and local subcircuit delineation between delta and nested theta rhythms. We show that spike timing during delta-nested theta rhythms controls an iterative, reciprocal interaction between deep and superficial cortical layers resembling the unsupervised learning processes proposed for laminar neural networks by Hinton and colleagues ... and mimicking the alternating cortical dynamics of sensory and memory processing during wakefulness."
43. Dynamic dopamine modulation in the basal ganglia: Learning in Parkinson (Frank et al 2004,2005)
See README file for all info on how to run models under different tasks and simulated Parkinson's and medication conditions.
44. Effects of increasing CREB on storage and recall processes in a CA1 network (Bianchi et al. 2014)
Several recent results suggest that boosting the CREB pathway improves hippocampal-dependent memory in healthy rodents and restores this type of memory in an AD mouse model. However, not much is known about how CREB-dependent neuronal alterations in synaptic strength, excitability and LTP can boost memory formation in the complex architecture of a neuronal network. Using a model of a CA1 microcircuit, we investigate whether hippocampal CA1 pyramidal neuron properties altered by increasing CREB activity may contribute to improve memory storage and recall. With a set of patterns presented to a network, we find that the pattern recall quality under AD-like conditions is significantly better when boosting CREB function with respect to control. The results are robust and consistent upon increasing the synaptic damage expected by AD progression, supporting the idea that the use of CREB-based therapies could provide a new approach to treat AD.
45. Electrically-coupled Retzius neurons (Vazquez et al. 2009)
"Dendritic electrical coupling increases the number of effective synaptic inputs onto neurons by allowing the direct spread of synaptic potentials from one neuron to another. Here we studied the summation of excitatory postsynaptic potentials (EPSPs) produced locally and arriving from the coupled neuron (transjunctional) in pairs of electrically-coupled Retzius neurons of the leech. We combined paired recordings of EPSPs, the production of artificial EPSPs (APSPs) in neuron pairs with different coupling coefficients and simulations of EPSPs produced in the coupled dendrites. ..."
46. ELL pyramidal neuron (Simmonds and Chacron 2014)
network model of ELL pyramidal neurons receiving both feedforward and feedback inputs
47. Engaging distinct oscillatory neocortical circuits (Vierling-Claassen et al. 2010)
"Selective optogenetic drive of fast-spiking (FS) interneurons (INs) leads to enhanced local field potential (LFP) power across the traditional “gamma” frequency band (20–80 Hz; Cardin et al., 2009). In contrast, drive to regular-spiking (RS) pyramidal cells enhances power at lower frequencies, with a peak at 8 Hz. The first result is consistent with previous computational studies emphasizing the role of FS and the time constant of GABAA synaptic inhibition in gamma rhythmicity. However, the same theoretical models do not typically predict low-frequency LFP enhancement with RS drive. To develop hypotheses as to how the same network can support these contrasting behaviors, we constructed a biophysically principled network model of primary somatosensory neocortex containing FS, RS, and low-threshold spiking (LTS) INs. ..."
48. Epilepsy may be caused by very small functional changes in ion channels (Thomas et al. 2009)
We used a previously published model of the dentate gyrus with varying degrees of mossy fibre sprouting.We preformed a sensitivity analysis where we systematically varied individual properties of ion channels. The results predict that genetic variations in the properties of sodium channels are likely to have the biggest impact on network excitability. Furthermore, these changes may be as small as 1mV, which is currently undetectable using standard experimental practices.
49. Escape response latency in the Giant Fiber System of Drosophila melanogastor (Augustin et al 2019)
"The Giant Fiber System (GFS) is a multi-component neuronal pathway mediating rapid escape response in the adult fruit-fly Drosophila melanogaster, usually in the face of a threatening visual stimulus. Two branches of the circuit promote the response by stimulating an escape jump followed by flight initiation. Our recent work demonstrated an age-associated decline in the speed of signal propagation through the circuit, measured as the stimulus-to-muscle depolarization response latency. The decline is likely due to the diminishing number of inter-neuronal gap junctions in the GFS of ageing flies. In this work, we presented a realistic conductance-based, computational model of the GFS that recapitulates our experimental results and identifies some of the critical anatomical and physiological components governing the circuit's response latency. According to our model, anatomical properties of the GFS neurons have a stronger impact on the transmission than neuronal membrane conductance densities. The model provides testable predictions for the effect of experimental interventions on the circuit's performance in young and ageing flies."
50. Failure of Deep Brain Stimulation in a basal ganglia neuronal network model (Dovzhenok et al. 2013)
"… Recently, a lot of interest has been devoted to desynchronizing delayed feedback deep brain stimulation (DBS). ... This study explores the action of delayed feedback stimulation on partially synchronized oscillatory dynamics, similar to what one observes experimentally in parkinsonian patients. …" Implemented by Andrey Dovzhenok, to whom questions should be addressed.
51. Fast oscillations in inhibitory networks (Maex, De Schutter 2003)
We observed a new phenomenon of resonant synchronization in computer-simulated networks of inhibitory neurons in which the synaptic current has a delayed onset, reflecting finite spike propagation and synaptic transmission times. At the resonant level of network excitation, all neurons fire synchronously and rhythmically with a period approximately four times the mean delay of the onset of the inhibitory synaptic current. ... By varying the axonal delay of the inhibitory connections, networks with a realistic synaptic kinetics can be tuned to frequencies from 40 to >200 Hz. ... We conclude that the delay of the synaptic current is the primary parameter controlling the oscillation frequency of inhibitory networks and propose that delay-induced synchronization is a mechanism for fast brain rhythms that depend on intact inhibitory synaptic transmission.
52. Feedforward heteroassociative network with HH dynamics (Lytton 1998)
Using the original McCulloch-Pitts notion of simple on and off spike coding in lieu of rate coding, an Anderson-Kohonen artificial neural network (ANN) associative memory model was ported to a neuronal network with Hodgkin-Huxley dynamics.
53. Gamma and theta rythms in biophysical models of hippocampus circuits (Kopell et al. 2011)
" ... the main rhythms displayed by the hippocampus, the gamma (30–90 Hz) and theta (4–12 Hz) rhythms. We concentrate on modeling in vitro experiments, but with an eye toward possible in vivo implications. ... We use simpler biophysical models; all cells have a single compartment only, and the interneurons are restricted to two types: fast-spiking (FS) basket cells and oriens lacunosum-moleculare (O-LM) cells. ... , we aim not so much at reproducing dynamics in great detail, but at clarifying the essential mechanisms underlying the production of the rhythms and their interactions (Kopell, 2005). ..."
54. Gamma genesis in the basolateral amygdala (Feng et al 2019)
Using in vitro and in vivo data we develop the first large-scale biophysically and anatomically realistic model of the basolateral amygdala nucleus (BL), which reproduces the dynamics of the in vivo local field potential (LFP). Significantly, it predicts that BL intrinsically generates the transient gamma oscillations observed in vivo. The model permitted exploration of the poorly understood synaptic mechanisms underlying gamma genesis in BL, and the model's ability to compute LFPs at arbitrary numbers of recording sites provided insights into the characteristics of the spatial properties of gamma bursts. Furthermore, we show how gamma synchronizes principal cells to overcome their low firing rates while simultaneously promoting competition, potentially impacting their afferent selectivity and efferent drive, and thus emotional behavior.
55. Gamma oscillations in hippocampal interneuron networks (Bartos et al 2002)
To examine whether an interneuron network with fast inhibitory synapses can act as a gamma frequency oscillator, we developed an interneuron network model based on experimentally determined properties. In comparison to previous interneuron network models, our model was able to generate oscillatory activity with higher coherence over a broad range of frequencies (20-110 Hz). In this model, high coherence and flexibility in frequency control emerge from the combination of synaptic properties, network structure, and electrical coupling.
56. Gamma oscillations in hippocampal interneuron networks (Wang, Buzsaki 1996)
The authors investigated the hypothesis that 20-80Hz neuronal (gamma) oscillations can emerge in sparsely connected network models of GABAergic fast-spiking interneurons. They explore model NN synchronization and compare their results to anatomical and electrophysiological data from hippocampal fast spiking interneurons.
57. Gamma-beta alternation in the olfactory bulb (David, Fourcaud-Trocmé et al., 2015)
This model, a simplified olfactory bulb network with mitral and granule cells, proposes a framework for two regimes of oscillation in the olfactory bulb: 1 - a weak inhibition regime (with no granule spike) where the network oscillates in the gamma (40-90Hz) band 2 - a strong inhibition regime (with granule spikes) where the network oscillates in the beta (15-30Hz) band. Slow modulations of sensory and centrifugal inputs, phase shifted by a quarter of cycle, possibly combined with short term depression of the mitral to granule AMPA synapse, allows the network to alternate between the two regimes as observed in anesthetized animals.
58. Gap junction coupled network of striatal fast spiking interneurons (Hjorth et al. 2009)
Gap junctions between striatal FS neurons has very weak ability to synchronise spiking. Input uncorrelated between neighbouring neurons is shunted, while correlated input is not.
59. Gap-junction coupled network activity depends on coupled dendrites diameter (Gansert et al. 2007)
"... We have previously shown that the amplitude of electrical signals propagating across gap-junctionally coupled passive cables is maximized at a unique diameter. This suggests that threshold-dependent signals may propagate through gap junctions for a finite range of diameters around this optimal value. Here we examine the diameter dependence of action potential propagation across model networks of dendro-dendritically coupled neurons. The neurons in these models have passive soma and dendrites and an action potential-generating axon. We show that propagation of action potentials across gap junctions occurs only over a finite range of dendritic diameters and that propagation delay depends on this diameter. ...". See paper for more and details.
60. Gating of steering signals through phasic modulation of reticulospinal neurons (Kozlov et al. 2014)
" ... We use the lamprey as a model for investigating the role of this phasic modulation of the reticulospinal activity, because the brainstem–spinal cord networks are known down to the cellular level in this phylogenetically oldest extant vertebrate. We describe how the phasic modulation of reticulospinal activity from the spinal CPG ensures reliable steering/turning commands without the need for a very precise timing of on- or offset, by using a biophysically detailed large-scale (19,600 model neurons and 646,800 synapses) computational model of the lamprey brainstem–spinal cord network. To verify that the simulated neural network can control body movements, including turning, the spinal activity is fed to a mechanical model of lamprey swimming. ..."
61. Grid cell oscillatory interference with noisy network oscillators (Zilli and Hasselmo 2010)
To examine whether an oscillatory interference model of grid cell activity could work if the oscillators were noisy neurons, we implemented these simulations. Here the oscillators are networks (either synaptically- or gap-junction--coupled) of one or more noisy neurons (either Izhikevich's simple model or a Hodgkin-Huxley--type biophysical model) which drive a postsynaptic cell (which may be integrate-and-fire, resonate-and-fire, or the simple model) which should fire spatially as a grid cell if the simulation is successful.
62. H-currents effect on the fluctuation of gamma/beta oscillations (Avella-Gonzalez et al., 2015)
This model was designed to study the impact of H-currents on the dynamics of cortical oscillations, and in paticular on the occurrence of high and low amplitude episodes (HAE, LAE) in network oscillations. The H-current is a slow, hyperpolarization-activated, depolarizing current that contributes to neuronal resonance and membrane potential. We characterized amplitude fluctuations in network oscillations by measuring the average durations of HAEs and LAEs, and explored how these were modulated by trains of external spikes, both in the presence and absence of H-channels. We looked at HAE duration, the frequency and power of network oscillations, and the effect of H-channels on the temporal voltage profile in single cells. We found that H-currents increased the oscillation frequency and, in combination with external spikes, representing input from areas outside the network, strongly decreased the synchrony of firing. As a consequence, the oscillation power and the duration of episodes during which the network exhibited high-amplitude oscillations were greatly reduced in the presence of H-channels.
63. Half-center oscillator database of leech heart interneuron model (Doloc-Mihu & Calabrese 2011)
We have created a database (HCO-db) of instances of a half-center oscillator computational model [Hill et al., 2001] for analyzing how neuronal parameters influence network activity. We systematically explored the parameter space of about 10.4 million simulated HCO instances and corresponding isolated neuron model simulations obtained by varying a set of selected parameters (maximal conductance of intrinsic and synaptic currents) in all combinations using a brute-force approach. We classified these HCO instances by their activity characteristics into identifiable groups. We built an efficient relational database table (HCO-db) with the resulting instances characteristics.
64. High frequency oscillations in a hippocampal computational model (Stacey et al. 2009)
"... Using a physiological computer model of hippocampus, we investigate random synaptic activity (noise) as a potential initiator of HFOs (high-frequency oscillations). We explore parameters necessary to produce these oscillations and quantify the response using the tools of stochastic resonance (SR) and coherence resonance (CR). ... Our results show that, under normal coupling conditions, synaptic noise was able to produce gamma (30–100 Hz) frequency oscillations. Synaptic noise generated HFOs in the ripple range (100–200 Hz) when the network had parameters similar to pathological findings in epilepsy: increased gap junctions or recurrent synaptic connections, loss of inhibitory interneurons such as basket cells, and increased synaptic noise. ... We propose that increased synaptic noise and physiological coupling mechanisms are sufficient to generate gamma oscillations and that pathologic changes in noise and coupling similar to those in epilepsy can produce abnormal ripples."
65. High frequency stimulation of the Subthalamic Nucleus (Rubin and Terman 2004)
" ... Using a computational model, this paper considers the hypothesis that DBS works by replacing pathologically rhythmic basal ganglia output with tonic, high frequency firing. In our simulations of parkinsonian conditions, rhythmic inhibition from GPi to the thalamus compromises the ability of thalamocortical relay (TC) cells to respond to depolarizing inputs, such as sensorimotor signals. High frequency stimulation of STN regularizes GPi firing, and this restores TC responsiveness, despite the increased frequency and amplitude of GPi inhibition to thalamus that result. We provide a mathematical phase plane analysis of the mechanisms that determine TC relay capabilities in normal, parkinsonian, and DBS states in a reduced model. This analysis highlights the differences in deinactivation of the low-threshold calcium T -current that we observe in TC cells in these different conditions. ..."
66. Hippocampal basket cell gap junction network dynamics (Saraga et al. 2006)
2 cell network of hippocampal basket cells connected by gap junctions. Paper explores how distal gap junctions and active dendrites can tune network dynamics.
67. Hippocampal CA1 NN with spontaneous theta, gamma: full scale & network clamp (Bezaire et al 2016)
This model is a full-scale, biologically constrained rodent hippocampal CA1 network model that includes 9 cells types (pyramidal cells and 8 interneurons) with realistic proportions of each and realistic connectivity between the cells. In addition, the model receives realistic numbers of afferents from artificial cells representing hippocampal CA3 and entorhinal cortical layer III. The model is fully scaleable and parallelized so that it can be run at small scale on a personal computer or large scale on a supercomputer. The model network exhibits spontaneous theta and gamma rhythms without any rhythmic input. The model network can be perturbed in a variety of ways to better study the mechanisms of CA1 network dynamics. Also see online code at and further information at
68. Hippocampal CA3 network and circadian regulation (Stanley et al. 2013)
This model produces the hippocampal CA3 neural network model used in the paper below. It has two modes of operation, a default mode and a circadian mode. In the circadian mode, parameters are swept through a range of values. This model can be quite easily adapted to produce theta and gamma oscillations, as certain parameter sweeps will reveal (see Figures). BASH scripts interact with GENESIS 2.3 to implement parameter sweeps. The model contains four cell types derived from prior papers. CA3 pyramidal are derived from Traub et al (1991); Basket, stratum oriens (O-LM), and Medial Septal GABAergic (MSG) interneurons are taken from Hajos et al (2004).
69. Hippocampus temporo-septal engram shift model (Lytton 1999)
Temporo-septal engram shift model of hippocampal memory. The model posits that memories gradually move along the hippocampus from a temporal encoding site to ever more septal sites from which they are recalled. We propose that the sense of time is encoded by the location of the engram along the temporo-septal axis.
70. Homosynaptic plasticity in the tail withdrawal circuit (TWC) of Aplysia (Baxter and Byrne 2006)
The tail-withdrawal circuit of Aplysia provides a useful model system for investigating synaptic dynamics. Sensory neurons within the circuit manifest several forms of synaptic plasticity. Here, we developed a model of the circuit and investigated the ways in which depression (DEP) and potentiation (POT) contributed to information processing. DEP limited the amount of motor neuron activity that could be elicited by the monosynaptic pathway alone. POT within the monosynaptic pathway did not compensate for DEP. There was, however, a synergistic interaction between POT and the polysynaptic pathway. This synergism extended the dynamic range of the network, and the interplay between DEP and POT made the circuit respond preferentially to long-duration, low-frequency inputs.
71. Hopfield and Brody model (Hopfield, Brody 2000)
NEURON implementation of the Hopfield and Brody model from the papers: JJ Hopfield and CD Brody (2000) JJ Hopfield and CD Brody (2001). Instructions are provided in the below readme.txt file.
72. Ih tunes oscillations in an In Silico CA3 model (Neymotin et al. 2013)
" ... We investigated oscillatory control using a multiscale computer model of hippocampal CA3, where each cell class (pyramidal, basket, and oriens-lacunosum moleculare cells), contained type-appropriate isoforms of Ih. Our model demonstrated that modulation of pyramidal and basket Ih allows tuning theta and gamma oscillation frequency and amplitude. Pyramidal Ih also controlled cross-frequency coupling (CFC) and allowed shifting gamma generation towards particular phases of the theta cycle, effected via Ih’s ability to set pyramidal excitability. ..."
73. Information-processing in lamina-specific cortical microcircuits (Haeusler and Maass 2006)
A major challenge for computational neuroscience is to understand the computational function of lamina-specific synaptic connection patterns in stereotypical cortical microcircuits.We approach this problem by studying ... the dynamical system defined by more realistic cortical microcircuit models as a whole and by investigating the influence that its laminar structure has on the transmission and fusion of information within this dynamical system. The circuit models that we examine consist of Hodgkin--Huxley neurons with dynamic synapses... We investigate to what extent this cortical microcircuit template supports the accumulation and fusion of information contained in generic spike inputs into layer 4 and layers 2/3 and how well it makes this information accessible to projection neurons in layers 2/3 and layer 5. ... We conclude that computer simulations of detailed lamina-specific cortical microcircuit models provide new insight into computational consequences of anatomical and physiological data. See paper for more and details.
74. Inhibition and glial-K+ interaction leads to diverse seizure transition modes (Ho & Truccolo 2016)
"How focal seizures initiate and evolve in human neocortex remains a fundamental problem in neuroscience. Here, we use biophysical neuronal network models of neocortical patches to study how the interaction between inhibition and extracellular potassium ([K+]o) dynamics may contribute to different types of focal seizures. Three main types of propagated focal seizures observed in recent intracortical microelectrode recordings in humans were modelled ..."
75. Investigation of different targets in deep brain stimulation for Parkinson`s (Pirini et al. 2009)
"We investigated by a computational model of the basal ganglia the different network effects of deep brain stimulation (DBS) for Parkinson’s disease (PD) in different target sites in the subthalamic nucleus (STN), the globus pallidus pars interna (GPi), and the globus pallidus pars externa (GPe). A cellular-based model of the basal ganglia system (BGS), based on the model proposed by Rubin and Terman (J Comput Neurosci 16:211–235, 2004), was developed. ... Our results suggest that DBS in the STN could functionally restore the TC relay activity, while DBS in the GPe and in the GPi could functionally over-activate and inhibit it, respectively. Our results are consistent with the experimental and the clinical evidences on the network effects of DBS."
76. Irregular spiking in NMDA-driven prefrontal cortex neurons (Durstewitz and Gabriel 2006)
Slow N-Methyl-D-aspartic acid (NMDA) synaptic currents are assumed to strongly contribute to the persistently elevated firing rates observed in prefrontal cortex (PFC) during working memory. During persistent activity, spiking of many neurons is highly irregular. ... The highest interspike-interval (ISI) variability occurred in a transition regime where the subthreshold membrane potential distribution shifts from mono- to bimodality, ... Predictability within irregular ISI series was significantly higher than expected from a noise-driven linear process, indicating that it might best be described through complex (potentially chaotic) nonlinear deterministic processes. Accordingly, the phenomena observed in vitro could be reproduced in purely deterministic biophysical model neurons. High spiking irregularity in these models emerged within a chaotic, close-to-bifurcation regime characterized by a shift of the membrane potential distribution from mono- to bimodality and by similar ISI return maps as observed in vitro. ... NMDA-induced irregular dynamics may have important implications for computational processes during working memory and neural coding.
77. Ketamine disrupts theta modulation of gamma in a computer model of hippocampus (Neymotin et al 2011)
"Abnormalities in oscillations have been suggested to play a role in schizophrenia. We studied theta-modulated gamma oscillations in a computer model of hippocampal CA3 in vivo with and without simulated application of ketamine, an NMDA receptor antagonist and psychotomimetic. Networks of 1200 multi-compartment neurons (pyramidal, basket and oriens-lacunosum moleculare, OLM, cells) generated theta and gamma oscillations from intrinsic network dynamics: basket cells primarily generated gamma and amplified theta, while OLM cells strongly contributed to theta. ..."
78. KInNeSS : a modular framework for computational neuroscience (Versace et al. 2008)
The xml files provided here implement a network of excitatory and inhibitory spiking neurons, governed by either Hodgkin-Huxley or quadratic integrate-and-fire dynamical equations. The code is used to demonstrate the capabilities of the KInNeSS software package for simulation of networks of spiking neurons. The simulation protocol used here is meant to facilitate the comparison of KInNeSS with other simulators reviewed in <a href="">Brette et al. (2007)</a>. See the associated paper "Versace et al. (2008) KInNeSS : a modular framework for computational neuroscience." for an extensive description of KInNeSS .
79. Knox implementation of Destexhe 1998 spike and wave oscillation model (Knox et al 2018)
" ...The aim of this study was to use an established thalamocortical computer model to determine how T-type calcium channels work in concert with cortical excitability to contribute to pathogenesis and treatment response in CAE. METHODS: The model is comprised of cortical pyramidal, cortical inhibitory, thalamocortical relay, and thalamic reticular single-compartment neurons, implemented with Hodgkin-Huxley model ion channels and connected by AMPA, GABAA , and GABAB synapses. Network behavior was simulated for different combinations of T-type calcium channel conductance, inactivation time, steady state activation/inactivation shift, and cortical GABAA conductance. RESULTS: Decreasing cortical GABAA conductance and increasing T-type calcium channel conductance converted spindle to spike and wave oscillations; smaller changes were required if both were changed in concert. In contrast, left shift of steady state voltage activation/inactivation did not lead to spike and wave oscillations, whereas right shift reduced network propensity for oscillations of any type...."
80. L4 cortical barrel NN model receiving thalamic input during whisking or touch (Gutnisky et al. 2017)
Excitatory neurons in layer 4 (L4) in the barrel cortex respond relatively strongly to touch but not to whisker movement (Yu et al., Nat. Neurosci. 2016). The model explains the mechanism underlying this effect. The network is settled to filter out most stationary inputs. Brief touch input passes through because it takes time until feed-forward inhibition silences excitatory neurons receiving brief and strong thalamic excitation.
81. L5 PFC microcircuit used to study persistent activity (Papoutsi et al. 2014, 2013)
Using a heavily constrained biophysical model of a L5 PFC microcircuit we investigate the mechanisms that underlie persistent activity emergence (ON) and termination (OFF) and search for the minimum network size required for expressing these states within physiological regimes.
82. Large scale model of the olfactory bulb (Yu et al., 2013)
The readme file currently contains links to the results for all the 72 odors investigated in the paper, and the movie showing the network activity during learning of odor k3-3 (an aliphatic ketone).
83. Large-scale model of neocortical slice in vitro exhibiting persistent gamma (Tomsett et al. 2014)
This model contains 15 neuron populations (8 excitatory, 7 inhibitory) arranged into 4 cortical layers (layer 1 empty, layers 2/3 combined). It produces a persistent gamma oscillation driven by layer 2/3. It runs using the VERTEX simulator, which is written in Matlab and is available from
84. Lateral dendrodenditic inhibition in the Olfactory Bulb (David et al. 2008)
Mitral cells, the principal output neurons of the olfactory bulb, receive direct synaptic activation from primary sensory neurons. Shunting inhibitory inputs delivered by granule cell interneurons onto mitral cell lateral dendrites are believed to influence spike timing and underlie coordinated field potential oscillations. Lateral dendritic shunt conductances delayed spiking to a degree dependent on both their electrotonic distance and phase of onset. Recurrent inhibition significantly narrowed the distribution of mitral cell spike times, illustrating a tendency towards coordinated synchronous activity. This result suggests an essential role for early mechanisms of temporal coordination in olfaction. The model was adapted from Davison et al, 2003, but include additional noise mechanisms, long lateral dendrite, and specific synaptic point processes.
85. Leech Heart (HE) Motor Neuron conductances contributions to NN activity (Lamb & Calabrese 2013)
"... To explore the relationship between conductances, and in particular how they influence the activity of motor neurons in the well characterized leech heartbeat system, we developed a new multi-compartmental Hodgkin-Huxley style leech heart motor neuron model. To do so, we evolved a population of model instances, which differed in the density of specific conductances, capable of achieving specific output activity targets given an associated input pattern. ... We found that the strengths of many conductances, including those with differing dynamics, had strong partial correlations and that these relationships appeared to be linked by their influence on heart motor neuron activity. Conductances that had positive correlations opposed one another and had the opposite effects on activity metrics when perturbed whereas conductances that had negative correlations could compensate for one another and had similar effects on activity metrics. "
86. Leech heart interneuron network model (Hill et al 2001, 2002)
We have created a computational model of the timing network that paces the heartbeat of the medicinal leech, Hirudo medicinalis. In the intact nerve cord, segmental oscillators are mutually entrained to the same cycle period. Although experiments have shown that the segmental oscillators are coupled by inhibitory coordinating interneurons, the underlying mechanisms of intersegmental coordination have not yet been elucidated. To help understand this coordination, we have created a simple computational model with two variants: symmetric and asymmetric. See references for more details. Biologically realistic network models with two, six, and eight cells and a tutorial are available at the links to Calabrese's web site below.
87. Lobster STG pyloric network model with calcium sensor (Gunay & Prinz 2010) (Prinz et al. 2004)
This pyloric network model simulator is a C/C++ program that saves 384 different calcium sensor values that are candidates for activity sensors (Gunay and Prinz, 2010). The simulator was used to scan all of the 20 million pyloric network models that were previously collected in a database (Prinz et al, 2004).
88. Long time windows from theta modulated inhib. in entorhinal–hippo. loop (Cutsuridis & Poirazi 2015)
"A recent experimental study (Mizuseki et al., 2009) has shown that the temporal delays between population activities in successive entorhinal and hippocampal anatomical stages are longer (about 70–80 ms) than expected from axon conduction velocities and passive synaptic integration of feed-forward excitatory inputs. We investigate via computer simulations the mechanisms that give rise to such long temporal delays in the hippocampus structures. ... The model shows that the experimentally reported long temporal delays in the DG, CA3 and CA1 hippocampal regions are due to theta modulated somatic and axonic inhibition..."
89. Mechanisms of very fast oscillations in axon networks coupled by gap junctions (Munro, Borgers 2010)
Axons connected by gap junctions can produce very fast oscillations (VFOs, > 80 Hz) when stimulated randomly at a low rate. The models here explore the mechanisms of VFOs that can be seen in an axonal plexus, (Munro & Borgers, 2009): a large network model of an axonal plexus, small network models of axons connected by gap junctions, and an implementation of the model underlying figure 12 in Traub et al. (1999) . The large network model consists of 3,072 5-compartment axons connected in a random network. The 5-compartment axons are the 5 axonal compartments from the CA3 pyramidal cell model in Traub et al. (1994) with a fixed somatic voltage. The random network has the same parameters as the random network in Traub et al. (1999), and axons are stimulated randomly via a Poisson process with a rate of 2/s/axon. The small network models simulate waves propagating through small networks of axons connected by gap junctions to study how local connectivity affects the refractory period.
90. MEG of Somatosensory Neocortex (Jones et al. 2007)
"... To make a direct and principled connection between the SI (somatosensory primary neocortex magnetoencephalography) waveform and underlying neural dynamics, we developed a biophysically realistic computational SI model that contained excitatory and inhibitory neurons in supragranular and infragranular layers. ... our model provides a biophysically realistic solution to the MEG signal and can predict the electrophysiological correlates of human perception."
91. Microcircuits of L5 thick tufted pyramidal cells (Hay & Segev 2015)
"... We simulated detailed conductance-based models of TTCs (Layer 5 thick tufted pyramidal cells) forming recurrent microcircuits that were interconnected as found experimentally; the network was embedded in a realistic background synaptic activity. ... Our findings indicate that dendritic nonlinearities are pivotal in controlling the gain and the computational functions of TTCs microcircuits, which serve as a dominant output source for the neocortex. "
92. Model of arrhythmias in a cardiac cells network (Casaleggio et al. 2014)
" ... Here we explore the possible processes leading to the occasional onset and termination of the (usually) non-fatal arrhythmias widely observed in the heart. Using a computational model of a two-dimensional network of cardiac cells, we tested the hypothesis that an ischemia alters the properties of the gap junctions inside the ischemic area. ... In conclusion, our model strongly supports the hypothesis that non-fatal arrhythmias can develop from post-ischemic alteration of the electrical connectivity in a relatively small area of the cardiac cell network, and suggests experimentally testable predictions on their possible treatments."
93. Model of eupnea and sigh generation in respiratory network (Toporikova et al 2015)
Based on recent in vitro data obtained in the mouse embryo, we have built a computational model consisting of two compartments, interconnected through appropriate synapses. One compartment generates sighs and the other produces eupneic bursts. The model reproduces basic features of simultaneous sigh and eupnea generation (two types of bursts differing in terms of shape, amplitude, and frequency of occurrence) and mimics the effect of blocking glycinergic synapses
94. Model of long range transmission of gamma oscillation (Murray 2007)
"... A minimal mathematical model was developed for a preliminary study of long-range neural transmission of gamma oscillation from the CA3 to the entorhinal cortex via the CAI region of the hippocampus, a subset within a larger complex set of pathways. A module was created for each local population of neurons with common intrinsic properties and connectivity to simplify the connection process and make the model more flexible. Three modules were created using MATLAB Simulink® and tested to confirm that they transmit gamma through the system. The model also revealed that a portion of the signal from CAI to the entorhinal cortex may be lost in transmission under certain conditions."
95. Model of the cerebellar granular network (Sudhakar et al 2017)
"The granular layer, which mainly consists of granule and Golgi cells, is the first stage of the cerebellar cortex and processes spatiotemporal information transmitted by mossy fiber inputs with a wide variety of firing patterns. To study its dynamics at multiple time scales in response to inputs approximating real spatiotemporal patterns, we constructed a large-scale 3D network model of the granular layer. ..."
96. Modelling platform of the cochlear nucleus and other auditory circuits (Manis & Compagnola 2018)
"Models of the auditory brainstem have been an invaluable tool for testing hypotheses about auditory information processing and for highlighting the most important gaps in the experimental literature. Due to the complexity of the auditory brainstem, and indeed most brain circuits, the dynamic behavior of the system may be difficult to predict without a detailed, biologically realistic computational model. Despite the sensitivity of models to their exact construction and parameters, most prior models of the cochlear nucleus have incorporated only a small subset of the known biological properties. This confounds the interpretation of modelling results and also limits the potential future uses of these models, which require a large effort to develop. To address these issues, we have developed a general purpose, bio-physically detailed model of the cochlear nucleus for use both in testing hypotheses about cochlear nucleus function and also as an input to models of downstream auditory nuclei. The model implements conductance-based Hodgkin-Huxley representations of cells using a Python-based interface to the NEURON simulator. ..."
97. Modulation of septo-hippocampal theta activity by GABAA receptors (Hajos et al. 2004)
Theta frequency oscillation of the septo-hippocampal system has been considered as a prominent activity associated with cognitive function and affective processes. ... In the present experiments we applied a combination of computational and physiological techniques to explore the functional role of GABAA receptors in theta oscillation. ... In parallel to these experimental observations, a computational model has been constructed by implementing a septal GABA neuron model with a CA1 hippocampal model containing three types of neurons (including oriens and basket interneurons and pyramidal cells; latter modeled by multicompartmental techniques; for detailed model description with network parameters see online addendum: This connectivity made the network capable of simulating the responses of the septo-hippocampal circuitry to the modulation of GABAA transmission, and the presently described computational model proved suitable to reveal several aspects of pharmacological modulation of GABAA receptors. In addition, computational findings indicated different roles of distinctively located GABAA receptors in theta generation.
98. Multiplication by NMDA receptors in Direction Selective Ganglion cells (Poleg-Polsky & Diamond 2016)
The model demonstrates how signal amplification with NMDARs depends on the synaptic environment. When direction selectivity (DS) detection is mediated by DS inhibition, NMDARs multiply other synaptic conductances. In the case of DS tuned excitation, NMDARs contribute additively.
99. Multitarget pharmacology for Dystonia in M1 (Neymotin et al 2016)
" ... We developed a multiscale model of primary motor cortex, ranging from molecular, up to cellular, and network levels, containing 1715 compartmental model neurons with multiple ion channels and intracellular molecular dynamics. We wired the model based on electrophysiological data obtained from mouse motor cortex circuit mapping experiments. We used the model to reproduce patterns of heightened activity seen in dystonia by applying independent random variations in parameters to identify pathological parameter sets. ..."
100. Muscle spindle feedback circuit (Moraud et al, 2016)
Here, we developed a computational model of the muscle spindle feedback circuits of the rat ankle that predicts the interactions between Epidural Stimulation and spinal circuit dynamics during gait.
101. Na channel mutations in the dentate gyrus (Thomas et al. 2009)
These are source files to generate the data in Figure 6 from "Mossy fiber sprouting interacts with sodium channel mutations to increase dentate gyrus excitability" Thomas EA, Reid CA, Petrou S, Epilepsia (2009)
102. Neocort. pyramidal cells subthreshold somatic voltage controls spike propagation (Munro Kopell 2012)
There is suggestive evidence that pyramidal cell axons in neocortex may be coupled by gap junctions into an ``axonal plexus" capable of generating Very Fast Oscillations (VFOs) with frequencies exceeding 80 Hz. It is not obvious, however, how a pyramidal cell in such a network could control its output when action potentials are free to propagate from the axons of other pyramidal cells into its own axon. We address this problem by means of simulations based on 3D reconstructions of pyramidal cells from rat somatosensory cortex. We show that somatic depolarization enables propagation via gap junctions into the initial segment and main axon, while somatic hyperpolarization disables it. We show further that somatic voltage cannot effectively control action potential propagation through gap junctions on minor collaterals; action potentials may therefore propagate freely from such collaterals regardless of somatic voltage. In previous work, VFOs are all but abolished during the hyperpolarization phase of slow-oscillations induced by anesthesia in vivo. This finding constrains the density of gap junctions on collaterals in our model and suggests that axonal sprouting due to cortical lesions may result in abnormally high gap junction density on collaterals, leading in turn to excessive VFO activity and hence to epilepsy via kindling.
103. Network bursts in cultured NN result from different adaptive mechanisms (Masquelier & Deco 2013)
It is now well established that cultured neuron networks are spontaneously active, and tend to synchronize. Synchronous events typically involve the whole network, and have thus been termed “network spikes” (NS). Using experimental recordings and numerical simulations, we show here that the inter-NS interval statistics are complex, and allow inferring the neural mechanisms at work, in particular the adaptive ones, and estimating a number of parameters to which we cannot access experimentally.
104. Network model of the granular layer of the cerebellar cortex (Maex, De Schutter 1998)
We computed the steady-state activity of a large-scale model of the granular layer of the rat cerebellum. Within a few tens of milliseconds after the start of random mossy fiber input, the populations of Golgi and granule cells became entrained in a single synchronous oscillation, the basic frequency of which ranged from 10 to 40 Hz depending on the average rate of firing in the mossy fiber population. ... The synchronous, rhythmic firing pattern was robust over a broad range of biologically realistic parameter values and to parameter randomization. Three conditions, however, made the oscillations more transient and could desynchronize the entire network in the end: a very low mossy fiber activity, a very dominant excitation of Golgi cells through mossy fiber synapses (rather than through parallel fiber synapses), and a tonic activation of granule cell GABAA receptors (with an almost complete absence of synaptically induced inhibitory postsynaptic currents). The model predicts that, under conditions of strong mossy fiber input to the cerebellum, Golgi cells do not only control the strength of parallel fiber activity but also the timing of the individual spikes. Provided that their parallel fiber synapses constitute an important source of excitation, Golgi cells fire rhythmically and synchronized with granule cells over large distances along the parallel fiber axis. See paper for more and details.
105. Network model with neocortical architecture (Anderson et al 2007,2012; Azhar et al 2012)
Architecturally realistic neocortical model using seven classes of excitatory and inhibitory single compartment Hodgkin-Huxley cells. This is an addendum to ModelDB Accession # 98902, Studies of stimulus parameters for seizure disruption (Anderson et al. 2007). Wiring is adapted from the minicolumn hypothesis and incorporates visual and neocortical wiring data. Simulation demonstrates spontaneous bursting onset and cessation. This activity can be induced by random fluctuations in the surrounding background input.
106. Network recruitment to coherent oscillations in a hippocampal model (Stacey et al. 2011)
"... Here we demonstrate, via a detailed computational model, a mechanism whereby physiological noise and coupling initiate oscillations and then recruit neighboring tissue, in a manner well described by a combination of Stochastic Resonance and Coherence Resonance. We develop a novel statistical method to quantify recruitment using several measures of network synchrony. This measurement demonstrates that oscillations spread via preexisting network connections such as interneuronal connections, recurrent synapses, and gap junctions, provided that neighboring cells also receive sufficient inputs in the form of random synaptic noise. ..."
107. NMDAR & GABAB/KIR Give Bistable Dendrites: Working Memory & Sequence Readout (Sanders et al., 2013)
" ...Here, we show that the voltage dependence of the inwardly rectifying potassium (KIR) conductance activated by GABA(B) receptors adds substantial robustness to network simulations of bistability and the persistent firing that it underlies. The hyperpolarized state is robust because, at hyperpolarized potentials, the GABA(B)/KIR conductance is high and the NMDA conductance is low; the depolarized state is robust because, at depolarized potentials, the NMDA conductance is high and the GABA(B)/KIR conductance is low. Our results suggest that this complementary voltage dependence of GABA(B)/KIR and NMDA conductances makes them a "perfect couple" for producing voltage bistability."
108. Nonlinear dendritic processing in barrel cortex spiny stellate neurons (Lavzin et al. 2012)
This is a multi-compartmental simulation of a spiny stellate neuron which is stimulated by a thalamocortical (TC) and cortico-cortical (CC) inputs. No other cells are explicitly modeled; the presynaptic network activation is represented by the number of active synapses. Preferred and non –preferred thalamic directions thus correspond to larder/smaller number of TC synapses. This simulation revealed that randomly activated synapses can cooperatively trigger global NMDA spikes, which involve participation of most of the dendritic tree. Surprisingly, we found that although the voltage profile of the cell was uniform, the calcium influx was restricted to ‘hot spots’ which correspond to synaptic clusters or large conductance synapses
109. Normal ripples, abnormal ripples, and fast ripples in a hippocampal model (Fink et al. 2015)
"...We use a computational model of hippocampus to investigate possible network mechanisms underpinning normal ripples, pathological ripples, and fast ripples. Our results unify several prior findings regarding HFO mechanisms, and also make several new predictions regarding abnormal HFOs. We show that HFOs are generic, emergent phenomena whose characteristics reflect a wide range of connectivity and network input. Although produced by different mechanisms, both normal and abnormal HFOs generate similar ripple frequencies, underscoring that peak frequency is unable to distinguish the two. Abnormal ripples are generic phenomena that arise when input to pyramidal cells overcomes network inhibition, resulting in high-frequency, uncoordinated firing. In addition, fast ripples transiently and sporadically arise from the precise conditions that produce abnormal ripples. Lastly, we show that such abnormal conditions do not require any specific network structure to produce coherent HFOs, as even completely asynchronous activity is capable of producing abnormal ripples and fast ripples in this manner. These results provide a generic, network-based explanation for the link between pathological ripples and fast ripples, and a unifying description for the entire spectrum from normal ripples to pathological fast ripples."
110. Olfactory bulb cluster formation (Migliore et al. 2010)
Functional roles of distributed synaptic clusters in the mitral-granule cell network of the olfactory bulb.
111. Olfactory bulb microcircuits model with dual-layer inhibition (Gilra & Bhalla 2015)
A detailed network model of the dual-layer dendro-dendritic inhibitory microcircuits in the rat olfactory bulb comprising compartmental mitral, granule and PG cells developed by Aditya Gilra, Upinder S. Bhalla (2015). All cell morphologies and network connections are in NeuroML v1.8.0. PG and granule cell channels and synapses are also in NeuroML v1.8.0. Mitral cell channels and synapses are in native python.
112. Olfactory bulb mitral and granule cell column formation (Migliore et al. 2007)
In the olfactory bulb, the processing units for odor discrimination are believed to involve dendrodendritic synaptic interactions between mitral and granule cells. There is increasing anatomical evidence that these cells are organized in columns, and that the columns processing a given odor are arranged in widely distributed arrays. Experimental evidence is lacking on the underlying learning mechanisms for how these columns and arrays are formed. We have used a simplified realistic circuit model to test the hypothesis that distributed connectivity can self-organize through an activity-dependent dendrodendritic synaptic mechanism. The results point to action potentials propagating in the mitral cell lateral dendrites as playing a critical role in this mechanism, and suggest a novel and robust learning mechanism for the development of distributed processing units in a cortical structure.
113. Olfactory bulb mitral and granule cell: dendrodendritic microcircuits (Migliore and Shepherd 2008)
This model shows how backpropagating action potentials in the long lateral dendrites of mitral cells, together with granule cell actions on mitral cells within narrow columns forming glomerular units, can provide a mechanism to activate strong local inhibition between arbitrarily distant mitral cells. The simulations predict a new role for the dendrodendritic synapses in the multicolumnar organization of the granule cells.
114. Olfactory bulb mitral cell gap junction NN model: burst firing and synchrony (O`Connor et al. 2012)
In a network of 6 mitral cells connected by gap junction in the apical dendrite tuft, continuous current injections of 0.06 nA are injected into 20 locations in the apical tufts of two of the mitral cells. The current injections into one of the cells starts 10 ms after the other to generate asynchronous firing in the cells (Migliore et al. 2005 protocol). Firing of the cells is asynchronous for the first 120 ms. However after the burst firing phase is completed the firing in all cells becomes synchronous.
115. Olfactory bulb mitral cell: synchronization by gap junctions (Migliore et al 2005)
In a realistic model of two electrically connected mitral cells, the paper shows that the somatically-measured experimental properties of Gap Junctions (GJs) may correspond to a variety of different local coupling strengths and dendritic distributions of GJs in the tuft. The model suggests that the propagation of the GJ-induced local tuft depolarization is a major mechanim for intraglomerular synchronization of mitral cells.
116. Olfactory Bulb Network (Davison et al 2003)
A biologically-detailed model of the mammalian olfactory bulb, incorporating the mitral and granule cells and the dendrodendritic synapses between them. The results of simulation experiments with electrical stimulation agree closely in most details with published experimental data. The model predicts that the time course of dendrodendritic inhibition is dependent on the network connectivity as well as on the intrinsic parameters of the synapses. In response to simulated odor stimulation, strongly activated mitral cells tend to suppress neighboring cells, the mitral cells readily synchronize their firing, and increasing the stimulus intensity increases the degree of synchronization. For more details, see the reference below.
117. Olfactory bulb network model of gamma oscillations (Bathellier et al. 2006; Lagier et al. 2007)
This model implements a network of 100 mitral cells connected with asynchronous inhibitory "synapses" that is meant to reproduce the GABAergic transmission of ensembles of connected granule cells. For appropriate parameters of this special synapse the model generates gamma oscillations with properties very similar to what is observed in olfactory bulb slices (See Bathellier et al. 2006, Lagier et al. 2007). Mitral cells are modeled as single compartment neurons with a small number of different voltage gated channels. Parameters were tuned to reproduce the fast subthreshold oscillation of the membrane potential observed experimentally (see Desmaisons et al. 1999).
118. Olfactory Computations in Mitral-Granule cell circuits (Migliore & McTavish 2013)
Model files for the entry "Olfactory Computations in Mitral-Granule Cell Circuits" of the Springer Encyclopedia of Computational Neuroscience by Michele Migliore and Tom Mctavish. The simulations illustrate two typical Mitral-Granule cell circuits in the olfactory bulb of vertebrates: distance-independent lateral inhibition and gating effects.
119. Optimal deep brain stimulation of the subthalamic nucleus-a computational study (Feng et al. 2007)
Here, we use a biophysically-based model of spiking cells in the basal ganglia (Terman et al., Journal of Neuroscience, 22, 2963-2976, 2002; Rubin and Terman, Journal of Computational Neuroscience, 16, 211-235, 2004) to provide computational evidence that alternative temporal patterns of DBS inputs might be equally effective as the standard high-frequency waveforms, but require lower amplitudes. Within this model, DBS performance is assessed in two ways. First, we determine the extent to which DBS causes Gpi (globus pallidus pars interna) synaptic outputs, which are burstlike and synchronized in the unstimulated Parkinsonian state, to cease their pathological modulation of simulated thalamocortical cells. Second, we evaluate how DBS affects the GPi cells' auto- and cross-correlograms.
120. Parallel odor processing by mitral and middle tufted cells in the OB (Cavarretta et al 2016, 2018)
"[...] experimental findings suggest that MC and mTC may encode parallel and complementary odor representations. We have analyzed the functional roles of these pathways by using a morphologically and physiologically realistic three-dimensional model to explore the MC and mTC microcircuits in the glomerular layer and deeper plexiform layers. [...]"
121. Parametric computation and persistent gamma in a cortical model (Chambers et al. 2012)
Using the Traub et al (2005) model of the cortex we determined how 33 synaptic strength parameters control gamma oscillations. We used fractional factorial design to reduce the number of runs required to 4096. We found an expected multiplicative interaction between parameters.
122. Parvalbumin-positive basket cells differentiate among hippocampal pyramidal cells (Lee et al. 2014)
This detailed microcircuit model explores the network level effects of sublayer specific connectivity in the mouse CA1. The differences in strengths and numbers of synapses between PV+ basket cells and either superficial sublayer or deep sublayer pyramidal cells enables a routing of inhibition from superficial to deep pyramidal cells. At the network level of this model, the effects become quite prominent when one compares the effect on firing rates when either the deep or superficial pyramidal cells receive a selective increase in excitation.
123. Persistent synchronized bursting activity in cortical tissues (Golomb et al 2005)
The program simulates a one-dimensional model of a cortical tissue with excitatory and inhibitory populations.
124. Principles of Computational Modelling in Neuroscience (Book) (Sterratt et al. 2011)
"... This book provides a step-by-step account of how to model the neuron and neural circuitry to understand the nervous system at all levels, from ion channels to networks. Starting with a simple model of the neuron as an electrical circuit, gradually more details are added to include the effects of neuronal morphology, synapses, ion channels and intracellular signaling. The principle of abstraction is explained through chapters on simplifying models, and how simplified models can be used in networks. This theme is continued in a final chapter on modeling the development of the nervous system. Requiring an elementary background in neuroscience and some high school mathematics, this textbook is an ideal basis for a course on computational neuroscience."
125. Prosthetic electrostimulation for information flow repair in a neocortical simulation (Kerr 2012)
This model is an extension of a model (<a href="">138379</a>) recently published in Frontiers in Computational Neuroscience. This model consists of 4700 event-driven, rule-based neurons, wired according to anatomical data, and driven by both white-noise synaptic inputs and a sensory signal recorded from a rat thalamus. Its purpose is to explore the effects of cortical damage, along with the repair of this damage via a neuroprosthesis.
126. Pyramidal neuron, fast, regular, and irregular spiking interneurons (Konstantoudaki et al 2014)
This is a model network of prefrontal cortical microcircuit based primarily on rodent data. It includes 16 pyramidal model neurons, 2 fast spiking interneuron models, 1 regular spiking interneuron model and 1 irregular spiking interneuron model. The goal of the paper was to use this model network to determine the role of specific interneuron subtypes in persistent activity
127. Rapid desynchronization of an electrically coupled Golgi cell network (Vervaeke et al. 2010)
Electrical synapses between interneurons contribute to synchronized firing and network oscillations in the brain. However, little is known about how such networks respond to excitatory synaptic input. In addition to detailed electrophysiological recordings and histological investigations of electrically coupled Golgi cells in the cerebellum, a detailed network model of these cells was created. The cell models are based on reconstructed Golgi cell morphologies and the active conductances are taken from an earlier abstract Golgi cell model (Solinas et al 2007, accession no. 112685). Our results show that gap junction coupling can sometimes be inhibitory and either promote network synchronization or trigger rapid network desynchronization depending on the synaptic input. The model is available as a neuroConstruct project and can executable scripts can be generated for the NEURON simulator.
128. Regulation of a slow STG rhythm (Nadim et al 1998)
Frequency regulation of a slow rhythm by a fast periodic input. Nadim, F., Manor, Y., Nusbaum, M. P., Marder, E. (1998) J. Neurosci. 18: 5053-5067
129. Response properties of neocort. neurons to temporally modulated noisy inputs (Koendgen et al. 2008)
Neocortical neurons are classified by current–frequency relationship. This is a static description and it may be inadequate to interpret neuronal responses to time-varying stimuli. Theoretical studies (Brunel et al., 2001; Fourcaud-Trocmé et al. 2003; Fourcaud-Trocmé and Brunel 2005; Naundorf et al. 2005) suggested that single-cell dynamical response properties are necessary to interpret ensemble responses to fast input transients. Further, it was shown that input-noise linearizes and boosts the response bandwidth, and that the interplay between the barrage of noisy synaptic currents and the spike-initiation mechanisms determine the dynamical properties of the firing rate. In order to allow a reader to explore such simulations, we prepared a simple NEURON implementation of the experiments performed in Köndgen et al., 2008 (see also Fourcaud-Trocmé al. 2003; Fourcaud-Trocmé and Brunel 2005). In addition, we provide sample MATLAB routines for exploring the sandwich model proposed in Köndgen et al., 2008, employing a simple frequdency-domain filtering. The simulations and the MATLAB routines are based on the linear response properties of layer 5 pyramidal cells estimated by injecting a superposition of a small-amplitude sinusoidal wave and a background noise, as in Köndgen et al., 2008.
130. Reverberatory bursts propagation and synchronization in developing cultured NNs (Huang et al 2016)
"Developing networks of neural systems can exhibit spontaneous, synchronous activities called neural bursts, which can be important in the organization of functional neural circuits. ... Using a propagation model we infer the spreading speed of the spiking activity, which increases as the culture ages. We perform computer simulations of the system using a physiological model of spiking networks in two spatial dimensions and find the parameters that reproduce the observed resynchronization of spiking in the bursts. An analysis of the simulated dynamics suggests that the depletion of synaptic resources causes the resynchronization. The spatial propagation dynamics of the simulations match well with observations over the course of a burst and point to an interplay of the synaptic efficacy and the noisy neural self-activation in producing the morphology of the bursts."
131. S cell network (Moss et al 2005)
Excerpts from the abstract: S cells form a chain of electrically coupled neurons that extends the length of the leech CNS and plays a critical role in sensitization during whole-body shortening. ... Serotonin ... increasedAP latency across the electrical synapse, suggesting that serotonin reduced coupling between S cells. ... Serotonin modulated instantaneous AP frequency when APs were initiated in separate S cells and in a computational model of S cell activity following mechanosensory input. Thus, serotonergic modulation of S cell electrical synapses may contribute to changes in the pattern of activity in the S cell network. See paper for more.
132. Self-organized olfactory pattern recognition (Kaplan & Lansner 2014)
" ... We present a large-scale network model with single and multi-compartmental Hodgkin–Huxley type model neurons representing olfactory receptor neurons (ORNs) in the epithelium, periglomerular cells, mitral/tufted cells and granule cells in the olfactory bulb (OB), and three types of cortical cells in the piriform cortex (PC). Odor patterns are calculated based on affinities between ORNs and odor stimuli derived from physico-chemical descriptors of behaviorally relevant real-world odorants. ... The PC was implemented as a modular attractor network with a recurrent connectivity that was likewise organized through Hebbian–Bayesian learning. We demonstrate the functionality of the model in a one-sniff-learning and recognition task on a set of 50 odorants. Furthermore, we study its robustness against noise on the receptor level and its ability to perform concentration invariant odor recognition. Moreover, we investigate the pattern completion capabilities of the system and rivalry dynamics for odor mixtures."
133. Simulated cortical color opponent receptive fields self-organize via STDP (Eguchi et al., 2014)
"... In this work, we address the problem of understanding the cortical processing of color information with a possible mechanism of the development of the patchy distribution of color selectivity via computational modeling. ... Our model of the early visual system consists of multiple topographically-arranged layers of excitatory and inhibitory neurons, with sparse intra-layer connectivity and feed-forward connectivity between layers. Layers are arranged based on anatomy of early visual pathways, and include a retina, lateral geniculate nucleus, and layered neocortex. ... After training with natural images, the neurons display heightened sensitivity to specific colors. ..."
134. Simulation studies on mechanisms of levetiracetam-mediated inhibition of IK(DR) (Huang et al. 2009)
Levetiracetam (LEV) is an S-enantiomer pyrrolidone derivative with established antiepileptic efficacy in generalized epilepsy and partial epilepsy. However, its effects on ion currents and membrane potential remain largely unclear. In this study, we investigated the effect of LEV on differentiated NG108-15 neurons. ... Simulation studies in a modified Hodgkin-Huxley neuron and network unraveled that the reduction of slowly inactivating IK(DR) resulted in membrane depolarization accompanied by termination of the firing of action potentials in a stochastic manner. Therefore, the inhibitory effects on slowly inactivating IK(DR) (Kv3.1-encoded current) may constitute one of the underlying mechanisms through which LEV affects neuronal activity in vivo.
135. Simulations of oscillations in piriform cortex (Wilson & Bower 1992)
"1. A large-scale computer model of the piriform cortex was constructed on the basis of the known anatomic and physiological organization of this region. 2. The oscillatory field potential and electroencephalographic (EEG) activity generated by the model was compared with actual physiological results. The model was able to produce patterns of activity similar to those recorded physiologically in response to both weak and strong electrical shocks to the afferent input. The model also generated activity patterns similar to EEGs recorded in behaving animals. 3. ..."
136. Single compartment Dorsal Lateral Medium Spiny Neuron w/ NMDA and AMPA (Biddell and Johnson 2013)
A biophysical single compartment model of the dorsal lateral striatum medium spiny neuron is presented here. The model is an implementation then adaptation of a previously described model (Mahon et al. 2002). The model has been adapted to include NMDA and AMPA receptor models that have been fit to dorsal lateral striatal neurons. The receptor models allow for excitation by other neuron models.
137. Single E-I oscillating network with amplitude modulation (Avella Gonzalez et al. 2012)
"... Intriguingly, the amplitude of ongoing oscillations, such as measured in EEG recordings, fluctuates irregularly, with episodes of high amplitude (HAE) alternating with episodes of low amplitude (LAE). ... Here, we show that transitions between HAE and LAE in the alpha/beta frequency band occur in a generic neuronal network model consisting of interconnected inhibitory (I) and excitatory (E) cells that are externally driven by sustained depolarizing currents(cholinergic input) and trains of action potentials that activate excitatory synapses. In the model, action potentials onto inhibitory cells represent input from other brain areas and desynchronize network activity, being crucial for the emergence of amplitude fluctuations. ..."
138. Sleep-wake transitions in corticothalamic system (Bazhenov et al 2002)
The authors investigate the transition between sleep and awake states with intracellular recordings in cats and computational models. The model describes many essential features of slow wave sleep and activated states as well as the transition between them.
139. Small world networks of Type I and Type II Excitable Neurons (Bogaard et al. 2009)
Implemented with NEURON 5.9, four model neurons with varying excitability properties affect the spatiotemporal patterning of small world networks of homogeneous and heterogeneous cell population.
140. Software for teaching neurophysiology of neuronal circuits (Grisham et al. 2008)
"To circumvent the many problems in teaching neurophysiology as a “wet lab,” we developed SWIMMY, a virtual fish that swims by moving its virtual tail by means of a virtual neural circuit. ... Using SWIMMY, students (1) review the basics of neurophysiology, (2) identify the neurons in the circuit, (3) ascertain the neurons’ synaptic interconnections, (4) discover which cells generate the motor pattern of swimming, (5) discover how the rhythm is generated, and finally (6) use an animation that corresponds to the activity of the motoneurons to discover the behavioral effects produced by various lesions and explain them in terms of their neural underpinnings. ..."
141. Spikes,synchrony,and attentive learning by laminar thalamocort. circuits (Grossberg & Versace 2007)
"... The model hereby clarifies, for the first time, how the following levels of brain organization coexist to realize cognitive processing properties that regulate fast learning and stable memory of brain representations: single cell properties, such as spiking dynamics, spike-timing-dependent plasticity (STDP), and acetylcholine modulation; detailed laminar thalamic and cortical circuit designs and their interactions; aggregate cell recordings, such as current-source densities and local field potentials; and single cell and large-scale inter-areal oscillations in the gamma and beta frequency domains. ..."
142. State dependent drug binding to sodium channels in the dentate gyrus (Thomas & Petrou 2013)
A Markov model of sodium channels was developed that includes drug binding to fast inactivated states. This was incorporated into a model of the dentate gyrus to investigate the effects of anti-epileptic drugs on neuron and network properties.
143. State-dependent rhythmogenesis in a half-center locomotor CPG (Ausburn et al 2017)
"The spinal locomotor central pattern generator (CPG) generates rhythmic activity with alternating flexion and extension phases. This rhythmic pattern is likely to result from inhibitory interactions between neural populations representing flexor and extensor half-centers. However, it is unclear whether the flexor-extensor CPG has a quasi-symmetric organization with both half-centers critically involved in rhythm generation, features an asymmetric organization with flexor-driven rhythmogenesis, or comprises a pair of intrinsically rhythmic half-centers. There are experimental data that support each of the above concepts but appear to be inconsistent with the others. In this theoretical/modeling study, we present and analyze a CPG model architecture that can operate in different regimes consistent with the above three concepts depending on conditions, which are defined by external excitatory drives to CPG half-centers. We show that control of frequency and phase durations within each regime depends on network dynamics, defined by the regime-dependent expression of the half-centers' intrinsic rhythmic capabilities and the operating phase transition mechanisms (escape vs. release). Our study suggests state dependency in locomotor CPG operation and proposes explanations for seemingly contradictory experimental data."
144. STDP promotes synchrony of inhibitory networks in the presence of heterogeneity (Talathi et al 2008)
"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 tpre and the postsynaptic excitation at time tpost, the synapse is strengthened (delta_t = t(post) - t(pre) > 0) or weakened (delta_t < 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. ..."
145. Striatal NN model of MSNs and FSIs investigated effects of dopamine depletion (Damodaran et al 2015)
This study investigates the mechanisms that are affected in the striatal network after dopamine depletion and identifies potential therapeutic targets to restore normal activity.
146. Structure-dynamics relationships in bursting neuronal networks revealed (Mäki-Marttunen et al. 2013)
This entry includes tools for generating and analyzing network structure, and for running the neuronal network simulations on them.
147. Studies of stimulus parameters for seizure disruption using NN simulations (Anderson et al. 2007)
Architecturally realistic neocortical model using seven classes of excitatory and inhibitory single compartment Hodgkin-Huxley cells. Wiring is adapted to minicolumn hypothesis and incorporates visual and neocortical data. Simulation demonstrates spontaneous bursting onset and cessation, and activity can be altered with external electric field.
148. Study of augmented Rubin and Terman 2004 deep brain stim. model in Parkinsons (Pascual et al. 2006)
" ... The model by Rubin and Terman [31] represents one of the most comprehensive and biologically plausible models of DBS published recently. We examined the validity of the model, replicated its simulations and tested its robustness. While our simulations partially reproduced the results presented by Rubin and Terman [31], several issues were raised including the high complexity of the model in its non simplified form, the lack of robustness of the model with respect to small perturbations, the nonrealistic representation of the thalamus and the absence of time delays. Computational models are indeed necessary, but they may not be sufficient in their current forms to explain the effect of chronic electrical stimulation on the activity of the basal ganglia (BG) network in PD."
149. Subiculum network model with dynamic chloride/potassium homeostasis (Buchin et al 2016)
This is the code implementing the single neuron and spiking neural network dynamics. The network has the dynamic ion concentrations of extracellular potassium and intracellular chloride. The code contains multiple parameter variations to study various mechanisms of the neural excitability in the context of chloride homeostasis.
150. Synaptic gating at axonal branches, and sharp-wave ripples with replay (Vladimirov et al. 2013)
The computational model of in vivo sharp-wave ripples with place cell replay. Excitatory post-synaptic potentials at dendrites gate antidromic spikes arriving from the axonal collateral, and thus determine when the soma and the main axon fire. The model allows synchronous replay of pyramidal cells during sharp-wave ripple event, and the replay is possible in both forward and reverse directions.
151. Synaptic information transfer in computer models of neocortical columns (Neymotin et al. 2010)
"... We sought to measure how the activity of the network alters information flow from inputs to output patterns. Information handling by the network reflected the degree of internal connectivity. ... With greater connectivity strength, the recurrent network translated activity and information due to contribution of activity from intrinsic network dynamics. ... At still higher internal synaptic strength, the network corrupted the external information, producing a state where little external information came through. The association of increased information retrieved from the network with increased gamma power supports the notion of gamma oscillations playing a role in information processing."
152. Synchronization by D4 dopamine receptor-mediated phospholipid methylation (Kuznetsova, Deth 2008)
"We describe a new molecular mechanism of dopamine-induced membrane protein modulation that can tune neuronal oscillation frequency to attention related gamma rhythm. This mechanism is based on the unique ability of D4 dopamine receptors (D4R) to carry out phospholipid methylation (PLM) that may affect the kinetics of ion channels. We show that by deceasing the inertia of the delayed rectifier potassium channel, a transition to 40 Hz oscillations can be achieved. ..."
153. Synchrony by synapse location (McTavish et al. 2012)
This model considers synchrony between mitral cells induced via shared granule cell interneurons while taking into account the spatial constraints of the system. In particular, since inhibitory inputs decay passively along the lateral dendrites, this model demonstrates that an optimal arrangement of the inhibitory synapses will be near the cell bodies of the relevant mitral cells.
154. Thalamic transformation of pallidal input (Hadipour-Niktarash 2006)
"In Parkinson’s disease, neurons of the internal segment of the globus pallidus (GPi) display the low-frequency tremor-related oscillations. These oscillatory activities are transmitted to the thalamic relay nuclei. Computer models of the interacting thalamocortical (TC) and thalamic reticular (RE) neurons were used to explore how the TC-RE network processes the low-frequency oscillations of the GPi neurons. ..."
155. Thalamocortical and Thalamic Reticular Network (Destexhe et al 1996)
NEURON model of oscillations in networks of thalamocortical and thalamic reticular neurons in the ferret. (more applications for a model quantitatively identical to previous DLGN model; updated for NEURON v4 and above)
156. Thalamocortical augmenting response (Bazhenov et al 1998)
In the cortical model, augmenting responses were more powerful in the "input" layer compared with those in the "output" layer. Cortical stimulation of the network model produced augmenting responses in cortical neurons in distant cortical areas through corticothalamocortical loops and low-threshold intrathalamic augmentation. ... The predictions of the model were compared with in vivo recordings from neurons in cortical area 4 and thalamic ventrolateral nucleus of anesthetized cats. The known intrinsic properties of thalamic cells and thalamocortical interconnections can account for the basic properties of cortical augmenting responses. See reference for details. NEURON implementation note: cortical SU cells are getting slightly too little stimulation - reason unknown.
157. The activity phase of postsynaptic neurons (Bose et al 2004)
We show, in a simplified network consisting of an oscillator inhibiting a follower neuron, how the interaction between synaptic depression and a transient potassium current in the follower neuron determines the activity phase of this neuron. We derive a mathematical expression to determine at what phase of the oscillation the follower neuron becomes active. This expression can be used to understand which parameters determine the phase of activity of the follower as the frequency of the oscillator is changed. See paper for more.
158. The origin of different spike and wave-like events (Hall et al 2017)
Acute In vitro models have revealed a great deal of information about mechanisms underlying many types of epileptiform activity. However, few examples exist that shed light on spike and wave (SpW) patterns of pathological activity. SpW are seen in many epilepsy syndromes, both generalised and focal, and manifest across the entire age spectrum. They are heterogeneous in terms of their severity, symptom burden and apparent anatomical origin (thalamic, neocortical or both), but any relationship between this heterogeneity and underlying pathology remains elusive. Here we demonstrate that physiological delta frequency rhythms act as an effective substrate to permit modelling of SpW of cortical origin and may help to address this issue. ..."
159. Turtle visual cortex model (Nenadic et al. 2003, Wang et al. 2005, Wang et al. 2006)
This is a model of the visual cortex of freshwater turtles that is based upon the known anatomy and physiology of individual neurons. The model was published in three papers (Nenadic et al., 2003; Wang et al., 2005; Wang et al., 2006), which should be consulted for full details on its construction. The model has also been used in several papers (Robbins and Senseman, 2004; Du et al., 2005; Du et al., 2006). It is implemented in GENESIS (Bower and Beeman, 1998).
160. Unbalanced peptidergic inhibition in superficial cortex underlies seizure activity (Hall et al 2015)
" ...Loss of tonic neuromodulatory excitation, mediated by nicotinic acetylcholine or serotonin (5HT3A) receptors, of 5HT3-immunopositive interneurons caused an increase in amplitude and slowing of the delta rhythm until each period became the "wave" component of the spike and wave discharge. As with the normal delta rhythm, the wave of a spike and wave discharge originated in cortical layer 5. In contrast, the "spike" component of the spike and wave discharge originated from a relative failure of fast inhibition in layers 2/3-switching pyramidal cell action potential outputs from single, sparse spiking during delta rhythms to brief, intense burst spiking, phase-locked to the field spike. The mechanisms underlying this loss of superficial layer fast inhibition, and a concomitant increase in slow inhibition, appeared to be precipitated by a loss of neuropeptide Y (NPY)-mediated local circuit inhibition and a subsequent increase in vasoactive intestinal peptide (VIP)-mediated disinhibition. Blockade of NPY Y1 receptors was sufficient to generate spike and wave discharges, whereas blockade of VIP receptors almost completely abolished this form of epileptiform activity. These data suggest that aberrant, activity-dependent neuropeptide corelease can have catastrophic effects on neocortical dynamics."
161. Updated Tritonia Swim CPG (Calin-Jagemann et al. 2007)
Model of the 3-cell core CPG (DSI, C2, and VSI-B) mediating escape swimming in Tritonia diomedea. Cells use a hybrid integrate-and-fire scheme pioneered by Peter Getting. Each model cell is reconstructed from extensive physiological measurements to precisely mimic I-F curves, synaptic waveforms, and functional connectivity.
162. Visual physiology of the layer 4 cortical circuit in silico (Arkhipov et al 2018)
"Despite advances in experimental techniques and accumulation of large datasets concerning the composition and properties of the cortex, quantitative modeling of cortical circuits under in-vivo-like conditions remains challenging. Here we report and publicly release a biophysically detailed circuit model of layer 4 in the mouse primary visual cortex, receiving thalamo- cortical visual inputs. The 45,000-neuron model was subjected to a battery of visual stimuli, and results were compared to published work and new in vivo experiments. ..."

Re-display model names without descriptions