Population Pharmacokinetic and Pharmacodynamic of Atorvastatin in Chinese Lung Transplant Recipients.

Atorvastatin is a widely prescribed cholesterol-lowering drug that inhibits 3-hydroxy-3-methylglutaryl coenzyme A reductase to reduce lipid levels. This study aimed to establish a population pharmacokinetic and pharmacodynamic model for atorvastatin in Chinese lung transplant recipients (LTRs), particularly those using voriconazole (VOR) and with different genotypes. It evaluated precise dosing regimens and analyzed the correlation between atorvastatin exposure and clinical outcomes. A nonlinear mixed-effects model was used for the population pharmacokinetic/pharmacodynamic (PK/PD) analysis. A one-compartment population PK model was developed, incorporating VOR, SLCO2A1 rs76906503, and SLC22A8 rs2187383 to assess apparent clearance and volume of distribution. LDL-C was modeled as a biomarker to evaluate atorvastatin efficacy. A Monte Carlo simulation was conducted to assess various dosing schemes and the effects of different covariates on achieving the target LDL concentration. The correlation between atorvastatin exposure and clinical outcomes was also evaluated. Results indicated that the average probability of target attainment for optimal dosing regimens across various covariate results exceeded 45.8%. Dosages of 10, 20, and 40 mg were deemed suitable for LTRs. A lower dose was recommended for LTRs taking VOR or with mutant-type genotypes to avoid overexposure and adverse reactions. The population PK/PD model offers valuable guidance for evaluating atorvastatin dosing regimens in clinical settings, particularly for LTRs using VOR and those with different genotypes.

Toward Model-Informed Precision Dosing for Remimazolam: A Population Pharmacokinetic-Pharmacodynamic Analysis.

Remimazolam, widely used for procedural sedation and general anesthesia, is a new ultra short-acting benzodiazepine for intravenous sedation and anesthesia. We aim to characterize the pharmacokinetics/pharmacodynamics (PK/PD) of remimazolam and its metabolite CNS 7054 in healthy Chinese volunteers using population analysis and suggest an optimal dosing regimen for sedation therapy. Data were collected from a single-center, placebo-controlled, randomized, and dose-escalation clinical pharmacology study. Forty-six healthy volunteers received a single infusion dose of remimazolam, while nine healthy subjects received a continuous infusion of remimazolam. A population PK/PD model was established and RxODE and Shiny in R were used to design the remimazolam dosing regimens. A three-compartment model best described the PK of remimazolam and a two-compartment model with one transit compartment was adopted for CNS 7054. The relationship between exposure and the bispectral index was best described using an effect compartment model with an inhibitory sigmoid model. Additionally, a web-based dashboard was developed to provide individualized dosing regimens, complemented by a graphical illustration of the PK/PD profiles of the proposed dosing regimen. The established population PK/PD model characterized the dose-exposure-response relationship of remimazolam well, which could be applied to optimize individual dosing regimens.

Application of machine learning techniques in population pharmacokinetics/pharmacodynamics modeling.

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

The history and future of population pharmacokinetic analysis in drug development.

1. The analysis of pharmacokinetic data has been in a constant state of evolution since the introduction of the term pharmacokinetics. Early work focused on mechanistic understanding of the absorption, distribution, metabolism and excretion of drug products.2. The introduction of non-linear mixed effects models to perform population pharmacokinetic analysis initiated a paradigm shift. The application of these models represented a major shift in evaluating variability in pharmacokinetic parameters across a population of subjects.3. While technological advancements in computing power have feuled the growth of population pharmacokinetics in drug development efforts, there remain many challenges in reducing the time required to incorporate these learnings into a model-informed development process. These challenges exist because of expanding datasets, increased number of diagnostics, and more complex mathematical models.4. New machine learning tools may be potential solutions for these challenges. These new methodologies include genetic algorithms for model selection, machine learning algorithms for covariate selection, and deep learning models for pharmacokinetic and pharmacodynamic data. These new methods promise the potential for less bias, faster analysis times, and the ability to integrate more data.5. While questions remain regarding the ability of these models to extrapolate accurately, continued research in this area is expected to address these questions.

