Power-law statistics and universal scaling in the absence of criticality

Critical states are sometimes identified experimentally through power-law statistics or universal scaling functions. We show here that such features naturally emerge from networks in self-sustained irregular regimes away from criticality. In these regimes, statistical physics theory of large interacting systems predict a regime where the nodes have independent and identically distributed dynamics. We thus investigated the statistics of a system in which units are replaced by independent stochastic surrogates, and found the same power-law statistics, indicating that these are not sufficient to establish criticality. We rather suggest that these are universal features of large-scale networks when considered macroscopically. These results put caution on the interpretation of scaling laws found in nature.Comment: in press in Phys. Rev.

Estimating input parameters from intracellular recordings in the Feller neuronal model

We study the estimation of the input parameters in a Feller neuronal model from a trajectory of the membrane potential sampled at discrete times. These input parameters are identified with the drift and the infinitesimal variance of the underlying stochastic diffusion process with multiplicative noise. The state space of the process is restricted from below by an inaccessible boundary. Further, the model is characterized by the presence of an absorbing threshold, the first hitting of which determines the length of each trajectory and which constrains the state space from above. We compare, both in the presence and in the absence of the absorbing threshold, the efficiency of different known estimators. In addition, we propose an estimator for the drift term, which is proved to be more efficient than the others, at least in the explored range of the parameters. The presence of the threshold makes the estimates of the drift term biased, and two methods to correct it are proposed

Sample Path Analysis of Integrate-and-Fire Neurons

Computational neuroscience is concerned with answering two intertwined questions that are based on the assumption that spatio-temporal patterns of spikes form the universal language of the nervous system. First, what function does a specific neural circuitry perform in the elaboration of a behavior? Second, how do neural circuits process behaviorally-relevant information? Non-linear system analysis has proven instrumental in understanding the coding strategies of early neural processing in various sensory modalities. Yet, at higher levels of integration, it fails to help in deciphering the response of assemblies of neurons to complex naturalistic stimuli. If neural activity can be assumed to be primarily driven by the stimulus at early stages of processing, the intrinsic activity of neural circuits interacts with their high-dimensional input to transform it in a stochastic non-linear fashion at the cortical level. As a consequence, any attempt to fully understand the brain through a system analysis approach becomes illusory. However, it is increasingly advocated that neural noise plays a constructive role in neural processing, facilitating information transmission. This prompts to gain insight into the neural code by studying the stochasticity of neuronal activity, which is viewed as biologically relevant. Such an endeavor requires the design of guiding theoretical principles to assess the potential benefits of neural noise. In this context, meeting the requirements of biological relevance and computational tractability, while providing a stochastic description of neural activity, prescribes the adoption of the integrate-and-fire model. In this thesis, founding ourselves on the path-wise description of neuronal activity, we propose to further the stochastic analysis of the integrate-and fire model through a combination of numerical and theoretical techniques. To begin, we expand upon the path-wise construction of linear diffusions, which offers a natural setting to describe leaky integrate-and-fire neurons, as inhomogeneous Markov chains. Based on the theoretical analysis of the first-passage problem, we then explore the interplay between the internal neuronal noise and the statistics of injected perturbations at the single unit level, and examine its implications on the neural coding. At the population level, we also develop an exact event-driven implementation of a Markov network of perfect integrate-and-fire neurons with both time delayed instantaneous interactions and arbitrary topology. We hope our approach will provide new paradigms to understand how sensory inputs perturb neural intrinsic activity and accomplish the goal of developing a new technique for identifying relevant patterns of population activity. From a perturbative perspective, our study shows how injecting frozen noise in different flavors can help characterize internal neuronal noise, which is presumably functionally relevant to information processing. From a simulation perspective, our event-driven framework is amenable to scrutinize the stochastic behavior of simple recurrent motifs as well as temporal dynamics of large scale networks under spike-timing-dependent plasticity

Informační procesy v neuronech

