Population Pharmacokinetic-pharmacodynamic Model Research Articles

Nitric Oxide Donor Sodium Nitroprusside Reduces Racemic Ketamine-But Not Esketamine-Induced Pain Relief.

The anesthetic, analgesic and antidepressant drug ketamine produces dissociation with symptoms of psychosis and anxiety, an effect attributed to neuronal nitric oxide depletion following N-methyl-d-aspartate blockade. There is evidence that dissociation induced by racemic ketamine, containing both ketamine enantiomers (S- and R-ketamine) but not esketamine (the S-isomer) is inhibited by nitric oxide (NO) donor sodium nitroprusside (SNP). We tested whether a similar intervention would reduce racemic and esketamine-induced analgesia in a randomized double-blind placebo-controlled trial. Seventeen healthy volunteers were treated with 0.5 μg.kg-1.min-1 SNP or placebo during a 3-h infusion of escalating doses of racemic ketamine (total dose 140 mg) or esketamine (70 mg). Pain pressure threshold (PPT) and arterial blood samples for measurement of S- and R-ketamine and their metabolites, S- and R-norketamine, were obtained. The data were analyzed with a population pharmacokinetic-pharmacodynamic model that incorporated the measured S- and R- ketamine and S- and R-norketamine isomers as input and PPT as output to the model. The potency of the 2 formulations in increasing PPT from baseline by 100% was 0.47 ± 0.12 (median ± standard error of the estimate) nmol/mL for esketamine and 0.62 ± 0.19 nmol/mL for racemic ketamine, reflecting the 52 ± 27% lower analgesic potency of R-ketamine versus S-ketamine. Modeling showed that SNP had no effect on S-ketamine potency but abolished the R-ketamine analgesic effect. Similar observations were made for S- and R-norketamine. Since SNP had no effect on S-ketamine analgesia, we conclude that SNP interacts on R-ketamine nociceptive pathways, possibly similar to its effects on R-ketamine activated dissociation pathways.

Dose-exposure-efficacy response of intravenous immunoglobulin G 10% in multifocal motor neuropathy.

Multifocal motor neuropathy is a rare chronic immune-mediated neuropathy with impaired grip strength representing a common symptom. While intravenous immunoglobulin G is an effective treatment for the disease, significant variation in treatment response has been observed but not well understood. This analysis characterized dose-exposure-response relationships in multifocal motor neuropathy, using grip strength as a clinical efficacy measure. Serum immunoglobulin G trough concentrations and grip strength data for the more affected hand from a Phase 3, randomized, double-blind, placebo-controlled, crossover trial of intravenous immunoglobulin 10% in 44 patients with multifocal motor neuropathy (NCT00666263) were used to develop a population pharmacokinetic-pharmacodynamic model. The model adequately described the observed pharmacokinetic and pharmacodynamic data and relationships between intravenous immunoglobulin 10% dose, serum immunoglobulin G trough levels, grip strength, and inter-patient variabilities in multifocal motor neuropathy. Model-based simulations for various dosing regimens (0.4-2.0 g/kg every 2-4 weeks) indicated that ≥1.6 g/kg/month would achieve clinically meaningful improvements in grip strength (≥4 kg) in ≥70% of patients. More frequent dosing at an equivalent monthly dose led to a more consistent response in grip strength. Furthermore, splitting the dose over multiple days for high doses (>1 g/kg) did not impact grip strength. These findings suggest that the majority of patients with multifocal motor neuropathy would respond rapidly to intravenous immunoglobulin 10% with a range of dosing regimens. Shorter dosing intervals may avoid the diminishing response seen with longer dosing intervals. Dose-splitting provided similar outcomes while offering flexibility and convenience.

A pharmacokinetic randomised interventional study to optimise dihydroartemisinin-piperaquine dosing for malaria preventive treatment in Malawian infants: A protocol for the OPTIMAL study

