Scaling clearance in paediatric pharmacokinetics: All models are wrong, which are useful?

Covariate modeling in pharmacometrics: General points for consideration

Modeling the relationships between covariates and pharmacometric model parameters is a central feature of pharmacometric analyses. The information obtained from covariate modeling may be used for dose selection, dose individualization, or the planning of clinical studies in different population subgroups. The pharmacometric literature has amassed a diverse, complex, and evolving collection of methodologies and interpretive guidance related to covariate modeling. With the number and complexity of technologies increasing, a need for an overview of the state of the art has emerged. In this article the International Society of Pharmacometrics (ISoP) Standards and Best Practices Committee presents perspectives on best practices for planning, executing, reporting, and interpreting covariate analyses to guide pharmacometrics decision making in academic, industry, and regulatory settings.

Spontaneous recovery from neuromuscular block after a single dose of a muscle relaxant in pediatric patients: A systematic review using a network meta-analytic and meta-regression approach

BACKGROUND

Age-related differences in the pharmacokinetics and pharmacodynamics of neuromuscular blocking agents (NMBAs) and the short duration of many surgical procedures put pediatric patients at risk of postoperative residual curarization (PORC). To date, the duration of neuromuscular blocking agent effect in children has not been analyzed in a quantitative review. The current meta-analysis aimed to compare spontaneous recovery following administration of various types and doses of neuromuscular blocking agents and to quantify the effect of prognostic variables associated with the recovery time in pediatric patients.

METHOD

We searched for randomized controlled trials (RCTs) and controlled clinical trials (CCTs) that compared the time to 25% T1 (t25), from 25% to 75% T1 (RI25-75), and to ≥90% train-of-four (tTOF90) neuromuscular recovery between common neuromuscular blocking agent treatments administered as a single bolus to healthy pediatric participants. We compared spontaneous t25, RI25-75, and tTOF90 between (1) neuromuscular blocking agent treatments and (2) age groups receiving a given neuromuscular blocking agent intervention and anesthesia technique. Bayesian random-effects network and pairwise meta-analyses along with meta-regression were used to evaluate the results.

RESULTS

We used data from 71 randomized controlled trials/controlled clinical trials including 4319 participants. Network meta-analysis allowed for the juxtaposition and ranking of spontaneous t25, RI25-75, and tTOF90 following common neuromuscular blocking agent interventions. For all neuromuscular blocking agents a log-linear relationship between dose and duration of action was found. With the neuromuscular blocking agent treatments studied, the average tTOF90 (mean[CrI95]) in children (>2-11 y) was 41.96 [14.35, 69.50] and 17.06 [5.99, 28.30] min shorter than in neonates (<28 d) and infants (28 d-12 M), respectively. We found a negative log-linear correlation between age and duration of neuromuscular blocking agent effect. The difference in the tTOF90 (mean[CrI95]) between children and other age groups increased by 21.66 [8.82, 34.53] min with the use of aminosteroid neuromuscular blocking agents and by 24.73 [7.92, 41.43] min with the addition of sevoflurane/isoflurane for anesthesia maintenance.

CONCLUSIONS

The times to neuromuscular recovery are highly variable. These can decrease significantly with age and are prolonged when volatile anesthetics are administered. This variability, combined with the short duration of many pediatric surgical procedures, makes quantitative neuromuscular monitoring mandatory even after a single dose of neuromuscular blocking agent.

High-dose methotrexate pharmacokinetics and its impact on prognosis of paediatric acute lymphoblastic leukaemia patients: A population pharmacokinetic study

This study delivers a comprehensive evaluation of the efficacy and pharmacokinetics of high-dose methotrexate (HDMTX) in a large cohort of Chinese paediatric acute lymphoblastic leukaemia patients. A total of 533 patients were included in the prognostic analysis. An association was observed between lower steady-state MTX concentrations (<56 μmol/L) and poorer outcomes in intermediate-/high-risk (IR/HR) patients. Subgroup analysis further revealed that this relationship between concentrations and prognosis was even more pronounced in patients with MLL rearrangements. In contrast, such an association did not emerge within the low-risk patient group. Additionally, utilizing population pharmacokinetic modelling (6051 concentrations from 815 patients), we identified the significant impact of physiological maturation, estimated glomerular filtration rate, sex and concurrent dasatinib administration on MTX pharmacokinetics. Simulation-based recommendations include a reduced dosage regimen for those with renal insufficiency and a specific 200 mg/kg dosage for infants under 1 year. The findings underscore the critical role of HDMTX in treating IR/HR populations and call for a reassessment of its application in lower-risk groups. An individualized pharmacokinetic dosage regimen could achieve the most optimal results, ensuring the largest proportion of steady-state concentrations within the optimal range.

Population pharmacokinetics of vancomycin in term neonates with perinatal asphyxia treated with therapeutic hypothermia

AIMS

Little is known about the population pharmacokinetics (PPK) of vancomycin in neonates with perinatal asphyxia treated with therapeutic hypothermia (TH). We aimed to describe the PPK of vancomycin and propose an initial dosing regimen for the first 48 h of treatment with pharmacokinetic/pharmacodynamic target attainment.

METHODS

Neonates with perinatal asphyxia treated with TH were included from birth until Day 6 in a multicentre prospective cohort study. A vancomycin PPK model was constructed using nonlinear mixed-effects modelling. The model was used to evaluate published dosing guidelines with regard to pharmacokinetic/pharmacodynamic target attainment. The area under the curve/minimal inhibitory concentration ratio of 400-600 mg*h/L was used as target range.

RESULTS

Sixteen patients received vancomycin (median gestational age: 41 [range: 38-42] weeks, postnatal age: 4.4 [2.5-5.5] days, birth weight: 3.5 [2.3-4.7] kg), and 112 vancomycin plasma concentrations were available. Most samples (79%) were collected during the rewarming and normothermic phase, as vancomycin was rarely initiated during the hypothermic phase due to its nonempirical use. An allometrically scaled 1-compartment model showed the best fit. Vancomycin clearance was 0.17 L/h, lower than literature values for term neonates of 3.5 kg without perinatal asphyxia (range: 0.20-0.32 L/h). Volume of distribution was similar. Published dosing regimens led to overexposure within 24 h of treatment. A loading dose of 10 mg/kg followed by 24 mg/kg/day in 4 doses resulted in target attainment.

CONCLUSION

Results of this study suggest that vancomycin clearance is reduced in term neonates with perinatal asphyxia treated with TH. Lower dosing regimens should be considered followed by model-informed precision dosing.

Training the next generation of pharmacometric modelers: a multisector perspective

The current demand for pharmacometricians outmatches the supply provided by academic institutions and considerable investments are made to develop the competencies of these scientists on-the-job. Even with the observed increase in academic programs related to pharmacometrics, this need is unlikely to change in the foreseeable future, as the demand and scope of pharmacometrics applications keep expanding. Further, the field of pharmacometrics is changing. The field largely started when Lewis Sheiner and Stuart Beal published their seminal papers on population pharmacokinetics in the late 1970's and early 1980's and has continued to grow in impact and use since its inception. Physiological-based pharmacokinetics and systems pharmacology have grown rapidly in scope and impact in the last decade and machine learning is just on the horizon. While all these methodologies are categorized as pharmacometrics, no one person can be an expert in everything. So how do you train future pharmacometricians? Leading experts in academia, industry, contract research organizations, clinical medicine, and regulatory gave their opinions on how to best train future pharmacometricians. Their opinions were collected and synthesized to create some general recommendations.

Population pharmacokinetics of lisinopril in hypertensive children and adolescents with normal to mildly reduced kidney function

AIMS

Lisinopril, an angiotensin-converting enzyme inhibitor, is a frequently prescribed antihypertensive drug in the paediatric population, while being used off-label under the age of 6 years in the USA and for all paediatric patients globally. The SAFEPEDRUG project (IWT-130033) investigated lisinopril pharmacokinetics in hypertensive paediatric patients corresponding with the day-to-day clinical population.

METHODS

The dose-escalation pilot study included 13 children with primary and secondary hypertension who received oral lisinopril once daily in the morning; doses ranged from 0.05 to 0.2 mg kg-1 . Patients were aged between 1.9 and 17.9 years (median 13.5 years) and weight ranged between 9.62 and 97.2 kg (median 53.2 kg). All data were analysed using Monolix version 2020R1 (Lixoft, France) and R version 3.6.2.

RESULTS

A 1-compartment model with first-order absorption and first-order elimination optimally describes the data. Parameter estimates of absorption rate constant (0.075 h-1 [0.062, 0.088], typical value [95% confidence interval]), volume of distribution (31.38 L 70 kg-1 [10.5, 52.3]) and elimination clearance (24.2 L h-1 70 kg-1 [19.5, 28.9]) show good predictive ability. Significant covariate effects include total body weight on elimination clearance, and distribution volume and estimated glomerular filtration rate (eGFR) on elimination clearance. The effects of eGFR on the elimination clearance are optimally described by a linear effect centred around 105 mL min-1  1.73 m-2 . The effects of body weight were implemented using fixed allometric exponents centred around an adult weight of 70 kg.

CONCLUSION

Lisinopril dose and regimen adjustments for paediatric patients should include eGFR on top of weight adjustments. An expanded model characterizing the pharmacodynamic effect is required to identify the optimal dose and dosing regimen.

Diamorphine pharmacokinetics and conversion factor estimates for intranasal diamorphine in paediatric breakthrough pain:systematic review

BACKGROUND

Intranasal diamorphine is a potential treatment for breakthrough pain but few paediatric data are available to assist dose estimation.

AIM

To determine an intranasal diamorphine dose in children through an understanding of pharmacokinetics.

DESIGN

