Risk Factors for Voriconazole-Associated Hepatotoxicity in Patients in the Intensive Care Unit

Development of a therapeutic drug-monitoring algorithm for outpatients receiving voriconazole: A multicentre retrospective study

AIMS

Although therapeutic drug monitoring (TDM) of voriconazole is performed in outpatients to prevent treatment failure and toxicity, whether TDM should be performed in all or only selected patients remains controversial. This study evaluated the association between voriconazole trough concentrations and clinical events.

METHODS

We investigated the aggravation of clinical symptoms, incidence of hepatotoxicity and visual disturbances, change in co-medications and interaction between voriconazole and co-medications in outpatients receiving voriconazole between 2017 and 2021 in three facilities. Abnormal trough concentrations were defined as <1.0 mg/L (low group) and >4.0 mg/L (high group).

RESULTS

A total of 141 outpatients (578 concentration measurements) met the inclusion criteria (treatment, 37 patients, 131 values; prophylaxis, 104 patients, 447 values). The percentages of patients with abnormal concentrations were 29.0% and 31.5% in the treatment and prophylaxis groups, respectively. Abnormal concentrations showed 50% of the concentrations at the first measurement in both therapies. Aggravation of clinical symptoms was most frequently observed in the low treatment group (18.2%). Adverse events were most common in the high group for both therapies (treatment, hepatotoxicity 6.3%, visual disturbance 18.8%; prophylaxis, hepatotoxicity 27.9%). No differences were found in changes to co-medications and drug interactions. In the prophylaxis group, prescription duration in the presence of clinical events tended to be longer than in their absence (47.4 ± 23.4 days vs 39.7 ± 21.9 days, P = .1132).

CONCLUSIONS

We developed an algorithm based on clinical events for appropriate implementation of TDM in outpatients. However, future interventions based on this algorithm should be validated.

Hypokalemia and Hyponatremia in Adult Patients Receiving Voriconazole Therapeutic Drug Monitoring

Hypokalemia and hyponatremia are common but easily ignored adverse events in treatment with voriconazole (VCZ) that can lead to serious consequences. We intend to investigate the incidence of VCZ-induced hypokalemia and hyponatremia and their risk factors based on real-world data. A prospective study was conducted. A total of 272 patients with 414 VCZ plasma trough concentrations (C0) and VCZ N-oxide concentrations (CN) were included. The incidence of hypokalemia was 18.0% (48/266). A total of 81.2% (39/48) of patients developed hypokalemia within 14 days, whereas 56.2% (27/48) of patients developed hypokalemia within 1 week. The proportion of female patients in the hypokalemia group was higher than that in the nonhypokalemia group, as was the proportion of patients receiving intravenous VCZ. In the multivariate analysis, the independent risk factors for hypokalemia were sex, combined use of antibiotics, and VCZ CN/C0. The incidence of hyponatremia was 7.9% (21/266). The proportion of patients over 47 years of age in the hyponatremia group was 71.4% (15/21). The number of days of VCZ use in the hyponatremia group was greater than that in the nonhyponatremia group. A total of 47.6% (10/21) of patients in the hyponatremia group had supratherapeutic VCZ C0 (>5.0 µg/mL). In conclusion, hypokalemia is more likely to occur in females, in patients receiving intravenous VCZ, and in patients with the combined use of antibiotics. Hyponatremia is more likely to occur in patients older than 47 years who have been using VCZ for a long time and have higher VCZ C0 values.

Metabolic characteristics of voriconazole - Induced liver injury in rats

Voriconazole (VOR) - induced liver injury is a common adverse reaction, and can lead to serious clinical outcomes. It is of great significance to describe the metabolic characteristics of VOR - induced liver injury and to elucidate the potential mechanisms. This study investigated the changes of plasma metabolic profiles in a rat model of VOR - induced liver injury by non - targeted metabolomics. Correlation analysis was performed between differentially expressed metabolites and plasma liver function indexes. The metabolites with strong correlation were determined for their predictive performance for liver injury using receiver operating characteristic (ROC) curve analysis. Potential biomarkers were then screened combined with liver pathological scores. Finally, the expression level of genes that involved in lipid metabolism were determined in rat liver to verify the mechanism of VOR - induced liver injury we proposed. VOR - induced liver injury in rats was characterized by plasma alanine aminotransferase (ALT) and aspartate aminotransferase (AST) elevation, the lipid droplets accumulation in liver, as well as inflammation and fibrosis. Significant changes of plasma metabolites were observed, with a decrease in lipid metabolites accounting for over 50% of all changed metabolites, and alterations of cholesterol and bile acids metabolites. The decrease of 3 phosphatidylcholine (PC) in plasma could indicate the occurrence of VOR - induced liver injury. Decreased fatty acids (FA) oxidation and bile acid excretion might be the potential mechanisms of VOR - induced liver injury. This study provided new insights into the molecular characterization of VOR - induced liver injury.

