Pharmacokinetic–pharmacodynamic analysis of voriconazole in Japanese patients with hematological malignancies

Evaluation of Pharmacokinetics and Safety With Bioequivalence of Voriconazole Injection of 2 Formulations in Chinese Healthy Volunteers: Bioequivalence Study

Voriconazole is a first-line medicine for treating invasive aspergillosis. We aimed to evaluate the bioequivalence (BE) of voriconazole injection in Chinese healthy volunteers (HVs). In this single-center, randomized, single-dose, 2-cycle, fasting-dose BE study, HVs (n = 24) were 1:1 divided into 2 groups (test [T]-reference [R] and R-T) and received 6 mg/kg of voriconazole intravenously with a 7-day washout. The plasma was collected for up to 72 hours at the time point after dosing on day 1/day 8. The plasma concentration of voriconazole was measured by liquid chromatography-tandem mass spectrometry. Pharmacokinetic parameters were ascertained on the basis of a noncompartmental model. In the BE study, the geometric mean ratios of the maximum concentration, area under the concentration-time curve from time 0 to the last measurable plasma concentration, and area under the concentration-time curve from time 0 to infinity were 101.1%, 105.6%, and 105.5%, respectively, and the 90%CI fell within 80%-125%. Adverse events were observed in 26.1% of subjects in the T formulation stage and 17.4% in the R formulation stage. Under the BE study, voriconazole values from T and R formulations were bioequivalent.

Voriconazole exposure is influenced by inflammation: A population pharmacokinetic model

BACKGROUND

Voriconazole is an antifungal drug used for the treatment of invasive fungal infections. Due to highly variable drug exposure, therapeutic drug monitoring (TDM) has been recommended. TDM may be helpful to predict exposure accurately, but covariates, such as severe inflammation, that influence the metabolism of voriconazole have not been included in the population pharmacokinetic (popPK) models suitable for routine TDM.

OBJECTIVES

To investigate whether the effect of inflammation, reflected by C-reactive protein (CRP), could improve a popPK model that can be applied in clinical care.

PATIENTS AND METHODS

Data from two previous studies were included in the popPK modelling. PopPK modelling was performed using Edsim++. Different popPK models were compared using Akaike Information Criterion and goodness-of-fit plots.

RESULTS

In total, 1060 voriconazole serum concentrations from 54 patients were included in this study. The final model was a one-compartment model with non-linear elimination. Only CRP was a significant covariate, and was included in the final model and found to affect the maximum rate of enzyme activity (V_max). For the final popPK model, the mean volume of distribution was 145 L [coefficient of variation percentage (CV%)=61%], mean Michaelis-Menten constant was 5.7 mg/L (CV%=119%), mean V_max was 86.4 mg/h (CV%=99%) and mean bioavailability was 0.83 (CV%=143%). Internal validation using bootstrapping resulted in median values close to the population parameter estimates.

CONCLUSIONS

This one-compartment model with non-linear elimination and CRP as a covariate described the pharmacokinetics of voriconazole adequately.

Population pharmacokinetics of voriconazole and initial dosage optimization in patients with talaromycosis

The high variability and unpredictability of the plasma concentration of voriconazole (VRC) pose a major challenge for clinical administration. The aim of this study was to develop a population pharmacokinetics (PPK) model of VRC and identify the factors influencing VRC PPK in patients with talaromycosis. Medical records and VRC medication history of patients with talaromycosis who were treated with VRC as initial therapy were collected. A total of 233 blood samples from 69 patients were included in the study. A PPK model was developed using the nonlinear mixed-effects models (NONMEM). Monte Carlo simulation was applied to optimize the initial dosage regimens with a therapeutic range of 1.0–5.5 mg/L as the target plasma trough concentration. A one-compartment model with first-order absorption and elimination adequately described the data. The typical voriconazole clearance was 4.34 L/h, the volume of distribution was 97.4 L, the absorption rate constant was set at 1.1 h^-1, and the bioavailability was 95.1%. Clearance was found to be significantly associated with C-reactive protein (CRP). CYP2C19 polymorphisms had no effect on voriconazole pharmacokinetic parameters. ‏Monte Carlo simulation based on CRP levels showed that a loading dose of 250 mg/12 h and a maintenance dose of 100 mg/12 h are recommended for patients with CRP ≤ 96 mg/L, whereas a loading dose of 200 mg/12 h and a maintenance dose of 75 mg/12 h are recommended for patients with CRP > 96 mg/L. The average probability of target attainment of the optimal dosage regimen in CRP ≤ 96 mg/L and CRP > 96 mg/L groups were 61.3% and 13.6% higher than with empirical medication, and the proportion of C_min > 5.5 mg/L decreased by 28.9%. In conclusion, the VRC PPK model for talaromycosis patients shows good robustness and predictive performance, which can provide a reference for the clinical individualization of VRC. Adjusting initial dosage regimens based on CRP may promote the rational use of VRC.

