Bursting Reverberation as a Multiscale Neuronal Network Process Driven by Synaptic Depression-Facilitation

Neuronal networks can generate complex patterns of activity that depend on membrane properties of individual neurons as well as on functional synapses. To decipher the impact of synaptic properties and connectivity on neuronal network behavior, we investigate the responses of neuronal ensembles from small (5–30 cells in a restricted sphere) and large (acute hippocampal slice) networks to single electrical stimulation: in both cases, a single stimulus generated a synchronous long-lasting bursting activity. While an initial spike triggered a reverberating network activity that lasted 2–5 seconds for small networks, we found here that it lasted only up to 300 milliseconds in slices. To explain this phenomena present at different scales, we generalize the depression-facilitation model and extracted the network time constants. The model predicts that the reverberation time has a bell shaped relation with the synaptic density, revealing that the bursting time cannot exceed a maximum value. Furthermore, before reaching its maximum, the reverberation time increases sub-linearly with the synaptic density of the network. We conclude that synaptic dynamics and connectivity shape the mean burst duration, a property present at various scales of the networks. Thus bursting reverberation is a property of sufficiently connected neural networks, and can be generated by collective depression and facilitation of underlying functional synapses.     

Structure-Dynamics Relationships in Bursting Neuronal Networks Revealed Using a Prediction Framework

The question of how the structure of a neuronal network affects its functionality has gained a lot of attention in neuroscience. However, the vast majority of the studies on structure-dynamics relationships consider few types of network structures and assess limited numbers of structural measures. In this in silico study, we employ a wide diversity of network topologies and search among many possibilities the aspects of structure that have the greatest effect on the network excitability. The network activity is simulated using two point-neuron models, where the neurons are activated by noisy fluctuation of the membrane potential and their connections are described by chemical synapse models, and statistics on the number and quality of the emergent network bursts are collected for each network type. We apply a prediction framework to the obtained data in order to find out the most relevant aspects of network structure. In this framework, predictors that use different sets of graph-theoretic measures are trained to estimate the activity properties, such as burst count or burst length, of the networks. The performances of these predictors are compared with each other. We show that the best performance in prediction of activity properties for networks with sharp in-degree distribution is obtained when the prediction is based on clustering coefficient. By contrast, for networks with broad in-degree distribution, the maximum eigenvalue of the connectivity graph gives the most accurate prediction. The results shown for small () networks hold with few exceptions when different neuron models, different choices of neuron population and different average degrees are applied. We confirm our conclusions using larger () networks as well. Our findings reveal the relevance of different aspects of network structure from the viewpoint of network excitability, and our integrative method could serve as a general framework for structure-dynamics studies in biosciences.     

Calcium Involvement in Regulation of Neuronal Bursting in Disinhibited Neuronal Networks: Insights from Calcium Studies in a Spherical Cell Model

Cytosolic calcium is involved in the regulation of many intracellular processes. Intracellular calcium may therefore potentially affect the behavior of both single neurons and synaptically connected neuronal assemblies. In computer model studies, we investigated calcium dynamics in spherical neurons during periods of recurrent neuronal bursting that were simulated in a disinhibited neuronal network. The model takes into account calcium influx via voltage-gated calcium channels, extrusion through the cell membrane, and binding to two different buffers representing fixed and mobile endogenous calcium buffers. Throughout the duration of the simulated recurrent neuronal bursting, the concentration of free fixed buffers shows a hyperbolic decrease in time at a rate that is not uniform inside a neuron. Recurrent calcium influxes associated with bursting lead to the formation of gradients in the concentration of the fixed buffer in the radial direction, and are accompanied by the redistribution of mobile buffers acting to compensate for these gradients. Simulated intracellular calcium transients have a slow component characterized by a gradual increase in the calcium baseline level that reaches a plateau 120-200 s after the onset of recurrent bursting. Using this model, we demonstrate what we believe is a novel mechanism of regulation of network excitability that occurs in conditions of prolonged and recurrent neuronal bursting in disinhibited networks. This mechanism is expressed via interaction of calcium clearance systems inside neurons with calcium-dependent potassium regulation of neuronal excitability in membranes. This is a network phenomenon because it arises largely by synaptic interactions. Therefore, it can serve as a network safety mechanism to prevent excessive and uncontrolled neuronal firing resulting from the lack of inhibition or after acute suppression of the inhibitory drive.     

Adaptive Robotic Control Driven by a Versatile Spiking Cerebellar Network