Model Informed Development of SIM0295 in Patients with Gout and Hyperuricemia and Healthy Volunteers Using a Population Pharmacokinetics/ Pharmacodynamics Approach.

BACKGROUND: SIM0295, a novel inhibitor of human uric acid transporter 1 (hURAT1), is used to treat patients with gout and hyperuricemia. This study aimed to develop population pharmacokinetics and pharmacodynamics (popPK/PD) models of SIM0295 and explore potential covariates to inform clinical drug development. RESEARCH DESIGN AND METHODS: Data were obtained from four phase I studies conducted in healthy Korean and Chinese subjects and two phase II studies conducted in Korean patients with gout and hyperuricemia. The popPK/PD model of SIM0295 was developed using nonlinear mixed effects modeling. RESULTS: SIM0295 pharmacokinetics was described using a two-compartment model with the absorption of four transit compartments and first-order elimination. PK parameters were normalized to weight via allometric scaling. Food was identified as a factor significantly affecting the absorption rate, with no clinical relevance. The sigmoid Emax model with a semi-mechanism of inhibition of serum uric acid (sUA) reabsorption was used to describe the exposure-response relationship. Additionally, Monte Carlo simulations demonstrated that approimately 9 mg/day of SIM0295 for 7 days could achieve the maximum decrease in sUA. CONCLUSION: The established popPK/PD model characterized the dose-exposure-response relationship for SIM0295 in healthy subjects and patients with gout and hyperuricemia and could be used to inform the drug development.

Developmental Population Pharmacokinetics-Pharmacodynamics of Meropenem in Chinese Neonates and Young Infants: Dosing Recommendations for Late-Onset Sepsis.

The pharmacokinetic (PK) studies of meropenem in Chinese newborns with late-onset sepsis (LOS) are still lacking. Causative pathogens of LOS and their susceptibility patterns in China differ from the data abroad. We, therefore, conducted a developmental population pharmacokinetic−pharmacodynamic analysis in Chinese newborns with the goal to optimize meropenem dosing regimens for LOS therapy. An opportunistic sampling strategy was used to collect meropenem samples, followed by model building and validation. A Monte Carlo simulation was performed to show the probability of target attainment (PTA) for various dosages. The information from 78 newborns (postmenstrual age: 27.4−46.1 weeks) was compiled and had a good fit to a 1-compartment model that had first order elimination. The median (range) values of estimated weight−normalized volume of distribution (V)and clearance (CL) were 0.60 (0.51−0.69) L/kg and 0.16 (0.04−0.51) L/h/kg, respectively. Covariate analysis revealed that postnatal age (PNA), gestational age (GA) and current weight (CW) were the most important factors in describing meropenem PK. Simulation results showed for LOS with a minimal inhibitory concentration (MIC) of 8 mg/L, the doses of 30 mg/kg 3 times daily (TID) as a 1-h infusion for newborns with GA ≤ 37 weeks and 40 mg/kg TID as a 3-h infusion for those with GA > 37 weeks were optimal, with PTA of 71.71% and 75.08%, respectively. In conclusion, we proposed an evidence-based dosing regimen of meropenem for LOS in Chinese newborns by using the population pharmacokinetic−pharmacodynamic analysis, based on domestic common pathogens and their susceptibility patterns.

Rivaroxaban population pharmacokinetic and pharmacodynamic modeling in Iranian patients.

