Cefepime-induced encephalopathy: Neural mass modeling of triphasic wave-like generalized periodic discharges with a high negative component (Tri-HNC).

AIM: Cefepime, a fourth-generation cephalosporin, acts as a GABAA receptor antagonist. Cefepime-induced encephalopathy (CIE) is frequently overlooked. We aimed to clarify the clinical features, characteristic electroencephalography (EEG), and mechanisms of CIE to aid in its early recognition.
METHODS: CIE cases documented by a single-center consultation-liaison team between April 2015 and March 2017 were retrospectively reviewed. For further investigation, neural mass modeling was performed in silico.
RESULTS: Three patients with CIE refused medication/examination and showed overt pain, palilalia, and much greater deterioration of eye and verbal response than the motor response, which was possibly related to GABAergic dysfunction. Triphasic wave-like generalized periodic discharges with a high negative component (Tri-HNC) were identified on the EEG of all three cases. The simulation reproduced the characteristic feature of 2-3 Hz Tri-HNC and recovery course on EEG, and a possible involvement of individual differences in pharmacological intervention. It also suggested that auto-inhibition (synaptic inputs from interneuron to interneuron) dysregulation contributed to generating Tri-HNC in CIE.
CONCLUSION: As CIE is iatrogenic and continues unless cefepime is stopped, early recognition is crucial. Physicians should be vigilant about altered mental status, pain, and verbal changes in patients taking cefepime. Tri-HNC on EEG can expedite the diagnosis of CIE, and the association between Tri-HNC and CIE suggests that an excitatory and inhibitory imbalance due to the dysfunction of GABAergic interneurons is the underlying mechanism. This modeling may offer a new method of investigating disorders related to GABAergic dysfunction.

Identifying the primary cause of altered mental status is difficult since various confounders are associated with the patient's condition. Although antibiotics are known to cause delirium, this adverse effect is largely under‐recognized.1 Cefepime, a fourth‐generation cephalosporin, acts as a GABAA receptor antagonist2 and often induces antibiotic‐associated encephalopathy.1 In the intensive care unit, 15% of patients treated with intravenous cefepime for at least 3 days suffer cefepime‐induced encephalopathy (CIE), but this condition is frequently overlooked.3 To aid in its early recognition, we reviewed three cases of CIE presenting unique triphasic wave‐like generalized periodic discharges with a high negative component (Tri‐HNC).

Generalized periodic discharges (GPD) are commonly encountered in metabolic encephalopathy and cerebral hypoxia4 and can be simulated using neural mass modeling.5, 6 Synaptic failure or neuronal damage of inhibitory interneurons presumably plays a critical role.4, 7 To improve our understanding of the mechanism underlying CIE, which is presumably related to GABAergic interneuron dysfunction, we performed neural mass modeling to reproduce the Tri‐HNC in silico.

Methods

Clinical data collection and ethical considerations

Three patients with CIE were referred to the consultation–liaison (C‐L) team between April 2015 and March 2017 at Tokyo Metropolitan Tama Medical Center, a 789‐bed tertiary care teaching hospital in Tokyo, Japan. The C‐L service, which is available around‐the‐clock, receives approximately 850 patients annually. EEG were recorded for all three patients. Twenty‐one electrodes (13 out of 21 for a simplified recording) were placed on the scalp in accordance with the International 10–20 system. The time constant and the high‐cut filter were set to 0.3 s and 60 Hz, respectively. We obtained written informed consent from Patients 1 and 2, and the study was approved by the institutional review board at Tokyo Metropolitan Tama Medical Center with a waiver of informed consent from Patient 3 (due to death).

Neural mass modeling