The cerebellum is involved in a large number of different neural processes, especially in associative learning and in fine motor control. To develop a comprehensive theory of sensorimotor learning and control, it is crucial to determine the neural basis of coding and plasticity embedded into the cerebellar neural circuit and how they are translated into behavioral outcomes in learning paradigms. Learning has to be inferred from the interaction of an embodied system with its real environment, and the same cerebellar principles derived from cell physiology have to be able to drive a variety of tasks of different nature, calling for complex timing and movement patterns. We have coupled a realistic cerebellar spiking neural network (SNN) with a real robot and challenged it in multiple diverse sensorimotor tasks. Encoding and decoding strategies based on neuronal firing rates were applied. Adaptive motor control protocols with acquisition and extinction phases have been designed and tested, including an associative Pavlovian task (Eye blinking classical conditioning), a vestibulo-ocular task and a perturbed arm reaching task operating in closed-loop. The SNN processed in real-time mossy fiber inputs as arbitrary contextual signals, irrespective of whether they conveyed a tone, a vestibular stimulus or the position of a limb. A bidirectional long-term plasticity rule implemented at parallel fibers-Purkinje cell synapses modulated the output activity in the deep cerebellar nuclei. In all tasks, the neurorobot learned to adjust timing and gain of the motor responses by tuning its output discharge. It succeeded in reproducing how human biological systems acquire, extinguish and express knowledge of a noisy and changing world. By varying stimuli and perturbations patterns, real-time control robustness and generalizability were validated. The implicit spiking dynamics of the cerebellar model fulfill timing, prediction and learning functions.     

Shaping Neuronal Network Activity by Presynaptic Mechanisms

Neuronal microcircuits generate oscillatory activity, which has been linked to basic functions such as sleep, learning and sensorimotor gating. Although synaptic release processes are well known for their ability to shape the interaction between neurons in microcircuits, most computational models do not simulate the synaptic transmission process directly and hence cannot explain how changes in synaptic parameters alter neuronal network activity. In this paper, we present a novel neuronal network model that incorporates presynaptic release mechanisms, such as vesicle pool dynamics and calcium-dependent release probability, to model the spontaneous activity of neuronal networks. The model, which is based on modified leaky integrate-and-fire neurons, generates spontaneous network activity patterns, which are similar to experimental data and robust under changes in the model''s primary gain parameters such as excitatory postsynaptic potential and connectivity ratio. Furthermore, it reliably recreates experimental findings and provides mechanistic explanations for data obtained from microelectrode array recordings, such as network burst termination and the effects of pharmacological and genetic manipulations. The model demonstrates how elevated asynchronous release, but not spontaneous release, synchronizes neuronal network activity and reveals that asynchronous release enhances utilization of the recycling vesicle pool to induce the network effect. The model further predicts a positive correlation between vesicle priming at the single-neuron level and burst frequency at the network level; this prediction is supported by experimental findings. Thus, the model is utilized to reveal how synaptic release processes at the neuronal level govern activity patterns and synchronization at the network level.     

Noise-Induced Coherence and Network Oscillations in a Reduced Bursting Model

The dynamics of the Hindmarsh-Rose (HR) model of bursting thalamic neurons is reduced to a system of two linear differential equations that retains the subthreshold resonance properties of the HR model. Introducing a reset mechanism after a threshold crossing, we turn this system into a resonant integrate-and-fire (RIF) model. Using Monte-Carlo simulations and mathematical analysis, we examine the effects of noise and the subthreshold dynamic properties of the RIF model on the occurrence of coherence resonance (CR). Synchronized burst firing occurs in a network of such model neurons with excitatory pulse-coupling. The coherence level of the network oscillations shows a stochastic resonance-like dependence on the noise level. Stochastic analysis of the equations shows that the slow recovery from the spike-induced inhibition is crucial in determining the frequencies of the CR and the subthreshold resonance in the original HR model. In this particular type of CR, the oscillation frequency strongly depends on the intrinsic time scales but changes little with the noise intensity. We give analytical quantities to describe this CR mechanism and illustrate its influence on the emerging network oscillations. We discuss the profound physiological roles this kind of CR may have in information processing in neurons possessing a subthreshold resonant frequency and in generating synchronized network oscillations with a frequency that is determined by intrinsic properties of the neurons. PACS 05.45.-a, 05.40.Ca, 87.18.Sn, 87.19     

