Pharmacokinetics of morphine in encephalopathic neonates treated with therapeutic hypothermia

Abstract

Objective

Morphine is a commonly used drug in encephalopathic neonates treated with therapeutic hypothermia after perinatal asphyxia. Pharmacokinetics and optimal dosing of morphine in this population are largely unknown. The objective of this study was to describe pharmacokinetics of morphine and its metabolites morphine-3-glucuronide and morphine-6-glucuronide in encephalopathic neonates treated with therapeutic hypothermia and to develop pharmacokinetics based dosing guidelines for this population.

Study design

Term and near-term encephalopathic neonates treated with therapeutic hypothermia and receiving morphine were included in two multicenter cohort studies between 2008–2010 (SHIVER) and 2010–2014 (PharmaCool). Data were collected during hypothermia and rewarming, including blood samples for quantification of morphine and its metabolites. Parental informed consent was obtained for all participants.

Results

244 patients (GA mean (sd) 39.8 (1.6) weeks, BW mean (sd) 3,428 (613) g, male 61.5%) were included. Morphine clearance was reduced under hypothermia (33.5°C) by 6.89%/°C (95% CI 5.37%/°C– 8.41%/°C, p<0.001) and metabolite clearance by 4.91%/°C (95% CI 3.53%/°C– 6.22%/°C, p<0.001) compared to normothermia (36.5°C). Simulations showed that a loading dose of 50 μg/kg followed by continuous infusion of 5 μg/kg/h resulted in morphine plasma concentrations in the desired range (between 10 and 40 μg/L) during hypothermia.

Conclusions

Clearance of morphine and its metabolites in neonates is affected by therapeutic hypothermia. The regimen suggested by the simulations will be sufficient in the majority of patients. However, due to the large interpatient variability a higher dose might be necessary in individual patients to achieve the desired effect.

Trial registration

www.trialregister.nl NTR2529.

About the Authors:

Laurent M. A. Favié

Roles Data curation, Formal analysis, Project administration, Resources, Software, Validation, Visualization, Writing – original draft, Writing – review & editing

* E-mail: [email protected]

Affiliations Department of Clinical Pharmacy, University Medical Center Utrecht, Utrecht, the Netherlands, Department of Neonatology, Wilhelmina Children's Hospital, University Medical Center Utrecht, Utrecht, the Netherlands

ORCID logo http://orcid.org/0000-0002-4894-5236

Floris Groenendaal

Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Resources, Supervision, Writing – review & editing

Affiliations Department of Neonatology, Wilhelmina Children's Hospital, University Medical Center Utrecht, Utrecht, the Netherlands, Brain Center Rudolf Magnus, University Medical Center Utrecht, Utrecht, the Netherlands

Marcel P. H. van den Broek

Roles Conceptualization, Funding acquisition, Investigation, Project administration, Resources, Writing – review & editing

Affiliations Department of Clinical Pharmacy, University Medical Center Utrecht, Utrecht, the Netherlands, Department of Clinical Pharmacy, St. Antonius Hospital, Nieuwegein, the Netherlands

Carin M. A. Rademaker

Roles Conceptualization, Funding acquisition, Methodology, Resources, Supervision, Writing – review & editing

Affiliation: Department of Clinical Pharmacy, University Medical Center Utrecht, Utrecht, the Netherlands

Timo R. de Haan

Roles Conceptualization, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Writing – review & editing

Affiliation: Department of Neonatology, Emma Children’s Hospital, Academic Medical Center, Amsterdam, the Netherlands

Henrica L. M. van Straaten

Roles Conceptualization, Investigation, Resources, Supervision, Writing – review & editing

Affiliation: Department of Neonatology, Isala Clinics, Zwolle, the Netherlands

Peter H. Dijk

Roles Investigation, Resources, Supervision, Writing – review & editing

Affiliation: Department of Neonatology, Groningen University Medical Centre, Groningen, the Netherlands

Arno van Heijst

Roles Investigation, Resources, Supervision, Writing – review & editing

Affiliation: Department of Neonatology, Radboud university medical center-Amalia Children’s Hospital, Nijmegen, the Netherlands

Jeroen Dudink

Roles Investigation, Resources, Supervision, Writing – review & editing

Affiliation: Department of Pediatrics, Division of Neonatology, Erasmus Medical Centre-Sophia Children’s Hospital, Rotterdam, the Netherlands

Koen P. Dijkman

Roles Investigation, Resources, Supervision, Writing – review & editing

Affiliation: Department of Neonatology, Máxima Medical Center Veldhoven, Veldhoven, the Netherlands

Monique Rijken

Roles Investigation, Resources, Supervision, Writing – review & editing

Affiliation: Department of Neonatology, Leiden University Medical Center, Leiden, the Netherlands

Inge A. Zonnenberg

Roles Investigation, Resources, Supervision, Writing – review & editing

Affiliation: Department of Neonatology, VU University Medical Center, Amsterdam, the Netherlands

Filip Cools

Roles Investigation, Resources, Supervision, Writing – review & editing

Affiliation: Department of Neonatology, UZ Brussel—Vrije Universiteit Brussel, Brussels, Belgium

Alexandra Zecic

Roles Investigation, Resources, Supervision, Writing – review & editing

Affiliation: Department of Neonatology, University Hospital Gent, Gent, Belgium

Johanna H. van der Lee

Roles Conceptualization, Funding acquisition, Methodology, Writing – review & editing

Affiliation: Paediatric Clinical Research Office, Emma Children’s Hospital, Academic Medical Center, University of Amsterdam, Amsterdam, The Netherlands

Debbie H. G. M. Nuytemans

Roles Project administration, Resources, Software, Writing – review & editing

Affiliation: Clinical Research Coordinator PharmaCool Study, Amsterdam, the Netherlands

Frank van Bel

Roles Conceptualization, Investigation, Resources, Supervision, Writing – review & editing

Affiliations Department of Neonatology, Wilhelmina Children's Hospital, University Medical Center Utrecht, Utrecht, the Netherlands, Brain Center Rudolf Magnus, University Medical Center Utrecht, Utrecht, the Netherlands

Toine C. G. Egberts

Roles Conceptualization, Funding acquisition, Resources, Supervision, Writing – review & editing

Affiliations Department of Clinical Pharmacy, University Medical Center Utrecht, Utrecht, the Netherlands, Department of Pharmacoepidemiology and Clinical Pharmacology, Faculty of Science, Utrecht University, Utrecht, the Netherlands

Alwin D. R. Huitema

Roles Data curation, Formal analysis, Investigation, Resources, Software, Supervision, Validation, Visualization, Writing – review & editing

Affiliations Department of Clinical Pharmacy, University Medical Center Utrecht, Utrecht, the Netherlands, Department of Pharmacy & Pharmacology, Netherlands Cancer Institute, Amsterdam, the Netherlands

on behalf of the PharmaCool study group

¶A complete list of non-author contributors appears in the Acknowledgments section.

