Population based models of cortical drug response: insights from anaesthesia

Efficient and Accurate Computational Model of Neuron with Spike Frequency Adaptation

Simplified models of neurons are widely used in computational investigations of large networks. One of the most important performance metrics of simplified models is their accuracy in reproducing action potential (spike) timing. In this article, we developed a simple, computationally efficient neuron model by modifying the adaptive exponential integrate and fire (AdEx) model [1] with sigmoid afterhyperpolarization current (Sigmoid AHP). Our model can precisely match the spike times and spike frequency adaptation of cortical pyramidal neurons. The accuracy was similar to a more complex two compartment biophysically realistic model of the same neurons. This work provides a simplified neuronal model with improved spike timing accuracy for use in modeling of large neural networks.Clinical Relevance- Accurate and computationally efficient single neuron model will enable large network modeling of brain regions involved in neurological and psychiatric disorders and may lead to a better understanding of the disorder mechanisms.

Transient neocortical gamma oscillations induced by neuronal response modulation

In this paper a mean field model of spatio-temporal electroencephalographic activity in the neocortex is used to computationally study the emergence of neocortical gamma oscillations as a result of neuronal response modulation. It is shown using a numerical bifurcation analysis that gamma oscillations emerge robustly in the solutions of the model and transition to beta oscillations through coordinated modulation of the responsiveness of inhibitory and excitatory neuronal populations. The spatio-temporal pattern of the propagation of these oscillations across the neocortex is illustrated by solving the equations of the model using a finite element software package. Thereby, it is shown that the gamma oscillations remain localized to the regions of neuronal modulation. Moreover, it is discussed that the inherent spatial averaging effect of commonly used electrocortical measurement techniques can significantly alter the amplitude and pattern of fast oscillations in neocortical recordings, and hence can potentially affect physiological interpretations of these recordings.

A Computational Study of a Spatiotemporal Mean Field Model Capturing the Emergence of Alpha and Gamma Rhythmic Activity in the Neocortex

In this paper, we analyze the spatiotemporal mean field model developed by Liley et al. in order to advance our understanding of the wide effects of pharmacological agents and anesthetics. Specifically, we use the spatiotemporal mean field model for capturing the electrical activity in the neocortex to computationally study the emergence of α - and γ -band rhythmic activity in the brain. We show that α oscillations in the solutions of the model appear globally across the neocortex, whereas γ oscillations can emerge locally as a result of a bifurcation in the dynamics of the model. We solve the dynamic equations of the model using a finite element solver package and show that our results verify the predictions made by bifurcation analysis.

Anesthesia modifies subthreshold critical slowing down in a stochastic Hodgkin-Huxley-like model with inhibitory synaptic input

The dynamics of a stochastic type-I Hodgkin-Huxley-like point neuron model exposed to inhibitory synaptic noise are investigated as a function of distance from spiking threshold and the inhibitory influence of the general anesthetic agent propofol. The model is biologically motivated and includes the effects of intrinsic ion-channel noise via a stochastic differential equation description as well as inhibitory synaptic noise modeled as multiple Poisson-distributed impulse trains with saturating response functions. The effect of propofol on these synapses is incorporated through this drug's principal influence on fast inhibitory neurotransmission mediated by γ-aminobutyric acid (GABA) type-A receptors via reduction of the synaptic response decay rate. As the neuron model approaches spiking threshold from below, we track membrane voltage fluctuation statistics of numerically simulated stochastic trajectories. We find that for a given distance from spiking threshold, increasing the magnitude of anesthetic-induced inhibition is associated with augmented signatures of critical slowing: fluctuation amplitudes and correlation times grow as spectral power is increasingly focused at 0 Hz. Furthermore, as a function of distance from threshold, anesthesia significantly modifies the power-law exponents for variance and correlation time divergences observable in stochastic trajectories. Compared to the inverse square root power-law scaling of these quantities anticipated for the saddle-node bifurcation of type-I neurons in the absence of anesthesia, increasing anesthetic-induced inhibition results in an observable exponent <-0.5 for variance and >-0.5 for correlation time divergences. However, these behaviors eventually break down as distance from threshold goes to zero with both the variance and correlation time converging to common values independent of anesthesia. Compared to the case of no synaptic input, linearization of an approximating multivariate Ornstein-Uhlenbeck model reveals these effects to be the consequence of an additional slow eigenvalue associated with synaptic activity that competes with those of the underlying point neuron in a manner that depends on distance from spiking threshold.

