A neural mass model for disturbance of alpha rhythm in the minimal hepatic encephalopathy.

One of the early markers of minimal hepatic encephalopathy (MHE) is the disruption of alpha rhythm observed in electroencephalogram (EEG) signals. However, the underlying mechanisms responsible for this occurrence remain poorly understood. To address this gap, we develop a novel biophysical model MHE-AWD-NCM, encompassing the communication dynamics between a cortical neuron population (CNP) and an astrocyte population (AP), aimed at investigating the relationship between alpha wave disturbance (AWD) and mechanistical principles, specifically concerning astrocyte-neuronal communication in the context of MHE. In addition, we introduce the concepts of peak power density and peak frequency within the alpha band as quantitative measures of AWD. Our model faithfully reproduces the characteristic EEG phenomenology during MHE and shows how impairments of communication between CNP and AP could promote AWD. The results suggest that the disruptions in feedback neurotransmission from AP to CNP, along with the inhibition of GABA uptake by AP from the extracellular space, contribute to the observed AWD. Moreover, we found that the variation of external excitatory stimuli on CNP may play a key role in AWD in MHE. Finally, the sensitivity analysis is also performed to assess the relative significance of above factors in influencing AWD. Our findings align with the physiological observations and provide a more comprehensive understanding of the complex interplay of astrocyte-neuronal communication that underlies the AWD observed in MHE, which potentially may help to explore the targeted therapeutic interventions for the early stage of hepatic encephalopathy.

Cortical Synchrony and Information Flow during Transition from Wakefulness to Light Non-Rapid Eye Movement Sleep.

Sleep is a highly stereotyped phenomenon, requiring robust spatiotemporal coordination of neural activity. Understanding how the brain coordinates neural activity with sleep onset can provide insights into the physiological functions subserved by sleep and the pathologic phenomena associated with sleep onset. We quantified whole-brain network changes in synchrony and information flow during the transition from wakefulness to light non-rapid eye movement (NREM) sleep, using MEG imaging in a convenient sample of 14 healthy human participants (11 female; mean 63.4 years [SD 11.8 years]). We furthermore performed computational modeling to infer excitatory and inhibitory properties of local neural activity. The transition from wakefulness to light NREM was identified to be encoded in spatially and temporally specific patterns of long-range synchrony. Within the delta band, there was a global increase in connectivity from wakefulness to light NREM, which was highest in frontoparietal regions. Within the theta band, there was an increase in connectivity in fronto-parieto-occipital regions and a decrease in temporal regions from wakefulness to Stage 1 sleep. Patterns of information flow revealed that mesial frontal regions receive hierarchically organized inputs from broad cortical regions upon sleep onset, including direct inflow from occipital regions and indirect inflow via parieto-temporal regions within the delta frequency band. Finally, biophysical neural mass modeling demonstrated changes in the anterior-to-posterior distribution of cortical excitation-to-inhibition with increased excitation-to-inhibition model parameters in anterior regions in light NREM compared with wakefulness. Together, these findings uncover whole-brain corticocortical structure and the orchestration of local and long-range, frequency-specific cortical interactions in the sleep-wake transition.SIGNIFICANCE STATEMENT Our work uncovers spatiotemporal cortical structure of neural synchrony and information flow upon the transition from wakefulness to light non-rapid eye movement sleep. Mesial frontal regions were identified to receive hierarchically organized inputs from broad cortical regions, including both direct inputs from occipital regions and indirect inputs via the parieto-temporal regions within the delta frequency range. Biophysical neural mass modeling revealed a spatially heterogeneous, anterior-posterior distribution of cortical excitation-to-inhibition. Our findings shed light on the orchestration of local and long-range cortical neural structure that is fundamental to sleep onset, and support an emerging view of cortically driven regulation of sleep homeostasis.

Incorporating slow NMDA-type receptors with nonlinear voltage-dependent magnesium block in a next generation neural mass model: derivation and dynamics.