A systematic review of the literature was undertaken to seek diamorphine pharmacokinetic parameters in neonates, children and adults. Parenteral and enteral diamorphine bioavailability were reviewed with respect to formation of the major metabolite, morphine. Clinical data quantifying equianalgesic effects of diamorphine and morphine were reviewed.

REVIEW SOURCES

PubMed (1960-2020); EMBASE (1980-2020); IPA (1973-2020) and original human research studies that reported diacetylmorphine and metabolite after any dose or route of administration.

RESULTS

The systematic review identified 19 studies: 16 in adults and 1 in children and 2 neonatal reports. Details of study participants were extracted. Age ranged from premature neonates to 67 years and weight 1.4-88 kg. Intranasal diamorphine bioavailability was predicted as 50%. The equianalgesic intravenous conversion ratio of morphine:diamorphine was 2:1. There was heterogeneity between pharmacokinetic parameter estimates attributed to routes of administration, lack of size standardisation, methodology and pharmacokinetic analysis. Estimates of the pharmacokinetic parameters clearance and volume of distribution were reduced in neonates. There were insufficient paediatric data to characterise clearance or volume maturation of either diamorphine or its metabolites.

CONCLUSIONS

We estimate equianalgesic ratios of intravenous morphine:diamorphine 2:1, intravenous morphine:intranasal diamorphine 1:1 and oral morphine:intranasal diamorphine of 1:3. These ratios are based on adult literature, but are reasonable for deciding on an initial dose of 0.1 mg/kg in children 4-13 years.

Pharmacokinetics and safety of prolonged paracetamol treatment in neonates: An interventional cohort study

AIMS

To investigate the pharmacokinetics and safety of prolonged paracetamol use (>72 h) for neonatal pain.

METHODS

Neonates were included if they received paracetamol orally or intravenously for pain treatment. A total of 126 samples were collected. Alanine aminotransferase and bilirubin were measured as surrogate liver safety markers. Paracetamol and metabolites were measured in plasma. Pharmacokinetic parameters for the parent compound were estimated with a nonlinear mixed-effects model.

RESULTS

Forty-eight neonates were enrolled (38 received paracetamol for >72 h). Median gestational age was 38 weeks (range 25-42), and bodyweight at inclusion was 2954 g (range 713-4750). Neonates received 16 doses (range 4-55) over 4.1 days (range 1-13.8). The median (range) dose was 10.1 mg/kg (2.9-20.3). The median oxidative metabolite concentration was 14.6 μmol/L (range 0.12-113.5) and measurable >30 h after dose. There was no significant difference (P > .05) between alanine aminotransferase and bilirubin measures at <72 h or >72 h of paracetamol treatment or the start and end of the study. Volume of distribution and paracetamol clearance for a 2.81-kg neonate were 2.99 L (% residual standard error = 8, 95% confidence interval 2.44-3.55) and 0.497 L/h (% residual standard error = 7, 95% confidence interval 0.425-0.570), respectively. Median steady-state concentration from the parent model was 50.3 μmol/L (range 30.6-92.5), and the half-life was 3.55 h (range 2.41-5.65).

CONCLUSION

Our study did not provide evidence of paracetamol-induced liver injury nor changes in metabolism in prolonged paracetamol administration in neonates.

Predicting Maternal and Infant Tetrahydrocannabinol Exposure in Lactating Cannabis Users: A Physiologically Based Pharmacokinetic Modeling Approach

A knowledge gap exists in infant tetrahydrocannabinol (THC) data to guide breastfeeding recommendations for mothers who use cannabis. In the present study, a paired lactation and infant physiologically based pharmacokinetic (PBPK) model was developed and verified. The verified model was used to simulate one hundred virtual lactating mothers (mean age: 28 years, body weight: 78 kg) who smoked 0.32 g of cannabis containing 14.14% THC, either once or multiple times. The simulated breastfeeding conditions included one-hour post smoking and subsequently every three hours. The mean peak concentration (Cmax) and area under the concentration-time curve (AUC(0-24 h)) for breastmilk were higher than in plasma (Cmax: 155 vs. 69.9 ng/mL; AUC(0-24 h): 924.9 vs. 273.4 ng·hr/mL) with a milk-to-plasma AUC ratio of 3.3. The predicted relative infant dose ranged from 0.34% to 0.88% for infants consuming THC-containing breastmilk between birth and 12 months. However, the mother-to-infant plasma AUC(0-24 h) ratio increased up to three-fold (3.4-3.6) with increased maternal cannabis smoking up to six times. Our study demonstrated the successful development and application of a lactation and infant PBPK model for exploring THC exposure in infants, and the results can potentially inform breastfeeding recommendations.

Pharmacokinetic pharmacodynamic modeling of analgesics and sedatives in children

Pharmacokinetic pharmacodynamic modeling is an important tool which uses statistical methodology to provide a better understanding of the relationship between concentration and effect of drugs such as analgesics and sedatives. Pharmacokinetic pharmacodynamic models also describe between-subject variability that allows identification of subgroups and dose adjustment for optimal pain management in individual patients. This approach is particularly useful in the pediatric population, where most drugs have received limited evaluation and dosing is extrapolated from adult practice. In children, the covariates of weight and age are used to describe size- and maturation-related changes in pharmacokinetics. It is important to consider both size and maturation in order to develop an accurate model and determine the optimal dose for different age groups. An adequate assessment of analgesic and sedative effect using pain scales or brain activity measures is essential to build reliable pharmacokinetic pharmacodynamic models. This is often challenging in children due to the multidimensional nature of pain and the limited sensitivity and specificity of some measurement tools. This review provides a summary of the pharmacokinetic and pharmacodynamic methodology used to describe the dose-concentration-effect relationship of analgesics and sedation in children, with a focus on the different pharmacodynamic endpoints and the challenges of pharmacodynamic modeling.

The Neonatal and Paediatric Pharmacokinetics of Antimicrobials study (NAPPA): investigating amoxicillin, benzylpenicillin, flucloxacillin and piperacillin pharmacokinetics from birth to adolescence

BACKGROUND

Pharmacokinetic (PK) data underlying paediatric penicillin dosing remain limited, especially in critical care.

OBJECTIVES

The primary objective of the Neonatal and Paediatric Pharmacokinetics of Antimicrobials study (NAPPA) was to characterize PK profiles of commonly used penicillins using data obtained during routine care, to further understanding of PK variability and inform future evidence-based dosing.

METHODS

NAPPA was a multicentre study of amoxicillin, co-amoxiclav, benzylpenicillin, flucloxacillin and piperacillin/tazobactam. Patients were recruited with informed consent. Antibiotic dosing followed standard of care. PK samples were obtained opportunistically or at optimal times, frozen and analysed using UPLC with tandem MS. Pharmacometric analysis was undertaken using NONMEM software (v7.3). Model-based simulations (n = 10 000) tested PTA with British National Formulary for Children (BNFC) and WHO dosing. The study had ethical approval.

RESULTS

For the combined IV PK model, 963 PK samples from 370 participants were analysed simultaneously incorporating amoxicillin, benzylpenicillin, flucloxacillin and piperacillin data. BNFC high-dose regimen simulations gave these PTA results (median fT>MIC at breakpoints of specified pathogens): amoxicillin 100% (Streptococcus pneumoniae); benzylpenicillin 100% (Group B Streptococcus); flucloxacillin 48% (MSSA); and piperacillin 100% (Pseudomonas aeruginosa). Oral population PK models for flucloxacillin and amoxicillin enabled estimation of first-order absorption rate constants (1.16 h-1 and 1.3 h-1) and bioavailability terms (62.7% and 58.7%, respectively).

CONCLUSIONS

NAPPA represents, to our knowledge, the largest prospective combined paediatric penicillin PK study undertaken to date, and the first paediatric flucloxacillin oral PK model. The PTA results provide evidence supportive of BNFC high-dose IV regimens for amoxicillin, benzylpenicillin and piperacillin.

Vancomycin population pharmacokinetics analysis in Chinese paediatric patients with varying degrees of renal function and ages: development of new practical dosing recommendations

OBJECTIVES

To describe the pharmacokinetics of vancomycin in a large Chinese paediatric cohort with varying degrees of renal function and ages and to develop practical dosing guidelines.

PATIENTS AND METHODS

We conducted a retrospective population pharmacokinetic study using data from paediatric patients who received vancomycin between June 2013 and June 2022. A non-linear mixed-effect modelling approach with a one-compartment model structure was applied. Monte Carlo simulations were used to stimulate an optimal dosage regimen to achieve the target of AUC24/MIC between 400 and 650.

RESULTS

We analysed a total of 673 paediatric patients and 1547 vancomycin serum concentrations. Covariate analysis revealed that physiological maturation, renal function, albumin and cardiothoracic surgery (CTS) significantly affected vancomycin pharmacokinetics. The typical clearance and volume of distribution, standardized to 70 kg, were 7.75 L/h (2.3% relative standard error, RSE) and 36.2 L (1.7% RSE), respectively. Based on the model, we proposed an optimal dosing regimen that considers the patient's age and estimate glomerular filtration rate (eGFR) to achieve a target AUC24/MIC for CTS and non-CTS patients. We also found that a loading dose of 20 mg/kg can help patients with an eGFR of <60 mL/min/1.73 m2 achieve the target AUC on the first day of treatment.

CONCLUSIONS

We established vancomycin pharmacokinetic parameters in Chinese paediatric patients and proposed a dosing guideline integrating eGFR, age and CTS status, potentially improving clinical outcomes and reducing nephrotoxicity risk.

Population pharmacokinetics of benznidazole in neonates, infants and children using a new pediatric formulation

BACKGROUND