Introduction

Hypoxic-ischemic encephalopathy (HIE) following perinatal asphyxia is one of the leading causes of death and disability in term and near term neonates. Therapeutic hypothermia (TH, lowering the core temperature to 33–34° C for 72h) is an established neuroprotective strategy and has become standard of care for these patients in developed countries.[1,2] In the Netherlands, approximately 150–200 neonates receive this treatment annually using whole-body cooling.[3]

Morphine is a commonly used drug in hypothermic neonates to provide analgesia and sedation, and is considered an important drug since stress may reduce the neuroprotective effects of TH.[4] Morphine undergoes extensive hepatic metabolism and its predominant metabolite is morphine-3-glucuronide (M3G) which is non-sedative. The less abundant metabolite morphine-6-glucuronide (M6G) is pharmacologically active with similar or greater sedative and analgesic effects compared to the parent compound.[5,6] Both glucuronide metabolites are formed by UDP glucuronosyltransferase 2B7 (UGT2B7).[5] The UGT2B7 enzyme activity in neonates is less than 10% of that in adults, but increases rapidly during the first days after birth.[7,8] Both metabolites are eliminated through the kidneys.[5] At birth, renal function is underdeveloped compared to older children and adults. In the first few weeks of life, a steady increase in renal function can be seen.[8] Thus, maturation of kidney function might influence metabolite elimination.[8,9]

Hypothermia might influence numerous physiological processes involved in drug metabolism. Hypothermia reduces cardiac output and increases vascular resistance, which leads to decreased liver perfusion. Decreased liver perfusion might result in decreased drug clearance, especially in drugs with a high hepatic extraction ratio. Furthermore, the activity of liver enzymes such as cytochrome P450 and UGT27B can be declined during TH resulting in impaired clearance. Likewise, TH can decrease renal drug clearance by reducing kidney perfusion and subsequent glomerular filtration but also through changes in tubular secretion and reabsorption.[10–12] Additionally, pharmacokinetics (PK) of drugs administered to these neonates may be altered due to hypoxia-ischemia related multi-organ failure.[13,14] In recent years, studies have been conducted that investigated the PK of frequently used drugs in neonates undergoing TH. The findings of these studies have led to dose recommendations for several antibiotics and anticonvulsive drugs.[15–21]

Morphine PK in normothermic neonates has been investigated in several studies, mostly involving preterm and term neonates following major thoracic and abdominal surgery.[22] Neonates with a postnatal age (PNA) below 10 days had a markedly reduced morphine clearance compared to older children which has been attributed to impaired glucuronidation. This effect was independent of birth weight (BW) or gestational age (GA). Maintenance dose in this group was reduced by 50% compared to older children to achieve morphine plasma concentrations between 10 and 40 μg/L.[9] This dosing algorithm has been prospectively validated and body weight has been shown to accurately predict morphine clearance across the entire pediatric population.[23–25]

Morphine PK in neonates with HIE undergoing TH has only sparsely been investigated. Róka et al. (2008) found elevated morphine plasma concentrations in neonates treated with TH (N = 10) compared to non-asphyxiated normothermic controls (N = 6) with similar infusion rates and cumulative doses.[26] Frymoyer et al. (2016) developed a population PK model for morphine, M3G and M6G during TH using data from 20 neonates. They concluded that morphine clearance during TH was lower compared to previous studies in normothermic asphyxiated neonates and advised a loading dose of 50 μg/kg followed by 5 μg/kg continuous infusion.[27] Both studies did not include data during and after rewarming. Additional characterization of morphine PK using a larger dataset is imperative to guide clinicians in the application of this widely used and important drug in this critically ill population.

The objective of the present study was to describe the PK of morphine and its metabolites in neonates with HIE both during and after TH using nonlinear mixed effect modelling and to develop pharmacokinetics based dosing guidelines based on a large dataset obtained from two multicenter studies conducted in the Netherlands and Belgium.

Patients and methods

Setting, study design and study population

The open label prospective SHIVER study was performed in the tertiary neonatal intensive care units (NICU) of the University Medical Center Utrecht, Utrecht and Isala Clinics, Zwolle. The open label prospective PharmaCool study was conducted in twelve tertiary NICUs in the Netherlands and Belgium.[28] In both studies, term neonates undergoing TH for HIE were eligible for inclusion. According to national treatment protocol, neonates with a GA between 36.0 and 42.0 weeks were cooled within 6 hours after birth to a core temperature of 33.5°C (accepted range 33.0–34.0°C) for 72 hours. Thereafter, patients were slowly (0.4°C/hour) rewarmed to normothermia (36.5°C). After rewarming, body temperature was stabilized at 36.5°C for 24 hours.[3] Exclusion criteria were severe congenital malformations, encephalopathy due to other causes than perinatal asphyxia and the absence of central venous or arterial access for non-invasive blood sampling. From each included patient, written parental informed consent was obtained. Inclusion took place between 2008–2010 (SHIVER) and 2010–2014 (PharmaCool). In total, 339 patients were screened and 277 included. For the present study, neonates participating in either study and receiving intravenous morphine were selected. Data analysis was completed in 2018. The SHIVER study was approved by the Institutional Review Board (IRB) of the University Medical Center Utrecht (no. 08/404) and subsequently approved by the IRB of the Isala Clinics, Zwolle. The PharmaCool study was approved by the IRB of the Academic Medical Center Amsterdam (no. 10/255) and subsequently approved by the IRBs of the VU Medical Center Amsterdam, University Medical Center Utrecht, Leiden University Medical Center, Erasmus Medical Center Rotterdam, Maxima Medical Center Veldhoven, Maastricht University Medical Center, Radboud University Medical Center Nijmegen, Isala Clincs Zwolle, University Medical Center Groningen, University Hospital Gent and University Hospital Brussels.

Morphine dosing and administration

In both studies, morphine was administrated as morphine hydrochloride according to local protocols and/or the attending physician’s discretion as an intravenous continuous infusion, often preceded by a loading dose. Morphine was generally started at the onset of TH or shortly before. Dose adjustments, including administration of any additional loading dose, were based on each patient’s clinical condition and were not influenced by the study protocol.

Pharmacokinetic sampling and bioanalysis

From all patients, 1 ml blood samples were obtained from an indwelling catheter on four consecutive days, both during hypothermia and rewarming/normothermia. Sampling was scheduled at designated time points at 24 hours intervals. This limited sampling strategy was designed to minimize patient risk while still obtaining sufficient information to achieve the study objective. In the SHIVER study, residual material from blood samples taken for clinical care were available for some patients. Plasma concentrations of morphine, M3G and M6G were determined using liquid chromatography-tandem mass spectrometry (LC-MS/MS). The lower limit of quantification (LLQ) was 10 μg/L for morphine and M3G and 5 μg/L for M6G. The calibration curves were linear from 10 to 1000 μg/L for morphine, 10 to 600 μg/L for M3G and 5 to 200 μg/L for M6G. Between-run and within-run coefficients of variation were <5% for morphine and M3G and <8% for M6G. Samples were stored at -80°C until analyses at the Clinical Pharmaceutical and Toxicological Laboratory of the Department of Clinical Pharmacy of the University Medical Center Utrecht, the Netherlands.

