A Physiological-Based Pharmacokinetic Model Embedded with a Target-Mediated Drug Disposition Mechanism Can Characterize Single-Dose Warfarin Pharmacokinetic Profiles in Subjects with Various CYP2C9 Genotypes under Different Cotreatments

Warfarin-A natural anticoagulant: A review of research trends for precision medication

BACKGROUND

Warfarin is a widely prescribed anticoagulant in the clinic. It has a more considerable individual variability, and many factors affect its variability. Mathematical models can quantify the quantitative impact of these factors on individual variability.

PURPOSE

The aim is to comprehensively analyze the advanced warfarin dosing algorithm based on pharmacometrics and machine learning models of personalized warfarin dosage.

METHODS

A bibliometric analysis of the literature retrieved from PubMed and Scopus was performed using VOSviewer. The relevant literature that reported the precise dosage of warfarin calculation was retrieved from the database. The multiple linear regression (MLR) algorithm was excluded because a recent systematic review that mainly reviewed this algorithm has been reported. The following terms of quantitative systems pharmacology, mechanistic model, physiologically based pharmacokinetic model, artificial intelligence, machine learning, pharmacokinetic, pharmacodynamic, pharmacokinetics, pharmacodynamics, and warfarin were added as MeSH Terms or appearing in Title/Abstract into query box of PubMed, then humans and English as filter were added to retrieve the literature.

RESULTS

Bibliometric analysis revealed important co-occuring MeShH and index keywords. Further, the United States, China, and the United Kingdom were among the top countries contributing in this domain. Some studies have established personalized warfarin dosage models using pharmacometrics and machine learning-based algorithms. There were 54 related studies, including 14 pharmacometric models, 31 artificial intelligence models, and 9 model evaluations. Each model has its advantages and disadvantages. The pharmacometric model contains biological or pharmacological mechanisms in structure. The process of pharmacometric model development is very time- and labor-intensive. Machine learning is a purely data-driven approach; its parameters are more mathematical and have less biological interpretation. However, it is faster, more efficient, and less time-consuming. Most published models of machine learning algorithms were established based on cross-sectional data sourced from the database.

CONCLUSION

Future research on personalized warfarin medication should focus on combining the advantages of machine learning and pharmacometrics algorithms to establish a more robust warfarin dosage algorithm. Randomized controlled trials should be performed to evaluate the established algorithm of warfarin dosage. Moreover, a more user-friendly and accessible warfarin precision medicine platform should be developed.

A physiologically-based pharmacokinetic/pharmacodynamic modeling approach for drug-drug-gene interaction evaluation of S-warfarin with fluconazole

Warfarin is a widely used anticoagulant, and its S-enantiomer has higher potency compared to the R-enantiomer. S-warfarin is mainly metabolized by cytochrome P450 (CYP) 2C9, and its pharmacological target is vitamin K epoxide reductase complex subunit 1 (VKORC1). Both CYP2C9 and VKORC1 have genetic polymorphisms, leading to large variations in the pharmacokinetics (PKs) and pharmacodynamics (PDs) of warfarin in the population. This makes dosage management of warfarin difficult, especially in the case of drug-drug interactions (DDIs). This study provides a whole-body physiologically-based pharmacokinetic/PD (PBPK/PD) model of S-warfarin for predicting the effects of drug-drug-gene interactions on S-warfarin PKs and PDs. The PBPK/PD model of S-warfarin was developed in PK-Sim and MoBi. Drug-dependent parameters were obtained from the literature or optimized. Of the 34 S-warfarin plasma concentration-time profiles used, 96% predicted plasma concentrations within twofold range compared to observed data. For S-warfarin plasma concentration-time profiles with CYP2C9 genotype, 364 of 386 predicted plasma concentration values (~94%) fell within the twofold of the observed values. This model was tested in DDI predictions with fluconazole as CYP2C9 perpetrators, with all predicted DDI area under the plasma concentration-time curve to the last measurable timepoint (AUClast) ratio within twofold of the observed values. The anticoagulant effect of S-warfarin was described using an indirect response model, with all predicted international normalized ratio (INR) within twofold of the observed values. This model also incorporates a dose-adjustment method that can be used for dose adjustment and predict INR when warfarin is used in combination with CYP2C9 perpetrators.

