Mechanism-based concepts of size and maturity in pharmacokinetics

Population pharmacokinetic analyses of methotrexate in pediatric patients: a systematic review

BACKGROUND AND OBJECTIVES

Methotrexate is widely utilized in the chemotherapy of malignant tumors and autoimmune diseases in the pediatric population, but dosing can be challenging. Several population pharmacokinetic models were developed to characterize factors influencing variability and improve individualization of dosing regimens. However, significant covariates included varied across studies. The primary objective of this review was to summarize and discuss population pharmacokinetic models of methotrexate and covariates that influence pharmacokinetic variability in pediatric patients.

METHODS

Systematic searches were conducted in the PubMed and EMBASE databases from inception to 7 July 2023. Reporting Quality was evaluated based on a checklist with 31 items. The characteristics of studies and information for model construction and validation were extracted, summarized, and discussed.

RESULTS

Eighteen studies (four prospective studies and fourteen retrospective studies with sample sizes of 14 to 772 patients and 2.7 to 93.1 samples per patient) were included in this study. Two-compartment models were the commonly used structural models for methotrexate, and the clearance range of methotrexate ranged from 2.32 to 19.03 L/h (median: 6.86 L/h). Body size and renal function were found to significantly affect the clearance of methotrexate for pediatric patients. There were limited reports on the role of other covariates, such as gene polymorphisms and co-medications, in the pharmacokinetic parameters of methotrexate pediatric patients. Internal and external evaluations were used to assess the performance of the population pharmacokinetic models.

CONCLUSION

A more rigorous external evaluation needs to be performed before routine clinical use to select the appropriate PopPK model. Further research is necessary to incorporate larger cohorts or pool analyses in specific susceptible pediatric populations to improve the understanding of predicted exposure profiles and covariate identification.

Absorption pharmacokinetics and feasibility of intranasal dexmedetomidine in patients under general anaesthesia

BACKGROUND

The use of intranasal dexmedetomidine is hampered by a limited understanding of its absorption pharmacokinetics.

METHODS

We examined the pharmacokinetics and feasibility of intranasal dexmedetomidine administered in the supine position to adult patients undergoing general anaesthesia. Twenty-eight patients between 35 and 80 years of age, ASA 1-3 and weight between 50 and 100 kg, who underwent elective unilateral total hip or knee arthroplasty under general anaesthesia were recruited. All patients received 100 μg of intranasal dexmedetomidine after anaesthesia induction. Six venous blood samples (at 0, 5, 15, 45, 60, 240 min timepoints from dexmedetomidine administration) were collected from each patient and dexmedetomidine plasma concentrations were measured. Concentration-time profiles after nasal administration were pooled with earlier data from a population analysis of intravenous dexmedetomidine (n = 202) in order to estimate absorption parameters using nonlinear mixed effects. Peak concentration (CMAX) and time (TMAX) were estimated using simulation (n = 1000) with parameter estimates and their associated variability.

RESULTS

There were 28 adult patients with a mean (SD) age of 66 (8) years and weight of 83 (10) kg. The mean weight-adjusted dose of dexmedetomidine was 1.22 (0.15) μg kg-1. CMAX 0.273 μg L-1 was achieved at 98 min after intranasal administration (TMAX). The relative bioavailability of dexmedetomidine was 80% (95% CI 75-91%). The absorption half-time (TABS = 120 min; 95% CI 90-147 min) was slower than that in previous pharmacokinetic studies on adult patients. Perioperative haemodynamics of all patients remained stable.

CONCLUSIONS

Administration of intranasal dexmedetomidine in the supine position during general anaesthesia is feasible with good bioavailability. This administration method has slower absorption when compared to awake patients in upright position, with consequent concentrations attained after TMAX for several hours.

Model-Informed Development of a Cost-Saving Dosing Regimen for Sacituzumab Govitecan

BACKGROUND

The antibody-drug conjugate sacituzumab govitecan is approved for metastatic triple-negative breast cancer and has shown promising results in various other types of cancer. Its costs may limit patient access to this novel effective treatment modality.

OBJECTIVE

The purpose of this study was to develop an evidence-based rational dosing regimen that results in targeted drug exposure within the therapeutic range while minimizing financial toxicity, to improve treatment access.

PATIENTS AND METHODS

Exposure equivalent dosing strategies were developed based on pharmacokinetic modeling and simulation by using the published pharmacokinetic model developed by the license holder. The alternative dose was based on the principle of using complete vials to prevent spillage and on the established non-linear relationship between body weight and systemic exposure. Equivalent exposure compared to the approved dosing regimen of 10 mg/kg was aimed for. Equivalent exposure was conservatively defined as calculated geometric mean ratios within the 0.9-1.11 boundaries for area under the concentration-time curve (AUC), trough concentration (Ctrough) and maximum concentration (Cmax) of the alternative dosing regimen compared to the approved dosing regimen. Since different vial sizes are available for the European Union (EU) and United States (US) market, because body weight distributions differ between these populations, we performed our analysis for both scenarios.

RESULTS

Dosing regimens of sacituzumab govitecan for the EU (< 50 kg: 400 mg, 50-80 kg: 600 mg, and > 80 kg: 800 mg) and US population (< 40 kg: 360 mg, 40-65 kg: 540 mg, 65-90 kg: 720 mg, and > 90 kg: 900 mg) were developed, based on weight bands. The geometric mean ratios for all pharmacokinetic outcomes were within the predefined equivalence boundaries, while the quantity of drug used was 21.5% and 19.0% lower for the EU and US scenarios, respectively.

CONCLUSIONS

With the alternative dosing proposal, an approximately 20% reduction in drug expenses for sacituzumab govitecan can be realized while maintaining an equivalent and more evenly distributed exposure throughout the body weight range, without notable increases in pharmacokinetic variability.

Population pharmacokinetics and CD20 binding dynamics for mosunetuzumab in relapsed/refractory B-cell non-Hodgkin lymphoma

Mosunetuzumab (Mosun) is a CD20xCD3 T-cell engaging bispecific antibody that redirects T cells to eliminate malignant B cells. The approved step-up dose regimen of 1/2/60/30 mg IV is designed to mitigate cytokine release syndrome (CRS) and maximize efficacy in early cycles. A population pharmacokinetic (popPK) model was developed from 439 patients with relapsed/refractory B-Cell Non-Hodgkin lymphoma receiving Mosun IV monotherapy, including fixed dosing (0.05-2.8 mg IV every 3 weeks (q3w)) and Cycle 1 step-up dosing groups (0.4/1/2.8-1/2/60/30 mg IV q3w). Prior to Mosun treatment, ~50% of patients had residual levels of anti-CD20 drugs (e.g., rituximab or obinutuzumab) from prior treatment. CD20 receptor binding dynamics and rituximab/obinutuzumab PK were incorporated into the model to calculate the Mosun CD20 receptor occupancy percentage (RO%) over time. A two-compartment model with time-dependent clearance (CL) best described the data. The typical patient had an initial CL of 1.08 L/day, transitioning to a steady-state CL of 0.584 L/day. Statistically relevant covariates on PK parameters included body weight, albumin, sex, tumor burden, and baseline anti-CD20 drug concentration; no covariate was found to have a clinically relevant impact on exposure at the approved dose. Mosun CD20 RO% was highly variable, attributed to the large variability in residual baseline anti-CD20 drug concentration (median = 10 μg/mL). The 60 mg loading doses increased Mosun CD20 RO% in Cycle 1, providing efficacious exposures in the presence of the competing anti-CD20 drugs. PopPK model simulations, investigating Mosun dose delays, informed treatment resumption protocols to ensure CRS mitigation.

Population pharmacokinetics of ravidasvir in adults with chronic hepatitis C virus infection and impact of antiretroviral treatment

Ravidasvir (RDV) is a novel NS5A inhibitor that exhibits potent pan-genotypic inhibition of hepatitis C virus (HCV) replication. Sofosbuvir (SOF) plus RDV was demonstrated to be efficacious and safe in adults with active HCV infection, including those living with HIV (LWHIV), in the STORM-C-1 trial. We assessed the population pharmacokinetics (PK) of RDV in a sub-study nested within STORM-C-1 conducted in Thailand and Malaysia. SOF (400 mg) plus RDV (200 mg) was administered orally once daily for 12 weeks to adults with chronic HCV infection, but without cirrhosis and for 24 weeks to those with compensated cirrhosis. Intensive and sparse PK samples were collected at 4, 8, and 12 weeks after treatment initiation. Population PK parameters of RDV and the impact of covariates were evaluated using nonlinear mixed-effects modeling. Five hundred ninety-four participants were included, 235 (40%) had compensated cirrhosis, and 189 (32%) were LWHIV. RDV plasma concentrations were best described by a two-compartment model with first-order elimination. Oral clearance (CL/F) and volume of distribution (Vd/F) parameters were allometrically scaled on fat-free mass. Concomitant antiretroviral treatment (ART) increased RDV CL/F by 30%-60%, with efavirenz-based ART having the largest impact. Females had 16% lower RDV CL/F than males, and higher albumin levels reduced RDV central volume of distribution. While several covariates impact RDV CL/F and Vd/F, the effect on RDV exposures was not clinically relevant based on the efficacy data reported in this diverse Asian adult population. There were no meaningful drug-drug interactions in adults LWHIV on ART.

Population pharmacokinetic modelling of primaquine exposures in lactating women and breastfed infants

Current guidelines advise against primaquine treatment for breastfeeding mothers to avoid the potential for haemolysis in infants with G6PD deficiency. To predict the haemolytic risk, the amount of drug received from the breast milk and the resulting infant drug exposure need to be characterised. Here, we develop a pharmacokinetic model to describe the drug concentrations in breastfeeding women using venous, capillary, and breast milk data. A mother-to-infant model is developed to mimic the infant feeding pattern and used to predict their drug exposures. Primaquine and carboxyprimaquine exposures in infants are <1% of the exposure in mothers. Therefore, even in infants with the most severe G6PD deficiency variants, it is highly unlikely that standard doses of primaquine (0.25-1 mg base/kg once daily given to the mother for 1-14 days) would cause significant haemolysis. After the neonatal period, primaquine should not be restricted for breastfeeding women (Clinical Trials Registration: NCT01780753).