We derive a next generation neural mass model of a population of quadratic-integrate-and-fire neurons, with slow adaptation, and conductance-based AMPAR, GABAR and nonlinear NMDAR synapses. We show that the Lorentzian ansatz assumption can be satisfied by introducing a piece-wise polynomial approximation of the nonlinear voltage-dependent magnesium block of NMDAR current. We study the dynamics of the resulting system for two example cases of excitatory cortical neurons and inhibitory striatal neurons. Bifurcation diagrams are presented comparing the different dynamical regimes as compared to the case of linear NMDAR currents, along with sample comparison simulation time series demonstrating different possible oscillatory solutions. The omission of the nonlinearity of NMDAR currents results in a shift in the range (and possible disappearance) of the constant high firing rate regime, along with a modulation in the amplitude and frequency power spectrum of oscillations. Moreover, nonlinear NMDAR action is seen to be state-dependent and can have opposite effects depending on the type of neurons involved and the level of input firing rate received. The presented model can serve as a computationally efficient building block in whole brain network models for investigating the differential modulation of different types of synapses under neuromodulatory influence or receptor specific malfunction.

Transition to seizure in focal epilepsy: From SEEG phenomenology to underlying mechanisms.

OBJECTIVE: For the pre-surgical evaluation of patients with drug-resistant focal epilepsy, stereo-electroencephalographic (SEEG) signals are routinely recorded to identify the epileptogenic zone network (EZN). This network consists of remote brain regions involved in seizure initiation. However, the pathophysiological mechanisms underlying typical SEEG patterns that occur during the transition from interictal to ictal activity in distant brain nodes of the EZN remain poorly understood. The primary aim is to identify and explain these mechanisms using a novel physiologically-plausible model of the EZN. METHODS: We analyzed SEEG signals recorded from the EZN in 10 patients during the transition from interictal to ictal activity. This transition consisted of a sequence of periods during which SEEG signals from distant neocortical regions showed stereotypical patterns of activity: sustained preictal spiking activity preceding a fast activity occurring at seizure onset, followed by the ictal activity. Spectral content and non-linear correlation of SEEG signals were analyzed. In addition, we developed a novel neuro-inspired computational model consisting of bidirectionally coupled neuronal populations. RESULTS: The proposed model captured the essential characteristics of the patient signals, including the quasi-synchronous onset of rapid discharges in distant interconnected epileptogenic zones. Statistical analysis confirmed the dynamic correlation/de-decorrelation pattern observed in the patient signals and accurately reproduced in the simulated signals. SIGNIFICANCE: This study provides insight into the abnormal dynamic changes in glutamatergic and γ-aminobutyric acid (GABA)ergic synaptic transmission that occur during the transition to seizures. The results strongly support the hypothesis that bidirectional connections between distant neuronal populations of the EZN (from pyramidal cells to vaso-intestinal peptide-positive interneurons) play a key role in this transition, while parvalbumin-positive interneurons intervene in the emergence of rapid discharges at seizure onset.

Laminar neural dynamics of auditory evoked responses: Computational modeling of local field potentials in auditory cortex of non-human primates.

Evoked neural responses to sensory stimuli have been extensively investigated in humans and animal models both to enhance our understanding of brain function and to aid in clinical diagnosis of neurological and neuropsychiatric conditions. Recording and imaging techniques such as electroencephalography (EEG), magnetoencephalography (MEG), local field potentials (LFPs), and calcium imaging provide complementary information about different aspects of brain activity at different spatial and temporal scales. Modeling and simulations provide a way to integrate these different types of information to clarify underlying neural mechanisms. In this study, we aimed to shed light on the neural dynamics underlying auditory evoked responses by fitting a rate-based model to LFPs recorded via multi-contact electrodes which simultaneously sampled neural activity across cortical laminae. Recordings included neural population responses to best-frequency (BF) and non-BF tones at four representative sites in primary auditory cortex (A1) of awake monkeys. The model considered major neural populations of excitatory, parvalbumin-expressing (PV), and somatostatin-expressing (SOM) neurons across layers 2/3, 4, and 5/6. Unknown parameters, including the connection strength between the populations, were fitted to the data. Our results revealed similar population dynamics, fitted model parameters, predicted equivalent current dipoles (ECD), tuning curves, and lateral inhibition profiles across recording sites and animals, in spite of quite different extracellular current distributions. We found that PV firing rates were higher in BF than in non-BF responses, mainly due to different strengths of tonotopic thalamic input, whereas SOM firing rates were higher in non-BF than in BF responses due to lateral inhibition. In conclusion, we demonstrate the feasibility of the model-fitting approach in identifying the contributions of cell-type specific population activity to stimulus-evoked LFPs across cortical laminae, providing a foundation for further investigations into the dynamics of neural circuits underlying cortical sensory processing.

