Assessing Prolongation of the Corrected QT Interval with Bedaquiline and Delamanid Coadministration to Predict the Cardiac Safety of Simplified Dosing Regimens.

Delamanid and bedaquiline are two drugs approved to treat drug-resistant tuberculosis, and each have been associated with corrected QT interval (QTc) prolongation. We aimed to investigate the relationships between the drugs' plasma concentrations and the prolongation of observed QT interval corrected using Fridericia's formula (QTcF) and to evaluate their combined effects on QTcF, using a model-based population approach. Furthermore, we predicted the safety profiles of once daily regimens. Data were obtained from a trial where participants were randomized 1:1:1 to receive delamanid, bedaquiline, or delamanid + bedaquiline. The effect on QTcF of delamanid and/or its metabolite (DM-6705) and the pharmacodynamic interactions under coadministration were explored based on a published model between bedaquiline's metabolite (M2) and QTcF. The metabolites of each drug were found to be responsible for the drug-related QTcF prolongation. The final drug-effect model included a competitive interaction between M2 and DM-6705 acting on the same cardiac receptor and thereby reducing each other's apparent potency, by 28% (95% confidence interval (CI), 22-40%) for M2 and 33% (95% CI, 24-54%) for DM-6705. The generated combined effect was not greater but close to "additivity" in the analyzed concentration range. Predictions with the final model suggested a similar QT prolonging potential with simplified, once-daily dosing regimens compared with the approved regimens, with a maximum median change from baseline QTcF increase of 20 milliseconds in both regimens. The concentrations-QTcF relationship of the combination of bedaquiline and delamanid was best described by a competitive binding model involving the two main metabolites. Model predictions demonstrated that QTcF prolongation with simplified once daily regimens would be comparable to currently used dosing regimens.

WHAT IS THE CURRENT KNOWLEDGE ON THE TOPIC?

Bedaquiline and delamanid use has been associated with corrected QT interval (QTc) prolongation. Combining the two drugs led to a no more than additive effect on the QTc interval.

WHAT QUESTION DID THIS STUDY ADDRESS?

What is the relation between drugs' concentrations and QTc prolongation?

WHAT DOES THIS STUDY ADD TO OUR KNOWLEDGE?

The respective metabolites were found to be responsible for the QTcF prolongation. The pharmacodynamic interaction was described by a competitive agonist model, indicating no synergic toxicity for the combination of bedaquiline and delamanid.

HOW MIGHT THIS CHANGE CLINICAL PHARMACOLOGY OR TRANSLATIONAL SCIENCE?

Alternative dosing regimens can be explored with this quantitative model‐based approach. The newly more‐practical‐to‐use once‐daily regimen for these drugs was explored through the developed model, and was predicted not to lead to a higher risk of safety events.

Mycobacterium tuberculosis is one of the leading causes of mortality globally, with about 1.5 million reported deaths in 2019. Drug‐resistant tuberculosis (TB) has poor treatment outcomes and, until recently, required prolonged therapy, including injectables, with considerable toxicity. Fortunately, since 2020, all‐oral shorter regimens (9–12 months) have been recommended by the World Health Organization (WHO) for the treatment of drug‐resistant TB. To further refine treatment, mismatching dosing schedules could be aligned within the multidrug regimens and once‐daily dosing is preferred to simplify use and improve adherence.

Bedaquiline and delamanid are two registered antituberculosis drugs with new and distinct mechanisms of action, making them good candidates for inclusion in shorter, all‐oral, once‐daily regimens. With bedaquiline targeting adenosine triphosphate synthase and delamanid acting on mycolic acid synthesis, the combination shows promising results in terms of efficacy to enhance culture conversion. , , , However, both drugs have been associated with the cardiotoxic potential of prolonging the QTc interval, a risk factor for sudden death. , While no pharmacokinetic interaction between bedaquiline and delamanid is anticipated, the pharmacodynamic interaction requires further study. Prior observational studies have been reassuring with regard to QTc prolongation when bedaquiline and delamanid are given together under compassionate use. , , , , The DELIBERATE (Delamanid Bedaquiline for Resistant Tuberculosis) clinical trial was a phase II randomized, controlled study that aimed to precisely quantify the effect of bedaquiline, delamanid, or both on the QTcF interval. Given the long terminal half‐lives of the studied TB drugs and their metabolites, leading to slow accumulation over time, the effect on QTcF interval was assessed over 6 months of multidrug treatment. The trial included 84 participants, of whom 80 proceeded to receive study treatment, and none of whom experienced absolute QTcF values over 500 milliseconds. The change from baseline QTcF of the combined arm vs. the sum of changes on bedaquiline alone plus delamanid alone arms had a mean QTcF of −0.1 milliseconds (95% CI, −8.0 to 7.7). The authors concluded that combining bedaquiline and delamanid should have a no more than additive effect on the QTc interval.

