Pharmacometric in silico studies used to facilitate a national dose standardisation process in neonatology - application to amikacin.

BACKGROUND AND AIMS: Pharmacometric in silico approaches are frequently applied to guide decisions concerning dosage regimes during the development of new medicines. We aimed to demonstrate how such pharmacometric modelling and simulation can provide a scientific rationale for optimising drug doses in the context of the Swiss national dose standardisation project in paediatrics using amikacin as a case study. METHODS: Amikacin neonatal dosage is stratified by post-menstrual age (PMA) and post-natal age (PNA) in Switzerland and many other countries. Clinical concerns have been raised for the subpopulation of neonates with a post-menstrual age of 30-35 weeks and a post-natal age of 0-14 days ("subpopulation of clinical concern"), as potentially oto-/nephrotoxic trough concentrations (Ctrough >5 mg/l) were observed with a once-daily dose of 15 mg/kg. We applied a two-compartmental population pharmacokinetic model (amikacin clearance depending on birth weight and post-natal age) to real-world demographic data from 1563 neonates receiving anti-infectives (median birth weight 2.3 kg, median post-natal age six days) and performed pharmacometric dose-exposure simulations to identify extended dosing intervals that would ensure non-toxic Ctrough (Ctrough <5 mg/l) dosages in most neonates. RESULTS: In the subpopulation of clinical concern, Ctrough <5 mg/l was predicted in 59% versus 79-99% of cases in all other subpopulations following the current recommendations. Elevated Ctrough values were associated with a post-natal age of less than seven days. Simulations showed that extending the dosing interval to ≥36 h in the subpopulation of clinical concern increased the frequency of a desirable Ctrough below 5 mg/l to >80%. CONCLUSION: Pharmacometric in silico studies using high-quality real-world demographic data can provide a scientific rationale for national paediatric dose optimisation. This may increase clinical acceptance of fine-tuned standardised dosing recommendations and support their implementation, including in vulnerable subpopulations.

Pharmacokinetics-Based Pediatric Dose Evaluation and Optimization Using Saliva - A Case Study.

Understanding pharmacokinetics (PK) in children is a prerequisite to determine optimal pediatric dosing. As plasma sampling in children is challenging, alternative PK sampling strategies are needed. In this case study we evaluated the suitability of saliva as alternative PK matrix to simplify studies in infants, investigating metamizole, an analgesic used off-label in infants. Six plasma and 6 saliva PK sample collections were scheduled after a single intravenous dose of 10 mg/kg metamizole. Plasma/saliva pharmacometric (PMX) modeling of the active metabolites 4-methylaminoantipyrine (4-MAA) and 4-aminoantipyrine (4-AA) was performed. Various reduced plasma sampling scenarios were evaluated by PMX simulations. Saliva and plasma samples from 25 children were included (age range, 5-70 months; weight range, 8.7-24.8 kg). Distribution of metamizole metabolites between plasma and saliva was without delay. Estimated mean (individual range) saliva/plasma fractions of 4-MAA and 4-AA were 0.32 (0.05-0.57) and 0.57 (0.25-0.70), respectively. Residual variability of 4-MAA (4-AA) in saliva was 47% (28%) versus 17% (11%) in plasma. A simplified sampling scenario with up to 6 saliva samples combined with 1 plasma sample was associated with similar PK parameter estimates as the full plasma sampling scenario. This case study with metamizole shows increased PK variability in saliva compared to plasma, compromising its suitability as single matrix for PK studies in infants. Nonetheless, rich saliva sampling can reduce the number of plasma samples required for PK characterization, thereby facilitating the conduct of PK studies to optimize dosing in pediatric patients.

Prospective study to characterize adalimumab exposure in pediatric patients with rheumatic diseases.