Recording Brain Electromagnetic Activity During the Administration of the Gaseous Anesthetic Agents Xenon and Nitrous Oxide in Healthy Volunteers

Anesthesia arguably provides one of the only systematic ways to study the neural correlates of global consciousness/unconsciousness. However to date most neuroimaging or neurophysiological investigations in humans have been confined to the study of γ-Amino-Butyric-Acid-(GABA)-receptor-agonist-based anesthetics, while the effects of dissociative N-Methyl-D-Aspartate-(NMDA)-receptor-antagonist-based anesthetics ketamine, nitrous oxide (N2O) and xenon (Xe) are largely unknown. This paper describes the methods underlying the simultaneous recording of magnetoencephalography (MEG) and electroencephalography (EEG) from healthy males during inhalation of the gaseous anesthetic agents N2O and Xe. Combining MEG and EEG data enables the assessment of electromagnetic brain activity during anesthesia at high temporal, and moderate spatial, resolution. Here we describe a detailed protocol, refined over multiple recording sessions, that includes subject recruitment, anesthesia equipment setup in the MEG scanner room, data collection and basic data analysis. In this protocol each participant is exposed to varying levels of Xe and N2O in a repeated measures cross-over design. Following relevant baseline recordings participants are exposed to step-wise increasing inspired concentrations of Xe and N2O of 8, 16, 24 and 42%, and 16, 32 and 47% respectively, during which their level of responsiveness is tracked with an auditory continuous performance task (aCPT). Results are presented for a number of recordings to highlight the sensor-level properties of the raw data, the spectral topography, the minimization of head movements, and the unequivocal level dependent effects on the auditory evoked responses. This paradigm describes a general approach to the recording of electromagnetic signals associated with the action of different kinds of gaseous anesthetics, which can be readily adapted to be used with volatile and intravenous anesthetic agents. It is expected that the method outlined can contribute to the understanding of the macro-scale mechanisms of anesthesia by enabling methodological extensions involving source space imaging and functional network analysis.

Transitions between asynchronous and synchronous states: a theory of correlations in small neural circuits

The study of correlations in neural circuits of different size, from the small size of cortical microcolumns to the large-scale organization of distributed networks studied with functional imaging, is a topic of central importance to systems neuroscience. However, a theory that explains how the parameters of mesoscopic networks composed of a few tens of neurons affect the underlying correlation structure is still missing. Here we consider a theory that can be applied to networks of arbitrary size with multiple populations of homogeneous fully-connected neurons, and we focus its analysis to a case of two populations of small size. We combine the analysis of local bifurcations of the dynamics of these networks with the analytical calculation of their cross-correlations. We study the correlation structure in different regimes, showing that a variation of the external stimuli causes the network to switch from asynchronous states, characterized by weak correlation and low variability, to synchronous states characterized by strong correlations and wide temporal fluctuations. We show that asynchronous states are generated by strong stimuli, while synchronous states occur through critical slowing down when the stimulus moves the network close to a local bifurcation. In particular, strongly positive correlations occur at the saddle-node and Andronov-Hopf bifurcations of the network, while strongly negative correlations occur when the network undergoes a spontaneous symmetry-breaking at the branching-point bifurcations. These results show how the correlation structure of firing-rate network models is strongly modulated by the external stimuli, even keeping the anatomical connections fixed. These results also suggest an effective mechanism through which biological networks may dynamically modulate the encoding and integration of sensory information.

Human seizures couple across spatial scales through travelling wave dynamics

