When will the Glomerular Filtration Rate in Former Preterm Neonates Catch up with Their Term Peers?

AIMS: Whether and when glomerular filtration rate (GFR) in preterms catches up with term peers is unknown. This study aims to develop a GFR maturation model for (pre)term-born individuals from birth to 18 years of age. Secondarily, the function is applied to data of different renally excreted drugs. METHODS: We combined published inulin clearance values and serum creatinine (Scr) concentrations in (pre)term born individuals throughout childhood. Inulin clearance was assumed to be equal to GFR, and Scr to reflect creatinine synthesis rate/GFR. We developed a GFR function consisting of GFRbirth (GFR at birth), and an Emax model dependent on PNA (with GFRmax, PNA50 (PNA at which half of GFR max is reached) and Hill coefficient). The final GFR model was applied to predict gentamicin, tobramycin and vancomycin concentrations. RESULT: In the GFR model, GFRbirth varied with birthweight linearly while in the PNA-based Emax equation, GA was the best covariate for PNA50, and current weight for GFRmax. The final model showed that for a child born at 26 weeks GA, absolute GFR is 18%, 63%, 80%, 92% and 96% of the GFR of a child born at 40 weeks GA at 1 month, 6 months, 1 year, 3 years and 12 years, respectively. PopPK models with the GFR maturation equations predicted concentrations of renally cleared antibiotics across (pre)term-born neonates until 18 years well. CONCLUSIONS: GFR of preterm individuals catches up with term peers at around three years of age, implying reduced dosages of renally cleared drugs should be considered below this age.

Maturation of Paracetamol Elimination Routes in Preterm Neonates Born Below 32 Weeks of Gestation.