Our objective in this analysis was to characterize the relationships between plasma concentrations of bedaquiline, its M2 metabolite, delamanid, its DM‐6705 metabolite, and QTcF interval, using a population pharmacokinetic–pharmacodynamic model. Then we simulated the effect of daily dosing regimens of bedaquiline and delamanid on QTcF interval prolongation, at doses proposed for ongoing and future clinical trials.

MATERIALS AND METHODS

Data

Data were obtained from the DELIBERATE clinical trial, which was a phase II, open‐label, randomized controlled study in adults with multidrug‐resistant TB (ClinicalTrials.gov number NCT 02583048). The trial was approved by local ethics committees at each site, and participants gave written informed consent.

Participants were randomized 1:1:1 to either bedaquiline, delamanid, or bedaquiline plus delamanid, together with a multidrug background treatment (28 participants in each group). Clofazimine and moxifloxacin, both of which cause significant QT prolongation, were not permitted in the background regimen. Delamanid was dosed 100 mg twice daily with food for 24 weeks, and bedaquiline was administered 1 hour after the meal together with other anti‐TB drugs at 400 mg daily for 2 weeks then 200 mg three times a week for 22 weeks. QT intervals were assessed by electrocardiogram (ECG) and corrected for the effect of heart rate using Fridericia's formula (QTcF). ECGs were performed in triplicate at baseline and every 2 weeks until Week 24, ~ 4–6 hours post dose, before lunch. Rich pharmacokinetic (PK) sampling was performed at Weeks 2, 8, and 24, (according to the schedule displayed in Table ), and sparse sampling was drawn, concurrent to ECG measurements, every 2 weeks.

Model selection and evaluation

Model selection was based on significant drops in objective function value (OFV) or Bayesian information criteria (BIC), goodness‐of‐fit plots, precision in parameter estimates, and scientific plausibility. For nested models, a significant drop in OFV for one added parameter (i.e., 1 degree of freedom) was 3.84 points, corresponding to a P value <0.05. Non‐nested models were compared with BIC for nonlinear mixed‐effect models, accounting for the numbers of subjects, observations, fixed effects, and random effects in the calculation of the OFV's penalty. Parameter precision was obtained with the R covariance (Hessian) matrix. Models were qualified using visual predictive checks (1,000 replicates), stratified by study arms. All triplicate ECG measurements were included (i.e., no averaging of data was carried out). Missing covariate values (see Table ) were imputed by replacing them with the population median for continuous covariates, and with the most common category for categorical covariates.

Table 1

Demographic characteristics at baseline and model‐predicted PK data at ECG timepoints for each analyte

	Bedaquiline alone arm (n = 28)	Delamanid alone arm (n = 26)	Delamanid + bedaquiline arm (n = 26)	Total (n = 80)