Background A newer malaria preventive treatment, dihydroartemisinin-piperaquine (DP), has been identified as an effective alternative to sulfadoxine-pyrimethamine, to which malaria parasites are increasingly becoming resistant. However, how best to dose DP to safely prevent malaria in infants when aligned with routine health facility visits remains unresolved. As infants are usually excluded from participating in early dose optimisation clinical trials, the present study seeks to shift the paradigm and develop optimised DP dosing strategies for malaria preventive treatment in infants. Methods A randomised, single-blind, placebo-controlled, two-arm, interventional study will be conducted in southern Malawi. At 10 weeks (2.5 months) of age, 220 eligible infants will be randomised to receive DP (intervention group, n=110) or placebo (control group, n=110) with routine vaccines. They will be followed until 12 months of age and receive three further DP or placebo treatment courses at 14 weeks, six- and nine months. Infants in the intervention group will contribute capillary samples for piperaquine concentrations pre-dose and at three-, seven-, 14- and 28-days post-DP dosing as well as capillary samples pre-dose and on day 28 post-DP to quantify malaria parasitaemia using microscopy and quantitative PCR. In the control group, infants will contribute capillary blood samples for malaria parasitaemia at the same time points as the intervention group. Malaria incidence and adverse events will be compared between the two groups. Population pharmacokinetic-pharmacodynamic modelling techniques will be applied to derive feasible, optimised, efficacious, and safe DP dosing strategies for malaria preventive treatment in infancy. Conclusions The findings will provide the much-needed evidence to inform DP dosing for malaria preventive treatment in infants when administered with routine health facility visits. Additionally, they will help inform optimal DP dosing for malaria treatment in infants. The trial was registered with the Pan African Clinical Trials Registry; (PACTR202211575727659) on 8 November 2022. Protocol version 3.1, dated 29 September 2022.

Factors Influencing Quetiapine Pharmacokinetic Variability: A Review of Population Pharmacokinetics

Abstract: Atypical antipsychotic quetiapine (QTP) exhibits high pharmacokinetic variability and population pharmacokinetic (PopPK) analysis is one of the approaches used to characterize factors influencing QTP pharmacokinetic variability. Though QTP is not regarded as a narrow therapeutic index drug, knowledge of this area is of importance. Thus, this review was conducted to summarize significant predictors for QTP pharmacokinetic variability identified using a PopPK analysis and to explore any knowledge gaps to be investigated. PubMed, Scopus, and CINAHL Complete databases were searched for eligible studies, and 75 articles were identified. Of these, only five studies were included as they were conducted using a nonlinear mixed-effects approach. This review found that only limited predictors for QTP pharmacokinetics were identified, with body weight being a predictor for the volume of distribution and age and γ-glutamyl transpeptidase being predictors for QTP clearance. None of the studies included elderly patients aged >65 years, and thus factors associated with aging were not investigated. Also, most of the participants in the PopPK analyses were from clinical trials which might not reflect real-world patients e.g., the impacts of polypharmacy may not be available. Moreover, while a population pharmacokinetic-pharmacodynamic model explaining QTP exposure and clinical response using the scores of the Brief Psychiatric Rating Scale is available, knowledge relevant to the relationship between exposure and QTP side effects has not been explored. Based on this limited information, future PopPK research encompassing a wide range of patient characteristics is required.

Leishmania blood parasite dynamics during and after treatment of visceral leishmaniasis in Eastern Africa: A pharmacokinetic-pharmacodynamic model.

With the current treatment options for visceral leishmaniasis (VL), recrudescence of the parasite is seen in a proportion of patients. Understanding parasite dynamics is crucial to improving treatment efficacy and predicting patient relapse in cases of VL. This study aimed to characterize the kinetics of circulating Leishmania parasites in the blood, during and after different antileishmanial therapies, and to find predictors for clinical relapse of disease. Data from three clinical trials, in which Eastern African VL patients received various antileishmanial regimens, were combined in this study. Leishmania kinetoplast DNA was quantified in whole blood with real-time quantitative PCR (qPCR) before, during, and up to six months after treatment. An integrated population pharmacokinetic-pharmacodynamic model was developed using non-linear mixed effects modelling. Parasite proliferation was best described by an exponential growth model, with an in vivo parasite doubling time of 7.8 days (RSE 12%). Parasite killing by fexinidazole, liposomal amphotericin B, sodium stibogluconate, and miltefosine was best described by linear models directly relating drug concentrations to the parasite elimination rate. After treatment, parasite growth was assumed to be suppressed by the host immune system, described by an Emax model driven by the time after treatment. No predictors for the high variability in onset and magnitude of the immune response could be identified. Model-based individual predictions of blood parasite load on Day 28 and Day 56 after start of treatment were predictive for clinical relapse of disease. This semi-mechanistic pharmacokinetic-pharmacodynamic model adequately captured the blood parasite dynamics during and after treatment, and revealed that high blood parasite loads on Day 28 and Day 56 after start of treatment are an early indication for VL relapse, which could be a useful biomarker to assess treatment efficacy of a treatment regimen in a clinical trial setting.

