Circuits that contain the Current : I Na,t

("Transient"; rapidly inactivating)
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 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.
6. 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.
7. A network of AOB mitral cells that produces infra-slow bursting (Zylbertal et al. 2017)
Infra-slow rhythmic neuronal activity with very long (> 10 s) period duration was described in many brain areas but little is known about the role of this activity and the mechanisms that produce it. Here we combine experimental and computational methods to show that synchronous infra-slow bursting activity in mitral cells of the mouse accessory olfactory bulb (AOB) emerges from interplay between intracellular dynamics and network connectivity. In this novel mechanism, slow intracellular Na+ dynamics endow AOB mitral cells with a weak tendency to burst, which is further enhanced and stabilized by chemical and electrical synapses between them. Combined with the unique topology of the AOB network, infra-slow bursting enables integration and binding of multiple chemosensory stimuli over prolonged time scale. The example protocol simulates a two-glomeruli network with a single shared cell. Although each glomerulus is stimulated at a different time point, the activity of the entire population becomes synchronous (see paper Fig. 8)
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 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.
11. 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.
12. Alpha rhythm in vitro visual cortex (Traub et al 2020)
The paper describes an experimental model of the alpha rhythm generated by layer 4 pyramidal neurons in a visual cortex slice. The simulation model is derived from that of Traub et al. (2005) J Neurophysiol, developed for thalamocortical oscillations.
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. 2012)
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 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.
20. 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.
21. 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. ..."
22. 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."
23. 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."
24. 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.
25. 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.
26. 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.
27. 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.
28. Current Dipole in Laminar Neocortex (Lee et al. 2013)
Laminar neocortical model in NEURON/Python, adapted from Jones et al 2009. https://bitbucket.org/jonescompneurolab/corticaldipole
29. 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.
30. 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."
31. 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.
32. 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.
33. 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. ..."
34. Electrodecrements in in vitro model of infantile spasms (Traub et al 2020)
The code is an extension of the thalamocortical model of Traub et al. (2005) J Neurophysiol. It is here applied to an in vitro model of the electrodecremental response seen in the EEG of children with infantile spasms (West syndrome)
35. 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. ..."
36. 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."
37. 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.
38. 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.
39. 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.
40. 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). ..."
41. 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.
42. 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.
43. 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.
44. 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.
45. 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. ..."
46. 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.
47. 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.
48. 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.
49. 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."
50. 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. ..."
51. 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.
52. 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 http://bitbucket.org/mbezaire/ca1 and further information at http://mariannebezaire.com/models/ca1
53. 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).
54. 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.
55. Homeostatic mechanisms may shape oscillatory modulations (Peterson & Voytek 2020)
"Neural oscillations are observed ubiquitously in the mammalian brain, but their stability is known to be rather variable. Some oscillations are tonic and last for seconds or even minutes. Other oscillations appear as unstable bursts. Likewise, some oscillations rely on excitatory AMPAergic synapses, but others are GABAergic and inhibitory. Why this diversity exists is not clear. We hypothesized Ca2+-dependent homeostasis could be important in finding an explanation. We tested this hypothesis in a highly simplified model of hippocampal neurons. In this model homeostasis profoundly alters the modulatory effect of neural oscillations. Under homeostasis, tonic AMPAergic oscillations actually decrease excitability and desynchronize firing. Tonic oscillations that are synaptically GABAergic-like those in real hippocampus-don't provoke a homeostatic response, however. If our simple model is correct, homeostasis can explain why the theta rhythm in the hippocampus is synaptically inhibitory: GABA has little to no intrinsic homeostatic response, and so can preserve the pyramidal cell's natural dynamic range. Based on these results we can also speculate that homeostasis may explain why AMPAergic oscillations in cortex, and hippocampus, often appear as bursts. Bursts do not interact with the slow homeostatic time constant, and so retain their normal excitatory effect."
56. 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.
57. 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. ..."
58. 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.
59. 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."
60. 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.
61. 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="http://dx.doi.org/10.1007/s10827-007-0038-6">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 .
62. 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...."
63. 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.
64. 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.
65. 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).
66. 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.
67. 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.
68. 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).
69. 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..."
70. 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.
71. 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. "
72. 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."
73. 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. ..."
74. 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. ..."
75. 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: http://geza.kzoo.edu/theta). 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.
76. 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.
77. 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. ..."
78. 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)
79. 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.
80. 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.
81. 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. ..."
82. 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."
83. Olfactory bulb cluster formation (Migliore et al. 2010)
Functional roles of distributed synaptic clusters in the mitral-granule cell network of the olfactory bulb.
84. 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.
85. 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.
86. 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.
87. 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.
88. 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.
89. 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).
90. 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.
91. 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. [...]"
92. 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.
93. 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."
94. 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
95. 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.
96. 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
97. 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.
98. 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."
99. Sensory-evoked responses of L5 pyramidal tract neurons (Egger et al 2020)
This is the L5 pyramidal tract neuron (L5PT) model from Egger, Narayanan et al., Neuron 2020. It allows investigating how synaptic inputs evoked by different sensory stimuli are integrated by the complex intrinsic properties of L5PTs. The model is constrained by anatomical measurements of the subcellular synaptic input patterns to L5PT neurons, in vivo measurements of sensory-evoked responses of different populations of neurons providing these synaptic inputs, and in vitro measurements constraining the biophysical properties of the soma, dendrites and axon (note: the biophysical model is based on the work by Hay et al., Plos Comp Biol 2011). The model files provided here allow performing simulations and analyses presented in Figures 3, 4 and 5.
100. 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. ..."
101. 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.
102. 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.
103. 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. ..."
104. Spike burst-pause dynamics of Purkinje cells regulate sensorimotor adaptation (Luque et al 2019)
"Cerebellar Purkinje cells mediate accurate eye movement coordination. However, it remains unclear how oculomotor adaptation depends on the interplay between the characteristic Purkinje cell response patterns, namely tonic, bursting, and spike pauses. Here, a spiking cerebellar model assesses the role of Purkinje cell firing patterns in vestibular ocular reflex (VOR) adaptation. The model captures the cerebellar microcircuit properties and it incorporates spike-based synaptic plasticity at multiple cerebellar sites. ..."
105. 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. ..."
106. 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.
107. State-dependent rhythmogenesis in a half-center locomotor CPG (Ausborn 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."
108. 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. ..."
109. 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.
110. 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."
111. 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.
112. 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."
113. 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. ..."
114. 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.
115. Systematic integration of data into multi-scale models of mouse primary V1 (Billeh et al 2020)
"Highlights • Two network models of the mouse primary visual cortex are developed and released • One uses compartmental-neuron models and the other point-neuron models • The models recapitulate observations from in vivo experimental data • Simulations identify experimentally testable predictions about cortex circuitry"
116. 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. ..."
117. 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)
118. 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.
119. Thalamocortical control of propofol phase-amplitude coupling (Soplata et al 2017)
"The anesthetic propofol elicits many different spectral properties on the EEG, including alpha oscillations (8-12 Hz), Slow Wave Oscillations (SWO, 0.1-1.5 Hz), and dose-dependent phase-amplitude coupling (PAC) between alpha and SWO. Propofol is known to increase GABAA inhibition and decrease H-current strength, but how it generates these rhythms and their interactions is still unknown. To investigate both generation of the alpha rhythm and its PAC to SWO, we simulate a Hodgkin-Huxley network model of a hyperpolarized thalamus and corticothalamic inputs. ..."
120. The microcircuits of striatum in silico (Hjorth et al 2020)
"Our aim is to reconstruct a full-scale mouse striatal cellular level model to provide a framework to integrate and interpret striatal data. We represent the main striatal neuronal subtypes, the two types of projection neurons (dSPNs and iSPNs) giving rise to the direct and indirect pathways, the fast-spiking interneurons, the low threshold spiking interneurons, and the cholinergic interneurons as detailed compartmental models, with properties close to their biological counterparts. Both intrastriatal and afferent synaptic inputs (cortex, thalamus, dopamine system) are optimized against existing data, including short-term plasticity. This model platform will be used to generate new hypotheses on striatal function or network dynamic phenomena."
121. 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. ..."
122. 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).
123. 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."
124. 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