**Êtes-vous un étudiant de l'EPFL à la recherche d'un projet de semestre?**

Travaillez avec nous sur des projets en science des données et en visualisation, et déployez votre projet sous forme d'application sur GraphSearch.

Publication# Estimation of the depth of anaesthesia from interacting physiological oscillators

Résumé

Nowadays, physiological monitoring is imperative for the safety of medical operations. However, systems which monitor the depth of anaesthesia are still far from reliable, such that still some patients may experience the trauma of remaining conscious under general anaesthesia during surgery. The long term goal of our interdisciplinary project "BRACCIA" was to develop a device to measure the depth of anaesthesia. And, in view of this important goal, the main objective of research was to establish how the couplings between the cardiac, respiratory and cortical oscillations change in anaesthesia. Under the framework of this project, our objectives were: 1) The detection of the deep-light change of anaesthesia from experimental recordings on rats, and furthermore, the investigation of the interdependencies among three physiological activities, namely, the cardiac activity (H), respiration (R) and cortical activities (B) from experimental recordings of rats and humans, for each state of the depth of anaesthesia. 2) The modelling of the slow brain waves, and to consider the effect of anaesthesia on this realized model. The analysis of the recordings were carried out with five methods. First method is the "S-estimator", which indirectly quantifies the amount of synchronization within a data set measuring the contraction of the embedding dimension of the state space. Second method is the "new S-estimator". In this new one, a linear transformation of the reconstructed state space trajectory orthonormalizes the state variables within each model such that global state space volume reduction becomes a measure of synchronization exclusively between the different models. Third method is the "embedding dimension analysis", which examines the time evolution of the embedding dimension obtained with false nearest neighbors method on each windowed time series. The last two methods consist in a more detailed analysis of the dependencies among three systems. Fourth method is the "coupling matrix". This calculates the coupling matrix, CM, which infers linear interactions between multivariate time series after constructing separately the self model of each signal from reconstructed states. The last method is the "nonparametric Granger causality". This method calculates the Granger causality, GC, which measures bivariate causal influence in frequency domain. Here, a nonparametric estimation approach was used to remove the difficulties such as uncertainty in model parameters. For this method, permutation tests are added to figure out the real causality. From the obtained results, in most of the groups of experimental recordings, a change in the synchronization inside of the whole system between deep and light anaesthesia or between resting state and anaesthesia was found. For the results on Ketamine-Xylazine (KX) anaesthetized rats, a decrease of this synchronization was so clear that we succeeded to detect the deep-light transition of the anaesthesia, and automatically in one group of the recordings. The changes in the synchronization exclusively between three sub-systems as measured with the "new S-estimator" was only found on the results on Pentobarbital anaesthetized rats. Interestingly, a change in the time evolution of the embedding dimension was found on B and H in most of the results of the recordings. Furthermore, some changes in the couplings were found according to the depth of anaesthesia. For KX anaesthetized rats and human, coupling direction change was also observed. The nonparametric GC agreed for some causality change obtained with coupling matrix, but still some disagreements remained. From our literature study about the rhythmic activity of the brain, we deduce that it is generally accepted in the neuroscience community that the rhythms of the EEG are created by the interaction between the cortex and the thalamus. Consequently, we chose one of the Hodgkin-Huxley based model which considers the thalamocortical assemblies, the model of Bazhenov et al. This model includes four layers of neurons where two of them are for the thalamus (RE: thalamic reticular neuron, TC: thalamic relay neuron) and the other two are for the cortex (PY: pyramidal neuron, IN: interneuron). Not without having difficulties to set parameters and equations, we succeeded to reproduce this model and simulated it. Furthermore, we performed simulations on this model varying the parameter of maximal conductance of Ca2+ channels or the closing/opening rate of GABAA receptors from one simulation to another, and determined how the collective behavior of PY neurons, which should be interpreted as a field potential analogous to the EEG can be modified according to these parameter changes. Finally, to understand the model itself, some studies on the reproduced model were carried out. As a results, from observing the collective spiking activity of the PY cells on the reproduced Bazhenov model, we confirmed the presence of slow oscillations, namely the δ-waves. The reaction of the model to changes of the parameters related to anesthetics showed that reinforcement of the inhibition that suppresses the spiking activities was effectively induced by all changes of the parameters, but the degree of this reaction was different in each case. From the studies of the size of the model, we found that the local synchrony and the mean field frequency is not influenced by the size of network, but the global synchrony is not preserved at larger sizes. In conclusion, our data analysis studies showed a clear possibility for the detection of the depth of anaesthesia and a clear change in the interdependency changes depending on the anaesthesia. Moreover, our modelling study took an essential step for the investigation of the effect of anaesthesia on a brain model. As future work, all the unsolved problems which appear in the analysis chapter of this thesis should be treated, but especially the embedding dimension analysis should be followed up in detail. For modelling, as we just made the first step, we should continue with the understanding of the model and the closing of the gap between the mathematical model and the reality.

