Scientific and Technological Research

Randomized cross-over clinical trial comparing two pharmacokinetic models of propofol using entropy indices^*

Ensayo clínico cruzado y aleatorizado para comparar dos modelos farmacocinéticos de propofol usando índices de entropía

Oscar Leonardo Mosquera-Dussán^a, Andrés Cárdenas^a, Daniel Alfonso Botero-Rosas^a,**, Andrés Yepes^a, Henry Oliveros^a, Rubén Henao^b, Fernando Ríos^a

^a Grupo PROSE1M, Facultad de Medicina, Doctorado en Biociencias, Universidad de La Sabana, Chía, Colombia
^b Grupo Calidad, Seguridad, y Educación en Salud, Sociedad Colombiana de Anestesiología y Reanimación (S.C.A.R.E.), Bogotá D.C., Colombia

^* Please cite this article as: Mosquera-Dussán OL, Cárdenas A, Botero-Rosas DA, Yepes A, Oliveros H, Henao R, et al. Ensayo clínico cruzado y aleatorizado para comparar dos modelos farmacocinéticos de propofol usando índices de entropía. Rev Colomb Anestesiol. 2016;44:193-200.
^** Corresponding author at: Facultad de Medicina, Doctorado en Biociencias, Universidad de La Sabana, Campus del Puente del Común, Km. 7, Autopista Norte de Bogotá, Chía, Colombia.
E-mail address daniel.botero@unisabana.edu.co (D.A. Botero-Rosas).

Abstract

Introduction: There are two different pharmacokinetic models (Marsh and Schnider) for the administration of total intravenous anesthesia with propofol, the parameter differences could have some impact over the depth of anesthesia.

Objective: To determine if there is a significant difference in the variability of depth of anesthesia suggesting that one model is superior in achieving a more stable and predictable depth of anesthesia during surgery.

Methods: A cross-over clinical trial was conducted on 16 healthy patients programmed for upper or lower limb ambulatory orthopedic surgery. Patients were randomly assigned to (i) effect site target controlled infusion of propofol with Marsh model at a target concentration of 2.5 μg/ml for 20 min, a 20 min washout, then propofol was administered with Schnider model at the same effect site target for the reminder of the surgery, or (ii) the reverse sequence. Differences in variability of depth of anesthesia, were assessed by comparing records of spectral entropy indices during surgery through an unpaired t-test.

Results: There was no evidence of significant difference in the mean variances of either spectral entropy indices between the two models (p-value: 0.57 for State Entropy, p-value: 0.51 for Response Entropy).

Conclusion: The study suggests that both pharmacokinetic models are equivalent in terms of stability of depth of anesthesia. It is important to keep testing determinants of the efficacy of the models in different types of population because their behavior according to individual characteristics of patients or variables such as cost-effectiveness could end up tilting the scale.

Keywords: Pharmacokinetics, Entropy, Deep sedation, Anesthesia, Propofol.

Resumen

Introducción: Hay dos modelos farmacocinéticos diferentes para la administración de la anestesia total intravenosa con propofol (Marsh y Schnider), las diferencias entre los parámetros podrían tener algún impacto sobre la profundidad anestésica.

Métodos: Estudio clínico cruzado, controlado y aleatorizado llevado a cabo en 16 pacientes programados para cirugía ambulatoria de ortopedia. Los pacientes fueron asignados aleatoriamente a i) infusión controlada por objetivo de propofol con el modelo de Marsh a una concentración objetivo en sitio de efecto de 2.5μg/ml durante 20 minutos, 20 minutos de periodo de lavado, seguido de infusión de propofol con modelo de Schnider a la misma concentración objetivo; o ii) la secuencia inversa. La diferencia en variabilidad de profundidad anestésica fue evaluada mediante la comparación de registros de índices de entropía con una prueba t no pareada.

Resultados: No se encontró evidencia de diferencias significativas de la varianza media en los índices de entropía espectral asociada a los modelos (valor-p: 0.57 para entropía de estado, valor-p: 0.51 para entropía de respuesta).