Population pharmacokinetic analysis

A population pharmacokinetic model was developed from morphine, M3G, and M6G concentration–time data using the nonlinear mixed effect modelling program NONMEM (version 7.3, Icon Development Solutions) with R (version 3.4.1), Xpose (version 4) for data visualization and Piraña for run management.[29] Morphine hydrochloride (molecular weight (MW) 321.8 g/mol) doses were converted to morphine base (MW 285.3 g/mol) and consecutively, all units of dose and concentration for morphine, M3G and M6G (MW 461.5 g/mol) were converted to μmol and μmol/L, respectively for the purpose of the pharmacokinetic analysis. BW was used as a descriptor for body size in our population and was related to pharmacokinetic parameters using allometric relationships. The exponent defining the relationship of BW and clearance (Cl) was fixed to 0.75 and the exponent defining the relationship of BW and volume of distribution (V) was fixed to 1. The fractions of morphine converted to the metabolite M3G and M6G in neonates under hypothermia were unknown (FM3G and FM6G, respectively); therefore, parameters relative to F were estimated (e.g. ClM3G/FM3G and VM3G/FM3G). Based on previously published pharmacokinetic models of morphine in neonates[27,30], one- and two-compartment models for morphine and subsequent one-compartment models for both metabolites were tested as a structural model. Morphine and metabolite data were fitted simultaneously.

To study the effect of hypothermia on pharmacokinetics, a dynamic model of temperature over time was included, which allowed prediction of the actual body temperature at each moment of sampling. For all patients, the reported start and end times of TH were used to determine the period of TH treatment. Body temperature during hypothermia was set at 33.5°C, with consecutive rewarming at 0.4°C/hour (i.e. rewarming time 7.5h) until 36.5°C after which body temperature was set to 36.5°C for the remainder of the study time. Calculated body temperature for each plasma sample was subsequently included in the PK model.

As renal function may be an important determinant for metabolite clearance and given that renal function cannot be estimated from a single serum creatinine (SCr) measurement in neonates, a model for SCr was developed using daily SCr values taken for clinical care from all patients. In this model, the elimination rate of SCr was used as a surrogate marker for renal function. PNA and GA were tested as covariates on both morphine and metabolite clearance. Inclusion of covariates was guided by effect size, biological plausibility and statistical significance (using the likelihood ratio test which assesses the difference in the NONMEM objective function value (OFV), which is equal to minus twice the log likelihood, with a p-value of <0.05 as cut-off for significance).

Interindividual variability (IIV) was modelled using a proportional model and tested on all parameters. Covariance between IIV components was included based on physiological plausibility and graphical exploration. A proportional error model was used to model residual unexplained variability. For each compound, separate error models were used. Parameter precision was assessed with sampling importance resampling (SIR).[31] Internal validation of the final model was evaluated by computing the normalized prediction distribution errors (NPDE, 1000 simulations).[32] Both graphical (e.g. goodness-of-fit plots, visual predictive check) and statistical model evaluation procedures were used to assess model adequacy.

Dosing regimen development

Simulations were conducted to test four different dosing regimens using the parameter estimates from the final pharmacokinetic model. To create the simulation dataset, the patient characteristics of each neonate included in this study were replicated five times. Morphine loading dose was simulated at PNA 4 hours, immediately followed by continuous infusion. The following dosing regimens were evaluated, based on the current clinical practice: 1. loading dose of 50 μg/kg followed by continuous infusion of 5 μg/kg/h; 2. loading dose of 50 μg/kg, continuous infusion of 10 μg/kg/h; 3. loading dose of 100 μg/kg, continuous infusion of 5 μg/kg/h; 4. loading dose of 100 μg/kg, continuous infusion of 10 μg/kg/h. The dynamic temperature model was used to introduce TH. For each neonate in the simulation dataset, TH (body temperature of 33.5°C for 72 hours) was simulated to start at PNA 5 hours, after which rewarming commenced at 0.4°C/hour. After rewarming, body temperature was fixed to 36.5°C for the remainder of the simulations. Hourly plasma concentrations were predicted until PNA 120 hours. Morphine plasma concentrations between 10 and 40 μg/L were considered effective and safe.

Results

Patient characteristics

For 244 neonates morphine dosing information and at least one morphine plasma concentration was available for analysis (Table 1). In general, loading doses between 50 and 100 μg/kg were given, followed by continuous infusion with doses varying between 5 and 25 μg/kg/h.

[Figure omitted. See PDF.]

Table 1. Patient characteristics.

https://doi.org/10.1371/journal.pone.0211910.t001

A total of 853 blood samples were analyzed (median 4 samples per patient, range 1–11). Samples with measurements below LLQ for all compounds (n = 23) were excluded from further analyses, leaving 830 viable samples. Of these, 550 (66.3%) were drawn during the hypothermic phase. For 18 patients (7.4%), only one sample was available. Plasma concentrations for morphine varied between 10.0 and 371.2 μg/L (Fig 1); for M3G between 11.0 and 930.6 μg/L and for M6G between 5.1 and 211.2 μg/L (S1 Fig). In one patient (0.41%), morphine plasma concentrations exceeding 300 μg/L were reached.

[Figure omitted. See PDF.]

Fig 1. Observed morphine plasma concentrations (μg/L).

Dotted lines indicate the proposed therapeutic window of 10–40 μg/L; solid line indicates the potentially toxic limit of 300 μg/L.

https://doi.org/10.1371/journal.pone.0211910.g001

Population pharmacokinetic analysis

A one-compartment model for morphine and subsequent one-compartment models for both metabolites provided the best fit for the data. Pharmacokinetic parameter estimates of the final model are shown in Table 2.

[Figure omitted. See PDF.]

Table 2. Final model pharmacokinetic parameter estimates and SIR results.

https://doi.org/10.1371/journal.pone.0211910.t002

Introduction of a peripheral compartment for morphine resulted in an unstable model with unrealistic intercompartmental clearance. GA and PNA were identified as covariates on morphine clearance (GA: p<0.001, PNA: p<0.001), but not on metabolite clearance. Morphine clearance was increased by 50.4% at PNA 5 days, compared to birth (increase of 0.42%/h, 95%CI 0.297%/h– 0.582%/h); at birth, morphine clearance in a neonate with GA 36 weeks was 46% lower compared to GA 40 weeks, while clearance in a neonate with GA 42 weeks is 23% higher (difference of 1.66%/d, 95%CI 1.30%– 1.94%). The elimination rate of SCr was introduced as a covariate on the clearance of the metabolites as a measure of renal function. The influence of this covariate was non-significant and therefore excluded from the final model.