PURPOSE: Despite being off-label, intravenous paracetamol (PCM) is increasingly used to control mild-to-moderate pain in preterm neonates. Here we aim to quantify the maturation of paracetamol elimination pathways in preterm neonates born below 32 weeks of gestation. METHODS: Datasets after single dose (rich data) or multiple doses (sparse data) of intravenous PCM dose (median (range)) 9 (3-25) mg/kg were pooled, containing 534 plasma and 44 urine samples of PCM and metabolites (PCM-glucuronide, PCM-sulfate, PCM-cysteine, and PCM-mercapturate) from 143 preterm neonates (gestational age 27.7 (24.0-31.9) weeks, birthweight 985 (462-1,925) g, postnatal age (PNA) 5 (0-30) days, current weight 1,012 (462-1,959) g. Population pharmacokinetic analysis was performed using NONMEM® 7.4. RESULTS: For a typical preterm neonate (birthweight 985 g; PNA 5 days), PCM clearance was 0.137 L/h, with glucuronidation, sulfation, oxidation and unchanged renal clearance accounting for 5.3%, 73.7%, 16.3% and 4.6%, respectively. Maturational changes in total PCM clearance and its elimination pathways were best described by birthweight and PNA. Between 500-1,500 g birthweight, total PCM clearance increases by 169%, with glucuronidation, sulfation and oxidation clearance increasing by 347%, 164% and 164%. From 1-30 days PNA for 985 g birthweight neonate, total PCM clearance increases by 167%, with clearance via glucuronidation and oxidation increasing by 551%, and sulfation by 69%. CONCLUSION: Birthweight and PNA are the most important predictors for maturational changes in paracetamol clearance and its glucuronidation, sulfation and oxidation. As a result, dosing based on bodyweight alone will not lead to consistent paracetamol concentrations among preterm neonates.

Total bodyweight and sex both drive pharmacokinetic variability of fluconazole in obese adults.

BACKGROUND: Fluconazole is commonly used to treat or prevent fungal infections. It is typically used orally but in critical situations, IV administration is needed. Obesity may influence the pharmacokinetics and therapeutic efficacy of a drug. In this study, we aim to assess the impact of obesity on fluconazole pharmacokinetics given orally or IV to guide dose adjustments for the obese population. METHODS: We performed a prospective pharmacokinetic study with intensive sampling in obese subjects undergoing bariatric surgery (n = 17, BMI ≥ 35 kg/m2) and non-obese healthy controls (n = 8, 18.5 ≤ BMI < 30.0 kg/m2). Participants received a semi-simultaneous oral dose of 400 mg fluconazole capsules, followed after 2 h by 400 mg IV. Population pharmacokinetic modelling and simulation were performed using NONMEM 7.3. RESULTS: A total of 421 fluconazole concentrations in 25 participants (total bodyweight 61.0-174 kg) until 48 h after dosing were obtained. An estimated bioavailability of 87.5% was found for both obese and non-obese subjects, with a 95% distribution interval of 43.9%-98.4%. With increasing total bodyweight, both higher CL and Vd were found. Sex also significantly impacted Vd, being 27% larger in male compared with female participants. CONCLUSIONS: In our population of obese but otherwise healthy individuals, obesity clearly alters the pharmacokinetics of fluconazole, which puts severely obese adults, particularly if male, at risk of suboptimal exposure, for which adjusted doses are proposed.

Modelling inflammatory biomarker dynamics in a human lipopolysaccharide (LPS) challenge study using delay differential equations.

Clinical studies in healthy volunteers challenged with lipopolysaccharide (LPS), a constituent of the cell wall of Gram-negative bacteria, represent a key model to characterize the Toll-like receptor 4 (TLR4)-mediated inflammatory response. Here, we developed a mathematical modelling framework to quantitatively characterize the dynamics and inter-individual variability of multiple inflammatory biomarkers in healthy volunteer LPS challenge studies. Data from previously reported LPS challenge studies were used, which included individual-level time-course data for tumour necrosis factor α (TNF-α), interleukin 6 (IL-6), interleukin 8 (IL-8) and C-reactive protein (CRP). A one-compartment model with first-order elimination was used to capture the LPS kinetics. The relationships between LPS and inflammatory markers was characterized using indirect response (IDR) models. Delay differential equations were applied to quantify the delays in biomarker response profiles. For LPS kinetics, our estimates of clearance and volume of distribution were 35.7 L h-1 and 6.35 L, respectively. Our model adequately captured the dynamics of multiple inflammatory biomarkers. The time delay for the secretion of TNF-α, IL-6 and IL-8 were estimated to be 0.924, 1.46 and 1.48 h, respectively. A second IDR model was used to describe the induced changes of CRP in relation to IL-6, with a delayed time of 4.2 h. The quantitative models developed in this study can be used to inform design of clinical LPS challenge studies and may help to translate preclinical LPS challenge studies to humans.

Population pharmacokinetics-pharmacodynamics of fondaparinux in dialysis-dependent chronic kidney disease patients undergoing chronic renal replacement therapy.

PURPOSE: Data on the anti-Xa efficacy of fondaparinux in dialysis-dependent chronic kidney disease (DD-CKD) patients are scarce. This study characterizes the pharmacokinetics (PK) and pharmacodynamics (PD) of fondaparinux in DD-CKD patients undergoing renal replacement therapy (RRT), to assess dosing strategies. METHODS: A retrospective, observational study was conducted using data on anti-Xa activity (112 samples) from 12 (3 male and 9 female) DD-CKD patients (median (IQR) age 71 years (63-88), weight 73 kg (59-98.5)). Eleven patients underwent high-flux or low-flux hemodialysis (HD) and one patient underwent peritoneal dialysis. Three patients were also treated with therapeutic plasma exchange (TPE). A non-linear mixed effects analysis was performed using NONMEM 7.3.0. RESULTS: The lab-specific slope of the relationship between fondaparinux concentration and anti-Xa levels was 1.18 IU/µg. In a one-compartment model, clearance (CL) and volume of distribution (Vd) were 0.05289 L/h and 5.55 L, respectively. High-flux HD was found to increase the CL of fondaparinux 2.26 times. TPE also considerably increased CL, but the fold-change could not be accurately estimated. Low-flux HD and peritoneal dialysis did not impact PK parameters. CONCLUSIONS: Model-based simulations showed that standard dosing (2.5 mg three times weekly before HD) results in a median anti-Xa activity of 0.55 IU/mL and 0.98 IU/mL, pre- and post-low-flux HD, respectively. In patients undergoing high-flux HD, these values are approximately 27% lower. Additional caution is warranted with TPE, as this treatment can reduce anti-Xa activity even further.

Effect of Different Exercise Training Modalities on Fasting Levels of Oxylipins and Endocannabinoids in Middle-Aged Sedentary Adults: A Randomized Controlled Trial.

This study aimed to investigate the effects of different exercise training programs on fasting plasma levels of oxylipins, endocannabinoids (eCBs), and eCBs-like molecules in middle-aged sedentary adults. A 12-week randomized controlled trial was conducted using a parallel group design. Sixty-five middle-aged adults (40-65 years old) were randomly assigned to: (a) no exercise (control group), (b) concurrent training based on international physical activity recommendations (PAR group), (c) high-intensity interval training (HIIT group), and (d) HIIT together with whole-body electromyostimulation (HIIT + EMS group). Plasma levels of oxylipins, eCBs, and eCBs-like molecules were determined in plasma samples before and after the intervention using targeted lipidomics. Body composition was assessed through dual-energy X-ray absorptiometry, and dietary intake through a food frequency questionnaire and three nonconsecutive 24-hr recalls. The physical activity recommendations, HIIT, and HIIT-EMS groups showed decreased plasma levels of omega-6 and omega-3-derived oxylipins, and eCBs and eCBs-like molecules after 12 weeks (all Δ ≤ -0.12; all p < .05). Importantly, after Bonferroni post hoc corrections, the differences in plasma levels of omega-6 and omega-3 oxylipins were not statistically significant compared with the control group (all p > .05). However, after post hoc corrections, plasma levels of anandamide and oleoylethanolamide were increased in the physical activity recommendations group compared with the control group (anandamide: Δ = 0.05 vs. -0.09; oleoylethanolamide: Δ = -0.12 vs. 0.013, all p ≤ .049). In conclusion, this study reports that a 12-week exercise training intervention, independent of the modality applied, does not modify fasting plasma levels of omega-6 and omega-3 oxylipins, eCBs, and eCBs-like molecules in middle-aged sedentary adults.

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.

Midazolam Infusion and Disease Severity Affect the Level of Sedation in Children: A Parametric Time-to-Event Analysis.

AIM: In critically ill mechanically ventilated children, midazolam is used first line for sedation, however its exact sedative effects have been difficult to quantify. In this analysis, we use parametric time-to-event (PTTE) analysis to quantify the effects of midazolam in critically ill children. METHODS: In the PTTE analysis, data was analyzed from a published study in mechanically ventilated children in which blinded midazolam or placebo infusions were administered during a sedation interruption phase until, based on COMFORT-B and NISS scores, patients became undersedated and unblinded midazolam was restarted. Using NONMEM® v.7.4.3., restart of unblinded midazolam was analysed as event. Patients in the trial were divided into internal and external validation cohorts prior to analysis. RESULTS: Data contained 138 events from 79 individuals (37 blinded midazolam; 42 blinded placebo). In the PTTE model, the baseline hazard was best described by a constant function. Midazolam reduced the hazard for restart of unblinded midazolam due to undersedation by 51%. In the blinded midazolam group, time to midazolam restart was 26 h versus 58 h in patients with low versus high disease severity upon admission (PRISM II < 10 versus > 21), respectively. For blinded placebo, these times were 14 h and 33 h, respectively. The model performed well in an external validation with 42 individuals. CONCLUSION: The PTTE analysis effectively quantified the effect of midazolam in prolonging sedation and also the influence of disease severity on sedation in mechanically ventilated critically ill children, and provides a valuable tool to quantify the effect of sedatives. Clinical trial number and registry URL: Netherlands Trial Register, Trial NL1913 (NTR2030), date registered 28 September 2009 https://www.trialregister.nl/trial/1913 .

Predicting Unacceptable Pain in Cardiac Surgery Patients Receiving Morphine Maintenance and Rescue Doses: A Model-Based Pharmacokinetic-Pharmacodynamic Analysis.

BACKGROUND: Optimal analgesic treatment following cardiac surgery is crucial for both patient comfort and successful postoperative recovery. While knowledge of both the pharmacokinetics and pharmacodynamics of analgesics is required to predict optimal drug dosing, models quantifying the pharmacodynamics are scarce. Here, we quantify the pharmacodynamics of morphine by modeling the need for rescue morphine to treat unacceptable pain in 118 patients after cardiac surgery. METHODS: The rescue morphine event data were analyzed with repeated time-to-event (RTTE) modeling using NONMEM. Postoperative pain titration protocol consisted of continuous morphine infusions (median duration 20.5 hours) with paracetamol 4 times daily and rescue morphine in case of unacceptable pain (numerical rating scale ≥4). RESULTS: Patients had a median age of 73 years (interquartile range [IQR]: 63-77) and median bodyweight of 80 kg (IQR: 72-90 kg). Most patients (55%) required at least 1 rescue morphine dose. The hazard for rescue morphine following cardiac surgery was found to be significantly influenced by time after surgery, a day/night cycle with a peak at 23:00 (95% confidence interval [CI], 19:35-02:03) each day, and an effect of morphine concentration with 50% hazard reduction at 9.3 ng·mL-1 (95% CI, 6.7-16). CONCLUSIONS: The pharmacodynamics of morphine after cardiac surgery was successfully quantified using RTTE modeling. Future studies can be used to expand the model to better predict morphine's pharmacodynamics on the individual level and to include the pharmacodynamics of other analgesics so that improved postoperative pain treatment protocols can be developed.

Rapid Increase in Clearance of Phenobarbital in Neonates on Extracorporeal Membrane Oxygenation: A Pilot Retrospective Population Pharmacokinetic Analysis.

OBJECTIVES: This study characterizes the changes in the pharmacokinetics of phenobarbital associated with extracorporeal membrane oxygenation treatment in neonates, to illustrate our findings and provide guidance on dosing. DESIGN: Retrospective pilot population pharmacokinetic analysis. SETTING: Neonatal ICU. PATIENTS: Thirteen critically ill neonates (birth body weight, 3.21 kg [2.65-3.72 kg]; postnatal age at start of treatment: 2 d [0-7 d]; gestational age: 38 wk [38-41 wk]) receiving venovenous or venoarterial extracorporeal membrane oxygenation. INTERVENTIONS: Phenobarbital administered in a loading dose of 7.5 mg/kg (8.5-16 mg/kg) and maintenance dose of 6.9 mg/kg/d (4.5-8.5 mg/kg/d). MEASUREMENTS AND MAIN RESULTS: Therapeutic drug monitoring data were available, yielding 5, 31, and 19 phenobarbital concentrations before, during, and after extracorporeal membrane oxygenation, respectively. Population pharmacokinetic analysis was performed using NONMEM 7.3.0 (ICON Development Solutions, Ellicott City, MD). Maturation functions for clearance and volume of distribution were obtained from literature. In a one-compartment model, clearance and volume of distribution for a typical neonate off extracorporeal membrane oxygenation and with a median birth body weight (3.21 kg) at median postnatal age (2 d) were 0.0096 L/hr (relative SE = 11%)) and 2.72 L (16%), respectively. During extracorporeal membrane oxygenation, clearance was found to linearly increase with time. Upon decannulation, phenobarbital clearance initially decreased and subsequently increased slowly driven by maturation. Extracorporeal membrane oxygenation-related changes in volume of distribution could not be identified, possibly due to sparse data collection shortly after extracorporeal membrane oxygenation start. According to the model, target attainment is achieved in the first 12 days of extracorporeal membrane oxygenation with a regimen of a loading dose of 20 mg/kg and a maintenance dose of 4 mg/kg/d divided in two doses with an increase of 0.25 mg/kg every 12 hours during extracorporeal membrane oxygenation treatment. CONCLUSIONS: We found a time-dependent increase in phenobarbital clearance during the first 12 days of extracorporeal membrane oxygenation treatment in neonates, which results in continuously decreasing phenobarbital exposure and increases the risk of therapeutic failure over time. Due to high unexplained variability, frequent and repeated therapeutic drug monitoring should be considered even with the model-derived regimen.

