Population Pharmacokinetic Model Research Articles

A Population Pharmacokinetic Model of Tocilizumab in Kidney Transplant Patients Treated for Chronic Active Antibody-Mediated Rejection: Comparison of Plasma Exposure Between Intravenous and Subcutaneous Administration Schemes.

Tocilizumab prevents the clinical worsening of chronic active antibody-mediated rejection (CAAMR) in kidney transplant recipients. Following a global shortage of the intravenous pharmaceutical form in 2022, patients were switched from monthly intravenous administration of 8 mg/kg to weekly subcutaneous injection of 162 mg, raising the question of bioequivalence between these schemes of administration. We aimed to compare the areas under the curve (AUC) of tocilizumab in virtual simulations of populations treated with the two administration schemes and to identify the covariates that could contribute to pharmacokinetic variability of tocilizumab in kidney transplant patients with CAAMR who received tocilizumab as salvage treatment. This retrospective monocentric study included 43 kidney transplant patients (202 tocilizumab concentrations) with CAAMR treated with intravenous or subcutaneous tocilizumab between December 2020 and January 2023. We developed a population pharmacokinetic model using nonlinear mixed effects modeling and identified the covariates that could contribute to tocilizumab AUC variability. Monte Carlo simulations were then performed to assess the subcutaneous and intravenous tocilizumab AUC for 0-28 days (M1), 56-84 days (M3), 140-168 days (M6), and 308-336 days (M12). Bioequivalence was defined by SC/IV AUC geometric mean ratios (GMRs) between 0.80 and 1.25. A two-compartment model with parallel linear and nonlinear elimination best described the concentration-time data. Significant covariates for tocilizumab clearance were body weight, urinary albumin-to-creatinine ratio (ACR), and inflammation status [C-reactive protein (CRP) ≥5 mg/L]. The GMR values and their 90% confidence intervals at M3, M6, and M12 were within the 0.8-1.25 margin for equivalence. Conversely, the 90% prediction intervals of the GMR were much wider than the 90% confidence intervals and did not fall within 0.8 and 1.25. From month 3 of treatment, the subcutaneous and intravenous tocilizumab administration schemes provided average bioequivalent pharmacokinetic exposure at the population level but not at the individual level. Body weight, inflammation, ACR, and administration scheme should be considered to personalize the dose of tocilizumab for patients with CAAMR. Further studies are required to determine the target of tocilizumab exposure in kidney transplant patients with CAAMR.

Population Pharmacokinetics of Loxoprofen and its alcoholic metabolites in healthy Korean men.

Loxoprofen has been actively used clinically to relieve musculoskeletal pain and inflammatory symptoms. However, there are few reports on quantitative pharmacokinetic (PK) prediction tools and diversity analyzes for loxoprofen within populations. The aim of this study was to identify effective covariates associated with explaining inter-individual PK variability through a population pharmacokinetic (Pop-PK) modeling approach for loxoprofen, and to provide a starting point for establishing scientific dosing regimens. The bioequivalence PK results of loxoprofen performed on 52 healthy Korean men and the physiological and biochemical parameters derived from each individual were used as base data for the development of a Pop-PK model of loxoprofen. In order to simultaneously predict the PKs of the active form according to loxoprofen exposure, previously reported PK results of trans-alcohol loxoprofen, an active metabolite of loxoprofen, were used to expand the model. The Pop-PK profiles of loxoprofen were described in terms of the basic structure of a non-sequential two absorption with 2-disposition compartment, and for inter-individual PK variations, peripheral compartment volume of distribution could be correlated with body surface area (BSA), and central compartment clearance with creatinine clearance (CrCL) and albumin levels. As a result of the model simulation, the concentrations of loxoprofen and its alcoholic metabolites in plasma significantly decreased as CrCL and albumin levels increased and decreased, respectively. On the other hand, it was confirmed that the higher the BSA, the greater the distribution of loxoprofen to the periphery, and the minimum concentrations of loxoprofen and alcoholic metabolites in plasma in steady-state increased by approximately 1.78-2 times, while the fluctuation between maximum and minimum concentrations decreased. The results suggest that patients with large BSA, impaired renal function, and high serum albumin levels may have significantly higher plasma exposure to loxoprofen and trans-alcohol loxoprofen. It was also suggested that the potential side effects in the gastrointestinal system and various tissues and the level of exposure in plasma due to long-term application of loxoprofen in this patient group could be causally explained. This study provides a very useful starting point for a scientific precision medicine approach to loxoprofen by discovering effective covariates and establishing a quantitative model that can explain the diversity of loxoprofen PKs within the population. The clinical study protocol used in this study was thoroughly reviewed and approved by the Institutional Review Board of the Institute of Bioequivalence and Bridging Study, Chonnam National University, Gwangju, Republic of Korea. The bioequivalence study permit numbers are as follows: 041113; 10.15.2004.

