Across vertebrate species, sleep consists of repeating cycles of NREM followed by REM. However, their respective functions, and their stereotypic cycling pattern are not well understood. Using a simplified biophysical network model, we investigate the potential role of cholinergic modulation, acting via the muscarinic receptors, on network dynamics and memory consolidation. We show that low and high cholinergic levels associated with NREM and REM sleep, respectively, may play critical, sequential roles in memory consolidation. The network dynamics that facilitate these roles arise through alteration of neural excitability and changes to network-wide excitatory/inhibitory balance. At low acetylcholine (ACh) levels, reduced activation of inhibitory neurons leads to network-wide disinhibition and bursts of synchronized activity led by engram neurons, driving recruitment of additional excitatory neurons into the engram. In contrast, at high ACh levels, increased network inhibition suppresses firing in all but the most strongly recruited excitatory neurons, pruning the expanded engram population. Together, these results provide a testable hypothesis regarding the role of sleep state-specific cholinergic modulation in the process of memory consolidation.

General anesthetics work through a variety of molecular mechanisms while resulting in the common end point of sedation and loss of consciousness. Generally, the administration of common anesthetics induces reduction in synaptic excitation while promoting synaptic inhibition. Exogenous modulation of the anesthetics' synaptic effects can help determine the neuronal pathways involved in anesthesia. For example, both animal and human studies have shown that exogenously induced increases in acetylcholine in the brain can elicit wakeful-like behavior despite the continued presence of the anesthetic. However, the underlying mechanisms of anesthesia reversal at the cellular level have not been investigated. Here we apply a computational model of a network of excitatory and inhibitory neurons to simulate the network-wide effects of anesthesia, due to changes in synaptic inhibition and excitation, and their reversal by cholinergic activation through muscarinic receptors. We use a differential evolution algorithm to fit model parameters to match measures of spiking activity, neuronal connectivity, and network dynamics recorded in the visual cortex of rodents during anesthesia with desflurane in vivo. We find that facilitating muscarinic receptor effects of acetylcholine on top of anesthetic-induced synaptic changes predicts the reversal of anesthetic suppression of neurons' spiking activity, functional connectivity, as well as pairwise and population interactions. Thus, our model predicts a specific neuronal mechanism for the cholinergic reversal of anesthesia consistent with experimental behavioral observations.

Hidden hearing loss (HHL) is an auditory neuropathy characterized by normal hearing thresholds but reduced amplitudes of the sound-evoked auditory nerve compound action potential (CAP). In animal models, HHL can be caused by moderate noise exposure or aging, which induces loss of inner hair cell (IHC) synapses. In contrast, recent evidence has shown that transient loss of cochlear Schwann cells also causes permanent auditory deficits in mice with similarities to HHL. Histological analysis of the cochlea after auditory nerve remyelination showed a permanent disruption of the myelination patterns at the heminode of type I spiral ganglion neuron (SGN) peripheral terminals, suggesting that this defect could be contributing to HHL. To shed light on the mechanisms of different HHL scenarios observed in animals and to test their impact on type I SGN activity, we constructed a reduced biophysical model for a population of SGN peripheral axons whose activity is driven by a well-accepted model of cochlear sound processing. We found that the amplitudes of simulated sound-evoked SGN CAPs are lower and have greater latencies when heminodes are disorganized, i.e. they occur at different distances from the hair cell rather than at the same distance as in the normal cochlea. These results confirm that disruption of heminode positions causes desynchronization of SGN spikes leading to a loss of temporal resolution and reduction of the sound-evoked SGN CAP. Another mechanism resulting in HHL is loss of IHC synapses, i.e., synaptopathy. For comparison, we simulated synaptopathy by removing high threshold IHC-SGN synapses and found that the amplitude of simulated sound-evoked SGN CAPs decreases while latencies remain unchanged, as has been observed in noise exposed animals. Thus, model results illuminate diverse disruptions caused by synaptopathy and demyelination on neural activity in auditory processing that contribute to HHL as observed in animal models and that can contribute to perceptual deficits induced by nerve damage in humans.

[This corrects the article DOI: 10.1371/journal.pcbi.1008499.].