Source officielle

Cette page est générée automatiquement et peut contenir des informations qui ne sont pas correctes, complètes, à jour ou pertinentes par rapport à votre recherche. Il en va de même pour toutes les autres pages de ce site. Veillez à vérifier les informations auprès des sources officielles de l'EPFL.

Concepts associés (43)

MOOCs associés (32)

Publications associées (42)

Rythme cérébral

Un rythme cérébral (appelé aussi activité neuro-électrique) désigne l'oscillation électromagnétique émise par le cerveau des êtres humains, mais également de tout être vivant. Le cortex frontal qui permet la cognition, la logique et le raisonnement est composé de neurones qui sont reliés entre eux par des synapses permettant la neurotransmission. Mesurables en volt et en hertz, ces ondes sont de très faible amplitude : de l'ordre du microvolt (chez l'être humain), elles ne suivent pas toujours une sinusoïde régulière.

Gamma wave

A gamma wave or gamma rhythm is a pattern of neural oscillation in humans with a frequency between 25 and 140 Hz, the 40 Hz point being of particular interest. Gamma rhythms are correlated with large scale brain network activity and cognitive phenomena such as working memory, attention, and perceptual grouping, and can be increased in amplitude via meditation or neurostimulation. Altered gamma activity has been observed in many mood and cognitive disorders such as Alzheimer's disease, epilepsy, and schizophrenia.

Thalamocortical radiations

In neuroanatomy, thalamocortical radiations also known as thalamocortical fibres, are the efferent fibres that project from the thalamus to distinct areas of the cerebral cortex. They form fibre bundles that emerge from the lateral surface of the thalamus. Thalamocortical fibers (TC fibres) have been referred to as one of the two constituents of the isothalamus, the other being microneurons. Thalamocortical fibers have a bush or tree-like appearance as they extend into the internal capsule and project to the layers of the cortex.

Neuronal Dynamics - Computational Neuroscience of Single Neurons

The activity of neurons in the brain and the code used by these neurons is described by mathematical neuron models at different levels of detail.

Neuronal Dynamics - Computational Neuroscience of Single Neurons

The activity of neurons in the brain and the code used by these neurons is described by mathematical neuron models at different levels of detail.

Simulation Neurocience

Learn how to digitally reconstruct a single neuron to better study the biological mechanisms of brain function, behaviour and disease.

Daniel Keller, Andrii Stepaniuk

Neuroscience incorporates knowledge from a range of scales, from single molecules to brain wide neural networks. Modeling is a valuable tool in understanding processes at a single scale or the interactions between two adjacent scales and researchers use a ...

Wulfram Gerstner, Johanni Michael Brea, Alireza Modirshanechi, Shuqi Wang

Fitting network models to neural activity is an important tool in neuroscience. A popular approach is to model a brain area with a probabilistic recurrent spiking network whose parameters maximize the likelihood of the recorded activity. Although this is w ...

2021Felix Schürmann, Armando Romani, Michele Migliore, Luca Leonardo Bologna

To build biophysically detailed models of brain cells, circuits, and regions, a data-driven approach is increasingly being adopted. This helps to obtain a simulated activity that reproduces the experimentally recorded neural dynamics as faithfully as possi ...

2023