Population Pharmacokinetics of Nirsevimab in Preterm and Term Infants

Nirsevimab, a monoclonal antibody with an extended half-life, is approved for the prevention of respiratory syncytial virus (RSV) disease in all infants in Canada, the EU, Great Britain, and the USA. A population pharmacokinetics (PK) model was built to describe the PK of nirsevimab in preterm and term infants, and to evaluate the influence of covariates, including body weight and age, in infants. Nirsevimab PK was characterized by a 2-compartment model with first-order clearance (CL) and first-order absorption following intramuscular (IM) administration. The typical CL in a 5 kg infant was 3.4 mL/day. Body weight and postmenstrual age were the primary covariates on CL, with minor effects for race, second RSV season, and antidrug antibody status (deemed not clinically relevant). Congenital heart disease (CHD) and chronic lung disease (CLD) did not significantly impact nirsevimab PK. The final population PK model, based on 8987 PK observations from 2683 participants across 5 clinical trials, successfully predicted PK in an additional cohort of 967 healthy infants. Weight-banded dosing (50 mg in infants <5 kg; 100 mg in infants ≥5 kg) was predicted to be appropriate for infants ≥1 kg in their first RSV season. Together, these data support weight-banded dosing of nirsevimab in all infants in their first RSV season, including in healthy infants, infants with CHD or CLD, and in infants born prematurely.

Population Pharmacokinetics of Pediatric Lopinavir/Ritonavir Oral Pellets in Children Living with HIV in Africa

Antiretroviral therapy for children living with HIV (CLHIV) under 3 years of age commonly includes lopinavir/ritonavir (LPV/r). However, the original liquid LPV/r formulation has taste and cold storage difficulties. To address these challenges, LPV/r oral pellets have been developed. These pellets can be mixed with milk or food for administration and do not require refrigeration. We developed the population pharmacokinetic (PK) model and assessed drug exposure of LPV/r oral pellets administered twice daily to CLHIV per World Health Organization (WHO) weight bands. The PK analysis included Kenyan and Ugandan children participating in the LIVING studies (NCT02346487) receiving LPV/r pellets (40/10 mg) and ABC/3TC (60/30 mg) dispersible tablets. Population PK models were developed for lopinavir (LPV) and ritonavir (RTV) to evaluate the impact of RTV on the oral clearance (CL/F) of LPV. The data obtained from the study were analyzed using nonlinear mixed-effects modeling approach. Data from 514 children, comprising a total of 2,998 plasma concentrations of LPV/r were included in the analysis. The LPV and RTV concentrations were accurately represented by a one-compartment model with first-order absorption (incorporating a lag-time) and elimination. Body weight influenced LPV and RTV PK parameters. The impact of RTV concentrations on the CL/F of LPV was characterized using a maximum effect model. Simulation-predicted target LPV exposures were achieved in children with this pellet formulation across the WHO weight bands. The LPV/r pellets dosed in accordance with WHO weight bands provide adequate LPV exposures in Kenyan and Ugandan children weighing 3.0 to 24.9 kg.

Covariate modeling in pharmacometrics: General points for consideration

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

Time-Varying Clearance in Milrinone Pharmacokinetics from Premature Neonates to Adolescents

BACKGROUND AND OBJECTIVES

Milrinone is an inotrope and vasodilator used for prophylaxis or treatment of low cardiac output syndrome after weaning from cardiopulmonary bypass (CPB). It is renally eliminated and has an acceptable therapeutic range of 100-300 μg/L, but weight-based dosing alone is associated with poor target attainment. We aimed to develop a population pharmacokinetic model for milrinone from premature neonates to adolescents, and to evaluate how age, renal function and recovery from CPB may impact dose selection.

METHODS