Pharmacokinetics and Pharmacodynamics of Drugs in Obese Pediatric Patients: How to Map Uncharted Clinical Territories.

Clinicians are increasingly faced with challenges regarding the pharmacological treatment of obese pediatric patients. To provide guidance for these treatments, a better understanding of the impact of obesity on pharmacological processes in children is needed. Results on pharmacological studies in adults show however ambiguous patterns regarding the impact of obesity on ADME processes or on drug pharmacodynamics. Additionally, based on the limited research performed in obese pediatric patients, it becomes clear that findings from obese adults cannot be expected to always translate directly to similar findings in obese children. To improve knowledge on drug pharmacology in obese pediatric patients, studies should focus on quantifying the impact of maturation, obesity, and other relevant variables on primary pharmacological parameters and on disentangling systemic (renal and/or hepatic) and presystemic (gut and/or first-pass hepatic) clearance. For this, data is required from well-designed clinical trials that include patients with not only a wide range in age but also a range in excess body weight, upon oral and intravenous dosing. Population modelling approaches are ideally suitable for this purpose and can also be used to link the pharmacokinetics to pharmacodynamics and to derive drug dosing regimens. Generalizability of research findings can be achieved by including mechanistic aspects in the data analysis, for instance, using either extrapolation approaches in population modelling or by applying physiologically based modelling principles. It is imperative that more and smarter studies are performed in obese pediatric patients to provide safe and effective treatment for this special patient population.