Systematic Evaluation of Voriconazole Pharmacokinetic Models without Pharmacogenetic Information for Bayesian Forecasting in Critically Ill Patients

Voriconazole (VRC) is used as first line antifungal agent against invasive aspergillosis. Model-based approaches might optimize VRC therapy. This study aimed to investigate the predictive performance of pharmacokinetic models of VRC without pharmacogenetic information for their suitability for model-informed precision dosing. Seven PopPK models were selected from a systematic literature review. A total of 66 measured VRC plasma concentrations from 33 critically ill patients was employed for analysis. The second measurement per patient was used to calculate relative Bias (rBias), mean error (ME), relative root mean squared error (rRMSE) and mean absolute error (MAE) (i) only based on patient characteristics and dosing history (a priori) and (ii) integrating the first measured concentration to predict the second concentration (Bayesian forecasting). The a priori rBias/ME and rRMSE/MAE varied substantially between the models, ranging from −15.4 to 124.6%/−0.70 to 8.01 mg/L and from 89.3 to 139.1%/1.45 to 8.11 mg/L, respectively. The integration of the first TDM sample improved the predictive performance of all models, with the model by Chen (85.0%) showing the best predictive performance (rRMSE: 85.0%; rBias: 4.0%). Our study revealed a certain degree of imprecision for all investigated models, so their sole use is not recommendable. Models with a higher performance would be necessary for clinical use.

External evaluation of population pharmacokinetic models for voriconazole in Chinese adult patients with hematological malignancy

OBJECTIVES

Patients with hematological malignancies are prone to invasive fungal disease due to long-term chemotherapy or radiotherapy. Voriconazole is a second-generation triazole broad-spectrum antibiotic used to prevent or treat invasive fungal infections. Many population pharmacokinetic (pop PK) models have been published for voriconazole, and various diagnostic methods are available to validate the performance of these pop PK models. However, most of the published models have not been strictly evaluated externally. The purpose of this study is to evaluate these models externally and assess their predictive capabilities.

METHODS

The external dataset consists of adults receiving voriconazole treatment at Fujian Medical University Union Hospital. We re-established the published models based on their final estimated values in the literature and used our external dataset for initial screening. Each model was evaluated based on the following outcomes: prediction-based diagnostics, prediction- and variability-corrected visual predictive check (pvcVPC), normalized prediction distribution errors (NPDE), and Bayesian simulation results with one to two prior observations.

RESULTS

A total of 237 samples from 166 patients were collected as an external dataset. After screening, six candidate models suitable for the external dataset were finally obtained for comparison. Among the models, none demonstrated excellent predictive performance. Bayesian simulation shows that all models' prediction precision and accuracy were significantly improved when one or two prior concentrations were given.

CONCLUSIONS

The published pop PK models of voriconazole have significant differences in prediction performance, and none of the models could perfectly predict the concentrations of voriconazole for our data. Therefore, extensive evaluation should precede the adoption of any model in clinical practice.

Using Child-Pugh Class to Optimize Voriconazole Dosage Regimens and Improve Safety in Patients with Liver Cirrhosis: Insights from a Population Pharmacokinetic Model-based Analysis

BACKGROUND

Cirrhotic patients are at a high risk of fungal infections. Voriconazole is widely used as prophylaxis and in the treatment of invasive fungal disease. However, the safety, pharmacokinetics, and optimal regimens of voriconazole are currently not well defined in cirrhotic patients.