Epilepsy—the propensity toward recurrent, unprovoked seizures—is a devastating disease affecting 65 million people worldwide. Understanding and treating this disease remains a challenge, as seizures manifest through mechanisms and features that span spatial and temporal scales. Here we address this challenge through the analysis and modelling of human brain voltage activity recorded simultaneously across microscopic and macroscopic spatial scales. We show that during seizure large-scale neural populations spanning centimetres of cortex coordinate with small neural groups spanning cortical columns, and provide evidence that rapidly propagating waves of activity underlie this increased inter-scale coupling. We develop a corresponding computational model to propose specific mechanisms—namely, the effects of an increased extracellular potassium concentration diffusing in space—that support the observed spatiotemporal dynamics. Understanding the multi-scale, spatiotemporal dynamics of human seizures—and connecting these dynamics to specific biological mechanisms—promises new insights to treat this devastating disease.

The authors record both local and long-range neural activity during human epileptic seizures to study the underlying multi-scale dynamics. They find that coupling of activity across spatial scales increases during seizures through propagating waves that are fit by a model that combines neural activity and potassium concentration dynamics.

A neural mass model of spontaneous burst suppression and epileptic seizures

The paper presents a neural mass model that is capable of simulating the transition to and from various forms of paroxysmal activity such as burst suppression and epileptic seizure-like waveforms. These events occur without changing parameters in the model. The model is based on existing neural mass models, with the addition of feedback of fast dynamics to create slowly time varying parameters, or slow states. The goal of this research is to establish a link between system properties that modulate neural activity and the fast changing dynamics, such as membrane potentials and firing rates that can be manipulated using electrical stimulation. Establishing this link is likely to be a necessary component of a closed-loop system for feedback control of pathological neural activity.

How the cortico-thalamic feedback affects the EEG power spectrum over frontal and occipital regions during propofol-induced sedation

Increasing concentrations of the anaesthetic agent propofol initially induces sedation before achieving full general anaesthesia. During this state of anaesthesia, the observed specific changes in electroencephalographic (EEG) rhythms comprise increased activity in the δ- (0.5-4 Hz) and α- (8-13 Hz) frequency bands over the frontal region, but increased δ- and decreased α-activity over the occipital region. It is known that the cortex, the thalamus, and the thalamo-cortical feedback loop contribute to some degree to the propofol-induced changes in the EEG power spectrum. However the precise role of each structure to the dynamics of the EEG is unknown. In this paper we apply a thalamo-cortical neuronal population model to reproduce the power spectrum changes in EEG during propofol-induced anaesthesia sedation. The model reproduces the power spectrum features observed experimentally both in frontal and occipital electrodes. Moreover, a detailed analysis of the model indicates the importance of multiple resting states in brain activity. The work suggests that the α-activity originates from the cortico-thalamic relay interaction, whereas the emergence of δ-activity results from the full cortico-reticular-relay-cortical feedback loop with a prominent enforced thalamic reticular-relay interaction. This model suggests an important role for synaptic GABAergic receptors at relay neurons and, more generally, for the thalamus in the generation of both the δ- and the α- EEG patterns that are seen during propofol anaesthesia sedation.

Emergence of spatially heterogeneous burst suppression in a neural field model of electrocortical activity

Burst suppression in the electroencephalogram (EEG) is a well-described phenomenon that occurs during deep anesthesia, as well as in a variety of congenital and acquired brain insults. Classically it is thought of as spatially synchronous, quasi-periodic bursts of high amplitude EEG separated by low amplitude activity. However, its characterization as a “global brain state” has been challenged by recent results obtained with intracranial electrocortigraphy. Not only does it appear that burst suppression activity is highly asynchronous across cortex, but also that it may occur in isolated regions of circumscribed spatial extent. Here we outline a realistic neural field model for burst suppression by adding a slow process of synaptic resource depletion and recovery, which is able to reproduce qualitatively the empirically observed features during general anesthesia at the whole cortex level. Simulations reveal heterogeneous bursting over the model cortex and complex spatiotemporal dynamics during simulated anesthetic action, and provide forward predictions of neuroimaging signals for subsequent empirical comparisons and more detailed characterization. Because burst suppression corresponds to a dynamical end-point of brain activity, theoretically accounting for its spatiotemporal emergence will vitally contribute to efforts aimed at clarifying whether a common physiological trajectory is induced by the actions of general anesthetic agents. We have taken a first step in this direction by showing that a neural field model can qualitatively match recent experimental data that indicate spatial differentiation of burst suppression activity across cortex.

