RESUMEN
The electroencephalogram (EEG) patterns recorded during general anesthetic-induced coma are closely similar to those seen during slow-wave sleep, the deepest stage of natural sleep; both states show patterns dominated by large amplitude slow waves. Slow oscillations are believed to be important for memory consolidation during natural sleep. Tracking the emergence of slow-wave oscillations during transition to unconsciousness may help us to identify drug-induced alterations of the underlying brain state, and provide insight into the mechanisms of general anesthesia. Although cellular-based mechanisms have been proposed, the origin of the slow oscillation has not yet been unambiguously established. A recent theoretical study by Steyn-Ross et al. (2013) proposes that the slow oscillation is a network, rather than cellular phenomenon. Modeling anesthesia as a moderate reduction in gap-junction interneuronal coupling, they predict an unconscious state signposted by emergent low-frequency oscillations with chaotic dynamics in space and time. They suggest that anesthetic slow-waves arise from a competitive interaction between symmetry-breaking instabilities in space (Turing) and time (Hopf), modulated by gap-junction coupling strength. A significant prediction of their model is that EEG phase coherence will decrease as the cortex transits from Turing-Hopf balance (wake) to Hopf-dominated chaotic slow-waves (unconsciousness). Here, we investigate changes in phase coherence during induction of general anesthesia. After examining 128-channel EEG traces recorded from five volunteers undergoing propofol anesthesia, we report a significant drop in sub-delta band (0.05-1.5 Hz) slow-wave coherence between frontal, occipital, and frontal-occipital electrode pairs, with the most pronounced wake-vs.-unconscious coherence changes occurring at the frontal cortex.
RESUMEN
BACKGROUND: Investigation of the nonlinear pattern dynamics of a reaction-diffusion system almost always requires numerical solution of the system's set of defining differential equations. Traditionally, this would be done by selecting an appropriate differential equation solver from a library of such solvers, then writing computer codes (in a programming language such as C or Matlab) to access the selected solver and display the integrated results as a function of space and time. This "code-based" approach is flexible and powerful, but requires a certain level of programming sophistication. A modern alternative is to use a graphical programming interface such as Simulink to construct a data-flow diagram by assembling and linking appropriate code blocks drawn from a library. The result is a visual representation of the inter-relationships between the state variables whose output can be made completely equivalent to the code-based solution. RESULTS: As a tutorial introduction, we first demonstrate application of the Simulink data-flow technique to the classical van der Pol nonlinear oscillator, and compare Matlab and Simulink coding approaches to solving the van der Pol ordinary differential equations. We then show how to introduce space (in one and two dimensions) by solving numerically the partial differential equations for two different reaction-diffusion systems: the well-known Brusselator chemical reactor, and a continuum model for a two-dimensional sheet of human cortex whose neurons are linked by both chemical and electrical (diffusive) synapses. We compare the relative performances of the Matlab and Simulink implementations. CONCLUSIONS: The pattern simulations by Simulink are in good agreement with theoretical predictions. Compared with traditional coding approaches, the Simulink block-diagram paradigm reduces the time and programming burden required to implement a solution for reaction-diffusion systems of equations. Construction of the block-diagram does not require high-level programming skills, and the graphical interface lends itself to easy modification and use by non-experts.