Population pharmacokinetic-pharmacodynamic modeling of serum biomarkers as predictors of tumor dynamics following lenvatinib treatment in patients with radioiodine-refractory differentiated thyroid cancer (RR-DTC).

Lenvatinib is a receptor tyrosine kinase (RTK) inhibitor targeting vascular endothelial growth factor (VEGF) receptors 1-3, fibroblast growth factor (FGF) receptors 1-4, platelet-derived growth factor receptor-α (PDGFRα), KIT, and RET that have been implicated in pathogenic angiogenesis, tumor growth, and cancer. The primary objective of this work was to evaluate, by establishing quantitative relationships, whether lenvatinib exposure and longitudinal serum biomarker data (VEGF, Ang-2, Tie-2, and FGF-23) are predictors for change in longitudinal tumor size which was assessed based on data from 558 patients with radioiodine-refractory differentiated thyroid cancer (RR-DTC) receiving either lenvatinib or placebo treatment. Lenvatinib PK was best described by a 3-compartment model with simultaneous first- and zero-order absorption and linear elimination from the central compartment with significant covariates (body weight, albumin <30 g/dL, ALP>ULN, RR-DTC, RCC, HCC subjects, and concomitant CYP3A inhibitors). Except for body weight, none of the covariates have any clinically meaningful effect on exposure to lenvatinib. Longitudinal biomarker measurements over time were reasonably well defined by a PK/PD model with common EC50, Emax, and a slope for disease progression for all biomarkers. Longitudinal tumor measurements over time were reasonably well defined by a tumor growth inhibition Emax model, which in addition to lenvatinib exposure, included model-predicted relative changes from baseline over time for Tie-2 and Ang-2 as having significant association with tumor response. The developed PK/PD models pave the way for dose optimization and potential prediction of clinical response.

Mathematical modelling of clonal reduction therapeutic strategies in acute myeloid leukemia

Over the years, the overall survival of older patients diagnosed with acute myeloid leukemia (AML) has not significantly increased. Although standard cytotoxic therapies that rapidly eliminate dividing myeloblasts are used to induce remission, relapse can occur due to surviving therapy-resistant leukemic stem cells (LSCs). Hence, anti-LSC strategies have become a key target to cure AML. We have recently shown that previously approved cardiac glycosides and glucocorticoids target LSC-enriched CD34+ cells in the primary human AML 8227 model with more efficacy than normal hematopoietic stem cells (HSCs). To translate these in vitro findings into humans, we developed a mathematical model of stem cell dynamics that describes the stochastic evolution of LSCs in AML post-standard-of-care. To this, we integrated population pharmacokinetic-pharmacodynamic (PKPD) models to investigate the clonal reduction potential of several promising candidate drugs in comparison to cytarabine, which is commonly used in high doses for consolidation therapy in AML patients. Our results suggest that cardiac glycosides (proscillaridin A, digoxin and ouabain) and glucocorticoids (budesonide and mometasone) reduce the expansion of LSCs through a decrease in their viability. While our model predicts that effective doses of cardiac glycosides are potentially too toxic to use in patients, simulations show the possibility of mometasone to prevent relapse through the glucocorticoid’s ability to drastically reduce LSC population size. This work therefore highlights the prospect of these treatments for anti-LSC strategies and underlines the use of quantitative approaches to preclinical drug translation in AML.

New light on ω-3 polyunsaturated fatty acids and diabetes debate: a population pharmacokinetic-pharmacodynamic modelling and intake threshold study

Objectiveω-3 polyunsaturated fatty acids (PUFA) are a key modifiable factor in the intervention of type 2 diabetes, yet recommendations for dietary consumption of ω-3 PUFA in type 2 diabetes remain ambiguous and controversial. Here, we revisit the subject in the light of population pharmacokinetic-pharmacodynamic (PPK-PD) modeling and propose a threshold for intake.Research design and methodsPlasma levels of ω-3 PUFA and glycosylated hemoglobin (HbA1c) were measured as pharmacokinetic and pharmacodynamic indicator, respectively. The nonlinear mixed effect analysis was used to construct a PPK-PD model for ω-3 PUFA and to quantify the effects of FADS gene polymorphism, age, liver and kidney function, and other covariables.ResultsData from 161 patients with type 2 diabetes in the community were modeled in a two-compartment model with primary elimination, and HDL was a statistically significant covariate. The simulation results showed that HbA1c showed a dose-dependent decrease of ω-3 PUFA plasma level. A daily intake of ω-3 PUFA at 0.4 g was sufficient to achieve an HbA1c level of 7% in more than 95% of patients.ConclusionsPPK/PD modeling was proposed as a multilevel analytical framework to quantitatively investigate finer aspects of the complex relationship between ω-3 PUFA and type 2 diabetes on genetic and non-genetic influence factors. The results support a beneficial role for ω-3 PUFA in type 2 diabetes and suggested the intake threshold. This new approach may provide insights into the interaction of the two and an understanding of the context in which changes occur.