Anesthetic action on extra-synaptic receptors: effects in neural population models of EEG activity

The role of extra-synaptic receptors in the regulation of excitation and inhibition in the brain has attracted increasing attention. Because activity in the extra-synaptic receptors plays a role in regulating the level of excitation and inhibition in the brain, they may be important in determining the level of consciousness. This paper reviews briefly the literature on extra-synaptic GABA and NMDA receptors and their affinity to anesthetic drugs. We propose a neural population model that illustrates how the effect of the anesthetic drug propofol on GABAergic extra-synaptic receptors results in changes in neural population activity and the electroencephalogram (EEG). Our results show that increased tonic inhibition in inhibitory cortical neurons cause a dramatic increase in the power of both δ− and α− bands. Conversely, the effects of increased tonic inhibition in cortical excitatory neurons and thalamic relay neurons have the opposite effect and decrease the power in these bands. The increased δ-activity is in accord with observed data for deepening propofol anesthesia; but is absolutely dependent on the inclusion of extrasynaptic (tonic) GABA action in the model.

Emergence from general anesthesia and the sleep-manifold

The electroencephalogram (EEG) during the re-establishment of consciousness after general anesthesia and surgery varies starkly between patients. Can the EEG during this emergence period provide a means of estimating the underlying biological processes underpinning the return of consciousness? Can we use a model to infer these biological processes from the EEG patterns? A frontal EEG was recorded from 84 patients. Ten patients were chosen for state-space analysis. Five showed archetypal emergences; which consisted of a progressive decrease in alpha power and increase peak alpha frequency before return of responsiveness. The five non-archetypal emergences showed almost no spectral EEG changes (even as the volatile general anesthetic decreased) and then an abrupt return of responsiveness. We used Bayesian methods to estimate the likelihood of an EEG pattern corresponding to the position of the patient on a 2-dimensional manifold in a state space of excitatory connection strength vs. change in intrinsic resting neuronal membrane conductivity. We could thus visualize the trajectory of each patient in the state-space during their emergence period. The patients who followed an archetypal emergence displayed a very consistent pattern; consisting of progressive increase in conductivity, and a temporary period of increased connection strength before return of responsiveness. The non-archetypal emergence trajectories remained fixed in a region of phase space characterized by a relatively high conductivity and low connection strength throughout emergence. This unexpected progressive increase in conductivity during archetypal emergence may be due to an abating of the surgical stimulus during this period. Periods of high connection strength could represent forays into dissociated consciousness, but the model suggests all patients reposition near the fold in the state space to take advantage of bi-stable cortical dynamics before transitioning to consciousness.

Study of GABAergic extra-synaptic tonic inhibition in single neurons and neural populations by traversing neural scales: application to propofol-induced anaesthesia

Anaesthetic agents are known to affect extra-synaptic GABAergic receptors, which induce tonic inhibitory currents. Since these receptors are very sensitive to small concentrations of agents, they are supposed to play an important role in the underlying neural mechanism of general anaesthesia. Moreover anaesthetic agents modulate the encephalographic activity (EEG) of subjects and hence show an effect on neural populations. To understand better the tonic inhibition effect in single neurons on neural populations and hence how it affects the EEG, the work considers single neurons and neural populations in a steady-state and studies numerically and analytically the modulation of their firing rate and nonlinear gain with respect to different levels of tonic inhibition. We consider populations of both type-I (Leaky Integrate-and-Fire model) and type-II (Morris-Lecar model) neurons. To bridge the single neuron description to the population description analytically, a recently proposed statistical approach is employed which allows to derive new analytical expressions for the population firing rate for type-I neurons. In addition, the work shows the derivation of a novel transfer function for type-I neurons as considered in neural mass models and studies briefly the interaction of synaptic and extra-synaptic inhibition. We reveal a strong subtractive and divisive effect of tonic inhibition in type-I neurons, i.e. a shift of the firing rate to higher excitation levels accompanied by a change of the nonlinear gain. Tonic inhibition shortens the excitation window of type-II neurons and their populations while maintaining the nonlinear gain. The gained results are interpreted in the context of recent experimental findings under propofol-induced anaesthesia.