Mechanistic and Quantitative Understanding of Pharmacokinetics in Zebrafish Larvae through Nanoscale Blood Sampling and Metabolite Modeling of Paracetamol.

Zebrafish larvae are increasingly used for pharmacological research, but internal drug exposure is often not measured. Understanding pharmacokinetics is necessary for reliable translation of pharmacological results to higher vertebrates, including humans. Quantification of drug clearance and distribution requires measurements of blood concentrations. Additionally, measuring drug metabolites is of importance to understand clearance in this model organism mechanistically. We therefore mechanistically studied and quantified pharmacokinetics in zebrafish larvae, and compared this to higher vertebrates, using paracetamol (acetaminophen) as a paradigm compound. A method was developed to sample blood from zebrafish larvae 5 days post fertilization. Blood concentrations of paracetamol and its major metabolites, paracetamol-glucuronide and paracetamol-sulfate, were measured. Blood concentration data were combined with measured amounts in larval homogenates and excreted amounts and simultaneously analyzed through nonlinear mixed-effects modeling, quantifying absolute clearance and distribution volume. Blood sampling from zebrafish larvae was most successful from the posterior cardinal vein, with a median volume (interquartile range) of 1.12 nl (0.676-1.66 nl) per blood sample. Samples were pooled (n = 15-35) to reach measurable levels. Paracetamol blood concentrations at steady state were only 10% of the external paracetamol concentration. Paracetamol-sulfate was the major metabolite, and its formation was quantified using a time-dependent metabolic formation rate. Absolute clearance and distribution volume correlated well with reported values in higher vertebrates, including humans. Based on blood concentrations and advanced data analysis, the mechanistic and quantitative understanding of paracetamol pharmacokinetics in zebrafish larvae has been established. This will improve the translational value of this vertebrate model organism in drug discovery and development. SIGNIFICANCE STATEMENT: In early phases of drug development, new compounds are increasingly screened in zebrafish larvae, but the internal drug exposure is often not taken into consideration. We developed innovative experimental and computational methods, including a blood-sampling technique, to measure the paradigm drug paracetamol (acetaminophen) and its major metabolites and quantify pharmacokinetics (absorption, distribution, elimination) in zebrafish larvae of 5 days post fertilization with a total volume of only 300 nl. These parameter values were scaled to higher vertebrates, including humans.