Age (years)
Median (min, max)	34.5 (20, 58)	35 (18, 73)	35.5 (18, 55)	35 (18, 73)
Sex
Male	22 (79%)	19 (73%)	20 (77%)	61 (76%)
Female	6 (21%)	7 (27%)	6 (23%)	19 (24%)
Race
Black African	13 (46%)	9 (35%)	11 (42%)	33 (41%)
Mixed race	11 (40%)	11 (42%)	14 (54%)	36 (45%)
White	0 (0%)	1 (4%)	0 (0%)	1 (1%)
Othera	0 (0%)	1 (4%)	0 (0%)	1 (1%)
Missing	4 (14%)	4 (15%)	1 (4%)	9 (11%)
HIV‐1
Negative	18 (64%)	14 (54%)	15 (58%)	47 (59%)
Positive	10 (36%)	10 (38%)	10 (38%)	30 (38%)
Missing	0 (0%)	2 (8%)	1 (4%)	3 (4%)
Baseline weight (kg)
Median (min, max) [missing n (%)]	54.00 (35, 80) [1 (4%)]	54.00 (38, 83) [2 (8%)]	52.35 (41, 71.5)	54 (35, 83) [3 (4%)]
Baseline potassium (mmol/L)
Median (min, max) [missing n (%)]	4.1 (3, 5.1) [1 (4%)]	4.2 (3.6, 5.5) [2 (8%)]	4.3 (3.3, 5.5) [1 (4%)]	4.15 (3, 5.5) [4 (5%)]
Baseline QTcF (ms) ‐ Mean of triplicates
Median (min, max)	394.3 (360.7, 461)	409.2 (364.3, 445)	389.2 (368.3, 422.3)	397.5 (360.7, 461)
Bedaquiline concentrations (ng/mL) – model‐predicted individual concentrations at all ECG timepoints
Median (10th perc., 90th perc.)	1,034 (673.8, 1759)	—	1,185 (530.6, 1,690)	1,100 (627.2, 1734)
M2 concentrations (ng/mL) – model‐predicted individual concentrations at all ECG timepoints
Median (10th perc., 90th perc.)	204.0 (119.7, 331.9)	—	190.2 (109.1, 365.0)	196.4 (113.9, 344.6)
Delamanid (ng/mL) – model‐predicted individual concentrations at all ECG timepoints
Median (10th perc., 90th perc.)	—	263.4 (174.5, 336.6)	254.2 (165.6, 357.7)	256.3 (172.7, 355.2)
DM‐6705 concentrations (ng/mL) – model‐predicted individual concentrations at all ECG timepoints
Median (10th perc., 90th perc.)	—	65.6 (25.39, 118.3)	70.0 (28.32, 125.2)	67.4 (26.27, 121.6)

"Other" indicates as being neither Black African, nor Mixed race, nor White.

ECG, electrocardiogram; max, maximum; min, minimum; perc., percentile; PK, pharmacokinetic; QTcF, QT interval corrected using Fridericia’s formula; —, not applicable.

Pharmacokinetic models

Parent and metabolite drug concentrations at the time of the ECG measurements for each individual were obtained by predictions with previously developed models, based on individual dosing information, demographics, and observed drug concentrations. , Diagnostic plots between observed concentrations and individual predicted concentrations were generated to assess the predictability of the models. Similar to the previously published evaluation of bedaquiline and/or M2 effects on delamanid/DM‐6705 PK, the effect of delamanid and/or DM‐6705 on bedaquiline/M2 PK was investigated by evaluating a possible arm effect on disposition parameters.

Pharmacokinetic–pharmacodynamic model

A previously developed model describing the relation between M2 concentration and QTcF interval was used as the base model. In addition to the drug effect (DE), this model included several components that influenced the QTcF interval, such as participants' characteristics (age, sex, race, electrolyte levels) on baseline QTcF, the circadian rhythm, and the effect of time on treatment (see Figure ). No covariate search was performed on top of the base model, given that no new potential covariates were available in the analysis data set, compared with the base model.

Our analysis consisted of three steps. First, with data from the bedaquiline arm only, the base model was fitted using priors. Given the similarity of the participants' population, the PRIOR functionality in NONMEM uses information from the base model (developed with 430 participants) to stabilize the estimation process of a smaller data set, such as the current data set with 28 participants in the bedaquiline arm. The drug QTcF effect was also confirmed to be driven only by M2 concentrations, and not bedaquiline concentrations. Secondly, the impact of delamanid and/or DM‐6705 concentrations on QTcF prolongation was analyzed on top of the base model, via step, linear, or maximal drug effect (Emax) functions. Data from the delamanid arm only were first fitted to define delamanid/DM‐6705 drug effect, then applied to the entire data set. Thirdly, pharmacodynamic drug interactions between bedaquiline and delamanid, or their metabolites, using the entire data set, were explored via an empirical approach by estimating an interaction parameter β ( ) and a more mechanistic approach using a competitive interaction model ( ). , DE stands for drug effect. A represents bedaquiline or M2, and B represents delamanid or DM‐6705.