In the simulation, we used a spatially homogeneous version of neural mass modeling, which is basically a mean field model of the cortex.4, 5, 6, 8 Figure 1a,b present an overview of the model. The model describes the dynamics of the average membrane potentials of a cortical macrocolumn comprised of excitatory (pyramidal) neurons and GABAergic inhibitory (inter‐)neurons. In short, both excitatory and inhibitory neurons receive input via intracortical synaptic projections as well as nonspecific excitatory input from regions not explicitly incorporated into the model, such as the thalamus. The excitatory neurons excite both themselves and the inhibitory interneurons through glutamate‐mediated synapses. The interneurons inhibit both themselves and the excitatory neurons through GABA‐mediated synapses.4, 7

Figure 1

An overview of the mean field model, simulation, and analysis. (a) From the cortical circuitry comprised of excitatory and inhibitory neurons, (b) only one excitatory neuron (Ex) and one inhibitory neuron (Inh) are extracted for the mean field model. Both excitatory and inhibitory neurons receive input via intracortical synaptic projections (Ex‐>Ex, Ex‐>Inh, Inh‐>Ex and Inh‐>Inh) as well as nonspecific excitatory input from regions not explicitly incorporated into the model, such as the thalamus (p ± σ). Nkx denotes the number of synaptic inputs from k to x. The indices e and i refer to the excitatory, and inhibitory neurons, respectively. Using the values of Nkx, the model generates the mean voltage of the excitatory and inhibitory neurons. (c) The mean voltage of the excitatory neuron (Nex = 3000, Nix = 380). (d) To illustrate the simulated scalp electroencephalography (EEG), the plus and minus ends were simply reversed and arbitrary units were used for the y‐axis. (e) Using normalized EEG, (f) negative/positive ratio, complex frequency, and complex duration were calculated. (g) The power of the waves and (h) spectrogram are also presented.

Model parameters and equations

The indices e and i used below refer to the excitatory and inhibitory neurons, respectively. The membrane potentials are denoted by Vk(t) for k = e,i. Their dynamics are governed by the following set of differential equations: where and denote the membrane time‐constants and resting‐potentials, respectively, and Ikx is proportional to the current flowing into population x due to the activity of population k. The currents Ine and Ini model the afferent, non‐specific input to the cortical column, and are modeled as uncorrelated white‐noise processes with the mean values pne and pni and the standard deviations and , respectively. The currents Ikx are given bywhere h is the response function of AMPA/GABAergic receptors located on the dendrites of neurons within population x. Nkx denotes the number of synaptic inputs from k to x, which is a variable in this study. The function Sk relates the membrane potential of population k to its firing‐rate, and is given by

The variables are dimensionless and model the dependence of AMPA/GABAergic synaptic conductance on the membrane potential of the post‐synaptic neural population k. They are given by

The baseline values for all model parameters were modified from previous reports5 and are listed in Table 1. For further details, see Tjepkema‐Cloostermans et al.6

Table 1

Model parameters, their symbols, and baseline values

Parameter	Symbol	Baseline value
Maximum spike‐rate	Qkmax	500 s−1
Spike‐threshold	Vkspike	−60 mV
SD of spike‐threshold	σk	5 mV
Synaptic efficacy	H GABA, H AMPA	0.71 mV
Reversal potential	E GABA, E AMPA	−90, 40 mV
Number of synaptic inputs from k to x	Nei, Nee, Nie, Nii	3000, 3000, 500, 500 (variable)
Membrane time‐constant	τe, τi	0.094, 0.042 s
Resting potential	Vkrest	−70 mV
Synaptic rate‐constant	γGABA, γAMPA	65, 300 s−1
Non‐specific firing rate	pne, pni	3460, 5070 s−1
SD of non‐specific fluctuation	σ_p_ne,σ_p_ni	1000, 150 s−1

The indices k and x refer to neural populations of type k and x for e (excitatory‐) and i (inhibitory‐) neurons

Illustrating EEG, in silico analysis, and statistics

With the values of Nkx, the model generated the mean voltage of the excitatory and inhibitory neurons. To illustrate the simulated scalp EEG, the plus and minus ends were simply reversed, and arbitrary units were used for the y‐axis (Fig. 1c,d).