The development and evaluation of dose-prediction tools for allopurinol therapy (Easy-Allo tools).

Dose escalation at the initiation of allopurinol therapy can be protracted and resource intensive. Tools to predict the allopurinol doses required to achieve target serum urate concentrations would facilitate the implementation of more efficient dose-escalation strategies. The aim of this research was to develop and externally evaluate allopurinol dosing tools, one for use when the pre-urate-lowering therapy serum urate is known (Easy-Allo1) and one for when it is not known (Easy-Allo2). A revised population pharmacokinetic-pharmacodynamic model was developed using data from 653 people with gout. Maintenance doses to achieve the serum urate target of <0.36 mmol L-1 in >80% of individuals were simulated and evaluated against external data. The predicted and observed allopurinol doses were compared using the mean prediction error (MPE) and root mean square error (RMSE). The proportion of Easy-Allo predicted doses within 100 mg of the observed was quantified. Allopurinol doses were predicted by total body weight, baseline urate, ethnicity and creatinine clearance. Easy-Allo1 produced unbiased and suitably precise dose predictions (MPE 2mg day-1 95% confidence interval [CI] -13-17, RMSE 91%, 90% within 100 mg of the observed dose). Easy-Allo2 was positively biased by about 70 mg day-1 and slightly less precise (MPE 70 mg day-1 95% CI 52-88, RMSE 131%, 71% within 100 mg of the observed dose). The Easy-Allo tools provide a guide to the allopurinol maintenance dose requirement to achieve the serum urate target of <0.36 mmol L-1 and will aid in the development of novel dose-escalation strategies for allopurinol therapy.

Is Gender an Important Factor in the Precision Medicine Approach to Levocetirizine?

Currently, there is insufficient information on the variability in levocetirizine pharmacometrics among individuals, a crucial aspect for establishing its clinical use. The gender differences in pharmacokinetics and the extent of variation in pharmacodynamics have not been definitively identified. The primary goal of this study was to investigate gender differences in levocetirizine pharmacokinetics and quantitatively predict and compare how these gender-related pharmacokinetic differences impact pharmacodynamics, using population pharmacokinetic-pharmacodynamic modeling. Bioequivalence results for levocetirizine (only from the control formulation) were obtained from both healthy Korean men and women. Physiological and biochemical parameters for each individual were utilized as pharmacokinetic comparison and modeling data between genders. Pharmacodynamic modeling was performed using reported data on antihistamine responses following levocetirizine exposure. Gender, weight, body surface area, peripheral distribution volume, albumin, central-peripheral inter-compartmental clearance, and the fifth sequential absorption rate constant were explored as effective covariates. A comparison of the model simulation results showed a higher maximum concentration and faster plasma loss in females than in males, resulting in a faster recovery to baseline of the antihistamine effect; however, the absolute differences between genders in the mean values were not large within 10 ng/mL (for plasma concentrations) or % (wheal and flare size changes). Regarding the pharmacokinetics and pharmacodynamics of levocetirizine, the gender effect may not be significant when applying the usual dosage (5 mg/day). This study will be useful for bridging the knowledge gap in scientific precision medicine by introducing previously unconfirmed information regarding gender differences in levocetirizine pharmacometrics.

Using real-world data to inform dosing strategies of rituximab for pediatric patients with frequent-relapsing or steroid-dependent nephrotic syndrome: a prospective pharmacokinetic-pharmacodynamic study.

