Skip Navigation LinksHome > September 2007 - Volume 107 - Issue 3 > Influence of Administration Rate on Propofol Plasma–Effect S...
Anesthesiology:
doi: 10.1097/01.anes.0000278902.15505.f8
Clinical Investigations

Influence of Administration Rate on Propofol Plasma–Effect Site Equilibration

Struys, Michel M. R. F. M.D., Ph.D.*; Coppens, Marc J. M.D.†; De Neve, Nikolaas M.D.‡; Mortier, Eric P. M.D., D.Sc.§; Doufas, Anthony G. M.D., Ph.D.∥; Van Bocxlaer, Jan F. P. Pha., Ph.D.#; Shafer, Steven L. M.D.**

Free Access
Supplemental Author Material
Article Outline
Collapse Box

Author Information

Collapse Box

Abstract

Background: The authors hypothesized a difference in plasma–effect site equilibration, depicted by a first-order constant ke0, depending on the injection rate of propofol.
Methods: Sixty-one patients received 2.5 mg/kg propofol given as a bolus or as a 1-, 2-, or 3-min infusion. The Bispectral Index was used to monitor drug effect. Propofol predicted plasma concentration was calculated using a three-compartment model and the effect site concentration over time as the convolution between the predicted plasma concentration and the disposition function of the effect site concentration. The authors evaluated the influence of the infusion rate on the ke0 by comparing the model with one ke0 for all groups with models estimating different ke0 values for each group. The authors also assessed the accuracy of two pharmacokinetic models after bolus injection.
Results: The best model based was a fixed (Bispectral Index ≥ 90) plus sigmoidal model (Bispectral Index < 90) with two values of ke0, one for the bolus (t½ ke0 = 1.2 min) and one for the infusions (t½ ke0 = 2.2 min). However, the tested pharmacokinetic models poorly predicted the arterial concentrations in the first minutes after bolus injection. Simulations showed the requirement for two ke0 values for bolus and infusion was mostly a compensation for the inaccurate prediction of arterial concentrations after a bolus.
Conclusion: Propofol plasma–effect site equilibration occurs more rapidly after a bolus than after rapid infusion, based on the electroencephalogram as a drug effect measure, mostly because of misspecification of the pharmacokinetic model in the first minutes after bolus.
Back to Top | Article Outline

ArticlePlus

Click on the links below to access all the ArticlePlus for this article.
Please note that ArticlePlus files may launch a viewer application outside of your web browser.
* http://links.lww.com/ALN/A484
* http://links.lww.com/ALN/A485
* http://links.lww.com/ALN/A486
* http://links.lww.com/ALN/A487
* http://links.lww.com/ALN/A488
PROPOFOL transfer between the plasma and effect site can be modeled as a first-order process characterized by ke0.1,2 The standard model of ke0 assumes that the rate of equilibration between the plasma and the site of drug effect is independent of the rate of drug administration. However, there are conflicting data on the rate of equilibration between the plasma and the site of propofol drug effect. In a study involving both bolus injections and intravenous infusions, Schnider et al.3 found that the rate of equilibration was rapid, with a half-time of equilibration, t½ ke0, of 1.5 min, and a peak effect, tpeak, of 1.7 min. Schnider's finding of rapid equilibration was subsequently validated by Struys et al.4 However, using continuous infusions of propofol, Doufas et al.5 found a much slower rate of plasma–effect site equilibration, with a t½ ke0 of 4.1 min, and a tpeak of 2.7 min. They also found that infusion rate had no influence on ke0.5
The maximum propofol infusion rate in the study of Doufas et al.5 was 60 mg/min, far lower than the maximum rate of approximately 500 mg/min required for Schnider et al. to give a 2.5-mg bolus over 20 s. Doufas et al.5 proposed that there could be a fundamental difference in plasma–effect site equilibration depending on whether propofol was given as a bolus or continuous infusion. We investigated this hypothesis.
Back to Top | Article Outline

Materials and Methods