One representative recording of real or simulated EEG was normalized to its amplitude (Fig. 1e). We defined the negative/positive ratio as the amplitude ratio of negative component over positive component of each GPD. We also defined the complex frequency and the complex duration as the number of GPD shown per second and the duration of each GPD, respectively, based on peak‐to‐peak interval times using semi‐automatic peak detection (Fig. 1f). We characterized spectral components based on fast Fourier transform (FFT). In detail, we computed the power spectrum of the whole waveform, to visualize overall characteristics (Fig. 1g), and spectrogram using half‐overlapping 0.4‐s time windows with a hamming window to visualize temporal variation (Fig. 1h) on MATLAB (R2018a, MathWorks, Inc., Natick, MA, USA), with a signal processing toolbox. Negative/positive ratio, complex frequencies, and complex duration were tested using Welch's analysis of variance on MATLAB. P < 0.05 was considered statistically significant.

Results

Case presentations

Case 1 was that of a 64‐year‐old woman on dialysis. Four days after cefepime was started, she suddenly began refusing medication and repeatedly complained about feeling pain with ‘ouch, ouch, ouch…’. Delirium was suspected, and the patient was referred to our C‐L team. Her Glasgow Coma Scale was E2V3M4, and her EEG showed 2–3 Hz Tri‐HNC (Day 1, Fig. 2a). Other diseases causing altered mental status were excluded as possible by close examinations. We recommended that cefepime be stopped and that the patient be observed closely without any additional interventions, such as benzodiazepines or anti‐epileptic drugs. Four days after the consultation, her mental status normalized, and her EEG showed an alpha rhythm with waves of various frequencies (Day 5, Fig. 2b). After she became fully alert, a follow‐up EEG showed no epileptic discharges (Days 12 and 31; Fig. 2c,d, respectively). Figure 2e shows the spectrogram throughout the clinical course. Alpha‐band oscillation recovered clearly by Day 31.

Figure 2

Monopolar electroencephalography (EEG) during the clinical course of cefepime‐induced encephalopathy. (a–d) Case 1, a 64‐year‐old woman on dialysis. (a) Unique 2–3 Hz triphasic wave‐like generalized periodic discharges with a high negative component (Tri‐HNC) are shown (Consultation Day 1). (b) After cessation of cefepime, her mental status recovered, and her EEG showed an alpha rhythm with waves of various frequencies (Day 5). Follow‐up EEG showed no epileptic discharges on (c) Day 12 or (d) Day 31. Scale: 1 s, 50 μV. (e) The spectrograms throughout the clinical course of Case 1. Alpha‐band oscillation had recovered clearly by Day 31. (f–i) Case 2, a 44‐year‐old man. EEG taken on consultation Days 1, 3, 15, and 288 (after discharge) are shown. (j) Case 3, a 94‐year‐old woman. EEG taken on consultation Day 1. All three cases represent Tri‐HNC. (k,l) The spectrograms throughout the clinical courses of Cases 2 and 3.

Case 2 was that of a 44‐year‐old man who experienced E1V2M4 and incontinence and refused to be examined. Case 3 involved a 94‐year‐old woman who experienced E1V2M4 and pain and refused medication. All three patients were cared for in a non‐intensive care unit setting and showed palilalia. Myoclonus was not noticed by patients themselves or staff. Their clinical course is summarized in Table 2, and the EEG in Cases 2 and 3 are shown in Figure 2f–i and Figure 2j, respectively. The spectrograms taken throughout the clinical course in Cases 2 and 3 are shown in Figure 2k,l.

Table 2

Summary of three cases of cefepime‐induced encephalopathy

Age (years)	Sex	Weight (kg)	eGFR	Adjusted for renal function	Indication	Clinical findings	Latency (days)	Intervention	Days to improvement
64	Female	52.7	21 (on dialysis)	No	Infectious endocarditis	AMS (E2V3M4), pain, palilalia, refusal of medication	4	Observation	4
44	Male	40.0	33	No	Pneumonia	AMS (E1V2M4), palilalia, refusal of examination	5	Observation	4
94	Female	38.3	14	Yes	Cholangitis, cellulitis	AMS (E1V2M4), pain, palilalia, refusal of medication	4	DZP, PHT	4