Objectives: Rituximab is frequently used off-label for the treatment of frequent-relapsing/steroid-dependent nephrotic syndrome (FRNS/SDNS). However, the optimal dosing schedules remain undetermined. The objective of this study was to establish a population pharmacokinetic-pharmacodynamic (PK-PD) model in pediatric patients with FRNS/SDNS, and to investigate dosing regimens that provide adequate suppression of B lymphocytes. Methods: A prospective, open-label, single-center study was conducted in Nephrology Department at Children's Hospital of Fudan University, and a two-compartment PK model of rituximab in pediatric FRNS/SDNS has been developed previously by our group. CD19+ lymphocyte count profiles were obtained from these patients. The presence of anti-rituximab antibodies was assessed prior to medication in children who had previously received rituximab or during follow-up at the last sampling point for PK analysis. PK-PD analyses were performed to describe the changes of CD19+ lymphocytes, with rituximab assumed to increase their death rate. Monte Carlo simulation was conducted to evaluate different dosing regimens. Results: In total, 102 measurements of CD19+ lymphocyte counts were available for PK-PD analysis. No detectable levels of anti-rituximab antibodies were observed during the PK follow-up period. A turnover model with saturable stimulatory action of rituximab on the removal of lymphocytes best characterized the relationship between rituximab concentration and CD19+ lymphocyte counts, where the Emax and EC50 were estimated to be 99.6*106/L and 5.87μg/mL, respectively. Simulations indicated that a single infusion of 750mg/m2 and 2 infusions of 375mg/m2 both yielded a 10-week suppression of CD19+ lymphocytes. Conclusion: This study represents a first attempt to quantitatively describe the PK-PD relationship of rituximab in pediatric patients with FRNS/SDNS, and provide a potential pathway for future precision dosing strategy for rituximab therapy. Further clinical studies are warranted to evaluate the efficacy and safety of different dosing schemes.

Population Pharmacokinetic-Pharmacodynamic Modeling of Neutrophil and Platelet Count for Lower-Intensity Therapy of CPX-351 Combined with Venetoclax in Acute Myeloid Leukemia

Background: Improved survival outcomes and a comparable safety profile with CPX-351 (a dual-drug liposomal encapsulation of daunorubicin and cytarabine in a synergistic 1:5 molar ratio) vs conventional 7+3 chemotherapy in a pivotal phase 3 trial led to the approval of CPX-351 for newly diagnosed therapy-related acute myeloid leukemia (tAML) or AML with myelodysplasia-related changes (AML MRC) in adults and pediatric (aged ≥1 year) patients in the United States and in adults in Europe. The approved induction dose is 1-2 cycles of CPX-351 100 units/m 2 (daunorubicin 44 mg/m 2 and cytarabine 100 mg/m 2) as a 90-minute intravenous infusion on days 1, 3, and 5 (days 1 and 3 for second induction). In this study, we developed a population pharmacokinetic-pharmacodynamic (PK-PD) model of absolute neutrophil count (ANC) and platelet count to determine the optimal dose for a lower-intensity therapy (LiT) of CPX-351 combined with venetoclax. Methods: The CPX-351 population PK-PD analysis was based on data collected from 3 clinical studies (phase 1 study [NCT00389428; J Clin Oncol 2011; 29(8):979-85];phase 2 study [NCT02238925; Cancer Chemother Pharmacol 2019; 84(1):163-73];phase 3 study [NCT01696084; J Clin Oncol 2018; 36(26):2684-92]) in patients with newly diagnosed high-risk/secondary AML. Venetoclax population PK-PD analysis was based on data from a comprehensive safety analysis of venetoclax monotherapy ( Clin Cancer Res 2018; 24(18):4371-9). ANC and platelet data used in the analysis were from 94 patients for CPX-351 and 272 patients for venetoclax. Generalized semi-mechanistic population PK-PD analyses were performed using non-linear mixed-effects modeling to characterize the time-course of ANC and platelet count following administration of CPX-351 and venetoclax monotherapy. The majority of system parameters were shared across the 2 therapies and a limited number of treatment-specific parameters were estimated by fitting CPX-351 and venetoclax PK time-courses simultaneously. An R Shiny app was developed for the simulation of different dosing regimens with favorable safety profiles to guide optimal dose selection for LiT of CPX-351 combined with venetoclax, assuming the 2 drug effects on ANC and platelet count were additive. Results: The generalized PK-PD model included multi-compartment transit models for ANC and platelet count with separate semi-mechanistic compartments to represent proliferating, maturing, and circulating cells. The model fitted CPX-351 and venetoclax PK-PD data with shared system parameters except for the drug-specific parameter IC 50 (concentration producing 50% of maximal inhibition of cell proliferation). For ANC, the estimated IC 50 values for CPX-351 and venetoclax were 295 μM (relative standard error [RSE]: 8.96%) and 120 μM (RSE: 18.0%), respectively; the shared mean transit time (MTT) estimate for ANC maturation was 66.6 h (RSE: 4.95%), and shared maximal inhibition of cell proliferation (I max) was fixed to 1. For platelet count, the estimated IC 50 values for CPX-351 and venetoclax were 0.0109 μM (RSE: 50.2%) and 16.5 μM (RSE: 8.73%); the shared MTT estimate for platelet maturation was 87.2 h (RSE: 9.01%) and shared I max estimate was 0.313 (RSE: 11.2%). Simulations with different doses of CPX-351 in combination with 400 mg venetoclax were performed. The simulated mean nadir ANC were 2.1, 1.5, 1.1, and 0.6 × 10 9/L for CPX-351 doses of 20, 40, 60, and 100 units/m 2, respectively, when administered on days 1 and 3 and combined with venetoclax 400 mg daily for 28 days. The simulated mean nadir platelet counts were 31.0, 30.3, 30.0, and 29.7 × 10 9/L for CPX-351 doses of 20, 40, 60, and 100 units/m 2, respectively, when administered on days 1 and 3 and combined with venetoclax 400 mg daily for 28 days. Model-based simulated outcomes of ANC and platelet count were consistent with the observed trial data. C onclusions: A generalized semi-mechanistic population PK-PD model was developed to adequately describe the time-courses of ANC and platelet counts in patients with AML following the administration of CPX-351 and venetoclax monotherapy and was used to infer the impact on ANC and platelet count when given in combination. Model-based simulations supported initial dose selection in the phase 1b study of combined CPX-351 LiT plus venetoclax as first-line treatment for patients with AML who are unfit for intensive chemotherapy (NCT04038437).

