Semi-Mechanistic Modelling and Simulation of Inhibition of Platelet Aggregation by Antiplatelet Agents

Physiologically-Based Pharmacokinetic-Pharmacodynamics Model Characterizing CYP2C19 Polymorphisms to Predict Clopidogrel Pharmacokinetics and Its Anti-Platelet Aggregation Effect Following Oral Administration to Coronary Artery Disease Patients With or Without Diabetes

Background and Objective: Clopidogrel (CLOP) is commonly used in coronary artery disease (CAD) patients with or without diabetes (DM), but these patients often suffer CLOP resistance, especially those with diabetes. This study was aimed to develop a physiologically-based pharmacokinetic-pharmacodynamic (PBPK-PD) model to describe the pharmacokinetics and pharmacodynamics of clopidogrel active metabolite (CLOP-AM) in CAD patients with or without DM.

Methods: The PBPK-PD model was first established and validated in healthy subjects and then in CAD patients with or without DM. The influences of CYP2C19, CYP2C9, CYP3A4, carboxylesterase 1 (CES1), gastrointestinal transit rates (K_t,i) and platelets response to CLOP-AM (k_irre) on predicted pharmacokinetics and pharmacodynamics were investigated, followed with their individual and integrated effects on CLOP-AM pharmacokinetics due to changes in DM status.

Results: Most predictions fell within 0.5–2.0 folds of observations, indicating successful predictions. Sensitivity analysis showed that contributions of interested factors to pharmacodynamics were CES1> k_irre> K_t,i> CYP2C19 > CYP3A4> CYP2C9. Mimicked analysis showed that the decreased exposure of CLOP-AM by DM was mainly attributed to increased CES1 activity, followed by decreased CYP2C19 activity.

Conclusion: The pharmacokinetics and pharmacodynamics of CLOP-AM were successfully predicted using the developed PBPK-PD model. Clopidogrel resistance by DM was the integrated effects of altered K_t,i, CYP2C19, CYP3A4, CES1 and k_irre.

Semi-Mechanistic Modeling of HY-021068 Based on Irreversible Inhibition of Thromboxane Synthetase

Background: HY-021068 [4-(2-(1H-imidazol-1-yl) ethoxy)-3-methoxybenzoate], developed by Hefei Industrial Pharmaceutical Institute Co., Ltd. (Anhui, China), is a potential thromboxane synthetase inhibitor under development as an anti-platelet agent for the treatment of stroke. A semi-mechanistic pharmacokinetic/pharmacodynamic (PK/PD) model was developed to characterize the PK of HY-021068 and its platelet aggregation inhibitory effect in beagle dogs. Method: Beagle dogs received single oral administration of 2.5 mg/kg HY-021068 or consecutively oral administration of 5 mg/kg HY-021068 once daily for 7 days. The plasma concentration of HY-021068 and the platelet aggregation rate (PAR) were determined by liquid chromatography tandem-mass spectrometry (LC-MS/MS) assay and a photometric method, respectively. The PK/PD data was sequentially fitted by Phoenix NLME. The PK/PD parameters of HY-021068 in beagle dogs were estimated by 2.5 and 5 mg/kg dosing on the 1st day, and then used to simulate the PAR of HY-021068 on the 7th day after 5 mg/kg dosing daily. Result: A one-compartment model with saturable Michaelis-Menten elimination was best fitted to the PK of HY-021068. A mechanistic PD model based on irreversible inhibition of thromboxane synthetase was constructed to describe the relationship between plasma concentration of HY-021068 and PAR. Diagnostic plots showed no obvious bias. Visual predictive check confirmed the stability and reliability of the model. Most of PK/PD observed data on the 7th day after 5 mg/kg dosing fell in the 90% prediction interval. Conclusion: We established a semi-mechanistic PK/PD model for characterizing the PK of HY-021068 and its anti-platelet effect in beagle dogs. The model can be used to predict the concentration and PAR under different dosage regimen of HY-021068, and might be served as a reference for dose design in the future clinical studies.

Population pharmacodynamics of cilostazol in healthy Korean subjects

Cilostazol is used for the treatment of intermittent claudication, ulceration and pain. This study was conducted to develop a population pharmacodynamic (PD) model for cilostazol's closure time (CT) prolongation effect in healthy Korean subjects based on a pharmacokinetic (PK) model previously developed. PD data were obtained from 29 healthy subjects who participated in a study conducted in 2009 at Severance Hospital. The PK model used was a two-compartment model with first order absorption. CT data were best described by a turnover model with a fractional turnover rate constant (K_out) inhibited by drug effects (Eff), which were represented by a sigmoid E_max model [Eff = E_max · C^γ / (EC₅₀^γ+Cγ)] with E_max being maximum drug effect, EC₅₀ drug plasma concentration at 50% of E_max, C drug plasma concentrations, and γ the Hill coefficient. For the selected PD model, parameter estimates were 0.613 hr⁻¹ for K_out, 0.192 for E_max, 730 ng/ml for EC₅₀ and 5.137 for γ. Sex and caffeine drinking status significantly influenced the baseline CT, which was 85.36 seconds in male non-caffeine drinkers and increased by 15.5% and 16.4% in females and caffeine drinkers, respectively. The model adequately described the time course of CT. This was the first population PD study for cilostazol's CT prolongation effect in a Korean population.

