Establishing best practices and guidance in population modeling: an experience with an internal population pharmacokinetic analysis guidance

Exploring the Influence of Obesity and CYP2 C19 Genotype on Lansoprazole Pharmacokinetics in Pediatric Patients: A Population Modeling Approach

BACKGROUND AND OBJECTIVE

Pediatric pharmacokinetics pose unique challenges, particularly concerning drug dosing in the presence of obesity and genetic variability in drug-metabolizing enzymes such as cytochrome P450 (CYP) 2C19. This study aimed to develop a lansoprazole pediatric population pharmacokinetic (popPK) model considering obesity, obesity-associated changes in inflammatory cytokines and the liver, and CYP2C19 genotype, and to assess implications for dosing strategies.

METHODS

Pediatric subjects with and without obesity (6-21 years, n = 47) received one oral dose of 1.2 mg/kg fat-free mass (FFM), not exceeding 60 mg. Plasma concentrations were measured at multiple time points, and a popPK analysis using NONMEM with stepwise covariate modeling was performed. Simulations compared exposures between children without and with obesity (BMI percentile ≥ 95th) after two dosing strategies: U.S. Food and Drug Administration (FDA)-approved total body weight-tiered and FFM-based dosing.

RESULTS

A total of 537 lansoprazole concentrations were modeled using a two-compartment model with Weibull absorption and linear elimination. Parameters were allometrically scaled to FFM with fixed exponents of 0.75 for clearances and 1 for volumes. CYP2C19 was identified as a significant covariate for CL/F. No significant differences in lansoprazole exposure were observed between children with and without obesity, with both dosing approaches. Higher exposure was noted in poor/intermediate metabolizers of CYP2C19. FFM-based dosing led to similar levels of exposure between children (aged 6-11 years) and adolescents (aged 12-17 years).

CONCLUSIONS

Obesity was not associated with differences in lansoprazole pharmacokinetics in children. A FFM-based dosing approach could result in comparable exposure between children and adolescents. Dose adjustments are supported for poor/intermediate metabolizers of CYP2C19.

Pharmacokinetic Modeling and Model-Based Hypothesis Generation for Dose Optimization of Clonidine in Neonates With Neonatal Opioid Withdrawal Syndrome

The No-POPPY study (NCT03396588), a double-blind, randomized trial compared morphine with clonidine therapy for neonatal opioid withdrawal syndrome (NOWS) and found that the duration of treatment was similar across groups. This is significant because perinatal use of morphine has the potential for neurodevelopmental consequences. Still, the clonidine group reached symptom stabilization (Finnegan score (FS) < 8) later than the morphine group and had a more significant number of patients who required adjunct therapy. However, the mean FS was consistently lower in the clonidine group after day 6. This prompted us to use pharmacokinetic (PK) and parametric time-to-event (TTE) modeling to simulate dosage schedules that may decrease the time to stabilization and reduce the need for adjunct therapy. Population PK (popPK) analysis was conducted, and the final model was a one-compartment model with first-order absorption and elimination, incorporating allometric scaling and age effect on apparent clearance (CL/F) and apparent volume (V/F). The population estimates for CL/F and V/F were 13.6 L/h/70 kg and 416 L/70 kg, respectively, similar to the reported values. A Weibull model described the TTE data best, followed by incorporating predicted average concentrations to yield the final Weibull accelerated failure time model. Simulations of dosing strategies showed that increasing both the starting and maximum doses could potentially shorten the time to stabilization, and thus, length of treatment and hospital stay. Given the hypothesis-generating nature of this analysis, the recommended dosing regimens should be tested prospectively to evaluate their benefits.

Recommended approaches for integration of population pharmacokinetic modelling with precision dosing in clinical practice

Current methods of dose determination have contributed to suboptimal and inequitable health outcomes in underrepresented patient populations. The persistent demand to individualise patient treatment, alongside increasing technological feasibility, is leading to a growing adoption of model-informed precision dosing (MIPD) at point of care. Population pharmacokinetic (popPK) modelling is a technique that supports treatment personalisation by characterising drug exposure in diverse patient groups. This publication addresses this important shift in clinical approach, by collating and summarising recommendations from literature. It seeks to provide standardised guidelines on best practices for the development of popPK models and their use in MIPD software tools, ensuring the safeguarding and optimisation of patient outcomes. Moreover, it consolidates guidance from key regulatory and advisory bodies on MIPD software deployment, as well as technical requirements for electronic health record integration. It also considers the future application and clinical impact of machine learning algorithms in popPK and MIPD. Ultimately, this publication aims to facilitate the incorporation of high-quality precision-dosing solutions into standard clinical workflows, thereby enhancing the effectiveness of individualised dose selection at point of care.

Model-informed dose optimization for prophylactic piperacillin-tazobactam in perioperative pediatric critically ill patients

Piperacillin/tazobactam (PTZ) is frequently prescribed during the perioperative period as prophylaxis in critically ill patients. Current international guidelines recommend that the pediatric intraoperative dosing regimen for PTZ be 90-112.5 mg/kg (80-100 mg/kg as piperacillin [PIP]) administered every 2 hours (Q2H). Concerns have been raised not only about the risk of nephrotoxicity due to elevated PIP exposure but also regarding the practicality of adhering to a 2-h dosing interval in clinical settings. To address these concerns, we employed population pharmacokinetic (PK) modeling and simulation approaches to optimize PTZ dosing regimens in pediatric intraoperative patients. PIP plasma concentration data were obtained from 34 patients using an opportunistic sampling strategy. A two-compartment model was found to adequately describe the PK data. Creatinine clearance was identified as a significant covariate on clearance. The inclusion of inter-occasion variability significantly improved model fit. Simulations across body weights of 10-70 kg and creatinine clearance of 20-130 mL/min/1.73 m2 demonstrated that 6-15 mg/kg Q2H, or a 10 mg/kg loading dose followed by 1.0-2.75 mg/kg/h continuous infusion would achieve free PIP concentrations being above the minimum inhibitory concentration (MIC) for 100% of the dosing interval (100% fT >1× MIC). For achieving 100% fT >4× MIC, 25-55 mg/kg Q2H or a 20 mg/kg loading dose followed by 3.25-9.25 mg/kg/h continuous infusion was derived. The model-informed simulations indicated that both lower Q2H doses and continuous infusion methods are clinically viable options and potentially resolve current clinical challenges during intraoperative dosing.

Semimechanistic Population PK/PD Modeling of Axatilimab in Healthy Participants and Patients With Solid Tumors or Chronic Graft-Versus-Host Disease

Axatilimab, a high-affinity humanized immunoglobulin G4 monoclonal antibody against colony-stimulating factor 1 receptor (CSF-1R), is approved for the treatment of chronic graft-versus-host disease (cGVHD), and under investigation for idiopathic pulmonary fibrosis and solid tumors. The population pharmacokinetics (PK) and pharmacodynamics (PD) of axatilimab were characterized in healthy participants and patients with solid tumors or cGVHD using data from four clinical studies with 325 participants, including 278 patients with cGVHD. The model structure reflected the mechanism of action of axatilimab: blocking CSF-1R signaling with axatilimab reduces the circulating levels of cells in the mononuclear phagocytic cell lineage (including nonclassical monocytic cells (NCMCs) and Kupffer cells), resulting in increases in circulating enzymes owing to reduced clearance by Kupffer cells. The structural model consisted of a two-compartment axatilimab PK model and turnover PD models for CSF-1, NCMCs, aspartate transaminase (AST), and creatine phosphokinase (CPK). Axatilimab PK and CSF-1 equations also included saturable clearance components to reflect the competitive binding of axatilimab and CSF-1 to CSF-1R. Covariate search was conducted with the conditional sampling use for the stepwise approach based on correlation tests (COSSAC) approach. Covariate effects on model parameters, steady-state axatilimab exposure, and NCMC concentrations were assessed. The final population PK/PD model was mathematically described with 6 ordinary differential equations and 39 model parameters. Among the 11 statistically significant covariates, one (body weight) and two (participant population type and baseline CPK) covariates affected axatilimab steady-state exposure and steady-state NCMC levels by > 20%, respectively. These results informed the axatilimab dosing strategy in patients with cGVHD.

Population pharmacokinetics and dosing optimization of imipenem in Chinese elderly patients

Objectives

To assess the pharmacokinetics and pharmacodynamics of imipenem in a retrospective cohort of hospitalized Chinese older patients.

Methods

A population pharmacokinetic (PPK) model was constructed utilizing a nonlinear mixed-effects modeling approach. The final model underwent evaluation through bootstrap resampling and visual predictive checks. Additionally, a population pharmacokinetic and pharmacodynamic analysis was conducted employing Monte Carlo simulations to investigate the impact of commonly used dosing regimens (0.25 g every 6 h, 0.5 g every 6 h, 0.5 g every 8 h, 1 g every 6 h, 1 g every 8 h, and 1 g every 12 h) on the likelihood of achieving the target therapeutic outcomes.

Results

A total of 370 observations available from 142 patients were incorporated in the PPK model. A two-compartment PPK model with linear elimination best predicted the imipenem plasma concentrations, with the creatinine clearance as a significant covariate of clearance. Typical estimates for clearance, inter-compartmental clearance, central and peripheral volume were 13.1 L·h-1, 11.9 L·h-1, 11.7 L, 29.3 L, respectively.

Conclusion

The pharmacokinetics of imipenem in elderly patients were effectively characterized by the established PPK model, which includes creatinine clearance as a key covariate. This research will enhance our understanding of imipenem elimination and support precision dosing in this patient demographic.

Population pharmacokinetic modelling of cetirizine concentrations in human breast milk-A contribution from the ConcePTION project