Effect of High Altitude Environment on Pharmacokinetic and Pharmacodynamic of Warfarin in Rats

BACKGROUND

High altitude environment affects the pharmacokinetic (PK) parameters of drugs and the PK parameters are an important theoretical basis for guiding the rational clinical use of drugs. Warfarin is an oral anticoagulant of the coumarin class commonly used in clinical practice, but it has a narrow therapeutic window and wide individual variation. However, the effect of high altitude environment on PK and pharmacodynamic (PD) of warfarin is unclear.

OBJECTIVE

The objective of this study is to investigate the effect of a high altitude environment on PK and PD of warfarin in rats.

METHOD

Rats were randomly divided into plain group and high altitude group and blood samples were collected through the orbital venous plexus after administration of 2 mg/kg warfarin. Warfarin concentrations in plasma samples were determined by liquid chromatography-tandem mass spectrometry (LC-MS/MS) and PK parameters were calculated by the non-compartment model using WinNonlin 8.1 software. Meanwhile, the expression of PXR, P-gp and CYP2C9 in liver tissues was also determined by western blotting. The effect of high altitude environment on PD of warfarin was explored by measuring activated partial thromboplastin time (APTT) and prothrombin time (PT) values and then calculated international normalized ratio (INR) values based on PT.

RESULTS

Significant changes in PK behaviors and PD of warfarin in high altitude-rats were observed. Compared with the plain-rats, the peak concentration (Cmax) and the area under the plasma concentration-time curve (AUC) increased significantly by 50.9% and 107.46%, respectively. At the same time, high altitude environment significantly inhibited the expression of PXR, P-gp and CYP2C9 in liver tissues. The results of the PD study showed that high altitude environments significantly prolonged PT, APTT and INR values.

CONCLUSION

High altitude environment inhibited the metabolism and increased the absorption of warfarin in rats and increased the effect of anticoagulant effect, suggesting that the optimal dose of warfarin for patients at high altitude should be reassessed.

Predicting In Vivo Target Occupancy (TO) Profiles via Physiologically Based Pharmacokinetic-TO Modeling of Warfarin Pharmacokinetics in Blood: Importance of Low Dose Data and Prediction of Stereoselective Target Interactions

Warfarin is well recognized for its high-affinity and capacity-limited binding to the pharmacological target and undergoes target-mediated drug disposition. Here, we developed a physiologically based pharmacokinetic (PBPK) model that incorporated saturable target binding and other reported hepatic disposition components of warfarin. The PBPK model parameters were optimized by fitting to the reported blood pharmacokinetic (PK) profiles of warfarin with no stereoisomeric separation after oral dosing of racemic warfarin (0.1, 2, 5, or 10 mg) using the Cluster Gauss-Newton method (CGNM). The CGNM-based analysis yielded multiple "accepted" sets for six optimized parameters, which were then used to simulate the warfarin blood PK and in vivo target occupancy (TO) profiles. When further analyses examined the impact of dose selection on uncertainty in parameter estimation by the PBPK modeling, the PK data from 0.1 mg dose (well below target saturation) was important in practically identifying the target binding-related parameters in vivo. When stereoselective differences were incorporated for both hepatic disposition and target interactions, our PBPK modeling predicted that R-warfarin (of slower clearance and lower target affinity than S-warfarin) contributes to TO prolongation after oral dosing of racemic warfarin. Our results extend the validity of the approach by which the PBPK-TO modeling of blood PK profiles can yield TO prediction in vivo (applicable to the drugs with targets of high affinity and abundance and limited distribution volume via nontarget interactions). Our findings support that model-informed dose selection and PBPK-TO modeling may aid in TO and efficacy assessment in preclinical and clinical phase 1 studies. SIGNIFICANCE STATEMENT: The current physiologically based pharmacokinetic modeling incorporated the reported hepatic disposition components and target binding of warfarin and analyzed the blood pharmacokinetic (PK) profiles from varying warfarin doses, practically identifying target binding-related parameters in vivo. By implementing the stereoselective differences between R- and S-warfarin, our analysis predicted the role of R-warfarin in prolonging overall target occupancy. Our results extend the validity of analyzing blood PK profiles to predict target occupancy in vivo, which may guide efficacy assessment.