Exploring Differences in Pharmacometrics of Rabeprazole between Genders via Population Pharmacokinetic-Pharmacodynamic Modeling.

Rabeprazole is a proton pump inhibitor that inhibits gastric acid production and increases gastric pH; it is widely used clinically as a treatment option for gastritis and gastric ulcers. However, information on the inter-individual variability of rabeprazole pharmacometrics, which is a key element in establishing its scientific clinical use, is still lacking. Particularly, the differences in pharmacokinetics between genders and the degree of variation in pharmacodynamics have not been clearly identified. Thus, the main purpose of this study was to explore any differences in rabeprazole pharmacokinetics between genders and to quantitatively predict and compare the effects of any differences in pharmacokinetics between genders on known pharmacodynamics using population pharmacokinetic-pharmacodynamic modeling. To compare pharmacokinetics and modeling data between genders, bioequivalence results were used simultaneously on healthy Korean men and women using the physiological and biochemical parameters derived from each individual. Pharmacodynamic modeling was performed based on the data of previously reported gastric pH changes in response to rabeprazole plasma concentrations, which was co-linked to the central compartmental bioavailable concentration in the population pharmacokinetic model. There was no significant difference in the level of rabeprazole exposure and elimination of plasma between genders following oral administration of 10 mg enteric-coated rabeprazole tablets; however, there was a clear delay in absorption in women compared to men. Additionally, a comparison of pharmacokinetic parameters normalized to body weight between genders showed that the maximum plasma concentrations were significantly higher in women than in men, again suggesting gender differences in rabeprazole absorption. The population pharmacokinetic profiles for rabeprazole were described using a three-sequential multi-absorption with lag time (Tlag) two-compartment model, whereas body surface area and gender were explored as effective covariates for absorption rate constant and Tlag, respectively. The effect of increased gastric pH due to plasma exposure to rabeprazole was explained using the Sigmoid Emax model, with the baseline as a direct response. The significantly longer rabeprazole Tlag in females delayed the onset of an effect by an average of 1.58 times (2.02-3.20 h), yet the overall and maximum effects did not cause a significant difference within 15%. In the relative comparison of the overall efficacy of rabeprazole enteric-coated tablet administration between genders, it was predicted based on the model that males would have higher efficacy. This study will be very useful in broadening the perspective of interpreting drug diversity between individuals and narrowing the gap in knowledge related to scientific precision medicine by presenting new information on gender differences in rabeprazole pharmacometrics that had not been previously identified.

Population pharmacokinetic-pharmacodynamic modeling of clopidogrel for dose regimen optimization based on CYP2C19 phenotypes: A proof of concept study.