WHAT IS KNOWN AND OBJECTIVE: Although predictable pharmacokinetic and pharmacodynamic of rivaroxaban allow fixed dosing regimens without routine coagulation monitoring, there is still the necessity to monitor and predict the effects of rivaroxaban in specific conditions and different populations. The current study was designed and conducted to analyze the rivaroxaban population pharmacokinetics in Iranian patients and establish a pharmacokinetic/pharmacodynamic model to predict the relationship between rivaroxaban concentration and its anticoagulant activity. METHODS: A sequential nonlinear mixed effect pharmacokinetic/pharmacodynamic modeling method was used to establish the relation between rivaroxaban concentration and anti-factor Xa activity, prothrombin time, and activated partial thromboplastin time (aPTT) as pharmacodynamic biomarkers in a population of sixty-nine Iranian patients under treatment with oral rivaroxaban. Rivaroxaban plasma concentration was quantified by a validated high-performance liquid chromatography-tandem mass spectrometry. RESULTS AND DISCUSSION: The typical population values (inter-individual variability%) of the oral volume of distribution and clearance for a one-compartment model were 61.2 L (21%) and 3.68 L·h-1 (61%), respectively. Creatinine clearance and Child-Turcotte-Pugh score were found to affect the clearance. A direct link linear structural model best fitted the data for both prothrombin time and aPTT. The baseline estimates of aPTT and prothrombin time in the population were 35.0 (15%) and 12.6 (2%) seconds, respectively. The slope of the relationship between apTT, prothrombin time, and rivaroxaban concentration was 0.033 (28%) and 0.018 (54%) s·ml·ng-1 , respectively. The selected model for anti-factor Xa activity consisted of a direct link inhibitory Emax model with Hill coefficient. The maximum level of inhibition (Emax ) was 4 IU·ml-1 . The concentration of rivaroxaban producing 50% of the maximum inhibitory effect (EC50 ) was 180 (24%) ng·ml-1 , and Hill coefficient (γ) was 1.44 (108%). No covariates showed a statistically significant effect on PT and activated partial thromboplastin time prolonging properties and anti-factor Xa activity. WHAT IS NEW AND CONCLUSION: Our results confirmed that pharmacokinetic/pharmacodynamic models similar to those of the other studies describe the relationship between the rivaroxaban concentration and its anticoagulant effect in Iranian patients. However, considerable differences were observed in the parameters of the pharmacodynamics-pharmacokinetic models with the results of other reports that can explain the unpredictable effects of rivaroxaban in some patients.

Population pharmacokinetic and pharmacodynamic analysis of rivaroxaban in Chinese patients with non-valvular atrial fibrillation.

Rivaroxaban, a direct factor Xa inhibitor, is widely used for stroke prevention in patients with non-valvular atrial fibrillation (NVAF). The aim of this study was to conduct a population pharmacokinetic-pharmacodynamic (PK-PD) analysis of rivaroxaban in Chinese patients with NVAF to assess ethnic differences and provide model-based precision dosing. A total of 256 rivaroxaban plasma concentrations and 244 prothrombin time (PT) measurements were obtained from 195 Chinese NVAF patients from a prospective clinical trial. The population PK-PD model was developed using nonlinear mixed effects modeling (NONMEM) software. The PK of rivaroxaban was adequately described using a one-compartment model with first-order adsorption and elimination. Estimated glomerular filtration rate (eGFR) was identified as a major covariate for apparent clearance. No single nucleotide polymorphism was identified as a significant covariate. PT exhibited a linear relationship with rivaroxaban concentration. Total bilirubin (TBIL) and eGFR were identified as significant covariates for baseline PT. According to the Monte Carlo simulation, 15 mg for Chinese patients with eGFR ≥50 mL/min and normal liver function yielded an exposure comparable to 20 mg for Caucasian patients. Patients with moderately impaired renal function may require a lower dose of rivaroxaban to avoid overexposure. Moreover, there was an approximate 26% increase in PT levels in patients with TBIL of 34 µmol/L and eGFR of 30 mL/min, which could increase the risk of major bleeding. The established population PK-PD model could inform individualized dosing for Chinese NVAF patients who are administered rivaroxaban.

Population pharmacokinetics-pharmacodynamics of fondaparinux in dialysis-dependent chronic kidney disease patients undergoing chronic renal replacement therapy.