The competitive interaction model relies on the assumption that the analytes bind to a common receptor and have a common mechanism of action to induce QT prolongation. DE stands for drug effect, Emax for maximal drug effect, Conc for concentration, and EC50 for the concentration needed to achieve half of the Emax. A represents bedaquiline or M2, and B represents delamanid or DM‐6705.

Tested interindividual variability (IIV) on added parameters was assumed to be log‐normally distributed, with the exception of the interaction parameter β where a normal distribution was tested, so that positive and negative values were permitted.

To account for a possible relation between albumin concentrations and potency, through parent drug or metabolite‐binding, a directly proportional relation between albumin concentration and potency was tested. This relation can be motivated as all four entities explored have a high protein binding. ,

Simulations

The final model was used to simulate the drug‐induced QTcF interval prolongation expected if bedaquiline and delamanid were to be administered together with once‐daily regimens currently tested in clinical trials (200 mg once daily (q.d.) for 8 weeks then 100 mg q.d. for bedaquiline, and 300 mg q.d. for delamanid). This analysis focused only on the drug effect component of the model (different from the total change from baseline QTcF interval, as already reported in the clinical reference publication ), to more precisely address the change in dosing regimens. One thousand profiles of the drug‐induced QTcF prolongation were generated, considering IIV in the pharmacokinetic and pharmacodynamic parameters. The participants’ characteristics for the simula were 35 years old, 54 kg, non‐Black race, corresponding to the typical individual in the DELIBERATE trial. The PK model of delamanid/DM‐6705 includes an adherence component and a morning/evening effect on relative bioavailabity that were not included in the simulations (centering the simulations around the daily average relative bioavailability). The nonlinear bioavailability of delamanid with increasing doses was accounted for by a fixed effect on bioavailability extracted from the literature.

Software

Data analysis was performed with the nonlinear mixed‐effects modeling software NONMEM version 7.4.4 (ICON plc, Gaithersburg, MD), with the first‐order conditional estimation method with interaction. Model diagnostics were supported by Perl‐speaks‐NONMEM (PsN, version 5.2.6) , and R software version 4.1.0 (R Core Team, Vienna, Austria).

RESULTS

A total of 80 participants with drug‐resistant TB were included in the analysis, out of the 84 planned per protocol. Four participants were excluded from the analysis as they did not receive any doses of study treatment. The demographic characteristics at baseline are summarized in Table .

The previously published PK models displayed good performance at predicting observed concentrations (Figure ). The summary of the model‐predicted PK data at ECG timepoints for each analyte can be found in Table .

Figure 1

Diagnostic plots of the pharmacokinetic models showing observed concentrations vs. individual predicted concentrations for metabolites M2 and DM‐6705. The gray dashed line represents the trendline across the data and the black full line represents the line of identity.

No PK drug–drug interactions were identified between bedaquiline and delamanid: neither bedaquiline on delamanid PK, nor delamanid on bedaquiline PK (see Table ).

As a first step, the evaluation of the base model with bedaquiline arm data, using the PRIOR functionality, resulted in parameter estimates and precision in line with published results (see Table and Ref. ). The exploration of bedaquiline's contribution to the drug effect resulted in a nonsignificant improvement of the data fit, either as a single agent or as a combination with M2 (see Table ). Those results confirmed the previously published relationship.

Table 2

Parameters estimates and uncertainty of the final model

Submodel	Parameters (unit)	Value (RSE%)	IIV %CV (RSE%)
Baseline	QTcF0 (ms)	401 (0.312)	3.73 (7.52)
Drug effect	Emax a (ms)	25.9 (14.4)
	EC50,M2 (ng/mL)	695 (30.1)	155 (13.6)
	EC50,DM‐6705 (ng/mL)	205 (40.9)
Time effect	QTmax (ms)	7.09 (8.72)	166b (10.3)
	T1/2 (weeks)	7.52 (13)
