No Arabic abstract
Gamma frequency oscillations (25-140 Hz), observed in the neural activities within many brain regions, have long been regarded as a physiological basis underlying many brain functions, such as memory and attention. Among numerous theoretical and computational modeling studies, gamma oscillations have been found in biologically realistic spiking network models of the primary visual cortex. However, due to its high dimensionality and strong nonlinearity, it is generally difficult to perform detailed theoretical analysis of the emergent gamma dynamics. Here we propose a suite of Markovian model reduction methods with varying levels of complexity and applied it to spiking network models exhibiting heterogeneous dynamical regimes, ranging from homogeneous firing to strong synchrony in the gamma band. The reduced models not only successfully reproduce gamma band oscillations in the full model, but also exhibit the same dynamical features as we vary parameters. Most remarkably, the invariant measure of the coarse-grained Markov process reveals a two-dimensional surface in state space upon which the gamma dynamics mainly resides. Our results suggest that the statistical features of gamma oscillations strongly depend on the subthreshold neuronal distributions. Because of the generality of the Markovian assumptions, our dimensional reduction methods offer a powerful toolbox for theoretical examinations of many other complex cortical spatio-temporal behaviors observed in both neurophysiological experiments and numerical simulations.
The random transitions of ion channels between conducting and non-conducting states generate a source of internal fluctuations in a neuron, known as channel noise. The standard method for modeling fluctuations in the states of ion channels uses continuous-time Markov chains nonlinearly coupled to a differential equation for voltage. Beginning with the work of Fox and Lu, there have been attempts to generate simpler models that use stochastic differential equation (SDEs) to approximate the stochastic spiking activity produced by Markov chain models. Recent numerical investigations, however, have raised doubts that SDE models can preserve the stochastic dynamics of Markov chain models. We analyze three SDE models that have been proposed as approximations to the Markov chain model: one that describes the states of the ion channels and two that describe the states of the ion channel subunits. We show that the former channel-based approach can capture the distribution of channel noise and its effect on spiking in a Hodgkin-Huxley neuron model to a degree not previously demonstrated, but the latter two subunit-based approaches cannot. Our analysis provides intuitive and mathematical explanations for why this is the case: the temporal correlation in the channel noise is determined by the combinatorics of bundling subunits into channels, and the subunit-based approaches do not correctly account for this structure. Our study therefore confirms and elucidates the findings of previous numerical investigations of subunit-based SDE models. Moreover, it presents the first evidence that Markov chain models of the nonlinear, stochastic dynamics of neural membranes can be accurately approximated by SDEs. This finding opens a door to future modeling work using SDE techniques to further illuminate the effects of ion channel fluctuations on electrically active cells.
The oxytocin effects on large-scale brain networks such as Default Mode Network (DMN) and Frontoparietal Network (FPN) have been largely studied using fMRI data. However, these studies are mainly based on the statistical correlation or Bayesian causality inference, lacking interpretability at physical and neuroscience level. Here, we propose a physics-based framework of Kuramoto model to investigate oxytocin effects on the phase dynamic neural coupling in DMN and FPN. Testing on fMRI data of 59 participants administrated with either oxytocin or placebo, we demonstrate that oxytocin changes the topology of brain communities in DMN and FPN, leading to higher synchronization in the FPN and lower synchronization in the DMN, as well as a higher variance of the coupling strength within the DMN and more flexible coupling patterns across time. These results together indicate that oxytocin may increase the ability to overcome the corresponding internal oscillation dispersion and support the flexibility in neural synchrony in various social contexts, providing new evidence for explaining the oxytocin modulated social behaviors. Our proposed Kuramoto model-based framework can be a potential tool in network neuroscience and offers physical and neural insights into phase dynamics of the brain.
We consider a pair of stochastic integrate and fire neurons receiving correlated stochastic inputs. The evolution of this system can be described by the corresponding Fokker-Planck equation with non-trivial boundary conditions resulting from the refractory period and firing threshold. We propose a finite volume method that is orders of magnitude faster than the Monte Carlo methods traditionally used to model such systems. The resulting numerical approximations are proved to be accurate, nonnegative and integrate to 1. We also approximate the transient evolution of the system using an Ornstein--Uhlenbeck process, and use the result to examine the properties of the joint output of cell pairs. The results suggests that the joint output of a cell pair is most sensitive to changes in input variance, and less sensitive to changes in input mean and correlation.
In this paper we propose an $infty-$dimensional cerebellar model of neural controller for realistic human biodynamics. The model is developed using Feynmans action-amplitude (partition function) formalism. The cerebellum controller is acting as a supervisor for an autogenetic servo control of human musculo-skeletal dynamics, which is presented in (dissipative, driven) Hamiltonian form. The $infty-$dimensional cerebellar controller is closely related to entropic motor control. Keywords: realistic human biodynamics, cerebellum motion control, $infty-$dimensional neural network
Anaerobic glycolysis in yeast perturbed by the reduction of xenobiotic ketones is studied numerically in two models which possess the same topology but different levels of complexity. By comparing both models predictions for concentrations and fluxes as well as steady or oscillatory temporal behavior we answer the question what phenomena require what kind of minimum model abstraction. While mean concentrations and fluxes are predicted in agreement by both models we observe different domains of oscillatory behavior in parameter space. Generic properties of the glycolytic response to ketones are discussed.