Prediction of Morphine Clearance in the Paediatric Population: How Accurate are the Available Pharmacokinetic Models?

Scoping Review of Paediatric Population Pharmacokinetic Models of Morphine

OBJECTIVE AND PURPOSE

This scoping review aimed to summarise all available population pharmacokinetic models of morphine and its metabolites (morphine-3-glucoronide [M3G], morphine-6-glucoronide [M6G]) in children and describe how morphine exposure varies across paediatric age groups and settings. Identifying the factors that contribute to pharmacokinetic variability may improve our understanding of a patient's pharmacodynamic response to morphine.

METHODS

We searched Embase and MEDLINE databases from inception to 8 March 2024 for paediatric population pharmacokinetic models of morphine and its metabolites. Two reviewers independently screened abstracts and full texts and extracted the data. The review was conducted according to Preferred Reporting Items for Systematic Reviews and Meta-Analysis (PRISMA) guidelines.

RESULTS

In total, 21 paediatric population pharmacokinetic models of morphine were identified; 12 studies also included morphine metabolites (M3G and/or M6G). Neonates and young children (< 6 years) were the most studied age groups (18/21; 86%), whereas older children (> 6 years) and adolescents (> 10 years) were included in only 6 of the 21 (29%) models. Morphine pharmacokinetics were most commonly described with two-compartment (52%) and one-compartment (38%) structure with first-order elimination. Several model covariates were identified: bodyweight, post-natal age for neonates, body temperature, therapeutic cooling, duration of mechanical ventilation, and genetic variation in drug transporters that mediate the uptake of morphine (e.g. OCT1).

CONCLUSION

Several population pharmacokinetic models of morphine and its metabolites in paediatrics have been published across diverse patient groups. Bodyweight and age-related covariates emerged as the most common factors affecting clearance and distribution; other covariates, including mechanical ventilation, therapeutic cooling, and genetic variation, also impacted morphine pharmacokinetics. Further research should focus on validating the predictive accuracy of paediatric morphine models in different patient populations and the combined effect of covariates, such as those related to critical illness and genetic variation, on morphine pharmacokinetics.

Postoperative breakthrough pain in paediatric cardiac surgery not reduced by increased morphine concentrations

BACKGROUND

Morphine is commonly used for postoperative analgesia in children. Here we studied the pharmacodynamics of morphine in children after cardiac surgery receiving protocolized morphine.

METHODS

Data on morphine rescue requirements guided by validated pain scores in children (n = 35, 3-36 months) after cardiac surgery receiving morphine as loading dose (100 μg kg^-1) with continuous infusion (40 μg kg^-1 h^-1) from a previous study on morphine pharmacokinetics were analysed using repeated time-to-event (RTTE) modelling.

RESULTS

During the postoperative period (38 h (IQR 23-46)), 130 morphine rescue events (4 (IQR 1-5) per patient) mainly occurred in the first 24 h (107/130) at a median morphine concentration of 29.5 ng ml^-1 (range 7-180 ng ml^-1). In the RTTE model, the hazard of rescue morphine decreased over time (half-life 18 h; P < 0.001), while the hazard for rescue morphine (21.9% at 29.5 ng ml^-1) increased at higher morphine concentrations (P < 0.001).

CONCLUSIONS

In this study on protocolized morphine analgesia in children, rescue morphine was required at a wide range of morphine concentrations and further increase of the morphine concentration did not lead to a decrease in hazard. Future studies should focus on a multimodal approach using other opioids or other analgesics to treat breakthrough pain in children.

IMPACT

In children receiving continuous morphine infusion, administration of rescue morphine is an indicator for insufficient effect or an event. Morphine rescue events were identified at a wide range of morphine concentrations upon a standardized pain protocol consisting of continuous morphine infusion and morphine as rescue boluses. The expected number of rescue morphine events was found to increase at higher morphine concentrations. Instead of exploring more aggressive morphine dosing, future research should focus on a multimodal approach to treat breakthrough pain in children.

Modelling Tools to Characterize Acetaminophen Pharmacokinetics in the Pregnant Population