Sleep is critical for memory consolidation, although the exact mechanisms mediating this process are unknown. Combining reduced network models and analysis of in vivo recordings, we tested the hypothesis that neuromodulatory changes in acetylcholine (ACh) levels during non-rapid eye movement (NREM) sleep mediate stabilization of network-wide firing patterns, with temporal order of neurons' firing dependent on their mean firing rate during wake. In both reduced models and in vivo recordings from mouse hippocampus, we find that the relative order of firing among neurons during NREM sleep reflects their relative firing rates during prior wake. Our modeling results show that this remapping of wake-associated, firing frequency-based representations is based on NREM-associated changes in neuronal excitability mediated by ACh-gated potassium current. We also show that learning-dependent reordering of sequential firing during NREM sleep, together with spike timing-dependent plasticity (STDP), reconfigures neuronal firing rates across the network. This rescaling of firing rates has been reported in multiple brain circuits across periods of sleep. Our model and experimental data both suggest that this effect is amplified in neural circuits following learning. Together our data suggest that sleep may bias neural networks from firing rate-based towards phase-based information encoding to consolidate memories.

Theta and gamma rhythms and their cross-frequency coupling play critical roles in perception, attention, learning, and memory. Available data suggest that forebrain acetylcholine (ACh) signaling promotes theta-gamma coupling, although the mechanism has not been identified. Recent evidence suggests that cholinergic signaling is both temporally and spatially constrained, in contrast to the traditional notion of slow, spatially homogeneous, and diffuse neuromodulation. Here, we find that spatially constrained cholinergic stimulation can generate theta-modulated gamma rhythms. Using biophysically-based excitatory-inhibitory (E-I) neural network models, we simulate the effects of ACh on neural excitability by varying the conductance of a muscarinic receptor-regulated K+ current. In E-I networks with local excitatory connectivity and global inhibitory connectivity, we demonstrate that theta-gamma-coupled firing patterns emerge in ACh modulated network regions. Stable gamma-modulated firing arises within regions with high ACh signaling, while theta or mixed theta-gamma activity occurs at the peripheries of these regions. High gamma activity also alternates between different high-ACh regions, at theta frequency. Our results are the first to indicate a causal role for spatially heterogenous ACh signaling in the emergence of localized theta-gamma rhythmicity. Our findings also provide novel insights into mechanisms by which ACh signaling supports the brain region-specific attentional processing of sensory information.

Acetylcholine (ACh) is a regulator of neural excitability and one of the neurochemical substrates of sleep. Amongst the cellular effects induced by cholinergic modulation are a reduction in spike-frequency adaptation (SFA) and a shift in the phase response curve (PRC). We demonstrate in a biophysical model how changes in neural excitability and network structure interact to create three distinct functional regimes: localized asynchronous, traveling asynchronous, and traveling synchronous. Our results qualitatively match those observed experimentally. Cortical activity during slow wave sleep (SWS) differs from that during REM sleep or waking states. During SWS there are traveling patterns of activity in the cortex; in other states stationary patterns occur. Our model is a network composed of Hodgkin-Huxley type neurons with a M-current regulated by ACh. Regulation of ACh level can account for dynamical changes between functional regimes. Reduction of the magnitude of this current recreates the reduction in SFA the shift from a type 2 to a type 1 PRC observed in the presence of ACh. When SFA is minimal (in waking or REM sleep state, high ACh) patterns of activity are localized and easily pinned by network inhomogeneities. When SFA is present (decreasing ACh), traveling waves of activity naturally arise. A further decrease in ACh leads to a high degree of synchrony within traveling waves. We also show that the level of ACh determines how sensitive network activity is to synaptic heterogeneity. These regimes may have a profound functional significance as stationary patterns may play a role in the proper encoding of external input as memory and traveling waves could lead to synaptic regularization, giving unique insights into the role and significance of ACh in determining patterns of cortical activity and functional differences arising from the patterns.

Although sleep is a fundamental behavior observed in virtually all animal species, its functions remain unclear. One leading proposal, known as the synaptic renormalization hypothesis, suggests that sleep is necessary to counteract a global strengthening of synapses that occurs during wakefulness. Evidence for sleep-dependent synaptic downscaling (or synaptic renormalization) has been observed experimentally, but the physiological mechanisms which generate this phenomenon are unknown. In this study, we propose that changes in neuronal membrane excitability induced by acetylcholine may provide a dynamical mechanism for both wake-dependent synaptic upscaling and sleep-dependent downscaling. We show in silico that cholinergically-induced changes in network firing patterns alter overall network synaptic potentiation when synaptic strengths evolve through spike-timing dependent plasticity mechanisms. Specifically, network synaptic potentiation increases dramatically with high cholinergic concentration and decreases dramatically with low levels of acetylcholine. We demonstrate that this phenomenon is robust across variation of many different network parameters.