Accurate and Efficient Simulation of Very High-Dimensional Neural Mass Models with Distributed-Delay Connectome Tensors.

This paper introduces methods and a novel toolbox that efficiently integrates high-dimensional Neural Mass Models (NMMs) specified by two essential components. The first is the set of nonlinear Random Differential Equations (RDEs) of the dynamics of each neural mass. The second is the highly sparse three-dimensional Connectome Tensor (CT) that encodes the strength of the connections and the delays of information transfer along the axons of each connection. To date, simplistic assumptions prevail about delays in the CT, often assumed to be Dirac-delta functions. In reality, delays are distributed due to heterogeneous conduction velocities of the axons connecting neural masses. These distributed-delay CTs are challenging to model. Our approach implements these models by leveraging several innovations. Semi-analytical integration of RDEs is done with the Local Linearization (LL) scheme for each neural mass, ensuring dynamical fidelity to the original continuous-time nonlinear dynamic. This semi-analytic LL integration is highly computationally-efficient. In addition, a tensor representation of the CT facilitates parallel computation. It also seamlessly allows modeling distributed delays CT with any level of complexity or realism. This ease of implementation includes models with distributed-delay CTs. Consequently, our algorithm scales linearly with the number of neural masses and the number of equations they are represented with, contrasting with more traditional methods that scale quadratically at best. To illustrate the toolbox's usefulness, we simulate a single Zetterberg-Jansen and Rit (ZJR) cortical column, a single thalmo-cortical unit, and a toy example comprising 1000 interconnected ZJR columns. These simulations demonstrate the consequences of modifying the CT, especially by introducing distributed delays. The examples illustrate the complexity of explaining EEG oscillations, e.g., split alpha peaks, since they only appear for distinct neural masses. We provide an open-source Script for the toolbox.

Effect of channel density, inverse solutions and connectivity measures on EEG resting-state networks reconstruction: A simulation study.

Along with the study of brain activity evoked by external stimuli, the past two decades witnessed an increased interest in characterizing the spontaneous brain activity occurring during resting conditions. The identification of connectivity patterns in this so-called "resting-state" has been the subject of a great number of electrophysiology-based studies, using the Electro/Magneto-Encephalography (EEG/MEG) source connectivity method. However, no consensus has been reached yet regarding a unified (if possible) analysis pipeline, and several involved parameters and methods require cautious tuning. This is particularly challenging when different analytical choices induce significant discrepancies in results and drawn conclusions, thereby hindering the reproducibility of neuroimaging research. Hence, our objective in this study was to shed light on the effect of analytical variability on outcome consistency by evaluating the implications of parameters involved in the EEG source connectivity analysis on the accuracy of resting-state networks (RSNs) reconstruction. We simulated, using neural mass models, EEG data corresponding to two RSNs, namely the default mode network (DMN) and dorsal attentional network (DAN). We investigated the impact of five channel densities (19, 32, 64, 128, 256), three inverse solutions (weighted minimum norm estimate (wMNE), exact low-resolution brain electromagnetic tomography (eLORETA), and linearly constrained minimum variance (LCMV) beamforming) and four functional connectivity measures (phase-locking value (PLV), phase-lag index (PLI), and amplitude envelope correlation (AEC) with and without source leakage correction), on the correspondence between reconstructed and reference networks. We showed that, with different analytical choices related to the number of electrodes, source reconstruction algorithm, and functional connectivity measure, high variability is present in the results. More specifically, our results show that a higher number of EEG channels significantly increased the accuracy of the reconstructed networks. Additionally, our results showed significant variability in the performance of the tested inverse solutions and connectivity measures. Such methodological variability and absence of analysis standardization represent a critical issue for neuroimaging studies that should be prioritized. We believe that this work could be useful for the field of electrophysiology connectomics, by increasing awareness regarding the challenge of variability in methodological approaches and its implications on reported results.

Deep learning based source imaging provides strong sublobar localization of epileptogenic zone from MEG interictal spikes.