PURPOSE: Data on the anti-Xa efficacy of fondaparinux in dialysis-dependent chronic kidney disease (DD-CKD) patients are scarce. This study characterizes the pharmacokinetics (PK) and pharmacodynamics (PD) of fondaparinux in DD-CKD patients undergoing renal replacement therapy (RRT), to assess dosing strategies. METHODS: A retrospective, observational study was conducted using data on anti-Xa activity (112 samples) from 12 (3 male and 9 female) DD-CKD patients (median (IQR) age 71 years (63-88), weight 73 kg (59-98.5)). Eleven patients underwent high-flux or low-flux hemodialysis (HD) and one patient underwent peritoneal dialysis. Three patients were also treated with therapeutic plasma exchange (TPE). A non-linear mixed effects analysis was performed using NONMEM 7.3.0. RESULTS: The lab-specific slope of the relationship between fondaparinux concentration and anti-Xa levels was 1.18 IU/µg. In a one-compartment model, clearance (CL) and volume of distribution (Vd) were 0.05289 L/h and 5.55 L, respectively. High-flux HD was found to increase the CL of fondaparinux 2.26 times. TPE also considerably increased CL, but the fold-change could not be accurately estimated. Low-flux HD and peritoneal dialysis did not impact PK parameters. CONCLUSIONS: Model-based simulations showed that standard dosing (2.5 mg three times weekly before HD) results in a median anti-Xa activity of 0.55 IU/mL and 0.98 IU/mL, pre- and post-low-flux HD, respectively. In patients undergoing high-flux HD, these values are approximately 27% lower. Additional caution is warranted with TPE, as this treatment can reduce anti-Xa activity even further.

Population pharmacokinetic-pharmacodynamic model-based exploration of alternative ustekinumab dosage regimens for patients with Crohn's disease.

AIMS: In the UNITI endoscopy sub-study, only 17.4% of patients with Crohn's disease (CD) on ustekinumab achieved endoscopic response and 10.9% of patients achieved endoscopic remission at week (w)44. We aimed to evaluate the impact of alternative ustekinumab dosage regimens on endoscopic outcomes based on population pharmacokinetic-pharmacodynamic (popPK-PD) modelling and simulation analysis. METHODS: Real-world data were obtained from 83 patients with moderate-to-severe CD (95% biological-refractory) enrolled in a prospective cohort study receiving intravenous ustekinumab (~6 mg/kg) followed by every eight-week (q8w) subcutaneous maintenance therapy (90 mg). Three sequential models were developed: a two-compartment popPK model linking ustekinumab dose to ustekinumab exposure, an indirect response popPK-PD model describing the effect of ustekinumab exposure on fecal calprotectin (fCal), and a logistic regression outcome model linking fCal to endoscopic outcomes. RESULTS: Ustekinumab clearance increased with decreasing serum albumin and increasing bodyweight. fCal decreased with increasing ustekinumab exposure. The probability of endoscopic response at w24 increased from 10.0% to 17.9% with fCal at w8 decreasing from 1800 µg/g to 694 µg/g (EC50 ). The probability of endoscopic remission at w24 increased from 2.1% to 10.0% with fCal at w8 decreasing from 1800 µg/g to 214 µg/g (EC50 ). Simulation-based comparison of q8w and q4w maintenance dosing regimens predicted 16.7% and 22.2% endoscopic response rates, respectively. Endoscopic remission rates were estimated to be 4.2% on q8w dosing and 6.7% on q4w dosing. CONCLUSIONS: The developed models can guide clinical trial design and support model-informed dose optimization (stratified or individualized dosing) to improve endoscopic outcomes.

Model-Informed Precision Dosing during Infliximab Induction Therapy Reduces Variability in Exposure and Endoscopic Improvement between Patients with Ulcerative Colitis.

