| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
BACKGROUND: There are two recognized methods of estimating the brain concentrations of IV anesthetic drugs: (i) use of pharmacokinetic/pharmacodynamic (PK/PD) modeling of drug effect, from arterial concentrations and electroencephalogram changes, and (ii) direct measurement of the uptake of drug in the brain, by simultaneously measuring arterial and jugular concentrations and cerebral blood flow (mass-balance method). These two methods have not been directly compared. Because an accurate estimate of the time taken for transfer of anesthetic drug from arterial blood to its effect-compartment in the brain is critical for accurate effect-compartment dosing in IV anesthesia, we compared the PK/PD and mass-balance methods for propofol, methohexital, and ketamine in a sheep model. METHODS: After instrumentation with arterial and sagittal-sinus cannulae, electrocorticogram, and sagittal sinus Doppler flow measurement seven adult sheep were given a random sequence of short anesthetic infusions with methohexital, ketamine, and propofol. Multiple blood samples were taken for measurement of the time course of the drug concentrations, and the electrocorticogram processed (approximate entropy, for propofol and methohexital and percentage high frequency time, for ketamine) to numerically quantify drug effect.
RESULTS: Using the PK/PD method the t CONCLUSIONS: Both methods of calculating the delay in transfer of drug from arterial blood to brain give similar values. Methohexital crosses into the brain much faster than either propofol or ketamine.
The pharmacokinetic (PK) profiles of the IV anesthetics propofol, ketamine, and methohexital drugs are well described (1) and are used to formulate dose regimens to optimize anesthetic depth. Because of the equilibration delay between plasma and brain compartments, some of the currently available pharmacologic models directly target brain concentration by incorporating plasma effect-compartment kinetics (2). The utility of these models for targeting effect-compartment concentrations and for optimizing anesthetic depth is dependent upon the accuracy of the incorporated PK/pharmacodynamic (PK/PD) data.
The kinetics of propofol uptake into the brain has been a matter of some debate. Animal and human studies using mass-balance principles have suggested that propofol exhibits a large and prolonged (t The primary aim of the present study was to directly compare the brain uptake kinetics of propofol using these two different methods: mass-balance and PK/PD modeling. The former method is based on calculation of the net flux of the drug into the brain from the difference in arteriovenous concentrations and the cerebral blood flow. PK/PD analysis is a well-established method for linking the measured blood concentration profile with drug effect, from which an estimate of effect-compartment (brain) concentration can be derived. The PD characteristics of ketamine and methohexital have been less well studied than those of propofol. In the only previous article that has investigated the PK/PD of induction of anesthesia with ketamine, Keo and effect-compartment concentrations were not calculated (9). Similarly, we are not aware of any previous studies that have modeled the PK/PD characteristics of methohexital induction. The second aim of this study, therefore, was to compare the PD profiles of ketamine and methohexital with that of propofol.
After approval by the Animal Ethics Committee of the University of Adelaide 7 adult Merino sheep, each weighing approximately 50 kg, were studied during anesthesia with propofol, ketamine, and methohexital. Under halothane anesthesia, the sheep were instrumented as previously described (10). This included insertion of a sagittal-sinus Doppler flowprobe for cerebral blood flow measurements, a femoral arterial intravascular cannula for arterial blood sampling, a sagittal-sinus catheter for sampling of cerebral effluent blood and a parasagittal linear array of stainless steel electrodes for quantifying the approximate entropy of the electrocorticogram (ECoG). The electrodes were spaced 2 mm apart, and penetrated 1–2 mm into the outer layers of the cortical gray matter. The ECoG was obtained using two electroencephalogram (EEG) monitors (A1000, Aspect Medical Systems, Natick, MA) with the output voltage from the cortical surface reduced electronically, so that the input into the EEG monitors was similar in magnitude to that observed with scalp electrodes. The signal was digitally sampled at 256/s, and downloaded without further filtering to a computer for further analysis. After allowing at least 24 h for the animal to recover and to check the stability of the recording apparatus, the sheep were given a series of short IV anesthetics. Propofol (200 mg), ketamine (300 mg) and methohexital (150 mg) were administered over 2 min in random order and the sheep were then allowed to recover fully (at least 2 h) before the administration of the next drug. The relatively short recovery time allowed all anesthetics to be delivered on the same day, eliminating deterioration in Doppler flow data quality over successive days. All variables were recorded for 20 min after induction of anesthesia. Arterial and sagittal sinus blood samples (1 mL) were taken at intervals of 15 s for 3 min, and then every minute up to 20 min, heparinized (25 IU) and stored at –5°C. The concentration of each drug was later assayed using high performance liquid chromatography (HPLC) and the fluorescence detection method. Sample sizes were 500 µL for each drug. Propofol concentrations in whole blood were assayed using a modification of a technique described by Mather et al. (11) based on solvent (n-hexane) extraction of samples and HPLC. A silica column was used, and the buffer was acetonitrile: 0.5% acetic acid (55:45). Fluorescence detection was used with excitation and emission wavelengths of 274 and 300 nm, respectively. Thymol was the internal standard and the limit of sensitivity was approximately 0.02 µg/mL. Ketamine was assayed using solvent (n-hexane) extraction of samples and HPLC with a C18 column and using an acetonitrile-phosphate buffer (42:58), pH-8, with UV detection (210 nm). Lidocaine was used as an internal standard, and the limit of sensitivity was approximately 0.1 µg/mL. Methohexital was assayed using solvent (n-hexane) extraction and HPLC with a C18 column and acetonitrile-phosphate buffer (38:62), pH-5.5, with UV detection (230 nm). Pentobarbital was used as an internal standard, and the limit of sensitivity was approximately 0.15 µg/mL. The linear calibration ranges were 1.25–40 µg/mL for propofol, 0.5–4 µg/mL for ketamine and 10–50 µg/mL for methohexital.
MA-Balance Measurement of Brain Concentration
Quantification of Anesthetic Effect on the ECoG The approximate entropy of the ECoG was calculated from non-overlapping 5 s ECoG segments in the standard manner, as described by Bruhn et al. (12), using a sample rate of 128/s, an embedding of m = 2, and noise threshold r = 0.2 x SD. Because the algorithm incorporates a lower threshold for noise, approximate entropy has been shown to seamlessly decrease as burst suppression increases. There is no previously published EEG method of quantifying ketamines effect. On observation of the ECoG signals it was apparent that anesthesia with ketamine is characterized by an ECoG pattern in which short episodes of low frequency, high amplitude oscillations alternate with the preexisting high frequency, low amplitude activity (see Fig. A1 in Appendix). With increasing depth of anesthesia, the proportion and length of the high frequency episodes decrease. The reduction in the percentage of high frequency component of the ECoG was quantified for ketamine using wavelet analysis (see Appendix for a more detailed explanation). Because these episodes of low frequency activity are localized in time, a wavelet transform was applied to the ECoG using the MatLab® function cwt (continuous wavelet transform, Morlet wavelet). A scale (approximately "1/frequency") of 48 was heuristically chosen as the high-low frequency boundary. The proportion of high frequency in the ECoG was quantified as the percentage of time that the wavelet coefficients for scales <48 were above an arbitrary coefficient threshold of 13, calculated as a running 10 s average. This method of quantification of ketamines effect was found to track concentrations well and was consistent across different animals (Fig. 3). A Matlab function for this method may be found in the Appendix.
(PK/PD) Modeling of Effect-Compartment Brain Concentration
where Cart is the arterial concentration of the anesthetic, Ceff is the effect-compartment concentration of the drug, and Keo is the first order rate constant for efflux from the effect-compartment. We estimated the effect-compartment concentration of each drug by iteratively running this model with a series of Keo steps. For each iteration, a nonlinear inhibitory sigmoid Emax curve was fitted to the data using the MatLab "nlinfit" function (nonlinear least-squares data fitting), defined as follows:
where effect is the approximate entropy or high frequency component of the ECoG, the Emax and Emin are the maximum and minimum ECoG effects recorded for each animal, the EC50 is the drug concentration at which the ECoG effect is midway between this maximum and minimum, the Ceff is the concentration at the effect-compartment, and
Data Handling and Statistical Analysis
Brain Concentration Time-Course The summary PK/PD parameters for each drug are shown in Tables 1–3. The time-course of concentration changes were similar for propofol and ketamine, as reflected in the similar lag times to peak brain concentration for both analysis methods (Table 2 and Fig. 2). This implies similar t Keo values for each method, 2.7 ± 1.1 min for propofol and 2.0 ± 0.4 min for ketamine. For methohexital, peak concentration was reached slightly more rapidly for the PK/PD analysis (129 ± 3 s compared with 156 ± 7 s, P < 0.05). The t Keo for methohexital was rapid (0.3 ± 0.1 min) and significantly faster than both propofol and ketamine (P < 0.05). This was reflected in a relatively short latency to peak brain concentration (Table 2 and Fig. 2).
Peak Brain Concentrations—Comparison of Analysis Methods
Anesthetic Effects on the ECoG
Brain Concentration Time-Course The brain uptake kinetics for the two measurement methods were compared in the present study by analyzing the brain concentration time-course for each method.
MA-balance analysis of brain concentration and PK/PD modeling gave similar brain uptake kinetics for propofol, with a t
The similarity in the time-course for changes in brain concentration between the analysis methods for propofol suggests that the disequilibrium between propofol blood concentrations and effect is due largely to the physicochemical diffusion properties of the drug, rather than to drug-receptor interaction kinetics. If there were a significant lag due to drug-receptor interactions, one would have expected the PK/PD effect-compartment data to lag behind the mass-balance-measured concentration profile. The time resolution of the ECoG analysis in our study was 10 s (around the same order of magnitude as the drug sampling time); thus the analysis lag did not contribute significantly to total t Methohexital had a very short lag time and rapid elimination compared with the other two drugs. These results indicate that methohexital has faster induction and emergence characteristics than propofol, the current drug of choice for rapid IV anesthesia.
Peak Brain Concentrations The calculated partition coefficients in our study differ from published values (15). This is because the latter are based solely on octanol:water partition coefficients, whereas our values are based on real brain concentrations, which are dependent on drug lipophilicity and the extent of protein binding. The slightly higher partition coefficient for ketamine compared with propofol is therefore explained on the basis of the low protein binding of ketamine (15). Similarly, the very low partition coefficient for methohexital may be explained by its comparatively high protein binding and modest lipophilicity. An alternative explanation for the difference in calculated brain concentrations between the methods of analysis is that the PK/PD-derived effect compartment concentrations were not representative of the overall brain concentration. Total and regional brain concentrations of propofol have been measured previously in rats using HPLC (16,17). However, these studies showed a uniform regional distribution of propofol throughout the brain, with the possible exception of the midbrain, which showed increased concentrations at low infusion rates (17). Overall, these findings indicate that the differences observed in the present study are unlikely to be attributable to a regional distribution effect. This conclusion is further supported on grounds that both measurement methods were effectively sampling from the cerebral cortex, not the entire cerebrum.
As shown in Figure A1, at low levels of ketamine (left side), the "awake" (low-amplitude, mixed frequency) electrocorticogram (ECoG) pattern is interrupted by episodes of large amplitude slow waves (1–5 Hz), lasting about 1 s. At the peak of the ketamine effect these slow-wave episodes have increased to take up most of the ECoG (right side). This can be easily quantified in frequency and time using the Morelet wavelet transform - shown in the lower figures. A Matlab function to achieve this is: function [perc] = ket_wave_funct(x,Fs,thres);
% Does cwt and then determines percentage time that the ECoG ... % signal is in "slow wave" mode (vs "awake/REM" mode) % INPUTS: % (1) ECoG signal (x) % (2) Sampling frequency (set for Fs = 256/s) % (3) Threshold for determining presence of slow wave episodes % OUTPUTS percentage time (perc) that the ECoG is in slow-wave mode % ——————————————————– % Transforms to continuous wavelet - morlet - low frequencies only ... % (i.e., scales 48 to 64) coefs = cwt(x,[48:8:64],morl); zin = abs(coefs(3,:)); % Absolute magnitude of wavelets len = length(zin); % Make a smooth envelope of wavelet coeffs zins = medfilt1(zin,40); % Separate series into slow segments and fast segments % The threshold may need optimization depending on the ECoG amplitude a = find(zins > thres); %find the bits with lots of low frequency power b = find(zins <= thres); %find the bits with no low frequency power hi = x + 1; lo = x; hi(b) = 0; % set all no LF bits to zero lo(a) = 0; % set all lots of LF bits to zero (just for testing graph) % Go through whole file and identify low sections and %time low> 0 llo = length(lo); perc = []; seg = Fs*5; %five second running averages kk = 1; for i = 1:seg:seg*fix(llo/seg)-seg; zzz = lo(i:i+seg); % 5 s segment aa = find(zzz = 0); perc(kk) = length(aa)/seg; kk = kk + 1; end %% Do figure to check visually OK % figure (9); plot(perc,k); % xlabel(Time (s)); % ylabel(Percentage slow waves); return
Supported by the Neurological Foundation of New Zealand. Reprints will not be available from the authors.
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|