Prediction of plasma trough concentration of voriconazole in adult patients using machine learning

OBJECTIVE

Plasma trough concentration of voriconazole (VCZ) was associated with its toxicity and efficacy. However, the nonlinear pharmacokinetic characteristics of VCZ make it difficult to determine the relationship between clinical characteristics and its concentration. We intended to present a machine learning (ML)-based method to predict toxic plasma trough concentration of VCZ (>5 μg/mL).

METHODS

A single center retrospective study was conducted. Three ML algorithms were used to estimate the concentration in adult patients, including random forest (RF), gradient boosting (GB), and extreme gradient boosting (XGBoost). The importance of variables was recognized by the SHapley Additive exPlanations (SHAP) method. In addition, an external validation set was used to validate the robustness of models.

RESULTS

A total of 1318 VCZ plasma concentration were included, with 33 variables enrolled in the model. Nine classification models were developed using the RF, GB, and XGBoost algorithms. Most models performed well for both the training set and test set, with an average balanced accuracy (BA) of 0.704 and an average accuracy (ACC) of 0.788. In addition, the average Matthews correlation coefficient value reached 0.484, which indicated the predicted values are meaningful. Based on the average BA and ACC values, the predictive ability of the models can be ranked from best to worst as follows: younger adult models > mixed models > elderly models, and XGBoost models > GBT models > RF models. The SHAP results showed that the top five influencing factors in younger adult patients (<60 years) were albumin, total bile acid (TBA), platelets count, age, and inflammation, while the top five influencing factors in elderly patients were albumin, TBA, aspartate aminotransferase, creatinine, and alanine aminotransferase. Furthermore, the prediction of external validation set for VCZ concentrations verified the high reliability of the models, for the ACC value of 0.822 by the best model.

CONCLUSIONS

The ML models can be reliable tools for predicting toxic concentration exposure of VCZ. The SHAP results may provide useful guidelines for dosage adjustment of VCZ.

The Development and Validation of a Predictive Model for Voriconazole-Related Liver Injury in Hospitalized Patients in China

Voriconazole is widely used in the treatment and prevention of invasive fungal diseases. Common drug-induced liver injuries increase the economic burdens and the risks of premature drug withdrawal and disease recurrence. This study estimated the disposal cost of voriconazole-related liver injury, explored the risk factors of voriconazole-related liver injury in hospitalized patients, and established a predictive model of liver injury to assist clinicians and pharmacists in estimating the probability or risk of liver injury after voriconazole administration to allow for early identification and intervention in patients at high risk of liver injury. A retrospective study was conducted on the selected inpatients whose blood concentration of voriconazole was measured in the West China Hospital of Sichuan University from September 2016 to June 2020. The incidence and disposal cost of voriconazole-related liver injuries were calculated. The incidence of voriconazole-related liver injury was 15.82% (217/1372). The disposal cost has been converted to 2023 at a discount rate of 5%. The median (P₂₅, P₇₅) disposal cost of severe liver injury (n = 42), general liver injury (n = 175), and non-liver injury (n = 1155) was 993.59 (361.70, 1451.76) Chinese yuan, 0.00 (0.00, 410.48) yuan, and 0.00 (0.00, 0.00) yuan, respectively, with a statistically significant difference (p < 0.001). Single factor analysis and multiple factor logistic regression were used to analyze the risk factors of voriconazole-related liver injury. The voriconazole-related liver injury was related to the trough concentration (C_min, OR 1.099, 95% CI 1.058-1.140), hypoproteinemia (OR 1.723, 95% CI 1.126-2.636), and transplantation status (OR 0.555, 95% CI 0.325-0.948). The prediction model of liver injury was Logit (P)= -2.219 + 0.094 × C_min + 0.544 × H_{ydroproteinemia} - 0.589 × Transplantation, and the prediction model nomogram was established. The model validation results showed that the C-index of the derivation set and validation set was 0.706 and 0.733, respectively. The area under the curve (AUC) of the receiver operating characteristic (ROC) curve was 0.705 and 0.733, respectively, indicating that the model had good prediction ability. The prediction model will be helpful to develop clinical individualized medication of voriconazole and to identify and intervene in the cases of patients at high risk of voriconazole-related liver injury early on, in order to reduce the incidence of voriconazole-related liver injuries and the cost of treatment.