DESIGN

Retrospective pharmacokinetics study.

SETTING

Two large, academic, tertiary-care medical center.

PATIENTS

Two hundred nineteen plasma trough concentrations (C_min ) from 120 cirrhotic patients and 83 plasma concentrations from 11 non-cirrhotic patients were included.

METHODS

Data pertaining to voriconazole were collected retrospectively. A population pharmacokinetics analysis was performed and model-based simulation was used to optimize voriconazole dosage regimens.

RESULTS

Voriconazole-related adverse events (AEs) developed in 29 cirrhotic patients, and the threshold C_min for AE was 5.12 mg/L. A two-compartment model with first-order elimination adequately described the data. The Child-Pugh class and body weight were the significant covariates in the final model. Voriconazole clearance in non-cirrhotic, Child-Pugh class A and B cirrhotic (CP-A/B) and Child-Pugh class C cirrhotic (CP-C) patients was 7.59, 1.86, and 0.93 L/hour, respectively. The central distribution volume and peripheral distribution volume was 100.8 and 55.2 L, respectively. The oral bioavailability was 91.6%. Model-based simulations showed that a loading dose regimen of 200 mg/12 hours intravenously or orally led to 65.0-75.7% of voriconazole C_min in therapeutic range on day 1, and the appropriate maintenance dosage regimens were 75 mg/12 hours and 150 mg/24 hours intravenously or orally for CP-A/B patients, and 50 mg/12 hours and 100 mg/24 hours intravenously or orally for CP-C patients. The predicted probability of achieving the therapeutic target concentration for optimized regimens at steady-state was 66.8-72.3% for CP-A/B patients and 70.3-74.0% for CP-C patients.

CONCLUSIONS

These results recommended that the halved loading dose regimens should be used, and voriconazole maintenance doses in cirrhotic patients should be reduced to one-fourth for CP-C patients and to one-third for CP-A/B patients compared to that for patients with normal liver function.

Population Pharmacokinetics of Voriconazole in Patients With Invasive Aspergillosis: Serum Albumin Level as a Novel Marker for Clearance and Dosage Optimization

BACKGROUND

Voriconazole (VRCZ) is an antifungal triazole recommended as an effective first-line agent for treating invasive aspergillosis.

OBJECTIVES

To develop a population pharmacokinetic model of VRCZ and trough concentration-based dosing simulation for dynamic patient conditions.

METHODS

The authors combined plasma VRCZ data from intensive sampling, and retrospective trough concentration monitoring for analysis. Nonlinear mixed-effects modeling with subsequent model validation was performed. The recommended dosage regimens were simulated based on the developed model.

RESULTS

The study participants included 106 patients taking oral VRCZ. A linear one-compartment model with first-order elimination and absorption best described the observed data. The CYP2C19 phenotypes did not influence the pharmacokinetic parameters. Serum albumin (SA) levels and gamma-glutamyl transferase significantly correlated with the VRCZ clearance rate, whereas the actual body weight influenced the volume. A visual predictive check showed good consistency with the observed data, whereas SA levels across the treatment course correlated with linear clearance, irrespective of the CYP2C19 phenotype. Patients with SA levels ≤30 g/L had lower linear clearance than that in patients with SA levels >30 g/L. Dosing simulation based on the developed model indicated that patients with SA levels of ≤30 g/L required a lower daily maintenance dose to attain the therapeutic trough level.

CONCLUSIONS

SA level was identified as a novel marker associated with VRCZ clearance. This marker may be a practical choice for physicians to perform therapeutic drug monitoring and optimize VRCZ dosage.

Impact of Albumin and Omeprazole on Steady-State Population Pharmacokinetics of Voriconazole and Development of a Voriconazole Dosing Optimization Model in Thai Patients with Hematologic Diseases