BACKGROUND: In pediatric rheumatic diseases (PRD), adalimumab is dosed using fixed weight-based bands irrespective of methotrexate co-treatment, disease activity (DA) or other factors that might influence adalimumab pharmacokinetics (PK). In rheumatoid arthritis (RA) adalimumab exposure between 2-8 mg/L is associated with clinical response. PRD data on adalimumab is scarce. Therefore, this study aimed to analyze adalimumab PK and its variability in PRD treated with/without methotrexate. METHODS: A two-center prospective study in PRD patients aged 2-18 years treated with adalimumab and methotrexate (GA-M) or adalimumab alone (GA) for ≥ 12 weeks was performed. Adalimumab concentrations were collected 1-9 (maximum concentration; Cmax), and 10-14 days (minimum concentration; Cmin) during ≥ 12 weeks following adalimumab start. Concentrations were analyzed with enzyme-linked immunosorbent assay (lower limit of quantification: 0.5 mg/L). Log-normalized Cmin were compared between GA-M and GA using a standard t-test. RESULTS: Twenty-eight patients (14 per group), diagnosed with juvenile idiopathic arthritis (71.4%), non-infectious uveitis (25%) or chronic recurrent multifocal osteomyelitis (3.6%) completed the study. GA-M included more females (71.4%; GA 35.7%, p = 0.13). At first study visit, children in GA-M had a slightly longer exposure to adalimumab (17.8 months [IQR 9.6, 21.6]) compared to GA (15.8 months [IQR 8.5, 30.8], p = 0.8). Adalimumab dosing was similar between both groups (median dose 40 mg every 14 days) and observed DA was low. Children in GA-M had a 27% higher median overall exposure compared to GA, although median Cmin adalimumab values were statistically not different (p = 0.3). Cmin values ≥ 8 mg/L (upper limit RA) were more frequently observed in GA-M versus GA (79% versus 64%). Overall, a wide range of Cmin values was observed in PRD (0.5 to 26 mg/L). CONCLUSION: This study revealed a high heterogeneity in adalimumab exposure in PRD. Adalimumab exposure tended to be higher with methotrexate co-treatment compared to adalimumab monotherapy although differences were not statistically significant. Most children showed adalimumab exposure exceeding those reported for RA with clinical response, particularly with methotrexate co-treatment. This highlights the need of further investigations to establish model-based personalized treatment strategies in PRD to avoid under- and overexposure. TRIAL REGISTRATION: NCT04042792 , registered 02.08.2019.

Characterizing Associations of QTc Interval with Nocturnal Glycemic Control in Children with Type 1 Diabetes.