AMS, altered mental status; eGFR, estimated glomerular filtration rate (mL/min/1.73m2); DZP, diazepam 5 mg; PHT, phenytoin 250 mg.

To examine the underlying mechanism, we performed a simulation using neural mass modeling.

Simulation of Tri‐HNC and the recovery course in CIE

We first used the same number for Nii and Nie (represented as Nix) in accordance with the original Liley model.5 To simulate cefepime‐induced neurotoxicity as GABAA receptor antagonism, Nix was decreased from 500 (100% input). As expected, Nix = 380 (76% input) partly explained the characteristic feature of Tri‐HNC (Fig. 3a) while the FFT indicated the dominant frequency to be 6–7 Hz (Fig. 3b). Figure 3c shows a spectrogram corresponding to Figure 3a.

Figure 3

Simulation of triphasic wave‐like generalized periodic discharges with a high negative component (Tri‐HNC) and recovery course on electroencephalography (EEG). (a) Simulated EEG explains the characteristic feature of 2–3 Hz Tri‐HNC (on Nix = 380, 76% input), while (b) fast Fourier transform indicates that the dominant frequency was 6–7 Hz. (c) A spectrogram corresponding to (a). (d) The simulated EEG at Nix = 500 (100% input). (e) Illustration of the power of each frequency component (from left to right: Nix = 380, 390, 400, 420, 440, 460, 480, and 500). The dominant background activity changed from 7 Hz (Nix = 380, 76% input) to 12 Hz (Nix = 500, 100% input). (f) The spectrograms near the phase‐transition point (Nix = 380, 381, 382, 390) qualitatively reproduced the clinical recovery course.

To assess the recovery course, Nix was changed from 380 in a stepwise fashion to 500. Figure 3d shows the simulated EEG at Nix = 500 (100% input). The power for each frequency (from left to right) is illustrated in Figure 3e. The dominant background activity changed from 7 Hz (Nix = 380, 76% input) to 12 Hz (Nix = 500, 100% input). The spectrograms near the phase‐transition point (Fig. 3f) qualitatively reproduced the clinical recovery course (Fig. 2e).

Simulation of pharmacological intervention

Benzodiazepines, such as clonazepam and diazepam, are frequently used for alleviating CIE; however, their effectiveness has not been proven and may vary in individual cases.9, 10, 11, 12, 13 To investigate the effect of pharmacological intervention, we fixed Nix at 380 as a constant and moved Nex (=Nee = Nei) simultaneously from 1500 to 3000, since some interneurons are insensitive to benzodiazepines.14, 15 At Nex = 1500 or 1800 (50% or 60% input, respectively), the simulated EEG showed a low voltage pattern (Fig. 4a,b). Although they showed alpha‐band oscillations at Nex = 2100, 2400, and 2700 (70%, 80%, and 90% input in Fig. 4c–e, respectively), strikingly, they became unstable again at Nex = 3000 (100% input, Fig. 4f), suggesting that pharmacological intervention had transient effects and that the permanent recovery of Nix played a critical role in stability.

Figure 4

Simulation of pharmacological intervention. (a,b) With Nix = 380 as a constant, the simulated electroencephalography (EEG) showed a low‐voltage pattern at Nex = 1500 or 1800 (50% and 60% input, respectively), while (c–e) they showed alpha‐band oscillations at Nex = 2100, 2400, and 2700 (70%, 80%, and 90% input, respectively). (f) Notably, the EEG pattern became unstable again at Nex = 3000 (100% input). These findings indicate individual pharmacological differences.