Clopidogrel is an antiplatelet drug used to reduce the risk of acute coronary syndrome and stroke. It is converted by CYP2C19 to its active metabolite; therefore, poor metabolizers (PMs) of CYP2C19 exhibit diminished antiplatelet effects. Herein, we conducted a proof-of-concept study for using population pharmacokinetic-pharmacodynamic (PK-PD) modeling to recommend a personalized clopidogrel dosing regimen for individuals with varying CYP2C19 phenotypes and baseline P2Y12 reaction unit (PRU) levels. Data from a prospective phase I clinical trial involving 36 healthy male participants were used to develop the population PK-PD model predicting the concentrations of clopidogrel, clopidogrel H4, and clopidogrel carboxylic acid, and linking clopidogrel H4 concentrations to changes in PRU levels. A two-compartment model effectively described the PKs of both clopidogrel and clopidogrel carboxylic acid, and a one-compartment model of those of clopidogrel H4. The CYP2C19 phenotype was identified as a significant covariate influencing the metabolic conversion of the parent drug to its metabolites. A PD submodel of clopidogrel H4 that stimulated the fractional turnover rate of PRU levels showed the best performance. Monte Carlo simulations suggested that PMs require three to four times higher doses than extensive metabolizers to reach the target PRU level. Individuals within the top 20th percentile of baseline PRU levels were shown to require 2.5-3 times higher doses than those in the bottom 20th percentile. We successfully developed a population PK-PD model for clopidogrel considering the impact of CYP2C19 phenotypes and baseline PRU levels. Further studies are necessary to confirm actual dosing recommendations for clopidogrel.

2311P Population pharmacokinetic-pharmacodynamic modeling to inform optimal dosing strategies for GI-101, a novel fusion protein, targeting IL2βγR and CTLA4

2311P Population pharmacokinetic-pharmacodynamic modeling to inform optimal dosing strategies for GI-101, a novel fusion protein, targeting IL2βγR and CTLA4

Methotrexate Polyglutamates Exposure - Response Modeling in a Large Cohort of Rheumatoid Arthritis Patients Starting Methotrexate.

Methotrexate polyglutamates (MTX-PG) concentrations in red blood cells (RBCs) have been suggested as a biomarker of response in patients with rheumatoid arthritis (RA) receiving low-dose MTX therapy. We investigated the association and interpatient variability between RBC-MTX-PG3-5 -exposure and response in patients with RA starting MTX. Data of three prospective cohorts were available. The relationship between exposure and Disease Activity Score in 28 joints (DAS28) was analyzed using a population pharmacokinetic-pharmacodynamic model. Relevant covariates were tested using full covariate modeling and backward elimination. From 395 patients, 3,401 MTX-PG concentrations and 1,337 DAS28 measurements were available between 0 and 300 days after MTX treatment onset. The developed model adequately described the time course of MTX-PG3-5 and DAS28. The median MTX-PG3-5 level at month 1 was 30.9 nmol/L (interquartile range (IQR): 23.6-43.7; n = 41) and at month 3: 69.3 nmol/L (IQR: 17.9-41.2; n = 351). Clearance of MTX-PG3-5 from RBCs was 28% lower (95% confidence interval (CI): 23.6-32.8%) in a woman and 10% lower (95% CI: 7.7-12.4%) in a 65-year-old compared with a 35-year-old patient. MTX-PG3-5 concentrations associated with DAS28: half-maximal effective concentration (EC50 ) was 9.14 nmol/L (95% CI: 4.2 nmol/L-14.1 nmol/L). EF at 80% (EC80 ) above 47 nmol/L was regarded as the optimal response. Independent of the MTX-PG 3-5 - response association, co-administration of disease-modifying antirheumatic drugs and corticosteroids improved response (additive effect on maximum effect (Emax )), whereas smoking, high body mass index and low albumin decreased Emax . In patients with RA starting MTX, RBC-MTX-PG3-5 was associated with clinical response. A dose increase is suggested when MTX-PG3-5 at month 1 is below 9.15 nmol/L, continued with the same dose when the concentration is above 47 nmol/L, and consider other treatment options above 78 nmol/L from 3months onwards.

Population pharmacokinetic-pharmacodynamic model of nirogacestat effects on B-cell maturation antigen in healthy subjects.