Population Pharmacokinetics of Trimethoprim/Sulfamethoxazole: Dosage Optimization for Patients with Renal Insufficiency or Receiving Continuous Renal Replacement Therapy.

The goal of the study was to describe the population pharmacokinetics of trimethoprim, sulfamethoxazole, and N-acetyl sulfamethoxazole in hospitalized patients. Furthermore, this study used the model to optimize dosing regimens of cotrimoxazole for Pneumocystis jirovecii pneumonia and in patients with renal insufficiency or with continuous renal replacement therapy (CRRT). This was a retrospective multicenter observational cohort study based on therapeutic drug monitoring (TDM) data from hospitalized patients treated with cotrimoxazole. We developed two population pharmacokinetic (POPPK) models: a model of trimethoprim and an integrated model with both sulfamethoxazole and N-acetyl sulfamethoxazole concentrations. Monte Carlo simulations were performed to determine the optimal dosing regimen. A total of 348 measurements from 168 patients were available. The estimated glomerular filtration rate (eGFR) and CRRT were included as covariates on the clearance of all three compounds. Cotrimoxazole TID 1,920 mg and b.i.d. 2,400 mg led to sufficient exposure for infections with P. jirovecii in patients without renal insufficiency. To reach equivalent exposure, a dose reduction of 33.3% is needed in patients with an eGFR of 10 mL/minute/1.73 m2 and of 16.7% for an eGFR of 30 mL/minute/1.73 m2. N-acetyl sulfamethoxazole accumulates in patients with a reduced eGFR. CRRT increased the clearance of sulfamethoxazole, but not trimethoprim or N-acetyl sulfamethoxazole, compared with the median clearance in the population. Doubling the sulfamethoxazole dose is needed for patients on CRRT to reach equivalent exposure.

Can edoxaban be used at extremes of bodyweight and in patients with a creatinine clearance ≥95 ml/min? – A population pharmacokinetic analysis

BackgroundClinical evidence surrounding edoxaban use in patients weighing <50 kg and >120 kg is lacking. The International Society of Thrombosis and Haemostasis Scientific and Standardisation Committee suggests avoiding edoxaban in patients >120 kg. Additionally, concerns exist regarding decreased efficacy in patients prescribed edoxaban for atrial fibrillation with a creatinine clearance (CrCl) >95 ml/min, a finding of the ENGAGE AF-TIMI 48 trial when edoxaban was compared to warfarin. ObjectiveTo derive a population pharmacokinetic (PopPK) model using clinical practice data, to understand the impact of bodyweight and renal function on edoxaban pharmacokinetics. MethodEdoxaban plasma concentrations and patient characteristics were collated from King's College Hospital anticoagulation clinics between 11/2016 and 08/2022. A PopPK model was developed using non-linear mixed effects modelling and used to simulate edoxaban concentrations at the extremes of bodyweight and with varying renal function. ResultsData from 409 patients (46 < 50 kg, 34 > 120 kg and 123 with a CrCl > 95 ml/min) provided 455 edoxaban plasma concentrations. A one-compartment model with between-subject variability on clearance with a proportional error model best described the data. The most significant covariates impacting on edoxaban exposure were CrCl and bodyweight. Our work suggests that edoxaban exposure in patients weighing up to 140 kg is comparable to those weighing 75 kg. Edoxaban exposure is reduced in patients weighing <50 kg due to the recommended dose reductions. There is also a reduction in AUCss when CrCl > 95 ml/min compared to CrCl 80 ml/min. ConclusionsOur population PK model for edoxaban suggests that renal function is a key driver for overall edoxaban exposure. Further clinical outcome data is required to understand clinical effectiveness and adverse outcomes.

Prediction of the Intra-T Lymphocyte Tacrolimus Concentration after Kidney Transplantation with Population Pharmacokinetic Modeling.

The intracellular tacrolimus concentration in CD3+ T lymphocytes is proposed to be a better representative of the active component of tacrolimus than the whole blood concentration. However, intracellular measurements are complicated. Therefore, the aim of this study was to describe the relationship between intracellular and whole blood tacrolimus concentrations in a population pharmacokinetic model. Twenty-eight de novo kidney transplant recipients, treated with a once-daily oral extended-release tacrolimus formulation, were followed during the first-month post-transplantation. Additional whole blood and intracellular tacrolimus concentrations were measured at day 6 ± 1 (pre-dose, 4 and 8 hours post-dose) and day 14 ± 3 (pre-dose) post-transplantation. Pharmacokinetic analysis was performed using nonlinear mixed effects modeling software (NONMEM). The ratio between intracellular (n = 109) and whole blood (n = 248) concentrations was best described by a two-compartment whole blood model with an additional intracellular compartment without mass transfer from the central compartment. The ratio remained stable over time. Prednisolone dose influenced the absorption rate of tacrolimus, while hemoglobin, CYP3A4*22 allele carrier, and CYP3A5 expresser status were associated with the oral clearance of tacrolimus (P-value < 0.001). Furthermore, the intracellular tacrolimus concentrations were correlated with the intracellular production of interleukin-2 (P-value 0.015). The intracellular tacrolimus concentration can be predicted from a measured whole blood concentration using this model, without the need for repeated intracellular measurements. This knowledge is particularly important when the intracellular concentration is ready to be implemented into clinical practice, to overcome the complexities of cell isolation and analytical methods.

