This article has Open Peer Review reports available.
Population pharmacokinetics of remifentanil in infants and children undergoing cardiac surgery
© Sam et al; licensee BioMed Central Ltd. 2009
Received: 16 September 2008
Accepted: 27 July 2009
Published: 27 July 2009
The aim of this study was to provide a model-based analysis of the pharmacokinetics of remifentanil in infants and children undergoing cardiac surgery with cardiopulmonary bypass (CPB).
We studied nine patients aged 0.5 to 4 years who received a continuous remifentanil infusion via a computer-controlled infusion pump during cardiac surgery with mildly hypothermic CPB were studied. Arterial blood samples taken prior to, during and after CPB were analyzed for remifentanil concentrations using a validated gas-chromatographic mass-spectrophotometric assay. We used population mixed-effects modeling to characterize remifentanil pharmacokinetics. The final model was evaluated by its predictive performance.
The pharmacokinetics of remifentanil was described by a 1-compartment model with adjustments for CPB. Population mean parameter estimates were 1.41 L for volume of distribution (V) and 0.244 L/min for clearance. V was increased during CPB and post-CPB to 2.41 times the pre-CPB value. The median prediction error and the median of individual median absolute prediction error were 2.44% and 21.6%, respectively.
Remifentanil dosage adjustments are required during and after CPB due to marked changes in the V of the drug. Simulations indicate that a targeted blood concentration of 14 ng/mL is achieved and maintained in 50% of typical patients by administration of an initial dose of 18 μg remifentanil followed by an infusion of 3.7 μg/min before, during and post-CPB, supplemented with a bolus dose of 25 μg given at the start of CPB.
Remifentanil is a selective μ-opioid receptor agonist that produces intense analgesia of rapid onset and ultra short duration . Accordingly, remifentanil provides the benefits of a high-dose opioid anesthetic during cardiac surgery (e.g. minimal sympathetic response to tracheal intubation and sternotomy) while allowing tracheal extubation in the operating room immediately following surgery.
CPB may alter the pharmacokinetics of anesthetic agents due to hypothermia, hemodilution, exclusion of the lungs from the circulation and a decrease in plasma protein concentration . However, the impact of CPB on the pharmacokinetics of remifentanil in infants and children has not been fully investigated. Davis et al  studied the pharmacokinetics of remifentanil before and after CPB in 12 pediatric patients undergoing repair of an atrial septal defect. They found that after CPB clearance values increased by 20% (from 38.7 ± 9.6 to 46.8 ± 14 mL/min/kg), without a meaningful change in the coefficient of variation. The authors concluded that the pharmacokinetics of remifentanil is predictable following CPB. However, this study failed to investigate the pharmacokinetics of a remifentanil infusion during CPB.
An unanticipated decrease in remifentanil plasma concentration during CPB may be detrimental. Significant intra-operative stress responses, postoperative complications and increased mortality have been reported in infants receiving inadequate analgesia during cardiac surgery [4, 5]. Therefore, the aim of our investigation was to determine the effect of CPB on the pharmacokinetics of remifentanil in infants and children during all phases of cardiac surgery using mixed-effect modeling in order to inform accurate dosage titration.
Nine infants and children scheduled for open-heart surgery requiring CBP were studied after Institutional Review Board approval and written informed parental consent were obtained. These patients were part of a prospective, randomized, controlled clinical trial to define the opioid analgesic requirement after a remifentanil-based anesthetic with or without spinal anesthetic blockade (SAB). Inclusion criteria were age 3 months to 6 years and planned tracheal extubation in the operating room after surgery (e.g. absence of severe pulmonary hypertension or heart failure). Exclusion criteria were contraindication to SAB and failure to obtain informed consent. Details of the study have been described elsewhere .
Premedication was given to patients over the age of 1 year (midazolam 0.5–0.75 mg/kg by mouth). Following placement of standard monitors, anesthesia was induced with sevoflurane and tracheal intubation was performed after administration of rocuronium. Anesthesia was maintained with isoflurane 0.3% and a continuous infusion of remifentanil with or without SAB with tetracaine (0.5 – 2.0 mg/kg) and morphine (0.007 mg/kg). Remifentanil was infused with a target-controlled infusion system programmed with parameters obtained from Minto et al  to maintain a constant predicted blood target concentration of 1–8 ng/mL. Remifentanil infusion rates were adjusted according to the judgment of the investigator in order to achieve a target mean arterial blood pressure (i.e. 50 to 60 mm Hg). Arterial blood pressure, heart rate, oxygen saturation and dosing history of remifentanil were continuously recorded during the surgery. At the completion of surgery, tracheal extubation was performed in the operating room and supplemental analgesia was achieved with fentanyl 0.003 mg/kg IV every 10 minutes as needed.
Prior to CPB, anticoagulation was established with an initial bovine heparin dose of 400 U/kg and additional heparin was administered during CPB to maintain celite activated clotting time (ACT) greater than 480 seconds. Non-pulsatile CPB was performed with a hollow fiber membrane oxygenator (Terumo CapioxC RX05, Terumo Cardiovascular Systems, Ann Arbor, MI), uncoated polyvinyl chloride bypass tubing and cannulae, and non-occlusive roller pump. The circuit was primed with normal saline, 25% albumin, mannitol, sodium bicarbonate, calcium chloride, methylprednisolone (30 mg/kg) and heparin. CPB circuit volumes were 450 mL for patients < 10 kg, 800 mL for patients 10–15 kg and 1,000–1,200 mL for patients > 15 kg. Banked packed red blood cells and fresh frozen plasma were added to achieve a hematocrit of about 30% during initiation of CPB. CPB flow rates were 200 mL/kg for infants with body weight less than 5 kilograms (kg), 150 mL/kg for those between 5 and 9 kg and 125 mL/kg for those between 10 and 17 kg; flows of 2.4 L/m2 were used in children over 17 kg. An initial dose of cardioplegia of 30 mL/kg was given, followed by 10 mL every 10–20 minutes thereafter. Hypothermia (28 – 32°C) was induced in all patients and blood gases were regulated according to alpha-stat regimen. Myocardial preservation was achieved using cold crystalloid cardioplegia. Target post-CPB hematocrit values varied from 35% to 50% depending upon the patient's cardiac and respiratory status. Antifibrinolytic agents were not administered.
Conventional ultra-filtration (CUF) was performed throughout CPB to achieve a filtrate volume of at least 120 mL/kg. Fluids (crystalloid, red blood cells or fresh frozen plasma) were added when necessary to provide sufficient volume in the CPB circuit to permit ultrafiltration. The polysulfone hemofilter used (MinntechC HPH 400; Minntech Corporation, Minneapolis, MN) employs hollow fiber technology and is rated to have a filtration cut-off to particles greater than 65,000 Daltons (Da). A transmembrane pressure gradient of at least 200 mm Hg was applied during ultrafiltration. After the addition of blood products, hemofiltration of the CPB circuit prime was performed before CPB to adjust pH and electrolyte concentrations and to remove inflammatory mediators. Filtrate volume from pre-CPB filtration ranged from 100 mL to 200 mL.
Arterio-venous modified ultra-filtration (MUF) was initiated on selected patients after separation from CPB according to the surgeons' preferences. Blood from the aortic cannula and from the CPB circuit venous reservoir was pumped through the hemofilter and was then warmed by a coiled heat exchanger (Medtronic MYOtherm XPR cardioplegia delivery system, Medtronic Inc., Minneapolis, MN) and returned via the cardioplegia circuit to the venous cannulae. Infusion rates were adjusted to maintain appropriate central venous and/or left atrial pressures. MUF was terminated when red cell salvage of circuit contents was judged by the perfusionist to be complete.
Whole blood arterial samples for assay of remifentanil concentrations were obtained before and 5 min after each adjustment of the infusion rate. No samples were taken prior to the initiation of remifentanil infusion. Samples were initially aspirated into heparinized syringes, immediately transferred to tubes containing 50% citric acid to inactivate plasma esterases and kept frozen at -20°C. Remifentanil concentrations were determined using a validated gas-chromatographic mass-spectrophotometric (GC-MS) assay with inter-, intraassay coefficients of variation and lower limit of quantification of 4.6%, 4.0% and 0.5 ng/mL, respectively .
where TV(P) refers to the typical value of the pharmacokinetic parameter for a patient with the reference covariate value and θcovariaterefers to the estimated fractional change in the typical value of the pharmacokinetic parameter for the investigated covariate.
where TVCL is the model predicted value for CL given the value for temperature. θCL represents the population central tendency for CL. θCL, TEMP represents a parameter quantifying the effect of temperature on CL, and TEMP is the body temperature in °C.
Using the basic pharmacokinetic model, each potential covariate was separately incorporated and tested for statistical significance by use of the NONMEM objective function and the standard errors of the parameters. The likelihood ratio test at the significance level α = 0.01 was used to discriminate between alternative hierarchical models. The α level of 0.01 corresponds to a reduction of 6.64 (χ2, p < 0.01; 1 degree of freedom) in the minimum objective function when 1 parameter is added to the model and was used to examine significance. In addition to the minimum objective function, diagnostic goodness-of-fit plots were used for model building and selection. If more than one significant covariate was found, the covariate model causing the largest decrease in objective function was chosen as a basis to explore the influence of additional covariates sequentially with the use of the same criteria.
where pi refers to the individual value of the respective pharmacokinetic parameter in the ith individual, θ is the typical value of the parameter, and exp(ρi) expresses the random difference between θ and pi. Values of ηi are assumed to be independently multi-variate and normally distributed, with mean zero and diagonal variance-covariance matrix Ω with diagonal elements (ω1 2,..., ωm 2).
where Cij is the jth remifentanil concentration of the ith individual predicted by the pharmacokinetic model, and Yij is the measured remifentanil concentration. The residual departure of Yij from Cij is represented by εij. Values of εij are assumed to be independently and normally distributed, with mean zero and variance σ2.
The likelihood ratio test at the significance level α = 0.01 was used to discriminate between hierarchical models. This corresponds to a decrease of ≥ 6.635 (one parameter difference) in the minimum objective function (-2 × logarithm of the likelihood of the results), as the difference in objective function between hierarchical models is approximately χ2distributed. When more than 1 parameter was added, a decrease in the objective function of ≥ 9.210, 11.345 and 13.277 was needed at the significance level of α = 0.01 for 2, 3 and 4 degrees of freedom, respectively.
where Ni is the number of samples obtained for the ith patient and j is the sample number for the ith patient. Model goodness-of-fit was checked by using diagnostic plots such as the PE over time plot and observed versus predicted plot. The final optimal model was selected based on the objective function, the predictive accuracy and goodness-of-fit plots.
The predictive performance of the population pharmacokinetic parameters of remifentanil by Minto et al  used to program STANPUMP was assessed by the prediction error (PE) and residual error plot. The intra-subject bias and precision of the predictions were assessed by quantifying the MDPE and MDAPE, respectively. The residual error plot is the percentage ratio of the measured concentration over the predicted concentration versus time.
We assessed the stability and robustness of the final models via bootstrapping . We obtained a new replication of original dataset (a bootstrap sample) by 1000 random draws of individual subject's data (with replacement) from the original dataset and we fitted the final model to each new dataset. We compared the final model parameter estimates to the mean of the parameter estimates obtained from 1000 bootstrap replicates of the final model.
Simulation was also performed using parameters of our final pharmacokinetic model to determine the optimal infusion and loading doses of remifentanil needed during the different intervals of cardiac surgery needed to maintain average concentration of 14 ng/mL. A typical patient of weight 11.4 kg and 60 min of surgery were assumed before and after 60 min of CPB for simulation purposes.
Atrial septal defect repair
Supraventricular aortic stenosis repair
Resection of pulmonary valve leaflets
Reconstruction of pulmonary valve
Reduction plasty of aneurysm of main pulmonary artery
Division of ductus arteriosus
Duration of CPB (min)*
75.0 ± 33.3
Population pharmacokinetic parameters and the stability of the parameters using the bootstrap resampling procedure.
1000 bootstrap replicates
Volume of distribution during pre-bypass, V (L)
Clearance, CL (L/min)
Effect of bypass and post-bypass on VPRE
Inter-individual variability ω CL (%)
Residual unexplained variability
Proportional residual error, coefficient of variation (%)
Pharmacokinetic parameter estimates of remifentanil in pediatric patients undergoing cardiac surgery.
Remifentanil is a widely used analgesic agent during cardiac surgery but little is known about its pharmacokinetics in infants and children undergoing CPB. We found that a 1-compartment pharmacokinetic model adjusted for CPB performed well. The relatively sparse sampling scheme we adopted in this study precludes the development of a multi-compartment model. Although a limited number of blood samples were collected during CPB, use of a V-adjusted CPB pharmacokinetic model produced a statistically significant improvement (p < 0.01) over the basic 1-compartment model.
The results of our study may be compared with those of other investigators studying the pharmacokinetics of remifentanil in pediatric and adult patients undergoing CPB [3, 12, 13] (Table 3). The population mean volume of distribution at steady-state (Vss) at pre-CPB, 124 mL/kg, estimated by our final CPB-adjusted pharmacokinetic model in this study is smaller than those previously reported in pediatric patients of 235 mL/kg [3, 12, 13]. However, the population mean Vss at CPB and post-CPB, 298 mL/kg is similar to the post-CPB value of 235 mL/kg reported by Davis et al . The relatively small population mean CL estimated in this study (21.4 mL/kg/min) maybe due to differences in patient characteristics from studies previously reported such as cardiac output, which has been reported to be a factor influencing the CL of remifentanil . The elimination half-lives of remifentanil were similar to those determined by previous investigators [3, 12, 13].
The results of our analysis indicate that V of remifentanil increases by 141% with the institution of CPB and remained increased during post-CPB. This is similar to the findings of Michelsen et al who found that the volume of distribution of remifentanil increased by 86% with institution of CPB, and remained increased during post-CPB  in adult patients. Hemodilution due to the CPB prime is the most likely reason for the increase in V. The 141% increase in the V of remifentanil in infants and children during and post- CPB compared to pre-CPB is higher than the expected range as the average CPB prime volume for patients between 10–15 kg is about 800 mL. The discrepancy may be due to drug adsorption to CPB circuit or additional sources of extracellular fluid secondary to rennin and atrial natriuretic peptide release as a result of non-pulsatile renal blood flow during CPB . Our results are also similar to the findings of Davis et al  who estimated similar volumes of distribution of remifentanil during post-CPB in children undergoing cardiac surgery.
The CL of remifentanil was found to be unaffected by CPB which is similar to that reported by Michelsen et al . Contrary to that reported by Michelsen et al  and Russell et al , temperature is not found to be a significant factor affecting the CL of remifentanil in this study. This may be due to a limited range of temperature (33.7–38.0°C) and mild hypothermia encountered in the majority of the patients during CPB. Moreover, the relatively short duration of CPB may have hindered the ability to identify changes in CL during bypass. The changes in the elimination half-lives (t1/2β) of remifentanil during the CPB and non-CPB follow the same trend and magnitude as the V (Table 3), suggesting that the pharmacokinetic advantage of short recovery from anesthesia post-CPB is modified in pediatric cardiac surgical patients.
MUF, first described by Naik et al , reduces total body water, induces hemoconcentration and improves postoperative recovery in children after open heart surgery [16, 17]. However, information is lacking regarding how MUF affects the pharmacokinetics of remifentanil in pediatric patients. Our study shows that MUF has little effect on the pharmacokinetics of remifentanil in pediatric cardiac patients.
The population model is able to predict blood concentrations accurately and with good precision as evidenced by a small MDPE and median individual MDAPE of 2.44% and 21.6%, respectively. These, together with the results of the diagnostic plots which showed random uniform scatter around the identity line (Figure 1) and an absence of bias, indicates the suitability of the final pharmacokinetic model as a basis for pharmacokinetic modeling. On the other hand, the poor predictive performance of the remifentanil target-infusion system incorporating the population pharmacokinetic parameters for remifentanil of Minto et al  developed using data from adult patients does not meet the criteria for clinical utility  and thus has a limited utility in pediatric patients.
One reason for pharmacokinetic analysis is to determine an effective dosing regimen. Therapeutic blood concentrations of remifentanil range from 1.5 ng/mL for sedation to 50 ng/mL for anesthetic effects in both humans and dogs [1, 12, 19, 20]. We selected an average concentration of 14 ng/mL is selected as the therapeutic concentration for the design of an effective dosing regimen for remifentanil. This concentration was chosen as a result of pharmacokinetic-pharmacodynamic analysis of the hypotensive effect of remifentanil in pediatric patients undergoing cranioplasty by the authors, which shows that this concentration leads to a 30% reduction in mean arterial pressure from the baseline. For a typical patient, an initial loading bolus dose of 18 μg remifentanil and a 3.7 μg/min continuous infusion before, during and post-CPB, supplemented with a bolus dose of 25 μg given at the start of CPB, achieves and maintains a targeted blood concentration of 14 ng/mL in 50% of the patients. The supplemental bolus dose is needed to compensate for the increase in V during CPB.
Because the population model was based on pediatric patients with mild hypothermia during CPB, this dosing recommendation may not apply to conditions of moderate to deep hypothermia during CPB as low body temperature has been reported to be a factor reducing CL during CPB [12, 13].
One limitation of our pharmacokinetic modeling is our model's failure to account for drug removed from the CPB circuit post-CPB. This problem can be solved by treating the CPB process as a separate compartment from the central compartment and the CPB compartment turn on/off when CPB starts/stops. Unfortunately, the relatively sparse sampling schedule in this study precludes using a multi-compartment model and the use of this approach to correct for the removal of drug from the system after CPB. Thus, we adopted an approach commonly reported in the literature of allowing the pharmacokinetic parameters to change during and after CPB for modeling.
We developed a CPB-adjusted pharmacokinetic model for remifentanil dosing in infants and children undergoing cardiac surgery with mild hypothermic CPB. Because the V of remifentanil was markedly increased during CPB and remained elevated post-CPB, a supplemental bolus dose of remifentanil is required during CPB in pediatric cardiac anesthesia. Studies are warranted to determine the pharmacokinetics of remifentanil in infants and children during moderate and deep hypothermia to determine how further temperature decreases affect pharmacokinetic parameters. Such future studies should lead to safer more consistently effective dosing of remifentanil in pediatric patients undergoing cardiac surgery with CPB.
Funding was partially paid by a grant from Physiometrix Inc, which is now owned by Hospira Inc. Also, partial funding was provided by a grant from the Department of Anesthesia at Stanford.
- Glass PSA, Hardman D, Kamiyama Y, Quill TJ, Marton G, Donn KH, Grosse CM, Hermann D: Preliminary pharmacokinetics and pharmacodynamics of an ultra-short-acting opioid – remifentanil (G187084B). Anesth Analg. 1993, 77: 1031-1040. 10.1213/00000539-199311000-00028.View ArticlePubMedGoogle Scholar
- Holley FO, Ponganis KV, Stanski DR: Effect of cardiopulmonary bypass on the pharmacokinetics of drugs. Clin Pharmacokinet. 1982, 7: 234-251. 10.2165/00003088-198207030-00004.View ArticlePubMedGoogle Scholar
- Davis PJ, Wilson AS, Siewers RD, Pigula FA, Landsman IS: The effects of cardiopulmonary bypass on remifentanil kinetics in children undergoing atrial septal defect repair. Anesth Analg. 1999, 89: 904-908. 10.1097/00000539-199910000-00016.PubMedGoogle Scholar
- Anand KJ, Hansen DD, Hickey PR: Hormonal-metabolic stress responses in neonates undergoing cardiac surgery. Anesthesiology. 1990, 73: 661-670. 10.1097/00000542-199010000-00012.View ArticlePubMedGoogle Scholar
- Anand KJ, Hickey PR: Halothane-morphine compared with high-dose sufentanil for anesthesia and postoperative analgesia in neonatal cardiac surgery. N Engl J Med. 1992, 326: 1-9.View ArticlePubMedGoogle Scholar
- Hammer GB, Ramamoorthy C, Cao H, Williams GD, Boltz MG, Kamra K, Drover DR: Postoperative analgesia after spinal blockade in infants and children undergoing cardiac surgery. Anesth Analg. 2005, 100: 1283-1288. 10.1213/01.ANE.0000148698.84881.10.View ArticlePubMedGoogle Scholar
- Minto CF, Schnider TW, Egan TD, Youngs E, Lemmens HJ, Gambus PL, Billard V, Hoke JF, Moore KH, Hermann DJ, et al: Influence of age and gender on the pharmacokinetics and pharmacodynamics of remifentanil: I. Model development. Anesthesiology. 1997, 86: 10-23. 10.1097/00000542-199701000-00004.View ArticlePubMedGoogle Scholar
- Beal SL, Sheiner LB, Boeckmann AJ, (eds): NONMEM user's guide. 1989, Ellicott City, MD: Icon Development SolutionsGoogle Scholar
- Anderson BJ, Meakin GH: Scaling for size: some implications for paediatric anaesthesia dosing. Paediatr Anaesth. 2002, 12: 205-219. 10.1046/j.1460-9592.2002.00616.x.View ArticlePubMedGoogle Scholar
- 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-122. 10.1097/00000542-199401000-00018.View ArticlePubMedGoogle Scholar
- Ette EI: Stability and performance of a population pharmacokinetic model. J Clin Pharmacol. 1997, 37: 486-495.View ArticlePubMedGoogle Scholar
- Michelsen LG, Holford NHG, Lu W, Hoke JF, Hug CC, Bailey JM: The pharmacokinetics of remifentanil in patients undergoing coronary artery bypass grafting with cardiopulmonary bypass. Anesth Analg. 2001, 93: 1100-1105. 10.1097/00000539-200111000-00006.View ArticlePubMedGoogle Scholar
- Russell D, Royston D, Rees PH, Gupta SK, Kenny GNC: Effect of temperature and cardiopulmonary bypass on the pharmacokinetics of remifentanil. Br J Anaesth. 1997, 79: 456-459.View ArticlePubMedGoogle Scholar
- Duthie DJ, Stevens JJ, Doyle AR, Baddoo HH, Gupta SK, Muir KT, Kirkham AJ: Remifentanil and pulmonary extraction during and after cardiac anesthesia. Anesth Analg. 1997, 84: 740-744. 10.1097/00000539-199704000-00007.PubMedGoogle Scholar
- Bartlett RH: Extracorporeal life support for cardiopulmonary failure. Curr Probl Surg. 1990, 27: 621-705. 10.1016/0011-3840(90)90015-W.View ArticlePubMedGoogle Scholar
- Naik SK, Knight A, Elliot MJ: A prospective randomized study of modified technique of ultrafiltration during pediatric open-heart surgery. Circulation. 1991, 84: 422-431.Google Scholar
- Journois D, Israel-Biet D, Pouard P, Rolland B, Silvester W, Vouhe P, Safran D: High-volume, zero-balanced hemofiltration to reduce delayed inflammatory response to cardiopulmonary bypass in children. Anesthesiology. 1996, 85: 965-976. 10.1097/00000542-199611000-00003.View ArticlePubMedGoogle Scholar
- Glass PS, Jacobs JR, Smith LR, Ginsberg B, Quill TJ, Bai SA, Reves JG: Pharmacokinetic model-driven infusion of fentanyl: assessment of accuracy. Anesthesiology. 1990, 73: 1082-1090.View ArticlePubMedGoogle Scholar
- Lang E, Kapila A, Shlugman D, Hoke JF, Sebel PS, Glass PS: Reduction of isoflurane minimal alveolar concentration by remifentanil. Anesthesiology. 1996, 85: 721-728. 10.1097/00000542-199610000-00006.View ArticlePubMedGoogle Scholar
- Randel GI, Fragen RJ, Librojo ES, Jamerson BD, Gupta S: Remifentanil blood concentration effect relationship at intubation and skin incision in surgical patients compared to alfentanil. Anesthesiology. 1994, 81: A375-10.1097/00000542-199409001-00375.View ArticleGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2253/9/5/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.