There is a major need for information on pharmacokinetics (PK) of benznidazole (BNZ) in children with Chagas disease (CD). We conducted a multicentre population PK, safety and efficacy study in children, infants and neonates with CD treated with BNZ (formulated in 100 mg tablets or 12.5 mg dispersible tablets, developed by the pharmaceutical company LAFEPE, in a collaboration with DNDi).

METHODS

81 children 0-12 years old were enrolled at 5 pediatric centers in Argentina. Diagnosis of T. cruzi infection was confirmed by direct microscopic examination, or at least two positive conventional serological tests. Subject enrolment was stratified by age: newborns to 2 years (minimum of 10 newborns) and >2-12 years. BNZ 7.5 mg/kg/d was administered in two daily doses for 60 days. Five blood samples per child were obtained at random times within pre-defined time windows at Day 0 at 2-5 h post-dose; during steady state, one sample at Day 7 and at Day 30; and two samples at 12-24 h after final BNZ dose at Day 60. The primary efficacy endpoint was parasitological clearance by qualitative PCR at the end of treatment.

RESULTS

Forty-one (51%) patients were under 2 years of age (including 14 newborns <1 month of age). Median age at enrolment was 22 months (mean: 43.2; interquartile range (IQR) 7-72 months). The median measured BNZ Cmax was 8.32 mg/L (IQR 5.95-11.8; range 1.79-19.38). Median observed BNZ Cmin (trough) concentration was 2 mg/L (IQR 1.25-3.77; range 0.14-7.08). Overall median simulated Css was 6.3 mg/L (IQR 4.7-8.5 mg/L). CL/F increased quickly during the first month of postnatal life and reached adult levels after approximately 10 years of age. Negative qPCR was observed at the end of treatment in all 76 patients who completed the treatment. Five patients discontinued treatment (3 due to AEs and 2 due to lack of compliance).

CONCLUSION

We observed lower BNZ plasma concentrations in infants and children than those previously reported in adults treated with comparable mg/kg doses. Despite these lower concentrations, pediatric treatment was well tolerated and universally effective, with a high response rate and infrequent, mild AEs.

TRIAL REGISTRATION

Registered in clinicaltrials.gov #NCT01549236.

Population Pharmacokinetics and Dosing Optimization of Ceftazidime in Term Asphyxiated Neonates during Controlled Therapeutic Hypothermia

Ceftazidime is an antibiotic commonly used to treat bacterial infections in term neonates undergoing controlled therapeutic hypothermia (TH) for hypoxic-ischemic encephalopathy after perinatal asphyxia. We aimed to describe the population pharmacokinetics (PK) of ceftazidime in asphyxiated neonates during hypothermia, rewarming, and normothermia and propose a population-based rational dosing regimen with optimal PK/pharmacodynamic (PD) target attainment. Data were collected in the PharmaCool prospective observational multicenter study. A population PK model was constructed, and the probability of target attainment (PTA) was assessed during all phases of controlled TH using targets of 100% of the time that the concentration in the blood exceeds the MIC (T>MIC) (for efficacy purposes and 100% T>4×MIC and 100% T>5×MIC to prevent resistance). A total of 35 patients with 338 ceftazidime concentrations were included. An allometrically scaled one-compartment model with postnatal age and body temperature as covariates on clearance was constructed. For a typical patient receiving the current dose of 100 mg/kg of body weight/day in 2 doses and assuming a worst-case MIC of 8 mg/L for Pseudomonas aeruginosa, the PTA was 99.7% for 100% T>MIC during hypothermia (33.7°C; postnatal age [PNA] of 2 days). The PTA decreased to 87.7% for 100% T>MIC during normothermia (36.7°C; PNA of 5 days). Therefore, a dosing regimen of 100 mg/kg/day in 2 doses during hypothermia and rewarming and 150 mg/kg/day in 3 doses during the following normothermic phase is advised. Higher-dosing regimens (150 mg/kg/day in 3 doses during hypothermia and 200 mg/kg/day in 4 doses during normothermia) could be considered when achievements of 100% T>4×MIC and 100% T>5×MIC are desired.

Allometric Scaling of Testosterone Enanthate Pharmacokinetics to Adolescent Hypogonadal Males (IM and SC Administration)

Context

Intramuscular (IM) testosterone enanthate (TE) and testosterone pellets were US Food and Drug Administration approved before 1962 for pediatric use but not studied in controlled trials in adolescents.

Objective

An analysis using nonlinear mixed effect (NLME) modeling was designed to evaluate the adult pharmacokinetics (PK) of subcutaneous (SC) and IM TE. This model was used to simulate SC and IM TE administration in adolescents of different weight groups.

Methods

Data from adult male patients in a phase 2 trial were used to characterize the PK of TE using population PK modeling for SC and IM administration: Allometry was used to scale PK parameters from the adult model to simulate adolescent (aged 12 to < 18 years) serum testosterone levels at body weights of 30, 40, 50, and 60 kg after weekly, every-other-week (EOW), and monthly SC and IM administration of 12.5, 25, 50, 75, and 100 mg TE regimens.

Results

The final data set included 714 samples from 15 patients receiving 100 mg SC TE and 123 samples from 10 patients receiving 200 mg IM TE. In simulated populations, average serum concentration SC:IM ratios were 0.783, 0.776, and 0.757 at steady state for weekly, EOW, and monthly dosing groups, respectively. Simulated regimens of 12.5 mg SC TE monthly produced serum testosterone levels representative of early puberty and simulated pubertal stage progression following multiple subsequent testosterone dose increases.

Conclusion

SC TE administration achieved a testosterone exposure-response relationship similar to IM TE in simulated adolescent hypogonadal males, which may reduce size of fluctuations in serum T and related symptoms.

Evaluation of treosulfan cumulative exposure in paediatric patients through population pharmacokinetics and dosing simulations

AIM

To investigate the pharmacokinetics (PK) of intravenous treosulfan in paediatric patients undergoing haematopoietic stem cell transplantation (HSCT) for a broad range of diseases and to explore the impact of different dosing regimens on treosulfan exposure (area under the concentration-time curve, AUC_0→∞ ) through dosing simulations.

METHODS

A prospective multicentre PK study was conducted using treosulfan concentration data (n = 423) collected from 53 children (median age 3.5, range 0.2-17.0 years) receiving three daily age-guided doses (10-14 g/m² ). Population PK modelling was performed using NONMEM software, utilising a stepwise forward selection backward elimination method and likelihood-ratio test for screening covariates to describe PK variability. Monte Carlo simulation was used to generate patient PK data for 10 000 virtual paediatric patients and cumulative AUC_0→∞ values were evaluated using age, body surface area (BSA) and model-based dosing regimens, targeting 4800 mg*h/L.

RESULTS

Treosulfan concentration data were described using a one-compartment PK model with first-order elimination. Population mean (95% CI) estimates for clearance (CL) and volume of distribution (V) were 16.3 (14.9-18.1) L/h and 41.9 (38.8-45.1) L, respectively. Allometrically scaled body weight was the best covariate descriptor for CL and V, and maturational age further explained variability in CL. Dosing simulations indicated that in young patient groups (<2 years), a model-based dosing regimen more accurately achieved the target AUC_0→∞ (58.3%) over the age (42.6%) and BSA-based (51.3%) regimens.

CONCLUSION

Treosulfan disposition was described through allometric body weight and maturational age descriptors. Model-informed dosing is recommended for patients under 2 years. Treosulfan PK parameters and AUC_0→∞ were not influenced by patient disease.

Optimising intravenous salbutamol in children: a phase 2 study

OBJECTIVE

The β2-agonists such as salbutamol are the mainstay of asthma management. Pharmacokinetic-pharmacodynamic (PKPD) models to guide paediatric dosing are lacking. We explored the relationship between salbutamol dose, serum concentration, effectiveness and adverse effects in children by developing a PKPD model.

DESIGN

A prospective cohort study of children admitted to hospital with acute asthma, who received intravenous salbutamol.

SETTING

Children were recruited in two cohorts: the emergency departments of two London hospitals or those retrieved by the Children's Acute Transport Service to three London paediatric intensive care units.

PATIENTS

Patients were eligible if aged 1-15 years, admitted for acute asthma and about to receive or receiving intravenous salbutamol.

INTERVENTIONS

Treatment was according to local policy. Serial salbutamol plasma levels were taken. Effectiveness measurements were recorded using the Paediatric Asthma Severity Score (PASS). Toxicity measurements included lactate, pH, glucose, heart rate, blood pressure and arrhythmias. PKPD modelling was performed with non-linear mixed-effect models.

MAIN OUTCOMES

Fifty-eight children were recruited with 221 salbutamol concentration measurements from 54 children. Median (range) age was 2.9 (1.1-15.2) years, and weight was 13.6 (8-57.3) kg. Ninety-five PASS measurements and 2078 toxicity measurements were obtained.

RESULTS

A two-compartment PK model adequately described the time course of salbutamol-plasma concentrations. An E_MAX (maximum drug effect) concentration-effect relationship described PASS and toxicity measures. PKPD simulations showed an infusion of 0.5 µg/kg/min (maximum 20 µg/min) for 4 hours after bolus achieves >90% maximal bronchodilation for 12 hours.

CONCLUSIONS

A paediatric PKPD model for salbutamol is described. An infusion of 0.5 µg/kg/min after bolus achieves effective bronchodilation. Higher rates are associated with greater tachycardia and hyperglycaemia.

Population pharmacokinetics of sublingually administered tacrolimus in infants and young children with liver transplantation

AIMS

Pharmacokinetics of tacrolimus after sublingual administration is not characterized in paediatric liver transplant patients. Therefore, we aimed to develop a population pharmacokinetic model of sublingually administered tacrolimus in patients who cannot swallow the capsules due to their age, sedation status and/or mechanical ventilation during the first weeks post-transplantation.