Enteral Acetaminophen Bioavailability in Pediatric Intensive Care Patients Determined With an Oral Microtracer and Pharmacokinetic Modeling to Optimize Dosing.

OBJECTIVES: Decreasing morbidity and mortality by rationalizing drug treatment in the critically ill is of paramount importance but challenging as the underlying clinical condition may lead to large variation in drug disposition and response. New microtracer methodology is now available to gain knowledge on drug disposition in the intensive care. On the basis of studies in healthy adults, physicians tend to assume that oral doses of acetaminophen will be completely absorbed and therefore prescribe the same dose per kilogram for oral and IV administration. As the oral bioavailability of acetaminophen in critically ill children is unknown, we designed a microtracer study to shed a light on this issue. DESIGN: An innovative microtracer study design with population pharmacokinetics. SETTING: A tertiary referral PICU. PATIENTS: Stable critically ill children, 0-6 years old, and already receiving IV acetaminophen. INTERVENTIONS: Concomitant administration of an oral C radiolabeled acetaminophen microtracer (3 ng/kg) with IV acetaminophen treatment (15 mg/kg every 6 hr). MEASUREMENTS: Blood was drawn from an indwelling arterial or central venous catheter up to 24 hours after C acetaminophen microtracer administration. Acetaminophen concentrations were measured by liquid chromatography-mass spectrometry and C concentrations by accelerated mass spectrometry. MAIN RESULTS: In 47 patients (median age of 6.1 mo; Q1-Q3, 1.8-20 mo) the mean enteral bioavailability was 72% (range, 11-91%). With a standard dose (15 mg/kg 4 times daily), therapeutic steady-state concentrations were 2.5 times more likely to be reached with IV than with oral administration. CONCLUSIONS: Microtracer studies present a new opportunity to gain knowledge on drug disposition in the intensive care. Using this modality in children in the pediatric intensive care, we showed that enteral administration of acetaminophen results in less predictable exposure and higher likelihood of subtherapeutic blood concentration than does IV administration. IV dosing may be preferable to ensure adequate pain relief.