Electromagnetic source imaging (ESI) offers unique capability of imaging brain dynamics for studying brain functions and aiding the clinical management of brain disorders. Challenges exist in ESI due to the ill-posedness of the inverse problem and thus the need of modeling the underlying brain dynamics for regularizations. Advances in generative models provide opportunities for more accurate and realistic source modeling that could offer an alternative approach to ESI for modeling the underlying brain dynamics beyond equivalent physical source models. However, it is not straightforward to explicitly formulate the knowledge arising from these generative models within the conventional ESI framework. Here we investigate a novel source imaging framework based on mesoscale neuronal modeling and deep learning (DL) that can learn the sensor-source mapping relationship directly from MEG data for ESI. Two DL-based ESI models were trained based on data generated by neural mass models and either generic or personalized head models. The robustness of the two DL models was evaluated by systematic computer simulations and clinical validation in a cohort of 29 drug-resistant focal epilepsy patients who underwent intracranial EEG (iEEG) evaluation or surgical resection. Results estimated from pre-operative MEG interictal spikes were quantified using the overlap with resection regions and the distance to the seizure-onset zone (SOZ) defined by iEEG recordings. The DL-based ESI provided robust results when no personalized head geometry is considered, reaching a spatial dispersion of 21.90 ± 19.03 mm, sublobar concordance of 83 %, and sublobar sensitivity and specificity of 66 and 97 % respectively. When using personalized head geometry derived from individual patients' MRI in the training data, personalized DL-based ESI model can further improve the performance and reached a spatial dispersion of 8.19 ± 8.14 mm, sublobar concordance of 93 %, and sublobar sensitivity and specificity of 77 and 99 % respectively. When compared to the SOZ, the localization error of the personalized approach is 15.78 ± 5.54 mm, outperforming the conventional benchmarks. This work demonstrates that combining generative models and deep learning enables an accurate and robust imaging of epileptogenic zone from MEG recordings with strong sublobar precision, suggesting its added value to enhancing MEG source localization and imaging, and to epilepsy source localization and other clinical applications.

A brain atlas of axonal and synaptic delays based on modelling of cortico-cortical evoked potentials.

Epilepsy presurgical investigation may include focal intracortical single-pulse electrical stimulations with depth electrodes, which induce cortico-cortical evoked potentials at distant sites because of white matter connectivity. Cortico-cortical evoked potentials provide a unique window on functional brain networks because they contain sufficient information to infer dynamical properties of large-scale brain connectivity, such as preferred directionality and propagation latencies. Here, we developed a biologically informed modelling approach to estimate the neural physiological parameters of brain functional networks from the cortico-cortical evoked potentials recorded in a large multicentric database. Specifically, we considered each cortico-cortical evoked potential as the output of a transient stimulus entering the stimulated region, which directly propagated to the recording region. Both regions were modelled as coupled neural mass models, the parameters of which were estimated from the first cortico-cortical evoked potential component, occurring before 80 ms, using dynamic causal modelling and Bayesian model inversion. This methodology was applied to the data of 780 patients with epilepsy from the F-TRACT database, providing a total of 34 354 bipolar stimulations and 774 445 cortico-cortical evoked potentials. The cortical mapping of the local excitatory and inhibitory synaptic time constants and of the axonal conduction delays between cortical regions was obtained at the population level using anatomy-based averaging procedures, based on the Lausanne2008 and the HCP-MMP1 parcellation schemes, containing 130 and 360 parcels, respectively. To rule out brain maturation effects, a separate analysis was performed for older (>15 years) and younger patients (<15 years). In the group of older subjects, we found that the cortico-cortical axonal conduction delays between parcels were globally short (median = 10.2 ms) and only 16% were larger than 20 ms. This was associated to a median velocity of 3.9 m/s. Although a general lengthening of these delays with the distance between the stimulating and recording contacts was observed across the cortex, some regions were less affected by this rule, such as the insula for which almost all efferent and afferent connections were faster than 10 ms. Synaptic time constants were found to be shorter in the sensorimotor, medial occipital and latero-temporal regions, than in other cortical areas. Finally, we found that axonal conduction delays were significantly larger in the group of subjects younger than 15 years, which corroborates that brain maturation increases the speed of brain dynamics. To our knowledge, this study is the first to provide a local estimation of axonal conduction delays and synaptic time constants across the whole human cortex in vivo, based on intracerebral electrophysiological recordings.

Ion gradient-driven bifurcations of a multi-scale neuronal model.