METHODS

Demographic, clinical and pharmacological variables, including tacrolimus whole blood concentrations obtained from therapeutic drug monitoring and data from dense-sampling pharmacokinetic profiles, were recorded in 26 paediatric patients with biliary atresia who underwent liver transplantation between 2016 and 2021. Population pharmacokinetic analysis was performed with NONMEM v7.4.

RESULTS

Disposition of tacrolimus was best characterized by a 2-compartment model with clearance achieving half of the maximum elimination capacity (CL_MAX  = 4.1 L/h) at 4.6 days post-transplantation (T₅₀ ). Compared to sedated patients, nonsedated status showed an increased first-order absorption rate constant (1.1 vs. 0.1 h^-1 ) and a 24% reduction in bioavailability (F_NS ) at 14 days post-transplant. The model was able to explain the oral absorption pattern in nonsedated patients as the result of gut bioavailability (0.9) and hepatic extraction ratio, with the latter being responsible for first-pass effects. Estimates of interindividual variability remained moderate (25.9% for the gut bioavailability) to high (79.8% for the apparent volume of distribution of the central compartment, and 101% for T₅₀ ).

CONCLUSION

A population pharmacokinetic model of sublingually administered tacrolimus in paediatric patients was developed to characterize different absorption mechanisms. Once the model is externally validated, the effect of post-transplant time on clearance and the sedation status may be considered in routine dosing management.

Integrative population pharmacokinetic/pharmacodynamic analysis of nemonoxacin capsule in Chinese patients with community-acquired pneumonia

Introduction: Nemonoxacin is an innovative quinolone antibiotic for treatment of community-acquired pneumonia (CAP). As more data are available from clinical studies, it is necessary to perform an integrative pharmacokinetic/pharmacodynamic (PK/PD) analysis to support and justify the optimal dosing regimen of nemonoxacin in clinical practice. Methods and Results: We developed a population PK model using non-linear mixed effect model based on the data of 195 Chinese subjects receiving nemonoxacin in phase I to III clinical trials. The base model was a standard two-compartment PK model defined by clearance (12 L/h) and central volume of distribution (86 L). Covariates included creatinine clearance (CL_cr), body weight (BW), sex, disease status and food. Compared to the subject with BW 60 kg, C_max and A U C 0 - 24 ,   ss reduced by 24% and 19% in the subject with BW 80 kg, respectively. Compared to the subject with CL_cr 150 ml/min, A U C 0 - 24 ,   ss and T_1/2 increased by 28% and 24%, respectively in the subject with CL_cr 30 ml/min. Compared to the fasted status, T_max of nemonoxacin increased by 1.2 h in the subject with fed status. Effects of sex and disease status on PK parameters were small (change of PK parameters ≤19%). AUC_0-24/MIC and %T > MIC were identified as the optimal PK/PD indices for predicting clinical efficacy. The AUC_0-24/MIC target was 63.3, 97.8, and 115.7 against Streptococcus pneumoniae, Staphylococcus aureus, and Haemophilus influenzae, respectively. The %T > MIC target was 7.96% against Klebsiella pneumoniae. Monte Carlo simulation showed that treatment with nemonoxacin 500 mg q24 h could attain a PK/PD cutoff value higher than the MIC₉₀ against S. pneumoniae and S. aureus. The corresponding cumulative fraction of response (CFR) was greater than 93%, while nemonoxacin 750 mg q24 h would provide higher PK/PD cutoff value against Haemophilus parainfluenzae, and higher CFR (83%) than 500 mg q24 h. Conclusion: Integrative PK/PD analysis justifies the reliable clinical and microbiological efficacy of nemonoxacin 500 mg q24 h in treating CAP caused by S. pneumoniae, S. aureus, and K. pneumoniae, irrespective of patient sex, mild renal impairment, empty stomach or not. However, nemonoxacin 750 mg q24 h would provide better efficacy than 500 mg q24 h for the CAP caused by H. parainfluenzae in terms of CFR.

External evaluation of published population pharmacokinetic models of posaconazole

Population pharmacokinetic (PopPK) models of posaconazole have been established to promote the precision dosing. However, the performance of these models extrapolated to other centers has not been evaluated. This study aimed to conduct an external evaluation of published posaconazole PopPK models to evaluate their predictive performance. Posaconazole PopPK models screened from the PubMed and MEDLINE databases were evaluated using an external dataset of 213 trough concentration samples collected from 97 patients. Their predictive performance was evaluated by prediction-based diagnosis (prediction error), simulation-based diagnosis (visual predictive check), and Bayesian forecasting. In addition, external cohorts with and without proton pump inhibitor were used to evaluate the models respectively. Ten models suitable for the external dataset were finally included into the study. In prediction-based diagnostics, none of the models met pre-determined criteria for predictive indexes. Only M4, M6, and M10 demonstrated favorable simulations in visual predictive check. The prediction performance of M5, M7, M8, and M9 evaluated using the cohort without proton pump inhibitor showed a significant improvement compared to that evaluated using the whole cohort. Consistent with our expectations, Bayesian forecasting significantly improved the predictive per-formance of the models with two or three prior observations. In general, the applicability of these published posaconazole PopPK models extrapolated to our center was unsatisfactory. Prospective studies combined with therapeutic drug monitoring are needed to establish a PopPK model for posaconazole in the Chinese population to promote individualized dosing.

Population pharmacokinetics of posaconazole in Chinese pediatric patients with acute leukaemia: effect of food on bioavailability and dose optimization

This study aimed to investigate the effect of food on the pharmacokinetics of posaconazole suspension in pediatric patients with acute leukaemia and to recommend optimal dosing strategies. This single-site, prospective, open-label, observational study was conducted in 42 patients and included 186 plasma concentrations of posaconazole. Sparse data were analyzed using population pharmacokinetic modeling. Monte Carlo simulations were conducted to predict the morning trough concentrations at steady-state with the proposed dose of 2-7 mg/kg three times daily (tid) or four times daily (qid) for bodyweights of 10-36 kg. The target concentrations were 700 ng/mL for prophylaxis and 1000 ng/mL for treatment. Dosage regimens with percentage of target attainment (PTA) ≥70% were recommended. A one-compartment model with allometric scaling adequately described the pharmacokinetic profile. The apparent clearance was 9.05 L/h (95% confidence interval [CI] 7.14-11.09) and the apparent volume of distribution was 283 L (95% CI 168-491) for a typical individual of 17.5 kg. The relative bioavailability with high-fat diet was as high as 1.95-fold compared with regular food. Following the intake of regular meals, 4 mg/kg qid was adequate with a PTA ≥ 71.8% for prophylaxis. A dosage of 6 mg/kg qid under a regular diet reached a PTA ≥ 73.4% for treatment. The recommended dosage of posaconazole for prophylaxis and treatment could be predicted by the pharmacokinetic model based on bodyweight and diet type in pediatric patients.

Population Pharmacokinetics of Intranasal Dexmedetomidine in Infants and Young Children

BACKGROUND

Intranasal dexmedetomidine provides noninvasive, effective procedural sedation for pediatric patients, and has been widely used in clinical practice. However, the dosage applied has varied fourfold in pediatric clinical studies. To validate an appropriate dosing regimen, this study investigated the pharmacokinetics of intranasal dexmedetomidine in Chinese children under 3 yr old.

METHODS

Intranasal dexmedetomidine 2 µg · kg-1 was administered to children with simple vascular malformations undergoing interventional radiological procedures. A population pharmacokinetic analysis with data from an optimized sparse-sampling design was performed using nonlinear mixed-effects modeling. Clearance was modeled using allometric scaling and a sigmoid postmenstrual age maturation model. Monte Carlo simulations were performed to assess the different dosing regimens.

RESULTS

A total of 586 samples from 137 children aged 3 to 36 months were included in the trial. The data were adequately described by a two-compartment model with first-order elimination. Body weight with allometric scaling and maturation function were significant covariates of dexmedetomidine clearance. The pharmacokinetic parameters for the median subjects (weight 10 kg and postmenstrual age 101 weeks) in the authors' study were apparent central volume of distribution 7.55 l, apparent clearance of central compartment 9.92 l · h-1, apparent peripheral volume of distribution 7.80 l, and apparent intercompartmental clearance 61.7 l · h-1. The simulation indicated that at the dose of 2 µg · kg-1, 95% of simulated individuals could achieve a target therapeutic concentration of 0.3 ng · ml-1 within 20 min, and the average peak concentration of 0.563 ng · ml-1 could be attained at 61 min.

CONCLUSIONS

The pharmacokinetic characteristics of intranasal dexmedetomidine were evaluated in Chinese pediatric patients aged between 3 and 36 months. An evidence-based dosing regimen at 2 µg · kg-1 could achieve a preset therapeutic threshold of mild to moderate sedation that lasted for up to 2 h.

EDITOR’S PERSPECTIVE

Population pharmacokinetics of ethambutol in African children: a pooled analysis

OBJECTIVES

Ethambutol protects against the development of resistance to co-administered drugs in the intensive phase of first-line anti-TB treatment in children. It is especially relevant in settings with a high prevalence of HIV or isoniazid resistance. We describe the population pharmacokinetics of ethambutol in children with TB to guide dosing in this population.

METHODS

We pooled data from 188 intensively sampled children from the DATiC, DNDi and SHINE studies, who received 15-25 mg/kg ethambutol daily according to WHO guidelines. The median (range) age and weight of the cohort were 1.9 (0.3-12.6) years and 9.6 (3.9-34.5) kg, respectively. Children with HIV (HIV+; n = 103) received ART (lopinavir/ritonavir in 92%).