Model-informed precision dosing (MIPD) may be a solution to therapeutic failure of infliximab for patients with ulcerative colitis (UC), as underexposure could be avoided, and the probability of endoscopic improvement (pEI; Mayo endoscopic subscore ≤ 1) could be optimized. To investigate in silico whether this claim has merit, four induction dosing regimens were simulated: 5 mg/kg (label dosing), 10 mg/kg, covariate-based MIPD (fat-free mass, corticosteroid use, and presence of extensive colitis at baseline), and concentration-based MIPD (based on the trough concentration at day 14). Covariate- and concentration-based MIPD were chosen to target the same median area under the infliximab concentration-time curve up to endoscopy at day 84 (AUCd84), as was predicted from 10 mg/kg dosing. Dosing at 5 mg/kg resulted in a mean ± standard deviation pEI of 55.7 ± 9.0%. Increasing the dose to 10 mg/kg was predicted to improve pEI to 65.1 ± 6.1%. Covariate-based MIPD reduced variability in exposure and pEI (65.1 ± 5.5%). Concentration-based MIPD decreased variability further (66.0 ± 3.9%) but did so at an increased average dose of 2293 mg per patient, as compared to 2168 mg for 10 mg/kg dosing. Mean pEI remained unchanged between 10 mg/kg dosing and MIPD, since the same median AUCd84 was targeted. In conclusion, quantitative simulations predict MIPD will reduce variability in exposure and pEI between patients with UC during infliximab induction therapy.

Optimal empiric treatment for KPC-2-producing Klebsiella pneumoniae infections in critically ill patients with normal or decreased renal function using Monte Carlo simulation.

BACKGROUND: Limited clinical studies describe the pharmacodynamics of fosfomycin (FOS), tigecycline (TGC) and colistin methanesulfonate (CMS) in combination against KPC-producing Klebsiella pneumoniae (KPC-Kp). Population pharmacokinetic models were used in our study. Monte Carlo simulation was conducted to calculate probability of target attainment (PTA) and cumulative fraction of response (CFR) of each agent alone and in combination against KPC-Kp in patients with normal or decreased renal function. RESULTS: The simulated regimen of FOS 6 g q8h reached ≥90% PTA against a MIC of 64 mg/L in patients with normal renal function. For patients with renal impairment, FOS 4 g q8h could provide sufficient antimicrobial coverage against a MIC of 128 mg/L. And increasing the daily dose could result to the cut-off value to 256 mg/L in decreased renal function. For TGC, conventional dosing regimens failed to reach 90% PTA against a MIC of 2 mg/L. Higher loading and daily doses (TGC 200/400 mg loading doses followed by 100 mg q12h/200 mg q24h) were needed. For CMS, none achieved 90% PTA against a MIC of 2 mg/L in normal renal function. Against KPC-Kp, the regimens of 200/400 mg loading dose followed by 100 q12h /200 mg q24h achieved > 80% CFRs regardless of renal function, followed by CMS 9 million IU loading dose followed by 4.5/3 million IU q12h in combination with FOS 8 g q8h (CFR 75-91%). CONCLUSIONS: The use of a loading dose and high daily dose of TGC and CMS in combination with FOS can provide sufficient antimicrobial coverage against critically ill patients infected with KPC-Kp.

Modelling of the relationship between infliximab exposure, faecal calprotectin and endoscopic remission in patients with Crohn's disease.

AIMS: Evidence for the benefits of pharmacokinetic (PK) and pharmacodynamic (PD) monitoring of infliximab in patients with Crohn's disease (CD) remains scarce. We aimed to develop a population (pop)PK/PD model to characterise the infliximab dose-exposure-biomarker-response (faecal calprotectin [fCal] and endoscopic remission [ER]) relationship. METHODS: Data were obtained from 116 patients with CD in a phase 4 dose-escalation study. Three sequential models were developed: a 2-compartment popPK model linking infliximab dose to exposure; an indirect response popPK/PD model describing the inhibitory effect of infliximab exposure on fCal; and a first-order Markov popPD model linking fCal to transitions between states of ER, no ER and dropout. RESULTS: Infliximab clearance increased with increasing fCal, decreasing albumin, increasing CD activity index and presence of anti-drug antibodies. Baseline fCal increased with increasing C-reactive protein and decreasing platelet count. Lower fCal increased the probability of attaining ER and decreased the probability of losing ER. Probability of dropping out given an earlier state of absence of ER increased with time. Large interpatient PK and PD variability resulted in a flat dose-response curve. Predicted fraction of patients achieving ER was 45% [30-61] (median [interquartile range], n = 50 000) when on 5 mg/kg infliximab (~46% observed in data). Simulations with 10 mg/kg induction doses predicted an increase to 48% [32-63]. This minor benefit at the population level argues against systematic 10 mg/kg induction dosing in all patients. CONCLUSION: Model-informed infliximab dose optimisation towards a predefined fCal concentration (while accounting for PK and PD variability) may improve effectiveness of infliximab therapy.