Metabolic limitations within the brain frequently arise in the context of aging and disease. As the largest consumers of energy within the brain, ion pumps that maintain the neuronal membrane potential are the most affected when energy supply becomes limited. To characterize the effects of such limitations, we analyze the ion gradients present in a conductance-based (Morris-Lecar) neural mass model. We show the existence and locations of Neimark-Sacker and period-doubling bifurcations in the sodium, calcium, and potassium reversal potentials and demonstrate that these bifurcations form physiologically relevant bounds of ion gradient variability. Within these bounds, we show how depolarization of the gradients causes decreased neural activity. We also show that the depolarization of ion gradients decreases inter-regional coherence, causing a shift in the critical point at which the coupling occurs and thereby inducing loss of synchrony between regions. In this way, we show that the Larter-Breakspear model captures ion gradient variability present at the microscale level and propagates these changes to the macroscale effects such as those observed in human neuroimaging studies.

Suppression of beta oscillations by delayed feedback in a cortex-basal ganglia-thalamus-pedunculopontine nucleus neural loop model.

Excessive neural synchronization of neural populations in the beta (ß) frequency range (12-35 Hz) is intimately related to the symptoms of hypokinesia in Parkinson's disease (PD). Studies have shown that delayed feedback stimulation strategies can interrupt excessive neural synchronization and effectively alleviate symptoms associated with PD dyskinesia. Work on optimizing delayed feedback algorithms continues to progress, yet it remains challenging to further improve the inhibitory effect with reduced energy expenditure. Therefore, we first established a neural mass model of the cortex-basal ganglia-thalamus-pedunculopontine nucleus (CBGTh-PPN) closed-loop system, which can reflect the internal properties of cortical and basal ganglia neurons and their intrinsic connections with thalamic and pedunculopontine nucleus neurons. Second, the inhibitory effects of three delayed feedback schemes based on the external globus pallidum (GPe) on ß oscillations were investigated separately and compared with those based on the subthalamic nucleus (STN) only. Our results show that all four delayed feedback schemes achieve effective suppression of pathological ß oscillations when using the linear delayed feedback algorithm. The comparison revealed that the three GPe-based delayed feedback stimulation strategies were able to have a greater range of oscillation suppression with reduced energy consumption, thus improving control performance effectively, suggesting that they may be more effective for the relief of Parkinson's motor symptoms in practical applications.

Design of Closed-Loop Control Schemes Based on the GA-PID and GA-RBF-PID Algorithms for Brain Dynamic Modulation.

Neurostimulation can be used to modulate brain dynamics of patients with neuropsychiatric disorders to make abnormal neural oscillations restore to normal. The control schemes proposed on the bases of neural computational models can predict the mechanism of neural oscillations induced by neurostimulation, and then make clinical decisions that are suitable for the patient's condition to ensure better treatment outcomes. The present work proposes two closed-loop control schemes based on the improved incremental proportional integral derivative (PID) algorithms to modulate brain dynamics simulated by Wendling-type coupled neural mass models. The introduction of the genetic algorithm (GA) in traditional incremental PID algorithm aims to overcome the disadvantage that the selection of control parameters depends on the designer's experience, so as to ensure control accuracy. The introduction of the radial basis function (RBF) neural network aims to improve the dynamic performance and stability of the control scheme by adaptively adjusting control parameters. The simulation results show the high accuracy of the closed-loop control schemes based on GA-PID and GA-RBF-PID algorithms for modulation of brain dynamics, and also confirm the superiority of the scheme based on the GA-RBF-PID algorithm in terms of the dynamic performance and stability. This research of making hypotheses and predictions according to model data is expected to improve and perfect the equipment of early intervention and rehabilitation treatment for neuropsychiatric disorders in the biomedical engineering field.

Space-time resolved inference-based neurophysiological process imaging: Application to resting-state alpha rhythm.