Pharmacological Characterization of Cultivated Neuronal Networks: Relevance to Synaptogenesis and Synaptic Connectivity

Mental disorders, such as schizophrenia or Alzheimer’s disease, are associated with impaired synaptogenesis and/or synaptic communication. During development, neurons assemble into neuronal networks, the primary supracellular mediators of information processing. In addition to the orchestrated activation of genetic programs, spontaneous electrical activity and associated calcium signaling have been shown to be critically involved in the maturation of such neuronal networks. We established an in vitro model that recapitulates the maturation of neuronal networks, including spontaneous electrical activity. Upon plating, mouse primary hippocampal neurons grow neurites and interconnect via synapses to form a dish-wide neuronal network. Via live cell calcium imaging, we identified a limited period of time in which the spontaneous activity synchronizes across neurons, indicative of the formation of a functional network. After establishment of network activity, the neurons grow dendritic spines, the density of which was used as a morphological readout for neuronal maturity and connectivity. Hence, quantification of neurite outgrowth, synapse density, spontaneous neuronal activity, and dendritic spine density allowed to study neuronal network maturation from the day of plating until the presence of mature neuronal networks. Via acute pharmacological intervention, we show that synchronized network activity is mediated by the NMDA-R. The balance between kynurenic and quinolinic acid, both neuro-active intermediates in the tryptophan/kynurenine pathway, was shown to be decisive for the maintenance of network activity. Chronic modulation of the neurotrophic support influenced the network formation and revealed the extreme sensitivity of calcium imaging to detect subtle alterations in neuronal physiology. Given the reproducible cultivation in a 96-well setup in combination with fully automated analysis of the calcium recordings, this approach can be used to build a high-content screening assay usable for neurotoxicity screening, target identification/validation, or phenotypic drug screening.     

Beyond Two-Cell Networks: Experimental Measurement of Neuronal Responses to Multiple Synaptic Inputs

Oscillations of large populations of neurons are thought to be important in the normal functioning of the brain. We have used phase response curve (PRC) methods to characterize the dynamics of single neurons and predict population dynamics. Our past experimental work was limited to special circumstances (e.g., 2-cell networks of periodically firing neurons). Here, we explore the feasibility of extending our methods to predict the synchronization properties of stellate cells (SCs) in the rat entorhinal cortex under broader conditions. In particular, we test the hypothesis that PRCs in SCs scale linearly with changes in synaptic amplitude, and measure how well responses to Poisson process-driven inputs can be predicted in terms of PRCs. Although we see nonlinear responses to excitatory and inhibitory inputs, we find that models based on weak coupling account for scaling and Poisson process-driven inputs reasonably accurately.     

Cancer Development and Progression: A Non-Adaptive Process Driven by Genetic Drift

The current mainstream in cancer research favours the idea that malignant tumour initiation is the result of a genetic mutation. Tumour development and progression is then explained as a sort of micro-evolutionary process, whereby an initial genetic alteration leads to abnormal proliferation of a single cell that leads to a population of clonally derived cells. It is widely claimed that tumour progression is driven by natural selection, based on the assumption that the initial tumour cells acquire some properties that endow such cells with a selective advantage over the normal cells from which the tumour cells are derived. The standard view assumes that the transformed bodily cell somehow acquires "responsiveness" to natural selection independently of the whole organism to which the cell belongs. Yet, it is never explained where such an acquired capacity to respond to natural selection by the individual bodily cell comes from. This situation poses many difficult questions that so far have been left unanswered. For example, there is no explanation why some cells belonging to an organised whole and as such having no independent capacity for survival, apparently become 'independent' entities, able to respond to selective pressures in an autonomous fashion and then to be evaluated by natural selection. Hereunder it is argued that such a qualitative change cannot be the consequence of specific genetic mutations. Moreover, it is shown that natural selection is unlikely to be acting within the organism during tumour development and progression and that tumour evolution is a random, non-adaptive process, driven by no fundamental biological principle. Thus, mutations in the so-called oncogenes and tumour suppressor genes observed in epithelial cancers (that constitute more than 90% of all cancers) are not the result of selection for better cellular growth or survival under restrictive conditions. Instead, here it is suggested that they are the consequence of genetic drift acting upon gene functions that become non-relevant, either for the individual or the species fitness, once the organism is past its reproductive prime and as such, they also become superfluous for cell survival in the short term. It is proposed that the origin of cancer is epigenetic and it is a consequence of the need for a continued turnover of the individuals that constitute a species.