RESULTS

Ethambutol pharmacokinetics were best described by a two-compartment model with first-order elimination and absorption transit compartments. Clearance was estimated to reach 50% of its mature value by 2 months after birth and 99% by 3 years. Typical steady-state apparent clearance in a 10 kg child was 15.9 L/h. In HIV+ children on lopinavir/ritonavir, bioavailability was reduced by 32% [median (IQR) steady-state Cmax = 0.882 (0.669-1.28) versus 1.66 (1.21-2.15) mg/L). In young children, bioavailability correlated with age. At birth, bioavailability was 73.1% of that in children 3.16 years or older.

CONCLUSIONS

To obtain exposure within the 2-6 mg/L recommended range for Cmax, the current doses must be doubled (or tripled with HIV+ children on lopinavir/ritonavir) for paediatric patients. This raises concerns regarding the potential for ocular toxicity, which would require evaluation.

Pharmacometrics in tuberculosis: progress and opportunities

Tuberculosis remains one of the leading causes of death by a communicable agent, infecting up to one-quarter of the world's population, predominantly in disadvantaged communities. Pharmacometrics employs quantitative mathematical models to describe the relationships between pharmacokinetics and pharmacodynamics, and to predict drug doses, exposures, and responses. Pharmacometric approaches have provided a scientific basis for improved dosing of antituberculosis drugs and concomitantly administered antiretrovirals at the population level. The development of modelling frameworks including physiologically-based pharmacokinetics, quantitative systems pharmacology and machine learning provides an opportunity to extend the role of pharmacometrics to in silico quantification of drug-drug interactions, prediction of doses for special populations, dose optimization and individualization, and understanding the complex exposure-response relationships of multidrug regimens in terms of both efficacy and safety, informing regimen design for future study. In this short clinically-focused review, we explore what has been done, and what opportunities exist for pharmacometrics to impact tuberculosis pharmacotherapy.

Pharmacokinetic modeling and simulation to understand diamorphine dose-response in neonates, children, and adolescents

Pharmacokinetic-pharmacodynamic modeling and simulation can facilitate understanding and prediction of exposure-response relationships in children with acute or chronic pain. The pharmacokinetics of diamorphine (diacetylmorphine, heroin), a strong opioid, remain poorly quantified in children and dose is often guided by clinical acumen. This tutorial demonstrates how a model to describe intranasal and intravenous diamorphine pharmacokinetics can be fashioned from a model for diamorphine disposition in adults and a model describing morphine disposition in children. Allometric scaling and maturation models were applied to clearances and volumes to account for differences in size and age between children and adults. The utility of modeling and simulation to gain insight into the analgesic exposure-response relationship is demonstrated. The model explains reported observations, can be used for interrogation, interpolated to determine equianalgesia and inform future clinical studies. Simulation was used to illustrate how diamorphine is rapidly metabolized to morphine via its active metabolite 6-monoacetylmorphine, which mediates an early dopaminergic response accountable for early euphoria. Morphine formation is then responsible for the slower, prolonged analgesic response. Time-concentration profiles of diamorphine and its metabolites reflected disposition changes with age and were used to describe intravenous and intranasal dosing regimens. These indicated that morphine exposure in children after intranasal diamorphine 0.1 mg.kg-1 was similar to that after intranasal diamorphine 5 mg in adults. A target concentration of morphine 30 μg.L-1 can be achieved by a diamorphine intravenous infusion in neonates 14 μg.kg-1 .h-1 , in a 5-year-old child 42 μg.kg-1 .h-1 and in an 15 year-old-adolescent 33 μg.kg-1 .h-1 .

A General Biphasic Bodyweight Model for Scaling Basal Metabolic Rate, Glomerular Filtration Rate, and Drug Clearance from Birth to Adulthood

The objective of this study is to propose a unified, continuous, and bodyweight-only equation to quantify the changes of human basal metabolic rate (BMR), glomerular filtration rate (GFR), and drug clearance (CL) from infancy to adulthood. The BMR datasets were retrieved from a comprehensive historical database of male and female subjects (0.02 to 64 years). The CL datasets for 17 drugs and the GFR dataset were generated from published maturation and growth models with reported parameter values. A statistical approach was used to simulate the model-generated CL and GFR data for a hypothetical population with 26 age groups (from 0 to 20 years). A biphasic equation with two power-law functions of bodyweight was proposed and evaluated as a general model using nonlinear regression and dimensionless analysis. All datasets universally reveal biphasic curves with two distinct linear segments on log-log plots. The biphasic equation consists of two reciprocal allometric terms that asymptotically determine the overall curvature. The fitting results show a superlinear scaling phase (asymptotic exponent >1; ca. 1.5-3.5) and a sublinear scaling phase (asymptotic exponent <1; ca. 0.5-0.7), which are separated at the phase transition bodyweight ranging from 5 to 20 kg with a mean value of 10 kg (corresponding to 1 year of age). The dimensionless analysis generalizes and offers quantitative realization of the maturation and growth process. In conclusion, the proposed mixed-allometry equation is a generic model that quantitatively describes the phase transition in the human maturation process of diverse human functions.

Allometric Scaling in Pharmacokinetic Studies in Anesthesiology

A clinical review is presented of basic allometric scaling theory and its application to pharmacokinetic models in anesthesia and other fields in the biologic sciences.

An update on the use of allometric and other scaling methods to scale drug clearance in children: towards decision tables

INTRODUCTION

When pediatric data are not available for a drug, allometric and other methods are applied to scale drug clearance across the pediatric age-range from adult values. This is applied when designing first-in-child studies, but also for off-label drug prescription.

AREAS COVERED

This review provides an overview of the systematic accuracy of allometric and other pediatric clearance scaling methods compared to gold-standard PBPK predictions. The findings are summarized in decision tables to provide a priori guidance on the selection of appropriate pediatric clearance scaling methods for both novel drugs for which no pediatric data are available and existing drugs in clinical practice.

EXPERT OPINION

While allometric scaling principles are commonly used to scale pediatric clearance, there is no universal allometric exponent (i.e., 1, 0.75 or 0.67) that can accurately scale clearance for all drugs from adults to children of all ages. Therefore, pediatric scaling decision tables based on age, drug elimination route, binding plasma protein, fraction unbound, extraction ratio, and/or isoenzyme maturation are proposed to a priori select the appropriate (allometric) clearance scaling method, thereby reducing the need for full PBPK-based clearance predictions. Guidance on allometric scaling when estimating pediatric clearance values is provided as well.

Pediatric pharmacokinetic considerations and implications for drug dosing

Optimizing the dosing of medicines for pediatric patients in routine clinical practice and determining the dose for clinical trials is still a challenging task. Children differ from adults in their response to drugs due to inherent differences in pharmacokinetics and/or pharmacodynamics, and responses may also vary among pediatric patients of different ages. However, the greatest disparities compared to adult pharmacokinetic profiles are observed in children below 2 years of age. The maturation of the liver and the kidneys, as well as the variation in body composition, are considered to be the main sources of pharmacokinetic variability. Hence, besides specific pharmacodynamic features, understanding age-related changes in drug absorption, distribution, and elimination is fundamental for optimizing drug efficacy and avoiding toxicity. This paper summarizes the pharmacokinetic changes throughout the childhood, along with the effect of developmental changes on drug dosage calculation. In clinical practice, age and body weight-based dosing regimens are usually used. In spite of dosing recommendations based on age and/or body weight, variabilities in pharmacokinetics and pharmacodynamic response remain, implying a need to monitor patients and optimize the dosing regimen according to physiological characteristics, disease characteristics and therapy.

Wide size dispersion and use of body composition and maturation improves the reliability of allometric exponent estimates

To evaluate study designs and the influence of dispersion of body size, body composition and maturation of clearance or reliable estimation of allometric exponents. Non-linear mixed effects modeling and parametric bootstrap were employed to assess how the study sample size, number of observations per subject, between subject variability (BSV) and dispersion of size distribution affected estimation bias and uncertainty of allometric exponents. The role of covariate model misspecification was investigated using a large data set ranging from neonates to adults. A decrease in study sample size, number of observations per subject, an increase in BSV and a decrease in dispersion of size distribution, increased the uncertainty of allometric exponent estimates. Studies conducted only in adults with drugs exhibiting normal (30%) BSV in clearance may need to include at least 1000 subjects to be able to distinguish between allometric exponents of 2/3 and 1. Nevertheless, studies including both children and adults can distinguish these exponents with only 100 subjects. A marked bias of 45% (95%CI 41-49%) in the estimate of the allometric exponent of clearance was obtained when maturation and body composition were ignored in infants. A wide dispersion of body size (e.g. infants, children and adults) is required to reliably estimate allometric exponents. Ignoring differences in body composition and maturation of clearance may bias the exponent for clearance. Therefore, pharmacometricians should avoid estimating allometric exponent parameters without suitable designs and covariate models. Instead, they are encouraged to rely on the well-developed theory and evidence that clearance and volume parameters in humans scale with theory-based exponents.

Development of weight and age-based dosing of daily primaquine for radical cure of vivax malaria

BACKGROUND

In many endemic areas, Plasmodium vivax malaria is predominantly a disease of young adults and children. International recommendations for radical cure recommend fixed target doses of 0.25 or 0.5 mg/kg/day of primaquine for 14 days in glucose-6-phosphate dehydrogenase normal patients of all ages. However, for many anti-malarial drugs, including primaquine, there is evidence that children have lower exposures than adults for the same weight-adjusted dose. The aim of the study was to develop 14-day weight-based and age-based primaquine regimens against high-frequency relapsing tropical P. vivax.

METHODS

