Analysis of the Hindriks and van Putten model for propofol anesthesia: Limitations and extensions

We present a detailed analysis of the Hindriks and van Putten thalamocortical mean-field model for propofol anesthesia [NeuroImage 60(23), 2012]. The Hindriks and van Putten (HvP) model predicts increases in delta and alpha power for moderate (up to 130%) prolongation of GABA_A inhibitory response, corresponding to light anesthetic sedation. Our analysis reveals that, for deeper anesthetic effect, the model exhibits an unexpected abrupt jump in cortical activity from a low-firing state to an extremely high-firing stable state (∼250 spikes/s), and remains locked there even at GABA_A prolongations as high as 300% which would be expected to induce full comatose suppression of all firing activity. We demonstrate that this unphysiological behavior can be completely suppressed with appropriate tuning of the parameters controlling the sigmoidal functions that map soma voltage to firing rate for the excitatory and inhibitory neural populations, coupled with elimination of the putative population-dependent anesthetic efficacies introduced in the HvP model. The modifications reported here constrain the anesthetized brain activity into a biologically plausible range in which the cortex now has access to a moderate-firing state ("awake") and a low-firing ("anesthetized") state such that the brain can transition from "awake" to "anesthetized" states at a critical level of drug concentration. The modified HvP model predicts a drug-effect hysteresis in which the drug concentration required for induction is larger than that at emergence. In addition, the revised model shows a decrease in the intensity and frequency of alpha-band fluctuations, transitioning to delta-band dominance, with deepening anesthesia. These predicted drug concentration-dependent changes in EEG dynamics are consistent with clinical reports.