Characterization of mavacamten pharmacokinetics in patients with hypertrophic cardiomyopathy to inform dose titration.

Mavacamten is a selective, allosteric, reversible cardiac myosin inhibitor that has been developed for the treatment of adults with symptomatic obstructive hypertrophic cardiomyopathy (HCM). A population pharmacokinetic (PopPK) model was developed to characterize mavacamten pharmacokinetics (PK) and the variation in mavacamten exposure associated with intrinsic and extrinsic factors. Data from 12 clinical studies (phases 1, 2, and 3) were used. Evaluable participants were those who had at least one mavacamten concentration measurement with associated sampling time and dosing information. The base model included key covariates: body weight, cytochrome P450 isozyme 2C19 (CYP2C19) phenotype with respect to PK, and formulation. The final model was generated using stepwise covariate testing and refinement processes. Simulations were performed to evaluate PK: apparent clearance (CL/F); apparent central and peripheral volumes of distribution; and steady-state average, trough, and maximum concentrations. Overall, 9244 measurable PK observations from 497 participants were included. A two-compartment model structure was selected. After stepwise covariate model building and refinement, additional covariates included were: specified mavacamten dose, omeprazole or esomeprazole administration, health/disease status, estimated glomerular filtration rate, fed status, and sex. The final PopPK model accurately characterized mavacamten concentrations. At any given dose, CYP2C19 phenotype was the most influential covariate on exposure parameters (e.g., median CL/F was reduced by 72% in CYP2C19:poor metabolizers compared with the reference participant [CYP2C19:normal metabolizer]). CL/F was also approximately 16% higher in women than in men but lower in participants receiving concomitant omeprazole or esomeprazole (by 33% and 42%, respectively) than in participants not receiving such concomitant therapy.

Model-informed Evidence for Clinical Non-inferiority of Every-2-Weeks Versus Standard Weekly Dosing Schedule of Cetuximab in Metastatic Colorectal Cancer.

Cetuximab was initially developed and approved as a first-line treatment in patients with unresectable metastatic colorectal cancer (mCRC) for weekly administration (250 mg/m2 Q1W with 400 mg/m2 loading dose). An every-2-weeks schedule (500 mg/m2 Q2W) was approved recently by several health authorities. Being synchronized with chemotherapy, Q2W administration should improve patients' convenience and healthcare resource utilization. Herein, we present evidence of non-inferiority of Q2W cetuximab, compared with Q1W dosing using pharmacometrics modeling and clinical trial simulation (CTS). Pooled data from five phase I-III clinical trials in 852 patients with KRAS wild-type mCRC treated with Q1W or Q2W cetuximab were modeled using a population exposure-tumor size (TS) model linked to overall survival (OS); exposure was derived from a previously established population pharmacokinetic model. A semi-mechanistic TS model adapted from the Claret model incorporated killing rate proportional to cetuximab area under the concentration-time curve over 2 weeks (AUC) with Eastern Cooperative Oncology Group (ECOG) status as covariate on baseline TS. The OS was modeled with Weibull hazard using ECOG, baseline TS, primary tumor location, and predicted percent change in TS at 8 weeks as covariates. Model-based simulations revealed indistinguishable early tumor shrinkage and survival between Q2W vs. Q1W cetuximab. CTS evaluated OS non-inferiority (predefined margin of 1.25) in 1,000 trials, each with 2,000 virtual patients receiving Q2W or Q1W cetuximab (1:1), and demonstrated non-inferiority in 94% of cases. Taken together, these analyses provide model-based evidence for clinical non-inferiority of Q2W vs. Q1W cetuximab in mCRC with potential benefits to patients and healthcare providers.

Pharmacokinetics of Dexamethasone in Children and Adolescents with Obesity.