Dose-linearity of the pharmacokinetics of an intravenous [14 C]midazolam microdose in children.

AIMS: Drug disposition in children may vary from adults due to age-related variation in drug metabolism. Microdose studies present an innovation to study pharmacokinetics (PK) in paediatrics; however, they should be used only when the PK is dose linear. We aimed to assess dose linearity of a [14 C]midazolam microdose, by comparing the PK of an intravenous (IV) microtracer (a microdose given simultaneously with a therapeutic midazolam dose), with the PK of a single isolated microdose. METHODS: Preterm to 2-year-old infants admitted to the intensive care unit received [14 C]midazolam IV as a microtracer or microdose, followed by dense blood sampling up to 36 hours. Plasma concentrations of [14 C]midazolam and [14 C]1-hydroxy-midazolam were determined by accelerator mass spectrometry. Noncompartmental PK analysis was performed and a population PK model was developed. RESULTS: Of 15 infants (median gestational age 39.4 [range 23.9-41.4] weeks, postnatal age 11.4 [0.6-49.1] weeks), 6 received a microtracer and 9 a microdose of [14 C]midazolam (111 Bq kg-1 ; 37.6 ng kg-1 ). In a 2-compartment PK model, bodyweight was the most significant covariate for volume of distribution. There was no statistically significant difference in any PK parameter between the microdose and microtracer, nor in the area under curve ratio [14 C]1-OH-midazolam/[14 C]midazolam, showing the PK of midazolam to be linear within the range of the therapeutic and microdoses. CONCLUSION: Our data support the dose linearity of the PK of an IV [14 C]midazolam microdose in children. Hence, a [14 C]midazolam microdosing approach may be used as an alternative to a therapeutic dose of midazolam to study developmental changes in hepatic CYP3A activity in young children.