Population Pharmacokinetics of Isavuconazole in Critical Care Patients with COVID-19-Associated Pulmonary Aspergillosis and Monte Carlo Simulations of High Off-Label Doses

Isavuconazole is a triazole antifungal agent recently recommended as first-line therapy for invasive pulmonary aspergillosis. With the COVID-19 pandemic, cases of COVID-19-associated pulmonary aspergillosis (CAPA) have been described with a prevalence ranging from 5 to 30%. We developed and validated a population pharmacokinetic (PKpop) model of isavuconazole plasma concentrations in intensive care unit patients with CAPA. Nonlinear mixed-effect modeling Monolix software were used for PK analysis of 65 plasma trough concentrations from 18 patients. PK parameters were best estimated with a one-compartment model. The mean of ISA plasma concentrations was 1.87 [1.29-2.25] mg/L despite prolonged loading dose (72 h for one-third) and a mean maintenance dose of 300 mg per day. Pharmacokinetics (PK) modeling showed that renal replacement therapy (RRT) was significantly associated with under exposure, explaining a part of clearance variability. The Monte Carlo simulations suggested that the recommended dosing regimen did not achieve the trough target of 2 mg/L in a timely manner (72 h). This is the first isavuconazole PKpop model developed for CAPA critical care patients underlying the need of therapeutic drug monitoring, especially for patients under RRT.

Therapeutic drug monitoring and safety evaluation of voriconazole in the treatment of pulmonary fungal diseases

Aims:

The gene polymorphism of voriconazole metabolism–related liver enzyme is notable in East Asia population. It casts a significant influence on the rational use of voriconazole. We conducted this study to investigate the relationship between steady-state voriconazole trough concentration (C_trough) and adverse effects (AEs), especially hepatotoxicity.

Methods:

We conducted a real-world study in the Jinling Hospital from January 2015 to June 2020. A total of 140 patients receiving voriconazole were enrolled in this study. The determination and scoring of voriconazole-associated hepatotoxicity were performed according to the Roussel Uclaf Causality Assessment Method scoring scale and the severity of hepatotoxicity was graded according to the Common Terminology Criteria for Adverse Events (CTCAE).

Results:

Elevated steady-state voriconazole C_trough with concomitant AEs are the most common reason for dose adjustments during treatment. Compared with the group without any AEs, voriconazole C_trough was significantly higher in the hepatotoxicity and neurotoxicity groups, and the incidence of both events showed an overall increasing trend with increasing voriconazole C_trough. Hepatotoxicity occurred in 66.7% of patients within 7 days of the first dose of voriconazole and 94.4% within 15 days of the dose. Steady-state voriconazole C_trough >3.61 mg/l was associated with an increased incidence of hepatotoxicity (area under the curve = 0.645, p = 0.047). Logistic regression analysis showed that timely voriconazole dose adjustment was a predictor of attenuated hepatotoxicity after adjustment for confounders, but hepatotoxicity was not associated with voriconazole C_trough measured at a single time point.

Conclusion:

Hepatotoxicity and neurotoxicity correlate with voriconazole C_trough, and dose reduction in patients with elevated steady-state voriconazole C_trough may prevent hepatotoxicity. In patients with early occurrence of hepatotoxicity, initial therapeutic drug monitoring (TDM) might predict the risk of hepatotoxicity. Follow-up TDM may be necessary to predict late onset hepatotoxicity.

Plain Language Summary

Safety of voriconazole for the treatment of pulmonary fungal diseases