Subsequently, the dynamic model of temperature over time was included as covariate on Cl. The influence of body temperature on clearance was separated into an effect on ClMORPHINE (a combination of hepatic and renal clearance) and CLMETABOLITES (renal clearance). As the effect of body temperature on M3G and M6G clearance were similar and separate effects for both metabolites did not improve model performance, this was estimated as a single effect in the final model. Morphine clearance during hypothermia was decreased by 20.7% (p<0.001, 6.89%/°C, 95% CI 5.37%/°C– 8.41%/°C) compared to normothermia. Metabolite clearance during hypothermia was decreased by 14.7% (p<0.001, 4.91%/°C, 95% CI 3.53%/°C– 6.22%/°C).

The influence of BW, GA, PNA and temperature on the average morphine clearance is predicted by the final model are depicted in Fig 2.

[Figure omitted. See PDF.]

Fig 2.

Average predicted morphine clearance over time before, during and after TH for neonates with BW 3.5 and GA 36, 38, 40 and 42 weeks, respectively (left) and for neonates with GA 40 weeks and BW 2.5, 3.0, 3.5 and 4.0 kg, respectively (right). Solid vertical lines represent the start and end of TH (33.5°C) simulated between 5h and 77h after birth; dashed vertical line indicates the return to normothermia (36.5°C) with rewarming simulated at. 0.4°C/h; TH = therapeutic hypothermia, BW = birth weight GA = gestational age.

https://doi.org/10.1371/journal.pone.0211910.g002

Model evaluation demonstrated that the final model was adequate in describing the data. Goodness-of-fit plots of observed versus population and individual predicted concentrations showed no systematic deviation and the weighted residuals were homogeneously scattered for both parent and metabolites (S2–S4 Figs). NPDE plots for morphine, M3G and M6G indicate that the NPDE follows the normal distribution and that the model does not contain major bias (S5–S7 Figs).

Dosing regimen

Morphine plasma concentrations after various dosing regimens were predicted using a simulation dataset of 1220 patients and the final PK parameter estimates. In all simulations, morphine clearance was markedly influenced by PNA and TH. Immediately after rewarming, average morphine clearance was increased by 63.4% compared to clearance at the start of TH. Of this increase, 29.6% could be attributed to an effect of PNA. A maintenance dose of 5 μg/kg/h preceded by a loading dose of 50 μg/kg resulted in plasma concentrations between 10 and 40 μg/L at PNA 12h in 88.2% of patients, while 7.8% of patients were below 10 μg/L and 4.0% above 40 μg/L. At PNA 48h, morphine plasma concentration exceeded 40 μg/L in 6.8% of patients. As clearance is not constant but increased over time, no steady state in morphine plasma concentration was reached in the first five days of life. At PNA 77 hours, TH was stopped resulting in an additional increase in clearance and drop in plasma concentration (Fig 3). Plasma concentrations for both metabolites accumulated during TH but reached steady state once clearance increased under normothermic conditions (S8 Fig). Morphine plasma concentrations for the other simulated dosing regimens are included as a supplement (S9 Fig).

[Figure omitted. See PDF.]

Fig 3. Simulated morphine plasma concentrations of the proposed dosing regimen of 5 μg/kg/h after loading dose of 50 μg/kg.

Solid line indicates the mean morphine plasma concentration; gray area represents the 95% prediction interval. Dotted horizontal lines indicate the proposed therapeutic window of 10–40 μg/L. Solid vertical lines indicate the start and end of TH (33.5°C) simulated between 5h and 77h after birth; dashed vertical line indicates the return to normothermia (36.5°C) with rewarming simulated at 0.4°C/h; TH = therapeutic hypothermia.

https://doi.org/10.1371/journal.pone.0211910.g003

Discussion

This study combining data from two large multicenter studies shows that clearance of both morphine and its metabolites is reduced during hypothermia in neonates with HIE compared to normothermia. Furthermore, the impact of BW, GA and PNA on the PK of morphine and its metabolites has been quantified. Reduction in clearance during TH is most likely caused by a decrease in perfusion of the liver and kidneys. Additionally, an effect of TH on activity of UGT2B7, the enzyme responsible for metabolizing morphine into M3G and M6G, might explain why morphine clearance is more strongly affected than metabolite clearance.

Although a therapeutic window for morphine plasma concentrations has not been firmly established, especially in neonates undergoing TH, the best available evidence suggests a preferred range between 10 and 40 μg/L, while levels above 300 μg/L have been associated with respiratory depression and prolonged mechanical ventilation.[25,33–38] Based on the simulations performed in this study, a starting dose of 50 μg/kg followed by 5 μg/kg/h is recommended to achieve morphine plasma concentrations between 10 and 40 μg/L, although the large interpatient variability (47.3% for Cl and 68.1% for V, Table 2) might lead to higher (>40 μg/l) plasma concentrations in individual patients. Contrarily, a higher morphine dose may be needed in some patients to ensure effective treatment. Clinicians should not be reluctant to increase the maintenance dose if the starting dose proves inadequate for a patient’s clinical condition both during and, if morphine is not stopped simultaneously with TH, after hypothermia.

The current practice rarely leads to plasma concentrations exceeding the potentially toxic upper limit of 300 μg/L, but may lead to unnecessary high morphine exposure. Neonatal opioid use has been associated with impaired cognitive and behavioral development in animal studies.[39] Long term follow-up studies in humans suggest a possible negative effect in early childhood that does not persist later in life.[40–43] Conversely, adequate management of pain and discomfort is needed to improve recovery, to ensure the effectiveness of TH and to prevent adverse physiological responses such as changes in intrathoracic or arterial pressure and vasoconstriction of vital organs.[4,38,39]

Morphine PK was best described using a one-compartment model for morphine and subsequent one-compartment models for each of the glucuronide metabolites. Previously, PK of the parent compound has been adequately described using both one-compartment[44,45] and two-compartment models.[9,23,27] A recently published meta-model combining data obtained in neonates and older children from five separate studies proposed a one-compartment model for morphine PK.[30] Metabolite PK was adequately described using one-compartment models for each metabolite in all studies.[9,23,27,44] Parameter estimates from this model extrapolated to a neonate with GA 40 weeks and BW 3.5 kg result in a higher morphine clearance (1.54 l/h) and a lower volume of distribution (5.25 l) compared to our findings.[30] This might be explained by the fact that our patient population consisted of critically ill term neonates admitted to a NICU. The meta-model incorporated data form both term and preterm neonates and from older children and adults as well. In all included studies, morphine was administered for post-operative pain. Morphine parameter estimates reported in a small dataset by Frymoyer et al. in the same population are in accordance with our findings, despite the differences in the underlying PK model (Cl 0.765 l/h, V 8.02 l).[27] Also, the impact of TH on morphine clearance in our study is similar to the effect found by Róka et al., who compared hypothermic neonates to non-asphyxiated normothermic controls using a non-parametric approach (Cl 0.69 l/h vs 0.89 l/h, decrease of 22.5%).[26]