Pregnancy affects the pharmacokinetics of sildenafil and its metabolite in the rabbit.

1. There is growing interest in the use of sildenafil during pregnancy for various maternal and fetal conditions. This study aims to investigate the effect of pregnancy on the maternal pharmacokinetics (PK) of sildenafil and its main metabolite desmethylsildenafil in rabbits. Using NONMEM, population PK modeling was performed based on plasma samples from 31 rabbits of whom 15 were pregnant and 16 were not. All received a single subcutaneous sildenafil dose of 10 mg/kg. One sample was obtained per rabbit at either 30, 60, 120, 360, 720 or 1320 min after sildenafil administration. 2. A two- and one-compartment PK-model best described the data for sildenafil and desmethylsildenafil, respectively. Compared to non-pregnant rabbits, the central and peripheral volume of distribution and inter-compartmental clearance of sildenafil were lower in pregnant rabbits [32.1 versus 12.2 L, 110 versus 44.4 L and 25.5 versus 12.1 L/h; all p < 0.05]. The formation clearance from sildenafil to desmethylsildenafil was also reduced during pregnancy [13.3 versus 7.8 L/h; p < 0.05]. 3. In contrast, the elimination clearance of desmethylsildenafil, was higher in pregnancy [73.5 versus 116. 9; p < 0.05]. In rabbits, pregnancy impacts PK parameters of sildenafil and its metabolite, leading to an increased peak concentration and 24 h exposure for sildenafil and a decreased 24 h exposure for desmethylsildenafil.

Predicting CYP3A-mediated midazolam metabolism in critically ill neonates, infants, children and adults with inflammation and organ failure.

AIMS: Inflammation and organ failure have been reported to have an impact on cytochrome P450 (CYP) 3A-mediated clearance of midazolam in critically ill children. Our aim was to evaluate a previously developed population pharmacokinetic model both in critically ill children and other populations, in order to allow the model to be used to guide dosing in clinical practice. METHODS: The model was evaluated externally in 136 individuals, including (pre)term neonates, infants, children and adults (body weight 0.77-90 kg, C-reactive protein level 0.1-341 mg l-1 and 0-4 failing organs) using graphical and numerical diagnostics. RESULTS: The pharmacokinetic model predicted midazolam clearance and plasma concentrations without bias in postoperative or critically ill paediatric patients and term neonates [median prediction error (MPE) <30%]. Using the model for extrapolation resulted in well-predicted clearance values in critically ill and healthy adults (MPE <30%), while clearance in preterm neonates was over predicted (MPE >180%). CONCLUSION: The recently published pharmacokinetic model for midazolam, quantifying the influence of maturation, inflammation and organ failure in children, yields unbiased clearance predictions and can therefore be used for dosing instructions in term neonates, children and adults with varying levels of critical illness, including healthy adults, but not for extrapolation to preterm neonates.

Characterization of Intestinal and Hepatic CYP3A-Mediated Metabolism of Midazolam in Children Using a Physiological Population Pharmacokinetic Modelling Approach.