Transition strategies from cangrelor to oral platelet P2Y12 receptor antagonists

Cangrelor is the first parenteral antagonist of the platelet P2Y12 receptor. This direct-acting antagonist of the platelet P2Y12 receptor should be considered an adjunct to a percutaneous coronary intervention in patients who have not been adequately pretreated with platelet P2Y12 receptor antagonists at the time of the procedure. The use of cangrelor requires transition to an oral platelet P2Y12 receptor antagonist. Transition strategies have been developed on the basis of pharmacologic characteristics of platelet P2Y12 receptor antagonists, results of pharmacodynamic studies, and results from clinical trials. Cangrelor blocks the binding to the platelet P2Y12 receptor of the active metabolite of the thienopyridines, clopidogrel and prasugrel. The active metabolite of thienopyridines is present in blood for a short interval after administration. For this reason, clopidogrel should be administered after cangrelor is stopped. Prasugrel can be administered at the end of the cangrelor infusion or up to 30 min before cangrelor is stopped. Ticagrelor is also a reversible direct-acting antagonist of the platelet P2Y12 receptor. Because there is no interaction between ticagrelor and cangrelor, ticagrelor can be administered before or during the infusion of cangrelor.

Development of a physiology-directed population pharmacokinetic and pharmacodynamic model for characterizing the impact of genetic and demographic factors on clopidogrel response in healthy adults

Clopidogrel (Plavix®), is a widely used antiplatelet agent, which shows high inter-individual variability in treatment response in patients following the standard dosing regimen. In this study, a physiology-directed population pharmacokinetic/pharmacodynamic (PK/PD) model was developed based on clopidogrel and clopidogrel active metabolite (clop-AM) data from the PAPI and the PGXB2B studies using a step-wise approach in NONMEM (version 7.2). The developed model characterized the in vivo disposition of clopidogrel, its bioactivation into clop-AM in the liver and subsequent platelet aggregation inhibition in the systemic circulation reasonably well. It further allowed the identification of covariates that significantly impact clopidogrel's dose-concentration-response relationship. In particular, CYP2C19 intermediate and poor metabolizers converted 26.2% and 39.5% less clopidogrel to clop-AM, respectively, compared to extensive metabolizers. In addition, CES1 G143E mutation carriers have a reduced CES1 activity (82.9%) compared to wild-type subjects, which results in a significant increase in clop-AM formation. An increase in BMI was found to significantly decrease clopidogrel's bioactivation, whereas increased age was associated with increased platelet reactivity. Our PK/PD model analysis suggests that, in order to optimize clopidogrel dosing on a patient-by-patient basis, all of these factors have to be considered simultaneously, e.g. by using quantitative clinical pharmacology tools.

Population pharmacokinetic and pharmacodynamic modeling of transformed binary effect data of triflusal in healthy Korean male volunteers: a randomized, open-label, multiple dose, crossover study

Background

Triflusal is a drug that inhibits platelet aggregation. In this study we investigated the dose-exposure-response relationship of a triflusal formulation by population pharmacokinetic (PK) and pharmacodynamic (PD) modeling of its main active metabolite, hydroxy-4-(trifluoromethyl) benzoic acid (HTB).

Methods

This study was a randomized, open-label, multiple-dose, two-period, two-treatment, comparative crossover design. All volunteers received a single oral loading dose of 900 mg of triflusal on Day 1, followed by a dose of 600 mg/day from Day 2 to 9. Using data from 34 healthy volunteers, 476 HTB plasma concentration data points and 340 platelet aggregation data points were used to construct PK and PD models respectively using NONMEM (version 6.2). As the PD endpoint was qualitative, we implemented binary analysis of ‘inhibition’ and ‘non-inhibition’ rather than using the actual value of the test. The final PK-PD model was evaluated using a visual predictive check (VPC) and bootstrap.

Results

The time-concentration profile of HTB over the entire dosing period was described by a one-compartment model with a first-order formation rate constant for HTB. Weight was selected as a covariate for clearance and volume of triflusal, respectively. The structure and the population estimates for triflusal PK were as follows: oral clearance (CL/F) = 0.2 · (weight/71.65)^0.845 L/h, oral volume of distribution (V/F) = 8.3 · (weight/71.65) L, and k_f = 0.341 h^-1. A sigmoid relationship between triflusal concentration and the probability of significant inhibition with shape factor was chosen as the final PD model. No time delay between concentration and response was identified. The final structure between predicted concentration and the probability of inhibition of platelet aggregation (IPA) relationship was as follows: Probability of . Thus, we concluded this relationship is more like quantal concentration-response relationship. The current dosing regimen was considered to be efficacious based on the EC₅₀ estimate of 84.9 μg/mL obtained in this study.

Conclusions

A PK and binary probability PD model of triflusal was successfully developed for Korean healthy volunteers. The model may be used to further prediction inhibition of platelet aggregation by triflusal.

Trial registration

Clinical Research Information Service (CRIS), KCT0001299 (Registered December 5, 2014)