This simulation was consistent with clinical observations and led us to form the hypothesis that the underlying mechanism involved the induction of low‐voltage EEG activity in the CIE patient by high‐dose benzodiazepine while low‐dose (or an appropriate concentration of) benzodiazepine transiently improved awareness without maintaining full‐recovery until sufficient cefepime clearance.

Excitatory and inhibitory imbalance explains the characteristic features seen on the EEG

We hypothesized that excitatory and inhibitory (E/I) imbalance was essential to explaining the features seen on the EEG. To determine which parameter was chiefly responsible for these characteristic features, we conducted stepwise parameter studies in silico.

First, we set Nix at 500 (100% input, constant) and experimentally changed Nee and Nei individually from 1200 to 6000 (40% to 200% input; Fig. 5a). This simulation suggested that alpha‐band oscillation was robust in the condition that both Nee and Nei were under 100% input and that Nee basically provided the dominant frequency. It also suggested that Nee and Nei balance was a key factor of phase‐transition, as previously reported.6, 7 As expected, the Tri‐HNC patterns can be generated at Nee = 4650 and Nei = 5850 (Fig. 6a), which indicated that increased E/I was essential to generate the Tri‐HNC.

Figure 5

Dominant frequencies and phase‐transition experimentally simulated using each excitatory and inhibitory parameter. (a) Nei from 1200 to 6000 (40% to 200% input) and Nee from 1800 to 4950 (60% to 165% input) under Nix constant (500, 100% input), suggesting that Nee provided the dominant frequency and that Nee and Nei balance was a key factor of phase‐transition. (b) Nie from 340 to 500 (68% to 100% input) and Nii from 290 to 500 (58% to 100% input) under Nex constant (3000, 100% input), suggesting that Nie provided the dominant frequency and that Nie and Nii balance was another key factor of phase‐transition. (c) Nei from 1800 to 6000 (60% to 200% input) and Nii from 430 to 500 (86% to 100% input) under Nee and Nie constant (100% input), indicating that Nii and Nei were insufficient to explain the generalized periodic discharge generation, and that alpha‐band oscillation was robust enough to overcome the external disturbance. Low Vol., low voltage EEG. *Phase‐transition point.

Figure 6

Simulation of other characteristic features and statistics. (a,b) A unique 2–3 Hz triphasic wave‐like generalized periodic discharges with a high negative component (Tri‐HNC) pattern similar to that seen in cefepime‐induced encephalopathy was reproduced with certain parameters. (a) Nee = 4650, Nei = 5850, and Nix = 500. (b) Nex = 3000, Nie = 350, and Nii = 292. (c) Tri‐HNC with another frequency (Nee = 4500, Nei = 5100, Nix = 500) and (d) spindle (Nee = 2700, Nei = 4500, Nix = 500) was also generated. (e) The negative/positive ratios (P = 0.10) did not differ statistically between the clinical cases and simulations, whereas (f) complex frequencies and (g) complex durations differed statistically (P < 0.05). Error bars indicate SD.

Cefepime theoretically antagonizes GABAergic inhibitory inputs. To determine whether Nie or Nii contributed more to generating GPD, we set Nex = 3000 (100% input, constant) and changed Nie and Nii individually from 280 to 500 (56% to 100% input, Fig. 5b). This simulation suggested that Nie provided the dominant frequency and that Nie and Nii balance was another key factor of phase‐transition, indicating that only Nie or Nii was not enough to generate GPD. Moreover, the 2–3 Hz Tri‐HNC pattern was most closely reproduced at Nie = 350 and Nii = 292 (Fig. 6b).

To investigate another possibility that inputs to interneuron is the underlying mechanism, we set Nee = 3000 (100% input, constant) and Nie = 500 (100% input, constant) and changed Nei and Nii individually as shown in Figure 5c; however, the dominant frequencies were restricted to the alpha‐band and GPD were not able to be reproduced with any of the parameters. This finding is similar in Nie and Nii stepwise study under Nei and Nii constant (data not shown). These findings demonstrated that only one parameter was insufficient to explain the GPD and that alpha‐band oscillation was robust enough to overcome the external disturbance.