Neurony spolu komunikují pomocí posloupností akčních potenciálů. Celý tento proces může být popsán detailními biochemickými modely membrány a iontových kanálů na neuronu nebo jednoduššími fenomenologickými mo- dely (typickým představitelem jsou tzv. "integrate-and-fire" modely) nebo případně ještě více abstraktními modely sledu akčních potenciálů bez při- hlédnutí k dynamice membrány neuronu. Vybrali jsme konkrétní variantu stochastického "leaky integrate-and-fire" modelu a porovnali jí s aktivitou biologického neuronu (nitrobuněčný zá- znam pořízený in-vivo). Provedli jsme statistický odhad parametrů modelu a na základě počítačových simulací úspěšně srovnali modelovaný záznam se záznamem z reálného neuronu. Při abstraktnější úrovni popisu je sled akčních potenciálů analyzován pouze jako množina bodových událostí v čase a základní otázka zní, jakým způsobem je vnější podnět kódován v zaznamenané posloupnosti akčních po- tenciálů. Bylo navrženo mnoho odlišných kódů pro řešení rozmanitých úloh v neuronových sítích. My jsme se zaměřili na otevřený problém neuronál- ního kódu v úloze prostorového slyšení u savců. V současnosti je zvažováno několik teorií vysvětlujících experimentální nálezy. V naší práci navrhujeme specifickou variantu modelu založeného na frekvenčním kódu. Zkonstruovaný neuronový obvod,...Neurons communicate by action potentials. This process can be described by very detailed biochemical models of neuronal membrane and its channels, or by simpler phenomenological models of membrane potential (integrate-and- fire models) or even by very abstract models when only time of spikes are considered. We took one particular description - stochastic leaky integrate-and-fire model - and compared it with recorded in-vivo intracellular activity of the neuron. We estimated parameters of this model, compared how the model simulation corresponds with a real neuron. It can be concluded that the data are generally consistent with the model. At a more abstract level of description, the spike trains are analyzed without considering exact membrane voltage and one asks how the external stimulus is encoded in the spike train emitted by neurons. There are many neuronal codes described in literature and we focused on the open problem of neural code responsible for spatial hearing in mammals. Several theories explaining the experimental findings have been proposed and we suggest a specific variant of so called slope-encoding model. Neuronal circuit mimick- ing auditory pathway up to the first binaural neuron was constructed and experimental results were reproduced. Finally, we estimated the minimal number of such...First Faculty of Medicine1. lékařská fakult

Coupled stochastic systems of Skorokhod type: well-posedness of a mathematical model and its applications

Population dynamics with complex biological interactions, accounting for uncertainty quantification, is critical for many application areas. However, due to the complexity of biological systems, the mathematical formulation of the corresponding problems faces the challenge that the corresponding stochastic processes should, in most cases, be considered in bounded domains. We propose a model based on a coupled system of reflecting Skorokhod-type stochastic differential equations with jump-like exit from a boundary. The setting describes the population dynamics of active and passive populations. As main working techniques, we use compactness methods and Skorokhod’s representation of solutions to SDEs posed in bounded domains to prove the well-posedness of the system. This functional setting is a new point of view in the field of modelling and simulation of population dynamics. We provide the details of the model, as well as representative numerical examples, and discuss the applications of a Wilson-Cowan-type system, modelling the dynamics of two interacting populations of excitatory and inhibitory neurons. Furthermore, the presence of random input current, reflecting factors together with Poisson jumps, increases firing activity in neuronal systems

A statistical model for in vivo neuronal dynamics

Single neuron models have a long tradition in computational neuroscience. Detailed biophysical models such as the Hodgkin-Huxley model as well as simplified neuron models such as the class of integrate-and-fire models relate the input current to the membrane potential of the neuron. Those types of models have been extensively fitted to in vitro data where the input current is controlled. Those models are however of little use when it comes to characterize intracellular in vivo recordings since the input to the neuron is not known. Here we propose a novel single neuron model that characterizes the statistical properties of in vivo recordings. More specifically, we propose a stochastic process where the subthreshold membrane potential follows a Gaussian process and the spike emission intensity depends nonlinearly on the membrane potential as well as the spiking history. We first show that the model has a rich dynamical repertoire since it can capture arbitrary subthreshold autocovariance functions, firing-rate adaptations as well as arbitrary shapes of the action potential. We then show that this model can be efficiently fitted to data without overfitting. Finally, we show that this model can be used to characterize and therefore precisely compare various intracellular in vivo recordings from different animals and experimental conditions.Comment: 31 pages, 10 figure

Contributions to statistical analysis methods for neural spiking activity

With the technical advances in neuroscience experiments in the past few decades, we have seen a massive expansion in our ability to record neural activity. These advances enable neuroscientists to analyze more complex neural coding and communication properties, and at the same time, raise new challenges for analyzing neural spiking data, which keeps growing in scale, dimension, and complexity. This thesis proposes several new statistical methods that advance statistical analysis approaches for neural spiking data, including sequential Monte Carlo (SMC) methods for efficient estimation of neural dynamics from membrane potential threshold crossings, state-space models using multimodal observation processes, and goodness-of-fit analysis methods for neural marked point process models. In a first project, we derive a set of iterative formulas that enable us to simulate trajectories from stochastic, dynamic neural spiking models that are consistent with a set of spike time observations. We develop a SMC method to simultaneously estimate the parameters of the model and the unobserved dynamic variables from spike train data. We investigate the performance of this approach on a leaky integrate-and-fire model. In another project, we define a semi-latent state-space model to estimate information related to the phenomenon of hippocampal replay. Replay is a recently discovered phenomenon where patterns of hippocampal spiking activity that typically occur during exploration of an environment are reactivated when an animal is at rest. This reactivation is accompanied by high frequency oscillations in hippocampal local field potentials. However, methods to define replay mathematically remain undeveloped. In this project, we construct a novel state-space model that enables us to identify whether replay is occurring, and if so to estimate the movement trajectories consistent with the observed neural activity, and to categorize the content of each event. The state-space model integrates information from the spiking activity from the hippocampal population, the rhythms in the local field potential, and the rat's movement behavior. Finally, we develop a new, general time-rescaling theorem for marked point processes, and use this to develop a general goodness-of-fit framework for neural population spiking models. We investigate this approach through simulation and a real data application