The recommended adult target dose of 0.5 mg/kg/day (30 mg in a 60 kg patient) is highly efficacious against tropical P. vivax and was assumed to produce optimal drug exposure. Primaquine doses were calculated using allometric scaling to derive a weight-based primaquine regimen over a weight range from 5 to 100 kg. Growth curves were constructed from an anthropometric database of 53,467 individuals from the Greater Mekong Subregion (GMS) to define weight-for-age relationships. The median age associated with each weight was used to derive an age-based dosing regimen from the weight-based regimen.

RESULTS

The proposed weight-based regimen has 5 dosing bands: (i) 5-7 kg, 5 mg, resulting in 0.71-1.0 mg/kg/day; (ii) 8-16 kg, 7.5 mg, 0.47-0.94 mg/kg/day; (iii) 17-40 kg, 15 mg, 0.38-0.88 mg/kg/day; (iv) 41-80 kg, 30 mg, 0.37-0.73 mg/kg/day; and (v) 81-100 kg, 45 mg, 0.45-0.56 mg/kg/day. The corresponding age-based regimen had 4 dosing bands: 6-11 months, 5 mg, 0.43-1.0 mg/kg/day; (ii) 1-5 years, 7.5 mg, 0.35-1.25 mg/kg/day; (iii) 6-14 years, 15 mg, 0.30-1.36 mg/kg/day; and (iv) ≥ 15 years, 30 mg, 0.35-1.07 mg/kg/day.

CONCLUSION

The proposed weight-based regimen showed less variability around the primaquine dose within each dosing band compared to the age-based regimen and is preferred. Increased dose accuracy could be achieved by additional dosing bands for both regimens. The age-based regimen might not be applicable to regions outside the GMS, which must be based on local anthropometric data. Pharmacokinetic data in small children are needed urgently to inform the proposed regimens.

A Population Pharmacokinetic Modeling and Simulation Study of Posaconazole Oral Suspension in Immunocompromised Pediatric Patients: A Short Communication

BACKGROUND

Posaconazole oral suspension emerged as a promising candidate for prophylaxis of invasive fungal infections in immunocompromised children. Its pharmacodynamic advantages include a broad-spectrum activity and a favorable safety profile; however, they are overshadowed by its large pharmacokinetic (PK) variability, which might cause subtherapeutic exposure. The aim of this study was to develop a population (pop) PK model based on rich sampling data to better understand the PK of posaconazole oral suspension in pediatric patients.

METHODS

Data were obtained from a prospective interventional study involving hospitalized pediatric patients with a hematologic malignancy and prophylactically treated with posaconazole oral suspension. After constructing the popPK model, the probability of target attainment (PTA; 100% T ≥ 0.7 mg/L) for prophylaxis under fixed, body weight-based, and body surface area-based dosing was evaluated using Monte Carlo simulation.

RESULTS

Fourteen patients contributed 112 posaconazole plasma concentrations. The PK of posaconazole was adequately described by a 1-compartment model with lag time 2.71 hours [13%]; nonlinear bioavailability ED50 99.1 mg/m2 (fixed); first-order absorption rate constant 0.325 hour-1 [27%]; apparent volume of distribution 1150 L [34%]; and apparent clearance 15.4 L/h [24%] (∼70-kg individual). The bioavailability decreased in the presence of diarrhea and co-treatment with a proton pump inhibitor (PPI). The unexplained interindividual variability in posaconazole PK remained large. The PTA was <85%, irrespective of the simulated dosing strategy. Patients without diarrhea and not administered a PPI had the highest PTA (85% under the fixed 300-mg dosing 4 times per day).

CONCLUSIONS

Therapeutic drug monitoring is recommended during prophylactic posaconazole therapy in immunocompromised pediatric patients. Large-scale comparative studies are needed to characterize the PK variability between different posaconazole formulations in this cohort.

Pharmacometrics meets statistics-A synergy for modern drug development

Modern drug development problems are very complex and require integration of various scientific fields. Traditionally, statistical methods have been the primary tool for design and analysis of clinical trials. Increasingly, pharmacometric approaches using physiology-based drug and disease models are applied in this context. In this paper, we show that statistics and pharmacometrics have more in common than what keeps them apart, and collectively, the synergy from these two quantitative disciplines can provide greater advances in clinical research and development, resulting in novel and more effective medicines to patients with medical need.

Abacavir pharmacokinetics in African children living with HIV: A pooled analysis describing the effects of age, malnutrition and common concomitant medications

Aims

Abacavir is part of WHO‐recommended regimens to treat HIV in children under 15 years of age. In a pooled analysis across four studies, we describe abacavir population pharmacokinetics to investigate the influence of age, concomitant medications, malnutrition and formulation.

Methods

A total of 230 HIV‐infected African children were included, with median (range) age of 2.1 (0.1–12.8) years and weight of 9.8 (2.5–30.0) kg. The population pharmacokinetics of abacavir was described using nonlinear mixed‐effects modelling.

Results

Abacavir pharmacokinetics was best described by a two‐compartment model with first‐order elimination, and absorption described by transit compartments. Clearance was predicted around 54% of its mature value at birth and 90% at 10 months. The estimated typical clearance at steady state was 10.7 L/h in a child weighing 9.8 kg co‐treated with lopinavir/ritonavir, and was 12% higher in children receiving efavirenz. During coadministration of rifampicin‐based antituberculosis treatment and super‐boosted lopinavir in a 1:1 ratio with ritonavir, abacavir exposure decreased by 29.4%. Malnourished children living with HIV had higher abacavir exposure initially, but this effect waned with nutritional rehabilitation. An additional 18.4% reduction in clearance after the first abacavir dose was described, suggesting induction of clearance with time on lopinavir/ritonavir‐based therapy. Finally, absorption of the fixed dose combination tablet was 24% slower than the abacavir liquid formulation.

Conclusion

In this pooled analysis we found that children on lopinavir/ritonavir or efavirenz had similar abacavir exposures, while concomitant TB treatment and super‐boosted lopinavir gave significantly reduced abacavir concentrations.

Nonparametric goodness-of-fit testing for parametric covariate models in pharmacometric analyses

The characterization of covariate effects on model parameters is a crucial step during pharmacokinetic/pharmacodynamic analyses. Although covariate selection criteria have been studied extensively, the choice of the functional relationship between covariates and parameters, however, has received much less attention. Often, a simple particular class of covariate‐to‐parameter relationships (linear, exponential, etc.) is chosen ad hoc or based on domain knowledge, and a statistical evaluation is limited to the comparison of a small number of such classes. Goodness‐of‐fit testing against a nonparametric alternative provides a more rigorous approach to covariate model evaluation, but no such test has been proposed so far. In this manuscript, we derive and evaluate nonparametric goodness‐of‐fit tests for parametric covariate models, the null hypothesis, against a kernelized Tikhonov regularized alternative, transferring concepts from statistical learning to the pharmacological setting. The approach is evaluated in a simulation study on the estimation of the age‐dependent maturation effect on the clearance of a monoclonal antibody. Scenarios of varying data sparsity and residual error are considered. The goodness‐of‐fit test correctly identified misspecified parametric models with high power for relevant scenarios. The case study provides proof‐of‐concept of the feasibility of the proposed approach, which is envisioned to be beneficial for applications that lack well‐founded covariate models.

Optimizing Meropenem Therapy for Severe Nosocomial Infections in Neonates

Meropenem pharmacokinetics in neonates exhibits large interindividual variability due to developmental changes occurring during the first month of life. The objective was to characterize meropenem pharmacokinetics through a population approach to determine effective dosing recommendations in neonates with severe nosocomial infections. Three blood samples from forty neonates were obtained once steady-state blood levels were achieved and plasma concentrations were determined with a validated chromatographic method. Data were used to develop and validate the one-compartment with first-order elimination population pharmacokinetic model obtained by non-linear mixed effect modeling. The final model was Clearance (L/h) = 2.23 × Creatinine Clearance (L/h) and Volume of distribution(L) = 6.06 × Body Surface Area(m²) × (1 + 0.60 if Fluticasone comedication). Doses should be adjusted based on said covariates to increase the likelihood of achieving therapeutic targets. This model explains 12.9% of the interindividual variability for meropenem clearance and 19.1% for volume of distribution. Stochastic simulations to establish initial dosing regimens to maximize the time above the MIC showed that the mean probabilities to achieve the PK/PD target (PTA) for microorganisms with a MIC of 2 and 8 µg/mL were 0.8 and 0.7 following i.v. bolus of 250 and 500 mg/m²/dose q8h, respectively. Meropenem extended 4h infusion would improve PTA in neonates with augmented creatinine clearance.

A Semi-Mechanistic Population Pharmacokinetic/Pharmacodynamic Model of Bortezomib in Pediatric Patients with Relapsed/Refractory Acute Lymphoblastic Leukemia

INTRODUCTION

The pharmacokinetics (PK) of the 20S proteasome inhibitor bortezomib are characterized by a large volume of distribution and a rapid decline in plasma concentrations within the first hour after administration. An increase in exposure was observed in the second week of treatment, which has previously been explained by extensive binding of bortezomib to proteasome in erythrocytes and peripheral tissues. We characterized the nonlinear population PK and pharmacodynamics (PD) of bortezomib in children with acute lymphoblastic leukemia.

METHODS

Overall, 323 samples from 28 patients were available from a pediatric clinical study investigating bortezomib at an intravenous dose of 1.3 mg/m² twice weekly (Dutch Trial Registry number 1881/ITCC021). A semi-physiological PK model for bortezomib was first developed; the PK were linked to the decrease in 20S proteasome activity in the final PK/PD model.

RESULTS