This review describes acetaminophen pharmacokinetics (PK) throughout pregnancy, as analyzed by three methods (non-compartmental analyses (NCA), population PK, and physiologically based PK (PBPK) modelling). Eighteen studies using NCA were reported in the scientific literature. These studies reported an increase in the volume of distribution (3.5-60.7%) and an increase in the clearance (36.8-84.4%) of acetaminophen in pregnant women compared to non-pregnant women. Only two studies using population PK modelling as a technique were available in the literature. The largest difference in acetaminophen clearance (203%) was observed in women at delivery compared to non-pregnant women. One study using the PBPK technique was found in the literature. This study focused on the formation of metabolites, and the toxic metabolite N-acetyl-p-benzoquinone imine was the highest in the first trimester, followed by the second and third trimester, compared with non-pregnant women. In conclusion, this review gave an overview on acetaminophen PK changes in pregnancy. Also, knowledge gaps, such as fetal and placenta PK parameters, have been identified, which should be explored further before dosing adjustments can be suggested on an evidence-based basis.

Model-Informed Bayesian Estimation Improves the Prediction of Morphine Exposure in Neonates and Infants

BACKGROUND

Pain control in infants is an important clinical concern, with potential long-term adverse neurodevelopmental effects. Intravenous morphine is routinely administered for postoperative pain management; however, its dose-concentration-response relationship in neonates and infants has not been well characterized. Although the current literature provides dosing guidelines for the average infant, it fails to control for the large unexplained variability in morphine clearance and response in individual patients. Bayesian estimation can be used to control for some of this variability. The authors aimed to evaluate morphine pharmacokinetics (PKs) and exposure in critically ill neonates and infants receiving standard-of-care morphine therapy and compare a population-based approach to the model-informed Bayesian techniques.

METHODS

The PKs and exposure of morphine and its active metabolites were evaluated in a prospective opportunistic PK study using 221 discarded blood samples from 57 critically ill neonates and infants in the neonatal intensive care unit. Thereafter, a population-based PK model was compared with a Bayesian adaptive control strategy to predict an individual's PK profile and morphine exposure over time.

RESULTS

Among the critically ill neonates and infants, morphine clearance showed substantial variability with a 40-fold range (ie, 2.2 to 87.1, mean 23.7 L/h/70 kg). Compared with the observed morphine concentrations, the population-model based predictions had an R of 0.13, whereas the model-based Bayesian predictions had an R of 0.61.

CONCLUSIONS

Model-informed Bayesian estimation is a better predictor of morphine exposure than PK models alone in critically ill neonates and infants. A large variability was also identified in morphine clearance. A further study is warranted to elucidate the predictive covariates and precision dosing strategies that use morphine concentration and pain scores as feedbacks.

Clinical pharmacology and dosing regimen optimization of neonatal opioid withdrawal syndrome treatments

In this paper, we review the management of neonatal opioid withdrawal syndrome (NOWS) and clinical pharmacology of primary treatment agents in NOWS, including morphine, methadone, buprenorphine, clonidine, and phenobarbital. Pharmacologic treatment strategies in NOWS have been mostly empirical, and heterogeneity in dosing regimens adds to the difficulty of extrapolating study results to broader patient populations. As population pharmacokinetics (PKs) of pharmacologic agents in NOWS become more well‐defined and knowledge of patient‐specific factors affecting treatment outcomes continue to accumulate, PK/pharmacodynamic modeling and simulation will be powerful tools to aid the design of optimal dosing regimens at the patient level. Although there is an increasing number of clinical trials on the comparative efficacy of treatment agents in NOWS, here, we also draw attention to the importance of optimizing the dosing regimen, which can be arguably equally important at identifying the optimal treatment agent.

Severity parameters for asphyxia or hypoxic-ischemic encephalopathy do not explain inter-individual variability in the pharmacokinetics of phenobarbital in newborns treated with therapeutic hypothermia

BACKGROUND

The current study uses a population modeling approach to evaluate and quantify the impact of severity of asphyxia and hypoxic-ischemic encephalopathy (HIE) on the pharmacokinetics of phenobarbital in asphyxiated newborns treated with therapeutic hypothermia.

METHODS

Included newborns received phenobarbital (the TOBY trial protocol). 120 plasma samples were available from 50 newborns, median (IQR) weight 3.3 (2.8-3.5) kg and gestational age 39 (39-40) weeks. NONMEM^® version 7.2 was used for the data analysis. Age, body weight, sex, concomitant medications, kidney and liver function markers, as well as severity parameters of asphyxia and HIE were tested as potential covariates of pharmacokinetics of phenobarbital. Severe asphyxia was defined as pH of arterial umbilical cord blood ≤7.1 and Apgar 5 ≤5, and severe HIE was defined as time to normalization of amplitude-integrated electroencephalography (aEEG) >24 h.

RESULTS

Weight was found to be the only statistically significant covariate for the volume of distribution. At weight of 1 kg volume of distribution was 0.91 L and for every additional kg it increased in 0.91 L. Clearance was 0.00563 L/h. No covariates were statistically significant for the clearance of phenobarbital.