PURPOSE: Changes in drug absorption and first-pass metabolism have been reported throughout the pediatric age range. Our aim is to characterize both intestinal and hepatic CYP3A-mediated metabolism of midazolam in children in order to predict first-pass and systemic metabolism of CYP3A substrates. METHODS: Pharmacokinetic (PK) data of midazolam and 1-OH-midazolam from 264 post-operative children 1-18 years of age after oral administration were analyzed using a physiological population PK modelling approach. In the model, consisting of physiological compartments representing the gastro-intestinal tract and liver,intrinsic intestinal and hepatic clearances were estimated to derive values for bioavailability and plasma clearance. RESULTS: The whole-organ intrinsic clearance in the gut wall and liver were found to increase with body weight, with a 105 (95% confidence interval (CI): 5-405) times lower intrinsic gut wall clearance than the intrinsic hepatic clearance (i.e. 5.08 L/h (relative standard error (RSE) 10%) versus 527 L/h (RSE 7%) for a 16 kg individual, respectively). When expressed per gram of organ, intrinsic clearance increases with increasing body weight in the gut wall, but decreases in the liver, indicating that CYP3A-mediated intrinsic clearance and local bioavailability in the gut wall and liver do not change with age in parallel. The resulting total bioavailability was found to be age-independent with a median of 20.8% in children (95%CI: 3.8-50.0%). CONCLUSION: In conclusion, the intrinsic CYP3A-mediated gut wall clearance is substantially lower than the intrinsic hepatic CYP3A-mediated clearance in children from 1 to 18 years of age, and contributes less to the overall first-pass metabolism compared to adults.

Can Population Modelling Principles be Used to Identify Key PBPK Parameters for Paediatric Clearance Predictions? An Innovative Application of Optimal Design Theory.

PURPOSE: Physiologically-based pharmacokinetic (PBPK) models are essential in drug development, but require parameters that are not always obtainable. We developed a methodology to investigate the feasibility and requirements for precise and accurate estimation of PBPK parameters using population modelling of clinical data and illustrate this for two key PBPK parameters for hepatic metabolic clearance, namely whole liver unbound intrinsic clearance (CLint,u,WL) and hepatic blood flow (Qh) in children. METHODS: First, structural identifiability was enabled through re-parametrization and the definition of essential trial design components. Subsequently, requirements for the trial components to yield precise estimation of the PBPK parameters and their inter-individual variability were established using a novel application of population optimal design theory. Finally, the performance of the proposed trial design was assessed using stochastic simulation and estimation. RESULTS: Precise estimation of CLint,u,WL and Qh and their inter-individual variability was found to require a trial with two drugs, of which one has an extraction ratio (ER) ≤ 0.27 and the other has an ER ≥ 0.93. The proposed clinical trial design was found to lead to precise and accurate parameter estimates and was robust to parameter uncertainty. CONCLUSION: The proposed framework can be applied to other PBPK parameters and facilitate the development of PBPK models.

Amikacin Pharmacokinetics To Optimize Dosing in Neonates with Perinatal Asphyxia Treated with Hypothermia.

Aminoglycoside pharmacokinetics (PK) is expected to change in neonates with perinatal asphyxia treated with therapeutic hypothermia (PATH). Several amikacin dosing guidelines have been proposed for treating neonates with (suspected) septicemia; however, none provide adjustments for cases of PATH. Therefore, we aimed to quantify the differences in amikacin PK between neonates with and without PATH to propose suitable dosing recommendations. Based on amikacin therapeutic drug monitoring data collected retrospectively from neonates with PATH, combined with a published data set, we assessed the impact of PATH on amikacin PK by using population modeling. Monte Carlo and stochastic simulations were performed to establish amikacin exposures in neonates with PATH after dosing according to the current guidelines and according to proposed model-derived dosing guidelines. Amikacin clearance was decreased 40.6% in neonates with PATH, with no changes in volume of distribution. Simulations showed that increasing the dosing interval by 12 h results in a decrease in the percentage of neonates reaching toxic trough levels (>5 mg/liter), from 40 to 76% to 14 to 25%, while still reaching efficacy targets compared to the results of current dosing regimens. Based on this study, a 12-h increase in the amikacin dosing interval in neonates with PATH is proposed to correct for the reduced clearance, yielding safe and effective exposures. As amikacin is renally excreted, further studies into other renally excreted drugs may be required, as their clearance may also be impaired.