Cetirizine is an antihistamine commonly used to treat allergic rhinitis and other allergic conditions. Cetirizine is often prescribed to breastfeeding mothers although there is limited information on infant exposure via breast milk. The aim of this study was to develop a popPK model based on data from a lactation study to predict cetirizine breast milk concentrations and estimate the relative infant dose (RID) in a breastfed infant. A popPK model was developed in NONMEM on data from a human lactation study including 35 women using cetirizine or levocetirizine while breastfeeding. Serial samples of breast milk were collected (n = 205) and the cetirizine concentrations quantified using a validated LC-MS/MS method. A one-compartment model of cetirizine in breast milk was developed and employed to calculate the relative infant dose (RID). Covariates related to the maternal characteristics and breastfeeding patterns were evaluated in the model; only milk sampling pumping duration was found to be a significant covariate, with an increasing pumping duration leading to an increased apparent milk volume of distribution (Vm). The mean RID was 1.99% with the highest RID being 3.36% at Cmax. PopPK modelling could be used to estimate infant exposure to cetirizine via breast milk. The low predicted exposure in infants supports that cetirizine is compatible with breastfeeding.

Model-informed drug discovery and development approaches to inform clinical trial design and regulatory decisions: A primer for the MENA region

Model-Informed Drug Discovery and Development (MID3) represents a transformative approach in pharmaceutical research, integrating quantitative models to inform and optimize decision-making throughout the drug development process. This review explores the current applications, challenges, and future prospects of MID3 within the Middle East and North Africa (MENA) region. By leveraging local data and advanced computational techniques, MID3 has the potential to significantly enhance the efficiency and success rates of drug development tailored to regional health priorities. We discussed successful case studies of applying MID3 at different phases of drug development and clinical trials. Furthermore, we emphasized the critical need for MENA countries to embrace MID3 by investing in workforce training, aligning regulatory frameworks, and fostering collaborative research initiatives. This call to action underscores the importance of a robust MID3 ecosystem, urging policymakers, academic institutions, and industry stakeholders to prioritize and support its integration into the MENA region's healthcare.

Effect of poly(ε-caprolactone) microspheres on population pharmacokinetic/pharmacodynamic model of a simple coumarin

This study aims to evaluated the impact of poly(ε-caprolactone) (PCL) microspheres on the pharmacokinetics and pharmacodynamics (PopPK/PD) of 6-methylcoumarin (6MC). For this, PCL microspheres loaded with 6MC were prepared using the emulsification-evaporation method. Particle size, zeta potential, drug loading, and entrapment efficiency were characterised by dynamic light scattering and UV spectrophotometry. In vitro release and pharmacokinetics in Wistar rats were assessed for free and encapsulated 6MC. Anti-inflammatory activity was evaluated using the carrageenan-induced paw edoema model, with PopPK and PopPK/PD models developed. Microspheres showed diameters between 2.9 and 7.1 µm, zeta potentials of -10 to -15 mV, and drug loading of 0.24 mg/mg. Encapsulation efficiency was 45.5% to 75.9%. PopPK models showed enhanced absorption and distribution, with increased anti-inflammatory potency of encapsulated 6MC. PCL microspheres significantly improved the pharmacokinetic and pharmacodynamic profiles of 6MC, enhancing its therapeutic potential for lipophilic drugs.

An integrated population pharmacokinetic model of febuxostat in pediatric patients with hyperuricemia including gout and adult population of healthy subjects and patients with renal dysfunction

The study objective was to validate febuxostat dosage and administration in pediatric patients with hyperuricemia including gout, using an integrated population pharmacokinetic (PopPK) analysis in the Japanese population. Integrated PopPK analysis of febuxostat used a nonlinear mixed-effects modeling (NONMEM) program on plasma febuxostat concentration data for 2611 samples from Japanese pediatric patients with hyperuricemia including gout (n = 29) and from adult subjects who are healthy or have renal dysfunction (n = 113). We described febuxostat pharmacokinetics using an integrated PopPK model applicable both to pediatric patients and to the adult population. The covariates of body weight and eGFR were identified for CL/F and the covariate of fasted/fed status for bioavailability. The range of steady-state exposures (Cmax,ss and AUCτ,ss) for 5, 10, 20, and 30 mg of febuxostat in fed pediatric patients weighing 20 to 40 kg was within that for 10, 20, 40, and 60 mg of febuxostat in fed pediatric patients and adults weighing 40 to 120 kg. Post hoc estimates of CL/F, adjusted by body weight, differed little between pediatric patients and the adult population in the renal function categories of normal, mild dysfunction, and moderate dysfunction. We successfully validated the febuxostat dose that provided the same level of exposure in pediatric patients as in the adult population: half the adult dose for pediatric patients weighing <40 kg and the full adult dose for pediatric patients weighing ≥40 kg. As in adults, the results support the use of febuxostat without dose adjustment in pediatric patients who have mild to moderate renal dysfunction.

Population pharmacokinetics of oral fosfomycin calcium in healthy women

BACKGROUND

Fosfomycin is an antibiotic extensively used to treat uncomplicated urinary tract infections in women, and it is available in different salts and formulations. The European Medicines Agency (EMA) recommends further studies to characterize the pharmacokinetics of fosfomycin calcium for oral administration and to justify its dosage recommendation.

OBJECTIVES

A population pharmacokinetic model of fosfomycin calcium was developed after oral administration to healthy women.

METHODS

A clinical trial (a randomized, open-label, bioavailability study of single and multiple doses of 1000 mg capsules, single dose of 500 mg capsule and single dose of 250 mg/5 mL suspension of oral fosfomycin calcium under fasted conditions in healthy women volunteers, Code: PD7522.22, EudraCT: 2020-001664-28) was carried out at the Clinical Trial Unit, Araba University Hospital (Vitoria-Gasteiz, Spain). Twenty-four healthy women were included in the study, and plasma samples were collected at different times over a period of 24 h. The concentration-time data of fosfomycin in plasma were modelled by a population approach using a nonlinear mixed-effects modelling implemented by NONMEM 7.4 (ICON Clinical Research LLC, North Wales, PA, USA).

RESULTS

The pharmacokinetics of fosfomycin was best described by a two-compartment model. Creatinine clearance and body weight were identified as covariates for fosfomycin clearance and volume of distribution, respectively.

CONCLUSIONS

This study provides relevant information on the pharmacokinetic profile of fosfomycin in women after oral administration as calcium salt. This population model may be very useful for establishing dosage recommendations of fosfomycin calcium to treat urinary tract infections in women.

A sub-pharmacological test dose does not predict individual docetaxel exposure in prostate cancer patients

PURPOSE

Docetaxel is a cytotoxic drug used for first-line treatment of various malignancies. It has a narrow therapeutic index and shows wide interpatient variability in clearance and toxicity. Tools for individual dose optimization are needed to maximize efficacy and avoid toxicity.

METHODS

We performed a proof-of-concept study (EudraCT 2016-003785-77) to evaluate whether pharmacokinetics after a sub-pharmacological test dose of 1000 µg docetaxel (millidose) could be used to predict therapeutic dose exposure. Thirty prostate cancer patients eligible for treatment with docetaxel as part of routine clinical care were included. An intravenous docetaxel millidose was administered 1-7 days prior to therapeutic docetaxel. After both doses plasma docetaxel concentrations were measured by ultra- high performance liquid chromatography-tandem mass spectrometry. The docetaxel clearance was estimated with non-linear mixed effects modeling.

RESULTS