This study aimed to identify factors that significantly influence the pharmacokinetics of voriconazole in Thai adults with hematologic diseases, and to determine optimal voriconazole dosing regimens. Blood samples were collected at steady state in 65 patients (237 concentrations) who were taking voriconazole to prevent or treat invasive aspergillosis. The data were analyzed using a nonlinear mixed-effects modeling approach. Monte Carlo simulation was applied to optimize dosage regimens. Data were fitted with the one-compartment model with first-order absorption and elimination. The apparent oral clearance (CL/F) was 3.43 L/h, the apparent volume of distribution (V/F) was 47.6 L, and the absorption rate constant (Ka) was fixed at 1.1 h⁻¹. Albumin and omeprazole ≥ 40 mg/day were found to significantly influence CL/F. The simulation produced the following recommended maintenance doses of voriconazole: 50, 100, and 200 mg every 12 h for albumin levels of 1.5–3, 3.01–4, and 4.01–4.5 g/dL, respectively, in patients who receive omeprazole ≤ 20 mg/day. Patients who receive omeprazole ≥ 40 mg/day and who have serum albumin level 1.5–3 and 3.01–4.5 g/dL should receive voriconazole 50 and 100 mg, every 12 h, respectively. Albumin level and omeprazole dosage should be carefully considered when determining the appropriate dosage of voriconazole in Thai patients.

Saliva for Precision Dosing of Antifungal Drugs: Saliva Population PK Model for Voriconazole Based on a Systematic Review

Precision dosing for many antifungal drugs is now recommended. Saliva sampling is considered as a non-invasive alternative to plasma sampling for therapeutic drug monitoring (TDM). However, there are currently no clinically validated saliva models available. The aim of this study is firstly, to conduct a systematic review to evaluate the evidence supporting saliva-based TDM for azoles, echinocandins, amphotericin B, and flucytosine. The second aim is to develop a saliva population pharmacokinetic (PK) model for eligible drugs, based on the evidence. Databases were searched up to July 2019 on PubMed® and Embase®, and 14 studies were included in the systematic review for fluconazole, voriconazole, itraconazole, and ketoconazole. No studies were identified for isavuconazole, posaconazole, flucytosine, amphotericin B, caspofungin, micafungin, or anidulafungin. Fluconazole and voriconazole demonstrated a good saliva penetration with an average S/P ratio of 1.21 (± 0.31) for fluconazole and 0.56 (± 0.18) for voriconazole, both with strong correlation (r = 0.89-0.98). Based on the evidence for TDM and available data, population PK analysis was performed on voriconazole using Nonlinear Mixed Effects Modeling (NONMEM 7.4). 137 voriconazole plasma and saliva concentrations from 11 patients (10 adults, 1 child) were obtained from the authors of the included study. Voriconazole pharmacokinetics was best described by one-compartment PK model with first-order absorption, parameterized by clearance of 4.56 L/h (36.9% CV), volume of distribution of 60.7 L, absorption rate constant of 0.858 (fixed), and bioavailability of 0.849. Kinetics of the voriconazole distribution from plasma to saliva was identical to the plasma kinetics, but the extent of distribution was lower, modeled by a scale factor of 0.5 (4% CV). A proportional error model best accounted for the residual variability. The visual and simulation-based model diagnostics confirmed a good predictive performance of the saliva model. The developed saliva model provides a promising framework to facilitate saliva-based precision dosing of voriconazole.

Population Pharmacokinetics of Voriconazole in Chinese Patients with Hematopoietic Stem Cell Transplantation

BACKGROUND AND OBJECTIVE

Voriconazole is widely recommended for the prevention and treatment of invasive fungal infections in hematopoietic stem cell transplantation patients. However, its use is limited by a narrow therapeutic range and large inter-individual variability. This study aimed to characterize the pharmacokinetics of voriconazole in Chinese hematopoietic stem cell transplantation patients, to explore factors affecting its pharmacokinetic parameters, and to provide recommendations for its optimal dosing regimens.

METHODS

A total of 121 serum concentration samples from 23 patients were retrospectively included. Voriconazole concentrations were detected, and patient clinical data were recorded. Population pharmacokinetic analysis was performed by a non-linear, mixed-effect modeling approach. Goodness-of-fit plots, bootstrap method, prediction-corrected visual predictive check and external validation by an independent group of seven patients were performed to evaluate the final model.

RESULTS

