Perturbation of whole-brain dynamics in silico reveals mechanistic differences between brain states.

Human neuroimaging research has revealed that wakefulness and sleep involve very different activity patterns. Yet, it is not clear why brain states differ in their dynamical complexity, e.g. in the level of integration and segregation across brain networks over time. Here, we investigate the mechanisms underlying the dynamical stability of brain states using a novel off-line in silico perturbation protocol. We first adjust a whole-brain computational model to the basal dynamics of wakefulness and deep sleep recorded with fMRI in two independent human fMRI datasets. Then, the models of sleep and awake brain states are perturbed using two distinct multifocal protocols either promoting or disrupting synchronization in randomly selected brain areas. Once perturbation is halted, we use a novel measure, the Perturbative Integration Latency Index (PILI), to evaluate the recovery back to baseline. We find a clear distinction between models, consistently showing larger PILI in wakefulness than in deep sleep, corroborating previous experimental findings. In the models, larger recoveries are associated to a critical slowing down induced by a shift in the model's operation point, indicating that the awake brain operates further from a stable equilibrium than deep sleep. This novel approach opens up for a new level of artificial perturbative studies unconstrained by ethical limitations allowing for a deeper investigation of the dynamical properties of different brain states.