An association between QT prolongation (Bazett's corrected QT interval, QTcB) of 7 milliseconds and nocturnal hypoglycemia, compared with euglycemia, has been observed in children with type 1 diabetes (T1D). The objective of this pharmacometric analysis was to understand this association and other sources of QTc variability quantitatively. Data originate from a prospective observational study (25 cardiac healthy children with T1D, aged 8.1-17.6 years) with continuous subcutaneous glucose and electrocardiogram measurements for 5 consecutive nights. Mixed-effect modeling was used to compare QTcB with individual heart-rate correction (QTcI). Covariate models accounting for circadian variation, age, and sex were evaluated, followed by an investigation of glucose-QTc relationships (with univariable and combined adjusted analysis). Factors potentially modifying sensitivity to QTc lengthening were explored. Random inter-individual variability was reduced in the QTcI versus QTcB model (±12.6 vs 14.1 milliseconds), and was further reduced in the adjusted covariate model (±9.7 milliseconds), accounting for the significantly (P < .01) shortened QTc in adolescent boys (-14.6 milliseconds), circadian variation (amplitude, 19.2 milliseconds; shift, 2.9 hours), and linear glucose-QTc relationship (delay rate, 0.56-h ; slope, 0.76 milliseconds [95%CI 0.67- 0.85 milliseconds] per 1 mmol/L decrease in glucose). Differing sensitivity was suggested to depend upon hemoglobin A1c (HbA1c), T1D duration, and time spent in nocturnal hypoglycemia. In conclusion, a clinically mild association of QTc prolongation with nocturnal hypoglycemia was confirmed and quantified in this pharmacometric analysis, and the longest QTc interval was around 03:00 a.m. The characterized delayed association with glucose highlights the relevance of both the extent and the duration of hypoglycemia. Further clinical studies are warranted to investigate whether these factors contribute to increased risk of hypoglycemia-associated cardiac arrhythmia in children with T1D.

Clinically practical pharmacometrics computer model to evaluate and personalize pharmacotherapy in pediatric rare diseases: application to Graves' disease.

Objectives: Graves' disease (GD) with onset in childhood or adolescence is a rare disease (ORPHA:525731). Current pharmacotherapeutic approaches use antithyroid drugs, such as carbimazole, as monotherapy or in combination with thyroxine hormone substitutes, such as levothyroxine, as block-and-replace therapy to normalize thyroid function and improve patients' quality of life. However, in the context of fluctuating disease activity, especially during puberty, a considerable proportion of pediatric patients with GD is suffering from thyroid hormone concentrations outside the therapeutic reference ranges. Our main goal was to develop a clinically practical pharmacometrics computer model that characterizes and predicts individual disease activity in children with various severity of GD under pharmacotherapy. Methods: Retrospectively collected clinical data from children and adolescents with GD under up to two years of treatment at four different pediatric hospitals in Switzerland were analyzed. Development of the pharmacometrics computer model is based on the non-linear mixed effects approach accounting for inter-individual variability and incorporating individual patient characteristics. Disease severity groups were defined based on free thyroxine (FT4) measurements at diagnosis. Results: Data from 44 children with GD (75% female, median age 11 years, 62% receiving monotherapy) were analyzed. FT4 measurements were collected in 13, 15, and 16 pediatric patients with mild, moderate, or severe GD, with a median FT4 at diagnosis of 59.9 pmol/l (IQR 48.4, 76.8), and a total of 494 FT4 measurements during a median follow-up of 1.89 years (IQR 1.69, 1.97). We observed no notable difference between severity groups in terms of patient characteristics, daily carbimazole starting doses, and patient years. The final pharmacometrics computer model was developed based on FT4 measurements and on carbimazole or on carbimazole and levothyroxine doses involving two clinically relevant covariate effects: age at diagnosis and disease severity. Discussion: We present a tailored pharmacometrics computer model that is able to describe individual FT4 dynamics under both, carbimazole monotherapy and carbimazole/levothyroxine block-and-replace therapy accounting for inter-individual disease progression and treatment response in children and adolescents with GD. Such clinically practical and predictive computer model has the potential to facilitate and enhance personalized pharmacotherapy in pediatric GD, reducing over- and underdosing and avoiding negative short- and long-term consequences. Prospective randomized validation trials are warranted to further validate and fine-tune computer-supported personalized dosing in pediatric GD and other rare pediatric diseases.

From Personalized to Precision Medicine in Oncology: A Model-Based Dosing Approach to Optimize Achievement of Imatinib Target Exposure.

Imatinib is a targeted cancer therapy that has significantly improved the care of patients with chronic myeloid leukemia (CML) and gastrointestinal stromal tumor (GIST). However, it has been shown that the recommended dosages of imatinib are associated with trough plasma concentration (Cmin) lower than the target value in many patients. The aims of this study were to design a novel model-based dosing approach for imatinib and to compare the performance of this method with that of other dosing methods. Three target interval dosing (TID) methods were developed based on a previously published PK model to optimize the achievement of a target Cmin interval or minimize underexposure. We compared the performance of those methods to that of traditional model-based target concentration dosing (TCD) as well as fixed-dose regimen using simulated patients (n = 800) as well as real patients' data (n = 85). Both TID and TCD model-based approaches were effective with about 65% of Cmin achieving the target imatinib Cmin interval of 1000-2000 ng/mL in 800 simulated patients and more than 75% using real data. The TID approach could also minimize underexposure. The standard 400 mg/24 h dosage of imatinib was associated with only 29% and 16.5% of target attainment in simulated and real conditions, respectively. Some other fixed-dose regimens performed better but could not minimize over- or underexposure. Model-based, goal-oriented methods can improve initial dosing of imatinib. Combined with subsequent TDM, these approaches are a rational basis for precision dosing of imatinib and other drugs with exposure-response relationships in oncology.

Semimechanistic modeling of copeptin and aldosterone kinetics and dynamics in response to rehydration treatment for diabetic ketoacidosis in children.

Diabetic ketoacidosis (DKA), a frequent complication of type 1 diabetes (T1D), is characterized by hyperosmolar hypovolemia. The response of water-regulating hormones arginine vasopressin (AVP; antidiuretic hormone) and aldosterone to DKA treatment in children is not well understood, although they may have potential as future diagnostic, prognostic, and/or treatment monitoring markers in diabetic patients. We aimed to characterize the dynamics of the response in copeptin (marker for AVP) and aldosterone secretion to rehydration treatment in pediatric patients with DKA. Data originated from a prospective, observational, multicenter study including 28 pediatric T1D patients treated for DKA (median age, 11.5 years; weight, 35 kg). Serial measurements of hormone levels were obtained during 72 h following rehydration start. Semimechanistic pharmacometric modeling was used to analyze the kinetic/dynamic relationship of copeptin and aldosterone secretion in response to the correction of hyperosmolality and hypovolemia, respectively. Modeling revealed different sensitivities for osmolality-dependent copeptin secretion during the first 72 h of rehydration, possibly explained by an osmotic shift introduced by hypovolemia. Response in aldosterone secretion to the correction of hypovolemia seemed to be delayed, which was well described by an extra upstream turnover compartment, possibly representing chronic upregulation of aldosterone synthase (cytochrome P450 11B2). In conclusion, semimechanistic modeling provided novel physiological insights in hormonal water regulation in pediatric patients during DKA treatment, providing rationale to further evaluate the potential of monitoring copeptin, but not aldosterone due to its delayed response, for future optimization of rehydration treatment to reduce the risk of acute complications such as cerebral edema.

Pharmacokinetics and tolerability of intranasal or intravenous administration of nalbuphine in infants.

OBJECTIVES: Intranasal nalbuphine could be a safe, efficacious and non-invasive alternative to parenteral pain medication in infants. We aimed to assess pharmacokinetics (PK) and tolerability of intranasal and intravenous nalbuphine administration in infants. METHODS: Prospective open-label study including infants 1-3 months of age admitted to the emergency department, receiving nalbuphine for procedural pain management. Patients were alternately allocated to a single nalbuphine dose of 0.05 mg/kg intravenously or 0.1 mg/kg intranasally. Nalbuphine PK samples were collected 15, 30 and 120-180 min after dosing. Area under the concentration time curve (AUC0-Tlast) was calculated by non-compartmental analysis (NCA) and compared by Wilcoxon test. Neonatal Infant Pain Score was assessed during nalbuphine administration and the following interventions: venous access, urinary catheterisation, lumbar puncture. RESULTS: Out of 52 study subjects receiving nalbuphine, 31 were eligible for NCA (11 intravenous, 20 intranasal). Median AUC0-Tlast after 0.05 mg/kg intravenously was 8.7 (IQR: 8.0-18.6) µg×L/hour vs 7.6 (5.4-10.4) µg×L/hour after intranasal administration of 0.1 mg/kg (p=0.091). Maximum serum concentration (Cmax) was observed 30 min after intranasal administration (3.5-5.6 µg/L). During intravenous and intranasal nalbuphine administration, mild to no pain was recorded in 71% and 67% of study subjects, respectively. CONCLUSION: This is the first study investigating intranasal administration of nalbuphine in infants suggesting an intranasal bioavailability close to 50%. Non-invasive intranasal application was well tolerated. Additional studies are warranted to optimise dosing and timing of interventions as Cmax is delayed by half an hour after intranasal administration. TRIAL REGISTRATION NUMBER: NCT03059511.

Therapeutic Drug Monitoring of Orally Administered Letermovir Prophylaxis in Allogeneic Hematopoietic Stem Cell Transplant Recipients.

With balanced safety-efficacy profile, letermovir anti-cytomegalovirus (CMV) prophylaxis is used in hematopoietic stem cell transplant recipients (HSCTR). We assessed feasibility and usefulness of letermovir therapeutic drug monitoring (TDM) in HSCTR. We performed a prospective observational study on letermovir-TDM including 40 consecutive adult CMV-seropositive allogeneic-HSCTR who received orally (PO) administered letermovir. Minimal blood concentrations of letermovir (Ctrough) were measured on days 3 and 7 postletermovir initiation and weekly thereafter. Letermovir-Ctrough remained stable during the first 70 days post-HSCT at a median of 286 µg/L (interquartile range, 131 to 591 µg/L), with large interpatient/intrapatient variability. No associations between breakthrough clinically significant CMV infection or detectable CMV DNAemia and letermovir-Ctrough were observed. Patients with letermovir-associated adverse events had higher letermovir-Ctrough than patients without (400 versus 266 µg/L, P = 0.02). Letermovir-Ctrough was similar in patients with or without gastrointestinal symptoms (280 versus 300 µg/L, P = 0.49). Acute grade ≥2 GvHD was associated with higher letermovir-Ctrough (479 versus 248 µg/L, P = 0.001), including gastrointestinal GvHD (499 versus 263 µg/L, P = 0.004). Concomitantly administered posaconazole and cyclosporine were associated with higher letermovir-Ctrough (707 versus 259 µg/L, P < 0.001 and 437 versus 248 µg/L, P = 0.01, respectively). In multivariable analysis, both posaconazole (odds ratio [OR], 4.9; 95% confidence interval [CI], 2.4 to 9.7; P < 0.0001) and cyclosporine-adjusted letermovir dose at 240 mg daily (OR, 3.5; 95% CI, 1.4 to 9.0; P = 0.01) were independently associated with higher letermovir-Ctrough. In conclusion, administration of PO letermovir led to measurable and relatively stable letermovir-Ctrough, without noticeable associations with clinical efficacy. Letermovir exposure was not affected by gastrointestinal symptoms, but with posaconazole and cyclosporine administration. Associations between letermovir and concomitantly administered agents and adverse events warrant additional clinical studies.

Copeptin Release in Arterial Hypotension and Its Association with Severity of Disease in Critically Ill Children.

Low copeptin levels may indicate inadequate arginine-vasopressin release promoting arterial hypotension, whereas high copeptin concentrations may reflect disease severity. This single-center prospective non-randomized clinical trial analyzed the course of blood copeptin in critically ill normo- and hypotensive children and its association with disease severity. In 164 patients (median age 0.5 years (interquartile range 0.1, 2.9)), the mean copeptin concentration at baseline was 43.5 pmol/L. Though not significantly different after 61 h (primary outcome, mean individual change: −12%, p = 0.36, paired t-test), we detected 1.47-fold higher copeptin concentrations during arterial hypotension when compared to normotension (mixed-effect ANOVA, p = 0.01). In total, 8 out of 34 patients (23.5%) with low copeptin concentrations <10 pmol/L were hypotensive. Copeptin was highest in the adjusted mixed-effect regression analysis within the first day (+20% at 14 h) and decreased significantly at 108 h (−27%) compared to baseline (p = 0.002). Moreover, we found a significant association with vasopressor-inotrope treatment intensity, infancy (1−12 months) and cardiopulmonary bypass (all p ≤ 0.001). In conclusion, high copeptin values were associated with arterial hypotension and severity of disease in critically ill children. This study does not support the hypothesis that low copeptin values might be indicative of arginine-vasopressin deficiency.

Pharmacometric Analysis of Intranasal and Intravenous Nalbuphine to Optimize Pain Management in Infants.

Objectives: The objective of this pharmacometric (PMX) study was to (i) characterize population pharmacokinetics (PPK) and exposure-pain response associations following intranasal (0.1 mg/kg) or intravenous (IV, 0.05 mg/kg) administration of nalbuphine, with the goal to (ii) evaluate strategies for optimized dosing and timing of painful interventions in infants 1-3 months old. Methods: PPK analysis of nalbuphine serum concentrations, prospectively collected 15, 30, and between 120 and 180 min post-dose, utilizing the software package Monolix. The final PPK model was applied to derive individual time-matched concentration predictions for each pain assessment (Neonatal Infant Pain Score, NIPS) after establishment of venous access and urinary catheterization or lumbar puncture. Drug exposure-pain response simulations were performed to evaluate potential benefits of higher doses with respect to a previously proposed target concentration of 12 mcg/L (efficacy threshold). Results: Thirty-eight of 52 study subjects receiving nalbuphine had at least one concentration measurement and were included in the pharmacometric analysis. A two-compartment model with allometric scaling was applied to describe population PK data, with intranasal bioavailability estimated to be 41% (95%CI: 26-56%). Model-based simulations showed that the proposed efficacy threshold (12 mcg/L) is expected to be exceeded with an IV dose of 0.05 mg/kg for 6 min, with 0.1 mg/kg for 30 min and with 0.2 mg/kg for 80 min. This efficacy threshold is not achieved with intranasal doses of 0.1 and 0.2 mg/kg, whereas an intranasal dose of 0.4 mg/kg is expected to exceed such threshold for 30 to 100 min. Conclusion: This PMX study confirmed that bioavailability of intranasal nalbuphine is close to 50%. Exposure-pain response simulations indicated that an intranasal dose of 0.4 mg/kg is required to provide a comparable pain control as achieved with an IV dose of 0.1-0.2 mg/kg. The optimal time window for painful procedures appears to be within the first 30 min after IV administration of 0.1 mg/kg nalbuphine, whereas such procedures should be scheduled 30 min after an intranasal dose of 0.4 mg/kg nalbuphine. Additional clinical studies are warranted to confirm these PMX based recommendations and to further optimize pain management in this vulnerable infant population.

The kinetic profiles of copeptin and mid regional proadrenomedullin (MR-proADM) in pediatric lower respiratory tract infections.

BACKGROUND: Kinetics of copeptin and mid regional proadrenomedullin (MR-proADM) during febrile pediatric lower respiratory tract infections (LRTI) are unknown. We aimed to analyze kinetic profiles of copeptin and MR-proADM and the impact of clinical and laboratory factors on those biomarkers. METHODS: This is a retrospective post-hoc analysis of a randomized controlled trial, evaluating procalcitonin guidance for antibiotic treatment of LRTI (ProPAED-study). In 175 pediatric patients presenting to the emergency department plasma copeptin and MR-proADM concentrations were determined on day 1, 3, and 5. Their association with clinical characteristics and other inflammatory biomarkers were tested by non-linear mixed effect modelling. RESULTS: Median copeptin and MR-proADM values were elevated on day 1 and decreased during on day 3 and 5 (-26%; -34%, respectively). The initial concentrations of MR-proADM at inclusion were higher in patients receiving antibiotics intravenously compared to oral administration (difference 0.62 pmol/L, 95%CI 0.44;1.42, p<0.001). Intensive care unit (ICU) admission was associated with a daily increase of MR-proADM (increase/day 1.03 pmol/L, 95%CI 0.43;1.50, p<0.001). Positive blood culture in patients with antibiotic treatment and negative results on nasopharyngeal aspirates, or negative blood culture were associated with a decreasing MR-proADM (decrease/day -0.85 pmol/L, 95%CI -0.45;-1.44), p<0.001). CONCLUSION: Elevated MR-proADM and increases thereof were associated with ICU admission suggesting the potential as a prognostic factor for severe pediatric LRTI. MR-proADM might only bear limited value for decision making on stopping antibiotics due to its slow decrease. Copeptin had no added value in our setting.