Neural processes are complex and difficult to image. This paper presents a new space-time resolved brain imaging framework, called Neurophysiological Process Imaging (NPI), that identifies neurophysiological processes within cerebral cortex at the macroscopic scale. By fitting uncoupled neural mass models to each electromagnetic source time-series using a novel nonlinear inference method, population averaged membrane potentials and synaptic connection strengths are efficiently and accurately inferred and imaged across the whole cerebral cortex at a resolution afforded by source imaging. The efficiency of the framework enables return of the augmented source imaging results overnight using high performance computing. This suggests it can be used as a practical and novel imaging tool. To demonstrate the framework, it has been applied to resting-state magnetoencephalographic source estimates. The results suggest that endogenous inputs to cingulate, occipital, and inferior frontal cortex are essential modulators of resting-state alpha power. Moreover, endogenous input and inhibitory and excitatory neural populations play varied roles in mediating alpha power in different resting-state sub-networks. The framework can be applied to arbitrary neural mass models and has broad applicability to image neural processes of different brain states.

Mean-Field Models for EEG/MEG: From Oscillations to Waves.

Neural mass models have been used since the 1970s to model the coarse-grained activity of large populations of neurons. They have proven especially fruitful for understanding brain rhythms. However, although motivated by neurobiological considerations they are phenomenological in nature, and cannot hope to recreate some of the rich repertoire of responses seen in real neuronal tissue. Here we consider a simple spiking neuron network model that has recently been shown to admit an exact mean-field description for both synaptic and gap-junction interactions. The mean-field model takes a similar form to a standard neural mass model, with an additional dynamical equation to describe the evolution of within-population synchrony. As well as reviewing the origins of this next generation mass model we discuss its extension to describe an idealised spatially extended planar cortex. To emphasise the usefulness of this model for EEG/MEG modelling we show how it can be used to uncover the role of local gap-junction coupling in shaping large scale synaptic waves.

Mean-Field Modeling of Brain-Scale Dynamics for the Evaluation of EEG Source-Space Networks.

Understanding the dynamics of brain-scale functional networks at rest and during cognitive tasks is the subject of intense research efforts to unveil fundamental principles of brain functions. To estimate these large-scale brain networks, the emergent method called "electroencephalography (EEG) source connectivity" has generated increasing interest in the network neuroscience community, due to its ability to identify cortical brain networks with satisfactory spatio-temporal resolution, while reducing mixing and volume conduction effects. However, no consensus has been reached yet regarding a unified EEG source connectivity pipeline, and several methodological issues have to be carefully accounted to avoid pitfalls. Thus, a validation toolbox that provides flexible "ground truth" models is needed for an objective methods/parameters evaluation and, thereby an optimization of the EEG source connectivity pipeline. In this paper, we show how a recently developed large-scale model of brain-scale activity, named COALIA, can provide to some extent such ground truth by providing realistic simulations of source-level and scalp-level activity. Using a bottom-up approach, the model bridges cortical micro-circuitry and large-scale network dynamics. Here, we provide an example of the potential use of COALIA to analyze, in the context of epileptiform activity, the effect of three key factors involved in the "EEG source connectivity" pipeline: (i) EEG sensors density, (ii) algorithm used to solve the inverse problem, and (iii) functional connectivity measure. Results showed that a high electrode density (at least 64 channels) is required to accurately estimate cortical networks. Regarding the inverse solution/connectivity measure combination, the best performance at high electrode density was obtained using the weighted minimum norm estimate (wMNE) combined with the weighted phase lag index (wPLI). Although those results are specific to the considered aforementioned context (epileptiform activity), we believe that this model-based approach can be successfully applied to other experimental questions/contexts. We aim at presenting a proof-of-concept of the interest of COALIA in the network neuroscience field, and its potential use in optimizing the EEG source-space network estimation pipeline.

Whole-Brain Modelling: Past, Present, and Future.

Whole-Brain Modelling is a scientific field with a short history and a long past. Its various disciplinary roots and conceptual ingredients extend back to as early as the 1940s. It was not until the late 2000s, however, that a nascent paradigm emerged in roughly its current form-concurrently, and in many ways joined at the hip, with its sister field of macro-connectomics. This period saw a handful of seminal papers authored by a certain motley crew of notable theoretical and cognitive neuroscientists, which have served to define much of the landscape of whole-brain modelling as it stands at the start of the 2020s. At the same time, the field has over the past decade expanded in a dozen or more fascinating new methodological, theoretical, and clinical directions. In this chapter we offer a potted Past, Present, and Future of whole-brain modelling, noting what we take to be some of its greatest successes, hardest challenges, and most exciting opportunities.

The essential role of hippocampo-cortical connections in temporal coordination of spindles and ripples.