Circadian rhythm	A24 (ms)	2.96 (40.9)
	φ_24 (hours)	4.76 (27.3)
	A12 (ms)	1.51 (25.6)
	φ_12 (hours)	4.32 (24.1)
Covariates	Effect of potassium levels (ms per IU/L)c	−1.25 (36.2)
	Effect of being a female (ms)	6.67 (21)
	Effect of being black (ms)	−7.14 (18.6)
	Effect of age (ms per year)d	0.366 (15)
Residual error model	Additive RUV (ms)	8.77 (3.98)	20.0 (8.9)
	Box‐Cox IIVe RUV	3.89 (21.5)
	Additive RUVrepl (ms)	5.29 (3.13)	23.5 (4.96)
	Box‐Cox IIVe RUVrepl	0.874 (34.7)

CV is reported as the square root of the variance. RSE of IIV and RUV is reported on the approximate standard deviation scale (standard error/variance estimate)/2.

A12, amplitude for the 12‐hour circadian rhythm cycles; A24, amplitude for the 24‐hour circadian rhythm cycles; CV, coefficient of variation; EC50, concentration needed to achieve half of Emax; EC50,DM‐6705, EC50 of DM‐6705, delamanid’s metabolite; EC50,M2, EC50 of M2, bedaquiline’s metabolite; Emax, maximal drug effect; IIV, interindividual variability; ms, milliseconds; QTcF0, baseline QT interval corrected using Fridericia’s formula; QTmax, maximal effect of time on treatment; RSE, relative standard error; RUV, residual unexplained variability; RUVrepl, replicate‐specific residual unexplained variability; T1/2, time needed to achieve half of QTmax; φ12, acrophase for the 12‐hour circadian rhythm cycles; φ24, acrophase for the 24‐hour circadian rhythm cycles.

Same maximal effect parameter (Emax) for M2 and DM‐6705.

IIV coded with a proportional model, whereas the others are coded with an exponential model.

Absolute change in QTcF0 (ms) per IU/L, different from the population median, 4.150 IU/L.

Absolute change in QTcF0 (ms) per year, different from the population median, 35 years.

Parameter estimate of the Box–Cox transformed distribution of IIV on ε components.

Secondly, analysis of the contribution of delamanid and DM‐6705 concentrations to the drug effect, with delamanid arm data only, demonstrated that only DM‐6705 (and not delamanid) played a significant role in QTcF prolongation. A model accounting for the effects of both delamanid and DM‐6705 concentrations was not better at describing the data. The function that described best the drug‐induced QTcF prolongation by DM‐6705 concentrations in the delamanid arm data was a linear function (see Table ), with a slope coefficient estimated to 7.91 milliseconds per 100 ng/mL of DM‐6705 (90% CI, 4.2 to 11.6). Those results are in line with a published linear model, defined with pediatric data, where the slope was estimated to 6.13 milliseconds per 100 ng/mL (90% CI, 1.6 to 10.7).

Thirdly, the characterization of the interaction between M2 and DM‐6705 concentrations to the drug‐induced QTcF prolongation was explored for the entire data set (80 individuals). Beforehand, the previously defined Emax function for M2 concentrations, and the linear function for DM‐6705 concentrations, were turned into two Emax functions sharing the same maximal effect parameter (no difference in OFV and same number of parameters). This allowed us to normalize the drug effect (by dividing by the maximal effect) and to compare the empirical approach with the Bliss independence theorem. The latter says that the interaction parameter β = 1 indicates Bliss independence (i.e., additive interaction), β > 1 antagonism and β < 1 synergy. In our analysis, β was estimated to 2.72 (relative standard error: 32%), indicating antagonism. Also, a competitive interaction model was tested, with and without estimation of two separated maximal effects. Finally, the competitive interaction model with a common estimated maximal effect (Eq. 3) was selected, over the empirical more flexible model, to be the best to describe the data (BICcompetitive = 15,491; BICempirical = 15,497). Addition of IIV on Emax or EC50,DM‐6705 did not improve the fit of the data. Further investigations to account for a possible relation between albumin concentrations and potency were not significant.