Accuracy of a Novel Transcutaneous PCO2 and PO2 Sensor with Optical PO2 Measurement in Neonatal Intensive Care: A Single-Centre Prospective Clinical Trial.

BACKGROUND AND OBJECTIVES: Transcutaneous PCO2 and PO2 measurement systems offer non-invasive blood gas trend monitoring. The aim of this prospective study was to assess bias and precision of a transcutaneous PCO2 and PO2 measurement system incorporating a novel pO2 sensor (Sentec OxiVenT™) in neonates ≥34 weeks of gestational age (GA) admitted to intensive care. METHODS: Transcutaneous PCO2 and PO2 were compared to arterial and capillary blood gas measurements. Bias and precision were calculated by fitting linear mixed models to account for repeated measurements, and influence of clinical covariates on bias and precision was assessed. RESULTS: We obtained 611 paired transcutaneous and blood gas measurements in 110 patients (median GA 38.3 [interquartile range 36.1-39.7] weeks; age 9 [4-15] days; weight 3,000 [2,500-3,500] g). Transcutaneous PCO2 showed significant bias to arterial PCO2 (+0.61; 95% confidence interval 0.46, 0.76 kPa), but not to capillary PCO2 (-0.23; -0.46, 0.002 kPa). Bias of transcutaneous PO2 was significant to arterial PO2 (-2.50; -2.94, -2.06 kPa), while no significant bias compared to capillary PO2 was observed (+0.17; -0.30, 0.64 kPa). Precision intervals were ±1.8/2.0 kPa for arterial versus capillary PCO2 and ±4.9/3.3 kPa for arterial versus capillary PO2 comparisons, respectively. Further, sensor operating temperature (43°C vs. 42°C), soft tissue oedema, vasoactive drugs, weight, and GA significantly altered bias (p < 0.05). CONCLUSIONS: The tested transcutaneous blood gas measurement system showed no significant bias compared to capillary PCO2 and PO2, acceptable bias to arterial PCO2, and limited agreement with arterial PO2. Precision intervals were wide for all comparisons.