We also examined a possible application of neural mass modeling. With appropriate parameters, neural mass modeling generated characteristic features, such as Tri‐HNC at another frequency (Fig. 6c) and spindle (Fig. 6d). These simulations suggested that E/I balance accounted for some of the characteristic features seen on the EEG.

Finally, we quantitatively compared the negative/positive ratio, complex frequencies, and complex durations of real and simulated EEG. Negative/positive ratios (Fig. 6e) did not differ statistically, whereas the complex frequencies (Fig. 6f) and the complex durations (Fig. 6g) differed statistically. The findings were based on a small number of cases and await validation by further studies.

Discussion

Clinical presentation

We noticed that 2–3 Hz GPD with a high negative component (above baseline) seemed to characterize cefepime‐induced encephalopathy9, 10, 11, 12, 16, 17, 18 in contrast to the triphasic waves induced by other illnesses, such as hepatic encephalitis10 or Creutzfeldt–Jakob disease.19 Moreover, the triphasic morphology with negative polarity in the dominant phase was perceived as differing from the ‘typical’ triphasic wave by clinicians (odds ratio: 0.09, 0.02–0.22),20 suggesting that the clinicians differentiated the Tri‐HNC from ‘typical’ triphasic wave; however, only a few studies have focused on this characteristic feature.10, 21

We identified Tri‐HNC in all three C‐L cases of CIE with chronic kidney disease (Fig. 2 and Table 2). Latency and days to improvement were almost the same as those previously reported,1, 22 and other diseases causing altered mental status were carefully excluded in all the cases. All the patients showed a similar, characteristic clinical course (Table 2), including refusal of medication/examination, complaints of overt pain and palilalia, and much greater deterioration of eye and verbal response than the motor response. Previous reports showed that patients with CIE presented altered mental status22 more often than seizures23, 24 or psychotic symptoms.1 Although it is still unclear when antibiotic‐associated altered mental status with GPD should be diagnosed as non‐convulsive status epilepticus,25 the characteristic clinical presentations in the present cases led the physicians to assume a psychological cause. Interestingly, previous reports suggested that aphasia was a diagnostic feature of CIE,17, 18, 23 and we hypothesized that palilalia was a mild sign of aphasia in patients with CIE. A previous study showed that memory formation, sensory processing, and pain were cognitive functions associated with GABAergic inhibition.26 Together with the evidence that cefepime acts as a GABAA receptor antagonist, the symptoms that the patients experienced can be summarized to be dysfunctions of GABAergic inhibition.

Neural mass modeling and the relevance of Tri‐HNC detection

Our neural mass modeling reproduced the Tri‐HNC and a characteristic feature seen on EEG in the course of recovery (Fig. 3) and suggested a potential mechanism for individual differences in patients’ responses to pharmacological interventions (Fig. 4). In silico analysis revealed that generating GPD was a circuit‐level phenomenon and offered some hints regarding the relation between alpha‐band oscillation and awareness (Fig. 5). Together with the biological insight, it also suggested that the number of synaptic inputs from interneuron to interneuron (Nii, auto‐inhibition) dysregulation contributed to generating Tri‐HNC in CIE, although auto‐inhibition dysregulation alone cannot explain the whole phenomena.

Based on the simulation of pharmacological intervention (Fig. 5), low‐dose (or an appropriate concentration of) benzodiazepine can transiently improve awareness; however, sufficient cefepime clearance was suggested to be a requirement for full‐recovery. Comparison with a large clinical information dataset is necessary for clarifying whether further pharmacological treatment has long‐term benefit and/or side‐effect.

A previous case report showed that flumazenil, a GABAA receptor antagonist, improved the triphasic waves caused by hepatic encephalopathy.27 We hypothesized that these different types of triphasic wave reflected essentially different activities of the brain. Although our neural mass modeling cannot provide a hypothesis as to which subtype of interneuron, such as PV, SST, or 5HT3aR,28 plays a key role, it emphasized the importance of the auto‐inhibitory system associated with interneuron microcircuits.29