Conclusión: El estudio sugiere que ambos modelos son equivalentes en términos de estabilidad de profundidad anestésica. Es importante continuar estudiando la eficacia de los modelos en diferentes tipos de población, dado que su comportamiento según características individuales de los pacientes o variables cómo costo-efectividad podrían inclinar la balanza.

Palabras clave: Farmacocinética, Entropía, Sedación profunda, Anestesia, Propofol.

Introduction

Currently, there are multiple efficacious anesthetic agents, inhaled as well as intravenous, and both types allow detailed titration and fast recovery with a good safety profile. However, total intravenous anesthesia has been demonstrating possible advantages over inhaled techniques, regarding not only safety during surgery and post-operative well being but also in terms of convenience of the administration technique in some specific circumstances (airway intervention, neuroanesthesia) and other issues like environmental impact. This is why total intravenous anesthesia has been gaining popularity in clinical practice, especially since the introduction of propofol, its use is increasingly widespread.^1-4

There are two pharmacokinetic models for the administration of total intravenous anesthesia (TIVA) with propofol, Marsh and Schnider models, which take into account interactions between body compartments to modify infusion rate and, in theory, maintain a constant plasmatic concentration.⁵ Although there is no evidence of clinically important differences between the two models, neither is it clear which of the two is able to predict more accurately the plasma concentrations, it is evident that they differ in the calculations of compartment volumes as well as diffusion velocity between them, which results in significant differences in total propofol dose administered, infusion velocity, and therefore estimated plasma and effect site concentration,⁶ differences that could have some impact over the depth of anesthesia. Table 1, presents respective equations for each model implemented in the B.Braun Space Infusion Pumps,⁷ corresponding parameters are: V1, distribution volume for central compartment. V2, V3, distribution volume for fast and slow peripheral compartments. k10, velocity constant for elimination rate. k12, k13, velocity constants from central compartment to peripheral compartments. k21, k31, velocity constants from peripheral compartments to central compartment.

As of today, there is no gold standard to quantitatively measure the state of consciousness and the depth of anesthesia, ordinarily their monitoring is based on judgment of the anesthesiologist, based on variables of autonomic activity, respiratory cycles and pupil size.¹⁰ However, there has been important progress in the analysis of electroencephalographic signals by techniques like bispectral index (BIS)¹¹ and spectral entropy indices (M-Entropy)¹² which have demonstrated good correlation with sedation levels evaluated clinically and different steps of anesthesia.^13,14 The M-Entropy module, particularly the Response Entropy (RE) index, was considered a better predictor of patient response to painful stimuli than BIS.¹⁵ This is why entropy measurement can be considered an indirect sign of depth of anesthesia, and allows a type of monitoring able to detect variations even within the same plane of anesthesia which makes it an ideal tool to evaluate dynamically and quantitatively the real repercussions of the differences observed between the two pharmacokinetic models for administration of anesthesia with propofol.

Methods

An uniform cross-over, controlled, clinical trial was designed and approved by the ethics committee of Clínica Universidad de la Sabana (C.U.S.), pilot trials were conducted with the surgical team involved in the procedures looking for homogeneity on criteria and management from all personal involved. The study took place from February to July 2015, it was registered in http://www.clinicaltrials.gov with study id number: ECC001-2015. Sixteen patients aged 18-65 years were included in the study, all with at least 8 h of fasting, who had not previously been taking medication with action on the central or autonomous nervous system (benzodiazepines, beta blockers, calcium channel blockers, alpha 2 agonists), with no communicative or auditory difficulties, with no propofol allergy, with pre-anesthetic evaluation (medical and nurse) by the Anesthesiology department of C.U.S., with American Society of Anesthesiology classification (ASA) I, programmed for upper or lower extremity ambulatory orthopedic surgery, who required regional anesthesia by sonogram guided nerve block combined with general anesthesia. For all patients, informed consent (surgical, anesthetic, and research inclusion) was signed and verified, check list was performed, adequate oxygen pressure, connected anesthetic vaporizers, secured and correctly assembled breathing system, emergency oxygen tank, adequate state of fluxometers, circuit, unidirectional valves, pressure relief valves, CO2 absorbent canister and finally adequate function of ventilator and monitors were noted. Surgical Safety Checklist: Safe Surgery Saves Lives was implemented. Monitoring with electrocardiography, non-invasive blood pressure, pulse oximetry, capnography and entropy was initiated. Adequate electrode impedance and reliable reading from the monitor (Datex/Ohmeda) was ensured.