Dexamethasone is a synthetic glucocorticoid approved for treating disorders of various organ systems in both adult and pediatric populations. Currently, approved pediatric dosing recommendations are weight-based, but it is unknown whether differences in dexamethasone drug disposition and exposure exist for children with obesity. This study aimed to develop a population pharmacokinetic (PopPK) model for dexamethasone with data collected from children with obesity. Dexamethasone was given as either IV or oral/enteral administration, and a salt factor correction was used for dexamethasone sodium phosphate injection. A PopPK analysis using dexamethasone plasma concentration versus time was performed using the software NONMEM. A virtual population of 1000 children with obesity across three age groups was generated for dosing simulations. Data from 59 study participants with 82 PK plasma samples were used in the PopPK analysis. A one-compartment model with first-order absorption and the inclusion of total body weight as a covariate characterized the data. No other covariates were included in the PopPK model. Single and multiple IV dose(s) of 0.5 and 1 mg/kg every 8 h resulted in 68% or more of virtual children with obesity attaining simulated exposures that were within exposure ranges previously reported in adult studies. In conclusion, this was the first study to characterize dexamethasone's PopPK in children with obesity. Simulation results suggest that virtual children with obesity receiving oral doses of 0.5 and 1 mg/kg had generally comparable dexamethasone exposures as adult estimates. Additional studies are needed to characterize the dexamethasone's target exposure in children.

Pharmacokinetic modeling of sulfamethoxazole-trimethoprim and sulfadiazine-trimethoprim combinations in broilers

Sulfonamides (S) are old bacteriostatic antibiotics which are widely prescribed in combination with trimethoprim (TMP) for the treatment of various diseases in food-producing animals such as poultry. Nowadays, the 1:5 dose ratio of TMP/S used in broilers is a direct transposition of the ratio determined in Human decades ago for TMP/sulfamethoxazole (SMX), aiming to obtain a supposed synergistic plasma concentration ratio of 1:19. However, major pharmacokinetics (PK) differences exist according to the sulfonamide used in the combination. Here, we generated new PK data in broilers after a cross-over design with IV and the oral administration of 2 major sulfonamides, sulfadiazine (SDZ) and SMX, in combination with TMP, and analyzed the data via a population pharmacokinetic (popPK) modeling approach. Results showed that TMP has a greater plasma to tissue distribution than both sulfonamides with a higher volume of distribution (0.51 L/kg for SDZ, 0.62 L/kg for SMX and 3.14 L/kg for TMP). SMX has the highest elimination half-life (2.83 h) followed by SDZ and TMP (2.01 h and 1.49 h, respectively). The oral bioavailability of the 3 molecules was approximately 100%. Bodyweight could explain some of the inter-individual variability in the volume of distribution of SDZ and SMX and the clearance of SDZ and TMP, as heavier broilers have higher typical values. Monte Carlo simulations of a large virtual broiler population (n = 1,000) showed that the targeted plasma ratio of TMP:S of 1:19 was rarely or never reached at the individual level for both combinations at the marketed doses and greatly varies over time and between individuals, questioning the relevance of the 1:5 dose ratio for current formulations of TMP/S.

Nedosiran population pharmacokinetic and pharmacodynamic modelling and simulation to guide clinical development and dose selection in patients with primary hyperoxaluria type 1.

To characterize pharmacokinetic and pharmacodynamic profiles of nedosiran in patients with primary hyperoxaluria type 1 (PH1), identify influential covariates and confirm therapeutic doses. A population pharmacokinetic (PK)/pharmacodynamic (PD) (POP-PKPD) model was developed to characterize the concentration-time course of nedosiran and the corresponding effect on 24-h urinary oxalate (Uox). Simulations of dosing to achieve clinically meaningful reduction in Uox in children, adolescents and adults with PH1 were performed. Analyses included PK data from 143 healthy participants and PH1/PH2 patients, and PD data from 46 PH1 patients. Nedosiran PK was described by a two-compartment model with dual n-transit absorption and parallel linear and nonlinear elimination. The relationship between nedosiran exposure and Uox was described by an indirect response model. Body weight, estimated glomerular filtration rate (eGFR) and disease status were identified as influential covariates for the POP-PK model. The simulation results supported a weight-banded dosing regimen of nedosiran sodium in adolescents and adults (≥12 years) with PH1 of 170 mg (weight ≥50 kg) and 136 mg (weight <50 kg), in children (6-11 years) with PH1 of 3.5mg/kg, and no dose adjustments for PH1 patients with relatively preserved kidney function (eGFR ≥ 30 mL/min/1.73m2). Following the proposed dosing regimens, the simulated median fold-changes in PK AUC0-τ,ss were acceptable (≤1.51 fold-change) and ~71% of PH1 patients across all age groups achieved near-normal Uox (<0.6mmol) by week 52. The final POP-PKPD model characterizes observed nedosiran PK and Uox data. Simulations support nedosiran dosing regimens in PH1 patients aged ≥6years with relatively preserved kidney function.

The use of quantitative clinical pharmacology approaches to support moxidectin dosing recommendations in lactation.