Introduction: Several studies have suggested an association between the concentration of voriconazole in the blood and liver damage, but the evidence is weak. This study aimed to investigate relationships between voriconazole drug concentration and side effects and to analyze the factors affecting liver damage caused by voriconazole.

Methods: We conducted a study at the Jinling Hospital from January 2015 to June 2020, in which a total of 140 patients were finally enrolled.

Results: Voriconazole doses were adjusted in 44 patients due to abnormal voriconazole drug concentration or side effects, 32 patients reduced the dose and 8 patients increased the dose. An elevated liver enzyme level was the most common cause for dose adjustment. After the first dose adjustment, most patients achieved the target drug concentration. A total of 18 patients were determined as probable or highly probable to have drug-induced liver injury from voriconazole. Voriconazole drug concentration was significantly higher in the liver damage and nervous system damage groups as compared with the group without any side effects, and most liver damage events occurred within 14 days of the first dose. Voriconazole drug concentration >3.61 mg/l was associated with an increased incidence of liver damage.

Conclusion: In this study, approximately one-third of patients with pulmonary fungal disease needed to adjust their dose after the standard dose of voriconazole treatment. The incidence of liver damage and nervous system damage showed an overall increasing trend with increasing voriconazole baseline concentrations. Initial therapeutic drug monitoring may be predictive of liver damage. Follow-up monitoring of liver enzymes may be needed.

Impact of cytochrome P450 2C19 polymorphisms on the clinical efficacy and safety of voriconazole: an update systematic review and meta-analysis

OBJECTIVE

To assess the impact of cytochrome P450 (CYP) 2C19 polymorphisms on the clinical efficacy and safety of voriconazole.

METHODS

We systematically searched PubMed, EMBASE, CENTRAL, ClinicalTrials.gov, and three Chinese databases from their inception to 18 March 2021 using a predefined search algorithm to identify relevant studies. Studies that reported voriconazole-treated patients and information on CYP2C19 polymorphisms were included. The efficacy outcome was success rate. The safety outcomes included overall adverse events, hepatotoxicity, and neurotoxicity.

RESULTS

A total of 20 studies were included. Intermediate metabolizers (IMs) and poor metabolizers (PMs) were associated with increased success rates compared with normal metabolizers (NMs) [risk ratio (RR), 1.18; 95% confidence interval (CI), 1.03-1.34; I2 = 0%; P = 0.02; RR, 1.28; 95% CI, 1.06-1.54; I2 = 0%; P = 0.01]. PMs were at increased risk of overall adverse events in comparison with NMs and IMs (RR, 2.18; 95% CI, 1.35-3.53; I2 = 0%; P = 0.001; RR, 1.80; 95% CI, 1.23-2.64; I2 = 0%; P = 0.003). PMs demonstrated a trend towards an increased incidence of hepatotoxicity when compared with NMs (RR, 1.60; 95% CI, 0.94-2.74; I2 = 27%; P = 0.08), although there was no statistically significant difference. In addition, there was no significant association between CYP2C19 polymorphisms and neurotoxicity.

CONCLUSION

IMs and PMs were at a significant higher success rate in comparison with NMs. PMs were significantly associated with an increased incidence of all adverse events compared with NMs and IMs. Researches are expected to further confirm these findings. Additionally, the relationship between hepatotoxicity and CYP2C19 polymorphisms deserves clinical attention.

Association between voriconazole exposure and Sequential Organ Failure Assessment (SOFA) score in critically ill patients