All patients were premedicated with 0.03 mg/kg IV midazolam and 2 μg/kg fentanyl, area was prepared with chlorhexidine in sterile fashion. Nerve was located using stimulating needle, nerve stimulator and extension line (ref. 5194.103 30°-21G-L 100mm, Locoplex®). For lower limb, advance current of 1.2 mA was used until nerve response was achieved at 0.5 mA to localize femoral and obturator nerves by anterior approach and sciatic nerve via posterior approach. For upper limb, a similar procedure was performed with sup-raclavicular or axillar approaches for nerve blockade. 20 ml of local anesthetic was infused over each nerve, using 0.375% levobupivacaine without epinephrine and 1% lidocaine with weight adjusted dose without reaching a maximal dose of 150 mg for levobupivacaine and 5 mg/kg for lidocaine. Absence of pain, paresthesia, dysesthesia or difficult infusion was verified during the administration of local anesthetic. Additionally, position of the needle was assessed by ultrasound guidance. To ensure randomization, an equilibrated list of values 0 and 1 was generated using randomization function of Matlab (The Mathworks Inc., Massachusetts, U.S.A.) version 2015a¹⁶, where 0 meant starting with Marsh model and 1 starting with Schnider. Then digital data collection of SE and RE was started every 5 s using the Datex/Ohmeda collect^TM serial communication acquisition software for later processing using Matlab. M-Entropy indices were registered during a 4min baseline period without propofol, then, the anesthesiologist was informed with which model to initiate the procedure according to the order generated by the list. A technique of TIVA through target controlled infusion (TCI) with B Braun laboratory infusion pump was used.⁷ Induction was conducted with 5 ng/ml remifentanil (Ultiva, Minto model) and after 3min, 1% propofol-lipuro infusion started at an effect site target concentration of 2.5 μg/ml with either Marsh or Schnider model. An additional 1 mg/kg propofol bolus was administered for placement of the laryngeal mask. TCI with propofol was administered for a period of 20 min from the beginning of surgery (period A) after which the infusion was suspended for 20 min of washout period (W). Next, the propofol infusion was configured and started again with the remaining pharmacokinetic model until the end of the surgery (period B). W was established by observation of plasma levels estimated by the infusion pump during pilot trials (Fig. 1).

Sample size was determined using the parameters shown in Table 2. Using Eq. (1) for integration of variances.¹⁸

Considering type I error (α = 0.05), and 80% of power (β = 0.2) to detect a expected difference in entropy values between groups of E = 10 units, the parameters from the normal standard probability distribution A (Eq. (2)) and B (Eq. (3)) were determined using Matlab 2015a.¹⁶

Then calculation of the sample size is (Eq. (4))

For statistical analysis, output variances for each patient from the periods associated to pharmacokinetic models were averaged, an unpaired t-test was performed and implemented (Eq. (5)) to confirm the absence of carryover according to Wellek and Blettner.¹⁷

Where:

m: number of patients in Marsh-Schnider group.
n: number of patients in Schnider-Marsh group.
N: total number of patients.
(X): mean of the sum of variances in Marsh-Schnider group. ]]> (Y): mean of the sum of variances in Schnider-Marsh group.

Then, the differences between models were evaluated. Statistical significance was estimated using an unpaired t-test, this time comparing the mean of the differences of variances of the entropy indices for each model in every patient as shown in Eq. (6):

Where:(X): mean of the difference of variances in Marsh-Schnider group