Clinical Protocol
After institutional ethics committee approval, written informed consent was obtained from 61 female patients with American Society of Anesthesiologists physical status I, aged 18–45 yr, scheduled to undergo ambulatory gynecologic surgery. Exclusion criteria included weight less than 70% or more than 130% of ideal body weight, neurologic disorder, and recent use of psychoactive medication, including alcohol.
All patients were randomly assigned to one of four groups to receive 2.5 mg/kg propofol (Diprivan 1%; AstraZeneca, London, United Kingdom) given as a bolus (within 10 s) (group 1) or given as a continuous infusion over 1, 2, or 3 min (groups 2, 3, and 4, respectively). Bolus administration was performed manually. The continuous infusions were administered using a Fresenius Modular DPS Infusion Pump connected to a Fresenius Base A (Fresenius Vial Infusion Systems, Brézins, France). To ensure synchronized data recording, all monitor and infusion data were continuously captured by a computer running RUGLOOP II (Demed, Temse, Belgium) via multiple RS 232 interfaces. By tracking the infused propofol volume continuously in groups 2, 3, and 4, RUGLOOP II calculated the corresponding plasma concentration using the three-compartment model previously published by Schnider et al.6 This model was selected because of its optimal performance in previous studies.5 For group 1, the plasma concentration was calculated post hoc using RUGLOOP II simulation mode and the Schnider propofol pharmacokinetic model.6
Propofol was infused via a large left forearm vein. Every patient received approximately 100 ml crystalloid fluid during the study period. No fluid load was given before induction. No patient received preanesthetic medication. No other drugs were given. All patients maintained spontaneous ventilation via a facemask delivering 100% O2. Before starting the drug administration, all patients were asked to close their eyes and relax for 2 min. Thereafter, baseline measures were taken. The operating room was kept silent to avoid noise-related stimulation and artifact.
Propofol drug effect was continuously monitored using the Bispectral Index (BIS; version 4.0; Aspect Medical Systems, Inc., Newton, MA). The BIS was derived from the frontal electroencephalogram and calculated by the A-2000 BIS® monitor using the four BIS®-Sensor electrodes. Electrode impedance was less than 5 kΩ. The smoothening time of the BIS® monitor was set at 15 s. The BIS data were captured in real time on a laptop computer using RUGLOOP II. Heart rate, noninvasive blood pressure, oxygen saturation measured by pulse oximetry, and capnography were recorded at 1-min time intervals using an S-5 monitor (Datex-Ohmeda, Helsinki, Finland) and were also captured electronically using RUGLOOP II. Averaging of the data was performed using 10-s intervals. All patients were monitored until return of consciousness after propofol administration, defined as spontaneous eye opening (without a stimulus).
Back to Top | Article Outline
Pharmacodynamic Modeling and Estimation of ke0
Equation (Uncited)
Equation (Uncited)
Image Tools
In our initial approach, the effect site was assumed to be linked to the plasma by a compartment of trivial volume with a first-order equilibrium constant of ke0. The relation between propofol effect site concentration (Ce) and the electroencephalographic measures of anesthetic drug effect was described using a classic sigmoid Emax model:
Equation (Uncited)
Equation (Uncited)
Image Tools
where Effect is the electroencephalographic effect (e.g., the measured BIS), E0 is the baseline measurement when no drug is present, Emax is the maximum possible drug effect, Ce is the calculated effect site concentration of propofol, Ce50 is the Ce associated with 50% maximal drug effect, and γ is the steepness of the concentration-versus-response relation. The model parameters were estimated using NONMEM V (GloboMax LLC, Hanover, MD). For Ce50 and ke0, interindividual variability was permitted using a log-normal distribution:
where Pi is the parameter value in the ith patient, PTV is the typical value of the parameter in the population, and η is a random variable with a mean of 0 and a variance of ω2. Individual variability is reported as ω, the SD of η in the log domain, which is approximately the coefficient of variation in the standard domain. Residual intraindividual variability was modeled using a standard additive error model.
After visual inspection of the BIS data above 90, and initial attempts to find a value of ke0 that fit all of the observations to a single pharmacodynamic model, we concluded that no model could be fit to the data. The raw BIS data showed an abrupt (< 15 s) transition from a BIS greater than 90 to a very low BIS, suggesting a nearly instantaneous state change. Therefore, we investigated the relation between propofol Ce and BIS using separate models for periods before and after the state change. The combined model described all BIS data above 90 by a single value model (“one size fits all”) and all data of 90 or less using the classic sigmoid Emax model.
For both model approaches (classic sigmoidal model or the combined fixed plus sigmoidal model), we assumed that the underlying sigmoidal model which describes the relation between the effect site concentration and the drug effect is not influenced by the method of drug administration. Indeed, this is a fundamental assumption underlying the standard model of the effect site. Therefore, we concurrently estimated the model parameters for all four groups, only permitting the value of ke0 to differ between groups. We evaluated the influence of the administration rate on the ke0 by comparing the log likelihood between a model with one ke0 for all administration rates with models estimating different ke0 values for each group. The addition of infusion rate specific values of ke0 was considered statistically significant when the log likelihood decreased by at least 6.63 (P < 0.01, chi-square test with 1 degree of freedom).7
NONMEM had difficulty simultaneously estimating the parameters of the structural model (e.g., the θ parameters) and the variance model (the ω parameters). Specifically, estimating both the structural and variance models simultaneously produced highly biased estimates of the parameters, as demonstrated by the observation that the post hoc estimates of the η parameter were uniformly either positive or negative. Therefore, we chose the best model using the naive pooled data method as described by Kataria et al.8 After identifying the best structural model, we fixed the parameters of the structural model at the naive pooled data estimates and estimated the parameters of the intersubject and intrasubject variability models. This approach comes very close to the generalized least-squares method. The generalized least-squares approach involves three steps. In the first step, the intersubject variability (i.e., ω) is fixed at zero, and the typical values of the structural model (i.e., θ) are estimated. In the second step, the typical values of the structural model (i.e., θ) are fixed at the estimates from the first step, and the intersubject variability (i.e., ω) is estimated. In the third step, the intersubject variability (i.e., ω) is fixed at the value from step 2, and the typical values of the structural model (i.e., θ) are estimated for the last time. Our approach was to implement the first two steps of the generalized least-squares method.
Equation (Uncited)
Equation (Uncited)
Image Tools
The effect site concentrations over time were calculated as the convolution of the predicted plasma concentrations over time with the disposition function of the effect site, ke0 e−ke0t. The convolution was based on a “connect-the-dots” approach previously used by Schnider et al.3 In brief, increasing plasma concentrations were modeled using a linear interpretation between adjacent plasma concentrations. As such, when concentrations are increasing, the slope in the concentration from time t1 to time t2 can be modeled as
Equation (Uncited)
Equation (Uncited)
Image Tools
where Cp is the plasma concentration, and thus the plasma concentrations from time t1 to time t2 are described by the formula
Equation (Uncited)
Equation (Uncited)
Image Tools
The convolution of this with ke0 e−ke0t calculates the effect site concentration at time t2 as a function of the effect site concentration at time t1 and the plasma concentrations over the interval
Equation (Uncited)
Equation (Uncited)
Image Tools
Similarly, when concentrations are decreasing, the slope of the log of the concentrations from time t1 to time t2 can be modeled as
Equation (Uncited)
Equation (Uncited)
Image Tools
and thus the plasma concentrations from time t1 to time t2 are described by the formula
Equation (Uncited)
Equation (Uncited)
Image Tools
The convolution of this with ke0 e−ke0t calculates the effect site concentration at time t2 as a function of the effect site concentration at time t1 and the plasma concentrations over the interval
The observed BIS value has a time delay for the measurement, which we fixed at 10 s and defined as lag time. This lag time is approximately the sum of half of the smoothing interval within the BIS® monitor and the time for the automated data collection.
Equation (Uncited)
Equation (Uncited)
Image Tools
We assessed the model performance by calculating the prediction error (PE) between measured (BISmeas) and post hoc Bayesian predicted BIS (BISpred) values, as
We also calculated median prediction error (MDPE) and absolute median prediction error, (MDAPE), for each patient.9 Differences among groups for MDPE and MDAPE were tested using a Student t test with Bonferroni correction for multiple testing.
To observe whether the model reflects reality, we compared the observed time of maximum BIS response (tmax,BIS) with the time to reach maximum Ce (tmax,Ce) after the specific administration of propofol in each group.
Back to Top | Article Outline
Bolus Validation Study
The three-compartment pharmacokinetic model published by Schnider et al.6 was used to predict the time course of the propofol plasma concentration in our study. For slowly administered continuous infusion of propofol, Doufas et al.5 found that the pharmacokinetic model published by Schnider et al.6 accurately predicted the propofol plasma concentrations in arterial blood. However, the accuracy of this model in the first 3 min after bolus injection is not described in the literature.
We therefore performed a study to validate the pharmacokinetics after bolus injection. In this study, we collected propofol arterial blood samples from 10 additional patients to evaluate the accuracy of the pharmacokinetic models published by Schnider et al.6 and by Marsh et al.10 in the first 5 min after intravenous propofol bolus. After additional ethics committee approval and written informed consent, 10 patients received a bolus dose of propofol (2.5 mg/kg) within 10 s in a large forearm vein. Propofol arterial blood samples were collected (contralateral from the injection of propofol) at 0, 30, 60, 120, 180, 240, and 300 s after injection. Propofol (bound and free) plasma concentrations were analyzed using a validated gas chromatographic–mass spectrometric method with solid-phase micro extraction. The time course of the propofol plasma concentration as predicted by both pharmacokinetic models (Schnider and Marsh) was simulated using RUGLOOP II. We calculated the PE of the predicted concentration and the MDAPE for both models for the first 5 min after injection.
The arterial concentrations in the first 5 min after bolus injection were compared with the predictions based on the pharmacokinetics reported by Schnider et al.6 and by Marsh et al.10 As described in the Results, neither pharmacokinetic parameter set accurately predicted the concentrations in the first 5 min after bolus injection. To see whether this misspecification might affect the estimation of ke0, we calculated the effect site concentrations over time after bolus injection as a convolution of the median concentrations in these 10 individuals with ke0 e−ke0t, the disposition function of the effect site. We then calculated the value of ke0 that predicted the observed time of peak BIS response.
Back to Top | Article Outline
Validation of the Applied Lag Time in the BIS Measurement
We fixed the lag time of the BIS at 10 s. Based on literature reports,11 a longer lag time might be more appropriate. Therefore, we applied a post hoc analysis on our final selected model in an attempt to explore whether longer lag times would result in a better model for both the sigmoidal model and the fixed BIS plus sigmoidal model. Lag times between 10 and 25 s were evaluated, and the NONMEM objective functions were compared.
Back to Top | Article Outline
Validation of the Modeling Methods Using Another Data Set
We validated our final model results against the data previously published by Doufas et al.5 In brief, Doufas et al. analyzed the pharmacokinetics and pharmacodynamics of propofol in 18 healthy volunteers receiving five consecutive target-controlled propofol infusions. During each infusion, predicted Ce increased linearly at a rate of 0.1, 0.3, 0.5, 0.7, or 0.9 μg · ml−1 · min−1 based on the Schnider pharmacokinetic–pharmacodynamic model. BIS was collected continuously during the infusions. Doufas et al.5 fit a combined pharmacokinetic–dynamic model to their data. To describe the pharmacodynamics of propofol, a classic sigmoid Emax model was used. No lag time in the BIS was included in the model. The study results can be seen in the original publication.5
We validated four models: (1) a sigmoid Emax model without a lag time (same as the model used in the Doufas article; (2) a sigmoidal Emax model with a 10-s lag time (equivalent to our sigmoidal Emax model); (3) a model consisting of a fixed estimate of BIS values of 90 or greater and a sigmoidal Emax model for BIS values less than 90, with no lag time; and (4) a model consisting of a fixed estimate of BIS values of 90 or greater and a sigmoidal Emax model for BIS values less than 90, with a 10-s lag time (equivalent to our best model). We encountered identical problems with concurrent identification of the structural and variance models that we described when fitting our own data, and so we used the naive pooled data approach to identify the best model among the four structural models considered. After identification of the best structural model, we fixed the structural parameters and estimated the interindividual variability.
We assessed the performance of the optimal model by calculating the PE between BISmeas and BISpred values, as well as the MDPE and MDAPE for each patient.9
Back to Top | Article Outline

Results

Model Estimation and Validation
Table 1
Table 1
Image Tools
Table 2
Table 2
Image Tools
All recorded data were used. No patients experienced hemodynamic or respiratory instability during the study. The demographics for patient in the four groups are shown in table 1. The amount of propofol given and the infusion rates are shown in table 2.
Fig. 1
Fig. 1
Image Tools
Figures 1A–D show the relation between the measured BIS and calculated plasma concentrations, Cp, versus time for the four groups, respectively. The dashed line shows the cutoff value of 90 used in the combined fixed plus sigmoidal model approach. Figures 1E–H show the relation between BIS and plasma concentration, revealing the hysteresis in the plasma concentration–versus–response relation.
Table 3
Table 3
Image Tools
Table 3 shows the NONMEM objective function (−2 log likelihood) for all investigated models. The NONMEM objective function improved by 46 points (P ≪ 0.001) when the bolus group (group 1) was given a ke0 distinct from the ke0 used for the three infusion groups (groups 2, 3, and 4). There was no significant benefit to adding distinct values of ke0 for the three infusion groups. Therefore, we selected the model with two values of ke0, one for the bolus (t½ ke0 = 1.2 min) and one for the infusions (t½ ke0 = 2.2 min).
Our final model included estimates of interindividual variability on ke0 and Ce50. In our final model, BIS values of 90 or greater were typically 95.6. BIS values less than 90 were described by a sigmoidal curve with ranging from a maximum value of 74 to a minimum of 25. The Ce50 for the sigmoidal portion was 5.1 μg/ml, with a γ of 3.4. The SDs in the log domain of Ce50 and ke0 were 0.27 and 0.83, respectively, which approximately correspond to the coefficient of variation in the standard domain. Additive residual intraindividual error was 7.2.
Fig. 2
Fig. 2
Image Tools
Figures 2A–D show the raw BIS values and the calculated effect site concentration, Ce, versus time, showing the parallel between the time course of BIS and Ce. Figures 2E–H show the BIS versus Ce. BIS values of 90 or greater are depicted as black circles. BIS values less than 90 are shown as connected lines. Two individuals receiving the bolus injection had baseline BIS values below 90 (fig. 2E), which escaped our efforts to censor the initial unresponsive portion of the BIS response. We did not decrease our censoring value (90) to censor these data because we did not wish to lose informative BIS values in other subjects.
Table 4
Table 4
Image Tools
Table 4 shows the median (range) results for the post hoc Bayesian estimates of Ce50 and t½ ke0 for each individual in the four groups. Table 4 also shows the time of maximum BIS change (tmax,BIS), the time of maximum Ce (tmax,Ce) after the specific propofol infusion, and the absolute value of the difference between these times (terror). The difference in the time of maximum BIS change and the time of the highest Ce can be used to assess the performance of the model. The median terror was 0.6 min or less for all groups. For all patients, the individual results for the post hoc Bayesian estimates of Ce50 and t½ ke0 in the four groups, tmax,BIS and tmax,Ce after the specific propofol infusion, and the absolute value of the difference between these times (tabs,error) can be found in a supplement on the Anesthesiology Web site (http://www.anesthesiology.org).
Fig. 3
Fig. 3
Image Tools
The performance accuracy for the final model was assessed by calculating the PE. Figures 3A–D show the individual PE versus time for the four groups, respectively. Median (range) MDPE and MDAPE for each group was calculated as shown in table 4. The MDAPE was 18% or less for all groups. No group showed significant bias (MDPE). All individual MDPE and MDAPE calculations can be found on the Anesthesiology Web site (http://www.anesthesiology.org).
Back to Top | Article Outline
Bolus Validation Study
Fig. 4
Fig. 4
Image Tools
Ten female patients were included, and all blood samples were analyzed. Demographics were similar in study groups 1 through 4 (age, 34 ± 8 yr; weight, 63 ± 10 kg; height, 163 ± 7 cm). For both models, measured versus predicted propofol plasma concentrations and individual PE% versus time are plotted in figure 4. The mean (SD) PEs were −40.50 (53.08) and −25.90 (36.80)% for the Marsh and Schnider models, respectively. The mean MDAPEs were 60.22 (28.24) and 39.45 (21.38)% for the Marsh and Schnider models, respectively.
Fig. 5
Fig. 5
Image Tools
Figure 5 shows the median arterial concentrations from the bolus. The effect site concentrations were computed based on a t½ ke0 of 2.0 min and predict a peak effect site concentration of 1.9 min, consistent with the time of peak BIS effect in our bolus study (table 4, group 1). All data can be found on the Anesthesiology Web site (http://www.anesthesiology.org).
Back to Top | Article Outline
Validation of the Applied Lag Time in the BIS Measurement
Table 5
Table 5
Image Tools
Table 5 shows the −2 log likelihood values (the NONMEM objective function) when applying different lag times. The best result is obtained when implementing a 10-s lag time in the combined fixed plus sigmoidal model approach.
Back to Top | Article Outline
Method Validation
Fig. 6
Fig. 6
Image Tools
Table 6
Table 6
Image Tools
All data from the 18 volunteers in the original publication by Doufas et al. were included in the model validation. Figure 6A shows the relation between BIS and Cp for the Doufas validation data set. Table 6 shows the −2 log likelihood values (the NONMEM objective function) and the typical values for all investigated models. Doufas' data were best described by a single estimate for those BIS observations of 90 or greater, and a sigmoidal Emax model for BIS values less than 90, with a 10-s delay, exactly as was the case for our data. Most critically, this model estimated a t½ ke0 of 2.1 min, very close to our estimate of 2.2 min for infusions. The SDs in the log domain of Ce50 and ke0 were 0.31 and 1.29, respectively. Additive residual intraindividual error was 9 (i.e., residual variance = 81).
Figure 6B shows the relation between BIS and Ce for each individual, based on the best model and the post hoc Bayesian estimates of Ce50 and ke0. The individual PEs are shown in figure 6C.
Table 7
Table 7
Image Tools
The median (range) post hoc Bayesian parameter estimates for all subjects in the validation data set are shown in table 7, as well as the individual MDPE and MDAPE. The model performed well, with an MDAPE less than 9% and minimal bias (5%). The Ce50 in the validation set is less than in our data set (3.6 vs. 5.1 μg/ml) but is similar to that seen in our bolus group (3.7 μg/ml; table 4). The estimate of γ in the validation set is less than in our data set (1.3 vs. 3.4). For all patients, individual data can be found on the Anesthesiology Web site (http://www.anesthesiology.org).
Back to Top | Article Outline

Discussion

Our goal was to determine whether the equilibration rate between plasma and effect site was influenced by rate of propofol administration. If so, this must be considered when designing drug infusions, and is particularly relevant for target-controlled infusion systems. We tested this by examining the time course of electroencephalographic (BIS) response to a bolus and three infusion rates. The apparent plasma–effect site equilibration for bolus injections is faster than for infusions. This is consistent with the hypothesis of Doufas et al.,5 who speculated that this might be the explanation for the difference between plasma–effect site equilibration in their study versus the study by Schnider et al.3 The accuracy of the pharmacodynamic models is demonstrated by the reasonable values of MDAPE and the agreement between the observed time to maximum BIS changes and the calculated time to the highest effect site concentration after a specific propofol administration (depending on the group).
Because this study aimed at answering the question of why different ke0s are found in the literature, even when using the same pharmacokinetic model, we wanted to explore bolus injections (as used by Schnider et al.3) and various infusion rates covering published work and clinical practice. Our bolus was 2.5 mg/kg given over 10 s. Our slowest administration was an infusion of 2.5 mg/kg given over 3 min, which is 50 mg · kg−1 · h−1. This is lower than the critical infusion rate of 80 mg · kg−1 · h−1 that Kazama et al.12 postulate causes incomplete mixing in the central compartment.
Classically, the sigmoid Emax model has been applied successfully for pharmacodynamic modeling of propofol. However, when observing our raw BIS data, we found that the BIS values showed no change during the initial increase in effect site propofol concentration, until there was an abrupt decrease. This may reflect an acute change in BIS with loss of consciousness, although the article by Doufas et al. examined BIS at loss of consciousness and did not detect any abrupt state change. This might be due to the fact that they used an Observers' Assessment of Anesthesia/Sedation score of less than 2 (= loss of responsiveness to shaking and shouting) as their endpoint.5 Alternatively, this delay may reflect the combination of the averaging algorithm to calculate BIS (i.e., smoothing rate) and the delay in adaptation of one of the artifact rejection preprocessing steps, that rejects large changes in the electroencephalogram as artifact until those changes persist for approximately 5 s (Scott Greenwald, Ph.D., Vice President of Research, Aspect Medical Systems, Newton, MA, verbal personal communication, April 2005). We found that the inclusion of a lag time of 10 s into our model resulted in the best fit. Longer lag times resulted in worse fits of the model to the data. Regardless, the initial observations when the BIS was unambiguously unresponsive to changing propofol concentrations (figs. 1E–H) were censored from the sigmoidal model by introducing a separate model for BIS values of 90 or greater and BIS values less than 90. This model significantly improved the quality of the fit with our data (P ≪ 0.001 for the 2 ke0 model; table 3) and with the validation data set from Doufas et al.5 (P ≪ 0.001 for the models with 10-s delays; table 6).
If propofol administration rate is influencing the time course of drug effect, a model incorporating different estimates of ke0 for different infusion rates would yield a better model fit than a model whereby one ke0 was used to describe the time course of effect site concentration independent of drug administration rate. We found that the model discriminating the bolus injection group from the three continuous infusion groups yielded a significantly better fit than a model with a single ke0 for all four groups. In our bolus group, the value of t½ ke0 was 1.2 min, very close to the value of 1.5 min reported by Schnider et al.,3 whose value of ke0 was primarily based on bolus propofol administration. In the Schnider study, the bolus was administered in a mean time of 18 s (range, 13–24 s), which is somewhat slower than in our study (table 2). In addition, our calculated time of maximum BIS effect in the bolus group was 1.7 min (table 4), exactly the same as reported by Schnider et al.3 Our t½ ke0 of 2.2 min for propofol administration by infusion is in agreement with the value computed from the validation data set of Doufas et al.5 of 2.1 min, both of which reflect fast infusions, but obviously less rapid than bolus injections.
Our findings are consistent with propofol modeling studies using a conventional continuous infusion, which have reported propofol t½ ke0s between 2.3 and 3.5 min. In the first commercial target-controlled infusion device (Diprifusor; AstraZeneca), the kinetic model described by Marsh et al.10 was linked to a t½ ke0 of 2.65 min as described by Schwilden et al.13 and based on slow continuous infusion data. Billard et al.14 administered propofol (continuous infusion of 0.5 mg · kg−1 · min−1 until burst suppression) and measured the electroencephalographic effect using three different electroencephalographic measures. The estimates of t½ ke0 were significantly higher (mean, 2.6 min) for delta power than those for spectral edge 95% (mean, 3.5 min) and Bispectral Index version 1.1 (mean, 3.5 min). White et al.15 used midlatency auditory evoked potentials to measure propofol drug effect, when administered as a continuous infusion (0.5 mg · kg−1 · min−1), and calculated a ke0 of 3.5 min (mean value). All of these effect measures were modeled without including any delay in the electroencephalographic measure.
Ke0 is dependent on the accuracy of the underlying pharmacokinetic model, and if the pharmacokinetic model is biased, the estimate of ke0 will biased by that error. One limitation of our study is that we did not include measurement of propofol plasma concentrations over time in our four groups. Similar to previously published work,5,16–18 we applied the three-compartment model published by Schnider et al.6 to predict the time course of the propofol plasma concentration. In the study by Doufas et al.,5 the pharmacokinetic performance of the propofol pharmacokinetics reported by Schnider et al.3 was validated with rapid arterial samples and found to be accurate. Because our study was designed specifically to test the hypothesis in the Doufas article that the difference in values of ke0 might depend on the difference between a bolus and a continuous infusion of propofol, we did not believe it was necessary to collect arterial blood samples in all patients during the bolus (similar to the infusion scheme used by Schnider in deriving the pharmacokinetics6) or the infusions (where the pharmacokinetics have been validated by both Doufas and Schnider).
However, the literature has not documented the accuracy of the Schnider pharmacokinetic model in the first minutes after bolus injection. Therefore, we added the bolus validation study and documented that neither model performed well in the first 3 min (fig. 4). The failure of both models in the first 3 min is an expected consequence of a flawed fundamental assumption with mammillary compartment models: instantaneous mixing in the central compartment. Conventional two- and three-compartment mammillary compartment models assume that drug added to the central compartment is instantaneously completely mixed, and that this mixed plasma instantaneously appears in the arterial circulation. This is not the case, as extensively reported by Henthorn et al.19–21 and others.22,23 The reason that more complex models have not been integrated into target-controlled infusion systems is that incorporating this into target-controlled infusion algorithms introduces considerable mathematical complexity. Because the problem with model misspecification is limited to the first few minutes after bolus injection, there is limited incentive to develop hybrid models to correct the fundamentally flawed assumption of instantaneous mixing.
The initial error in the pharmacokinetic models is evident for the full 5 min with the Marsh model but seems to be resolving by 3 min with the Schnider model. Therefore, the Schnider model comes closest to fitting the bolus validation data and would be the better of the two models to use for this purpose.
The simulation results from the bolus validation study (fig. 5) show that a t½ ke0 of 2.0 min, coupled to a more accurate pharmacokinetic model, predicts a peak effect site concentration of 1.9 min, as observed with our bolus data (table 4). This is fairly close to the t½ ke0 of 2.2 min that we estimated for the three infusion groups, and demonstrates that most of the difference in ke0 between bolus and infusion administration is an artifact caused by the inability of mammillary models to accurately describe the concentrations in the first few minutes after bolus injection. The actual rate of blood–brain equilibration may not differ between bolus and infusion methods of administration.
This does not preclude the possibility that there are physiologic reasons the rate of plasma–effect site equilibration might change with infusion rate. Upton and Ludbrook24,25 demonstrated that propofol decreases cerebral blood flow in a dose-dependent manner. This might explain why the plasma–effect site equilibration time of propofol could change between a bolus, whose high early concentrations acutely decrease cerebral blood flow, and an infusion, where the changes in cerebral blood flow would be small.
The rate of propofol administration influences the apparent rate of plasma–effect site equilibration. This is partly, and perhaps mostly, a result of the inability of conventional pharmacokinetic models to accurately describe the first few minutes after bolus injection. Target-controlled infusion devices need to select the value of ke0 that will most accurately represent the effect site concentrations over time. The error in the pharmacokinetic model with rapid administration can be partly compensated for by the selection of ke0. If the maximum infusion rate is between 300 and 900 ml/h, this corresponds to the infusion rates in our three infusion groups, and the t½ ke0 should be approximately 2.2 min (ke0 = 0.32). However, if the infusion is closer to our bolus rate (2.5 mg/kg over 10 s ≈ 6,300 ml/h), the faster t½ ke0 of 1.2 min will better approximate the time course of drug effect. These values of ke0 are for use with the Schnider pharmacokinetics. If these data are to be implemented with other pharmacokinetics, then they should be implemented to achieve the predicted time of peak effect of 1.5 min if the device delivers an induction dose of propofol over a minute or less. If the pump is unable to infusion propofol that quickly, selecting ke0 based on a peak effect of 1.8 min is appropriate.2
We conclude that the observation of Doufas et al. is correct, although the explanation is surprising. There is a difference in the apparent rate of plasma–effect site equilibration between propofol boluses and propofol infusions. The difference is mostly caused by misspecification in the pharmacokinetic model over the first few minutes after bolus injection. A pharmacokinetic model that accurately predicted propofol concentration in the first few minutes after bolus injection, and for many hours after infusion, might improve the accuracy of target-controlled infusion administration when targeting the site of drug effect. If propofol infusion rate affects the apparent ke0 even when the pharmacokinetic model is unbiased, there may also be a physiologic basis for dependence of ke0 on propofol infusion rate.
The authors thank Dorine Dyzers, C.N. (Ghent University Hospital, Gent, Belgium), and Scott Greenwald, Ph.D. (Vice-President of Research, Aspect Medical, Newton, Massachusetts), for their help during the preparation of the manuscript. The authors also thank Dennis Fisher, M.D. (President, P Less Than Company, San Francisco, California), for his many helpful comments during the review process.
Back to Top | Article Outline

References

1. Sheiner LB, Stanski DR, Vozeh S, Miller RD, Ham J: Simultaneous modeling of pharmacokinetics and pharmacodynamics: Application to d-tubocurarine. Clin Pharmacol Ther 1979; 25:358–71

2. Minto CF, Schnider TW, Gregg KM, Henthorn TK, Shafer SL: Using the time of maximum effect site concentration to combine pharmacokinetics and pharmacodynamics. Anesthesiology 2003; 99:324–33

3. Schnider TW, Minto CF, Shafer SL, Gambus PL, Andresen C, Goodale DB, Youngs EJ: The influence of age on propofol pharmacodynamics. Anesthesiology 1999; 90:1502–16

4. Struys MM, De Smet T, Depoorter B, Versichelen LF, Mortier EP, Dumortier FJ, Shafer SL, Rolly G: Comparison of plasma compartment versus two methods for effect compartment–controlled target-controlled infusion for propofol. Anesthesiology 2000; 92:399–406

5. Doufas AG, Bakhshandeh M, Bjorksten AR, Shafer SL, Sessler DI: Induction speed is not a determinant of propofol pharmacodynamics. Anesthesiology 2004; 101:1112–21

6. Schnider TW, Minto CF, Gambus PL, Andresen C, Goodale DB, Shafer SL, Youngs EJ: The influence of method of administration and covariates on the pharmacokinetics of propofol in adult volunteers. Anesthesiology 1998; 88:1170–82

7. Beal S, Sheiner L: NONMEM User's Guide. San Francisco, NONMEM Project Group, University of California, 1992

8. Kataria BK, Ved SA, Nicodemus HF, Hoy GR, Lea D, Dubois MY, Mandema JW, Shafer SL: The pharmacokinetics of propofol in children using three different data analysis approaches. Anesthesiology 1994; 80:104–22

9. Varvel JR, Donoho DL, Shafer SL: Measuring the predictive performance of computer-controlled infusion pumps. J Pharmacokinet Biopharm 1992; 20:63–94

10. Marsh B, White M, Morton N, Kenny GN: Pharmacokinetic model driven infusion of propofol in children. Br J Anaesth 1991; 67:41–8

11. Pilge S, Zanner R, Schneider G, Blum J, Kreuzer M, Kochs EF: Time delay of index calculation: Analysis of cerebral state, Bispectral, and Narcotrend indices. Anesthesiology 2006; 104:488–94

12. Kazama T, Ikeda K, Morita K, Kikura M, Ikeda T, Kurita T, Sato S: Investigation of effective anesthesia induction doses using a wide range of infusion rates with undiluted and diluted propofol. Anesthesiology 2000; 92:1017–28

13. Schwilden H, Stoeckel H, Schuttler J: Closed-loop feedback control of propofol anaesthesia by quantitative EEG analysis in humans. Br J Anaesth 1989; 62:290–6

14. Billard V, Gambus PL, Chamoun N, Stanski DR, Shafer SL: A comparison of spectral edge, delta power, and bispectral index as EEG measures of alfentanil, propofol, and midazolam drug effect. Clin Pharmacol Ther 1997; 61:45–58

15. White M, Schenkels MJ, Engbers FH, Vletter A, Burm AG, Bovill JG, Kenny GN: Effect-site modelling of propofol using auditory evoked potentials. Br J Anaesth 1999; 82:333–9

16. Vanluchene AL, Vereecke H, Thas O, Mortier EP, Shafer SL, Struys MM: Spectral entropy as an electroencephalographic measure of anesthetic drug effect: A comparison with Bispectral Index and processed midlatency auditory evoked response. Anesthesiology 2004; 101:34–42

17. Vereecke HE, Struys MM, Mortier EP: A comparison of bispectral index and ARX-derived auditory evoked potential index in measuring the clinical interaction between ketamine and propofol anaesthesia. Anaesthesia 2003; 58:957–61

18. Struys MM, Vereecke H, Moerman A, Jensen EW, Verhaeghen D, De Neve N, Dumortier FJ, Mortier EP: Ability of the Bispectral Index, autoregressive modelling with exogenous input-derived auditory evoked potentials, and predicted propofol concentrations to measure patient responsiveness during anesthesia with propofol and remifentanil. Anesthesiology 2003; 99:802–12

19. Henthorn TK, Avram MJ, Krejcie TC: Intravascular mixing and drug distribution: The concurrent disposition of thiopental and indocyanine green. Clin Pharmacol Ther 1989; 45:56–65

20. Henthorn TK, Krejcie TC, Shanks CA, Avram MJ: Time-dependent distribution volume and kinetics of the pharmacodynamic effector site (letter). J Pharm Sci 1992; 81:1136–8

21. Krejcie TC, Henthorn TK, Niemann CU, Klein C, Gupta DK, Gentry WB, Shanks CA, Avram MJ: Recirculatory pharmacokinetic models of markers of blood, extracellular fluid and total body water administered concomitantly. J Pharmacol Exp Ther 1996; 278:1050–7

22. Upton RN, Ludbrook GL: A physiological model of induction of anaesthesia with propofol in sheep: I. Structure and estimation of variables. Br J Anaesth 1997; 79:497–504

23. Upton RN, Grant C, Martinez AM, Ludbrook GL: Recirculatory model of fentanyl disposition with the brain as the target organ. Br J Anaesth 2004; 93:687–97

24. Upton RN, Ludbrook GL, Grant C, Doolette DJ: The effect of altered cerebral blood flow on the cerebral kinetics of thiopental and propofol in sheep. Anesthesiology 2000; 93:1085–94

25. Ludbrook GL, Visco E, Lam AM: Propofol: Relation between brain concentrations, electroencephalogram, middle cerebral artery blood flow velocity, and cerebral oxygen extraction during induction of anesthesia. Anesthesiology 2002; 97:1363–70

Cited By:

This article has been cited 23 time(s).

Acta Anaesthesiologica Scandinavica
Prospective evaluation of the time to peak effect of propofol to target the effect site in children
Munoz, HR; Leon, PJ; Fuentes, RS; Echevarria, GC; Cortinez, LI
Acta Anaesthesiologica Scandinavica, 53(7): 883-890.
10.1111/j.1399-6576.2009.01995.x
CrossRef
British Journal of Anaesthesia
Pharmacokinetic models for propofol: defining and illuminating the devil in the detail
Absalom, AR; Mani, V; de Smet, T; Struys, MMRF
British Journal of Anaesthesia, 104(2): 262-264.

Pediatric Anesthesia
Pediatric models for adult target-controlled infusion pumps
Anderson, BJ
Pediatric Anesthesia, 20(3): 223-232.
10.1111/j.1460-9592.2009.03072.x
CrossRef
Clinical Pharmacology & Therapeutics
Contributions of PK/PD modeling to intravenous anesthesia
Minto, CF; Schnider, TW
Clinical Pharmacology & Therapeutics, 84(1): 27-38.
10.1038/clpt.2008.100
CrossRef
Biopharmaceutics & Drug Disposition
Simple Pharmacometric Tools for Oral Anti-Diabetic Drug Development: Competitive Landscape For Oral Non-Insulin Therapies in Type 2 Diabetes
Samtani, MN
Biopharmaceutics & Drug Disposition, 31(): 162-177.
10.1002/bdd.700
CrossRef
Anaesthesist
Target-controlled infusion Clinical relevance and special features when using pharmacokinetic models
Ihmsen, H; Schraag, S; Kreuer, S; Bruhn, J; Albrecht, S
Anaesthesist, 58(7): 708-715.
10.1007/s00101-009-1575-3
CrossRef
British Journal of Anaesthesia
Study of the time course of the clinical effect of propofol compared with the time course of the predicted effect-site concentration: performance of three pharmacokinetic-dynamic models
Coppens, M; Van Limmen, JGM; Schnider, T; Wyler, B; Bonte, S; Dewaele, F; Struys, MMRF; Vereecke, HEM
British Journal of Anaesthesia, 104(4): 452-458.
10.1093/bja/aeq028
CrossRef
British Journal of Anaesthesia
Dose-dependency of pharmacokinetic/pharmacodynamic parameters after intravenous bolus doses of cisatracurium
Chen, C; Yamaguchi, N; Varin, F
British Journal of Anaesthesia, 101(6): 788-797.
10.1093/bja/aen308
CrossRef
International Journal of Adaptive Control and Signal Processing
Adaptive disturbance rejection control. for compartmental. systems with application to intraoperative anesthesia influenced by hemorrhage and hemodilution effects
Volyanskyy, K; Haddad, WM; Bailey, JM
International Journal of Adaptive Control and Signal Processing, 23(1): 1-29.
10.1002/acs.1029
CrossRef
British Journal of Anaesthesia
Evaluation of the predictive performance of four pharmacokinetic models for propofol
Glen, JB; Servin, F
British Journal of Anaesthesia, 102(5): 626-632.
10.1093/bja/aep043
CrossRef
Journal of Chromatography B-Analytical Technologies in the Biomedical and Life Sciences
Formation of the N-methylpyridinium ether derivative of propofol to improve sensitivity, specificity and reproducibility of its detection in blood by liquid chromatography-mass spectrometry
Thieme, D; Sachs, H; Schelling, G; Hornuss, C
Journal of Chromatography B-Analytical Technologies in the Biomedical and Life Sciences, 877(): 4055-4058.
10.1016/j.jchromb.2009.10.010
CrossRef
Anesthesia and Analgesia
The Dynamic Relationship Between End-Tidal Sevoflurane Concentrations, Bispectral Index, and Cerebral State Index in Children
Fuentes, R; Cortinez, LI; Struys, MMRF; Delfino, A; Munoz, H
Anesthesia and Analgesia, 107(5): 1573-1578.
10.1213/ane.0b013e318181ef88
CrossRef
British Journal of Anaesthesia
Pharmacodynamic modelling of the bispectral index response to propofol-based anaesthesia during general surgery in children
Jeleazcov, C; Ihmsen, H; Schmidt, J; Ammon, C; Schwilden, H; Schuttler, J; Fechner, J
British Journal of Anaesthesia, 100(4): 509-516.
10.1093/bja/aem408
CrossRef
Acta Anaesthesiologica Scandinavica
Predictive ability of propofol effect-site concentrations during fast and slow infusion rates
Sepulveda, PO; Cortinez, LI; Recart, A; Munoz, HR
Acta Anaesthesiologica Scandinavica, 54(4): 447-452.
10.1111/j.1399-6576.2009.02183.x
CrossRef
Journal of Pharmacokinetics and Pharmacodynamics
A two-compartment effect site model describes the bispectral index after different rates of propofol infusion
Bjornsson, MA; Norberg, A; Kalman, S; Karlsson, MO; Simonsson, USH
Journal of Pharmacokinetics and Pharmacodynamics, 37(3): 243-255.
10.1007/s10928-010-9157-1
CrossRef
Academic Emergency Medicine
Subdissociative-dose Ketamine versus Fentanyl for Analgesia during Propofol Procedural Sedation: A Randomized Clinical Trial
Messenger, DW; Murray, HE; Dungey, PE; van Vlymen, J; Sivilotti, MLA
Academic Emergency Medicine, 15(): 877-886.
10.1111/j.1553-2712.2008.00219.x
CrossRef
Pediatric Anesthesia
Which model for propofol TCI in children
Constant, I; Rigouzzo, A
Pediatric Anesthesia, 20(3): 233-239.
10.1111/j.1460-9592.2010.03269.x
CrossRef
British Journal of Anaesthesia
Bispectral index-guided induction of general anaesthesia in patients undergoing major abdominal surgery using propofol or etomidate: a double-blind, randomized, clinical trial
Petrun, AM; Kamenik, M
British Journal of Anaesthesia, 110(3): 388-396.
10.1093/bja/aes416
CrossRef
Anesthesiology
Response Surface Modeling of the Interaction between Propofol and Sevoflurane
Schumacher, PM; Dossche, J; Mortier, EP; Luginbuehl, M; Bouillon, TW; Struys, MM
Anesthesiology, 111(4): 790-804.
10.1097/ALN.0b013e3181b799ef
PDF (2382) | CrossRef
Anesthesiology
Early Phase Pharmacokinetics but Not Pharmacodynamics Are Influenced by Propofol Infusion Rate
Masui, K; Kira, M; Kazama, T; Hagihira, S; Mortier, EP; Struys, MM
Anesthesiology, 111(4): 805-817.
10.1097/ALN.0b013e3181b799c1
PDF (1607) | CrossRef
Anesthesiology
Take It to the Limit1 (One More Time)
Fisher, DM
Anesthesiology, 107(3): 367-368.
10.1097/01.anes.0000279485.30953.a9
PDF (233) | CrossRef
Current Opinion in Anesthesiology
Pharmacokinetic–pharmacodynamic modeling in anesthesia, intensive care and pain medicine
Sadean, MR; Glass, PS
Current Opinion in Anesthesiology, 22(4): 463-468.
10.1097/ACO.0b013e32832c3c6c
PDF (115) | CrossRef
European Journal of Anaesthesiology (EJA)
Individual titration of propofol plasma target improves anaesthetic stability in patients undergoing major abdominal surgery: a comparison with manually controlled infusion: Retracted
Mayer, J; Boldt, J; Triem, JG; Schellhaa, A; Mengistu, AM; Suttner, S
European Journal of Anaesthesiology (EJA), 25(9): 741-747.
10.1017/S0265021508004328
PDF (530) | CrossRef
Back to Top | Article Outline

Supplemental Digital Content

Back to Top | Article Outline

© 2007 American Society of Anesthesiologists, Inc.

Publication of an advertisement in Anesthesiology Online does not constitute endorsement by the American Society of Anesthesiologists, Inc. or Lippincott Williams & Wilkins, Inc. of the product or service being advertised.
Login

Article Tools

Images

Share