2022 Volume 45 Issue 7 Pages 948-954
Some population pharmacokinetic models for amiodarone (AMD) did not incorporate N-desethylamiodarone (DEA) concentration. Glucocorticoids activate CYP3A4 activity, metabolizing AMD. In contrast, CYP3A4 activity may decrease under inflammation conditions. However, direct evidence for the role of glucocorticoid or inflammation on the pharmacokinetics of AMD and DEA is lacking. The pilot study aimed to address this gap using a population pharmacokinetic analysis of AMD and DEA. A retrospective cohort observational study in adult patients who underwent AMD treatment with trough concentration measurement was conducted at Tokyo Women’s Medical University, Medical Center East from June 2015 to March 2019. Both structural models of AMD and DEA applied 1-compartment models, which included significant covariates using a stepwise forward selection and backward elimination method. The eligible 81 patients (C-reactive protein level: 0.26 [interquartile range; 0.09–1.92] mg/dL) had a total of 408 trough concentrations for both AMD and DEA. The median trough concentrations were 0.49 [0.31–0.81] µg/mL for AMD and 0.43 [0.28–0.71] µg/mL for DEA during a median follow-up period of 446 [147–1059] d. Three patients received low-dose oral glucocorticoid. The final model identified that AMD clearance was 7.9 L/h, and the apparent DEA clearance was 10.3 L/h. Co-administered glucocorticoids lowered apparent DEA clearance by 35%. These results indicate that co-administered glucocorticoids may increase DEA concentrations in patients without severe inflammation.
Amiodarone (AMD) and N-desethylamiodarone (DEA), an active metabolite of AMD, block the potassium channels responsible for cardiac repolarization, thereby suppressing atrial fibrillation and ventricular tachycardia.1,2) AMD is effective to control arrhythmia (e.g., atrial fibrillation and ventricular/ventricular fibrillation).3,4) The pharmacokinetics of AMD depends on age and body size.5–7) A trough concentration of AMD >0.5 µg/mL is required to obtain an antiarrhythmic effect.8) Clinicians recently have prescribed low-dose AMD targeting lower concentrations.9) AMD frequently causes adverse events in patients with high AMD trough concentrations.10) For instance, DEA concentration >0.6 µg/mL is a risk factor for interstitial pneumonia.11)
Previous population pharmacokinetic analyses focused on covariates for AMD, although these have not examined covariates of DEA that cause lung damage.5–7,12) Because AMD is metabolized primarily via CYP3A4, glucocorticoids induce CYP3A4 producing DEA from AMD.13,14) However, drug-drug interaction between AMD (or DEA) and glucocorticoids remains unclear. Furthermore, there is a concern that inflammation inhibits CYP3A4 activity and increases exposure to AMD as with midazolam and perampanel, which positively correlates with C-reactive protein (CRP) levels.15,16)
Accordingly, this pilot study aimed to conduct a comprehensive assessment of glucocorticoid and inflammation using a population pharmacokinetic analysis of AMD and DEA.
A retrospective cohort observational study was conducted at Tokyo Women’s Medical University, Medical Center East (450 bed, branch hospital of academic medical center) from June 2015 to March 2019. We screened study patients who underwent AMD treatment orally or intravenously with monitored blood concentrations of AMD and DEA. Patients were excluded if they were aged <20 years or had not undergone CRP measurement. The study protocol was approved by the institutional review board of Tokyo Women’s Medical University and followed the Declaration of Helsinki (Approval No. 5151).
Data CollectionWe extracted clinical data using electronic medical records from the commencement of AMD treatment to March 2019 or the discontinuation of AMD. All blood samples for trough concentrations were taken before the next administration of AMD. AMD and DEA trough concentrations were quantified by performing LC coupled with ion trap mass spectrometry. The limit of quantification was 20.0 ng/mL for AMD and 50.0 ng/mL for DEA.
Population Pharmacokinetic ModelingPharmacokinetic models for AMD and DEA were constructed using non-linear mixed effect analysis with Phoenix® NLME Version 8.1 software (Certara USA, Inc., Princeton, NJ, U.S.A.). The Laplace approximation computational algorithm was used. A model including both AMD and DEA was simultaneously used and sequentially fitted to the data, that is, AMD parameters were first estimated, then, fixing those parameters estimated, DEA parameters were estimated. Since only AMD and DEA trough concentrations were available, we could not construct a multi-compartment model and estimate the absorption rate constant (ka) of AMD. Therefore, a pharmacokinetic model for AMD was constructed as a 1-compartment model with a first-order absorption process. Oral bioavailability (Foral) was fixed as 54%.17) As ka data was not available, we arbitrarily considered ka as 1.0 h−1. The structural model was constructed considering amiodaron clearance (CLAMD) and volume distribution (VdAMD) as the pharmacokinetic parameters; VdAMD was fixed as 14.0 L/kg because we could not obtain peak concentration data.17) Similarly, the structural model for DEA was constructed as a 1-compartment model with clearance (CLDEA) and volume distribution (VdDEA) in addition to a metabolic fraction of AMD to DEA (fm). Since we could not obtain detailed DEA pharmacokinetics data, we assumed that VdDEA was compatible with VdAMD and fixed VdDEA as 14.0 L/kg following AMD data. We tested an additive or a proportional error model to describe residual variability, defined as follows:
![]() |
where Cobs and Cpred denote observed and predicted trough concentration (AMD or DEA), and ε denotes measurement error including intraindividual variability, analytical error, and dosing error.
Interindividual variability of AMD and DEA clearance was described using an exponential random-effects model defined as follows:
![]() |
where CLAMD and CLDEA denote the parameter for AMD and DEA clearance, tv CLAMD, and tv CLDEA are the typical values for CLAMD and CLDEA, respectively, and each η represents interindividual variability. We explored covariates for CLAMD and CLDEA. The continuous covariate was modelled as follows:
![]() |
where Ptv denotes the typical value of the parameter, θtv denotes the parameter for the covariate equal median value, and θ denotes a covariate scale factor for the covariate. In addition, the categorical covariate was modelled as follows:
![]() |
where n is an indicator of the covariate (0 or 1).
The candidate covariates comprised sex, age, body mass index, cardiac comorbidity, CRP, and glucocorticoids. Time-varying CRP values were used in the covariate modeling. The last observation carried forward was used to impute missing values for clinical laboratory data. We checked the presence of multicollinearity of covariates. The stepwise forward selection and backward elimination method determined covariates according to the objective function value (OFV). The significance level for forward selection and backward elimination was set at p < 0.05 (ΔOFV 3.84) and p < 0.01 (ΔOFV 6.63).
Pharmacokinetic Model ValidationGoodness-of-fit plots showed the correlations of observed vs. predicted concentrations, observed vs. individual predicted concentrations, conditional weighted residuals vs. predicted concentrations, and conditional weighted residuals vs. time after dose. We visually compared observed with predicted values based on prediction corrected visual predictive check (pc-VPC) using 1000 hypothetical data sets. Using operating bootstrap methods, we obtained mean estimates with 95% confidence intervals (95% CI) of respective parameters from the generated 1000 resampling datasets. The final model was considered reliable if the parameter estimates were within the 95% CI. The results of bootstrap methods were compared with those obtained from the final model.
Pharmacokinetic Model ApplicationWe performed simulations to confirm an oral dosing schedule for AMD treatment. The initial oral dose was fixed at 200 mg twice daily for a week. Maintenance oral dose was assigned as follows: 1) 50 mg once daily; 2) 50 mg twice daily; 3) 100 mg once daily; 4) 100 mg twice daily. Body weight was fixed as 50, 70, and 90 kg. The results were summarized according to glucocorticoid use. We simulated steady-state trough concentrations at day 150 after AMD treatment.
Outcome AnalysisWe performed an exploratory analysis examining the association between trough concentration and the typical outcome of overt hypothyroidism associated with AMD using univariate logistic regression analysis. Moreover, the ratio of DEA to AMD concentration and the sum were investigated similarly. We defined hypothyroidism as a thyroid-stimulating hormone level of >10 mU/L.18)
Statistical AnalysisStatistical analysis was conducted using JMP® pro 14 (SAS Institute Inc., Cary, NC, U.S.A.). Normally distributed continuous data were expressed as the mean ± standard deviation, whereas non-normally distributed data were expressed as median values and interquartile range (IQR). Categorical data were expressed as numbers (%). When we conducted a logistic regression analysis, a two-sided p-value <0.05 was considered statistically significant.
We enrolled 266 patients (>20 years) who underwent AMD therapy and therapeutic drug monitoring at our hospital from June 2015 to March 2019. There were no patients who were younger than 20 years. We excluded 185 patients who had no CRP data. Therefore, we identified 81 patients in this study. Baseline clinical data are summarized in Table 1. The median initial daily dose of AMD [IQR] was 200 [100–400] mg. There were 10 (12%) patients who orally administered AMD every other day. The median CRP level at the baseline was 0.26 [0.09–1.92] mg/dL, which corresponded to the low range. All three patients received low-dose oral prednisolone. The median trough concentrations were 0.49 [0.31–0.81] µg/mL for AMD and 0.43 [0.28–0.71] µg/mL for DEA during a follow-up period of 446 [147–1059] days.
Variables | Data, N = 81 |
---|---|
Demographical data | |
Male, N (%) | 55 (68) |
Age, years | 64 ± 15 |
Height, cm | 164 ± 11 |
Body weight, kg | 61 ± 15 |
Body mass index, kg/m2 | 23 ± 4 |
Cardiac comorbidity | |
Congestive heart failure, N (%) | 48 (59) |
Atrial fibrillation, N (%) | 27 (33) |
VT or VF, N (%) | 39 (48) |
Detail of amiodarone exposure | |
Initial daily dose, mg | 200 [100–400] |
Oral routea), N (%) | 61 (75) |
Intravenous route, N (%) | 20 (25) |
Clinical laboratory data | |
CRP, mg/dL | 0.26 [0.09–1.92] |
Concurrent medications | |
Glucocorticoidb), N (%) | 3 (4) |
N, number; VT, ventricular tachycardia; VF, ventricular fibrillation; CRP, C-reactive protein. a) 10 patients received amiodarone every other day. b) All patients received oral prednisolone. Normally distributed continuous data are expressed as mean ± standard deviation, whereas non-normally distributed data are expressed as median values and interquartile range (IQR). Categorical data are expressed as numbers and percentages.
A total of 408 AMD and 408 DEA measurement points were available for the population pharmacokinetic modeling. The median number of measurement points was 6 [2–14] per patient. There was a total of 26 AMD and 26 DEA measurement points in patients who received prednisolone. The final model parameters are summarized in Table 2. The proportional error model was incorporated in the structured model. The stepwise forward selection method determined that co-administered glucocorticoids lowered apparent CLDEA by 35% (Table 2). We did not conduct the stepwise backward elimination method since the final models included no covariates for AMD and only glucocorticoid for DEA. The overall η-shrinkage was 21.5 and 23.5% for CLAMD and CLDEA, respectively. Additionally, the overall ε-shrinkage was 8.7 and 8.9% for AMD concentration and DEA concentration, respectively.
Parameter | Estimated mean | RSE% | Bootstrap mean | 95% CI |
---|---|---|---|---|
AMD | ||||
Population value | ||||
CLAMD, L/h | 7.9 | 10.9 | 7.7 | 7.6 to 7.8 |
VdAMD, L/kg (fixed) | 14.0 | NA | 14.0 | NA |
Foral, % (fixed) | 54.0 | NA | 54.0 | NA |
ka, /h (fixed) | 1.0 | NA | 1.0 | NA |
Interindividual variability | ||||
ωCLAMD, % | 43.3 | 23.3 | 46.9 | 46.1 to 47.8 |
Intraindividual variability | ||||
σ, % | 31.3 | 10.8 | 30.9 | 30.7 to 31.1 |
DEA | ||||
Population value | ||||
CLDEA/fm, L/h | 10.3 | 11.0 | 10.1 | 10.0 to 10.2 |
Glucocorticoids on CLDEA/fm | −0.35 | 42.5 | −0.37 | −0.38 to −0.36 |
VdDEA, L/kg (fixed) | 14.0 | NA | 14.0 | NA |
Interindividual variability | ||||
ωCLDEA, % | 39.4 | 21.1 | 42.6 | 41.7 to 43.5 |
Intraindividual variability | ||||
σ, % | 23.9 | 11.8 | 23.3 | 23.1 to 23.5 |
RSE%, relative standard error (%); 95% CI, 95% confidence interval; AMD, amiodarone; CLAMD, amiodaron clearance; VdAMD, amiodarone volume distribution; Foral, oral bioavailability; ka, absorption constant; DEA; N-desethylamiodarone, CLDEA, N-desethylamiodarone clearance; fm, fraction of metabolism; VdDEA, N-desethylamiodarone volume distribution; NA, not applicable. The final model equations were as follows. CLAMD = 7.9 L/h Vd_AMD = 14.0 × Body weight CLDEA/fm = 10.3 × (1−0.35 × [if glucocorticoid]) L/h VdDEA = 14.0 × Body weight The difference in objective function value of final model was 5.0. The overall η-shrinkage was 21.5 and 23.5% for CLAMD and CLDEA, respectively. The overall ε-shrinkage was 8.7 and 8.9% for AMD concentration and DEA concentration, respectively.
Goodness-of-fit plots are shown in Figs. 1 and 2. Although some tendency was observed, the model has shown some correlation and margin of conditional weighted residual of ±3, indicating that the final model had the acceptable capacity to predict AMD and DEA trough concentrations. The results of pc-VPC also confirmed an acceptable agreement between the observed and predicted values although there were outliers as shown in Fig. 3. Bootstrap methods were 100% successful and appropriately estimated the final model parameters with mean values of the bootstrap parameters closely equivalent to those of the model parameter estimates (Table 2). A summary of three patients who received glucocorticoid has shown in Table 3.
A: Observed concentrations and population predicted concentrations. Y-axis and X-axis represent observed concentrations and population predicted concentrations, respectively. Black circles represent each point, while the solid line represents the reference line. B: Observed concentrations and individual predicted concentrations. Y-axis and X-axis represent observed concentrations and individual predicted concentrations, respectively. Black circles represent each point, while the solid line represents the reference line. C: Conditional weighted residual and time after the first dose. Y-axis and X-axis represent conditional weighted residual and time after the first dose, respectively. The solid horizontal line represents the reference line, while black circles represent each point. D: Conditional weighted residual and predicted concentrations. Y-axis and X-axis represent conditional weighted residual and predicted concentrations, respectively. Black circles represent each point. CWRES; conditional weighted residual.
A: Observed concentrations and population predicted concentrations. Y-axis and X-axis represent observed concentrations and population predicted concentrations, respectively. Black circles represent each point, while the solid line represents the reference line. B: Observed concentrations and individual predicted concentrations. Y-axis and X-axis represent observed concentrations and individual predicted concentrations, respectively. Black circles represent each point, while the solid line represents the reference line. C: Conditional weighted residual and time after the first dose. Y-axis and X-axis represent conditional weighted residual and time after the first dose, respectively. The solid horizontal line represents the reference line, while black circles represent each point. D: Conditional weighted residual and predicted concentrations. Y-axis and X-axis represent conditional weighted residual and predicted concentrations, respectively. Black circles represent each point. CWRES; conditional weighted residual.
A: X-axis and Y-axis represent the dosing period and prediction-corrected observed and predicted values of amiodarone (logarithmic scale). Opened circles represent each point. The solid and dotted lines represent predicted median and 95% confidence interval. The lower limit of predicted concentration is below 0. B: X-axis and Y-axis represent the dosing period and prediction-corrected observed and predicted values of N-desethylamiodarone (logarithmic scale). Opened circles represent each point. The solid and dotted lines represent predicted median and 95% confidence interval. The lower limit of predicted concentration is below 0.
Sex | Age year | BW kg | BMI kg/m2 | Daily dose mg | Treatment duration day | CRP mg/dL | Mean troughAMD µg/mL | Mean troughDEA µg/mL | Measurement pointsa) | |
---|---|---|---|---|---|---|---|---|---|---|
1 | M | 74 | 46 | 19 | 200 | 683 | 0.17 | 1.08 | 0.90 | 34 |
2 | M | 88 | 64 | 24 | 200 | 390 | 1.43 | 0.98 | 1.18 | 4 |
3 | F | 43 | 45 | 18 | 50 | 1225 | 0.02 | 0.16 | 0.22 | 14 |
BW, body weight; BMI, body mass index; CRP, C-reactive protein; troughAMD, amiodarone trough concentration; troughDEA, N-desethylamiodarone trough concentration; M, male; F, female. Mean troughAMD and troughDEA values were calculated by each value during the study period. a) Measurement points were available for AMD and DEA.
The summary of simulations revealed that glucocorticoid co-administration led to approximately 50% higher steady-state DEA trough concentration regardless of body weight and dosing schedule (Table 4).
Daily dose | With glucocorticoid | Without glucocorticoid | ||
---|---|---|---|---|
AMD, µg/mL | DEA, µg/mL | AMD, µg/mL | DEA, µg/mL | |
Body weight, 50 kg | ||||
50 mg once | 0.08 | 0.10 | 0.08 | 0.07 |
50 mg twice | 0.16 | 0.20 | 0.16 | 0.13 |
100 mg once | 0.15 | 0.20 | 0.15 | 0.13 |
100 mg twice | 0.33 | 0.41 | 0.33 | 0.27 |
Body weight, 70 kg | ||||
50 mg once | 0.08 | 0.10 | 0.08 | 0.07 |
50 mg twice | 0.17 | 0.20 | 0.17 | 0.13 |
100 mg once | 0.16 | 0.20 | 0.16 | 0.13 |
100 mg twice | 0.33 | 0.41 | 0.33 | 0.25 |
Body weight, 90 kg | ||||
50 mg once | 0.08 | 0.10 | 0.08 | 0.07 |
50 mg twice | 0.17 | 0.20 | 0.17 | 0.13 |
100 mg once | 0.16 | 0.20 | 0.16 | 0.13 |
100 mg twice | 0.34 | 0.41 | 0.34 | 0.27 |
AMD; amiodarone, DEA; N-desethylamiodarone.
Out of all patients, 13 (16%) experienced overt hypothyroidism associated with AMD. The significant relationship between DEA to AMD ratio and the development of overt hypothyroidism was present in this population (odds ratio = 23.9, 95% CI: 2.19 to 261.2, p < 0.01) (Supplementary Table 1).
Although there were only 3 patients concomitantly administered glucocorticoids, we used 26 blood samples of the patients and showed that glucocorticoids lead to the accumulation of DEA concentration in patients without severe CRP elevation.
Several population pharmacokinetic models have identified significant covariates, such as age and body size, in Japanese patients.5,7) In contrast, no covariates were found in a certain clinical study with a limited sample size.19) Thus, it is controversial about covariates for AMD. Our final model parameters (e.g., half-lives of about 80 h) were similar to the values reported in a previous study.7)
Glucocorticoids penetrate cells and bind to nuclear receptors, thereby exerting a 1.5-fold higher CYP3A4 activity than the control.20) Indeed, high dose tacrolimus was required to reach the target trough concentration in patients with low-dose prednisolone.21) Thus, we considered that glucocorticoids induce the conversion from AMD to DEA although we could not conclude whether our findings resulted from a decrease in CLDEA or an increase in fm. Since Fukuchi et al. suggested that elderly patients had high trough concentrations of AMD7) (Table 3), concomitant use of glucocorticoid may have a significant impact on AMD metabolism in elderly patients. Interstitial pneumonia is one of the adverse events associated with AMD.22,23) The risk of interstitial pneumonia was reported to increase by 1.3 times with the elevation of DEA trough concentration 0.2 µg/mL independent of AMD trough concentration.11) Thus, the effect of glucocorticoids on AMD pharmacokinetics may contribute to the increased risk of interstitial pneumonia. Although we did not collect the data on the onset of interstitial pneumonia, AMD dose adjustment may be an efficient strategy to minimize the risk of interstitial pneumonia in patients who co-administered AMD and glucocorticoids.
The previous report has shown that the inhibition constant of AMD for CYP3A4 is 271.6 µM and DEA for CYP3A4 is 12.1 µM.24) The present study confirmed that the median trough concentrations were 0.49 [0.31–0.81] µg/mL (0.76 µM) for AMD and 0.43 [0.28–0.71] µg/mL (0.66 µM) for N-DEA, respectively. Therefore, AMD and DEA could not exceed the inhibition constant for CYP3A4. Accordingly, we consider that the effect of AMD and DEA on glucocorticoids might be limited.
CRP had little or no effect on AMD pharmacokinetics because no patients had elevated CRP of more than 10.0 mg/dL in our study. A significant correlation between CRP level and medications (e.g., voriconazole and tramadol) has been reported in a wide range of CRP levels.25–27) Notably, AMD and DEA have long half-lives with high interindividual variability in long-term treatment, which supports our findings (half-lives: about 80 h).28) Conversely, the half-life of CRP was relatively short (about 6 to 12 h).29) Therefore, the transient and mild inflammatory reaction has limited effects on AMD treatment.
Preliminarily, our finding indicated a significant relationship between DEA to AMD ratio and the onset of overt hypothyroidism, which has been a controversial issue.30–32) The fundamental experimental data suggest that DEA plays an essential role in exerting cytotoxicity.33) On the contrary, AMD tends to accumulate in thyroid tissue.34) These controversial data should be resolved by additional research.
This study has several limitations. First, this study was a single-center retrospective analysis, thus, unknown confounding factors were inevitable. Second, the sample size was small; therefore, we could not assess the dose–response of glucocorticoids. Additionally, no patients developed hyperthyroidism (thyroid-stimulating hormone <0.35 mU/L).18) Third, we could not investigate the role of di-N-desethylamiodarone, a metabolite of DEA, exerts hepatotoxicity.35) Finally, inflammation status was only evaluated according to the CRP levels because there was no data on the other inflammatory marker (e.g., interleukin 6).
This pilot study suggested that glucocorticoid increases DEA exposure in the absence of severe inflammation. Our population pharmacokinetic model might represent a step toward adequate dose adjustment of AMD in clinical practice. Further study is required to examine the effect of glucocorticoid on AMD pharmacokinetics.
The authors declare no conflict of interest.
This article contains supplementary materials.