Moxidectin is approved by the US Food and Drug Administration (US FDA) for the treatment of onchocerciasis (river-blindness) due to Onchocerca volvulus in patients aged 12 years and older. In onchocerciasis-endemic areas, mass drug administration (MDA) programs with ivermectin, with or without vector control, aim to control the disease, reduce morbidity, interrupt transmission, and more recently, achieve elimination. Moxidectin has the potential to be used in MDA programs. In countries where onchocerciasis is endemic, infants are often breastfed up to the age of 2 years, suggesting that some women are likely to be lactating during such periodic MDA programs. Quantitative analyses of non-clinical and clinical data using non-compartmental analysis and population based pharmacokinetic (popPK) modeling as well as physiologically based pharmacokinetic modeling (PBPK) were performed to determine the amount of moxidectin excreted in breast milk and subsequent exposures in the infant. The results of the analyses were similar. Concentrations of moxidectin in breast milk followed a similar pattern to those in plasma, with maximum concentrations occurring approximately 4 hours after dosing followed by a rapid decline in both breast milk and plasma. As early as two days after dosing, concentrations of moxidectin in breast milk were below the threshold for acceptable daily intake levels established by the European Medicines Agency (EMA) and FDA for secondary exposures from veterinary use, and below the WHO recommended relative infant dose (RID) safety threshold. The analyses were conducted to support prescribers and policy makers on dosing recommendations for moxidectin in lactation.

Lamivudine dosing for preterm infants exposed to HIV: a population pharmacokinetic modelling and simulation study.

To develop a pragmatic twice daily lamivudine dosing strategy for preterm infants from 24 to 37 completed weeks of gestation. Data were combined from eight pharmacokinetic studies in neonates and infants receiving lamivudine oral solution. A population pharmacokinetic model was developed using non-linear mixed effects regression. Different lamivudine dosing strategies, stratified by gestational age at birth (GA) bands, were simulated in a virtual population of preterm infants, aimed at maintaining lamivudine drug exposures (AUC0-12) within a reference target range of 2.95 to 13.25 µg·h/mL, prior to switching to WHO-weight band doses when ≥4 weeks of age and weighing ≥3 kg. A total of 154 infants (59% female) contributed 858 lamivudine plasma concentrations. Median (range) GA at birth was 38 (27-41) weeks. At the time of first pharmacokinetic sampling infants were older with median postnatal age (PNA) of 6.3 (0.52-26.6) weeks. Lamivudine concentrations were described by a one-compartment model, with CL/F and V/F allometrically scaled to weight. Maturation of CL/F was described using an Emax model based on PNA. CL/F was also adjusted on GA to allow extrapolation for extreme prematurity. Simulations predicted an optimal lamivudine dosing for infants GA ≥24 to <30 weeks of 2 mg/kg twice daily from birth until weighing 3 kg; and for GA ≥30 to <37 weeks, 2 mg/kg twice daily for the first 4 weeks of life, followed by 4 mg/kg twice daily until weighing 3 kg. Model-based predictions support twice daily pragmatic GA band dosing of lamivudine for preterm infants, but clinical validation is warranted.

Identifying optimal dosing strategies for meropenem in the paediatric intensive care unit through modelling and simulation.

Meropenem, a β-lactam antibiotic commonly prescribed for severe infections, poses dosing challenges in critically ill patients due to highly variable pharmacokinetics. We sought to develop a population pharmacokinetic model of meropenem for critically ill paediatric and young adult patients. Paediatric intensive care unit patients receiving meropenem 20-40 mg/kg every 8 h as a 30 min infusion were prospectively followed for clinical data collection and scavenged opportunistic plasma sampling. Nonlinear mixed effects modelling was conducted using Monolix®. Monte Carlo simulations were performed to provide dosing recommendations against susceptible pathogens (MIC ≤ 2 mg/L). Data from 48 patients, aged 1 month to 30 years, with 296 samples, were described using a two-compartment model with first-order elimination. Allometric body weight scaling accounted for body size differences. Creatinine clearance and percentage of fluid balance were identified as covariates on clearance and central volume of distribution, respectively. A maturation function for renal clearance was included. Monte Carlo simulations suggested that for a target of 40% fT > MIC, the most effective dosing regimen is 20 mg/kg every 8 h with a 3 h infusion. If higher PD targets are considered, only continuous infusion regimens ensure target attainment against susceptible pathogens, ranging from 60 mg/kg/day to 120 mg/kg/day. We successfully developed a population pharmacokinetic model of meropenem using real-world data from critically ill paediatric and young adult patients with an opportunistic sampling strategy and provided dosing recommendations based on the patients' renal function and fluid status.

Population pharmacokinetic modeling of sotatercept in healthy participants and patients with pulmonary arterial hypertension.

