| ||||||||||||||
|
|
|||||||||||||

*Department of Anesthesiology and Intensive Care Medicine, University of Saarland, Homburg/Saar, Germany, and the
Department of Anesthesiology and Intensive Care Medicine, University of Bonn, Bonn, Germany
Address correspondence and reprint requests to Wolfram Wilhelm, MD, DEAA, Department of Anesthesiology and Intensive Care Medicine, University of Saarland, 66421 Homburg/Saar, Germany. Address email to wolfram.wilhelm{at}uniklinik-saarland.de
| Abstract |
|---|
|
|
|---|
IMPLICATIONS: This study in 18 adult patients undergoing radical prostatectomy describes the relationship between Narcotrend index and bispectral index versus predicted propofol effect compartment concentrations. In terms of prediction probability, the performance of the Narcotrend index and the bispectral index to predict propofol effect site concentrations was comparable.
| Introduction |
|---|
|
|
|---|
In contrast to the well-established and intensively evaluated BIS, published data regarding the Narcotrend monitor are still sparse. Although first studies indicate the favorable use of the Narcotrend to titrate anesthesia in clinical practice (5,6), we studied the performance of Narcotrend index compared to BIS as a measure of anesthetic drug effect during propofol anesthesia. It was the hypothesis of this investigation that Narcotrend and BIS are equivalent in detecting differences in EEG dynamics.
| Methods |
|---|
|
|
|---|
6 k
for the Narcotrend and
7.5 k
for the BIS as required by the manufacturers. Induction of anesthesia was started with a fentanyl bolus of 0.1 mg; 5 min later 36 mg · kg-1 · h-1 propofol was given for hypnosis using an infusion system (Graseby 3400; Graseby Medical Limited, Watford, UK). After loss of consciousness oxygen was given by face mask ventilation, the propofol infusion was reduced to 16 mg · kg-1 · h-1, patients received 0.5 mg/kg of atracurium, the trachea was intubated 3 min later, and the lungs were ventilated to an end-tidal carbon dioxide concentration of 35 mm Hg. Immediately after intubation, patients received 15 mL bupivacaine 0.5% epidurally, whereas the propofol infusion was adjusted according to clinical needs.
To ensure conditions of constant surgical stimulation, only the EEG data obtained during prostata preparation were used for further analysis. In all cases, a waiting period of at least 45 min after induction of anesthesia was maintained to minimize the effects of the fentanyl administered for induction. For study measurements the propofol infusion rate was increased to 26 mg · kg-1 · h-1 until a substantial burst suppression pattern was recognized in the raw EEG. Subsequently, the propofol infusion rate was reduced to 0.5 mg · kg-1 · h-1 until a BIS value of 50 was reached. Two of these cycles were performed in each patient; complete neuromuscular blockade during the investigation period was ensured by repetitive injections of 0.25 mg/kg atracurium.
Study data were automatically recorded in intervals of 5 s during the 2 concentration curves. BIS values were recorded using the software program Hyperterminal (Microsoft, Redmond, WA). The Narcotrend index was recorded using a modified research laptop version of the Narcotrend software 4.0. The propofol infusion pump (Graseby 3400) was connected via an RS-232 interface to a computer running the software Rugloop (simulation version, Department of Anesthesia, Ghent University, Ghent, Belgium) which automatically calculated and recorded propofol plasma and effect site concentrations using the Diprifusor data set as published by Marsh et al. (7). The synchronization of the data was performed automatically by the software Access 2000 (Microsoft).
The correlation between Narcotrend or BIS and the predicted propofol effect compartment concentration was calculated using prediction probability (PK). Smith et al. (8) elucidated the prediction probability as follows: Given two randomly selected data points with distinct anesthetic drug concentration, the PK value describes the probability that the EEG correctly predicts which of the data points is the one with the larger (or smaller) anesthetic drug concentration. As a nonparametric measure, PK is independent of scale units and does not require knowledge of underlying distributions or efforts to linearize or to otherwise transform scales. Furthermore, PK can be computed for any degree of coarseness or fineness of the scales. Thus, PK fully uses the available data without imposing additional arbitrary constraints.
The relationship between ordinal variables x (indicator value, i.e., in our study BIS value and Narcotrend index value, respectively) and y (i.e., in our study propofol effect-site concentration) is described in terms of the rank ordering of the x and y values for pairs of data points. A concordance occurs in our study when the x values and y values for a pair of data points are rank ordered in the way that the data point with the larger propofol effect site concentration exhibits the lower BIS or Narcotrend index value. A discordance occurs when the x and y values are rank ordered in the opposite direction, i.e., when the data point with the larger propofol effect-site concentration exhibits the higher BIS or Narcotrend index value. An indicator-only, or x-only, tie (tie in indicator value but not in propofol effect-site concentration) occurs when the x values are tied but the y values are not (i.e., the same bispectral or Narcotrend index value but different propofol effect-site concentration). Similarly, a y-only tie occurs when the x values are not tied but the y values are (i.e., different BIS or Narcotrend index value but the same propofol effect-site concentration). Finally, a joint tie occurs when there are ties in both x and in y (i.e., the same BIS or Narcotrend index value and the same propofol effect site concentration). Concordances are desired; discordances and indicator-only ties are undesirable. Ties in y (both y-only and joint ties) will not be considered according to the PK algorithm.
PK has been defined (8) as PK = (Pc + 0.5*Ptx)/(Pc + Pd + Ptx), where Pc, Pd, and Ptx are the respective probabilities that two data points drawn at random, independently and with replacement, from the population are a concordance, a discordance, or an x-only tie.
When comparing indicators, however, it is necessary to gather data using the same procedure and over the same distribution of e.g., propofol concentrations; therefore it is recommended to measure the EEG index values simultaneously for the same subjects.
Data Analysis and Statistics
First, we calculated the prediction probability for Narcotrend index and BIS using the propofol effect compartment concentrations as calculated by the Diprifusor data set with a fixed ke0 value which is integrated in the Diprifusor algorithm.
In a second separate step we estimated an individual ke0 value for each patient using the fractional sigmoid Emax model for the relationship between effect compartment concentration and the respective EEG parameter (Hill equation) (9) with input from a first-order effect-site model: equation
|
|
and E = E0 * (1-ceff
/(C50
+ ceff
))
Ceff is the apparent effect site concentration, Cplasma is the predicted propofol plasma concentration, keo is the first order rate constant determining the efflux from the effect compartment, E0 is the measured baseline effect of each individual, C50 is the concentration that causes 50% of the maximum effect, and
describes the slope of the concentration-response relation.
Data from each person were fitted separately. The parameters of the above models were estimated using nonlinear regression with ordinary least squares. The computations were performed on a spreadsheet using the Excel 2000 software program (Microsoft), and the parameters were optimized with the Solver tool within Excel.
We now calculated the prediction probability for Narcotrend index and BIS using the propofol effect compartment concentrations calculated from the predicted plasma concentrations (as obtained from the Diprifusor data set) and the individually optimized ke0 values. Differences in PK values between Narcotrend index and BIS were tested with Wilcoxons test. Statistical significance was assumed at P < 0.05.
Data are presented as mean and SD. Statistical calculations were performed using SigmaStat 2.03 and SigmaPlot 2000 computer software (SPSS GmbH, Erkrath, Germany).
| Results |
|---|
|
|
|---|
The predicted propofol effect site concentrations reached from a smallest concentration of 2.0 ± 0.4 µg/mL to a largest concentration of 6.3 ± 1.3 µg/mL. With increasing propofol plasma concentrations, Narcotrend and BIS decreased (and vice versa) with a time delay (Fig. 1). Plotting the Narcotrend index versus propofol plasma concentration or the BIS versus propofol plasma concentration revealed hysteresis, which collapsed by introduction of the propofol effect-site concentration (Fig. 2). The performance of the Narcotrend index (PK = 0.88 ± 0.03) as an index of propofol effect-site concentrations, as predicted by the Diprifusor data set, was not significantly different from that of the BIS (PK = 0.85 ± 0.04). Figure 3 shows the individual PK values for each patient.
|
|
|
The mean ke0 individually optimized while using Narcotrend index as a measure of drug effect was 0.20 ± 0.05 min-1, whereas it was 0.16 ± 0.07 min-1 when the BIS was used. The estimates of ke0 for Narcotrend and BIS are summarized in Figure 4.
|
| Discussion |
|---|
|
|
|---|
For the data analysis in this study we preferred the PK value as the objective function instead of the R2 value that had been used in previous studies investigating EEG parameters as measures of anesthetic drug effect (1012). Calculating R2 necessitates the use of a specific pharmacodynamic model. A fractional sigmoid Emax model based on the Hill equation is usually used to that purpose. However, in the present study a burst suppression pattern was observed in all patients. The appearance of burst suppression leads to a switch in the BIS algorithm with a possible model misspecification while using a single fractional sigmoid Emax model (13). Therefore, we decided to use the model-independent calculation of the prediction probability PK. Given two randomly selected data points with distinct anesthetic drug concentrations, the prediction probability PK is the probability that the EEG parameter correctly predicts which of the data points is the one with the larger (or smaller) anesthetic drug concentration. This method was previously shown to be an adequate alternative to the R2 value for estimating the performance of an EEG parameter as a measure of anesthetic drug effect (14) and to predict propofol target concentration (15).
Nonetheless, it should be noted that we used the fractional sigmoid Emax model for estimating individual values for ke0 because no model-independent method is available to estimate ke0 for increasing and decreasing anesthetic drug concentrations with a nonzero starting point.
We determined that a ke0 of 0.25 min-1 is incorporated in the Diprifusor data set. This is the same ke0 value that is used in the computer program STANPUMP, based on preliminary work of Dyck et al.,1 and that has also been used experimentally (16). This ke0 value is very close to the ke0 of 0.239 min-1 of Schüttler et al.,2 which has been used by Vuyk et al. (17) and others for their simulations. Our individually optimized ke0 values are lower. This is in accordance with a mean ke0 value of 0.20 min-1 for propofol as reported by Billard et al. (18) with an older BIS version. Interestingly, when investigating volatile anesthetics Olofsen and Dahan (12) found lower ke0 values for BIS than for spectral edge frequency. We also found lower ke0 values for BIS than for Narcotrend index. This may be explained by different calculation intervals; i.e., 30 s for BIS XP (19) and 20 s for the Narcotrend (6). With very short calculation intervals of 2 s Schnider et al. (20) reported even larger ke0 values of 0.456 min-1 for propofol. Summarizing these results, it appears that the ke0 value predominantly depends on the length of the calculation interval, whereas other explanations, e.g., a different site of action for BIS and for the other EEG parameters, are less probable.
We did not measure propofol plasma concentrations but relied on the predicted propofol concentrations. The prediction of propofol concentrations by the Diprifusor is based on a three-compartment model and the calculated pharmacokinetic variable set by Marsh et al. (7). In general, pharmacokinetic variable sets are determined from collectives of patients or subjects. However, the pharmacokinetic variable set of a single individual patient can certainly be different, resulting in a prediction error, i.e., the difference between predicted and measured concentrations, of up to 30% (21). Nevertheless, the prediction of propofol concentrations based on pharmacokinetic principles still makes sense. Predicted and measured concentrations are almost always parallel, and the percentage of deviation remains relatively constant during the complete duration of anesthesia (22). Thus, even though the absolute value of the predicted concentrations may be incorrect, the duration and changes undertaken are in relation correctly reflected. Eventually, this is sufficient for our study design, as only the performance of the EEG parameter to reflect the dynamic (the relative and not the absolute) changes in the single patient was the subject of the study. The excellent PK values obtained in this investigation further support the adequacy of our approach.
Although the algorithms of Narcotrend index (6) and BIS (23) are completely different, both have a common composition based on an optimization of subparameters using extensive databases. The first step in the development of the automatic classification algorithm of the Narcotrend was to build a database of typical examples for all substages from A to F. For this purpose, artifact-free EEG epochs from a large database of one-channel EEG recordings collected from the electrode positions C3-P3 during anesthesia with thiopental-enflurane or propofol were selected. Initially, more than 1000 artifact-free EEG epochs with a length of 20 s were visually classified to form the basis for the development of automatic classification algorithms. From these EEG epochs numerous quantitative EEG features from the time and the frequency domain were extracted, e.g., spectral parameters, entropy measures, and autoregressive parameters (24). The extracted EEG parameters were statistically analyzed to identify a subset of EEG parameters that were best suited to discriminate between the different visually determined EEG substages. Age-related changes in the EEG were incorporated. The resulting parameters were combined in multivariate discriminant functions to classify an EEG epoch into one of the EEG substages between A and E. The discriminant analysis yields probabilities for the degree of similarity of an EEG epoch with the typical EEG stages A to E during anesthesia. In addition, algorithms for the classification of stage F were developed based on the proportion and intensity of very flat EEG segments. The algorithms were revised continuously including extensions for the standard use of frontal electrode positions. Obviously these efforts were successful. First studies indicate the favorable use of the Narcotrend to titrate anesthesia in clinical practice with reduction of drug consumption and shortening of recovery times (5,6). In addition, studies using an older version of the Narcotrend software, before the introduction of the Narcotrend index from 0 to 100, showed dose-dependent changes of the Narcotrend stages A to F, during emergence from propofol anesthesia (15) and during emergence from desflurane anesthesia (25). The Narcotrend index appeared to read larger at small propofol concentrations and smaller at larger propofol concentrations than the BIS. This phenomenon underlines that the parameterization of the raw EEG by Narcotrend index and BIS is different and that a simple 1:1 translation of the index values is obviously inadequate. Further investigations are needed to clarify the comparability of Narcotrend index and BIS values and to guide the clinician in interpreting various Narcotrend index values.
In summary, we were able to demonstrate a close quantitative correlation between Narcotrend index and predicted propofol effect compartment concentrations. In this study, the Narcotrend index detected differences in EEG dynamics in the observed concentration range as well as BIS. However, the small size of the study group and its narrow composition, as well as the particular design of the anesthetic, might limit the ability of these data to be extrapolated to other or larger populations.
| Acknowledgments |
|---|
| Footnotes |
|---|
2 Schüttler J, Schwilden H, Stoeckel H. Pharmakokinetic-dynamic modeling of Diprivan [abstract]. Anesthesiology 1986;65:A-549. ![]()
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
S. Kreuer, J. Bruhn, E. Walter, R. Larsen, C. C. Apfel, U. Grundmann, A. Biedler, and W. Wilhelm Comparative Pharmacodynamic Modeling Using Bispectral and Narcotrend-Index With and Without a Pharmacodynamic Plateau During Sevoflurane Anesthesia Anesth. Analg., April 1, 2008; 106(4): 1171 - 1181. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. N. Sessler and K. Varney Patient-Focused Sedation and Analgesia in the ICU Chest, February 1, 2008; 133(2): 552 - 565. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. K. Ellerkmann, M. Soehle, T. M. Alves, V.-M. Liermann, I. Wenningmann, H. Roepcke, S. Kreuer, A. Hoeft, and J. Bruhn Spectral entropy and bispectral index as measures of the electroencephalographic effects of propofol. Anesth. Analg., May 1, 2006; 102(5): 1456 - 1462. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. F. Russell The Narcotrend 'depth of anaesthesia' monitor cannot reliably detect consciousness during general anaesthesia: an investigation using the isolated forearm technique Br. J. Anaesth., March 1, 2006; 96(3): 346 - 352. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Kreuer, J. Bruhn, C. Stracke, L. Aniset, M. Silomon, R. Larsen, and W. Wilhelm Narcotrend or Bispectral Index Monitoring During Desflurane-Remifentanil Anesthesia: A Comparison with a Standard Practice Protocol Anesth. Analg., August 1, 2005; 101(2): 427 - 434. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|