Qualitatively and quantitatively, the Tri‐HNC had a blunt ascent and descent with a high negative component, relatively long complex duration (indicating triphasic waves rather than a generalized spike–wave complex25), and 2–3 Hz complex frequencies (Fig. 6e–g). Of note, we succeeded in simulating the most prominent feature, high negative component, quantitatively (Fig. 6e). The psychological clinical presentations (verbal changes, response to stimuli, and fluctuating attention) and the lack of a need for additional interventions showed that Tri‐HNC differed from epileptic discharges and may require a new designation25 although the matter remains controversial.

To explore other possible applications of neural mass modeling, we simulated the EEG of other disorders related to the dysfunction of GABAergic interneurons. We were able to generate a 0.5–1 Hz Tri‐HNC similar to that seen in some patients with glutamic acid decarboxylase antibody‐associated encephalitis, who showed 0.3–0.5 Hz Tri‐HNC30, 31 (Fig. 6c). It is possible, at any rate, that Tri‐HNC is common to various cases of GABAergic interneuron instability. Unfortunately, our neural mass modeling was not able to reproduce ‘typical’ triphasic waves observed in hepatic encephalopathy, which suggested that some metabolic encephalopathies might not be approximated by the simplified mean field model of synaptic networks. Nonetheless, computational modeling may contribute to clarifying its pathophysiology and offer a new method for investigating disorders related to GABAergic dysfunction.32

Limitations

We are aware of several limitations to this study. First, not all patients referred to the C‐L team for consultation underwent an EEG examination. Although all three patients presented Tri‐HNC, it is possible that other patients with CIE were overlooked as often happens in cases of delirium. Second, the number of cases was too small to draw any firm conclusions about the clinical course. Further, comprehensive, prospective studies are needed for clinical investigation.

Next, neural mass modeling was aimed for an approximation of spatially and temporally homogenous EEG. In the model, the interactions between excitatory and inhibitory neurons are averaged over the whole brain and only one excitatory neuron and one inhibitory neuron with their interactions represent the networks. Relatively complicated concepts, such as the dipole moment or source and sink, are not incorporated into the model; however, Tri‐HNC may be one of the best‐matched patterns in this simplified model given its spatially and temporally homogenous pattern in the real EEG. A more sophisticated model could reproduce the ‘typical’ triphasic wave, while our simplified model is just enough to reproduce Tri‐HNC.

This study identified a connection between the features of real EEG and neural mass modeling‐simulated EEG, suggesting that E/I imbalance due to the dysfunction of GABAergic interneurons was the underlying mechanism of CIE (Fig. S1). As theoretical molecular mechanisms remain speculative at best, further biological experiments are needed to verify our hypothesis.

Conclusion

As CIE is iatrogenic and potentially self‐limiting only when cefepime is stopped, early recognition of the symptoms is crucial. In addition to the US Food and Drug Administration recommendation to be aware of the risk of seizures and to recognize the need for dose adjustments when using cefepime,33 physicians also need to be more vigilant in regards to altered mental status, pain, and verbal changes in patients taking cefepime.

Neural mass modeling reproduced Tri‐HNC, a characteristic feature of EEG in patients with CIE. Together with the clinical presentations, Tri‐HNC can expedite the diagnosis of CIE. The association between CIE and Tri‐HNC suggested that E/I imbalance due to the dysfunction of GABAergic interneurons is the underlying mechanism of CIE. This modeling may offer a new method for investigating disorders related to GABAergic dysfunction.

Disclosure statement

The authors have no conflicts of interest to disclose.

Author contributions

H.T. conceived the study. H.T., Y.H., and N.Y. collected the data. H.T. and N.A. analyzed the data. H.T. drafted the first manuscript. All the authors critically revised the manuscript for intellectual content and approved the final version.

Figure S1. Graphical abstract.

Click here for additional data file.