Pharmacokinetic-pharmacodynamic Model Research Articles

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

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. 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 (~6mg/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. 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. The developed models can guide clinical trial design and support model-informed dose optimization (stratified or individualized dosing) to improve endoscopic outcomes.

Predicting Inter-individual Variability During Lipid Resuscitation of Bupivacaine Cardiotoxicity in Rats: A Virtual Population Modeling Study.

IntroductionIntravenous lipid emulsions (ILE) have been credited for successful resuscitation in drug intoxication cases where other cardiac life-support methods have failed. However, inter-individual variability can function as a confounder that challenges our ability to define the scope of efficacy for lipid interventions, particularly as relevant data are scarce. To address this challenge, we developed a quantitative systems pharmacology model to predict outcome variability and shed light on causal mechanisms in a virtual population of rats subjected to bupivacaine toxicity and ILE intervention.Materials and MethodsWe combined a physiologically based pharmacokinetic–pharmacodynamic model with data from a small study in Sprague-Dawley rats to characterize individual-specific cardiac responses to lipid infusion. We used the resulting individual parameter estimates to posit a population distribution of responses to lipid infusion. On that basis, we constructed a large virtual population of rats (N = 10,000) undergoing lipid therapy following bupivacaine cardiotoxicity.ResultsUsing unsupervised clustering to assign resuscitation endpoints, our simulations predicted that treatment with a 30% lipid emulsion increases bupivacaine median lethal dose (LD50) by 46% when compared with a simulated control fluid. Prior experimental findings indicated an LD50 increase of 48%. Causal analysis of the population data suggested that muscle accumulation rather than liver accumulation of bupivacaine drives survival outcomes.ConclusionOur results represent a successful prediction of complex, dynamic physiological outcomes over a virtual population. Despite being informed by very limited data, our mechanistic model predicted a plausible range of treatment outcomes that accurately predicts changes in LD50 when extrapolated to putatively toxic doses of bupivacaine. Furthermore, causal analysis of the predicted survival outcomes indicated a critical synergy between scavenging and direct cardiotonic mechanisms of ILE action.Supplementary InformationThe online version contains supplementary material available at 10.1007/s40268-021-00353-4.

Optimizing Oral Targeted Anticancer Therapies Study for Patients With Solid Cancer: Protocol for a Randomized Controlled Medication Adherence Program Along With Systematic Collection and Modeling of Pharmacokinetic and Pharmacodynamic Data.

BackgroundThe strengthening or substitution of intravenous cytotoxic chemotherapy cycles by oral targeted anticancer therapies, such as protein kinase inhibitors (PKIs), has provided impressive clinical benefits and autonomy as well as a better quality of life for patients with cancer. Despite these advances, adverse event management at home and medication adherence remain challenging. In addition, PKI plasma concentrations vary significantly among patients with cancer receiving the same dosage, which could explain part of the observed variability in the therapeutic response.ObjectiveThe aim of this optimizing oral targeted anticancer therapies (OpTAT) study is to optimize and individualize targeted anticancer treatments to improve patient care and self-monitoring through an interprofessional medication adherence program (IMAP) combined with measurement PKI plasma concentrations.MethodsThe OpTAT study has two parts: (1) a 1:1 randomized medication adherence program, in which the intervention consists of regular motivational interviewing sessions between the patient and the pharmacist, along with the delivery of PKIs in electronic monitors, and (2) a systematic collection of blood samples and clinical and biological data for combined pharmacokinetic and pharmacodynamic analysis. On the basis of the electronic monitor data, medication adherence will be compared between groups following the three operational definitions: implementation of treatment during the persistent period, persistence with treatment and longitudinal adherence. The implementation will be described using generalized estimating equation models. The persistence of PKI use will be represented using a Kaplan-Meier survival curve. Longitudinal adherence is defined as the product of persistence and implementation. PKI pharmacokinetics will be studied using a population approach. The relationship between drug exposure and efficacy outcomes will be explored using Cox regression analysis of progression-free survival. The relationship between drug exposure and toxicity will be analyzed using a pharmacokinetic-pharmacodynamic model and by logistic regression analysis. Receiver operating characteristic analyses will be applied to evaluate the best exposure threshold associated with clinical benefits.ResultsThe first patient was included in May 2015. As of June 2021, 262 patients had participated in at least one part of the study: 250 patients gave at least one blood sample, and 130 participated in the adherence study. Data collection is in process, and the final data analysis is planned to be performed in 2022.ConclusionsThe OpTAT study will inform us about the effectiveness of the IMAP program in patients with solid cancers treated with PKIs. It will also shed light on PKI pharmacokinetic and pharmacodynamic properties, with the aim of learning how to adapt the PKI dosage at the individual patient level to increase PKI clinical suitability. The IMAP program will enable interprofessional teams to learn about patients’ needs and to consider their concerns about their PKI self-management, considering the patient as an active partner.Trial RegistrationClinicalTrials.gov NCT04484064; https://clinicaltrials.gov/ct2/show/NCT04484064.International Registered Report Identifier (IRRID)DERR1-10.2196/30090

Sustained Minimal Residual Disease Negativity Predicted in Chronic Lymphocytic Leukemia Patients Treated with Venetoclax Combination Therapy for 2 Years: an Integrated Mechanistic Analysis of Multiple Phase I and II Studies

Sustained Minimal Residual Disease Negativity Predicted in Chronic Lymphocytic Leukemia Patients Treated with Venetoclax Combination Therapy for 2 Years: an Integrated Mechanistic Analysis of Multiple Phase I and II Studies

Developing pharmacokinetics-pharmacodynamics model of valproic acid syrup based on prediction of population pharmacokinetics parameter and seizure frequency in Indonesian pediatric epilepsy outpatients.

Valproic acid (VPA) is a broad-spectrum antiepileptic drug with known efficacy profile in pediatric patients, despite of its narrow therapeutic index. There is lack of VPA's pharmacokinetics profile in Indonesian pediatric subjects, partly due to limited pediatric blood volume taken for conducting therapeutic drug monitoring. This study aimed to determine the correlation between VPA pharmacokinetics parameters based on population data and seizure frequency in pediatric epilepsy outpatients. This observational study was conducted at Sanglah General Hospital during June-December 2019. The subjects of this research were 38 pediatric epilepsy patients who adhered to VPA syrup monotherapy for at least 3weeks. Five subjects randomly selected for blood sample collection. Thus, VPA concentration level in the blood being analysed as a comparison to its concentration predicted from Yukawa's steady state equation. Monolix2019R2® software was used to identify VPA population pharmacokinetics-pharmacodynamics (PK-PD) parameters at steady state level. Population PK-PD of VPA syrup at steady state levelwere ka_pop = 6.25/h, Vd_pop = 3.36L, Cl_pop = 3.17·e-11mL/min, IC50_pop = 1.85·e-6, correlation of Vd_pop and Cl_pop = 0.966. Kendall Tau Correlation of predicted VPA steady state concentration and frequency of seizure was-0.66. Mean prediction error between predicted steady state concentration of five subjects and their related blood levels was≤25% and considered as within clinically acceptable limit. It needs further study to develop best matched PK-PD model of VPA syrup at steady state condition in pediatric epilepsy.

Pharmacokinetic-Pharmacodynamic Model of Vedolizumab for Targeting Endoscopic Remission in Patients With Crohn Disease: Posthoc Analysis of the LOVE-CD Study.

BackgroundHigher serum concentrations of vedolizumab have been associated with improved outcomes in inflammatory bowel disease. It is unclear how vedolizumab exposure is linked to endoscopic remission in Crohn disease (CD). We aimed to develop a pharmacokinetic-pharmacodynamic model linking vedolizumab exposure to endoscopic remission in CD.MethodsData were obtained from the first 110 patients participating in a phase 4 prospective multicenter trial (LOVE-CD; ClinicalTrials.gov identifier: NCT02646683), where vedolizumab was dosed at 300 mg every 8 weeks and serum concentrations and antibodies to vedolizumab were measured before each infusion. Concentration-time profiles were described by a 2-compartment model with parallel linear and nonlinear elimination. A first-order discrete-time Markov model was used to describe the relationship between pharmacokinetic exposure metrics and the probability of endoscopic remission (Simple Endoscopic Score for CD < 4).ResultsLinear clearance was 0.215 L/d, and the volume of distribution of the central compartment was 4.92 L. Linear clearance was higher and vedolizumab exposure was lower in patients with lower serum albumin concentrations, in the presence of antibodies to vedolizumab, and in patients with previous exposure to other biologic therapy. A week 22 vedolizumab concentration of 20.0 mg/L was predicted to yield a 35% probability of achieving endoscopic remission at week 26. Model-based simulations suggested that endoscopic remission rates of 46.5% or 40.0% could be reached with every-4-weeks dosing in patients who were naïve or previously exposed to biologic therapy, respectively.ConclusionsModel-informed dosing of vedolizumab in CD provides a foundation for future research aiming to maximize endoscopic remission rates.

Application of Physiologically Based Pharmacokinetic-Pharmacodynamic Modeling in Preterm Neonates to Guide Gentamicin Dosing Decisions and Predict Antibacterial Effect.

Clinical studies in preterm neonates are rarely performed due to ethical concerns and difficulties associated with trials and recruitment. Consequently, dose selection in this population is primarily empirical. Scaling neonatal doses from adult doses does not account for developmental changes and may not accurately predict drug kinetics. This is especially important for gentamicin, a narrow therapeutic index aminoglycoside antibiotic. While gentamicin's bactericidal effect is associated with its peak plasma concentration, keeping trough concentrations below 1 µg/mL prevents toxicity and also helps to counteract adaptive resistance in bacteria such as Escherichia coli. In this study, physiologically based pharmacokinetic-pharmacodynamic (PBPK-PD) modeling was used to support and/or guide dosing decisions and to predict the antibacterial effect in preterm neonates. A gentamicin PBPK model was successfully verified in healthy adults and preterm neonates across all gestational ages. Clinical data from a neonatal intensive care unit at NYU Langone Hospital-Long Island was used to identify dosing regimens associated with increased incidence of elevated gentamicin trough concentrations in different preterm patient cohorts. Model predictions demonstrated that a higher dose with an extended-dosing interval (every 36 hours) in neonates with a postmenstrual age of 30 to 34 weeks and ≥35 weeks, with postnatal age 8 to 28 days and 0 to 7 days, respectively, were more likely to have a trough <1 µg/mL when compared with once-daily (every 24 hours) dosing. PBPK-PD modeling suggested that a higher dose administered every 36 hours may provide effective antibacterial therapy.

P307 Modelling of the relationship between ustekinumab exposure, faecal calprotectin and endoscopic outcomes in patients with Crohn’s disease

Abstract Background In the UNITI endoscopy sub-study, only 17.4% of patients with Crohn’s disease (CD) on ustekinumab achieved endoscopic response and 10.9% achieved endoscopic remission at week (w)44. We aimed to investigate if improved endoscopic outcomes can be achieved through dose optimisation based on a population pharmacokinetic-pharmacodynamic (popPK-PD) modelling and simulation analysis. Methods Real-world data were obtained from 83 patients with moderate-to-severe CD (94% multi-refractory) enrolled in a prospective cohort study receiving ustekinumab 6 mg/kg induction and every eight-week (q8w) 90 mg maintenance therapy. Ustekinumab serum concentrations were measured at mid-dose (w4) and trough (w8, w16, w24). Faecal calprotectin (fCal) was measured at baseline and at w4, w8, w16, w24. Endoscopic response (≥50% decrease in simple endoscopic score for CD [SES-CD]) and endoscopic remission (SES-CD ≤2) were assessed at w24. Modelling and simulation were performed using NONMEM 7.4. Results 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 fCal, and a logistic regression popPD model linking fCal at w8 to endoscopic outcomes at w24 (Figure 1). Ustekinumab clearance increased with decreasing serum albumin and increasing bodyweight. The terminal half-life of ustekinumab in a median patient (bodyweight 65 kg, serum albumin 42.7 g/L) was 20.4 days. 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 1,800 μg/g to 694 μg/g (Figure 2a) The probability of endoscopic remission at w24 increased from 2.1% to 10.0% with fCal at w8 decreasing from 1,800 μg/g to 214 μg/g (Figure 2b).The results from the simulation-based comparison of q8w and q4w maintenance dosing are shown in Table 1. Dose doubling (180 mg q8w), as opposed to interval halving (90 mg q4w), was predicted to result in a ustekinumab trough concentration of 2.4 μg/mL instead of 4.8 μg/mL. Conclusion The developed model can guide clinical trial design and support model-informed dose optimisation to improve endoscopic outcome rates. Although our analyses showed that q4w dosing resulted in higher ustekinumab and lower fCal concentrations, the proportion of patients achieving endoscopic remission was limited.

Pharmacokinetic-pharmacodynamic modelling of atazanavir in hair among adolescents on antiretroviral treatment in Zimbabwe

BackgroundDrug potency is a pharmacological parameter defining dose or concentration of drug required to obtain 50% of the drug’s maximal effect. Pharmacokinetic-pharmacodynamic modelling and simulation allows estimation of potency and evaluate strategies improving treatment outcome. The objective of our study is to determine potency of atazanavir in hair, defined as atazanavir level in hair associated with 50% probability of failing to achieve viral load below 1000 copies/ml among adolescents, and explore the effect of participant specific variables on potency.MethodsA secondary analysis was performed on data from a previous study conducted in HIV-infected adolescents failing 2nd line ART from Harare central hospital, Zimbabwe, between 2015 and 2016. We simulated atazanavir concentrations in hair using NONMEM (version 7.3) ADVAN 13, based on a previously established pharmacokinetic model. Logistic regression methods were used for PKPD analysis. Simulations utilising PKPD model focused on estimation of potency and exploring the effect of covariates.ResultsThe potency of atazanavir in hair was found to be 4.5 ng/mg hair before adjusting for covariate effects. Participants at three months follow-up, reporting adequate adherence, having normal BMI-for-age, and cared for by mature guardians had increased potency of atazanavir in hair of 2.6 ng/mg, however the follow-up event was the only statistically significant factor at 5% level.ConclusionAtazanavir in hair in the range 2.6 to 4.5 ng/mg is associated with above 50% probability of early viral load suppression. Adherence monitoring to adolescents with lower potency of atazanavir is recommended. The effect self-reported adherence level, BMI-for-age, and caregiver status require further evaluation.

Towards Evidence-Based Weaning: a Mechanism-Based Pharmacometric Model to Characterize Iatrogenic Withdrawal Syndrome in Critically Ill Children

For the management of iatrogenic withdrawal syndrome (IWS) in children, a quantitative understanding of the dynamics of IWS of commonly used opioids and sedatives is lacking. Here, we introduce a new mechanism-based pharmacokinetic-pharmacodynamic (PKPD) modeling approach for studying IWS in pediatric clinical datasets. One thousand seven hundred eighty-two NRSwithdrawal scores of IWS severity were analyzed, which were collected from 81 children (age range: 1 month–18 years) that received opioids or sedatives by continuous infusion for 5 days or more. These data were successfully fitted with a PKPD model consisting of a plasma and a dependence compartment that well characterized the dynamics of IWS from morphine, fentanyl, and ketamine. The results suggest that (1) instead of decreasing the infusion rate by a set percentage at set intervals, it would be better to lengthen the weaning period when higher infusion rates are administered prior to weaning; (2) for fentanyl specifically, the risk of IWS might be lower when weaning with smaller dose reductions every 12 h instead of weaning with greater dose reductions every 48 h. The developed PKPD model can be used to evaluate the risk of IWS over time and the extent to which it is affected by different weaning strategies. The results yield hypotheses that could guide future clinical research on optimal weaning strategies. The mechanism-based PKPD modeling approach can be applied in other datasets to characterize the IWS dynamics of other drugs used in pediatric intensive care.Graphical abstract

The Combination of Cell Cultured Technology and In Silico Model to Inform the Drug Development.

Human-derived in vitro models can provide high-throughput efficacy and toxicity data without a species gap in drug development. Challenges are still encountered regarding the full utilisation of massive data in clinical settings. The lack of translated methods hinders the reliable prediction of clinical outcomes. Therefore, in this study, in silico models were proposed to tackle these obstacles from in vitro to in vivo translation, and the current major cell culture methods were introduced, such as human-induced pluripotent stem cells (hiPSCs), 3D cells, organoids, and microphysiological systems (MPS). Furthermore, the role and applications of several in silico models were summarised, including the physiologically based pharmacokinetic model (PBPK), pharmacokinetic/pharmacodynamic model (PK/PD), quantitative systems pharmacology model (QSP), and virtual clinical trials. These credible translation cases will provide templates for subsequent in vitro to in vivo translation. We believe that synergising high-quality in vitro data with existing models can better guide drug development and clinical use.

Pharmacokinetic-Pharmacodynamic Modelling of Systemic IL13 Blockade by Monoclonal Antibody Therapy: A Free Assay Disguised as Total.

A sequential pharmacokinetic (PK) and pharmacodynamic (PD) model was built with Nonlinear Mixed Effects Modelling based on data from a first-in-human trial of a novel biologic, MEDI7836. MEDI7836 is a human immunoglobulin G1 lambda (IgG1λ-YTE) monoclonal antibody, with an Fc modification to reduce metabolic clearance. MEDI7836 specifically binds to, and functionally neutralizes interleukin-13. Thirty-two healthy male adults were enrolled into a dose-escalation clinical trial. Four active doses were tested (30, 105, 300, and 600 mg) with 6 volunteers enrolled per cohort. Eight volunteers received placebo as control. Following single subcutaneous administration (SC), individual time courses of serum MEDI7836 concentrations, and the resulting serum IL13 modulation in vivo, were quantified. A binding pharmacokinetic-pharmacodynamic (PK-PD) indirect response model was built to characterize the exposure-driven modulation of the target over time by MEDI7836. While the validated bioanalytical assay specification quantified the level of free target (i.e., a free IL13 assay), emerging clinical data suggested dose-dependent increase in systemic IL13 concentration over time, indicative of a total IL13 assay. The target time course was modelled as a linear combination of free target and a percentage of the drug-target complex to fit the clinical data. This novel PK-PD modelling approach integrates independent knowledge about the assay characteristics to successfully elucidate apparently complex observations.

Pharmacokinetic-pharmacodynamic modeling for hepatic delivery and efficacy of antisense oligonucleotides with lipophilic ligands in mice.

Conjugation with lipophilic ligands such as cholesterol and α-tocopherol dramatically improves the delivery and efficacy of antisense oligonucleotides (ASOs) in the liver. To estimate the hepatic ASO concentration and the efficacy of ASOs conjugated with lipophilic ligands in mice, we constructed a pharmacokinetic-pharmacodynamic (PK-PD) model that consisted of a two-linear compartment model for the plasma and the hepatic ASO concentration, and two indirect response models for the hepatic apolipoprotein B (Apo-B) mRNA and plasma total cholesterol. The model provided a good fit of the hepatic ASO concentration although it showed an overprediction of Apo-B mRNA and an underprediction of the plasma total cholesterol within 2-fold at a later time after single intravenous administration of ASOs conjugated with lipophilic ligands. In addition, the model simulations indicated that the efficacy at a dose regimen of ASOs conjugated with lipophilic ligands (0.2mg/kg, once a week) in mice was comparable to that at an effective dose of unchanged ASO (2.5mg/kg, once a week). Although further studies are required to refine the parameters of the PK-PD model, this approach could be used to guide dose-ranging pharmacological studies for ASOs conjugated with lipophilic ligands in mice.

Pharmacokinetic-pharmacodynamic modeling approach for dose prediction of the optimal long-acting injectable formulation of finasteride

A controlled drug release formulation based on the subcutaneous injection of poly (lactic-co-glycolic acid) (PLGA) microspheres loaded with finasteride was prepared and evaluated for monthly delivery. After selection of biodegradable polymer and polymer-to-finasteride ratio, the formulation was characterized. Scanning electron microscopy (SEM) and laser-light particle size analysis were used to examine the morphology, surface structure, and particle size. High‑performance liquid chromatography (HPLC) was used to determine the drug loading, while liquid chromatography with tandem mass spectrometry (LC-MS/MS) was employed to analyze plasma finasteride concentrations. Results showed that the PLGA microspheres were spherical and of an appropriate size. The formulation stably releases the drug from the microspheres and the release sustained for a month without burst release, which was the desired duration. In vivo pharmacokinetic-pharmacodynamic (PK-PD) studies were conducted in beagle dogs through the administration of PROPECIA® (as a reference drug) per oral and subcutaneous injection of the long-acting injectable microsphere formulation (LAIF) loaded with five different doses of finasteride. From the acquired plasma data, PK-PD models for both PROPECIA®-administered group and LAIFs-injected groups were developed and validated. PK-PD profiles of both groups were predicted for up to one month. The predicted PK-PD profile of all LAIFs showed the achievability of monthly delivery and pharmacological effects without burst release, compared to the simulated PK-PD profile of PROPECIA®. According to the predicted PK-PD profiles, the formulation loaded with 16.8mg of finasteride was determined to be the optimal dose. The data obtained from the PK-PD model could be used as the basis for the estimation of a first-in-human dose of the formulation.

A Novel Integrated Pharmacokinetic-Pharmacodynamic Model to Evaluate Combination Therapy and Determine In Vivo Synergism.

Understanding pharmacokinetic (PK)-pharmacodynamic (PD) relationships is essential in translational research. Existing PK-PD models for combination therapy lack consideration of quantitative contributions from individual drugs, whereas interaction factor is always assigned arbitrarily to one drug and overstretched for the determination of in vivo pharmacologic synergism. Herein, we report a novel generic PK-PD model for combination therapy by considering apparent contributions from individual drugs coadministered. Doxorubicin (Dox) and sorafenib (Sor) were used as model drugs whose PK data were obtained in mice and fit to two-compartment model. Xenograft tumor growth was biphasic in mice, and PD responses were described by three-compartment transit models. This PK-PD model revealed that Sor (contribution factor = 1.62) had much greater influence on overall tumor-growth inhibition than coadministered Dox (contribution factor = 0.644), which explains the mysterious clinical findings on remarkable benefits for patients with cancer when adding Sor to Dox treatment, whereas there were none when adding Dox to Sor therapy. Furthermore, the combination index method was integrated into this predictive PK-PD model for critical determination of in vivo pharmacologic synergism that cannot be correctly defined by the interaction factor in conventional models. In addition, this new PK-PD model was able to identify optimal dosage combination (e.g., doubling experimental Sor dose and reducing Dox dose by 50%) toward much greater degree of tumor-growth inhibition (>90%), which was consistent with stronger synergy (combination index = 0.298). These findings demonstrated the utilities of this new PK-PD model and reiterated the use of valid method for the assessment of in vivo synergism.Significance StatementA novel pharmacokinetic (PK)-pharmacodynamic (PD) model was developed for the assessment of combination treatment by considering contributions from individual drugs, and combination index method was incorporated to critically define in vivo synergism. A greater contribution from sorafenib to tumor-growth inhibition than that of coadministered doxorubicin was identified, offering explanation for previously inexplicable clinical observations. This PK-PD model and strategy shall have broad applications to translational research on identifying optimal dosage combinations with stronger synergy toward improved therapeutic outcomes.

Pharmacokinetic-pharmacodynamic modeling analysis and anti-inflammatory effect of Wangbi capsule in the treatment of adjuvant-induced arthritis.

Clinically, Wangbi Capsule (WBC) is widely used in the treatment of Rheumatoid arthritis (RA) because of its remarkable therapeutic effect. To reveal the mechanism, a pharmacokinetic-pharmacodynamic (PK-PD) model was developed for the first time to assess the relationship between time-concentration (dose)-effect. Freund's Complete Adjuvant was used to induce the adjuvant-induced arthritis model. Multi-indices were used to evaluate the therapeutic effect and an S-Imax PK-PD model was established based on the concentrations of osthole, 5-O-methylvisamminoside, cimifugin, albiflorin, paeoniflorin and icariin and the levels of interleukin-1β and prostaglandin E2 using a two-compartment PK model together with a PD model with an effect-site compartment. The results suggest that WBC can treat RA by regulating the levels of prostaglandin E2 and interleukin-1β. For the PK-PD model, the parameters indicated that WBC had a large safety margin and all six bioactive ingredients of WBC have therapeutic effects on RA. Among them icariin, osthole and 5-O-methylvisamminoside may be the main effective substances. This study provided a scientific basis for further study of population pharmacokinetics / population pharmacodynamics (PPK/PPD), to develop a reasonable administration plan and improve individualized drug therapy.

Proposal of a prediction model describing the metabolic interaction between warfarin and sorafenib using a population pharmacokinetic approach.

We attempted to examine the applicability of a population pharmacokinetic-pharmacodynamic (PK-PD) model describing the metabolic interaction between warfarin and sorafenib due to CYP2C9 inhibition and to predict the plasma concentrations of sorafenib and S-warfarin, the international normalized ratio (INR), and the optimal maintenance dose of warfarin in the presence of sorafenib in vivo. The sorafenib inhibition constant for S-warfarin metabolism was determined in vitro, and the unbound fraction in the liver was estimated using the published equations. A population PK-PD model describing the interaction between warfarin and sorafenib assuming competitive metabolic inhibition of S-warfarin by sorafenib was developed using NONMEM. The model was evaluated using clinical data and INR collected from the literature. The observed time courses of INR retrieved from Japanese and Caucasian patients given warfarin and sorafenib were mostly within the 90% range of the predicted values. Then, we predicted the plasma sorafenib and S-warfarin concentrations and INR after administration of warfarin (3mg/day) alone and warfarin + sorafenib (800mg/day). The predicted mean plasma S-warfarin concentration and INR at steady state were almost 6 and 2times greater, respectively, in the presence of sorafenib than those for warfarin alone. The predicted S-warfarin concentrations and INR after reduction of the warfarin dose (0.5mg/day) in the presence of sorafenib were comparable to those after 3mg/day warfarin alone. The proposed population PK-PD model has the potential to predict an increase in INR quantitatively after concurrent administration of warfarin and sorafenib.

A Physiologically Based Pharmacokinetic-Pharmacodynamic Model for Capecitabine in Colorectal Cancer Rats: Simulation of Antitumor Efficacy at Various Administration Schedules.

Capecitabine is an oral prodrug of 5-fluorouracil and is widely used for colorectal cancer (CRC) treatment. However, knowledge of its antitumor efficacy after modification of the dosing schedule is insufficient. The aim of this study was to predict the antitumor efficacy of capecitabine using a physiologically based pharmacokinetic-pharmacodynamic (PBPK-PD) model based on metabolic enzyme activities. CRC model rats were administrated 180mg/kg of capecitabine for 2weeks. Blood samples were collected at 0, 0.25, 0.5, 0.75, 1, 1.5, 2, 3, 4, 6, and 8h following capecitabine administration. Plasma concentrations of capecitabine and its metabolites were measured on days 1, 7, and 14. Metabolic enzyme activities were determined in vitro using the liver and small intestine of the CRC model rats. A PBPK-PD model was developed based on metabolic enzyme activities. The antitumor efficacy of capecitabine after regimen modification was simulated using the PBPK-PD model. Capecitabine antitumor efficacy was dose-dependent. A dose of >500μmol/kg was needed to inhibit tumor growth. After capecitabine regimen modification, a 1-week postponement of capecitabine administration was more efficacious than a reduction in the dosage to 80%. The PBPK-PD model could simulate the antitumor efficacy at various capecitabine administration schedules. PBPK-PD models can contribute to the development of an appropriate CRC chemotherapy regimen with capecitabine.

Antinociceptive effects of levomethadone in standing horses sedated with romifidine

To evaluate the antinociceptive effect of a bolus of intravenous levomethadone administered to horses during romifidine constant rate infusion (CRI). Prospective, randomized, masked, crossover experimental study. A group of eight adult Warmblood horses (seven geldings, one mare) aged 6.6 ± 4.4 years, weighing 548 ± 52 kg [mean ± standard deviation (SD)]. Levomethadone 0.1 mg kg-1 or an equivalent volume of saline (control) was administered intravenously to standing horses 60 minutes after starting a romifidine CRI. Blood samples to quantify romifidine and levomethadone plasma concentrations by capillary electrophoresis were collected up to 150 minutes after levomethadone administration. The nociceptive withdrawal reflex threshold (NWRT) was determined continuously using an automated threshold tracking device. Sedation and cardiopulmonary variables were assessed at regular intervals. A pharmacokinetic-pharmacodynamic (PK-PD) model was elaborated. Data are presented as mean ± SD or median (interquartile range, 25%-75%) where appropriate. Differences between groups were considered statistically significant for p < 0.05. Horses exhibited higher NWRTs after levomethadone administration than after saline (123 ± 9% versus 101 ± 9% relative to baseline, p < 0.05). The PK-PD model identified a contribution of levomethadone to the NWRT increase. Effect size was variable among individuals. No adverse reactions to levomethadone administration were observed. A slight effect of levomethadone on sedation scores was evident for the 60 minutes following its administration. A single injection of levomethadone has the potential to increase the NWRT during romifidine CRI in horses and can be administered in combination with α2-adrencoceptor agonists to enhance antinociception in horses. However, individual variation is marked.

Population pharmacokinetics of mycophenolic acid in paediatric patients.

Mycophenolic acid (MPA) is widely used in paediatric kidney transplant patients and sometimes prescribed for additional indications. Population pharmacokinetic or pharmacodynamic modelling has been frequently used to characterize the fixed, random and covariate effects of MPA in adult patients. However, MPA population pharmacokinetic data in the paediatric population have not been systematically summarized. The objective of this narrative review was to provide an up-to-date critique of currently available paediatric MPA population pharmacokinetic models, with emphases on modelling techniques, pharmacological findings and clinical relevance. PubMed and EMBASE were searched from inception of database to May 2020, where a total of 11 studies have been identified representing kidney transplant (n = 4), liver transplant (n = 1), haematopoietic stem cell transplant (n = 1), idiopathic nephrotic syndrome (n = 2), systemic lupus erythematosus (n = 2), and a combined population consisted of kidney, liver and haematopoietic stem cell transplant patients (n = 1). Critical analyses were provided in the context of MPA absorption, distribution, metabolism, excretion and bioavailability in this paediatric database. Comparisons to adult patients were also provided. With respect to clinical utility, Bayesian estimation models (n = 6) with acceptable accuracy and precision for MPA exposure determination have also been identified and systematically evaluated. Overall, our analyses have identified unique features of MPA clinical pharmacology in the paediatric population, while recognizing several gaps that still warrant further investigations. This review can be used by pharmacologists and clinicians for improving MPA pharmacokinetic-pharmacodynamic modelling and patient care.