Modeling of levothyroxine in newborns and infants with congenital hypothyroidism: challenges and opportunities of a rare disease multi-center study.

Modeling of retrospectively collected multi-center data of a rare disease in pediatrics is challenging because laboratory data can stem from several decades measured with different assays. Here we present a retrospective pharmacometrics (PMX) based data analysis of the rare disease congenital hypothyroidism (CH) in newborns and infants. Our overall aim is to develop a model that can be applied to optimize dosing in this pediatric patient population since suboptimal treatment of CH during the first 2 years of life is associated with a reduced intelligence quotient between 10 and 14 years. The first goal is to describe a retrospectively collected dataset consisting of 61 newborns and infants with CH up to 2 years of age. Overall, 505 measurements of free thyroxine (FT4) and 510 measurements of thyrotropin or thyroid-stimulating hormone were available from patients receiving substitution treatment with levothyroxine (LT4). The second goal is to introduce a scale/location-scale normalization method to merge available FT4 measurements since 34 different postnatal age- and assay-specific laboratory reference ranges were applied. This method takes into account the change of the distribution of FT4 values over time, i.e. a transformation from right-skewed towards normality during LT4 treatment. The third goal is to develop a practical and useful PMX model for LT4 treatment to characterize FT4 measurements, which is applicable within a clinical setting. In summary, a time-dependent normalization method and a practical PMX model are presented. Since there is no on-going or planned development of new pharmacological approaches for CH, PMX based modeling and simulation can be leveraged to personalize dosing with the goal to enhance longer-term neurological outcome in children with the rare disease CH.