The predominant activity of slow wave sleep is cortical slow oscillations (SOs), thalamic spindles and hippocampal sharp wave ripples. While the precise temporal nesting of these rhythms was shown to be essential for memory consolidation, the coordination mechanism is poorly understood. Here we develop a minimal hippocampo-cortico-thalamic network that can explain the mechanism underlying the SO-spindle-ripple coupling indicating of the succession of regional neuronal interactions. Further we verify the model predictions experimentally in naturally sleeping rodents showing our simple model provides a quantitative match to several experimental observations including the nesting of ripples in the spindle troughs and larger duration but lower amplitude of the ripples co-occurring with spindles or SOs compared to the isolated ripples. The model also predicts that the coupling of ripples to SOs and spindles monotonically enhances by increasing the strength of hippocampo-cortical connections while it is stronger at intermediate values of the cortico-hippocampal projections.

NMDA-receptor antibodies alter cortical microcircuit dynamics.

NMDA-receptor antibodies (NMDAR-Abs) cause an autoimmune encephalitis with a diverse range of EEG abnormalities. NMDAR-Abs are believed to disrupt receptor function, but how blocking this excitatory synaptic receptor can lead to paroxysmal EEG abnormalities-or even seizures-is poorly understood. Here we show that NMDAR-Abs change intrinsic cortical connections and neuronal population dynamics to alter the spectral composition of spontaneous EEG activity and predispose brain dynamics to paroxysmal abnormalities. Based on local field potential recordings in a mouse model, we first validate a dynamic causal model of NMDAR-Ab effects on cortical microcircuitry. Using this model, we then identify the key synaptic parameters that best explain EEG paroxysms in pediatric patients with NMDAR-Ab encephalitis. Finally, we use the mouse model to show that NMDAR-Ab-related changes render microcircuitry critically susceptible to overt EEG paroxysms when these key parameters are changed, even though the same parameter fluctuations are tolerated in the in silico model of the control condition. These findings offer mechanistic insights into circuit-level dysfunction induced by NMDAR-Ab.

Generative modelling of the thalamo-cortical circuit mechanisms underlying the neurophysiological effects of ketamine.

Cortical recordings of task-induced oscillations following subanaesthetic ketamine administration demonstrate alterations in amplitude, including increases at high-frequencies (gamma) and reductions at low frequencies (theta, alpha). To investigate the population-level interactions underlying these changes, we implemented a thalamo-cortical model (TCM) capable of recapitulating broadband spectral responses. Compared with an existing cortex-only 4-population model, Bayesian Model Selection preferred the TCM. The model was able to accurately and significantly recapitulate ketamine-induced reductions in alpha amplitude and increases in gamma amplitude. Parameter analysis revealed no change in receptor time-constants but significant increases in select synaptic connectivity with ketamine. Significantly increased connections included both AMPA and NMDA mediated connections from layer 2/3 superficial pyramidal cells to inhibitory interneurons and both GABAA and NMDA mediated within-population gain control of layer 5 pyramidal cells. These results support the use of extended generative models for explaining oscillatory data and provide in silico support for ketamine's ability to alter local coupling mediated by NMDA, AMPA and GABA-A.

Comparing dynamic causal models of neurovascular coupling with fMRI and EEG/MEG.

This technical note presents a dynamic causal modelling (DCM) procedure for evaluating different models of neurovascular coupling in the human brain - using combined electromagnetic (M/EEG) and functional magnetic resonance imaging (fMRI) data. This procedure compares the evidence for biologically informed models of neurovascular coupling using Bayesian model comparison. First, fMRI data are used to localise regionally specific neuronal responses. The coordinates of these responses are then used as the location priors in a DCM of electrophysiological responses elicited by the same paradigm. The ensuing estimates of model parameters are then used to generate neuronal drive functions, which model pre- or post-synaptic activity for each experimental condition. These functions form the input to a model of neurovascular coupling, whose parameters are estimated from the fMRI data. Crucially, this enables one to evaluate different models of neurovascular coupling, using Bayesian model comparison - asking, for example, whether instantaneous or delayed, pre- or post-synaptic signals mediate haemodynamic responses. We provide an illustrative application of the procedure using a single-subject auditory fMRI and MEG dataset. The code and exemplar data accompanying this technical note are available through the statistical parametric mapping (SPM) software.