Fifty paediatric patients (aged 4 days to 16 years) were studied after undergoing cardiac surgery supported by CPB. Data from 29 premature neonates (23-28 weeks' postmenstrual age) treated for prophylaxis of low systemic blood flow were available for a pooled pharmacokinetic analysis. Population parameters were estimated using non-linear mixed effects modelling (NONMEM 7.5.1).

RESULTS

There were 369 milrinone measurements available for analysis. A one-compartment model with zero-order input and first-order elimination was used to describe milrinone disposition. Population parameters were clearance 17.8 L/70 kg [95% CI 15.8-19.9] and volume 20.4 L/h/70 kg [95% CI 17.8-22.1]. Covariates included size, postmenstrual age and renal function for clearance, and size and postnatal age for volume. Milrinone clearance is reduced by 39.5% [95% CI 24.0-53.7] immediately after bypass, and recovers to baseline clearance with a half-time of 12.0 h [95% CI 9.7-15.2]. Milrinone volume was 2.07 [95% CI 1.87-2.27] times greater at birth than the population standard and decreased over the first days of life with a half-time of 0.977 days [95% CI 0.833-1.12].

CONCLUSION

Milrinone is predominately renally eliminated and so renal function is an important covariate describing variability in clearance. Increasing clearance over time likely reflects increasing cardiac output and renal perfusion due to milrinone and return to baseline following CPB.

Eculizumab dose tapering should take into account the nonlinearity of its pharmacokinetics

AIMS

Eculizumab is a monoclonal antibody targeting complement protein C5 used in renal diseases. As recommended dosing regimen leads to unnecessarily high concentrations in some patients, tailored dosing therapeutic drug monitoring was proposed to reduce treatment cost. The objectives of the present work were (i) to investigate the target-mediated elimination of eculizumab and (ii) whether a pharmacokinetic model integrating a nonlinear elimination allows a better prediction of eculizumab concentrations than a linear model.

METHODS

We analysed 377 eculizumab serum concentrations from 44 patients treated for atypical haemolytic uraemic syndrome and C3 glomerulopathy with a population pharmacokinetic approach. Critical concentrations (below which a non-log-linear decline of concentration over time is evidenced) were computed to estimate the relevance of the target-mediated elimination. Simulations of dosing regimens were then performed to predict probabilities of target attainment (i.e. trough >100 mg/L).

RESULTS

Pharmacokinetics of eculizumab was nonlinear and followed a mixture of first-order (CL = 1.318 mL/day/kg) and Michaelis-Menten elimination (Vmax = 26.07 mg/day, Km = 24.06 mg/L). Volume of distribution (72.39 mL/kg) and clearance were weight-dependent. Critical concentrations (Vmax/CL) ranged from 144.7 to 759.7 mg/L and were inversely related to body weight (P = .013). Nonlinearity was thus noticeable at therapeutic concentrations. Simulations predicted that 1200 mg of eculizumab every 21 days would allow 85% and 76% of patients to maintain a therapeutic exposure, for 50 or 90 kg body weight, respectively.

CONCLUSIONS

Our study investigates the nonlinear elimination of eculizumab and discusses the importance of accounting for eculizumab target-mediated elimination in therapeutic drug monitoring.

Population pharmacokinetics of cyclosporine A in pediatric patients with thalassemia undergoing allogeneic hematopoietic stem cell transplantation

PURPOSE

To establish the population pharmacokinetics (PPK) model of cyclosporine A(CsA) in pediatric patients with thalassemia undergoing allogeneic hematopoietic stem cell transplantation (HSCT), aiming at providing a reference for clinical dose individualization of CsA.

METHODS

Children with thalassemia who underwent allogeneic HSCT were enrolled retrospectively. The PPK structural model and the random variable model of CsA were established on NONMEN. And goodness of fit plots (GOFs), visual predictive check (VPC), and bootstrap and normalized prediction distribution errors (NPDE) were used to evaluate the final model.

RESULTS

A one-compartment model with first-order absorption was employed to fit the base model. A total of 74 pediatric patients and 600 observations of whole blood concentration were included. The final model included weight (WT) in clearance (CL), alongside post-operative day (POD), fluconazole (FLUC), voriconazole (VORI), posaconazole (POSA), and red blood cell count (RBC) significantly. All the model evaluations were passed.

CONCLUSION

In the PPK model based on the pediatric cohort on CsA with thalassemia undergoing allogeneic HSCT, WT, POD, FLUC, VORI, POSA, and RBC were found to be the significant factors influencing CL of CsA. The reliability and robustness of the final model were excellent. It is expected that the PPK model can assist in individualizing dosing strategy clinically.

Assessing Pharmacogenomic loci Associated with the Pharmacokinetics of Vamorolone in Boys with Duchenne Muscular Dystrophy

Human genetic variation (polymorphisms) in genes coding proteins involved in the absorption, distribution, metabolism, and elimination (ADME) of drugs can have a strong effect on drug exposure and downstream efficacy and safety outcomes. Vamorolone, a dissociative steroidal anti-inflammatory drug for treating Duchenne muscular dystrophy (DMD), primarily undergoes oxidation by CYP3A4 and CYP3A5 and glucuronidation by UDP-glucuronosyltransferases. This work assesses the pharmacokinetics (PKs) of vamorolone and sources of interindividual variability (IIV) in 81 steroid-naïve boys with DMD aged 4 to <7 years old considering the genetic polymorphisms of CYPS3A4 (CYP3A4*22, CYP3A4*1B), CYP3A5 (CYP3A5*3), and UGT1A1 (UGT1A1*60) utilizing population PK modeling. A one-compartment model with zero-order absorption (Tk0, duration of absorption), linear clearance (CL/F), and volume (V/F) describes the plasma PK data for boys with DMD receiving a wide range of vamorolone doses (0.25-6 mg/kg/day). The typical CL/F and V/F values of vamorolone were 35.8 L/h and 119 L, with modest IIV. The population Tk0 was 3.14 h yielding an average zero-order absorption rate (k0) of 1.16 mg/kg/h with similar absorption kinetics across subjects at the same vamorolone dose (i.e., no IIV on Tk0). The covariate analysis showed that none of the genetic covariates had any significant impact on the PKs of vamorolone in boys with DMD. Thus, the PKs of vamorolone is very consistent in these young boys with DMD.

Pediatric anesthesia in Australia and New Zealand and health inequity among First Nations and Māori children

Australia and New Zealand are two countries in the Southern Pacific region. They share many pediatric anesthesia similarities in terms of medical organizational systems, education, training, and research, however there are important differences between the two nations in relation to geography, the First Nations populations and the history of colonization. While the standards for pediatric anesthesia and the specialty training requirements are set by the Australian and New Zealand College of Anesthetists and the Society for Pediatric Anesthesia in New Zealand and Australia, colonization has created distinct challenges that each nation now faces in order to improve the anesthetic care of its pediatric population. Australia generally has a high standard of living and good access to health care; disparities exist for First Nations People and for those living in rural or remote areas. Two influences have shaped training within New Zealand over the past 40 years; establishment of a national children's hospital in 1990 and, more importantly, acknowledgement that the First Nations people of New Zealand (Māori) have suffered because of failure to recognize their rights consequent to establishing a partnership treaty between Māori and the British Crown in 1840. Health inequities among Māori in New Zealand and First Nations People in Australia have implications for the health system, culturally appropriate approaches to treatment, and the importance of having an appreciation of First Nations people's history and culture, language, family structure, and cultural safety. Trainees in both countries need to be adequately supported in these areas in order for the sub-specialty of pediatric anesthesia to develop further and improve the anesthetic and surgical outcomes of our children.

Population pharmacokinetic rationale for intravenous contezolid acefosamil followed by oral contezolid dosage regimens

Contezolid is a novel oxazolidinone antibiotic with a promising safety profile. Oral contezolid and its intravenous (IV) prodrug contezolid acefosamil (CZA) are in development for treatment of diabetic foot and acute bacterial skin and skin structure infections (ABSSSI). The prodrug CZA is converted to active contezolid via intermediate MRX-1352. This study aimed to provide the pharmacokinetic rationale for safe, effective, and flexible dosage regimens with initial IV CZA followed by oral contezolid. We simultaneously modeled plasma concentrations from 110 healthy volunteers and 74 phase 2 patients with ABSSSI via population pharmacokinetics (using the importance sampling estimation algorithm), and optimized dosage regimens by Monte Carlo simulations. This included data on MRX-1352, contezolid, and its metabolite MRX-1320 from 66 healthy volunteers receiving intravenous CZA (150-2400 mg) for up to 28 days, and 74 patients receiving oral contezolid [800 mg every 12 h (q12h)] for 10 days. The apparent total clearance for 800 mg oral contezolid with food was 16.0 L/h (23.4% coefficient of variation) in healthy volunteers and 17.7 L/h (53.8%) in patients. CZA was rapidly converted to MRX-1352, which subsequently transformed to contezolid. The proposed dosage regimen used an IV CZA 2000 mg loading dose with 1000 mg IV CZA q12h as maintenance dose(s), followed by 800 mg oral contezolid q12h (with food). During each 24-h period, Monte Carlo simulations predicted this regimen to achieve consistent areas under the curve of 91.9 mg·h/L (range: 76.3-106 mg·h/L) under all scenarios. Thus, this regimen was predicted to reliably achieve efficacious contezolid exposures independent of timing of switch from IV CZA to oral contezolid.IMPORTANCEThis study provides the population pharmacokinetic rationale for the dosage regimen of the intravenous (IV) prodrug contezolid acefosamil (CZA) followed by oral contezolid. We developed the first integrated population model for the pharmacokinetics of the MRX-1352 intermediate prodrug, active contezolid, and its main metabolite MRX-1320 based on data from three clinical studies in healthy volunteers and phase 2 patients. The proposed regimen was predicted to reliably achieve efficacious contezolid exposures independent of timing of switch from IV CZA to oral contezolid.

Learning pharmacometric covariate model structures with symbolic regression networks

Efficiently finding covariate model structures that minimize the need for random effects to describe pharmacological data is challenging. The standard approach focuses on identification of relevant covariates, and present methodology lacks tools for automatic identification of covariate model structures. Although neural networks could potentially be used to approximate covariate-parameter relationships, such approximations are not human-readable and come at the risk of poor generalizability due to high model complexity.In the present study, a novel methodology for the simultaneous selection of covariate model structure and optimization of its parameters is proposed. It is based on symbolic regression, posed as an optimization problem with a smooth loss function. This enables training of the model through back-propagation using efficient gradient computations.Feasibility and effectiveness are demonstrated by application to a clinical pharmacokinetic data set for propofol, containing infusion and blood sample time series from 1031 individuals. The resulting model is compared to a published state-of-the-art model for the same data set. Our methodology finds a covariate model structure and corresponding parameter values with a slightly better fit, while relying on notably fewer covariates than the state-of-the-art model. Unlike contemporary practice, finding the covariate model structure is achieved without an iterative procedure involving manual interactions.

Quantifying the Effect of Methotrexate on Adalimumab Response in Psoriasis by Pharmacokinetic-Pharmacodynamic Modeling

Previously, we showed that the combination of methotrexate and adalimumab treatment leads to less antidrug antibody development. In this study, we quantify the pharmacokinetics/pharmacodynamics (PK/PD) of adalimumab and evaluate the influence of methotrexate cotreatment. A population PK-PD model was developed using prospective data from 59 patients with psoriasis (baseline PASI = 12.6) receiving adalimumab over 49 weeks. Typical PK and PD parameters and their corresponding interpatient variability were estimated. We performed a covariate analysis to assess whether interpatient variability could be explained by addition of methotrexate and other covariates. In total, 330 PASIs, 252 adalimumab serum concentrations, and 247 antidrug antibody titers were available. Presence of antidrug antibodies (adalimumab group = 46.7%, adalimumab + methotrexate group = 38.7%; P = .031) was correlated with increased adalimumab apparent clearance (P < .001). In the PD model, the use of concomitant methotrexate was borderline to significantly correlated with a decreased half-maximal inhibitory concentration (adalimumab concentration for which clinical response score is reduced by half; P < .10). On the basis of our PK-PD model, concomitant use of methotrexate indirectly increases adalimumab concentration, partially through less antidrug antibodies formation, which may result in better efficacy.

Pharmacokinetic Profile of Brepocitinib with Topical Administration in Atopic Dermatitis and Psoriasis Populations: Strategy to Inform Clinical Trial Design in Adult and Pediatric Populations

INTRODUCTION

Topical brepocitinib, a tyrosine kinase (TYK)2/Janus kinase (JAK)1 inhibitor, is in development for psoriasis (PsO) and atopic dermatitis (AD). Quantitative analyses of prior clinical trial data were used to inform future clinical trial designs.

METHODS

Two phase 2b studies in patients with AD and PsO were used to characterize the amount of topical brepocitinib and the resultant systemic trough concentration (CTrough) using a linear mixed-effects regression (LMER). This model was used to predict brepocitinib systemic CTrough for higher treated body surface areas (BSAs) in adults and children. Information from non-clinical and clinical trials with oral brepocitinib was leveraged to set safety thresholds. This combined approach was used to inform future dose-strength selection and treated BSA limits.

RESULTS

Data from 256 patients were analyzed. Patient type, dose strength, and frequency had significant impacts on the dose-exposure relationship. Systemic concentration in patients with PsO was predicted to be 45% lower than in patients with AD from the same dose. When topically applied to the same percentage BSA, brepocitinib systemic exposures are expected to be comparable between adults and children. The systemic steady-state exposure after 3% once daily and twice daily (2 mg/cm2) cream applied to less than 50% BSA in patients with AD and PsO, respectively, maintains at least a threefold margin to non-clinical safety findings and clinical hematologic markers.

CONCLUSION

The relationship between the amount of active drug applied and brepocitinib systemic CTrough, described by LMER, may inform the development strategy for dose optimization in the brepocitinib topical program.

Tofacitinib pharmacokinetics in children and adolescents with juvenile idiopathic arthritis

These analyses characterized tofacitinib pharmacokinetics (PKs) in children and adolescents with juvenile idiopathic arthritis (JIA). Data were pooled from phase I (NCT01513902), phase III (NCT02592434), and open-label, long-term extension (NCT01500551) studies of tofacitinib tablet/solution (weight-based doses administered twice daily [b.i.d.]) in patients with JIA aged 2 to less than 18 years. Population PK modeling used a nonlinear mixed-effects approach, with covariates identified using stepwise forward-inclusion backward-deletion procedures. Simulations were performed to derive dosing recommendations for children and adolescents with JIA. Two hundred forty-six pediatric patients were included in the population PK model. A one-compartment model with first-order elimination and absorption with body weight as a covariate for oral clearance and apparent volume of distribution sufficiently described the data. Oral solution was associated with comparable average concentration (Cavg) and slightly higher (113.9%) maximum concentration (Cmax) versus tablet, which was confirmed by a subsequent randomized, open-label, bioavailability study conducted in healthy adult participants (n = 12) by demonstrating adjusted geometric mean ratios (90% confidence interval) between oral solution and tablet of 1.04 (1.00-1.09) and 1.10 (1.00-1.21) for area under the curve extrapolated to infinity and Cmax, respectively (NCT04111614). A dosing regimen of 3.2 mg b.i.d. solution in patients 10 to less than 20 kg, 4 mg b.i.d. solution in patients 20 to less than 40 kg, and 5 mg b.i.d. tablet/solution in patients greater than or equal to 40 kg, irrespective of age, was proposed to achieve constant Cavg across weight groups. In summary, population PK characterization informed a simplified tofacitinib dosing regimen that has been implemented in pediatric patients with JIA.

Effects of Aripiprazole on Olanzapine Population Pharmacokinetics and Initial Dosage Optimization in Schizophrenia Patients

Objective

Olanzapine has already been used to treat schizophrenia patients; however, the initial dosage recommendation when multiple drugs are used in combination, remains unclear. The purpose of this study was to explore the drug-drug interaction (DDI) of multiple drugs combined with olanzapine and to recommend the optimal administration of olanzapine in schizophrenia patients.

Methods

In this study, we obtained olanzapine concentrations from therapeutic drug monitoring (TDM) database. In addition, related medical information, such as physiological, biochemical indexes, and concomitant drugs was acquired using medical log. Sixty-five schizophrenia patients were enrollmented for analysis using population pharmacokinetic model by means of nonlinear mixed effect (NONMEM).

Results

Weight and combined use of aripiprazole significantly affected olanzapine clearance. Without aripiprazole, for once-daily olanzapine administration dosages, 0.6, 0.5 mg/kg/day were recommended for 40-70, and 70-100 kg schizophrenia patients, respectively; for twice-daily olanzapine administration dosages, 0.6, 0.5 mg/kg/day were recommended for 40-60, and 60-100 kg schizophrenia patients, respectively. With aripiprazole, for once-daily olanzapine administration dosages, 0.4, 0.3 mg/kg/day were recommended for 40-53, and 53-100 kg schizophrenia patients, respectively; for twice-daily olanzapine administration dosages, 0.4 mg/kg/day was recommended for 40-100 kg schizophrenia patients, respectively.

Conclusion

Aripiprazole significantly affected olanzapine clearance, and when schizophrenia patients use aripiprazole, the olanzapine dosages need adjust. Meanwhile, we firstly recommended the optimal initial dosages of olanzapine in schizophrenia patients.

Pharmacokinetics of antitubercular drugs in patients hospitalized with HIV-associated tuberculosis: a population modeling analysis

Background

Early mortality among hospitalized HIV-associated tuberculosis (TB/HIV) patients is high despite treatment. The pharmacokinetics of rifampicin, isoniazid, and pyrazinamide were investigated in hospitalized TB/HIV patients and a cohort of outpatients with TB (with or without HIV) to determine whether drug exposures differed between groups.

Methods

Standard first-line TB treatment was given daily as per national guidelines, which consisted of oral 4-drug fixed-dose combination tablets containing 150 mg rifampicin, 75 mg isoniazid, 400 mg pyrazinamide, and 275 mg ethambutol. Plasma samples were drawn on the 3rd day of treatment over eight hours post-dose. Rifampicin, isoniazid, and pyrazinamide in plasma were quantified and NONMEM ® was used to analyze the data.

Results

Data from 60 hospitalized patients (11 of whom died within 12 weeks of starting treatment) and 48 outpatients were available. Median (range) weight and age were 56 (35 - 88) kg, and 37 (19 - 77) years, respectively. Bioavailability and clearance of the three drugs were similar between TB/HIV hospitalized and TB outpatients. However, rifampicin's absorption was slower in hospitalized patients than in outpatients; mean absorption time was 49.9% and 154% more in hospitalized survivors and hospitalized deaths, respectively, than in outpatients. Higher levels of conjugated bilirubin correlated with lower rifampicin clearance. Isoniazid's clearance estimates were 25.5 L/h for fast metabolizers and 9.76 L/h for slow metabolizers. Pyrazinamide's clearance was more variable among hospitalized patients. The variability in clearance among patients was 1.70 and 3.56 times more for hospitalized survivors and hospitalized deaths, respectively, than outpatients.

Conclusions

We showed that the pharmacokinetics of first-line TB drugs are not substantially different between hospitalized TB/HIV patients and TB (with or without HIV) outpatients. Hospitalized patients do not seem to be underexposed compared to their outpatient counterparts, as well as hospitalized patients who survived vs who died within 12 weeks of hospitalization.

A Comprehensive Mixed-Method Approach to Characterize the Source of Diurnal Tacrolimus Exposure Variability in Children: Systematic Review, Meta-analysis, and Application to an Existing Data Set

Tacrolimus is widely reported to display diurnal variation in pharmacokinetic parameters with twice-daily dosing. However, the contribution of chronopharmacokinetics versus food intake is unclear, with even less evidence in the pediatric population. The objectives of this study were to summarize the existing literature by meta-analysis and evaluate the impact of food composition on 24-hour pharmacokinetics in pediatric kidney transplant recipients. For the meta-analysis, 10 studies involving 253 individuals were included. The pooled effect sizes demonstrated significant differences in area under the concentration-time curve from time 0 to 12 hours (standardized mean difference [SMD], 0.27; 95% confidence interval [CI], 0.03-0.52) and maximum concentration (SMD, 0.75; 95% CI, 0.35-1.15) between morning and evening dose administration. However, there was significant between-study heterogeneity that was explained by food exposure. The effect size for minimum concentration was not significantly different overall (SMD, -0.09; 95% CI, -0.27 to 0.09) or across the food exposure subgroups. A 2-compartment model with a lag time, linear clearance, and first-order absorption best characterized the tacrolimus pharmacokinetics in pediatric participants. As expected, adding the time of administration and food composition covariates reduced the unexplained within-subject variability for the first-order absorption rate constant, but only caloric composition significantly reduced variability for lag time. The available data suggest food intake is the major driver of diurnal variation in tacrolimus exposure, but the associated changes are not reflected by trough concentrations alone.

A population pharmacokinetic model of methotrexate in Korean patients with haematologic malignancy

AIMS

This study was conducted to develop a population pharmacokinetic (PK) model of methotrexate in Korean patients with haematologic malignancy, identify factors affecting methotrexate PK, and propose an optimal dosage regimen for the Korean population.

METHODS

Data were retrospectively collected from 188 patients with acute leukaemia or non-Hodgkin's lymphoma who were admitted to Severance Hospital during the period from November 2005 to January 2016. Using demographic factors and laboratory results as potential covariates for PK parameters, model development was performed using NONMEM and optimal dosing regimens were developed using the final PK model.

RESULTS

A two-compartment model incorporating body weight via allometry best described the data, yielding typical parameter values of 25.09 L for central volume of distribution ( V 1 ), 17.65 L for peripheral volume of distribution ( V 2 ), 12.89 L/h for clearance (CL) and 0.655 L/h for inter-compartmental clearance in a 50 kg patient. Covariate analyses showed that, at the weight of 50 kg, CL decreased by 0.11 L/h for each 1-year increase in age above 14 years old and decreased 0.8-fold when serum creatinine level doubled, indicating the importance of age-specific dose individualization in methotrexate treatment. Volume of distribution at steady state derived from PK parameters (=  V 1  +  V 2 ) was 0.85 L/kg, which was similar to those in the Western or Chinese populations. Optimal doses simulated from the final model successfully produced the PK measures close to the target chosen.

CONCLUSIONS

The population PK model and optimal dosage regimens developed in this study can be used as a basis to achieve precision dosing in Korean patients with haematologic malignancy.

Model-Informed Vancomycin Dosing Optimization to Address Delayed Renal Maturation in Infants and Young Children with Critical Congenital Heart Disease

Ensuring safe and effective drug therapy in infants and young children often requires accounting for growth and organ development; however, data on organ function maturation are scarce for special populations, such as infants with congenital diseases. Children with critical congenital heart disease (CCHD) often require multiple staged surgeries depending on their age and disease severity. Vancomycin (VCM) is used to treat postoperative infections; however, the standard pediatric dose (60-80 mg/kg/day) frequently results in overexposure in children with CCHD. In this study, we characterized the maturation of VCM clearance in pediatric patients with CCHD and determined the appropriate dosing regimen using population pharmacokinetic (PK) modeling and simulations. We analyzed 1,254 VCM serum concentrations from 152 postoperative patients (3 days-13 years old) for population PK analysis. The PK model was developed using a two-compartment model with allometrically scaled body weight, estimated glomerular filtration rate (eGFR), and postmenstrual age as covariates. The observed clearance in patients aged ≤ 1 year and 1-2 years was 33% and 40% lower compared with that of non-CCHD patients, respectively, indicating delayed renal maturation in patients with CCHD. Simulation analyses suggested VCM doses of 25 mg/kg/day (age ≤ 3 months, eGFR 40 mL/min/1.73 m2 ) and 35 mg/kg/day (3 months < age ≤ 3 years, eGFR 60 mL/min/1.73 m2 ). In conclusion, this study revealed delayed renal maturation in children with CCHD, could be due to cyanosis and low cardiac output. Model-informed simulations identified the lower VCM doses for children with CCHD compared with standard pediatric guidelines.

Pharmacokinetics, Safety, and Tolerability of the Novel Tetrameric, High-Relaxivity, Macrocyclic Gadolinium-Based Contrast Agent Gadoquatrane in Healthy Adults

OBJECTIVES

Gadolinium (Gd)-based contrast agents are well established in clinical routine and have been proven safe and effective. However, there is a need for "next-generation" Gd-based contrast agents that would allow lowering the Gd dose used for routine contrast-enhanced magnetic resonance imaging procedures. The objective of this first-in-human study was to investigate the pharmacokinetic profile, safety, and tolerability of gadoquatrane, a novel high-relaxivity Gd-based contrast agent.

MATERIALS AND METHODS

This study was conducted in 2018/2019 as a prospective, randomized, single-blind, single-dose, placebo-controlled, escalating-dose study. Healthy volunteers were randomly assigned (6:2) to intravenous administration of gadoquatrane (0.025 to 0.2 mmol Gd/kg body weight) or placebo. Study procedures included collection of blood samples and excreta for pharmacokinetic analyses and safety assessments.

RESULTS

Forty-nine healthy study participants (mean age ± SD, 35 ± 6.3 years; 24 female) were evaluated. The effective half-life of gadoquatrane in plasma was short and similar in all dose groups (1.4-1.7 hours). Plasma concentrations around the lower quantitation limit (0.0318 μmol Gd/L) were reached 15-72 hours after administration. The volume of distribution at steady state was ~0.2 L/kg in all dose groups. The clearance (total and renal) was ~0.1 L/h per kilogram in all groups. Across dose groups, the exposure of gadoquatrane increased dose-proportionally. Metabolite profiling revealed no hint of degradation in vivo or release of free Gd. Seven of 36 participants (19.4%) receiving gadoquatrane and 4 of 13 participants (30.8%) receiving placebo experienced mild or moderate treatment-emergent adverse events. No serious adverse events occurred. The analysis of the Gd concentration-QTc interval relationship indicated no risk of QT/QTc prolongation (>10 milliseconds) with gadoquatrane at clinical dose levels.

CONCLUSIONS

Gadoquatrane with its high-relaxivity, pharmacokinetic similarity to established Gd-based contrast agents and high tolerability is a promising "next-generation" contrast agent for magnetic resonance imaging.

Is the current therapeutic dosage of nadroparin adequate for neonates and infants under 8 months with thromboembolic disease? a population pharmacokinetic study from a national children's medical center

Objectives: Nadroparin, a low-molecular-weight-heparin is commonly used off-label in neonates and infants for thromboembolic events prevention. However, the recommended dosing regimen often fails to achieve therapeutic target ranges. This study aimed to develop a population pharmacokinetic (PK) model of nadroparin to determine an appropriate dosing regimen for neonates and infants less than 8 months. Methods: A retrospective chart review was conducted on patients treated with nadroparin at Children's Hospital of Fudan University between July 2021 and December 2023. A population PK model was developed using anti-Xa levels, and its predictive performance was evaluated internally. Monte Carlo simulations were performed to design an initial dosing schedule targeting anti-Xa levels between 0.5 and 1 IU/mL. Results: A total of 40 neonates and infants aged less than 8 months with gestational age ranging from 25 to 41 weeks treated with nadroparin were enrolled in the study for analysis. A one-compartment PK model with first order absorption and elimination was adequately fitted to the data. Creatinine clearance was identified as a significant factor contributing to inter-individual variability in clearance. The typical population parameter estimates of clearance, distribution volume and absorption rate in this population were 0.211 L/h, 1.55 L and 0.495 h-1, respectively. Our findings suggest that current therapeutic doses of nadroparin (150-200 IU/kg q12 h) may result in subtherapeutic exposure, thus higher doses might be required. Conclusion: The present study offers the first estimation of PK parameters for nadroparin in preterm or term neonates and infants less than 8 months utilizing the model. Our findings have potential implications for recommending initial personalized dosages, particularly among patient populations exhibiting similar characteristics.

Fixed-Dose Combination Formulations in Solid Oral Drug Therapy: Advantages, Limitations, and Design Features

Whilst monotherapy is traditionally the preferred treatment starting point for chronic conditions such as hypertension and diabetes, other diseases require the use of multiple drugs (polytherapy) from the onset of treatment (e.g., human immunodeficiency virus acquired immunodeficiency syndrome, tuberculosis, and malaria). Successful treatment of these chronic conditions is sometimes hampered by patient non-adherence to polytherapy. The options available for polytherapy are either the sequential addition of individual drug products to deliver an effective multi-drug regimen or the use of a single fixed-dose combination (FDC) therapy product. This article intends to critically review the use of FDC drug therapy and provide an insight into FDC products which are already commercially available. Shortcomings of FDC formulations are discussed from multiple perspectives and research gaps are identified. Moreover, an overview of fundamental formulation considerations is provided to aid formulation scientists in the design and development of new FDC products.

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

BACKGROUND

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

AIM

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

DESIGN

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

REVIEW SOURCES

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

RESULTS

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

CONCLUSIONS

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

Dabigatran pharmacokinetic-pharmacodynamic in sheep: Informing dose for anticoagulation during cardiopulmonary bypass

BACKGROUND

The effect of the anticoagulant, dabigatran, and its antagonist, idarucizumab, on coagulation remains poorly quantified. There are few pharmacokinetic-pharmacodynamic data available to determine dabigatran dose in humans or animals undergoing cardiopulmonary bypass.

METHODS

Five sheep were given intravenous dabigatran 4 mg/kg. Blood samples were collected for thromboelastometric reaction time (R-time) and drug assay at 5, 15, 30, 60, 120, 240, 480 min, and 24 h. Plasma dabigatran concentrations and R-times were analyzed using an integrated pharmacokinetic-pharmacodynamic model using non-linear mixed effects. The impact of idarucizumab 15 mg/kg administered 120 min after dabigatran 4 mg/kg and its effect on R-time was observed.

RESULTS

A 2-compartment model described dabigatran pharmacokinetics with a clearance (CL 0.0453 L/min/70 kg), intercompartment clearance (Q 0.268 L/min/70 kg), central volume of distribution (V1 2.94 L/70 kg), peripheral volume of distribution (V2 9.51 L/70 kg). The effect compartment model estimates for a sigmoid EMAX model using Reaction time had an effect site concentration (Ce50 64.2 mg/L) eliciting half of the maximal effect (EMAX 180 min). The plasma-effect compartment equilibration half time (T1/2keo) was 1.04 min. Idarucizumab 15 mg/kg reduced R-time by approximately 5 min.

CONCLUSIONS

Dabigatran reversibly binds to the active site on the thrombin molecule, preventing activation of coagulation factors. The pharmacologic target concentration strategy uses pharmacokinetic-pharmacodynamic information to inform dose. A loading dose of dabigatran 0.25 mg/kg followed by a maintenance infusion of dabigatran 0.0175 mg/kg/min for 30 min and a subsequent infusion dabigatran 0.0075 mg/kg/min achieves a steady state target concentration of 5 mg/L in a sheep model.

A new population pharmacokinetic model for recombinant factor IX-Fc fusion concentrate including young children with haemophilia B

AIMS

Recombinant factor IX Fc fusion protein (rFIX-Fc) is an extended half-life factor concentrate administered to haemophilia B patients. So far, a population pharmacokinetic (PK) model has only been published for patients aged ≥12 years. The aim was to externally evaluate the predictive performance of the published rFIX-Fc population PK model for patients of all ages and develop a model that describes rFIX-Fc PK using real-world data.

METHODS

We collected prospective and retrospective data from patients with haemophilia B treated with rFIX-Fc and included in the OPTI-CLOT TARGET study (NTR7523) or United Kindom (UK)-EHL Outcomes Registry (NCT02938156). Predictive performance was assessed by comparing predicted with observed FIX activity levels. A new population PK model was constructed using nonlinear mixed-effects modelling.

RESULTS

Real-world data were obtained from 37 patients (median age: 16 years, range 2-71) of whom 14 were aged <12 years. Observed FIX activity levels were significantly higher than levels predicted using the published model, with a median prediction error of -48.8%. The new model showed a lower median prediction error (3.4%) and better described rFIX-Fc PK, especially for children aged <12 years. In the new model, an increase in age was correlated with a decrease in clearance (P < .01).

CONCLUSIONS

The published population PK model significantly underpredicted FIX activity levels. The new model better describes rFIX-Fc PK, especially for children aged <12 years. This study underlines the necessity to strive for representative population PK models, thereby avoiding extrapolation outside the studied population.

Preliminary pharmacokinetics and patient experience of jet-injected dexmedetomidine in healthy adults

Jet injection is a drug delivery system without a needle. A compressed liquid drug formulation pierces the skin, depositing the drug into the subcutaneous or intramuscular tissues. We investigated the pharmacokinetics and patient experience of dexmedetomidine administered using jet injection in six healthy adult study participants. This needleless jet injection device was used to administer dexmedetomidine 0.5 μg/kg to the subcutaneous tissues overlying the deltoid muscle. Serum concentrations of dexmedetomidine were assayed at approximately 5 minutes, 15 minutes, 30 minutes, 1 hour and 4 hours after administration. Pharmacokinetic interrogation of concentration time profiles estimated an absorption half time for jet-injected dexmedetomidine of 21 minutes (coefficient of variation 69.4%) with a relative bioavailability assumed unity. In our samples the measured median peak (range) concentration was 0.164 μg/l (0.011-0.325 μg/l), observed in the sample taken at a median (range) of 13.5 minutes (11-30 minutes). The Richmond agitation sedation scale was used to assess the sedative effect, and scored 0 (alert and calm) or -1 (drowsy) in all participants. Five of the six participants stated they would prefer jet injection to needle injection in the future and one had no preference. The findings suggest that the use of a larger dose (>2 μg/kg) would be required to achieve the clinically relevant target concentration of 1 μg/l necessary to achieve deeper sedation (Richmond agitation sedation scale ≤3).

Optimizing dosing of the cycloserine pro-drug terizidone in children with rifampicin-resistant tuberculosis

There are no pharmacokinetic data in children on terizidone, a pro-drug of cycloserine and a World Health Organization (WHO)-recommended group B drug for rifampicin-resistant tuberculosis (RR-TB) treatment. We collected pharmacokinetic data in children <15 years routinely receiving 15-20 mg/kg of daily terizidone for RR-TB treatment. We developed a population pharmacokinetic model of cycloserine assuming a 2-to-1 molecular ratio between terizidone and cycloserine. We included 107 children with median (interquartile range) age and weight of 3.33 (1.55, 5.07) years and 13.0 (10.1, 17.0) kg, respectively. The pharmacokinetics of cycloserine was described with a one-compartment model with first-order elimination and parallel transit compartment absorption. Allometric scaling using fat-free mass best accounted for the effect of body size, and clearance displayed maturation with age. The clearance in a typical 13 kg child was estimated at 0.474 L/h. The mean absorption transit time when capsules were opened and administered as powder was significantly faster compared to when capsules were swallowed whole (10.1 vs 72.6 min) but with no effect on bioavailability. Lower bioavailability (-16%) was observed in children with weight-for-age z-score below -2. Compared to adults given 500 mg daily terizidone, 2022 WHO-recommended pediatric doses result in lower exposures in weight bands 3-10 kg and 36-46 kg. We developed a population pharmacokinetic model in children for cycloserine dosed as terizidone and characterized the effects of body size, age, formulation manipulation, and underweight-for-age. With current terizidone dosing, pediatric cycloserine exposures are lower than adult values for several weight groups. New optimized dosing is suggested for prospective evaluation.

Semi-mechanistic population pharmacokinetic/pharmacodynamic modeling of a Plasmodium elongation factor 2 inhibitor cabamiquine for prevention and cure of malaria

Cabamiquine is a novel antimalarial agent that demonstrates the potential for chemoprevention and treatment of malaria. In this article, the dose-exposure-response relationship of cabamiquine was characterized using a population pharmacokinetic (PK)/pharmacodynamic (PD) model, incorporating the effects of cabamiquine on parasite dynamics at the liver and blood stages of malaria infection. Modeling was performed sequentially. First, a three-compartmental population PK model was developed, comprising linear elimination, a transit absorption model in combination with first-order absorption, and a recirculation model. Second, this model was expanded into a PK/PD model using parasitemia data from an induced blood stage malaria (IBSM) human challenge model. To describe the parasite growth and killing in the blood, a turnover model was used. Finally, the liver stage parasite dynamics were characterized using data from a sporozoite challenge model (SpzCh), and system parameters were fixed based on biological plausibility. Cabamiquine concentration in the central compartment was used to drive parasite killing at the blood and liver stages. Blood stage minimum inhibitory concentrations (MICb) were estimated at 7.12 ng/mL [95% confidence interval (CI95%): 6.26-7.88 ng/mL] and 1.28 ng/mL (CI95%: 1.12-1.43 ng/mL) for IBSM and SpzCh populations, respectively, while liver stage MICl was lower (0.61 ng/mL; CI95%: 0.24-0.96 ng/mL). In conclusion, a population PK/PD model was developed by incorporating parasite dynamics and drug activity at the blood and liver stages based on clinical data and biological knowledge. This model can potentially facilitate antimalarial agent development by supporting the efficient selection of the optimal dosing regimen.

Pharmacokinetics and Proposed Dosing of Levetiracetam in Children With Obesity

OBJECTIVE

Characterize levetiracetam pharmacokinetics (PK) in children with obesity to inform dosing.

METHODS

Children 2 to <21 years old receiving standard of care oral levetiracetam across two opportunistic studies provided blood samples. Levetiracetam plasma PK data were analyzed with a nonlinear mixed-effects modeling approach. Indirect measures for body size and covariates were tested for model inclusion. Individual empirical Bayesian estimates using the final model parameters were compared by obesity status. Monte Carlo simulation using total body weight was performed in children with normal estimated glomerular filtration rate to identify dosing for children with obesity that resulted in comparable exposures to normal weight adults and children after receiving label dosing.

RESULTS

The population PK model was developed from 341 plasma concentrations from 169 children. A 1-compartment model best fit the data with fat-free mass as a significant covariate. Compared with children with normal weight, children with obesity had significantly lower body weight-normalized clearance (median [range], 4.77 [1.49-10.44] and 3.71 [0.86-13.55] L/h/70 kg, respectively). After label dosing with the oral formulation in children with obesity 4 to <16 years old, maximum and minimum steady-state concentrations were higher (25% and 41%, respectively [oral solution] and 27% and 19%, respectively [tablet]) compared with children with normal weight. Comparable exposures between children with and without obesity were achieved with weight-tiered dosing regimens of <75 kg or ≥75 kg.

CONCLUSIONS

Weight-tiered dosing for levetiracetam oral solution and tablets for children with obesity 4 to <16 years old results in more comparable exposures to children of normal weight.

Ceftriaxone Pharmacokinetics and Pharmacodynamics in 2 Pediatric Patients on Extracorporeal Membrane Oxygenation Therapy

BACKGROUND

Critically ill patients with cardiac or respiratory failure may require extracorporeal membrane oxygenation (ECMO). Antibiotics are frequently administered when the suspected cause of organ failure is an infection. Ceftriaxone, a β-lactam antibiotic, is commonly used in patients who are critically ill. Although studies in adults on ECMO have suggested minimal impact on ceftriaxone pharmacokinetics, limited research exists on ceftriaxone pharmacokinetics/pharmacodynamics (PK/PD) in pediatric ECMO patients. We report the PK profiles and target attainment of 2 pediatric patients on ECMO who received ceftriaxone.

METHODS

Ceftriaxone concentrations were measured in 2 pediatric patients on ECMO using scavenged opportunistic sampling. PK profiles were generated and individual PK parameters were estimated using measured free ceftriaxone concentrations and a published population PK model in children who are critically ill, using Bayesian estimation.

RESULTS

Patient 1, an 11-year-old boy on venovenous ECMO for respiratory failure received 2 doses of 52 mg/kg ceftriaxone 12 hours apart while on ECMO and additional doses every 12 hours off ECMO. On ECMO, ceftriaxone clearance was 13.0 L/h/70 kg compared with 7.6 L/h/70 kg off ECMO, whereas the model-predicted mean clearance in children who are critically ill without ECMO support was 6.54 L/h/70 kg. Patient 2, a 2-year-old boy on venoarterial ECMO due to cardiac arrest received 50 mg/kg ceftriaxone every 12 hours while on ECMO for >7 days. Only clearance while on ECMO could be estimated (9.1 L/h/70 kg). Trough concentrations in both patients were >1 mg/L (the breakpoint for Streptococcus pneumoniae ) while on ECMO.

CONCLUSIONS

ECMO increased ceftriaxone clearance above the model-predicted clearances in the 2 pediatric patients studied. Twelve-hour dosing allowed concentrations to remain above the breakpoint for commonly targeted bacteria but not 4 times the breakpoint in one patient, suggesting that precision dosing may be beneficial to ensure target attainment in children on ECMO.

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.

Abacavir Drug Exposures in African Children Under 14 kg Using Pediatric Solid Fixed Dose Combinations According to World Health Organization Weight Bands

BACKGROUND

The pharmacokinetics of abacavir (ABC) in African children living with HIV (CLHIV) weighing <14 kg and receiving pediatric fixed dose combinations (FDC) according to WHO weight bands dosing are limited. An ABC population pharmacokinetic model was developed to evaluate ABC exposure across different World Health Organization (WHO) weight bands.

METHODS

Children enrolled in the LIVING study in Kenya and Uganda receiving ABC/lamivudine (3TC) dispersible tablets (60/30 mg) according to WHO weight bands. A population approach was used to determine the pharmacokinetic parameters. Monte Carlo simulations were conducted using an in silico population with demographic characteristics associated with African CLHIV. ABC exposures (AUC0-24) of 6.4-50.4 mg h/L were used as targets.

RESULTS

Plasma samples were obtained from 387 children. A 1-compartment model with allometric scaling of clearance (CL/F) and volume of distribution (V/F) according to body weight best characterized the pharmacokinetic data of ABC. The maturation of ABC CL/F was characterized using a sigmoidal Emax model dependent on postnatal age (50% of adult CL/F reached by 0.48 years of age). Exposures to ABC were within the target range for children weighing 6.0-24.9 kg, but children weighing 3-5.9 kg were predicted to be overexposed.

CONCLUSIONS

Lowering the ABC dosage to 30 mg twice daily or 60 mg once daily for children weighing 3-5.9 kg increased the proportion of children within the target and provided comparable exposures. Further clinical study is required to investigate clinical implications and safety of the proposed alternative ABC doses.

Dolutegravir pharmacokinetics in Ugandan patients with TB and HIV receiving standard- versus high-dose rifampicin

Higher rifampicin doses may improve tuberculosis treatment outcomes. This could however exacerbate the existing drug interaction with dolutegravir. Moreover, the metabolism of dolutegravir may also be affected by polymorphism of UGT1A1, a gene that codes for uridine diphosphate glucuronosyltransferase. We used population pharmacokinetic modeling to compare the pharmacokinetics of dolutegravir when coadministered with standard- versus high-dose rifampicin in adults with tuberculosis and HIV, and investigated the effect of genetic polymorphisms. Data from the SAEFRIF trial, where participants were randomized to receive first-line tuberculosis treatment with either standard- 10 mg/kg or high-dose 35 mg/kg rifampicin alongside antiretroviral therapy, were used. The dolutegravir model was developed with 211 plasma concentrations from 44 participants. The median (interquartile range) rifampicin area under the curve (AUC) in the standard- and high-dose arms were 32.3 (28.7-36.7) and 153 (138-175) mg·h/L, respectively. A one-compartment model with first-order elimination and absorption through transit compartments best described dolutegravir pharmacokinetics. For a typical 56 kg participant, we estimated a clearance, absorption rate constant, and volume of distribution of 1.87 L/h, 1.42 h-1, and 12.4 L, respectively. Each 10 mg·h/L increase in the AUC of coadministered rifampicin from 32.3 mg·h/L led to a 2.3 (3.1-1.4) % decrease in dolutegravir bioavailability. Genetic polymorphism of UGT1A1 did not significantly affect dolutegravir pharmacokinetics. Simulations of trough dolutegravir concentrations show that the 50 mg twice-daily regimen attains both the primary and secondary therapeutic targets of 0.064 and 0.3 mg/L, respectively, regardless of the dose of coadministered rifampicin, unlike the once-daily regimen.

Population pharmacokinetic model of cefepime for critically ill adults: a comparative assessment of eGFR equations

Cefepime exhibits highly variable pharmacokinetics in critically ill patients. The purpose of this study was to develop and qualify a population pharmacokinetic model for use in the critically ill and investigate the impact of various estimated glomerular filtration rate (eGFR) equations using creatinine, cystatin C, or both on model parameters. This was a prospective study of critically ill adults hospitalized at an academic medical center treated with intravenous cefepime. Individuals with acute kidney injury or on kidney replacement therapy or extracorporeal membrane oxygenation were excluded. A nonlinear mixed-effects population pharmacokinetic model was developed using data collected from 2018 to 2022. The 120 included individuals contributed 379 serum samples for analysis. A two-compartment pharmacokinetic model with first-order elimination best described the data. The population mean parameters (standard error) in the final model were 7.84 (0.24) L/h for CL1 and 15.6 (1.45) L for V1. Q was fixed at 7.09 L/h and V2 was fixed at 10.6 L, due to low observed interindividual variation in these parameters. The final model included weight as a covariate for volume of distribution and the eGFRcr-cysC (mL/min) as a predictor of drug clearance. In summary, a population pharmacokinetic model for cefepime was created for critically ill adults. The study demonstrated the importance of cystatin C to prediction of cefepime clearance. Cefepime dosing models which use an eGFR equation inclusive of cystatin C are likely to exhibit improved accuracy and precision compared to dosing models which incorporate an eGFR equation with only creatinine.

Pharmacokinetics and optimized dosing of dispersible and non-dispersible levofloxacin formulations in young children

BACKGROUND

Levofloxacin is used for treatment and prevention of rifampicin-resistant (RR)-TB in children. Recent data showed higher exposures with 100 mg dispersible compared with non-dispersible tablet formulations with potentially important dosing implications in children. We aimed to verify and better characterize this finding.

METHODS

We conducted a crossover pharmacokinetic trial in children aged ≤5 years receiving levofloxacin RR-TB preventive therapy. Pharmacokinetic sampling was done after 15-20 mg/kg doses of levofloxacin with 100 mg dispersible and crushed 250 mg non-dispersible levofloxacin formulations. A population pharmacokinetic model was developed.

RESULTS

Twenty-five children were included, median (IQR) weight and age 12.2 (10.7-15.0) kg and 2.56 (1.58-4.03) years, respectively. A two-compartment model with first-order elimination and transit compartment absorption best described levofloxacin pharmacokinetics. Allometric scaling adjusted for body size, and maturation of clearance with age was characterized. Typical clearance in a 12 kg child was estimated at 4.17 L/h. Non-dispersible tablets had 21.5% reduced bioavailability compared with the dispersible formulation, with no significant differences in other absorption parameters.Dosing simulations showed that current recommended dosing for both formulations result in median exposures below adult-equivalent exposures at a 750 mg daily dose, mainly in children >6 months. Higher levofloxacin doses of 16-30 mg/kg for dispersible and 20-38 mg/kg for crushed non-dispersible tablets may be required in children >6 months.

CONCLUSIONS

The dispersible paediatric levofloxacin formulation has improved bioavailability compared with the crushed non-dispersible adult formulation, but exposures remain below those in adults. We propose optimized age- and weight-based dosing for levofloxacin, which require further evaluation.

Population pharmacokinetic and pharmacokinetic-pharmacodynamic modeling of bempedoic acid and low-density lipoprotein cholesterol in healthy subjects and patients with dyslipidemia

Population pharmacokinetics (popPK) of bempedoic acid and the popPK/pharmacodynamic (popPK/PD) relationship between bempedoic acid concentrations and serum low-density lipoprotein cholesterol (LDL-C) from baseline were characterized. A two-compartment disposition model with a transit absorption compartment and linear elimination best described bempedoic acid oral pharmacokinetics (PK). Multiple covariates, including renal function, sex, and weight, had statistically significant effects on the predicted steady-state area under the curve. Mild (estimated glomerular filtration rate (eGFR) 60 to < 90 mL/min vs. ≥ 90 mL/min) and moderate (eGFR 30 to < 60 mL/min vs. ≥ 90 mL/min) renal impairment, female sex, low (< 70 kg vs. 70-100 kg) and high (> 100 kg vs. 70-100 kg) body weight were predicted to have a 1.36-fold (90% confidence interval (CI) 1.32, 1.41), 1.85-fold (90% CI 1.74, 2.00), 1.39-fold (90% CI 1.34, 1.47), 1.35-fold (90% CI 1.30, 1.41), and 0.75-fold (90% CI 0.72, 0.79) exposure difference relative to their reference populations, respectively. An indirect response model described changes in serum LDL-C with a model-predicted 35% maximal reduction and bempedoic acid IC50 of 3.17 µg/mL. A 28% reduction from LDL-C baseline was predicted for a steady-state average concentration of 12.5 µg/mL after bempedoic acid (180 mg/day) dosing, accounting for approximately 80% of the predicted maximal LDL-C reduction. Concurrent statin therapy, regardless of intensity, reduced the maximal effect of bempedoic acid but resulted in similar steady-state LDL-C levels. While multiple covariates had statistically significant effects on PK and LDL-C lowering, none were predicted to warrant bempedoic acid dose adjustment.

Population pharmacokinetics of esaxerenone, a novel non-steroidal mineralocorticoid receptor blocker, in patients with essential hypertension, patients with diabetic nephropathy, and healthy volunteers

OBJECTIVES

Esaxerenone is a novel, non-steroidal mineralocorticoid receptor (MR) blocker with improved selectivity and affinity for MR. The objectives of this study were to model the population pharmacokinetics of esaxerenone in a diverse population and to evaluate the effect of covariates on pharmacokinetics parameters.

METHODS

A total of 8263 plasma esaxerenone concentrations from 166 healthy volunteers, 1097 hypertensive patients and 360 patients with diabetic nephropathy were pooled. A three-compartment model with sequential zero- and first-order absorption was used to describe the time-courses of plasma esaxerenone following single and multiple doses once daily for up to 12 weeks. Covariate effects were estimated using the full covariate modeling approach. Clinical relevance of covariates was ascertained using tornado plots.

RESULTS

Esaxerenone was estimated to have high bioavailability (85.3%), low clearance (3.28 L/h) and relatively large distribution volume at steady state (94.8 L). Body weight (-26 to +36%) and coadministration of itraconazole (+64%) or rifampicin (-68%) were associated with a greater influence on esaxerenone exposure.

CONCLUSIONS

The most influential covariates on esaxerenone exposure were coadministrations of itraconazole and rifampicin, followed by body weight. The clinical relevance of effects of renal impairment, mild to moderate hepatic impairment, and age is limited.

Cefepime pharmacokinetics in critically ill children and young adults undergoing continuous kidney replacement therapy

OBJECTIVES

Cefepime is an antibiotic commonly used to treat sepsis and is cleared by renal excretion. Cefepime dosing requires adjustment in patients with decreased kidney function and in those receiving continuous kidney replacement therapy (CKRT). We aimed to characterize cefepime PK in a diverse cohort of critically ill paediatric patients on CKRT.

METHODS

Patients were identified from an ongoing pharmacokinetic/pharmacodynamic (PK/PD) study of beta-lactam antibiotics, and were included if they had received at least two cefepime doses in the ICU and were on CKRT for at least 24 h. PK parameters were estimated using MwPharm++ with Bayesian estimation and a paediatric population PK model. Target attainment was assessed as time of free cefepime concentrations above minimum inhibitory concentration (fT > 1× or 4 × MIC).

RESULTS

Seven patients were included in the study (ages 2 to 20 years). CKRT indications included liver failure (n = 1), renal failure (n = 4) and fluid overload (n = 2). Total effluent flow rates ranged from 1833 to 3115 (mean 2603) mL/1.73 m2/h, while clearance was 2.11-3.70 (mean 3.0) L/h/70 kg. Effluent flows were lower, but clearance and fT > MIC were similar to paediatric data published previously. Using Pseudomonas aeruginosa MIC breakpoints, all patients had 100% of dosing interval above MIC, but only one had 100% of dosing interval above 4× MIC.

CONCLUSIONS

Since most patients failed to attain stringent targets of 100% fT > 4×  MIC, model-informed precision dosing may benefit such patients.

Nonlinear Mixed Effects Modeling of Glucagon Kinetics in Healthy Subjects

OBJECTIVE

To date, the lack of a model of glucagon kinetics precluded the possibility of estimating and studying glucagon secretion in vivo, e.g., using deconvolution, as done for other hormones like insulin and C-peptide. Here, we used a nonlinear mixed effects technique to develop a robust population model of glucagon kinetics, able to describe both the typical population kinetics (TPK) and the between-subject variability (BSV), and relate this last to easily measurable subject characteristics.

METHODS

Thirty-four models of increasing complexity (variably including covariates and correlations among random effects) were identified on glucagon profiles obtained from 53 healthy subjects, who received a constant infusion of somatostatin to suppress endogenous glucagon production, followed by a continuous infusion of glucagon (65 ng/kg/min). Model selection was performed based on its ability to fit the data, provide precise parameter estimates, and parsimony criteria.

RESULTS

A two-compartment model was the most parsimonious. The model was able to accurately describe both the TPK and the BSV of model parameters as function of body mass and body surface area. Parameters were precisely estimated, with central volume of distribution V1 = 5.46 L and peripheral volume of distribution V2 = 5.51 L. The introduction of covariates resulted in a significant shrinkage of the unexplained BSV and considerably improved the model fit.

CONCLUSION

We developed a robust population model of glucagon kinetics.

SIGNIFICANCE

This model provides a deeper understanding of glucagon kinetics and is usable to estimate glucagon secretion in vivo by deconvolution of plasma glucagon concentration data.

Model-Informed Support of Dose Selection for Prophylactic Treatment with Dalcinonacog Alfa in Adult and Paediatric Hemophilia B Patients

INTRODUCTION

Dalcinonacog alfa (DalcA), a novel subcutaneously administered recombinant human factor IX (FIX) variant is being developed for adult and paediatric patients with hemophilia B (HB). DalcA has been shown to raise FIX to clinically meaningful levels in adults with HB. This work aimed to support dosing regimen selection in adults and perform first-in-paediatric dose extrapolations using a model-based pharmacokinetic (PK) approach.

METHODS

A population PK model was built using adult data from two clinical trials (NCT03186677, NCT03995784). With allometry in the model, clinical trial simulations were performed to study alternative dosing regimens in adults and children. Steady-state trough levels and the time-to-reach target were derived to inform dose selection.

RESULTS

Almost 90% of the adults were predicted to achieve desirable FIX levels, i.e. 10% FIX activity, following daily 100 IU/kg dosing, with 90% of the subjects reaching target within 1.6-7.1 days. No every-other-day regimen met the target. A dose of 125 IU/kg resulted in adequate FIX levels down to 6 years, whereas a 150 IU/kg dose was needed below 6 down to 2 years of age. For subjects down to 6 years that did not reach target with 125 IU/kg, a dose escalation to 150 IU/kg was appropriate. The children below 6 to 2 years were shown to need a dose escalation to 200 IU/kg if 150 IU/kg given daily was insufficient.

CONCLUSION

This study supported the adult dose selection for DalcA in the presence of sparse data and enabled first-in-paediatric dose selection to achieve FIX levels that reduce risk of spontaneous bleeds.

Initial dosage optimisation of cyclosporine in Chinese paediatric patients undergoing allogeneic haematopoietic stem cell transplantation based on population pharmacokinetics: a retrospective study

OBJECTIVE

Improved understanding of cyclosporine A (CsA) pharmacokinetics in children undergoing allogeneic haematopoietic stem cell transplantation (allo-HSCT) is crucial for effective prevention of acute graft-versus-host disease and medication safety. The aim of this study was to establish a population pharmacokinetic (Pop-PK) model that could be used for individualised therapy to paediatric patients undergoing allo-HSCT in China.

DESIGN, SETTING AND PARTICIPANTS

A retrospective analysis of 251 paediatric HSCT patients who received CsA intravenously in the early post transplantation period at Women and Children's Medical Center in Guangzhou was conducted.

ANALYSIS MEASURES

The model building dataset from 176 children was used to develop and analyse the CsA Pop-Pk model by using the nonlinear mixed effect model method. The basic information was collected by the electronic medical record system. Genotype was analysed by matrix-assisted time-of-flight mass spectrometry. The stability and predictability of the final model were verified internally, and a validation dataset of 75 children was used for external validation. Monte Carlo simulation is used to adjust and optimise the initial dose of CsA in paediatric allo-HSCT patients.

RESULTS

The typical values for clearance (CL) and volume of distribution ([Formula: see text]) were 14.47 L/hour and 2033.53 L, respectively. The body weight and haematocrit were identified as significant variables for V, while only body weight had an impact on CL. The simulation based on the final model suggests that paediatrics with HSCT required an appropriate intravenous dose of 5 mg/kg/day to reach the therapeutic trough concentration.

CONCLUSIONS

The CsA Pop-PK model established in this study can quantitatively describe the factors influencing pharmacokinetic parameters and precisely predict the intrinsic exposure to CsA in children. In addition, our dosage simulation results can provide evidence for the personalised medications TRIAL REGISTRATION NUMBER: ChiCTR2000040561.

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.

Model-informed pediatric dose selection of marzeptacog alfa (activated): An exposure matching strategy

Marzeptacog alfa (activated) (MarzAA) is an activated recombinant human rFVII variant intended for subcutaneous (s.c.) administration to treat or prevent bleeding in individuals with hemophilia A (HA) or B (HB) with inhibitors, and other rare bleeding disorders. The s.c. administration provides benefits over i.v. injections. The objective of the study was to support the first-in-pediatric dose selection for s.c. MarzAA to treat episodic bleeding episodes in children up through 11 years in a registrational phase III trial. Assuming the same exposure-response relationship as in adults, an exposure matching strategy was used with a population pharmacokinetics model. A sensitivity analysis evaluating the impact of doubling in absorption rate and age-dependent allometric exponents on dose selection was performed. Subsequently, the probability of trial success, defined as the number of successful trials for a given pediatric dose divided by the number of simulated trials (n = 1000) was studied. A successful trial was defined as outcome where four, three, or two out of 24 pediatric subjects per trial were allowed to fall outside the adult exposures after s.c. administration of 60 μg/kg. A dose of 60 μg/kg in children with HA/HB was supported by the clinical trial simulations to match exposures in adults. The sensitivity analyses further supported selection of the 60 μg/kg dose level in all age groups. Moreover, the probability of trial success evaluations given a plausible design confirmed the potential of a 60 μg/kg dose level. Taken together, this work demonstrates the utility of model-informed drug development and could be helpful for other pediatric development programs for rare diseases.

Population pharmacokinetic analyses of tacrolimus in non-transplant patients: a systematic review

BACKGROUND AND OBJECTIVES

Tacrolimus (TAC) has been increasingly used in patients with non-transplant settings. Because of its large between-subject variability, several population pharmacokinetic (PPK) studies have been performed to facilitate individualized therapy. This review summarized published PPK models of TAC in non-transplant patients, aiming to clarify factors affecting PKs of TAC and identify the knowledge gap that may require further research.

METHODS

The PubMed, Embase databases, and Cochrane Library, as well as related references, were searched from the time of inception of the databases to February 2023, to identify TAC population pharmacokinetic studies modeled in non-transplant patients using a non-linear mixed-effects modeling approach.

RESULTS

Sixteen studies, all from Asian countries (China and Korea), were included in this study. Of these studies, eleven and four were carried out in pediatric and adult patients, respectively. One-compartment models were the commonly used structural models for TAC. The apparent clearance (CL/F) of TAC ranged from 2.05 to 30.9 L·h^-1 (median of 14.9 L·h^-1). Coadministered medication, genetic factors, and weight were the most common covariates affecting TAC-CL/F, and variability in the apparent volume of distribution (V/F) was largely explained by weight. Coadministration with Wuzhi capsules reduced CL/F by about 19 to 43%. For patients with CYP3A5*1*1 and *1*3 genotypes, the CL/F was 39-149% higher CL/F than patients with CYP3A5*1*1.

CONCLUSION

The optimal TAC dosage should be adjusted based on the patient's co-administration, body weight, and genetic information (especially CYP3A5 genotype). Further studies are needed to assess the generalizability of the published models to other ethnic groups. Moreover, external validation should be frequently performed to improve the clinical practicality of the models.

Estimation of FMO3 Ontogeny by Mechanistic Population Pharmacokinetic Modelling of Risdiplam and Its Impact on Drug-Drug Interactions in Children

BACKGROUND AND OBJECTIVE

Spinal muscular atrophy (SMA) is a progressive neuromuscular disease caused by insufficient levels of survival motor neuron (SMN) protein. Risdiplam (EvrysdiTM) increases SMN protein and is approved for the treatment of SMA. Risdiplam has high oral bioavailability and is primarily eliminated through hepatic metabolism by flavin-containing monooxygenase3 (FMO3) and cytochrome P450 (CYP) 3A, by 75% and 20%, respectively. While the FMO3 ontogeny is critical input data for the prediction of risdiplam pharmacokinetics (PK) in children, it was mostly studied in vitro, and robust in vivo FMO3 ontogeny is currently lacking. We derived in vivo FMO3 ontogeny by mechanistic population PK modelling of risdiplam and investigated its impact on drug-drug interactions in children.

METHODS

Population and physiologically based PK (PPK and PBPK) modelling conducted during the development of risdiplam were integrated into a mechanistic PPK (Mech-PPK) model to estimate in vivo FMO3 ontogeny. A total of 10,205 risdiplam plasma concentration-time data from 525 subjects aged 2 months-61 years were included. Six different structural models were examined to describe the in vivo FMO3 ontogeny. Impact of the newly estimated FMO3 ontogeny on predictions of drug-drug interaction (DDI) in children was investigated by simulations for dual CYP3A-FMO3 substrates including risdiplam and theoretical substrates covering a range of metabolic fractions (fm) of CYP3A and FMO3 (fmCYP3A:fmFMO3 = 10%:90%, 50%:50%, 90%:10%).

RESULTS

All six models consistently predicted higher FMO3 expression/activity in children, reaching a maximum at the age of 2 years with an approximately threefold difference compared with adults. Different trajectories of FMO3 ontogeny in infants < 4 months of age were predicted by the six models, likely due to limited observations for this age range. Use of this  in vivo FMO3 ontogeny function improved prediction of risdiplam PK in children compared to in vitro FMO3 ontogeny functions. The simulations of theoretical dual CYP3A-FMO3 substrates predicted comparable or decreased CYP3A-victim DDI propensity in children compared to adults across the range of fm values. Refinement of FMO3 ontogeny in the risdiplam model had no impact on the previously predicted low CYP3A-victim or -perpetrator DDI risk of risdiplam in children.

CONCLUSION

Mech-PPK modelling successfully estimated in vivo FMO3 ontogeny from risdiplam data collected from 525 subjects aged 2 months-61 years. To our knowledge, this is the first investigation of in vivo FMO3 ontogeny by population approach using comprehensive data covering a wide age range. Derivation of a robust in vivo FMO3 ontogeny function has significant implications on the prospective prediction of PK and DDI in children for other FMO3 substrates in the future, as illustrated in the current study for FMO3 and/or dual CYP3A-FMO3 substrates.

CLINICAL TRIAL REGISTRY NUMBERS

NCT02633709, NCT03032172, NCT02908685, NCT02913482, NCT03988907.