Computer simulations of synchrony and oscillations evoked by two coherent inputs

Coherent oscillations have been reported in multiple cortical areas. This study examines the characteristics of output spikes through computer simulations when the neural network model receives periodic/aperiodic spatiotemporal spikes with modulated/constant populational activity from two pathways. Synchronous oscillations which have the same period as the input are observed in response to periodic input patterns regardless of populational activity. The results confirm that the output frequency of synchrony is essentially determined by the period of the repeated input patterns. On the other hand, weak periodic outputs are observed when aperiodic spikes are input with modulated populational activity. In this case, higher firing rates are necessary to input for higher frequency oscillations. The spike-timing-dependent plasticity suppresses the spikes which do not contribute to the synchrony for periodic inputs. This effect corresponds to the experimental reports that learning sharpens the synchrony in the motor cortex. These results suggest that spatiotemporal spike patterns should be entrained on modulated populational activity to transmit oscillatory information effectively in the convergent pathway.

Ketamine, Propofol, and the EEG: A Neural Field Analysis of HCN1-Mediated Interactions

Ketamine and propofol are two well-known, powerful anesthetic agents, yet at first sight this appears to be their only commonality. Ketamine is a dissociative anesthetic agent, whose main mechanism of action is considered to be N-methyl-d-aspartate (NMDA) antagonism; whereas propofol is a general anesthetic agent, which is assumed to primarily potentiate currents gated by γ-aminobutyric acid type A (GABA_A) receptors. However, several experimental observations suggest a closer relationship. First, the effect of ketamine on the electroencephalogram (EEG) is markedly changed in the presence of propofol: on its own ketamine increases θ (4–8 Hz) and decreases α (8–13 Hz) oscillations, whereas ketamine induces a significant shift to beta band frequencies (13–30 Hz) in the presence of propofol. Second, both ketamine and propofol cause inhibition of the inward pacemaker current I_h, by binding to the corresponding hyperpolarization-activated cyclic nucleotide-gated potassium channel 1 (HCN1) subunit. The resulting effect is a hyperpolarization of the neuron’s resting membrane potential. Third, the ability of both ketamine and propofol to induce hypnosis is reduced in HCN1-knockout mice. Here we show that one can theoretically understand the observed spectral changes of the EEG based on HCN1-mediated hyperpolarizations alone, without involving the supposed main mechanisms of action of these drugs through NMDA and GABA_A, respectively. On the basis of our successful EEG model we conclude that ketamine and propofol should be antagonistic to each other in their interaction at HCN1 subunits. Such a prediction is in accord with the results of clinical experiment in which it is found that ketamine and propofol interact in an infra-additive manner with respect to the endpoints of hypnosis and immobility.

Genuine and spurious phase synchronization strengths during consciousness and general anesthesia

Spectral content in a physiological dataset of finite size has the potential to produce spurious measures of coherence. This is especially true for electroencephalography (EEG) during general anesthesia because of the significant alteration of the power spectrum. In this study we quantitatively evaluated the genuine and spurious phase synchronization strength (PSS) of EEG during consciousness, general anesthesia, and recovery. A computational approach based on the randomized data method was used for evaluating genuine and spurious PSS. The validity of the method was tested with a simulated dataset. We applied this method to the EEG of normal subjects undergoing general anesthesia and investigated the finite size effects of EEG references, data length and spectral content on phase synchronization. The most influential factor for genuine PSS was the type of EEG reference; the most influential factor for spurious PSS was the spectral content. Genuine and spurious PSS showed characteristic temporal patterns for each frequency band across consciousness and anesthesia. Simultaneous measurement of both genuine and spurious PSS during general anesthesia is necessary in order to avoid incorrect interpretations regarding states of consciousness.

Cortico-cardio-respiratory network interactions during anesthesia