8027 Background: B-cell maturation antigen (BCMA) is expressed on the membrane of normal plasma cells (PCs) and multiple myeloma (MM) cells and is the target of several investigational agents and approved products for MM treatment. Mechanistically, BCMA is cleaved from the cell surface by the enzyme gamma secretase, resulting in reduced levels of membrane-bound BCMA (mbBCMA) and release of soluble BCMA. Gamma secretase inhibitors (GSIs) have been shown to increase BCMA receptor density and potentiate the activity of several BCMA-targeted therapies in both clinical trials and in vitro studies. Although the effect of GSIs on BCMA receptor density has been well characterized in vitro, their effect on BCMA dynamics is incompletely understood in humans. This study sought to evaluate the pharmacodynamics (PD) of the GSI nirogacestat on mbBCMA in healthy subjects. Methods: A total of 23 healthy subjects were administered single doses of 50, 150, or 300 mg nirogacestat or multiple doses of 100 mg every 12 hours. Density of mbBCMA on PCs isolated from whole blood and bone marrow was measured by flow cytometry at baseline and 1, 2, 4, 8, 24, and 48 hours after single doses and 24 and 48 hours after the first dose in subjects receiving multiple doses. Nirogacestat serum concentrations were also measured at multiple time points to fully characterize the pharmacokinetics (PK). A population PK-PD model of the effect of nirogacestat on BCMA dynamics in humans was developed and used to simulate various administration dose levels and schedules to find an optimal nirogacestat regimen that continually increases mbBCMA over the dosing interval. Results: A 2-compartment model with linear absorption and elimination described serum nirogacestat in healthy subjects. A sequential PK-PD model that included an effect compartment and drug effect on mbBCMA cleavage rate described mbBCMA values. Simulations indicated that a 100-mg dose administered twice daily sufficiently maintained elevated and clinically relevant mbBCMA levels. Exposure and efficacy benchmarks from in vitro experiments and published clinical studies were used to contextualize simulations and hypothesize various dose and schedule scenarios to minimize dose of the BCMA-targeting therapy while maximizing BCMA potentiation. Conclusions: These results suggest that nirogacestat dose and schedule may have direct implications on the potentiation of BCMA-targeting therapies. Simulations suggest that the 100-mg, twice-daily nirogacestat dose results in elevated mbBCMA levels that could even allow for use of lower doses of BCMA-targeted agents while maintaining or improving efficacy. Study results along with this PK-PD model will be used to support nirogacestat dose and schedule in ongoing and planned studies in patients with MM.

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 pharmacokinetic-pharmacodynamic model of dexmedetomidine in elderly patients undergoing sedation after abdominal aortic surgery

Background. Dexmedetomidine (DEX) is a widely used sedative agent for treating post-surgery patients. It also acts on hemodynamic parameters like heart rate or cardiac output. This study aimed to develop a pharmacokinetic-pharmacodynamic (PK/PD) model of DEX using bispectral index (BIS) and cardiac output (CO) as a response. Methodology and results. 21 mechanically ventilated elderly cardiac patients undergoing abdominal aortic surgery were enrolled in the study. DEX was given to maintain moderate or deep sedation. Genotypes of ADR2A*55 were identified using real-time PCR-HRM. Data were analyzed using nonlinear mixed-effect modelling. A two-compartment model described DEX pharmacokinetics. The sigmoid Emax and linear models were used to describe BIS and CO measurements. The typical value of EC50 for DEX effects on BIS was 3.62 ng/ml, and the slope between CO and DEX concentrations was 0.819 (L/min)/(ng/ml). We were unable to show the effects of considered covariates on DEX pharmacodynamics. Conclusions. WE proposed the PK/PD model of DEX to understand better the BIS and CO changes observed after surgery. The measured CI values were in the reference range showing that the used doses of DEX ensured stable cardiac function in the studied patients.

Predicting treatment response to vancomycin using bacterial DNA load as a pharmacodynamic marker in premature and very low birth weight neonates: A population PKPD study.

Background: While positive blood cultures are the gold standard for late-onset sepsis (LOS) diagnosis in premature and very low birth weight (VLBW) newborns, these results can take days, and early markers of possible treatment efficacy are lacking. The objective of the present study was to investigate whether the response to vancomycin could be quantified using bacterial DNA loads (BDLs) determined by real-time quantitative polymerase chain reaction (RT-qPCR). Methods: VLBW and premature neonates with suspected LOS were included in a prospective observational study. Serial blood samples were collected to measure BDL and vancomycin concentrations. BDLs were measured with RT-qPCR, whereas vancomycin concentrations were measured by LC-MS/MS. Population pharmacokinetic-pharmacodynamic modeling was performed with NONMEM. Results: Twenty-eight patients with LOS treated with vancomycin were included. A one-compartment model with post-menstrual age (PMA) and weight as covariates was used to describe the time PK profile of vancomycin concentrations. In 16 of these patients, time profiles of BDL could be described with a pharmacodynamic turnover model. The relationship between vancomycin concentration and first-order BDL elimination was described with a linear-effect model. Slope S increased with increasing PMA. In 12 patients, no decrease in BDL over time was observed, which corresponded with clinical non-response. Discussion: BDLs determined through RT-qPCR were adequately described with the developed population PKPD model, and treatment response to vancomycin using BDL in LOS can be assessed as early as 8h after treatment initiation.