Therapeutic drug monitoring (TDM) is essential for voriconazole to ensure optimal drug exposure, mainly in critically ill patients for whom voriconazole demonstrated a large variability. The study aimed at describing factors associated with trough voriconazole concentrations in critically ill patients and evaluating the impact of voriconazole concentrations on adverse effects. A 2-year retrospective multicenter cohort study (NCT04502771) was conducted in six intensive care units. Adult patients who had at least one voriconazole TDM were included. Univariable and multivariable linear regression analyses were performed to identify predictors of voriconazole concentrations, and univariable logistic regression analysis, to study the relationship between voriconazole concentrations and adverse effects. During the 2-year study period, 70 patients were included. Optimal trough voriconazole concentrations were reported in 37 patients (52.8%), subtherapeutic in 20 (28.6%), and supratherapeutic in 13 (18.6%). Adverse effects were reported in six (8.6%) patients. SOFA score was identified as a factor associated with an increase in voriconazole concentration (p = 0.025), mainly in the group of patients who had SOFA score ≥ 10. Moreover, an increase in voriconazole concentration was shown to be a risk factor for occurrence of adverse effects (p = 0.011). In that respect, critically ill patients who received voriconazole treatment must benefit from a TDM, particularly if they have a SOFA score ≥ 10. Indeed, identifying patients who are overdosed will help to prevent voriconazole related adverse effects. This result is of utmost importance given the recognized COVID-19-associated pulmonary aspergillosis in ICU patients for whom voriconazole is among the recommended first-line treatment.

Model-Oriented Dose Optimization of Voriconazole in Critically Ill Children

This study aimed to employ a population pharmacokinetic (PK) model to optimize the dosing regimen of voriconazole (VRC) in children with a critical illness. A total of 99 children aged from 0.44 to 13.58 years were included in this study. The stability and predictive performance of the final model were evaluated by statistical and graphical methods. The optimal dosing regimen was proposed for children with different body weights, CYP2C19 phenotypes, and coadministrations with omeprazole. The PK of VRC was described by a two-compartment model with nonlinear Michaelis-Menten elimination. Body weight, CYP2C19 phenotype, and omeprazole were significant covariates on the maximum velocity of elimination (V_max), which had an estimated typical value of 18.13 mg · h^-1. Bayesian estimation suggested that the dose-normalized concentration and total exposure (peak concentration [C_max]/D, trough concentration [C_min]/D, and area under the concentration-time curve over 24 h [AUC₂₄]/D) were significantly different between extensive metabolizer (EM) patients and poor metabolizer (PM) patients. To achieve the target concentration early, two loading doses of 9 mg · kg^-1 of body weight every 12 h (q12h) were reliable for most children, whereas three loading doses of 6 to 7.5 mg · kg^-1 q8h were warranted for young children weighing ≤18 kg (except for PM patients). The maintenance doses decreased about 30 to 40% in PM patients compared to that in EM patients. For children aged <2 years, in EM patients, the maintenance dose could be as high as 9 mg · kg^-1. The maintenance dose of VRC was supposed to decrease slightly when coadministered with omeprazole. A population PK model of intravenous VRC for critically ill children has been successfully developed. It is necessary to adjust dosing regimens according to the CYP2C19 genotype. Optimal dosing regimens have been recommended based on the final model.

Application of a Physiologically Based Pharmacokinetic Model to Characterize Time-dependent Metabolism of Voriconazole in Children and Support Dose Optimization

Background: Voriconazole is a potent antifungal drug with complex pharmacokinetics caused by time-dependent inhibition and polymorphisms of metabolizing enzymes. It also exhibits different pharmacokinetic characteristics between adults and children. An understanding of these alterations in pharmacokinetics is essential for pediatric dose optimization.

Objective: To determine voriconazole plasma exposure in the pediatric population and further investigate optimal dosage regimens.

Methods: An adult and pediatric physiologically based pharmacokinetic (PBPK) model of voriconazole, integrating auto-inhibition of cytochrome P450 3A4 (CYP3A4) and CYP2C19 gene polymorphisms, was developed. The model was evaluated with visual predictive checks and quantitative measures of the predicted/observed ratio of the area under the plasma concentration-time curve (AUC) and maximum concentration (C_max). The validated pediatric PBPK model was used in simulations to optimize pediatric dosage regimens. The probability of reaching a ratio of free drug (unbound drug concentration) AUC during a 24-h period to minimum inhibitory concentration greater than or equal to 25 (fAUC_24h/MIC ≥ 25) was assessed as the pharmacokinetic/pharmacodynamic index.

Results: The developed PBPK model well represented voriconazole's pharmacokinetic characteristics in adults; 78% of predicted/observed AUC ratios and 85% of C_max ratios were within the 1.25-fold range. The model maintained satisfactory prediction performance for intravenous administration in pediatric populations after incorporating developmental changes in anatomy/physiology and metabolic enzymes, with all predicted AUC values within 2-fold and 73% of the predicted C_max within 1.25-fold of the observed values. The simulation results of the PBPK model suggested that different dosage regimens should be administered to children according to their age, CYP2C19 genotype, and infectious fungal genera.