General anesthetics are used during medical and surgical procedures to reversibly induce a state of total unconsciousness in patients. Here, we investigate, from a dynamic network perspective, how the cortical and cardiovascular systems behave during anesthesia by applying nonparametric spectral techniques to cortical electroencephalography, electrocardiogram and respiratory signals recorded from anesthetized rats under two drugs, ketamine-xylazine (KX) and pentobarbital (PB). We find that the patterns of low-frequency cortico-cardio-respiratory network interactions may undergo significant changes in network activity strengths and in number of network links at different depths of anesthesia dependent upon anesthetics used.

Distributed control in a mean-field cortical network model: implications for seizure suppression

Brain electrical stimulation (BES) has long been suggested as a means of controlling pathological brain activity. In epilepsy, control of a spatially localized source, the seizure focus, may normalize neuronal dynamics. Consequently, most BES research has been directed at controlling small, local, neuronal populations. At a higher level, pathological seizure activity can be viewed as a network event that may begin without a clear spatial focus or in multiple sites and spread rapidly through a distributed cortical network. In this paper, we begin to address the implications of local control in a network scenario. To do so, we explore the efficacy of local BES when deployed over a larger-scale neuronal network, for instance, using a grid of stimulating electrodes on the cortex. By introducing a mean-field model of neuronal interactions we are able to identify limitations in network controllability based on physiological constraints that suggest the need for more nuanced network control strategies.

Modeling effect of GABAergic current in a basal ganglia computational model

Electrical high frequency stimulation (HFS) of deep brain regions is a method shown to be clinically effective in different types of movement and neurological disorders. In order to shed light on its mode of action a computational model of the basal ganglia network coupled the HFS as injection current into the cells of the subthalamic nucleus (STN). Its overall increased activity rendered a faithful transmission of sensorimotor input through thalamo-cortical relay cells possible. Our contribution uses this model by Rubin and Terman (J Comput Neurosci, 16, 211-223, 2004) as a starting point and integrates recent findings on the importance of the extracellular concentrations of the inhibiting neurotransmitter GABA. We are able to show in this computational study that besides electrical stimulation a high concentration of GABA and its resulting conductivity in STN cells is able to re-establish faithful thalamocortical relaying, which otherwise broke down in the simulated parkinsonian state.

Understanding the effects of anesthetic agents on the EEG through neural field theory

Anesthetic and analgesic agents act through a diverse range of pharmacological mechanisms. Existing empirical data clearly shows that such "microscopic" pharmacological diversity is reflected in their "macroscopic" effects on the human electroencephalogram (EEG). Based on a detailed mesoscopic neural field model we theoretically posit that anesthetic induced EEG activity is due to selective parametric changes in synaptic efficacy and dynamics. Specifically, on the basis of physiologically constrained modeling, it is speculated that the selective modification of inhibitory or excitatory synaptic activity may differentially effect the EEG spectrum. Such results emphasize the importance of neural field theories of brain electrical activity for elucidating the principles whereby pharmacological agents effect the EEG. Such insights will contribute to improved methods for monitoring depth of anesthesia using the EEG.

Dynamical changes in neurological diseases and anesthesia