The plasma PK data were adequately described using a two-compartment model with linear elimination. Increased concentrations were observed in week 2 compared with week 1, which was described using a Langmuir binding model. The decrease in 20S proteasome activity was best described by a direct effect model with a sigmoidal maximal inhibitory effect, representing the relationship between plasma concentrations and effect. The maximal inhibitory effect was 0.696 pmol AMC/s/mg protein (95% confidence interval 0.664-0.728) after administration.

CONCLUSION

The semi-physiological model adequately described the nonlinear PK and PD of bortezomib in plasma. This model can be used to further optimize dosing of bortezomib.

Ampicillin Pharmacokinetics During First Week of Life in Preterm and Term Neonates

BACKGROUND AND AIMS

Ampicillin is 1 of the most commonly used antibiotics for treatment of early onset sepsis, but its pharmacokinetics (PK) is poorly characterized. We aimed to define the dose of ampicillin for late preterm and term neonates by evaluating its PK in serum, cerebrospinal (CSF), and epithelial lining fluid.

METHODS

A prospective study included neonates receiving ampicillin for suspected or proven early onset sepsis and pneumonia. PK samples were collected at steady state, at predose and 5 minutes, 1 hour, 3 hours, 8 hours, and 12 hours after ampicillin 3-minute infusion. Ampicillin concentrations were measured by ultra-high-performance liquid chromatography. Noncompartmental anaysis (NCA) and population pharmacokinetic (pop-PK) modeling were performed and probability of therapeutic target attainment was simulated.

RESULTS

In 14 neonates (GA of 32-42 wks; mean BW 2873 g), PK parameters (mean ± SD) in NCA were the following: half-life 7.21 ± 7.97 hours; volume of distribution (Vd) 1.07 ± 0.51 L; clearance (CL) 0.20 ± 0.13 L/h; 24-hour area under the concentration-time curve 348.92 ± 114.86 mg*h/L. In pop-PK analysis, a 2-compartmental model described the data most adequately with the final parameter estimates of CL 15.15 (CV 40.47%) L/h/70kg; central Vd 24.87 (CV 37.91%) L/70kg; intercompartmental CL 0.39 (CV 868.56) L/h and peripheral Vd 1.039 (CV 69.32%) L. Peutic target attainment simulations demonstrated that a dosage of 50 mg/kg q 12 hours attained 100% fT > MIC 0.25 mg/L, group B streptococcal breakpoint.

CONCLUSIONS

We recommend ampicillin dosage 50 mg/kg q 12 hours for neonates with gestational age ≥32 weeks during the first week of life.

Variation in Target Attainment of Beta-Lactam Antibiotic Dosing Between International Pediatric Formularies

As antimicrobial susceptibility of common bacterial pathogens decreases, ensuring optimal dosing may preserve the use of older antibiotics in order to limit the spread of resistance to newer agents. Beta-lactams represent the most widely prescribed antibiotic class, yet most were licensed prior to legislation changes mandating their study in children. As a result, significant heterogeneity persists in the pediatric doses used globally, along with quality of evidence used to inform dosing. This review summarizes dosing recommendations from the major pediatric reference sources and tries to answer the questions: Does beta-lactam dose heterogeneity matter? Does it impact pharmacodynamic target attainment? For three important severe clinical infections-pneumonia, sepsis, and meningitis-pharmacokinetic models were identified for common for beta-lactam antibiotics. Real-world demographics were derived from three multicenter point prevalence surveys. Simulation results were compared with minimum inhibitory concentration distributions to inform appropriateness of recommended doses in targeted and empiric treatment. While cephalosporin dosing regimens are largely adequate for target attainment, they also pose the most risk of neurotoxicity. Our review highlights aminopenicillin, piperacillin, and meropenem doses as potentially requiring review/optimization in order to preserve the use of these agents in future.

Population Pharmacokinetics of Intraventricular Vancomycin in Neonatal Ventriculitis, A Preterm Pilot Study

•

Pharmacokinetics modelling of intraventricular vancomycin in a preterm pilot study.

•

Intraventricular vancomycin follows a one compartment model in neonatal ventriculitis treatment.

•

Ventricular Index, a dosing parameter, does not influence cerebrospinal fluid vancomycin levels.

Aim

Intraventricular vancomycin is an effective treatment for neonatal ventriculitis, as the cerebrospinal fluid (CSF) vancomycin levels reach adequate concentrations to achieve microbiological cure. There is no robust data on intraventricular vancomycin pharmacokinetics in the preterm population. This pilot population pharmacokinetic modelling study examines the pharmacokinetic behaviour of intraventricular vancomycin in the preterm population of < 28 weeks gestation, to inform the feasibility of future prospective studies.

Methods

The study comprised 8 preterm infants with neonatal ventriculitis (median gestation age 25.3 weeks; range 23.9 - 27.7). Population pharmacokinetics (non-linear mixed effects modelling) were described with one- and two-compartment models to fit plasma concentrations of vancomycin. A CSF compartment was added to the plasma modelling and mass transfer examined. Three covariates (serum creatinine, ventricular index (VI) and CSF protein) were tested on the final model. Area under the curve (AUC) and average CSF concentration (C average) predictions were generated from the final model and compared with time to microbiological cure.

Results

A one-compartment model provided the best fit to the data. There was no appreciable transfer between plasma and CSF. None of the covariates provided a significant reduction in the objective function value (OFV). Generally, time to sterilisation with higher CSF AUC (0-24) and C average tends to be shorter, however this should be interpreted with caution as data is erratic.

Conclusion

This pilot population pharmacokinetic analysis provides important information to warrant changes in the management of intraventricular vancomycin treatment in the preterm population, such as the current use of VI as a dosing parameter. Further study with a larger data pool is necessary to investigate the influence of VI on CSF vancomycin and ascertain dosing strategies.

Early prognostic performance of miR155-5p monitoring for the risk of rejection: Logistic regression with a population pharmacokinetic approach in adult kidney transplant patients

Previous results from our group and others have shown that urinary pellet expression of miR155-5p and urinary CXCL-10 production could play a key role in the prognosis and diagnosis of acute rejection (AR) in kidney transplantation patients. Here, a logistic regression model was developed using NONMEM to quantify the relationships of miR155-5p urinary expression, CXCL-10 urinary concentration and tacrolimus and mycophenolic acid (MPA) exposure with the probability of AR in adult kidney transplant patients during the early post-transplant period. Owing to the contribution of therapeutic drug monitoring to achieving target exposure, neither tacrolimus nor MPA cumulative exposure was identified as a predictor of AR in the studied population. Even though CXCL-10 urinary concentration showed a trend, its effect on AR was not significant. In contrast, urinary miR155-5p expression was prognostic of clinical outcome. Monitoring miR155-5p urinary pellet expression together with immunosuppressive drug exposure could be very useful during routine clinical practice to identify patients with a potential high risk of rejection at the early stages of the post-transplant period. This early risk assessment would allow for the optimization of treatment and improved prevention of AR.

Allometric scaling of therapeutic monoclonal antibodies in preclinical and clinical settings

Constant technological advancement enabled the production of therapeutic monoclonal antibodies (mAbs) and will continue to contribute to their rapid expansion. Compared to small-molecule drugs, mAbs have favorable characteristics, but also more complex pharmacokinetics (PK), e.g., target-mediated nonlinear elimination and recycling by neonatal Fc-receptor. This review briefly discusses mAb biology, similarities and differences in PK processes across species and within human, and provides a detailed overview of allometric scaling approaches for translating mAb PK from preclinical species to human and extrapolating from adults to children. The approaches described here will remain vital in mAb drug development, although more data are needed, for example, from very young patients and mAbs with nonlinear PK, to allow for more confident conclusions and contribute to further growth of this field. Improving mAb PK predictions will facilitate better planning of (pediatric) clinical studies and enable progression toward the ultimate goal of expediting drug development.

Concept and utility of population pharmacokinetic and pharmacokinetic/pharmacodynamic models in drug development and clinical practice

Due to frequent clinical trial failures and consequently fewer new drug approvals, the need for improvement in drug development has, to a certain extent, been met using model-based drug development. Pharmacometrics is a part of pharmacology that quantifies drug behaviour, treatment response and disease progression based on different models (pharmacokinetic - PK, pharmacodynamic - PD, PK/PD models, etc.) and simulations. Regulatory bodies (European Medicines Agency, Food and Drug Administration) encourage the use of modelling and simulations to facilitate decision-making throughout all drug development phases. Moreover, the identification of factors that contribute to variability provides a basis for dose individualisation in routine clinical practice. This review summarises current knowledge regarding the application of pharmacometrics in drug development and clinical practice with emphasis on the population modelling approach.

The Neonatal and Juvenile Pig in Pediatric Drug Discovery and Development

Pharmacotherapy in pediatric patients is challenging in view of the maturation of organ systems and processes that affect pharmacokinetics and pharmacodynamics. Especially for the youngest age groups and for pediatric-only indications, neonatal and juvenile animal models can be useful to assess drug safety and to better understand the mechanisms of diseases or conditions. In this respect, the use of neonatal and juvenile pigs in the field of pediatric drug discovery and development is promising, although still limited at this point. This review summarizes the comparative postnatal development of pigs and humans and discusses the advantages of the juvenile pig in view of developmental pharmacology, pediatric diseases, drug discovery and drug safety testing. Furthermore, limitations and unexplored aspects of this large animal model are covered. At this point in time, the potential of the neonatal and juvenile pig as nonclinical safety models for pediatric drug development is underexplored.

Population Pharmacokinetics of Tapentadol in Children from Birth to <18 Years Old

OBJECTIVE

The main aim of this analysis was to characterize the pharmacokinetics (PK) of tapentadol in pediatric patients from birth to <18 years old who experience acute pain, requiring treatment with an opioid analgesic.

