Effect of vitamin K2 on the anticoagulant activity of warfarin during the perioperative period of catheter ablation: Population analysis of retrospective clinical data
© The Author(s). 2016
Received: 31 March 2016
Accepted: 5 July 2016
Published: 4 August 2016
Catheter ablation is a non-medication therapy for atrial fibrillation, and during the procedure, warfarin is withdrawn in the preoperative period to prevent the risk of bleeding. In case of emergency, vitamin K2 can be intravenously administered to antagonize the anticoagulant activity of warfarin. The aims of this study were to conduct population pharmacokinetic/pharmacodynamic modeling for retrospective clinical data and to investigate the effect of vitamin K2 on the anticoagulant activity of warfarin in the perioperative period of catheter ablation.
A total of 579 international normalized ratio (INR) values of prothrombin time from 100 patients were analyzed using the nonlinear mixed-effects modeling program NONMEM. A 1-compartment model was adapted to the pharmacokinetics of warfarin and vitamin K2, and the indirect response model was used to investigate the relationship between plasma concentration and the pharmacodynamic response of warfarin and vitamin K2. Since no plasma concentration data for warfarin and vitamin K2 were available, 3 literally available pharmacokinetic parameters were used to simultaneously estimate 1 pharmacokinetic parameter and 5 pharmacodynamic parameters.
The population parameters obtained not only successfully explained the observed INR values, but also indicated an increase in sensitivity to warfarin in patients with reduced renal function. Simulations using these parameters indicated that vitamin K2 administration of more than 20 mg caused a slight dose-dependent decrease in INR on the day of catheter ablation and a delayed INR elevation after warfarin re-initiation.
A pharmacokinetic/pharmacodynamic model was successfully built to explain the retrospective INR data during catheter ablation. Simulation studies suggest that vitamin K2 should be administered with care and that more than 20 mg is unnecessary in the preoperative period of catheter ablation.
KeywordsWarfarin Vitamin K2 Pharmacodynamics INR NONMEM Catheter ablation
Atrial fibrillation is the most common sustained cardiac arrhythmia and a major cause of stroke [1, 2]. In order to prevent stroke, an anticoagulant drug, warfarin, is usually used since aspirin was proven ineffective in retrospective analyses . The anticoagulant effect of warfarin does not always correlate with its dose, and polymorphisms in cytochrome P450 (CYP) 2C9 and vitamin K epoxide reductase complex subunit 1 (VKORC1) genes have been proven to influence interindividual variability in the optimal doses, in addition to patients’ primary diseases and characteristics such as age or ethnicity [4, 5]. In Japanese patients, warfarin dose adjustments based on their prothrombin time, an international normalized ratio (INR) of 1.6–2.6 (age ≥ 70 years) or 2.0-3.0 (age < 70 years), are recommended for effective therapy to avoid life-threatening bleeding [6, 7]. When hemorrhagic complications occur, warfarin withdrawal is required and vitamin K2 or fresh frozen plasma administration is recommended [8–10].
In atrial fibrillation treatment, antiarrythmic agents are often used, while catheter ablation is also an available option as a non-medication therapy . When catheter ablation, an invasive procedure for complete cure of atrial fibrillation, is selected, anticoagulant therapy with warfarin is withdrawn in the preoperative period to prevent the risk of bleeding, although catheter ablation is sometimes performed in periprocedural therapeutic anticoagulation with warfarin if possible. In some patients, discontinuation of warfarin is not sufficient to lower the INR to the required level before catheter ablation. In such cases, vitamin K2 is intravenously administered to antagonize the anticoagulant activity of warfarin resulting in prompt recovery of INR to a safe level. Some reports have mentioned the use of pharmacokinetic/pharmacodynamic models for an anticoagulant drug and have conducted population analyses; however, only warfarin was investigated using these models [11, 12]. The effect of vitamin K2 dose on controlling the anticoagulant activity of warfarin during the perioperative period of catheter ablation has not yet been reported. The aims of this study are to build a population pharmacokinetic/pharmacodynamic model not only for warfarin, but also for vitamin K2, by using routine clinical data of patients who had been diagnosed with atrial fibrillation and received a catheter ablation, and to obtain information on the optimal vitamin K2 dose in the preoperative period before catheter ablation.
Patients and data studied
We retrospectively collected data from patients who have had a catheter ablation for atrial fibrillation at the Department of Cardiovascular Medicine, Kyoto University Hospital from January to December in 2008. During this period, 126 Japanese patients underwent catheter ablation, and 111 of these patients were treated with warfarin on the day of admission. A total of 100 patients whose INR values were between 1.0 and 3.0 in the hospitalization period were included in this study. We used 579 INR values obtained from 100 patients during the perioperative period. Clinical laboratory data and medication history for the patients studied were collected from electrical medical records. No patients were taking any medications that may have clinically significantly altered the pharmacokinetics of warfarin, except 4 patients with amiodaron and 1 patient with bucolome [13, 14].
Pharmacokinetic/pharmacodynamic model building
where Cp 1 and Cp 3 represent the plasma concentration of warfarin and vitamin K2, respectively; and Vd 1 and Vd 3 represent the distribution volume; and k 10 and k 30 represent the elimination rate constant for each drug, respectively. Since no plasma concentration data were available for warfarin and vitamin K2, and INR values were the available data for this study, reported pharmacokinetic parameters for warfarin in Japanese patients  and the distribution volume for vitamin K2 in the product information (Eisai Co., Ltd., Tokyo, Japan) were used in the analysis: k 10 = 0.0129 (1/h), Vd 1 = 0.183 (L/kg) and Vd 3 = 0.051 (L/kg). Therefore, k 30 was the only pharmacokinetic parameter to be estimated in this analysis.
Since the predicted values were outputted by the nonlinear mixed-effects modeling program (NONMEM)  using TT values, these were then converted into INR values when necessary by solving the quadratic equation obtained from Equation 4.
where P ij is the i-th individual pharmacokinetic or pharmacodynamic parameter for patient j; P pop,i is the i-th population mean parameter; and η ij is the individual random perturbation from the population mean parameter that is distributed with a mean of zero and variance ωι 2 . TT jk is the observed TT value at time k for patient j; TT * jk is the corresponding predicted TT value; and ε jk represents the independent identically distributed error with a mean of zero and variance of σ 2 for the TT value.
where OBJ is the objective function values calculated using the NONMEM and M is the number of independently adjusted parameters within the model.
where RF = 1 if serum creatinine was higher than our in-hospital reference value, namely 1.1 mg/dL or higher for men, and 0.8 mg/dL or higher for women, otherwise RF = 0. P * pop,i is the i-th population mean parameter in the patient whose serum creatinine is within our in-hospital reference value. The parameter set that had the smallest objective function value was selected, and the null hypothesis that θ was not statistically different from unity was examined using the likelihood ratio test. A difference of 7.88 in OBJ with 1 degree of freedom was used to measure statistical significance (P < 0.005 by the chi-squared distribution).
Simulation for INR transition
(A) Effect of vitamin K2 dose
Simulations were carried out using the obtained population mean parameters based on a typical patient whose body weight was 50 kg with/without renal failure. The maintenance dose of warfarin was set to 3 mg/day (7 PM) and was stopped on day −1 (the day prior to the operation), and 5 mg/day was administered for 2 days after the operation as a loading dose, followed by a maintenance dose of 3 mg/day. Vitamin K2 was administered at 20 mg 0, 1, 2, or 3-times every 4 hours after 4 PM on day −1 with the total dose administered ranging from 0 mg to 60 mg.
For quantitative evaluation, we obtained 4 parameters, namely ∆INR, 1st loading, 95 % recovery, and INR/day. The ∆INR represents the difference in INR values between before warfarin withdrawal and before the loading dose; the 1st loading represents an INR increase after the first warfarin loading dose; and the 95 % recovery represents the time needed for INR elevation in the postoperative period up to 95 % of the preoperative steady state INR value. In addition, INR/day was calculated by dividing ∆INR by 95 % recovery (day).
(B) Effect of warfarin dose
Simulations with various warfarin maintenance doses were conducted. As a maintenance dose, 3 to 6 mg of warfarin was administered and it was stopped on day −1 without vitamin K2 administration. Warfarin (2 mg) was added to each maintenance dose as a loading dose, and it was administered for 2 days after the operation, followed by each maintenance dose. Cases where 20 to 60 mg of vitamin K2 was administered were also simulated.
(C) Effect of interindividual variability
Simulations were also conducted using several parameter sets in which 1 of the mean parameters was altered using the interindividual variability (+ or – ω) from the population mean value. Warfarin and vitamin K2 doses were set to 3 and 20 mg, respectively, in each simulation.
Patients’ characteristics and INR transitions
Number or median (min-max)
Total number of patients (M/F)
Body weight (kg)
Warfarin maintenance dose (mg)
Warfarin loading dose (mg)
Number of patients treated with vitamin K2
Total dosage of vitamin K2 (mg)
Total bilirubin concentration (mg/dL)
Serum albumin (g/dL)
Serum creatinine concentration (mg/dL)
Estimated glomerular filtration rate (mL/min/1.73 m2)
When interindividual variability was considered for all population pharmacokinetic/pharmacodynamic mean parameters (η = 6), AIC was 3398. To simplify the model in which only η ks and η IC50 were included (η = 2), AIC was 3394, and was decreased to 3393 when another η for k 30 was included in the model (η = 3). Thus, the model with the minimum AIC value was adopted, which reflected the interindividual variability of k s , IC 50, and k 30 .
There were only 4 patients out of 100 patients whose total bilirubin concentration exceeded our in-hospital reference value, and those values were not remarkably high. Therefore, the effect of hepatic function on population mean parameters was not further examined. The anticoagulant effect of warfarin is generally considered to be associated with its unbound plasma concentration . We examined the effect of serum albumin concentration on the IC 50 or k 10 , but we could not obtain any significant effects.
Final population pharmacokinetic and pharmacodynamic parameters
k s (%/h)
k d (1/h)
IC 50 (μg/mL)
EC 50 (μg/mL)
k 30 (1/h)
ω ks 2
ω IC50 2
ω k30 2
Residual variability (%)
Validity of population mean parameters
Effect of renal function on INR transition
(A) Effect of renal function on INR transitions.
Decreased renal function
Vitamin K2 (mg)
1st Loading (×10−1)
95 % Recovery (h)
(B) Effects of combinations of various warfarin maintenance doses and vitamin K2 doses on INR transitions.
Vitamin K2 (mg)
1st Loading (×10−1)
95 % Recovery (h)
Vitamin K2 (mg)
1st Loading (×10−1)
95 % Recovery (h)
(C) Effect of interindividual variability on INR transitions.
1st Loading (×10−1)
95 % Recovery (h)
Effect of warfarin dose on INR transition
Effect of interindividual variability on INR transition
Figure 6b shows the effects of interindividual variability on INR transition. The simulated curves suggested that the interindividual variability of k 30 had a relatively small effect on INR variability, while k s and IC 50 had greater effects although they varied by 26.5 % or 37.9 %, respectively, from each population mean value. The INR values under a warfarin maintenance dose of 3 mg ranged from 1.47 to 1.98, and INR values after warfarin withdrawal ranged from 1.23 to 1.55, depending on k s , IC 50, and k 30 values. Table 3C shows quantitative indices of the results of Fig. 6b. The ∆INR values ranged from 0.237 to 0.416, from 0.220 to 0.504, and from 0.294 to 0.310, when k s , IC 50 , and k 30 were increased or decreased by the interindividual variability from the population mean value, respectively. The interindividual variability of k s , IC 50 , and k 30 had similar effects on the 1st loading. Unlike the ∆INR values, the interindividual variability of k s and IC 50 had a small effect on the 95 % recovery, while the k 30 value strongly affected the 95 % recovery.
It is widely known that the warfarin dose suitable for a patient varies among individuals and that careful monitoring of its anticoagulant activity is necessary for preventing excessive anticoagulation or hemorrhagic events [6, 7]. Vitamin K2 can effectively antagonize warfarin, for example, in the preoperative period and when life-threatening bleeding occurs . Although the recommended dose of vitamin K2 was under 5 mg , 20–70 mg of vitamin K2 was administered to decrease the INR value in the preoperative period (Table 1). Thus, caution must be exercised to find a balance between over- and under-coagulation. The pharmacokinetics and pharmacodynamics of warfarin have been studied since 1960’s [11, 12, 14, 21, 22], while combined pharmacokinetic/pharmacodynamic analyses of both warfarin and vitamin K formulations have not yet been reported. In the present study, we built a model that describes the pharmacokinetics/pharmacodynamics of these drugs for the first time. However, because this is a retrospective study wherein only patients’ pharmacodynamic data were used and because we converted the INR values to TT values while calculating the pharmacokinetic/pharmacodynamic parameters, special attention should be paid when drawing conclusions from the results obtained herein. Additionally, the obtained pharmacokinetic and pharmacodynamic parameters should be carefully treated, since these values greatly depended on the fixed pharmacokinetic parameters of warfarin and vitamin K2 in the model.
Final population pharmacokinetic/pharmacodynamic parameters had reasonably small relative standard errors except ω k30 2 (Table 3), and both individual predicted TT and INR values were well correlated with the observed values (Figs. 3 and 4), indicating that reliable population mean parameters were obtained in this study. Some patients had the INR values between 1.0-1.5 on the day of admission (Fig. 2). We could not check drug compliance in the patients before the hospitalization, but good compliance was expected in the hospital. Since the prediction bias of the TT was not observed against the time (data not shown), effects of non-compliance on the present results were considered to be small. Since coadministration of amiodarone or bucolome was reported to inhibit the warfarin metabolism mediated by CYP2C9 [13, 14], we examined the effect of these drugs on the k 10 . Although the coadministration of these drugs decreased k 10 , this effect did not reach a statistical significance level (−2LLD = 7.61 < 7.88). Therefore, we did not include the effect of amiodarone and bucolome in the final model. The estimated population mean parameters for k s , k d, and IC 50 were similar to those in a previous report , and interindividual variability for k s , IC 50 , and k 30 was minimal, although the intraindividual variability was quite significant.
The several simulations of INR transition by the obtained population pharmacokinetic/pharmacodynamic parameters showed that vitamin K2 could antagonize the anticoagulant activity of warfarin in a dose-dependent manner. While more than 20 mg of vitamin K2 showed only a small effect on the extent of INR decreases in the preoperative period, the time required for warfarin to exert its anticoagulation activity again in the postoperative period depended on the total dose of vitamin K2. An inability to anticoagulate promptly after the operation may lead to prolonged hospitalization and consequently decrease patients’ quality of life, as well as increase medical costs. Although it is important to examine the effect of less than 20 mg vitamin K2 on INR, we could not obtain clinical data using less than 20 mg vitamin K2. Effects of lower dose of vitamin K2 on INR remains to be examined in a future study.
In this study, we clarified the enhanced anticoagulant activity of warfarin in patients with decreased renal function. Warfarin is well known to inhibit the vitamin K-dependent synthesis pathway of coagulation factors in the liver and to be degraded in the liver . Thus, great caution is required while using warfarin in patients with hepatic disorders . According to the package insert of warfarin, caution is also required while use in those with renal dysfunction. Recent studies reported that renal function influences warfarin responsiveness and hemorrhagic complications [23, 24]. The maintenance warfarin dose was positively correlated with kidney function in Japanese patients . Precise mechanisms for the enhanced sensitivity to warfarin in patients with decreased renal function should be investigated further in future studies.
We built and analyzed a pharmacokinetic/pharmacodynamic model of both warfarin and vitamin K2 by using retrospective clinical data during the catheter ablation. Simulations using the obtained population pharmacokinetic/pharmacodynamic parameters indicated that vitamin K2 should be administered with care and that more than 20 mg is unnecessary in the preoperative period of catheter ablation. Low-dose (5 mg or less) of vitamin K is recommended in the guideline .
AIC, Akaike information criterion; Cp, plasma concentration; CYP, cytochrome P450; EC 50 , 50 % effective concentration; E max , maximum effect; IC 50 , 50 % inhibitory concentration; INR, international normalized ratio; k d , degradation rate constant; k s , synthesis rate constant; LLD, log likelihood difference; OBJ, objective function; TT, thrombotest; Vd, distribution volume, k, elimination rate constant; VKORC1, vitamin K epoxide reductase complex subunit 1.
The authors would like to thank Ms. Mio Kikuchi for her kind assistance of data entry and helpful suggestions.
Availability of data and materials
The data will not be shared because of human data.
ZZ, IY, SO, and SS conceived the study, designed the protocol. ZZ, IY, SO, and YM carried out the study and drafted the manuscript. SS, MH, TK, AA, KI, and KM participated in interpretation of the data and contributed the discussions. All authors read and approved the final manuscript.
The present address of IY is the Department of Pharmacy, Kobe University Hospital.
The authors have no competing interests to declare for this study.
Consent for publication
Ethics approval and consent to participate
This study was conducted in accordance with the Declaration of Helsinki and its amendments. The study protocol was approved by the Ethics Committee of the Kyoto University Graduate School of Medicine and Kyoto University Hospital (R0264).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Haïssaguerre M, Jaïs P, Shah D, Takahashi A, Hocini M, Quiniou G, Garrigue S, Le Mouroux A, Le Métayer P, Clémenty J. Spontaneous initiation of atrial fibrillation by ectopic beats originating in the pulmonary veins. N Engl J Med. 1998;339:659–66.View ArticlePubMedGoogle Scholar
- Wellens H. Atrial fibrillation--the last big hurdle in treating supraventricular tachycardia. N Engl J Med. 1994;331:944–5.View ArticlePubMedGoogle Scholar
- trial Fibrillation Investigators: Atrial Fibrillation, Aspirin, Anticoagulation Study; Boston Area Anticoagulation Trial for Atrial Fibrillation Study; Canadian Atrial Fibrillation Anticoagulation Study; Stroke Prevention in Atrial Fibrillation Study; Veterans Affairs Stroke Prevention in Nonrheumatic Atrial Fibrillation Study. Risk factors for stroke and efficacy of antithrombotic therapy in atrial fibrillation: Analysis of pooled data from five randomized controlled trials. Arch Intern Med. 1994;154:1449–57.View ArticleGoogle Scholar
- Higashi M, Veenstra D, Kondo L, Wittkowsky A, Srinouanprachanh S, Farin F, Rettie A. Association between CYP2C9 genetic variants and anticoagulation-related outcomes during warfarin therapy. J Am Med Assoc. 2002;287:1690–8.View ArticleGoogle Scholar
- Obayashi K, Nakamura K, Kawana J, Ogata H, Hanada K, Kurabayashi M, Hasegawa A, Yamamoto K, Horiuchi R. VKORC1 gene variations are the major contributors of variation in warfarin dose in Japanese patients. Clin Pharmacol Ther. 2006;80:169–78.View ArticlePubMedGoogle Scholar
- Yasaka M, Minematsu K, Yamaguchi T. Optimal intensity of international normalized ratio in warfarin therapy for secondary prevention of stroke in patients with non-valvular atrial fibrillation. Intern Med. 2001;40:1183–8.View ArticlePubMedGoogle Scholar
- Joint Working Groups. The Japanese Circulation Society, The Japanese College of Cardiology, The Japanese Society of Electrocardiology, and The Japanese Heart Rhythm Society. Guidelines for pharmacotherapy of atrial fibrillation (JCS2013). 2013. http://www.j-circ.or.jp/guideline/pdf/JCS2013_inoue_h.pdf. Accessed 30 March 2016.
- Joint working Groups. Guidelines for management of anticoagulant and antiplatelet therapy in cardiovascular disease (JCS2004). Circ J. 2004;68 Suppl 4:1153–219.Google Scholar
- The Haemostasis and Thrombosis Task Force for the British Committee for Standards in Haematology. Guidelines on oral anticoagulation: third edition. Br J Haematol. 1998;101:374–87.Google Scholar
- Shetty H, Fennerty A, Routledge P. Clinical pharmacokinetic considerations in the control of oral anticoagulant therapy. Clin Pharmacokinet. 1989;16:238–53.View ArticlePubMedGoogle Scholar
- Sato S, Toda T, Yamazaki M, Hongo F, Kurosawa N, Owada E, Hatta E, Nakamura M, Nakanishi K, Okamoto F, Sakai K, Hirokami M, Hanawa N, Tanaka S. Warfarin maintenance dose estimation program “WfTDM” based on thrombo-test value. Jpn J Ther Drug Monit. 2006;23:10–6.Google Scholar
- Hamberg AK, Dahl ML, Barban M, Scordo MG, Wadelius M, Pengo V, Padrini R, Jonsson EN. A PK-PD model for predicting the impact of age, CYP2C9, and VKORC1 genotype on individualization of warfarin therapy. Clin Pharmacol Ther. 2007;81:529–38.View ArticlePubMedGoogle Scholar
- Holbrook AM, Pereira JA, Labiris R, McDonald H, Douketis JD, Crowther M, Wells PS. Systematic overview of warfarin and its drug and food interactions. Arch Intern Med. 2005;165:1095–106.View ArticlePubMedGoogle Scholar
- Osawa M, Hada N, Matsumoto K, Hasegawa T, Kobayashi D, Morimoto Y, Yamaguchi M, Kanamoto I, Nakagawa T, Sugibayashi K. Usefulness of coadministration of bucolome in warfarin therapy: pharmacokinetic and pharmacodynamic analysis using outpatient prescriptions. Int J Pharm. 2005;293:43–9.View ArticlePubMedGoogle Scholar
- Jusko WJ, Ko HC. Physiological indirect response models characterize diverse types of pharmacodynamic effects. Clin Pharmacol Ther. 1994;56:406–19.View ArticlePubMedGoogle Scholar
- Dayneka NL, Garg V, Jusko WJ. Comparison of 4 basic models of indirect pharmacodynamic responses. J Pharmacokinet Biopharm. 1993;21:457–78.View ArticlePubMedPubMed CentralGoogle Scholar
- Tie J, Stafford D. Structure and function of vitamin K epoxide reductase. Vitam Horm. 2008;78:103–30.View ArticlePubMedGoogle Scholar
- Gogstad G, Wadt J, Smith P, Brynildsrud T. Utility of a modified calibration model for reliable conversion of thromboplastin times to international normalized ratios. Thromb Haemost. 1986;56:178–82.PubMedGoogle Scholar
- Beal SL, Boeckmann AJ, Sheiner LB. NONMEM Users Guides. NONMEM Project Group. San Francisco: University of California; 1992.Google Scholar
- Akaike H. A new look at the statistical model identification. IEEE Trans Autom Control. 1974;19:716–23.View ArticleGoogle Scholar
- Sheiner LB. Computer-aided long-term anticoagulation therapy. Comput Biomed Res. 1969;2:507–18.View ArticlePubMedGoogle Scholar
- Nagashima RA, O’reilly RA, Levy G. Kinetics of pharmacologic effects in man: the anticoagulant action of warfarin. Clin Pharmacol Ther. 1969;10:22–35.View ArticlePubMedGoogle Scholar
- Limdi NA, Beasley TM, Baird MF, Goldstein JA, McGwin G, Arnett DK, Acton RT, Allon M. Kidney function influences warfarin responsiveness and hemorrhagic complications. J Am Soc Nephrol. 2009;20:912–21.View ArticlePubMedPubMed CentralGoogle Scholar
- Limdi NA, Limdi MA, Cavallari L, Anderson AM, Crowley MR, Baird MF, Allon M, Beasley TM. Warfarin dosing in patients with impaired kidney function. Am J Kidney Dis. 2010;56:823–31.View ArticlePubMedPubMed CentralGoogle Scholar
- Ichihara N, Ishigami T, Umemura S. Effect of impaired renal function on the maintenance dose of warfarin in Japanese patients. J Cardiol. 2015;65:178–84.View ArticlePubMedGoogle Scholar