Serum Creatinine and Serum Cystatin C are Both Relevant Renal Markers to Estimate Vancomycin Clearance in Critically Ill Neonates.

Purpose: Serum creatinine (SCr) is used as a marker of kidney function to guide dosing of renally eliminated drugs. Serum Cystatin C (S-CysC) has been suggested as a more reliable kidney marker than SCr in adults and children. Purpose of this study was to investigate S-CysC as alternative renal marker to SCr for estimating vancomycin clearance in neonates undergoing intensive care. Methods: Vancomycin pharmacokinetics (PK), SCr and S-CysC data were collected in patients undergoing vancomycin treatment in the neonatal intensive care unit of Robert Debré Hospital - Paris. A population PK analysis was performed utilizing routine therapeutic drug monitoring samples. S-CysC and SCr were compared as covariates on vancomycin clearance using stepwise covariate modeling (forward inclusion [p < 0.05] and backward elimination [p < 0.01]). Model performance was evaluated by graphical and statistical criteria. Results: A total of 108 vancomycin concentrations from 66 patients (postmenstrual age [PMA] of 26-46 weeks) were modeled with an allometric one-compartment model. The median (range) values for SCr and S-CysC were 41 (12-153) µmol/l and 1.43 (0.95-2.83) mg/l, respectively. Following stepwise covariate model building, SCr was retained as single marker of kidney function (after accounting for weight and PMA) in the final model. Compared to the final model based on SCr, the alternative model based on S-CysC showed very similar performance (e.g. BIC of 578.3 vs. 576.4) but included one additional covariate: impact of mechanical ventilation on vancomycin clearance, in addition to the effects of size and maturation. Conclusion: ill neonates. However, if using S-CysC for this purpose mechanical ventilation needs to be taken into account.