Sotatercept is a breakthrough, first-in-class biologic, that is FDA-approved for the treatment of pulmonary arterial hypertension (PAH). A population pharmacokinetic (PopPK) model was developed using data from two phase 1 studies in healthy participants, and two phase 2 studiesand one phase 3 study in participants with PAH. The pooled sotatercept PK data encompassed single intravenous (IV) or subcutaneous (SC) doses ranging from 0.01 to 3.0 mg/kg, as well as multiple SC doses ranging from 0.03 to 1.0 mg/kg, with PK samples collected up to a maximum of ~150 weeks following Q3W and Q4W dosing regimens. The final PopPK analysis included 350 participants, with 30 and 320 participants receiving sotatercept IV and SC, respectively. A two-compartment model with a first-order absorption rate constant and a linear disposition from central compartment well-described sotatercept PK. The estimated bioavailability is ~66%; bioavailability, clearance (CL), and central volume (VC) have low to moderate inter-individual variability. Time-varying body weight and baseline albumin concentration were statistically significant predictors of PK; CL and VC were predicted to increase with increasing body weight, while CL was predicted to decrease with increasing baseline albumin concentration. However, the magnitude of covariate effects is not predicted to meaningfully alter the disposition of sotatercept. Altogether, the PopPK modeling results demonstrate favorable PK characteristics (low to moderate variability and typical bioavailability), supporting sotatercept as a SC biological agent for the treatment of patients with PAH.

A Model-Informed Drug Development Approach to Design a Phase 3 Trial of Teverelix Drug Product in Advanced Prostate Cancer Patients with Increased Cardiovascular Risk.

Teverelix drug product (DP) is a parenteral gonadotropin-releasing hormone (GnRH) antagonist that has been successfully tested in phase 2 trials for hormone-sensitive advanced prostate cancer (APC) and benign prostatic hyperplasia (BPH). In previous APC trials, teverelix DP was administered as intramuscular (IM) and subcutaneous (SC) injections, using a loading dose and (in a single trial) a maintenance dose. Our objective was to derive an optimal dosing regimen for phase 3 clinical development, using a pharmacometrics modeling approach. Data from 9 phase 2 studies (229 patients) was utilized to develop a population pharmacokinetic (PK) model that described the concentration profile accommodating both IM and SC routes of administration. A 2-compartment model with sequential first-order absorption (slow and fast) and lag times best described the PK profiles of teverelix following SC and IM administration. An indirect response model with inhibition of production rate was fit to describe testosterone (T) concentrations based on physiological relevance. The final population PK-pharmacodynamic model was used to conduct simulations of various candidate dosing regimens to select the optimal dosing regimen to achieve clinical castration (T<0.5ng/mL by day 28) and to sustain clinical castration for 26 weeks. Model simulation showed that a loading dose of 360mg SC and 180mg IM with a maintenance dose of 360mg SC 6-weekly (Q6W) starting at day 28 can achieve a ≥95% castration rate up to 52 weeks. This dose regimen was selected for phase 3 clinical development, which includes cardiovascular safety assessment in comparison to a GnRH agonist.

Nicotine Gum in Thai Smokers with Different CYP2A6 Enzymes: A Population Pharmacokinetic Analysis

Objective: Despite the popularity of nicotine gum in Thailand, population pharmacokinetics of nicotine gum in the Thai population has not been investigated. This study aimed to develop a population pharmacokinetic (POPPK) model of nicotine and to quantify the effects of genetic and non-genetic factors to nicotine pharmacokinetics. Materials and Methods: Secondary data collected from a previous clinical trial assessing cytochrome P450 2A6 (CYP2A6) genotypes in Thai smokers was investigated. Eighteen participants who had received a single dose of 2 mg nicotine gum were included. Blood samples were collected before, at 0.25, 0.5, 1, 1.5, 2, 2.5, 3, 4.5 and 6 hours after nicotine administration. POPPK analysis was performed using nonlinear mixed effect modelling. Results: One-compartment with 1st order elimination and absorption with 6-transit compartments best described the data. CYP2A6 enzyme activity was a significant covariate on the nicotine clearance. Apparent elimination clearance (CL/F) for a person with 100% CYP2A6 activity was 266.0 L/h. CL/F would be 36.0 L/h in a subject with 0% CYP2A6 activity. However, the impact of non-genetic factors (monthly alcohol consumption, Fagerstrom Test for Nicotine Dependence score and the number of cigarettes per day) on pharmacokinetics of nicotine were not found. Conclusion: This first report on population pharmacokinetics of nicotine gum in Thai smokers provided the pharmacokinetic model and quantified CL/F for smokers with different CYP2A6 genotypes. A markedly lower exposure to nicotine in the Thai population compared to others highlights the need for more studies to ensure the efficacy of nicotine gum in the Thai population.