Geometric mean docetaxel clearance was 57.9 L/h (GCV 78.6%) after admission of a millidose and 40.3 L/h (GCV 60.7%) after admission of a therapeutic dose. The millidose and therapeutic dose in a single patient were not significantly correlated (Spearman's rho R = 0.02, P = 0.92).

CONCLUSION

Docetaxel pharmacokinetics at milli- and therapeutic dose level showed insufficient correlation for individual dose optimization. However, the clearance of a docetaxel millidose and full dose are within the same order of magnitude. Therefore, docetaxel millidose pharmacokinetics could potentially facilitate prediction of docetaxel pharmacokinetics at a population level in situations where therapeutic dose levels are impractical, such as pharmacokinetic drug-drug interaction studies or pediatric studies.

Clofazimine serum concentration and safety/efficacy in nontuberculous mycobacterial pulmonary disease treatment

BACKGROUND

Clofazimine (CFZ) has shown promising effects against Mycobacterium avium-intracellulare complex pulmonary disease (MAC-PD) and Mycobacterium abscessus species pulmonary disease (MABS-PD). However, the optimal CFZ dose remains unknown. We aimed to explore the relationship between steady-state CFZ concentration and its safety and efficacy in MAC-PD and MABS-PD.

METHODS

This prospective observational study focused on patients with MAC-PD and MABS-PD treated with CFZ (UMIN 000041053). To understand the safety and efficacy profile of CFZ and elucidate its optimal concentration, we analyzed CFZ-induced pigmentation grade, QTc interval, and culture conversion outcomes in relation to serum CFZ concentration using Student's t-test, a concentration-QTc model, and multivariable logistic regression analysis, respectively. In total, 64 patients (34 with MAC-PD; 30 with MABS-PD) were included.

RESULTS

The steady-state concentration of CFZ was higher in the moderate-to-severe pigmentation group than in the none-to-light pigmentation group (P < 0.001). At a CFZ concentration of 1 mg/L, the QTc interval was prolonged by 17.3 ms (95 % confidence interval [CI], 3.9-25.4) from baseline. Culture conversion was achieved in 33 (51.6 %) patients. The only significant predictor of culture conversion was surgery (adjusted odds ratio, 5.4; 95 % CI, 1.3-38.0). CFZ concentration and MIC of CFZ less than 0.25 mg/L were not associated with culture conversion in this study.

CONCLUSION

CFZ-induced pigmentation and QT interval prolongation are associated with serum CFZ concentrations. CFZ dosage may be optimized by monitoring serum CFZ concentration.

Pharmacokinetic Models of Tafenoquine: Insights for Optimal Malaria Treatment Strategies

Tafenoquine (TQ) is a new 8-aminoquinoline antimalarial drug developed by the US Army for Plasmodium vivax malaria treatment. Modeling and simulation are essential tools for drug development and improving rationality in pharmacotherapy, and different modeling approaches are used. This study aims to summarize and explore the pharmacokinetic (PK) models available for tafenoquine in the literature. An integrative methodology was used to collect and review published data. Fifteen articles were identified using three modeling approaches: non-compartmental analysis (NCA), population pharmacokinetic analysis (popPK), and pharmacokinetic/pharmacodynamic analysis (PK/PD). An NCA was mainly used to describe the PK profile of TQ and to compare its PK profile alone to those obtained in association with other drugs. PopPK was used to assess TQ population PK parameters, covariates' impact, and dose selection. PK/PD helped understand the relationship between TQ concentrations, some adverse events common for 8-aminoquilones, and the efficacy assessment for Plasmodium falciparum. In summary, pharmacokinetic models were widely used during TQ development. However, there is still a need for different modeling approaches to support further therapeutic questions, such as treatment for special populations and potential drug-drug interactions.

More Than a Gut Feeling─A Combination of Physiologically Driven Dissolution and Pharmacokinetic Modeling as a Tool for Understanding Human Gastric Motility

In vivo studies of formulation performance with in vitro and/or in silico simulations are often limited by significant gaps in our knowledge of the interaction between administered dosage forms and the human gastrointestinal tract. This work presents a novel approach for the investigation of gastric motility influence on dosage form performance, by combining biopredictive dissolution tests in an innovative PhysioCell apparatus with mechanistic physiology-based pharmacokinetic modeling. The methodology was based on the pharmacokinetic data from a large (n = 118) cohort of healthy volunteers who ingested a capsule containing a highly soluble and rapidly absorbed drug under fasted conditions. The developed dissolution tests included biorelevant media, varied fluid flows, and mechanical stress events of physiological timing and intensity. The dissolution results were used as inputs for pharmacokinetic modeling that led to the deduction of five patterns of gastric motility and their prevalence in the studied population. As these patterns significantly influenced the observed pharmacokinetic profiles, the proposed methodology is potentially useful to other in vitro-in vivo predictions involving immediate-release oral dosage forms.

Population pharmacokinetics of daptomycin in critically ill patients receiving extracorporeal membrane oxygenation

BACKGROUND

Daptomycin is widely used in critically ill patients for Gram-positive bacterial infections. Extracorporeal membrane oxygenation (ECMO) is increasingly used in this population and can potentially alter the pharmacokinetic (PK) behaviour of antibiotics. However, the effect of ECMO has not been evaluated in daptomycin. Our study aims to explore the effect of ECMO on daptomycin in critically ill patients through population pharmacokinetic (PopPK) analysis and to determine optimal dosage regimens based on both efficacy and safety considerations.

METHODS

A prospective, open-label PK study was carried out in critically ill patients with or without ECMO. The total concentration of daptomycin was determined by UPLC-MS/MS. NONMEM was used for PopPK analysis and Monte Carlo simulations.

RESULTS

Two hundred and ninety-three plasma samples were collected from 36 critically ill patients, 24 of whom received ECMO support. A two-compartment model with first-order elimination can best describe the PK of daptomycin. Creatinine clearance (CLCR) significantly affects the clearance of daptomycin while ECMO has no significant effect on the PK parameters. Monte Carlo simulations showed that, when the MICs for bacteria are  ≥1 mg/L, the currently recommended dosage regimen is insufficient for critically ill patients with CLCR > 30 mL/min. Our simulations suggest 10 mg/kg for patients with CLCR between 30 and 90 mL/min, and 12 mg/kg for patients with CLCR higher than 90 mL/min.

CONCLUSIONS

This is the first PopPK model of daptomycin in ECMO patients. Optimal dosage regimens considering efficacy, safety, and pathogens were provided for critical patients based on pharmacokinetic-pharmacodynamic analysis.

Application of machine learning techniques in population pharmacokinetics/pharmacodynamics modeling

Population pharmacokinetics/pharmacodynamics (pop-PK/PD) consolidates pharmacokinetic and pharmacodynamic data from many subjects to understand inter- and intra-individual variability due to patient backgrounds, including disease state and genetics. The typical workflow in pop-PK/PD analysis involves the determination of the structure model, selection of the error model, analysis based on the base model, covariate modeling, and validation of the final model. Machine learning is gaining considerable attention in the medical and various fields because, in contrast to traditional modeling, which often assumes linear or predefined relationships, machine learning modeling learns directly from data and accommodates complex patterns. Machine learning has demonstrated excellent capabilities for prescreening covariates and developing predictive models. This review introduces various applications of machine learning techniques in pop-PK/PD research.

Population pharmacokinetics of everolimus in renal transplant recipients receiving long-term multiple immunosuppressive therapy

Everolimus is used for immunosuppression after renal transplantation. This study aimed to develop a population pharmacokinetic (PopPK) model of everolimus using therapeutic drug monitoring (TDM) data of patients under long-term multiple immunosuppressive therapy, including tacrolimus. To develop the model, 185 renal transplant recipients with 3358 everolimus blood concentrations during a median postoperative period of 35.3 months were included. The PopPK model is described as a one-compartment model with first-order absorption. The population mean of apparent clearance is 8.92 L/h (relative standard error = 3.6%), and this negatively correlated with the dose-normalized concentration (C/D) of tacrolimus and hematocrit value, and positively correlated with a daily dose of everolimus (i.e. TDM effect). The usefulness of dose adjustment using the final popPK model was assessed by a simulation study. The ratio of the first trough measurement within the therapeutic range of 3-8 ng/mL increased from 69.8% in the original dose to 87.9% in the individual dose calculated by the final PopPK model. The tacrolimus C/D ratio before initiating everolimus therapy and the hematocrit value were useful to estimate the initial dose of everolimus and can improve the safety and effectiveness of immunosuppressive therapy involving everolimus.

Population pharmacokinetics of oxcarbazepine active metabolite in Chinese paediatric patients with epilepsy: Model-based dose optimization

The pharmacological activity of oxcarbazepine (OXC) is primarily exerted through its active 10-monohydroxy metabolite (MHD). Nonetheless, there is limited pharmacokinetic information available regarding paediatric patients with epilepsy treated with OXC, especially in infants and toddlers. Concurrently, this drug exhibits substantial variability in pharmacokinetics and therapeutic response across different individuals. We aimed to develop a model to quantitatively investigate factors that affect MHD pharmacokinetics to formulate a dosage guideline for OXC in Chinese paediatric patients. A total of 297 MHD trough concentrations were obtained from 287 epileptic children. Six body weight (BW)-based allometric models were used for population pharmacokinetic modelling, while investigating the impact of other covariates on the apparent clearance. The one-compartment model and age cut-off model for the apparent clearance (CL/F) were established to describe the pharmacokinetics of MHD. The probability to obtain target trough concentration ranges (TTCRs) of MHD between 3 and 35 mg/L was determined by Monte Carlo simulations for doses ranging from 8 to 90 mg/kg/day. A new dose optimization strategy combining the dosage guidelines and Bayesian method provides a tailored approach for Chinese paediatric epileptic patients based on their individual BW and desired TTCRs of MHD, and also supports current dose recommendations, with the exception of children weighing ≤5 kg.

Pharmacokinetic Modeling and Simulation with Pharmacogenetic Insights Support the Relevance of Therapeutic Drug Monitoring for Myeloablative Busulfan Dosing in Adult HSCT

Therapeutic drug monitoring (TDM) of busulfan (Bu) is well-established in pediatric hematopoietic stem cell transplantation (HSCT), but its use in adults is limited due to a lack of clear recommendations and scarcity of evidence regarding its utility. GSTA1 promoter variants are reported to affect Bu clearance in both adults and pediatric patients. This study aimed to evaluate the value of preemptive genotyping GSTA1 and body composition (obesity) in individualizing Bu dosing in adults, through pharmacokinetic (PK) modeling and simulations. A population pharmacokinetic (PopPK) model was developed and validated with data from 60 adults who underwent HSCT. Simulations assessed different dosing scenarios based on body size metrics and GSTA1 genotypes. Due to the limited number of obese patients in the cohort, the effect of obesity on Bu pharmacokinetics (PK) was evaluated in silico using a physiologically-based pharmacokinetic (PBPK) model and relevant virtual populations from Simcyp software. Patients with at least 1 GSTA1*B haplotype had 17% lower clearance on average. PopPK simulations indicated that adjusting doses based on genotype increased the probability of achieving the target exposure (3.7 to 5.5 mg.h/L) from 53% to 60 % in GSTA1*A homozygous patients, and from 50% to 61% in *B carriers. Still, Approximately 40% of patients would not achieve this therapeutic window without TDM. A 2-sample optimal design was validated for routine model-based Bu first dose AUC0-∞ estimation, and the model was implemented in the Tucuxi user-friendly TDM software. PBPK simulations confirmed body surface area-based doses of 29 to 31 mg/m2/6h as the most appropriate, regardless of obesity status. This study emphasizes the importance of individualized Bu dosing strategies in adults to achieve therapeutic targets. Preemptive genotyping alone may not have a significant clinical impact, and routine TDM may be necessary for optimal transplantation outcomes.

Pharmacokinetics and pharmacodynamics of cannabis-based medicine in a patient population included in a randomized, placebo-controlled, clinical trial

Information on the pharmacokinetics (PK) and pharmacodynamics (PD) of orally administered cannabis-based medicine (CBM) in capsule formulation in patient populations is sparse. In this exploratory study, we aimed to evaluate the PK and PD in a probable steady state of CBM in neuropathic pain and spasticity in a population of patients with multiple sclerosis (MS). Of 134 patients participating in a randomized, double-blinded, placebo-controlled, trial, 23 patients with MS (17 female) mean age 52 years (range 21-67) were enrolled in this substudy. They received oral capsules containing Δ9 -tetrahydrocannabinol (THC, n = 4), cannabidiol (CBD, n = 6), a combination (THC&CBD, n = 4), or placebo (n = 9). Maximum doses were 22.5 mg (THC) and 45 mg (CBD) a day divided into three administrations. PD parameters were evaluated for pain and spasticity. Blood samples were analyzed using an ultra-high-performance liquid chromatography-tandem mass spectrometer after protein precipitation and phospholipid removal. PK parameters were estimated using computerized modeling. The variation in daily dose and PK between individuals was considerable in a steady state, yet comparable with previous reports from healthy controls. Based on a simulation of the best model, the estimated PK parameters (mean) for THC (5 mg) were Cmax 1.21 ng/mL, Tmax 2.68 h, and half-life 2.75 h, and for CBD (10 mg) were Cmax 2.67 ng/mL, Tmax 0.10 h, and half-life 4.95 h, respectively. No effect was found on the PD parameters, but the placebo response was considerable. More immediate adverse events were registered in the active treatment groups compared with the placebo group.

Population pharmacokinetics and model-based dosing optimization of teicoplanin in elderly critically ill patients with pneumonia

PURPOSE

To evaluate the population pharmacokinetics and pharmacodynamics of teicoplanin in elderly critically ill patients with pneumonia for optimal dosages.

METHODS

Fifteen critically ill patients (9 men) ≥ 60 years received teicoplanin 6 mg/kg for three doses followed by standard maintenance doses (6 mg/kg q24h) with renal dosing adjustment. Serial plasma samples from all patients were analyzed simultaneously by population pharmacokinetic modeling using NONMEM. Probability of target attainment (PTA) was calculated through Monte Carlo simulations for various dosing regimens to achieve adequate systemic exposures.

RESULTS

The median (interquartile range, IQR) age, body mass index, and creatinine clearance (CrCl) was 75 (64-78) years, 22.5 (20.8-25.4) kg/m2, and 64 (47-106) mL/min, respectively. The median (IQR) peak and trough concentration was 46.5 (42.7-51.0) and 8.7 (7.2-9.5) mg/L. The population pharmacokinetic model showed slower clearance (CL) and larger peripheral volume of distribution (V2) in patients with reduced CrCl: CL (L/h) = 0.629 × (CrCl/64)0.656, V2 (L) = 55.7 × (CrCl/64)-0.665. Model-based simulations showed PTAs ≥85% only for higher-dose regimens (12 mg/kg) up to an MIC of 0.5 mg/L.

CONCLUSIONS

Standard teicoplanin dosages for pneumonia may provide inadequate systemic exposures in elderly critically ill patients. High-dose regimens should be considered as empiric therapy or for less susceptible pathogens.

Tutorial on model selection and validation of model input into precision dosing software for model-informed precision dosing

There has been rising interest in using model-informed precision dosing to provide personalized medicine to patients at the bedside. This methodology utilizes population pharmacokinetic models, measured drug concentrations from individual patients, pharmacodynamic biomarkers, and Bayesian estimation to estimate pharmacokinetic parameters and predict concentration-time profiles in individual patients. Using these individualized parameter estimates and simulated drug exposure, dosing recommendations can be generated to maximize target attainment to improve beneficial effect and minimize toxicity. However, the accuracy of the output from this evaluation is highly dependent on the population pharmacokinetic model selected. This tutorial provides a comprehensive approach to evaluating, selecting, and validating a model for input and implementation into a model-informed precision dosing program. A step-by-step outline to validate successful implementation into a precision dosing tool is described using the clinical software platforms Edsim++ and MwPharm++ as examples.

Dosage optimization of tacrolimus based on the glucocorticoid dose and pharmacogenetics in adult patients with systemic lupus erythematosus

BACKGROUND

The purpose of the study was to develop a genotype-incorporated population pharmacokinetic (PPK) model of tacrolimus (TAC) in adults with systemic lupus erythematosus (SLE) to investigate the factors influencing TAC pharmacokinetics and to develop an individualized dosing regimen based on the model. In addition, a non-genotype-incorporated model was also established to assess its predictive performance compared to the genotype-incorporated model.

METHODS

A total of 365 trough concentrations from 133 adult SLE patients treated with TAC were collected to develop a genotype-incorporated PPK model and a non-genotype-incorporated PPK model of TAC using a nonlinear mixed-effects model (NONMEM). External validation of the two models was performed using data from an additional 29 patients. Goodness-of-fit diagnostic plots, bootstrap method, and normalized predictive distribution error test were used to validate the predictive performance and stability of the final models. The goodness-of-fit of the two final models was compared using the Akaike information criterion (AIC). The dosing regimen was optimized using Monte Carlo simulations based on the developed optimal model.

RESULTS

The typical value of the apparent clearance (CL/F) of TAC estimated in the final genotype-incorporated model was 14.3 L h-1 with inter-individual variability of 27.6%. CYP3A5 polymorphism and coadministered medication were significant factors affecting TAC-CL/F. CYP3A5 rs776746 GG genotype carriers had only 77.3% of the TAC-CL/F of AA or AG genotype carriers. Omeprazole reduced TAC-CL/F by 3.7 L h-1 when combined with TAC, while TAC-CL/F increased nonlinearly as glucocorticoid dose increased. Similar findings were demonstrated in the non-genotype-incorporated PPK model. Comparing these two models, the genotype-incorporated PPK model was superior to the non-genotype-incorporated PPK model (AIC = 643.19 vs. 657.425). Monte Carlo simulation based on the genotype-incorporated PPK model indicated that CYP3A5 rs776746 AA or AG genotype carriers required a 1/2-1 fold higher dose of TAC than GG genotype carriers to achieve the target concentration. And as the daily dose of prednisone increases, the dose of TAC required to reach the target concentration increases appropriately.

CONCLUSIONS

We developed the first pharmacogenetic-based PPK model of TAC in adult patients with SLE and proposed a dosing regimen based on glucocorticoid dose and CYP3A5 genotype according to the model, which could facilitate individualized dosing for TAC.

Quantitative risk-benefit profiles of oral contraceptives, insulin sensitizers and antiandrogens for women with polycystic ovary syndrome: A model-based meta-analysis

Oral contraceptives (OCs), insulin sensitizers, and antiandrogens (AAs), alone or in combination, are commonly used for treating non-fertility indications in polycystic ovary syndrome (PCOS). However, unclear risk-benefit profiles jeopardize their appropriate clinical applications. This study aimed to quantitatively evaluate the effects of the aforementioned medications and to compare their risk-benefit profiles. Randomized controlled trials published until 14th March 2022 were searched in PubMed and Embase. A model-based meta-analysis was developed to examine the time-effect profiles of each medication. The maximal percentage change of the effect (Emax) and time to achieve half of Emax (T50) were estimated. Primary outcomes included menstruation, hirsutism score, free androgen index (FAI), body mass index (BMI), insulin sensitivity, and lipid profiles. Overall, 200 studies (9,685 patients and 385 arms) were identified for modeling. OCs performed exceptionally well in improving menstruation (Emax: 149%; T50: 7.44 weeks), hirsutism score (Emax: 66.2%; T50: 26.2 weeks), and FAI (Emax: 75.7%; T50: 0.51 weeks). However, OCs elevated the triglyceride (TG) level (Emax: 12.6%; T50:1.19 weeks). After 12-week OC treatment, the TG level of approximately 30% of patients, whose baselines were normal, exceeded the reference limit. This suggested that OC-induced dyslipidemia should be routinely monitored. The maximal BMI-lowering effect of metformin was similar to that of placebo (Emax: 3.80%); however, metformin had a shorter T50 (6.67 weeks versus 12.9 weeks). Further, active lifestyle intervention plus placebo significantly decreased BMI (Emax: 8.78%). Adding metformin to active lifestyle intervention accelerated the BMI-lowering effect within 24 weeks, whereas with the extension of this addition beyond 24 weeks, BMI did not reduce further, which indicated that benefits were limited from this prolonged addition. AAs were less potent in reducing hirsutism score (Emax: 40.2% versus 66.2%) and FAI (Emax: 34.5% versus 75.7%) compared to OCs. OC plus metformin combined OC-derived androgen-suppressing effects and metformin-derived insulin-sensitizing effects, and partially relieved the OC-induced TG increase (Emax: 9.76%). Baseline dependency was found in most clinical responses, implying that pharmacotherapies tailored based on baselines achieved more clinical improvements. This study presents new quantitative evidence on pharmacotherapies for PCOS. Currently, long-term risk-benefit profiles and emerging therapies are inadequately reported and require more further research.

Pharmacokinetic analysis across studies to drive knowledge-integration: A tutorial on individual patient data meta-analysis (IPDMA)

Answering challenging questions in drug development sometimes requires pharmacokinetic (PK) data analysis across different studies, for example, to characterize PKs across diverse regions or populations, or to increase statistical power for subpopulations by combining smaller size trials. Given the growing interest in data sharing and advanced computational methods, knowledge integration based on multiple data sources is increasingly applied in the context of model-informed drug discovery and development. A powerful analysis method is the individual patient data meta-analysis (IPDMA), leveraging systematic review of databases and literature, with the most detailed data type of the individual patient, and quantitative modeling of the PK processes, including capturing heterogeneity of variance between studies. The methodology that should be used in IPDMA in the context of population PK analysis is summarized in this tutorial, highlighting areas of special attention compared to standard PK modeling, including hierarchical nested variability terms for interstudy variability, and handling between-assay differences in limits of quantification within a single analysis. This tutorial is intended for any pharmacological modeler who is interested in performing an integrated analysis of PK data across different studies in a systematic and thorough manner, to answer questions that transcend individual primary studies.

Model-Informed Precision Dosing Guidance of Ethosuximide Developed from a Randomized Controlled Clinical Trial of Childhood Absence Epilepsy

Ethosuximide was identified as the optimal option for new-onset childhood absence epilepsy (CAE) in a randomized, two-phase dose escalation comparative effectiveness trial of ethosuximide, lamotrigine, and valproic acid. However, 47% of ethosuximide initial monotherapy participants experienced short-term treatment failure. This study aimed to characterize the initial monotherapy ethosuximide exposure-response relationship and to propose model-informed precision dosing guidance. Dose titration occurred over a 16-20-week period until patients experienced seizure freedom or intolerable side effects. Subjects with initial monotherapy failure were randomized to one of the other two medications and dose escalation was repeated. A population pharmacokinetic model was created using plasma concentration data (n = 1,320), collected at 4-week intervals from 211 unique participants during both the initial and second monotherapy phases. A logistic regression analysis was performed on the initial monotherapy cohort (n = 103) with complete exposure-response data. Eighty-four participants achieved seizure freedom with a wide range of ethosuximide area under the curves (AUC) ranging from 420 to 2,420 μg·h/mL. AUC exposure estimates for achieving a 50% and 75% probability of seizure freedom were 1,027 and 1,489 μg·h/mL, respectively, whereas the corresponding cumulative frequency of intolerable adverse events was 11% and 16%. Monte Carlo Simulation indicated a daily dose of 40 and 55 mg/kg to achieve 50% and 75% probability of seizure freedom in the overall population, respectively. We identified the need for adjusted mg/kg dosing in different body weight cohorts. This ethosuximide proposed model-informed precision dosing guidance to achieve seizure freedom carries promise to optimize initial monotherapy success for patients with CAE.

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 subcutaneous alemtuzumab in kidney transplantation

AIM

Alemtuzumab is a monoclonal antibody used as induction immunosuppressive therapy in kidney transplantation. It targets CD52 on lymphocytes, inducing profound immune cell depletion upon administration. Owing to its off-label status in kidney transplantation, its pharmacokinetic characteristics are largely unknown in this setting, and its current fixed dosing algorithm originates from other populations. We developed a population pharmacokinetic model for alemtuzumab in kidney transplant recipients and investigated the potential of personalized alemtuzumab therapy.

METHODS

In total, 362 pharmacokinetic observations drawn 0-165 days after transplantation were available from 61 adult kidney transplant recipients who received two consecutive doses of 15 mg alemtuzumab subcutaneously. A population pharmacokinetic model was developed using nonlinear mixed-effects modelling and applied to simulate various dosing regimens.

RESULTS

The alemtuzumab concentration-time data were best described by a two-compartmental model with first-order absorption and parallel first-order and time-varying concentration-dependent elimination, with between-subject variability on the first-order elimination (39.6%) and central distribution volume (39.6%). Alemtuzumab pharmacokinetics varied with body size, rendering lighter individuals exposed to lympholytic alemtuzumab concentrations (>0.1 mg/L) for prolonged durations as compared to their heavier peers. This between-subject variability could be reduced through lean bodyweight-adjusted dosing, showing a twofold to threefold reduction in the slope of the median alemtuzumab exposure over the bodyweight range.

CONCLUSION

Alemtuzumab displays substantial pharmacokinetic variability in kidney transplant recipients, which may warrant a personalized treatment strategy. Lean bodyweight-adjusted dosing poses an option for individualized dosing, but further evaluation of its potential clinical benefit is warranted.

Rationale for the selection of dual primary endpoints in prevention studies of cognitively unimpaired individuals at genetic risk for developing symptoms of Alzheimer's disease

BACKGROUND

There is a critical need for novel primary endpoints designed to detect early and subtle changes in cognition in clinical trials targeting the asymptomatic (preclinical) phase of Alzheimer's disease (AD). The Alzheimer's Prevention Initiative (API) Generation Program, conducted in cognitively unimpaired individuals at risk of developing AD (e.g., enriched by the apolipoprotein E (APOE) genotype), used a novel dual primary endpoints approach, whereby demonstration of treatment effect in one of the two endpoints is sufficient for trial success. The two primary endpoints were (1) time to event (TTE)-with an event defined as a diagnosis of mild cognitive impairment (MCI) due to AD and/or dementia due to AD-and (2) change from baseline to month 60 in the API Preclinical Composite Cognitive (APCC) test score.

METHODS

Historical observational data from three sources were used to fit models to describe the TTE and the longitudinal APCC decline, both in people who do and do not progress to MCI or dementia due to AD. Clinical endpoints were simulated based on the TTE and APCC models to assess the performance of the dual endpoints versus each of the two single endpoints, with the selected treatment effect ranging from a hazard ratio (HR) of 0.60 (40% risk reduction) to 1 (no effect).

RESULTS

A Weibull model was selected for TTE, and power and linear models were selected to describe the APCC score for progressors and non-progressors, respectively. Derived effect sizes in terms of reduction of the APCC change from baseline to year 5 were low (0.186 for HR = 0.67). The power for the APCC alone was consistently lower compared to the power of TTE alone (58% [APCC] vs 84% [TTE] for HR = 0.67). Also, the overall power was higher for the 80%/20% distribution (82%) of the family-wise type 1 error rate (alpha) between TTE and APCC compared to 20%/80% (74%).

CONCLUSIONS

Dual endpoints including TTE and a measure of cognitive decline perform better than the cognitive decline measure as a single primary endpoint in a cognitively unimpaired population at risk of AD (based on the APOE genotype). Clinical trials in this population, however, need to be large, include older age, and have a long follow-up period of at least 5 years to be able to detect treatment effects.

Population pharmacokinetic analysis reveals no impact of aprepitant on the pharmacokinetics of ifosfamide, 2-dechloroifosfamide, and 3-dechloroifosfamide

PURPOSE

Several case reports and retrospective series have clearly pointed to the role of aprepitant, an antiemetic drug, in the development of encephalopathy when used with ifosfamide. Described as an inhibitor of several CYP metabolic pathways, aprepitant is suspected of drug-drug-interaction on ifosfamide pharmacokinetics. The pharmacokinetics of ifosfamide and two of its metabolites (2-dechloroifosfamide and 3-dechloroifosfamide) was studied in patients with soft tissue sarcomas to evaluate the impact of aprepitant administration.

METHODS

A population pharmacokinetic approach was applied to analyze data obtained in 42 patients at cycle 1 (without aprepitant) and cycle 2 (with aprepitant for 34 of them).

RESULTS

A previously published pharmacokinetic model including a time-dependency process well fit the data. Aprepitant had no impact on ifosfamide or its two metabolite pharmacokinetic parameters.

CONCLUSION

This study suggests that aprepitant does not lead to a significant modification of ifosfamide metabolization, even though other metabolites such as 4 hydroxyifosfamide and chloroacetaldehyde were not monitored in this study.

Real world population pharmacokinetic study in children and young adults with inflammatory bowel disease discovers novel blood and stool microbial predictors of vedolizumab clearance

BACKGROUND

Vedolizumab for inflammatory bowel disease (IBD) is often intensified based on distinct pharmacokinetics in children. Prior adult-specific population pharmacokinetic models have identified limited covariates of drug clearance.

AIMS

To establish a population pharmacokinetic model for children and young adults to identify novel covariates of drug clearance to better account for paediatric-specific inter-patient variability in vedolizumab pharmacokinetics; a key secondary exploratory aim was to identify microbial signatures of pharmacokinetic outcomes in a subset of patients.

METHODS

The study included data from 463 observed vedolizumab concentrations (59 peaks and 404 troughs) from 74 patients with IBD (52 with Crohn's disease and 22 with ulcerative colitis or unclassified IBD, median age 16 years). Pharmacokinetic analysis was conducted with non-linear mixed effects modelling. For the evaluation of the exposure-response relationship, clinical outcomes were evaluated by trough levels, clearance and vedolizumab exposure. Whole-genome metagenomic sequencing was conducted at baseline and week 2.

RESULTS

A two-compartment population pharmacokinetic model was identified with a clear correlation between CL and weight, erythrocyte sedimentation rate, and hypoalbuminemia. Trough concentrations before infusion 3 (37 μg/ml) and before infusion 4 (20 μg/ml) best predicted steroid-free clinical remission at infusion 4. Using faecal metagenomics, we identified an early (baseline and week 2) abundance of butyrate-producing species and pathways that were associated with an infusion 4 trough concentration >20 μg/ml.

CONCLUSIONS

This novel paediatric vedolizumab pharmacokinetic model could inform precision dosing. While additional studies are needed, an abundance of faecal butyrate producers is associated with early response to vedolizumab, suggesting that microbial analysis may be beneficial to biological selection.

Population pharmacokinetics and dosing regimen optimization of levetiracetam in epilepsy during pregnancy

AIMS

The pharmacokinetics of levetiracetam (LEV) significantly changed during pregnancy. It is a great challenge to predict the adjusted doses of LEV to reach the preconception target concentrations. This study aimed to establish a population pharmacokinetic model of LEV in women with epilepsy (WWE) during pregnancy to analyse the factors of pharmacokinetic variability and to develop a model-based individualized dosing regimen.

METHODS

A total of 166 concentration-time points from 37 WWE during pregnancy treated with LEV were collected to analyse LEV pharmacokinetics with nonlinear mixed-effects modelling. The dosing regimen was optimized by Monte Carlo simulations based on the final model.

RESULTS

The LEV pharmacokinetics in pregnant WWE were best described by a 1-compartment model of first-order absorption and elimination. The population typical value of apparent clearance (CL/F) in the final model was estimated to be 3.82 L/h (95% confidence interval 3.283-4.357 L/h) with a relative standard error of 7.2%. Both total body weight (TBW) and trimester of pregnancy were significantly associated with LEV-CL/F during pregnancy; LEV-CL/F increased by 42.72% when TBW increased from 55 to 65 kg from the first trimester to the second trimester. Monte Carlo simulations showed that dosing regimens for LEV should be individualized based on the patient's TBW and trimester of pregnancy to maximize the likelihood of achieving the therapeutic range.

CONCLUSION

This first population pharmacokinetic study of LEV in WWE during pregnancy supports the use of a weight-based and pregnancy-based dosing regimen and can lay a foundation for further optimizing the individualized dosing regimens.

Population Pharmacokinetics of CMAB008 (an Infliximab Biosimilar) and Remicade® in Healthy Subjects and Patients with Moderately to Severely Active Rheumatoid Arthritis

INTRODUCTION

CMAB008 is a monoclonal antibody developed as a biosimilar to infliximab (Remicade^®, Janssen). The pharmacokinetic characteristics of CMAB008 and Remicade^® in healthy subjects and patients with moderately to severely active rheumatoid arthritis (RA) were investigated using a population modeling approach, and the pharmacokinetic similarity of CMAB008 to Remicade^® was assessed.

METHODS

The population pharmacokinetic model was developed on the basis of intensive pharmacokinetic data from a phase 1 study in healthy male subjects and combined intensive and sparse pharmacokinetic data from a phase 3 study in patients with RA.

RESULTS

A two-compartment model with first-order elimination adequately described CMAB008 and Remicade^® concentration data in healthy subjects and patients with RA. The analysis of covariates identified anti-drug antibody (ADA), neutralizing antibody (NAB), real-time body weight (BWT), and real-time albumin (ALB) as significant covariates on clearance, and BWT was also a significant covariate for the central volume of distribution. The treatment type (CMAB008 versus Remicade^®) and the study population (healthy subjects versus patients with RA) were not identified as significant covariates on the pharmacokinetics of infliximab, demonstrating pharmacokinetic similarity between CMAB008 and Remicade^® in both study populations. The effect of BWT and ALB changes on exposures to infliximab was within the acceptable range, suggesting that the 3 mg/kg regimen is appropriate in clinical practice for patients with RA and BTW and ALB distribution within the range evaluated in the current analysis.

CONCLUSIONS

The pharmacokinetic characteristics were similar between CMAB008 and Remicade^® in healthy subjects and patients with RA. CMAB008 can be considered bioequivalent to Remicade^®.

CLINICAL TRIAL REGISTRATION

ClinicalTrials.gov identifiers NCT04779892, NCT03478111.

Determination of vancomycin exposure target and individualized dosing recommendations for critically ill patients undergoing continuous renal replacement therapy

STUDY OBJECTIVE

Few studies have been conducted to quantify the exposure target of vancomycin in intensive care unit (ICU) patients undergoing continuous renal replacement therapy (CRRT) and provide optimized dosage regimens. We aimed to determine vancomycin exposure target and dosing recommendations using data from an open database in critically ill patients undergoing CRRT.

DESIGN

A retrospective observational cohort study.

DATA SOURCE

A large public database.

PATIENTS

The adult patients who received intravenous vancomycin and CRRT treatment in the database between 2017 and 2019 were reviewed to determine eligibility. A total of 180 patients with 1186 observations were included in the population pharmacokinetic (PPK) model development. The clinical efficacy of vancomycin was analyzed in 159 eligible patients.

METHODS

A PPK model was developed to estimate individual pharmacokinetic (PK) parameters. The area under the concentration-time curve (AUC) was estimated by a Bayesian approach based on individual vancomycin concentrations. Multivariate logistic regression analyses were performed to identify the factors of clinical outcomes. Threshold of vancomycin exposure in predicting efficacy was identified via receiver operating characteristic (ROC) curve. Dosing recommendations were designed using Monte Carlo Simulations (MCS) based on the optimized exposure target.

MEASUREMENTS AND MAIN RESULTS

On covariate analysis, CRRT intensity significantly affected vancomycin PK. The AUC above 427 mg*h/L was the only significant predictor of clinical efficacy (adjusted odds ratio (aOR): 1.008, 95% confidence interval (CI): 1.004-1.011, p = 0.000). MCS indicated that vancomycin dosage regimens of 5 mg/kg q12h or 7.5 mg/kg q12h were recommended for patients with CRRT intensities of 20-25 mL/kg/h or 25.1-45 mL/kg/h, respectively.

CONCLUSIONS

An AUC threshold of 427 mg*h/L (assuming the minimal inhibitory concentration (MIC) = 1 mg/L) was a recommended efficacy exposure target of vancomycin for critically ill patients undergoing CRRT. Vancomycin 5-7.5 mg/kg q12h is recommended as the initial dosage regimens for ICU patients undergoing CRRT.

Population PK and Semimechanistic PK/PD Modeling and Simulation of Relugolix Effects on Testosterone Suppression in Men with Prostate Cancer

Relugolix, the first orally active, nonpeptide gonadotropin-releasing hormone receptor antagonist, is approved in the United States and the European Union for the treatment of adult patients with advanced prostate cancer. The recommended dosing regimen is a 360-mg loading dose followed by a 120-mg daily dose. Relugolix and testosterone concentration data and clinical information from two phase I studies, two phase II studies, and the phase III safety and efficacy study (HERO) were used to develop a population pharmacokinetic (PopPK) model and a semimechanistic population pharmacokinetic/pharmacodynamic (PopPK/PD) model that characterized relugolix exposure and its relationship to testosterone concentrations. Age, body weight, and Black/African American race had at most minimal effects on relugolix exposure or testosterone concentrations with no clinical relevance. Simulations using the PopPK/PD model confirmed the recommended dosing regimen of relugolix, with the median simulated testosterone concentrations predicted to achieve castration levels (< 50 ng/dL) and profound castration levels (< 20 ng/dL) by day 2 and day 9, respectively, and demonstrated that 97.3% and 85.5% of the patients remained at castration levels (< 50 ng/dL) upon temporary interruption of treatment for 7 days and 14 days, respectively. Collectively, simulations based on the PopPK and PopPK/PD models were consistent with actual data from clinical studies, reflecting the high predictiveness of the models and supporting the reliability of model-based simulations. These models can be used to provide guidance regarding dosing recommendations under various circumstances (e.g., temporary interruption of treatment, if needed) for relugolix.

A Survey of Population Pharmacokinetic Reports Submitted to the USFDA: An Analysis of Common Issues in NDA and BLA from 2012 to 2021

BACKGROUND AND OBJECTIVE

Population pharmacokinetic (PopPK) analysis is one of the important components of regulatory submission. The purpose of this study was to survey PopPK analysis reports in new drug applications (NDAs) and biological licensing applications (BLAs) submitted to the US Food and Drug Administration (FDA) and to retrieve the information regarding PopPK-related issues in the US FDA review reports.

METHODS

We surveyed NDAs and BLAs over the period from 2012 to 2022 from Drug@FDA databases and extracted the review reports from the website. We explored the issues identified in these reports and sorted them into four categories: (i) data used for PopPK analysis, (ii) base model development, (iii) final model development, and (iv) model evaluation. The percentage of the issues in each category and the total issues were calculated.

RESULTS

Data from a total of 317 applications were analyzed. Of these, 122 applications had at least one PopPK-related issue, and the count went up to 168 issues. The proportion of issues for each category was ranked as follows: model evaluation (39%, 65 of 168), final model development (34%, 57 of 168), data used for PopPK analysis (14%, 24 of 168), and base model development (13%, 22 of 168). The most common issues were related to the goodness-of-fit plots, covariate selection, and high shrinkage of ETA, accounting for 17.9%, 15.5%, and 11.3% of the total issues, respectively.

CONCLUSION

The findings from this study may help the applicant understand the FDA's thinking on the PopPK analysis.

Studying Telmisartan Plasma Exposure, Kidney Distribution, Receptor Occupancy, and Response in Patients With Type 2 Diabetes Using [11 C]Telmisartan

The angiotensin receptor blocker telmisartan slows progression of kidney disease in patients with type 2 diabetes (T2D), yet many patients remain at high risk for progressive kidney function loss. The underlying mechanisms for this response variation might be attributed to differences in angiotensin-1 receptor occupancy (RO), resulting from individual variation in plasma drug exposure, tissue drug exposure, and receptor availability. Therefore, we first assessed the relationship between plasma telmisartan exposure and urinary-albumin-to-creatinine-ratio (UACR) in 10 patients with T2D and albuminuria (mean age 66 years, median UACR 297 mg/g) after 4 weeks treatment with 80 mg telmisartan once daily. Increasing telmisartan exposure associated with a larger reduction in UACR (Pearson correlation coefficient (PCC) = -0.64, P = 0.046, median change UACR: -40.1%, 95% confidence interval (CI): -22.9 to -77.4%, mean telmisartan area under the curve (AUC) = 2927.1 ng·hour/mL, 95% CI: 723.0 to 6501.6 ng·hour/mL). Subsequently, we assessed the relation among plasma telmisartan exposure, kidney distribution, and angiotensin-1 RO in five patients with T2D (mean age 60 years, median UACR 72 mg/g) in a separate positron emission tomography imaging study with [¹¹ C]Telmisartan. Individual plasma telmisartan exposure correlated with telmisartan distribution to the kidneys (PCC = 0.976, P = 0.024). A meaningful RO could be calculated in three patients receiving 120 mg oral telmisartan, and although high exposure seems related to higher RO, with AUC_0-last of 31, 840, and 274 ng·hour/mL and corresponding RO values 5.5%, 44%, and 59%, this was not significant (P = 0.64). Together these results indicate, for the first time, a relationship among interindividual differences in plasma exposure, kidney tissue distribution, RO, and ultimately UACR response after telmisartan administration.

Optimal exposure to docetaxel in adjuvant chemotherapy for early-stage breast cancer

WHAT IS KNOWN AND OBJECTIVE

Drug-induced neutropenia is the main reason for the dose limitation of docetaxel in patients with breast cancer. The area under the drug concentration-time curve (AUC) of docetaxel is associated with neutropenia. However, the optimal exposure to docetaxel for receiving postoperative adjuvant chemotherapy remains unclear. Therefore, we aimed to evaluate the relationship between the docetaxel AUC and neutropenia, identify potential influencing factors, and explore the best monitoring target for docetaxel when treating patients with early-stage breast cancer using a population pharmacokinetic (PopPK) model.

METHODS

Docetaxel plasma concentration, demographics, clinical data, and related laboratory data were collected. PopPK analyses were performed using a nonlinear mixed-effect modelling program. The docetaxel AUC was determined using the maximum a posteriori Bayesian (MAPB) method. The docetaxel exposure-toxicity threshold measured from the AUC for neutropenia was determined using the receiver operating characteristic (ROC) curve. The correlation between docetaxel exposure and neutropenia was analysed using multivariable logistic regression.

RESULTS

Among the 70 participants, 47 (67.1%) developed severe neutropenia. The PopPK analysis showed that the typical drug clearance (CL) rate was 37.4 L/h. Age was a significant covariate of CL rate, and aspartate aminotransferase and albumin levels were covariables of the volume of distribution. The multivariable regression analysis showed that AUC >3.0 mg.h/L (odds ratio [OR], 5.940; 95% confidence interval [CI], 1.693-20.843; P = 0.005), platinum use (OR, 0.156; 95% CI, 0.043-0.562; P = 0.005) and baseline haemoglobin level (OR, 0.938; 95% CI, 0.887-0.993; P = 0.027) were significant factors influencing the occurrence of grade 3/4 neutropenia. The AUC of first cycle may not predict the occurrence rates of grade 3/4 neutropenia in later cycles.

WHAT IS NEW AND CONCLUSION

We developed a docetaxel PopPK model for patients with early-stage breast cancer. Age and AST and ALB levels were significant covariates. AUC estimated using the MAPB method can predict the toxicity of docetaxel in patients with breast cancer. Docetaxel AUC >3.0 mg.h/L, absence of platinum use and low baseline haemoglobin level were risk factors for docetaxel-induced grade 3/4 neutropenia.

STUDY REGISTRATION

Chinese Clinical Trial Center Registry (ChiCTR2200056460).

Pharmacogenetics-based population pharmacokinetic analysis and dose optimization of valproic acid in Chinese southern children with epilepsy: Effect of ABCB1 gene polymorphism

Objective: The aim of this study was to establish a population pharmacokinetic (PPK) model of valproic acid (VPA) in pediatric patients with epilepsy in southern China, and provide guidance for individualized medication of VPA therapy. Methods: A total of 376 VPA steady-state trough concentrations were collected from 103 epileptic pediatric patients. The PPK parameter values for VPA were calculated by using the nonlinear mixed-effects modeling (NONMEM) method, and a one-compartment model with first-order absorption and elimination processes was applied. Covariates included demographic information, concomitant medications and selected gene polymorphisms. Goodness-of-fit (GOF), bootstrap analysis, and visual predictive check (VPC) were used for model evaluation. In addition, we used Monte Carlo simulations to propose dose recommendations for different subgroup patients. Results: A significant effect of the patient age and ABCB1 genotypes was observed on the VPA oral clearance (CL/F) in the final PPK model. Compared with patients with the ABCB1 rs3789243 AA genotype, CL/F in patients with GG and AG genotypes was increased by 8% and reduced by 4.7%, respectively. The GOF plots indicated the satisfactory predictive performance of the final model, and the evaluation by bootstrap and VPC showed that a stable model had been developed. A table of individualized dosing regimens involving age and ABCB1 genotype was constructed based on the final PPK model. Conclusion: This study quantitatively investigated the effects of patient age and ABCB1 rs3789243 variants on the pharmacokinetic variability of VPA. The PPK models could be beneficial to individual dose optimization in epileptic children on VPA therapy.

Population Pharmacokinetic Modeling and Probability of Target Attainment of Ceftaroline in Brain and Soft Tissues

Ceftaroline, approved to treat skin infections and pneumonia due to methicillin-resistant Staphylococcus aureus (MRSA), has been considered for the treatment of central nervous system (CNS) infections. A population pharmacokinetic (popPK) model was developed to describe ceftaroline soft tissue and cerebrospinal fluid (CSF) distributions and investigate the probability of target attainment (PTA) of the percentage of the dosing interval that the unbound drug concentration exceeded the MIC (%fT_>MIC) to treat MRSA infections. Healthy subjects' plasma and microdialysate concentrations from muscle and subcutaneous tissue following 600 mg every 12 h (q12h) and q8h and neurosurgical patients' plasma and CSF concentrations following single 600-mg dosing were used. Plasma concentrations were described by a two-compartment model, and tissue concentrations were incorporated as three independent compartments linked to the central compartment by bidirectional transport (clearance in [CL_in] and CL_out). Apparent volumes were fixed to physiological interstitial values. Healthy status and body weight were identified as covariates for the volume of the central compartment, and creatinine clearance was identified for clearance. The CSF glucose concentration (GLUC) was inversely correlated with CL_in,CSF. Simulations showed a PTA of >90% in plasma and soft tissues for both regimens assuming an MIC of 1 mg/L and a %fT_>MIC of 28.8%. Using the same target, patients with inflamed meninges (0.5 < GLUC ≤ 2 mmol/L) would reach PTAs of 99.8% and 97.2% for 600 mg q8h and q12h, respectively. For brain infection with mild inflammation (2 < GLUC ≤ 3.5 mmol/L), the PTAs would be reduced to 34.3% and 9.1%, respectively. Ceftaroline's penetration enhanced by meningeal inflammation suggests that the drug could be a candidate to treat MRSA CNS infections.

Pharmacokinetics and Adverse Effects of Clofazimine in the Treatment of Pulmonary Non-Tuberculous Mycobacterial Infection

Clofazimine (CFZ) is used to treat pulmonary non-tuberculous mycobacterial (NTM) infection; however, its pharmacokinetics remain unexplored in patients with pulmonary NTM, and the relationship between CFZ serum concentration and adverse effects has not been investigated. The objectives of this study were to characterize the pharmacokinetics of CFZ in pulmonary NTM disease treatment and to investigate the relationship between the steady-state CFZ serum concentration and adverse effects. A prospective observational study was conducted on 45 patients with pulmonary NTM treated with CFZ (UMIN000041053). A maximum of five serum samples per patient were taken at the CFZ trough, and serum concentration was measured using high-performance liquid chromatography-mass spectrometry (HPLC-MS). The pharmacokinetics of CFZ were analyzed using a nonlinear mixed effect model. The relationships among steady-state CFZ serum concentration and adverse effects, pigmentation, and heart rate-corrected QT (QTc) interval were investigated. Twenty-six patients had M. avium or M. intracellulare infection and nineteen had M. abscessus infection. The primary CFZ dosage was 50 mg/day. The estimated apparent CFZ clearance, apparent volume of distribution, and half-life were 2.4 L/h, 2,960 L, and 36 days, respectively. The combined use of rifampicin and CFZ significantly reduced CFZ exposure by 22%. Although there was no relationship between CFZ serum concentration and pigmentation intensity, the QTc interval was significantly correlated with CFZ serum concentration. The estimation of accurate pharmacokinetics for CFZ required approximately 5 months of monitoring. The relationship between the serum concentration and specific adverse effects of CFZ confirmed that CFZ serum concentration was not associated with pigmentation but did affect the QTc interval.

Population pharmacokinetic of paracetamol and atorvastatin with co-administration of semaglutide

Semaglutide is a glucagon-like-peptide-1 (GLP-1) analogue marketed for once-weekly subcutaneous administration for type 2 diabetes mellitus. Like other long-acting GLP-1 analogues, semaglutide reduces gastric emptying and, potentially, alters the rate of absorption of orally co-administered drugs. The objective of the current analysis was to evaluate the effects on the gastric emptying rate caused by semaglutide on pharmacokinetic model parameters of paracetamol and atorvastatin in healthy subjects. Non-linear mixed effect modeling was used to estimate population pharmacokinetic model parameters of paracetamol and atorvastatin after single doses with or without semaglutide. The absorption rate (ka) of paracetamol decreased by 53% when co-administered with semaglutide. For atorvastatin, ka and transit compartment rate (ktr) decreased by 72% and 91%, respectively. Thus, gastric emptying, measured as T50, i.e., drug disappearance from the absorption compartments, showed an additional 5-min delay for paracetamol and a 67-min delay for atorvastatin when co-administered with semaglutide. Semaglutide affected pharmacokinetic model parameters of paracetamol and atorvastatin, and minor quantitative differences in gastric emptying between placebo vs. semaglutide administration were observed. However, these effects of semaglutide were considered not to be of clinical relevance.

Population pharmacokinetic modelling of imatinib in healthy subjects receiving a single dose of 400 mg

Purpose

Imatinib is indicated for treatment of CML, GIST, etc. The population pharmacokinetics (popPK) of imatinib in patients under long-term treatment are reported in literature. Data obtained from bioequivalence trials for healthy subjects were used to evaluate the influence of demographic and pharmacogenetic factors on imatinib pharmacokinetics (PK) in a collective without concurrent drugs, organ dysfunction, inflammation etc. In addition, the differences in PK between the healthy subjects and a patient cohort was examined to identify possible disease effects.

Methods

26 volunteers were administered orally with single dose of 400 mg imatinib. 16–19 plasma samples per volunteer were collected from 0.5 up to 72 h post-dose. The popPK was built and post hoc estimates were compared with previously published PK parameters evaluated by non-compartmental analysis in the same cohort. The predictivity of the model for data collected from 40 patients with gastrointestinal stromal tumors at steady state was evaluated.

Results

The popPK was best described by a two-compartment transit model with first-order elimination. No significant covariates were identified, probably due to the small cohort and the narrow range of demographic covariates; CYP3A5 phenotypes appeared to have some influence on the clearance of imatinib. Good agreement between non-compartment and popPK analyses was observed with the differences of the geometric means/ median of PK estimates below 10%. The model indicated lower clearance for patients compared to healthy volunteers (p value < 0.01).

Conclusion

The two-compartment transit model adequately describes the absorption and distribution of imatinib in healthy volunteers. For patients, a lower clearance of imatinib compared to healthy volunteer was estimated by the model. The model can be applied for dose individualization based on trough concentrations assuming no significant differences in absorption between patients and healthy volunteers.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00280-022-04454-y.

Population Pharmacokinetics and Dosing Regimen of Lithium in Chinese Patients With Bipolar Disorder

Background: Lithium is an effective medication approved for the treatment of bipolar disorder (BD). It has a narrow therapeutic index (TI) and requires therapeutic drug monitoring. This study aimed to conduct a population pharmacokinetics (PPK) analysis of lithium and investigate the appropriateness of the dosing regimen according to different patient characteristics.

Methods: A total of 476 lithium concentrations from 268 patients with bipolar disorder were analyzed using nonlinear mixed-effects modeling. Monte Carlo simulations were employed to investigate the influence of covariates, such as weight, creatinine clearance, and daily doses of lithium concentrations, and to determine the individualized dosing regimens for patients.

Results: Lithium PK was described by a one-compartment model with first-order absorption and elimination processes. The typical estimated apparent clearance was 0.909 L/h⁻¹ with 16.4% between-subject variability in the 62 kg patients with 116 ml/min creatinine clearance and 600 mg daily doses. To achieve a target trough concentration (0.4–0.8 mmol/L) in the maintenance phase, the regimen of 500 mg than 750 mg daily dose was recommended for patients with renal insufficiency and weighing 100 kg.

Conclusion: A PPK model for lithium was developed to determine the influence of patient characteristics on lithium pharmacokinetics. Weight, creatinine clearance, and total daily dose of lithium can affect the drug’s clearance. These results demonstrate the nonlinear renal excretion of lithium; hence, dosage adjustments are recommended for patients with renal insufficiency.

Introduction of an artificial neural network–based method for concentration‐time predictions

Pharmacometrics and the application of population pharmacokinetic (PK) modeling play a crucial role in clinical pharmacology. These methods, which describe data with well‐defined equations and estimate physiologically interpretable parameters, have not changed substantially during the past decades. Although the methods have proven their usefulness, they are often resource intensive and require a high level of expertise. We investigated whether a method based on artificial neural networks (ANNs) may provide an alternative approach for the prediction of concentration‐time curve to supplement the gold standard methods. In this work, we used simulated data to overcome the requirement for a large clinical training data set, implemented a pharmacologically reasonable network architecture to improve extrapolation to different dosing schemes, and used transfer learning to quickly adapt the predictions to new patient groups. We demonstrate that ANNs are able to learn the shape of concentration‐time curves and make individual predictions based on a short sequence of PK measurements. Furthermore, an ANN trained on simulated data was applied to real clinical data and was demonstrated to extrapolate to different dosing schemes. We also adapted the ANN trained on simulated healthy subjects to simulated hepatic impaired patients through transfer learning. In summary, we demonstrate how ANNs could be leveraged in a PK workflow to efficiently make individual concentration‐time predictions, and we discuss the current limitations and advantages of such an ANN‐based method.

Population Pharmacokinetic Modeling and Simulation of TV-46000: A Long-Acting Injectable Formulation of Risperidone

TV‐46000 is a long‐acting subcutaneous antipsychotic that uses a novel copolymer drug delivery technology in combination with a well‐characterized molecule, risperidone, that is in clinical development as a treatment for schizophrenia. A population pharmacokinetic (PPK) modeling and simulation approach was implemented to identify TV‐46000 doses and dosing schedules for clinical development that would provide the best balance between clinical efficacy and safety. The PPK model was created by applying pharmacokinetic data from a phase 1 study of 97 patients with a diagnosis of schizophrenia or schizoaffective disorder who received either single or repeated doses of TV‐46000. The PPK model was used to characterize the complex release profile of the total active moiety (TAM; the sum of the risperidone and 9‐OH risperidone concentrations) concentration following subcutaneous injections of TV‐46000. The PK profile was best described by a double Weibull function of the in vivo release rate and by a 2‐compartment disposition and elimination model. Simulations were performed to determine TV‐46000 doses and dosing schedules that maintained a median profile of TAM concentrations similar to published TAM exposure following oral risperidone doses that have been correlated to a 40% to 80% dopamine‐D2 receptor occupancy therapeutic window. The simulations showed that therapeutic dose ranges for TV‐46000 are 50 to 125 mg for once‐monthly and 100 to 250 mg for the once every 2 months regimens. This PPK model provided a basis for prediction of patient‐specific exposure and dopamine‐D2 receptor occupancy estimates to support further clinical development and dose selection for the phase 3 studies.

R-praziquantel integrated population pharmacokinetics in preschool- and school-aged African children infected with Schistosoma mansoni and S. haematobium and Lao adults infected with Opisthorchis viverrini

Racemic praziquantel (PZQ) is the standard treatment for schistosomiasis and liver fluke infections (opisthorchiasis and clonorchiasis). The development of an optimal pediatric formulation and dose selection would benefit from a population pharmacokinetic (popPK) model. A popPK model was developed for R-PZQ, the active enantiomer of PZQ, in 664 subjects, 493 African children (2-15 years) infected with Schistosoma mansoni and S. haematobium, and 171 Lao adults (15-78 years) infected with Opisthorchis viverrini. Racemate tablets were administered as single doses of 20, 40 and 60 mg/kg in children and 30, 40 and 50 mg/kg in 129 adults, and as 3 × 25 mg/kg apart in 42 adults. Samples collected by the dried-blood-spot technique were assayed by LC-MS/MS. A two-compartment disposition model, with allometric scaling and dual first-order and transit absorption, was developed using Phoenix™ software. Inversely parallel functions of age described the apparent oral bioavailability (BA) and clearance maturation in children and ageing in adults. BA decreased slightly in children with dose increase, and by 35% in adults with multiple dosing. Crushing tablets for preschool-aged children increased the first-order absorption rate by 64%. The mean transit absorption time was 70% higher in children. A popPK model for R-PZQ integrated African children over 2 years of age with schistosomiasis and Lao adults with opisthorchiasis, and should be useful to support dose optimization in children. In vitro hepatic and intestinal metabolism data would help refining and validating the model in younger children as well as in target ethnic pediatric and adult groups.

Eculizumab precision dosing algorithm for thrombotic microangiopathy in children and young adults undergoing HSCT

This large study examined eculizumab PK/PD in bleeding and nonbleeding patients with TA-TMA.

PK/PD model-based eculizumab dosing is needed for bleeding in TA-TMA, whereas fixed-dose regimens are effective in nonbleeding patients.

Transplant-associated thrombotic microangiopathy (TA-TMA) is a fatal posttransplant complication of hematopoietic stem cell transplantation. We recently reported that survival for TA-TMA has been improved by early intervention with eculizumab, a complement C5 inhibitor, guided by pharmacokinetic/pharmacodynamic (PK/PD) model-informed precision dosing. However, patients with gastrointestinal bleeding showed poor survival, even when treated with more frequent doses. Our objective was to develop separate models in bleeding and nonbleeding patients with TA-TMA and to propose precision dosing algorithms. Eculizumab PK/PD was analyzed in 19 bleeding and 38 nonbleeding patients (0.5-29.9 years of age). A complement activation biomarker (sC5b-9) and body weight were identified as significant determinants of eculizumab clearance regardless of bleeding. Eculizumab clearance after the first dose was higher in bleeding than in nonbleeding patients (83.8 vs 61.3 mL/h per 70 kg; P = .07). The high clearance was maintained over treatment doses in bleeding patients, whereas nonbleeding patients showed a time-dependent decrease in clearance. sC5b-9 levels were highest before the first dose and decreased over time, regardless of bleeding complications. A Monte Carlo Simulation analysis showed that the current dosing protocols recommended for atypical hemolytic uremic syndrome had <15% probability of attaining the target concentration of >100 μg/mL eculizumab in nonbleeding patients. We identified an intensified loading protocol to reach 80% target attainment. Our data clearly showed the need for individualized dosing for patients with significant bleeding and for ongoing dose adjustments to optimize outcomes. The developed models will be incorporated into a clinical decision guideline for precision dosing to improve outcomes in children and young adults with TA-TMA.