Conclusion: The PBPK model integrating CYP3A4 auto-inhibition and CYP2C19 gene polymorphisms successfully predicted voriconazole pharmacokinetics during intravenous administration in children and could further be used to optimize dose strategies. The infectious fungal genera should be considered in clinical settings, and further research with large sample sizes is required to confirm the current findings.

Model-based Voriconazole Dose Optimization in Chinese Adult Patients With Hematologic Malignancies

PURPOSE

The objective of this study was to characterize the population pharmacokinetics of voriconazole and to identify factors that significantly affect pharmacokinetic parameters and to further investigate optimal dosage regimens in Chinese adult patients with hematologic malignancies.

METHODS

A prospective population pharmacokinetic analysis was performed on 186 concentration measurements obtained from 41 adult patients with hematologic malignancies. All enrolled patients were treated with voriconazole for diagnosed or suspected invasive fungal diseases. Oral voriconazole was routinely administered at a maintenance dose of 200 mg q12h. Serial blood samples were collected after steady-state of each patient. Monte Carlo simulation was applied to optimize dosage strategies.

FINDINGS

A one-compartment model with first-order absorption and elimination adequately described the data. The typical voriconazole clearance was 4.18 L/h, the volume of distribution was 88.9 L, and the absorption rate constant was 0.729 h^-1. Clearance and steady-state exposure (AUC_0-12) were found to be significantly associated with age and CYP2C19 phenotype. The average AUC_0-12 of elderly patients (aged 60-90 years) was 2.1 times higher than that of relative younger patients (aged 18-59 years). The average AUC_0-12 of poor metabolizers (PMs) was approximately 2.5 and 1.8 times higher than that of extensive and intermediate metabolizers (IMs), respectively. Considering both efficacy and tolerability, dosage regimens of 100 and 50 mg orally administered every 12 hours were recommended for elderly IMs and PMs, respectively.

IMPLICATIONS

A population pharmacokinetic model for voriconazole in Chinese adult patients with hematologic malignancies was successfully developed and could well capture voriconazole's pharmacokinetic characteristics. Age and CYP2C19 phenotype were found to significantly influence voriconazole clearance and should be taken into consideration clinically for dose optimization. The optimal dosage strategies in specific clinical scenarios were proposed in this study based on model simulation. Because of the high incidence of mutant CYP2C19*2 and *3 alleles, genetic testing seems to be necessary for Asian elderly patients when voriconazole treatment is initiated.

Impact of voriconazole plasma concentrations on treatment response in critically ill patients

WHAT IS KNOWN AND OBJECTIVE

Several authors have demonstrated the relationship between voriconazole concentrations and the risk of therapeutic failure and adverse events However, the information about voriconazole concentrations in the critically ill patient is scarce. The aim of this study was to analyse the plasma concentrations and pharmacokinetic behaviour of voriconazole in critically ill patients and their association with the treatment response and development of toxicity.

METHODS

A prospective, observational study was conducted. Patients admitted to an intensive care unit and on treatment with intravenous voriconazole were included. Plasma concentrations were measured between days 4 and 7 from the start of the treatment. The pharmacokinetic analysis was performed using the NONMEM^® software. A regression model was used to evaluate the variables associated with the values outside the therapeutic range, as well as the relationship between the plasma concentrations and the treatment response and the development of hepatotoxicity.

RESULTS AND DISCUSSION

A total of 33 patients were included. Plasma concentrations outside the therapeutic range (1-5.5 mg/L) were observed in 15 patients, being above the established range in 9 (27.3%) cases, and below it in 6 (18.2%) cases. The presence of a bilirubin value of >1.5 mg/dL and a C-reactive protein >100 mg/dL was associated with supra-therapeutic concentrations. Voriconazole concentrations greater than 5.5 mg/dL were associated with the development of hepatotoxicity.

WHAT IS NEW AND CONCLUSIONS

There is a wide variation in voriconazole concentrations in critically ill patients, being associated with a high frequency of adverse events. Close monitoring of these values is required in order to decrease the risk of therapeutic failure and toxicity.