Population Pharmacokinetics and Target Attainment of Allopurinol and Oxypurinol Before, During, and After Cardiac Surgery with Cardiopulmonary Bypass in Neonates with Critical Congenital Heart Disease.

The CRUCIAL trial (NCT04217421) is investigating the effect of postnatal and perioperative administration of allopurinol on postoperative brain injury in neonates with critical congenital heart disease (CCHD) undergoing cardiac surgery with cardiopulmonary bypass (CPB) shortly after birth. This study aimed to characterize the pharmacokinetics (PK) of allopurinol and oxypurinol during the preoperative, intraoperative, and postoperative phases in this population, and to evaluate target attainment of the current dosing strategy. Nonlinear mixed-effects modeling was used to develop population PK models in 14 neonates from the CRUCIAL trial who received up to five intravenous allopurinol administrations throughout the postnatal and perioperative periods. Target attainment was defined as achieving an allopurinol concentration >2 mg/L in at least two-thirds of the patients during the first 24h after birth and between the start and 36h after cardiac surgery with CPB. A two-compartment model for allopurinol was connected to a one-compartment model for oxypurinol with an auto-inhibition effect on the conversion, which best described the PK. In a typical neonate weighing 3.5kg who underwent cardiac surgery at a postnatal age (PNA) of 5.6 days, the clearance (CL) of allopurinol and oxypurinol at birth was 0.95 L/h (95% confidence interval 0.75-1.2) and 0.21 L/h (0.17-0.27), respectively, which subsequently increased with PNA to 2.97 L/h and 0.41 L/h, respectively, before CPB. During CPB, allopurinol and oxypurinol CL decreased to 1.38 L/h (0.9-1.87) and 0.12 L/h (0.05-0.22), respectively. Post-CPB, allopurinol CL increased to 2.21 L/h (1.74-2.83), while oxypurinol CL dropped to 0.05 L/h (0.01-0.1). Target attainment was 100%, 53.8%, and 100% at 24h postnatally, 24h after the start of CPB, and 36h after the end of cardiac surgery, respectively. The combined concentrations of allopurinol and oxypurinol maintained ≥ 90% inhibition of xanthine oxidase (IC90XO) throughout the postnatal and perioperative period. The minimal target concentration of allopurinol was not achieved at every predefined time interval in the CRUCIAL trial; however, the dosing strategy used was deemed adequate, since it yielded concentrations well exceeding the IC90XO. The decreased CL of both compounds during CPB suggests influence of the hypothermia, hemofiltration, and the potential sequestration of allopurinol in the circuit. The reduced CL of oxypurinol after CPB is likely attributable to impaired kidney function.

Effect of Blood Product Resuscitation on Cefazolin Pharmacokinetics in Trauma Patients.

Background: Antibiotics are frequently administered prophylactically to trauma patients with various injury patterns to prevent infectious complications. Trauma patients may also require large volume resuscitation with blood products. Limited data are available to support antibiotic dosing recommendations in this population. We hypothesized that we would be able to develop a population pharmacokinetic model of cefazolin, a frequently used antibiotic in the trauma scenario, from remnant blood samples by pharmacokinetic analysis of trauma patients. Methods: Remnant plasma from standard of care chemistry/hematology assessments was retrieved within 48 h of collection and assayed to determine cefazolin concentrations. Population pharmacokinetic analyses were conducted in Pmetrics using R. Linear regression was conducted to assess the effect of blood product resuscitation volume on cefazolin pharmacokinetic parameters. Results: Cefazolin concentrations best fitted a two-compartment model (Akaike information criterion: 443.9). The mean ± standard deviation parameters were total body clearance (4.3 ± 1.9L), volume of the central compartment (Vc: 7.7 ± 6.9L), and intercompartment transfer constants (k12: 1.3 ± 0.98 h-1, k21: 0.6 ± 0.45 h-1). No statistical relationships were observed between blood products, volume of blood products, and cefazolin clearance or Vc (R2: 0.0004-0.21, p = 0.08-0.95). Using a 5,000-patient Monte Carlo simulation, 2 g with repeated dosing every 2 h until end of surgery was required to achieve 93.2% probability of 100% free time above the minimum inhibitory concentration (MIC) (fT > MIC) at the ECOFF value for Staphylococcus aureus (2 mg/L). Conclusions: In these 15 trauma patients receiving blood transfusion, no relationship with blood volume resuscitation and cefazolin pharmacokinetics was observed. On the basis of this pharmacokinetic model, frequent cefazolin doses are required to maintain 100% fT > MIC.

Population Pharmacokinetics of Capivasertib in Patients with Advanced or Metastatic Solid Tumours.

