Emergent complex neural dynamics

A large repertoire of spatiotemporal activity patterns in the brain is the basis for adaptive behaviour. Understanding the mechanism by which the brain's hundred billion neurons and hundred trillion synapses manage to produce such a range of cortical configurations in a flexible manner remains a fundamental problem in neuroscience. One plausible solution is the involvement of universal mechanisms of emergent complex phenomena evident in dynamical systems poised near a critical point of a second-order phase transition. We review recent theoretical and empirical results supporting the notion that the brain is naturally poised near criticality, as well as its implications for better understanding of the brain

Suppression of growth by multiplicative white noise in a parametric resonant system

The author studied the growth of the amplitude in a Mathieu-like equation with multiplicative white noise. The approximate value of the exponent at the extremum on parametric resonance regions was obtained theoretically by introducing the width of time interval, and the exponents were calculated numerically by solving the stochastic differential equations by a symplectic numerical method. The Mathieu-like equation contains a parameter $\alpha$ that is determined by the intensity of noise and the strength of the coupling between the variable and the noise. The value of $\alpha$ was restricted not to be negative without loss of generality. It was shown that the exponent decreases with $\alpha$, reaches a minimum and increases after that. It was also found that the exponent as a function of $\alpha$ has only one minimum at $\alpha \neq 0$ on parametric resonance regions of $\alpha = 0$. This minimum value is obtained theoretically and numerically. The existence of the minimum at $\alpha \neq 0$ indicates the suppression of the growth by multiplicative white noise.Comment: The title and the description in the manuscript are change

Signal Propagation in Feedforward Neuronal Networks with Unreliable Synapses

