Sleep onset is characterized by a departure from arousal, and can be separated into well-differentiated stages: NREM (which encompasses three substages: N1, N2 and N3) and REM (Rapid Eye Movement). Awake brain dynamics are maintained by various wake-promoting mechanisms, particularly the neuromodulators Acetylcholine (ACh) and Noradrenaline (NA), whose levels naturally decrease during the transition to sleep. The combined influence of these neurotransmitters on brain connectivity during sleep remains unclear, as previous models have examined them mostly in isolation or only in deep sleep. In this study, we analyze fMRI data obtained from healthy individuals and employ a whole-brain model to investigate how changes in brain neurochemistry during NREM sleep, specifically involving ACh and NA, affect the Functional Connectivity (FC) of the brain. FC analysis reveals distinct connectivity changes: a decrease in Locus Coeruleus (LC) connectivity with the cortex during N2 and N3, and a decrease in Basal Forebrain (BF) connectivity with the cortex during N3. Additionally, compared to Wakefulness (W), there is a transition to a more integrated state in N1 and a more segregated state in N3. Using a Wilson-Cowan whole-brain model, informed by an empirical connectome and a heterogeneous receptivity map of neuromodulators, we explored possible mechanisms underlying these dynamics. We fit the model adjusting the coupling and input-output slope of the whole-brain model to account for ACh and NA, respectively, and show that region-specific neurotransmitter effect is key to explain their effects on FC. This work enhances our understanding of neurotransmitters' roles in modulating sleep stages and their significant contribution to brain state transitions between different states of consciousness, both in health and disease.

The human brain generates a rich repertoire of spatio-temporal activity patterns, which support a wide variety of motor and cognitive functions. These patterns of activity change with age in a multi-factorial manner. One of these factors is the variations in the brain's connectomics that occurs along the lifespan. However, the precise relationship between high-order functional interactions and connnectomics, as well as their variations with age are largely unknown, in part due to the absence of mechanistic models that can efficiently map brain connnectomics to functional connectivity in aging. To investigate this issue, we have built a neurobiologically-realistic whole-brain computational model using both anatomical and functional MRI data from 161 participants ranging from 10 to 80 years old. We show that the differences in high-order functional interactions between age groups can be largely explained by variations in the connectome. Based on this finding, we propose a simple neurodegeneration model that is representative of normal physiological aging. As such, when applied to connectomes of young participant it reproduces the age-variations that occur in the high-order structure of the functional data. Overall, these results begin to disentangle the mechanisms by which structural changes in the connectome lead to functional differences in the ageing brain. Our model can also serve as a starting point for modeling more complex forms of pathological ageing or cognitive deficits.

Segregation and integration are two fundamental principles of brain structural and functional organization. Neuroimaging studies have shown that the brain transits between different functionally segregated and integrated states, and neuromodulatory systems have been proposed as key to facilitate these transitions. Although whole-brain computational models have reproduced this neuromodulatory effect, the role of local inhibitory circuits and their cholinergic modulation has not been studied. In this article, we consider a Jansen & Rit whole-brain model in a network interconnected using a human connectome, and study the influence of the cholinergic and noradrenergic neuromodulatory systems on the segregation/integration balance. In our model, we introduce a local inhibitory feedback as a plausible biophysical mechanism that enables the integration of whole-brain activity, and that interacts with the other neuromodulatory influences to facilitate the transition between different functional segregation/integration regimes in the brain.

Neuromodulators, such as neuropeptides, can regulate and reconfigure neural circuits to alter their output, affecting in this way animal physiology and behavior. The interplay between the activity of neuronal circuits, their modulation by neuropeptides, and the resulting behavior, is still poorly understood. Here, we present a quantitative framework to study the relationships between the temporal pattern of activity of peptidergic neurons and of motoneurons during Drosophila ecdysis behavior, a highly stereotyped motor sequence that is critical for insect growth. We analyzed, in the time and frequency domains, simultaneous intracellular calcium recordings of peptidergic CCAP (crustacean cardioactive peptide) neurons and motoneurons obtained from isolated central nervous systems throughout fictive ecdysis behavior induced ex vivo by Ecdysis triggering hormone. We found that the activity of both neuronal populations is tightly coupled in a cross-frequency manner, suggesting that CCAP neurons modulate the frequency of motoneuron firing. To explore this idea further, we used a probabilistic logistic model to show that calcium dynamics in CCAP neurons can predict the oscillation of motoneurons, both in a simple model and in a conductance-based model capable of simulating many features of the observed neural dynamics. Finally, we developed an algorithm to quantify the motor behavior observed in videos of pupal ecdysis, and compared their features to the patterns of neuronal calcium activity recorded ex vivo. We found that the motor activity of the intact animal is more regular than the motoneuronal activity recorded from ex vivo preparations during fictive ecdysis behavior; the analysis of the patterns of movement also allowed us to identify a new post-ecdysis phase.

Slow-wave sleep cortical brain activity, conformed by slow-oscillations and sleep spindles, plays a key role in memory consolidation. The increase of the power of the slow-wave events, obtained by auditory sensory stimulation, positively correlates with memory consolidation performance. However, little is known about the experimental protocol maximizing this effect, which could be induced by the power of slow-oscillation, the number of sleep spindles, or the timing of both events' co-occurrence. Using a mean-field model of thalamocortical activity, we studied the effect of several stimulation protocols, varying the pulse shape, duration, amplitude, and frequency, as well as a target-phase using a closed-loop approach. We evaluated the effect of these parameters on slow-oscillations (SO) and sleep-spindles (SP), considering: (i) the power at the frequency bands of interest, (ii) the number of SO and SP, (iii) co-occurrences between SO and SP, and (iv) synchronization of SP with the up-peak of the SO. The first three targets are maximized using a decreasing ramp pulse with a pulse duration of 50 ms. Also, we observed a reduction in the number of SO when increasing the stimulus energy by rising its amplitude. To assess the target-phase parameter, we applied closed-loop stimulation at 0°, 45°, and 90° of the phase of the narrow-band filtered ongoing activity, at 0.85 Hz as central frequency. The 0° stimulation produces better results in the power and number of SO and SP than the rhythmic or random stimulation. On the other hand, stimulating at 45° or 90° change the timing distribution of spindles centers but with fewer co-occurrences than rhythmic and 0° phase. Finally, we propose the application of closed-loop stimulation at the rising zero-cross point using pulses with a decreasing ramp shape and 50 ms of duration for future experimental work.

The capability of cortical regions to flexibly sustain an "ignited" state of activity has been discussed in relation to conscious perception or hierarchical information processing. Here, we investigate how the intrinsic propensity of different regions to get ignited is determined by the specific topological organisation of the structural connectome. More specifically, we simulated the resting-state dynamics of mean-field whole-brain models and assessed how dynamic multistability and ignition differ between a reference model embedding a realistic human connectome, and alternative models based on a variety of randomised connectome ensembles. We found that the strength of global excitation needed to first trigger ignition in a subset of regions is substantially smaller for the model embedding the empirical human connectome. Furthermore, when increasing the strength of excitation, the propagation of ignition outside of this initial core-which is able to self-sustain its high activity-is way more gradual than for any of the randomised connectomes, allowing for graded control of the number of ignited regions. We explain both these assets in terms of the exceptional weighted core-shell organisation of the empirical connectome, speculating that this topology of human structural connectivity may be attuned to support enhanced ignition dynamics.