Overactivation of the PI3K/AKT pathway can occur in many cancers. Capivasertib is a potent, selective pan-AKT inhibitor. The objectives of this analysis were to develop a population pharmacokinetic model for capivasertib and to quantitatively assess the impact of intrinsic and extrinsic factors on the pharmacokinetics of capivasertib. Pharmacokinetic data from four phase I and II studies were combined. Capivasertib was administered orally at a dose range of 80-800mg twice daily over 28-day and 21-day cycles as monotherapy or in combination with paclitaxel or fulvestrant, using continuous dosing or one of two intermittent dosing schedules: either 4 days on, 3 days off (4/3) or 2days on, 5days off (2/5). Several models and approaches were tested for their ability to describe capivasertib disposition. The covariates assessed included dose, schedule, age, body weight, race, sex, creatinine clearance, hepatic function, renal function, smoking status, food effect, formulation, and concomitant use with paclitaxel, fulvestrant, cytochrome P450, family 3, subfamily A (CYP3A)inducers, CYP3A inhibitors and acid-reducing agents. A total of 3963 capivasertib plasma concentrations from 441 patients were included. Capivasertib pharmacokinetics was adequately described by a three-compartment model where the apparent clearance (CL/F) presented a moderate time-dependent and dose-dependent clearance. Following oral administration of multiple doses of capivasertib (400 mg twice daily; [4/3]), the initial CL/F was 62.2 L/h (between-subject variability 39.3%), and after approximately 120 hours, CL/F decreased by 18%. The effective half-life was 8.34 h. Steady state was predicted to be reached on every third and fourth dosing day each week from the second week with exposure levels that produced robust inhibition of AKT but not of other related kinases. The area under the plasma concentration-time curve and maximum plasma concentration of capivasertib were proportional between the dose levels of 80-480mg after multiple doses but more than proportional beyond 480mg. Schedule, age, race, sex, creatinine clearance, hepatic function, renal function, smoking status and concomitant use with fulvestrant, CYP3A inducers, CYP3A inhibitors or acid-reducing agents were not significant covariates for capivasertib pharmacokinetics. Concomitant use of paclitaxel, food effect and formulation statistically significantly affected capivasertib pharmacokinetics, but the effect was low. Body weight was statistically significantly related to capivasertib CL/F, with a 12% reduction in CL/F at steady state and a 14% increase in the area under the curve for 12 hours at steady state and maximum concentration at steady state at a lower body weight (47 kg vs 67 kg reference). Capivasertib pharmacokinetics showed moderate between-subject variability, and most covariates assessed had no significant impact. Body weight, dose, concomitant use of paclitaxel, food effect and formulation showed statistically significant effects. However, these were predicted to impact exposure to capivasertib by <20% and were not expected to be clinically relevant. Based on the population pharmacokinetics, no a priori dose adjustment is needed for intrinsic and extrinsic factors.

Prospective Trial on the Pharmacokinetics of Clopidogrel in Hemodialysis Patients

IntroductionPatients undergoing hemodialysis (HDP) exhibit extensive cardiovascular risk. The widely prescribed anti-platelet agent clopidogrel is metabolically activated by cytochrome enzymes, which may be impaired by uremia and chronic low-grade inflammation, typically present in HDP. We conducted a prospective multicenter study to investigate the pharmacokinetics and pharmacodynamics of clopidogrel in HDP and healthy volunteers (HV). MethodsWe enrolled HDP receiving long-term clopidogrel (75mg) and pantoprazole treatment (40 mg). HV received a loading dose of 300 mg clopidogrel, followed by 75 mg once daily. Pantoprazole, a substrate and probe drug of CYP2C19, was administered intravenously (40 mg). Plasma concentrations were quantified by mass spectrometry. Pharmacokinetics were calculated, and a population pharmacokinetic model was developed. The primary endpoint was the maximum concentration of clopidogrel’s active metabolite (Cmax_CAM). Platelet aggregation was measured using ADP-induced whole-blood aggregometry. ResultsSeventeen HDP and 16 HV were included. Cmax_CAM was significantly lower in HDP compared to HV (median [IQR] 12.2 [4.6-23.4] vs. 24.7 [17.8-36.5] ng/mL, p=0.02). The Cmax ratio of CAM to prodrug was 8.5-fold lower in HDP, and an 82.7% reduced clopidogrel clearance, including CAM formation, was found using population pharmacokinetic modeling. ADP-induced platelet aggregation at 120 min was significantly higher in HDP than in HV (median [IQR]: 26 [14-43] vs. 12 [11-18] U, p=0.004. Pantoprazole terminal half-life was ∼1.7-fold higher in HDPs compared to HVs. ConclusionOur data demonstrate an altered metabolism of clopidogrel in HDP in the context of lower CYP2C19 activity, with potential implications for other substances metabolized by this enzyme.