An integrated PK-PD model for cortisol and the 17-hydroxyprogesterone and androstenedione biomarkers in children with congenital adrenal hyperplasia.

AIMS: The aim of this study was to characterize the pharmacokinetic/pharmacodynamic relationships of cortisol and the adrenal biomarkers 17-hydroxyprogesterone and androstenedione in children with congenital adrenal hyperplasia (CAH). METHODS: A nonlinear mixed-effect modelling approach was used to analyse cortisol, 17-hydroxyprogesterone and androstenedione concentrations obtained over 6 hours from children with CAH (n = 50). A circadian rhythm was evident and the model leveraged literature information on circadian rhythm in untreated children with CAH. Indirect response models were applied in which cortisol inhibited the production rate of all three compounds using an Imax model. RESULTS: Cortisol was characterized by a one-compartment model with apparent clearance and volume of distribution estimated at 22.9 L/h/70 kg and 41.1 L/70 kg, respectively. The IC50 values of cortisol concentrations for cortisol, 17-hydroxyprogesterone and androstenedione were estimated to be 1.36, 0.45 and 0.75 µg/dL, respectively. The inhibitory effect was found to be more potent on 17OHP than D4A, and the IC50 values were higher in salt-wasting subjects than simple virilizers. Production rates of cortisol, 17-hydroxyprogesterone and androstenedione were higher in simple-virilizer subjects. Half-lives of cortisol, 17-hydroxyprogesterone and androstenedione were 60, 47 and 77 minutes, respectively. CONCLUSION: Rapidly changing biomarker responses to cortisol concentrations highlight that single measurements provide volatile information about a child's disease control. Our model closely captured observed cortisol, 17-hydroxyprogesterone and androstenedione concentrations. It can be used to predict concentrations over 24 hours and allows many novel exposure metrics to be calculated, e.g., AUC, AUC-above-threshold, time-within-range, etc. Our long-range goal is to uncover dose-exposure-outcome relationships that clinicians can use in adjusting hydrocortisone dose and timing.

Model-Based Simulation of Maintenance Therapy of Childhood Acute Lymphoblastic Leukemia.

Acute lymphoblastic leukemia is the most common malignancy in childhood. Successful treatment requires initial high-intensity chemotherapy, followed by low-intensity oral maintenance therapy with oral 6-mercaptopurine (6MP) and methotrexate (MTX) until 2-3 years after disease onset. However, intra- and inter-individual variability in the pharmacokinetics (PK) and pharmacodynamics (PD) of 6MP and MTX make it challenging to balance the desired antileukemic effects with undesired excessive myelosuppression during maintenance therapy. A model to simulate the dynamics of different cell types, especially neutrophils, would be a valuable contribution to improving treatment protocols (6MP and MTX dosing regimens) and a further step to understanding the heterogeneity in treatment efficacy and toxicity. We applied and modified a recently developed semi-mechanistic PK/PD model to neutrophils and analyzed their behavior using a non-linear mixed-effects modeling approach and clinical data obtained from 116 patients. The PK model of 6MP influenced the accuracy of absolute neutrophil count (ANC) predictions, whereas the PD effect of MTX did not. Predictions based on ANC were more accurate than those based on white blood cell counts. Using the new cross-validated mathematical model, simulations of different treatment protocols showed a linear dose-effect relationship and reduced ANC variability for constant dosages. Advanced modeling allows the identification of optimized control criteria and the weighting of specific influencing factors for protocol design and individually adapted therapy to exploit the optimal effect of maintenance therapy on survival.

A Population Pharmacokinetic and Exposure-Response Model of Golimumab for Targeting Endoscopic Remission in Patients With Ulcerative Colitis.

