ABSTRACT
Therapeutic drug monitoring (TDM) of moxifloxacin is recommended to improve the response to tuberculosis treatment and reduce acquired drug resistance. Limited sampling strategies (LSSs) are able to reduce the burden of TDM by using a small number of appropriately timed samples to estimate the parameter of interest, the area under the concentration-time curve. This study aimed to develop LSSs for moxifloxacin alone (MFX) and together with rifampin (MFX+RIF) in tuberculosis (TB) patients. Population pharmacokinetic (popPK) models were developed for MFX (n = 77) and MFX+RIF (n = 24). In addition, LSSs using Bayesian approach and multiple linear regression were developed. Jackknife analysis was used for internal validation of the popPK models and multiple linear regression LSSs. Clinically feasible LSSs (one to three samples, 6-h timespan postdose, and 1-h interval) were tested. Moxifloxacin exposure was slightly underestimated in the one-compartment models of MFX (mean –5.1%, standard error [SE] 0.8%) and MFX+RIF (mean –10%, SE 2.5%). The Bayesian LSSs for MFX and MFX+RIF (both 0 and 6 h) slightly underestimated drug exposure (MFX mean –4.8%, SE 1.3%; MFX+RIF mean –5.5%, SE 3.1%). The multiple linear regression LSS for MFX (0 and 4 h) and MFX+RIF (1 and 6 h), showed mean overestimations of 0.2% (SE 1.3%) and 0.9% (SE 2.1%), respectively. LSSs were successfully developed using the Bayesian approach (MFX and MFX+RIF; 0 and 6 h) and multiple linear regression (MFX, 0 and 4 h; MFX+RIF, 1 and 6 h). These LSSs can be implemented in clinical practice to facilitate TDM of moxifloxacin in TB patients.
INTRODUCTION
Each year, the global tuberculosis (TB) incidence declines with approximately 2%, while by 2020 an annual 4 to 5% decline is strived for by the World Health Organization (WHO) (1). Multidrug-resistant TB (MDR-TB) remains a major problem, with an estimated number of 458,000 cases in 2017 (1). Currently, the worldwide success rate of MDR-TB treatment is 55%, and this is considered low compared to a success rate of 85% for drug-susceptible TB (DS-TB) (1).
Moxifloxacin (MFX), a fluoroquinolone, is one of the most important drugs for the treatment of MDR-TB (2), but it has also been used as an alternative to first-line anti-TB drugs if not well tolerated or suggested to include in case of isoniazid resistance (3–5). In general, the toxicity profile of moxifloxacin is rather mild, though it includes concentration-dependent corrected QT interval prolongation and, rarely, tendinopathy (6–9). A clinically relevant drug-drug interaction is the combination of moxifloxacin with rifampin, since these two drugs can be used concomitantly in TB treatment. Rifampin (RIF) lowers the moxifloxacin area under the concentration-time curve of 0 to 24 h (AUC0–24) by approximately 30% by inducing phase II metabolizing enzymes (glucuronosyltransferase and sulfotransferase) (10–12).
The efficacy of fluoroquinolones is related to the ratio of AUC0–24 to MIC (AUC0–24/MIC) (13, 14). The fluoroquinolone exposure is effective against Gram-negative bacteria at an AUC0–24/MIC ratio of >100 to 125 and against Gram-positive species at an AUC0–24/MIC ratio of >25 to 30 (13, 15, 16). An in vitro moxifloxacin exposure of unbound (f)AUC0–24/MIC ratio of >53 was able to substantially decrease the total population of Mycobacterium tuberculosis with by >3 log10 CFU/ml, as well as suppress the emergence of drug resistance, while an fAUC0–24/MIC ratio of >102 completely killed the fluoroquinolone-sensitive population of M. tuberculosis without the development of drug resistance (17). Approximately 50% of moxifloxacin is assumed to be protein bound, although protein binding is highly variable between individuals and might be concentration dependent (13, 16, 18, 19). Corresponding with an fAUC0–24/MIC ratio of >53 and a fraction unbound of 0.5, the target total (bound and unbound) AUC0–24/MIC ratio of >100 to 125 is regularly used in TB, because individual data of protein binding is often lacking.(18, 20, 21) In case of a proven susceptibility for moxifloxacin while lacking a MIC value of the strain, the target AUC0–24 is generally set at >50 to 65 mg ⋅ h/liter based on a critical concentration of 0.5 mg/liter (22, 23).
Therapeutic drug monitoring (TDM) is recommended by the American Thoracic Society for all second-line drugs, including moxifloxacin (24, 25). It is important to monitor the moxifloxacin exposure in TB patients to determine an individualized dose, because of substantial interindividual pharmacokinetic variability and relevant drug-drug interactions with the risk of treatment failure and developing drug resistance (18, 26–28). However, routine TDM to estimate the AUC0–24 requiring frequent blood sampling is time-consuming, a burden for patients and health care professionals, and expensive. Optimizing the sampling schedule by developing a limited sampling strategy (LSS) could overcome these difficulties with TDM in TB treatment (29).
There are two main methods to develop an LSS: the Bayesian approach and multiple linear regression (30). The advantages of the Bayesian approach are the flexible timing of samples as the population pharmacokinetic model can correct for deviations and that it takes a number of parameters into account for example sex, age, and kidney function, leading to a more accurate estimation of AUC0–24. The advantage of multiple linear regression-based LSSs is that these do not require modeling software and AUC0–24 can be easily estimated using only an equation and the measurement of drug concentrations. The disadvantage is that samples must be taken exactly according to the predefined schedule and the population of interest should be comparable because patient characteristics are not included in the equations to estimate drug exposure (30).
Pranger et al. described a LSS for moxifloxacin for the first time using t = 4 and 14 h postdose samples (21). This sampling strategy can be considered unpractical to be used in daily practice. Magis-Escurra et al. described LSSs to simultaneously estimate the AUC0–24 of all first-line drugs, together with moxifloxacin (t = 1, 4, and 6 h or t = 2, 4, and 6 h) but did not differentiate between patients using moxifloxacin alone and moxifloxacin in combination with rifampin (20). Therefore, the influence of the drug-drug interaction between moxifloxacin and rifampin, namely, an increased moxifloxacin clearance, was not taken into account in these LSSs.
Therefore, the aim of this study was to develop and validate two population pharmacokinetic models of moxifloxacin (alone and with rifampin), along with clinically feasible LSSs using the Bayesian approach, as well as multiple linear regression, for the purpose of TDM of moxifloxacin in TB patients.
RESULTS
Study population.The group treated with moxifloxacin alone (MFX) included pharmacokinetic profiles of 77 TB patients, and the group with moxifloxacin together with rifampin (MFX+RIF) included profiles of 24 TB patients (Fig. 1). The baseline characteristics sex, age, and height were significantly different (P < 0.05) between these two groups (Table 1). In addition, the AUC0–24 calculated with the trapezoidal rule (AUC0–24,ref) was significantly lower, and the time of peak concentration (Tmax) was significantly earlier in the MFX+RIF group (P < 0.05, Table 2). Several abnormal pharmacokinetic curves (e.g., delayed absorption or single aberrant data point) were observed in both the MFX and the MFX+RIF group.
Moxifloxacin concentrations of the pharmacokinetic curves of MFX (n = 77) and MFX+RIF (n = 24).
Patient characteristics of the study populationa
Noncompartmental parameters of MFX and MFX+RIF
Population pharmacokinetic model.For both MFX and MFX+RIF, an one-compartment model with lag time resulted in the lowest Akaike Information Criterion (AIC) values and described the data best (Table 3). Two-compartment models were not favorable for either MFX or MFX+RIF. A statistical comparison of the pharmacokinetic parameters of the MFX versus MFX+RIF model is provided in Table 4. The total body clearance (CL) was higher, and the lag time (Tlag) was shorter in the MFX+RIF model (P < 0.05). Internal validation of the two models resulted in a mean underestimation of AUC0–24 of −5.1% (standard error [SE] 0.8%) in the MFX model and a mean underestimation of −10% (SE 2.5%) in the MFX+RIF model (Fig. 2A and Fig. 3A). In the validation of the MFX model, an r2 of 0.98, a y-axis intercept of –0.3 (95% confidence interval [CI] = –1.1 to 0.5), and a slope of 0.96 (95% CI = 0.94 to 0.98) were found in the Passing Bablok regression (Fig. 2B). For the MFX+RIF model, an r2 of 0.94, y-axis intercept of –1.0 (95% CI = −4.1 to 0.9), and a slope of 0.98 (95% CI = 0.92 to 1.07) were found in the Passing Bablok regression (Fig. 3B).
Starting parameters of the default one-compartment and two-compartment models of MFX and MFX+RIF, together with the parameters of the final models based on AIC
Comparison of pharmacokinetic parameters of the population pharmacokinetic model of MFX versus MFX+RIF
Bland-Altman plot (A) and Passing Bablok regression (B) of internal validation (n − 7) of population pharmacokinetic model of MFX (n = 77).
Bland-Altman plot (A) and Passing Bablok regression (B) of internal validation (n − 2) of population pharmacokinetic model of MFX+RIF (n = 24).
LSS using the Bayesian approach.The best performing LSSs of MFX and MFX+RIF are shown in Table 5 and Table 6, including the mean prediction error (MPE), the root mean-squared error (RMSE), and r2 to evaluate the performance of the LSSs. The performance of the LSS using t = 2 and 6 h samples was evaluated as well because this strategy is currently used in many health facilities for TDM of anti-TB drugs (31). Not all strategies met the preset acceptance criteria (RMSE < 15%, MPE < 5%).(21) Low r2 values were observed that were caused by high interindividual variability in performance of the LSSs.
LSSs of moxifloxacin without rifampin using the Bayesian approach
LSSs of moxifloxacin with RIF using the Bayesian approach
For the MFX model, an LSS using t = 0 and 6 h samples was chosen for further evaluation (RSME = 15.17%, MPE = 2.42%, r2 = 0.874) because it required one sample less than the three-sample strategies, while the RMSE was only slightly above 15%. The internal validation showed a mean underestimation of −4.8% (SE 1.3%). However, low AUC0–24 values were more frequently overestimated in contrast to an AUC0–24 of >40 mg ⋅ h/liter mainly being underestimated by the LSS (Fig. 4A). The Passing Bablok regression showed an r2 of 0.94, a y-axis intercept of 3.4 (95% CI = 1.6 to 4.9), and a slope of 0.85 (95% CI = 0.80 to 0.91) (Fig. 4B).
Bland-Altman plot (A) and Passing Bablok regression (B) of internal validation of Bayesian LSS (t = 0 and 6 h) of MFX (n = 77).
For the MFX+RIF model, an LSS using t = 0 h and 6 h samples was chosen for further evaluation (RSME = 15.81%, MPE = 2.35%, r2 = 0.885) because of the benefit of requiring only two samples while performance in terms of RSME and MPE remained acceptable. The internal validation showed a mean underestimation of −5.5% (SE 3.1%) in the Bland-Altman plot and an r2 of 0.90, a y-axis intercept of –1.3 (95% CI = –4.4 to 2.8), and a slope of 1.0 (95% CI = 0.88 to 1.10) in the Passing Bablok regression (Fig. 5).
Bland-Altman plot (A) and Passing Bablok regression (B) of internal validation of Bayesian LSS (t = 0 and 6 h) of MFX+RIF (n = 24).
LSS using multiple linear regression.Tables 7 and 8 show the best-performing LSSs for MFX and MFX+RIF. The performance of the frequently used LSS using t = 2 and 6 h samples was evaluated as well and included in the tables. None of the MFX LSSs met the acceptance criteria (RMSE < 15%, MPE < 5%) as bias was above 5% for all combinations. For MFX+RIF, the two three-sample strategies and LSS using t = 1 and 6 h samples met the acceptance criteria.
LSSs of MFX without RIF determined using linear regression
LSSs of MFX+RIF determined using multiple linear regression
The MFX LSS using t = 0 and 4 h samples (RSME = 9.25%, MPE = 6.85%, r2 =0.957) had a performance comparable to the three-sample strategies while being more clinically feasible and therefore was chosen for further evaluation. In contrast to the Bayesian LSSs for MFX and MFX+RIF, a t = 0 and 6 h strategy was not feasible using a multiple linear regression approach since its performance was substantially worse (RMSE = 12.01, MPE = 9.43, r2 = 0.905) than the LSS using t = 0 and 4 h samples. Internal validation of this t = 0 and 4 h LSS for MFX showed a mean overestimation of 0.2% (SE 1.3%) in the Bland-Altman plot and an r2 of 0.95, a y-axis intercept of 0.1 (95% CI = –2.1 to 1.6), and a slope of 0.99 (95% CI = 0.95 to 1.06) in the Passing Bablok regression (Fig. 6).
Bland-Altman plot (A) and Passing Bablok regression (B) of internal validation (n − 6) of LSS using multiple linear regression (t = 0 and 4 h) of MFX (n = 66).
For MFX+RIF, the LSS using t = 1 and 6 h samples (RSME = 6.09%, MPE = 4.83%, r2 = 0.971) was chosen for further evaluation, because of clinical suitability in addition to good performance (RMSE < 15%, MPE < 5%). Internal validation showed a mean overestimation of 0.9% (SE 2.1%) in the Bland-Altman plot and an r2 of 0.96, a y-axis intercept of −0.2 (95% CI = –4.9 to 2.3), and a slope of 1.02 (95% CI = 0.88 to 1.15) in the Passing Bablok regression (Fig. 7).
Bland-Altman plot (A) and Passing Bablok regression (B) of internal validation (n − 1) of LSS using multiple linear regression (t = 1 and 6 h) of MFX+RIF (n = 14).
DISCUSSION
In this study, we successfully developed a population pharmacokinetic model for moxifloxacin alone and in combination with rifampin. Furthermore, we developed and validated sampling strategies using the Bayesian approach (MFX and MFX+RIF t = 0 and 6 h) and multiple linear regression (MFX t = 0 and 4 h; MFX+RIF t = 1 and 6 h) for both groups as well.
It was decided to develop two separate population pharmacokinetic models and therefore also separate LSSs, for moxifloxacin alone and in combination with rifampin, after observing a significant effect of rifampin on the pharmacokinetics of moxifloxacin. The population pharmacokinetic model of MFX+RIF showed an approximately 35% higher total body clearance of moxifloxacin compared to the MFX pharmacokinetic model (Table 4). This was to be expected since rifampin enhances the metabolism of moxifloxacin and increases in total body clearance of 45 to 50% have been reported by others (10, 32). As a result of this drug-drug interaction, the pharmacokinetic profiles of MFX+RIF showed reduced moxifloxacin concentrations and 25% lower median moxifloxacin AUC0–24 values after administration of a similar dose (Fig. 1, Table 2). The latter is confirmed by a significant −17% difference in dose-corrected AUC0–24,ref between the MFX and MFX+RIF group (Table 2). The decrease in moxifloxacin exposure by rifampin was estimated at 30% in previous studies (10, 12, 32), although others found nonsignificant or smaller decreases in moxifloxacin AUC0–24 (21, 33). In this study, we observed only a slightly smaller effect of rifampin on the total body clearance and exposure than previously reported. This might be explained by the possibility that maximal enzyme induction was not yet achieved at the moment of sampling in a few cases, since it generally takes around 10 to 14 days of rifampin treatment to reach maximal induction (34). Furthermore, we encountered a significant, but small, difference in lag time between the MFX and MFX+RIF models and in the Tmax of the included pharmacokinetic profiles. Faster absorption of moxifloxacin in combination with rifampin was found in other studies as well; however, some reported the opposite effect. This could suggest that lag time and Tmax was not influenced by rifampin but more likely by other differences between the MFX and MFX+RIF group, such as concomitantly taken TB drugs or interindividual differences in absorption due to disease state.
In addition to the population pharmacokinetic models, we developed and validated LSSs using the Bayesian approach, as well as multiple linear regression for MFX and MFX+RIF. LSSs of moxifloxacin have been described before. Pranger et al. found a Bayesian LSS with a comparable performance (RMSE = 15%, MPE = –1.5%, r2 = 0.90) compared to our LSSs for MFX and MFX+RIF (21). The LSS of Magis-Escurra et al. performed better (RMSE = 1.45%, MPE = 0.58%, r2 = 0.9935) than the multiple linear regression LSSs proposed in this study (20). However, a smaller sample size (n = 12) was used to establish the equation, and this was not externally validated. Further, we provided suitable sampling strategies for multiple situations, in patients using moxifloxacin alone or together with rifampin, and for centers that either do or do not have pharmacokinetic modeling software available. Health care professionals may select the LSS that is the most applicable to the circumstances.
The Bayesian LSS for MFX (t = 0 and 6 h) showed a slight downward trend between the bias of the estimated AUC0–24 and the mean of the estimated and actual AUC0–24 (Fig. 4). Low AUC0–24 values were more frequently overestimated in comparison to higher AUC0–24 values. A possible cause might be that we could not differentiate between metabolic clearance and renal clearance in both population pharmacokinetic models due to a small range of creatinine clearance in the study population. A relatively high exposure of moxifloxacin in patients with renal insufficiency could be underestimated since renal function may be overestimated and the other way around for patients with normal renal function and relatively low exposures. The pharmacokinetic modeling software will fit a curve, with the greatest likelihood of being the actual pharmacokinetic curve based on drug concentrations at 0 and 6 h, together with patient characteristics and data of the entire population. However, when the influence of creatinine clearance is not available, the software will pick a fit with average parameters, causing overestimation in low AUC0–24 and underestimation in high AUC0–24 ranges. We decided not to validate one of the better performing three-sample strategies from Table 5, since we focused on developing a clinically feasible LSS with a strong preference for only two samples. Furthermore, we aimed to provide a simple and well-performing alternative LSS for MFX using multiple linear regression (t = 0 and 4 h). We recommend using this LSS instead of the Bayesian LSS for MFX, particularly when low drug exposure is suspected, because overestimation of AUC0–24 can lead to sub therapeutic dosing with treatment failure and acquired drug resistance as possible harmful consequence (26, 35, 36).
In this study we decided to validate one LSS for each situation (Bayesian or multiple linear regression; MFX or MFX+RIF), due to the significant influence of rifampin on the pharmacokinetics of moxifloxacin and so there would be a suitable LSS for every patient in each health care center. The LSSs using multiple linear regression performed rather well in our study population but are less flexible in patients with different characteristics. A Bayesian LSS is therefore preferred for patients who are not comparable to our study populations since the population pharmacokinetic model is able to include some patient characteristics. Clinicians are guided to the best option for TDM of moxifloxacin by following the decision tree in Fig. 8. For implementation of moxifloxacin TDM using LSSs in daily practice, it would be convenient to be able to use one sampling strategy for both MFX and MFX+RIF. This study showed that it is possible to use t = 0 and 6 h samples in a Bayesian LSS for both MFX and MFX+RIF, and probably even in a multiple linear regression LSS for MFX+RIF, after successful validation. Unfortunately, a multiple linear regression strategy for MFX alone using t = 0 and 6 h samples was not feasible because of inferior performance. Considering that TB patients are treated with a combination of multiple anti-TB drugs, one single LSS suitable for all drugs of interest is the ideal situation but, unfortunately, also rather challenging due to the various pharmacokinetic properties of the different drugs. Others did succeed in developing a LSS using multiple linear regression for simultaneously estimating exposure of all first-line drugs and moxifloxacin in a small population of TB patients (20). A 2 and 6 h postdose sampling strategy is frequently used for TDM of anti-TB drugs since it is believed to be able to estimate Cmax, as well as to detect delayed absorption (31). However, better performances were found for the LSSs proposed in this study, although the 2 and 6 h LSSs performed within acceptable limits as well in the Bayesian approach and the multiple linear regression.
Clinical guide for choosing the best LSS for TDM of moxifloxacin alone or in combination with rifampin.
In general, we noticed large interindividual pharmacokinetic variation in terms of moxifloxacin concentrations (Fig. 1), Cmax, and AUC0–24 (Table 2) as described earlier (18), but also in Ka and CL/F (Table 4). Patients received 400, 600, or 800 mg of moxifloxacin; this obviously influenced the drug concentration, Cmax, and the AUC0–24, but not all variation could be explained by different dosage regimes. For MFX the AUC0–24 corrected to a 400-mg standard dose ranged from 10.2 to 79.1 mg ⋅ h/liter, and for MFX+RIF the AUC0–24 corrected to a 400-mg standard dose ranged from 10.0 to 47.4 mg ⋅ h/liter. This substantial interindividual variation is the reason why TDM of moxifloxacin is helpful to ensure optimal drug exposure and thus minimize the risk of treatment failure and developing acquired drug resistance (26, 27). The estimated AUC0–24 using one of the LSSs proposed, together with the MIC of the M. tuberculosis strain, will provide valuable information on the optimal moxifloxacin dose to be used in an individual patient.
A limitation of the study is the exclusion of the creatinine clearance from the population pharmacokinetic model. As discussed earlier, this could have led to the observed bias in the MFX LSS using 0 and 6 h samples since approximately 20% of moxifloxacin is eliminated unchanged in the urine. On the contrary, a well-performing LSS using multiple linear regression (t = 0 and 4 h) is a suitable alternative for MFX. The lack of prospective or external validation of the population pharmacokinetic model and LSSs could be considered another limitation. However, we were able to collect a large data set to develop the model and clinically feasible LSSs using a sufficient number of pharmacokinetic profiles. A strength of our study is that a large part of our dataset consisted of drug concentrations which were collected as part of daily routine TDM. During visual check of the data we noticed several abnormal curves (both MFX and MFX+RIF) that, for instance, showed delayed absorption with Tmax values of 4 to 6 h. These curves were not excluded from the study. The models and LSSs appeared to be able to adapt to this delayed absorption. In most cases, the subsequent decision to either increase the dose or not was similar. For these reasons, we expect the results reported here to represent the clinical practice of TDM using these LSSs very closely. The small sample size of the MFX+RIF group can be considered a limitation as well, although comparable to previously published LSS studies (21, 37–40). We consider this sample size to be sufficient for exploratory objectives, since this is the first study that developed separate LSSs for moxifloxacin alone and in combination with rifampin. Future research can build on the results described in this study.
In conclusion, we developed and validated two separate pharmacokinetic models for moxifloxacin alone and in combination with rifampin in TB patients. We provided data to show significant differences in drug clearance and drug exposure between these groups. Furthermore, we developed and validated LSS based on the Bayesian approach (MFX and MFX+RIF, 0 and 6 h) and multiple linear regression (MFX, 0 and 4 h; MFX+RIF, 1 and 6 h) that can be used to perform TDM on moxifloxacin in TB patients.
MATERIALS AND METHODS
Study population.This study used three databases. Database 1 consisted of retrospective data of routine TDM in 67 tuberculosis patients treated at Tuberculosis Center Beatrixoord, University Medical Center Groningen, Groningen, The Netherlands, and was collected between January 2006 and May 2017, partly published earlier (18). All patients received moxifloxacin (with or without rifampin) as part of their daily TB treatment and pharmacokinetic curves were obtained as part of routine TDM care. Each patient was only included once. Various sampling schedules were used, but most profiles included t = 0, and 1, 2, 3, 4, and 8 h postdose samples. Pharmacokinetic profiles consisting of less than three data points were excluded. The second database included data of 25 TB patients participating in a clinical study in Thessaloniki, Greece (33). After at least 12 days of treatment with moxifloxacin with or without rifampin, blood samples were collected at t = 0 and at 1, 1.5, 2, 3, 4, 6, 9, 12, and 24 h after drug intake. The third database consisted of pharmacokinetic data of nine Brazilian TB patients receiving 400 mg of moxifloxacin (no rifampin) daily in an early bactericidal activity study.(14) At day 5, blood samples were collected at t = 0 and at 1, 2, 4, 8, 12, 18 and 24 h after drug intake.
As steady state is reached within 3 to 5 days of treatment with moxifloxacin, all data were collected during steady-state conditions (11). In general, no informed consent was required, due to the retrospective nature of the study.
The total study population was split in two groups—patients that received moxifloxacin alone (MFX) and patients that received moxifloxacin together with rifampin (MFX+RIF)—because of the pharmacokinetic drug-drug interaction between rifampin and moxifloxacin (10). Since sample collection in the MFX+RIF group was performed after a median number of days on rifampin treatment of 35 (interquartile range [IQR], 13 to 87 days), maximum enzyme induction by rifampin was expected to be reached in most patients (41).
Patient characteristics of both groups were tested for significant differences, median (IQR) using the Mann-Whitney U test, and number (%) using the Fisher exact test in IBM SPS Statistics (version 23; IBM Corp., Armonk, NY). P values of <0.05 were considered significant.
Population pharmacokinetic model.For each group, MFX and MFX+RIF, a population pharmacokinetic model was developed using the iterative two-stage Bayesian procedure of the KinPop module of MWPharm (version 3.82; Mediware, The Netherlands). Since the pharmacokinetics of moxifloxacin have been described with one compartment (14, 21). as well as two-compartment models (42, 43), both types were evaluated. The population pharmacokinetic parameters of the models were assumed to be log normally distributed, with a residual error and concentration-dependent standard deviation (SD; SD = 0.1 + 0.1 × C, where C is the moxifloxacin concentration in mg/liter). Because the bioavailability (F) of moxifloxacin is almost complete (11) and pharmacokinetic data following intravenous administration were not available, F was fixed at 1 in the analysis and pharmacokinetic parameters are presented relative to F. Moxifloxacin is mainly metabolized in the liver by glucuronosyltransferase and sulfotransferase (ca. 80%) (11). Only total body clearance (CL), the sum of metabolic and renal clearance, was included in the model development because it was not possible to determine renal clearance due to a small range of creatinine clearance values in our data set.
We started the analysis with a single default one-compartment model for both MFX and MFX+RIF developed by Pranger et al. using a very similar methodology (21). This study found comparable pharmacokinetic parameters of MFX and MFX+RIF, although likely due to a small sample size. Two default two-compartment models were used, one for MFX and one for MFX+RIF (42, 44). Modeling was started with all parameters fixed, and the Akaike Information Criterion (AIC) was used to evaluate the model (45). Subsequently, one by one parameters were Bayesian estimated, and each step was evaluated by calculation of the AIC. A reduction in the AIC with at least three points was regarded as a significant improvement of the model (46). One-compartment models included the parameters CL, volume of distribution (V), and absorption rate constant (Ka). Two-compartment models included the parameters Ka, CL, the intercompartmental clearance (CL12), the central volume of distribution (V1), the volume of distribution of the second compartment (V2), and the lag time for absorption (Tlag). Afterward, Tlag was added to the best performing one compartment model and evaluated for goodness of fit as well because of oral intake of moxifloxacin. The default two-compartment models already included Tlag. The final models of MFX and MFX+RIF were chosen based on AIC values.
The final models were internally validated using 11 different (n – 7) submodels for MFX and 12 (n – 2) submodels for MFX+RIF, each leaving out randomly chosen pharmacokinetic curves. All pharmacokinetic curves were excluded once (jackknife analysis). The Bayesian fitted AUC0–24 of each left out curve (AUC0–24,fit) was compared to the AUC0–24 calculated with the trapezoidal rule (AUC0–24,ref) using a Bland-Altman plot and Passing Bablok regression (Analyze-it 4.81; Analyze-it Software Ltd., Leeds, United Kingdom). In the calculation of AUC0–24,ref, moxifloxacin concentrations at t = 0 and 24 h after drug intake were assumed to be equal due to steady-state conditions. The Cmax (mg/liter) was defined as the highest observed moxifloxacin concentration and Tmax (h) as the time at which Cmax occurred. Noncompartmental parameters (i.e., AUC0–24,ref, dose-corrected AUC0–24,ref to the standard dose of 400 mg, Cmax, and Tmax) and population pharmacokinetic model parameters of the MFX and MFX+RIF group were compared and tested for significant differences using the Mann-Whitney U test.
LSS using Bayesian approach.Using the Bayesian approach, we performed two separate analyses to develop LSSs: one for MFX and one for MFX+RIF. Using Monte Carlo simulation in MWPharm, 1,000 virtual pharmacokinetic profiles were created to represent the pharmacokinetic data used in the development of the LSS. The reference patient for the Monte Carlo simulation was selected based on representative pharmacokinetic data and patient characteristics. For MFX, a 36-year-old male with a bodyweight of 57 kg, a height of 1.60 m, a body mass index (BMI) of 22.2 kg/m2, a serum creatinine level of 74 μmol/liter, and a moxifloxacin dose of 7.0 mg/kg was chosen. For MFX+RIF, a 56-year-old male with a body weight of 56 kg, a height of 1.63 m, a BMI of 21.1 kg/m2, a serum creatinine level of 80 μmol/liter, and a moxifloxacin dose of 7.1 mg/kg was selected. The LSSs were optimized using the steady-state AUC0–24. Only clinically feasible LSSs using one to three samples between 0 and 6 h postdose and a sample interval of 1 h were tested. The LSSs were evaluated using acceptance criteria for precision and bias (RMSE < 15%, MPE < 5%) (18). For both MFX and MFX+RIF, one LSS was chosen for internal validation based on performance, as well as clinical feasibility. The AUC0–24 estimated with the chosen LSS (AUC0–24,est) was compared to the AUC0–24,ref using a Bland-Altman plot and Passing Bablok regression. In addition, the performance of a LSS using 2 and 6 h postdose samples was evaluated because this is LSS is frequently used for TDM of anti-TB drugs (31).
LSS using multiple linear regression.Two separate analyses (MFX and MFX+RIF) using multiple linear regression were performed.
Only clinically suitable LSSs (one to three samples, 0 to 6 h postdose, and sample interval of 1 h) were included in the analysis. Each analysis excluded the pharmacokinetic curves without data at the selected time points of the LSS, resulting in a variable number of included curves (N). Multiple linear regression in Microsoft Office Excel 2010 was used to evaluate the correlation of moxifloxacin concentrations at the chosen time points of the LSS and AUC0–24,ref. The acceptance criteria (RMSE < 15%, MPE < 5%) were applied to each LSS (18). Internal validation using 11 different (n – 6) subanalyses for MFX and 14 (n – 1) subanalyses for MFX+RIF was used to evaluate the performance of the LSSs. Each subanalysis excluded randomly chosen profiles, and all profiles were excluded once (jackknife analysis). Agreement of AUC0–24,est and AUC0–24,ref was tested using a Bland-Altman plot and Passing Bablok regression.
ACKNOWLEDGMENT
No funding was received for the present study. The Brazilian TBRU moxifloxacin study was funded by the U.S. National Institutes of Health (NO1-AI95383 and HHSN266200700022C).
FOOTNOTES
- Received 20 February 2019.
- Returned for modification 26 March 2019.
- Accepted 15 April 2019.
- Accepted manuscript posted online 22 April 2019.
- Copyright © 2019 American Society for Microbiology.