Within the brain, the interplay between connectivity patterns of neurons and their spatiotemporal dynamics is believed to be intricately linked to the bases of behavior, such as the process of storing, consolidating, and retrieving memory traces. Memory is believed to be stored in the synaptic patterns of anatomical circuitry in the form of increased connectivity densities within subpopulations of neurons. At the same time, memory recall is thought to correspond to activation of discrete areas of the brain corresponding to those memories. Such regional subpopulations can selectively activate during memory recall or retrieval, signifying the process of accessing a single memory or concept. It has been shown previously that recovery of single memory activity patterns is mediated by global neuromodulation signifying transition into different cognitive states such as sleep or awake exploration. We examine how underlying topology can affect memory awake activation and sleep reactivation when such memories share increasing proportions of neurons. The results show that while single memory activation is diminished with increased overlap, pattern separation can be recovered by offsetting excitatory associations between two memories with targeted and heterogeneous inhibitory feedback. Such findings point to the importance of excitatory-to-inhibitory current balance at both the global and local levels in the context of memory retrieval and replay, and highlight the role of network topology in memory management processes.

Spatiotemporal pattern formation in neuronal networks depends on the interplay between cellular and network synchronization properties. The neuronal phase response curve (PRC) is an experimentally obtainable measure that characterizes the cellular response to small perturbations, and can serve as an indicator of cellular propensity for synchronization. Two broad classes of PRCs have been identified for neurons: Type I, in which small excitatory perturbations induce only advances in firing, and Type II, in which small excitatory perturbations can induce both advances and delays in firing. Interestingly, neuronal PRCs are usually attenuated with increased spiking frequency, and Type II PRCs typically exhibit a greater attenuation of the phase delay region than of the phase advance region. We found that this phenomenon arises from an interplay between the time constants of active ionic currents and the interspike interval. As a result, excitatory networks consisting of neurons with Type I PRCs responded very differently to frequency modulation compared to excitatory networks composed of neurons with Type II PRCs. Specifically, increased frequency induced a sharp decrease in synchrony of networks of Type II neurons, while frequency increases only minimally affected synchrony in networks of Type I neurons. These results are demonstrated in networks in which both types of neurons were modeled generically with the Morris-Lecar model, as well as in networks consisting of Hodgkin-Huxley-based model cortical pyramidal cells in which simulated effects of acetylcholine changed PRC type. These results are robust to different network structures, synaptic strengths and modes of driving neuronal activity, and they indicate that Type I and Type II excitatory networks may display two distinct modes of processing information.

It has been hypothesized that the brain uses combinatorial as well as temporal coding strategies to represent stimulus properties. The mechanisms and properties of the temporal coding remain undetermined, although it has been postulated that oscillations can mediate formation of this type of code. Here we use a generic model of the vertebrate olfactory bulb to explore the possible role of oscillatory behavior in temporal coding. We show that three mechanisms—synaptic inhibition, slow self-inhibition and input properties—mediate formation of a temporal sequence of simultaneous activations of glomerular modules associated with specific odorants within the oscillatory response. The sequence formed depends on the relative properties of odorant features and thus may mediate discrimination of odorants activating overlapping sets of glomeruli. We suggest that period-doubling transitions may be driven through excitatory feedback from a portion of the olfactory network acting as a coincidence modulator. Furthermore, we hypothesize that the period-doubling transition transforms the temporal code from a roster of odorant components to a signal of odorant identity and facilitates discrimination of individual odorants within mixtures.

We modify neural networks models of the Hopfield type so that they can recognize the degree of novelty of the input stimuli on a local level. The networks control themselves the quality of recognition and can also recognize locally the bits of information in the input patterns which do not agree with known patterns, i.e. stored memories. This task is achieved by introducing local variations of the noise level beta in the network. Noise level in a given location depends on the flip frequency of the neurons close to that location.