(Y): mean of the difference of variances in Schnider-Marsh group

Finally, post hoc power (1 - β) was calculated as a function of different degrees of mean variance difference (δ).¹⁹ Data processing and statistical tests were done with a function developed on Matlab (analisis_pk.m) which takes as input arguments the times registered during data recollection and M-entropy series to calculate the variance of each period, then statistics test were applied.

Results

Eight patients were allocated initially to each group and underwent the described intervention. Data from two patients was excluded because the surgical time was not enough to complete the study design. The demographics of the population studied are presented in Table 3, organized according to the groups associated to the pharmacokinetic models.

Carryover effect pre-test

With the sum of the variances of each group, the pre-test for carryover effect between periods was implemented. The cumulative effect is reported in the form of median ± standard deviation, for SE and sequence Marsh-Schnider, SE-C(X) = 42.80 ±57.14, sequence Schnider-Marsh SE-C(Y) = 53.70 ±48.52. Cumulative effect for RE-C(X) = 84.89 ±39.63, and RE-C(Y) = 91.93 ±48.90. Fig. 3 shows the box-plot of the cumulative effect for each group in the RE and SE series (Cx-RE, Cx-SE). There was no evidence of significant carryover effect between periods, two sided unpaired t-test for SE (p-value = 0.54, T=0.63) and RE p-value = 0.92, T = -0.10).

Statistical test of differences

The difference of variance for each period in both groups was obtained (Table 4), Fig. 3 also displays box-plot of the differences between periods related to each Datex-Ohmeda indice (Dx-SE, Dx-RE).

The statistical Td (Eq. (2)) was calculated to detect significant differences in the SE an RE series between the two pharmacokinetic models, no significant difference was found in the mean variance of entropy indices between Marsh and Schnider pharmacokinetic models, two sided unpaired t-test for SE (p-value = 0.57, T = -0.58) and RE (p-value = 0.51, T = -0.68).

Lastly, post hoc power of the study to detect differences between the mean variance of the models with SE and RE entropy was calculated. A vector differences from -125 to 125 with 2.5 step increments in variance units was used. Fig. 4 shows the corresponding power curve for differences in mean variance and mean standard deviation. The particular case for a difference of 100 variance units, therefore 10 standard deviation units corresponds to a power of 84.66% for SE and 85.99% for RE.

Discussion

To the best knowledge of the authors, no studies were found that evaluated the mean variability of the general anesthesia in relation to the pharmacokinetic model used during total intravenous anesthesia with propofol. This was an initial limiting factor to determine the sample size and it was necessary to use estimated values for the parameters from the results of pilot trials. Cross-over design allows a reduced sample size because every subject is its own control, which lowers the inter-individual variability and the effect of covariates such as weight, age, and height. It was not considered necessary to maintain a double-blinding since none of the models is a "standard of care" which means there is no clear expectancy from the anesthesiologist about the results. Additionally the process was standardized so there was no supplementary management implemented, this way, the dosage of the infusions was controlled by the infusion pump according to the pharmacokinetic models and standard analgesic management was provided. In relation to possible bias from co-interventions, the initial pharmacokinetic model was revealed to the anesthesiologist right before onset of induction and after nerve blockade for regional anesthesia. Finally, response to intervention was measured and registered in an automated fashion without mediation from the anesthesiologist, and the M-entropy indices are an endpoint obtained from objective data which makes it less susceptible to bias.

There was one period of intervention with atypical values of variance in patient 3 (induction with Marsh model, 60 kg, 1.72 m, 22 ages, male). Further inspection of the entropy graphs series (SE, RE) revealed a temporal gap in the indices in relation to the clinical assessment (airway placement), therefore, definition of the times for period A included entropy values (SE, RE) associated with an awake patient (Fig. 5) which produced in those cases the atypical values of variance. Data that exceeded the statistical threshold defined by quartiles (Q3 + 3*(Q3 - Q1)) were initially considered atypical and excluded from the statistical analysis, but it was decided to not discard the information that outlier values may provide, data marked as atypical was included for analysis.