In this paper, we systematically investigate both the synfire propagation and firing rate propagation in feedforward neuronal network coupled in an all-to-all fashion. In contrast to most earlier work, where only reliable synaptic connections are considered, we mainly examine the effects of unreliable synapses on both types of neural activity propagation in this work. We first study networks composed of purely excitatory neurons. Our results show that both the successful transmission probability and excitatory synaptic strength largely influence the propagation of these two types of neural activities, and better tuning of these synaptic parameters makes the considered network support stable signal propagation. It is also found that noise has significant but different impacts on these two types of propagation. The additive Gaussian white noise has the tendency to reduce the precision of the synfire activity, whereas noise with appropriate intensity can enhance the performance of firing rate propagation. Further simulations indicate that the propagation dynamics of the considered neuronal network is not simply determined by the average amount of received neurotransmitter for each neuron in a time instant, but also largely influenced by the stochastic effect of neurotransmitter release. Second, we compare our results with those obtained in corresponding feedforward neuronal networks connected with reliable synapses but in a random coupling fashion. We confirm that some differences can be observed in these two different feedforward neuronal network models. Finally, we study the signal propagation in feedforward neuronal networks consisting of both excitatory and inhibitory neurons, and demonstrate that inhibition also plays an important role in signal propagation in the considered networks.Comment: 33pages, 16 figures; Journal of Computational Neuroscience (published

Uncovering the Dynamics of Cardiac Systems Using Stochastic Pacing and Frequency Domain Analyses

Alternans of cardiac action potential duration (APD) is a well-known arrhythmogenic mechanism which results from dynamical instabilities. The propensity to alternans is classically investigated by examining APD restitution and by deriving APD restitution slopes as predictive markers. However, experiments have shown that such markers are not always accurate for the prediction of alternans. Using a mathematical ventricular cell model known to exhibit unstable dynamics of both membrane potential and Ca2+ cycling, we demonstrate that an accurate marker can be obtained by pacing at cycle lengths (CLs) varying randomly around a basic CL (BCL) and by evaluating the transfer function between the time series of CLs and APDs using an autoregressive-moving-average (ARMA) model. The first pole of this transfer function corresponds to the eigenvalue (λalt) of the dominant eigenmode of the cardiac system, which predicts that alternans occurs when λalt≤−1. For different BCLs, control values of λalt were obtained using eigenmode analysis and compared to the first pole of the transfer function estimated using ARMA model fitting in simulations of random pacing protocols. In all versions of the cell model, this pole provided an accurate estimation of λalt. Furthermore, during slow ramp decreases of BCL or simulated drug application, this approach predicted the onset of alternans by extrapolating the time course of the estimated λalt. In conclusion, stochastic pacing and ARMA model identification represents a novel approach to predict alternans without making any assumptions about its ionic mechanisms. It should therefore be applicable experimentally for any type of myocardial cell

A Simple Artificial Life Model Explains Irrational Behavior in Human Decision-Making

Although praised for their rationality, humans often make poor decisions, even in simple situations. In the repeated binary choice experiment, an individual has to choose repeatedly between the same two alternatives, where a reward is assigned to one of them with fixed probability. The optimal strategy is to perseverate with choosing the alternative with the best expected return. Whereas many species perseverate, humans tend to match the frequencies of their choices to the frequencies of the alternatives, a sub-optimal strategy known as probability matching. Our goal was to find the primary cognitive constraints under which a set of simple evolutionary rules can lead to such contrasting behaviors. We simulated the evolution of artificial populations, wherein the fitness of each animat (artificial animal) depended on its ability to predict the next element of a sequence made up of a repeating binary string of varying size. When the string was short relative to the animats’ neural capacity, they could learn it and correctly predict the next element of the sequence. When it was long, they could not learn it, turning to the next best option: to perseverate. Animats from the last generation then performed the task of predicting the next element of a non-periodical binary sequence. We found that, whereas animats with smaller neural capacity kept perseverating with the best alternative as before, animats with larger neural capacity, which had previously been able to learn the pattern of repeating strings, adopted probability matching, being outperformed by the perseverating animats. Our results demonstrate how the ability to make predictions in an environment endowed with regular patterns may lead to probability matching under less structured conditions. They point to probability matching as a likely by-product of adaptive cognitive strategies that were crucial in human evolution, but may lead to sub-optimal performances in other environments

Mitochondrial chaotic dynamics: Redox-energetic behavior at the edge of stability

Mitochondria serve multiple key cellular functions, including energy generation, redox balance, and regulation of apoptotic cell death, thus making a major impact on healthy and diseased states. Increasingly recognized is that biological network stability/instability can play critical roles in determining health and disease. We report for the first-time mitochondrial chaotic dynamics, characterizing the conditions leading from stability to chaos in this organelle. Using an experimentally validated computational model of mitochondrial function, we show that complex oscillatory dynamics in key metabolic variables, arising at the “edge” between fully functional and pathological behavior, sets the stage for chaos. Under these conditions, a mild, regular sinusoidal redox forcing perturbation triggers chaotic dynamics with main signature traits such as sensitivity to initial conditions, positive Lyapunov exponents, and strange attractors. At the “edge” mitochondrial chaos is exquisitely sensitive to the antioxidant capacity of matrix Mn superoxide dismutase as well as to the amplitude and frequency of the redox perturbation. These results have potential implications both for mitochondrial signaling determining health maintenance, and pathological transformation, including abnormal cardiac rhythms.publishedVersionKembro, Jackelyn Melissa. Universidad Nacional de Córdoba. Facultad de Ciencias Exactas, Físicas y Naturales; Argentina.Kembro, Jackelyn Melissa. Consejo Nacional de Investigaciones Científicas y Técnicas. Instituto de Investigaciones Biológicas y Tecnológicas; Argentina.Cortassa, Sonia. National Institutes of Health. NIH · NIA Intramural Research Program; Estados Unidos.Lloyd, David. Cardiff University. School of Biosciences 1; Inglaterra.Sollot, Steven. Johns Hopkins University. Laboratory of Cardiovascular Science; Estados Unidos.Sollot, Steven. Johns Hopkins University. Laboratory of Cardiovascular Science; Estados Unidos

Critical synchronization dynamics of the Kuramoto model on connectome and small world graphs

The hypothesis, that cortical dynamics operates near criticality also suggests, that it exhibits universal critical exponents which marks the Kuramoto equation, a fundamental model for synchronization, as a prime candidate for an underlying universal model. Here, we determined the synchronization behavior of this model by solving it numerically on a large, weighted human connectome network, containing 804092 nodes, in an assumed homeostatic state. Since this graph has a topological dimension $d < 4$, a real synchronization phase transition is not possible in the thermodynamic limit, still we could locate a transition between partially synchronized and desynchronized states. At this crossover point we observe power-law--tailed synchronization durations, with $\tau_t \simeq 1.2(1)$, away from experimental values for the brain. For comparison, on a large two-dimensional lattice, having additional random, long-range links, we obtain a mean-field value: $\tau_t \simeq 1.6(1)$. However, below the transition of the connectome we found global coupling control-parameter dependent exponents $1 < \tau_t \le 2$, overlapping with the range of human brain experiments. We also studied the effects of random flipping of a small portion of link weights, mimicking a network with inhibitory interactions, and found similar results. The control-parameter dependent exponent suggests extended dynamical criticality below the transition point.Comment: 12 pages, 9 figures + Supplemenraty material pdf 2 pages 4 figs, 1 table, accepted version in Scientific Report

Failure of adaptive self-organized criticality during epileptic seizure attacks

Critical dynamics are assumed to be an attractive mode for normal brain functioning as information processing and computational capabilities are found to be optimized there. Recent experimental observations of neuronal activity patterns following power-law distributions, a hallmark of systems at a critical state, have led to the hypothesis that human brain dynamics could be poised at a phase transition between ordered and disordered activity. A so far unresolved question concerns the medical significance of critical brain activity and how it relates to pathological conditions. Using data from invasive electroencephalogram recordings from humans we show that during epileptic seizure attacks neuronal activity patterns deviate from the normally observed power-law distribution characterizing critical dynamics. The comparison of these observations to results from a computational model exhibiting self-organized criticality (SOC) based on adaptive networks allows further insights into the underlying dynamics. Together these results suggest that brain dynamics deviates from criticality during seizures caused by the failure of adaptive SOC.Comment: 7 pages, 5 figure

Brain Performance versus Phase Transitions

We here illustrate how a well-founded study of the brain may originate in assuming analogies with phase-transition phenomena. Analyzing to what extent a weak signal endures in noisy environments, we identify the underlying mechanisms, and it results a description of how the excitability associated to (non-equilibrium) phase changes and criticality optimizes the processing of the signal. Our setting is a network of integrate-and-fire nodes in which connections are heterogeneous with rapid time-varying intensities mimicking fatigue and potentiation. Emergence then becomes quite robust against wiring topology modification—in fact, we considered from a fully connected network to the Homo sapiens connectome—showing the essential role of synaptic flickering on computations. We also suggest how to experimentally disclose significant changes during actual brain operation.The authors acknowledge support from the Spanish Ministry of Economy and Competitiveness under the project FIS2013-43201-P