BACKGROUND: Unlike other anti-tumor necrosis factor alpha antibodies, golimumab does not deliver on its promise of effectiveness for treating patients with ulcerative colitis. We investigated the value of therapeutic drug monitoring for optimizing golimumab therapy. METHODS: We analyzed the golimumab pharmacokinetics data of 56 patients with moderate to severe ulcerative colitis. Induction and maintenance golimumab concentrations (296 venipuncture, 414 serum) were used to develop a population pharmacokinetic model. Exposure-response relationships were analyzed using the data of 40/56 patients with available endoscopy data. Receiver operating characteristic curve analysis was performed, and an exposure-response Markov model was developed, linking golimumab exposure to probabilities of transitioning between Mayo endoscopic subscore (MES) states from baseline to week (w)14. RESULTS: Golimumab pharmacokinetics was best described by a 2-compartment model with linear absorption and elimination. Antibodies to golimumab and previous biological therapy reduced golimumab exposure. Still, interindividual pharmacokinetic variability (IIVPK) remained largely unexplained. Endoscopic remission (ER; MESw14 ≤ 1) was achieved in 14/40 (35%) patients. Golimumab serum trough concentration thresholds of 7.4 mg/L (w6) and 3.2 mg/L (w14) predicted ER at w14 (positive predictive values [pv+] 83% and 91%, pv- 82% and 67%, respectively). The 3.2-mg/L target predicted 38% and 44% chances of achieving ER in patients with MESbaseline of 3 and 2, respectively. CONCLUSIONS: Personalized, model-based induction dosing aiming at here-established target concentrations may account for IIVPK and thus provide patients with more equal chances of achieving ER. As <50% of patients attained the exposure targets, higher golimumab induction dosing requires investigation to secure its future in clinical practice.

Translational Approach to Predicting the Efficacy of Maraviroc-Based Regimens as HIV Preexposure Prophylaxis.

Maraviroc-based regimens have been explored as preexposure prophylaxis (PrEP) against human immunodeficiency virus (HIV). In this study, we utilized mucosal tissue drug exposure data, combined with target concentrations generated in vitro, in a pharmacokinetic-pharmacodynamic analysis to predict the effects of drug combinations and adherence on PrEP efficacy. Mucosal tissue concentrations of maraviroc were measured in 24 healthy women. The 90% effective concentrations (EC90) of maraviroc (alone and combined with tenofovir and emtricitabine) for protection against HIV were identified in CD4+ T cells. Monte Carlo simulations were performed to identify dosing strategies to protect colorectal and female genital tract (FGT) tissues from HIV infection. Colorectal maraviroc concentrations were 350-fold higher than in the FGT. Under steady-state conditions, our model predicted that one 300-mg dose/week was sufficient to protect colorectal tissue from HIV in 99% of the population, while 300 mg daily would protect the FGT in only 63% of the population. FGT protection increased to >90% when maraviroc was used in combination with tenofovir (5 doses/week) or emtricitabine (3 doses/week). Poor adherence resulted in a drastic decrease in efficacy in the FGT but not colorectal tissue. However, greater forgiveness was seen when maraviroc was combined with tenofovir or emtricitabine, suggesting that maraviroc should not be used alone as PrEP.

Downregulation of ABCB1 gene in patients with total hip or knee arthroplasty influences pharmacokinetics of rivaroxaban: a population pharmacokinetic-pharmacodynamic study.

PURPOSE: Rivaroxaban is a substrate for ABCB1 transporter and is commonly used in patients undergoing hip or knee replacement surgery for thromboprophylaxis. The objective of this study was to develop a population pharmacokinetic-pharmacodynamic (PK-PD) model to investigate the influence of ABCB1 gene expression and polymorphism on rivaroxaban exposure and anticoagulation effects. METHODS: Five blood samples per patient were collected during 5 days after the surgery for the determination of rivaroxaban concentration in plasma and for determination of prothrombin time and partial thromboplastin time. Non-linear mixed effects model was used for a population PK-PD analysis and for testing covariate effects. RESULTS: A one-compartment PK model with first-order absorption adequately described the pharmacokinetic data. The typical oral clearance (CL/F) was 6.12 L/h (relative standard error, 15.8%) and was associated with ABCB1 expression. Compared to base line before the surgery, a significant ABCB1 downregulation was observed 5 days after the surgery (p < 0.001). Prothrombin time and partial thromboplastin time were both linearly associated to the logarithm of the rivaroxaban plasma concentration. CONCLUSIONS: We confirmed that variable rivaroxaban CL/F is associated with ABCB1 expression, which is in accordance with previous studies on P-glycoprotein involvement in rivaroxaban PK. Furthermore, we observed the downregulation of ABCB1 expression after the surgery. The cause remains unclear and further research is needed to explain the underlying mechanisms.