GA was identified as a significant covariate on morphine clearance despite the relatively narrow range of GA (36–42 weeks) in this population. Previous reports investigating the PK of other drugs in the PharmaCool study population have reported similar effects.[18–20] In the present study, this might be explained by a lower baseline UGT2B7 activity in neonates with a lower GA. However, due to the large interpatient variability, this finding did not translate into a dosing advise differentiated by GA. For each of the situations presented in Fig 2, the simulated dosing regimen yields average morphine plasma concentrations between 10 and 40 μg/L. The majority of patients included in this study had a GA between 38 and 41 weeks (202/244, 82.8%); a lower mg/kg dose for neonates with a GA of 36 and 37 weeks only marginally improved the fraction of patients within the therapeutic window at PNA 12 hours. Furthermore, a differentiated dosing regimen based on GA within this relatively small subpopulation of NICU patients is deemed undesirable as this will be error-prone. Therefore, the proposed dosing regimen is advised for all neonates treated with TH after HIE. Currently, TH is explored in preterm neonates (GA 34–35 weeks) as well.[46] Extrapolation of our results to these neonates should be done with caution due to the absence of this patient population in the current dataset.

Our data show an increase in morphine clearance during the first five days after birth. This effect could be identified independently of the effect of body temperature. The increase of clearance over time can be attributed to maturation of UGT2B7. Maturation of this enzyme in normothermic neonates has been described for morphine but also for other drugs predominantly glucuronidated by UGT2B7.[9,47–50] Recovery of organ function after perinatal asphyxia might also play a role in this observed increase in clearance during the first days of life. During asphyxia, the liver is deprived of oxygen, resulting in hepatocyte damage. Alanine aminotransferase (ALAT) and aspartate aminotransferase (ASAT) levels are commonly used to indicate hepatocyte damage and both are frequently elevated in neonates with HIE. Peak levels are reached within 72 hours and normalize within the first two weeks of life, indicating recovery of liver function in this timeframe.[51–53] Unfortunately, ALAT and ASAT are poor predictors for hepatic drug metabolism and cannot be reliably used to predict hepatic clearance. Therefore, it was not possible to distinguish between enzyme maturation and recovery of liver damage.

Introducing temperature as a continuous covariate allowed for a more precise estimate of the effect of body temperature on the pharmacokinetics of morphine and its metabolites since also samples during rewarming were available. Additionally, including temperature as a dichotomous covariate proved to be a less adequate fit for the data. Some assumptions had to be made for this dynamic model. Firstly, we assumed that the average body temperature of each neonate during TH was 33.5°C and that possible fluctuations between 33.0 and 34.0°C would have a negligible effect. Body temperature during TH was therefore fixed to 33.5°C. Secondly, we assumed that rewarming for each neonate occurred according to national protocol at 0.4°C/hour. In clinical practice, rewarming is sometimes slowed or halted if seizures occur during rewarming. As this information was not available in our dataset, rewarming for each neonate was set at 0.4°C/hour. Thirdly, as body temperature after TH is stabilized at 36.5°C for 24 hours, we opted for a fixed body temperature of 36.5°C for each neonate after rewarming. As morphine is primarily administrated to prevent stress during hypothermia and is often stopped with or shortly after TH, we believe that this assumption is an accurate representation of clinical practice. Lastly, we assume a linear effect between body temperature and clearance and therefore report and effect per°C. Although we have no evidence for non-linearity, the limited sampling strategy was insufficient to exclude this. However, in event of non-linearity, the reported effect is an average effect per°C between 33.5 and 36.5°C and does not alter the overall effect of TH on clearance.

SCr is a specific marker for renal function and is widely used in adults to predict reduced clearance of renally excreted drugs. In neonates, SCr levels in the first few days of life are confounded by maternal SCr levels due to maternal transfer. We considered a dynamic model of SCr over time a better predictor of changes in renal function. No relationship between serum creatinine and renal clearance of the metabolites could be identified. Additionally, increased metabolite clearance over time in the first five days of life could not be observed. This effect was found in the same population for amoxicillin and benzylpenicillin, drugs that are predominantly excreted renally in unmetabolized form.[18,20] Data collection up to five days after birth might have been too short to detect maturation of renal function since steady state M3G and M6G plasma concentrations are not reached during hypothermia (S8 Fig). Also, as maturation of renal function in the first few days after birth occurs simultaneously with TH, the effect of maturation on metabolite clearance might not be distinguishable from the effect of hypothermia.

Pharmacodynamics (PD) end points of morphine were not incorporated in the final model. Although the COMFORT-B score, as indicator for pain and stress, was routinely recorded in this population, the timing of this score in relation to morphine dosing was often unclear. Pain expression in hypothermic neonates differs from normothermic neonates, making it uncertain whether the COMFORT-B score is suitable for treatment evaluation in this population.[38,54,55] Additionally, this scale has not been developed to distinguish between adequate treatment effect (eg. adequate sedation) and supratherapeutic effects (eg. oversedation).

Strengths of this study are the large number of included patients recruited from twelve tertiary NICUs in two countries, making this study population representative for all neonates treated with TH after HIE. The dosing regimen advised in this study corresponds with the dosing advise postulated by Frymoyer et al.[27] Reconstructing the full profile of hypothermia and rewarming enabled us to accurately assess the influence of body temperature on clearance. In all participating centers, TH is applied using a uniform cooling device and with a joint treatment protocol, thereby decreasing the chance of treatment variation.[3]

Limitations are the lack of PD end points associated with morphine. Clinicians are insufficiently supported by robust tools to facilitate morphine dose adjustments. Development and validation of such an instrument should be the focus of future research. Subsequently, this tool could be used to prospectively validate our dosing regimen. Ideally, future PD studies will also incorporate M6G plasma concentrations for determining exposure-response as the contribution of this active metabolite to the effects attributed to morphine in neonates are largely unknown.

For ethical reasons, it was not possible to answer this research question using a prospective randomized controlled trial comparing hypothermic with non-hypothermic patients. Also, comparison to an adequate historical control group is not feasible as morphine nor metabolite plasma concentrations are available from before 2008.

Conclusion

Clearance of morphine and its metabolites is reduced in neonates treated with TH for HIE. Even though the current clinical practice only very rarely leads to morphine plasma concentrations exceeding 300 μg/L, a relatively low starting dose of 50 μg/kg followed by continuous infusion of 5 μg/kg/h is recommended in critically ill neonates treated with TH for HIE. However, due to the large interpatient variability, the uncertainty regarding the supposed therapeutic window and the undesirable effect of discomfort in this population, a higher maintenance dose may be required if the starting dose proves inadequate for the clinical condition of the individual patient.

Supporting information

[Figure omitted. See PDF.]

S1 Fig.

observed plasma concentrations for M3G (left) and M3G (right). M3G = morphine-3-glucuronde, M6G = morphine-6-glucuronide.