CONCLUSIONS

Phenobarbital dose adjustments are not indicated in the studied population, irrespective of the severity of asphyxia or HIE.

Electronic Health Record-Embedded Decision Support Platform for Morphine Precision Dosing in Neonates

Morphine is the opioid most commonly used for neonatal pain management. In intravenous form, it is administered as continuous infusions and intermittent injections, mostly based on empirically established protocols. Inadequate pain control in neonates can cause long-term adverse consequences; however, providing appropriate individualized morphine dosing is particularly challenging due to the interplay of rapid natural physiological changes and multiple life-sustaining procedures in patients who cannot describe their symptoms. At most institutions, morphine dosing in neonates is largely carried out as an iterative process using a wide range of starting doses and then titrating to effect based on clinical response and side effects using pain scores and levels of sedation. Our background data show that neonates exhibit large variability in morphine clearance resulting in a wide range of exposures, which are poorly predicted by dose alone. Here, we describe the development and implementation of an electronic health record-integrated, model-informed decision support platform for the precision dosing of morphine in the management of neonatal pain. The platform supports pharmacokinetic model-informed dosing guidance and has functionality to incorporate real-time drug concentration information. The feedback is inserted directly into prescribers' workflows so that they can make data-informed decisions. The expected outcomes are better clinical efficacy and safety with fewer side effects in the neonatal population.

Influence of OCT1 Ontogeny and Genetic Variation on Morphine Disposition in Critically Ill Neonates: Lessons From PBPK Modeling and Clinical Study

Morphine is commonly used for analgesia in the neonatal intensive care unit (NICU) despite having highly variable pharmacokinetics (PKs) between individual patients. The pharmacogenetic (PG) effect of variants at the loci of organic cation transporter 1 (OCT1) and UDP-glucuronosyltransferase 2B7 (UGT2B7) on age-dependent morphine clearance were evaluated in a cohort of critically ill neonatal patients using an opportunistic sampling design. Our primary results demonstrate the significant influence of OCT1 genotype (P < 0.05) and gestational age (P ≤ 0.005) on morphine PKs. A physiologically based pharmacokinetic (PBPK) model for morphine that accounted for OCT1 ontogeny and PG effect in post-term neonates adequately described the clinically observed variability in morphine PKs. This study serves as a proof of concept for genotype-dependent drug transporter ontogeny in neonates.

PBPK Model of Morphine Incorporating Developmental Changes in Hepatic OCT1 and UGT2B7 Proteins to Explain the Variability in Clearances in Neonates and Small Infants

Morphine has large pharmacokinetic variability, which is further complicated by developmental changes in neonates and small infants. The impacts of organic cation transporter 1 (OCT1) genotype and changes in blood-flow on morphine clearance (CL) were previously demonstrated in children, whereas changes in UDP-glucuronosyltransferase 2B7 (UGT2B7) activity showed a small effect. This study, targeting neonates and small infants, was designed to assess the influence of developmental changes in OCT1 and UGT2B7 protein expression and modified blood-flow on morphine CL using physiologically based pharmacokinetic (PBPK) modeling. The implementation of these three age-dependent factors into the pediatric system platform resulted in reasonable prediction for an age-dependent increase in morphine CL in these populations. Sensitivity of morphine CL to changes in cardiac output increased with age up to 3 years, whereas sensitivity to changes in UGT2B7 activity decreased. This study suggests that morphine exhibits age-dependent extraction, likely due to the developmental increase in OCT1 and UGT2B7 protein expression/activity and hepatic blood-flow.

Update of the Scientific Opinion on opium alkaloids in poppy seeds

Poppy seeds are obtained from the opium poppy (Papaver somniferum L.). They are used as food and to produce edible oil. The opium poppy plant contains narcotic alkaloids such as morphine and codeine. Poppy seeds do not contain the opium alkaloids, but can become contaminated with alkaloids as a result of pest damage and during harvesting. The European Commission asked EFSA to provide an update of the Scientific Opinion on opium alkaloids in poppy seeds. The assessment is based on data on morphine, codeine, thebaine, oripavine, noscapine and papaverine in poppy seed samples. The CONTAM Panel confirms the acute reference dose (ARfD) of 10 μg morphine/kg body weight (bw) and concluded that the concentration of codeine in the poppy seed samples should be taken into account by converting codeine to morphine equivalents, using a factor of 0.2. The ARfD is therefore a group ARfD for morphine and codeine, expressed in morphine equivalents. Mean and high levels of dietary exposure to morphine equivalents from poppy seeds considered to have high levels of opium alkaloids (i.e. poppy seeds from varieties primarily grown for pharmaceutical use) exceed the ARfD in most age groups. For poppy seeds considered to have relatively low concentrations of opium alkaloids (i.e. primarily varieties for food use), some exceedance of the ARfD is also seen at high levels of dietary exposure in most surveys. For noscapine and papaverine, the available data do not allow making a hazard characterisation. However, comparison of the dietary exposure to the recommended therapeutical doses does not suggest a health concern for these alkaloids. For thebaine and oripavine, no risk characterisation was done due to insufficient data. However, for thebaine, limited evidence indicates a higher acute lethality than for morphine and the estimated exposure could present a health risk.