Visual predictive checks of QTcF vs. time and change from baseline QTcF vs. time, showing adequate model fit to the data are depicted in Figure . The final model parameter estimates with their precision are summarized in Table . To note, although IIV parameters on maximal effect of time on treatment and EC50,M2 were large (166% and 155%, respectively), the uncertainty of these parameters was low (10.3% and 13.6%, respectively).

Figure 2

Visual predictive checks of the final models. Panel (a) represents QTcF interval over time after start of treatment per arm, and panel (b) represents change from baseline QTcF interval over time after start of treatment per arm. The solid and dashed lines represent the median, the 2.5th, and 97.5th percentiles of the observed data (black circles), respectively, and the shaded areas the simulation‐based 95% confidence intervals for the corresponding percentiles. QTcF, QT interval corrected using Fridericia’s formula.

Simulation of once‐daily regimens

The final model predicted similar QTcF prolongation with the approved regimens and the once‐daily regimens, as depicted in Figure . It predicted median drug‐induced QTcF increases at 2, 8, and 24 weeks of 8.4 milliseconds (95% CI, 4.2–22.0), 12.5 milliseconds (95%, CI, 6.1–24.5), and 11.9 milliseconds (95% CI, 6.1–24.0) with once‐daily regimens compared with 10.9 milliseconds (95% CI, 5–24.4), 12.1 milliseconds (95% CI, 6.2–24.5), and 12.4 milliseconds (95% CI, 6.2–24.2) with the approved regimens. At Week 24, 15.2 and 18.3% of the participants were simulated to have a change from baseline QTcF interval over 30 milliseconds with the once‐daily and approved regimens, respectively. Similarly, 4.1% of the participants with both regimens were simulated to have an absolute QTcF interval over 450 milliseconds at Week 24.

Figure 3

Simulated drug‐induced QTcF increase with approved regimens (400 mg daily for 14 days, then 200 mg thrice‐weekly for bedaquiline, and 100 mg twice‐daily for delamanid) and once‐daily regimens (200 mg daily for 8 weeks then 100 mg daily for bedaquiline, and 300 mg daily for delamanid) for a typical participant (35 years old, 54 kg, non‐Black). The solid line represents the median of the simulated data, and the limits of the shaded area represent the 2.5th and 97.5th percentiles of the simulated data. QTcF, QT interval corrected using Fridericia’s formula.

DISCUSSION

In this work, we used a pharmacometric approach to characterize the QTcF prolongation and its relation to parent drug/metabolite exposures when bedaquiline and delamanid are administered together. Indeed, QTcF interval prolongation has been reported as an electrophysiologic effect of both drugs separately. To perform this work, we analyzed data from a trial that was designed to assess the pharmacodynamic interaction when both drugs were administered alone or together (on top of multidrug background treatment). First, we confirmed the absence of a pharmacokinetic interaction between bedaquiline and delamanid, i.e., bedaquiline was not impacting delamanid's pharmacokinetics, and delamanid was not impacting bedaquiline's pharmacokinetics. Then, after characterizing the effect of each drug individually and confirming the metabolites to be the drivers of the QTcF prolongations, the pharmacodynamic interaction between the two metabolites was explored. A flexible model allowing the interaction to be either additive, antagonistic, or synergistic, when compared with the Bliss independence theorem, indicated an antagonistic relationship between the metabolites. Thus, a competitive binding model was tested and selected as being the best to describe the data. It denotes that metabolites M2 and DM‐6705 both act as agonists of the human ether‐a‐go‐go‐related gene (hERG) receptor (by binding and blocking the hERG potassium channel) , but will act as antagonists for the other metabolite by competing to bind to the same receptor. Therefore, each metabolite reduces the other metabolite's apparent potency (i.e., increasing apparent EC50, see Eq. 3). It results in a shift of the concentration–response curve to the right, meaning that a larger concentration of drug is required to give the same prolongation of the QTcF interval in the presence of the other metabolite. The ratios of median concentration to typical EC50 for each metabolite were of the same magnitude (0.28 for M2 and 0.33 for DM‐6705), depicting a similar extent of effect for both metabolites on each other's potency, in the range of analyzed concentrations (Figure ). Moreover, the low values of the ratios and the near‐parallel plot lines (Figure ) indicate that the combined effect was close to “additivity,” where the combined effect would be equal to the sum of two effects separately. However, the time course of each metabolite is different, so the contribution to the drug effect of each metabolite over time is illustrated in Figure .