PATIENTS AND METHODS

Data from four clinical trials and 148 pediatric patients who received a single dose of tapentadol oral or intravenous solution were included. Population PK analysis was performed to determine the contribution of size-related (bodyweight) and function-related (maturation) factors to the changes in oral bioavailability (F), volume of distribution (V), and clearance (CL) with age. Simulations were carried out to compare pediatric exposures to reference adult values.

RESULTS

A one-compartment model with allometric scaling on disposition parameters (using theoretical or estimated exponents) and maturation functions on CL and F best described tapentadol PK. The estimated allometric exponents for CL (0.603) and V (0.820) differed slightly from the theoretical values of 0.75 for CL and 1 for V. A maximum in CL/F was observed at about 2-3 years when expressed on a bodyweight basis. Results for younger children as well as the F estimate were sensitive to the scaling approach, but CL/F and V/F as a function of age for the two scaling approaches led to similar curves within the bioequivalence range except below 5 weeks of age. Model-based simulations indicated that the doses used in the included clinical trials lead to exposures within the lower half of the targeted adult exposure.

CONCLUSION

The development of tapentadol is one of the first examples following a systematic approach for analgesic drug development for children. Our analysis enabled a full characterization and robust understanding of tapentadol PK in children from birth to <18 years, including preterm infants, and showed the importance of evaluating the sensitivity of the inferences of the PK parameters to the selected scaling approach.

Population pharmacokinetics of propofol in neonates and infants: Gestational and postnatal age to determine clearance maturation

AIMS

Develop a population pharmacokinetic model describing propofol pharmacokinetics in (pre)term neonates and infants, that can be used for precision dosing (e.g. during target-controlled infusion) of propofol in this population.

METHODS

A nonlinear mixed effects pharmacokinetic analysis (Monolix 2018R2) was performed, based on a pooled study population in 107 (pre)term neonates and infants.

RESULTS

In total, 836 blood samples were collected from 66 (pre)term neonates and 41 infants originating from 3 studies. Body weight (BW) of the pooled study population was 3.050 (0.580-11.440) kg, postmenstrual age (PMA) was 36.56 (27.00-43.00) weeks and postnatal age (PNA) was 1.14 (0-104.00) weeks (median and min-max range). A 3-compartment structural model was identified and the effect of BW was modelled using fixed allometric exponents. Elimination clearance maturation was modelled accounting for the maturational effect on elimination clearance until birth (by gestational age [GA]) and postpartum (by PNA and GA). The extrapolated adult (70 kg) population propofol elimination clearance (1.64 L min^-1 , estimated relative standard error = 6.02%) is in line with estimates from previous population pharmacokinetic studies. Empirical scaling of BW on the central distribution volume in function of PNA improved the model fit.

CONCLUSIONS

It is recommended to describe elimination clearance maturation by GA and PNA instead of PMA on top of size effects when analyzing propofol pharmacokinetics in populations including preterm neonates. Changes in body composition in addition to weight changes or other physio-anatomical changes may explain the changes in central distribution volume. The developed model may serve as a prior for propofol dose finding and target-controlled infusion in (preterm) neonates.

Perioperative pharmacokinetics and pharmacodynamics of meloxicam in emus (Dromaius novaehollandiae) of different age groups using nonlinear mixed effect modelling

Meloxicam is a widely used nonsteroidal anti-inflammatory drug in avian species. However, variability in pharmacokinetic (PK) and pharmacodynamic (PD) parameters in birds warrants species-specific studies for dose and dosing interval optimization. We performed a perioperative PK study of meloxicam (0.5 mg/kg, intravenously) on emus of three different age groups: 3 chicks (5 weeks old, 3.5 kg), 4 juveniles (26 weeks old, 18.8 kg) and 6 adults (66 weeks old, 38.8 kg). A two-compartment population PK model including weight as a significant covariate on clearance and central volume of distribution (V1) best fitted the data. The typical values (20 kg bird) for clearance and V1 were 0.54 L/kg/h and 0.095 L/kg. Both parameters significantly decreased with increasing weight/age. Meloxicam potency and selectivity for COX-1 and COX-2 were measured in whole blood assays (TxB₂ production endpoint). Meloxicam was partially selective in emus (IC₅₀ COX-1:COX-2 = 9.1:1). At the current empirical dose (0.5 mg/kg/24 hr), plasma meloxicam concentration is above IC₅₀ of COX-2 for only 2 hr. PK/PD predicted dose required for 80% COX-2 inhibition over 24 hr were 3.4, 1.4 and 0.95 L/kg/day in chicks, juveniles and adult emus, respectively. The safety, therapeutic efficacy and practicality of modifying the daily dose or dose interval should be considered for dose recommendations in emus.

Population Pharmacokinetics of Unbound and Total Teicoplanin in Critically Ill Pediatric Patients

BACKGROUND AND OBJECTIVES

Teicoplanin is a highly protein-bound antibiotic, increasingly used to treat serious Gram-positive infections in critically ill children. Maturational and pathophysiological intensive care unit-related changes often lead to altered pharmacokinetics. In this study, the objectives were to develop a pediatric population-pharmacokinetic model of unbound and total teicoplanin concentrations, to investigate the impact of plasma albumin levels and renal function on teicoplanin pharmacokinetics, and to evaluate the efficacy of the current weight-based dosing regimen.

METHODS

An observational pharmacokinetic study was performed and blood samples were collected for quantification of unbound and total concentrations of teicoplanin after the first dose and in assumed steady-state conditions. A population-pharmacokinetic analysis was conducted using a standard sequential approach and Monte Carlo simulations were performed for a probability of target attainment analysis using previously published pharmacokinetic-pharmacodynamic targets.

RESULTS

A two-compartment model with allometric scaling of pharmacokinetic parameters and non-linear plasma protein binding best described the data. Neither the inclusion of albumin nor the renal function significantly improved the model and no other covariates were supported for inclusion in the final model. The probability of target attainment analysis showed that the standard dosing regimen does not satisfactory attain the majority of the proposed targets.

CONCLUSIONS

We successfully characterized the pharmacokinetics of unbound and total teicoplanin in critically ill pediatric patients. The highly variable unbound fraction of teicoplanin could not be predicted using albumin levels, which may support the use of therapeutic drug monitoring of unbound concentrations. Poor target attainment was shown for the most commonly used dosing regimen, regardless of the pharmacokinetic-pharmacodynamic target evaluated.

Effects of intraoperative dexmedetomidine on the incidence of acute kidney injury in pediatric cardiac surgery patients: A randomized controlled trial

BACKGROUND

Perioperative dexmedetomidine use has been reported to reduce the incidence of postoperative acute kidney injury after adult cardiac surgery. However, large-scale randomized controlled trials evaluating the effect of dexmedetomidine use on acute kidney injury in pediatric patients are lacking.

AIMS

We investigated whether intraoperative dexmedetomidine could reduce the incidence of acute kidney injury in pediatric cardiac surgery patients.

METHODS

In total, 141 pediatric patients were randomly assigned to dexmedetomidine or control groups. After anesthetic induction, patients in the dexmedetomidine group were administered 1 µg/kg of dexmedetomidine over 10 minutes and an additional 0.5 µg/kg/h of dexmedetomidine during surgery. Additionally, 1 µg/kg of dexmedetomidine was infused immediately after cardiopulmonary bypass was initiated. The incidence of acute kidney injury was defined following Kidney Disease Improving Global Outcomes guidelines.

RESULTS

The final analysis included 139 patients. The incidence of acute kidney injury did not differ between dexmedetomidine and control groups (16.9% vs 23.5%; odds ratio 0.661; 95% CI 0.285 to 1.525; P = .33). Similarly, neither the incidence of abnormal postoperative estimated glomerular filtration rate values (P = .96) nor the incidence of arrhythmia, mechanical ventilation duration, length of stay in the intensive care unit, and hospitalization differed between the two groups.

CONCLUSIONS

Intraoperative dexmedetomidine did not reduce acute kidney injury incidence in pediatric cardiac surgery patients.

β-Lactam antimicrobial pharmacokinetics and target attainment in critically ill patients aged 1 day to 90 years: the ABDose study

Background

The pharmacokinetics of β-lactam antibiotics in critical illness remain poorly characterized, particularly in neonates, children and the elderly. We undertook a pharmacokinetic study of commonly used β-lactam antibiotics in critically ill patients of all ages. The aims were to produce a whole-life β-lactam pharmacokinetic model and describe the extent to which standard doses achieve pharmacokinetic/pharmacodynamic targets associated with clinical cure.

Patients and methods

A total of 212 critically ill participants with an age range from 1 day (gestational age 24 weeks) to 90 years were recruited from a UK hospital, providing 1339 pharmacokinetic samples. Population pharmacokinetic analysis was undertaken using non-linear mixed-effects modelling (NONMEM) for each drug. Pooled data were used to estimate maturation and decline of β-lactam pharmacokinetics throughout life.

Results

Pharmacokinetic models for eight drugs were described, including what is thought to be the first benzylpenicillin model in critically ill adults. We estimate that 50% of adult β-lactam clearance is achieved by 43 weeks post-menstrual age (chronological plus gestational age). Fifty percent of decline from peak adult clearance occurs by 71 years. Paediatric participants were significantly less likely than adults to achieve pharmacokinetic/pharmacodynamic targets with standard antibiotic doses (P &lt; 0.01).

Conclusions

We believe this to be the first prospective whole-life antibiotic pharmacokinetic study in the critically ill. The study provides further evidence that standard antibiotic doses fail to achieve pharmacokinetic/pharmacodynamic targets associated with clinical success in adults, children and neonates. Maturation and decline parameters estimated from this study could be adopted as a standard for future prospective studies.