The Cauchy problem for one-dimensional spiking neuron models.
Ontology highlight
ABSTRACT:
SUBMITTER: Brette R
PROVIDER: S-EPMC2289251 | biostudies-other | 2008 Mar
REPOSITORIES: biostudies-other
Similar Datasets
Project description:An increasing amount of behavioral and neurophysiological data suggests that the brain performs optimal (or near-optimal) probabilistic inference and learning during perception and other tasks. Although many machine learning algorithms exist that perform inference and learning in an optimal way, the complete description of how one of those algorithms (or a novel algorithm) can be implemented in the brain is currently incomplete. There have been many proposed solutions that address how neurons can perform optimal inference but the question of how synaptic plasticity can implement optimal learning is rarely addressed. This paper aims to unify the two fields of probabilistic inference and synaptic plasticity by using a neuronal network of realistic model spiking neurons to implement a well-studied computational model called the Helmholtz Machine. The Helmholtz Machine is amenable to neural implementation as the algorithm it uses to learn its parameters, called the wake-sleep algorithm, uses a local delta learning rule. Our spiking-neuron network implements both the delta rule and a small example of a Helmholtz machine. This neuronal network can learn an internal model of continuous-valued training data sets without supervision. The network can also perform inference on the learned internal models. We show how various biophysical features of the neural implementation constrain the parameters of the wake-sleep algorithm, such as the duration of the wake and sleep phases of learning and the minimal sample duration. We examine the deviations from optimal performance and tie them to the properties of the synaptic plasticity rule.
| S-EPMC4405618 | biostudies-literature
Project description:Realistic modeling of neurons are quite successful in complementing traditional experimental techniques. However, their networks require a computational power beyond the capabilities of current supercomputers, and the methods used so far to reduce their complexity do not take into account the key features of the cells nor critical physiological properties. Here we introduce a new, automatic and fast method to map realistic neurons into equivalent reduced models running up to > 40 times faster while maintaining a very high accuracy of the membrane potential dynamics during synaptic inputs, and a direct link with experimental observables. The mapping of arbitrary sets of synaptic inputs, without additional fine tuning, would also allow the convenient and efficient implementation of a new generation of large-scale simulations of brain regions reproducing the biological variability observed in real neurons, with unprecedented advances to understand higher brain functions.
| S-EPMC3514644 | biostudies-literature
Project description:Point process generalized linear models (PP-GLMs) provide an important statistical framework for modeling spiking activity in single-neurons and neuronal networks. Stochastic stability is essential when sampling from these models, as done in computational neuroscience to analyze statistical properties of neuronal dynamics and in neuro-engineering to implement closed-loop applications. Here we show, however, that despite passing common goodness-of-fit tests, PP-GLMs estimated from data are often unstable, leading to divergent firing rates. The inclusion of absolute refractory periods is not a satisfactory solution since the activity then typically settles into unphysiological rates. To address these issues, we derive a framework for determining the existence and stability of fixed points of the expected conditional intensity function (CIF) for general PP-GLMs. Specifically, in nonlinear Hawkes PP-GLMs, the CIF is expressed as a function of the previous spike history and exogenous inputs. We use a mean-field quasi-renewal (QR) approximation that decomposes spike history effects into the contribution of the last spike and an average of the CIF over all spike histories prior to the last spike. Fixed points for stationary rates are derived as self-consistent solutions of integral equations. Bifurcation analysis and the number of fixed points predict that the original models can show stable, divergent, and metastable (fragile) dynamics. For fragile models, fluctuations of the single-neuron dynamics predict expected divergence times after which rates approach unphysiologically high values. This metric can be used to estimate the probability of rates to remain physiological for given time periods, e.g., for simulation purposes. We demonstrate the use of the stability framework using simulated single-neuron examples and neurophysiological recordings. Finally, we show how to adapt PP-GLM estimation procedures to guarantee model stability. Overall, our results provide a stability framework for data-driven PP-GLMs and shed new light on the stochastic dynamics of state-of-the-art statistical models of neuronal spiking activity.
| S-EPMC5325182 | biostudies-literature
Project description:We present a computational multiscale model for the efficient simulation of vascularized tissues, composed of an elastic three-dimensional matrix and a vascular network. The effect of blood vessel pressure on the elastic tissue is surrogated via hyper-singular forcing terms in the elasticity equations, which depend on the fluid pressure. In turn, the blood flow in vessels is treated as a one-dimensional network. Intravascular pressure and velocity are simulated using a high-order finite volume scheme, while the elasticity equations for the tissue are solved using a finite element method. This work addresses the feasibility and the potential of the proposed coupled multiscale model. In particular, we assess whether the multiscale model is able to reproduce the tissue response at the effective scale (of the order of millimeters) while modeling the vasculature at the microscale. We validate the multiscale method against a full scale (three-dimensional) model, where the fluid/tissue interface is fully discretized and treated as a Neumann boundary for the elasticity equation. Next, we present simulation results obtained with the proposed approach in a realistic scenario, demonstrating that the method can robustly and efficiently handle the one-way coupling between complex fluid microstructures and the elastic matrix.
| S-EPMC8671283 | biostudies-literature
Project description:As computational models of the cardiovascular system are applied in modern personalized medicine, maximizing certainty of model input becomes crucial. A model with a high number of arterial segments results in a more realistic description of the system, but also requires a high number of parameters with associated uncertainties. In this paper, we present a method to optimize/reduce the number of arterial segments included in one-dimensional blood flow models, while preserving key features of flow and pressure waveforms. We quantify the preservation of key flow features for the optimal network with respect to the baseline networks (a 96-artery and a patient-specific coronary network) by various metrics and quantities like average relative error, pulse pressure and augmentation pressure. Furthermore, various physiological and pathological states are considered. For the aortic root and larger systemic artery pressure waveforms a network with minimal description of lower and upper limb arteries and no cerebral arteries, sufficiently captures important features such as pressure augmentation and pulse pressure. Discrepancies in carotid and middle cerebral artery flow waveforms that are introduced by describing the arterial system in a minimalistic manner are small compared with errors related to uncertainties in blood flow measurements obtained by ultrasound.
| S-EPMC6303799 | biostudies-literature
Project description:Low-dimensional attractive manifolds with flows prescribing the evolution of state variables are commonly used to capture the lawful behavior of behavioral and cognitive variables. Neural network dynamics underlie many of the mechanistic explanations of function and demonstrate the existence of such low-dimensional attractive manifolds. In this study, we focus on exploring the network mechanisms due to asymmetric couplings giving rise to the emergence of arbitrary flows in low dimensional spaces. Here we use a spiking neural network model, specifically the theta neuron model and simple synaptic dynamics, to show how a qualitatively identical set of basic behaviors arises from different combinations of couplings with broken symmetry, in fluctuations of both firing rate and spike timing. We further demonstrate how such network dynamics can be combined to create more complex processes. These results suggest that 1) asymmetric coupling is not always a variance to be averaged over, 2) different networks may produce the same dynamics by different dynamical routes and 3) complex dynamics may be formed by simpler dynamics through a combination of couplings.
| S-EPMC3656844 | biostudies-other
Project description:Different global patterns of brain activity are associated with distinct arousal and behavioral states of an animal, but how the brain rapidly switches between different states remains unclear. We here report that repetitive high-frequency burst spiking of a single rat cortical neuron could trigger a switch between the cortical states resembling slow-wave and rapid-eye-movement sleep. This is reflected in the switching of the membrane potential of the stimulated neuron from slow UP/DOWN oscillations to a persistent-UP state or vice versa, with concurrent changes in the temporal pattern of cortical local field potential (LFP) recorded several millimeters away. These results point to the power of single cortical neurons in modulating the behavioral state of an animal.
| S-EPMC2913066 | biostudies-literature
Project description:An unusual feature of the cerebellar cortex is that its output neurons, Purkinje cells, release GABA (?-aminobutyric acid). Their high intrinsic firing rates (50?Hz) and extensive convergence predict that their target neurons in the cerebellar nuclei would be largely inhibited unless Purkinje cells pause their spiking, yet Purkinje and nuclear neuron firing rates do not always vary inversely. One indication of how these synapses transmit information is that populations of Purkinje neurons synchronize their spikes during cerebellar behaviours. If nuclear neurons respond to Purkinje synchrony, they may encode signals from subsets of inhibitory inputs. Here we show in weanling and adult mice that nuclear neurons transmit the timing of synchronous Purkinje afferent spikes, owing to modest Purkinje-to-nuclear convergence ratios (?40:1), fast inhibitory postsynaptic current kinetics (?(decay) = 2.5?ms) and high intrinsic firing rates (?90?Hz). In vitro, dynamically clamped asynchronous inhibitory postsynaptic potentials mimicking Purkinje afferents suppress nuclear cell spiking, whereas synchronous inhibitory postsynaptic potentials entrain nuclear cell spiking. With partial synchrony, nuclear neurons time-lock their spikes to the synchronous subpopulation of inputs, even when only 2 out of 40 afferents synchronize. In vivo, nuclear neurons reliably phase-lock to regular trains of molecular layer stimulation. Thus, cerebellar nuclear neurons can preferentially relay the spike timing of synchronized Purkinje cells to downstream premotor areas.
| S-EPMC3268051 | biostudies-literature
Project description:Insects are able to solve basic numerical cognition tasks. We show that estimation of numerosity can be realized and learned by a single spiking neuron with an appropriate synaptic plasticity rule. This model can be efficiently trained to detect arbitrary spatiotemporal spike patterns on a noisy and dynamic background with high precision and low variance. When put to test in a task that requires counting of visual concepts in a static image it required considerably less training epochs than a convolutional neural network to achieve equal performance. When mimicking a behavioral task in free-flying bees that requires numerical cognition, the model reaches a similar success rate in making correct decisions. We propose that using action potentials to represent basic numerical concepts with a single spiking neuron is beneficial for organisms with small brains and limited neuronal resources.
| S-EPMC7005464 | biostudies-literature
Project description:We propose a novel biologically plausible computational model of working memory (WM) implemented by a spiking neuron network (SNN) interacting with a network of astrocytes. The SNN is modeled by synaptically coupled Izhikevich neurons with a non-specific architecture connection topology. Astrocytes generating calcium signals are connected by local gap junction diffusive couplings and interact with neurons via chemicals diffused in the extracellular space. Calcium elevations occur in response to the increased concentration of the neurotransmitter released by spiking neurons when a group of them fire coherently. In turn, gliotransmitters are released by activated astrocytes modulating the strength of the synaptic connections in the corresponding neuronal group. Input information is encoded as two-dimensional patterns of short applied current pulses stimulating neurons. The output is taken from frequencies of transient discharges of corresponding neurons. We show how a set of information patterns with quite significant overlapping areas can be uploaded into the neuron-astrocyte network and stored for several seconds. Information retrieval is organized by the application of a cue pattern representing one from the memory set distorted by noise. We found that successful retrieval with the level of the correlation between the recalled pattern and ideal pattern exceeding 90% is possible for the multi-item WM task. Having analyzed the dynamical mechanism of WM formation, we discovered that astrocytes operating at a time scale of a dozen of seconds can successfully store traces of neuronal activations corresponding to information patterns. In the retrieval stage, the astrocytic network selectively modulates synaptic connections in the SNN leading to successful recall. Information and dynamical characteristics of the proposed WM model agrees with classical concepts and other WM models.
| S-EPMC8044545 | biostudies-literature