Dynamics of neuronal networks can be altered in at least two ways: by changes in connectivity, that is, the physical architecture of the network, or changes in the amplitudes and kinetics of the intrinsic and synaptic currents within and between the elements making up a network. We argue that the latter changes are often overlooked as sources of alterations in network behavior when there are also structural (connectivity) abnormalities present; indeed, they may even give rise to the structural changes observed in these states. Here we look at two clinically relevant states (Parkinson's disease and schizophrenia) and argue that non-structural changes are important in the development of abnormal dynamics within the networks known to be relevant to each disorder. We also discuss anesthesia, since it is entirely acute, thus illustrating the potent effects of changes in synaptic and intrinsic membrane currents in the absence of structural alteration. In each of these, we focus on the role of changes in GABAergic function within microcircuits, stressing literature within the last few years.

Effects of anesthetic agents on brain blood oxygenation level revealed with ultra-high field MRI

During general anesthesia it is crucial to control systemic hemodynamics and oxygenation levels. However, anesthetic agents can affect cerebral hemodynamics and metabolism in a drug-dependent manner, while systemic hemodynamics is stable. Brain-wide monitoring of this effect remains highly challenging. Because T₂*-weighted imaging at ultra-high magnetic field strengths benefits from a dramatic increase in contrast to noise ratio, we hypothesized that it could monitor anesthesia effects on brain blood oxygenation. We scanned rat brains at 7T and 17.2T under general anesthesia using different anesthetics (isoflurane, ketamine-xylazine, medetomidine). We showed that the brain/vessels contrast in T₂*-weighted images at 17.2T varied directly according to the applied pharmacological anesthetic agent, a phenomenon that was visible, but to a much smaller extent at 7T. This variation is in agreement with the mechanism of action of these agents. These data demonstrate that preclinical ultra-high field MRI can monitor the effects of a given drug on brain blood oxygenation level in the absence of systemic blood oxygenation changes and of any neural stimulation.

Effects of transcranial Direct Current Stimulation (tDCS) on cortical activity: a computational modeling study

Although it is well-admitted that transcranial Direct Current Stimulation (tDCS) allows for interacting with brain endogenous rhythms, the exact mechanisms by which externally-applied fields modulate the activity of neurons remain elusive. In this study a novel computational model (a neural mass model including subpopulations of pyramidal cells and inhibitory interneurons mediating synaptic currents with either slow or fast kinetics) of the cerebral cortex was elaborated to investigate the local effects of tDCS on neuronal populations based on an in-vivo experimental study. Model parameters were adjusted to reproduce evoked potentials (EPs) recorded from the somatosensory cortex of the rabbit in response to air-puffs applied on the whiskers. EPs were simulated under control condition (no tDCS) as well as under anodal and cathodal tDCS fields. Results first revealed that a feed-forward inhibition mechanism must be included in the model for accurate simulation of actual EPs (peaks and latencies). Interestingly, results revealed that externally-applied fields are also likely to affect interneurons. Indeed, when interneurons get polarized then the characteristics of simulated EPs become closer to those of real EPs. In particular, under anodal tDCS condition, more realistic EPs could be obtained when pyramidal cells were depolarized and, simultaneously, slow (resp. fast) interneurons became de- (resp. hyper-) polarized. Geometrical characteristics of interneurons might provide some explanations for this effect.

Dynamic causal models and physiological inference: a validation study using isoflurane anaesthesia in rodents

Generative models of neuroimaging and electrophysiological data present new opportunities for accessing hidden or latent brain states. Dynamic causal modeling (DCM) uses Bayesian model inversion and selection to infer the synaptic mechanisms underlying empirically observed brain responses. DCM for electrophysiological data, in particular, aims to estimate the relative strength of synaptic transmission at different cell types and via specific neurotransmitters. Here, we report a DCM validation study concerning inference on excitatory and inhibitory synaptic transmission, using different doses of a volatile anaesthetic agent (isoflurane) to parametrically modify excitatory and inhibitory synaptic processing while recording local field potentials (LFPs) from primary auditory cortex (A1) and the posterior auditory field (PAF) in the auditory belt region in rodents. We test whether DCM can infer, from the LFP measurements, the expected drug-induced changes in synaptic transmission mediated via fast ionotropic receptors; i.e., excitatory (glutamatergic) AMPA and inhibitory GABA(A) receptors. Cross- and auto-spectra from the two regions were used to optimise three DCMs based on biologically plausible neural mass models and specific network architectures. Consistent with known extrinsic connectivity patterns in sensory hierarchies, we found that a model comprising forward connections from A1 to PAF and backward connections from PAF to A1 outperformed a model with forward connections from PAF to A1 and backward connections from A1 to PAF and a model with reciprocal lateral connections. The parameter estimates from the most plausible model indicated that the amplitude of fast glutamatergic excitatory postsynaptic potentials (EPSPs) and inhibitory postsynaptic potentials (IPSPs) behaved as predicted by previous neurophysiological studies. Specifically, with increasing levels of anaesthesia, glutamatergic EPSPs decreased linearly, whereas fast GABAergic IPSPs displayed a nonlinear (saturating) increase. The consistency of our model-based in vivo results with experimental in vitro results lends further validity to the capacity of DCM to infer on synaptic processes using macroscopic neurophysiological data.

Effects of the anesthetic agent propofol on neural populations

The neuronal mechanisms of general anesthesia are still poorly understood. Besides several characteristic features of anesthesia observed in experiments, a prominent effect is the bi-phasic change of power in the observed electroencephalogram (EEG), i.e. the initial increase and subsequent decrease of the EEG-power in several frequency bands while increasing the concentration of the anaesthetic agent. The present work aims to derive analytical conditions for this bi-phasic spectral behavior by the study of a neural population model. This model describes mathematically the effective membrane potential and involves excitatory and inhibitory synapses, excitatory and inhibitory cells, nonlocal spatial interactions and a finite axonal conduction speed. The work derives conditions for synaptic time constants based on experimental results and gives conditions on the resting state stability. Further the power spectrum of Local Field Potentials and EEG generated by the neural activity is derived analytically and allow for the detailed study of bi-spectral power changes. We find bi-phasic power changes both in monostable and bistable system regime, affirming the omnipresence of bi-spectral power changes in anesthesia. Further the work gives conditions for the strong increase of power in the δ-frequency band for large propofol concentrations as observed in experiments.

Large-scale neural dynamics: simple and complex

We review the use of neural field models for modelling the brain at the large scales necessary for interpreting EEG, fMRI, MEG and optical imaging data. Albeit a framework that is limited to coarse-grained or mean-field activity, neural field models provide a framework for unifying data from different imaging modalities. Starting with a description of neural mass models, we build to spatially extend cortical models of layered two-dimensional sheets with long range axonal connections mediating synaptic interactions. Reformulations of the fundamental non-local mathematical model in terms of more familiar local differential (brain wave) equations are described. Techniques for the analysis of such models, including how to determine the onset of spatio-temporal pattern forming instabilities, are reviewed. Extensions of the basic formalism to treat refractoriness, adaptive feedback and inhomogeneous connectivity are described along with open challenges for the development of multi-scale models that can integrate macroscopic models at large spatial scales with models at the microscopic scale.

The electrocortical effects of enflurane: experiment and theory

BACKGROUND

High concentrations of enflurane will induce a characteristic electroencephalogram pattern consisting of periods of suppression alternating with large short paroxysmal epileptiform discharges (PEDs). In this study, we compared a theoretical computer model of this activity with real local field potential (LFP) data obtained from anesthetized rats.

METHODS

After implantation of a high-density 8 x 8 electrode array in the visual cortex, the patterns of LFP and multiunit spike activity were recorded in rats during 0.5, 1.0, 1.5, and 2.0 minimum alveolar anesthetic concentration (MAC) enflurane anesthesia. These recordings were compared with computer simulations from a mean field model of neocortical dynamics. The neuronal effect of increasing enflurane concentration was simulated by prolonging the decay time constant of the inhibitory postsynaptic potential (IPSP). The amplitude of the excitatory postsynaptic potential (EPSP) was modulated, inverse to the neocortical firing rate.

RESULTS

In the anesthetized rats, increasing enflurane concentrations consistently caused the appearance of suppression pattern (>1.5 MAC) in the LFP recordings. The mean rate of multiunit spike activity decreased from 2.54/s (0.5 MAC) to 0.19/s (2.0 MAC). At high MAC, the majority of the multiunit action potential events became synchronous with the PED. In the theoretical model, prolongation of the IPSP decay time and activity-dependent EPSP modulation resulted in output that was similar in morphology to that obtained from the experimental data. The propensity for rhythmic seizure-like activity in the model could be determined by analysis of the eigenvalues of the equations.

CONCLUSION

It is possible to use a mean field theory of neocortical dynamics to replicate the PED pattern observed in LFPs in rats under enflurane anesthesia. This pattern requires a combination of a moderately increased total area under the IPSP, prolonged IPSP decay time, and also activity-dependent modulation of EPSP amplitude.