Median variance and respective standard deviation in entire dataset for SE in Marsh period was 27.93 ±78.97, in Schnider period was 20.26 ± 39.45; for RE in marsh period was 34.01 ±123.08, and in Schnider period was 36.03 ±60.59. An increment in the standard deviation of data point corresponding to the Marsh period were observed, although significant differences between pharmacokinetic models were not found (p >0.05), it was interesting to note that the outlier values produced a bigger impact on the groups that initiate the procedure with the Marsh model, this could be related to the difference in the pharmacodynamic parameter Keo, the lower value used in the marsh model (0.26 min^-1) in this study could provide an erroneous estimation of the propofol effect time, leading to an incorrect estimation of propofol dose. We used effect site target concentration for both models, it has been suggested that if the Marsh model is used in effect site targeting mode, it should be used with the faster keo of 1.2 min^-1 recommended by Struys and colleagues.^6,20

The estimated washout time can be considered enough for the drug administered during period A not to affect the results of period B. Even though the use of remifentanil can modify depth of anesthesia, the present study focuses on analyzing the variation in entropy index and not on its mean value, thus, it is considered a constant administration for the whole intervention and in every patient to not affect the results of the study.

No information was found on the value of variance that could be considered clinically significant, and no statistically significant differences between the models were found for either of the entropy indices studied. Nevertheless, considering the range of appropriate anesthetic plane for this indices (40-60)¹² it can be reasoned that a value of variance under 100 could be tolerable because in spite of the fluctuations they would remain in an adequate range. Data analysis showed a power over 80% to detect differences of 100 units of variance between the models for State Entropy index, which means that even if there is a difference, the actual mean variability of depth of anesthesia between the two models is probably not enough to define a clear superiority of one model in the clinical setting.

Absalom et al., maintain in their article that propofol infusion controlled by the Marsh model increases usage of propofol when compared to Schnider model, authors remark that the important question is not, which model delivers the largest or smallest dose of drug, but which one produces the more accurate predictions of plasma and effect site concentration. Fig. 6 presents a 3-D visualization of the multi-compartment models for propofol in a patient (70 kg, 170 cm, 25 years; male) programmed in Matlab software. Each compartment is represented by spheres according to respective volumes, a wider link between compartment represents a higher velocity constant, and thus a faster circulation of pro-profol between compartments, blue links represent velocity parameters from central compartment V1 to peripheral compartments V2, V3, and the elimination rate; red links are the velocity parameters into V1.

The development of the pharmacokinetic models^3,4 did not include severely obese patients, for the case of a 96 kg, 1.70 m, 22 years, male person, the only difference in the Schnider model would be an increment in the elimination rate constant; the Marsh model would considers bigger compartments but the same elimination rate, and velocityconstants between compartments. It is important to keep testing determinants of the efficacy of the models and use the TCI pumps with caution in different types of population because their behavior according to individual characteristics such as the severely obese patient, or variables such as cost-effectiveness could end up tilting the scale.

Ethical disclosures

Protection of human and animal subjects. The authors declare that the procedures followed were in accordance with the regulations of the relevant clinical research ethics committee and with those of the Code of Ethics of the World Medical Association (Declaration of Helsinki).

Confidentiality of data. The authors declare that they have followed the protocols of their work center on the publication of patient data.

Right to privacy and informed consent. The authors have obtained the written informed consent of the patients or subjects mentioned in the article. The corresponding author is in possession of this document.

Funding

This study was funded by Colombian Society of Anesthesiol-ogy and Reanimation, in the contest modality of investigation proposals.

Conflicts of interest

The authors have no conflicts of interest to declare.

1. Cole CD, Gottfried ON, Gupta DK. Total intravenous anesthesia: advantages for intracraneal surgery. Neurosurgery. 2007;61:369-78.         [ Links ]

2. Hans P, Bonhomme V. Why we still use intravenous drugs as the basic regimen for neurosurgical anaesthesia. Curr Opin Anaesthesiol. 2006;19:498-503.         [ Links ]