A one-compartment model with first-order elimination successfully described the data. The absorption rate constant was fixed at 1.1 h^-1 and bioavailability was fixed at 0.895. The typical values for voriconazole clearance and distribution volume were 9.52 L/h and 155 L, respectively. CYP2C19*2 genotype and mycophenolate mofetil combination presented a significant impact on the clearance. Compared with CYP2C19*2 carriers, voriconazole clearance was proven to be higher in CYP2C19*1/*1 patients.

CONCLUSIONS

A population pharmacokinetic model of voriconazole was successfully established in Chinese hematopoietic stem cell transplantation patients. Based on the final model, CYP2C19*2 genotyping coupled with therapeutic drug monitoring seems to be useful to guide voriconazole dosing and to explain subtherapeutic concentrations in clinical practice.

Therapeutic drug monitoring of voriconazole for treatment and prophylaxis of invasive fungal infection in children

Voriconazole therapeutic drug monitoring is not consistently recommended due to its high interpatient and intrapatient variability. Here, we aimed to describe our experience with voriconazole for treatment and prophylaxis of invasive fungal infections in paediatric patients. A fully validated high-performance liquid chromatography-mass spectrometry method was used to quantify voriconazole concentration in plasma, at the end of dosing interval. A high interindividual variability was shown. We enrolled 237 children, 83 receiving intravenous and 154 oral voriconazole. A positive correlation between drug dose and drug plasma exposure was observed. Considering intravenous route, patients with higher serum creatinine had higher voriconazole concentrations; a positive correlation was found among drug exposure and age. Sex significantly influenced drug levels: males had higher median drug concentrations than females (P < 0.001). Close voriconazole pharmacokinetics monitoring should help individualize antifungal therapy for children.

Comparative Evaluation of the Predictive Performances of Three Different Structural Population Pharmacokinetic Models To Predict Future Voriconazole Concentrations

Bayesian methods for voriconazole therapeutic drug monitoring (TDM) have been reported previously, but there are only sparse reports comparing the accuracy and precision of predictions of published models. Furthermore, the comparative accuracy of linear, mixed linear and nonlinear, or entirely nonlinear models may be of high clinical relevance. In this study, models were coded into individually designed optimum dosing strategies (ID-ODS) with voriconazole concentration data analyzed using inverse Bayesian modeling. The data used were from two independent data sets, patients with proven or suspected invasive fungal infections (n = 57) and hematopoietic stem cell transplant recipients (n = 10). Observed voriconazole concentrations were predicted whereby for each concentration value, the data available to that point were used to predict that value. The mean prediction error (ME) and mean squared prediction error (MSE) and their 95% confidence intervals (95% CI) were calculated to measure absolute bias and precision, while ΔME and ΔMSE and their 95% CI were used to measure relative bias and precision, respectively. A total of 519 voriconazole concentrations were analyzed using three models. MEs (95% CI) were 0.09 (-0.02, 0.22), 0.23 (0.04, 0.42), and 0.35 (0.16 to 0.54) while the MSEs (95% CI) were 2.1 (1.03, 3.17), 4.98 (0.90, 9.06), and 4.97 (-0.54 to 10.48) for the linear, mixed, and nonlinear models, respectively. In conclusion, while simulations with the linear model were found to be slightly more accurate and similarly precise, the small difference in accuracy is likely negligible from the clinical point of view, making all three approaches appropriate for use in a voriconazole TDM program.

Pharmacokinetic/pharmacodynamic analysis of voriconazole against Candida spp. and Aspergillus spp. in children, adolescents and adults by Monte Carlo simulation

The objective of this study was to investigate the cumulative fraction of response of various voriconazole dosing regimens against six Candida and six Aspergillus spp. in immunocompromised children, immunocompromised adolescents, and adults. Using pharmacokinetic parameters and pharmacodynamic data, 5000-subject Monte Carlo simulations (MCSs) were conducted to evaluate the ability of simulated dosing strategies in terms of fAUC/MIC targets of voriconazole. According to the results of the MCSs, current voriconazole dosage regimens were all effective for children, adolescents and adults against Candida albicans, Candida parapsilosis and Candida orthopsilosis. For adults, dosing regimens of 4 mg/kg intravenous every 12 h (q12h) and 300 mg orally q12h were sufficient to treat fungal infections by six Candida spp. (C. albicans, C. parapsilosis, Candida tropicalis, Candida glabrata, Candida krusei and C. orthopsilosis) and five Aspergillus spp. (Aspergillus fumigatus, Aspergillus flavus, Aspergillus terreus, Aspergillus niger and Aspergillus nidulans). However, high doses should be recommended for children and adolescents in order to achieve better clinical efficacy against A. fumigatus and A. nidulans. The current voriconazole dosage regimens were all ineffective against A. niger for children and adolescents. All voriconazole dosage regimens were not optimal against Aspergillus versicolor. This is the first study to evaluate clinical therapy of various voriconazole dosing regimens against Candida and Aspergillus spp. infections in children, adolescents and adults using MCS. The pharmacokinetic/pharmacodynamic-based dosing strategy provided a theoretical rationale for identifying optimal voriconazole dosage regimens in children, adolescents and adults in order to maximise clinical response and minimise the probability of exposure-related toxicity.

Population Pharmacokinetics in China: The Dynamics of Intravenous Voriconazole in Critically Ill Patients with Pulmonary Disease

Pharmacokinetic research in China on the use of voriconazole in critically ill adult patients with different pulmonary diseases remains to be explored. This study evaluated the population pharmacokinetics of the use of voriconazole (VRC) in critically ill patients to determine covariate effects on VRC pharmacokinetics by NONMEM, which could further optimize VRC dosing in this population. A one-compartment model with first-order absorption and elimination best fit the data, giving 4.28 L/h clearance and 93.4 L volume of distribution of VRC. The model variability, described as an approximate percentage coefficient of interindividual variability in clearance and volume of distribution, was 72.94% and 26.50%, respectively. A significant association between Cmin and drug response or grade 2 hepatotoxicity was observed (p=0.002, <0.001, respectively, 1.5-4.0 µg/mL) via logistic multivariate regression. Monte Carlo simulations at 100, 150, 200, and 250 mg dosage predicted effectiveness at 45.99%, 99.76%, 98.76%, and 67.75% within the 1.5-4.0 µg/mL range, suggesting that a 150 or 200 mg intravenous dose twice daily is best suited to achieve the target steady state trough concentration range in critically ill patients with pulmonary disease.

Bioavailability and population pharmacokinetics of voriconazole in lung transplant recipients

This study was undertaken to characterize the pharmacokinetics and bioavailability of voriconazole in adult lung transplant patients during the early postoperative period, identify factors significantly associated with various pharmacokinetic parameters, and make recommendations for adequate dosing regimens. Thirteen lung transplant patients received two intravenous infusions (6 mg/kg, twice daily [b.i.d.]) immediately posttransplant followed by oral doses (200 mg, b.i.d.) for prophylaxis. Blood samples (9/interval) were collected during one intravenous and one oral dosing interval from each patient. Voriconazole plasma concentrations were measured by high-pressure liquid chromatography (HPLC). NONMEM was used to develop pharmacokinetic models, evaluate covariate relationships, and perform Monte Carlo simulations. There was a good correlation (R(2) = 0.98) between the area under the concentration-time curve specific for the dose evaluated (AUC(0-∞)) and trough concentrations. A two-compartment model adequately described the data. Population estimates of bioavailability, clearance, V(c), and V(p) were 45.9%, 3.45 liters/h, 54.7 liters, and 143 liters. Patients with cystic fibrosis (CF) exhibited a significantly lower bioavailability (23.7%, n = 3) than non-CF patients (63.3%, n = 10). Bioavailability increased with postoperative time and reached steady levels in about 1 week. V(p) increased with body weight. Bioavailability of voriconazole is substantially lower in lung transplant patients than non-transplant subjects but significantly increases with postoperative time. CF patients exhibit significantly lower bioavailability and exposure of voriconazole and therefore need higher doses. Intravenous administration of voriconazole during the first postoperative day followed by oral doses of 200 mg or 400 mg appeared to be the optimal dosing regimen. However, voriconazole levels should be monitored, and the dose should be individualized based on trough concentrations as a good measure of drug exposure.