https://doi.org/10.1371/journal.pone.0211910.s001

(TIF)

S2 Fig. Morphine goodness-of-fit plots.

A = observed vs population predicted plasma concentrations; B = observed vs individual predicted plasma concentrations; C = population conditional weighted residuals vs population predicted plasma concentrations; D = population conditional weighted residuals vs time after birth; solid line indicates the linear regression line.

https://doi.org/10.1371/journal.pone.0211910.s002

(TIF)

S3 Fig. M3G goodness-of-fit plots.

A = observed vs population predicted plasma concentrations; B = observed vs individual predicted plasma concentrations; C = population conditional weighted residuals vs population predicted plasma concentrations; D = population conditional weighted residuals vs time after birth; M3G = morphine-3-glucuronide; solid line indicates the linear regression line.

https://doi.org/10.1371/journal.pone.0211910.s003

(TIF)

S4 Fig. M6G goodness-of-fit plots.

A = observed vs population predicted plasma concentrations; B = observed vs individual predicted plasma concentrations; C = population conditional weighted residuals vs population predicted plasma concentrations; D = population conditional weighted residuals vs time after birth; M6G = morphine-6-glucuronide; solid line indicates the linear regression line.

https://doi.org/10.1371/journal.pone.0211910.s004

(TIF)

S5 Fig. Normalized prediction distribution errors (NPDEs) of the final pharmacokinetic model for morphine.

A = kernel density plot of NPDE with a normal, Gaussian distribution overlaid for comparative purposes; B = Q-Q plot of theoretical quantiles vs sample quantiles; C = NPDE vs Time; D = NPDE vs predicted plasma concentrations; solid lines in figures C and D represent the observed median, 5th and 95th percentiles, red box represent the predicted 90% confidence interval around the median, blue boxes represent the predicted 90% confidence intervals around the 5th and 95th percentiles.

https://doi.org/10.1371/journal.pone.0211910.s005

(TIF)

S6 Fig. Normalized prediction distribution errors (NPDEs) of the final pharmacokinetic model for M3G.

A = kernel density plot of NPDE with a normal, Gaussian distribution overlaid for comparative purposes; B = Q-Q plot of theoretical quantiles vs sample quantiles; C = NPDE vs Time; D = NPDE vs predicted plasma concentrations; M3G = morphine-3-glucuronide; solid lines in figures C and D represent the observed median, 5th and 95th percentiles, red box represent the predicted 90% confidence interval around the median, blue boxes represent the predicted 90% confidence intervals around the 5th and 95th percentiles.

https://doi.org/10.1371/journal.pone.0211910.s006

(TIF)

S7 Fig. Normalized prediction distribution errors (NPDEs) of the final pharmacokinetic model for M6G.

A = kernel density plot of NPDE with a normal, Gaussian distribution overlaid for comparative purposes; B = Q-Q plot of theoretical quantiles vs sample quantiles; C = NPDE vs Time; D = NPDE vs predicted plasma concentrations; M6G = morphine-6-glucuronide; solid lines in figures C and D represent the observed median, 5th and 95th percentiles, red box represent the predicted 90% confidence interval around the median, blue boxes represent the predicted 90% confidence intervals around the 5th and 95th percentiles.

https://doi.org/10.1371/journal.pone.0211910.s007

(TIF)

S8 Fig.

simulated plasma concentration time profiles for M3G (left) and M6G (right) of the proposed morphine dosing regimen of 5 μg/kg/h after loading dose of 50 μg/kg. Solid line indicates the mean plasma concentration; gray area represents the 95% prediction interval. M3G = morphine-3-glucuronde, M6G = morphine-6-glucuronide.

https://doi.org/10.1371/journal.pone.0211910.s008

(TIF)

S9 Fig.

Simulated morphine plasma concentrations of the dosing regimens of 10 μg/kg/h after loading dose of 50 μg/kg (left), 5 μg/kg/h after loading dose of 100 μg/kg (center) and 10 μg/kg/h after loading dose of 100 μg/kg (right). Solid line indicates the mean morphine plasma concentration; gray area represents the 95% prediction interval. Dotted horizontal lines indicate the proposed therapeutic window of 10–40 μg/L. Solid vertical lines indicate the start and end of TH (33.5°C) simulated between 5h and 77h after birth; dashed vertical line indicates the return to normothermia (36.5°C) with rewarming simulated at 0.4°C/h; TH = therapeutic hypothermia.

https://doi.org/10.1371/journal.pone.0211910.s009

(TIF)

Acknowledgments

Lead author for the PharmaCool study group is Timo R. de Haan, [email protected].

The following PharmaCool study group members are non-author contributors: Mieke J. Brouwer PhD2, S.M. Mulder—de Tollenaer MD, PhD6, L.J.M. Groot Jebbink—Akkerman MD, PhD6, Djien Liem MD, PhD8, Katerina Steiner MD, PhD8, Sinno H.P. Simons MD, PhD9, Rogier C.J. de Jonge MD, PhD9, Annelies A. Bos MD, PhD9, Michel Sonnaert MD, PhD13, Fleur Anne Camfferman MD, PhD13. Affiliations correspond with author affiliations.

The authors would like to thank Mark de Groot for contributions with regard to data management.

Citation: Favié LMA, Groenendaal F, van den Broek MPH, Rademaker CMA, de Haan TR, van Straaten HLM, et al. (2019) Pharmacokinetics of morphine in encephalopathic neonates treated with therapeutic hypothermia. PLoS ONE 14(2): e0211910. https://doi.org/10.1371/journal.pone.0211910

References

1. Jacobs SE, Berg M, Hunt R, Tarnow-Mordi WO, Inder TE, Davis PG. Cooling for newborns with hypoxic ischaemic encephalopathy. In: Jacobs SE, editor. Cochrane Database of Systematic Reviews. Chichester, UK: John Wiley & Sons, Ltd; 2013. pp. 385–385. https://doi.org/10.1002/14651858.CD003311.pub3 pmid:23440789

2. Azzopardi D V, Strohm B, Edwards AD, Dyet L, Halliday HL, Juszczak E, et al. Moderate hypothermia to treat perinatal asphyxial encephalopathy. N Engl J Med. 2009;361: 1349–1358. pmid:19797281

3. Groenendaal F, Casaer A, Dijkman KP, Gavilanes AWD, De Haan TR, Ter Horst HJ, et al. Introduction of hypothermia for neonates with perinatal asphyxia in the Netherlands and flanders and the dutch-flemish working group on neonatal neurology. Neonatology. 2013;104: 15–21. pmid:23615314

4. Thoresen M, Satas S, LØberg EM, Whitelaw A, Acolet D, Lindgren C, et al. Twenty-four hours of mild hypothermia in unsedated newborn pigs starting after a severe global hypoxic-ischemic insult is not neuroprotective. Pediatr Res. 2001;50: 405–411. pmid:11518829