Optimising infliximab induction dosing for patients with ulcerative colitis.

AIMS: The therapeutic failure of infliximab therapy in patients with ulcerative colitis remains a challenge even 2 decades after its approval. Therapeutic drug monitoring (TDM) has shown value during maintenance therapy, but induction therapy has still not been explored. Patients may be primary nonresponders or underexposed with the standard dosing regimen. We aimed to: (i) develop a population pharmacokinetic-pharmacodynamic model; (ii) identify the best exposure metric that predicts mucosal healing; and (iii) build an exposure-response (ER) model to demonstrate model-based dose finding during induction therapy with infliximab. METHODS: Data were retrospectively collected from a clinical database. A total of 583 samples, from 204 patients, was used to develop a population pharmacokinetic model to generate exposure metrics for subsequent ER modelling. A subset of 159 patients was used to develop a logistic regression ER model, describing the relationship between infliximab exposure and ordered transitions between Mayo endoscopic subscore (MES) 3, 2 and ≤1 (baseline to post-induction). RESULTS: A 1-compartment population pharmacokinetic model with interindividual and interoccasion variability was found to fit the data best. Covariates influencing exposure were C-reactive protein, albumin, baseline MES, fat-free mass, concomitant corticosteroid use and pancolitis. The cumulative area under the infliximab concentration-time curve until endoscopy (CAUCendoscopy ) was found to be the best exposure metric for predicting mucosal healing (baseline MES >1 and post-induction MES ≤1). The model predicted that 70% of patients will attain mucosal healing with infliximab administered at days 0, 14 and 42 and a target CAUCendoscopy of 3752 mg/L*day at day 84. CONCLUSIONS: TDM-based dose individualisation targeting CAUCendoscopy has the potential to improve the effectiveness of infliximab during induction therapy.

Development of an Enantioselective and Biomarker-Informed Translational Population Pharmacokinetic/Pharmacodynamic Model for Etodolac.

Cyclooxygenase-2 (COX-2) isoform has a critical role in the development of pain. Inhibition of COX-2 in vitro serves as a biomarker for nonsteroidal anti-inflammatory drugs (NSAIDs). The NSAID concentrations yielding 80% COX-2 inhibition (IC80) correlate with therapeutic doses to achieve analgesia across multiple COX-2 inhibitors. However, there are no time-course models relating COX-2 inhibition with decreased pain. This study aimed to characterize the relationship between NSAID concentrations, in vitro COX-2 inhibition, and acute pain decrease in humans over time by a translational approach using clinical pharmacokinetic and literature reported in vitro and clinical pharmacodynamic data. In a two-way cross-over study, eight healthy volunteers received 300 and 400 mg racemic etodolac, a preferential COX-2 inhibitor. R- and S-etodolac were determined by LC-MS/MS and simultaneously modeled. Literature in vitro IC50 data for COX-2 inhibition by S-etodolac were used to fit adjusted pain score profiles from dental patients receiving etodolac. External model qualification was performed using published ibuprofen data. Etodolac absorption was highly variable due to gastric transit kinetics and low aqueous solubility. The disposition parameters differed substantially between enantiomers with a total clearance of 2.21 L/h for R-etodolac and 26.8 L/h for S-etodolac. Volume of distribution at steady-state was 14.6 L for R-etodolac and 45.8 L for S-etodolac. Inhibition of COX-2 by 78.1% caused a half-maximal pain decrease. The time-course of pain decrease following ibuprofen was successfully predicted via the developed translational model. This proposed enantioselective pharmacodynamic-informed approach presents the first quantitative time-course model for COX-2 induced pain inhibition in patients.