A review of perioperative anesthesia and analgesia for infants: updates and trends to watch

This review focuses on pharmacokinetics and pharmacodynamics of opioid and non-opioid analgesics in neonates and infants. The unique physiology of this population differs from that of adults and impacts drug handling. Morphine and remifentanil are described as examples of older versus recently developed opiates to compare and contrast pharmacokinetics and pharmacodynamics in infants. Exploration of genetics affecting both pharmacokinetics and pharmacodynamics of opiates is an area of active research, as is the investigation of a new class of mu-opiate-binding agents which seem selective for analgesic pathways while having less activity in pathways linked to side effects. The kinetics of acetaminophen and of ketorolac as examples of parenteral non-steroidal analgesics in infants are also discussed. The growth in regional anesthesia for peri-operative analgesia in infants can fill an important role minimizing intra-operative anesthetic exposure to opioids and transitioning to post-operative care. Use of multi-modal techniques is recommended to decrease undesirable opiate-related side effects in this vulnerable population.

Mechanistic Population Pharmacokinetics of Morphine in Neonates With Abstinence Syndrome After Oral Administration of Diluted Tincture of Opium

Conducting and analyzing clinical trials in vulnerable neonates are extremely challenging. The aim of this analysis is to develop a morphine population pharmacokinetics (PK) model using data collected during a randomized control trial in neonates with abstinence syndrome (NAS). A 3-compartment morphine structural PK model after intravenous (IV) administration from previously published work was utilized as prior, whereas an allometric scaling method with physiological consideration was used to extrapolate a PK profile from adults to pediatrics. The absorption rate constant and bioavailability were estimated in NAS after oral administration of diluted tincture of opium (DTO). Goodness-of-fit plots along with normalized prediction distribution error and bootstrap method were performed for model evaluation. We successfully extrapolated the PK profile from adults to pediatrics after IV administration. The estimated first-order absorption rate constant and bioavailability were 0.751 hour(-1) and 48.5%, respectively. Model evaluations showed that the model can accurately and precisely describe the observed data. The population pharmacokinetic model we derived for morphine after oral administration of DTO is reasonable and acceptable; therefore, it can be used to describe the PK and guide future studies. The integration of the previous population PK knowledge as prior information successfully overcomes the logistic and practical issue in vulnerable neonate population.

Developmental changes in morphine clearance across the entire paediatric age range are best described by a bodyweight-dependent exponent model

BACKGROUND AND OBJECTIVE

Morphine clearance has been successfully scaled from preterm neonates to 3-year-old children on the basis of a bodyweight-based exponential (BDE) function and age younger or older than 10 days. The aim of the current study was to characterize the developmental changes in morphine clearance across the entire paediatric age range.

METHODS

Morphine and morphine-3-glucuronide (M3G) concentration data from 358 (pre)term neonates, infants, children and adults, and morphine concentration data from 117 adolescents were analysed using NONMEM 7.2. Based on available data, two models were developed: I. using morphine data; II. using morphine and M3G data.

RESULTS

In model I, morphine clearance across the paediatric age range was very well described by a BDE function in which the allometric exponent decreased in a sigmoidal manner with bodyweight (BDE model) from 1.47 to 0.88, with half the decrease in exponent reached at 4.01 kg. In model II, the exponent for the formation and elimination clearance of M3G was found to decrease from 1.56 to 0.89 and from 1.06 to 0.61, with half the decrease reached at 3.89 and 4.87 kg, respectively. Using the BDE model, there was no need to use additional measures for size or age.

CONCLUSION

The BDE model was able to scale both total morphine clearance and glucuronidation clearance through the M3G pathway across all age ranges between (pre)term neonates and adults by allowing the allometric exponent to decrease across the paediatric age range from values higher than 1 for neonates to values lower than 1 for infants and children.