Hemodialysis (HD) dose and ultrafiltration rate are associated with survival in pediatric and adolescent patients on chronic HD-a large observational study with follow-up to young adult age.

BACKGROUND: Hemodialysis (HD) dose targets and ultrafiltration rate (UFR) limits for pediatric patients on chronic HD are not known and are derived from adults (spKt/V>1.4 and <13 ml/kg/h). We aimed to characterize how delivered HD dose and UFR are associated with survival in a large cohort of patients who started HD in childhood. METHODS: Retrospective analysis on a cohort of patients <30 years, on chronic HD since childhood (<19 years), having received thrice-weekly HD 2004-2016 in outpatient DaVita centers. OUTCOME: Survival while remaining on HD. PREDICTORS: (I) primary analysis: mean delivered dialysis dose stratified as spKt/V ≤1.4/1.4-1.6/>1.6 (Kaplan-Meier analysis), (II) secondary analyses: UFR and alternative dialysis adequacy measures [eKt/V, body-surface normalized Kt/BSA] on continuous scale (Weibull regression model). RESULTS: A total of 1780 patients were included (age at the start of HD: 0-12y: n=321, >12-18y: n=1459; median spKt/V=1.55, eKt/V=1.31, Kt/BSA=31.2 L/m2, UFR=10.6 mL/kg/h). (I) spKt/V<1.4 was associated with lower survival compared to spKt/V>1.4-1.6 (P<0.001, log-rank test), and spKt/V>1.6 (P<0.001), with 10-year survival of 69.3% (59.4-80.9%) versus 83.0% (76.8-89.8%) and 84.0% (79.6-88.5%), respectively. (II) Kt/BSA was a better predictor of survival than spKt/V or eKt/V. UFR was additionally associated with survival (P<0.001), with increased mortality <10/>18 mL/kg/h. Associations did not alter significantly following adjustment for demographic characteristics (age, etiology of kidney disease, and ethnicity). CONCLUSIONS: Our results suggest usefulness of targeting Kt/BSA>30 L/m2 for best long-term outcomes, corresponding to spKt/V>1.4 (>12 years) and >1.6 (<12 years). In contrast to adults, higher UFR of 10-18 ml/kg/h was not associated with greater mortality in this population.

Drug-drug and drug-food interactions in an infant with early-onset SCN2A epilepsy treated with carbamazepine, phenytoin and a ketogenic diet.