3. Yasny JS, White J. Environmental implications of anesthetic gases. Anesth Prog. 2012;59:154-8.         [ Links ]

4. Lauder GR. Total intravenous anesthesia will supercede inhalational anesthesia in pediatric anesthetic practice. Pediatr Anesth. 2015;25:52-64.         [ Links ]

5. Tafur LA, Lema E. Anestesia total intravenosa. Rev Colomb Anestesiol. 2010;38:215-31.         [ Links ]

6. Absalom AR, Man V, De Smet T, Struys MMRF. Pharmacokinetic models for propofol - defining and illuminating the devil in the detail. Br J Anaesth. 2009;103:26-37.         [ Links ]

7. B.Braun. B.Braun Space Infusion Pumps [Internet]. Space System - Technical Data [cited 2015 Oct 1] . p. 7. Available from: http://www.space.bbraun.com/documents/Space_System_TechnicaLData_l.pdf.         [ Links ]

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

9. Schnider T, Minto C, Pedro L, Dds DB, Steven L, Elizabeth J. The influence of method of administration and covariates on the pharmacokinetics of propofol in adult volunteers methods study design. Anesthesiology. 1998;88:1170-82.         [ Links ]

10. Landers R, Wen P, Pather S. Depth of anaesthesia: measuring or guessing? IEEE International Conference on Nano/Molecular Medicine and Engineering, IEEE NANOMED. 2010. p. 76-81.         [ Links ]

11. Rampil IJ. A primer for EEG signal processing in anesthesia. Anesthesiology. 1998;89:980-1002.         [ Links ]

12. Viertio-Oja H, Maja V, Sárkelá M, Talja P, Tenkanen N, Tolvanen-Laakso H, et al. Description of the EntropyTM algorithm as applied in the Datex-Ohmeda S/5TM Entropy Module. Acta Anaesthesiol Scand. 2004;48:154-61.         [ Links ]

13. Schmidt GN, Bischoff P, Standl T, Hellstern A, Teuber O, Schulte Esch J. Comparative evaluation of the Datex-Ohmeda S/5 Entropy Module and the Bispectral Index monitor during propofol-remifentanil anesthesia. Anesthesiology. 2004;101:1283-90.         [ Links ]

14. Hasak L, Wujtewicz M, Owczuk R. Assessment of the depth of anaesthesia during inhalational and intravenous induction of general anaesthesia. Anesthesiol Intensive Ther. 2014;46:274-9.         [ Links ]

15. Wheeler P, Hoffman WE, Baughman VL, Koenig H. Response entropy increases during painful stimulation. J Neurosurg Anesthesiol. 2005;17:86-90.         [ Links ]

16. The MathWorks Inc. Statistics and Machine Learning Toolbox TM User's Guide r2015a;2015:7916 p.         [ Links ]

17. Wellek S, Blettner M. On the proper use of the crossover design in clinical trials. Dtsch Arztebl Int. 2012;109:276-81.         [ Links ]

18. Siyasinghe NM, Sooriyarachchi MR. Guidelines for calculating sample size in 2 x 2 crossover trials: a simulation study. J Natl Sci Found Sri Lanka. 2011;39:77-89.         [ Links ]

19. Walpole RE, Myers RH, Myers SL. In: Roig PE, Bonilla J, Ávalos Ó, editors. Probabilidad y estadística para ingenieros. 6th ed. Prentice-Hall Hispanoamericana SA; 1999. p. 752.         [ Links ]

20. Struys MM, De Smet T, Depoorter B, Vesichelen LF, Mortier EP, Dumortier FJ, et al. Comparison of plasma compartment versus two methods for effect compartment-controlled target-controlled infusion for propofol. Anesthesiology. 2000;92:399-406.         [ Links ]

21. Swinhoe CF, Peacock JE, Glen JB, Reilly CS. Evaluation of the predictive performance of a "Diprifusor" TCI system. Anaesthesia. 1998;53 Suppl. 1:61-7.         [ Links ]