Figure 4

Drug‐induced QTcF increase vs. metabolite M2 or DM‐6705 concentrations, stratified by arm. For the bedaquiline + delamanid arm (black line), for each panel, while the concentration of one metabolite is increasing, the concentration of the other metabolite is constant (set to median of observed concentrations). QTcF, QT interval corrected using Fridericia’s formula.

Figure 5

Contribution of each metabolite (M2 or DM‐6705) to the drug‐induced QTcF increase over time after start of study. The solid line represents the median contribution among all participants and the dashed lines represent the individual contributions. QTcF, QT interval corrected using Fridericia’s formula.

The newly tested once‐daily regimens for bedaquiline and delamanid are of key importance to improve patients' quality of life as well as logistics for programs for provision of multidrug regimens. To assess the cardiotoxic risk of the newly offered once‐daily regimens to TB patients, we performed simulations with the developed model to compare the drug‐induced QTcF prolongation between the approved regimens and the newly proposed once‐daily regimens (200 mg q.d. for 8 weeks then 100 mg q.d. for bedaquiline, and 300 mg q.d. for delamanid). The simulations suggested that once‐daily regimens of bedaquiline and delamanid taken together are as safe as the approved regimens. Indeed, with the once‐daily regimens, the predicted maximal QTcF increase over 24 weeks of treatment, as well as the prevalence of patients with change from baseline QTcF interval over 30 milliseconds at 24 weeks, were not higher than the predicted values with the current approved regimens. From a safety perspective, these results provide assurance in the use of all‐oral once‐daily regimens.

There are several limitations to the model developed. One is that the estimation of two separate Emax parameters for M2 and DM‐6705 in the competitive interaction model was not supported by the data. It would mean that at high concentrations (higher than the analyzed concentrations range), the two metabolites would share the same maximal possible toxicity. So, extrapolation outside the range of analyzed exposures could be uncertain; however, in the analysis, the simulated concentrations were in the same range as the concentrations used for the model development. Another limitation is that the simulations only account for the bedaquiline and delamanid coadministration, as they reproduce the DELIBERATE trial design where there was no coadministration of drugs with drug–drug interaction with bedaquiline or delamanid (e.g., antiretroviral drugs such as lopinavir/ritonavir), and no coadministration of drugs with QT liability (e.g., clofazimine). So, to fully assess the drug‐induced QTcF prolongation in new clinical trials with once‐daily regimens, other drug–drug interactions should be accounted for in the simulations.

In conclusion, we established the PK‐QTcF relationship of the combination of bedaquiline and delamanid. The metabolites, M2 and DM‐6705, were found to be responsible for the drug‐induced QTcF prolongation, and the pharmacodynamic interaction was best described by a competitive agonist model, indicating no synergic toxicity for the combination of bedaquiline and delamanid. Predictions from the model support the use of these drugs together in once‐daily regimens currently being tested in clinical trials, as the predicted QTcF prolongation is comparable to what is seen with currently used dosing regimens.

FUNDING

This analysis was supported by AIDS Clinical Trials Group (ACTG), National Institutes of Health (NIH).

CONFLICT OF INTEREST

Study bedaquiline and delamanid provided to NIH by Janssen and Otsuka, respectively, for the parent study. All other authors declared no competing interests for this work.

AUTHOR CONTRIBUTIONS

All authors wrote the manuscript. L.T., M.O.K., S.L.R., Y.S.C., J.S., C.M.U., A.H.D., G.M., K.E.D., and E.M.S. designed the research. L.T. performed the research. L.T., M.O.K., J.M., and E.M.S. analyzed the data.

Figure S1

Click here for additional data file.

Figure S2

Click here for additional data file.

Table S1

Click here for additional data file.

Table S2

Click here for additional data file.

Table S3

Click here for additional data file.