Sedation is integral to the management of critically ill patients requiring mechanical ventilation. It is currently provided by either infusions or intermittent boluses of drugs with hypnotic and analgesic properties.Reference Park, Coursin and Ely 1 Oversedation during critical illness is associated with adverse outcomes, and clinical strategies designed to avoid this can decrease ventilation times and intensive care unit (ICU) length of stay.Reference Gehlbach and Kress 2 , Reference Tonner, Weiler, Paris and Scholz 3 These strategies can also improve clinical outcomes, decrease complications (e.g., ventilator-associated pneumonia), and decrease illness cost.Reference Kress, Pohlman, O’Connor and Hall 4
At present, sedation is usually managed by clinical assessments, often in conjunction with a protocol for adjusting drug doses. This approach has been shown to decrease ventilation times and costs of sedative drugs and is considered best practice in current guidelines.Reference Nasraway, Jacobi, Murray and Lumb 5 Several clinical scales have been proposed and validated for assessing sedation level reliably; these are integral parts of most sedation protocols.Reference Riker, Picard and Fraser 6 - Reference Ely, Truman and Shintani 8 Clinical sedation scales are, however, prone to interrater variability and provide only intermittent information of the patient’s state. The most commonly used assessment scale is the Ramsay scale,Reference Jackson, Proudfoot, Cann and Walsh 9 which controversially has not been validated in an ICU environment.
Electroencephalogram (EEG)-based monitoring devices such as the Bispectral Index (BIS) (Aspect Medical Systems, Needham MA) and the Entropy measurement (GE Healthcare, Helsinki Finland) are used in monitoring the depth of anesthesia in the operating room. The BIS monitor is also intended for use in the ICU, but although correlation between BIS and Ramsay has been established,Reference Mondello, Siliotti and Noto 10 studies have shown the confounding effect of frontal electromyographic (FEMG) activity on both the EntropyReference Walsh, Ramsay, Lapinlampi, Särkelä, Viertiö-Oja and Meriläinen 11 and BISReference Vivien, Di Maria, Ouattara, Langeron, Coriat and Riou 12 in the ICU setting.
In conscious, unparalyzed, and lightly anesthetized patients, painful and/or stressful stimuli are associated with increased FEMG amplitude.Reference Edmonds, Couture, Stolzy and Paloheimo 13 FEMG power may also increase in association with certain sounds in deeply anesthetized or comatose patients who are clinically unresponsive to surgical or electrical stimulation.Reference Edmonds, Couture, Stolzy and Paloheimo 13 Because FEMG is a major confounder for EEG-based algorithms, we hypothesized that measuring FEMG activity would provide a surrogate measure of sedation depth. Because deepening sedation tends to suppress naturally occurring arousals as well as responses to external stimuli,Reference Dutton, Smith, Bennett, Archer and Smith 14 measuring the number and magnitude of these responses should give a surrogate measure of the depth of sedation.
For quantifying FEMG activity in ICU patients, we developed a novel algorithm that analyzes the response patterns in the spontaneous FEMG signal measured from the forehead and outputs a single digit with arbitrary units, the Responsiveness Index (RI), ranging from 0, corresponding to no variability in FEMG activity, to 100, corresponding to strong and variable FEMG activity. The characteristics of the response patterns in the FEMG signal to predefined stimuli were studied as part of the algorithm development. The measurement itself does not use active stimulation because the patient is continuously subjected to several sources of external and internal stimuli in the ICU. External sources of stimulation can include lights, noise, and care-giving procedures, whereas internal stimuli can include pain, stress, and anxiety.
We have previously evaluated the RI performance in critically ill patients recovering from cardiac surgery and shown a better correlation to clinical sedation scores compared with state and response entropy.Reference Walsh, Lapinlampi, Ramsay, Särkelä, Uutela and Viertiö-Oja 15 The current article first describes a study of EMG responses and the physiological background that led to the approach used in RI and then describes the development of the algorithm details using the development data.
Methods
Physiological Background
There is abundant neuroanatomical and neurophysiological support for an interaction of the components of the reticular formation (RF) in the brainstem that are responsible for arousal and pain response and of the motor nuclei in the brainstem that are responsible for activation of skeletal muscles of the forehead.Reference Kandel, Schwartz and Jessell 16 The facial muscles are innervated by the seventh cranial nerve (facial nerve) that emerges from the brainstem between the pons and the medulla. The facial nerve nucleus is situated within the RF in the lower pons. The dorsal aspect of the facial nucleus contributes to innervation of the forehead muscles, which are monitored with surface electrodes used for RI calculation. The third cranial (oculomotor) nucleus, which is intimately associated with the RF in the dorsal midbrain, sends fibers that accompany the seventh nerve in its course. These supply the anterior part of the occipitalis-frontalis, the orbicularis palpebrarum, and the corrugator supercilii.
The rostral pontine RF is also known to regulate somatic motor tone,Reference Jouvet 17 - Reference Hendricks, Morrison and Mann 19 and the progression from wakefulness to rapid eye movement sleep is marked by a decrease of this tone.Reference Wills and Chase 20 The activity of the area is modulated by a number of neurotransmitters, notably gamma aminobutyric acid (GABA).Reference Morales, Engelhardt, Soja, Pereda and Chase 21 - Reference Xi, Morales and Chase 25 The evolution of wakefulness from isoflurane anesthesia has been characterized experimentally by increasing GABA levels in the RF, which increases EMG and decreases EEG activity.Reference Vanini, Watson, Lydic and Baghdoyan 26
Even anesthetized patients may show limb movements as a response to painful stimulation. These reflexes can be spinal, and it is suggested that immobility caused by anesthesia is primarily produced by effect on the spinal cord.Reference Rampil, Mason and Singh 27 , Reference Rampil 28 In contrast, the frontal muscles are controlled by the brainstem and require brain activation for even simple reflexes.
With increasing level of consciousness, cortical processing and the facial muscles expressing the emotional state of the patient are likely to generate additional responses. Startle-type reflexes can be induced by many types of stimulation.Reference Brown, Rothwell, Thompson, Britton, Day and Marsden 29 Finally, when the patient is fully conscious, FEMG activity resulting from speech and other voluntarily generated muscle activations occurs. Our hypothesis is that different levels of sedation are reflected at the threshold of generating frontal muscle movements and therefore the FEMG response amounts will be correlated with the levels of sedation.
Methods
The development of the RI algorithm was divided into three stages: characterizing ICU patients’ FEMG activity related to standardized stimuli, developing the algorithm by an iterative process, and evaluating the algorithm with clinical data. The purpose of the first stage was to test the hypothesis that a standardized vocal stimulus causes a consistent response pattern in the frontal EMG activity of ICU patients. In the second stage, the observed FEMG response characteristics would be used to develop an algorithm that quantifies FEMG activity. The third stage would be an initial test on the clinical usability of the algorithm.
FEMG Response to Vocal Stimulus
To study the characteristics of the frontal FEMG response to stimuli, we used a previously collected data set of 19 ICU patients (11 females, eight males) with a mean age (standard deviation [SD]) of 46.5 (17.1) years. Patients were admitted to the Helsinki Surgical Hospital (Helsinki, Finland) ICU because of hepatic encephalopathy. After receiving written informed consent from the patient or immediate family member, each subject was monitored for 2-7 days. Approval for the study was obtained from the local ethics committee. Monitoring was divided into active and passive periods. During active monitoring periods (~8 hours per day on average) the patients were stimulated once every 3 hours with a prerecorded vocal stimulus. The command “open your eyes” (in Finnish) was used and applied via headphones (Sony MDR-XD200) at an 85-dB intensity. A positive response was recorded if the patient opened his or her eyes within 10 seconds from the start of the stimulus. The stimuli’s time stamps were recorded automatically and the research nurse recorded the response (positive or negative) manually. During passive monitoring periods, the stimuli were not applied.
The FEMG signals were recorded with the GE Datex-Ohmeda Entropy module connected to a GE Datex-Ohmeda S/5 monitor. The primary Entropy sensor was placed on the patient’s forehead symmetrically across the hemispheres approximately 4 cm above the orbital bridge. A secondary Entropy sensor was placed below the first one, with one electrode on the temple and the other above the contralateral eye (Figure 1). Two sensors were used to study the effect of placing the electrodes on different facial muscles. The primary sensor covered the frontalis muscle on both hemispheres. The negative pole of the secondary electrode was placed at the temporalis muscle and the positive pole on the frontalis and the corrugator supercilii. The sampling frequency of the measurement was 400 Hz and the frequency band was 0.5-150 Hz.
Placement of primary and secondary Entropy sensors in the recordings. Numbers 1 and 2 denote the active electrodes; G represents the ground electrode.
Epochs starting 2 minutes before each stimulus and ending 2 minutes after the stimulus were extracted from the FEMG waveform. A visual inspection was performed to remove epochs that contained artifacts, including poor electrode contacts and impedance tests performed by either one of the sensors. All epochs containing artifacts or more than 20 seconds of missing data were removed before analysis. From a total of 415 assessments, the research nurse classified 244 as positive response and 171 as no response.
An estimate of the FEMG power signal was calculated by dividing each epoch to 0.5 s long windows. Root mean square power from the 65-95 Hz frequency band was calculated using the fast-Fourier transform. The frequency band was chosen conservatively to easily avoid the mains frequencies (and multiples) and to contain dominantly FEMG activity.
The median and 25th and 75th quartile FEMG powers were analyzed separately for patients who produced a clinically observable response (opened eyes within 10 seconds from the stimulus) and for those who did not respond to the stimulus.
By observing individual FEMG power curves, we discovered that the latency from the start of the stimulus to the FEMG response onset differed within and between patients. For comparing the FEMG responses, we defined a feature-aligned FEMG response by matching the individual responses by response onset rather than the stimulus onset. Response onset time was defined as the first time point after the stimulus when FEMG power increased to a value greater than the 20-second mean FEMG power plus 3 SDs before the stimulus. Only the cases in which the response onset occurred within 30 seconds or less after the stimulus were included. Feature-aligned median and 25th percentile and 75th percentile powers were calculated to describe the form of a typical FEMG response. Of the 244 assessments classified as positive response, 53 were rejected, leaving 191 assessments for the feature-aligned analysis. Wilcoxon signed-rank test was performed to a median power before and after stimulus to analyze the significance of differences in the power levels.
Responsiveness Algorithm Description
In the second stage, the RI algorithm was developed following an iterative process using real ICU patient data (described in the following section). The starting point was to detect patterns resembling the FEMG response patterns related to the vocal stimulus and quantify the amount and magnitude of those patterns over time. The performance of the algorithm was measured by comparing the output with the clinical sedation scores (the modified Ramsay score) in the development data set. Based on the performance metrics in the development data, the algorithm was tuned by iterating the algorithm parameters.
Development Data Set
The development data for the RI algorithm have been published in the team’s previous article (Walsh et al., 2008).Reference Walsh, Ramsay, Lapinlampi, Särkelä, Viertiö-Oja and Meriläinen 11 The inclusion criteria for the patients were (1) admission to the ICU because of a critical illness, (2) requirement for mechanical ventilation, and (3) written consent from the patient or relatives. The exclusion criteria were (1) drug overdose, (2) need for neuromuscular paralysis at the time of enrollment, (3) age <16 years, (4) pregnancy, and (5) status epilepticus.
After local ethics committee approval and written informed consent, each subject was monitored for a maximum of 72 hours. Monitoring was stopped earlier if mechanical ventilation was discontinued, the patient or a relative requested discontinuation, the patient died, or technical reasons made continuation impossible. Two Entropy sensors were placed on the patient’s forehead as described in Figure 2. As a reference measure for the depth of sedation, we used a self-developed modified Ramsay scoring system that used standardized vocal, loud vocal, and tetanic stimuli. The vocal stimulus was a verbal command “open your eyes” spoken to the patient. The loud vocal stimulus was a similar verbal command with louder intensity. The tetanic stimulus was a repeated electrical stimulus (50 Hz, 40 mA, 250 ms) applied to the ulnar nerve. The purpose of this modification was (1) to diminish the variation in the assessment by using machine-generated stimuli and (2) to enable discrimination of deep sedation levels in which the patient is unconscious by applying tetanic stimulus. The assessment started with the lightest stimulus and proceeded step by step to stronger stimuli. If the patient at any step showed an observable response, no further stimuli were applied. The correspondence of the standard stimuli, the modified Ramsay scores, and the expected responses is shown in Table 1. A single trained research nurse performed all the modified Ramsay assessments, minimizing variation between the assessments. The interval between the assessments was at least 30 minutes and depended on the clinical duties of the research nurse. The RI values were not shown to the clinical staff during data collection.
Flow chart of the responsiveness index algorithm. (a, b) The impulse response of the filter used to detect the response patterns from the FEMG data. (c) The weighting function g and (d) the scaling function S used in the RI algorithm.
Definition of the modified Ramsay scale
N.A.= not applicable.
Based on the clinician’s review of the patient record and patient state after regaining consciousness, a retrospective classification into two subgroups was made: group 1 had a low probability of encephalopathy (17 patients) and group 2 (13 patients) had a high probability of encephalopathy. This exploratory analysis attempted to identify patients with significant encephalopathy and/or delirium caused by critical illness, which would be a confounder to normal responses to vocal and possibly physical stimuli.
Algorithm Description
The RI algorithm (Figure 2) calculates the root mean square FEMG power in 5-second epochs using the frequency band between 50 and 150 Hz. The power calculation is performed in a time domain after applying a high-pass finite impulse response filter and a 10-Hz comb filter to the signal. The comb filter is used to attenuate the mains frequencies (50 and 60 Hz and multiples). The power value time series is then subjected to a filter that extracts the steep rises related to patient responses from the signal. The filter is implemented as a finite impulse response filter with a step-shaped impulse response depicted in Figure 2b. The shape of the filter was chosen based on the results of the median feature-aligned FEMG power response analysis. For constant FEMG activity, the filter output is zero regardless of the baseline level. The length of the filter is eight samples corresponding to 40 seconds in the FEMG power time series. If the filtered power value at time point t is denoted by P F (t), the RI value at time t is written as
$$RI(t)=S\left( {\mathop{\sum}\limits_{n=t{\minus}\Delta }^t {g(n)log\left( {{{P_{F} (n)} \over {1\mu V}}{\plus}1} \right)} } \right)$$
where S is a scaling function, g is a weighting function, and Δ is the number of P F values included in the summation. Hence, the RI is a weighted average of the logarithm of the P F values adding 1 µV in the given time window. The addition of 1 µV was introduced to set the minimum value of the logarithm to zero. The purpose of the weighting function g (Figure 2c) is to operate as a low-pass filter and to give more weight to recent FEMG changes compared with those occurring in the past. The value of Δ was set at 720 samples (corresponding to 1 hour). The weighting function g gives a 50% contribution to the latest approximate 17 minutes of the 1-hour time window in the RI calculation.
For enhanced usability, the RI value is scaled as a 0-100 range by a scaling function S shown in Figure 2d. The introduction of the scaling was driven by the clinical need to maximize discrimination within deeper sedation states where clinical assessments based on intermittent assessments lack discrimination. The scaling is linear in the range 0-50, which is the relevant range for detecting deeper sedation states, whereas values >50 are scaled nonlinearly. The maximum RI value is set at 100.
The RI algorithm is tolerant to typical EEG artifacts such as patient eye movements, blinks, and movements because these phenomena are present in the low end of the frequency spectrum and are discarded in the filtering steps of the algorithm. Short-lasting (<1 second) EMG bursts in this context are regarded as artifacts and also filtered out. The frequency band used in the calculation ranges from 50 to 150 Hz with removal of the main frequencies. FEMG, when measured with surface electrodes with relatively large area, has a wide frequency spectrum, and the choice of the frequency band is not critical as long as a wide enough band is used.
ROC Analysis and Cutoff Values
RI values were calculated from all the patients in the development data set. The RI values preceding each modified Ramsay level assessment by 1 minute were extracted from the data. The 1-minute advance was used to avoid the effect of the Ramsay assessment procedure on the RI value. The sensitivity and specificity and the related RI threshold value for detecting deeper sedation states (defined as modified Ramsay levels 5 and 6) were determined by receiver operating characteristic (ROC) analysis. Sensitivities and specificities were calculated by altering the threshold RI value. The cost of false negatives and false positives was presumed to be equal, so the optimal threshold was selected by maximizing the sum of sensitivity and specificity. This corresponds to choosing a point on the ROC curve maximizing the distance from the diagonal.
The total number of Ramsay scores, number of rejected data points (with reasoning), and distribution of data points to different Ramsay levels are shown in Table 2.
Number of observations across modified Ramsay scores in the development data
RI= Responsiveness Index.
Subgroup 1/2= low/high probability of encephalopathy.
Results
FEMG Response to Vocal Stimulus
The median and 25th and 75th quartile FEMG powers from the primary Entropy sensor are shown in Figure 3 for patients not responding to the vocal stimulus (a) and for the patients who did respond (b) to the stimulus. In the responding group, the poststimulus FEMG power was significantly higher than the prestimulus power, with a median difference of 0.33 µV (quartiles: −0.00, 1.38, p<0.001) between 30 seconds median power before and after the stimulus. In the nonresponding group, the pre- and poststimulus powers did not differ (median difference 0.00 µV, quartiles: −0.04, 0.05, p=0.73). For the secondary Entropy sensor, the results were similar. In the responding group, poststimulus FEMG power was significantly higher than the prestimulus power with a median difference of 0.23 µV (quartiles: −0.07, 0.86, p<0.001). In the nonresponding group, the pre- and poststimulus powers did not differ (median difference 0.00 µV, quartiles: −0.05, 0.05, p=0.96).
Median FEMG power (25th and 75th quartiles) related to (a) the patients not responding to the vocal stimulus and (b) the patients responding. (c) Feature-aligned power. Time t = 0 is the time of the stimulus onset in (a, b). In (c), the time t=0 indicates the response onset.
As shown in Figure 3c, the feature-aligned FEMG power shows a very sharp rise followed by an elevated baseline level after the response. The initial peak in the aligned response curve may be enhanced by the way the alignment was performed based on the power data itself and by transient muscle activity related to the eye-opening command. However, a fast rise followed by a sustained power increase was a typical finding also in the unaveraged EMG power; therefore, the stepwise increase pattern was used for the RI algorithm.
RI Algorithm Performance
The data set included 30 critically ill patients’ general ICU data from Edinburgh Royal Infirmary (22 males, 8 females) with a mean age (SD) of 59.0 (17.5) years. The algorithm performance was measured within the development data set with the purpose of guiding the iterative algorithm development. The reported results are calculated with the final algorithm settings described previously in the Algorithm Description section.
The RI value distributions across the modified Ramsay levels in all patients and subgroups 1 and 2 are shown in Figure 4. The optimal threshold values for detecting deep sedation in all patients and subgroups 1 and 2 are presented in Table 3. The corresponding ROC curves are shown in Figure 5. Figure 6 shows an example of a patient case in the development data set (subgroup 1) with EEG, FEMG power, filtered FEMG power, RI, and the modified Ramsay score.
Boxes and whiskers of RI values across different modified Ramsay levels. The white boxes indicate the primary Entropy sensor and the gray boxes the secondary Entropy sensor. The horizontal lines within the boxes denote median RI values and the box edges the 25th and 75th quartiles. Statistical outliers denoted by a + sign are defined as data points whose distance from the 25th or 75th quartile lines was greater than 1.5 times the interquartile range.
The ROC curves for the primary (black curves) and secondary (gray curves). Entropy sensors in all patients and the subgroups with a low and high probability of encephalopathy. The optimal RI threshold values are shown as filled circles.
Algorithm example from actual patient case in the development data set in sub-group 1. From top to bottom are shown (a) the EEG waveform, (b) FEMG power, (c) the Responsiveness Index, and (d) the modified Ramsay score.
The optimal Responsiveness Index threshold values and related sensitivities and specificities related to detecting deep sedation (Richmond Agitation-Sedation Scale >5) with the RI algorithm
Subgroup 1/2= low/high probability of encephalopathy.
Discussion
The RI algorithm has the potential to be a robust and objective method for quantifying the balance between the sedation state and the stimulation level of an ICU patient. It is not a direct measure of the depth of sedation and should not be treated as such. For example, a patient receiving little or no sedation and sleeping calmly may respond strongly to stimulation-related care procedures or noise. The same patient will show low RI values if he or she is sleeping calmly and no stimulation is applied.
There are some limitations related to this study. First, the patient group used in the first phase (characterization of ICU patients’ FEMG activity) was limited to a specific condition and does not represent general ICU patients. Furthermore, hepatic encephalopathy may cause neural alterations that in turn might affect patients’ responses to stimuli. Second, the development data set used in the RI algorithm development was quite small, containing only 30 patients. However, the algorithm is relatively simple and robust with few adjustable parameters. Good results were previously reported in another patient population,Reference Walsh, Lapinlampi, Ramsay, Särkelä, Uutela and Viertiö-Oja 15 which provides confidence in the method and warrants further investigations in larger patient populations.
Some of the potentially encephalopathic patients showed increased levels of EMG responses without purposeful movements as reaction to the auditory stimulus (Figure 3, modified Ramsay scores 5 and 6). Our hypothesis is that encephalopathy may alter the patients’ reactions to stimuli. In particular, responses to stimuli expecting a cognitive reaction (such as a command to open one’s eyes) might be compromised. This would make traditional sedation assessments (including the modified Ramsay score used in this study) unreliable in these patients because responses might be suppressed from the encephalopathy instead of the sedation state.
In the analysis, we defined Ramsay levels 5 and 6 as deeper sedation. The definition of deep sedation, however, depends on the clinical requirements of individual patients. Recent studies recommend a target Ramsay score of 3Reference Brook, Ahrens and Schaiff 30 or from 3 to 44 in general. Guidelines, however, recommend that the need for sedation should be assessed according to patient and situation.Reference Jacobi, Fraser and Coursin 31 The main purpose of the RI algorithm is to alert the hospital staff to deep sedation. The discrimination at lighter sedation levels was considered less important in the development phase because light sedation is usually easily observable from patient behavior. The time scale of the RI algorithm is rather long (60 minutes on average) and therefore it is not applicable in cases where the time scale of events is short, such as procedural sedation.
Measuring the performance of the RI algorithm by comparing it with clinical sedation scales such as Ramsay and the Richmond Agitation-Sedation Scale may understate the actual performance of the algorithm. These clinical scales are the current gold standard, but they measure a slightly different phenomenon. Our hypothesis is that the RI algorithm measures the balance between the sedative drug effect and the level of stimulation (i.e., if the sedation level remains constant, the RI value may change when the stimulation level changes). The clinical scales, on the other hand, measure patients’ response to certain stimuli as perceived by a caring nurse.
Using neuromuscular blocking agents (NMBA) is believed to prevent the use of a developed monitoring method because it is based on changes in FEMG activity. However, upper facial muscles are more resistant to the effect of NMBAs than the hand muscles, where the effect of the neuromuscular block is usually measured.Reference Paloheimo, Wilson, Edmonds, Lucas and Triantafillou 32 Today the use of NMBAs in the ICU is restricted to limited patient groups and are not recommended in the current guidelines.Reference Nasraway, Jacobi, Murray and Lumb 5
This study did not address the patients’ related outcomes, which could be looked at in a bigger trial.
Acknowledgements and Funding
The authors thank research nurses Pam Ramsay, Petra Peltola, and Kristiina Järvelä for their work in collecting the data in the Royal Infirmary of Edinburgh and in Helsinki University Hospital. We wish to show our appreciation for the late Dr. Markku Paloheimo’s guidance and pioneering research on the application of EMG in anesthesia and critical care.
T.P. Lapinlampi, H.E. Viertiö-Oja, K.H. Uutela, and M.O.K. Särkelä are employees of GE Healthcare Finland Oy, as was M. Helin was at the time of the study. A. Vakkuri, G.B. Young, and T.S. Walsh are or have been consultants for GE Healthcare. GE Healthcare owns patent rights related to the RI measurement. The study was designed in cooperation with GE Healthcare and Edinburgh Royal Infirmary (T.S.W.).
The study was funded by GE Healthcare and the SalWe Research Program for Mind and Body (Tekes - the Finnish Funding Agency for Technology and Innovation grant 1104/10).
Statement of Authorship
TL carried out the statistical analysis for RI algorithm performance, participated in the RI algorithm development, and was the main writer of the manuscript. HV participated in the design of the study and RI algorithm development. MH carried out the statistical analysis for FEMG responses to vocal stimuli. KU participated in the design of the study and RI algorithm development. MS participated in the RI algorithm development. AV participated in the design of the study and data collection from the hepatic encephalopathy patients. BY wrote the Physiological Background section. TW participated in the design of the study and data collection from the general ICU patients. All authors read and approved the final manuscript.