5. De Gregori S, De Gregori M, Ranzani GN, Allegri M, Minella C, Regazzi M. Morphine metabolism, transport and Brain disposition. Metabolic Brain Disease. 2012. pp. 1–5. pmid:22193538

6. Allegaert K, Van De Velde M, Van Den Anker J. Neonatal clinical pharmacology. Paediatr Anaesth. 2014;24: 30–38. pmid:23617305

7. De Wildt SN, Kearns GL, Leeder JS, Van Den Anker JN. Glucuronidation in humans. Pharmacogenetic and developmental aspects. Clinical Pharmacokinetics. 1999. pp. 439–452. pmid:10427468

8. Kearns GL, Abdel-Rahman SM, Alander SW, Blowey DL, Leeder JS, Kauffman RE. Developmental Pharmacology—Drug Disposition, Action, and Therapy in Infants and Children. N Engl J Med. 2003;349: 1157–1167. pmid:13679531

9. Knibbe CAJ, Krekels EHJ, Van Den Anker JN, De Jongh J, Santen GWE, Van Dijk M, et al. Morphine glucuronidation in preterm neonates, infants and children younger than 3 years. Clin Pharmacokinet. 2009;48: 371–385. pmid:19650676

10. Zanelli S, Buck M, Fairchild K, Zhou J, Poloyac SM. Physiologic and pharmacologic considerations for hypothermia therapy in neonates. Expert Opin Drug Metab Toxicol. Nature Publishing Group; 2011;7: 803–16.

11. Van den Broek MPH, Groenendaal F, Egberts ACG, Rademaker CMA. Effects of Hypothermia on Pharmacokinetics and Pharmacodynamics. Clin Pharmacokinet. 2010;49: 277–294. pmid:20384391

12. Zhou J, Poloyac SM. The effect of therapeutic hypothermia on drug metabolism and response: cellular mechanisms to organ function. Expert Opin Drug Metab Toxicol. 2011;7: 803–16. pmid:21473710

13. Tarcan A, Tiker F, Güvenir H, Gürakan B. Hepatic involvement in perinatal asphyxia. J Matern Neonatal Med. 2007;20: 407–410. pmid:17674246

14. Shah P, Riphagen S, Beyene J, Perlman M. Multiorgan dysfunction in infants with post-asphyxial hypoxic-ischaemic encephalopathy. Arch Dis Child Fetal Neonatal Ed. 2004;89: F152–5. pmid:14977901

15. Van Den Broek MPH, Van Straaten HLM, Huitema ADR, Egberts T, Toet MC, De Vries LS, et al. Anticonvulsant effectiveness and hemodynamic safety of midazolam in full-term infants treated with hypothermia. Neonatology. 2015;107: 150–156. pmid:25572061

16. Van den Broek MPH, Rademaker CMA, van Straaten HLM, Huitema ADR, Toet MC, de Vries LS, et al. Anticonvulsant treatment of asphyxiated newborns under hypothermia with lidocaine: efficacy, safety and dosing. Arch Dis Child Fetal Neonatal Ed. 2013;98: F341–5. pmid:23303304

17. Van Den Broek MPH, Groenendaal F, Toet MC, Van Straaten HLM, Van Hasselt JGC, Huitema ADR, et al. Pharmacokinetics and clinical efficacy of phenobarbital in asphyxiated newborns treated with hypothermia: A thermopharmacological approach. Clin Pharmacokinet. 2012;51: 671–679. pmid:23018530

18. Bijleveld YA, Mathôt RAA, van der Lee JH, Groenendaal F, Dijk PH, van Heijst A, et al. Population Pharmacokinetics of Amoxicillin in Term Neonates Undergoing Moderate Hypothermia. Clin Pharmacol Ther. 2017;00.

19. Bijleveld YA, De Haan TR, van der Lee JH, Groenendaal F, Dijk PH, van Heijst A, et al. Altered gentamicin pharmacokinetics in term neonates undergoing controlled hypothermia. Br J Clin Pharmacol. 2016;81: 1067–1077. pmid:26763684

20. Bijleveld YA, de Haan TR, van der Lee JH, Groenendaal F, Dijk PH, van Heijst A, et al. Evaluation of a System-specific Function to Describe the Pharmacokinetics of Benzylpenicillin in Term Neonates Undergoing Moderate Hypothermia. Antimicrob Agents Chemother. 2018; AAC.02311–17. pmid:29378710

21. Frymoyer A, Meng L, Bonifacio SL, Verotta D, Guglielmo BJ. Gentamicin pharmacokinetics and dosing in neonates with hypoxic ischemic encephalopathy receiving hypothermia. Pharmacotherapy. 2013; pmid:23553582

22. Baarslag MA, Allegaert K, Van Den Anker JN, Knibbe CAJ, Van Dijk M, Simons SHP, et al. Paracetamol and morphine for infant and neonatal pain; still a long way to go? Expert Review of Clinical Pharmacology. 2017. pp. 111–126. pmid:27785937

23. Wang C, Sadhavisvam S, Krekels EHJ, Dahan A, Tibboel D, Danhof M, et al. Developmental changes in morphine clearance across the entire paediatric age range are best described by a bodyweight-dependent exponent model. Clin Drug Investig. 2013;33: 523–534. pmid:23754691

24. Krekels EHJ, DeJongh J, Van Lingen RA, Van Der Marel CD, Choonara I, Lynn AM, et al. Predictive performance of a recently developed population pharmacokinetic model for morphine and its metabolites in new datasets of (preterm) neonates, infants and children. Clin Pharmacokinet. 2011;50: 51–63.

25. Krekels EHJ, Tibboel D, De Wildt SN, Ceelie I, Dahan A, Van Dijk M, et al. Evidence-based morphine dosing for postoperative neonates and infants. Clin Pharmacokinet. 2014;53: 553–563. pmid:24496960

26. Róka A, Melinda KT, Vásárhelyi B, Machay T, Azzopardi D, Szabó M. Elevated morphine concentrations in neonates treated with morphine and prolonged hypothermia for hypoxic ischemic encephalopathy. Pediatrics. 2008;121: e844–9. pmid:18381513

27. Frymoyer A, Bonifacio SL, Drover DR, Su F, Wustoff CJ, Van Meurs KP. Decreased Morphine Clearance in Neonates With Hypoxic Ischemic Encephalopathy Receiving Hypothermia. J Clin Pharmacol. 2017;57: 64–76. pmid:27225747

28. De Haan TR, Bijleveld YA, van der Lee JH, Groenendaal F, van den Broek MPH, Rademaker CMA, et al. Pharmacokinetics and pharmacodynamics of medication in asphyxiated newborns during controlled hypothermia. The PharmaCool multicenter study. BMC Pediatr. 2012;12: 45. pmid:22515424

29. Keizer RJ, Karlsson MO, Hooker A. Modeling and simulation workbench for NONMEM: Tutorial on Pirana, PsN, and Xpose. CPT Pharmacometrics Syst Pharmacol. 2013; pmid:23836189