Sodium channel 2 subunit α (SCN2A) mutations cause difficult-to-treat early-onset epilepsy. Effective treatment includes high-dose phenytoin or carbamazepine ± ketogenic diet (KD). We describe an infant with early-onset SCN2A-epilepsy with subtherapeutic carbamazepine concentration during transition from phenytoin treatment to avoid long-term neurotoxicity. The transition from high-dose phenytoin (20 mg kg-1 d-1 , concentration: ≥20 mg/L) with KD, to carbamazepine (50-75 mg kg-1 d-1 , concentration: 9-12 mg/L) lasted 85 days, which we suspected was due to significant drug-drug and/or drug-food interactions. Model-based analysis of carbamazepine pharmacokinetics quantified significant time- and dose-dependent phenytoin-mediated CYP3A4 induction and carbamazepine concentration-dependent auto-induction (apparent clearance increased up to 2.5/3-fold). Lower carbamazepine concentrations under KD were modelled as decreased relative bioavailability (44%), potentially related to decreased fraction absorbed (unexpected for this lipophilic drug), increased intestinal/hepatic metabolism and/or decreased protein-binding with KD. This suggests importance of carbamazepine-concentration monitoring during KD-introduction/removal and necessity of high carbamazepine doses to achieve therapeutic concentrations, especially in infants treated with high-dose phenytoin.

Identifying key predictors of mortality in young patients on chronic haemodialysis-a machine learning approach.

BACKGROUND: The mortality risk remains significant in paediatric and adult patients on chronic haemodialysis (HD) treatment. We aimed to identify factors associated with mortality in patients who started HD as children and continued HD as adults. METHODS: The data originated from a cohort of patients <30 years of age who started HD in childhood (≤19 years) on thrice-weekly HD in outpatient DaVita dialysis centres between 2004 and 2016. Patients with at least 5 years of follow-up since the initiation of HD or death within 5 years were included; 105 variables relating to demographics, HD treatment and laboratory measurements were evaluated as predictors of 5-year mortality utilizing a machine learning approach (random forest). RESULTS: A total of 363 patients were included in the analysis, with 84 patients having started HD at <12 years of age. Low albumin and elevated lactate dehydrogenase (LDH) were the two most important predictors of 5-year mortality. Other predictors included elevated red blood cell distribution width or blood pressure and decreased red blood cell count, haemoglobin, albumin:globulin ratio, ultrafiltration rate, z-score weight for age or single-pool Kt/V (below target). Mortality was predicted with an accuracy of 81%. CONCLUSIONS: Mortality in paediatric and young adult patients on chronic HD is associated with multifactorial markers of nutrition, inflammation, anaemia and dialysis dose. This highlights the importance of multimodal intervention strategies besides adequate HD treatment as determined by Kt/V alone. The association with elevated LDH was not previously reported and may indicate the relevance of blood-membrane interactions, organ malperfusion or haematologic and metabolic changes during maintenance HD in this population.

Characterizing dynamics of serum creatinine and creatinine clearance in extremely low birth weight neonates during the first 6 weeks of life.

BACKGROUND: Characterizing the dynamics of serum creatinine concentrations (Scr) and associated creatinine clearance (CLcr) as a measure of kidney function in extremely low birth weight (≤ 1000 g; ELBW) neonates remains challenging. METHODS: We performed a retrospective study that included longitudinal Scr (enzymatic assay) data from 148 ELBW neonates up to 6 weeks after birth. Change of Scr and inter-individual variability was characterized with nonlinear mixed-effect modeling. Key covariates such as gestational age (GA), mode of delivery (MOD), and treatment with ibuprofen or inotropic agents were investigated. RESULTS: A total of 2814 Scr concentrations were analyzed. GA was associated with Scr at birth (higher with advancing GA), and GA and MOD showed an association with postnatal maturation of CLcr (faster clearance increase with advancing GA and after C-section). Small CLcr decrease (≤ 5%) was quantified during ibuprofen treatment. For a GA of 27 weeks, mean Scr (estimated CLcr) at birth was 0.61 mg/dl (0.23 ml/min), increasing to 0.87 mg/dl (0.27 ml/min) at day three, and decreasing to 0.36 mg/dl (0.67 ml/min) at day 42 after birth. CONCLUSIONS: We report the first mathematical model able to characterize Scr and CLcr in ELBW neonates during the first 6 weeks of life in a quantitative manner as a function of GA, MOD, and ibuprofen treatment. This model allows the derivation of GA-adjusted reference ranges for ELBW neonates and provides a rationale for normative Scr concentrations, and as such will help clinicians to further optimize monitoring and treatment decisions in this vulnerable patient population.