A hidden Markov model reliably characterizes ketamine-induced spectral dynamics in macaque local field potentials and human electroencephalograms
Authors: Indie C Garwood, Sourish Chakravarty, Jacob Donoghue, Meredith Mahnke, Pegah Kahali, Shubham Chamadia, Oluwaseun Akeju, Earl K Miller, Emery N Brown
Journal: PLOS Computational Biology
Accepted: 2021-07-15T00:00:00Z


Ketamine is an NMDA receptor antagonist commonly used to maintain general anesthesia. At anesthetic doses, ketamine causes high power gamma (25-50 Hz) oscillations alternating with slow-delta (0.1-4 Hz) oscillations. These dynamics are readily observed in local field potentials (LFPs) of non-human primates (NHPs) and electroencephalogram (EEG) recordings from human subjects. However, a detailed statistical analysis of these dynamics has not been reported. We characterize ketamine’s neural dynamics using a hidden Markov model (HMM). The HMM observations are sequences of spectral power in seven canonical frequency bands between 0 to 50 Hz, where power is averaged within each band and scaled between 0 and 1. We model the observations as realizations of multivariate beta probability distributions that depend on a discrete-valued latent state process whose state transitions obey Markov dynamics. Using an expectation-maximization algorithm, we fit this beta-HMM to LFP recordings from 2 NHPs, and separately, to EEG recordings from 9 human subjects who received anesthetic doses of ketamine. Our beta-HMM framework provides a useful tool for experimental data analysis. Together, the estimated beta-HMM parameters and optimal state trajectory revealed an alternating pattern of states characterized primarily by gamma and slow-delta activities. The mean duration of the gamma activity was 2.2s([1.7,2.8]s) and 1.2s([0.9,1.5]s) for the two NHPs, and 2.5s([1.7,3.6]s) for the human subjects. The mean duration of the slow-delta activity was 1.6s([1.2,2.0]s) and 1.0s([0.8,1.2]s) for the two NHPs, and 1.8s([1.3,2.4]s) for the human subjects. Our characterizations of the alternating gamma slow-delta activities revealed five sub-states that show regular sequential transitions. These quantitative insights can inform the development of rhythm-generating neuronal circuit models that give mechanistic insights into this phenomenon and how ketamine produces altered states of arousal.
1 Author summary: Monitoring brain activity during anesthesia can provide insights into the underlying mechanisms of how anesthetics elicit altered states of consciousness. Ketamine, a commonly used anesthetic, is known to cause short duration bursts of high frequency electrophysiological activity in the brain, but the neural mechanisms underlying this activity are not well understood. A key limitation in developing accurate models of the underlying mechanism is a lack of detailed knowledge of the dynamic structure and spectral properties of ketamine-induced oscillations. In this work, we address this limitation by developing a statistical framework to quantify ketamine-induced neural activity. Our framework is based on a hidden Markov model, which assumes that the neural activity switches among discrete states, each of which has its own distinctive probabilistic spectral representation. By estimating this versatile statistical model from electrophysiology data, we generated detailed descriptions of the dynamic properties and oscillatory signatures associated with ketamine-induced neurophysiological states in non-human primates and in human patients. Furthermore, we identified additional ketamine-induced states that have not yet been reported. In summary, our detailed quantitative descriptions of ketamine-induced spectra can aid further developments of neurophysiological mechanistic models of ketamine as well as biomarker discovery for clinical anesthesia monitoring.