30. Knøsgaard KR, Foster DJR, Kreilgaard M, Sverrisdóttir E, Upton RN, van den Anker JN. Pharmacokinetic models of morphine and its metabolites in neonates: Systematic comparisons of models from the literature, and development of a new meta-model. Eur J Pharm Sci. Elsevier B.V.; 2016;92: 117–130. pmid:27373670

31. Dosne AG, Bergstrand M, Harling K, Karlsson MO. Improving the estimation of parameter uncertainty distributions in nonlinear mixed effects models using sampling importance resampling. J Pharmacokinet Pharmacodyn. 2016;43: 583–596. pmid:27730482

32. Brendel K, Comets E, Laffont C, Laveille C, Mentré F. Metrics for external model evaluation with an application to the population pharmacokinetics of gliclazide. Pharm Res. 2006; pmid:16906454

33. Olkkola KT, Maunuksela EL, Korpela R, Rosenberg PH. Kinetics and dynamics of postoperative intravenous morphine in children. Clin Pharmacol Ther. 1988;44: 128–136. pmid:3135138

34. Bouwmeester NJ, van den Anker JN, Hop WCJ, Anand KJS, Tibboel D. Age- and therapy-related effects on morphine requirements and plasma concentrations of morphine and its metabolites in postoperative infants. Br J Anaesth. 2003;90: 642–652. pmid:12697593

35. Kart T, Christrup LL, Rasmussen M. Recommended use of morphine in neonates, infants and children based on a literature review: Part 2—Clinical use. Paediatr Anaesth. 1997;7: 93–101. pmid:9188108

36. Pacifici G. Metabolism and pharmacokinetics of morphine in neonates: A review. Clinics. 2016;71: 474–480. pmid:27626479

37. Chay PCW, Duffy BJ, Walker JS. Pharmacokinetic-pharmacodynamic relationships of morphine in neonates. Clin Pharmacol Ther. 1992;51: 334–342. pmid:1544290

38. Zeilmaker GA, Pokorna P, Mian P, Wildschut ED, Knibbe CAJ, Krekels EHJ, et al. Pharmacokinetic considerations for pediatric patients receiving analgesia in the intensive care unit; targeting postoperative, ECMO and hypothermia patients. Expert Opin Drug Metab Toxicol. Taylor & Francis; 2018;00: 1–12.

39. Attarian S, Tran LC, Moore A, Stanton G, Meyer E, Moore RP. The neurodevelopmental impact of neonatal morphine administration. Brain Sciences. 2014. pp. 321–334. pmid:24961764

40. De Graaf J, Van Lingen RA, Simons SHP, Anand KJS, Duivenvoorden HJ, Weisglas-Kuperus N, et al. Long-term effects of routine morphine infusion in mechanically ventilated neonates on children’s functioning: Five-year follow-up of a randomized controlled trial. Pain. 2011;152: 1391–1397. pmid:21402444

41. Van Den Bosch GE, White T, El Marroun H, Simons SHP, Van Der Lugt A, Van Der Geest JN, et al. Prematurity, Opioid Exposure and Neonatal Pain: Do They Affect the Developing Brain? Neonatology. 2015;108: 8–15. pmid:25871803

42. Steinhorn R, McPherson C, Anderson PJ, Neil J, Doyle LW, Inder T. Neonatal morphine exposure in very preterm infants—Cerebral development and outcomes. J Pediatr. 2015;166: 1200–1207.e4. pmid:25919729

43. De Graaf J, Van Lingen RA, Valkenburg AJ, Weisglas-Kuperus N, Groot Jebbink L, Wijnberg-Williams B, et al. Does neonatal morphine use affect neuropsychological outcomes at 8 to 9 years of age? Pain. 2013;154: 449–458. pmid:23352760

44. Bouwmeester NJ, Anderson BJ, Tibboel D, Holford NHG. Developmental pharmacokinetics of morphine and its metabolites in neonates, infants and young children. Br J Anaesth. 2004;92: 208–217. pmid:14722170

45. Anand KJS, Anderson BJ, Holford NHG, Hall RW, Young T, Shephard B, et al. Morphine pharmacokinetics and pharmacodynamics in preterm and term neonates: Secondary results from the NEOPAIN trial. Br J Anaesth. 2008;101: 680–689. pmid:18723857

46. Rao R, Trivedi S, Vesoulis Z, Liao SM, Smyser CD, Mathur AM. Safety and Short-Term Outcomes of Therapeutic Hypothermia in Preterm Neonates 34–35 Weeks Gestational Age with Hypoxic-Ischemic Encephalopathy. J Pediatr. 2017;183: 37–42. pmid:27979578

47. Cook SF, Stockmann C, Samiee-Zafarghandy S, King AD, Deutsch N, Williams EF, et al. Neonatal Maturation of Paracetamol (Acetaminophen) Glucuronidation, Sulfation, and Oxidation Based on a Parent–Metabolite Population Pharmacokinetic Model. Clin Pharmacokinet. Springer International Publishing; 2016;55: 1395–1411. pmid:27209292

48. Krekels EHJ, Danhof M, Tibboel D, Knibbe CAJ. Ontogeny of hepatic glucuronidation; methods and results. Curr Drug Metab. 2012;13: 728–43. pmid:22452455

49. Van Donge T, Mian P, Tibboel D, Van Den Anker J, Allegaert K. Drug metabolism in early infancy: opioids as an illustration. Expert Opin Drug Metab Toxicol. Taylor & Francis; 2018;14: 287–301. pmid:29363349

50. Capparelli E V., Englund JA, Connor JD, Spector SA, McKinney RE, Palumbo P, et al. Population pharmacokinetics and pharmacodynamics of zidovudine in HIV-infected infants and children. J Clin Pharmacol. 2003;43: 133–140. pmid:12616665

51. Muniraman H, Gardner D, Skinner J, Paweletz A, Vayalakkad A, Chee YH, et al. Biomarkers of hepatic injury and function in neonatal hypoxic ischemic encephalopathy and with therapeutic hypothermia. Eur J Pediatr. European Journal of Pediatrics; 2017;176: 1295–1303. pmid:28741035

52. Choudhary M, Sharma D, Dabi D, Lamba M, Pandita A, Shastri S. Hepatic Dysfunction in Asphyxiated Neonates: Prospective Case-Controlled Study. Clin Med Insights Pediatr. 2015;9: CMPed.S21426. pmid:25674030

53. Chhavi N, Zutshi K, Singh NK, Awasthi A, Goel A. Serum liver enzyme pattern in birth asphyxia associated liver injury. Pediatr Gastroenterol Hepatol Nutr. 2014;17: 162–9. pmid:25349832

54. Foulkes T, Wood JN. Mechanisms of cold pain. Channels. 2007. pp. 154–160. pmid:18690033

55. Rutkove SB. Effects of temperature on neuromuscular electrophysiology. Muscle and Nerve. 2001. pp. 867–882. pmid:11410914